某水庫工程邊坡穩定分析

帕爾哈提·艾則孜

（新疆維吾爾自治區水利水電勘測設計研究院，新疆 烏魯木齊 830000）

克拉瑪依市位于我國西北部邊陲，常年干旱缺水、人民群眾引水困難，水資源短缺嚴重制約了地區經濟發展。為了緩解當地人民群眾缺水現狀，保障地區工業發展需水和進一步促進當地社會經濟的發展，在準格爾沙漠北部，開工建設了一座以農業灌溉和工業供水為主、兼顧防洪的具有綜合開發任務的小（1）型水利樞紐工程。水庫建成后，將年新增工業及城鎮生活供水水量1611萬m3，年保障農業灌溉供水量3964.72萬m3，調蓄洪水可將水庫下游的防洪標準由目前的不足5年一遇提高到10年一遇，減少洪水對下游的危害，同時將有效滿足新增工業供水的需求，為新型工業化建設提供水資源保障。該水庫兩岸山體陡峻，開挖難度大，易產生滑坡坍塌等邊坡地質災害。本文針對工程不同部位邊坡的特性及破壞型，采用彈塑性有限單元法對具有代表性的大壩左岸邊坡開挖施工過程進行仿真模擬。

1 邊坡地質條件

水庫壩址區為地質條件為粉砂質粘土巖，山體相對平緩，岸坡一般25°～35°，材料分區從上到下依次為P1m1-2、P1m1-1、P1q4、P1q3、P1q2、P1q1、P1L、S2h1-2、S2h1-1、S1sh2、S1sh1、S1L、O2sh+b、O1m3，除 S1sh 及 O2sh+b 層以中硬的灰巖、砂巖為主外，其余各層主要為頁巖、粘土巖及泥質粉砂巖，巖石軟弱，易風化，其中，兩岸強、弱風化層厚5 m～17 m，河床強、弱風化層厚6 m～10 m。在P1m、P1q層中，據地質測繪與調查，兩岸地表長大裂隙共計135條。長度一般為10 m～30 m，大于30 m的占極少數，其發育方向主要有兩組，即 270°～300°，301°～330°組。多為陡傾角裂隙，溶蝕、卸荷張開，裂隙優勢走向與谷坡走向近平行，卸荷作用明顯。

2 邊坡穩定性分析

邊坡穩定分析方法較多，一般根據邊坡類型、滑移機制、邊坡變形與穩定控制要求等進行選擇。針對工程不同部位邊坡的特性及破壞型式，大壩區邊坡巖體上硬下軟，采用彈塑性有限單元法對水庫大壩左岸邊坡開挖施工過程進行仿真模擬，對邊坡的應力狀態、變形及屈服情況進行分析；同時，采用離散單元法對可能的失穩破壞模式進行分析，并對其整體穩定性及局部穩定性進行評價。

2.1 彈塑性有限單元法邊坡穩定分析

2.1.1 荷載

a.邊坡巖體自重；

b.開挖過程中的荷載：由開挖卸荷產生。

水庫工程區域構造應力較低，初始應力主要由自重產生，所以初始應力場按自重進行模擬。區域天然情況地下水較低，對邊坡穩定分析時未考慮地下水的影響。

2.1.2 巖體力學參數

計算中所用到的巖體力學參數采用現場工程地質勘查報告中的數據，卸荷帶巖體力學參數按卸荷帶外巖體的0.6倍進行計算[1-2]。巖體材料的本構關系采用Drucker-Prager準則。計算斷面內材料分區從上到下依次為P1m1-2、P1m1-1、P1q4、P1q3、P1q2、P1q1、P1L、S2h1-2、S2h1-1、S1sh2、S1sh1、S1L、O2sh+b、O1m3，各巖層靠近臨空面附近均存在一定范圍的卸荷帶。

因為斜坡表面附近相對松散的塊基本上在挖掘范圍之內，這些松散塊的局部穩定性在計算模型中不被考慮，而被認為是均勻連續體。整個計算與邊坡的整體穩定性有關[3]。

2.1.3 計算工況

選擇大壩左岸邊坡斷面進行計算，按照天然情況（自重應力場）及邊坡開挖，不考慮系統錨固措施兩個工況進行應力模擬分析。

2.1.4 計算結論

a.邊坡巖體的位移場主要是由開挖應力釋放和調整引起的。開挖完成后最大水平位移3.8 cm，位于464.0 m高程；最大鉛直位移3.1 cm，位于428.00 m高程。位移矢量見圖1。

圖1 開挖后位移矢量

b.開挖完成后位移矢見圖1，垂直開挖附近的主應力（第一主應力）方向接近于垂直開挖方向，主壓力應力方向改為拉伸應力方向，主要的拉伸應力為0到0.90 MPa。二次主應力（第三主應力）方向接近平行開挖面方向，基本處于壓應力狀態。在邊坡開挖臺階的拐角處，應力集中程度不同，地表附近的最大壓力應力值約為3.4 MPa。

c.在自重作用下邊坡在坡腳部位具有一定范圍的壓剪屈服區（見圖2），開挖后屈服區分布見圖3。與自重情況相比，在坡腳部位屈服區往深度方向有一定程度的發展，但各級坡均未現貫通屈服區，因此，邊坡整體是穩定的且有一定安全余度。

圖2 自重工況屈服區

圖3 開挖后屈服區

2.2 離散單元法失穩破壞模式分析

2.2.1 模型概化

根據巖體的統計參數、層間位錯和裂隙，將邊坡巖體劃分為塊體，構成離散單元的計算網格。塊通過切線、正常Springs、切線和普通阻尼器連接。塊的力的大小和方向是由塊的相對位置和運動決定的[4]。

考慮到實際裂縫分布具有一定的連接率，離散單元法將邊坡劃分為連接率為100%的單元[5]。因此，應考慮單元的結構面參數的實際連接率，并選擇巖體的加權平均值和結構面參數。考慮到巖體一旦產生滑動粘結強度，其內聚度取結構面本身的值，通過加權公式減小強度指數f：

式中：f為加權后的結構面強度指標；f1為巖體強度指標；f2為結構面強度指標；Φ為連通率。

按上式取75%的連通率根據巖體結構面抗剪強度進行概化得到計算中使用的巖體結構面抗剪強度參數。

裂隙的線剛度無實測值，根據經驗選取塊體之間彈簧的法向和切向線剛度均為10 GPa。

2.2.2 開挖邊坡安全系數定義

分析開挖邊坡穩定安全性的方法如下：模擬邊坡在自重及開挖卸荷作用下穩定后，按比例降低結構面抗剪強度，直到邊坡有比較大的滑坡或特別大的位移發生，據此來評價邊坡的穩定安全性。可以定義邊坡的穩定安全系數如下：

式中，fcr、ccr為結構面的臨界摩擦系數和凝聚力；f、c為結構面的初始摩擦系數和凝聚力。

離散單元法尚無可遵循的規范，其安全系數的定義不同于剛體極限平衡法，故其對邊坡的穩定性評價亦不能套用極限平衡法允許的抗滑穩定安全系數。

2.2.3 計算網格

據該斷面所揭示的斷層、層間錯動、層面及相對發育的節理裂隙進行網格劃分。圖4網格較密，表示節理裂隙十分發育；圖5網格稀疏近一倍，表示裂隙節理不十分發育，或邊坡經加固后巖體完整性得到提高。

圖4 離散元密網格

圖5 離散元稀網格

2.2.4 計算成果

a.當節理裂隙十分發育時，左岸邊坡的穩定安全系數不小于1.25。對于節理裂隙不十分發育，或邊坡經過支護加固處理巖體完整性得到提高后，左岸邊坡的穩定安全系數達到2.5以上。

b.開挖邊坡與自然邊坡的穩定安全系數基本相同，均不小于1.25。在邊坡結構面抗剪強度降低到臨界值以下時，下部軟巖的開挖引起坡腳滑動。因此，對上部已開挖邊坡必須支護完畢后才能對下部軟巖進行開挖。

c.在邊坡結構面抗剪強度降低到臨界值以下時，邊坡發生淺層滑動失穩，深層巖體仍保持穩定。失穩破壞形態見圖6和圖7。

圖6 密網格結構面強度降至初始值0.6倍破壞形態

圖7 稀網格結構面強度降至初始值0.4倍破壞形態

3 結論

彈塑性有限單元法對邊坡開挖施工過程進行仿真模擬，對邊坡的應力狀態、變形及屈服情況進行分析：邊坡巖體的位移場主要是由開挖引起的應力釋放和調整引起的。開挖完成后，開挖面附近的主應力方向接近豎直開挖面的方向，主應力的局部方向由主應力方向向主拉應力方向變化。亞主應力方向接近平行開挖面方向，基本處于壓應力狀態。在邊坡開挖臺階的拐角處，有不同程度的應力集中。在自重作用下邊坡在坡腳部位具有一定范圍的壓剪屈服區，與自重情況相比，在坡腳部位屈服區往深度方向有一定程度的發展，但各級坡均未現貫通屈服區，因此，邊坡整體是穩定的且有一定安全余度。

采用離散單元法對可能的失穩破壞模式進行分析，并對其整體穩定性及局部穩定性進行評價：當節理裂隙十分發育時，左岸邊坡的穩定安全系數不小于1.25。開挖邊坡和自然邊坡的穩定安全系數相當，左岸自然邊坡的穩定安全系數不小于1.25。在邊坡結構面抗剪強度降低到臨界值以下時，下部軟巖的開挖引起坡腳滑動。因此，對上部已開挖邊坡必須支護完畢后才能對下部軟巖進行開挖。在邊坡結構面抗剪強度降低到臨界值以下時，邊坡發生淺層滑動失穩，深層巖體仍保持穩定。