江俊杰，劉東澤，盧應(yīng)發(fā)

（1. 湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院，湖北 武漢 430068； 2. 中冶武勘工程技術(shù)有限公司，湖北 武漢 430080）

0 引 言

自2003年三峽水庫(kù)開(kāi)始周期性蓄水以來(lái)，由于庫(kù)區(qū)的地質(zhì)條件被極大地改變，三峽庫(kù)區(qū)沿岸產(chǎn)生大量涉水滑坡［1］，這些滑坡一旦發(fā)生大規(guī)模失穩(wěn)，將對(duì)人民的生命財(cái)產(chǎn)安全造成巨大威脅［2］。

由于滑坡變形是一個(gè)復(fù)雜、作用機(jī)理不明確的過(guò)程，因此諸多學(xué)者對(duì)滑坡的變形情況與穩(wěn)定性影響因素進(jìn)行了研究，如唐軍峰等［3］利用SEEP/W 軟件對(duì)庫(kù)水位升降和降雨進(jìn)行模擬，并計(jì)算不同條件下的穩(wěn)定性系數(shù)，以此來(lái)研究不同庫(kù)水位波動(dòng)速率及降雨強(qiáng)度作用下滑坡體的滲流特征和穩(wěn)定性變化規(guī)律；曾剛等［4］采用Geo-studio 軟件得到滑坡的安全系數(shù)變化規(guī)律。分析了在庫(kù)水位變動(dòng)和最危險(xiǎn)庫(kù)水位變動(dòng)疊加降雨兩種工況下滑坡的穩(wěn)定性。謝林沖等［5］在滑坡的地表位移監(jiān)測(cè)數(shù)據(jù)、庫(kù)水位變動(dòng)及降雨資料的基礎(chǔ)上分析了滑坡的時(shí)空演化過(guò)程，并通過(guò)耦合分析確定了庫(kù)水位下降和降雨是滑坡變形的關(guān)鍵因素。王樂(lè)等［6］根據(jù)非飽和滲流原理利用有限元軟件對(duì)不同降雨類型聯(lián)合庫(kù)水位驟降情況下的滑坡穩(wěn)定性進(jìn)行分析。雷德鑫等［7］采用蒙特卡羅隨機(jī)抽樣法，分析了滑坡在不同環(huán)境下的穩(wěn)定性，同時(shí)分析了滑坡的穩(wěn)定性系數(shù)對(duì)滑體物理力學(xué)參數(shù)的敏感性。胡致遠(yuǎn)等［8］通過(guò)建立滑坡多滑帶地質(zhì)模型，來(lái)研究庫(kù)水位聯(lián)合降雨作用下滑坡的變形和穩(wěn)定性。張夏冉等［9］采用極限平衡的定量計(jì)算與破壞概率的定性分析，分析了不同庫(kù)水降速及不同降雨條件下的滑坡穩(wěn)定性。肖先煊等［10］通過(guò)建立滑坡的地質(zhì)力學(xué)模型，來(lái)分析滑坡體在降雨及庫(kù)水位升降作用下滑體內(nèi)水壓力、土壓力的變化規(guī)律及穩(wěn)定性。

諸多研究表明庫(kù)水位波動(dòng)和降雨是三峽庫(kù)區(qū)滑坡發(fā)生變形的主要誘因［11］。為分析庫(kù)水位降落聯(lián)合降雨耦合作用下的滑坡破壞機(jī)理，本研究以雞腦殼包滑坡為例，通過(guò)建立二維、三維庫(kù)水位聯(lián)合降雨的耦合數(shù)值分析模型，分析滑坡在不同環(huán)境下的變形演變規(guī)律，同時(shí)建立滑坡變形破壞的數(shù)值預(yù)測(cè)模型，預(yù)測(cè)最危險(xiǎn)情況下雞腦殼包滑坡的變形特征，再采用Morgenstern-Prince 法計(jì)算不同庫(kù)水位降速率、不同降雨強(qiáng)度及降雨歷時(shí)組成的降雨過(guò)程對(duì)應(yīng)的穩(wěn)定系數(shù)，為三峽庫(kù)區(qū)水位變化及降雨條件下沿岸滑坡的監(jiān)測(cè)防治提供參考。

1 數(shù)值分析模型

1.1 工程概況

雞腦殼包滑坡位于巫山縣曲尺鄉(xiāng)柑園村，長(zhǎng)江干流長(zhǎng)江左岸，滑坡前緣涉水，距三峽大壩141.30 km，滑坡體處于巫山縣向斜南冀，附近無(wú)大的斷層，裂隙不發(fā)育。

雞腦殼包滑坡體長(zhǎng)600 m，寬500 m，后緣高程355 m，滑坡前緣直抵長(zhǎng)江河邊，剪出口位于135 m 處，在庫(kù)區(qū)蓄水水位以下，滑坡全貌如圖1所示。整個(gè)滑坡以下伏基巖為滑床，左右兩側(cè)分別以溝谷、沖溝為界，呈橫長(zhǎng)形，滑坡中部和后緣稍陡，前緣較緩，平均坡度30°，坡向150°，后緣陡壁坡度35°，高10 m?；麦w厚50 m，滑體面積30 萬(wàn)m2，體積1 500 萬(wàn)m3。

圖1 雞腦殼包滑坡全貌Fig.1 General view of Jinaokebao landslide

雞腦殼包滑坡所處的巫山縣年降雨分布圖如圖2 所示，其中強(qiáng)降雨主要集中在4 月份至10 月份，其余月份降雨相對(duì)較少。

圖2 巫山縣年降雨量Fig.2 Annual rainfall in Wushan County

1.2 計(jì)算參數(shù)

根據(jù)雞腦殼包滑坡地勘資料確定計(jì)算相關(guān)參數(shù)，具體數(shù)值見(jiàn)表1，二維計(jì)算剖面的地質(zhì)剖面圖如圖3所示。

表1 雞腦殼包滑坡數(shù)值計(jì)算參數(shù)取值表Tab.1 The value table of numerical calculation parameters of Jinaokebao landslide

圖3 滑坡剖面圖Fig.3 Landslide profile map

1.3 邊界條件

降雨和庫(kù)水位變化時(shí)滑坡滲流的邊界條件如圖4 所示，其中滑坡前緣庫(kù)水浸沒(méi)部分為水頭邊界。斜坡表面處降雨引起的入滲為流量邊界，當(dāng)降雨強(qiáng)度大于巖土體入滲速度時(shí)，邊界流量值取巖土體入滲速度值，反之則取降雨強(qiáng)度值。模型底面和兩側(cè)為自由滲透邊界且基巖的滲透性較小，可視為不滲透邊界。

圖4 邊界條件示意圖Fig.4 Schematic diagram of boundary conditions

1.4 計(jì)算模型

通過(guò)地質(zhì)剖面建立二維計(jì)算模型，利用GEO-SLOPE 有限元軟件對(duì)雞腦殼包滑坡主縱剖面進(jìn)行四邊形單元網(wǎng)格剖分，節(jié)點(diǎn)數(shù)為8 130，單元數(shù)為8 035，網(wǎng)格圖如圖5所示。

圖5 二維計(jì)算網(wǎng)格模型圖Fig.5 2D computational grid model diagram

根據(jù)雞腦殼包滑坡的實(shí)際地貌特征建立思維三維數(shù)值計(jì)算模型，取值范圍為分別平行、垂直長(zhǎng)江水流方向886、850 m，模型底面高程為0 m，滑體、滑帶、基巖均包含在計(jì)算域，該計(jì)算域剖分了36 864 個(gè)六面體單元，共計(jì)42 250 個(gè)節(jié)點(diǎn)，三維計(jì)算模型與網(wǎng)格如圖6所示。

圖6 三維計(jì)算網(wǎng)格模型圖Fig.6 3D computational grid model diagram

1.5 工況設(shè)計(jì)

三峽工程全年水位調(diào)度圖如圖7 所示。其中6 月至9 月為汛期，此期間是強(qiáng)降雨的高發(fā)期，因此三峽庫(kù)區(qū)為滿足蓄洪要求，會(huì)將庫(kù)水位從162 m 快速下降至145 m。而在10月至次年5月為非汛期，在此期間內(nèi)降雨量相對(duì)較少，庫(kù)水位變成為緩慢消落方式，水庫(kù)水位從175 m緩慢下降至145 m。所以在計(jì)算中考慮快慢兩種消落方式。

圖7 三峽水庫(kù)運(yùn)行水位曲線圖Fig.7 Chart of operating water level of Three Gorges Reservoir

由于汛期與非汛期的降雨強(qiáng)度并不一致，所以依據(jù)降雨資料，總結(jié)不同時(shí)期的20 年一遇與50 年一遇的區(qū)域降雨強(qiáng)度，如表2 所示。結(jié)合水位調(diào)動(dòng)情況，在二維計(jì)算時(shí)設(shè)置7 種工況對(duì)滑坡的穩(wěn)定性進(jìn)行分析，具體工況如表3所示，在三維分析時(shí)設(shè)置5種工況，具體工況分析如表4所示。

表2 雞腦殼包滑坡所在區(qū)域降雨強(qiáng)度 mm/dTab.2 Rainfall intensity in the area where Jinaokebao landslide is located

表3 二維計(jì)算工況Tab.3 2D calculation conditions

表4 三維計(jì)算工況Tab.4 3D calculation conditions

2 滑坡二維穩(wěn)定性分析

2.1 滑坡穩(wěn)定性計(jì)算結(jié)果

二維計(jì)算采用Morgenstern-Prince 法計(jì)算雞腦殼包滑坡的穩(wěn)定性，其中靜水時(shí)的穩(wěn)定系數(shù)為1.091。按照“三峽庫(kù)區(qū)地質(zhì)災(zāi)害防治工程地質(zhì)勘查技術(shù)要求”，當(dāng)穩(wěn)定性系數(shù)f≥1.10 為“穩(wěn)定”，介于1.05≤f＜1.10之間為“基本穩(wěn)定”，介于1.00≤f＜1.05之間為“欠穩(wěn)定”，當(dāng)穩(wěn)定性系數(shù)f＜1.00 時(shí)為“不穩(wěn)定”，所以庫(kù)水位靜止時(shí)滑坡基本穩(wěn)定。在庫(kù)水降落與不同強(qiáng)度降雨耦合情況下的滑坡穩(wěn)定系數(shù)變化曲線如圖8所示。

2.2 二維穩(wěn)定性分析

由圖8 可知在僅有庫(kù)水位變化的情況下，不論快速下降還是緩慢下降，滑坡的穩(wěn)定性都會(huì)降低，但在高落差情況（175 m降至145 m）的穩(wěn)定性明顯要低于低落差情況（162 m 降至145 m），即在僅有水位調(diào)度時(shí)，供水期滑坡會(huì)更加危險(xiǎn)。

在疊加不同強(qiáng)度降雨時(shí)，滑坡穩(wěn)定性會(huì)出現(xiàn)明顯的降低，由圖8（a）可知，在50年一遇降雨強(qiáng)度下的穩(wěn)定系數(shù)直接由無(wú)降雨的1.091 降至1.086，說(shuō)明強(qiáng)降雨在滑坡體前期就作用明顯，結(jié)合圖8（b）可知，20年的降雨條件下滑坡前期也會(huì)明顯導(dǎo)致滑坡的穩(wěn)定性降低，但隨著時(shí)間的持續(xù)，開(kāi)始與無(wú)降雨條件下的穩(wěn)定性開(kāi)始接近，所以主要影響滑坡的穩(wěn)定性的因子是水位變化和強(qiáng)降雨。

根據(jù)穩(wěn)定系數(shù)變化曲線可知工況3（庫(kù)水位從175 m水位緩降至145 m+50 年一遇暴雨）為雞腦殼包滑坡最危險(xiǎn)工況，此時(shí)由于庫(kù)水位下降和降雨影響，滑坡的穩(wěn)定系數(shù)最小值為1.032，滑坡處于欠穩(wěn)定狀態(tài)。

3 滑坡三維穩(wěn)定性分析

采用Abaqus對(duì)雞腦殼包滑坡在表4中確定的工況1～工況5的滑坡變形破壞過(guò)程進(jìn)行模擬。

3.1 孔隙水壓力計(jì)算結(jié)果分析

當(dāng)水庫(kù)蓄水至175 m 時(shí)，工況1 條件下孔隙水壓力計(jì)算結(jié)果如圖9（a）所示，此時(shí)庫(kù)水位直接影響滑坡涉水區(qū)域的孔隙水壓力。

圖9 雞腦殼包滑坡三維模型不同工況下孔隙水壓力計(jì)算結(jié)果分布圖（單位：kPa）Fig.9 Distribution diagram of pore water pressure calculation results under different working conditions of 3D model of Jinaokebao landslidte

在庫(kù)水位下降且沒(méi)有降雨影響時(shí)，工況2 和工況4 孔隙水壓力計(jì)算結(jié)果如圖9（b）、（d）所示，此時(shí)滑坡前緣地下水分布受庫(kù)水位變化直接影響，當(dāng)庫(kù)水位下降時(shí)，滑坡前緣的地下水位相應(yīng)下降，但由于滑坡前緣孔隙水壓力消散需要時(shí)間，因此存在滯后性。由于滑坡后緣離庫(kù)水位變化帶的距離相對(duì)較遠(yuǎn)，所以滑坡后緣的孔隙水壓力基本未發(fā)生變化。

在庫(kù)水位下降且遭遇50 年一遇降雨影響時(shí)，工況3 與工況5 孔隙水壓力計(jì)算結(jié)果如圖9（c）、（e）所示，此時(shí)庫(kù)水位變化和降雨聯(lián)合作用直接影響滑坡前緣地下水分布，當(dāng)庫(kù)水位下降時(shí)，地下水位也有滯后性下降。由于滑坡后緣離庫(kù)水位變化帶的距離相對(duì)較遠(yuǎn)，所以滑坡后緣的孔隙水壓力主要受降雨影響，而滑坡后緣的孔隙水壓力受降雨影響時(shí)的變化對(duì)于滑坡整體而言并不大。

由圖9 可知，涉水庫(kù)岸的孔隙水壓力受庫(kù)水位下降的影響較大，降雨主要影響滑坡表層。

3.2 應(yīng)力計(jì)算結(jié)果分析

水庫(kù)蓄水至175 m水位后，大部分區(qū)域拉、壓應(yīng)力分布變化較小。工況1 應(yīng)力分布云圖如圖10 所示，此時(shí)拉應(yīng)力主要分布在滑體下部；滑體上部壓應(yīng)力值較大，中部壓應(yīng)力值較小。

在庫(kù)水位下降且沒(méi)有降雨影響時(shí)，結(jié)合應(yīng)力最值統(tǒng)計(jì)表（表5）可知，工況2 和工況4 拉應(yīng)力值減小，壓應(yīng)力值的絕對(duì)值增大，當(dāng)庫(kù)水位由175 m 緩降至145 m 時(shí)，拉應(yīng)力最大值由938.3 kPa 減小為551.9 kPa，壓應(yīng)力最大值的絕對(duì)值由1 459 kPa 增大為1 515 kPa。拉應(yīng)力主要分布在滑體下部，滑體上部壓應(yīng)力值較大，滑體中部壓應(yīng)力值較小。當(dāng)庫(kù)水位由162 m 水位驟降到145 m 時(shí)，拉應(yīng)力主要分布在滑體下部，最大值為558.8 kPa；滑體上部壓應(yīng)力值較大，其絕對(duì)值的最大值為1 516 kPa；滑體中部壓應(yīng)力值較小。由于庫(kù)水位下降使?jié)B流場(chǎng)發(fā)生改變，通過(guò)滲流場(chǎng)與應(yīng)力場(chǎng)的耦合作用，導(dǎo)致應(yīng)力場(chǎng)的改變，這是應(yīng)力場(chǎng)發(fā)生改變的重要原因。

在庫(kù)水位下降且遭遇50年一遇降雨影響時(shí)，結(jié)合應(yīng)力最值統(tǒng)計(jì)表（表5）可知，工況3 和工況5 拉應(yīng)力值有所增大，壓應(yīng)力值的絕對(duì)值有所減小，當(dāng)庫(kù)水位由175 m 緩降到145 m 且遭遇50 年一遇降雨影響時(shí)，拉應(yīng)力最大值由551.9 kPa 增大到603.9 kPa，壓應(yīng)力最大值的絕對(duì)值由1 515 kPa 減小到1 510 kPa。拉應(yīng)力主要分布在滑體下部，滑體上部壓應(yīng)力值較大，滑體中部壓應(yīng)力值較小。當(dāng)庫(kù)水位由162 m 突降到145 m 且遭遇50年一遇降雨影響時(shí)，拉應(yīng)力最大值由558.8 kPa 增大到635.2 kPa，壓應(yīng)力最大值的絕對(duì)值由1 516 kPa 減小到1 421 kPa，拉應(yīng)力主要分布在滑體下部，滑體上部壓應(yīng)力值較大，滑體中部壓應(yīng)力值較小。在滲流場(chǎng)與應(yīng)力場(chǎng)相互耦合作用下，庫(kù)水下降和降雨的共同作用引起滲流場(chǎng)的改變，是應(yīng)力場(chǎng)發(fā)生變化的主要原因。

結(jié)合孔隙水壓力變化可知，水位變化與降雨作用導(dǎo)致滲流場(chǎng)變化，而滲流場(chǎng)與應(yīng)力場(chǎng)相互影響，因此不同工況下應(yīng)力場(chǎng)的改變主要是滲流場(chǎng)的變化。

3.3 位移計(jì)算結(jié)果分析

圖11 為工況1 條件下位移分布圖，此時(shí)滑體上部發(fā)生較大水平位移，滑體中下部水平位移相對(duì)較小。這是由于該滑坡為推移式滑坡，而推移式滑坡是從滑坡后緣發(fā)生變形而產(chǎn)生破壞。在滑體中部發(fā)生較大的向下垂直位移，滑體上部和下部向下垂直位移相對(duì)較小。這是由于滑體中部較為陡峭，使其在自重作用下產(chǎn)生較大的向下垂直位移，而上部和下部相對(duì)較為平緩，垂直位移相對(duì)較小。

匯總滑坡位移最大值得到表6，對(duì)比工況2～5 的水平位移，4 種工況下較大水平位移均發(fā)生在滑體上部，在庫(kù)水位下降且沒(méi)有降雨影響時(shí)，工況2和工況4的水平位移均大于工況1時(shí)的水平位移，這是由于庫(kù)水下降，滑體內(nèi)產(chǎn)生向外的滲透力，使滑體下滑力增大，導(dǎo)致滑體水平位移增大，但水位下降速率對(duì)水平位移影響不大。在庫(kù)水位下降且遭遇50年一遇降雨影響時(shí)，由于滑坡的滲透力進(jìn)一步增大，滑動(dòng)力也隨之增大，導(dǎo)致工況3和工況5 的水平位移大于工況2 和工況4，其中工況3 條件下滑坡的水平位移最大。

表6 滑坡位移最大值統(tǒng)計(jì)表Tab.6 Landslide displacement maximum statistical table

對(duì)比工況2～工況5 垂直位移，滑坡中部均發(fā)生較大的向下垂直位移，而滑坡上部和下部向下垂直位移相對(duì)較小。在庫(kù)水位下降且遭遇50 年一遇降雨影響時(shí)，垂直位移最為明顯，其中工況3 時(shí)垂直位移最大，這是由于庫(kù)水位下降導(dǎo)致滑體向下滑移，同時(shí)在自重和降雨作用下，使滑體產(chǎn)生的向下垂直位移增大。當(dāng)庫(kù)水位下降且沒(méi)有降雨影響時(shí)，對(duì)比工況2 和工況4 的垂直位移場(chǎng)分布，此時(shí)工況2 滑體中部的垂直位移量為4.35～26.55 mm，滑體上部垂直位移相對(duì)較小，量值為0～19.81 mm，滑體下部垂直位移較小，量值為0～13.08 mm；而工況4 滑體中部的垂直位移量值為6.55～26.44 mm，滑體上部垂直位移相對(duì)較小，量值為0～19.73 mm，滑體下部垂直位移較小，量值為0～13.02 mm。結(jié)果表明水位下降速率對(duì)滑坡的垂直位移有影響。

4 數(shù)值預(yù)測(cè)模擬

由數(shù)值計(jì)算分析結(jié)果可知工況3（庫(kù)水位由175 m水位緩降到145 m 水位+50 年一遇降雨）為雞腦殼包滑坡最不利工況。提取滑坡主滑剖面上地表位移監(jiān)測(cè)點(diǎn)WS07-04、WS07-05 和WS07-06，如圖12 所示。在最不利工況（工況3）下各監(jiān)測(cè)點(diǎn)的位移變化如圖13所示。

圖13 監(jiān)測(cè)點(diǎn)WS07-04、WS07-05、WS07-06位移變化趨勢(shì)Fig.13 Change trend of displacement in monitoring points WS07-04， WS07-05 and WS07-06

根據(jù)監(jiān)測(cè)點(diǎn)位移變化趨勢(shì)可知工況3（庫(kù)水位由175 m水位緩降到145 m 水位+50 年一遇降雨）條件下，處于滑坡中后部的地表變形監(jiān)測(cè)點(diǎn)WS07-04 和WS07-05 的位移遠(yuǎn)大于滑坡前部的地表變形監(jiān)測(cè)點(diǎn)WS07-06位移，說(shuō)明滑坡的位移主要發(fā)生在中后部。

5 庫(kù)水位下降速率、降雨與滑坡穩(wěn)定性關(guān)系

5.1 庫(kù)水位下降速率與滑坡穩(wěn)定性關(guān)系

為進(jìn)一步分析庫(kù)水位下降速率與滑坡穩(wěn)定性關(guān)系，取不同庫(kù)水位下降速率進(jìn)行滑坡滲流計(jì)算，由于二維計(jì)算結(jié)果表明工況3（庫(kù)水位從175 m 水位緩降至145 m+50年一遇暴雨）時(shí)滑坡最危險(xiǎn)，所以基于此條件下僅改變庫(kù)水位下降速率作為變量對(duì)滑坡進(jìn)行分析，計(jì)算時(shí)以1.2m/d 為初始下降速率逐漸增加，每增加0.2 m/d 為一檔，采用Morgenstern-Prince 法計(jì)算滑坡在不同庫(kù)水降速下的穩(wěn)定系數(shù)，計(jì)算結(jié)果見(jiàn)表7。

表7 不同庫(kù)水降速下滑坡穩(wěn)定系數(shù)Tab.7 Stability coefficient of downhill slope of different reservoir water slowdowns

由表7可知，隨著庫(kù)水下降速率的增大，滑坡的穩(wěn)定系數(shù)逐漸減小，庫(kù)水下降速率從1.2 m/d 增加到4.2 m/d 時(shí)，穩(wěn)定系數(shù)從1.032降至0.997，處于失穩(wěn)狀態(tài)。

從圖14 可知，隨著庫(kù)水下降速率的增大，滑坡的穩(wěn)定性逐漸降低，主要是由于庫(kù)水位下降對(duì)滑坡產(chǎn)生滲流作用，滲流場(chǎng)受降速的影響很大，庫(kù)水位下降速率越大，滲透力也越大，滑坡穩(wěn)定性系數(shù)越低。所以應(yīng)加強(qiáng)滑坡地表位移的監(jiān)測(cè)，做好災(zāi)前預(yù)報(bào)和準(zhǔn)備工作，及時(shí)防治可以有效減少及避免滑坡的發(fā)生。

圖14 不同庫(kù)水降速下滑坡穩(wěn)定系數(shù)Fig.14 Stability coefficient of downhill slope of different reservoir water slowdowns

5.2 降雨與滑坡穩(wěn)定性關(guān)系

為了分析降雨與滑坡穩(wěn)定性關(guān)系，取100 年一遇暴雨和50年一遇暴雨對(duì)應(yīng)的降雨強(qiáng)度和3 天、4 天、5 天的降雨歷時(shí)進(jìn)行組合。采用Morgenstern-Prince 法，計(jì)算滑坡在不同降雨條件下的穩(wěn)定系數(shù)，計(jì)算結(jié)果見(jiàn)表8。

表8 不同降雨條件下滑坡穩(wěn)定系數(shù)Tab.8 Landslide stability coefficient under different rainfall conditions

由表8可知，隨著降雨時(shí)間和降雨強(qiáng)度的不斷增加，滑坡的穩(wěn)定系數(shù)不斷減??；當(dāng)降雨過(guò)程為重現(xiàn)期100 年、歷時(shí)3 天（降雨強(qiáng)度為78.315 mm/d）時(shí)，滑坡穩(wěn)定系數(shù)為0.975，處于失穩(wěn)狀態(tài)；當(dāng)降雨過(guò)程為重現(xiàn)期50 年、歷時(shí)4 天（降雨強(qiáng)度為64.535 mm/d）時(shí)，滑坡穩(wěn)定性系數(shù)降至0.968，處于失穩(wěn)狀態(tài)。

6 結(jié) 論

（1）根據(jù)二維穩(wěn)定性分析可知，主要影響滑坡的穩(wěn)定性的因子是水位變化和強(qiáng)降雨，其中庫(kù)水位從175 m 緩降至145 m且遭遇50 年一遇暴雨時(shí)，為雞腦殼包滑坡最危險(xiǎn)工況，此時(shí)滑坡的穩(wěn)定系數(shù)最小值為1.032，滑坡處于欠穩(wěn)定狀態(tài)。

（2）根據(jù)三維變形破壞分析可知，庫(kù)水位下降與強(qiáng)降雨造成滑體內(nèi)滲流場(chǎng)發(fā)生很大變化，使坡體內(nèi)的滲透力發(fā)生很大的改變，進(jìn)而引起位移發(fā)生很大的變化；雞腦殼包滑坡在庫(kù)水位由175 m 水位緩降至145 m 水位并且疊加50年一遇降雨條件下位移最大。因此，工況3為雞腦殼包滑坡最不利工況，此時(shí)滑坡位移主要發(fā)生在滑坡中后部。

（3）庫(kù)水位下降速率與滑坡穩(wěn)定性關(guān)系的分析表明，滑坡的穩(wěn)定系數(shù)隨著庫(kù)水降速的增大而減小，當(dāng)庫(kù)水下降速率為4.2 m/d 時(shí)，滑坡的穩(wěn)定性系數(shù)為0.997，此時(shí)滑坡處于失穩(wěn)狀態(tài)。降雨與滑坡穩(wěn)定性關(guān)系的分析表明，滑坡的穩(wěn)定系數(shù)隨著降雨強(qiáng)度和降雨時(shí)長(zhǎng)的增大而減小。因此，庫(kù)水位下降速率越大、降雨強(qiáng)度越大、降雨時(shí)長(zhǎng)越長(zhǎng)，滑坡越容易失穩(wěn)。

（4）綜合以上分析，庫(kù)水位下降和降雨對(duì)雞腦殼包滑坡的穩(wěn)定性影響較大。鑒于滑坡巖土體物理力學(xué)參數(shù)在經(jīng)歷多次水庫(kù)水位波動(dòng)和降雨過(guò)程會(huì)發(fā)生變化的現(xiàn)實(shí)，建議加強(qiáng)觀測(cè)，并根據(jù)觀測(cè)結(jié)果及時(shí)反演更新滑坡物理力學(xué)參數(shù)。以此為基礎(chǔ)，對(duì)滑坡在后續(xù)蓄水和降雨條件下的變形破壞過(guò)程進(jìn)行及時(shí)預(yù)測(cè)分析，為該滑坡的監(jiān)測(cè)預(yù)警提供合理可靠的依據(jù)。