某水庫滑坡在降雨作用下的穩(wěn)定性模擬分析

王 煥 黃海峰，3

(1. 三峽大學(xué) 湖北長江三峽滑坡國家野外科學(xué)觀測研究站， 湖北 宜昌 443002； 2. 三峽大學(xué) 三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心，湖北 宜昌 443002； 3. 三峽大學(xué) 湖北省地質(zhì)災(zāi)害防治工程技術(shù)研究中心， 湖北 宜昌 443002)

?

某水庫滑坡在降雨作用下的穩(wěn)定性模擬分析

王 煥1，2黃海峰1，2，3

(1. 三峽大學(xué) 湖北長江三峽滑坡國家野外科學(xué)觀測研究站， 湖北 宜昌 443002； 2. 三峽大學(xué) 三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心，湖北 宜昌 443002； 3. 三峽大學(xué) 湖北省地質(zhì)災(zāi)害防治工程技術(shù)研究中心， 湖北 宜昌 443002)

水庫型滑坡主要受庫水位升降變化而發(fā)生變形，因此對水庫滑坡的研究主要集中在庫水變動影響下的滑坡變形機理及穩(wěn)定性分析方面．然而通過對三峽庫區(qū)滑坡調(diào)查發(fā)現(xiàn)，有一類水庫滑坡變形主要受降雨影響，而研究這類滑坡在降雨和庫水共同作用下的變形及穩(wěn)定性同樣具有重要意義．本文以三峽庫區(qū)某水庫滑坡為例，在詳細調(diào)查該滑坡地質(zhì)條件及變形發(fā)展過程基礎(chǔ)上，運用數(shù)值模擬方法對其滲流場及穩(wěn)定性進行了定量分析．結(jié)果表明，該滑坡由于特殊的地形條件，其變形主要受降雨影響，受庫水位升降影響較小．

水庫滑坡； 降雨； 滲流場； 穩(wěn)定性分析

0 引 言

三峽大壩建成以后，三峽庫區(qū)汛期限制水位為145 m，以便攔蓄洪水，非汛期水位最高則達到175 m[1]．因此，庫水位的大幅上升和下降成為誘發(fā)水庫滑坡地質(zhì)災(zāi)害的主要影響因素．如世界著名的意大利瓦依昂滑坡[2-4]，我國湖南省柘溪水庫的塘巖光滑坡[5]，以及三峽庫區(qū)的千將坪滑坡[6]等，都是在水庫蓄水期間發(fā)生的災(zāi)難性滑坡．據(jù)不完全統(tǒng)計，三峽庫區(qū)在175 m水位范圍內(nèi)共有大小滑坡2 000余處，各類變形體更是廣泛分布[7-8]．庫水周期性波動對岸坡產(chǎn)生的浸泡軟化、流水沖刷、浪蝕、浮力減重和動水壓力作用[9]，是導(dǎo)致大部分水庫滑坡變形失穩(wěn)的主要因素．然而，調(diào)查發(fā)現(xiàn)，三峽庫區(qū)存在一類變形主要受降雨影響的水庫滑坡，研究這類滑坡受降雨影響發(fā)生變形的特征及機理同樣具有重要實踐指導(dǎo)意義．

關(guān)于降雨影響滑坡變形的作用機理，前人做了很多研究，文高原等[10]研究了降雨過程中模型邊坡的變形、位移以及破壞特征，對比分析了此類邊坡在有無降雨情況下邊坡的破壞模式．Lee Y S等[11]表明填土邊坡失穩(wěn)破壞是由于降雨在坡體內(nèi)部的滲流作用極大的降低了邊坡內(nèi)部土體的有效內(nèi)摩擦角而引起的；Huang C C等[12]研究了不同降雨強度對邊坡破壞模式的影響，得到了不同降雨強度下隨時間變化邊坡土體的含水量、滑體滑動距離的變化曲線．李邵軍等[13]發(fā)現(xiàn)庫區(qū)典型土質(zhì)邊坡的變形破壞呈現(xiàn)為牽引破壞模式，坡體變形由前緣向后緣逐漸發(fā)展．林鴻州等[14]人把降雨強度和累積降雨量聯(lián)系起來，進而研究了降雨條件下雨水入滲對邊坡穩(wěn)定性的影響以及兩種因素對滑坡發(fā)生發(fā)展的影響；巨昆侖[15]通過在邊坡土體中埋設(shè)土壓力盒和孔隙水壓力盒，得出隨著降雨的持續(xù)進行，坡內(nèi)土壓力和孔隙水壓力不斷增大，至降雨結(jié)束后一段時間趨于穩(wěn)定．吳仁銑[16]發(fā)現(xiàn)當(dāng)孔隙水壓力超過某一臨界值，斜坡上的巖土體就會沿潛在滑移面發(fā)生滑動；譚玲[17]分析了研究降雨誘發(fā)滑坡的機理，認為降雨誘發(fā)滑坡主要通過土體的吸力降低，孔隙水壓力增大；坡體自重增大；抗剪強度參數(shù)劣化；裂縫發(fā)育產(chǎn)生動靜水壓力以及雨對坡面的沖刷等作用．

綜上所述，降雨對滑坡變形影響主要表現(xiàn)在以下方面：一是滑坡體淺層土體迅速達到飽和，坡面形成地表徑流，對坡面造成沖刷；二是滲透到坡體內(nèi)部導(dǎo)致滲流場的變化，從而使滑坡土體的含水率增大及土體抗剪強度的降低；三是濕潤鋒下移導(dǎo)致坡體基質(zhì)吸力降低，進而導(dǎo)致抗滑力減小，下滑力增大，直至失穩(wěn)[18-20]．本文選擇三峽庫區(qū)某代表性滑坡，通過詳細調(diào)查其地質(zhì)條件及變形發(fā)展過程，并采用數(shù)值模擬方法針對稍高于研究區(qū)實際降雨量(300mm)疊加庫水位變動的6種工況，對該滑坡的滲流場及穩(wěn)定性進行了定量分析，得出了該滑坡的變形特征及變形機理，可以為三峽庫區(qū)降雨型水庫滑坡變形的相關(guān)研究提供參考．

1 滑坡特征及變形機制分析

1.1 基本地質(zhì)特征

滑坡位于長江左岸，平面形態(tài)總體呈舌形(圖1、圖2)，剖面形態(tài)呈凹形(圖3)．滑體后緣位于高程290 m基巖陡坎下部；前緣高程145 m，沒于長江；左側(cè)邊界以陡緩相接凹狀地形為界；右側(cè)以土體與基巖接觸面為界，以形成土體貫穿裂縫為依據(jù)．

圖1 滑坡全貌照片

圖2 滑坡工程地質(zhì)平面圖

圖3 滑坡剖面圖

滑體物質(zhì)為第四系松散堆積碎塊石土，平均厚度10 m，黃褐色、紫紅色，土石比5∶5，土體結(jié)構(gòu)松散，碎塊石成分主要為粘土質(zhì)粉細砂巖，塊石粒徑多為20～50 cm，碎石粒徑多為2～15 cm，碎石多呈次棱角形，不規(guī)則排列．滑帶物質(zhì)為碎石土，土石比8∶2．滑床基巖巖性為侏羅系下統(tǒng)香溪群(J1x)中厚層粘土質(zhì)粉細砂巖，巖層產(chǎn)狀為傾向17～45°，傾角23～32°．發(fā)育一組節(jié)理裂隙，傾向160～190°，傾角60～70°．

斜坡體前緣直抵長江，地下水類型主要為松散堆積層孔隙水和基巖裂隙水．由于滑坡區(qū)地勢低洼因此易匯集地表徑流，地表水入滲形成松散堆積層孔隙水．地下水除了接受一定量的大氣降水的面狀垂直滲透補給外，還受滑坡后緣西部高陡剝蝕低山區(qū)的基巖裂隙地下水補給，由于坡度較陡，排水條件通暢，主要沿坡面向水庫排泄；同時在蓄水期間斜坡前部地下水接受庫水補給．據(jù)現(xiàn)場滲水試驗結(jié)果表明，滑坡體滲透系數(shù)為7.67×10-4～4.08×10-2cm/s，屬弱-中等透水．

1.2 變形過程及特征

該滑坡從2007年至2014年間，不斷發(fā)生地表變形，表1總結(jié)了滑坡變形過程及特征，其變形裂縫分布如圖2所示．

表1 滑坡變形過程及特征統(tǒng)計表

圖4 中部公路下張拉裂縫

圖5 左側(cè)邊界裂縫照片

圖6 前緣左側(cè)塌岸

總體來看，裂縫變形具有以下規(guī)律：

1)從空間上來看，裂縫首先出現(xiàn)在滑坡東部前緣，隨時間推移，向西側(cè)中部延伸，中部和后部公路依次下座，接著滑坡后緣出現(xiàn)明顯張拉裂縫，并向兩側(cè)延伸擴展形成邊界．

2)從時間上來看，滑坡每次發(fā)生明顯變形時，庫水位多在150 m左右，即庫水淹沒滑體范圍極少，但之前均出現(xiàn)了強降雨過程．

3)從變形程度來看，裂縫水平位移一般由初始的5～8 cm后逐漸增大到10 cm以上，而后緣邊界裂縫甚至達到1 m；裂縫垂直位移一般由初始的2～10 cm增大到20 cm左右；而最長裂縫的長度由初始的40 m不斷延伸增長到400 m．

1.3 變形機制分析

滑坡物質(zhì)為周邊山體的崩坡積物長期堆積所形成的碎石土，結(jié)構(gòu)松散．當(dāng)堆積物超過一定程度，打破原有平衡狀態(tài)，加之坡體前緣受長江河谷下切以及庫水不斷沖刷作用，形成良好臨空條件，最終形成滑坡．從滑坡變形過程及特征來看，滑坡多次變形主要表現(xiàn)為前緣塌岸與滑坡中后部的大量裂縫．而每次滑坡的明顯變形均伴隨有前期出現(xiàn)的連續(xù)強降雨或極端降雨過程．例如，2007年8月2日滑坡體東部出現(xiàn)變形，雨量記錄顯示8月2日前后出現(xiàn)連續(xù)大幅降雨；2008年8月5日，緊鄰東部變形區(qū)的西側(cè)坡體出現(xiàn)變形，資料顯示連續(xù)出現(xiàn)多天大幅降雨，而當(dāng)月降雨量甚至接近220 mm，等等．同時，變形時間對應(yīng)的庫水位卻一般保持在150 m左右，由于滑坡前緣高程為145 m，可見此時庫水淹沒下的滑坡體積很小．綜合以上分析認為，該滑坡變形主要與降雨密切相關(guān)，而與庫水位變動關(guān)系不大．

結(jié)合前人研究成果，分析降雨影響滑坡變形的主要作用機理包括：1)雨水入滲使滑坡體由不飽和向飽和狀態(tài)發(fā)展，降低土體的抗剪強度參數(shù)c及φ值，進而降低抗滑力；2)增加坡體自重，增加下滑力；3)形成動、靜水壓力，對滑坡體產(chǎn)生向外的推力．

綜上，該滑坡受高水位浸泡和降雨作用影響，變形則包括前緣局部塌岸和中后部坡體裂縫．其中，坡體裂縫變形是該滑坡的主要宏觀變形跡象，而降雨是主要影響因素．

2 滑坡滲流場及穩(wěn)定性數(shù)值模擬

2.1 非飽和滲流與孔壓應(yīng)力耦合理論[21]

基于飽和、非飽和土的滲流理論，將滑坡簡化為各向同性材料的飽和、非飽和問題，滲流控制方程為

式中，kx、ky為x、y向的滲透系數(shù)；H為總水頭；Q為施加的邊界流量；mw為比水容量．

方程顯示了同區(qū)域同時段流入、流出的總凈流量與體積含水率增量相同．邊界條件需要滿足：

求解以上3個方程，得到滲流計算的孔壓變化量后，相應(yīng)于孔壓改變的體積變化接著被計算，涉及到的耦合方程為

式中，[L]為耦合矩陣；Δd為增量位移；Δu為孔壓增量．

2.2 數(shù)值模型建立

本文基于飽和、非飽和土的滲流理論采用巖土軟件GeoStudio中的SEEP/W和SLOPE/W進行滑坡滲流及穩(wěn)定性模擬分析．具體過程為：首先采用SEEP/W進行滲流分析，以得到流場的水頭分布；然后將滲流場的水頭值輸入到邊坡穩(wěn)定性計算軟件SLOPE/W中，再在該軟件中用條分法計算滑坡的穩(wěn)定性系數(shù)[20]．

根據(jù)野外現(xiàn)場地質(zhì)調(diào)查及收集資料分析，選取滑坡體的主剖面I-I'(圖3)為計算剖面，建立相應(yīng)的有限元計算模型．其中節(jié)點數(shù)為4 703，單元數(shù)為4 697，最終的網(wǎng)格模型如圖7所示．

圖7 計算模型

Seep/W中滲透系數(shù)曲線、土-水特征曲線可分別由飽和體積含水率、土體的滲透系數(shù)實現(xiàn)．由于未測試本滑坡土樣的土-水特征曲線，因此采用樣本函數(shù)方法進行估算．滲透系數(shù)方程則是根據(jù)滲透系數(shù)和土-水特征曲線，采用Fredlund法得出[13]．最終得到的滑坡體土-水特征曲線和滲透性函數(shù)(如圖8～9所示)．

圖8 滑體土水特征曲線

圖9 滑體的滲透函數(shù)

同時，根據(jù)野外現(xiàn)場調(diào)查、類比三峽庫區(qū)其它類似滑坡，并結(jié)合地下水位和變形實測資料反演，最終確定的巖土體物理力學(xué)參數(shù)見表2．

表2 滑體物理力學(xué)參數(shù)表

該滑坡在水庫蓄水和降雨條件下的變形破壞是其在既定的三峽水庫運行期水位波動的條件下，耦合致災(zāi)的結(jié)果．為了確定滑坡的穩(wěn)定性變化規(guī)律，必須找到使得該滑坡在水庫運行期水位波動工況條件或者是在水庫運行期水位波動條件下的成災(zāi)降雨過程．

據(jù)相關(guān)資料顯示，宜昌地區(qū)全年最大月平均降雨量為222.8 mm(如圖10所示)．為了分析降雨對滑坡的影響情況以及預(yù)測滑坡在增大降雨強度情況下的穩(wěn)定性變化，本次計算工況選取稍大于最大月降雨量，即300 mm的降雨強度進行疊加計算．按照3 d平均分配，每天降雨強度值為100 mm．降雨添加時間為庫水位變動155～151.4 m區(qū)間內(nèi)．結(jié)合該滑坡的工程地質(zhì)條件、水文地質(zhì)條件，設(shè)置見表3．

圖10 1981～2010年宜昌地區(qū)月平均降雨量(資料來源：據(jù)中國氣象數(shù)據(jù)網(wǎng))

工況水位從175m降至159m水位從159m降至145m降雨情況10.13m/d0.6m/d無20.13m/d0.8m/d無30.13m/d1.0m/d無40.13m/d0.6m/d降雨強度為300mm/3d50.13m/d0.8m/d降雨強度為300mm/3d60.13m/d1.0m/d降雨強度為300mm/3d

2.3 滑坡滲流場變化特征

1)175 m降至159 m水位(無降雨)滲流場模擬

1～3工況下滲流場邊界條件設(shè)置如下：首先以175 m水位作為穩(wěn)定水頭，右邊界設(shè)置190 m的穩(wěn)定水頭邊界，下邊界與左邊界考慮為不透水邊界，通過SEEP計算出175 m水位穩(wěn)定以后的地下水位線，如圖11所示．然后把該水位線設(shè)置為初始水位線，水位以下為水頭邊界，設(shè)置為不同工況下的庫水位變動．

圖11 以不同速率下降至145 m時的地下水位線

由上圖可以看出，庫水位分別以不同速率從175 m水位降至145 m水位時刻，某滑體內(nèi)中前部的地下水位線隨著庫水的下降而下降，庫水緩慢滲出滑體，滑坡體內(nèi)地下水位線的下降滯后于庫水位的下降，滑坡體內(nèi)外形成一定的水頭差，地下水位線呈現(xiàn)一定上凸趨勢．為了能直觀的分析，將前緣滲流場進行局部放大，分析發(fā)現(xiàn)，庫水位下降速率越大，地下水線上凸效果越明顯．水位下降過程中，滲透壓力指向坡外，將對滑坡穩(wěn)定性產(chǎn)生一定不利影響．

2)175 m降至145 m水位(降雨強度300 mm/3 d)滲流場模擬

由庫水位下降與降雨同步進行滲流場模擬可知，降雨作用對地下水滲流場作用明顯，尤其是在滑坡前緣部分(如圖12所示)．具體來看，降雨前地下水位在滑坡體表面以下3.2 m左右；降雨期間，地下水位線與滑坡體表面重合，滑坡體淺層部分達到暫飽和狀態(tài)．

圖12 降雨引起的地下水位變化圖

3)庫水位降速為1.0 m/d

連續(xù)降雨3 d之后的滑坡體孔隙水壓力分布情況如圖13所示．

圖13 孔隙水壓力分布圖

4)孔隙水壓力分布圖(kPa)

對于1～6工況，選取滑坡體的A點處進行孔隙水壓力分析，時間為降雨前一天至降雨后的6 d，結(jié)果如圖14所示．

圖14 孔隙水壓力變化曲線

由圖可知，對于6種工況，滑坡體A點在降雨前一天的孔隙水壓力均為-300 kPa左右．降雨發(fā)生以后，孔隙水壓力逐漸增大，即基質(zhì)吸力逐漸減小(負孔隙水壓力即為基質(zhì)吸力)．降雨第一天基質(zhì)吸力為127 kPa，至降雨第3 d孔隙水壓力值變?yōu)?，此時土體達到暫飽和狀態(tài)．對于工況1～3，曲線幾乎水平，說明庫水位以不同速率下降對于A點處的基質(zhì)吸力影響很小．對于4～6工況，曲線變化一致，都呈現(xiàn)先迅速下降，然后緩慢回升，降雨作用對滑坡體的基質(zhì)吸力影響顯著．

2.4 穩(wěn)定性分析

利用Geo-slope軟件得出戴家坪滑坡的滑定性系數(shù)變化如圖15～16所示，其中，1-6工況下的穩(wěn)定性系數(shù)最小值依次為1.001 1，1，0.999 1，0.993，0.991，0.989 1．

圖15 175 m降至159 m的穩(wěn)定性系數(shù)曲線圖

圖16 159 m降至145 m的穩(wěn)定性系數(shù)圖

從以上結(jié)果可以看出，各工況下庫水位以0.13 m/d從175 m下降到159 m的過程中，各工況下的穩(wěn)定性系數(shù)曲線重合，變化趨勢為穩(wěn)定性系數(shù)與庫水位呈正相關(guān)，隨著庫水位下降而下降，由1.027降至1.016，穩(wěn)定性系數(shù)下降0.011．主要是由于庫水位下降改變滑坡的滲流場，是動水壓力指向坡體外側(cè)，增大下滑力，進而降低滑坡的穩(wěn)定性．

庫水位從159 m下降至145 m的過程中，穩(wěn)定性系數(shù)曲線圖特點主要為：1)工況1、2、3的最小值，依次遞減，相差最大為0.002．2)降雨前，工況1、2、3分別與4、5、6曲線重合．降雨作用以后，工況1、2、3與4、5、6的曲線開始分開，隨著時間的推移，穩(wěn)定性系數(shù)差值越來越大．由表3可以看出，穩(wěn)定性系數(shù)最小值相差0.008 1．工況2和工況5的穩(wěn)定性系數(shù)最小值相差0.009．工況3和工況6的穩(wěn)定性系數(shù)最小值相差0.01．3)6種工況對應(yīng)的穩(wěn)定性曲線都表現(xiàn)出先下降后上升，最后趨于穩(wěn)定的特點．

綜合以上分析，可以得出：1)庫水位下降速率越小，穩(wěn)定性系數(shù)降幅越小．主要是由于庫水位下降時間越長，越有利于滑坡體中的地下水排出，使動水壓力降低，基質(zhì)吸力增大，從而有利于滑坡的穩(wěn)定性．2)降雨工況與無降雨工況對應(yīng)的滑坡穩(wěn)定性系數(shù)差值范圍為0.008～0.01．而不同庫水位下降速率對應(yīng)的滑坡穩(wěn)定性系數(shù)差值最大為0.002．并且無降雨工況對應(yīng)的穩(wěn)定性系數(shù)都大于1，而降雨之后穩(wěn)定性系數(shù)均小于1.由原來的穩(wěn)定轉(zhuǎn)為不穩(wěn)定狀態(tài)．可見降雨是影響該滑坡的主要因素．3)穩(wěn)定性系數(shù)先下降后上升主要由于隨著時間推移，經(jīng)過庫水浸泡過的土體水分得以蒸發(fā)或者流出土體，使得滑坡體中的孔隙水壓力降低，基質(zhì)吸力增大，穩(wěn)定性得到提高．

3 結(jié) 論

本文以三峽庫區(qū)某降雨型水庫滑坡為例，根據(jù)滑坡變形過程及特征，結(jié)合數(shù)值模擬分析方法，對該滑坡的變形機理及穩(wěn)定性進行了系統(tǒng)研究，主要獲得以下認識：1)滑坡變形過程及特征表明，該滑坡變形主要表現(xiàn)為中后部產(chǎn)生大量裂縫以及前緣部分塌岸，主要影響因素應(yīng)該是降雨．2)滲流場模擬表明，該滑坡體主要受降雨入滲作用引起滲流場變化，從而引發(fā)滑坡主要變形；庫水位變化對滑坡的變形影響范圍有限；3)穩(wěn)定性分析表明，該滑坡在某些降雨條件下，穩(wěn)定性明顯下降；而庫水對其穩(wěn)定性影響較小．4)綜合分析表明，該滑坡雖然是水庫滑坡，但由于滑體組成為周邊山體的崩坡積物長期堆積所形成的碎石土，結(jié)構(gòu)松散，加之滑坡區(qū)地勢低洼因此易匯集地表逕流，雨水入滲使滑坡體由不飽和向飽和狀態(tài)發(fā)展，加重土體自重，降低土體的抗剪強度參數(shù)c及φ值以及基質(zhì)吸力，促使滑坡失穩(wěn)向下滑移，產(chǎn)生微裂縫，裂縫的產(chǎn)生又加劇了雨水入滲，進而降低了滑坡的穩(wěn)定性．

本文研究表明，雖然庫水變化是誘發(fā)絕大多數(shù)水庫滑坡變形的主要因素，但降雨的影響也不容忽視．由于本文在數(shù)值模擬時沒有考慮滑坡體裂縫，庫水位的連續(xù)周期性變化，以及數(shù)值模擬本身的局限性，所得的滲流場變化及穩(wěn)定性變化可能較實際情況偏小，但是足以反映其變化規(guī)律．通過連續(xù)的地質(zhì)調(diào)查和數(shù)值模擬分析，有利于揭示和掌握該類滑坡變形機理，為其他同類滑坡的穩(wěn)定性分析起到重要的參考意義．

[1] 易 武，孟召平，易慶林．三峽庫區(qū)滑坡預(yù)測理論與方法[M]．北京：科學(xué)出版社，2011．

[2] HELLER V，HAGER W H．Wave Ttypes of Landslidegenerated Impulse Waves[J]．Ocean Engineering，2011，38(4)： 630-640．

[3] 劉傳正．意大利瓦依昂水庫滑坡五十年[J]．水文地質(zhì)工程地質(zhì)，2013，40(5)：卷首語．

[4] 王蘭生．意大利瓦依昂水庫滑坡考察[J]．中國地質(zhì)災(zāi)害與防治學(xué)報，2007，18(3)：145-148．

[5] 肖詩榮，劉德富，胡志宇．世界三大典型水庫型順層巖質(zhì)滑坡工程地質(zhì)比較研究[J]．工程地質(zhì)學(xué)報，2010，18(1)：52-59．

[6] 殷躍平，彭軒明．三峽庫區(qū)千將坪滑坡失穩(wěn)探討[J]．水文地質(zhì)工程地質(zhì)，2007，34(3)：5-54．

[7] 孫廣忠．中國典型滑坡[M]．北京：科學(xué)出版社，1998．

[8] 陳 飛．長江流域地質(zhì)災(zāi)害及防治[M]．武漢：長江出版社，2007．

[9] 黃潤秋，許 強，戚國慶．降雨及水庫誘發(fā)滑坡的評價與預(yù)測[M]．北京：科學(xué)出版社，2007．

[10] 文高原，姚鵬運．降雨前、后夯實填土邊坡破壞模式試驗研究[J]．巖石力學(xué)與工程學(xué)報，2005，24(5)：747-754．

[11] Lee Y S，Cheuk C Y，SassaK． Instability Caused by AseepageIm-pediment in Layered Fill Slopes[J]． Canada Geotech，2008，45：1410-1425．

[12] Huang C C，Lo C L，Jang J S，et al． Internal Soil Moisture Response to Rainfall-induced Slope Failures and Debris Discharge[J]． Engineering Geology，2008，101：134-145．

[13] 李邵軍，Knappett J A， 馮夏庭．庫水位升降條件下邊坡失穩(wěn)離心模型試驗研究[J]．巖石力學(xué)與工程學(xué)報， 2008， 27(8)： 1586-1593．

[14] 林鴻州，于玉貞，李廣信，等．降雨特性對土質(zhì)邊坡失穩(wěn)的影響[J]．巖石力學(xué)與工程學(xué)報，2009，28(1)：198-204．

[15] 巨昆侖．長期干旱背景下降雨型松散堆積體滑坡機理研究[D]．西安：長安大學(xué)，2015．

[16] 吳仁銑．降雨誘發(fā)的滑坡作用機制研究[D]．長沙：中南大學(xué)，2013．

[17] 譚 玲．降雨誘發(fā)松散土質(zhì)滑坡機理研究[D]．重慶：交通大學(xué)，2013．

[18] 張桂榮，程 偉．降雨及庫水位聯(lián)合作用下秭歸八字門滑坡穩(wěn)定性預(yù)測[J]．巖土力學(xué)，2011，32(增 1)：476-482．

[19] 文 星，許 強，童龍云．三峽庫區(qū)黃土坡滑坡降雨入滲模型研究[J]．巖土力學(xué)，2013，34(12)：3527-3533．

[20] Geo-slope International Ltd．Geo-slope Users'Manual[M]． Calgary：Alberta，2001．

[21] 魏寶龍．降雨條件下某滑坡堆積體穩(wěn)定性研究[J]．三峽大學(xué)學(xué)報(自然科學(xué)版)，2015，37(2)：21-25．

[責(zé)任編輯 周文凱]

Simulation Analysis of Reservoir Landslide Stability under Rainfall Condition

Wang Huan1,2Huang Haifeng1,2,3

(1. National Field Observation & Research Station of Landslide in Three Gorges Reservoir Area of Yangtze River，China Three Gorges Univ., Yichang 443002, China; 2.Collaborative Innovation Center for Geo-Hazards & Eco-Environment in Three Gorges Area, China Three Gorges Univ., Yichang 443002, China; 3. Hubei Engineering Research Center of Geological Hazards Prevention, China Three Gorges Univ., Yichang 443002, China)

The deformation of reservoir landslide is mainly affected by water level changes; so the study of reservoir landslide mainly focuses on the mechanism of landslide deformation and stability analysis under the influence of reservoir water level changes. However, by investigating the landslide in Three Gorges area we found that there is a certain type of reservoir landslide whose deformation is mainly influenced by rainfall; so researching the deformation and stability of this type under the influence of both rainfall and reservoir water level change will be of great importance.Taking a reservoir landslide in Three Gorges area for example, using numerical simulation method to analyze the seepage field and stability quantitatively based on detailed investigations of the geology and deformation process of the landslide. The results indicate that the landslide due to the special topographical conditions, the deformation is mainly affected by rainfall, less affected by the reservoir water level fluctuation.

reservoir landslide; rainfall; seepage field; stability analysis

10.13393/j.cnki.issn.1672-948X.2016.05.009

2016-04-24

國家自然科學(xué)基金項目(41302260)；湖北省科技支撐計劃項目(2015BCE070)；湖北省自然科學(xué)基金創(chuàng)新群體項目(2015CFA025)；湖北省水電工程智能視覺監(jiān)測重點實驗室開放基金(2014KLA11)資助

黃海峰(1978-)，男，副教授，博士，主要從事滑坡災(zāi)害監(jiān)測、3S技術(shù)應(yīng)用等方面的研究． E-mail： hhf@ctgu.edu.cn

P642.22

A

1672-948X(2016)05-0046-07