曹秋麗

(山東省濟寧市水利事業(yè)發(fā)展中心,山東 濟寧 272000)

1 概 述

許多地區(qū)在秋季大雨季節(jié),降雨引起的邊坡破壞較為常見。雨水滲透導(dǎo)致土壤非飽和區(qū)吸水能力下降,飽和區(qū)孔隙水壓力增加,最終導(dǎo)致邊坡破壞的發(fā)生[1]。同時,徑流沿坡面流動有一定的侵蝕作用,增加了入滲水量,進(jìn)一步降低了邊坡穩(wěn)定性。對于大面積集水區(qū)地形產(chǎn)生的徑流,從地勢較高的區(qū)域向凹地排水流動,在強降雨期間,徑流匯集的路堤邊坡發(fā)生破壞的風(fēng)險較高[2]。在分析雨水入滲引起的非飽和土邊坡失穩(wěn)時,為了簡化,通常忽略徑流的影響,無法完全再現(xiàn)降雨-徑流入滲的實際過程,尤其是在降雨強度遠(yuǎn)大于入滲能力的強降雨條件下[3]。因此,同時考慮地表徑流和土壤水流至關(guān)重要。

近年來,地表徑流和土壤水流的相互作用機制引起許多學(xué)者的關(guān)注[4]。為了降低計算復(fù)雜性和提高計算效率,基于淺水方程簡化的運動學(xué)公式嘗試用于徑流模擬。劉昱[5]采用具有時變結(jié)構(gòu)的ForecastNet徑流預(yù)測模型,有效提高了徑流模擬及多步預(yù)測精度。李根[6]采用Navier-Stokes方程和Richards方程,分別模擬了地表徑流和土壤水流,利用有限體積法求解,結(jié)果表明地表徑流對土體入滲有較大提高。

本文對地表徑流、土壤水流與產(chǎn)流時間的關(guān)系進(jìn)行分析,建立地表徑流和土壤水流的耦合模型,模擬強降雨條件下降雨、徑流和滲透之間的關(guān)系,評估降雨-徑流引起的邊坡不穩(wěn)定性,以確定邊坡易發(fā)生失穩(wěn)的風(fēng)險區(qū)域。

2 工程背景

2021年夏季,我國南部某省受到臺風(fēng)影響,并帶來強降雨,引發(fā)了強烈的山體滑坡、河堤坍塌。根據(jù)調(diào)查得到的相關(guān)數(shù)據(jù),該地區(qū)在3天內(nèi)的累計降雨量462mm,峰值降雨強度50mm/h。利用所提出的耦合模型,對該地區(qū)的邊坡穩(wěn)定性進(jìn)行分析,并基于激光剖面數(shù)字高程數(shù)據(jù),建立該地區(qū)天然山區(qū)的3D模型,對地表和土壤水流分析及邊坡穩(wěn)定性進(jìn)行評估。

3 模擬方法

在降雨情況下,地表徑流和土壤水流之間通常相互作用,雨水入滲包括前期的降雨入滲和后期的徑流入滲兩個階段。Green-Ampt模型用于估算地表的滲透能力,以確定地下水流分析的邊界條件。與在一定水頭壓差下的土壤導(dǎo)水率不同,滲透能力被定義為在給定條件下土壤吸收水的速率,在暴雨開始時滲透能力較強,后期隨著表層土壤飽水而迅速下降。在降雨事件發(fā)生的前期,降雨強度小于地表的滲透能力,Richards方程可直接用于地表處的邊界條件求解;當(dāng)降雨強度超過土壤的入滲能力時,地表就會產(chǎn)生徑流,通量邊界條件變?yōu)閴毫吔鐥l件。

本研究通過使用有限元軟件COMSOL Multiphysics,利用COMSOL的PDEs模塊,建立二維表面流動模型和理查茲方程,計算地表徑流深度和速度的分布以及土壤中孔隙水壓力和體積含水量的分布。利用COMSOL的固體力學(xué)模塊,建立三維土力學(xué)模型,計算有效應(yīng)力。由于該模型主要判斷邊坡是否發(fā)生破壞,以便及時進(jìn)行預(yù)警工作,沒有考慮大變形,因此該模型采用線性彈性模型。

地下水流模型和土壤力學(xué)模型以兩種方式耦合:①將與體積含水量相關(guān)的體荷載函數(shù)應(yīng)用于線彈性土壤力學(xué)模型,以表明水分變化對自重和應(yīng)力分布的影響;②基于Bishop有效應(yīng)力方程,考慮體積含水量變化對孔隙水壓力的影響來評估有效應(yīng)力。具體公式如下:

σ′=σ-u+χ(u2-uw)

(1)

式中:σ′為有效應(yīng)力,kPa;σ為總應(yīng)力,kPa;u為孔隙空氣壓力,kPa;uw為孔隙水壓力,kPa;χ為吸力系數(shù),根據(jù)飽和程度,在0～1范圍內(nèi)取值,與有效飽和度、殘余飽和度有關(guān)。

該模型由4部分組成,即風(fēng)化花崗巖、土壤、河堤路和瀝青層。土壤、路堤等邊坡的曼寧系數(shù)值取為0.3m/s1/3,瀝青層的曼寧系數(shù)為0.1m/s1/3。根據(jù)歷史資料測得的平均值,將地下水位的初始條件設(shè)置為距離地表5.5m。對于數(shù)值模型的邊界條件,底面為壓力邊界,側(cè)面為無滲流邊界,表面為降雨和徑流相結(jié)合的氣候邊界。模擬時間為2021年夏季該地區(qū)7月25日17:00至7日28日20:00,共計75h,時間步長1h。

4 結(jié)果與討論

4.1 徑流對邊坡穩(wěn)定性的影響

為了分析地表徑流和土壤水流對邊坡穩(wěn)定性的影響,研究兩種不同的情況,即考慮徑流和不考慮徑流的邊坡不穩(wěn)定性分析。根據(jù)數(shù)值模擬情況,該地區(qū)每小時降雨量達(dá)到最大值時,大量的水聚集在兩個位置,使更多的水滲入河堤,這正是河堤的破壞點。圖1為位置1和位置2的徑流深度隨時間的變化,以及每小時的降雨強度和入滲能力。由圖2可以看出,隨著降雨入滲量的增加,土壤的入滲能力急劇下降。在約23h時,降雨強度超過了滲透能力,這是沿坡面產(chǎn)生徑流的時間。因此,徑流深度較大會增加滲透,并導(dǎo)致沿斜坡坡面的侵蝕,從而導(dǎo)致邊坡失穩(wěn)。在考慮徑流的情況下,根據(jù)位置1和位置2的不穩(wěn)定區(qū)域,可以確定安全區(qū)域;不穩(wěn)定區(qū)域延伸至路堤深處,在位置1處達(dá)到3m,在位置2處達(dá)到1.5m。在不考慮徑流的情況下,邊坡破壞僅發(fā)生在位置1和位置2的周邊范圍較小的區(qū)域,無法滿足實際情況。

圖1 徑流深度隨時間的變化

圖2 降雨量與滲入量的關(guān)系

在考慮徑流的情況下,徑流對邊坡表面的水量起到重要影響作用,導(dǎo)致孔隙水壓力更顯著地增加。路堤表層孔隙水壓力由負(fù)值增加到正值,導(dǎo)致邊坡失穩(wěn)。相反,在不考慮徑流的情況下,孔隙水壓力的增量很小。從邊坡實際穩(wěn)定情況來看,考慮徑流的分析結(jié)果與實際破壞區(qū)域較好重合,而不考慮徑流的不穩(wěn)定區(qū)域小于實際情況。其中,位置1處的實際滑面達(dá)到13m,比模擬結(jié)果更深,其原因可能是忽略了速度較大的徑流引起的侵蝕效應(yīng),高速徑流會沖刷帶走地表土壤顆粒。結(jié)果表明,徑流對路堤邊坡的破壞有顯著影響,是強降雨期間大面積邊坡穩(wěn)定性分析中需要考慮的關(guān)鍵因素[7]。

為了進(jìn)一步研究徑流對邊坡穩(wěn)定性的影響,分析位置1和位置2的有效飽和度、孔隙水壓力和局部安全系數(shù)隨時間變化的規(guī)律及原因,結(jié)果見圖3、圖4和圖5。由圖3、圖4、圖5可知,滲透水會導(dǎo)致有效飽和度和孔隙水壓力的增加,從而導(dǎo)致非飽和土壤的吸力減少,最終降低局部安全系數(shù)。此外,考慮徑流時,有效飽和度和孔隙水壓力的增加以及局部安全系數(shù)的減少更為顯著。產(chǎn)流后,有效飽和度和孔隙水壓力顯著增加,局部安全系數(shù)急劇下降,表明來自上游的徑流允許更多的水滲入路堤,以減少路堤巖土體的吸水性。因此,在徑流匯集區(qū)域發(fā)生筑堤邊坡破壞的可能性遠(yuǎn)大于沿線的其他位置。

圖3 有效飽和度

圖4 孔隙水壓力

圖5 局部安全系數(shù)

徑流是導(dǎo)致河堤邊坡破壞的關(guān)鍵因素之一,當(dāng)流速超過臨界侵蝕速度時,坡面同時可能發(fā)生侵蝕。臨界侵蝕速度定義公式如下:

(2)

式中:V為臨界侵蝕速度,m/s;d50為沉積物的中值直徑,mm,本文試驗中取值為2.7mm;位置1和位置2的臨界侵蝕速度計算為0.53m/s。

進(jìn)一步分析徑流對河堤坡面侵蝕的影響,圖6為位置1和位置2路堤邊坡上的徑流速度。徑流速度超過臨界侵蝕速度時,會導(dǎo)致沿路堤邊坡的嚴(yán)重侵蝕。由圖6可知,徑流速度隨著降雨的前期而逐漸增加,在57h時達(dá)到峰值,在降雨后期又逐漸下降。

圖6 徑流速度

4.2 瀝青層對徑流和邊坡穩(wěn)定性的影響

目前,城市化發(fā)展使雨水流滲入地面的模式發(fā)生了巨大變化。道路瀝青層等不透水覆蓋系統(tǒng),減少了滲入土壤的水量;防滲層通過減少滲透,有效防止了邊坡的破壞失效。因瀝青層有限的滲透能力會導(dǎo)致更多的地表徑流量,本研究討論了瀝青層對徑流和邊坡不穩(wěn)定性評估的影響。為此,將數(shù)值模型中瀝青層的輸入?yún)?shù)替換為其下方材料土壤或路堤的參數(shù),比較考慮和不考慮瀝青層兩種情況下的徑流和邊坡穩(wěn)定性的模擬結(jié)果。在考慮瀝青層的情況下,位置1和位置2的徑流速度相較于不考慮瀝青層時顯著增加,在強降雨期間,瀝青層的存在增加了徑流速度,超過臨界侵蝕速度的高速徑流覆蓋了更多的區(qū)域,從而造成侵蝕和破壞。

圖7為道路上的徑流深度和徑流速度。由圖7(a)可知,由于瀝青層的滲透能力較低,徑流出現(xiàn)的時間更早。兩種情況下的最大徑流深度差異不大,而考慮瀝青層的最大徑流速度高于未考慮瀝青層的最大徑流速度,見圖7(b)。因此,在強降雨期間,瀝青層加速了徑流速度,并導(dǎo)致更多的水流沿著道路流動,使路堤加速破壞。另一方面,瀝青層的低滲透性防止了更多的水滲入土壤地面,反而會減少破壞區(qū)域,而在不考慮瀝青層的情況下,導(dǎo)致大量的水滲入。因此,在考慮瀝青層的情況下,道路區(qū)域周圍的局部安全系數(shù)相對較大,表明低滲透性的瀝青層能夠保持邊坡穩(wěn)定性,并降低道路周圍部分的破壞規(guī)模。在位置2處,實際的邊坡破壞僅發(fā)生在河堤路堤上,與考慮瀝青層的模擬結(jié)果吻合良好。

圖7 瀝青層的徑流深度與速度

5 結(jié) 論

本文通過分析降水引起的邊坡地表徑流和土壤徑流,對某邊坡案例的地表和土壤水流進(jìn)行了分析,并進(jìn)行了邊坡穩(wěn)定評估。結(jié)論如下:

1)地表徑流和土壤水流的耦合模型可以再現(xiàn)當(dāng)降雨強度超過集水區(qū)土壤地面的滲透能力時產(chǎn)生的徑流,可用于評估降雨-徑流引起的邊坡不穩(wěn)定性。

2)滲透水會導(dǎo)致有效飽和度和孔隙水壓力的增加,從而導(dǎo)致非飽和土壤的吸力減少,并最終降低局部安全系數(shù)。徑流對路堤邊坡的破壞有顯著影響,是強降雨期間大面積邊坡穩(wěn)定性分析中需要考慮的關(guān)鍵因素。徑流速度超過臨界侵蝕速度時,會導(dǎo)致沿河堤邊坡的嚴(yán)重侵蝕。

3)人工設(shè)施如鋪設(shè)瀝青層的坡面,會導(dǎo)致徑流速度的增加。同時,瀝青材料的弱滲透性可有效阻止降雨的滲透,以限制邊坡的破壞規(guī)模。