庫水位驟降條件下某土質(zhì)滑坡滲透穩(wěn)定性分析

江趙平 謝東亞

摘要：為研究庫水位驟降條件下某滑坡體深層與淺層的滲透特性以及相應的穩(wěn)定性，根據(jù)非飽和滲流原理，利用有限元軟件Geo-slope對該滑坡體在降雨和庫水位驟降工況下的孔壓變化以及安全系數(shù)變化規(guī)律進行了模擬。計算結(jié)果表明：庫水位下降速率越大，滑坡體內(nèi)部孔壓下降幅度與下降速率越大，但滑坡體上部的孔壓下降幅度較下部要小;滑坡體上部淺層安全系數(shù)隨著時間逐漸增大，下部淺層和深層滑坡安全系數(shù)呈現(xiàn)先減小后增大的趨勢;滑坡體下部淺層發(fā)生失穩(wěn)，然后引起深層整體滑動，而上部淺層較為穩(wěn)定。研究成果可為庫水位驟降條件下滑坡失穩(wěn)機理的認識提供一定參考。

關鍵詞：滑坡穩(wěn)定性;庫水位驟降;土質(zhì)滑坡;深層滑坡;淺層滑坡

中圖法分類號：TU43

文獻標志碼：A

庫水位變動是影響滑坡穩(wěn)定性的重要因素，滑坡在庫水位作用下的失穩(wěn)機理主要體現(xiàn)在以下幾個方面：①庫水位驟降導致滑坡上的附加水壓突然卸載，引發(fā)滑坡失穩(wěn);②庫水位驟降導致滑坡內(nèi)部水位線出現(xiàn)“滯后”效應，指向滑坡體外部的滲流力加劇滑坡失穩(wěn)幾率?；率Х€(wěn)后的危害主要體現(xiàn)在以下幾點：①對庫區(qū)周圍居民的生命財產(chǎn)安全造成危害;②減小水庫有效庫容，影響水庫正常穩(wěn)定運行;③對滑坡周圍的水工建筑物造成一定危害。因此，針對庫水位驟降下的滑坡滲透特性以及穩(wěn)定性研究顯得很有必要。

針對庫水位驟降下的滑坡穩(wěn)定性研究，國內(nèi)外學者進行了很多有益探討，如魯芃等對白家包滑坡不同庫水位日降幅條件下滑坡滲流場及應力應變進行了數(shù)值模擬;郭子正利用灰色模型對三峽庫區(qū)某堆積滑坡體進行了庫水位驟降下的安全性評價;鐘啟明基于非飽和理論計算了某黏土心墻壩在庫水位驟降情況下的上游壩坡穩(wěn)定系數(shù)。然而，這些研究均未涉及滑坡在庫水位驟降下的深層淺層滑動面穩(wěn)定性分析，事實上，對于滑坡體失穩(wěn)形式來說，在降雨工況下已有研究表明，滑坡過程表現(xiàn)為先發(fā)生前期的淺層失穩(wěn)然后是后期的整體滑動，然而對于庫水位驟降下深層、淺層滑坡的研究較少，對相應的規(guī)律還沒有一個清晰的認識。

本文以恩施州鶴峰縣境內(nèi)坪江電站滑坡體為研究對象，利用Geo-slope軟件的Seep/w和Slope/w對滑坡在庫水位驟降條件下深層淺層滑坡滲流特性以及穩(wěn)定性進行了研究，為相關工程治理提供了一定的參考。

1計算理論

1.1非飽和滲流理論

非飽和滲流的控制方程寫成張量的形式為式中，Kij為飽和滲透張量;Kr為相對透水率;hc為壓力水頭;Q為源匯項;C（hc）為容水度;0為與壓力水頭相關的函數(shù);n為土體內(nèi)部的孔隙率;Ss為單位貯水量。

土水特征曲線是衡量土體內(nèi)部滲透系數(shù)（體積含水量）與基質(zhì)吸力之間的關系，比較廣泛使用的是Fredlund&Xing模型，體積含水量函數(shù)可以表達為其中，修正函數(shù)Cw采用下式計算：式中，Ow為體積含水量;0s為飽和體積含水量，m與殘余含水量相關;n與體積含水量函數(shù)拐點斜率相關;a為與體積含水量拐點處的基質(zhì)吸力;w為基質(zhì)吸力;w為殘余含水量對應的基質(zhì)吸力。

得出體積含水量之后，根據(jù)滲透系數(shù)函數(shù)，便可得出滲透系數(shù)隨基質(zhì)吸力的變化曲線：式中，kw為土體內(nèi)部對應與不同基質(zhì)吸力下的實時滲透系數(shù);ks為土體內(nèi)部飽和滲透系數(shù);y代表負孔隙水壓力算法的虛擬變量;？為最小負孔隙水壓力;N為最終函數(shù)所描述的最大負孔隙水壓力。

1.2非飽和抗剪強度理論

土體強度理論采用Fredlund雙應力變量公式：式中，t為土體實時抗剪強度;c'與0'為有效強度參數(shù);0n為法向總應力與孔隙氣壓力的差值;ua為孔隙空氣壓力;uw為孔隙水壓力;0b為由負孔隙水壓力提高的強度。

2計算模型及參數(shù)

2.1計算模型

某滑坡體位于恩施州鶴峰縣坪江水電站庫岸，滑坡坡度約為200～30°，滑坡體長約83.7m，平均厚度11.3m。該庫區(qū)死水位高程1095m，正常蓄水位高程1125m，地質(zhì)剖面圖如圖1（a）所示，滑坡典型剖面圖如圖1（b）所示，相應的網(wǎng)格剖分圖如圖1（c）所示。全局單元尺寸約為8m，整個模型一共剖分為1988個節(jié)點，1901個單元。

2.2初始條件以及邊界條件

（1）初始條件。模型左側(cè)邊界為1140m定水頭邊界，右側(cè)為175m的定水邊界，計算所得的滲流場為整個計算工況的初始條件。

（2）邊界條件。ef為定水頭邊界;abc為水位變動邊界;如為降雨人滲的流量邊界;af，de為不透水邊界;邊界上流量Q為0。

2.3土體參數(shù)及計算工況

根據(jù)室內(nèi)試驗和工程地質(zhì)勘測結(jié)果綜合確定土體參數(shù)，滑體和滑床的物理力學參數(shù)見表1。

土水特征曲線如圖2所示。

為研究不同速率庫水位驟降下的滑坡體深層淺層滲透穩(wěn)定性規(guī)律，特取庫水位下降速率為0.5，1m/d和1.5m/d這3種工況進行數(shù)值模擬研究。

3計算結(jié)果分析

3.1孔壓變化規(guī)律

為定量研究庫水位驟降下滑坡體內(nèi)不同位置孔壓變化規(guī)律，特取如圖l（a）所示的兩個監(jiān)測斷面，即上部監(jiān)測斷面與下部監(jiān)測斷面，實時監(jiān)測庫水位驟降過程中的孔壓變化，不同監(jiān)測面的孔壓變化規(guī)律如圖3—5所示。

由圖可見：孔壓隨深度呈現(xiàn)線性變化。隨著庫水位的降低，坡體內(nèi)部孑L壓逐漸降低，工況l上部孔壓在計算時間內(nèi)平均降低0.06%，下部孔壓平均降低32.61%。滑坡體上部孔壓變化幅度要遠小于下部孔壓，這也就解釋了以往文獻中庫水位驟降情況下坡內(nèi)浸潤線出現(xiàn)“延遲”的現(xiàn)象，即距離岸邊浸潤線下降幅度較大而遠離岸邊浸潤線變化幅度較小?？v向來看，工況1～3上部監(jiān)測點孔壓在計算時間內(nèi)平均降低幅度分別為0.061%，0.068%和0.075%，而下部監(jiān)測點在計算時間內(nèi)的平均降幅分別為32.61%，35.87%和38.83%，可見庫水位降低速率越大，滑坡體內(nèi)部的孔壓下降幅度也越大。

3.2滑坡穩(wěn)定性規(guī)律

不同庫水位驟降速率下的安全系數(shù)變化規(guī)律見圖6。

由圖可見，①上部淺層滑坡隨著庫水位下降，安全系數(shù)呈現(xiàn)不斷升高趨勢。這是因為滑坡體上部孔壓變化幅度較小，隨著庫水位的降低，上部地下水位線逐漸降低，土體強度不斷增加，從而安全系數(shù)呈現(xiàn)逐漸增大趨勢。同時，庫水位下降速率越大，安全系數(shù)越小。②下部淺層滑動面安全系數(shù)變化規(guī)律與深層滑動面安全系數(shù)變化規(guī)律較為一致，庫水位下降速率越大，最小安全系數(shù)出現(xiàn)的時刻越早，同時，最小安全系數(shù)也越小。這是因為庫水位下降速率越快，上下孔壓形成的差值也就越大，從而出現(xiàn)浸潤線“延遲”現(xiàn)象，向下的滲流力加劇了滑坡的不穩(wěn)定，使最小安全系數(shù)更小。③庫水位下降后安全系數(shù)在計算時間內(nèi)最終大于初始安全系數(shù)。這是因為滑坡體內(nèi)部的水位降低，土體的強度增加，從而安全系數(shù)比初始狀態(tài)下更高。④下部淺層的最小安全系數(shù)要小于同工況下的深層滑動最小安全系數(shù)。這說明了在庫水位驟降的情況下，首先發(fā)生的是滑坡體下部淺層滑動，然后才會引發(fā)深層滑動，而滑坡體上部淺層滑動面較為穩(wěn)定，不會發(fā)生失穩(wěn)破壞。

4結(jié)論

（1）庫水位下降速率越大，滑坡體內(nèi)部孔壓下降幅度與下降速率越大，但是滑坡體上部的孔壓下降幅度較下部要小。

（2）滑坡體上部淺層安全系數(shù)隨時間逐漸增大，而下部淺層和深層滑坡安全系數(shù)呈現(xiàn)先減小后增大的趨勢。

（3）總體而言，最小安全系數(shù)從小到大的排序為：下部淺層最小安全系數(shù)、深層滑坡最小安全系數(shù)、上部淺層最小安全系數(shù)。說明滑坡體的下部淺層先發(fā)生失穩(wěn)，然后引起深層整體滑動，而上部淺層較為穩(wěn)定。因此，在實際工程中應該注意著重對滑坡體下部采取相應的治理措施。

參考文獻：

[1]中村浩之，王恭先，論水庫滑坡[J].水土保持通報，1990（1）：53-64.

[2]盧博，郭永成，趙二平，等，庫水位變化和降雨條件下滑坡滲流特性及穩(wěn)定性分析[J]三峽大學學報：自然科學版，2017，39（2）：54-59.

[3]魯芃，王世梅，郭振，等.庫水位日降幅對白家包滑坡穩(wěn)定性影響分析[J].人民長江，2017，48（9）：50-53.

[4]郭子正，殷坤龍，唐揚，等.庫水位下降及降雨作用下麻柳林滑坡穩(wěn)定性評價與預測[J].地質(zhì)科技情報，2017（4）：260-265.

[5]鐘啟明，霍家平，劉若星，庫水位驟降對非飽和壩坡穩(wěn)定性的影響[J].水利水電科技進展，2012，32（6）：84-86.

[6]付宏淵，曾鈴，王桂堯，等.降雨人滲條件下軟巖滑坡穩(wěn)定性分析[J].巖土力學，2012，33（8）：2359-2365.

[7]Fredlund D G，Rahardjo H.Soil mechanics for unsaturat-ed soils[M]. Canada：Wiley， 1993.

[8]劉釗，柴軍瑞，陳興周，等.庫水位驟降時壩體滲流場及壩坡穩(wěn)定性分析[J].西安理工大學學報，201 1，27（4）：466-470.

[9]王冬林，李宗利，張洪泉，庫水位驟降對均質(zhì)土壩壩坡穩(wěn)定的影響分析[J].人民黃河，2011，33 （4）：147-149.