王銳 孫杰 趙紅軍 孔俊

摘要：為了探究滲流水文過程對(duì)非飽和土質(zhì)邊坡穩(wěn)定性的影響機(jī)制，建立了能夠模擬地下水流和分析邊坡穩(wěn)定性的耦合模型?；?種表述邊坡水文過程的HSB模型開展建模對(duì)比分析，揭示邊坡內(nèi)非飽和滲流對(duì)水文過程的影響機(jī)制;通過耦合水文模型與邊坡穩(wěn)定極限平衡模型，比較了淺層邊坡在不同水力滲流條件下，因鄰近水位升降所造成穩(wěn)定性系數(shù)變化的差異，探討非飽和側(cè)向滲流過程對(duì)淺層邊坡穩(wěn)定性的影響機(jī)制。分析結(jié)果表明，土體非飽和區(qū)的儲(chǔ)水、釋水能力和滲透能力影響著非飽和區(qū)域流的滲流流速大小;在淺層邊坡中，非飽和側(cè)向滲流作用可以促進(jìn)水體在邊坡內(nèi)部的遷移而引起邊坡內(nèi)部水位的快速抬高，致使邊坡穩(wěn)定性降低，并導(dǎo)致邊坡失穩(wěn)時(shí)間提前。研究結(jié)果可以為臨近地域水位變動(dòng)區(qū)邊坡工程實(shí)施提供相關(guān)技術(shù)參考。

關(guān)鍵詞：水文地質(zhì)學(xué);滑坡;非飽和流;邊坡穩(wěn)定性;流域水文模型

中圖分類號(hào)：P641文獻(xiàn)標(biāo)識(shí)碼：ADOI： 10.7535/hbgykj.2021yx03012

Abstract： In order to explore the influence mechanism of seepage hydrological process on the stability of unsaturated soil slope， a coupling model was established to simulate groundwater flow and analyze slope stability. Based on three Hillslope-storage Boussinesq models describing the hydrological process of the slope， the modeling comparative analysis was carried out to reveal the influence mechanism of unsaturated seepage in the slope on the hydrological process. By coupling the hydrological model with the slope stability limit equilibrium model， the variations of stability coefficients of shallow slope caused by the rise and fall of adjacent water levels under different hydraulic seepage conditions were compared， and the influence mechanism of unsaturated lateral flow process on the stability of shallow slope was discussed. The results show that the water storage and release capacity and the permeability of soil in the unsaturated area affect the seepage velocity in unsaturated area. In the shallow slope， unsaturated lateral seepage action can promote the water movement in the slope， which will lead to the rapid rise of the internal water level of the slope and the reduction of the slope stability， resulting in the slope instability time. The research results can provide relevant technical reference for the implementation of slope engineering in the adjacent area with water level fluctuation.

Keywords：hydrogeology; landslide; unsaturated flow; slope stability; watershed hydrological model

中國幅員遼闊、人口眾多，是世界上受到地質(zhì)災(zāi)害影響較為嚴(yán)重的國家之一。隨著人類生產(chǎn)生活需求的不斷增長，開采資源的程度逐年遞增，改造自然的力度不斷加大，致使地質(zhì)災(zāi)害的發(fā)生頻率逐年增加?；率且环N受多因素影響的地貌改造過程[1]，也是最常見的一種地質(zhì)災(zāi)害類型[2]。據(jù)統(tǒng)計(jì)，全球每年有數(shù)千人因滑坡事故而死亡或受傷[3]，給生產(chǎn)生活帶來不可估量的損失。

國內(nèi)外研究成果表明邊坡內(nèi)部水文過程對(duì)邊坡穩(wěn)定性具有不可忽視的影響[4-5]，若在邊坡失穩(wěn)預(yù)警工作中忽略非飽和流的作用，邊坡內(nèi)部水位的計(jì)算將出現(xiàn)低估問題，進(jìn)而影響邊坡穩(wěn)定安全系數(shù)，使得邊坡預(yù)警系統(tǒng)難以準(zhǔn)確預(yù)測(cè)滑坡發(fā)生時(shí)間和滑塌面積的大小。選擇合適的流域水文模型有利于更加高效地描述淺層邊坡內(nèi)部水文過程。相較于二維/三維的變飽和模型（如Richards方程），即使面對(duì)小規(guī)模問題也需要進(jìn)行大量的計(jì)算[6]，Hillslope-storage Boussinesq（簡(jiǎn)稱HSB）模型由于擁有更高的計(jì)算效率和魯棒性而受到較多關(guān)注。目前滑坡的相關(guān)研究都側(cè)重于分析降雨情況下的邊坡穩(wěn)定性[7-10]，對(duì)以側(cè)向流為主要誘因的淺層滑坡研究較少。

本文創(chuàng)新性地引入HSB水文模型來模擬邊坡內(nèi)部的水位變化過程，引入3種對(duì)非飽和流處理方式不同的水文模型：Troch-HSB模型[11]、Hilberts-HSB模型[12]和Kong-HSB模型[13]。其中Troch-HSB模型為考慮飽和流影響的模型;Hilberts-HSB模型進(jìn)一步考慮了非飽和區(qū)的垂向流作用;Kong-HSB模型則在Hilberts-HSB模型的基礎(chǔ)上又考慮了非飽和區(qū)的側(cè)向流。通過比對(duì)上述3種HSB模型由于對(duì)非飽和流處理方式不同而在模型的理論參數(shù)上產(chǎn)生的差異，分析非飽和流在邊坡內(nèi)部水文過程中的影響。再將3種模型與邊坡穩(wěn)定性模型耦合到一起，模擬淺層邊坡中的水文過程和邊坡穩(wěn)定性變化過程，進(jìn)而分析非飽和側(cè)向滲流對(duì)邊坡穩(wěn)定計(jì)算結(jié)果的影響，從而為類似工程提供有價(jià)值的參考依據(jù)。

第3期王銳，等：非飽和流對(duì)邊坡內(nèi)水文過程及穩(wěn)定性的影響分析河北工業(yè)科技第38卷

1模型的比較與分析

1.1變飽和邊坡水文模型的比較

為將Boussinesq地下水方程用于描述邊坡內(nèi)部的水文過程，TROCH等[11]引入有效含水量S（x，t）來描述邊坡的幾何特征和內(nèi)部介質(zhì)的孔隙度，并假設(shè)垂直坡面方向上水壓滿足靜壓分布且垂直于x的水平方向無流動(dòng)，得到了一維的水動(dòng)力方程。

1.2特征參數(shù)分析

由式（4）可知，非飽和流的作用可分為2類：一是非飽和區(qū)垂向流，非飽和區(qū)中垂向流對(duì)邊坡儲(chǔ)水和排水都起著重要作用，在水位變化過程中非飽和區(qū)也會(huì)發(fā)生吸水或排水，從而影響坡內(nèi)部的水位變化和含水量的空間分布，這主要反映在參數(shù)有效孔隙度f中;二是非飽和區(qū)側(cè)向流，會(huì)促進(jìn)邊坡內(nèi)的水分運(yùn)移，這主要反映在參數(shù)滲透流量q中。通過比較式（4）可知：

圖2所示為Kong-HSB模型和Troch-HSB模型中有效孔隙度的比值隨水位的變化。從圖中可以看出，對(duì)于砂土，在坡內(nèi)水位較低時(shí)，Kong-HSB模型和Troch-HSB模型的有效孔隙度幾乎相等，在坡內(nèi)較高水位才會(huì)出現(xiàn)明顯差異;而對(duì)于壤土，Kong-HSB模型和Troch-HSB模型有效孔隙度的差異隨著水位的抬高不斷拉大，即使當(dāng)坡內(nèi)水位趨近于0時(shí)，Kong-HSB模型的有效孔隙度也僅相當(dāng)于Troch-HSB模型的0.9倍。

可見，在相同含水量變化時(shí)考慮非飽和流的影響下，水位的波動(dòng)更大，同時(shí)在相同水位變化下，考慮非飽和流情況下的側(cè)向滲透流量更大。

2非飽和區(qū)對(duì)水文變化的影響機(jī)制

為揭示模型有效孔隙度參數(shù)f對(duì)邊坡內(nèi)部水位波動(dòng)的影響機(jī)制，設(shè)定初始水位分別為0.4，0.8，1.0 m，含水量變化取值為0.05 m3，分別利用Kong-HSB模型和Troch-HSB模型模擬邊坡內(nèi)部含水量的空間分布。2種模型在各工況的水位變化量如表1所列，含水量的空間分布如圖3所示，虛線引向坐標(biāo)軸表示的是邊坡內(nèi)部水位高度。從表1和圖3可以發(fā)現(xiàn)，在保證初始水分分布和含水量變化量相同的情況下，Kong-HSB與Troch-HSB模型水位變化差異會(huì)受到土壤屬性和邊坡初始水位的影響，這與圖2中反映的有效孔隙度的相對(duì)關(guān)系一致。

根據(jù)圖3與式（6）可知：在砂土中非飽和區(qū)持水性弱，非飽和區(qū)的含水量相對(duì)較小，Kong-HSB模型與Troch-HSB模型的有效孔隙度幾乎相等，所以2種模型計(jì)算得到的水位變化幾乎相同，只有當(dāng)初始水位足夠高時(shí)，有效孔隙度差異所導(dǎo)致的水位波動(dòng)差異才會(huì)開始顯現(xiàn);在壤土中，由于非飽和區(qū)厚度大且含水量也大，即使在初始水位很低的情況下2個(gè)模型的有效孔隙度也會(huì)有顯著差異，這種差異造成了如表1所示的2種模型水位波動(dòng)的不同，且當(dāng)初始水位升高后，2種模型的有效孔隙度的差異增大導(dǎo)致水位波動(dòng)差異隨之增大。

3非飽和流對(duì)邊坡穩(wěn)定性的影響

3.1邊坡穩(wěn)定性分析模型

相比于邊坡的長度和寬度，所討論的淺層邊坡厚度足夠小，可以采用極限平衡法結(jié)合半無限邊坡假定來分析邊坡穩(wěn)定性[16-17]。假定土條兩側(cè)側(cè)向力的影響可以忽略，即可以認(rèn)為圖4中條塊的側(cè)面法向力E、剪切力X的增量ΔE，ΔX均為0。

3.2非飽和流對(duì)邊坡穩(wěn)定系數(shù)的影響

將3種HSB水文模型與邊坡穩(wěn)定性模型耦合，形成了Troch-HSB-SM模型、Hilberts-HSB-SM模型和Kong-HSB-SM模型。將非飽和流作用所造成邊坡內(nèi)部水位和含水量分布差異的相關(guān)參數(shù)引入到邊坡穩(wěn)定性計(jì)算，分析邊坡穩(wěn)定安全系數(shù)變化和失穩(wěn)時(shí)間差異。

為分析河流筑壩、流域匯流情況下邊坡鄰近地域水位升降所引起的側(cè)向流動(dòng)對(duì)邊坡穩(wěn)定性的影響，以壤土為例設(shè)計(jì)了2種坡長不同的工況：考慮坡長分別為10，50 m時(shí)，邊坡失穩(wěn)的擾動(dòng)因素均為坡頂邊界水位抬升，詳細(xì)信息見表2，土壤參數(shù)設(shè)定如表3，其中的土壤參數(shù)α，β，n1，n2，Ks來自TROCH的研究[14]。

模擬結(jié)果顯示，隨著時(shí)間的推移，水分滲入邊坡引起水位抬高，模型得到的邊坡穩(wěn)定系數(shù)逐漸減小。在初始階段3種模型穩(wěn)定系數(shù)變化非常接近，這是因?yàn)榇藭r(shí)水位抬升所造成的水壓力占主導(dǎo)地位。之后對(duì)于工況1中較短的10 m壤土邊坡，Kong-HSB-SM模型預(yù)測(cè)的邊坡穩(wěn)定性的下降最快，也最早出現(xiàn)邊坡失穩(wěn)現(xiàn)象，其次是Hilberts-HSB-SM模型和Troch-HSB-SM模型。工況2中壤土邊坡坡長為50 m時(shí)，Kong-HSB-SM模型和Hilberts-HSB-SM模型預(yù)測(cè)的邊坡穩(wěn)定系數(shù)F幾乎同時(shí)以相似的速度下降并接近邊坡失穩(wěn)，說明非飽和側(cè)向流對(duì)邊坡穩(wěn)定性的影響也會(huì)受到邊坡長度的影響，邊坡長度越短，非飽和流誘發(fā)的邊坡失穩(wěn)越明顯，失穩(wěn)時(shí)間也越提前。

3.3非飽和流對(duì)邊坡穩(wěn)定影響的動(dòng)力學(xué)機(jī)制

根據(jù)圖7邊坡邊界流量-時(shí)間過程曲線來分析坡頂與坡腳的滲水過程，發(fā)現(xiàn)在入滲初期時(shí)，由于坡頂邊界處的水位不斷抬高產(chǎn)生的水力梯度使得此時(shí)的坡頂入滲量較大;而當(dāng)坡頂邊界處水位抬高至恒定時(shí)，坡頂一側(cè)地下水位依舊升高，水力梯度減弱使得坡頂入滲量減少;當(dāng)坡頂一側(cè)地下水位升至與其邊界水位近乎相同時(shí)，此階段水力梯度僅受坡度影響，此時(shí)的坡頂入滲量保持恒定。同樣在滲透初期，由于上邊界入滲水未到達(dá)坡腳，坡腳一側(cè)水位和滲出量保持恒定;隨著上游入滲水到達(dá)坡腳，坡腳水位抬高產(chǎn)生的水力梯度加大了坡腳處的滲出量;當(dāng)坡腳滲出量增大到一定值時(shí)，會(huì)使得坡頂入滲量與坡腳滲出量達(dá)到平衡。

根據(jù)圖7與上述分析，非飽和區(qū)側(cè)向流的入滲過程可分為3個(gè)階段，如圖8所示：第1階段為滲透初期，與坡腳非飽和區(qū)的滲出量相比，坡頂一側(cè)較大的水力梯度使更多的水分從坡頂非飽和區(qū)滲入邊坡，這個(gè)階段非飽和區(qū)側(cè)向流促進(jìn)水分滲入邊坡，增加邊坡含水量，降低邊坡穩(wěn)定性;第2階段，升高的坡頂?shù)叵滤皇蛊渌μ荻葴p小，此時(shí)坡頂非飽和區(qū)入滲量減小，而坡腳非飽和區(qū)滲出量保持不變或略有增加，此階段非飽和區(qū)側(cè)向流作用微弱;第3階段，坡腳一側(cè)水力梯度增大使坡腳非飽和區(qū)滲出量增加，此時(shí)坡腳非飽和區(qū)滲出量大于坡頂非飽和區(qū)滲入量，非飽和區(qū)側(cè)向流協(xié)助邊坡的排水，降低邊坡內(nèi)的水位，緩解邊坡的失穩(wěn)性。

而對(duì)滑坡失穩(wěn)問題而言，非飽和區(qū)側(cè)向流入滲過程一般處于第1或第2階段。根據(jù)上述理論可解釋圖6：對(duì)于10 m邊坡工況，在失穩(wěn)前非飽和區(qū)側(cè)向流長時(shí)間處于第1階段，非飽和側(cè)向流作用明顯，所以Kong-HSB-SM模型預(yù)測(cè)的邊坡穩(wěn)定性下降最快;而對(duì)于長度50 m邊坡工況，在失穩(wěn)前非飽和區(qū)側(cè)向流的作用長時(shí)間處于第2階段，此階段非飽和側(cè)向流作用微弱，所以Kong-HSB-SM模型和Hilberts-HSB-SM模型預(yù)測(cè)結(jié)果接近。

4結(jié)語

通過分析計(jì)算，比對(duì)了3種表述邊坡水文過程的HSB模型，分析非飽和滲流對(duì)邊坡內(nèi)部水文過程的影響和作用機(jī)制，并耦合HSB模型與邊坡穩(wěn)定極限平衡模型，模擬了淺層邊坡在不同水動(dòng)力條件下的安全系數(shù)變化情況，探討非飽和側(cè)向滲流對(duì)淺層邊坡穩(wěn)定性的影響。研究主要結(jié)論如下。

1）對(duì)于非飽和流顯著的壤土，在相同含水量變化量的條件下，考慮非飽和流模型的水位波動(dòng)更大;在相同水位的條件下，考慮非飽和流的模型水分運(yùn)移更多。

2）在以側(cè)向流為主要滑坡誘發(fā)因素的淺層滑坡中，壤土邊坡中非飽和流作用可以促進(jìn)水分在邊坡內(nèi)部的運(yùn)移，導(dǎo)致邊坡穩(wěn)定性的降低，忽視非飽和流的邊坡失穩(wěn)時(shí)間預(yù)測(cè)會(huì)滯后。

3）非飽和區(qū)側(cè)向流的作用可劃分成3個(gè)階段，第1階段非飽和區(qū)側(cè)向流促進(jìn)水分滲入邊坡，降低邊坡穩(wěn)定性;第2階段非飽和區(qū)側(cè)向流作用減弱;第3階段非飽和區(qū)側(cè)向流能協(xié)助邊坡的排水，降低邊坡內(nèi)的水位，緩解邊坡的失穩(wěn)性。

本文假定研究對(duì)象為均質(zhì)各向同性邊坡，但在自然環(huán)境中非均質(zhì)性邊坡更為常見，且本文未能考慮更復(fù)雜的邊坡幾何形態(tài)（如起伏、發(fā)散、收斂等特征），這些可能會(huì)改變邊坡內(nèi)部側(cè)向流特征，對(duì)研究結(jié)果有一定的影響;今后可以通過引入更多的參數(shù)來描述邊坡地質(zhì)和形態(tài)，探究其非飽和流的作用。另外，本文主要分析了不受降雨入滲影響的淺層滑坡，在后續(xù)研究中可考慮降雨和側(cè)向流的共同作用促成的淺層滑坡。

參考文獻(xiàn)/References：

[1]宋玉延. 滑坡治理雙排抗滑樁受力特征及設(shè)計(jì)計(jì)算方法研究[D]. 武漢： 武漢理工大學(xué)， 2013.

SONG Yuyan. Study on Force Characteristics and Design and Calculation[D]. Wuhan： Wuhan University of Technology， 2013.

[2]SINGH T N， VERMA A K， SARKAR K. Static and dynamic analysis of a landslide[J]. Geomatics Natural Hazards and Risk， 2010， 1（4）： 323-338.

[3]PETLEY D. Global patterns of loss of life fromlandslides[J]. Geology， 2012， 40（10）： 927-930.

[4]CHEN X Z， CUI Y F. The formation ofthe Wulipo landslide and the resulting debris flow in Dujiangyan City， China[J]. Journal of Mountain Science， 2017， 14（6）： 1100-1112.

[5]CONTE E， DONATO A， TRONCONE A. A simplified method for predicting rainfall-induced mobility of activelandslides[J]. Landslides， 2017， 14（1）： 35-45.

[6]BRESCIANI E， DAVY P， DREUZY J R. Is the Dupuit assumption suitable for predicting the groundwater seepage area in hillslopes？[J]. Water Resources Research， 2014， 50（3）： 2394-2406.

[7]陳靈芝， 古振榮， 李揚(yáng). 降雨入滲條件下非飽和土高邊坡穩(wěn)定性有限元分析[J]. 齊齊哈爾大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2020， 36（6）： 77-82.

CHEN Lingzhi， GU Zhenrong， LI Yang. Finite element analysis of high slope stability of unsaturated soil under rainfall infiltration[J]. Journal of Qiqihar University（Natural Science Edition）， 2020， 36（6）： 77-82.

[8]劉鋒， 胡江洋， 劉軍勇， 等. 連續(xù)降雨入滲作用下非飽和土路基穩(wěn)定性研究[J]. 路基工程， 2020（1）： 55-60.

LIU Feng， HU Jiangyang， LIU Junyong， et al. Study on the stability of unsaturated soil subgrade under continuous rainfall infiltration[J]. Subgrade Engineering， 2020（1）： 55-60.

[9]毋遠(yuǎn)召， 馬文禮， 魏占璽， 等. 降雨及滲流條件下層狀砂質(zhì)板巖邊坡變形破壞模式研究[J]. 河北工業(yè)科技， 2020， 37（4）： 230-236.

WU Yuanzhao， MA Wenli， WEI Zhanxi， et al. Research on deformation and failure modes of layered sandy slate slopes under rainfall and seepage conditions[J]. Hebei Journal of Industrial Science and Technology， 2020， 37（4）： 230-236.

[10]周娟， 朱峰. 中西部地區(qū)黏性土質(zhì)邊坡降雨動(dòng)態(tài)滲流的非飽和特性分析[J]. 河北工業(yè)科技， 2020， 37（4）： 273-279.

ZHOU Juan， ZHU Feng. Analysis on unsaturated characteristics of dynamic rainfall seepage of clay soil slope in Central and Western China[J]. Hebei Journal of Industrial Science and Technology， 2020， 37（4）： 273-279.

[11]TROCH P A， PANICONI C， van LOON E E. Hillslope-storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes： 1. Formulation and characteristic response[J]. Water Resources Research， 2003， 39（11）： 1316-1328.

[12]HILBERTS A G J， TROCH P A， PANICONI C， et al. Low-dimensional modeling ofhillslope subsurface flow： Relationship between rainfall， recharge， and unsaturated storage dynamics[J]. Water Resources Research， 2007， 43（3）： 3445-3459.

[13]KONG J， SHEN C J， LUO Z Y， et al. Improvement of thehillslope-storage Boussinesq model by considering lateral flow in the unsaturated zone[J]. Water Resources Research， 2016， 52（4）： 2965-2984.

[14]TROCH P A. Conceptual Basin-Scale Runoff Process Models for Humid Catchments： Analysis， Synthesis and Applications[D]. Ghent： Ghent University， 1993.

[15]van GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturatedsoils[J]. Soil Science Society of America Journal， 1980， 44（5）： 892-898.

[16]李寧， 許建聰， 欽亞洲.降雨誘發(fā)淺層滑坡穩(wěn)定性的計(jì)算模型研究[J]. 巖土力學(xué)， 2012， 33（5）： 209-214.

LI Ning， XU Jiancong， QIN Yazhou，et al. Research on calculation model for stability evaluation of rainfall-induced shallow landslides[J]. Rock and Soil Mechanics， 2012，33（5）： 209-214.

[17]彭子茂， 黃震. 模擬非飽和流的新型相對(duì)滲透系數(shù)模型[J]. 長江科學(xué)院院報(bào)， 2020， 37（7）： 115-119.

PENG Zimao， HUANG Zhen. A relative permeability model for unsaturated flow in consideration of hysteresis effect[J]. Journal of Yangtze River Scientific Research Institute， 2020， 37（7）： 115-119.