王曉磊，王旭春，管曉明，張 鵬

(1.青島理工大學土木工程學院，山東 青島 266033;2.北京交通大學土木建筑工程學院，北京 100044;3.中國礦業(yè)大學力學與建筑工程學院，北京 100083)

降雨誘發(fā)邊坡失穩(wěn)占據(jù)了滑坡災害相當?shù)谋戎?，造成了大量的人員傷亡與財產(chǎn)損失?！吨袊卮蟮刭|(zhì)災害實例分析》［1］列舉了27個國家重大地質(zhì)災害，其中有15個就是由于降雨所引發(fā)?！吨袊湫突隆罚?］一書中實錄了90 多例滑坡，90%以上發(fā)生在雨中和雨后。降雨誘發(fā)邊坡失穩(wěn)在我國廣泛存在，造成了大量的人員傷亡與財產(chǎn)損失［3-4］。

以往的研究中，膨脹土、雜填土和裂隙黏土等含裂隙土質(zhì)邊坡的穩(wěn)定分析大多采用傳統(tǒng)的靜力平衡法和統(tǒng)計比較判別法等，坡體統(tǒng)一按飽和狀態(tài)進行計算，并且忽略了裂隙的作用［5-6］，這與現(xiàn)場實際情況差別較大。近年來，人們逐漸認識到土體的非飽和狀態(tài)以及邊坡裂隙對邊坡穩(wěn)定性的影響，并做了大量的研究工作，結(jié)果表明考慮裂隙和邊坡降雨入滲影響的穩(wěn)定性分析方法較為合理和實用［7-10］。殷宗澤等［11］提出了一種以條分法為基礎，近似反映裂隙影響的膨脹土邊坡穩(wěn)定性分析方法，將膨脹土坡劃分成3個不同強度指標的亞層。闕云等［12］對受降雨入滲影響的裂隙性黏土邊坡進行了研究，分析了降雨條件下不同裂隙深度和開口寬度等因素對裂隙性黏土邊坡滲流場的影響。以上研究對含裂隙土質(zhì)邊坡穩(wěn)定性分析理論的發(fā)展起到了極大的促進作用，為后續(xù)研究工作奠定了基礎。目前研究中多采用在整個邊坡均勻布置裂隙的方式，未考慮邊坡裂隙的出現(xiàn)形式，并且對降雨條件下坡體存在裂隙時的應變情況研究較少。

本文采用邊坡穩(wěn)定性分析及治理軟件Slide 對暴雨情形下含裂隙土質(zhì)邊坡的穩(wěn)定性進行分析，確定邊坡潛在滑動面的位置。根據(jù)分析結(jié)果布置促使邊坡發(fā)生推移式滑坡和牽引式滑坡的裂隙，最后應用MIDAS/GTS 軟件模擬分析降雨條件下的坡體應變規(guī)律和邊坡穩(wěn)定情況。計算過程中考慮降雨邊界條件、體積含水率、滲透系數(shù)與孔隙壓力之間的關系。

1 暴雨情形下裂隙對非飽和土質(zhì)邊坡穩(wěn)定性的影響

降雨尤其是暴雨是滑坡活動最重要的觸發(fā)因素，而邊坡裂隙的存在為降雨入滲提供了進入內(nèi)部的通道，改變了降雨的入滲路徑，使得降雨這一危險因素對災害的觸發(fā)能力得到充分的發(fā)揮，可以歸納為以下幾點:

a.裂隙增加了入滲能力。邊坡裂隙增大了雨水與坡體的接觸面積，在降雨過程中，地表徑流的雨水大量涌入裂隙之中，增加了雨水入滲量。隨著裂隙的發(fā)展，接觸范圍越來越大，并且在裂隙水頭的作用下，雨水的入滲能力大幅度增強，裂隙周圍土體的含水量迅速增大，邊坡土體基質(zhì)吸力下降，自重增加。基質(zhì)吸力的下降改變了非飽和土體的應力狀態(tài)，使土體產(chǎn)生應變，其外在表現(xiàn)為坡體的宏觀位移。而土體在飽和之后，自重的增加導致下滑力變大，同時巖土體黏聚力、內(nèi)摩擦角及基質(zhì)吸力減小，導致抗剪強度降低。

b.靜水壓力作用。強降雨下裂隙中充滿積水，對邊坡產(chǎn)生靜水壓力，壓力P可分解為平行和垂直于滑動面的分力，在此壓力作用下邊坡將產(chǎn)生向臨空面的位移，促使裂隙的發(fā)展。同時平行于滑動面的分力增大了邊坡的下滑力，而平行于滑移線的分力使坡體正壓力減小，降低邊坡局部穩(wěn)定性。裂隙開裂越深，裂隙積水的靜水壓力越大，靜水壓力作用越強。

c.裂隙積水滲流作用。裂隙積水使裂隙處存在水頭差，水體在裂隙周圍土體孔隙中流動，產(chǎn)生動水壓力，沿流線方向的單位體積滲透力為

式中:ρw為水的密度，kg/m3;g為重力加速度，m/s2;i為滲流坡降;L為滲流長度，m;h為水頭損失，m。動水壓力對土體骨架產(chǎn)生拖曳力，導致土體中的應力與變形發(fā)生變化。一旦邊坡滲流通道形成，土體在此滲透力的作用下，邊坡松散土體中的土顆粒被沖蝕，加劇邊坡的不穩(wěn)定性。

2 實例分析

2.1 計算模型、初始條件和邊界條件

2010年雨季期間，由于幾次強降雨，山西省安太堡露天煤礦東北角發(fā)生滑坡。本文以該煤礦周邊陡幫邊坡一勘探剖面為研究對象，具體地層尺寸及分布狀況見圖1。

圖1 安太堡露天煤礦剖面(單位:m)

坡體在計算范圍內(nèi)沒有地下水位，所以在計算過程中，初始狀態(tài)不考慮邊坡兩側(cè)土體滲流情況。Mein等［13］提出:當降雨量小于土體的入滲能力時，雨水全部入滲，當降雨量超過土體的入滲能力時產(chǎn)生地表徑流。降雨過程中，由于考慮的是暴雨情況，降雨強度超過土體的入滲能力，邊坡裂隙很快被地表徑流充滿，因此在降雨過程中，始終認為裂隙處的壓力水頭為裂隙內(nèi)水的深度，坡面不考慮積水情況。

本文的降雨入滲模型以土體體積含水率作為控制變量，通過在坡面施加節(jié)點流量的方式實現(xiàn)降雨，運用土壤水動力學的基本理論，求解給定降雨條件下的邊坡土體瞬態(tài)體積含水率的數(shù)值解，并將降雨計算結(jié)果應用于下一步邊坡穩(wěn)定性的計算之中。邊坡土體非穩(wěn)態(tài)滲流控制方程為

模型邊界條件:底部完全約束，兩側(cè)邊界水平位移約束，頂部為自由邊界。

2.2 巖土計算參數(shù)

計算中所用巖土計算參數(shù)來自現(xiàn)場的勘測資料，同一地層取試驗平均值，具體見表1。研究區(qū)域土體為非飽和土，其體積含水率的變化取決于應力狀態(tài)的變化和土體特性的變化。體積含水率是孔隙水壓力的函數(shù)，滲透系數(shù)又是體積含水率的函數(shù)，所以滲透系數(shù)也是孔隙水壓的函數(shù)。滲透系數(shù)和壓力水頭之間的關系用滲透系數(shù)曲線表示，體積含水率和孔隙水壓力之間的土水特征曲線采用Van Genuchten模型。降雨過程中體積含水率的變化對土體的強度影響采用畢肖普提出的非飽和土的強度公式［9］:

式中:σ'為有效應力;σ為總應力;ua為孔隙氣壓力;uw為孔隙水壓力;χ為反映土中氣相影響的系數(shù)。

表1 巖土計算參數(shù)

由于滑坡下部為巖質(zhì)邊坡，本文的研究對象主要為發(fā)生滑坡的上層土質(zhì)邊坡，限于篇幅只給出粉土的滲透系數(shù)曲線和土水特征曲線，見圖2和圖3。

圖2 粉土滲透系數(shù)和壓力水頭的關系

圖3 粉土體積含水率和孔隙水壓力的關系

2.3 數(shù)值模擬方案

針對安太堡礦區(qū)露天邊坡幾次不同程度滑坡均出現(xiàn)在暴雨之后的情況，本次模擬暴雨量為120 mm/d，分析不同降雨持續(xù)時間對邊坡應變和穩(wěn)定性的影響規(guī)律。為了使邊坡裂隙位置更符合實際狀態(tài)，首先采用Slide 軟件對邊坡進行分析，確定邊坡滑移線位置。然后根據(jù)模擬結(jié)論將裂隙布置如下:①潛在滑體的上部滑移線位置布置裂隙，模擬邊坡上部滑移線附近土體進水軟化，發(fā)生推移式滑坡;②邊坡中下部布置裂隙，模擬雨水進入下部土體，下部滑移線區(qū)域土體軟化，發(fā)生牽引式滑坡，詳細模擬方案見表2。

表2 數(shù)值模擬方案

2.4 計算結(jié)果分析

2.4.1 降雨過程中坡體應變規(guī)律

研究區(qū)域回填土邊坡結(jié)構(gòu)松散，孔隙性和透水性較強，使其更易受到降雨的影響，表層土體可以產(chǎn)生較大的塑性變形和顯著的流變變形。以下針對1號邊坡頂部平臺布置裂隙與頂部平臺和坡面中部同時布置裂隙兩種情況，對暴雨影響下的坡體應變進行分析:

a.1號邊坡頂部平臺布置裂隙坡體應變規(guī)律。以方案6為例，僅在1號邊坡頂部平臺出現(xiàn)張拉裂隙時，坡體內(nèi)部剪應變見圖4。隨著雨水入滲，降雨持續(xù)1 d 時，1號邊坡整個回填土區(qū)域出現(xiàn)較大的剪應變，最大值出現(xiàn)在1號邊坡坡腳位置。由于土體自重作用，坡腳位置區(qū)域存在應力集中現(xiàn)象。隨著雨水的不斷入滲，降雨持續(xù)2 d 時，1號邊坡回填土層和下層粉土接觸面附近區(qū)域也出現(xiàn)了較大的剪應變，且塑性區(qū)不斷向下部平臺擴展。到降雨持續(xù)3 d時，回填土與粉土接觸區(qū)域的剪應變范圍變大，這是由于雨水入滲使非飽和土含水量增加，基質(zhì)吸力下降，而基質(zhì)吸力的降低又使得邊坡非飽和物料產(chǎn)生應變，同時滑移線的輪廓變得非常明顯。

圖4 方案6 坡體剪應變

b.1號邊坡頂部平臺和坡面中部布置裂隙坡體應變規(guī)律。以方案8為例，1號邊坡頂部平臺和坡面中部同時布置裂隙時，坡體剪應變見圖5。降雨持續(xù)1 d 時，1號邊坡整個回填土層發(fā)生較大的剪應變，并且1號邊坡坡腳以及回填土與粉土接觸區(qū)域都出現(xiàn)剪應變極值。下方2號邊坡回填土與粉土接觸區(qū)域也形成了連貫的剪應變帶。隨著降雨的持續(xù)，剪應變區(qū)域縮小，極值只出現(xiàn)在1號邊坡坡腳下方回填土與粉土接觸區(qū)域。降雨持續(xù)3 d 時，隨著滲流的進行水分在坡體內(nèi)部運移，剪應變極值區(qū)幾乎消失。

圖5 方案8 坡體剪應變

對比以上兩種裂隙布置情況，后者降雨入滲影響的速度更快，且影響范圍較前者大。降雨持續(xù)1 d時，后者的塑性區(qū)已經(jīng)擴展到2號邊坡上方區(qū)域。而前者始終對1號邊坡存在較大的影響。降雨持續(xù)3 d后，后者的塑性區(qū)從2號邊坡一直貫通到1號邊坡的坡頂。

2.4.2 降雨過程中穩(wěn)定系數(shù)的變化規(guī)律

通過降雨滲流計算，將得到的瞬態(tài)孔隙水壓力作為荷載應用于邊坡的穩(wěn)定性分析中。采用有限元法的強度折減法，計算邊坡穩(wěn)定性。

a.1號邊坡頂部平臺布置裂隙穩(wěn)定性計算。在1號邊坡頂部布置深度為0.5 m的裂隙，隨著降雨的入滲，不同開口寬度的裂隙所對應的穩(wěn)定系數(shù)如圖6 所示。可以看出，裂隙的寬度對于安全系數(shù)的影響不大，這是因為開口寬度和裂隙充滿積水的速度有關，本文設定裂隙瞬間被雨水充滿，因而裂隙的開口寬度對邊坡穩(wěn)定性未有較強的影響;裂隙深度不同時，降雨影響下邊坡穩(wěn)定性變化見圖7。隨著裂隙深度的開展，穩(wěn)定系數(shù)逐步降低，在1～1.5 d之間穩(wěn)定系數(shù)降低速率最大，且在1.5 d 左右邊坡最不穩(wěn)定，隨后穩(wěn)定系數(shù)又出現(xiàn)反彈的趨勢。這是由于當雨水入滲到一定程度時土體達到飽和則入滲速度降到很低，此時隨著坡體內(nèi)部水體的運移，邊坡穩(wěn)定性出現(xiàn)回升。

圖6 邊坡頂部布置0.5 m 裂隙時降雨穩(wěn)定系數(shù)

圖7 邊坡頂部不同裂隙深度時降雨穩(wěn)定系數(shù)

b.1號邊坡頂部平臺和坡面中部同時布置裂隙時穩(wěn)定性計算。考慮1號邊坡頂部平臺和坡面中部同時布置裂隙，當坡面裂隙開裂較淺時，降雨產(chǎn)生的裂隙積水未能影響到重力荷載作用下坡體最危險滑移線的位置，此時隨著雨水的入滲，其穩(wěn)定系數(shù)的降低程度相比于只在坡頂布置裂隙時變化不大，如圖8 所示;但當裂隙開裂達到能夠影響到潛在滑動面位置深度時，邊坡潛在滑動面區(qū)域土體遇水軟化，強度降低，穩(wěn)定系數(shù)急劇下降。由圖5(c)可以看出坡體塑性變形向上延伸，坡體產(chǎn)生牽引式滑坡。

圖8 邊坡頂部和坡面中部布置裂隙時降雨穩(wěn)定系數(shù)

3 結(jié) 論

a.降雨過程中邊坡塑性區(qū)首先在坡腳位置出現(xiàn)，隨著降雨的不斷進行，塑性區(qū)自坡腳位置向上部不斷發(fā)展，直到形成貫通的塑性區(qū)。

b.降雨過程中穩(wěn)定系數(shù)并不是持續(xù)降低。當裂隙周圍土體完全飽和之后，土體的滲透系數(shù)變得很小，此時降雨入滲變得極為緩慢，坡體內(nèi)部的水體發(fā)生了運移，邊坡的穩(wěn)定性又得到了回升。

c.模擬結(jié)果表明在持續(xù)降雨1.5 d 左右邊坡的穩(wěn)定系數(shù)最小。當裂隙開裂較深時，邊坡穩(wěn)定系數(shù)降低到1 以下，存在發(fā)生滑坡危險。但分析結(jié)果表明，滑坡的發(fā)生存在于邊坡的表層土體。

d.本文提出的邊坡裂隙布置方法體現(xiàn)了邊坡實際情況下裂隙的產(chǎn)生以及發(fā)展，相比以往坡體均布裂隙的布置方法更接近邊坡裂隙的真實情況。

［1］鐘立勛.中國重大地質(zhì)災害實例分析［J］.中國地質(zhì)災害與防治學報，1999，10(3):1-7.

［2］殷躍平.中國典型滑坡［M］.北京:中國大地出版社，2007.

［3］張先發(fā)，李明華，張小剛.長江上游暴雨與滑坡崩塌關系［J］.地理，1995(3):102-106.

［4］楊文東.降雨型滑坡特征及其穩(wěn)定穩(wěn)定分析研究［D］.武漢:武漢理工大學，2006.

［5］廖世文.膨脹土與鐵路工程［M］.北京:中國鐵道出版社，1984.

［6］劉特洪.工程建設中的膨脹土問題［M］.北京:中國建筑工業(yè)出版社，1997.

［7］姚海林，鄭少河，李文斌，等.降雨入滲對非飽和膨脹土邊坡穩(wěn)定性影響的參數(shù)研究［J］.巖石力學與工程學報，2002，21(7):1034-1039.

［8］姚海林，鄭少河，陳守義，等.考慮裂隙及雨水入滲影響的膨脹土邊坡穩(wěn)定性分析［J］.巖土工程學報，2001，23(5):606-609.

［9］陳鐵林，鄧剛，陳生水，等.裂隙對非飽和土邊坡穩(wěn)定性的影響［J］.巖土工程學報，2006，28(2):210-215.

［10］李煥強，孫紅月，孫新民，等.降雨入滲對邊坡性狀影響的模型實驗研究［J］.巖土工程學報，2009，31(4):589-594.

［11］殷宗澤，徐彬.反映裂隙影響的膨脹土邊坡穩(wěn)定性分析［J］.巖土工程學報，2011，33(3):454-459.

［12］闕云，胡昌斌，姚曉琴.降雨入滲對裂隙性黏土邊坡穩(wěn)定性作用機理的分析［J］.福州大學學報，2009，37(3):423-429.

［13］MEIN R G，LARSON C L.Modelling infiltration during a steady rain［J］.Water Resource Research，1973，9(2):384-393.