考慮不同降雨概率和雨型的滑坡滲流場(chǎng)及穩(wěn)定性變化率特征分析

胡 濤， 樊 鑫， 王 碩

(1.中國(guó)地質(zhì)大學(xué)工程學(xué)院，武漢 430074；2.貴州省地質(zhì)工程勘查院，貴陽(yáng) 550002；3.貴州省銅仁市應(yīng)急管理局，銅仁 554300)

降雨是誘發(fā)滑坡的一個(gè)非常重要的外部因素。據(jù)統(tǒng)計(jì)，中國(guó)已發(fā)生的滑坡中大部分都是由降雨直接誘發(fā)或與降雨作用有關(guān)[1-4]。降雨誘發(fā)滑坡的機(jī)制主要是降雨導(dǎo)致滑坡地下水位上升、滲透壓力增大；且降雨增加了坡體自重并導(dǎo)致其抗剪強(qiáng)度降低，進(jìn)而誘發(fā)滑坡發(fā)生[5-7]。

在強(qiáng)降雨作用下，貴州省思南縣滑坡地質(zhì)災(zāi)害頻發(fā)，給當(dāng)?shù)厝罕娚a(chǎn)生活帶來較大危害，也引起了諸多研究人員的重視。董遠(yuǎn)峰等[8]采用連續(xù)介質(zhì)法對(duì)貴州省思南滑坡的危險(xiǎn)性進(jìn)行研究；張抽勇[9]對(duì)思南縣延安路某滑坡的誘發(fā)因素和地質(zhì)特征等進(jìn)行了分析；黃姍等[10]對(duì)思南縣的低緩松散堆積層滑坡的失穩(wěn)機(jī)制進(jìn)行了深入探討。由上述文獻(xiàn)可知，思南縣滑坡以中型松散土質(zhì)滑坡為主，其滑坡誘發(fā)因素主要是強(qiáng)降雨過程和不良人類工程活動(dòng)。可見，對(duì)思南縣降雨誘發(fā)滑坡的失穩(wěn)機(jī)制進(jìn)行深入研究具有重要意義。

各種降雨工況下(降雨強(qiáng)度、降雨歷時(shí)和降雨模式等工況的組合)的滑坡穩(wěn)定性變化規(guī)律研究得到了諸多學(xué)者的關(guān)注。丁志洋等[11]研究了不同雨強(qiáng)和降雨歷時(shí)等作用下的土質(zhì)邊坡地下水入滲規(guī)律；任佳等[12]研究了雨強(qiáng)對(duì)樹坪滑坡穩(wěn)定性的影響；吳宏偉等[13]針對(duì)香港地區(qū)降雨滑坡的研究表明，降雨強(qiáng)度、雨型及降雨持續(xù)時(shí)間對(duì)滑坡地下水浸潤(rùn)線和穩(wěn)定性系數(shù)變化具有非常重要的影響; Huang等[14]深入探究了雨強(qiáng)、坡體滲透性等對(duì)三峽庫(kù)區(qū)庫(kù)岸滑坡穩(wěn)定性變化的作用規(guī)律；馬紫娟等[15]對(duì)降雨強(qiáng)度、降雨歷時(shí)和地震等組合工況下的禮縣爛山滑坡穩(wěn)定性變化特征進(jìn)行了深入研究；楊背背等[16]研究了不同雨強(qiáng)對(duì)堆積層滑坡穩(wěn)定性變化規(guī)律的影響。與此同時(shí)，也有較多學(xué)者從降雨誘發(fā)滑坡的數(shù)值模擬角度出發(fā)，探討降雨和滑坡物理力學(xué)機(jī)制變化之間的響應(yīng)關(guān)系。李德心[17]基于無限邊坡模型，得到了前期降雨對(duì)邊坡失穩(wěn)的作用規(guī)律。蔡征龍等[18]基于簡(jiǎn)化Bishop積分法，對(duì)降雨作用下的滑坡失效概率模型進(jìn)行了研究；王雪冰等[19]專門研究淺層土質(zhì)滑坡在強(qiáng)降雨作用下的穩(wěn)定性分析的改進(jìn)模型；譚洋洋等[20]針對(duì)強(qiáng)降雨作用下的電塔邊坡穩(wěn)定性進(jìn)行深入研究。汪丁建等[21]對(duì)強(qiáng)降雨過程中邊坡穩(wěn)定性變化規(guī)律開展了數(shù)值模擬研究;吳云等[22]在降雨過程影響下采用Geo-studio軟件對(duì)滑坡穩(wěn)定性系數(shù)進(jìn)行計(jì)算。

另外，已有研究表明連續(xù)5 d降雨與滑坡的發(fā)生有著密切的關(guān)系。唐棟等[23]研究了土水特征曲線中當(dāng)土體飽和滲透系數(shù)保持不變，而其他參數(shù)變化的情況下不同連續(xù)5 d降雨對(duì)滑坡穩(wěn)定性系數(shù)變化規(guī)律的影響；張玉成等[24]在分析滑坡地質(zhì)資料和氣象資料的基礎(chǔ)上深入探討了滑坡與降雨的分布特點(diǎn)、滑坡與降雨頻次之間的關(guān)系、滑坡與降雨在時(shí)間上的關(guān)系等。因此本文也采用連續(xù)5 d降雨進(jìn)行研究。

大部分學(xué)者對(duì)累積降雨量及不同雨型工況下的滑坡穩(wěn)定性變化規(guī)律的研究不夠深入全面，且沒有統(tǒng)計(jì)出累積降雨量及雨型的發(fā)生概率。因此擬以貴州省思南縣官寨滑坡為例，在思南縣近57 a內(nèi)的日降雨數(shù)據(jù)基礎(chǔ)上，運(yùn)用泊松分布統(tǒng)計(jì)出不同累積降雨量和不同雨型的發(fā)生概率，再結(jié)合飽和-非飽和滲流理論和非線性有限元分析法，揭示滑坡在不同累積降雨量和雨型工況下的地下滲流場(chǎng)和穩(wěn)定性的變化規(guī)律。

1 研究方法

主要對(duì)思南縣官寨滑坡在降雨作用下的穩(wěn)定性變化機(jī)制進(jìn)行研究。首先依據(jù)思南縣歷史降雨資料，統(tǒng)計(jì)近50 a來思南縣發(fā)生的連續(xù)5 d降雨過程的次數(shù)；再用泊松分布統(tǒng)計(jì)連續(xù)5 d降雨產(chǎn)生的累積降雨量的概率分布特征，并計(jì)算不同累積降雨量工況下的各雨型的發(fā)生概率，以此建立降雨工況的組合模式。之后在Geo-Studio數(shù)值模擬軟件中建立官寨滑坡的物理模型，并將各降雨工況添加到該物理模型上；最后采用飽和-非飽和滲流理論計(jì)算得到官寨滑坡的地下水滲流規(guī)律并采用極限平衡法計(jì)算其穩(wěn)定性變化規(guī)律。具體研究思路如圖1所示。

圖1 降雨工況下官寨滑坡穩(wěn)定性變化研究思路Fig.1 Research approach of stability changes rule of Guanzhai landslide under rainfall conditions

1.1 連續(xù)5 d降雨概率統(tǒng)計(jì)

通過對(duì)思南縣57 a(1960—2017年)日降雨量數(shù)據(jù)進(jìn)行分析，統(tǒng)計(jì)連續(xù)5 d降雨的次數(shù)。再根據(jù)降雨強(qiáng)度變化(事件j)[25]將其劃分為上升型(j=A)、下降型(j=B)、先升后降型(j=C)、先降后升型(j=D)等4種雨型。其中各雨型統(tǒng)計(jì)概率分別為PA=0.18、PB=0.23、PC=0.42、PD=0.17。之后，將其所有雨型類別的連續(xù)5 d累積降雨量劃分為20個(gè)等級(jí)，每個(gè)等級(jí)的雨量相差10 mm，其中累積降雨量為10 mm以內(nèi)的忽略不計(jì)，如雨量為10.1～20 mm的記為P10～20。根據(jù)泊松分布定義，將隨機(jī)變量(X)的所有取值定義為0,1,2，…，且得到每個(gè)取值的具體概率值分別為

(1)

式(1)中：K=0，1,2，…；λ>0設(shè)定為常量，則認(rèn)為X服從泊松分布且其參數(shù)是λ，記為X～π(λ)。每個(gè)等級(jí)雨量所占總雨量的發(fā)生概率服從泊松分布，其中λ為各年平均各等級(jí)雨量所占總雨量的頻率，由此可以求出一次降雨落在某個(gè)等級(jí)的概率，記為Pi={i=20,30，…，210}，結(jié)果如圖2、表1所示，其中根據(jù)泊松分布得到的預(yù)測(cè)值與原真實(shí)值均方誤差為

(2)

圖2 不同累積降雨量的概率分布Fig.2 Probability distribution of different cumulative rainfall

雨量iP70～80P90～100P120～130P140～150P170～180P200～210概率Pj0.0310.0280.0180.0110.0030.004

1.2 滑坡飽和非飽和滲流理論

隨著降雨的入滲，滑坡坡體的地下水位線也會(huì)隨之變化，水位線之下形成飽和區(qū)，反之水位線之上形成非飽和區(qū)。滑坡非飽和區(qū)與飽和區(qū)的巖土體中的地下水處于相互連通的狀態(tài)，這種相互連通的運(yùn)動(dòng)特征將引起坡體地下水位的變化，也就是飽和與非飽和滲流場(chǎng)需要解決的問題[20-22]。基于飽和與非飽和滲流場(chǎng)的地下水滲流各向異性的特征，可以根據(jù)質(zhì)量守恒及達(dá)西定律計(jì)算出二維平面中的飽和與非飽和滲流控制表達(dá)式為

(3)

式(3)中:kx、ky為x、y方向上的滲透系數(shù) ，m/s；γw為水的容重 ，106 g/m3；h為水頭，m；h=u/γw+z，u為孔隙水壓力，kPa，z為水頭邊界所在的位置，m；g是地球引力常數(shù)，N/kg；Q是滑坡表面流量的邊界特征；mw定義為比水容重；θw是坡體單位體積內(nèi)的含水率數(shù)值，定義為對(duì)基質(zhì)吸力函數(shù)求偏導(dǎo)后得到的負(fù)值。

(4)

式(4)中：mw定義為對(duì)土水特征曲線的斜率取絕對(duì)值，μw和μa分別表達(dá)為孔隙中的水壓力和氣壓力，kPa。通過計(jì)算或試驗(yàn)得出坡體土水特征值并繪制相關(guān)的滲透曲線，然后將用于模型計(jì)算的初始和邊界等條件與Geo-Studio軟件中的滲流控制函數(shù)相耦合，最終可計(jì)算出坡體地下水的暫態(tài)浸潤(rùn)或滲流特征。將坡體初始特征定義為

h(x,y,0)=h0(x,y,0)

(5)

水頭邊界條件為

(6)

流量邊界條件為

(7)

式中：k為滲透系數(shù)；n為單位邊界法向量。

1.3 穩(wěn)定性計(jì)算理論

降雨過程會(huì)影響巖土體的非飽水和飽水區(qū)域互相發(fā)生變化，因此，巖土體安全系數(shù)計(jì)算不但需基于飽水區(qū)地下水影響，也要基于非飽水區(qū)的基質(zhì)吸力的作用。現(xiàn)利用Fredlund提出的非飽和土中考慮負(fù)孔隙水壓力的抗剪強(qiáng)度計(jì)算方案，對(duì)邊坡穩(wěn)定性進(jìn)行計(jì)算[25-27]，即

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

(8)

式(8)中：σ為法向總應(yīng)力，N；c′、φ′為土體有效應(yīng)力參數(shù)；ua為孔隙氣壓力(假定為常規(guī)大氣壓力, kPa)；ua-uw定義為基底吸力值，N；uw定義為巖體的孔隙水壓力值，kPa；φb定義為常量，表達(dá)抗剪強(qiáng)度值在基底吸力增加的工況下表現(xiàn)出的遞增規(guī)律。

2 官寨滑坡及其地質(zhì)模型

2.1 官寨滑坡概述

思南縣位于中亞熱帶氣候區(qū)，近17 a平均降水量為1 162.81 mm，最大年降水量為1998年的1 529.45 mm，最大日降水量為237 mm(1998年7月22日)，時(shí)空差異明顯，地域性強(qiáng)，降水主要集中在5—8月。思南縣官寨滑坡位于官寨村斜坡地帶，地貌類型以斜坡侵蝕和剝蝕等為主，滑坡區(qū)總體地形坡度為20°～40°，呈階梯狀。整體地勢(shì)西高東低，滑坡后緣最大高程為800.00 m，前緣剪出口最低高程為667.00 m，該滑坡體縱長(zhǎng)395 m，前緣寬208 m，體積約80萬m3，屬于中型中層堆積層滑坡。從地形上看，該滑坡所在的位置受兩側(cè)小山脊的控制，地勢(shì)比較低洼，滑坡前緣出水點(diǎn)、濕地出露較多，從西向東、從上至下呈“溝槽”狀,易于匯集地表水。官寨滑坡的平剖面圖如圖3和圖4所示。

該滑坡滑體為斜坡后緣崩落的大量灰?guī)r塊體、泥巖塊體及部分粉質(zhì)黏土組成，厚度變化較大，厚2～20 m，其中兩側(cè)緣較薄，中軸線一帶厚度較大。該滑坡巖土結(jié)構(gòu)松散。根據(jù)鉆探取樣揭露及現(xiàn)場(chǎng)調(diào)查，并考慮到滑體中灰?guī)r、泥巖塊體的含量約占50%的實(shí)際情況，滑坡總體重度γ應(yīng)有所增大，本次研究綜合選用的數(shù)值模擬參數(shù)值為：天然狀態(tài)下γ取21.0 kN/m3，飽水狀態(tài)下γ′取22.0 kN/m3，滑坡滑床為志留系秀山組泥巖，呈強(qiáng)-中風(fēng)化，巖層產(chǎn)狀220°∠10°。

官寨滑坡最早的變形始于1995年，當(dāng)時(shí)的變形主要集中在滑坡的中后緣。此后，每年都有一定的變形跡象，特別是每年的雨季，變形更加明顯。逐步擴(kuò)展到了滑坡的全范圍，在斜坡上的耕地出現(xiàn)了明顯的分級(jí)錯(cuò)落、建筑物室外地坪出現(xiàn)大量開裂、房屋出現(xiàn)傾斜等現(xiàn)象。可見，該滑坡的穩(wěn)定性情況堪憂，當(dāng)?shù)鼐用癜踩咽艿絿?yán)重威脅。

2.2 滑坡地質(zhì)模型

根據(jù)官寨滑坡工程地質(zhì)特征，確定最危險(xiǎn)剖面作為滑坡滲流場(chǎng)和穩(wěn)定性計(jì)算的地質(zhì)剖面。由于滑帶比較薄，建立地質(zhì)模型設(shè)定滑帶與滑坡體為同一材料。組成該滑坡的巖土體結(jié)構(gòu)為斜坡后緣崩落的大量灰?guī)r塊體、泥巖塊體及部分粉質(zhì)黏土組成，厚度變化較大，大約厚2～20 m。二維有限網(wǎng)格劃分共剖分為493個(gè)單元，619個(gè)節(jié)點(diǎn)，如圖5所示，滑體下面滑床為基巖，定義為不透水層根據(jù)勘察資料和室內(nèi)實(shí)驗(yàn)，確定官寨滑坡的滑體和滑床參數(shù)取值如表2所示，采用SEEP/W模塊進(jìn)行滑坡浸潤(rùn)性變化模擬時(shí)，視滑坡巖土體材料具有非飽和土的性質(zhì)，其巖土體積含水率及滲透系數(shù)設(shè)定為巖土體孔隙水壓力的函數(shù)。再利用Fredlund and Xing 模型方程對(duì)數(shù)據(jù)點(diǎn)進(jìn)行曲線擬合確定這種函數(shù)關(guān)系[17]。采用SEEP/W模塊中的Van Genuchten函數(shù)曲線和土體飽和時(shí)的參數(shù)，來確定滑體和滑床在非飽和特征下的滲透系數(shù)以及體積含水量與基質(zhì)吸力之間的關(guān)系，如圖6所示。

圖3 官寨滑坡勘察平面圖Fig.3 Topographical map of the Guanzhai landslide

圖4 官寨滑坡剖面圖Fig.4 Geological cross-section of Guanzhai landslide

圖5 官寨滑坡數(shù)值模擬的主剖面Fig.5 Main sliding section of Guanzhai landslide for numerical modeling

圖6 滑體水土特征曲線和滲透函數(shù)曲線圖Fig.6 Soil-water characteristic curve and permeability function of the slide mass

部位土體重度γ/(kN·m-3)黏聚力c/kPa內(nèi)摩擦角φ/(°)飽和體積含水率/(m3·m-3)滲透系數(shù)K/(m·s-1)滑體21.523.212.280.453.5×10-5

3 各工況下滑坡滲流場(chǎng)與穩(wěn)定性分析

3.1 滑坡計(jì)算工況及邊界條件

根據(jù)官寨滑坡當(dāng)?shù)氐慕涤昵闆r，現(xiàn)將降雨以8 d為一個(gè)計(jì)算階段，其中前5 d為連續(xù)降雨，后3 d停雨。通過設(shè)計(jì)5 d持續(xù)降雨的雨量和相關(guān)雨型，來探究降雨對(duì)坡體地下水浸潤(rùn)線及安全系數(shù)的作用特征。研究中依次設(shè)定了80、130、180 mm 3種累積雨量，且每種雨量對(duì)應(yīng)著上升型、下降型、先升后降型、先降后升型共4種降雨類型(表3)。各個(gè)工況下模型計(jì)算邊界條件為：本文模型采用Geo-studio中SEEP模塊中的瞬態(tài)分析，整個(gè)滑坡表面設(shè)定為水體入滲邊界。在雨強(qiáng)小于坡體入滲的強(qiáng)度時(shí)，可取滑坡表面為流量邊界，當(dāng)降雨強(qiáng)度大于入滲強(qiáng)度時(shí)，多余雨量隨坡體流下。

3.2 不同雨型及累積降雨量的滲流場(chǎng)

官寨滑坡在不同工況下地下滲流場(chǎng)如圖7～圖9所示。圖7表示各種雨型在工況為1-1、2-1、3-1、4-1下的地下浸潤(rùn)線分布，圖8表示為1-2、2-2、3-2、4-2工況下的地下浸潤(rùn)線分布，圖9表示各種雨型在工況為1-3、2-3、3-3、4-3下的地下浸潤(rùn)線分布。由圖7～圖9可知，隨著降雨的持續(xù)，雨水持續(xù)入滲坡體內(nèi)部引起地下水浸潤(rùn)線位置持續(xù)升高；且當(dāng)?shù)谝惶旖涤炅枯^大時(shí)其地下浸潤(rùn)線初始位置也相對(duì)較高。

表3 用于穩(wěn)定性和滲流計(jì)算的各種工況Table 3 Different combinations of operating conditions for landslide stability calculations

從圖7～圖9中可以看出，對(duì)于上升型降雨其浸潤(rùn)線位置隨著降雨的持續(xù)而均勻升高；對(duì)于下降型降雨其浸潤(rùn)線位置隨著降雨強(qiáng)度的遞減其上升的趨勢(shì)有著明顯的減弱；對(duì)于先升后降型降雨，其浸潤(rùn)線位置變化與降雨強(qiáng)度的變化規(guī)律性類似，呈現(xiàn)出浸潤(rùn)線前3 d大幅上升后2 d變化較小的特征；對(duì)于先降后升型降雨，其浸潤(rùn)線位置在前3 d上升較慢而在后2 d上升加快。一般而言，在第5天停止降雨后，坡體浸潤(rùn)線在第6天有一定幅度的下降，而在第7、第8天只有較小幅度下降。

3.3 不同雨型及降雨量穩(wěn)定性分析

使用M-P法計(jì)算官寨滑坡穩(wěn)定性，M-P法是基于非飽和土的Mohr-Coulomb強(qiáng)度規(guī)則而構(gòu)建的，它首先在Seep/W模塊中計(jì)算滑坡在不同工況條件下的飽和與非飽和滲流場(chǎng)，然后將滲流場(chǎng)數(shù)值模擬得到的滑坡暫態(tài)孔隙水壓力導(dǎo)入Slope/W模塊中，再根據(jù)地勘資料確定滑動(dòng)面，最后Slope/W模塊在地下水滲流場(chǎng)模擬的基礎(chǔ)上計(jì)算出穩(wěn)定性系數(shù)。由此，把動(dòng)態(tài)的滲流場(chǎng)變化過程與穩(wěn)定性計(jì)算相結(jié)合，計(jì)算并繪制出不同累積降雨量和雨型工況下滑坡穩(wěn)定性系數(shù)變化規(guī)律。官寨滑坡在不同的連續(xù)5 d累積降雨量和不同雨型工況下的穩(wěn)系數(shù)變化如圖10所示。

圖7 連續(xù)5 d累積降雨80 mm各種雨型下浸潤(rùn)線分布Fig.7 Transient changes under operating conditions of cumulative rainfall of 80 mm

圖8 連續(xù)累積降雨為130 mm各種雨型浸潤(rùn)線分布Fig.8 Transient changes under operating conditions of cumulative rainfall of 130 mm

圖9 連續(xù)累積降雨為180 mm各種雨型浸潤(rùn)線分布Fig.9 Transient changes under operating conditions of cumulative rainfall of 180 mm

圖10 不同雨型不同累積降雨下穩(wěn)定系數(shù)變化Fig.10 Stability coefficient changes rule under different rainfall pattern and cumulative rainfall

由圖10可知：當(dāng)累積降雨量為80 mm時(shí)其最低安全系數(shù)為1.84；當(dāng)累積降雨量為130 mm時(shí)其最低安全系數(shù)為1.82；當(dāng)累積降雨量為180 mm時(shí)其最低安全系數(shù)為1.80。說明當(dāng)累積降雨量依次增加時(shí)其最低穩(wěn)定系數(shù)也在逐漸減小。這是因?yàn)殚_始時(shí)土體處于非飽和狀態(tài)，隨著累積降雨量的增加，土體的非飽和度在降低，c、φ值減小導(dǎo)致抗剪強(qiáng)度降低,進(jìn)而引起穩(wěn)定系數(shù)減小。其中，上升型降雨作用下滑坡前4 d穩(wěn)定性系數(shù)均比其他3類要大，下降型降雨作用下滑坡前4 d穩(wěn)定系數(shù)均比其他3類要小，先升后降型和先降后升型降雨作用下的滑坡穩(wěn)定性系數(shù)介于上升型和下降型降雨之間。這是因?yàn)樯仙徒涤瓿跏冀涤炅亢苄?dǎo)致初始穩(wěn)定系數(shù)較大，而下降型初始雨量很大故其穩(wěn)定系數(shù)較小。另外兩種降雨則介于兩者之間。

另外，累積降雨量一定時(shí)各種雨型在第5天時(shí)穩(wěn)定性系數(shù)均較為接近，這可能是因?yàn)楦饔晷颓? d累積降雨量都相同。從停雨后第1天開始，滑坡穩(wěn)定系數(shù)略微下降，之后2 d穩(wěn)定系數(shù)保持平穩(wěn)。因?yàn)橥Ｓ旰笃麦w內(nèi)滲流還沒有穩(wěn)定會(huì)導(dǎo)致穩(wěn)定系數(shù)繼續(xù)下降，當(dāng)滲流穩(wěn)定后其穩(wěn)定系數(shù)也保持平穩(wěn)。

在累積降雨量相同的工況下，隨著降雨量的增加上升型降雨曲線斜率越來越大，說明其穩(wěn)定系數(shù)變化率越來越大且在第5天達(dá)到最大；下降型降雨曲線斜率越來越小在第5天達(dá)到最小；另外兩種雨型變化率與其雨型保持一致。

綜合上述分析：連續(xù)5 d累積降雨量一定時(shí)，不同雨型對(duì)滑坡穩(wěn)定性系數(shù)有著不同的影響，前期累積降雨量和滑坡穩(wěn)定系數(shù)有著緊密的關(guān)系。前期累積降雨量越大各雨型的穩(wěn)定性系數(shù)也就越低，因?yàn)槌掷m(xù)降雨入滲作用下的坡體含水率會(huì)逐漸變大，濕潤(rùn)鋒不斷下移，降雨使得滑體重度持續(xù)增大且土體基質(zhì)吸力逐漸降低，導(dǎo)致滑體的抗剪強(qiáng)度減小，增大了其發(fā)生滑坡的風(fēng)險(xiǎn)。

4 結(jié)論

通過分析貴州省思南縣54 a降雨資料，采用泊松分布法來確定某雨型的連續(xù)5 d累積降雨發(fā)生的概率，并結(jié)合數(shù)值模擬來分析思南縣官寨滑坡在連續(xù)5 d降雨作用下的滑坡滲流場(chǎng)和穩(wěn)定性變化規(guī)律，得出以下結(jié)論。

(1)連續(xù)5 d累積降雨量的發(fā)生概率從0.155 (20 mm)逐步下降至0.004 (210 mm)；且在所有雨型中發(fā)生概率最大的為先升后降型降雨，概率最小的為先降后升型降雨。

(2)隨著降雨的持續(xù)，滑坡地下浸潤(rùn)線不斷地上升；對(duì)于上升型降雨其浸潤(rùn)線位置持續(xù)均勻地升高，下降型降雨其浸潤(rùn)線位置的上升的趨勢(shì)明顯減弱，對(duì)于先升后降型降雨其浸潤(rùn)線位置呈現(xiàn)出前3 d大幅上升后2 d逐步減弱的現(xiàn)象，而先降后升型降雨作用下其浸潤(rùn)線位置有著前3 d上升較慢后2 d上升加快的趨勢(shì)。

(3)上升型降雨作用下滑坡前4 d穩(wěn)定性系數(shù)均比其他3類雨型作用下的要大，下降型降雨作用下滑坡前4 d穩(wěn)定系數(shù)均比其他日類雨型作用下的要小。

(4)穩(wěn)定性變化率方面，上升型降雨作用下滑坡穩(wěn)定性變化率在逐漸上升，下降型降雨作用下滑坡穩(wěn)定性變化率在逐漸下降，其他兩種雨型作用下的滑坡穩(wěn)定性變化率介于兩者之間。