降雨入滲非飽和花崗巖殘積土邊坡穩定性分析

張淑寶

(漳州廈蓉高速擴建工程有限公司 福建漳州 363000)

降雨入滲非飽和花崗巖殘積土邊坡穩定性分析

張淑寶

(漳州廈蓉高速擴建工程有限公司 福建漳州 363000)

針對福建地區典型邊坡，基于非飽和滲流理論以及非飽和土力學理論，分析降雨入滲條件下邊坡形態(高度、坡率)、物理力學參數(密度、有效黏聚力、有效內摩擦角)等對邊坡滲流場與穩定性的影響機理。結果表明：降雨對矮邊坡穩定性的影響明顯更強。降雨過程中，天然含水量越大，邊坡安全系數下降越快。邊坡穩定性隨土體有效粘聚力和有效內摩擦角的增大而增大，變化率分別為0.0426/kPa、0.0319/°。粗粒土邊坡滑裂面基本為直線，易形成邊坡淺部失穩。增大有效粘聚力降低有效內摩擦角，邊坡滑裂面逐漸向邊坡深部發展，同時滑裂面的形狀逐漸由直線轉變為曲線型。

路基邊坡；降雨；非飽和滲流；力學響應

0 引言

福建地區地處山區，素有“八山一水一分田”之稱，境內存在大量降雨型滑坡，據統計，2006年“珍珠”“桑美”等臺風特大暴雨，引發了3000多處地質災害。2009年的“莫拉克”臺風導致福建省發生了滑坡共計389處。因此，降雨誘發邊坡失穩機理研究始終是工程界的熱點。自從考慮非飽和滲流和非飽和土力學理論后，國內外學者對降雨誘發非飽和邊坡失穩機理進行了大量研究，獲得了較為豐碩的成果[1-5]，這些成果對于進一步認識降雨型滑坡失穩機理具有重要意義。然而，上述成果對于花崗巖殘積土研究較少，實際上，福建省境內近30%面積分布廣泛花崗巖。基此，本文針對福建省典型花崗巖殘積土邊坡及降雨特點，采用有限元Geo-Studio軟件中的Seep/W和slope模塊模擬降雨入滲的過程，分析降雨入滲條件下邊坡形態、物理力學參數等對非飽和邊坡滲流場和穩定性的影響規律。

1 非飽和土抗剪強度理論

采用Fredlund等1978年提出非飽和強度準則[6]：

τf=c′+σn-ua)tanφ′+(ua-uw)tanθb

(1)

式中，

τf為非飽和土抗剪強度；

c′、φ′分別為有效粘聚力和有效內摩擦角；

φb為吸力內摩擦角；σn為總應力；

ua-uw為土體的基質吸力；

Δc′=(ua-uw)tanφb表示有效凝聚力隨(ua-uw)的增量。

2 數值分析算例與及計算工況

2.1 計算模型及參數

數值模擬網格劃分如圖1所示，邊坡坡率為1∶0.75，由上而下分兩層，基本物理力學參數如表1所示。土-水特征曲線以及滲透曲線均是含水率的函數，其中土-水特征曲線采用濾紙法測得，并采用V-G模型擬合得到，如圖2所示。降雨強度的選取結合臺風期降雨強度觀測資料，取50mm/d。

圖1 數值計算模型(單位：m)

土層γ(kN/m3)E(MPa)?(%)ks(m/s)c′(kPa)φ′(°)φb(°)表層181039．251．588e-6252010下層192039．251．588e-8352613

注：γ為土體的重度，E為土體的楊氏模量，?s為土體的飽和體積含水量，ks為土體的飽和滲透系數。

2.2 降雨邊界

降雨邊界分兩種工況：

(1)當降雨強度小于土體入滲能力時，雨水全部入滲，入滲速率取降雨強度，邊界條件為Neumann邊界[7-10]：

(2)

(2)當降雨強度大于土體入滲能力時，將形成積水模式，入滲速率取飽和滲透系數，邊界條件轉為Dirichlet邊界：

(a)土-水特征曲線

(b)上層土滲透曲線

(c)底層土滲透曲線圖2 土-水特征曲線和滲透系數曲線

(3)

3 不同條件邊坡的穩定性

3.1 邊坡形態的影響

3.1.1路基高度的影響

表2為不同路塹的邊坡高度在50mm/d降雨條件下穩定性計算結果。由表2可見，不考慮降雨條件下，邊坡高度由6m增加至14m時，其安全系數由1.777下降到1.121，其變化率(安全系數差值與邊坡高度差值比)為0.082/m，可見邊坡高度對邊坡穩定性的影響非常明顯。在降雨作用下，矮邊坡受到的影響大，如在降雨持續24h時，6m與12m高邊坡安全系數分別下降了0.04、0.008，原因是由于邊坡較高時，雨水下滲到坡底深處所需要的時間更長。

表2 不同路塹邊坡高度的安全系數

3.1.2路基坡率的影響

取0.5、0.75、1、1.5等4種不同坡率，如圖3所示，計算邊坡在50mm/d降雨強度持續24h時的安全系數，結果如表3所示。由表3可知，不考慮降雨條件下，邊坡坡率由0.5增加到1.5時，其安全系數由1.356上升到1.961。而在降雨24h后，各種坡率的邊坡安全系數分別下降0.021、0.014、0.026、0.025。由此說明，在降雨持時較短時，表層土體水分未滲入邊坡，因而邊坡穩定性與坡率的關系不大。

圖3 不同坡率模擬模型

坡率持時(h)安全系數0．50．751．01．501．3561．4511．6221．96161．3501．4441．6181．959121．3451．4391．6151．954241．3351．4371．5961．936

3.2 土體的物理力學參數的影響

3.2.1天然含水量的影響

表4為降雨強度為50mm/d時，不同天然含水量邊坡安全系數隨持時的變化規律。不考慮降雨條件時，邊坡的穩定性隨土體天然含水量的增大而減小，原因在于當含水量較低時，土坡體比較干燥，基質吸力比較大，在降雨條件下，雨水入滲大部分首先被上部土體吸收，下滲的雨量不多，且速度較慢，除了在接近坡腳的地方有所抬高，深部土體的含水量改變不多。如含水量從30.6%上升到39.2%，邊坡安全系數下降了0.254。在降雨過程中，天然含水量越大，邊坡安全系數下降越快。原因在于天然含水量越大，坡體的飽和程度越高，雨水的入滲速率快，在相同降雨條件下，水位上升得更快、更高，土體抗剪強度越低。如在降雨持續6h時，30.6%條件下僅下降了0.007，而含水量為39.2%的情況卻下降了0.296。

表4 不同天然含水量下邊坡的安全系數隨持時變化

3.2.2土體重度的影響

降雨強度為50mm/d時，不同土體重度的邊坡安全系數在降雨前與降雨后24h的計算結果如表5所示。由表5可知，不考慮降雨的情況下，重度從16kN/m3增至19.5kN/m3時，其安全系數由1.601降至1.387，降幅達0.214；在降雨持續24h后，不同重度對應的邊坡，其安全系數變化值不同，如重度為16kN/m3的安全系數下降了0.046，而重度為19.5kN/m3時的安全系數僅下降了0.016。因此，重度越大，其邊坡穩定性對降雨的影響越不敏感。

表5 不同重度時路塹邊坡的安全系數

3.2.3有效粘聚力的影響

降雨強度為50mm/d時，不同有效粘聚力的邊坡安全系數在降雨后24h的計算結果如圖4所示。在相同條件下，邊坡穩定性隨土體有效粘聚力的增大而增大，變化率為0.0426/kPa。黨靖[6]等通過試驗研究了土體強度參數與含水量的關系表明，強度參數粘聚力和內摩擦角隨含水量的增加而減小，當含水量達26%后，粘聚力隨含水量增加而減小的趨勢減緩。由此可見，雨水入滲，一方面使粘土中的粘土礦物產生水化作用，粘土中產生離子交換作用，即將易溶膠結物帶走，顆粒間的結合力減弱，有效粘聚力降低，同時，雨水的潤滑作用使得顆粒之間的摩擦力降低。另一方面，坡體含水量增大，離子的濃度逐漸降低，顆粒間擴散層厚度增大，強度降低，即降雨入滲使得邊坡土體的有效粘聚力降低。

圖5為不同有效粘聚力下的邊坡最危險滑裂面的形狀示意圖。可見，當有效粘聚力接近零時，土體由粗粒組成，如砂或礫石土等，其滑裂面基本為直線，并且滑裂面一般位于邊坡淺部，從而易形成邊坡淺部失穩；隨著有效粘聚力的增大，邊坡滑裂面逐漸向邊坡深部發展，同時滑裂面的形狀逐漸由直線轉變為曲線型。

圖4 安全系數-有效粘聚力關系曲線

3.2.4有效內摩擦角的影響

圖6為在降雨24h(降雨強度為50mm/d)后不同有效內摩擦角的邊坡安全系數計算結果。穩定性系數隨土體有效內摩擦角與有效粘聚力有相同的變化趨勢，其變化率為0.0319/°。圖7為不同有效內摩擦角下的邊坡最危險滑裂面的形狀示意圖。當土體有效內摩擦角較大時，其土體粗顆粒較多，因而滑裂面接近一直線，并且位于邊坡淺層。當土體有效內摩擦角逐漸較小時，邊坡的滑裂面逐漸向邊坡深層發展，同時，滑裂面的形狀也隨之變為曲線型，這與有效粘聚力的變化趨勢相反。

4 結論

(a)c′=10kPa

(b)c′=20kPa

(c)c′=35kPa圖5 不同有效粘聚力在50mm/d降雨強度作用24h時最危險滑裂面

圖6 安全系數-有效內摩擦角關系曲線

(1)在降雨作用下，邊坡高度對邊坡穩定性的影響非常明顯。由于邊坡較高時，雨水下滲到坡底深處所需要的時間更長，因而相對于高邊坡，矮邊坡穩定性受到的影響更大。在降雨持續較短時，表層土體水分未滲入邊坡，因而邊坡穩定性與坡率的關系不大。

(a)φ′=5°

(b)φ′=15°

(c)φ′=25°

(d)φ′=35°圖7 不同有效粘內摩擦角在50mm/d降雨強度作用24時最危險滑裂面

(2)不考慮降雨條件時，邊坡的穩定性隨土體天然含水量的增大而減小。在降雨過程中，天然含水量越大，邊坡安全系數下降越快。邊坡重度越大，其邊坡穩定性對降雨的影響越不敏感。

(3)有效粘聚力和有效內摩擦角對邊坡穩定性影響顯著，對滑裂面位置影響較大。在相同條件下，邊坡穩定性隨土體有效粘聚力和有效內摩擦角的增大而增大，變化率分別為0.0426/kPa、0.0319/°。

(4)土體由粗粒組成，如砂或礫石土等，其滑裂面基本為直線，并且滑裂面一般位于邊坡淺部，從而易形成邊坡淺部失穩；隨著有效粘聚力的增大或有效內摩擦角的降低，邊坡滑裂面逐漸向邊坡深部發展，同時滑裂面的形狀逐漸由直線轉變為曲線型。

[1] Sun Y，Nishigaki M，Kohno I.A stusy on stability analysis of shallow layer slope due to raining permeation[C]//In：Proc.of the 1st Int.Conf.on Unsaturated Soil.Rotterdam：Balkema，1995，1135-1141.

[2] Fredlund D G，Xing A.The relationship of the unsaturated soil shear strength to the soil-water characteristic curve[J].Can.Geotech.J.,1996，33：440-448.

[3] 徐晗，朱以文，蔡元奇，等.降雨入滲條件下非飽和土邊坡穩定分析[J].巖土力學，2005，26(12):1957-1962.

[4] 魏寧，茜平一，傅旭東.降雨和蒸發對土質邊坡穩定性的影響[J].巖土力學，2006，27(5):778-786.

[5] 崔溦，楊建,劉麗.連續降雨作用下非飽和土邊坡的穩定性分析[J].水文地質工程地質，2007(3)：55-57.

[6] 黨靖，胡李俐，南帥.含水量及天然密度對土體抗剪強度參數的影響研究[J].中國西部科技，2007，9 (5)：14-15.

[7] 唐棟，李慶典，周創.考慮前期降雨過程的邊坡穩定性分析[J].巖土力學，2013,34(11)：3240-3248.

[8] 鄧福根，金愛兵，李龍福，等.降雨和蒸發對粉質黏土邊坡穩定性影響分析[J].采礦工程，2012,5(33)：27-31.

[9] 劉子振，言志信.降雨入滲非飽和黏土邊坡穩定性的極限平衡條分法研究[J].巖土力學，2016，37(2)：351-356.

[10] 李庚，陳超，蔣斌，等.降雨入滲的紅黏土邊坡穩定性分析研究[J].山東交通科技，2017(1)：36-74.

StabilityAnalysisonUnsaturatedGraniteResidualSoilSlopeUndertheRainfallAction

ZHANGShubao

(Xiarong Expressway Expansion Project Limited Company,Zhangzhou 363000)

This paper indicates that the influence mechanism of some factors including slope shape (e.g.height,slope),physical and mechanical parameters (e.g.density,effective cohesion and effective internal friction angle) on slope seepage field and stability for typical slope in Fujian Province based on the unsaturated seepage theory and unsaturated soil theory.The results show that rainfall has greater influence on low slope.The slope with high water content is easy to be instable during the rainfall process.The safety factor of slope increases with the increase of the effective cohesion and effective internal friction angle of the soil,and the change rates are 0.0426 / kPa and 0.0319 /°,respectively.The slippery surface of the coarse-grained soil slope is straight and shallow slope will tend to be instable.With the increase of the effective cohesion or the decrease of the effective internal friction angle,the slope slip surface will locate in the deep of slope and the shape of the sliding surface will tend to curve gradually.

Rainfall; Slope; Unsaturated seepage; Mechanical response

張淑寶(1972- )，男，高級工程師。

E-mail:804298921@qq.com

2017-08-21

U416.1

A

1004-6135(2017)12-0074-05