楊 鵬,鄭 光,黃 金 成,劉 震 東,陳 明 浩

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

0 引 言

金沙江上游受構(gòu)造破碎、巖體結(jié)構(gòu)復(fù)雜、河谷深切、強(qiáng)卸荷改造等多種不良因素的疊合[1],區(qū)內(nèi)地質(zhì)災(zāi)害頻發(fā)[2-3],如2018年金沙江上游西藏江達(dá)縣白格村先后發(fā)生2次大型高位滑坡-堰塞堵江事件,造成巨大損失[4-5]。若修建水庫(kù),將促進(jìn)該區(qū)域地質(zhì)災(zāi)害的發(fā)育,威脅水庫(kù)正常運(yùn)行以及當(dāng)?shù)鼐用癜踩?。?jù)金沙江水電規(guī)劃報(bào)告,金沙江作為中國(guó)最大的水電基地和“西電東送”主力,其上游將進(jìn)行12級(jí)梯級(jí)水電站建設(shè)。旭龍水庫(kù)位于金沙江干流上游河段,水庫(kù)建成蓄水后可能發(fā)生庫(kù)岸變形破壞等的工程地質(zhì)問題[6-7],加之該區(qū)域地質(zhì)環(huán)境差異大、新構(gòu)造運(yùn)動(dòng)強(qiáng)烈、生態(tài)環(huán)境脆弱,若滑坡失穩(wěn)下滑,將造成巨大損失和危害。

在國(guó)外,由于水庫(kù)蓄水導(dǎo)致滑坡失穩(wěn)的典型實(shí)例如瓦依昂水庫(kù)滑坡[8],造成了災(zāi)難性后果,使大量學(xué)者開始關(guān)注此類滑坡。根據(jù)國(guó)內(nèi)關(guān)于三峽庫(kù)區(qū)典型滑坡的研究成果,水庫(kù)蓄水后對(duì)滑坡的影響可主要?dú)w結(jié)為滲流特性研究和滑坡穩(wěn)定性評(píng)價(jià)2個(gè)方面。比如大部分學(xué)者對(duì)三峽水庫(kù)正常運(yùn)行期間,不同的庫(kù)水升降速率、滑體滲透系數(shù)對(duì)滑坡穩(wěn)定性的影響進(jìn)行研究[9-10],得出了一些影響穩(wěn)定性系數(shù)的因素和規(guī)律。也有學(xué)者對(duì)滑坡進(jìn)行流-固耦合的數(shù)值模擬計(jì)算[11-14],研究分析了滑坡的滲流特征及穩(wěn)定性變化。在蓄水條件下,學(xué)者多關(guān)注某一水位或變化水位后的滑坡穩(wěn)定性,但是對(duì)于水庫(kù)蓄水正常運(yùn)行后滑坡穩(wěn)定性的動(dòng)態(tài)變化以及滑坡變形預(yù)警模型的研究較少。坡體變形是直接反映滑坡穩(wěn)定性和變化特征的重要參考。近年來,基于滑坡變形監(jiān)測(cè)技術(shù),采用滑坡變形切線角預(yù)警模型[15],對(duì)貴州省龍井村滑坡、甘肅黑方臺(tái)黃土滑坡[16]進(jìn)行了多次成功預(yù)警。滑坡變形切線角模型的實(shí)時(shí)過程化、無量綱化特點(diǎn)適用于幾乎所有的漸變型或突發(fā)型滑坡,并取得了良好的預(yù)警防災(zāi)效果。利用該模型的普遍適用性,結(jié)合坡體穩(wěn)定性系數(shù)的動(dòng)態(tài)變化,可對(duì)庫(kù)水位變化下的坡體變形階段進(jìn)行預(yù)測(cè)判別。

本文通過數(shù)值計(jì)算得到水庫(kù)蓄水后旭龍水電站庫(kù)區(qū)格亞頂滑坡堆積體的滲流特征、應(yīng)力位移及穩(wěn)定性系數(shù)動(dòng)態(tài)變化,結(jié)合滑坡變形切線角預(yù)警判據(jù),對(duì)水庫(kù)蓄水后該坡體的變形階段進(jìn)行預(yù)測(cè)判別,對(duì)于旭龍水庫(kù)地質(zhì)災(zāi)害防治具有重要意義。

1 滑坡概況

旭龍水電站水庫(kù)正常蓄水位2 302 m,總庫(kù)容8.47億m3,具有日調(diào)節(jié)能力。格亞頂堆積體位于金沙江干流上游河段右岸,屬高山峽谷地貌,距離壩址約3.8 km。如圖1所示,滑坡平面形態(tài)呈不規(guī)則形態(tài),主滑方向約53°,平均坡度約30°。后緣高程約2 750 m,平均坡度約為47°,為崩塌形成的陡壁,兩側(cè)山脊有斷層F6-1穿過,基巖裂隙發(fā)育,前緣直抵金沙江江面,高程約2 183 m。

格亞頂堆積體存在Ⅰ區(qū)滑坡堆積體和Ⅱ區(qū)崩塌堆積體兩部分。Ⅰ區(qū)堆積平臺(tái)位于坡體下部,為滑坡堆積體平臺(tái),高程約2 230 m。Ⅰ區(qū)組成物質(zhì)主要為灰黑色粗粒土,巖性為花崗巖、變質(zhì)巖碎石、塊石,巖質(zhì)新鮮、堅(jiān)硬,粒徑為0.1～0.5 m,見圖2(a)。Ⅱ區(qū)堆積平臺(tái)位于坡體中部、Ⅰ區(qū)平臺(tái)左后側(cè),為上部的崩滑堆積體組成,高程約2 300 m。Ⅱ區(qū)物質(zhì)組成為黃褐色塊、灰黑色碎石土,主要為碎石、角礫及礫砂等,碎石巖性為石英閃長(zhǎng)巖,間夾云母片巖,見圖2(b)。水庫(kù)建成蓄水后,正常蓄水位將達(dá)到坡體中下部,淹沒到坡體下滑段,如圖3所示?？紤]水庫(kù)30 m的波動(dòng)水位,坡體受動(dòng)水壓力作用,將產(chǎn)生變形甚至失穩(wěn)下滑堵江,威脅下游水庫(kù)、村落和城鎮(zhèn)。動(dòng)水壓力作用指水在土體孔隙中滲流時(shí),對(duì)其周圍骨架產(chǎn)生的滲透力,其大小與水力梯度成正比。如果動(dòng)水壓力指向坡外,將降低滑坡體的穩(wěn)定性,從而引起滑坡復(fù)活[17]。本文將利用GeoStudio軟件對(duì)格亞頂堆積體進(jìn)行流-固耦合分析,預(yù)測(cè)其在庫(kù)水位變動(dòng)條件下的變形及穩(wěn)定性。

圖3 格亞頂堆積體1-1′剖面Fig.3 1-1′ profile of Geyading accumulation body

2 庫(kù)水位變動(dòng)下滑坡數(shù)值模型

選用GeoStudio中的SEEP/W、SIGMA/W、SLOPE/W等3種模塊相互結(jié)合的方式進(jìn)行滲流-應(yīng)力耦合分析,對(duì)蓄水后的格亞頂堆積體進(jìn)行變形及穩(wěn)定性分析。

2.1 模型構(gòu)建及邊界條件

根據(jù)工程地質(zhì)資料以及水文資料,本文選取1-1′地質(zhì)剖面建立數(shù)值模型,如圖4所示。設(shè)計(jì)單元尺寸為10 m,共計(jì)3 699個(gè)節(jié)點(diǎn),3 500個(gè)單元,并設(shè)置合理的邊界條件。具體邊界設(shè)置如下:① 初始滲流邊界:將前緣坡面初始水位高程設(shè)置為天然江面2 183 m,底部為零流量邊界,兩側(cè)邊界為自由滲流邊界,計(jì)算初始滲流狀態(tài)。② 瞬態(tài)滲流邊界:將前緣坡面設(shè)置為動(dòng)態(tài)水頭邊界,水位以一定的速率升降,以此進(jìn)行瞬態(tài)滲流分析。③ 位移邊界:約束模型兩側(cè)的水平位移、模型底部的水平和豎向位移。④ 監(jiān)測(cè)點(diǎn):為了便于分析庫(kù)水位變動(dòng)過程中滑坡失穩(wěn)的內(nèi)在規(guī)律,設(shè)置監(jiān)測(cè)點(diǎn) A～D進(jìn)行分析,其水位高程分別為2 230,2 272,2 302 m和2 320 m。

圖4 計(jì)算模型Fig.4 Calculation model

2.2 參數(shù)選取以及計(jì)算工況

假定模型介質(zhì)均為理想彈塑性材料,服從莫爾-庫(kù)倫強(qiáng)度準(zhǔn)則,材料的物理力學(xué)參數(shù)是通過室內(nèi)試驗(yàn)、現(xiàn)場(chǎng)試驗(yàn)及工程類比等方法確定,物理力學(xué)性質(zhì)參數(shù)見表1。運(yùn)用SEEP/W模塊自帶的樣本函數(shù)擬合得出巖土體的土水特征曲線(見圖5);利用飽和滲透系數(shù)與土水特征曲線,結(jié)合Van Genuchten模型推導(dǎo)出滲透系數(shù)曲線(見圖6)?？紤]水庫(kù)運(yùn)行調(diào)度及庫(kù)水位變幅特征變動(dòng)范圍為30 m,即2 183 m→2 302 m→2 272 m的庫(kù)水位變動(dòng)對(duì)滑坡變形及穩(wěn)定性的影響。庫(kù)水位以1 m/d的速率升降,動(dòng)態(tài)邊界水頭函數(shù)如圖7所示。

表1 巖土體物理力學(xué)參數(shù)取值Tab.1 Physical and mechanical parameters of rock and soil mass

圖5 土水特征曲線Fig.5 Soil-water characteristic curve

圖6 滲透系數(shù)曲線Fig.6 Permeability coefficient curve

圖7 動(dòng)態(tài)邊界水頭函數(shù)Fig.7 Dynamic boundary head function

3 數(shù)值模擬結(jié)果分析

3.1 滲流場(chǎng)數(shù)值分析

由于水庫(kù)運(yùn)行,岸坡外的水位周期性變化較大,岸坡土體的飽和度也處于變化中,水分的相互補(bǔ)給改變了岸坡內(nèi)的滲流場(chǎng)。斜坡巖土體的不均勻性和庫(kù)水的滲透作用將直接影響岸坡的穩(wěn)定性,所以,對(duì)庫(kù)水作用下的滲流規(guī)律進(jìn)行分析具有重要意義。

利用 GeoStudio 軟件的后處理功能,得到瞬態(tài)滲流分析中地下水浸潤(rùn)線動(dòng)態(tài)變化特征,如圖8所示。浸潤(rùn)線指的是地下水位線中孔隙水壓力為0的線條。各階段浸潤(rùn)線變化特征如下:

1.1 材料 ①一次性無紡布，規(guī)格120 cm×120 cm、60 cm×60 cm各60塊;②雙層棉布，規(guī)格120 cm×120 cm、60 cm×60 cm各60塊(在包布一角訂一塊3 cm×3 cm的“補(bǔ)丁”，用于記錄使用次數(shù));③3M壓力蒸汽滅菌化學(xué)指示卡及指示膠帶。

圖8 水位升降過程中岸坡地下水浸潤(rùn)線動(dòng)態(tài)變化Fig.8 Dynamic changes of groundwater infiltration line on bank slope during the rise and fall of water level

(1) 在水位上升階段,如圖8(a)所示,水位在上升過程中,因?yàn)槠麦w滲透速率明顯小于庫(kù)水升降速率,導(dǎo)致坡體內(nèi)浸潤(rùn)線的位置明顯滯后于庫(kù)水位,上升時(shí)浸潤(rùn)線在坡面附近向下彎折,隨著庫(kù)水上升,這種滯后性在坡體中越明顯,說明在水位上升過程中,水由坡面向坡體內(nèi)逐漸補(bǔ)給。

(2) 在2 302 m水位維持階段,如圖8(b)所示,浸潤(rùn)線逐漸緩慢上升,明顯還存在一定滯后性,這種現(xiàn)象在巖土體的滲透系數(shù)足夠大時(shí),將不再明顯。與剛上升至2 302 m水位時(shí)相比,209 d浸潤(rùn)線上升明顯。這種上升幅度隨時(shí)間逐漸減小,并且在局部出現(xiàn)彎折重合的現(xiàn)象。

(3) 在水位下降階段,如圖8(c)所示,浸潤(rùn)線下降幅度隨時(shí)間逐漸減小,在209～239 d期間,坡體浸潤(rùn)線向下彎折幅度較大,說明浸潤(rùn)線在水位下降階段變化較大。

(4) 在2 272 m水位維持階段,如圖8(d)所示,浸潤(rùn)線在水位維持的90 d中逐漸緩慢下降,在坡體中浸潤(rùn)線的滯后現(xiàn)象已不明顯。水位變動(dòng)階段相比于水位維持階段,浸潤(rùn)線變化幅度大,這也使得應(yīng)力在水位變動(dòng)階段出現(xiàn)明顯變化(見圖9～10)。

圖9 監(jiān)測(cè)點(diǎn)最大剪應(yīng)力變化特征Fig.9 Variation characteristics of maximum shear stress at monitoring points

圖10 監(jiān)測(cè)點(diǎn)最大有效應(yīng)力變化特征Fig.10 Variation characteristics of maximum effective stress at monitoring points

3.2 坡體變形分析

庫(kù)水位升降引起坡體內(nèi)部應(yīng)力的變化,繼而產(chǎn)生位移的變化。根據(jù)設(shè)置的監(jiān)測(cè)點(diǎn),庫(kù)水位在升降過程中不同位置的最大剪應(yīng)力、最大有效應(yīng)力和位移變化如圖9～11所示。各階段應(yīng)力-位移變化特征如下:

(1) 在水位上升階段,各監(jiān)測(cè)點(diǎn)的最大剪應(yīng)力都呈急劇上升的趨勢(shì);最大有效應(yīng)力在庫(kù)水位未達(dá)到監(jiān)測(cè)點(diǎn)高程前變化不明顯,在水位上升至監(jiān)測(cè)點(diǎn)高程附近時(shí),最大有效應(yīng)力開始急劇下降,如A,B,C點(diǎn),而D點(diǎn)因高程位于庫(kù)水位2 302 m之上,受庫(kù)水作用相比于其他點(diǎn)較小,變化不明顯。這說明在蓄水過程中,水位上升對(duì)附近巖土體應(yīng)力場(chǎng)影響較大。

對(duì)于監(jiān)測(cè)點(diǎn)的位移,各監(jiān)測(cè)點(diǎn)位在庫(kù)水?dāng)D壓作用下首先表現(xiàn)出向坡體內(nèi)部的位移,隨后在庫(kù)水位升高至監(jiān)測(cè)點(diǎn)附近時(shí),開始出現(xiàn)向坡外的位移。這是因?yàn)殡S著庫(kù)水位升高,坡體受到指向坡外的動(dòng)水壓力作用,產(chǎn)生指向坡外方向的位移,也使豎直位移表現(xiàn)出持續(xù)增長(zhǎng)的趨勢(shì)。根據(jù)圖11所示,在庫(kù)水位上升階段,坡體位移發(fā)生明顯突變,在蓄水后可能在監(jiān)測(cè)點(diǎn)高程附近出現(xiàn)局部失穩(wěn)、滑塌。

注:X和Y方向的位移以朝向坡外為正。圖11 監(jiān)測(cè)點(diǎn)位移隨時(shí)間變化Fig.11 Change of monitoring points displacement with time

(3) 在水位維持階段,坡體內(nèi)應(yīng)力變化不明顯,逐漸趨于穩(wěn)定,而各點(diǎn)位移則隨著庫(kù)水位維持一定時(shí)間后出現(xiàn)下降,并逐漸穩(wěn)定在一定范圍。這說明位移主要發(fā)生在庫(kù)水位上升和下降階段,而在水位維持階段的變形較小。

3.3 坡體穩(wěn)定性系數(shù)變化

(1) 庫(kù)水上升并維持水位對(duì)坡體的影響具體可分為2個(gè)方面:① 水位上升階段,坡體內(nèi)的水位上升滯后于庫(kù)水位(見圖8(a))從而形成水位差,水向坡體內(nèi)入滲,產(chǎn)生沿坡內(nèi)方向的動(dòng)水壓力和孔隙靜水壓力,相當(dāng)于施加在坡體表面的一個(gè)正應(yīng)力,產(chǎn)生向坡內(nèi)的位移(見圖11),并且?guī)焖灰韵虏糠质艿剿母×ψ饔?下滑力減小,從而使穩(wěn)定性系數(shù)上升,成為坡體穩(wěn)定性的有利因素,如圖12所示;② 水位維持階段,水向坡體內(nèi)入滲,土顆粒間的膠結(jié)物質(zhì)被溶解,使土體黏聚力降低,土粒間的潤(rùn)滑作用增強(qiáng)使內(nèi)摩擦角減小,即水的軟化作用使土的抗剪強(qiáng)度降低,與基巖接觸面(滑面)有效應(yīng)力減小,抗滑阻力降低,成為對(duì)坡體的不利因素,使穩(wěn)定性系數(shù)下降。

圖12 滑坡穩(wěn)定性系數(shù)與時(shí)間的關(guān)系Fig.12 Relationship between slope stability coefficient and time

(2) 庫(kù)水下降并維持水位對(duì)坡體的影響。在水位下降階段,由于巖土體的持水性主要與顆粒大小有關(guān),顆粒越小的巖土體滯留的水越多,從而使浸潤(rùn)線的下降滯后于庫(kù)水位,坡體內(nèi)具有較高的孔隙水壓且不能及時(shí)消散,形成非穩(wěn)態(tài)滲流場(chǎng),對(duì)周圍土體骨架產(chǎn)生向臨空面及向下的滲流動(dòng)水壓力(動(dòng)態(tài)擴(kuò)張力),同時(shí),由水頭差產(chǎn)生的靜水壓力相當(dāng)于對(duì)前緣施加了一個(gè)推力,使坡體穩(wěn)定性急劇下降,水位下降至2 272 m時(shí)穩(wěn)定性系數(shù)減小至最低值0.951,此時(shí)最危險(xiǎn)滑移面位置如圖13所示。在2 272 m水位維持階段,穩(wěn)定性系數(shù)因?yàn)榻?rùn)線逐漸趨于穩(wěn)定,形成穩(wěn)態(tài)滲流場(chǎng),并且應(yīng)力位移變化較小,使穩(wěn)定系數(shù)出現(xiàn)一定上升。

圖13 最危險(xiǎn)滑移面位置(239 d)(孔隙水壓力單位:kPa)Fig.13 Location of the most dangerous slip surface (239 d)

在水庫(kù)蓄水運(yùn)行中,穩(wěn)定性系數(shù)急劇下降的階段并未出現(xiàn)在水位上升階段,而是發(fā)生在2 302 m水位維持階段和下降階段,且穩(wěn)定性系數(shù)最小的時(shí)間點(diǎn)出現(xiàn)在水位下降階段。因此,在旭龍水電站蓄水運(yùn)行后,應(yīng)該關(guān)注蓄水初期以及庫(kù)水下降過程中的滑坡變形情況,及時(shí)制定相應(yīng)的災(zāi)害防治措施。

3.4 位移切線角預(yù)警判據(jù)

坡體從變形產(chǎn)生到最終失穩(wěn)破壞,累計(jì)位移(S)-時(shí)間(t)曲線一般會(huì)經(jīng)歷初始(期)變形階段、等(勻)速變形階段和加速變形階段,即所謂的斜坡變形演化3階段規(guī)律[18]。位移切線角是指S-t曲線中某一時(shí)刻變形曲線的切線與橫坐標(biāo)軸之間的夾角,即S-t曲線上某一時(shí)刻用角度表示的曲線斜率。坡體在失穩(wěn)破壞前的位移切線角一般為89.0°～89.5°,并以此作為滑坡的預(yù)警預(yù)報(bào)判據(jù)之一,即位移切線角預(yù)警判據(jù)。許強(qiáng)等[15]通過對(duì)斜坡累計(jì)位移-時(shí)間曲線進(jìn)行坐標(biāo)變換,實(shí)現(xiàn)縱橫坐標(biāo)同量綱化得到改良后的累積位移曲線。為滿足預(yù)警需要,許強(qiáng)等對(duì)滑坡的加速變形階段進(jìn)行了細(xì)分,并給出了基于變形的滑坡4級(jí)綜合預(yù)警判據(jù),如表2所示[19]。

表2 滑坡預(yù)警級(jí)別的定量劃分標(biāo)準(zhǔn)Tab.2 Quantitative classification criteria of landslide warning levels

改良S-t曲線處理具體方法如下[15]:① 通過位移累積曲線并結(jié)合坡體變形跡象,確定坡體等速變形階段的變形速率,利用公式(1)得到改良后的位移-時(shí)間變形曲線;② 通過公式(2)計(jì)算出單位時(shí)間間隔的切線角。

(1)

式中:S(i)為某一單位時(shí)間段(一般采用1個(gè)監(jiān)測(cè)周期,如1 d、7 d等)內(nèi)坡體累計(jì)位移量;v為等速變形階段的位移速率;T(i)為改良后與時(shí)間相同量綱的縱坐標(biāo)值。

(2)

式中:αi為改進(jìn)的切線角;ti為某一監(jiān)測(cè)時(shí)刻;Δt為與計(jì)算S時(shí)對(duì)應(yīng)的單位時(shí)間段(一般采用1個(gè)監(jiān)測(cè)周期,如1 d、7 d等);ΔT為單位時(shí)間段內(nèi)T(i)的變化量。

根據(jù)數(shù)值模擬得到的庫(kù)水位變動(dòng)下坡體的穩(wěn)定性系數(shù)變化規(guī)律(見圖12),選取庫(kù)水作用帶的監(jiān)測(cè)點(diǎn)C的位移數(shù)據(jù),利用單位時(shí)間步長(zhǎng)記錄的X,Y方向位移量,計(jì)算出單位時(shí)間間隔的切線角,如圖14所示。

圖14 點(diǎn)C的累積位移與切線角曲線Fig.14 Cumulative displacement and tangent angle curve at point C

利用位移切線角預(yù)警判據(jù),根據(jù)滑坡預(yù)警級(jí)別的定量劃分標(biāo)準(zhǔn)(見表2),判別坡體在蓄水后不同庫(kù)水位作用階段的變形狀態(tài)。在水庫(kù)蓄水運(yùn)行后,切線角迅速增加,達(dá)到一定水位后開始下降,并且在庫(kù)水維持階段逐漸趨于平緩穩(wěn)定。這說明在庫(kù)水上升作用下,坡體變形速率突增,切線角明顯增加。在庫(kù)水下降階段,切線角也有明顯的突增,說明坡體在庫(kù)水升降過程中變形較強(qiáng)烈,而在水位維持階段則變形較小,與前文的應(yīng)力-位移分析相符。該坡體在蓄水初期變形強(qiáng)烈,預(yù)測(cè)變形階段為初加速變形,預(yù)警級(jí)為警示級(jí),見表3。

表3 滑坡變形階段預(yù)測(cè)Tab.3 Prediction of landslide deformation stage

4 結(jié) 論

(1) 庫(kù)水位變動(dòng)影響著滑坡體的滲流場(chǎng)。在庫(kù)水位上升階段,庫(kù)水向坡體內(nèi)入滲,地下水位滯后于庫(kù)水位,浸潤(rùn)線向坡內(nèi)彎折凹陷;在庫(kù)水位下降階段,水由滑坡體向坡外排泄,浸潤(rùn)線向下彎折,并且下降幅度逐漸減小;在水位維持階段,浸潤(rùn)線隨時(shí)間趨于平緩,變化幅度逐漸減小,并出現(xiàn)局部彎折重合。

(2) 水庫(kù)蓄水運(yùn)行過程中坡體的應(yīng)力-位移變化規(guī)律:最大剪應(yīng)力與最大有效應(yīng)力受庫(kù)水位升降影響,在坡體內(nèi)變化明顯,突變主要發(fā)生在庫(kù)水上升階段。隨著庫(kù)水位上升,水平位移先增加后減小,豎向位移則逐漸增加。對(duì)水平位移影響較大的主要是水位上升階段,而對(duì)豎直位移影響較大的主要是水位升降階段。

(3) 水庫(kù)蓄水運(yùn)行過程中,坡體的穩(wěn)定性變化規(guī)律:在庫(kù)水上升階段,因指向坡內(nèi)方向的水壓力作用和浮力作用,坡體穩(wěn)定性系數(shù)明顯上升;在水位維持階段,因?yàn)槭艿綆?kù)水的軟化作用,巖土體抗滑力下降,穩(wěn)定性系數(shù)減小;在庫(kù)水位下降階段穩(wěn)定性系數(shù)繼續(xù)減小至最低值0.951,而在2 272 m水位維持階段則出現(xiàn)一定回升。

(4) 根據(jù)穩(wěn)定性系數(shù)變化規(guī)律,結(jié)合數(shù)值計(jì)算獲取的位移曲線,通過累加法計(jì)算出累積位移曲線,并利用切線角預(yù)警判據(jù),對(duì)該坡體在庫(kù)水位變動(dòng)下的變形階段進(jìn)行預(yù)測(cè)判定。結(jié)果表明:坡體變形主要發(fā)生在蓄水初期,變形階段為初加速變形,預(yù)警級(jí)為警示級(jí)。