裂隙性質對含裂隙基坑土坡滲流特性影響數值模擬研究*

朱 余 佳

(安徽審計職業學院 工程管理系，合肥 230001)

0 引 言

邊坡土在降雨循環，自然風化作用下，不可避免地出現各種形狀，各種深度，各種角度的裂隙，這些裂隙的存在影響了改變了邊坡的均質度，使得降雨入滲呈現強各向異性的性質[1-3]。雨水將優先通過這些滲流通道滲透至邊坡內部，影響整個邊坡的滲流場，從而影響邊坡的穩定性[4]。因此，對裂隙邊坡在降雨條件隙的滲流特性及穩定性影響規律的研究將對邊坡的排水設計及滑坡災害的防治提供一定的參考作用。

國內外學者從試驗及數值模擬方面對裂隙降雨展開研究，試驗方面:David[5]對某實際邊坡進行了概化，通過室內等比例試驗對裂隙傾角對邊坡滲流特性的影響進行了研究；Chen[6]對含節理裂隙邊坡降雨條件下的滲流穩定性進行了室內試驗。數值模擬方面:王培清[7]考慮了紅黏土邊坡的裂隙各向異性，對降雨條件下的滲流特性展開數值模擬；闕云[8]利用Geo-studio研究了黏土邊坡裂隙滲流特性及穩定性規律；張英瑩[9]利用Matlab程序及圖像編寫技術對膨脹土在循環降雨下的邊坡失穩規律進行了研究。以上研究對裂隙邊坡滲流及邊坡穩定性質進行了一定有益的探討，并進行了一些規律性的總結，但是對于邊坡裂隙性質的全面研究較少，裂隙存在著不同各向異性特征，如裂隙傾角[10]、裂隙方向[11]、裂隙深度[12]及裂隙各向異性程度[13]等，對邊坡的滲流特性影響巨大，在降雨條件下可能會改變邊坡整個滲流場的分布。

采用數值方法，基于非飽和滲流原理，建立了不同裂隙性質的有限元數值模型，對降雨條件下含裂隙邊坡進行了非飽和滲流計算，重點研究了裂隙部位的滲流場變化規律，研究結果為全面認識不同裂隙性質對邊坡滲流場的影響規律提供了一定的參考。

1 數值模型及計算理論

1.1 數值模擬方法

數值模擬在對邊坡裂隙的模型上，主要有兩種考慮方法:

(1) 裂隙屬于邊坡內部的不連續邊界，當降雨入滲時，雨水將優先經過裂隙進行入滲，因此可以將裂隙作為定水位或者是孔壓邊界的施加部位，在有限元計算中根據實際情況施加相應的邊界條件[14]。但是由于裂隙尺寸較小，裂隙的間隙遠小于模型網格的最小尺寸，建模存在困難，尤其是對于邊坡表面裂隙發育充分的區域。

(2) 因為裂隙作為邊坡表層的優勢滲流通道，因此可以將裂隙均質化處理，將裂隙等效成為連續介質，單對其參數單獨賦予，主要考慮裂隙的各向異性及強滲透性[15]，由于這種方法建模簡便，因此采用這種建模方法進行裂隙性質的討論。

1.2 裂隙性質參數

對于實際裂隙發育的邊坡表層，裂隙可能存在不同的方位角α，定義為裂隙的走向與豎直方向的夾角；裂隙的豎直向滲透系數與水平向滲透系數的比值η=ky/kx，為裂隙的各向異性程度；裂隙的長度L，定義為裂隙底端點與裂隙與坡面交點連線長度；裂隙位置，為不失一般性，分別取裂隙位置位于上部，中部及下部的情況。

裂隙性質參數示意圖如圖1所示。裂隙單元的局部坐標與整體坐標的換算公式如下:

(1)

式(1)中:ky為裂隙的豎向滲透系數，kx為裂隙水平向滲透系數，而α則為方位角，是裂隙的走向與豎直方向的夾角。

圖1 裂隙性質參數示意圖Fig. 1 Schematic diagram of fracture property parameters

1.3 非飽和滲流方程

非飽和滲流微分方程為

(2)

式(2)中，kr是透水率；kij是滲透張量；hc是水頭；Q是源匯；C(hc)是容水度；θ是水頭函數；n是孔隙率；Ss是單位貯水量。

2 計算模型及參數

2.1 模型概況

模型選擇某高速公路裂隙邊坡為典型案例，邊坡高10 m，為一二級邊坡，在邊坡中部設置一1 m寬的馬道。為減小邊界的影響，分別向邊坡左右拓寬25 m建模，模型深度取為12 m。整個模型一共劃分為2 145個節點，2 213個單元。其中為研究不同裂隙位置的影響，特取如圖2所示的上部裂隙，中部裂隙與下部裂隙進行討論研究，裂隙采用第1.1節中的均質化方法進行建模。

圖2 模型及網格劃分Fig. 2 Model and grid division

2.2 邊界條件

降雨條件下的邊坡滲流計算需要指定邊坡的滲流初始條件，對于算例而言，取如圖藍色虛線所示的地下浸潤線計算的滲流場作為整個工況的初始滲流場。

邊界條件如下:ab為23 m定水頭邊界，df為12 m定水頭邊界，bcde為降雨入滲邊界，降雨強度設置為0.02 m/d，降雨時間設置為7 d，計算總時間設置為12 d，而af為不透水邊界。

2.3 土體及裂隙的非飽和滲流參數

(a) 滲透系數函數

3 計算與結果分析

3.1 計算工況

表1 計算工況Table 1 Calculation conditions

3.2 計算結果

對不同工況下的裂隙處孔壓及體積含水率隨深度的變化規律進行統計，結果列于本節，為對比分析，特計算無裂隙工況作為對比工況(工況E)。

3.2.1 完整邊坡降雨滲流特性(工況E)

工況E為不含裂隙邊坡降雨條件下的孔壓及體積含水率變化規律，如圖4所示。

(a) 孔壓變化

完整邊坡相對應的虛擬裂隙面孔壓與體積含水率變化呈現單一變化規律，在降雨時刻邊坡的表層孔壓迅速增大到飽和(孔壓到達0 kPa或者體積含水率達到0.5)，停雨后孔壓緩慢減小，濕潤鋒不斷向邊坡深度發展，這與文獻[18]的研究結果類似，說明采用的數值方法可信。

3.2.2 不同各向異性程度下的裂隙邊坡滲流規律

工況A不同各向異性程度下的裂隙面滲流特性如圖5所示。

(a) 孔壓變化

由圖5可見，不同滲透比下的邊坡裂隙面滲流特性差異較大。與不含裂隙邊坡工況(工況E)相比，降雨下裂隙面孔壓整體增大，濕潤鋒進展更深入，裂隙各向異性程度越大，相同時間下的孔壓上升越快，對于各向異性為1的情況而言，降雨條件下孔壓或者體積含水率只有在表層達到最大，而在裂隙內部孔壓或者體積含水率逐漸減小，在裂隙末端達到最小，對于各向異性大于1的情況，降雨情況下裂隙內部的孔壓或者體積含水率隨高程幾乎呈現線性增大的規律，停雨后裂隙內部的孔壓或者體積含水率迅速減小，雨水逐漸向邊坡深部滲透，值得注意的一點是，對于各向異性程度大于100的情況，不同各向異性程度對孔壓的影響幾乎一致，這是由于各向異性程度大于100時，垂直向的滲透系數已經足夠大，雨水將主要沿著裂隙走向進行滲透，因此水平向的滲透行為幾乎對裂隙滲流特性無影響。

3.2.3 不同裂隙深度對邊坡滲流特性的影響

工況B為不同裂隙深度下的邊坡裂隙面滲流規律，如圖6所示。

(a) 孔壓變化

在工況B的計算條件下，不同裂隙面的滲流特性較為類似，降雨條件下的孔壓隨裂隙深度呈現線性變化，在裂隙底部孔壓達到最大，隨后迅速減小到初始孔壓，停雨后孔壓逐漸減小，對于體積含水率變化而言，降雨條件下裂隙內部體積含水率達到最大，并在裂隙深度處保持飽和至裂隙底部，裂隙深度以下體積含水率迅速折減，停雨后表層體積含水率逐漸降低。不同深度下，裂隙的深度越深，裂隙內部的高孔壓及體積含水率區越大，同時最大孔壓也越大。

3.2.4 不同裂隙面下的邊坡滲流特性

裂隙位于邊坡上部，中部及下部(圖1)時，裂隙面的孔壓變化如圖7所示。

(a) 上部裂隙

由圖7可見，裂隙處于邊坡不同位置下的裂隙面滲流特性有較大差異，對于上部裂隙而言在裂隙面孔壓呈現線性增大的規律，在裂隙高程以下，孔壓則逐漸降低，停雨后裂隙面的孔壓逐漸降低，而裂隙高程以下的孔壓則逐漸增大，體現了濕潤鋒不斷下移的過程。對于中部裂隙而言，滲流特性與上部裂隙一致，但是孔壓整體較上部裂隙要大，對于下部裂隙而言，裂隙面對孔壓的變化無影響，因為邊坡下部不但受到降雨入滲的影響，同時受到上部雨水的排泄，使得邊坡下部的孔壓相對較大，因此孔壓整體呈現線性增大的趨勢。

3.2.5 不同裂隙傾角下的裂隙面滲流特性

工況D為不同裂隙傾角下的邊坡孔壓及體積含水率變化曲線:

由圖8可知:計算工況而言，裂隙傾角較小時(0°～30°)，孔壓與體積含水率變化則較為一致，而當裂隙傾角較大時(60°)，差異則較大，主要體現在孔壓及體積含水率變化幅度的降低，這是因為裂隙的傾角較小時，雨水的豎向滲透占據滲流行為的主要部分，因此接近于裂隙垂直的情況，孔壓及體積含水率的變化較為一致，即雨水更容易垂直入滲，裂隙內部的孔壓及體積含水率較大，當裂隙的傾角較大時，裂隙更傾向于水平分布，垂直于裂隙面的滲透系數較低，因此裂隙的垂向滲流行為被阻斷，因而孔壓及體積含水率整體上較小。

(a) 孔壓變化

4 結論與討論

基于均質化處理方法，將裂隙等效為不同各向異性，不同傾角，不同深度，不同位置的有限元模型，對各種工況下的降雨裂隙滲透特性進行了有限元計算，得結論:

(1) 降雨條件下含裂隙邊坡裂隙面孔壓及體積含水率上升幅度及范圍要大于完整邊坡，且濕潤鋒進展得更加徹底。

(2) 各向異性程度較小時，降雨條件下孔壓或者體積含水率只有在表層達到最大，而在裂隙內部孔壓或者體積含水率逐漸減小，在裂隙末端達到最小，各向異性程度較大時，裂隙內部的孔壓或者體積含水率隨高程幾乎呈現線性增大的規律，各向異性足夠大差異較小。

(3) 不同深度下，裂隙的深度越深，裂隙內部的高孔壓及體積含水率區越大，同時最大孔壓也越大。

(4) 裂隙位于邊坡上部及中部時，裂隙面孔壓呈現線性增大的規律，在裂隙高程以下，孔壓則逐漸降低，降雨對下部裂隙內部的孔壓變化無明顯影響，孔壓呈現線性增大趨勢。

(5) 裂隙角度與土體內部含水量與孔壓上升速率呈負相關，裂隙傾角較小時(0°～30°)，孔壓與體積含水率變化則較為一致，而當裂隙傾角較大時(60°)，差異則較大。

(6) 僅僅對裂隙邊坡的滲流特性進行了均質化處理，對于實際情況應當考慮裂隙擴展演化，同時考慮到復雜邊界的影響。