張萬清， 徐玉華， 鄭 銳， 吳澤鵬

(1.福建省水土保持工作站，福建 福州 350001；2.福州大學(xué)土木工程學(xué)院，福建 福州 350003)

0 引 言

隧道工程等建設(shè)項目在施工過程中因爆破等會有棄土棄石產(chǎn)生。這些渣土結(jié)構(gòu)松散，在堆積的不同階段中會形成具有分層現(xiàn)象的坡面，可能會產(chǎn)生失穩(wěn)破壞。有許多的令人觸目驚心的事故說明降雨是觸發(fā)滑坡問題的主要原因之一[1]。

有關(guān)渣土場邊坡穩(wěn)定已經(jīng)進行了大量的研究，不少學(xué)者分析了內(nèi)摩擦角等土壤參數(shù)因素對邊坡穩(wěn)定性意義[2,3]。劉建偉針對渣土場邊坡容易發(fā)生坑蝕的問題，收集了渣場在不同坡度位置的渣場數(shù)據(jù)，并測量了相關(guān)的物理指標，并進行了直接剪切試驗并分析原因[4]。劉暉等在分析受到暴雨影響的基礎(chǔ)上，采用GEO-SLOPE軟件對棄渣場整體穩(wěn)定性進行計算，應(yīng)用瑞典圓弧法進行堆渣穩(wěn)定性分析[5]。王光輝通過實驗室試驗，使用極限平衡法計算邊坡的穩(wěn)定性，并使用不同的方法計算安全系數(shù)，分析了渣土堆積層在降雨影響下的邊坡穩(wěn)定性[6]。

以一互層渣土堆積層邊坡模型為例，分析不同降雨峰值分布對孔隙水壓力分布和邊坡危險性的影響。初步探討了不同降雨條件對渣土堆積層邊坡穩(wěn)定性的影響，為生產(chǎn)建設(shè)項目棄渣場的風險管理提供了一定的理論參考。

1 考慮降雨影響邊坡穩(wěn)定分析方法

1.1 非飽和滲流及邊坡穩(wěn)定計算

非飽和滲流模擬中，由GeoStudio的SEEP/W模塊進行，滲流的控制方程式為：

(1)

滲透系數(shù)的曲線如下式[7]：

kw=ksΘp

(2)

式中：kw表示非飽和土的滲透系數(shù)；ks表示飽和滲透系數(shù)；Θ=(θw/θs)表示體積含水量與飽和體積含水量的比率，p為滲透系數(shù)曲線中的擬合參數(shù)。

計算中所采用的土水特征曲線如圖1表示：

圖1 土水特征曲線

邊坡穩(wěn)定性計算選取應(yīng)用廣泛的Morgenstern-Price法。非飽和土抗剪強度計算選擇Frelund等提出的抗剪強度公式：

τ=c′+(σn-ua)tanψ′+(ua-uw)tanψb

(3)

式中，ua為孔隙空氣壓力；uw為孔隙水壓力；ψb為由負孔隙水壓力的作用而提高的強度；ψ′為有效摩擦角；c′為有效粘聚力。

1.2 計算模型及邊界條件

通過實地調(diào)研，發(fā)現(xiàn)渣土堆積層組成主要是由巖石在生產(chǎn)爆破后的開挖、開挖后經(jīng)坡面下落、運輸汽車-推土機排放所決定的，考慮到渣土各次堆放時間跨度大的影響，其分層現(xiàn)象就會變得明顯，主要包括水平分層和傾斜分層。選取角度為19°的傾斜分層，建立滲流分析的有限元模型，根據(jù)施工階段分為5個傾斜分層。邊坡基礎(chǔ)計算模型如圖2所示，坡面總體高度為20m，坡角為26°。地下水位高程在右側(cè)為7m，左側(cè)為11m。有限元采用結(jié)構(gòu)化四邊形單元，單元尺寸為0.5m。

當降雨強度和邊坡飽和滲透系數(shù)相比較，前者較大時，此時的雨水是不完全入滲的形式，這時把他定義為定水頭邊界；而后者較大時，雨水完全入滲，此時將坡面視為為流量邊界。斜坡模型的底部以及左側(cè)與基巖接觸的其余部分都是不透水邊界。

圖2 土水特征曲線

由于通常渣土堆積層成分復(fù)雜常夾雜石塊，重度變化較大，因此計算時需考慮重度的變異性。本數(shù)值模擬黏聚力與內(nèi)摩擦角均視為正態(tài)分布，將重度視為對數(shù)正態(tài)分布，計算參數(shù)參考工程實際參數(shù)水平。

1.3 降雨工況設(shè)計

依據(jù)福建省2018年降雨特性分析，短時極端降雨最大雨量達274.0mm[8]。考慮到在實際生活中，降雨是個不均勻的過程，因此按照時間和降雨量之間的動態(tài)的關(guān)系，設(shè)置過程降雨總量為300mm，降雨歷時為5d的四種降雨工況：前鋒型、后鋒型、中鋒型、等強形。假設(shè)四類工況的降雨量和降雨時間均相同，探究不同雨量峰值出現(xiàn)順序?qū)υ炼逊e層的危險性的影響，具體降雨工況如圖3所示。

圖3 按降雨量峰值變化雨型工況圖

2 計算結(jié)果分析

2.1 滲流場分析

不同雨型工況以邊界條件的形式施加在邊坡上，模擬不同工況過程中保證降雨歷時、總降雨量等條件不變。為反應(yīng)孔隙水壓隨時間及空間的變化規(guī)律，提取四類工況降雨過程中邊坡三個分層界面處各節(jié)點土體孔隙水壓，孔隙水壓隨時間變化如圖4。

不同的降雨工況下分層界面處孔隙水壓力的變化規(guī)律有著明顯的差異性。對于工況2后鋒型及工況4中鋒型雨型來說，初期雨量較小，因此界面處剛開始孔壓增長緩慢，隨中后期雨量增大后迅速增大，降雨結(jié)束時，后鋒型雨型相較其他幾種雨型工況，孔壓增長較大；工況1削弱型和工況3等強型在初期，降雨較為集中，土體含水量和孔隙水壓的變化也隨之迅速增長。

當土體飽和滲透系數(shù)大于滲透系數(shù)時，初期降雨的入滲升讀內(nèi)的土體孔隙水壓力隨降雨量增大而增大。持續(xù)入滲過程中，削弱型工況出現(xiàn)孔隙水壓變小的現(xiàn)象，回落較為明顯。后鋒型雨型邊坡孔壓在降雨過程中持續(xù)增大，對于等強型降雨，降雨的過程中，淺層的孔隙水壓增大到一定的程度后逐漸趨于穩(wěn)定，而其之后的降雨量的增加，僅僅增加其入滲的深度，對于淺層的孔隙水壓并沒有明顯的影響。

圖4 四種工況下孔壓隨時間變化圖

2.2 危險性評價

由于不同降雨條件下堆積層內(nèi)孔壓分布不同，內(nèi)部的應(yīng)力狀態(tài)也隨之改變，進而會對堆積層滑坡的危險性產(chǎn)生不同的影響。但由于四種工況的降雨歷時較短，因此四種工況的失效概率較為接近。各工況下失效概率隨累計降雨量變化如圖5、6，四類雨型危險性評價表如表1所示。

圖5 不同雨量峰值工況下失效概率隨時間變化圖

圖6 不同工況下失效概率隨累計降雨量變化圖

表1 四類雨型危險性評價表

由表1可知，在降雨模擬結(jié)束之后，在后鋒型降雨條件下，邊坡失效概率達到了27.45%，相比于其余幾種工況，失穩(wěn)的可能性是最大的，危險程度屬于低危險，接近于中危險。前鋒型降雨工況的失效相比于后鋒型降雨工況低3.8%，且相比于其他幾種降雨條件，前鋒型降雨條件下的堆積層邊坡失效概率為最低。中鋒型和等強型降雨后，堆積層滑坡的危險性已經(jīng)接近于遞增型，說明等強型降雨相對更危險。

3 結(jié) 論

降雨的類型以降雨量的峰值的出現(xiàn)的時間先后順序來進行分類。以降雨量的峰值的出現(xiàn)的時間先后順序可以把降雨類型分為前鋒型、后鋒型、中鋒型和等強型。在相同的降雨的持續(xù)時間和相同的降雨總量的情形下，降雨量的峰值出現(xiàn)的時間先后對邊坡產(chǎn)生了明顯的影響，當降雨期的峰值集中在降雨的后期時，會出現(xiàn)的淺層孔壓回落的現(xiàn)象，危險性是最高，而當降雨期的峰值集中在降雨的前期時，危險性較低。而當中鋒型雨型位于前鋒和后鋒之間時，接近于后鋒型，但是中鋒型的降雨入滲的深度比強型的小，因此他的滑動面的范圍也相對而言比較小，其邊坡相對安全。