曾鈴，李光裕，史振寧，邱祥，卞漢兵，李董可

?

前期降雨對(duì)非飽和覆蓋層邊坡性狀及穩(wěn)定性的影響

曾鈴1，李光裕1，史振寧2，邱祥2，卞漢兵1，李董可3

(1. 長(zhǎng)沙理工大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙，410114；2. 長(zhǎng)沙理工大學(xué) 交通運(yùn)輸工程學(xué)院, 湖南 長(zhǎng)沙，410114；3. 長(zhǎng)沙理工大學(xué)城南學(xué)院，湖南 長(zhǎng)沙，410114)

基于非飽和滲流及非飽和抗剪強(qiáng)度理論，以前期降雨強(qiáng)度為變量，對(duì)考慮前期不同降雨強(qiáng)度下的邊坡體積含水率、有效應(yīng)力場(chǎng)、位移場(chǎng)以及邊坡穩(wěn)定性變化過(guò)程進(jìn)行數(shù)值模擬分析，得到前期降雨強(qiáng)度變化對(duì)覆蓋層邊坡性狀以及穩(wěn)定性的影響。研究結(jié)果表明：前期降雨強(qiáng)度越大，無(wú)差別降雨作用后邊坡各位置含水率越大，對(duì)邊坡最終含水狀態(tài)起決定性影響的是前期降雨強(qiáng)度；無(wú)差別降雨過(guò)后，邊坡相同位置的有效應(yīng)力與前期降雨強(qiáng)度成反比，邊坡深部有效應(yīng)力下降幅度較小而邊坡坡腳處有效應(yīng)力下降幅度較大，邊坡坡腳處位移與前期降雨強(qiáng)度呈正比；無(wú)差別降雨過(guò)后，邊坡安全系數(shù)隨著前期降雨強(qiáng)度的增大而減小，前期降雨對(duì)邊坡穩(wěn)定性的影響主要表現(xiàn)在前期降雨強(qiáng)度越大，邊坡初始安全系數(shù)越小。

前期降雨；非飽和土；覆蓋層邊坡；邊坡穩(wěn)定性

第四紀(jì)沉積物形成時(shí)間短，成巖作用不充分，多為覆蓋在堅(jiān)硬基巖上松散、多孔隙的軟弱土層以及強(qiáng)風(fēng)化巖層，以此種形式存在的邊坡稱為覆蓋層邊坡，廣泛分布于我國(guó)大部分地區(qū)。大量研究資料表明，當(dāng)覆蓋層邊坡遭遇長(zhǎng)期降雨或短期極端降雨時(shí)極易發(fā)生滑坡、塌方等災(zāi)害，而路塹的開(kāi)挖則會(huì)進(jìn)一步導(dǎo)致此類地區(qū)公路兩側(cè)形成大量不良地質(zhì)邊坡[1?4]，對(duì)交通安全產(chǎn)生極大威脅。降雨會(huì)提高邊坡內(nèi)部孔隙水壓力，并且會(huì)對(duì)邊坡穩(wěn)定性產(chǎn)生明顯影響。目前，針對(duì)降雨對(duì)邊坡穩(wěn)定性影響的研究主要集中于對(duì)土體內(nèi)部孔隙水壓力變化的分析，如：張玉等[5]對(duì)古水水電站大型滑坡進(jìn)行了分析，認(rèn)為雨水入滲至滑帶土層降低了滑面處的有效應(yīng)力，進(jìn)而導(dǎo)致邊坡失穩(wěn)；CAREY等[6]對(duì)英國(guó)懷特島上Lowtherville滑坡進(jìn)行了分析，認(rèn)為孔隙水壓力增加導(dǎo)致有效應(yīng)力降低，最終導(dǎo)致滑坡發(fā)生。但現(xiàn)有研究中較少對(duì)土體中孔隙水壓力以及體積含水率分布等進(jìn)行分析。一般地，降雨是引發(fā)邊坡失穩(wěn)最常見(jiàn)也是最重要的因素。眾多研究表明，邊坡受后續(xù)降雨的影響程度與其前期降雨?duì)顟B(tài)密切相關(guān)，而且前期降雨也會(huì)對(duì)邊坡穩(wěn)定狀態(tài)起決定性作用，如 DAHAL等[7]分析了由降雨引起的尼泊爾小喜馬拉雅山滑坡，認(rèn)為前期降雨所造成的土體吸力喪失所導(dǎo)致的黏土礦物黏聚力降低是造成滑坡的主要原因。唐棟等[8]針對(duì)前期降雨過(guò)程對(duì)三峽庫(kù)區(qū)某邊坡穩(wěn)定性的影響進(jìn)行了分析，指出累計(jì)前期降雨量可以作為判斷邊坡是否穩(wěn)定的指標(biāo)之一。IBSEN等[9]認(rèn)為意大利Porretta- Vergato地區(qū)滑坡的原因主要為前期降雨，并提出對(duì)滑坡的啟動(dòng)分析應(yīng)考慮該地區(qū)6月之內(nèi)的降雨?duì)顟B(tài)。以上研究表明前期降雨與邊坡穩(wěn)定性具有緊密聯(lián)系，前期降雨強(qiáng)度越大，則后續(xù)降雨過(guò)后邊坡失穩(wěn)風(fēng)險(xiǎn)越大。目前的研究多是直接分析前期降雨和邊坡穩(wěn)定性這兩者之間的關(guān)系，并沒(méi)有進(jìn)一步討論邊坡體積含水率分布、邊坡應(yīng)力場(chǎng)以及邊坡位移場(chǎng)等變化的影響，而這些因素正是邊坡是否穩(wěn)定的關(guān)鍵[10]。結(jié)合邊坡在前期降雨影響下的應(yīng)力變化與位移變化對(duì)邊坡穩(wěn)定性進(jìn)行詳細(xì)分析具有重要意義。為此，本文作者基于現(xiàn)有研究成果，根據(jù)浙西山區(qū)當(dāng)?shù)氐刭|(zhì)與氣候條件，采用有限元方法對(duì)不同前期降雨條件下杭新景(杭州—新安江—景德鎮(zhèn))高速公路太真隧道上部非飽和覆蓋層邊坡性狀變化進(jìn)行模擬，并分析前期降雨與邊坡體積含水率分布狀態(tài)、有效應(yīng)力場(chǎng)、位移場(chǎng)以及穩(wěn)定性的關(guān)系，以便為同類型滑坡的滑移機(jī)理分析以及工程處置提供參考。

1 非飽和滲流及穩(wěn)定性計(jì)算方法

在非飽和滲流過(guò)程的有限元計(jì)算方法中，RICHARDS等[11]考慮非飽和土體的滲透系數(shù)與土體的基質(zhì)吸力存在函數(shù)關(guān)系，建立了土體中液體瞬時(shí)流動(dòng)的二維控制方程：

式中：m為基質(zhì)吸力水頭，kPa；(m)為非飽和土滲透系數(shù)，m/s；(m)為比水容量。

在非飽和土強(qiáng)度理論方面，F(xiàn)REDLUND等[12]將非飽和土中的基質(zhì)吸力概念引入傳統(tǒng)的抗剪強(qiáng)度分析，針對(duì)非飽和土抗剪強(qiáng)度提出了擴(kuò)展莫爾?庫(kù)侖準(zhǔn)則，即

將上述擴(kuò)展摩爾?庫(kù)侖公式代入普遍條分的極限平衡法中，得到邊坡穩(wěn)定性計(jì)算公式[15]：

式中：為擬定滑動(dòng)面的半徑；為垂直土條滑動(dòng)面方向的作用力；為土條底面長(zhǎng)度；為土條自重力；為土條重心距滑動(dòng)圓心水平距離；為滑動(dòng)圓心距作用力延長(zhǎng)線的垂直距離；為動(dòng)荷載下土條橫向作用力系數(shù)；為土條重心水平線距滑動(dòng)圓心距離；為點(diǎn)荷載；為滑動(dòng)圓心距作用線的距離；為坡頂拉裂縫充水后的水平推力；為滑動(dòng)圓心距作用線的距離?；率芰顟B(tài)如圖1所示。

圖1 普遍條分的極限平衡法土條受力狀態(tài)

2 工程概況

杭新景高速公路7標(biāo)段太真隧道位于浙西侵蝕剝蝕中低山區(qū)，地層為第四系殘坡積粉砂質(zhì)泥巖及炭質(zhì)泥巖夾煤層，含節(jié)理裂隙發(fā)育的黏性土碎石，完整性較差，滲透性較好，強(qiáng)度及抗剪性能差異較大。2013?01—08，該地區(qū)各月降雨量較往年同期較小，邊坡覆蓋層整體呈非飽和狀態(tài)。2013?08底，臺(tái)風(fēng)“潭美”帶來(lái)連續(xù)強(qiáng)降雨，導(dǎo)致山體失穩(wěn)坍塌，洞口滑塌后地形改變強(qiáng)烈，原隧道左側(cè)的山脊形態(tài)已完全破壞。滑塌體主要由風(fēng)化巖組成，上部有部分殘坡積的含黏性土碎石，基巖主要為強(qiáng)—中風(fēng)化細(xì)砂巖、粉砂質(zhì)泥巖，局部夾炭質(zhì)泥巖和石煤，呈碎塊狀、大塊狀?；碌湫涂v斷面如圖2所示。風(fēng)化砂巖與中風(fēng)化炭質(zhì)泥巖如圖3所示。在現(xiàn)場(chǎng)勘查過(guò)程中，發(fā)現(xiàn)邊坡裸露位置砂巖處于完全風(fēng)化狀態(tài)，經(jīng)降雨浸泡后基本失去原有強(qiáng)度。觀察隧道中開(kāi)挖出的棄方可知邊坡內(nèi)部炭質(zhì)泥巖與砂巖處于中風(fēng)化狀態(tài)，具有一定強(qiáng)度。

圖2 太真隧道上覆邊坡滑動(dòng)方向剖面圖

(a) 強(qiáng)風(fēng)化砂巖；(b) 中風(fēng)化炭質(zhì)泥巖

3 數(shù)值計(jì)算模型

根據(jù)邊坡典型剖面圖，采用有限元計(jì)算軟件SEEP/W建立分析模型，網(wǎng)格劃分以及測(cè)點(diǎn)選取如圖4所示?？紤]到地下水位線位于基巖下部，且覆蓋層邊坡滑動(dòng)面多位于基巖與覆蓋層交界處，在交界處設(shè)定4個(gè)測(cè)點(diǎn)。測(cè)點(diǎn)1~4的坐標(biāo)分別為(140，22)，(105，27)，(86，40)和(69，51)。綜合考慮計(jì)算效率與計(jì)算精度，模型共劃分為1 201個(gè)節(jié)點(diǎn)，1 162個(gè)網(wǎng)格，在滲流、應(yīng)力以及位移計(jì)算過(guò)程中，最大迭代步數(shù)為100步，規(guī)定迭代結(jié)果誤差小于1%即視為計(jì)算收斂。通過(guò)對(duì)計(jì)算結(jié)果進(jìn)行觀察，每一步計(jì)算均在100次迭代內(nèi)達(dá)到收斂，說(shuō)明計(jì)算結(jié)果精確度達(dá)到要求。

為分析不同前期降雨強(qiáng)度對(duì)邊坡?tīng)顟B(tài)的影響，在真實(shí)降雨?duì)顟B(tài)基礎(chǔ)上額外設(shè)定3種前期降雨?duì)顟B(tài)，即1~7月份降雨量分別為真實(shí)降雨強(qiáng)度的80%，120%和150%，8月份降雨量保持不變，具體降雨方案如表1所示。為區(qū)別具體的降雨過(guò)程，將1~7月份降雨過(guò)程稱為前期降雨，將8月份降雨過(guò)程稱為無(wú)差別降雨。在滲流計(jì)算中，定義坡面為單位流量邊界以模擬降雨過(guò)程，模型兩側(cè)設(shè)置為自由透水邊界，設(shè)置模型底部水平方向、豎直方向固定，左側(cè)水平方向固定。根據(jù)對(duì)現(xiàn)場(chǎng)合適地區(qū)所采集的土樣在室內(nèi)進(jìn)行力學(xué)性能和滲透性試驗(yàn)所得的結(jié)果，并參考該滑坡地段工程地質(zhì)勘察報(bào)告所提供的資料，得到巖土質(zhì)參數(shù)取值如表2所示。

圖4 邊坡有限元計(jì)算模型及測(cè)點(diǎn)位置

圖5所示為非飽和土的體積含水率與基質(zhì)吸力關(guān)系的擬合關(guān)系，采用式(5)所示VAN GENUCTHEN模型[16]對(duì)土水特征曲線進(jìn)行擬合。該模型的表達(dá)式中包含了曲線的轉(zhuǎn)折點(diǎn)信息，可以更有效地表征在進(jìn)氣壓力和趨近殘余含水量狀態(tài)時(shí)的平滑過(guò)渡情況。

式中：e為土體有效飽和度；為基質(zhì)吸力；，和為擬合參數(shù)。根據(jù)數(shù)值計(jì)算軟件中模塊自帶的樣本函數(shù)確定土體的滲透系數(shù)與基質(zhì)吸力關(guān)系曲線，如圖6所示。

表1 不同前期降雨量計(jì)算方案

表2 太真隧道上覆邊坡巖石性質(zhì)及參數(shù)

1—強(qiáng)風(fēng)化砂巖；2—中風(fēng)化砂巖；3—中風(fēng)化炭質(zhì)泥巖。

1—強(qiáng)風(fēng)化砂巖；2—中風(fēng)化砂巖；3—中風(fēng)化炭質(zhì)泥巖。

4 數(shù)值計(jì)算結(jié)果分析

4.1 邊坡體積含水率分析

基于非飽和滲流理論，采用有限元法對(duì)覆蓋層邊坡典型剖面進(jìn)行體積含水率分布計(jì)算。圖7所示為以不同前期降雨強(qiáng)度為初始條件的邊坡在無(wú)差別降雨后的體積含水率分布狀態(tài)。通過(guò)對(duì)比觀察可以發(fā)現(xiàn)：前期降雨強(qiáng)度對(duì)邊坡在無(wú)差別降雨之后的體積含水率有較大影響；在工況1下，無(wú)差別降雨后，邊坡表面以下5 m內(nèi)體積含水率約為0.20，飽和度約為60%；5 m以下區(qū)域土體體積含水率為0.10~0.15，飽和度為30%~50%；在工況3下，邊坡表面下部5 m內(nèi)體積含水率明顯上升，達(dá)到0.25左右，比工況2下略有升高；而在工況4下，在無(wú)差別降雨過(guò)后，邊坡表面10 m以下土體的體積含水率均達(dá)到0.25以上，坡頂、坡腳等部分區(qū)域也基本達(dá)到飽和狀態(tài)。由此可見(jiàn)：前期降雨強(qiáng)度越大，無(wú)差別降雨過(guò)后邊坡表面含水率越高，雨水入滲深度越大。

根據(jù)工程經(jīng)驗(yàn)，邊坡潛在滑動(dòng)面多經(jīng)過(guò)邊坡坡腳處，因此，本文將著重分析測(cè)點(diǎn)1(見(jiàn)圖4)處的體積含水率變化。從圖7(a)中的前期降雨過(guò)程可以看出：前期降雨結(jié)束后，土體的體積含水率隨著前期降雨強(qiáng)度的增加而顯著增加，而對(duì)于無(wú)差別降雨而言，此時(shí)即為初始狀態(tài)，對(duì)無(wú)差別降雨的影響起著決定性作用。測(cè)點(diǎn)1處體積含水率變化趨勢(shì)見(jiàn)圖8。從圖8可以看出：測(cè)點(diǎn)1處的體積含水率在不同工況下均有上升，但上升幅度并無(wú)明顯差異。由此可以得出：在后續(xù)降雨強(qiáng)度相同時(shí)，最終決定土體體積含水率的是前期降雨?duì)顟B(tài)。

4.2 應(yīng)力場(chǎng)分析

由于邊坡內(nèi)部有效應(yīng)力的變化會(huì)導(dǎo)致覆蓋層邊坡穩(wěn)定性發(fā)生變化[17?19]，因此，有必要對(duì)降雨入滲過(guò)程中邊坡內(nèi)部有效應(yīng)力的變化規(guī)律進(jìn)行細(xì)致分析。對(duì)4個(gè)測(cè)點(diǎn)(如圖4所示)在計(jì)算過(guò)程中豎直方向的有效應(yīng)力變化進(jìn)行觀察，結(jié)果如圖9所示。對(duì)比4個(gè)測(cè)點(diǎn)的有效應(yīng)力變化趨勢(shì)發(fā)現(xiàn)：由于前期降雨的作用，各個(gè)測(cè)點(diǎn)在不同工況下的初始狀態(tài)均有著明顯差別?？傮w而言，前期降雨強(qiáng)度越大，在初始狀態(tài)下的有效應(yīng)力越低，而且在無(wú)差別降雨過(guò)程中不同工況下各個(gè)測(cè)點(diǎn)處的有效應(yīng)力下降量基本相同。由此可見(jiàn)：前期降雨強(qiáng)度對(duì)邊坡內(nèi)部有效應(yīng)力場(chǎng)的影響在前期降雨結(jié)束時(shí)各個(gè)測(cè)點(diǎn)的有效應(yīng)力上已經(jīng)有所體現(xiàn)，無(wú)差別降雨過(guò)程對(duì)有效應(yīng)力降低的影響程度無(wú)明顯差異。

(a) 工況1 ；(b) 工況2； (c) 工況3；(d) 工況4

1—工況1；2—工況2；3—工況3；4—工況4。

圖10所示為測(cè)點(diǎn)1~4在不同降雨工況下有效應(yīng)力下降幅度的變化。從圖10可見(jiàn)：有效應(yīng)力降低程度與測(cè)點(diǎn)所在位置有關(guān)，具體體現(xiàn)在邊坡內(nèi)部的測(cè)點(diǎn)4在不同工況下達(dá)到初始狀態(tài)時(shí)有效應(yīng)力基本不變；測(cè)點(diǎn)3在不同工況下達(dá)到初始狀態(tài)時(shí)的有效應(yīng)力差值略比測(cè)點(diǎn)4的大；而測(cè)點(diǎn)2在不同工況下達(dá)到初始狀態(tài)時(shí)的有效應(yīng)力差值則明顯增加；測(cè)點(diǎn)1在不同工況下達(dá)到初始狀態(tài)時(shí)，豎直方向有效應(yīng)力差值最大，達(dá)到12.2 kPa，明顯比其他測(cè)點(diǎn)的大，可見(jiàn)降雨對(duì)邊坡影響最大的位置是邊坡坡腳處。產(chǎn)生以上現(xiàn)象的原因主要是坡腳處上覆土層較薄，雨水入滲所用時(shí)間較短，導(dǎo)致豎向有效應(yīng)力下降最快，幅度最大，而邊坡內(nèi)部上覆土層較厚，降雨入滲速度較小。

4.3 位移場(chǎng)分析

在降雨過(guò)程中，隨著土體體積含水率增加，有效應(yīng)力逐漸減小，邊坡表面及坡腳處會(huì)產(chǎn)生一定程度變形。為分析不同前期降雨強(qiáng)度對(duì)邊坡位移的影響，列出不同工況下無(wú)差別降雨后的邊坡水平方向位移，如圖11所示。從圖11可以看出：邊坡在降雨作用下坡面與坡腳均產(chǎn)生一定位移，而坡腳處的位移較大，在工況1~4條件下的位移分別為18.7，20.0，20.9和22.8 mm，可見(jiàn)坡腳在無(wú)差別降雨之后的位移隨前期降雨強(qiáng)度的增大而增大。而從圖12可見(jiàn)：在前期降雨過(guò)程中，邊坡坡腳處水平方向位移隨著降雨持續(xù)而增大，前期降雨強(qiáng)度越大，初始狀態(tài)下的坡腳位移也越大；而在無(wú)差別降雨過(guò)程中，不同工況下坡腳處位移增加量基本相同。因此，可以認(rèn)為前期降雨對(duì)邊坡坡腳位移的影響主要體現(xiàn)在邊坡在初始狀態(tài)下存在一定位移，并且初始狀態(tài)下的位移隨前期降雨強(qiáng)度的增大而增大。

(a) 測(cè)點(diǎn)1； (b) 測(cè)點(diǎn)2； (c) 測(cè)點(diǎn)3；(d) 測(cè)點(diǎn)4

圖10 前期降雨引起邊坡豎向有效應(yīng)力下降趨勢(shì)

4.4 邊坡穩(wěn)定性分析

采用普遍條分極限平衡法對(duì)不同前期降雨過(guò)程中的邊坡穩(wěn)定性進(jìn)行分析，得到邊坡安全系數(shù)變化規(guī)律如圖13所示。從圖13可以看出：邊坡初始安全系數(shù)為2.2，處于穩(wěn)定狀態(tài)，但隨著前期降雨持續(xù)，邊坡安全系數(shù)逐漸下降；在工況1~4條件下，邊坡在經(jīng)歷前期降雨后的初始狀態(tài)安全系數(shù)分別為1.506，1.223，1.004和0.776，可見(jiàn)前期降雨強(qiáng)度越大，邊坡在初始狀態(tài)下的安全系數(shù)越小。值得注意的是：在工況3和工況4下，邊坡在無(wú)差別降雨結(jié)束之前安全系數(shù)已經(jīng)小于1.000，邊坡處于失穩(wěn)狀態(tài)，但為了對(duì)比分析邊坡在無(wú)差別降雨過(guò)程中的含水率、有效應(yīng)力以及位移變化，仍考慮邊坡處于極限平衡狀態(tài)，并未發(fā)生破壞。在無(wú)差別降雨過(guò)程中，邊坡安全系數(shù)進(jìn)一步下降，無(wú)差別降雨結(jié)束后，工況1~4下的邊坡最終安全系數(shù)分別為1.492，1.130，0.917和0.703，與初始狀態(tài)相比，下降量基本相同。據(jù)以上分析，可以認(rèn)為前期降雨對(duì)邊坡穩(wěn)定性的影響主要表現(xiàn)在前期降雨強(qiáng)度對(duì)初始狀態(tài)下邊坡安全系數(shù)的影響。

(a) 工況1；(b) 工況2；(c) 工況3；(d) 工況4

1—工況1；2—工況2；3—工況3；4—工況4。

1—工況1；2—工況2；3—工況3；4—工況4。

5 結(jié)論

1) 前期降雨強(qiáng)度越大，無(wú)差別降雨作用后的邊坡各位置含水率越大，無(wú)差別降雨過(guò)程中邊坡內(nèi)部含水率增加量基本相同，對(duì)邊坡最終含水狀態(tài)起決定性影響的是前期降雨強(qiáng)度。

2) 邊坡相同位置的有效應(yīng)力與前期降雨強(qiáng)度成反比。在相同降雨強(qiáng)度下，邊坡深部有效應(yīng)力下降幅度較小，而邊坡坡腳處有效應(yīng)力下降幅度較大。邊坡坡腳處位移與前期降雨強(qiáng)度密切相關(guān)，降雨對(duì)邊坡內(nèi)部位移影響較小。

3) 無(wú)差別降雨過(guò)后邊坡的穩(wěn)定性與前期降雨強(qiáng)度成反比。前期降雨對(duì)邊坡穩(wěn)定性的影響主要表現(xiàn)在前期降雨強(qiáng)度越大，邊坡初始安全系數(shù)越小。

[1] 李濤, 付宏淵, 周功科, 等. 降雨入滲條件下粗粒土路堤暫態(tài)飽和區(qū)發(fā)展規(guī)律及穩(wěn)定性研究[J]. 水文地質(zhì)工程地質(zhì), 2013, 40(5): 74?80. LI Tao, FU Hongyuan, ZHOU Gongke, et al. A study of development law and stability of transient saturated areas of coarse-grained soil embankment under rainfall infiltration[J]. Hydrogeology & Engineering Geology, 2013, 40(5): 74?80.

[2] 馬世國(guó), 韓同春, 徐日慶, 等. 強(qiáng)降雨條件下含傾斜基巖層的邊坡穩(wěn)定分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 46(7): 2673?2678. MA Shiguo, HAN Tongchun, XU Riqing, et al. Stability analysis of slope with inclined bedrock beneath under intense rainfall[J]. Journal of Central South University (Science and Technology), 2015, 46(7): 2673?2678.

[3] 張俊然, 許強(qiáng), 孫德安. 吸力歷史對(duì)非飽和土力學(xué)性質(zhì)的影響[J]. 巖土力學(xué), 2013, 34(10): 2810?2814. ZHANG Junran, XU Qiang, SUN Dean. Effect of suction history on mechanical behavior of unsaturated soils[J]. Rock and Soil Mechanics, 2013, 34(10): 2810?2814.

[4] 劉子振, 言志信, 凌松耀, 等. 非飽和土邊坡抗剪強(qiáng)度的力學(xué)參數(shù)影響及靈敏度分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 43(11): 4508?4513. LIU Zizhen, YAN Zhixin, LING Songyao, et al. Sensitivity analysis and influence of mechanical parameters on shear strength of unsaturated soil slope[J]. Journal of Central South University (Science and Technology), 2012, 43(11): 4508?4513.

[5] 張玉, 徐衛(wèi)亞, 鄒麗芳, 等. 降雨條件下大型滑坡體滲流穩(wěn)定性分析[J]. 巖土力學(xué), 2013, 34(3): 833?841. ZHANG Yu, XU Weiya, ZOU Lifang, et al. Analysis of seepage stability of large-scale landslide under rainfall condition[J]. Rock and Soil Mechanics, 2013, 34(3): 833?841.

[6] CAREY J M, PETLEY D N. Progressive shear-surface development in cohesive materials, implications for landslide behavior[J]. Engineering Geology, 2014, 177(14): 54?65.

[7] DAHAL R K, HASEGAWA S, YAMANAKA M. Comparative analysis of contributing parameters for rainfall-triggered landslides in the Lesser Himalaya of Nepal[J]. Environment Geology 2009, 58(3): 567?586.

[8] 唐棟, 李典慶, 周創(chuàng)兵, 等. 考慮前期降雨過(guò)程的邊坡穩(wěn)定性分析[J]. 巖土力學(xué), 2013, 34(11): 3239?3248. TANG Dong, LI Dianqing, ZHOU Chuangbing, et al. Slope stability analysis considering antecedent rainfall process[J]. Rock and Soil Mechanics, 2013, 34(11): 3239?3248.

[9] IBESN M L, CASAGLI N. Rainfall patterns and related landslide incidence in the Porretta—Vergato region, Italy[J]. Landslides, 2004, 1(2): 143?150.

[10] 陳立平, 張頂立, 房倩, 等. 一種基于數(shù)值應(yīng)力場(chǎng)的強(qiáng)度各向異性邊坡穩(wěn)定性分析方法[J]. 巖土力學(xué), 2014, 35(12): 3611?3618. CHEN Liping, ZHANG Xiangli, FANG Qian, et al. A new search algorithm for strength anisotropy slope stability analysis based on numerical stress field[J]. Rock and Soil Mechanics, 2014, 35(12): 3611?3618.

[11] RICHARDS L A. Capillary conduction of liquids through porous medium[J]. Journal of Physics, 1931, 1(5): 318?333.

[12] FREDLUND D G, MORGENSTERN N R, WIDGER R A. Shear strength of unsaturated soils[J]. Canadian Geotechnical Journal, 1978, 15(3): 313?321.

[13] VANAPALLI S K, FREDLUND D G. Comparison of different procedures to predict unsaturated soil shear strength[J]. Geotechnical Special Publication, 2000, 287(99): 195?209.

[14] ESCARIO V, JUCA J, COPPE M S. Strength and deformation of party saturated soils[C]// Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering. Riode Janerio, Brazil, 1989: 43?46.

[15] JOHN K. Geo-slope International Ltd. Stability modeling with slope/W2007 version[R]. Calgary, Canada: Geo-slope International Ltd, 2008: 315?331.

[16] VAN GENUCHTEN M T. A closed form equation for prediction the hydraulic conductivity of unsaturated soils[J]. Soil science society of American Journal, 1980, 44: 892?898.

[17] 陳正漢. 非飽和土與特殊土力學(xué)的基本理論研究[J]. 巖土工程學(xué)報(bào), 2014, 36(2): 201?272. CHEN Zhenghan. On basic theories of unsaturated soils and special soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(2): 201?272.

[18] 李濤. 考慮降雨及開(kāi)挖影響下的厚覆蓋層邊坡滲流特征及穩(wěn)定性[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 47(5): 1708?1714. LI Tao. Seepage characteristics and stability of overburden slope considering rainfall and excavation[J]. Journal of Central South University (Science and Technology), 2016, 47(5): 1708?1714.

[19] 繆海繽, 王建國(guó), 費(fèi)曉歐, 等. 基于孔隙水壓力消散的排土場(chǎng)邊坡動(dòng)態(tài)穩(wěn)定性研究[J]. 煤炭學(xué)報(bào), 2017, 42(9): 2302?2306. MIAO Haibin, WANG Jianguo, FEI Xiaoou, et al. Study on dynamic stability of dump slope based on the dissipation of pore pressure[J]. Journal of China Coal Society, 2017, 42(9): 2302?2306.

(編輯 陳燦華)

Effect of antecedent rainfall on characteristics and stability of unsaturated overburden slope

ZENG Ling1, LI Guangyu1, SHI Zhenning2, QIU Xiang2, BIAN Hanbing1, LI Dongke3

(1. School of Civil Engineering, Changsha University of Science & Technology, Changsha 410114, China; 2. School of Traffic and Transportation Engineering, Changsha University of Science & Technology, Changsha 410114, China;3. Chengnan College, Changsha University of Science & Technology, Changsha 410114, China)

Based on the unsaturated seepage and unsaturated soil shear strength theory, the rainfall intensity was used as an variables. Numerical simulation of the volumetric water content changing process, effective stress field, displacement field and slope stability were studied in different conditions of antecedent rainfall intensity. The impact of the antecedent rainfall intensity on the properties and stability of unsaturated overburden slope was obtained. The results show that larger antecedent rainfall intensity leads to greater water content of the slope after the discriminate rainfall. The antecedent rainfall intensity plays a decisive role in the final water content of the slope. After the discriminate rainfall, the effective stress in the same position of the slope is inversely proportional to the intensity of the antecedent rainfall. The decrease of effective stress is small inside the slope but great at the bottom of the slope. The displacement at the bottom of slope is proportional to the rainfall intensity. After the discriminate rainfall, the slope safety factor decreases with the increase of antecedent rainfall intensity, the effect of rainfall on slope stability is mainly manifested as greater antecedent rainfall intensity leads to smaller initial safety factor of the slope.

antecedent rainfall; unsaturated soil; overburden slope; slope stability

TU457

A

1672?7207(2018)01?0238?09

10.11817/j.issn.1672-7207.2018.01.030

2017?03?20；

2017?05?06

國(guó)家自然科學(xué)基金資助項(xiàng)目(56178073, 51678074, 51508040, 51508079)；湖南省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016SK2023)；湖南省教育廳優(yōu)秀青年基金資助項(xiàng)目(17B013) (Projects(56178073, 51678074, 51508040, 51508079) supported by the National Natural Science Foundation of China; Project(2016SK2023) supported by Key Research and Development Plan of Hunan Province,; Project(17B013) supported by Outstanding Youth Foundation of Education Department in Hunan Province)

曾鈴，博士，講師，從事巖土工程、道路工程研究；E-mail:zlbingqing3@126.com