邊坡群穩定性耦聯分析方法在三峽庫區的應用

華小軍,樊啟祥,尚毅梓,程 建,張玉柱,劉志武

(1.中國長江電力股份有限公司,湖北 宜昌 443002; 2.中國長江三峽集團公司,北京 100038)

我國建成運行的高壩大庫多位于岸坡陡峻、構造復雜、巖體破碎、軟弱地層發育、崩滑災害體廣布等地質環境極脆弱的地帶。由于防洪、發電等要求,水庫水位變動頻繁,改變了庫區邊坡的水文地質環境,加劇了滑坡等地質災害風險[1-3],反過來也制約了水庫綜合效益的發揮。以三峽為例,根據2009年7月湖北省、重慶市庫區26個區縣政府上報統計資料及國土資源部門調查,崩塌、滑坡共5 386處,其中涉水的崩塌、滑坡1 945處,不涉水的崩塌、滑坡3 441處[4]。因此,研究三峽庫區邊坡穩定性具有重要意義。

目前,對水庫邊坡穩定性問題的研究已經較多。胡鐵松等[5]提出了滑坡預測的人工神經網絡方法,并將該方法應用于清江隔河巖庫區兩個主要滑坡的預測,預測結果表明了該方法的有效性以及在滑坡預測中的良好應用前景;賈官偉等[6]利用圖像動態采集系統記錄了水位驟降引致滑坡的發生、發展過程和滑坡模式,為探討邊坡失穩原因和失穩模式提供了依據;蔣秀玲等[7]采用Morgenstern-Price法對滑坡穩定性進行計算,對三峽水庫水位變動下的庫岸滑坡穩定性進行了評價;向杰等[8]研究了三峽水庫蓄水誘發神女峰庫岸破壞機制,并簡述了在庫水位變化情況下,岸坡穩定性分析的計算方法;韓凌風等[9]研究了庫水位升降作用下庫岸邊坡滲流特征,分析了該庫岸邊坡的孔隙水壓力變化特征和滲流速度動態變化規律,為滲流與變形耦合分析提供了依據;鄧華鋒等[10]提出了對庫岸邊坡分段的分析方法,從力學機理上較好的解釋了庫岸邊坡在水位上升或下降過程中安全系數先減小后增大的原因;段永祥等[11]進行了庫岸滑坡穩定性物理模擬研究,試驗結果表明水庫水位升降對滑坡前緣涉水部分影響較大,對中部和后緣影響較小,降雨對滑坡前緣、中部、后緣影響相對顯著;張嘎[12]提出了巖土域的概念,認為相鄰邊坡體的特性和變化規律具有相關性和相似性,在此基礎上構建了邊坡群穩定性耦聯分析概念性模型。

以上研究大多針對整個邊坡的定性評價或特定邊坡的定量分析,缺乏對邊坡整體安全性的定量分析,且對庫區邊坡群地質資料、水文信息等數據缺乏時的研究較少。本文基于張嘎建立的邊坡群穩定性耦聯分析概念性模型,提出一種庫區邊坡群穩定性分析方法。以三峽庫區一個典型邊坡群為實例,基于“實時反饋響應、同步修正模型”的思路,對該方法的有效性進行驗證和分析。

1 三峽庫區滑坡體概況

所研究的邊坡為三峽庫區某滑坡體(圖1),該滑坡體三維幾何計算模型由滑坡體、軟弱帶和基巖3部分組成,在平面上取包圍滑坡體的大約1 240 m×1 000 m的區域。有研究表明[13],從2004年1月15日起該滑坡體開始蠕動變形,進入滑坡發展階段,在2007年7月期間,局部地區的滑坡形變達到最大30 cm。在2009年5月至8月期間由于降雨影響和庫水位下降，該滑坡體發生了規模較大的形變,一些觀測點的數據表明位移量在50～550 mm之間,滑坡體往北方向滑移,局部存在偏西或者偏東方向的調整,滑坡引起的裂縫在滑坡體上隨處可見。

圖1 三峽庫區某滑坡體30 m分辨率DEM圖

前期沒有對該滑坡體進行勘察,信息缺乏,地質資料、水文信息等相關計算參數都沒有提取到。這些資料需要聯系當地國土資源局獲取,有一定的難度,而三峽庫區滑坡體資料本身也屬于敏感信息,故根據衛星遙感信息圖,找到一幅滑坡前的數字高程模型。運用ArcGis遙感信息技術軟件對DEM原始數據進行處理,提取滑坡體的等高線和剖面信息,可以獲取滑坡體的形態信息見圖2。

圖2 三峽庫區某滑坡體地形地貌

2 庫區邊坡群穩定性耦聯分析方法

2.1 巖土域相似性原理

庫區邊坡群具有相似的地質條件和外部環境,雖然每個邊坡體的幾何特征、土層構成及特性等要素有所不同,但也存在著明顯的共同屬性。如各邊坡體的地質構造與形成過程相似,因此其力學特性、滲透特性具有相關性，應力歷史具有相似性;各邊坡體的邊界條件具有宏觀一致性,其外水位、外荷載(如地震等)等具有相等或者相似的性質;各邊坡體的變形破壞規律也有相似性。這些相關性使得庫區內一個邊坡體的數據信息對其他邊坡體具有參考意義,因此將處于同一地域、具有內在聯系的多個巖土體(如邊坡等)構成的系統稱為巖土域。

2.2 耦聯分析概念性模型

巖土域耦聯分析概念性模型結構見圖3。該模型包括輸入層、模型層、計算層和響應層4個層。其中,輸入層用于收集多源數據和確定邊界條件和參數等,并傳遞到模型層;模型層用于建立單個巖土體的分析模型和相互間的耦聯關系;計算層用于多源數據的挖掘和分析,以及巖土域和所含巖土體分析模型的求解;響應層用于匯總各巖土體的響應及其耦聯特性,并與輸入層進行交互。

圖3 巖土域耦聯分析概念性模型

巖土域耦聯分析概念性模型具有如下主要特點:①模型以數據驅動為主,一旦有新的數據進入響應層或者輸入層,分析過程就將啟動,而傳統方法是過程驅動;②模型具有開放性,各巖土體的模型和響應相互獨立;③模型結果是定量的,而不是概率和可靠度分析;④數據和參數是多源的和靈活的。

2.3 庫區邊坡群穩定性耦聯分析步驟

a. 輸入層。輸入滑坡體初始滑坡剖面的參數和水位。容重和強度參數采用取樣試驗與工程類比后的參數,水位根據水庫蓄水情況變化。

b. 模型層。采用簡化Bishop條分法建立邊坡穩定性的分析模型,計算公式見式(1)(2),將滑坡體初始滑坡剖面實際滑坡形態導入耦聯分析模型。

(1)

(2)

式中:Fs為邊坡的安全系數;下標i為土條序號;b為土條的寬度;W和α分別為土條的重量和傾角;c和φ分別為土的黏聚力和內摩擦角;Mc(γ)和Mφ(γ)分別為黏聚力折減系數和摩擦力折減系數;γ為土條的浮容重;mαi為抗滑力矩。

c. 計算層。運用數據挖掘、極限平衡法等分析方法進行迭代計算,反算得出滑坡體達到極限平衡時的參數。

d. 響應層 運用上一步得出的參數對滑坡面進行預測推算,如果得出的滑坡形態與實際滑坡形態不符合則自動調整參數重新迭代計算,直至符合為止,此時的參數作為最終的計算參數,運用這些參數對滑坡體其他可能滑坡剖面進行安全分析,得出每個剖面的安全系數。

3 計算結果與分析

根據工程經驗和現場蠕動變形痕跡,對滑坡體的各個可能滑坡方向進行標記,劃分9個斷面進行計算分析,如圖4所示。由于滑坡體的初始滑裂面和變形體位置在剖面8處,故選取剖面8為計算基準剖面。

圖4 三峽某滑坡體選取斷面標記

將剖面8實際滑坡形態導入耦聯分析模型,模型運用極限平衡分析方法進行迭代計算,反算得出滑坡體達到極限平衡時的參數。運用此參數對滑坡面進行預測推算,如果得出的滑坡形態與實際滑坡形態不符合則自動調整參數重新迭代計算,直至符合為止,得出主要計算參數:土層容重為10 kN/m3,黏聚強度為28 kPa,內摩擦角為25°。剖面8的安全系數為1.126,計算出的剖面8滑坡形態見圖5。

圖5 耦聯分析模型模擬的剖面8滑坡形態

將圖5與Wang等[14]的實測結果進行對比可知運用耦聯分析模型得出的滑坡形態與該滑坡體的實際滑坡形態基本一致,可得出最終的計算參數。運用最終的計算參數對其他可能滑坡剖面進行安全分析計算,剖面1～7和剖面9的安全系數計算結果分別為0.828、1.276、0.906、2.204、1.144、1.252、1.194、1.064,耦聯分析模型計算出的剖面1滑坡形態見圖6,其他類似。

圖6 耦聯分析模型模擬的剖面1滑坡形態

計算結果顯示,最危險的剖面為剖面1、剖面3、剖面5、剖面9,其中剖面1、剖面3為欠穩定狀態(本次計算將安全系數小于1.15視為潛在危險滑裂面),預測的主要滑坡趨勢集中在北北方向,如圖7所示。實際滑坡位移矢量圖見圖8。整體上看,滑坡體作近北方向滑移,局部地區移動方向有偏西或者偏東的調整,平均位移在30 cm左右。滑坡體的頭部區域是位移量相對較大的區域,箭頭方向顯示這個區域從東南向西北方向順山坡從高處向低處滑移,最大位移接近70 cm。坡體的運動方向在靠近長江的舌部順著地勢發生改變,西邊沿北北西滑動,東邊沿北北東滑動。這個形變場特征和該滑坡體的地形特征吻合甚好,位移大小和方向與李小凡等[13]對滑坡體的野外觀測結果基本吻合,預測的主要滑坡趨勢軌跡方向與實際趨勢形態吻合良好。

圖7 預測滑坡趨勢標記

圖8 實際滑坡位移矢量圖

4 結 語

基于提出的巖土域概念,將巖土體分析擴展到巖土域分析,通過利用其他巖土體數據,采用耦聯分析、大數據挖掘等手段提高了邊坡穩定的分析精度。

提出了一種庫區邊坡群穩定性分析方法,以三峽庫區一個典型邊坡群為實例,在邊坡群資料不全、信息缺乏的情況下,利用遙感信息技術提取相關特性,計算了該邊坡的穩定性。結果表明預測的滑坡體位移大小和方向與野外觀測結果基本吻合,預測的主要滑坡趨勢軌跡方向與實際趨勢形態吻合良好,驗證了本文方法的有效性。

本文只對水位固定狀態下邊坡的穩定性進行了定量分析,后期可以將此耦聯模型應用于水位變化條件下的邊坡穩定性耦聯過程分析。