上部覆蓋層下伏頁巖一類潛在滑坡演化變形機(jī)制分析

李瑞園　鄭博　梁軍林　容洪流

摘 要：為研究上部覆蓋層下伏頁巖這一類型的滑坡失穩(wěn)變形，需要確定其變形的潛在影響因素，從定量角度上分析這一類型滑坡災(zāi)害的變形演化問題。以容縣六王鎮(zhèn)西流河潛在滑坡為研究對象，借助鉆孔勘探、室內(nèi)試驗和地質(zhì)調(diào)繪，詳細(xì)分析了西流河滑坡的工程地質(zhì)概況，并推測其可能的變形演化機(jī)制，然后利用Geostudio數(shù)值軟件建立西流河滑坡的滲流模型，模擬其在不同降雨工況下的滲流規(guī)律。研究發(fā)現(xiàn)：西流河滑坡在持續(xù)降雨過程中穩(wěn)定性系數(shù)從1.224逐步降低至1.09，穩(wěn)定性系數(shù)不斷下降，滑坡體內(nèi)由垂直入滲逐步轉(zhuǎn)化沿滑動帶穩(wěn)定滲流;當(dāng)遭遇極端暴雨情況（50年重現(xiàn)期）時，滑坡體內(nèi)穩(wěn)定滲流場顯著增厚，飽和滑帶土增多，滑坡安全儲備達(dá)到極限。由此可以得出：上部覆蓋層下伏頁巖一類滑坡變形演化過程主要分滑坡體上部覆蓋層堆積階段，垂直入滲至頁巖富水軟化階段和穩(wěn)定滲流至潛在滑動帶軟化蠕變這3個階段;其中上部覆蓋層的高滲透性和下伏頁巖的隔水特征是導(dǎo)致頁巖風(fēng)化帶軟化蠕變和穩(wěn)定滲流場形成的關(guān)鍵因素。研究結(jié)果為上部覆蓋層下伏頁巖這一類滑坡影響因素的研究提供了一些證據(jù)，對研究滑坡演化變形機(jī)制具有借鑒意義。

關(guān)鍵詞：巖土力學(xué);滑坡;廣西容縣;頁巖;滲流模擬;失穩(wěn)機(jī)制

中圖分類號：P642.2

文獻(xiàn)標(biāo)識碼：A

Abstract： In order to study the landslide instability deformation of the shale under the upper overburden， it is necessary to determine the potential influencing factors of the deformation and analyze the deformation evolution of this type of landslide disaster from a quantitative perspective.Taking the potential landslide of Xiliu River in Liuwang Town， Rongxian County as the research object， ?with the aid of drilling exploration， indoor test， geological annotation， the general situation of engineering geology of Xiliu River landslide was analyzed in detail， and the possible mechanism of deformation evolution was speculated. Then， the seepage model was established by using Geostudio numerical software to simulate its seepage law under different rainfall conditions. The results show that the stability coefficient of Xiliu River landslide gradually decreases from 1.224 to 1.09 in the process of continuous rainfall and the landslide body gradually changes from vertical infiltration to the stable seepage flow along the slip zone; When extreme rainstorm is encountered （50year recurrence period）， the stable seepage field in the landslide body thickens significantly， the saturated slide zone soil increases， and the landslide safety reserve reaches the limit. It can be concluded that the evolution process of the landslide deformation of the underlying shale in the upper overburden is mainly divided into three stages： The accumulation stage of the overburden of the landslide body， the vertical infiltration stage to the waterrich softening stage of the shale， and the stable seepage stage to the softening creep stage of the potential sliding zone. The high permeability of the upper overburden and the waterproof characteristics of the underlying shale are the key factors leading to the formation of in the shale weathering zone. The results of this study provide some evidence for influencing factors and evolutionary deformation mechanism of the landslide such as the underlying shale in the upper overburden.

Keywords：geomechanics;landslide; Rongxian county， Guangxi; shale; seepage simulation; instability mechanism

廣西壯族自治區(qū)容縣近年來由于礦產(chǎn)開發(fā)和過度資源開發(fā)導(dǎo)致滑坡、崩塌地質(zhì)災(zāi)害頻發(fā)。據(jù)廣西地質(zhì)環(huán)境監(jiān)測總站的監(jiān)測數(shù)據(jù)統(tǒng)計，容縣地區(qū)20年來共發(fā)生了180余次規(guī)模大小不一的滑坡、崩塌地質(zhì)災(zāi)害，在廣西省各縣市地質(zhì)災(zāi)害中排名靠前，滑坡等地質(zhì)災(zāi)害造成的直接經(jīng)濟(jì)損失超過0.5億元。容縣滑坡地質(zhì)災(zāi)害缺乏系統(tǒng)性分區(qū)評價和預(yù)警系統(tǒng)，滑坡的防災(zāi)減災(zāi)工作依然面臨較為嚴(yán)峻的挑戰(zhàn)。

滲透破壞作用是滑坡滲流穩(wěn)定性分析的重要研究方向，國內(nèi)外相關(guān)學(xué)者對降雨入滲所誘發(fā)的滑坡滲透破壞研究成果較多。LUMB[1]研究了香港地區(qū)的氣象資料，建立了滑坡與降雨的耦合關(guān)系，提出了一種簡化的降雨垂直入滲模型;安關(guān)峰等[2]引入了有限元分析和有限差分法模擬非飽和土在不同降雨條件下的穩(wěn)定性分析規(guī)律;嚴(yán)紹軍等[3]結(jié)合室內(nèi)模擬試驗，模擬了滑坡體在降雨滲流場中隨時間變化的規(guī)律，認(rèn)為水在滑坡體內(nèi)的運(yùn)動分為入滲、水力重分布和排水3個階段。

涉及廣西地區(qū)滑坡災(zāi)害的研究成果主要包括：廖麗萍等[4]搜集了廣西壯族自治區(qū)容縣2004-2015年的滑坡崩塌調(diào)查圖及崩滑統(tǒng)計資料，統(tǒng)計發(fā)現(xiàn)僅容縣地區(qū)現(xiàn)有滑坡總數(shù)為772個，其中高發(fā)密度區(qū)（六王鎮(zhèn)）滑坡約占總數(shù)的73.3%，同時建立了容縣崩塌滑坡地質(zhì)災(zāi)害的易發(fā)性評價指標(biāo);黃希明等[5]以容縣平梨滑坡為研究對象，分析了降雨和天然工況下該類型滑坡的穩(wěn)定性問題，認(rèn)為該滑坡前緣存在應(yīng)力集中和小規(guī)模崩塌的可能;文海濤等[6]和邱恩露等[7]從統(tǒng)計學(xué)角度分析了容縣滑坡的空間分布特征和時間分布特征，認(rèn)為強(qiáng)降雨周期和泥頁巖地層地區(qū)是滑坡高發(fā)頻率地區(qū)。

分析現(xiàn)有研究資料發(fā)現(xiàn)，容縣地區(qū)滑坡分析研究主要集中在數(shù)據(jù)統(tǒng)計、定性分析和現(xiàn)象描述階段[812]，對于特定成因類型的滑坡尤其是高發(fā)密度區(qū)滑坡類型的定量化分析評價成果較少。同時大多數(shù)研究成果是基于土質(zhì)滑坡體的研究，關(guān)于土頁巖復(fù)合地層滑坡的研究相對較少。針對這一問題，筆者以容縣六王鎮(zhèn)典型滑坡體西流河滑坡（上部覆蓋層、下伏頁巖地層）為研究對象，結(jié)合現(xiàn)有的地質(zhì)勘探資料，建立該滑坡的地層結(jié)構(gòu)模型，從地層地質(zhì)、降雨和水文等方面重點(diǎn)剖析了該潛在滑坡的演化變形機(jī)制;同時，結(jié)合試驗獲取滑帶土的非飽和土水特性參數(shù)，開展西流河滑坡的滲流穩(wěn)定性計算，分析這類型滑坡在水巖相互作用下的變形劣化過程，從而驗證其演化變形機(jī)制，最后提出對這類型滑坡的預(yù)防加固措施建議。

1 西流河滑坡概況

西流河滑坡位于廣西壯族自治區(qū)容縣大王鎮(zhèn)，東經(jīng)：112°58′46″，北緯：25°30′37″，滑坡總體形態(tài)呈寬“U”字型，近東西向分布，面積約為2.11萬m2，滑動體厚度為12～15 m，體積為2.4～3.0萬m3，滑坡體滑動方式以推移式為主。根據(jù)現(xiàn)場地質(zhì)測繪，滑坡體現(xiàn)狀條件下存在較為明顯的滑動傾向，所在區(qū)域存在3處較為明顯的陡坎，滑體后緣公路出現(xiàn)較多微裂縫，滑體前緣部分土體堆積隆起，樹木向河道傾斜。依據(jù)目前的調(diào)查結(jié)果推斷滑坡體規(guī)模如下：滑體橫向?qū)挾?62～370 m，縱向長度610～619 m，厚度12～15 m，兩側(cè)地形較為陡峭，后緣以繞山公路為界，前緣為水稻田，潛在滑動方向為SE83°，其全貌如圖1所示。

1.1 地 層

根據(jù)西流河滑坡鉆孔勘探資料，該滑坡體主要地層分布情況如圖2所示。

從圖2可以看出，地層分為4層：第1層為耕植土，灰黑色、黃褐色，主要由粉土、砂巖碎塊石構(gòu)成，潮濕，稍密松散，夾雜大量植物根系，厚1.00～1.50 m;第2層粉質(zhì)黏土（Qel+dl3），黃褐色，中密，稍濕，主要由粉土及少量細(xì)砂構(gòu)成，厚約8.60 m;第3層為頁巖風(fēng)化帶，青黑色、青褐色，泥質(zhì)結(jié)構(gòu)，可塑狀態(tài)，含有較多的粗顆粒，為頁巖風(fēng)化而成，厚0.20～3.00 m，厚度變化不均勻;第4層為頁巖，青黑色、黃褐色，薄層狀結(jié)構(gòu)，大部分呈碎塊狀及薄片狀，其遇水很短時間內(nèi)發(fā)生崩解，耐久性很差，揭露厚度12.70 m，未揭穿。

1.2 地質(zhì)構(gòu)造

研究區(qū)內(nèi)斷裂以繼承性跳躍式運(yùn)動為特點(diǎn)，斷裂延續(xù)時間較長，斷流形成時代劃分為中生代末至新生代初較合適，同時期在第三紀(jì)末至第四紀(jì)初普遍復(fù)活。其中，褶皺、斷裂較為發(fā)育，多數(shù)褶皺構(gòu)造以向斜為主，區(qū)內(nèi)巖體破碎，有助于水的滲透侵蝕，從而劣化巖體。第四紀(jì)以來本區(qū)域斷裂無活動跡象，區(qū)域穩(wěn)定性較好。西流河滑坡位于一個向斜東南角方向，受構(gòu)造裂隙切割導(dǎo)致頁巖整體產(chǎn)狀不均勻，完整性較差，在降雨入滲過程中容易受雨水侵蝕形成軟弱滑動面。

1.3 氣象與水文

研究區(qū)屬于溫帶大陸性氣候，多年平均年降水量為1 329.8 mm，降水量年內(nèi)分配極不均勻，降水主要集中在6—9月份，汛期降水量可占年降水量的75%，是形成本區(qū)雨洪、內(nèi)澇的主要原因。最大年降水量發(fā)生在7月，占全年降水量的23.7%，3—5月的降水量僅占年降水量的14.8%，容易出現(xiàn)春旱。11月至次年3月降水量僅占年降水量的7.6%。降水量年際變化比較明顯，有明顯的豐枯交替性。歷年最大降水量2 124.7 mm（1984年）。50年內(nèi)極端暴雨記錄境內(nèi)24 h最大降雨量為699 mm。

勘察場區(qū)地下水類型為第四系松散巖類孔隙潛水及基巖風(fēng)化帶網(wǎng)狀孔隙裂隙水，補(bǔ)給方式受大氣降水補(bǔ)給，勘察時為豐水期，初見水位埋深2.10～3.60 m，靜止水位埋深3.20～4.10 m，靜止水位標(biāo)高494.92～517.20 m，水位年變幅1.50 m左右。

西流河滑坡所在區(qū)域雨季漫長，巖土體長期飽和，上部覆蓋層與頁巖交界處容易產(chǎn)生地下水富集。旱季高溫導(dǎo)致水體蒸發(fā)，地層變得干燥而脆裂，頁巖具有較強(qiáng)的崩解性，季節(jié)反復(fù)交替作用下滑帶土的力學(xué)性質(zhì)不斷劣化可能形成貫通的滑帶。

1.4 人類活動

人類活動對西流河滑坡的影響主要體現(xiàn)在如下3個方面。

1）植被開采導(dǎo)致的水土流失破壞，植被對于上部覆蓋層下伏頁巖類型滑坡的保護(hù)作用十分明顯。

2）工程改造活動，西流河滑坡體后緣的公路開挖邊坡以及前緣的開荒種地對于滑坡體的穩(wěn)定性均是不利因素。

3）生活用水不合理排放，附近居民用水不規(guī)范排放改變滑坡體原有的水文條件，為地下水增加了來源，將會加速軟化滑帶土。

2 上部覆蓋層下伏頁巖滑坡建模分析

由于西流河滑坡這一類型滑坡體受滲透破壞影響明顯，筆者采用具有成熟的滲流模擬能力的Geostudio軟件建立有限元模型[1314]，用以模擬西流河滑坡在不同降雨工況所引發(fā)的滲流場分布和穩(wěn)定性問題。

2.1 網(wǎng)格劃分及建模參數(shù)

1）建模網(wǎng)格

根據(jù)西流河滑坡實(shí)測斷面CAD圖，利用插值法生成模型尺寸，網(wǎng)格劃分精度為1 m，整個模型網(wǎng)格節(jié)點(diǎn)為1 211個，模型網(wǎng)格數(shù)為1 986個，如圖3所示。

2）地層劃分

滑坡體主要分為3層，第1層為粉質(zhì)黏土及耕植土層，由于耕植土性狀不均一且厚度較薄，筆者考慮將其作為粉質(zhì)黏土層考慮;第2層為頁巖風(fēng)化帶;第3層為基巖頁巖層。

3）降雨滲流場

采用軟件中的Seep/W，Sigma/W，Slope/W模塊進(jìn)行耦合分析，其中降雨工況依據(jù)容縣當(dāng)?shù)氐臍庀筚Y料，筆者將西流河滑坡體受降雨條件影響的強(qiáng)度值做如下規(guī)定。

降雨條件選取年平均降雨量為1 329.8 mm，根據(jù)年平均降雨日為164 d，得出降雨強(qiáng)度平均值為0.39 mm/h。極端暴雨條件，根據(jù)24 h極端降雨量（考慮50年內(nèi)記錄的極端值699 mm），得出極端暴雨條件下降雨強(qiáng)度平均值為29.12 mm/h。降雨過程采用振幅函數(shù)來進(jìn)行循環(huán)定義，其中降雨強(qiáng)度隨時間的變化周期設(shè)定為72 h，在0～24 h內(nèi)從0線性增大至最高平均值，24～48 h內(nèi)保持平均降雨強(qiáng)度，48～72 h內(nèi)降雨強(qiáng)度線性減小至0。

地下水位根據(jù)鉆孔揭露實(shí)際地下水位考慮。

4）物理力學(xué)參數(shù)

模型中各土層計算參數(shù)考慮室內(nèi)土工試驗和現(xiàn)場原位試驗綜合確定，如表1所示。

2.2 滲流分析成果

模擬結(jié)果如圖5所示。

未發(fā)生降雨條件時，如圖5 a）所示，滑坡體內(nèi)沒有滲流發(fā)生，其穩(wěn)定性系數(shù)為1.224，處于穩(wěn)定狀態(tài)。

當(dāng)降雨發(fā)生35 d后，如圖5 b）所示，坡體內(nèi)處于降雨垂直入滲過程，當(dāng)遇水滲透至頁巖層時滲透能力減弱，其上部土體逐漸進(jìn)入飽和狀態(tài)，此時滑坡穩(wěn)定性系數(shù)1.194，穩(wěn)定性有所降低。

降雨發(fā)生90 d后，如圖5 c）所示，降雨入滲過程由坡體垂直入滲轉(zhuǎn)變?yōu)檠仫L(fēng)化帶向下游入滲，形成穩(wěn)定的滲流面，穩(wěn)定性系數(shù)在這一階段大幅度降低至1.09，這一變化過程符合滲流模擬分析。

當(dāng)給予模型施加極端降雨工況后，滑坡體內(nèi)的飽和土層厚度加大，穩(wěn)定的滲流場變厚，這些變化均增加了滑坡的下滑力，根據(jù)自動搜索的滑動面計算出來此時滑坡體的穩(wěn)定性系數(shù)為1.078，處于極限平衡狀態(tài)，隨時有發(fā)生整體滑動的風(fēng)險。

從模型的模擬結(jié)果來看，上覆第四系堆積物下伏頁巖類型滑坡的變形破壞機(jī)制可以初步推斷為上部具有高滲透性的土體在充分吸收降雨或外部水源的垂直補(bǔ)給后，遇到透水性微弱的基巖頁巖層，逐漸在基覆界面富水飽和軟化，并形成穩(wěn)定滲流場;持續(xù)降雨過程中滑體飽和土層增加，滲流場增厚，導(dǎo)致下滑力陡增，但是滑動帶收到軟化蠕變導(dǎo)致抗滑力下降，當(dāng)總的下滑力大于抗滑力時，滑坡發(fā)生破壞。

3 上部覆蓋層下伏頁巖滑坡演化分析

通過對西流河滑坡地形地質(zhì)條件、氣象和水文以及人工活動等因素分析，結(jié)合

數(shù)值模擬成果，認(rèn)為西流河滑坡作為上部覆蓋層下伏頁巖一類滑坡的主要變形演化過程分為3個階段，如圖6所示。

1）上部覆蓋層堆積階段

如圖6 a）所示，該階段主要為滑坡形成過程，由于研究區(qū)地形剛剛形成，上部覆蓋層主要來源于殘坡積物堆積形成，同時頁巖作為基巖并未大規(guī)模風(fēng)化，該階段是滑坡形成的初始階段。

2）接觸帶富水軟化形成風(fēng)化帶

如圖6 b）所示，上部殘坡積物由于以粉質(zhì)黏土夾雜砂、碎塊石，其滲流系數(shù)較高，強(qiáng)降雨作用下地表水大規(guī)模下滲，再遇到滲透性很小的侏羅系頁巖層后下滲水在基覆界面富集，同時在坡內(nèi)形成向坡腳的穩(wěn)定滲流場。這一軟化過程使得頁巖逐步發(fā)生軟化泥化，風(fēng)化帶逐步形成。

3）風(fēng)化帶長期飽水蠕變形成滑動帶

如圖6 c）所示，隨著頁巖風(fēng)化帶內(nèi)土體長期飽水軟化或者干濕循環(huán)交替，在基覆界面逐漸形成貫通的滑動面，極端降雨條件下大流量地表水短時間進(jìn)入滑體并向下游滲流，增加滑體的下滑力;同時，地下水的滲入使得滑帶土飽水抗剪強(qiáng)度降低。當(dāng)滑體的下滑力增加到大于抗滑力時，將會在基巖與覆蓋層接觸面形成貫通的軟化滑動帶，從而導(dǎo)致滑坡發(fā)生。

4 結(jié) 論

通過研究廣西容縣西流河潛在滑坡的變形演化機(jī)制并結(jié)合數(shù)值模擬軟件進(jìn)行驗證，分析了上部覆蓋層下伏頁巖一類滑坡的滲流變形和破壞機(jī)制，得出如下結(jié)論：

1）容縣較廣西其他地區(qū)滑坡發(fā)生頻率和規(guī)模均顯著偏大，這與該地區(qū)的強(qiáng)降雨條件和特殊的地層地質(zhì)條件密不可分。

2）西流河滑坡變形演化主要有3個階段：第1階段為覆蓋層堆積階段，是滑坡體形成初始階段;第2階段為降雨垂直入滲軟化基覆界線，形成頁巖風(fēng)化帶階段;第3階段為持續(xù)降雨在覆蓋層底部形成穩(wěn)定滲流場，導(dǎo)致風(fēng)化帶富水軟化形成滑動帶，當(dāng)抗滑力小于總的下滑力時滑坡變形破壞發(fā)生。

3）數(shù)值模擬結(jié)果表明，在降雨持續(xù)過程中，西流河滑坡的穩(wěn)定性系數(shù)從1.224逐步降低至1.090，滑坡體內(nèi)的滲流變化規(guī)律與推測基本一致，當(dāng)發(fā)生極端降雨條件時，滑坡體處于極限平衡狀態(tài)。

4）針對這一地層類型滑坡體，建議在初期進(jìn)行干預(yù)，阻斷或降低地表垂直降水入滲路徑，做好排水截水措施，防止形成穩(wěn)定滲流場，加強(qiáng)植被保護(hù)，建立完善的滑坡活動監(jiān)測預(yù)警預(yù)報系統(tǒng)。

參考文獻(xiàn)/References：

[1] LUMB P. Effect of Rain Storms on Slope Stability[M].Hong Kong： Local Property & Printing Company， Limited， 1962.

[2] 安關(guān)峰，殷坤龍，唐輝明.黃土坡滑坡的離散元研究[J].地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報，2002，27（4）：464466.

AN Guanfeng， YIN Kunlong， TANG Huiming. Discrete element analysis of Huangtupo landslide[J]. Earth Science——Journal of China University of Geosciences， 2002，27（4）：464466.

[3] 嚴(yán)紹軍， 唐輝明，項偉. 降雨對滑坡穩(wěn)定性影響過程分析[J]. 水文地質(zhì)工程地質(zhì)， 2007（2）：3336.

YAN Shaojun， TANG Huiming， XIANG Wei. Effect of rainfall on the stability of landslides[J]. Hydrogeology & Engineering Geology， 2007（2）： 3336.

[4] 廖麗萍，于淼，文海濤，等.廣西東南部容縣崩塌滑坡的易發(fā)性評價[J].地球與環(huán)境，2019，47（4）：518526.

LIAO Liping， YU Miao， WEN Haitao， et al. Evaluation on the susceptibility of collapse and landslide in Rongxian county， southeastern Guangxi[J]. Earth and Environment，2019，47（4）：518526.

[5] 黃希明，江思義，李海良.桂東南容縣平梨滑坡的成因機(jī)制及穩(wěn)定性評價[J].資源信息與工程，2019，34（6）：112118.

[6] 文海濤，韋朝華，廖麗萍，等.桂東南容縣地質(zhì)災(zāi)害發(fā)育與時空分布特征[J].水土保持通報，2017，37（5）：182188.

WEN Haitao， WEI Chaohua，LIAO Liping， et al. Occurrence and temporalspatial distribution characteristics of geological hazards in Rongxian county in southeast Guangxi Zhuang autonomous region[J]. Bulletin of Soil and Water Conservation，2017，37（5）：182188.

[7] 邱恩露，黃志強(qiáng)，陸仕壹，等.降雨對滑坡穩(wěn)定性影響的分析及模擬——以融水苗族自治縣良寨鄉(xiāng)安全村為例[J].南方國土資源，2019（10）：2324.

[8] 李剛，寧波，楊勃.關(guān)中環(huán)線臨渭段滑坡成因分析與防治對策[J].河北工業(yè)科技，2016，33（1）：9094.

LI Gang， NING Bo， YANG Bo.Analysis and protection of landslide in linwei segment of Guanzhong Loop Highway[J]. Hebei Journal of Industrial Science and Technology，2016，33（1）：9094.

[9] 韓同春， 馬世國， 徐日慶. 強(qiáng)降雨條件下氣壓對滑坡延時效應(yīng)研究[J]. 巖土力學(xué)， 2013，34（5）：13601366.

HAN Tongchun， MA Shiguo， XU Riqing. Research on delayed effect of landslides caused by air pressure under heavy rainfall[J]. Rock and Soil Mechanics， 2013，34（5）：13601366.

[10]溫銘生， 王連俊， 李鐵峰， 等. 基于實(shí)時監(jiān)測的河口糖廠滑坡監(jiān)測系統(tǒng)研究[J]. 工程地質(zhì)學(xué)報， 2011， 19（1）：7682.

WEN Mingsheng，WANG Lianjun，LI Tiefeng，et al.Realtime monitoring based landslide monitoring systemat Hekou sugar refnery factory[J].Journal of Engineering Geology， 2011， 19（1）： 7682.

[11]申春妮，方祥位，王和文，等.吸力、含水率和干密度對重塑非飽和土抗剪強(qiáng)度影響研究[J].巖土力學(xué)，2009， 30（5） ： 13471351.

SHEN Chunni， FANG Xiangwei， WANG Hewen et al. Research on effects of suction， water content and dry density on shear strength of remolded unsatura ted soils[J].Rock and Soil Mechanics， 2009， 30（5）：13471351.

[12]樊玉敬，錢龍，劉碩，等.雷達(dá)崗滑坡形成機(jī)理及穩(wěn)定性分析[J].河北工業(yè)科技，2013，30（2）：6365.

FAN Yujing， QIAN Long， LIU Shuo，et al. Formation mechanism and stability of leidaganglandslide[J]. Hebei Journal of Industrial Science and Technology，2013，30（2）：6365.

[13]李汝成， 王復(fù)明. 降雨入滲對泥巖土混填路堤穩(wěn)定性的影響 [J]. 巖石力學(xué)與工程學(xué) 報 ， 2008， 27（11）：22602266.

LI Rucheng，WANG Fuming. Effect of rainfall infiltration on stability of mudstonesoil mixture embankment[J].Chinese Journal of Rock Mechanics and Engineering，2008， 27（11）： 22602266.

[14]祝艷波，余宏明，高建偉，等.巴東組非飽和紅土強(qiáng)度與變形特性試驗研究[J].工程地質(zhì)學(xué)報，2012，20（6）： 10501056.

ZHU Yanbo， YU Hongming， GAO Jianwei，et al.Experimental study on deformation and strength characteristics of Badong unsaturated red clay[J].Journal of Engineering Geology，2012，20（6）：10501056.

[15]鄭建浩， 張志斌， 袁俊平. 降雨過程殘積土高邊坡土壤吸力動態(tài)監(jiān)測[J]. 中外公路， 2006， 26（2）： 4750.

[16]翟會君，翟亞鋒，朱濤，等.基于回歸ELM神經(jīng)網(wǎng)絡(luò)模型的滑坡變形及失穩(wěn)預(yù)測模型[J].河北工業(yè)科技，2017，34（6）：440447.

ZHAI Huijun，ZHAI Yafeng，ZHU Tao，et al.Prediction model of landslide deformation and instability based on regressionELM neural network model[J].Hebei Journal of Industrial Science and Technology，2017，34（6）：440447.

[17]胡廷飛，王輝，胡傳旺 等. 礫石覆蓋厚度對斥水土壤入滲特性的影響及模型優(yōu)選[J]. 水土保持學(xué)報，2019，33（2）：1729.

HU Tingfei， WANG Hui， HU Chuanwang， et al. Effect of thickness of gravel cover on infiltration characteristics of water repellent soils and its model optimization[J].Journal of Soil and Water Conservation， 2019， 33（2）：1729.