李 俊，林 鋒，丁秀美，葛海龍

（成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點(diǎn)實(shí)驗室，四川 成都 610059）

0 引言

庫岸邊坡穩(wěn)定問題是國內(nèi)外相關(guān)學(xué)者長期關(guān)注的焦點(diǎn)。庫區(qū)兩岸斜坡的穩(wěn)定與水位的波動起伏有直接關(guān)系，關(guān)系到水上航運(yùn)交通及下游地區(qū)生命財產(chǎn)安全。Z.G.Niu等[1]研究了洪水滲透條件下堤邊坡的穩(wěn)定性并模擬了泄洪時二維非穩(wěn)定滲流場。近年來，已有大量有關(guān)三峽庫區(qū)滑坡穩(wěn)定性評價的成果發(fā)表[2-4]。劉雪梅[5]研究了三峽庫區(qū)萬州區(qū)地貌特征和滑坡演化過程，提出地貌演化受內(nèi)部營力子系統(tǒng)，外動力子系統(tǒng)及地質(zhì)體子系統(tǒng)的綜合作用。王佳佳[6]基于GIS柵格數(shù)據(jù)模型，應(yīng)用信息量理論對萬州滑坡開展了災(zāi)害易發(fā)性評價，建立了易發(fā)性評價指標(biāo)體系。

一般認(rèn)為，水庫滑坡主要是由于水庫蓄水或運(yùn)營導(dǎo)致岸坡水動力條件不利而演化造成的，包括庫水位的漲落以及暴雨入滲2個方面。從宏觀上來看，庫水位變動誘發(fā)滑坡機(jī)理主要表現(xiàn)為兩個方面：一是水庫水位變動及浸泡引起滑坡巖土體物理力學(xué)性質(zhì)的改變；另一個是水庫水位變動引起滑坡受力狀態(tài)的變化。降雨強(qiáng)度和持續(xù)時間影響雨水入滲和徑流排泄，會改變坡體內(nèi)的力學(xué)性質(zhì)?？赏ㄟ^計算機(jī)模擬解釋降雨引發(fā)滑坡問題，戚國慶等[7]數(shù)值模擬了降雨入滲過程中，邊坡非飽和區(qū)基質(zhì)吸力的變化，暫態(tài)飽和區(qū)的形成、發(fā)展變化及暫態(tài)水壓力的大小和分布。張培文[8]利用降雨條件下地表水入滲和產(chǎn)流藕合的計算程序，模擬了小降雨強(qiáng)度和大雨強(qiáng)度對地下水位變化的影響。由于斜坡的非穩(wěn)定滲流模型有飽和模型和非飽和模型，浸潤面以上為非飽和，以下為飽和，滲流計算時將整個幾何區(qū)域作為滲流區(qū)域，計算方法簡單。

五尺壩滑坡于2001年9月，在暴雨作用下，滑坡體位移變形。2008年汶川地震后，一級平臺和二級平臺出現(xiàn)裂縫，受蓄水影響，變形加劇。本文借鑒前人研究的成果，剔除了對該滑坡影響小的災(zāi)害因素，只考慮滑坡體自重、降雨和地下水三種營力作用。荷載組合工況能反映滑坡受力，評價結(jié)果符合實(shí)際情況。

1 滑坡區(qū)地理及地質(zhì)環(huán)境概況

1.1 滑坡區(qū)地理位置

五尺壩滑坡位于重慶市萬州區(qū)長江右 岸斜坡，緊鄰新鄉(xiāng)鎮(zhèn)和武陵鎮(zhèn)（圖1）。

圖1 五尺壩滑坡地理位置示意圖Fig.1 Location of Wuchiba landslide

1.2 滑坡區(qū)地質(zhì)環(huán)境條件

1.2.1 地形地貌

滑坡屬侵蝕堆積河谷低山丘陵地貌，主滑方向為305o，與斜坡坡向一致，坡度10～20o，滑坡邊界特征明顯，坡體平面呈橫長型，坡形外凸，剖面形態(tài)呈直線形，坡頂較平緩，后緣略呈弧型，該斜坡結(jié)構(gòu)類型為順向斜坡?；虑熬壣烊腴L江，剪出口位于長江江岸，現(xiàn)已沒入水下。后緣可見基巖出露，左右兩側(cè)均為山脊，中部發(fā)育一縱向沖溝，沖溝穿越整個滑體。

滑坡前緣高程120m，后緣高程255m，滑體長360m，寬1500m，坡度約20o?；麦w上平臺不發(fā)育。勘察揭露該滑坡體厚6～16m，平均厚11m，滑坡總面積54萬m2，總體積594萬m3?；w后緣出露基巖壁，高程在250～210m之間。中部高程210～180m之間坡面分布大片稻田和少量民房，坡面植被發(fā)育。前緣高程180m以下至江面。

1.2.2 地層巖性

研究區(qū)發(fā)育近水平地層，出露為侏羅系中統(tǒng)新田溝組（J2x）砂巖、泥質(zhì)砂巖互層及泥質(zhì)灰?guī)r；上覆較厚的第四系崩積、滑坡堆積成因的粉質(zhì)粘土夾碎塊石；滑體后緣出露基巖壁，巖性為紫紅色砂巖夾灰質(zhì)泥巖；中部高程210～180m之間巖性為粘土夾紫紅色砂巖碎塊石土。

1.2.3 水文地質(zhì)條件

長江水位高于滑坡前緣剪出口，屬于涉水滑坡，滑坡體內(nèi)地表水系不發(fā)育，匯集雨水通過沖溝迅速排入長江?；律辖ㄓ醒咛?，有蓄水。地下水類型以松散層孔隙潛水為主，主要由大氣降水補(bǔ)給，沿斜坡向下經(jīng)短途徑流，向長江排泄或在坡體前緣低洼處向外滲出，排泄途徑暢通。地下水動態(tài)主要受大氣降雨的影響，變化較大。

2 滑坡成因機(jī)制分析

五尺壩滑坡兩側(cè)寬大沖溝發(fā)育，后緣為基巖陡壁出露，有較好的發(fā)育滑坡的邊界條件?；虑熬壟R空涉江，江水常年侵蝕坡腳，削弱坡體前緣阻滑段強(qiáng)度，滑體中部發(fā)育縱向沖溝，地表水易匯集在沖溝中軟化坡體。自水庫運(yùn)行后，水位上升使坡體浸水面積增大，進(jìn)一步加劇了坡體強(qiáng)度的弱化，從滑坡變形的時間和變形現(xiàn)象上看，滑坡的發(fā)生與庫水的運(yùn)行在同一時期，滑坡變形跡象也集中在坡體前緣，中后部不明顯。所以，誘發(fā)滑坡的直接因素是庫水的周期性變化。

通過調(diào)查分析，推斷該滑坡屬于蠕滑-拉裂破壞模式，經(jīng)歷了如下三個變形破壞發(fā)展階段（圖2）。

（1）隨著長江的下蝕，斜坡前緣臨空面坡體沿下伏泥巖夾層產(chǎn)生剪應(yīng)力集中，向臨空方向發(fā)生蠕動變形，潛在滑移面產(chǎn)生（圖2a）。

（2）受降雨等影響，剪應(yīng)力集中帶附近泥巖軟化、抗剪強(qiáng)度下降，后緣裂縫產(chǎn)生，地表水下滲，產(chǎn)生靜水壓力，滑體飽水加載，穩(wěn)定性急劇下降（圖2b）。

（3）滑移面孔隙水壓力增加，上覆滑體沿滑動面產(chǎn)生滑動變形，失穩(wěn)、解體、運(yùn)動至坡腳停積，并在后期地質(zhì)營力、人類活動改造下，逐漸形成了現(xiàn)地貌特征（圖2c）。

圖2 五尺壩斜坡變形破壞發(fā)展階段Fig.2 Evolving stage of deformation and failure in Wuchiba slope

3 滑坡滲流及穩(wěn)定性分析

3.1 剖面及參數(shù)選取

五尺壩滑坡前期勘查布置有四條縱剖面，從滑坡左側(cè)的Ⅰ—Ⅰ′剖面至右側(cè)Ⅳ—Ⅳ′剖面（圖3）。根據(jù)前期鉆探資料和水文資料，確定滑帶的位置和地下水位情況，監(jiān)測發(fā)現(xiàn)Ⅱ—Ⅱ′剖面和Ⅲ—Ⅲ′剖面出現(xiàn)變形，中段、北東段未見明顯變形跡象，因此選取這兩個剖面進(jìn)行分析。計算參數(shù)取值見表1。

表1 五尺壩滑坡滲流、穩(wěn)定分析參數(shù)取值匯總

3.2 方案設(shè)計

計算采用Geostudio2007軟件中的SEEP/W模塊進(jìn)行滲流場模擬，計算根據(jù)勘察時期的地下水位155m作為初始水位建立穩(wěn)態(tài)工況，以它為初始狀態(tài)依次建立滲流工況。先進(jìn)行庫水位上升過程的模擬計算，然后再以上升后的穩(wěn)定水位作為初始水位進(jìn)行庫水位下降過程的模擬計算，庫水位每上升或下降到一個高度，將保持這一水位繼續(xù)計算十天內(nèi)的穩(wěn)定系數(shù)。滲流計算時，模型左右兩側(cè)采用定水頭邊界。模型剖分后有853個節(jié)點(diǎn)，1572個單元。

本文研究認(rèn)為滑體自重、庫水位升降、降雨是該滑坡失穩(wěn)的主要原因，綜合各種影響因素，計算工況如表2所示。

表2 穩(wěn)定性計算工況及荷載組合

圖3 五尺壩滑坡監(jiān)測工程平面布置Fig.3 Plane layout of monitoring project in Wuchiba landslide

考慮水庫蓄水對滑坡地下水的影響，地下水位線以上的堆積體采用天然重度進(jìn)行計算，以下的采用飽和重度，庫水位以下的堆積體采用浮重度進(jìn)行計算。水位升降速率為2m/d，當(dāng)庫水位從155m升到175m時，需要10天，在此基礎(chǔ)上再設(shè)置10天的從瞬態(tài)到穩(wěn)態(tài)的滲流期，則滲流計算共20天。從175m降至145m需要15天，加上10天穩(wěn)定滲流期則為25天；同理145m升至175m加上10天穩(wěn)定滲流期則為25天。于是前三個工況總天數(shù)共計70天。根據(jù)三峽庫區(qū)情況，選取降雨強(qiáng)度為220mm/d，降雨分兩次，即水位下降至165m和再次從145m上升時開始降雨，兩次降雨歷時均為3天。

3.3 穩(wěn)定性分析

數(shù)值模擬庫水位工況組合旨在獲得滑坡穩(wěn)定性系數(shù)在外界因素影響下的變化趨勢。在只考慮自重的天然斜坡，坡體內(nèi)的孔隙水壓力為負(fù)，負(fù)的孔隙水壓力可理解成基質(zhì)吸力，它能產(chǎn)生增大斜坡穩(wěn)定性的表觀粘聚力[9]，引起坡體內(nèi)塊石土抗剪強(qiáng)度的增加。對于完全處于飽和狀態(tài)的斜坡土體，基質(zhì)吸力的消減會極大地減小斜坡的穩(wěn)定性，從而可能引發(fā)坡體失穩(wěn)破壞。由圖4和圖5可以看出，庫水位經(jīng)歷了兩次上升和一次下降過程，且穩(wěn)定系數(shù)的增減與庫水位的升降呈現(xiàn)一致性變化規(guī)律。降雨明顯降低了穩(wěn)定系數(shù)，其中，Ⅱ-Ⅱ′剖面受降雨影響穩(wěn)定系數(shù)平均降低了0.05，Ⅲ-Ⅲ′剖面平均降低了0.06。總的來說，Ⅲ-Ⅲ′剖面受降雨影響程度要大于Ⅱ-Ⅱ′剖面?；路€(wěn)定特征具體分析如下：

圖4 剖面Ⅱ-Ⅱ′穩(wěn)定系數(shù)隨時間的變化Fig.4 Change of stability coefficient versus time for profile Ⅱ-Ⅱ′

圖5 剖面Ⅲ-Ⅲ′穩(wěn)定系數(shù)隨時間的變化Fig.5 Change of stability coefficient versus time for profile Ⅲ-Ⅲ′

（1）后緣基巖水位穩(wěn)定和無降雨的條件下，蓄水從145m至175m過程中，其位移、土壓力和孔隙水壓力的變化主要受庫水位控制?；略诠r一至工況三情況下，Ⅱ-Ⅱ′剖面的穩(wěn)定系數(shù)介于1.143～1.287之間，Ⅲ-Ⅲ′剖面的穩(wěn)定系數(shù)介于1.097～1.198，兩個剖面均處于穩(wěn)定狀態(tài)，只是Ⅲ-Ⅲ′剖面的穩(wěn)定性略差一些。庫水位的上升，滑坡前緣不斷被庫水浸沒，水面下的前緣部分相當(dāng)于受到了指向坡體的靜水壓力；此外，地下水從庫岸滲入斜坡，地下水位的上升滯后于庫水位的上升，所以穩(wěn)定系數(shù)增大。庫水位急劇降落，斜坡內(nèi)容易產(chǎn)生超孔隙水壓力，再加上斜坡趾部的坍塌，易引起斜坡大面積移動，降低穩(wěn)定系數(shù)。

（2）工況四從第5天開始連續(xù)降雨3天后，穩(wěn)定系數(shù)呈下降趨勢，兩個剖面降雨后的穩(wěn)定系數(shù)變化范圍分別為1.074～1.228和1.051～1.157，最后都趨于定值，且都沒有低于安全系數(shù)1.05。由于滑坡后緣臨近石質(zhì)陡坡，降雨在后緣聚水，抬升基巖裂隙水水位，降雨與庫水的耦合效應(yīng)表現(xiàn)在滑體內(nèi)地下水水位的變動，即飽水面積比變化。若浸潤線愈陡則動水壓力愈大，基巖水位抬升后與庫水連同，將對具有相對隔水地板滑坡產(chǎn)生承壓作用。這一系列的變動和起伏都會給滑坡的穩(wěn)定性造成不良影響。

（3）工況五第1天連續(xù)降雨3天后，穩(wěn)定系數(shù)也下降了，但當(dāng)水位回升到155m時，穩(wěn)定系數(shù)又開始增加，兩個剖面最終都分別趨于一個定值，且這個最終值以Ⅱ-Ⅱ′剖面居較大。現(xiàn)場勘查結(jié)果顯示，滑體自重+暴雨+庫水位175m降到145m時，滑坡變形最大，位移具有明顯的分區(qū)特征，主要分布滑坡后緣和一級平臺、二級平臺之間，出現(xiàn)在滑坡前緣坡肩處。短期內(nèi)的集中降雨，可以在較短時間內(nèi)，使坡體內(nèi)的地下水位迅速抬升，孔隙水壓力迅速增加，坡體抗剪強(qiáng)度降低，從而使滑坡的穩(wěn)定性條件迅速惡化。另一方面，降雨滲入使土體產(chǎn)生物理化學(xué)效應(yīng)，軟化滑面，不利于滑坡的穩(wěn)定。

3.4 最不利工況組合分析

考慮到五尺壩滑坡地質(zhì)條件，結(jié)合以上數(shù)值模擬分析，認(rèn)為觸發(fā)五尺壩滑坡的降雨強(qiáng)度在220mm/d以上。該區(qū)歷史統(tǒng)計的一次降雨過程累積降雨量從未超過660mm，因此單一降雨因素觸發(fā)五尺壩滑坡失穩(wěn)的可能性較小。為了評價極端條件下滑坡的穩(wěn)定情況，計算中設(shè)定降雨強(qiáng)度達(dá)到特大暴雨級別250mm/d，水位升降降速率為3m/d，降雨歷時62天，其它方案設(shè)計同工況四和工況五。

圖6 最不利工況下穩(wěn)定系數(shù)隨時間的變化Fig.6 Change of stability coefficient versus time under the worst condition

圖6顯示兩個剖面的穩(wěn)定系數(shù)與庫水位變化趨勢相同。可以看出，當(dāng)庫水位降落速率突然增大至3m/d時，穩(wěn)定系數(shù)在很短時間內(nèi)急劇減小。Ⅱ-Ⅱ′剖面當(dāng)庫水位降落至145m時，穩(wěn)定系數(shù)減至1.054，接近臨界安全系數(shù)1.05；而Ⅲ-Ⅲ′剖面的穩(wěn)定系數(shù)在水位還沒降落至最低時，已經(jīng)低于了臨界值，且水位回升到175m時才到達(dá)警戒值以上。庫水位升降速率的加快對滑體內(nèi)超孔隙水壓力有很大影響?；w內(nèi)超孔隙水壓力消散的時間決定了穩(wěn)定系數(shù)的大小，時間越長，則最低水位時的穩(wěn)定系數(shù)與計算初期的穩(wěn)定系數(shù)之差值越大。

所以可以認(rèn)為在極端條件下，五尺壩斜坡存在失穩(wěn)的危險。應(yīng)加強(qiáng)庫水位變化動態(tài)監(jiān)測和滑坡位移監(jiān)測，并做好地面的排水措施。

4 結(jié)論

（1）五尺壩滑坡地貌類型為侵蝕堆積河谷低山丘陵地貌，所處地層局部上覆較厚的第四系崩積、滑坡堆積成因的粉質(zhì)粘土夾碎塊石，由于庫水的浸泡和運(yùn)動，使滑坡體前緣處于周期性的干濕變換環(huán)境，土體中的土顆?；蚰z結(jié)物質(zhì)被帶走，使土體的力學(xué)參數(shù)發(fā)生變化。

（2）五尺壩滑坡為蠕滑-拉裂破壞模式，后緣基巖陡壁，前緣臨空涉江，在水的營力和后期人類改造下形成了現(xiàn)今地貌。

（3）瞬態(tài)滲流分析了兩個剖面在水位升降和降雨條件下滑坡的穩(wěn)定性，這兩種影響因素都會降低此滑坡的穩(wěn)定性，兩種因素共同作用更會加劇穩(wěn)定系數(shù)的喪失，但對整體穩(wěn)定性影響不大。降雨與庫水位升降導(dǎo)致的滑體內(nèi)孔隙水壓力的增減決定了穩(wěn)定系數(shù)的變化趨勢。

（4）為了應(yīng)對極端條件滑坡災(zāi)害，評價了人為設(shè)定的工況組合引起的滑坡失穩(wěn)條件，并提出了預(yù)防措施。

[1]NIU Zhi-guo，LU Jun，YOU Ri，et al.The stability analysis of the dike slope in Bdg reservoir underthe seepage of flood[J].Procedia Engineering，2011，15：5264-5267.

[2]尚柟吉，周勇，張元才，等.三峽庫區(qū)云陽青術(shù)滑坡成因機(jī)制及穩(wěn)定性評價研究[J].甘肅水利水電技術(shù)，2010，46（8）：12-16.

[3]秦凱旭，石豫川，劉漢超，等.三峽庫區(qū)某滑坡體成因機(jī)制分析與穩(wěn)定性評價[J].水土保持研究，2006，13（5）：84-86.

[4]盧書強(qiáng)，易慶林，易武，等.三峽庫區(qū)臥沙溪滑坡變形失穩(wěn)機(jī)制分析[J].中國地質(zhì)災(zāi)害與防治學(xué)報，2013，24（2）：21-24.

[5]劉雪梅.三峽庫區(qū)萬州區(qū)地貌特征及滑坡演化過程研究博士學(xué)位論文][D].武漢：中國地質(zhì)大學(xué)，2010.

[6]王佳佳，殷坤龍，肖莉麗.基于GIS和信息量的滑坡災(zāi)害易發(fā)性評價——以三峽庫區(qū)萬州區(qū)為例[J].巖石力學(xué)與工程學(xué)報，2014，3（4）：798-808.

[7]戚國慶，黃潤秋，速寶玉，等.巖質(zhì)邊坡降雨入滲過程的數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報，2003，22（4）：625-629.

[8]張培文.降雨條件下飽和—非飽和土徑流滲流耦合數(shù)值模擬研究[博士學(xué)位論文][D].大連：大連理工大學(xué)，2002.

[9]Vandamme J, Zou Q.Investigation of slope instability induced by seepage and erosion by a particle method[J].Computers and Geotechnics，2013，48：9-20.