基于Scoops3D模型的區(qū)域庫岸邊坡穩(wěn)定性分析

趙志陽 楊雪琪 宋揚

摘 要：庫岸邊坡失穩(wěn)滑移是工程中常見的地質災害，針對庫水位變動及降水影響下的區(qū)域性庫岸邊坡穩(wěn)定性問題，引入Scoops 3D模型，采用瞬態(tài)非飽和滲流理論計算降雨入滲，并利用Boussinesq理論計算庫水位勻速變動情況下的浸潤線，綜合分析區(qū)域性庫岸邊坡的穩(wěn)定性狀態(tài)。以三峽水庫萬州段為例，采用3種工況進行對比分析，表明區(qū)域內庫岸邊坡的穩(wěn)定性對降雨的響應更為明顯，計算結果一定程度上反映了庫岸邊坡的穩(wěn)定性時空分布。

關鍵詞：庫岸滑坡;穩(wěn)定性;降雨;庫水位;Scoops 3D模型

中圖分類號：TV698.1 ? 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.04.028

Abstract：Reservoir landslide is one of the main geological hazards in the project. Aiming at the stability of regional reservoir bank slope under the influence of reservoir water level fluctuation and precipitation， Scoops 3D model was introduced to calculate rainfall infiltration by using transient unsaturated seepage theory and Boussinesq theory to calculate the saturation line under the condition of uniform variation of reservoir water level. Based on this， the stability state of regional reservoir bank slope was comprehensively analyzed. Taking Wanzhou section of the Three Gorges Reservoir as an example， three sets of working conditions were used for comparative analysis. The results show that the stability of the reservoir bank slope in the region responds is more obviously to rainfall， and the calculation results reflect the temporal and spatial distribution of the stability of the reservoir bank slope to a certain extent.

Key words： reservoir landslide; stability; rainfall; reservoir water level; Scoops 3D model

長江三峽工程建成后對航運、發(fā)電、灌溉等起著非常重要的作用，但是三峽水庫庫水位變幅大，岸坡受水流沖刷持續(xù)崩退，對岸坡的穩(wěn)定性造成了一定的影響，據(jù)統(tǒng)計，自20世紀90 年代以來，僅長江中下游河道的崩岸就多達數(shù)百處，崩岸段總長達1 520 km，占兩岸岸線總長的35.7%[1]。針對岸坡失穩(wěn)機理，眾多學者進行了深入研究，如劉磊等[2]將區(qū)域穩(wěn)定性計算模型引入岸坡穩(wěn)定性分析中，深入分析了區(qū)域性岸坡穩(wěn)定性;黃達等[3]綜合分析了長江航道塔坪H2岸坡的變形規(guī)律與失穩(wěn)機理;張景昱等[4]基于水–巖循環(huán)作用，建立了巖體強度劣化模型，對岸坡的長期穩(wěn)定性進行了分析;李卓等[5]設計了岸坡滑坡試驗裝置，以龍江水電站近壩岸坡為例進行了針對性分析;李景娟等[6]對彎曲航道水流淘刷作用下的岸坡穩(wěn)定性進行了分析;吳生才[7]分析了混合型生態(tài)修復技術在岸坡防護中的應用。上述研究在分析區(qū)域性庫岸邊坡穩(wěn)定性時，很少同時考慮降雨與庫水位變動對岸坡的影響，眾多研究表明，降雨與庫水位變動對岸坡的穩(wěn)定性起著極其重要的作用，因此綜合分析庫水位變動與降雨共同影響下的區(qū)域性岸坡穩(wěn)定性具有重要意義。

基于此，筆者引入美國地質調查局Reid等[8]研發(fā)的Scoops 3D模型（三維斜坡穩(wěn)定性計算模型），采用瞬態(tài)非飽和滲流理論計算降雨入滲，并利用Boussinesq理論計算庫水位勻速變動情況下的浸潤線，綜合分析庫水位與降雨共同影響下的區(qū)域性庫岸邊坡穩(wěn)定性，最后以三峽水庫萬州段為例進行驗證。

1 庫岸邊坡穩(wěn)定性分析基本原理

1.1 Scoops 3D模型

Scoops 3D是采用Fortran語言編寫的一個開源計算程序[8]，可系統(tǒng)地搜索并計算整個數(shù)字高程模型（DEM）中數(shù)百萬計的三維潛在滑動面的安全系數(shù)。Scoops 3D采用剛體極限平衡分析中的簡化Bishop法或瑞典條分法作為理論基礎，以球面作為潛在滑面，可以模擬復雜地形并考慮各種荷載，如地下水、地震等。Scoops 3D在計算分析中，首先將潛在滑坡體網格化為眾多三維柱體，然后在DEM上生成球心，并據(jù)此生成半徑不等的大量球體，以球體與DEM相交得到的球面作為潛在滑面，最后對球面網格化的潛在滑坡體采用剛體極限平衡法求解得到安全系數(shù)。三維柱體中的潛在滑面示意見圖1。

2 工程概況

2.1 工程地質概況

以三峽水庫萬州段為例建模計算。該段地層巖性以沉積巖為主，主要包括三疊紀和中生代侏羅紀頁巖、泥巖、灰?guī)r和砂巖等，研究區(qū)內寬闊向斜和緊密背斜間隔分布，無大型斷裂構造，地貌單元主要為構造剝蝕低山丘陵、侵蝕堆積階地以及堆積河漫灘等，區(qū)域內存在第四系松散堆積層，厚度為0.5～50.0 m。根據(jù)地質災害調查統(tǒng)計發(fā)現(xiàn)，研究區(qū)內受庫水位變化影響的大型堆積滑坡有87個，厚度為10～20 m，基本處于穩(wěn)定狀態(tài)。

2.2 庫水位變動特征

三峽水庫水位單日最大降幅曾被嚴格限制在60 cm，2016年10月中旬經專家組評審，決定將庫水位單日最大降幅增大到120 cm。三峽大壩自建成蓄水后，庫水位的漲落呈現(xiàn)周期性變化，大體可分為高水位期、緩慢蓄水期、快速蓄水期、水位波動期、快速下降期、緩慢下降期6個階段（見表1）。

2.3 降雨數(shù)據(jù)

本次庫岸邊坡穩(wěn)定性分析主要針對庫水位變動的快速下降期，根據(jù)萬州區(qū)1964—2013年共50 a的降雨資料，統(tǒng)計得到三峽水庫處于快速下降期的日最大降雨量和持續(xù)降雨天數(shù)，并根據(jù)Gumbel概率密度分布模型擬合日最大降雨量的概率密度曲線（見圖4），從而得到概率為2%（重現(xiàn)期50 a）的日最大降雨量為134 mm，概率為10%（重現(xiàn)期10 a）的日最大降雨量為100 mm。

2.4 巖土物理力學參數(shù)

按照巖土滲透系數(shù)量級將萬州段庫岸邊坡劃分為1區(qū)（滲透系數(shù)為10-7量級）、2區(qū)（滲透系數(shù)為10-6量級）、3區(qū)（滲透系數(shù)為10-5量級）3個區(qū)，根據(jù)邊坡巖性，將2區(qū)進一步細分為2-A、2-B、2-C 3個亞區(qū)，各區(qū)的巖土物理力學參數(shù)見表2。根據(jù)表2各區(qū)滲透性參數(shù)及Boussinesq潛水一維非穩(wěn)定滲流理論，即可計算確定浸潤線高程。

2.5 計算工況的確定

為了綜合分析降雨與庫水位變動對邊坡穩(wěn)定性的影響，先設置庫水位從175 m以60 cm/d的速率等速下降到159 m（工況0），在此基礎上再設置3種工況：①工況1，庫水位從159 m以60 cm/d的速率等速下降到145 m，降雨采用重現(xiàn)期為10 a的單日最大降雨量100 mm，降雨歷時6 h，降雨強度16.7 mm/h;②工況2，庫水位從159 m以120 cm/d的速率等速下降到145 m，降雨采用重現(xiàn)期為10 a的單日最大降雨量100 mm，降雨歷時6 h，降雨強度16.7 mm/h;③工況3，庫水位從159 m以60 cm/d的速率等速下降到145 m，降雨采用重現(xiàn)期為50 a的單日最大降雨量134 mm，降雨歷時6 h，降雨強度22.3 mm/h。

2.6 模型參數(shù)

Scoops 3D模型在計算過程中，需要建立不同的球面與DEM相交形成眾多滑動面，在此基礎上尋求穩(wěn)定性安全系數(shù)最小的滑動面作為潛在滑動面，這些球面的建立首先需要在DEM上形成眾多球心，球心的搜索范圍為三維搜索矩陣。由于DEM的最大高程為330 m，因此設置三維搜索矩陣的高程范圍為200～400 m，間距為5 m×5 m，共形成1 915×2 793個柵格，其主要參數(shù)見表3。

3 庫岸邊坡穩(wěn)定性分析及驗證

為了便于庫岸邊坡穩(wěn)定性分析，根據(jù)邊坡穩(wěn)定性安全系數(shù)，將邊坡穩(wěn)定性等級劃分為低、較低、中、較高、高5類（見表4）。

3.1 175 m初始水位計算結果及分析

在進行上述工況計算分析前，為了驗證本文方法的可行性和準確性，針對坡體內有地下水和無地下水兩種情形進行敏感性分析，并將其計算結果（見表5）與庫岸邊坡的實際穩(wěn)定狀態(tài)進行對比分析。

從表5可以看出，兩種情形下三峽水庫萬州段的庫岸邊坡穩(wěn)定性都處于較高狀態(tài)。無地下水時，較高及高穩(wěn)定性狀態(tài)的邊坡面積共5 314萬m2（占總面積的94.79%），較低及低穩(wěn)定性狀態(tài)的邊坡面積共152萬m2（占總面積的2.71%）。考慮地下水時，較高及高穩(wěn)定性狀態(tài)的邊坡面積共5 022萬m2，占總面積的89.58%，比不考慮地下水的減小了5.21%;較低及低穩(wěn)定性狀態(tài)的邊坡面積共324萬m2，占總面積的5.78%，比不考慮地下水的增大了3.07%。上述分析表明，考慮地下水后庫岸邊坡的穩(wěn)定性下降明顯，其中較低以下穩(wěn)定狀態(tài)的邊坡主要分布于三峽水庫萬州段的上游，通過地質調查發(fā)現(xiàn)，該區(qū)域覆蓋層較厚，其巖土體的抗剪強度明顯低于其他區(qū)域的，因此本文的計算結果與實際情況是相符的。

3.2 不同工況對比分析及驗證

將上述計算工況中擬定的庫水位下降速率及降雨強度等資料輸入Scoops 3D模型中，分別計算3種工況下庫岸邊坡的穩(wěn)定性狀態(tài)，結果見圖5及表6。

由圖5可知，隨著庫水位下降速率與降雨強度增大，庫岸邊坡的穩(wěn)定性狀態(tài)趨于惡化，低穩(wěn)定性等級的庫岸邊坡在部分地區(qū)出現(xiàn)了貫通現(xiàn)象，主要分布于上游的溪口鄉(xiāng)、燕山鄉(xiāng)和瀼渡鎮(zhèn)，因此實際工作中上述地區(qū)應當重點關注，防止出現(xiàn)險情。由表6可知，工況0時庫岸邊坡的整體穩(wěn)定性比較好，較高及高穩(wěn)定性狀態(tài)的邊坡面積占總面積的81.92%，較低及低穩(wěn)定性狀態(tài)的邊坡面積僅占總面積的12.59%;隨著庫水位繼續(xù)下降及考慮降雨的影響（工況1），較高及高穩(wěn)定性狀態(tài)的邊坡面積減小到總面積的69.77%，較低及低穩(wěn)定性狀態(tài)的邊坡面積增大到總面積的21.09%;繼續(xù)加大庫水位的下降速率（工況2），較高及高穩(wěn)定性狀態(tài)的邊坡面積減小到總面積的63.06%（比工況1減小了6.71%），較低及低穩(wěn)定性狀態(tài)的邊坡面積增大到總面積的27.51%（比工況1增大了6.42%）;在工況1的基礎上加大降雨強度（工況3），較高及高穩(wěn)定性狀態(tài)的邊坡面積減小到總面積的55.89%（比工況1減小了13.88%），較低及低穩(wěn)定性狀態(tài)的邊坡面積增大到總面積的34.58%（比工況1增大了13.49%）。通過上述分析可知，相比于工況1，增大降雨強度（工況3），庫岸邊坡的穩(wěn)定性狀態(tài)惡化更為明顯，因此庫岸邊坡的穩(wěn)定性對降雨強度的變化更為敏感，這與近年來庫岸邊坡的地質災害事件常發(fā)生于降雨期是相吻合的。

4 結 論

引入Scoops 3D模型，采用瞬態(tài)非飽和滲流理論計算降雨入滲，并以Boussinesq理論計算庫水位勻速變動情況下的浸潤線，綜合分析庫水位與降雨共同影響下的區(qū)域性庫岸邊坡穩(wěn)定性。以三峽水庫萬州段為例，庫水位為175 m時，針對有無地下水影響兩種情形下的庫岸邊坡穩(wěn)定性計算結果，對該方法進行了驗證，在此基礎上設置了3組對比工況，分析了庫水位快速下降和降雨對庫岸邊坡穩(wěn)定性的影響，結果表明：庫岸邊坡的總體穩(wěn)定性較好，工況0情形下，較高及高穩(wěn)定性狀態(tài)的邊坡面積占總面積的81.92%，較低及低穩(wěn)定性狀態(tài)的邊坡面積僅占總面積的12.59%;低穩(wěn)定性等級的庫岸邊坡主要分布于上游的溪口鄉(xiāng)、燕山鄉(xiāng)和瀼渡鎮(zhèn)，實際工作中上述地區(qū)應當重點關注，防止出現(xiàn)險情;通過多工況對比分析，發(fā)現(xiàn)庫岸邊坡的穩(wěn)定性對降雨強度的變化更為敏感，這與近年來庫岸邊坡的地質災害事件常發(fā)生于降雨期相吻合。

參考文獻：

[1] 高清洋，李旺生，楊陽，等.長江中下游河道崩岸研究現(xiàn)狀及展望[J]. 水運工程，2016，8（12）：99-105.

[2] 劉磊，殷坤龍，王佳佳，等.降雨影響下的區(qū)域滑坡危險性動態(tài)評價研究：以三峽庫區(qū)萬州主城區(qū)為例[J].巖石力學與工程學報，2016，35（3）：558-569.

[3] 黃達，顧東明，陳智強，等. 三峽庫區(qū)塔坪H2古滑坡臺階狀復活變形的庫水-降雨耦合作用機制[J].巖土工程學報，2017，39（12）：2203-2211.

[4] 張景昱，宛良朋，潘洪月，等. 考慮水-巖作用特點的典型岸坡長期穩(wěn)定性分析[J]. 巖土工程學報，2017，39（10）：1851-1858.

[5] 李卓，何勇軍，李宏恩，等. 前期降雨作用下邊坡滑坡模型試驗[J]. 河海大學學報（自然科學版），2016，44（5）：400-405.

[6] 李景娟，段祥寶，謝羅峰. 彎曲航道水流特性及岸坡穩(wěn)定性研究[J]. 水運工程，2016，12（5）：137-147.

[7] 吳生才. 混合型生態(tài)修復技術在內河航道岸坡防護中的應用[J]. 水運工程，2011，8（9）：106-110.

[8] REID M E， CHRISTIAN S B， BRIEN D L， et al. Scoops 3D Software to Analyze 3D Slope Stability Throughout a Digital Landscape[EB/OL].（2015-04-10） [2017-12-03]. http：//dx.doi.org/10.3133/tm14A1.

[9] 劉磊.三峽水庫萬州區(qū)庫岸滑坡災害風險評價研究[D].武漢：中國地質大學（武漢），2016：22-41.

[10] SAVAGE W Z， GODT J W， BAUM R L. A Model for Spatially and Temporally Distributed Shallow Landslide Initiation by Rainfall Infiltration[G]//SMILLE T. Debris-Flow Hazards Mitigation： Mechanics， Prediction and Assessment. Rotterdam： [s. n.]， 2003： 179-187.

[11] 薛禹群.地下水動力學原理[M].北京：地質出版社，1986：46-49.

[12] 陳崇希，林敏.地下水動力學[M].武漢：中國地質大學出版社，1999：5-14.

【責任編輯 張華興】