武子桓

（遼寧西北發(fā)電有限責(zé)任公司，遼寧沈陽110116）

1 工程概況

三灣電站是遼寧省十三五重點水利工程項目，是一座以防洪和發(fā)電為主，兼有其他多種功能的綜合性水利工程。電站壩址位于丹東鳳城市東湯鎮(zhèn)三灣村境內(nèi)的愛河干流上[1]。電站工程主要由大壩、溢洪道、導(dǎo)流泄洪洞及發(fā)電引水工程組成。其中，電站的大壩為漿砌石拱壩壩型設(shè)計，最大壩高66.00 m，壩軸線長231.10 m；大壩上游水庫的設(shè)計庫容為 2 430 萬 m3，設(shè)計裝配 3×8 MW 貫流式水力發(fā)電機組，裝機總?cè)萘繛?4 MW。三灣水電站的導(dǎo)流泄洪洞主要承擔(dān)施工期間的導(dǎo)流和電站運行期間的泄洪任務(wù)[2]。導(dǎo)流泄洪洞位于大壩右岸，上部是典型的低山丘陵地貌，同時相對埋深僅為15 m 左右，且?guī)r體較為松散，其中存在較多的裂隙和斷裂破碎帶，較小的埋深和復(fù)雜的地質(zhì)條件給開挖施工工作造成了諸多困擾。穿越卵礫石層的地下洞室工程條件復(fù)雜，施工難度較大，對其施工工法進行研究具有重要的工程實踐價值。為了保證施工的順利進行，首先需要結(jié)合工程設(shè)計和地質(zhì)環(huán)境特點，選擇合適的施工工法。基于此，工程項目部初步擬定了三臺階法和七步開挖法，通過對比和優(yōu)選確定最佳的施工工法[4]。

2 計算模型的構(gòu)建

2.1 軟件選擇

Midas GTS NX 2019 是一款針對巖土領(lǐng)域研發(fā)的通用有限元分析軟件，作為一個全面的有限元分析軟件包，可處理各種巖土工程設(shè)計應(yīng)用，支持靜力分析、動力分析、滲流分析、應(yīng)力-滲流耦合分析、固結(jié)分析、施工階段分析、邊坡穩(wěn)定分析等多種分析類型，適用于各種實際工程的準(zhǔn)確建模與分析，并提供了多種專業(yè)化建模助手和數(shù)據(jù)庫[5]。同時，該軟件作為新一代巖土分析軟件，應(yīng)用了最尖端的圖形和分析技術(shù)，可以支持最新的OS 圖形用戶界面。直觀的界面將使新用戶能夠輕松地將軟件集成到他們的工作流程中。此外，該軟件還集成了混合網(wǎng)格生成功能，該功能可創(chuàng)建使用六面體和四面體元素最佳組合的網(wǎng)格集，在建模復(fù)雜幾何形狀較尖銳的曲線和拐角時更有效[6]。鑒于該軟件的優(yōu)勢和此次研究的實際需求，利用Midas GTS NX 2019 進行研究對象的有限元模型構(gòu)建，對不同施工工法下的典型部位位移進行模擬計算，以獲取最佳施工方案。

2.2 計算模型

根據(jù)相關(guān)研究成果，地下洞室工程在開挖完畢后，距離開挖斷面中心點5～7 倍開挖直徑的范圍內(nèi)存在圍巖應(yīng)力應(yīng)變的影響[7]。為了保證結(jié)果的準(zhǔn)確性，研究中的計算模型邊界取10 倍洞徑，模型的長度為50 m，最終確定模型的尺寸為80 m×50 m×80 m。計算模型以隧洞中線指向下游的方向為Y軸正方向，以垂直于Y軸指向右側(cè)的方向為X軸正方向，以豎直向上的方向為Z軸正方向。對模型的底部施加全位移約束，模型的左右兩側(cè)施加X方向位移約束，模型的前后施加Y方向位移約束，模型的上部為自由邊界條件[8]。對模型進行六面體八節(jié)點實體單元劃分。

2.3 參數(shù)選取

為了簡化模型，提高模型的計算效率，對研究洞段的圍巖采用摩爾-庫侖本構(gòu)模型，支護結(jié)構(gòu)混凝土按線彈性本構(gòu)模型處理；初支結(jié)構(gòu)和二次襯砌的厚度分別為27 cm 和45 cm，均采用板單元模擬；系統(tǒng)錨桿直徑22 mm，長4.0 m，采用植入式桁架單元定義；超前小導(dǎo)管的加固循環(huán)為6.0 m，在預(yù)加固完成之后，圍巖的粘聚力、內(nèi)摩擦角及彈性模量可以分別提升90%，50%和95%左右。結(jié)合上述分析及工程資料，確定如表1 所示的模型物理力學(xué)參數(shù)。

表1 模型物理力學(xué)參數(shù)

3 計算結(jié)果與分析

利用已構(gòu)建的計算模型，對采用三臺階開挖施工工法下的隧洞位移進行模擬計算，從結(jié)果中提取拱頂沉降和周邊收斂變形量進行統(tǒng)計計算，結(jié)果見表2。由表2 可以看出，從拱頂沉降占總位移的比例大小來看，上臺階最大、中臺階次之、下臺階最小；從周邊收斂占總位移的比例來看，其排列順序分別為中臺階、上臺階和下臺階。從位移計算結(jié)果來看，如果對研究洞段采用三臺階法開挖，下臺階開挖過程中的位移變形量相對較小，而上臺階和中臺階開挖過程中的位移變形量相對較大。

表2 三臺階法位移變形計算結(jié)果

利用已構(gòu)建的計算模型對采用七步法開挖施工工法下的隧洞位移進行模擬計算，從結(jié)果中提取拱頂沉降和周邊收斂變形量并進行統(tǒng)計計算，結(jié)果見表3。由表3 可以看出，上臺階開挖過程中的拱頂沉降量占比最大，其次是中臺階核心土，第三次是下臺階核心土；上臺階開挖過程中的周邊收斂變形量占比也最大，中臺階左，第三次中臺階右。總體來看，上臺階與核心土在開挖過程中會導(dǎo)致較大的位移變形。

表3 七步法位移變形計算結(jié)果

對比表2 和表3 可以看出，三臺階開挖法的拱頂沉降位移量的最大值為40.25 mm，七步法的最大值為35.22 mm；七步法的位移變形量相對較小，具有較為明顯的優(yōu)勢，推薦在工程施工中采用。

在確定采用七步法進行進口段開挖施工之后，需確定合理的開挖步長，為施工的順利進行提供必要的支持。研究中設(shè)計了0.5，1.0，2.0，3.0 m等4 種不同的開挖進尺，利用構(gòu)建的數(shù)值模型進行位移變形模擬計算。根據(jù)計算結(jié)果，繪制出拱頂沉降變形量和周邊收斂變形量隨開挖進尺的變化曲線，如圖1 所示。從圖1 看，隨著開挖進尺的增大，拱頂沉降和周邊收斂變形量均呈現(xiàn)出不斷增大的變化特征。同時，當(dāng)開挖進尺較?。ㄐ∮?.0 m）時，拱頂沉降和周邊收斂變形量的增加幅度相對較小，當(dāng)開挖進尺較大（大于1.0 m）時，拱頂沉降和周邊收斂變形量的增加幅度開始迅速增大。結(jié)合計算結(jié)果，同時考慮工程施工的經(jīng)濟性，建議采取1.0 m 的開挖進尺。

圖1 位移變形隨開挖進尺變化曲線

4 結(jié)論

此次研究是利用數(shù)值模擬的方法對遼寧省三灣水電站導(dǎo)流泄洪洞進口的開挖施工工法和施工參數(shù)進行比選和優(yōu)化研究，得出：1）施工如果采用三臺階法，拱頂沉降和周邊收斂變形主要出現(xiàn)在上臺階和中臺階開挖過程中；如果采用七步法，拱頂沉降和周邊收斂變形主要出現(xiàn)在上臺階和核心土的開挖過程中；2）兩種開挖施工工法位移量的對比顯示，七步法的位移變形量較小，建議在工程施工中選用；3）同時根據(jù)模擬計算結(jié)果，考慮工程施工成本，施工中應(yīng)選擇1.0 m 的開挖進尺。