張晉梅 張震 李家春

摘要：為了研究棄渣場邊坡在降雨條件下的失穩(wěn)機理，以貴州省道真縣至甕安市高速公路路邊某棄渣場邊坡為研究對象，利用GEO-Studio中的SEEP/W模塊，在室內(nèi)試驗基礎(chǔ)上研究了棄渣場邊坡在降雨條件下暫態(tài)飽和區(qū)的分布規(guī)律，并利用SLOPE/W模塊研究了暫態(tài)飽和區(qū)分布與邊坡穩(wěn)定性之間的聯(lián)系。結(jié)果表明：① 降雨對棄渣場邊坡穩(wěn)定性的影響在降雨停止之后依舊會持續(xù)很長一段時間，這段時間安全系數(shù)會持續(xù)降低并最終達到最小值。② 坡面暫態(tài)飽和區(qū)從開始形成到坡面坡頂全部飽和只需要很短的一段時間，這段時間也是邊坡安全系數(shù)下降速度最快的時間段;坡面暫態(tài)飽和區(qū)形成以后，降雨對邊坡穩(wěn)定性的影響顯著降低。

關(guān) 鍵 詞：

棄渣場邊坡; 暫態(tài)飽和區(qū); 邊坡穩(wěn)定性; 非飽和土; 基質(zhì)吸力

中圖法分類號： U417.2

文獻標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2021.06.026

隧道開挖等工程會產(chǎn)出大量的棄渣，其主要由松散土、碎石等組成。大孔隙、非飽和是棄渣的主要特點，其具有殘積層的滲透性質(zhì)，同時由于其壓實度不夠，導(dǎo)致其強度較低，容易發(fā)生滑坡泥石流等災(zāi)害。最為大眾所熟知的棄渣場滑坡災(zāi)害就是發(fā)生在2015年12月的深圳光明新區(qū)棄渣場滑坡災(zāi)害，滑坡共造成4人死亡，76人失聯(lián)。棄渣場滑坡的危害程度之大在學(xué)術(shù)界已經(jīng)形成共識，迄今為止已有不少學(xué)者對此開展了研究。吳謙等[1]對某棄渣場進行了研究，分析了邊坡穩(wěn)定系數(shù)對參數(shù)的敏感性，并借助Monte Carlo算法對其穩(wěn)定性可靠度進行了研究。李葉鑫[2]采用主成分分析方法，分析了不同堆積時間下的棄渣場入滲特征。吳偉東等[3]建立基于（PPC）模型的風(fēng)險評價體系，結(jié)合貴昆線鐵路10座棄渣場研究發(fā)現(xiàn)棄渣場邊坡條件和地形與地基條件兩個維度對棄渣場的綜合風(fēng)險影響最大。孫朝燚等[4]為了研究坡體走向?qū)吰路€(wěn)定性的影響，針對某核電廠棄渣場，通過現(xiàn)場地質(zhì)調(diào)查和區(qū)域劃分研究了棄渣場邊坡的走向夾角對其潛在破壞機制的影響。聶超等[5]結(jié)合前期勘察資料模擬研究了不同降雨條件下地下水位等值線分布狀態(tài)，并對不同降雨條件下高填方邊坡變形特征進行研究分析，得出不同降雨條件下的滲流場分布情況以及相對應(yīng)的變形規(guī)律。高飛等[6]以某棄渣場設(shè)計為背景，通過創(chuàng)建三角網(wǎng)體積曲面探討了AutoCAD Civil 3D在棄渣場水土保持工程設(shè)計中應(yīng)用的主要優(yōu)點，以及該軟件在水土保持工程設(shè)計中的應(yīng)用前景。

本文在前人研究的基礎(chǔ)上對棄渣場邊坡降雨條件下暫態(tài)飽和區(qū)發(fā)生發(fā)展規(guī)律進行了分析研究，并將其發(fā)展規(guī)律與棄渣場邊坡安全系數(shù)的變化結(jié)合起來進行了綜合分析。研究得到了該類型棄渣場邊坡暫態(tài)飽和區(qū)發(fā)展規(guī)律以及與邊坡穩(wěn)定性之間的聯(lián)系，可為棄渣場的選址、坡面排水設(shè)施的設(shè)計施工、邊坡穩(wěn)定性監(jiān)測等工程提供參考。

1 數(shù)值分析理論

基于達西定律的二維有限元方程為[7]

ΚH+MH=Q（1）

式中：[K]為單元特征矩陣;{H}為節(jié)點水頭向量;[M]為單元質(zhì)量矩陣;{Q}為節(jié)點流量向量。

1980年Van Genuchten[8]基于大量試驗提出了VG統(tǒng)計模型，如式（2）和式（3）所示。

θ=θr+θs-θr1+ψanm（2）

式中：θ為體積含水率;θs為飽和體積含水率;θr為殘余體積含水率;ψ為基質(zhì)吸力;a，m，n為擬合參數(shù)。

kr=1-ψan-11+ψan-m21+ψanm2（3）

式中：kr為相對滲透系數(shù);由于忽略了孔隙氣壓的影響即ua=0，所以-p=u即為基質(zhì)吸力;ψ為基質(zhì)吸力;a，m，n為擬合參數(shù)。

本文采用VG模型擬合試驗數(shù)據(jù)的方法得到滲流過程中基質(zhì)吸力和體積含水率以及滲透系數(shù)的關(guān)系。同時采用極限平衡法對邊坡的穩(wěn)定性進行分析，利用靜力力學(xué)平衡和土的Mohr-Coulomb破壞準(zhǔn)則對土體穩(wěn)定性進行分析判斷，即通過分析極限狀態(tài)下坡體的穩(wěn)定性得到坡體的安全系數(shù)[9]。采用經(jīng)典的Bishop法進行穩(wěn)定性分析。

2 項目概況

本文研究點地貌類型為構(gòu)造侵蝕溶蝕低中山地貌，坡度多在30°～50°，坡體上部坡度60°～80°，懸崖峭壁遍布并伴有崩塌堆積體（見圖1）。坡面多堆積，坡腳一般有沖刷能力較強的溪流流過。巖性主要為二疊系灰?guī)r及志留系泥巖及泥灰?guī)r等，少量奧陶系的泥巖、頁巖及灰?guī)r。其中位于二疊系地層段的巖溶地貌發(fā)育，溶溝、溶槽較發(fā)育。主體山脈大致呈東西走向。整體地勢西高東低，最高海拔高程807.5 m，最低海拔高程767.0 m，相對高差40.5 m。

棄渣體地層出露主要為揚子一級地層區(qū)，路線走廊帶出露的地層有第四系、二疊系、志留系、奧陶系、寒武系等。巖性以灰?guī)r、白云巖為主夾少量泥巖、頁巖。棄渣體所處屬亞熱帶季風(fēng)濕潤氣候區(qū)，年平均氣溫16.2 ℃，年平均降水量1 071.8 mm。

綜上所述，基巖為相對穩(wěn)定的灰?guī)r組成，且坡度較陡;同時，考慮到其降雨量充沛，應(yīng)重點考慮降雨對其穩(wěn)定性的影響，故而，在現(xiàn)場勘察的基礎(chǔ)上切取圖2所示棄渣體的剖面，建立二維數(shù)值分析模型。

3 計算參數(shù)與方案設(shè)計

為了對邊坡孔隙水壓力、暫態(tài)飽和區(qū)以及棄渣體穩(wěn)定性進行分析，利用簡單高效的地質(zhì)工程通用仿真軟件GEO-Studio中的SEEP/W模塊建立如圖2所示計算模型。SEEP/W模塊是一款分析多孔透水材料中的地下水滲流以及孔隙水壓力分布狀態(tài)等問題的數(shù)值分析軟件，對飽和非飽和問題中考慮歷時的研究尤其出色。為了分析不同時刻不同深度孔隙水壓力分布情況，取x=50 m處斷面孔壓進行分析。

模型邊界條件：兩側(cè)按零流量邊界處理，底邊按不透水邊界處理，地下水水頭設(shè)置為零。而上邊界由于降雨入滲的復(fù)雜性，需要將入滲分為2種情況[10]：① 降雨強度小于邊坡入滲能力，此時降落到坡面的雨水全部滲入，按流量邊界考慮，而當(dāng)降雨強度大于邊坡入滲能力時，坡面形成積水，此時應(yīng)按壓力邊界處理。實際情況下單位流量大小隨歷時改變，為了控制降雨雨型等因素對暫態(tài)飽和區(qū)分布的影響，本文參考當(dāng)?shù)?月降雨數(shù)據(jù)，設(shè)置強度大小為0.052 m/h的均勻降雨邊界，此數(shù)值大于棄渣飽和滲透系數(shù)0.022 m/h。降雨歷時5 h，計算總時長12 d。

棄渣體土水特征曲線采用試驗數(shù)據(jù)點進行擬合得到，具體方法是將壓力板儀產(chǎn)生的吸力數(shù)據(jù)導(dǎo)入SEEP/W模塊的材料定義界面中，并基于公式（2）、（3）提到的VG模型擬合得到。同理，滲透系數(shù)曲線采用VG模型在土水特征曲線的基礎(chǔ)上擬合生成如圖3所示曲線，橫坐標(biāo)為對數(shù)坐標(biāo);圖4橫縱坐標(biāo)都為對數(shù)坐標(biāo)。

在降雨過程中為了分析邊坡穩(wěn)定性與暫態(tài)飽和區(qū)之間的關(guān)系，需在暫態(tài)飽和區(qū)各個發(fā)展階段對其穩(wěn)定性進行計算分析。通過室內(nèi)大型直剪試驗確定材料力學(xué)參數(shù)如表1所列?？紤]到擋土墻的整體性較好，且正常工作，所以在分析過程中賦予其較大的強度。

4 數(shù)據(jù)分析

4.1 棄渣場滲流分析

將降雨邊界施加到模型表面后，計算得到x=50 m斷面處體積含水率隨時間的分布，如圖5所示。降雨發(fā)生后，邊坡表面體積含水率并不是線性增大，而是以增幅逐漸降低的方式增大。在3 h后坡面體積含水率維持在0.18不再變化，此時地表棄渣土體達到飽和狀態(tài)。在5 h降雨停止后，到6 h坡面體積含水率依舊維持在0.18，表明隨著降雨的停止，坡面體積含水率并不會立馬開始下降。在y=16 m位置可以看到，由于土體性質(zhì)的不同，在層面處體積含水率有一個較大的拐點。同時由于地下水的存在，使得在基巖部分y=18 m以下深度范圍，深度越大含水率越高，且是以幅度逐漸增大的趨勢增加。

當(dāng)土體體積含水量達到飽和狀態(tài)的0.9時，即可認為此時棄渣體達到飽和狀態(tài)[11]。這對應(yīng)于本文當(dāng)棄渣體體積含水率達到0.16時，即可認為棄渣達到飽和狀態(tài)。如圖5虛線所示，在降雨停止約7 h（即計算12 h）坡面仍處于飽和狀態(tài)。由此得出結(jié)論：坡面暫態(tài)飽和區(qū)并不會隨著降雨的停止立刻開始消散。

圖6為x=50 m處孔隙水壓力隨深度的變化。在坡頂淺層部分變化規(guī)律與體積含水率相似，在降雨3 h后坡面即達到飽和狀態(tài)且孔壓不再增加。直到5 h，孔壓依舊維持不變，說明坡面積水深度不再增加，在積水達到一定深度即在坡面形成徑流沿著坡面往坡底流動。此徑流易對坡面產(chǎn)生沖刷，影響邊坡的穩(wěn)定，所以應(yīng)在坡頂設(shè)置截水溝等排水設(shè)施避免坡面被沖刷發(fā)生失穩(wěn)。

5 h降雨停止后，孔壓開始逐漸下降，這是由于地表水的下滲所致。地表孔壓在降雨停止163 h后（7 d）依舊未恢復(fù)到降雨前水平，說明降雨對邊坡的影響是一個較漫長的過程。因此在對邊坡進行安全監(jiān)測的時候，工作重心不能僅僅放在降雨期間，還要適當(dāng)延長。

同時還可以發(fā)現(xiàn)隨著降雨的進行每個時間點的拐點也是在不斷向右邊偏移。這是因為深層土體逐漸受到降雨的影響，基質(zhì)吸力開始逐漸消散，孔壓開始變大。通過該現(xiàn)象可以判斷出土體受降雨影響的深度，對應(yīng)該棄渣場，降雨影響的深度約為20 m。

在坡頂棄渣填筑較淺部分，雨水很快入滲到層面位置，由于滲透系數(shù)的差異，雨水沿著層面往坡腳滲流。如果降雨時間延長到一定閾值，雨水會在層面位置匯集，降低層面位置土體強度并導(dǎo)致棄渣體沿著層面發(fā)生滑坡失穩(wěn)。所以在對棄渣場進行選址時，在充分考慮基巖強度的同時，基巖透水性能指標(biāo)也是一個值得考慮的因素。在該棄渣場已經(jīng)填筑形成的情況下，應(yīng)在擋土墻層面位置設(shè)置排水孔，疏導(dǎo)層間水流。

圖7（a）中，降雨1 h，此時坡面未達到飽和狀態(tài)，雨水全部入滲。圖7（b），降雨2 h，坡面飽和度開始逐漸增大，坡腳飽和度相對增大得更快，并在坡頂位置形成徑流（藍色虛線）。圖7（c），降雨3 h，坡面已經(jīng)全部達到飽和狀態(tài)。圖7（d），降雨5 h，相對3 h時坡面暫態(tài)飽和區(qū)狀態(tài)，除了雨水持續(xù)入滲，坡面暫態(tài)飽和區(qū)分布不再發(fā)生變化。圖7（e），降雨停止20.5 h（即計算25.5 h后）坡頂暫態(tài)飽和區(qū)消散，緊接著26 h坡面暫態(tài)飽和區(qū)消散。

綜合以上分析筆者認為：

① 降雨發(fā)生期間，影響該棄渣場暫態(tài)飽和區(qū)分布的只是前3 h。在暫態(tài)飽和區(qū)穩(wěn)定形成以后，降雨對暫態(tài)飽和區(qū)分布的影響顯著降低，此時因為坡面已經(jīng)飽和并形成徑流。

② 降雨期間，飽和區(qū)從坡腳逐漸發(fā)展到坡頂，而降雨停止之后暫態(tài)飽和區(qū)的消散是從坡頂往坡腳逐漸消散。

4.2 棄渣場邊坡穩(wěn)定性分析

棄渣場土體力學(xué)參數(shù)如表1所列，采用Bishop法對邊坡穩(wěn)定性進行分析。在本次計算中考慮了基質(zhì)吸力對邊坡抗剪強度τ的影響，如公式（5）所示[12]：

τ=c ′+σn-σatanφ′+ua-uwtanφb（5）

式中：τ為抗剪強度;c′有效凝聚力;φ′為有效內(nèi)摩擦角;ua-uw為基質(zhì)吸力，由于不考慮氣壓影響，所以-uw即為基質(zhì)吸力;φb為基質(zhì)吸力引起的內(nèi)摩擦角，本文中φb的大小由土水特征曲線控制。

在對棄渣場穩(wěn)定性進行分析的過程中，考慮到棄渣場基巖的穩(wěn)定性很好，所以在進行棄渣場穩(wěn)定性分析的過程中將基巖設(shè)置為不可滑動層。為了簡化計算，將擋土墻強度模型也設(shè)置為摩爾庫倫模型?？紤]到隨著降雨類型以及降雨強度的不同，邊坡會表現(xiàn)出淺層局部以及深層整體兩種滑動類型，對棄渣體邊坡采用指定滑動面位置的方法分別從其局部以及整體兩個方面進行分析研究。即指定局部以及整體兩個滑面的大致范圍，由軟件在此范圍內(nèi)分別自動搜索得到兩個最危險滑面，局部滑動面以及整體滑動面能夠代表降雨對淺層滑動以及深層滑動兩種滑動形態(tài)的影響程度。圖8為滑動面位置示意圖，為了研究暫態(tài)飽和區(qū)分布與邊坡穩(wěn)定性之間的聯(lián)系，將暫態(tài)飽和區(qū)的發(fā)展劃分為：坡腳開始飽和階段、坡頂全飽和階段、坡頂消散階段、坡腳暫態(tài)飽和區(qū)消散4個階段。4個階段的時間節(jié)點分別對應(yīng)2.0，3.0，25.5，26.0 h。

降雨對邊坡穩(wěn)定性影響的方式主要有3種：① 降雨導(dǎo)致土質(zhì)軟化，凝聚力c、內(nèi)摩擦角φ降低;② 隨著降雨入滲，雨水填充土體內(nèi)部孔隙，基質(zhì)吸力降低，導(dǎo)致邊坡抗剪強度降低;③ 雨水滲入坡體，導(dǎo)致其自重增大，增大了下滑力，導(dǎo)致其穩(wěn)定性降低[13]。

本文在前人研究的基礎(chǔ)上，將暫態(tài)飽和區(qū)發(fā)展分布規(guī)律與邊坡安全系數(shù)的變化建立聯(lián)系進行分析。

圖9的局部滑動面中：在降雨初期（0～2 h），邊坡安全系數(shù)幾乎以勻速降低。而在2～5 h，安全系數(shù)降低的幅度是逐漸減小的。在25.5，26.0 h分別為坡頂與坡面暫態(tài)飽和區(qū)消散的臨界時間點，從局部安全系數(shù)變化曲線來看，25.5，26.0 h時并未表現(xiàn)出明顯的數(shù)據(jù)突變。在192 h，在雨水長時間的入滲之后，棄渣體邊坡安全系數(shù)開始恢復(fù)。

分析認為：0～2 h時，由于坡面未飽和，所以此時的降雨全部滲入土體，未在坡面匯集形成徑流流失。且此時入滲雨水主要集中在坡頂位置，不會出現(xiàn)軟化滑動帶土體的情況。所以前2 h降雨入滲對棄渣體安全系數(shù)的影響主要體現(xiàn)在加大了其自重使下滑力增大。2～5 h時坡面已經(jīng)達到飽和狀態(tài)形成暫態(tài)飽和區(qū)，此時雨水的入滲就變成了有壓入滲，且會在坡面形成徑流，雨水并非完全滲入土體，此時下滑力的增大趨勢變緩，安全系數(shù)的下降速率降低。192 h時由于雨水逐漸從坡體內(nèi)排出（或蒸發(fā)），安全系數(shù)開始增大。此現(xiàn)象與前文棄渣體邊坡穩(wěn)定性的監(jiān)測應(yīng)持續(xù)一段時間的建議論相吻合。

圖10的整體滑動面中，對比局部安全系數(shù)的變化可以發(fā)現(xiàn)，局部安全系數(shù)、整體安全系數(shù)與降雨歷時的變化規(guī)律與趨勢幾乎相同，不同點是降雨對局部安全系數(shù)的影響更大，整體安全系數(shù)的增大點出現(xiàn)的時間更晚。

筆者分析認為：由于整體穩(wěn)定性的下滑力涉及的土方更多，而降雨對其重度（下滑力）的增大相對于原下滑力較小。由于滲入棄渣體的雨水流出整體滑體需要的時間更長，導(dǎo)致整體安全系數(shù)增大點出現(xiàn)的時間相對更晚。兩者的變化趨勢相同，說明雖然整體滑面和局部滑面涉及的深度不同，但是雨水對該棄渣體整體以及局部安全系數(shù)的影響方式是相同的。

綜上分析，局部滑動面更容易受影響，所以在對該棄渣場進行治理的過程中應(yīng)該更多地關(guān)注淺層滑動。

結(jié)合前面暫態(tài)飽和區(qū)4個分布階段的分析，在坡面暫態(tài)飽和區(qū)形成過程中，邊坡安全系數(shù)下降速度最快。其中局部安全系數(shù)在坡腳開始飽和階段（2 h）下降速度最快;整體安全系數(shù)在坡腳開始飽和階段（2 h）和坡頂全飽和階段（3 h）的下降速度相同。

5 結(jié) 論

（1） 坡面暫態(tài)飽和區(qū)并不會隨著降雨的停止立刻開始消散。當(dāng)兩層土質(zhì)滲透系數(shù)相差較大時，降雨會沿著層面流動，所以在對棄渣場進行選址時，應(yīng)考慮場地中不同巖層滲透系數(shù)的差異。在該棄渣場已經(jīng)填筑完成的情況下，應(yīng)在擋土墻層面位置設(shè)置排水孔，疏導(dǎo)層間水流。

（2） 降雨對邊坡的影響是一個漫長的過程，降雨停止后的一段時間邊坡孔壓仍未恢復(fù)至正常水平。

（3） 降雨發(fā)生期間，影響該棄渣場暫態(tài)飽和區(qū)分布的只是前3 h;在暫態(tài)飽和區(qū)穩(wěn)定形成以后，降雨對暫態(tài)飽和區(qū)分布的影響顯著降低。

降雨期間，飽和區(qū)從坡腳逐漸發(fā)展到坡頂，而降雨停止之后暫態(tài)飽和區(qū)的消散是從坡頂向坡腳逐漸消散。

（4） 局部滑動面更容易受降雨的影響，所以在對該棄渣場進行治理的過程中應(yīng)該更多地關(guān)注淺層滑動。整體安全系數(shù)恢復(fù)更慢，如果需判斷棄渣體降雨后的整體穩(wěn)定性則需延長監(jiān)測時長。

（5） 前3 h是暫態(tài)飽和區(qū)形成的過程，安全系數(shù)下降速度最快，而當(dāng)暫態(tài)飽和區(qū)達到相對穩(wěn)定狀態(tài)以后，降雨對邊坡穩(wěn)定性的影響逐漸降低。

參考文獻：

[1] 吳謙，毛雪松，劉龍旗，等.某棄渣場邊坡穩(wěn)定性的可靠度分析[J].桂林理工大學(xué)學(xué)報，2017，37（3）：475-480.

[2] 李葉鑫，史東梅，呂剛.不同恢復(fù)年限棄渣場入滲特征研究與評價[J].水土保持學(xué)報，2017，31（3）：91-95.

[3] 吳偉東，茍?zhí)魄?，許博浩，等.基于改進PPC模型的鐵路棄渣場綜合風(fēng)險評價系統(tǒng)[J].中國安全生產(chǎn)科學(xué)技術(shù)，2019，15（8）：181-186.

[4] 孫朝燚，陳從新，鄭允，等.基于空間效應(yīng)的棄渣場邊坡穩(wěn)定性方法探討[J].西南交通大學(xué)學(xué)報，2019，54（1）：97-105.

[5] 聶超，韓文喜，宋柏冰.西南某機場高填方邊坡于不同降雨條件下的變形影響[J].科學(xué)技術(shù)與工程，2018，18（3）：174-183.

[6] 高飛，何京拔，湛漢溪.AutoCAD Civil 3D軟件在棄渣場設(shè)計中的應(yīng)用[J].人民長江，2017，48（12）：12-14，41.

[7] 任青文，張林飛，沈雷，等.考慮非飽和滲流過程的巖體變形規(guī)律分析[J].巖石力學(xué)與工程學(xué)報，2018，37（增2）：4100-4107.

[8] 張海闊，姜翠玲，李亮，等.基于HYDRUS-1D模擬的變水頭入滲條件下VG模型參數(shù)敏感性分析[J].河海大學(xué)學(xué)報（自然科學(xué)版），2019，47（1）：32-40.

[9] 盧玉林，陳曉冉.地下水滲流作用下土坡穩(wěn)定性的簡化Bishop法解[J].應(yīng)用力學(xué)學(xué)報，2018，35（3）：524-529，687.

[10] 毛雪松，宋玉品，吳謙，等.滲流條件下棄渣場邊坡的穩(wěn)定性分析[J].蘭州理工大學(xué)學(xué)報.2018，44（2）：139-143.

[11] 曾鈴，史振寧，付宏淵，等.降雨入滲對邊坡暫態(tài)飽和區(qū)分布特征的影響[J].中國公路學(xué)報，2017，30（1）：25-34.

[12] 藺港，孔令剛，詹良通，等.基于太沙基土拱效應(yīng)考慮基質(zhì)吸力影響的松動土壓力計算模型[J].巖土力學(xué)，2015，36（7）：2095-2104.

[13] 李寧，劉冠麟，陳有亮，等.降雨誘發(fā)淺層滑坡影響因素的解析分析[J].上海理工大學(xué)學(xué)報，2018，40（1）：65-75.

（編輯：鄭 毅）

Study on distribution and stability of transient saturated zone of slope in

abandoned slag yard under rainfall condition

ZHANG Jinmei1，ZHANG Zhen2，LI Jiachun2

（1.Economic and Technological ResearchInstitute，State Grid Yubei Power Co.，Ltd，Beijing 10000，China; 2.Highway College，Chang′an University，Xi′an 710064，China）

Abstract：

In order to clarify the instability mechanism of the slope of the slag yard under rainfall conditions，the abandoned slag yard of a highway from Daozhento Weng′an City in Guizhou Province was taken as a research object.Based on indoor test parameters，we used the SEEP/W module in GEO-Studio to study the distribution of transient saturation zone of abandoned slag yard slope under rainfall conditions and also applied the SLOPE/W module to obtain relationship between slope stability and the distribution of transient saturation zone.It is concluded that：①The effect of rainfall on the stability of the slag slope will last for a long time even after the rain stops.During this period，the safety factor will continue to decrease and finally reach the minimum value.②It takes only a short period for the slope transient saturation zone to develop to the whole slope，during which the slope safety factor declines fastest.After the formation of the transient saturation zone on the slope，the influence of rainfall on the stability of slope is significantly reduced.

Key words：

abandoned slag yard;transient saturation zone;slope stability;unsaturated soil;matrix suction