強降雨作用下昔格達邊坡滲流特性及穩(wěn)定性分析
2020-09-02 06:41:35楊世豪蘇立君張崇磊李丞胡兵立
土木建筑與環(huán)境工程 2020年4期
楊世豪 蘇立君 張崇磊 李丞 胡兵立
摘要：昔格達地層分布于中國西南，其工程性質(zhì)極差。為分析上覆第四紀殘坡積物的昔格達組粉砂土邊坡在強降雨作用下的滲流特性及穩(wěn)定性，結(jié)合石棉縣莫家崗滑坡，用數(shù)值模擬方法研究坡體滲流規(guī)律;用MorgensternPrice方法計算各時刻坡體穩(wěn)定性系數(shù)，并分析其變化規(guī)律。結(jié)果表明：雨水在坡體中逐層入滲，雨停后滲流過程將繼續(xù)發(fā)展，由于存在土層分界面，坡體滲流特性不同于均質(zhì)坡體;降雨過程中，孔隙水壓力發(fā)生增長的范圍在界面附近加速拓展，雨停后，土層分界面附近第四紀殘坡積物仍處于近飽和狀態(tài);降雨入滲導(dǎo)致坡體邊坡穩(wěn)定性系數(shù)下降，且穩(wěn)定性系數(shù)變化對降雨的響應(yīng)存在滯后性。
關(guān)鍵詞：昔格達組;層狀邊坡;降雨入滲;邊坡穩(wěn)定
中圖分類號：P642.3 文獻標志碼：A 文章編號：20966717（2020）04001909
收稿日期：20191127
基金項目：國家自然科學基金（41761144077）
作者簡介：楊世豪（1995 ），主要從事地質(zhì)災(zāi)害啟動機理研究，Email：347295340@qq.com。
蘇立君（通信作者），男，研究員，博士生導(dǎo)師，Email：sulijun1976@163.com。
Received：20191127
Foundation item：National Natural Science Foundation of China （No. 41761144077）
Author brief：Yang Shihao （1995 ）， main research interest： slope reliability， Email：347295340@qq.com.
Su Lijun （corresponding author）， professor， doctoral supervisor， Email： sulijun1976@163.com.
Analysis of seepage characteristics and stability of Xigeda Formation slope under heavy rainfall
Yang Shihao1，3，Su Lijun1，2，3，Zhang Chonglei1，2，Li Cheng1，3，Hu Bingli1，3
（1.Key Laboratory of Mountain Hazards and Earth Surface Process; Institute of Mountain Hazards and Environment，Chinese Academy of Sciences， Chengdu 610041， P.R. China; 2. Center for Excellence in Tibetan Plateau Earth Sciences， Chinese Academy of Sciences， Beijing 100101， P.R. China; 3. University of Chinese Academy of Sciences， Beijing 100049， P.R. China）
Abstract： Xigeda Formation is distributed in southwest of China， and its engineering properties are very poor. In order to analyze seepage characteristics and stability of silty sand slope of Xigeda Formation overlying Quaternary residual slope sediments under heavy rainfall conditions， seepage law of slope is studied by numerical simulation based on Mojiagang landslide in Shimian County; MorgensternPrice method is used to calculate safety factor of landslide and to analyze their variation laws. The results show that： rainwater infiltrates into slope layer by layer， and seepage process will continue after rainfall. As the existence of soil layer interface， seepage characteristics are different from homogeneous slope; during the rainfall process， the range of pore water pressure growth accelerates near the interface， and the Quaternary residual slope sediments near the interface is still in a nearsaturated state after the rainfall; rainfall infiltration leads to reduction of safety factor of the slope， and stability has a lag in response to rainfall.
Keywords：Xigeda Formation; layered soil slope; rainfall infiltration; slope stability
昔格達地層為廣泛分布于中國西南地區(qū)的半成巖河湖相沉積地層，形成于晚更新世與上第三系之間，主要由粉砂土和黏土組成[12]。隨著工程活動在昔格達地層分布區(qū)廣泛開展，其工程性質(zhì)引起學者關(guān)注。周平等[3]、王志杰等[4]研究發(fā)現(xiàn)，昔格達地層巖土體水穩(wěn)性差，極易崩解。杜翔宇等[5]研究昔格達地層巖土體的微觀結(jié)構(gòu)特性，并對其強度指標與含水率關(guān)系進行闡述。因力學性質(zhì)特殊，昔格達地層分布區(qū)滑坡頻發(fā)，相關(guān)學者對其進行了研究。黃紹檳等[6]將昔格達滑坡歸為覆蓋層滑坡、昔格達基底滑坡和昔格達組層面滑坡3類。丁文富等[7]提出，昔格達地層與上覆松散堆積層的接觸界面、與下伏地層的界面以及層間軟弱沉積結(jié)構(gòu)面是該地層中主要的不良結(jié)構(gòu)面。綜上，昔格達地層邊坡多為層狀結(jié)構(gòu)且滑面多位于不同土層接觸面上。
降雨是常見滑坡誘發(fā)因素[8]，對昔格達地層巖土體遇水軟化的特殊土體而言，這種作用尤為顯著。降雨對邊坡穩(wěn)定性的主要影響在于土體孔隙水壓力增加造成抗剪強度降低[9]。此外，雨水入滲將在邊坡表層巖土體形成暫態(tài)飽和區(qū)，使該區(qū)域土體自重增加，下滑力增大，邊坡穩(wěn)定性降低[10]。White等[11]的研究采用降雨條件下土體體積含水率、孔隙水壓力的變化規(guī)律來表征邊坡的滲流特性。鑒于降雨入滲會對邊坡穩(wěn)定性造成嚴重影響，相關(guān)研究建立了邊坡滲流特性與穩(wěn)定性的關(guān)系[12]。張建等[13]對降雨觸發(fā)淺層坡體失穩(wěn)的遲滯現(xiàn)象及其與土質(zhì)參數(shù)的關(guān)聯(lián)性進行研究。層狀邊坡滲流規(guī)律較均質(zhì)邊坡更復(fù)雜。韓同春等[14]認為，濕潤鋒至接觸面時，引起接觸面孔隙水壓力上升是滑坡主要誘因。馬吉倩等[15]通過數(shù)值模擬方法，得到降雨條件及坡積土層厚度對含水率、孔隙水壓力沿高程分布的影響規(guī)律。石振明等[16]通過改進GreenAmpt模型，提出考慮降雨入滲的多層非飽和土邊坡穩(wěn)定性分析方法。上述研究闡述了層狀邊坡降雨滲流特征與穩(wěn)定性變化規(guī)律。
昔格達地層邊坡多為層狀邊坡，但目前針對降雨入滲作用下昔格達地層邊坡內(nèi)部滲流場分布規(guī)律研究較少。以四川省雅安市石棉縣莫家崗滑坡為例，利用有限元方法模擬經(jīng)歷降雨邊坡內(nèi)滲流過程，并以MorgensternPrice方法計算邊坡各時刻穩(wěn)定性系數(shù)。通過研究坡內(nèi)滲流特點及邊坡在降雨作用下穩(wěn)定性系數(shù)的變化規(guī)律，為昔格達地層滑坡治理提供參考。
1計算理論
降雨入滲實質(zhì)為坡內(nèi)土體含水率上升，由非飽和狀態(tài)向飽和狀態(tài)發(fā)展的過程。入滲過程中，含水率隨深度分布可分為4部分，即：淺層飽和帶、含水率變化較大的過渡帶、含水率均勻分布的傳導(dǎo)帶和濕潤度隨深度減小的濕潤區(qū)，濕潤區(qū)前緣為濕潤鋒[17]。
坡體中，水的運動基本方程為根據(jù)達西定律及質(zhì)量守恒定律推導(dǎo)得到的RICHARDS方程[18]。
式中：x為水平方向;z為豎直方向;hm為基質(zhì)吸力水頭;k（hm）為非飽和土滲透系數(shù)函數(shù);C（hm）為比水容量及水土特征曲線斜率;t為時間。其中，非飽和土滲透系數(shù)k（hm）有別于飽和土，達西定律中，k（hm）等于常數(shù)滲透系數(shù)，為與體積含水率相關(guān)的函數(shù)，反映非飽和土中的水力關(guān)系對滲流的影響。
MorgensternPrice法可對任意形狀滑動面進行求解，同時，滿足豎向力、水平力及力矩的平衡，是一種常用極限平衡求解方法?；居嬎阍頌?，假定兩相鄰?fù)翖l的法向條間力和切向條間力之間存在一對水平方向坐標的函數(shù)關(guān)系，根據(jù)整個滑動土體邊界條件進行迭代，求出問題的解[7]。作為典型的復(fù)雜條分法，其計算過程比一般條分法復(fù)雜，但計算結(jié)果更為準確，因而，適用于利用數(shù)值模擬軟件進行求解。
2計算模型及工況設(shè)計
2.1工程概況
莫家崗滑坡位于石棉縣美羅鄉(xiāng)獅子村4組，后緣地理坐標：N29°17′07.9″，E102°26′52.8″，主滑方向為184°，剖面呈上陡下緩形態(tài)，主滑方向坡高80 m，總體坡度16°，為中型淺層土質(zhì)滑坡。坡體為典型層狀結(jié)構(gòu)，上部為第四系全新統(tǒng)殘坡積物（Q4el+dl），黃褐色粉質(zhì)粘土，局部夾碎塊石，厚度1.0～4.0 m;下部為第三系中統(tǒng)昔格達組粉砂土，強度較差;基巖為第三系中統(tǒng)昔格達組粉砂巖。
2.2計算參數(shù)及邊界條件
由試驗測得，滑坡中第四紀殘坡積物及昔格達組粉砂土相關(guān)參數(shù)見表1。殘坡積物水土特征曲線由軟件內(nèi)置樣本函數(shù)擬合;昔格達粉砂土水土特征曲線由粒徑級配曲線擬合。兩類土滲透系數(shù)函數(shù)根據(jù)相應(yīng)土壤水分特征曲線，由VG模型擬合[19]，昔格達組粉砂土粒徑級配曲線如圖1所示，兩類土水土特征曲線及滲透系數(shù)函數(shù)曲線如圖2所示。
采用Geostudio中的seep/w模塊計算邊坡滲流規(guī)律，geostudio是適用于巖土工程數(shù)值模擬的計算軟件，seep/w為軟件內(nèi)專門用于非飽和土滲流計算的分析模塊。seep/w可進行穩(wěn)態(tài)以及瞬態(tài)滲流計算分析，可以獲得土體內(nèi)部孔隙水壓力、體積含水率等指標的空間及時間變化規(guī)律。此外，seep/w模塊的計算結(jié)果可與其他模塊進行耦合計算，如與SLPOE/W耦合，獲得考慮邊坡內(nèi)部基質(zhì)吸力空間以及時間分布條件下邊坡的穩(wěn)定性系數(shù)。
莫家崗滑坡平面圖利用無人機航拍圖像獲得，如圖3所示，剖面圖如圖4所示，其中，剖面圖參考文獻[20]。坡體地下水位以上網(wǎng)格尺寸為1 m，以下網(wǎng)格尺寸為5 m。坡頂、坡中及坡腳設(shè)置A、B、C三處觀測面。穩(wěn)態(tài)分析中，模型上表面為自由邊界，兩側(cè)及底部為不透水邊界。瞬態(tài)分析包括降雨和不降雨兩段，其中，降雨期間，坡體表面設(shè)為單位流量邊界，流量等于降雨量，雨停后為自由邊界。
結(jié)合石棉縣當?shù)厍闆r及中國氣象局關(guān)于降雨強度的劃分標準，設(shè)定模擬共歷時7 d，前4天為降雨過程，降雨強度分別設(shè)為工況1大暴雨（0.006 m/h）和工況2大雨（0.001 6 m/h）兩種，后3天無降雨作用。
3計算結(jié)果分析
3.1降雨作用下邊坡含水率變化情況
3.1.1坡頂監(jiān)測面
圖5、圖6為兩種工況下模擬過程監(jiān)測面A的含水率分布。在大暴雨工況下，t=48 h，坡體出現(xiàn)飽和區(qū)，隨降雨持續(xù)，飽和區(qū)擴大;t=96 h，最大入滲深度約3 m，飽和區(qū)深度達2 m;雨停后，飽和區(qū)消散，雨水繼續(xù)入滲。t=168 h，入滲深度達3.5 m。A處滲流僅發(fā)生在覆蓋層中，為均質(zhì)土體中滲流。降雨時，含水率沿深度分布曲線含飽和區(qū)、過渡區(qū)、傳導(dǎo)區(qū)、濕潤區(qū)及濕潤鋒，與張建等[13]的結(jié)論一致。雨停后，高程87.5 m以上同深度含水率，隨時間延長而減少;高程87.5 m以下同深度含水率，隨時間延長而增大，濕潤鋒深度增加。在大雨工況條件下，降雨入滲總深度明顯減小，t=96 h時，降雨入滲深度僅為2 m，但與大暴雨工況下的滲流場特征相對比，可以看出，在降雨結(jié)束后，大雨工況條件下，雨水的入滲深度及同深度處的含水率都在繼續(xù)增長。
3.1.2坡中監(jiān)測面含水率
圖7、圖8為兩種工況下模擬過程監(jiān)測面B的含水率分布。在大暴雨工況下，0 3.1.3坡腳監(jiān)測面含水率 圖9、圖10為兩種工況下模擬過程監(jiān)測面C的含水率分布。坡腳覆蓋層厚度1.04 m，t=48 h時，覆蓋層已飽和，且雨水滲入昔格達地層中;t=96 h時，雨水最終入滲至粉砂土層內(nèi)約3 m。因入滲深度范圍內(nèi)兩種土體滲透性質(zhì)不同，含水率分布曲線分層明顯。降雨過程中，覆蓋層迅速飽和，因而，從t=48 h開始，該層土體處于飽和狀態(tài);雨停后，層面附近，殘坡積物含水率約0.35，昔格達組粉砂土含水率約0.17，表明，雨停后層面附近土體含水率較高。在大雨工況下，由于坡腳覆蓋層厚度較薄，降雨期間，雨水已經(jīng)滲透至昔格達地層內(nèi)部，但由于入滲量較小，未形成明確分區(qū)。與大暴雨工況下的滲流場特征進行對比可發(fā)現(xiàn)，在降雨結(jié)束后，水分在坡體內(nèi)的最大入滲深度仍在不斷增加，且同深度處的土體含水率不斷增加。 3.2降雨作用下邊坡孔隙水壓力變化情況 3.2.1坡頂監(jiān)測面孔隙水壓力變化情況 圖11、圖12為兩種工況下模擬過程中監(jiān)測面A孔隙水壓力分布情況。兩種工況下，A處滲流均發(fā)生在覆蓋層中，為均質(zhì)土體中滲流。在大暴雨工況下，雨水入滲深度范圍內(nèi)，孔隙水壓力持續(xù)增長，土體發(fā)生孔隙水壓力增長范圍勻速拓展;由圖11、圖12可知，對于均質(zhì)土體，雨停后，孔隙水壓力隨高程分布曲線可分為兩個段，高程87.5 m以上土體中孔隙水壓力消散;87.5 m以下土體中濕潤鋒深度增加，同深度孔隙水壓力增加。在大雨工況下，水分滲透滲透深度更小，因此，全過程僅發(fā)生在上層覆蓋物中。可以看出，模擬全過程土體中，孔隙水壓力始終在逐漸增大，但未達到正值。 3.2.2坡中監(jiān)測面孔隙水壓力變化情況 圖13、圖14為兩種工況下模擬過程中監(jiān)測面B孔隙水壓力分布。在大暴雨工況下，0 3.2.3坡腳監(jiān)測面孔隙水壓力變化情況 圖15、圖16為兩種工況下模擬過程中監(jiān)測面C孔隙水壓力分布。在大暴雨工況下，因覆蓋層厚度僅1.04 m，t=48 h時，昔格達組粉砂土中出現(xiàn)飽和區(qū)，孔隙水壓力增長至正值，層面附近孔隙水壓力增長范圍加速拓展;由圖15、圖16可知，因兩種土質(zhì)滲流特性差異，雨停后，坡內(nèi)孔隙水壓力隨高程分布曲線存在分區(qū)，分區(qū)特征與坡中監(jiān)測面B類似，覆蓋層中，近坡體表面孔隙水壓力減少明顯;昔格達粉砂土地層中，t=96 h時，昔格達粉砂土中含水率曲線飽和區(qū)與過渡區(qū)交點為同時刻坡內(nèi)孔隙水壓力最大部位，該深度以下同深度孔隙水壓力隨時間增長而增長，濕潤鋒深度加大。在大雨工況下，孔隙水壓力隨降雨過程的發(fā)展而不斷增長，但在t=96 h時，土體內(nèi)部最大孔孔隙水壓力仍未達到正值，降雨結(jié)束后，由于土體中水分繼續(xù)入滲，孔隙水壓力逐漸增長為正值。 3.3昔格達地層層狀結(jié)構(gòu)邊坡滲流特征討論 參考對層狀結(jié)構(gòu)邊坡滲流特性的相關(guān)研究成果，并結(jié)合《工程地質(zhì)手冊》等資料，對普通巖土體相關(guān)參數(shù)的總結(jié)可以發(fā)現(xiàn)，與普通層狀結(jié)構(gòu)土質(zhì)邊坡相比，昔格達地層邊坡中，上層松散堆積物覆蓋層的滲透系數(shù)及飽和含水率的參數(shù)與下層昔格達地層巖土體的滲透系數(shù)及飽和含水率值之間的差距更大。這種層狀結(jié)構(gòu)邊坡中，上下層土體之間在滲透系數(shù)以及飽和含水率等參數(shù)相差懸殊的特點，使得昔格達地層層狀結(jié)構(gòu)邊坡在經(jīng)歷降雨過程時，在土層交界面處土體含水率以及孔隙水壓力的突變較普通層狀結(jié)構(gòu)土質(zhì)邊坡更為劇烈，突變量值更大。 昔格達地層層狀結(jié)構(gòu)邊坡所具有的相較于普通層狀結(jié)構(gòu)土質(zhì)邊坡更為劇烈的土層交界面處，土體含水率以及孔隙水壓力的突變表明，該類邊坡更易在土層交界面附近形成高含水率區(qū)域，且該區(qū)域內(nèi)含水率更高，持續(xù)時間更久，消散更緩慢，這一特點對邊坡穩(wěn)定性造成較大影響。首先，由于接觸面附近持續(xù)處于高含水率狀態(tài)，降雨結(jié)束后，邊坡內(nèi)部仍會持續(xù)發(fā)生雨水滲流過程，使得邊坡安全系數(shù)持續(xù)下降，其滯后性較普通層狀結(jié)構(gòu)土質(zhì)邊坡更為明顯。其次，由于地層接觸面附近的覆蓋層土體含水率較高，且不能較快消散，相較于普通層狀結(jié)構(gòu)土質(zhì)邊坡，昔格達地層層狀結(jié)構(gòu)邊坡土體內(nèi)部水分更易發(fā)生沿土層接觸面的流動，對邊坡內(nèi)部土顆粒產(chǎn)生較大滲流力，甚至對土體顆粒產(chǎn)生淘蝕。 3.4邊坡穩(wěn)定性系數(shù)變化情況 圖17為兩種工況條件下邊坡穩(wěn)定性系數(shù)隨降雨入滲的變化曲線。由圖17可知，兩種工況下，邊坡安全系數(shù)的變化規(guī)律基本一致，降雨期間，穩(wěn)定性系數(shù)持續(xù)下降，下降速率先增大后減小，反映出雨水入滲至坡體不同部位時，穩(wěn)定性系數(shù)變化規(guī)律不同;24 h 3.5邊坡破壞機理分析 根據(jù)相關(guān)研究結(jié)論，莫家崗滑坡為受歷時長的強降雨誘發(fā)的昔格達地層牽引式滑坡，滑面位于基覆界面部位，其變形過程與降雨過程密切相關(guān)。研究區(qū)滑坡最早發(fā)生于2008年汛期，降雨導(dǎo)致斜坡后緣巖土體發(fā)生強烈變形; 2012年雨季斜坡后緣再次產(chǎn)生張拉裂縫;受暴雨影響，2013年7月斜坡發(fā)生進一步強烈變形，形成局部錯臺，斜坡中下部變形十分明顯，附近鄉(xiāng)村公路受路基出現(xiàn)不均勻沉降;2014年7月，研究區(qū)突降暴雨，致使斜坡局部強烈變形;2014年8月，研究區(qū)持續(xù)強降雨，斜坡體多處發(fā)生變形，斜坡后緣產(chǎn)生張拉裂縫[20]。 結(jié)合滲流場數(shù)值模擬結(jié)果，可以將莫家崗滑坡破壞類型概述為兩種，即局部淺層滑塌以及深層牽引型滑坡。莫家崗滑坡為典型的層狀邊坡，雨水在邊坡中逐層滲入，降雨入滲，首先對邊坡上層第四紀殘坡積物覆蓋層產(chǎn)生影響，造成覆蓋層土體處于飽和狀態(tài)，加之表層土體受人類活動改造劇烈，易引起局部淺層滑塌。結(jié)合滲流場模擬結(jié)果以及邊坡結(jié)構(gòu)特點，可以對邊坡深層牽引型滑坡的破壞機理做出解釋：邊坡中下部第四紀殘坡積物覆蓋層厚度較薄，降雨可滲流至下層昔格達組地層中，且由于兩種土體的水土性質(zhì)差異懸殊，造成界面處長期處于高含水率狀態(tài)，由于昔格達地層土體力學性質(zhì)差，在長期受雨水作用下，易造成坡體下部巖土體發(fā)生滑面位于土層交界面的變形破壞，因而，邊坡強烈變形多發(fā)生于坡體中下部。坡體上部覆蓋層土體厚度較大，雨水不易滲入下層昔格達地層中，因此，上部坡體能保持較好完整性，但由于下部坡體力學性質(zhì)不斷減弱，上部坡體抗力減少，整體受牽引下滑，后緣產(chǎn)生拉裂縫。 4結(jié)論 1）昔格達地層邊坡多為層狀結(jié)構(gòu)，該工程為上覆第四紀坡殘積物的昔格達組粉砂土邊坡;降雨過程中，雨水在坡體中逐層入滲;雨停之后，入滲過程仍將持續(xù)，但殘坡積物與昔格達組粉砂土中含水率及孔隙水壓力分布特征有明顯差異。 2）土層交界面處含水率、孔隙水壓力變化明顯;界面上部土體長期處于高含水率狀態(tài)，使昔格達地層靠近界面處土體長期受水分軟化作用，易形成軟弱夾層，誘發(fā)滑面位于界面的滑坡;土體中孔隙水壓力發(fā)生增長的影響范圍在界面附近會加速拓展，使邊坡抗滑力減弱，導(dǎo)致穩(wěn)定性系數(shù)加速降低。 3）降雨入滲導(dǎo)致邊坡穩(wěn)定性降低，當雨水入滲至坡體不同部位時，穩(wěn)定性系數(shù)表現(xiàn)出不同的變化規(guī)律;降雨結(jié)束后，穩(wěn)定性系數(shù)不會立即回升，表明邊坡對降雨入滲的響應(yīng)存在滯后性。參考文獻： [1] 劉惠軍， 聶德新. 昔格達地層研究綜述[J]. 地球科學進展， 2004， 19（Sup1）： 8082. LIU H J， NIE D X. The Overview of the Xigeda strata study [J]. Advance in Earth Sciences， 2004， 19（Sup1）： 8082. （in Chinese） [2] 羅運利， 劉東生. 昔格達組沉積環(huán)境演化與旋回地層學研究[J]. 第四紀研究， 1998， 18（4）： 373. LUO Y L， LIU D S. Cyclostratigraphy of the Xigeda formation and the paleoenvironmental implications [J]. Quaternary Sciences， 1998， 18（4）： 373. （in Chinese） [3] 周平， 王志杰， 徐海巖， 等. 考慮含水率的昔格達地層隧道圍巖穩(wěn)定及亞級分級研究[J]. 土木工程學報， 2017， 50（12）： 97110. ZHOU P， WANG Z J， XU H Y， et al. Stability and subclassification study on the tunnel surrounding rock of Xigeda strata considering the influence of moisture content [J]. China Civil Engineering Journal， 2017， 50（12）： 97110. （in Chinese） [15] 馬吉倩， 付宏淵， 王桂堯， 等. 降雨條件下成層土質(zhì)邊坡的滲流特征[J]. 中南大學學報（自然科學版）， 2018， 49（2）： 464471. MA J Q， FU H Y， WANG G Y， et al. Seepage characteristics of layered soil slope under rainfall conditions [J]. Journal of Central South University（Science and Technology）， 2018， 49（2）： 464471. （in Chinese） [16] 石振明， 沈丹祎， 彭銘， 等. 考慮多層非飽和土降雨入滲的邊坡穩(wěn)定性分析[J]. 水利學報， 2016， 47（8）： 977985. SHI Z M， SHEN D Y， PENG M， et al. Slope stability analysis by considering rainfall infiltration in multilayered unsaturated soils [J]. Journal of Hydraulic Engineering， 2016， 47（8）： 977985. （in Chinese） [17] 雷志棟， 楊詩秀， 謝森傳. 土壤水動力學[M]. 北京： 清華大學出版社， 1988. LEI Z D， YANG S X， XIE S C. Soil Water Dynamics [M]. Beijing： Tsinghua University Press， 1988. [18] RICHARDS L A. Capillary conduction of liquids through porous mediums [J]. Physics， 1931， 1（5）： 318333. [19] VAN GENUCHTEN M T. A closedform equation for predicting the hydraulic conductivity of unsaturated soils [J]. Soil Science Society of America Journal， 1980， 44（5）： 892898. [20] 覃亮， 楊學之， 王海波. 石棉縣某昔格達地層滑坡防治對策探討[J]. 成都大學學報（自然科學版）， 2017， 36（1）： 102105. QIN L， YANG X Z， WANG H B. Discussion on prevention countermeasures of Xigeda strata landslide in Shimian County [J]. Journal of Chengdu University（Natural Science Edition）， 2017， 36（1）： 102105. （in Chinese） （編輯鄧云）
