王仁槐 寇學文

摘? 要：庫岸穩(wěn)定對水庫蓄水后的安全和長期運行有重要影響，將常見庫岸類型劃分為基巖庫岸、土質(zhì)庫岸、巖土混合庫岸，探討圖解法和卡丘金法計算水庫塌岸的結果差異，利用圖解法對某水庫塌岸進行預測。結果顯示，隨著距壩距離的增加，塌岸寬度先增加后減少，擬合曲線大致呈拋物線型，同時塌岸寬度大于100 m的嚴重塌岸貢獻塌岸總方量的89%。該文的數(shù)據(jù)分析方法及研究成果對類似工程勘察設計、施工具有一定的借鑒意義。

關鍵詞：庫岸；穩(wěn)定性預測；水利工程；工程勘察；圖解法

中圖分類號：TV522? ? ? 文獻標志碼：A? ? ? ? ? 文章編號：2095-2945（2024）19-0082-04

Abstract： The stability of the reservoir bank has an important influence on the safety and long-term operation of the reservoir after impoundment. The common reservoir bank types are divided into bedrock reservoir bank， soil reservoir bank and geotechnical mixed reservoir bank. The difference between the results of the graphic method and the Kaqiujin method for calculating reservoir bank collapse is discussed. The graphic method is used to predict the bank collapse of a reservoir. The results show that with the increase of the distance from the dam， the width of bank collapse increases at first and then decreases， and the fitting curve is roughly parabolic， and the serious bank collapse whose width is more than 100 m contributes 89% of the total amount of bank collapse. The data analysis method and research results of this paper can be used for reference to the survey， design and construction of similar projects.

Keywords： reservoir bank; stability prediction; hydraulic engineering; engineering survey; graphic method

庫岸崩塌通常發(fā)生在庫水位劇烈變化時，當庫水位迅速上升時，周圍的地下水位也隨之上升，使地下水和庫水共同作用于巖體和庫岸，對巖體和庫岸產(chǎn)生物理和力學作用[1]，表現(xiàn)為巖體強度及抗剪強度降低，再加上波浪的淘刷，容易發(fā)生崩塌，使庫岸由單一坡度變?yōu)樯习枷峦沟膹推麦w[2]；庫水位下降威脅邊坡穩(wěn)定的機理主要體現(xiàn)在孔隙水壓激增，水位下降后帶來的水頭差引發(fā)了滲流拖拽效應，甚至可能導致滲透破壞，在剪出口形成一個應力集中區(qū)域，該應力集中區(qū)域會向坡體內(nèi)部和上部擴張，最終誘發(fā)邊坡整體失穩(wěn)[3]。降雨是誘發(fā)庫岸崩塌的重要因素之一，降雨的入滲，一方面增加了不穩(wěn)定巖體的重量，從而增加了巖體的下滑力；另一方面對軟弱巖層產(chǎn)生較強的軟化作用，加劇了崩塌的形成[4]。

庫岸穩(wěn)定性是一項復雜的系統(tǒng)工程，庫岸崩塌的產(chǎn)生既有確定性也有不確定性。巖體參數(shù)內(nèi)摩擦角φ、黏聚力c、重度γ對崩塌是否發(fā)生具有決定性意義[5]，但影響庫岸穩(wěn)定性的因素多而復雜，且各因素的影響程度也不相同，其相互間又存在著一定的牽連，各種因素的影響很難做出經(jīng)典數(shù)學模型加以統(tǒng)一量度[6]。隨著數(shù)學理論的發(fā)展，神經(jīng)網(wǎng)絡方法、模糊綜合評判方法、灰色理論、模糊理論、層次分析法和信息量法等開始用于庫岸穩(wěn)定性評價，這些評價方法取得了重要的進展，但仍然存在一定的局限[7]，實際工程中使用較少。

目前，庫岸塌岸預測方法有以佐洛塔廖夫為代表的條件類比圖解法、以康德拉捷夫為代表的數(shù)學分析法、平衡剖面法、動力法、徐瑞春的塌岸預測圖解的若干修正和超前信息法等[8]。其中，卡丘金圖解法的實質(zhì)是依據(jù)實測的洪枯水位變幅帶各類巖性岸坡長期穩(wěn)定坡角，根據(jù)幾何關系用圖解法求解岸坡最終塌岸預測寬度，其精度取決于計算參數(shù)的選定。該方法在山區(qū)河道型水庫的塌岸預測中，適用于水庫的中游和上游地帶，而預測水庫下游的結果往往與實際相差甚遠，實際塌岸寬度要比預測的小得多，卡丘金圖解法的預測結果偏于安全[8]，因此，在實際工程中應用最為廣泛。

某水庫庫區(qū)山高谷深，溝壑縱橫，庫區(qū)人煙稀少、偶有耕地?；鶐r地層巖性為砂巖、泥巖、粉砂巖。河谷多呈寬“U”字型，庫岸出露地層巖性主要為三疊系上統(tǒng)胡家村組（T3h）、更新統(tǒng)風積黃土（Q3eol）、殘坡積（Q4dl+el）和崩坡積層（Q4dl+col）。本文利用卡丘金圖解法對水庫運行后的塌岸長度進行分析預測，根據(jù)塌岸長度分析預測結果進行塌岸嚴重程度分級，最后，依據(jù)庫區(qū)地形地貌、地質(zhì)特征、地質(zhì)構造等開展庫岸穩(wěn)定性評價。

1? 庫岸類型分布

根據(jù)水庫區(qū)地形地貌、地層結構與岸坡的組合關系，將庫岸分為基巖庫岸、土質(zhì)庫岸、巖土混合庫岸，其中巖土混合庫岸又可根據(jù)基巖坡度分為A、B 2個亞類，水庫蓄水后，A類向B類發(fā)展。圖1為庫區(qū)岸坡結構類型示意圖。

基巖庫岸基巖面高程高于正常蓄水位，基巖巖性為三疊系上統(tǒng)胡家村組（T3h）的砂巖、粉砂巖泥巖地層，地質(zhì)構造簡單，地層產(chǎn)狀近于水平。砂巖、粉砂巖夾泥巖地段岸坡較陡，約50°，局部近直立，形成高陡巖壁，天然岸坡穩(wěn)定性較好，臨溝未見大的滑坡和崩塌體。

黃土庫岸廣泛分布上更新統(tǒng)風積黃土（Q3eol），岸坡表層主要由黃土類土組成，上部岸坡較緩，坡度一般在10～30°，下部為三疊系的砂巖、粉砂巖夾泥巖地層。

巖土混合庫岸A基巖出露高程低于正常蓄水位，基巖巖性為砂巖、粉砂巖、泥巖，地層產(chǎn)狀近于水平。砂巖、粉砂巖夾泥巖地段岸坡較陡，約50°，局部近直立，形成高陡巖壁，岸坡穩(wěn)定性較好。上覆更新統(tǒng)風積黃土（Q3eol），受庫水位影響，存在塌岸問題。

巖土混合庫岸B岸坡下部基巖巖性以砂巖、粉砂巖夾泥巖為主，地層產(chǎn)狀近于水平。岸坡坡度在15～25°，岸坡坡腳附近為河流堆積物和黃土坡積物（主要是Q4dl+el和Q4dl+col），低于正常蓄水位而被淹沒。由于基巖出露高程低于正常蓄水位，上覆更新統(tǒng)風積黃土（Q3eol），受庫水位影響，存在塌岸問題。

庫區(qū)岸坡結構類型分布如圖2所示。庫區(qū)庫岸主要以基巖為主，占比接近70%，其余3類合計占比約30%，現(xiàn)狀條件下，庫岸總體保持穩(wěn)定。

2? 塌岸主要影響因素

2.1? 地層巖性

地層巖性是影響塌岸是否發(fā)生的最主要因素。由碎石土、砂礫石組成的庫岸及基巖庫岸抗沖刷能力強，塌岸發(fā)育緩慢。由黃土、砂土等組成的庫岸抗沖刷能力弱，一般隨著蓄水初期庫水位迅速上升，侵蝕、搬運、堆積劇烈，直至庫岸穩(wěn)定。本水庫易產(chǎn)生塌岸的是巖土混合庫岸A，其下部基巖具備一定抗沖刷能力，上部土質(zhì)容易發(fā)生塌岸，向巖土混合庫岸B發(fā)展。

2.2? 水庫運行方式

一般情況下，多年調(diào)節(jié)水庫水位變化較慢，對庫岸影響較小。年內(nèi)調(diào)節(jié)水庫水位高低頻繁波動，塌岸反復進行，持續(xù)時間長，水下坡角逐步減少至穩(wěn)定坡腳，根據(jù)月潭水庫[9]及新疆某水庫[10]的經(jīng)驗，該值通常小于20°。本水庫的主要功能為調(diào)蓄，對庫岸的影響可視為年內(nèi)調(diào)節(jié)水庫。

2.3? 波浪

波浪對塌岸的影響主要表現(xiàn)為波浪對庫岸的淘蝕和對塌落物質(zhì)進行搬運[11]。波浪的高度與風速、風在自由水面的吹程及水深等因素有關[12]。波浪越高，對岸坡的破壞力就越強。

2.4? 其他因素

影響塌岸發(fā)育的其他因素還有凍融、降雨等。岸坡土體孔隙和裂隙中的水結冰后體積膨脹產(chǎn)生冰劈作用，并破壞了土體結構，解凍后，土體強度降低，致使岸坡發(fā)生破壞。降雨會增加不穩(wěn)定土體的容重，同時降低土體黏聚力，導致庫岸穩(wěn)定系數(shù)迅速下降從而引起邊坡失穩(wěn)。

3? 塌岸預測

水庫塌岸預測的目的是對水庫運行后的最終塌岸寬度進行分析預測，評價塌岸對工程水庫運行、對塌岸區(qū)和周邊環(huán)境的影響，并提出防治措施。目前常用的方法有圖解法和卡丘金法，盧欣[13]曾采用2種方法分別計算同一庫區(qū)24個剖面的塌岸寬度，按照圖解法計算結果由小到大順序排列，結果如圖3所示，結果顯示當利用圖解法計算塌岸寬度達到29 m以上時，卡丘金法與圖解法計算結果相對一致。

根據(jù)水庫運行特點和庫岸結構類型，初步考慮利用圖解法進行岸坡塌岸的初步預測評價，預測參數(shù)的選取主要通過現(xiàn)狀調(diào)查，類比同類工程實測數(shù)據(jù)確定。浪高計算采用官廳水庫公式[14]，水下穩(wěn)定坡角α、水上穩(wěn)定坡角β按表1取值。本工程土質(zhì)庫岸、巖土混合庫岸A、巖土混合庫岸B容易產(chǎn)生塌岸的主要為粉砂土，計算塌岸分布情況如圖4所示。

根據(jù)計算結果，庫區(qū)塌岸寬度大多在35～150 m之間，在該區(qū)間圖解法與卡丘金法計算結果相對一致。隨著距壩距離的增加，塌岸寬度先增加后減少，擬合曲線大致呈拋物線型，曲線方程為y=-13.391x2+43.842x+110.24，塌岸寬度最大處發(fā)生在距壩距離3.2 km左右。

發(fā)生塌岸的庫岸類型主要是土質(zhì)庫岸和巖土混合庫岸A，計算14處庫岸中，土質(zhì)庫岸和巖土混合庫岸A各5處，其他4處為土質(zhì)庫岸與巖土混合庫岸A的“結合體”。巖石庫岸和巖土混合庫岸B相對穩(wěn)定，基本不發(fā)生或很少發(fā)生塌岸。巖土混合庫岸A發(fā)生塌岸，經(jīng)過侵蝕、搬運、堆積發(fā)育為巖土混合庫岸B。土質(zhì)庫岸會通過塌岸，逐步發(fā)育為基巖庫岸及巖土混合庫岸。

根據(jù)塌岸寬度將塌岸分為輕微塌岸（0～40 m）、中等塌岸（40～100 m）、嚴重塌岸（大于100 m），各類型塌岸占比為1∶3∶6，同時嚴重塌岸集中分布在距壩4.0 km以內(nèi)。

土質(zhì)庫岸與巖土混合庫岸A的“結合體”產(chǎn)生塌岸方量占比約為40%，巖土混合庫岸A產(chǎn)生塌岸方量占比約50%，土質(zhì)庫岸產(chǎn)生塌岸方量僅有10%。此外，嚴重塌岸（塌岸寬度大于100 m）貢獻塌岸總方量的89%。

4? 結論

本文介紹了水庫庫岸常見的4種類型，并就某水庫的庫岸類型分布進行了研究，數(shù)據(jù)顯示該水庫的基巖庫岸占比超過70%，水庫庫岸整體穩(wěn)定；同時探討了圖解法和卡丘金法計算水庫塌岸的結果差異，認為當計算塌岸寬度達到29 m以上時，卡丘金法與圖解法計算結果相對一致，本庫區(qū)塌岸寬度多為35～150 m之間；根據(jù)計算結果，嚴重塌岸（大于100 m）貢獻了絕大部分塌岸方量。本文的數(shù)據(jù)分析方法及研究成果可為類似工程勘察設計、施工提供參考，具有一定的借鑒意義。

參考文獻：

[1] 崔中興，柴軍瑞，仵彥卿.西北地區(qū)水庫黃土庫岸穩(wěn)定分析[J].水土保持通報，2004（1）：26-28.

[2] 程昌華，陳學兵.波浪對庫岸坍塌變形的試驗研究[J].重慶交通學院學報，1994（S1）：54-59.

[3] 樊任華，程建設，樊義永.庫水位驟降作用下庫岸滑坡形成機理研究[J].水利建設與管理，2018，38（5）：59-63.

[4] 王志旺，楊健，張保軍，等.水庫庫岸滑坡穩(wěn)定性研究[J].巖土力學，2004（11）：1837-1840.

[5] 李強，管昌生.庫岸滑坡穩(wěn)定可靠性分析中若干規(guī)律的探討[J].巖石力學與工程學報，2002（7）：999-1002.

[6] 羅軼，任光明，王志紅，等.某水庫庫岸穩(wěn)定性的模糊綜合評判[J].長江科學院院報，2011，28（6）：67-70.

[7] 李鵬岳，倪化勇，王春山，等.基于改進層次分析-可拓學模型的庫岸穩(wěn)定評價[J].人民長江，2018，49（5）：52-57.

[8] 張梁，王俊杰，閻宗嶺.山區(qū)庫岸塌岸預測方法綜述[J].重慶交通大學學報（自然科學版），2010，29（2）：227-232.

[9] 吳炳富.月潭水庫庫岸穩(wěn)定問題研究[J].治淮，2020（6）：38-41.

[10] 石永磊.水庫庫岸穩(wěn)定與塌岸預測研究[J].黑龍江水利科技，2020，48（11）：43-46，56.

[11] 宋岳，段世委，陳書文.官廳水庫塌岸影響因素分析[J].水利水電工程設計，2004（1）：34-37.

[12] 碾壓式土石壩設計規(guī)范：SL 274—2020[S].2020.

[13] 盧欣.圖解法和卡丘金法庫岸塌岸分析對比研究[J].江西水利科技，2017，43（2）：84-87.

[14] 吳亮，譚健敏，姜志水.水工設計中波高的計算方法探討[J].科技創(chuàng)新與應用，2016（11）：213.