大華橋水電站拉古滑坡穩定性分析研究

黃志義，張永輝，孫凱輝

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

大華橋水電站是瀾滄江上游河段一庫七級開發方案的第六級電站，壩型為碾壓混凝土重力壩，最大壩高106 m，庫容25 600萬m3[1]。目前，大華橋水電站已經建成并開始蓄水運行，庫區范圍內沿線地形陡峭，覆蓋層分布厚度不均，沿線主要發育3處大型滑坡體，對水庫安全運行和周邊居民生產生活造成重大潛在隱患，是影響大華橋水電站后期運行的重大危險因素[2]。根據現場地質調繪，3處滑坡體均存在不同程度的變形跡象，其中，拉古滑坡（見圖1）前緣局部岸坡已經失穩滑動，坡體分布大型深切沖溝5條，坡體上游中部分布有拉古中村、清口寨等居住點，對拉古滑坡的穩定性進行研究顯得十分必要[3～5]。針對這一問題，本文依據前期勘察測繪成果，對拉古滑坡的地形地質和水位地質條件進行分析的基礎上，利用現場監測手段對拉古滑坡的潛在滑動狀態進行了歷時1年的監測，對拉古滑坡在不同蓄水條件、不同工況下的安穩定性問題進行分析計算。

圖1 拉古滑坡全貌

1 滑坡概況

拉古滑坡位于大華橋水電站大壩壩址上游約12 km處，滑坡所在山體高程為1400～2000 m，滑坡平面呈舌形，前緣較寬，后緣窄小，整體滑動方向為近南北向，沿瀾滄江沿岸分布寬度為1600 m，受地形切割影響滑坡表面呈現典型的階梯斜坡地貌[6]。根據前期勘察測繪，滑坡堆積體方量約為5800萬m3，規模巨大，同時滑坡體上部存在居民村落和農田。在水庫蓄水至正常蓄水位后，滑坡體坡腳被淹沒深度約47 m[7]。根據滑坡體的滑動形態和物質構成，將拉古滑坡劃分為3個主要區，其中Ⅲ區位于滑坡后緣，Ⅰ區位于滑坡前緣上游，Ⅱ區位于滑坡前緣下游，現狀條件下Ⅰ區為典型的滑坡區，Ⅱ區、Ⅲ區目前相對穩定。

根據現場地質勘探成果，研究區上部覆蓋層為崩坡積碎塊石土，土體結構松散～稍密，主要由粉土及次棱角碎塊石組成，塊石直徑為2～15 cm，覆蓋層厚度一般為15～30 m；下伏基巖為全強風化紫紅色板巖夾少量砂巖，巖體結構較完整，其中板巖占比約為95%，巖層產狀一般為NE10°～15°NW∠65°～85°，均傾向岸內。根據鉆孔巖芯揭露，拉古滑坡未見明顯的貫通滑動帶，但局部發育次級滑動帶，滑帶土以礫石、碎屑及粘粒構成，礫石一般粒徑3～8 mm，呈次棱角狀，半固結狀態。

根據滑坡場地內長觀孔觀測結果，地下水位埋深一般為20～30 m，最大埋深60 m，地下水位主要以基巖裂隙水和覆蓋層潛水為主，受大氣降水補給和臨近山谷河流側向補給。長期水位監測結果發現，拉古滑坡體內實測地下水位高程在1 469.89～1 847.25 m之間，水位變動幅度為0.4～4.1 m，水位變動帶基本位于上部覆蓋層內。

2 滑坡體監測

2.1 監測方案布置

從前期勘測成果來看，滑坡體Ⅰ區滑動跡象最為明顯，Ⅱ區和Ⅲ區現狀條件下較為穩定，本文重點分析Ⅰ區滑坡體的變形監測過程。在Ⅰ區1-1剖面布設1套陣列式位移計和3處測斜儀（見表1）。

表1 拉古滑坡Ⅰ區監測儀器統計表

2.2 監測結果

依據監測數據，繪制深度位移變化曲線（見圖2～圖4）。

圖2 LG-IN1-1測斜孔孔深累計位移曲線

從圖2～圖4可見，深度位移隨著深度增加而減少，隨著時間增加而增大。在測斜孔1-1中存在一處滑移帶（距孔口11～13 m附近），在測斜孔1-2中存在一處滑移帶（距孔口24～26 m附近），在陣列式位移計孔1-1中存在一處滑移帶（距孔口42～45 m附近）。

圖3 LG-IN1-2測斜孔孔深累計位移曲線

圖4 LG-SAA1-1陣列式位移計累計位移曲線

2.3 滑坡體位移等值線分析

為充分研究拉古滑坡潛在滑移面，將監測斷面位移隨深度變化過程進行等值線成圖（見圖5），發現該滑坡體存在3處潛在滑移帶，其三維等值線云圖見圖6。

圖5 拉古滑坡Ⅰ區1-1剖面位移等值線圖

從拉古滑坡三維等值線云圖（圖6）來看，滑坡體深部滑坡主要發育在Ⅰ區中下部和Ⅱ區中部，滑移帶不連續，初步分析認為其主要滑移方式為局部次級滑動推動，其可以作為滑坡穩定性分析和評價的前提。

圖6 拉古滑坡三維等值線云圖

2.4 監測成果初步分析

Ⅰ區1-1剖面高程1480 m以上天然坡度為20°～25°；高程1480 m以下天然地形坡度為25°～30°，下伏基巖面坡度與地形基本一致，地下水位基本位于覆蓋層內。該斷面布置了4個監測點，從低到高分布大約高程為1534、1588、1692、1880 m。通過監測數據表明：該斷面上部變形以深層滑移為主，測斜孔LG-IN1-1在孔深12 m處最大位移為49.34 mm、LG-IN1-2在孔深25 m處最大位移為37.94 mm，可能在該處存在一滑動帶，屬于覆蓋層內部的滑移；該斷面下部變形為深層滑移，陣列式位移計LG-SAA1-1在孔深43 m處累計位移分別為87（X方向）、512（Y方向）、279 mm（Z方向），可能在基覆界線附近存在一滑動帶。該斷面主要位于Ⅰ區，距離上游側沖溝較近，沖溝切割較深，對滑坡的變形位移影響較大。

3 滑坡穩定分析

3.1 邊坡建模

3.1.1 模型選擇及巖體參數

依據拉古滑坡1-1剖面監測成果，中上部的潛在滑裂面基本處于基覆界線處，而中下部的潛在滑面位于滑坡體中部。穩定分析中復核的滑移模式的滑裂面位置見圖7。

圖7 拉古滑坡Ⅰ區潛在滑動帶示意圖

結合地質勘察的建議值，采用反演所得天然狀態的抗剪強度參數，根據水庫蓄水前岸坡的穩定情況，按持久、短暫、偶然等工況，計算各典型斷面蓄水前的安全系數，驗證反演參數的合理性，在此基礎上結合地質勘察成果確定計算所需分析參數（見表2）。

表2 滑坡體穩定性分析參數

3.1.2 邊坡設計

邊坡設計按下列兩類作用組合：①基本組合：自重+岸邊外水壓力+地下水壓力+加固力；②偶然組合：基本組合+地震作用。

邊坡工程應按下列3種設計工況進行設計：

（1）持久設計工況：主要為邊坡正常運用工況，此時應采用基本組合設計。

（2）短暫設計工況：包括施工期缺少或部分缺少加固力；缺少排水設施或施工用水形成地下水位增高；運行期暴雨或久雨、或可能的泄流霧化雨，以及地下排水失效形成的地下水位增高；水庫水位奠降等情況。此時應采用基本組合設計。

（3）偶然設計工況：主要為遭遇地震、水庫緊急放空等情況，此時應采用偶然組合設計。

根據《水電水利工程邊坡設計規范》（DL/T5353-2006），邊坡為B類水庫邊坡，級別為Ⅱ級，邊坡設計安全系數標準：持久1.15～1.05，短暫1.10～1.05，偶然1.05～1.00。

3.2 穩定性計算成果

采用極限平衡法中的下限解法，進行二維極限平衡法分析。其中，折線滑面采用考慮力與力矩平衡條件的摩根斯頓－普萊斯法（Morgenstern-Price），圓弧滑面采用簡化畢肖普法（Bishop）。根據前期監測成果可知，拉古滑坡主要的破壞模式為整體沿基覆界線剪出、中部至底部沿基覆界線剪出和中部沿覆蓋層局部剪出。沿滑坡體Ⅰ區選取一處順滑坡體剖面分別按蓄水前后持久、短暫、偶然3種工況進行安全系數計算，計算結果見表3。

表3 蓄水前后滑坡體安全系數計算成果

從表3可知，蓄水前后，各工況下I區滑坡體整體和局部的安全系數均滿足規范要求。其中，中部沿覆蓋層局部剪出滑動模式安全系數不高，僅高于規范要求下限值，安全裕度不大，處于臨界穩定狀態，存在局部蠕滑變形的可能。計算結果與實際監測情況相符。

4 結論

（1）滑坡體深部滑坡主要發育在Ⅰ區中下部和Ⅱ區中部，滑移帶不連續，其主要滑動方式為局部次級滑體牽引，其可以作為滑坡穩定性分析和評價的前提。

（2）Ⅰ區1-1剖面上部變形以深層滑移為主，下部變形為深層滑移，滑坡體變形位移受上游沖溝切割作用影響較大。

（3）蓄水前后，各工況下I區滑坡體整體和局部的安全系數均滿足規范要求。其中，中部沿覆蓋層局部剪出滑移模式安全系數不高，僅高于規范要求下限值，安全裕度不大，處于臨界穩定狀態，存在局部蠕滑變形的可能。