田壩村變形體的變形機理及穩定性分析

張小寶，李 坤

(水利部水利水電規劃設計總院，北京 100120)

1 工程概況

天花板水電站壩址位于云南省魯甸縣翠屏鄉與巧家縣鉛廠鄉境內，廠址位于魯甸縣翠屏鄉境內。壩型為碾壓混凝土雙曲拱壩，壩高107m，水庫河道狹窄，總庫容為7871萬m3。

田壩村斜坡位于天花板水電站庫區左岸岸坡上，與已建成的大壩最近距離為1.0km左右，考慮到坡體目前變形較大，對已建成的天花板水電站構成潛在威脅(如圖1所示)。

圖1 田壩村斜坡研究區與天花板水電站空間相對位置

2 田壩村變形岸坡基本地質條件及變形機理分析

2.1 岸坡基本地質條件

2.1.1地形地貌

田壩村變形體位于天花板水電站庫區左岸，近壩庫段1km處，為典型的古堆積體岸坡，沿河向寬度約為1.3km，變形體前緣位于河床，高程為980m，后緣拉裂縫最大高程為1662m。邊坡上分布有灣子社、窩當、楊家天井、歐家寨、炭山、黃家坪子等村落，有200余戶居民。邊坡自然坡度為20°～50°，邊坡上多為梯田，蓄水前邊坡穩定性較好。蓄水后受前緣塌岸、降雨及水位驟降影響，邊坡在2011年1月開始發生變形，經勘察粗略估算，古堆積體總方量約為8000萬m3，發生變形的方量約為3500萬m3，屬于超大型的滑坡變形體。

2.1.2地層巖性

根據鉆探、坑槽探資料，邊坡堆積物厚度較大，中部和前緣的堆積物厚度均在40m以上，前緣最大堆積厚度達到了150m，后緣厚度較薄，一般為10～20m。邊坡堆積體由崩積、崩坡積及坡殘積物組成，主要為土夾碎(塊)石，局部塊石含量較高，塊石直徑最大可達20余米。下伏基巖為寒武系和震旦系的白云巖、白云質灰巖及砂巖等。

2.1.3水文地質特征

在邊坡發生變形前地下水位較高，地表溝谷內有小溪等徑流分布，邊坡上的居民生活用水主要由邊坡上的地表徑流供給。

邊坡發生變形后，邊坡上的縱向、橫向裂縫發育，地表水下滲，地表徑流消失。根據鉆孔水位監測資料成果，變形體的地下水位主要分布在基巖與覆蓋層附近，大氣降水通過地表裂縫下滲，以暗流的形式排到牛欄江。

2.1.4變形體物理力學參數

本項研究主要從現場地質調查入手，通過試驗成果以及切割小剖面進行統計的手段來總結碎石土內部礫石粒徑分布的主要特點，根據現場情況及崩塌堆積體的成因和分布特征，提出物理力學參數的建議值，見表1。

2.2 變形岸坡分區

蓄水前岸坡基本穩定，2011年1月17日蓄水后，經當地政府有關部門反映，老店鎮灣子社、坪子社及陡地社部分村民住房地基及房屋墻體出現開裂變形。根據收集到的田壩村變形體地形地貌特征、裂縫分布形態、鉆孔資料等，并在分析變形體成因和變形特征的基礎上，將田壩村變形體大致劃分為三個區，Ⅰ區崩塌堆積變形體、Ⅱ區滑坡體和Ⅲ區崩塌堆積體，如圖2所示。

圖2 變形體分區圖

2.3 變形體的成因機理分析

田壩村變形體基巖與土石混合體接觸面從整體上控制著變形體的變形軌跡，但由于變形體前緣瀕臨牛欄江，隨著江水的逐漸淘蝕，塌岸現象日益嚴重，前緣阻力逐漸減小，從而促使臨江岸坡局部失穩滑動，當岸坡失穩滑動后，進一步促使后部坡體產生下滑趨勢，產生拉裂縫，導致連鎖反應，從這個角度看，該變形體具有牽引式的特點。

具體來說，田壩村變形體演化過程如下：

(1)崩塌堆積物的形成

崩積的間歇發生—形成崩積堆積物—物理化學風化—形成土石混合體—壓實—形成基巖上覆表層堆積體。

(2)變形體的形成

水位驟降、大氣降水及風浪影響—邊坡前緣發生塌岸—前緣的局部塊體變形失穩—牽引后緣堆積物發生變形—坡內裂縫及邊界裂縫形成—降雨及地表水入滲—加劇邊坡變形—變形體形成。

3 變形體穩定性定量評價

本章節擬借助Geo-Studio數值軟件在研究區選取5個典型計算剖面(如圖3所示)進行定量計算，以便能夠對其穩定性現狀和發展趨勢進行預估評價。

3.1 崩坡積碎石土強度參數反演

3.1.1參數敏感性分析

首先進行強度參數敏感性分析，了解黏聚力和內摩擦角二者對安全系數的影響程度，便于實施后續計算參數調整過程。

從圖4中可以看出，黏聚力和內摩擦角二者對安全系數的影響基本成線性，隨著黏聚力和摩擦角的增大，安全系數基本成線性增大。同時可以直觀地看出：黏聚力對邊坡安全系數不是很敏感，而內摩擦角異常敏感。

3.1.2強度參數反演

在參數敏感性分析基礎上，此次強度參數反演時內摩擦角取值范圍為20°～40°，黏聚力取值范圍為0～10kPa?；趨蹈鞣N組合，計算圖4中潛在滑動面的安全系數，結果見表2。顯然，若按照表1所示的當地參數經驗取值，當黏聚力按照0kPa考慮時，內摩擦角不論在26°～29°范圍內如何變化，其潛在滑動面安全系數在0.87～0.99變化，且均小于1，因此上述經驗參數取值可以理解為崩坡積碎石土的下限。考慮到該崩坡積碎石土仍具有一定的黏結性，從而將其黏聚力確定為5～8kPa，那么為了使潛在滑動面安全系數達到1.00，摩擦角至少為29。

表1 變形體土體物理力學參數建議值表

圖3 5個典型地質剖面空間位置(單位：m)

基于上述分析，結合地形特征，此次邊坡穩定性計算時，崩坡積碎石土的黏聚力取值范圍為5～8kPa，內摩擦角取值范圍為29°～32°。

3.1.3數值計算材料參數

基于上述強度參數反演取值，并結合經驗，此次數值計算相關材料參數見表3。

圖4 黏聚力、內摩擦角敏感性分析

表2 潛在滑動面的安全系數

表3 數值計算材料參數

3.2 數值計算模型

基于地質剖面PM1、PM2、PM3、PM4、PM5建立水文地質概化模型，進而建立相應的數值計算模型，如圖5所示。

3.3 失穩模式預測

基于上述所建立的模型，通過剛體極限平衡法依次對研究區Ⅰ-1區、Ⅰ-2區、Ⅰ-3區、Ⅱ區、Ⅲ區坡體在現狀自然條件下、降雨工況、庫水位變化以及地震作用下可能發生的失穩模式進行了預測，結果見表4。

表4 各剖面計算結果統計表

圖5 各剖面計算模型

3.4 穩定性評價

表5為DL/T 5353—2006《水電水利工程邊坡設計規范》[5]中水利水電工程邊坡設計安全系數，此次研究區坡體按照水庫邊坡B類考慮。

另外，根據DL 5180—2003《水電樞紐工程等級劃分及設計安全標準》[6]的規定，天花板水電站為三等中型工程，樞紐主要建筑物中碾壓混凝土雙曲拱壩、泄水建筑物、進水口、引水隧洞、發電廠房等主要建筑物為3級。

綜上分析，此次邊坡穩定性評價時，地震工況按照偶然狀況考慮，取1.00作為考慮地震工況時邊體穩定性判別依據；降雨、水位下降工況按照短暫狀況考慮，取1.05作為考慮降雨或者水位驟降工況時邊體穩定性判別依據；自然工況按照持久狀況考慮，取1.10作為自然工況時邊體穩定性判別依據。

4 結論

Ⅰ區PM3在現狀自然條件、降雨、水位下降工況下，其安全系數均接近規范規定值，PM4在各工況下，其安全系數均低于規范值，穩定性差。

Ⅱ區在非地震工況下，整體滑動的安全系數接近規范值；僅在現狀自然條件、降雨、水位下降工況下，前緣局部的安全系數低于規范值，穩定性差。

Ⅲ區在現狀自然條件、降雨、水位下降工況下，其安全系數均高于規范規定值，因此認為上述工況下Ⅲ區處于較穩定狀態；而在地震作用下，其穩定性差。

變形體前緣牛欄江河道狹窄，根據測量資料，蓄水后由于岸坡塌岸，變形體前緣河床已經比蓄水前抬高了將近20m，間接起到了壓腳作用，綜合分析認為田壩村變形體產生整體快速失穩的可能性不大，隨著前緣塌岸進一步變形，前緣堆積物厚度加大，有逐漸趨穩的可能，建議進一步加強變形監測，控制水位驟降工況下的穩定性。

2014年8月3日，云南省昭通市魯甸縣發生6.5級地震，根據監測資料，變形體的變形發生一次較大突變后趨穩，根據后期監測資料，再無變形發生。

天花板水電站田壩村變形體是水庫蓄水后，由于塌岸、暴雨及水位驟降引起的近壩庫岸土質岸坡失穩的典型案例。變形發生后，居民房屋開裂、稻田干涸、果樹傾倒及溪水斷流，直接影響了6個自然村居民的生產生活，造成了較大的財產損失，是一次深刻的教訓，這對今后水庫土質岸坡勘察、古堆積體的穩定性研究具有較強的借鑒意義。

表5 水利水電工程邊坡設計安全系數