水位漲落引起的庫岸再造對岸坡穩定性的影響研究

楊 廣，晏鄂川，李慶偉

(中國地質大學(武漢)工程學院，湖北 武漢 430074)

三峽工程建成后，上游水位抬升超過百米，庫區淹沒總面積近8萬km2，涉及移民117.15萬人。庫水位升高引起庫岸斜坡水文地質條件發生顯著的改變，產生了許多工程地質問題。據統計，90%以上的斜坡破壞與水的作用有關。已有文獻[1-4]對庫水位升降作用下的邊坡滲流場和穩定性進行了研究，認為水是導致坡體變形破壞的主控因素。在三峽庫區，庫水位大幅波動造成岸坡滲流場改變的同時，其強烈的侵蝕堆積作用也不斷地影響著庫岸再造的強度，宏觀上即表現為地形地貌的改變。孫冬梅等[5]對庫水位下降時某土質岸坡非穩定滲流進行了分析，認為孔隙氣壓力對邊坡的穩定十分不利；連志鵬等[6]分析了在庫水位和降雨聯合作用下的岸坡滲流特征，并在此基礎上探討了岸坡穩定性的變化特征；盧書強等[7]通過分析滑坡的變形監測數據，認為庫水位下降速率影響滑坡的動態變形特征；張旭等[8]基于非飽和土理論建立了流-固耦合模型，并通過有限元法揭示了水分在坡體內運移與土體變形相互影響的空間效應。

針對庫水位對岸坡整體穩定性的影響，目前的研究往往局限于地下水滲流，而忽略了地表水體對岸坡地形的改造作用。因此，本文在現有分析方法的基礎上，以巫山縣某岸坡為例，結合詳細的勘察資料，并考慮庫水位周期性作用對岸坡地形的改變，利用數值模擬方法對該岸坡地下水滲流特征和穩定性進行了研究。

1 岸坡工程地質概況

本文研究的庫岸斜坡位于重慶巫山縣大寧河與長江交匯處。岸坡北臨大寧河，南臨長江，屬于構造侵蝕、剝蝕河谷地貌，地勢東高西低，山脊總體呈北東向展布，如圖1所示。該岸坡地表沿山脊呈馬鞍型，形成了兩個寬緩平臺：第一平臺位于墳院子一帶，高程為230～280 m，為人工建設后的神女廟，兩側岸坡坡角為20°～40°，高程越高坡角越陡；第二平臺為大、小和尚包一帶，高程為170～190 m，兩側岸坡坡角為20°～30°。山脊平臺之間臨水一側各發育紅梁子溝和小狗架子溝，沖溝切割較深，達幾十米，溝長為80～150 m。岸坡最低點為長江、大寧河，原始河床高程為70～80 m，相對高差為100～200 m。

圖1 2017年岸坡地形圖Fig.1 Topographic map of the bank slope in 2017

該岸坡坡體物質主要為第四系松散堆積物和三疊系嘉陵江組第四段灰巖(見圖2)。其中，第四系堆積物厚度大，廣泛分布于山體斜坡、臺地、階地及河漫灘上，成因類型有沖洪積物、殘坡積物、崩積物和滑坡堆積物，成分由粉質黏土夾碎石、碎塊石土及砂卵礫石等組成，表層較松散，深部較密實；基巖未在研究區內出露，為鉆孔揭露，以灰黑色中厚層微—細晶灰巖為主，堅硬性脆，溶蝕孔隙發育，較破碎，有膠結。

圖2 2017年岸坡典型地質剖面圖Fig.2 Typical geological section of the bank slope in 2017

2 岸坡地形對比分析

1997年11月，長江首次截流，三峽水庫庫水位開始升高；到2010年，三峽水庫第一次實現175 m蓄水目標。至今，三峽工程已連續8年成功實現175 m試驗性蓄水目標。自三峽庫區蓄水以來，原先位于江面以上的岸坡巖土體遇水飽和，物理力學強度降低，然后不斷受到周期性的河流地質作用，產生沖刷和淤積。對于土質河岸，河流的持續沖刷使得岸坡變陡，造成塌岸，崩塌后的土體在坡下淤積，產生的淤積物又被水流沖走，開始新的沖刷和崩塌[9]。如Simons等認為河岸沖刷強度不僅與水流的拖曳力有關，還受到風浪和岸坡幾何形態的影響[10]；石盛玉等[11]對比不同時期長江某河段水下地形資料，分析了該段河槽的沖淤演變特征，并對其微地貌進行了高分辨率的觀測；張芳枝等[12]認為河流的沖刷作用加大了堤岸和河床的塑性區范圍，降低了堤岸的安全系數。為了分析三峽庫區開始蓄水至今岸坡地形的變化，本次研究選取了1997年長江三峽庫區巫山縣遷建新址地形圖和2017年實測庫岸勘查平面圖進行了對比，見圖3和圖4。

圖3 岸坡地形平面對比圖Fig.3 Planar comparison graph of the bank slope terrain

圖4 岸坡地形剖面對比圖Fig.4 Comparison graph of the bank slope terrain profile

由圖3和圖4可見，隨著河流地質作用高程不斷抬升，岸坡不同位置和高程處地形改變的方式和程度均有所不同:北側岸坡主要表現為大寧河的淤積作用，南側則表現為長江的側蝕作用，且這種河流的沖淤作用在低水位表現得更加明顯(見圖3)；175 m水位線以上的岸坡，蓄水前后的地形變化不大，而175 m水位線以下的岸坡，地形變化明顯，特別是在水下寬緩平臺的區域，最大淤積厚度達12 m。分析原因認為：大寧河屬典型的山區河流，與長江呈大角度Y型交匯，在北側岸坡附近，大寧河河水流速急劇降低，局部地段形成回水區，產生泥沙淤積，而南側岸坡屬長江干流的凹岸，受江水沖刷嚴重，側蝕現象明顯；兩條沖溝內地形受到地表暫時性水流的影響，淤積作用更突出；由于庫區水位調度的限制，河流沖淤作用在水位變動帶附近表現明顯。現場調查也證實，南側岸坡受到江水的沖刷掏蝕，上覆土體在自重作用下產生崩塌，塌岸現象更嚴重。

3 岸坡地下水滲流特征和穩定性的數值模擬分析

3.1 岸坡數值模型構建

本次研究以整個岸坡為研究對象，考慮地形的變化，并選取庫岸段兩條典型剖面，對蓄水前后不同地形岸坡的地下水滲流特征和穩定性進行數值模擬計算，即通過數值分析軟件GeoStudio的SEEP/W模塊獲取飽和-非飽和條件下岸坡瞬態滲流場，并將結果導入SLOPE/W模塊中進行水-土耦合作用下的岸坡穩定性計算。

本次研究選取1997年和2017年同一岸坡的兩條典型剖面進行計算，以蓄水后2017年實測岸坡A-A′剖面為例(即圖2典型地質剖面)，滲流計算模型見圖5。在模型中定義了3種實體材料，即粉質黏土夾碎石、碎塊石土和灰巖。飽和-非飽和模型的土-水特征曲線和滲透系數隨基質吸力的變化曲線，依土體類型由SEEP/W模塊自帶的模型進行擬合，巖土體物理力學及水力學參數根據土工試驗和工程類比獲得，具體賦值見表1。

圖5 岸坡滲流計算模型Fig.5 Seepage calculation model of the bank slope

為了分析岸坡在庫水位長期作用下，地形改變前后坡體穩定性變化特征，本次計算假設不同時期岸坡均處于同一工況，從而對水庫運行后期岸坡穩定性的變化做出合理預測。因此，岸坡滲流計算取降水位工況，模型邊界條件為：基巖頂面為零流量邊界，左右為給定水頭邊界。其中，給定的水頭值依據三峽工程正常運行時庫水位調度特征(見圖6)賦值，初始庫水位為175 m，庫水位下降速度為0.6 m/d，歷時50 d，并持續145 m低水位運行10 d。岸坡穩定性計算采用M-P法，最危險滑動面通過自動搜索確定。

表1 巖土體物理力學及水力學參數Table 1 Physico-mechanical and hydraulic parameters of the rock and soil

圖6 三峽工程正常運行期庫水位調度圖(資料來源：三峽庫區地質災害防治工程地質勘查技術要求)Fig.6 Reservoir water level scheduling chart during the normal operation of the Three Gorges Project 1.防破壞線；2.限制供水線；Ⅰ.防洪區；Ⅱ.裝機預想出力區；Ⅲ.保證出力區；Ⅳ.降低出力區

3.2 岸坡地下水滲流特征數值模擬計算結果與分析

岸坡在降水位工況下地下水滲流數值模擬計算結果見圖7和圖8。

圖7 岸坡地下水浸潤線及流速矢量分布圖Fig.7 Distribution of the saturation line and velocity vector of groundwater of the bank slope

圖8 不同庫水位下岸坡地下水總水頭云圖Fig.8 Gross head contour of groundwater in bank slope with different water level

由圖7和圖8可見，岸坡地下水滲流具有對稱性。圖7(a)顯示，地下水在岸坡中部位置(相對靠近右側)出現分水嶺，分別從兩側坡腳自由水面排出，左側岸坡地下水流速矢量更為密集;圖7(b)顯示，隨著庫水位的下降，岸坡內地下水平緩下降，水力坡降開始增大，至第50 d達到最大，之后隨著庫水位的穩定運行，地下水浸潤線逐漸趨于平緩，表現出明顯的滯后效應。對比分析圖8(a)至(c)可以看出，岸坡地下水總水頭最小值為當前庫水位，相應位置處總水頭值隨著庫水位的降低而減小，左側岸坡總水頭分布變化比右側岸坡大，即左側岸坡地下水滲流的水力坡降比右側大。

分析產生上述地下水滲流特征的原因主要包括兩個方面：一是岸坡所處的水環境，該岸坡左臨大寧河，右臨長江，且兩側水位受庫區統一調度，因此岸坡地下水具有類似河間地塊的滲流特征；二是兩側岸坡地形差異，左側岸坡地形較陡峭，而右側岸坡地形較平緩，造成左側岸坡地下水滲流水力梯度比右側大，流速矢量更為集中。

3.3 岸坡穩定性數值模擬計算結果與分析

岸坡在降水位工況下穩定性數值模擬計算結果見圖9。

圖9 岸坡在降水位工況下不同剖面穩定性系數變化曲線Fig.9 Curves of stability coefficient of different profiles of the bank slope under the reservoir drawdown condition注：2017、1997代表年份；A、B代表剖面；1、2分別代表大寧河和長江

由圖9可見，總體上看，隨著庫水位的下降，岸坡穩定性系數均有不同程度的降低，并在剛到達最低庫水位時出現最小值，且A-A′剖面岸坡穩定性變化幅度較B-B′剖面更大。對比不同地形條件下相同岸坡的穩定性計算結果，結果表明：2017年岸坡更趨于穩定，即水體對岸坡地形的改造有利于坡體的穩定。分析原因認為：隨著庫水位的下降，地下水產生的向坡外的滲透壓力不斷增大，當庫水位降至145 m時，滲透壓力達最大；隨著庫水位維持在低水位運行，地下水浸潤線開始變緩，水力坡降開始減小，滲透壓力降低，岸坡穩定性逐漸升高。

由上可見，自三峽庫區開始蓄水以來，岸坡受庫水的沖刷淤積，當局部或整體失穩時，發生庫岸再造，產生新的地形仍不斷受到河流地質作用，最終形成使坡體自穩能力更強、趨于更穩定的狀態。

4 結 論

本文以巫山縣某岸坡為例，對比分析了岸坡自庫區開始蓄水以來的地形變化，然后利用數值模擬方法研究了地下水滲流特征，并對不同庫水位下同一岸坡的穩定性進行了數值模擬計算，得出如下結論：

(1) 岸坡不同位置和高程處地形改變的方式和程度有所不同，主要表現為北側大寧河的淤積作用和南側長江的側蝕作用，且這種河流的沖淤作用在低水位表現得更加明顯。175 m水位線以上的岸坡，蓄水前后的地形變化不大；而175 m水位線以下的岸坡，蓄水前后的地形變化明顯。

(2) 對于兩側涉水岸坡，其地下水滲流具有對稱性。分水嶺大約出現在岸坡中部位置，地下水從兩側坡腳自由水面排出;隨著庫水位下降，岸坡地下水總水頭最小值同步減小，地下水浸潤線平緩下降，并表現出明顯的滯后效應；受地形影響，坡度較陡一側岸坡地下水滲流水力梯度更大。

(3) 隨著庫水位的下降，岸坡各剖面穩定性系數降低，在低水位時達到最小值。同一工況下，2017年實測岸坡穩定性較1997年有所提高。

通訊作者：晏鄂川(1969—)，男，博士，教授，主要從事巖土體穩定性評價與利用方面的研究，E-mail：yecyec6970@163.com