楊發(fā)棟

（中國電建集團華東勘測設計研究院有限公司，浙江 杭州 311122）

水利樞紐工程中，順向高邊坡及偏壓隧洞的開挖、支護，是根據(jù)相關地質資料，利用數(shù)字模擬軟件以及建模分析的方式，加載工程各項參數(shù)，呈現(xiàn)工程施工的發(fā)展趨勢，為水利工程的施工、管理提供基礎數(shù)據(jù)[1]。在大型順向高邊坡、偏壓隧洞的開挖施工中，缺乏偏壓特性、地質屬性以及施工的時效性等研究資料時，應在現(xiàn)有相關的資料基礎上，結合實際的施工經驗，建立有限元模型，進行工程巖土材料的時效性分析。

1 工程概況

某水庫為大型低壩水利樞紐工程，水庫主要以農業(yè)灌溉為主，兼顧發(fā)電、航運及防洪；工程區(qū)位于河道夾角的峽谷處，山勢陡峭，水庫總庫容為15.1 億m3，校核庫容為1.05 億m3，正常蓄水位為313 m；1 號泄洪洞出口段為一級邊坡，施工區(qū)域地質斷裂帶較多，1 號泄洪洞修筑開挖尺寸15 m×15 m，出口段與放空洞的凈距為19.1 m。工程區(qū)域巖層為泥盆系下統(tǒng)D1，整體為單斜巖層構造，其中1 號泄洪洞出口段320°∠26°，隧洞修筑線路與巖層的走向夾角為50°～70°。工程區(qū)域有三個與隧洞斜切的斷層，均為東北走向，破碎帶的寬度在1.0～1.5 m，類型為碎裂巖系的斷層角礫巖、未固結或弱固結的泥狀巖石等，其巖石結構極為松散，具有透水性強等特點；工程區(qū)域地下水含量較大，多為軟塑粘性土。

1 號泄洪洞出口段，主要為強風化泥質粉砂巖，右側邊坡深度在50 m 以上，巖層變化情況極為復雜，部分區(qū)域呈現(xiàn)較大的傾斜角度，不具備施工條件，1 號泄洪洞出口段工程地質情況如圖1 所示。

圖1 1 號泄洪洞出口段工程地質情況

根據(jù)地質勘探的結果顯示，本次工程區(qū)域的地質斷層主要為F6，F(xiàn)244，F(xiàn)243 等，此區(qū)域裂隙破碎帶軟弱夾層較多，屬于中等透水性能，地下水折減系數(shù)0.5～0.7。其中，D16主要由紫紅色厚石英砂巖層、砂巖夾泥巖、含礫石英砂巖、細砂巖和粉砂巖構成；D15主要由雜色厚層狀含礫粗粒石英砂巖與紫紅色、石英巖狀砂巖夾薄層構成。

2 模型分析

2.1 計算方法

本次研究采用連續(xù)介質力學計算模型以及彈塑性力學作為本次研究的基礎，對巖土材料變形程度進行分析，使用C-M 準則以及流動法則進行巖土材料的粘彈性應力、形變，以及工程開挖支護的建模分析。

C-M 準則對巖土材料的粘彈性應力表示為：

式中：τf——巖土材料的抗剪強度值；

c——巖土材料的粘聚力；

σn——巖土材料的正應力值；

φ——巖土材料內摩擦角度。

在忽略巖土材料的強化作用下，假設發(fā)生屈服產生無法變形時，開挖釋放荷載所產生的應力變形，均可作為粘性、非粘性變形，計算公式：

式中：δ——全部荷載釋放后形成的位移；

δuv——開挖形成的瞬間位移（非粘性變形）；

δv——開挖過程中因持續(xù)變形而產生的位移（粘性變形）；

F——巖土材料所有荷載；

Fuv——非粘性變形的荷載值；

Fv——粘性變形的荷載值；

α——荷載瞬間釋放的因子。

釋放荷載的計算過程，是將粘性、彈塑性荷載作為節(jié)點力，施加于開挖節(jié)點，進行迭代求解，計算出開挖瞬間變形的非粘性彈塑性荷載值；消除施加在開挖節(jié)點上的應力后，添加支護參數(shù)，求出巖土材料在支護作用下產生的粘性變形增量。

2.2 計算模型

根據(jù)工程地質條件及施工方案，建立施工區(qū)域的地形、巖層結構、邊坡、1 號泄洪洞以及放空洞等模型，網(wǎng)格模式如圖2 所示。

圖2 網(wǎng)格模式圖

2.3 計算參數(shù)

根據(jù)本工程的地質參數(shù)及施工經驗，完成數(shù)值模擬計算。由于彈塑性作用受地下水的影響較大，以往同類工程地下水的折減系數(shù)取0.7～0.9，參考到本次工程的地質參數(shù)，最終確定本次地下水折減系數(shù)為0.5～0.7[2]，本工程各基巖物理學參數(shù)如表1 所示。

表1 基巖物理學參數(shù)

3 研究結果分析

工程相鄰放空洞的貫通區(qū)，未發(fā)生明顯塑性屈服區(qū)。假設巖土材料為理想的彈塑性體，在施工中瞬間非粘性變形荷載將得到釋放，而粘性變形的荷載則是隨著施工的過程逐漸釋放。荷載瞬間釋放因子，根據(jù)相關研究結果及同類工程經驗，取0.8。

3.1 開挖施工中邊坡應力分析

結合基巖物理學參數(shù)，對邊坡初始應力進行模擬分布，發(fā)現(xiàn)工程邊坡中的應力值較低，模擬計算中的最大數(shù)值為2.8 MPa。以計算出的初始應力值為基礎，使用邊坡有限元抗剪強度折減法進行結構的穩(wěn)定性能評估，獲得邊坡的安全系數(shù)為1.58。根據(jù)邊坡剪切應變增量圖顯示，應力分布情況受地質斷層的顯著影響，其中較大的潛在滑動面處于施工區(qū)域高程的260 m與邊坡的底部。

在邊坡開挖施工中所發(fā)生的重力分布轉移，將會在表面引起拉應力，對施工區(qū)域局部的穩(wěn)定性造成波動。邊坡開挖施工時，臨近洞口部高程為230 m 的區(qū)域出現(xiàn)了大幅度的變形。因此，對該區(qū)域布置監(jiān)控點，以獲得地層位移的變化情況。放空洞、1 號泄洪洞監(jiān)測點位分布如圖3 所示。

模型中監(jiān)控點的監(jiān)測結果顯示，該區(qū)域中最大的位移變化為14.5 mm，施工中產生的荷載釋放，將導致邊坡前表層位移及塌陷，整體變形程度較為顯著，因此需要對邊坡進行有效的支護加固處理。監(jiān)控點模擬位移增量如圖4[3]所示。

3.2 偏壓隧洞開挖支護措施

由于1 號泄洪洞的地質條件較為復雜，因此使用管棚預支護進行處理，以最大限度地優(yōu)化支護段的支護強度，在開挖施工后對偏壓出口發(fā)生位移段進行監(jiān)測。根據(jù)監(jiān)測結果發(fā)現(xiàn)，各個監(jiān)控樁位的位移均朝向泄洪洞洞口方向，在使用管棚預支護后，發(fā)生的最大位移距離為4.9 mm。不同樁號的最大位移距離如表2所示。

圖3 放空洞、1 號泄洪洞監(jiān)測點位分布圖

圖4 監(jiān)控點模擬位移增量

表2 不同樁號的最大位移距離 單位：mm

1 號泄洪洞管棚施工的部分措施為開挖后實施，K0+565—K0+502 之間的區(qū)域為五類圍巖，與斷層F6，F(xiàn)243 相切，最大的破碎帶寬度在1.3～5.0 m，并且由斷層角礫巖、未固結或弱固結的泥狀巖石構成，結構極為松散具有強透水性，穩(wěn)定狀態(tài)較差，巖體的風化程度較強，施工中常會出現(xiàn)破損、巖石脫落的現(xiàn)象，成洞難度較大。在使用管棚預支護施工后，仍然出現(xiàn)了最大4.9 mm 的位移，導致泄洪洞的管棚、襯砌發(fā)生局部開裂、洞室承載力下降的情況，由于對結構的安全產生了較大的影響，因此，需要對結構的發(fā)展狀況進行穩(wěn)定性的評估。

本工程中1 號泄洪洞兩側山體嚴重地不對稱，因此產生了較大的偏壓。為了詳細地分析偏壓對泄洪洞支護的影響，在泄洪洞出口段的邊墻及拱頂位置設置了三處監(jiān)控點位。發(fā)現(xiàn)此區(qū)域開挖時的瞬時荷載釋放所產生的位移增加較大，模擬計算的結果與工程的實際監(jiān)測資料的發(fā)展趨勢相一致。說明運用本研究方法，能夠有效地模擬工程施工過程情況，避免應用彈塑性力學理論進行工程開挖模擬時，瞬時彈性變形增量，造成后期加固支護措施未能發(fā)揮作用的情況[4]。出口段各監(jiān)測點位移增量如圖5 所示。

圖5 出口段各監(jiān)測點位移增量

根據(jù)上圖顯示，邊墻兩側監(jiān)控點位的最大位移增量為2.0 mm、0.2 mm，數(shù)值相差顯著，原因是山體偏壓效應致使管棚出現(xiàn)受力不均，根據(jù)進一步監(jiān)測發(fā)現(xiàn)，在山體較高側的管棚拱腳最大應力為6.1 MPa，較低側的管棚拱腳最大應力為1.0 MPa，監(jiān)控點斷面的應力情況均呈現(xiàn)內大外小，具有明顯的偏壓狀態(tài)。導致圍巖由內向外偏移，內側最大位移增量為2.0 mm，說明構造偏壓對洞室穩(wěn)定會造成較大的影響[5]。放空洞、1 號泄洪洞夾層的層狀巖體結構，在四類圍巖地層中能夠產生較大的發(fā)展。

4 結論

根據(jù)對工程區(qū)域在增量、應力及塑性區(qū)三個方面的研究發(fā)現(xiàn)，使用工程中制定的支護方案實施隧洞施工，將會造成失穩(wěn)的可能性，并且在地下水因素的作用下，將會增加失穩(wěn)的程度，因此，需要改善、強化當前的支護方案，避免隧洞邊壁、邊坡表面出現(xiàn)局部的倒塌、破壞。