去學水電站色叢組庫岸邊坡整體穩定性分析

于立宏，鄭亞臣，范文進，韓立軍，張懷芝

(中國電建集團北京勘測設計研究院有限公司，北京 100024)

1 工程概況

去學水電站位于四川省甘孜藏族自治州得榮縣和云南迪慶藏族自治州香格里拉縣交界部位的金沙江二級支流碩曲河干流上。水庫正常蓄水位2330m，死水位2310m，總庫容1.326億m3，瀝青混凝土心墻壩，最大壩高164.2m，裝機246MW。地震基本烈度為Ⅷ度。

色叢組位于一第四系更新統沖洪積階地上，階地對岸為高聳近直立高度超過800m陡崖，階地庫岸邊坡高20～25m，地形坡度65°～75°，基巖巖性為板巖，覆蓋層厚20～30m，成份主要為混合土卵石，弱膠結，密實狀態，弱透水性，混合土卵石與基巖接觸部位為一厚10～20cm的粉土質礫層，粉土質礫層為該區軟弱夾層，如圖1所示。覆蓋層有沿粉土質礫層整體下滑趨勢[1]。水庫蓄水后，地下水位會被壅高，必然導致該區粉土質礫層部分位于地下水位以下，由于水文環境的改變，從而導致粉土質礫層性狀發生變化[2]，覆蓋層整體是否會沿著粉土質礫層發生滑動破壞，對于去學水電站水庫蓄水方案確定、庫岸防護方案選擇、水庫是否正常運營等意義重大。

圖1 色叢組典型工程地質剖面圖

2 色叢組庫岸邊坡整體穩定性分析

2.1 工程地質分析

該區為階地，所以色叢組庫岸邊坡形成時間要早于現代河床內覆蓋層形成的時間[3]。從圖1中可以看出，階地區基巖面已高于河流水位高程，高于河床內河底基巖面43m，說明碩曲河在下切階地下面基巖的漫長過程中，色叢組庫岸邊坡一直處于沿粉土質礫層下滑的狀態，在這段漫長的地質歷史時間內，邊坡地表不斷受到水流沖刷、風化剝蝕作用[4]，相當于邊坡不斷處于“卸載”過程中，最終形成了如今的邊坡地貌形態、地形剖面，邊坡在地質歷史時期中覆蓋層不斷經受著“成巖”作用[5]，邊坡坡體自身整體性在不斷加強。

中國地震局分析預報中心地震數據庫查明：自1515年有地震記錄以來，邊坡區可能遭受過17次Ⅴ度烈度和4次Ⅵ度烈度地震影響，其中1970年1月—2016年10月邊坡近場區就記錄到M2.0～2.9級地震46次，M3.0～3.9級地震31次，M4.0～4.9級地震12次。當地老鄉(尤其是長者)也未見到過該區覆蓋層邊坡整體發生過滑動破壞。具體代表性強震對該工程區影響程度見表1。

由上述分析及表1中成果可以看出，該區覆蓋層邊坡在漫長的地質歷史中經受住了地震烈度為Ⅵ度地震作用而未發生整體破壞，說明該區可能經受了地震動峰值加速度值至少為0.04g的地震而保持整體穩定[6]。如果此種情況下，該區碩曲河河水位高于去學水電站水庫正常蓄水位情況，那么可以確定去學水電站蓄水后，該區邊坡依然穩定；如果不是此種情況時，可以假定地震動峰值加速度值為0.04g地震發生時，邊坡覆蓋層整體處于不利受力狀態情況下的極限平衡狀態，即可以假定邊坡覆蓋層整體在地震動峰值加速度值為0.04g地震發生時穩定性安全系數為1，這樣就可以在該區選擇2個典型單寬工程地質剖面利用極限平衡公式組成方程組進行反演計算[7]，從而求出邊坡抗滑動破壞時相對保守的滑動面抗剪參數C值、φ值，受力分析原理如圖2所示。

表1 代表性強震對色叢組工程區影響程度

圖2 地震發生時物體受力原理

由圖2分析推導出的邊坡單位寬度極限平衡條件下穩定性計算公式為：

(1)

式中，F—安全系數；G—下滑物體重力；G′—水平地震力，G′=0.04G；α—邊坡坡角；l—滑面長度；c，φ—粉土質礫層抗剪強度參數值。

2個典型單寬剖面的具體情況見表2。

表2 2個典型單寬剖面的具體情況

剖面1-1計算過程及成果如下：

=[(3.103×107×cos 20°-0.04×3.103×107×sin 20°)·tgφ+130.8×106·c]/(3.103×107×sin 20°+0.04×3.103×107×cos 20°)

令F=1，則：

=2.439tgφ+11.104C

(2)

剖面2-2計算過程及成果如下：

=[(2.528×107×cos 20°-0.04×2.528×107×sin 20°)·tgφ+122.2×106·c]/(2.528×107×sin 20°+0.04×2.528×107×cos 20°)

令F=1，則：

=2.439tgφ+12.734C

(3)

方程(2)、(3)聯立成方程組，解得C=0MPa，tgφ=0.41。

2.2 計算分析

水庫蓄水后，因為混合土卵石透水性相對較強，庫水和地下水連通性相對較好，地下水位被壅高，一部分粉土質礫層將會位于地下水位以下，覆蓋層會受到地下水產生的力作用，受地下水影響粉土質礫層抗剪強度會降低。對地下水位以下粉土質礫層抗剪強度參數進行折減[8]，根據葛洲壩工程泥化夾層的試驗研究資料成果，將地下水位以下的粉土質礫層抗剪強度參數數值取2.1節算出值的0.8倍。

利用圖2分析推導出的計算公式為：

(4)

式中，F—安全系數；G1—地下水位以上下滑物體重力；G2—地下水位以下考慮地下水影響的下滑物體有效重力；α—邊坡坡角；l1—地下水位以上滑面長度；l2—地下水位以下滑面長度；C1、φ1—地下水位以上粉土質礫抗剪強度值，C1=0.0MPa，tgφ1=0.41；C2、φ2—地下水位以下粉土質礫抗剪強度值，C2=0.0MPa,tgφ2=0.41×0.8=0.328。

水庫正常蓄水位運行時，該區庫水位高程為2330(水庫正常蓄水位)+0.17(尾水超高值)=2330.17m，2個典型工程地質單寬剖面的具體情況見表3。

剖面1-1計算過程及成果如下：

=(2.04×107×cos20°×0.41+0.46×107×cos 20°×0.328)/[(2.04+0.46)×107×sin 20°]=1.08

表3 2個典型工程地質單寬剖面的具體情況

剖面2-2計算過程及成果如下：

=(0.58×107×cos 20°×0.41+0.50×107×cos 20°×0.328)/[(0.58+0.50)×107×sin 20°]=1.02

計算成果表明2個典型工程地質剖面覆蓋層安全系數均大于1，說明水庫正常蓄水位運行工況下，2個典型的工程地質剖面上覆蓋層整體是安全的，覆蓋層邊坡整體是穩定的。

2.3 綜合分析

大壩建成水庫下閘蓄水后，粉土質礫層抗剪參數軟化折減至水庫蓄水后計算所用到的參數值是一個漫長過程[9]，同時，隨著水庫的形成與運行，碩曲河河水變成庫水，水庫淤積的進行，時間的推移，淤積物會越來越多[10]，會致邊坡坡角覆蓋的淤積物沉降、密實，增加了邊坡整體穩定性。同時，由于水庫蓄水，水位高程抬高，邊坡一部分將位于水下，地下水和庫水作用令混合土卵石層局部穩定性受到影響而降低[11]，邊坡局部將發生變形以至破壞，邊坡局部穩定性差的地方由于環境改變將發生坍方，將產生庫岸再造現象[12]，塌坍下來的物質將堆積在坡角，這些物質增加了邊坡整體穩定性。這些因素共同作用下，該邊坡在水庫蓄水后，明顯遭受到“削坡減載”和“壓腳”作用[13]，有利于庫岸邊坡整體穩定。

綜上分析，該區庫岸邊坡在水庫蓄水后整體穩定性依然較好。

3 結語

如今，去學水電站水庫蓄水已二年，色叢組庫岸邊坡僅局部發生了崩塌破壞，邊坡整體沒有發生滑動破壞，也未見有發生整體滑動破壞的跡象，邊坡整體穩定性增強了，這和分析結論相一致，說明邊坡整體穩定性分析是正確的。

無論工程地質條件多么復雜，在信息和科技發達的今天，只要我們加大理論和實踐相結合的力度，認真研究分析，總會找到解決問題的途徑和方法，色叢組庫岸邊坡整體穩定性分析是在工程地質分析基礎上，充分利用地震的科研成果進行的，效果較好，為大家提供了一個解決問題的思路。