基于簡布法與強度折減法的滑坡穩定性研究

李青山，康小兵，許 模，李子隆，李 強

(成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059)

?

基于簡布法與強度折減法的滑坡穩定性研究

李青山，康小兵，許 模，李子隆，李 強

(成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059)

基于簡布法與強度折減法，對四川巴中黃果樹滑坡進行了穩定性研究。結果表明，簡布法所得的穩定系數隨著蓄水水位的升高而減小，在蓄水水位以1m/d的速度升降工況中，穩定系數隨著水位的上升和下降先減小后增大，如果其最小值達到了坡體滑動條件，坡體即滑動，如果未達到滑動條件則漸漸恢復到一定值；而強度折減法所計算出的穩定系數在工況1、2下比簡布法計算的結果稍大，在工況3、4下穩定系數誤差在中間時步較大，在其他時步相對較小。2種方法的分析結果穩定系數差異較小，都符合分析要求，均可為滑坡穩定性提供有效分析與評價。

穩定系數；簡布法；強度折減法；黃果樹滑坡

0 引 言

近年來，分析與評價邊坡穩定性的各種方法得到了廣泛的應用[1]，如簡布法、有限差分法、有限單元法、強度折減法等。徐衛亞等[2]分別對考慮加固力作用的邊坡穩定性分析的Morgens-Price和Sarma法推導了穩定系數的計算過程；連鎮營等[3]分析了強度折減有限元法在開挖邊坡的穩定性；楊有成等[4]闡明了在邊坡穩定分析中邊界范圍和網格的疏密對強度折減法結果的影響；鄭穎人，趙尚毅等[5]探討了有限元強度折減法研究進展。通常的強度折減法分析滑坡穩定性是以土體彈塑性有限元分析為基礎的，通過逐步折減土體的抗剪強度，改變土體的屈服準則，計算坡體的應力應變，考察土坡的位移發展情況，判斷土坡是否處于穩定狀態，并以破壞前土坡達到臨界狀態的折減系數為穩定系數，其穩定系數具有“強度儲備”的意義。眾多學者在強度折減計算及其工程應用方面進行了大量研究，如 Matsui[6]、Dawson[7]、趙尚毅[8]、唐春安[9]等進行了富有成效地探索，推動了強度折減理論在邊坡穩定性分析與評價方面的應用。

圖1 黃果樹滑坡工程地質平面

邊坡的實際失穩過程并非一次瞬間完成，而是由坡體局部應力集中或力學參數的減弱，導致坡體局部范圍出現損傷破壞，隨后破損范圍逐漸擴展，并最終形成整體滑面[10]。強度折減法的優點是能夠考慮坡體的應力-應變關系、流固耦合作用，分析坡體內一點的水平和豎直位移變化、應力變化。對比簡布法，強度折減法有以下優點：①能夠對具有復雜地質地貌的邊坡進行穩定性計算；②能夠考慮土體的應力-應變關系，假設條件無需太多簡化，形成了嚴密的理論體系；③能夠模擬出土坡失穩過程及滑動面形狀；④計算穩定系數時，無需假設滑動面位置及形狀，更不必進行條分計算；⑤可以計算滲流作用，得出孔隙水壓力，對于土體變形和滲流作用的耦合計算有很大的幫助[11]。本文采用簡布法和強度折減法，對四川巴中黃果樹滑坡進行計算分析，對比分析2種方法的準確性。

1 工程地質概況

巴中黃果樹滑坡位于恩陽河左岸，距水庫壩址約10.9 km。該滑坡為一古滑坡堆積體，總體上是一個經過多次滑移且滑坡體物質堆積差異性較大的古滑坡，表現為前緣和中部稍寬、后緣則逐步收斂，河流改道和巖層反翹，在平面上滑坡體整體上呈圈椅狀。通過現場勘查，結合滑坡地貌特征和變形特征，將滑坡分為3個區：古滑坡復活體(1區)、古滑坡坡體(2區)、覆蓋層變形體(3區)。1區縱長(近似垂直恩陽河方向)約390 m，橫寬(近似平行恩陽河方向)約280 m，面積約1.01×105m2，堆積體平均厚度約19 m，滑體方量約1.9×106m3；2區縱長約190 m，橫寬約350 m，面積4.8×104m2，覆蓋層平均厚度20 m，變形體方量約9.6×105m3；3區縱長約357 m，橫寬190 m，面積5.9×104m2，變形體平均厚度9 m，變形體方量約5.3×105m3。此次對該滑坡分析只針對1區和2區，3區地形起伏變化不大、覆蓋層較薄、大部分地區被殘坡積物覆蓋，且通過現場調查，該區未見明顯的地表變形跡象，故不做討論。黃果樹滑坡工程地質平面見圖1。C-C′工程地質縱剖面見圖2。

根據滑體及滑帶土體粒度成分特點，本次對滑坡體塊碎石土、滑帶土進行了原狀樣和室內重塑樣的土力學試驗，分別考慮在天然含水和飽水狀態下，對滑帶土、塊碎石土進行直剪試驗獲取巖土體強度參數。力學計算參數取值見表1。

2 穩定性對比分析

根據巴中黃果樹滑坡現場調查以及計算區域歷史發生過的大的滑動跡象和該水庫蓄水實際情況，選擇垂直于恩陽河方向的C-C′剖面進行分析。考慮了4種工況：蓄水水位353 m(工況1)、蓄水水位372 m(工況2)、水位從353 m以1 m/d的速度上漲到 372 m(工況3)、水位從372 m以1 m/d的速度下降到353 m(工況4)。運用簡布法、強度折減法分別計算了該滑坡的4種工況的穩定性，根據計算結果對該滑坡進行穩定性評價。

圖2 黃果樹滑坡C-C＇工程地質縱剖面

計算材料天然狀態抗剪強度飽和狀態抗剪強度粘聚力c/kPa內摩擦角φ/(°)粘聚力c/kPa內摩擦角φ/(°)天然密度ρd/kN·m-3飽和密度ρSat/kN·m-3變形模量E/kPa泊松比μ滲透系數K/10-6m·d-1滑體2118171518.52148000.321潛在滑帶10.57.587.152626.8460.30.1古滑帶91381127.2295000.331滑床15201217.523.22515000.20.1

2.1 簡布法

運用簡布法分別計算了該滑坡的4種工況的穩定性，工況1穩定系數為1.07；工況2為1.00；工況3為0.95；工況4為0.89。該滑坡在1、2工況時穩定系數分別為1.07和1.0，參照滑坡工程規范，該滑坡處于欠穩定狀態；而在工況3時，根據計算結果，大約在第2～3 d時，水位在355～356 m時，穩定性系數下降接近1，出現臨界狀態，大約在第7 d，水位升至360 m時，穩定性系數最低，為0.95，隨著計算的進行，穩定系數漸漸變大[12]；工況4時，根據計算結果，大約在第3 d初水位降至370 m時，穩定性系數下降接近1，出現臨界狀態，大約在第6 d，水位降至362 m時，穩定性系數最低，為0.89，隨著計算的進行，穩定系數漸漸變大。工況3、4穩定系數隨時間的變化見圖3。根據3、4工況的穩定系數，表明了該滑坡處于不穩定狀態。該方法計算的滑動面見圖4。

圖3 穩定系數隨時間的變化

圖4 簡布法計算的滑動面

2.2 強度折減法

用強度折減法對該滑坡的C-C′剖面進行穩定性分析。工況1穩定系數為1.11；工況2為1.01；工況3為0.97；工況4為0.92。與簡布法計算結果相比，該方法計算結果在353 m水位時較簡布法大，而水位上漲到372 m之后的計算結果比較接近，這也和實際工程狀況相吻合；工況3、4的穩定系數則與簡布法的計算不同，在穩定系數減小過程中，其速度變化平均值和簡布法保持一致，但最小值要比簡布法小，在增大過程中速度平均值變化同樣和水位下降過程中保持一致，說明強度折減法在滑坡體計算過程中穩定系數減小和恢復的速度和簡布法一致，且最大不平衡力曲線也更進一步證實了該算法的正確性，強度折減法計算的4種工況下的滑動面見圖5。

分析2種方法在不同工況下穩定系數值，工況1～4呈變下降趨勢，表明了該滑坡在水庫處于高水位時比低水位危險；水位下降過程的穩定系數值比水位上升時大，說明水位下降比水位上升時危險，這極有可能是由于地下水滲透力和揚壓力引起的[13]。

圖5 強度折減法計算的各工況下的滑動面

2.3 簡布法與強度折減法間的誤差分析

工況3、4中2種方法穩定系數的方差曲線見圖6。從圖6可知，在工況3中，2種方法穩定系數之間方差先增大后減小，說明兩者間的誤差也在中間時步最大，其他時步基本保持一致，這也表明了在地下水下滲和恢復穩定狀態其他時間段中，2種方法穩定系數計算結果是相近的；工況4變化規律和工況3的結果保持一致。經綜合分析，總體結果符合工程實際需要。

圖6 2種方法穩定系數方差曲線

2.4 位移、應力變化分析

在水位從353 m以1 m/d的速度上漲到372 m之后，又以1 m/d的速度下降到原初始水位353 m的過程中，設置3個水平位移監測點和3個豎直位移監測點。觀察監測點1、2、3的水平歷史位移曲線可知，1、2點的位移較大，分別為0.43、0.37 m，而3點的位移較小，為0.03 m。在相同位置的監測點4、5、6豎向位移分別為0.23、0.55、0.09 m，表明該滑坡中部和前緣變形較大。后緣水平和豎向位移較小的緣故是因為后面古滑坡體現在處于相對穩定狀態[14]，前緣和中部都是屬于古滑坡復活體。從XY方向剪應力云圖(見圖7)也可看出，前緣和中部應力最大，而后緣較小，處于古滑坡復活體的整體區域都處于剪應力增大區。經綜合分析，該復活體蓄水后處于不穩定狀態。相比簡布法而言，強度折減法分析滑坡在滿足穩定系數與簡布法一致時，更能分析具體點流固耦合作用下水平和豎向位移變化、剪應力變化，從宏觀到微觀監測到該滑坡的發展過程，而這也是傳統的簡布法所不能做到的，總體分析結果也符合要求。對比2種方法的分析結果，2種方法均可為滑坡穩定性提供有效分析與評價。

圖7 XY方向剪應力云圖

3 結 語

本文采用簡布法與強度折減法對巴中黃果樹滑坡4種工況進行穩定性分析，得如下結論：

(1)2種方法在不同工況下穩定系數的計算結果表明，工況1到工況4，該滑坡穩定性成下降趨勢，水庫處于高水位比低水位危險，水位下降比水位上升時危險。

(2)在水庫處于穩定水位時，強度折減法穩定系數較簡布法大；在水位升降過程中，兩者間的誤差在中間時步較大；在地下水下滲和恢復穩定狀態其他時間段中，計算結果相近。2種方法總體結果都符合工程實際需要。

(3)強度折減法分析滑坡在滿足穩定系數和簡布法一致時，更能分析具體點流固耦合作用下水平和豎向位移變化、應力變化，從宏觀到微觀監測到該滑坡的發展的過程。

(4)對比2種方法的分析結果，2種方法均可為滑坡穩定性提供有效分析與評價。

[1]王科， 王常明， 王彬， 等. 基于Morgenstern-Price法和強度折減法的邊坡穩定性對比分析[J]. 吉林大學報， 2013， 43(3)： 903- 905．

[2]徐衛亞， 周家文， 石崇. 極限平衡分析中加固力對巖質邊坡穩定性的影響[J]. 水利學報， 2007， 38(9)： 1056- 1064．

[3]連鎮營， 韓國城， 孔憲京. 強度折減有限元法研究開挖邊坡的穩定性[J]. 巖土工程學報， 2001， 23(4)： 403- 411．

[4]楊有成， 李群， 陳新澤， 等. 對強度折減法若干問題的討論[J]. 巖土力學， 2008， 29(4)： 1103- 1106．

[5]鄭穎人， 趙尚毅， 宋雅坤. 有限元強度折減法研究進展[J]. 后勤工程學院學報， 2005， 24(3)： 1- 6.

[6]MATSUI T， SAN K C. Finite element slope stability analysis by shear strength reduction technique[J]. Soils and Foundations， 1992， 32(1)： 59- 70．

[7]DAWSONE M， ROTHW H， DRESCHER A. Slope stability analysis by strength reduction[J]. Geotechnique， 1999， 49(6)： 835- 840．

[8]趙尚毅， 鄭穎人， 時衛民， 等. 用有限元強度折減法求邊坡穩定安全系數[J]. 巖土工程學報， 2002， 24(3)： 343- 346．

[9]唐春安， 李連崇， 李常文， 等. 巖土工程穩定性分析RFPA 強度折減法[J]. 巖石力學與工程學報， 2006， 25(8)： 1522- 1530．

[10]陳國慶， 黃潤秋， 周輝， 等. 邊坡漸進破壞的動態強度折減法研究[J]. 巖土力學， 2013， 34(4)： 1140- 1147．

[11]陳婷， 程圣國， 潘飛. 降雨模式對滑坡穩定性的影響[J]. 水力發電， 2014， 40(4)： 25- 27．

[12]張永輝， 郭松， 王永彬. 某滑坡成因機制分析及穩定性評價[J]. 水力發電， 2013， 39(1)： 30- 31．

[13]郭曉光， 黃潤秋， 鄧輝， 等. 平推式滑坡多級拉陷槽形成過程及成因機理分析[J]. 工程地質學報， 2013， 21(5)： 770- 779．

[14]程圣國， 陳高峰. 考慮參數相關性的滑坡穩定敏感性分析模式研究[J]. 水力發電， 2008， 34(6)： 38- 40．

(責任編輯 楊 健)

Research on Slope Stability Based on N.Janbu Law and Strength Reduction Method

LI Qingshan, KANG Xiaobing, XU Mo, LI Zilong, LI Qiang

(State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology,Chengdu 610059, Shichuan, China)

The stability of Huangguoshu Landslide is studied by using N.Janbu Law and strength reduction method. The results show that： (a) the N. Janbu law’s stability factor will decrease with the raising of water level; (b) when the water level is changed with the speed of 1m per day, the stability factor will decrease first and then increase; and (c) if stability factor reaches the minimum slope sliding condition, the slope will slide, and if it does not reach the sliding conditions, the stability factor will gradually recover to a certain situation. The results also show that： the safety factor resulted from strength reduction method under operation condition 1 and 2 is slightly larger than N. Janbu’s results, the safety factors of two methods under operation condition 3 and 4 have larger difference which is larger in middle time of operation condition change and smaller in other change time of operation condition. The analysis results’ difference of two methods are relative smaller and both methods can provide effective analysis and evaluation of landslide stability．

stability factor; N. Janbu Law; strength reduction method; Huangguoshu Landslide

2016- 03- 18

地質災害防治與地質環境保護國家重點實驗室自由探索課題資助項目(SKLGP2013Z003)

李青山(1990—)，男，重慶人，碩士研究生，研究方向為工程與環境水文地質；康小兵(通訊作者)．

TU457(271)

A

0559- 9342(2016)06- 0034- 05