基于考慮含水率對土體抗剪強度參數影響的邊坡穩定性分析

王存利 唐棟 李毅峰 蔣中明 王亞雄

摘? 要：降雨是誘發滑坡的關鍵因素，一般認為其致災機理主要在于導致土體重度增加、孔隙水壓力提高及土的抗剪強度降低。然而在水的作用下，土體抗剪強度參數還存在一定的弱化現象，這對邊坡的穩定性有一定的影響。該文以東嶺信滑坡為研究對象，基于對滑帶重塑土的直剪試驗，研究含水率對土體抗剪強度參數的影響，并基于GeoStudio軟件二次開發實現考慮含水率對土體抗剪強度參數影響的邊坡穩定性分析。結果表明，抗剪強度參數隨含水率的增加而逐漸減小，滑帶土的黏聚力和內摩擦角整體上均可與含水率擬合成線性關系。忽略水對巖土體抗剪強度參數的弱化作用會高估邊坡的穩定性。研究結果對于進一步認識降雨誘發滑坡機理具有一定理論意義。

關鍵詞：非飽和土；邊坡穩定；抗剪強度；含水率；數值模擬

中圖分類號：TU94+3.2? ? ? 文獻標志碼：A? ? ? ? ? 文章編號：2095-2945（2024）06-0061-04

Abstract： Rainfall is the key factor to induce landslides， and it is generally believed that the disaster mechanism mainly lies in the increase of soil weight， the increase of pore water pressure and the decrease of shear strength of soil. However， under the action of water， the shear strength parameters of soil still have a certain weakening phenomenon， which has a certain impact on the stability of the slope. Taking Donglingxin landslide as the research object， based on the direct shear test of remolded soil in sliding zone， the effect of water content on soil shear strength parameters is studied， and the slope stability analysis considering the effect of water content on soil shear strength parameters is realized based on the secondary development of GeoStudio software. The results show that the shear strength parameters decrease gradually with the increase of water content， and the cohesion and internal friction angle of sliding zone soil can be linearly related to water content as a whole. Ignoring the weakening effect of water on shear strength parameters of rock and soil will overestimate the stability of slope. The research results have certain theoretical significance for further understanding the mechanism of rainfall-induced landslides.

Keywords： unsaturated soil; slope stability; shear strength; water content; numerical simulation

受特殊的地質地貌和氣候條件影響，我國是世界上受泥石流、滑坡等地質災害危害最嚴重的國家之一[1]。大量的研究與工程實踐已證實，水對邊坡穩定性的影響是非常顯著的[2]，其致災機理主要在于邊坡內巖土體含水率發生變化導致土體重度增加、孔隙水壓力提高及土的抗剪強度降低等[3]。

為準確計算邊坡穩定性，眾多學者針對邊坡內部含水率受降雨入滲影響的分布變化規律，通過不同的途徑開展了大量研究。潘永亮等[4]基于改進Green-Ampt模型，推導了可適用不同降雨的入滲函數，理想地描述了邊坡降雨入滲情況。靳紅華等[5]基于Chen-Yong改進模型，對具有兩級臺階的高切邊坡進行研究，揭示了高切邊坡在不同降雨工況下的內部滲流場的演化規律。Tang等[6]基于三峽庫區實測降雨資料，研究了不同初始條件對不同土體邊坡穩定性影響，建議了能夠反映邊坡含水狀態的初始條件選取方法。史振寧等[7]通過室內一維降雨入滲模型試驗，研究了在考慮邊坡內部初始含水率分布時，坡體內部含水率在降雨入滲條件下的變化規律，為準確把握降雨入滲條件下坡體內部滲流情況及邊坡穩定性提供了理論依據。陳善雄等[8]基于極限平衡法，提出了一套考慮降雨入滲的非飽和土邊坡穩定性分析方法。李寧等[9]基于有限元強度折減法，在ABAQUS有限元軟件中將土體抗剪強度參數設置為場變量的函數，對邊坡穩定性進行分析。Tang等[10]基于簡化Bishop法和Janbu法，提出了考慮暫態飽和區的邊坡穩定性分析方法。這些研究為邊坡穩定性分析提供了重要參考，但忽略了水對巖土體抗剪強度參數的弱化作用，可能影響到邊坡穩定性評估的準確性。

滑帶土是指在滑坡的發生和發展過程中遭受擠壓、剪切、搓揉和研磨，而在滑帶內特定的物理化學條件下所形成的強度較低，呈可塑狀、軟塑狀甚至流塑狀的巖土體。其包含了多種黏土礦物成分，并沿滑坡面沉淀和聚集。滑帶土的含水率對其抗剪強度有著重要的影響[11]。

本文以東嶺信滑坡堆積體為研究對象，基于對滑帶重塑土的直剪試驗，研究含水率對土體抗剪強度參數的影響，并基于GeoStudio軟件二次開發實現了考慮含水率對土體抗剪強度參數影響的邊坡穩定性分析。

1? 研究區概況

東嶺信滑坡堆積體位于三板溪水電站庫區上游的右側，為清水江河流“Z”字形轉彎段，堆積體部位在河段右岸呈弧形凸出。此處原河床江面狹窄，寬50～60 m，原河床水面高程約427 m，河床高程423 m左右。堆積體的北東、東側為三板溪水庫臨庫岸坡，北西側和南側為山脊，約高程700 m。勘察資料顯示，東嶺信滑坡堆積體的地下水埋深較深。根據2019年4月至2022年4月連續3年DSW1、DSW2、DSW3和DSW4四個地下水位的監測數據，即使在雨水充沛的季節，受降雨影響庫區水位出現上升，該地區的地下水位仍表現穩定，受降雨影響很小（圖1）。

根據地質勘察資料，滑坡體滑帶處于邊坡體較深處，且整條滑帶沿地下水位分布，如圖2所示，滑帶的位置一般是邊坡土體抗剪強度較弱，抗滑力難以維持邊坡自身的穩定。從圖2中可以清楚看到，整條滑帶完全處于地下水位之下，而滑帶下方是基巖，所以可以理解為滑帶是浸在水中的土體，以至于其抗剪強度參數很容易被水弱化。

2? 試驗研究

為研究含水率對滑帶土體抗剪強度參數的影響規律，在實驗室內對滑帶土進行重塑制樣，通過直剪試驗測定不同含水率土體的抗剪強度參數，并用數學方法描述含水率與滑帶重塑土體抗剪強度參數的關系。

根據滑帶土的勘察分析報告，根據實測重塑滑帶土試樣的液限和塑限的結果，預設了5種含水率，分別為11%、13%、15%、17%和20%。按照GB/T 50123—2019《土工試驗方法標準》對不同含水率的各組滑帶重塑土樣進行直剪試驗，直剪儀的剪切速率設定為0.02 mm/min，最大剪切位移設定為4 mm。

試驗得到的不同含水率土體對應的抗剪強度參數如圖3所示。可見，隨著含水率的增加，土樣的黏聚力和內摩擦角均呈現減小趨勢，并且具有良好的線性關系。其中黏聚力c的擬合公式為c=-1.006 ω+39.033，內摩擦角?的擬合公式為?=-1.689 2ω+43.5，二者的擬合優度均在0.9以上，擬合效果較好。

3? 數值模擬分析

本節將設置一個降雨案例，利用GeoStudio軟件的SEEP/W和SLOPE/W 2個模塊，實現分析降雨條件下考慮含水率對土體抗剪強度參數影響的邊坡穩定性分析。

3.1? 數值模型

以東嶺信滑坡堆積體2-2剖面為基準建立數值分析模型，如圖4所示。模型包含穩態滲流、瞬態滲流和邊坡穩定分析3個部分，以穩態滲流分析結果作為瞬態滲流分析的初始條件，邊坡穩定性分析則基于瞬態滲流結果進行。

降雨工況選定為均勻型降雨，降雨強度為30 mm/d，屬于大雨等級，降雨歷時24 h，然后以模擬此工況的降雨為瞬態滲流分析的邊界條件，施加在模型地表面，得到此降雨工況作用下邊坡體內部巖土體積含水量分布情況如圖5所示。經分析，基巖以上的土體體積含水量等值線呈近似指數函數或直線函數形式分布，各等值線兩兩基本平行，并且隨著土體深度的增加相鄰2條曲線間的距離逐漸減小。

在本案例的降雨條件下，邊坡體基巖上部土體積含水量分布可以近似用數學的方式描述，具體數學關系式如式（1）所示

θ=0.000 45x+0.001 5y+0.835 5 ， （1）

式中：θ為體積含水量；x為橫坐標，y為縱坐標。

聯立含水率與滑帶重塑土體抗剪強度參數的數學模型和含水率與模型位置坐標的數學模型，建立滑帶土體抗剪強度參數與位置坐標的數學模型，如式（2）和式（3）所示，以確定該降雨工況下邊坡中某點土體抗剪強度參數。

c=-0.000 26x+0.000 88y+38.544 3 ， （2）

?=-0.000 44x+0.001 47y+42.679 4 ，（3）

式中：c為黏聚力，kPa；?為內摩擦角°。

然后利用C#語言將上述建立的數學關系式寫成動態鏈接庫，放入軟件接口里，供SLOPE/W執行程序在運行時調用。

3.2? 邊坡穩定性分析

根據以上考慮含水率對土體抗剪強度參數影響的方法，對邊坡穩定性進行計算得到的安全系數為1.019；而不考慮含水率對邊坡土軟化的影響，邊坡的安全系數為1.240。不考慮含水率對土抗剪強度參數的影響，將高估邊坡的穩定性。此外，對比圖6和圖7可以發現，考慮水對邊坡土的軟化作用時，邊坡的最危險滑移面是發生在滑帶位置處的；而未考慮水對邊坡土的軟化作用時，邊坡的最危險滑移面是發生在滑帶位置以上的土層的。可見，考慮水對土體的軟化作用時，邊坡的最危險滑移面發生位置較深，這也與東嶺信滑坡堆積體的勘察資料更一致。所以，考慮含水率對土體抗剪強度參數的影響，可以有效提高利用數值模擬分析非飽和土坡穩定性的準確性。

4? 結論

1）滑帶土室內慢剪試驗成果表明，隨著含水率的增加，土體的黏聚力和內摩擦角均呈現減小趨勢，并且具有良好的線性關系。

2）基于GeoStudio軟件二次開發實現了考慮含水率對土體抗剪強度參數影響的邊坡穩定性分析。結果表明，忽略水對巖土體抗剪強度參數的弱化作用會高估邊坡的穩定性。考慮水對土體的軟化作用時，邊坡的最危險滑移面發生位置更深，數值模擬的結果與勘察資料更接近。進行非飽和土邊坡穩定性分析時，考慮含水率對土體抗剪強度參數的影響是非常有必要的。

參考文獻：

[1] 唐棟，李典慶，周創兵，等.考慮前期降雨過程的邊坡穩定性分析[J].巖土力學，2013，34（11）：3239-3248.

[2] 蘇永華，李誠誠.強降雨下基于Green-Ampt模型的邊坡穩定性分析[J].巖土力學，2020，41（2）：389-398.

[3] 楊世豪，蘇立君，張崇磊，等.強降雨作用下昔格達邊坡滲流特性及穩定性分析[J].土木與環境工程學報（中英文），2020，42（4）：19-27.

[4] 潘永亮，簡文星，李林均，等.基于改進Green-Ampt模型的花崗巖殘積土邊坡降雨入滲規律研究[J].巖土力學，2020，41（8）：2685-2692.

[5] 靳紅華，王林峰，任青陽，等.基于Chen-Young模型及不平衡推力法的降雨型高切坡穩定性演化過程分析[J].重慶交通大學學報（自然科學版），2022，41（5）：107-117.

[6] TANG D，LI D， CAO Z. Slope stability analysis in the Three Gorges Reservoir Area considering effect of antecedent rainfall[J].Georisk： Assessment and Management of Risk for Engineered Systems and Geohazards，2017，11（2）：161-172.

[7] 史振寧，戚雙星，劉登生.考慮初始含水率分布的覆蓋層邊坡非飽和入滲過程與穩定性計算方法[J].工程地質學報，2022，30（4）：999-1009.

[8] 陳善雄，陳守義.考慮降雨的非飽和土邊坡穩定性分析方法[J].巖土力學，2001（4）：447-450.

[9] 李寧，許建聰.基于場變量的邊坡穩定分析有限元強度折減法[J].巖土力學，2012，33（1）：314-318.

[10] TANG D， JIANG Z， YUAN T， et al. Stability Analysis of Soil Slope Subjected to Perched Water Condition[J].Ksce Journal of Civil Engineering，2020，24（9）：2581-2590.

[11] 苗發盛，趙帆程，吳益平，等.基于滲透-環剪試驗的三峽庫區童家坪滑坡滑帶土強度特性研究[J].巖土工程學報，2023，45（7）：1480-1489.