樹坪滑坡切坡抗滑增穩優化治理數值模擬研究

吳麗麗，賀可強，全恩厚，全恩偉，朱裕祥，吳遠鋒

（1.青島理工大學 土木工程學院，山東 青島 266033；2.山東華通路橋工程有限公司，山東 臨沂 273400）

0 引 言

在地形地貌的演變過程中，滑坡作為一種典型的塊體運動形式，是一種重大的自然地質災害，給人們的生產、生活帶來嚴重威脅［1］。堆積層滑坡是滑坡的一種基本類型，具體指發生在第四系及近代松散堆積層的一類滑坡，該類滑坡規模大、突發性強，容易受多種因素影響，且滑移條件復雜，在我國滑坡分布中占比較大［2］，尤其是我國的西南、西北地區和長江三峽庫區。堆積層滑坡已對我國國民經濟和人民生活等方面造成了極大危害，滑坡的防治問題得到越來越多的重視。

滑坡治理工程中應考慮的影響因素主要有滑坡類型、滑坡體積大小、穩定性以及滑坡區域的地質構造條件、構（建）筑物類型及其分布范圍、施工實施用具和季節性等，需要采用多種治理形式進行綜合治理，如圍截排水、卸載壓腳、打入抗滑樁、高壓注漿、錨固預應力錨索等［3］。對于涉水條件下大型堆積層滑坡的治理，由于滑坡體積大、滑體厚，主要采用切坡卸載和相關工程措施（堆載、排水、支護等）綜合治理方法［4-9］。夏艷華等［10］在切坡卸載設計中采用傳遞系數法思想，引入剩余抗滑力概念，利用預期安全系數逆向計算各個條塊的剩余抗滑力，給出切坡優化方案，確定最優條件下的切坡位置和規模；高俊等［11］針對某滑坡治理時切坡壓腳土方量不平衡問題進行研究，為滑坡體設計出高壓腳少卸荷、低壓腳多卸荷兩種治理方案并進行對比，發現采用低壓腳多卸荷治理方案效果更好，且當堆載體積達到一定程度后，滑坡體的抗滑安全系數將保持不變?；谂R界壓腳體進行滑坡治理可解決切坡土方量與堆載土方量的不平衡現象，降低工程造價，為實際工程建設提供一定的參考價值。張帆宇等［12］采用物理模型試驗及工程地質分析的定性分析與FLAC3D數值模擬定量分析相結合的方法，全面分析了袁家灣高速公路修建過程中因開挖高陡切坡引發的順層巖質滑坡失穩過程和演化機理，發現人工改造邊坡會不同程度地影響邊坡地質環境，并導致其他條件對邊坡穩定性影響程度也發生改變，即具有伴生效應；周中等［13］選取貴州晴隆典型堆積層邊坡進行了現場切坡開挖試驗，并同時進行原位監測，結果表明，在切坡開挖的影響下堆積層邊坡偏于發生淺層牽引式破壞，且滑動變形范圍在坡面的0～4 m，坡面變形最大，從坡面至坡面以下變形逐漸變小。

本文基于上述研究成果和堆積層邊坡切坡卸載過程中的下滑動力變化與變形響應特點，選用MIDAS GTS NX有限元軟件及其內置的強度折減法對樹坪滑坡A-A剖面優化治理方案的抗滑增穩效應進行分析，以期為該滑坡的治理提供參考。

1 切坡卸載對堆積層邊坡穩定性的 影響

1.1 切坡卸載條件下堆積層邊坡下滑動力及變形響應演化特征與規律

堆積層邊坡的切坡卸載過程是一個穩定的單向增加過程，即隨著切坡進行，邊坡將由不穩定狀態（或欠穩定狀態）逐漸進入穩定狀態。為了闡明切坡卸載過程對堆積層邊坡穩定性狀態的影響，設堆積層邊坡為均質體，且該均質體呈各向同性，坡體厚度也呈均一變化。將邊坡潛在滑動條塊作為研究對象，堆積層邊坡滑動條塊受力情況如圖1所示。圖1中，Gi為滑動條塊重力；Ri為滑面上的抗滑力；Ni為滑動條塊有效應力；ΔPi為下滑動力變化量；Hi為滑動條塊高度；hi為切坡高度；θi為滑面傾角。

由圖1可知，切坡前滑動條塊的下滑力Ta、抗滑力Ra分別為

圖1 堆積層邊坡滑動條塊受力示意圖Fig.1 Stress diagram of sliding strip of accumulation slope

切坡后滑動條塊的下滑力Tb、抗滑力Rb為

式中：c為滑面黏聚力；φ為滑面內摩擦角；θ為滑面傾角；l為滑動條塊底面長度；γ為滑動條塊重度。

因此，切坡引起的下滑動力變化量為：

堆積層邊坡切坡卸載時坡體變形主要受切坡影響，宏觀上表現為卸荷回彈變形，變形發展主要經歷變形趨緩階段和變形穩定階段［14］。

變形趨緩階段：隨著切坡卸載進行，坡體變形逐漸得到控制，并表現為趨緩變形，坡體穩定性得以提高。

變形穩定階段：隨著切坡結束，坡體變形逐漸停止并趨于穩定。

另外，依據堆積層邊坡自身特點和構造特征，將邊坡坡體假設為理想彈塑性。對邊坡坡體進行計算時，其斷面受力與變形物理參數均取平均值［15］。依據彈塑性力學基本原理，堆積層邊坡切坡卸載引起的下滑動力與位移變化量關系式為

根據式（6）～（8）可知

式中：ΔSi、ε-i、σ-i分別為切坡卸載引起的滑動條塊位移變化量、應變均值、應力均值；E-i為切坡卸載時相應的滑坡體變形模量均值；Vi為滑動條塊體積（單位寬度條件下），Vi＝Hil cosθ。

由式（5）～（9）可知，在其他因素相對不變的情況下，堆積層邊坡切坡引起的邊坡下滑動力變化量與位移變化量主要取決于切坡過程中切坡高度的變化。隨著切坡過程中切坡高度增加，滑坡的下滑動力降低，抗滑力增加，滑坡位移變化量減小?？梢哉J為切坡卸載過程是堆積層邊坡下滑動力變化與位移變化的主要動力影響因素。

1.2 切坡卸載條件下堆積層邊坡穩定性演化規律

根據式（3）～（4）和邊坡穩定性極限平衡原理，確定切坡后邊坡的穩定性系數Fsi為

由式（10）可知，在切坡卸載過程中，只有切坡高度變化量Hi-hi是變化參數，因此邊坡穩定性系數Fsi與切坡卸載動態變化過程直接相關。

2 實例分析

2.1 樹坪滑坡工程概況

樹坪滑坡［16］與三峽大壩相距約47 km，位于長江南岸的百福坪背斜東翼，滑坡整體由南至北展布，地形起伏狀態為南高北低，且陡、緩相間?；麦w縱向（南北）長約800 m，橫向（東北）寬約700 m，面積約55×104m2，厚度約30～70 m，厚度均值約50 m，體積約2 750×104m3?；聳|側和中部變形較大，且為主滑區，體積約為1 575×104m3?；麦w發育于沙鎮溪背斜近軸部，屬于逆向坡，其地層產狀120°～173°，傾角9°～38°，該滑坡斜坡走向總體上與長江平行。從樹坪滑坡的物質組成上分析，其滑體主要組成物質為崩坡積碎塊石土，滑坡下部以粉質黏土為主；滑帶以碎石角礫的粉質黏土為主，其中角礫大部分為砂石、泥巖和泥灰巖?；孪路鶐r是三疊系中統巴東組地層，以紫紅、灰綠色的中厚層狀粉砂巖夾泥巖，以及灰、淺灰色的中厚層狀泥巖為主要成分，剖面圖如圖2所示。樹坪滑坡因與三峽庫區較近，屬于涉水、大型老滑坡，該滑坡區域的地表水系來自長江，水位145～175 m，主要受三峽水庫的蓄水控制?；伦冃问加?003年三峽水庫蓄水，因常年水位不斷升降，再加上降雨影響，導致樹坪滑坡產生了較大的累計變形。

圖2 樹坪滑坡工程地質剖面圖Fig.2 Engineering geological sectional drawing of Shuping landslide

為了有效降低滑坡變形速率，2014年對樹坪滑坡進行應急治理［17］，應急治理工程的設計工況為：自重+庫水位175～145 m，整體安全系數1.03，上層滑體、潛在滑體的安全系數均需達到1.10；校核工況為：自重+庫水位175～145 m+50 a一遇暴雨，整體安全系數需達到1.0，上層滑體與潛在滑體的安全系數則需達到1.05。該應急治理工程的總體方案為：削方+壓腳+地表排水溝。主要工程量為：每隔20 m高差設置一級馬道，削方工程土石方開挖體積575 774 m3，壓腳工程土石方回填體積373 962 m3，另外設計5條排水溝，總長2 850 m。同時根據滑坡規模和危害程度確定治理工程等級為Ⅰ級?；缕矫嫘螒B特征見圖3。

根據《滑坡防治工程勘察規范》（GB／T 32864-2016）要求，滑坡防治工程的穩定狀態評價標準為：F＜1.00時，滑坡屬于不穩定狀態；1.00≤F＜1.05時，滑坡屬于欠穩定狀態；1.05≤F＜1.15時，滑坡屬于基本穩定狀態；F≥1.15時，滑坡屬于穩定狀態。削坡減載工程中，堆積體或土質邊坡高度超過10 m時，應進行分級放坡。根據規范中對滑坡穩定狀態的評價標準及施工技術的要求可知，有必要對樹坪滑坡提出優化治理方案。

2.2 樹坪滑坡優化治理方案選取與數值建模

2.2.1 優化治理方案的選取

本文以減小滑坡后緣下滑力、增加滑坡前緣抗滑力為目的，根據樹坪滑坡的地形地貌條件以及《滑坡防治工程設計與施工技術規范》（DZ／T 0219—2006）中的要求，同時綜合考慮多種條件，如滑坡地區地質、施工季節條件等，選取切坡、堆載作為優化治理手段，堆載作為安全儲備可提高穩定性。該滑坡具體優化治理方案為：切坡位置為中部（滑移面呈凸出狀），堆載位置為坡腳，切坡坡形為階梯狀；根據《邊坡工程處治技術》［18］中堆積體邊坡坡率參考值并結合相關工程經驗，該滑坡切坡參數為：切坡自上而下分6級（層），單級坡高10 m，坡率為1∶1.75（第1～3層）、1∶2.25（第4～6層），平臺寬5～8 m。樹坪滑坡切坡堆載示意圖如圖4所示。

2.2.2 計算模型的建立及參數選取

采用巖土有限元軟件MIDAS GTS NX進行數值模擬，考慮巖體材料的非線性特征，采用強度折減法進行計算分析，可為巖土工程中的實際施工提供可靠的技術支持。以樹坪滑坡典型剖面A-A為原型建立坡體，長度880 m，前后緣高程分別為20.3 m、418 m的計算模型，概化模型如圖5所示，共劃分5 939個節點，5 832個單元。根據樹坪滑坡的野外和室內試驗，通過工程類比，最后確定滑坡的物理力學參數如表1所示。

表1 樹坪滑坡物理力學參數Tab.1 Physical and mechanical parameters of Shuping landslide

圖5 樹坪滑坡A-A剖面數值模擬概化模型Fig.5 Generalized model of numerical simulation of Shuping landslide section A-A

2.2.3 計算工況

為了驗證優化治理方案的適用性，根據三峽庫區水位的調度方案和實測降雨量，并考慮樹坪滑坡動態卸荷、自重應力等因素的影響，設置4種計算工況，見表2。

表2 計算工況Tab.2 Calculation conditions

3 結果分析

3.1 優化治理方案的應力場分析

優化治理方案實施后，坡體的水平應力、豎向應力如圖6～7所示。其中，坡體表面的水平應力、豎向應力均表現為壓應力，其大小隨著切坡進行有所減小，這是切坡卸荷作用造成的，但整體沒有出現應力集中現象。最大剪應力云圖見圖8，坡體各臺階處沒有出現應力集中現象，表明切坡后邊坡不會發生剪切破壞。

圖6 邊坡切坡后水平應力Fig.6 Horizontal stress of slope after cutting

圖7 邊坡切坡后豎向應力Fig.7 Vertical stress of slope after cutting

圖8 邊坡切坡后最大剪應力Fig.8 Maximum shear stress of slope after cutting

3.2 優化治理方案的位移場分析

從宏觀上看，樹坪滑坡發生的橫向變形大于豎向變形，因此模擬結果中主要分析坡體的水平方向位移。由于模擬結果中第1～3切坡層的最大水平位移變化較?。t色區域范圍為零），第4～6切坡層的最大水平位移變化比較明顯，因此以第4～6層切坡的最大水平位移為例分析其變化規律。圖9～11所示分別為第4～6切坡層的水平位移，可以看出，邊坡臨空面附近的最大水平位移表現出擴大趨勢，且第4、5、6切坡層的最大水平位移紅色區域分別為9.2%、11.4%、11.3%，第4～5、5～6切坡層最大水平位移變化量分別為2.2%、0.1%，即最大水平位移的變化量呈現出降低趨勢?？梢?，該模擬演化位移特征與切坡卸載條件下堆積層邊坡變形減小的規律相符。

圖9 第4層切坡水平位移Fig.9 Horizontal displacement of cutting slope in the fourth layer

圖10 第5層切坡水平位移Fig.10 Horizontal displacement of cut slope of the fifth layer

圖11 第6層切坡水平位移Fig.11 Horizontal displacement of cut slope of the sixth layer

3.3 優化治理方案的塑性區及穩定性分析

切坡后坡體的有效塑性區見圖12。由圖12可以看出，切坡坡體表面塑性區不明顯，沒有擴展和貫通趨勢。通過強度折減法計算得出的穩定性系數見圖13，可以看出，隨著切坡進行，穩定性系數逐步升高，與前述位移場變化結果吻合。

圖12 邊坡切坡后有效塑性區Fig.12 Effective plastic zone of slope after cutting

圖13 邊坡穩定性系數隨切坡變化時序曲線Fig.13 Time series curve of slope stability coefficient changing with cut slope

3.4 不同計算工況下穩定性模擬結果分析

通過MIDASGTSNX有限元軟件內置的強度折減法對邊坡的穩定性進行計算分析，得出各計算工況下應急治理工程與本文優化治理方案的整體穩定性系數，見表3。

表3 各計算工況穩定性系數與切方量Tab.3 Stability coefficient and cutting volume of each calculation condition

從表3可以看出，工況2、3條件下應急治理工程穩定性系數滿足要求（F≥1.15）；在175 m水位和實測工況條件下，穩定性系數下降，分別處于欠穩定（1.00≤F＜1.05）和基本穩定狀態（1.05≤F＜1.15）。本文提出的優化治理方案在各計算工況下穩定性系數均滿足要求（F≥1.15），這與堆積層邊坡在切坡卸載條件下坡體穩定性增加的規律相符。由前文可知，應急治理工程的切方量為575 774 m3，優化治理方案的切方量為531 647 m3，與應急治理工程相比切方量降低了8%。由上述結果可以看出，本文提出的優化治理方案對樹坪滑坡起到了良好的抗滑增穩效果，且具有較好的適用性。

4 結 論

（1）堆積層邊坡的切坡卸載過程是一個單向增加穩定性的過程，該類邊坡的抗滑增穩效應表現為：下滑動力減小，變形響應降低。其中，下滑動力變化量主要受切坡高度、邊坡傾角、坡體物理力學性質等因素影響，切坡過程中切坡高度變化是堆積層滑坡下滑動力最明顯的動態變化因素，說明切坡卸載過程是影響該類邊坡穩定性的主要因素。

（2）隨著樹坪滑坡切坡進行，坡體表面的水平應力、豎向應力均表現為壓應力，壓應力值因切坡卸載作用減小，切坡完成后的最大剪應力沒有出現應力集中現象，不會發生剪切破壞；隨著切坡進行，位移變化量開始降低，穩定性系數逐步升高，體現了切坡卸載對樹坪滑坡坡體產生的抗滑增穩效應。

（3）為驗證樹坪滑坡優化治理方案的適用性，對4種計算工況下的優化治理方案與應急治理工程穩定性進行數值模擬分析，由模擬結果可知，各計算工況下優化治理方案的穩定性能滿足要求，切坡方量也比應急治理工程降低了8%。