土石壩飽和—非飽和穩定滲流場影響分析

韓 雪

(通河縣水務技術服務中心，黑龍江 通河 150900)

?

土石壩飽和—非飽和穩定滲流場影響分析

韓 雪

(通河縣水務技術服務中心，黑龍江 通河 150900)

土石壩是當今最為常用的一種壩型，而防滲體系的構建對壩體穩定具有十分重要的作用，國內外的多起水利工程安全事故足以引起足夠的重視。進行土石壩的滲流分析需要將滲流場定義為飽和和非飽和區的共同體，過去的經驗已然證明二者的相輔相成。二者相比非飽和滲流的運動特性相對復雜，影響因子也較飽和滲流多。文章結合工程實例，運用有限元數值模擬出土石壩在飽和與非飽和滲流中滲流場的變化規律，得出正確合理的滲流場理論，研究成果在土石壩的防滲體系的構建上具有一定的參考價值。

土石壩；非飽和土；滲流；有限元

0 前 言

非飽和土是指經開挖、重塑和碾壓密實的土體。土石壩土體介質的非飽和特性至今依然是現代水利工程界的熱門課題[1]。因為其不僅廣泛存在于水利工程，農業領域也運用廣泛，如地下水位變化、土體水份滲透及蒸發等的變化對地下水的影響，其主要作用在土體的非飽和區域。在土木工程領域，雨水的滲透對地基穩定性的影響也要分析飽和-非飽和的滲流分布規律。

土石壩的滲流場由飽和區和非飽和區共同組成[1]，二者相互制約又相互聯系。影響因素諸多，降雨量、水位的變化等均能改變飽和區和非飽和區的作用區域，理論上單獨對非飽和區的滲流研究更能反映工程實際，假如單獨進行飽和區的滲流分析，直接加大了計算難度，伴隨著有時計算陷入死循環難以求解。

影響非飽和滲流運動的因素眾多，且形態也較復雜，比如固、液、氣三種不同形態、溫度、體積比、壓強和空氣壓力等，從力學角度范疇將這是多因子相互作用的共同結果，也是相互耦合的受力變化過程。目前國內隨著多年的研究發現，對于非飽和滲流參數的確定至關重要，參數的確定才能計算分析得出非飽和滲流的滲流場分布規律[2]。為使更加密切聯系工程實際，并為工程設計運行提供正確合理的理論依據，足以見土石壩飽和-非飽和滲流場的數值模擬分析的重要性。

1 飽和-非飽和滲流有限元解法

1.1 定解條件

對穩定滲流方程進行求解時，列出邊界條件；對非穩定滲流方程求解時，列出初始條件和全過程的邊界條件[3]。

基于流動的數學模型，可分為3類邊界條件：

1)第1類邊界條件是水頭邊界條件，假定已知邊界的水頭分布。表達式如下：

h|Γ1=H1(x,y,z,t)

(1)

2)第2類邊界條件是流量邊界條件，假定已知邊界的位勢函數，表達式如下：

(2)

3)第3類邊界條件是混合邊界條件，定義邊界外水頭差和水位能量呈線性關系，即：

(3)

式中：α、β為正常數。

1.2 飽和-非飽和滲流的基本微分方程

在非飽和土體中，通?？紫读黧w組成部分有水和氣，假定只考慮流水作用，結合達西定律和水體積變化值，建立非飽和滲流基本微分方程[4]，即：

(4)

假定總應力對時間導數為零，孔隙氣壓力同理，則上式表達為：

(5)

式中：m2w為壩基水的黏滯力變化下體積系數?？蓪⑹?4)變為

(6)

若假設是穩定滲流分析，即：

(7)

非飽和與飽和土體的在滲流方程表達式上雖相同，但非飽和滲透系數為變值，作用點變化，系數隨之改變。將上式展開，得出基本的非線性方程表達式如下：

(8)

1.3 飽和-非飽和滲流的有限元解法

國內目前對土體飽和與非飽和滲流的研究公式主要是運用Galerkin加權余量法，推導出二維有限元飽和-非飽和滲流問題公式，具體推導過程如下[3]：

首先已知近似解假定一個試探函數，列出加權余量方程，基于格林原理的方法消除二次微分，并將所得式代入假定函數及邊界條件，進而得出組合滲透矩陣，推導得到參數矩陣。以單元形式表達如下：

∫A[B]T[K][B]dA×{Hi}e+∫A[N]Tλ[N]dA

(9)

數值積分后得出：

(10)

在時間微分的處理中首先計算出特定時間的節點水頭，再利用有限差分法繼續求解[5]。得出的各單元數值代入式(5)，進而得出總體流動方程組。要求各單元各節點具有相似性，水頭必須保持一致。下面對各節點處的水頭{H}求解過程需依據總體流動方程組。各節點處的水頭值直接影響滲透系數，采用迭代法來實現。文章基于SEEP程序采用向后差分法進行迭代，在運算中設置容許誤差和迭代次數上限。該程序通過求解得出的節點收斂水頭值，進而計算得到滲流速度和孔隙水壓力。

2 工程實例

建于中國西南地區的某均質土壩，主要功能是調蓄和灌溉，兼以發電功能，大壩壩頂長度505.0m，壩頂高程854.0m，最大壩高68.0m。上、下游壩坡比為1∶3、1∶2.75。該壩的最大橫剖面圖如圖1所示[6]。大壩內排水為水平褥墊方式?；鶐r上部為強風化層，基巖下部巖體為相對隔水層。

圖1 大壩橫剖面圖

2.1 有限元模型與計算參數

坐標系定義如下：X軸為順河流方向上游指向下游為正；Y軸是豎直方向，指向向上為正，有限元網格剖分圖如圖2。

圖2 大壩有限元模型

模型的計算區域范圍選?。貉睾恿鞣较蚍謩e向上、下游各取68m；豎直方向從壩基向下取80m。

有限元網格劃分：該模型由顏色各異的5個分區組成，分別表示大壩、強弱透水地層、防滲帷幕和下游排水體五部分。計算軟件采用GeoStudio，模塊為SEEP/W。該模塊的特點是適用性強，能夠自動剖分單元，模型劃分單元選取三角形和四邊形單元，這樣劃分很大程度上能提高計算準確度。網格剖分后共生成單元2113個，結點2245個。

對大壩各分區材料進行了力學實驗，并結合類似工程實例檢測結果，該土石壩的各物理力學參數見表1[7]。

表1 土石壩各分區材料物理力學參數表

2.2 穩定滲流場對比分析

文章研究情況為正常工況下的穩定滲流場。定義邊界條件如下：壩基處認定為不透水邊界，即滲流量為零；在壩體上游面和左側強透水層高程定義總水頭145m，在壩體下游面定義總水頭94m，將上、下游的總水頭值定義為水頭邊界。

2.2.1 飽和穩定滲流場

定義材料類型：在運用SEEP/W軟件[8]分析時，將大壩的各分區材料滲透系數定義為飽和滲透的。

通過有限元計算分析，飽和穩定滲流場在正常工況下的總水頭、壓力水頭、孔隙水壓力分布圖及計算結果已得出工況下的浸潤線見圖3-5。

圖3 正常工況下飽和穩定滲流場總水頭等值線分布圖

圖4 正常工況下飽和穩定滲流場壓力水頭等值線分布圖

圖5 正常工況下飽和穩定滲流場孔隙水壓力等值線分布圖

2.2.2 飽和-非飽穩定滲流場

定義材料類型：根據以往浸潤線作為飽和-非飽和的分界線，結合工程實際定義大壩材料為非飽和狀態下的參數，其他4分區材料定義為飽和狀態下的材料參數。

通過有限元計算分析可得到飽和-非飽和穩定滲流場的總水頭、壓力水頭、孔隙水壓力分布圖見圖6-8。

2.2.3 飽和與飽和-非飽和穩定滲流場對比分析

對比分析圖3和圖6，大壩下游的浸潤線高度在非飽和穩定滲流場中和飽和時相比偏低。分析原因在飽和穩定滲流分析時，理論認為浸潤線位置處無水能的損失，僅為滲流場的分界線；在非飽和穩定滲流場分析影響時，水量從上向下由飽和區進入非飽和區，原因是非飽和區的基質吸引力造成的，這樣使得非飽和滲流的浸潤線高度降低，壓力水頭下降，由此帶來水頭損失的增加。

圖6 飽和-非飽和穩定滲流場總水頭等值線分布圖

圖7 飽和-非飽和穩定滲流場壓力水頭等值線分布圖

圖8 飽和-非飽和穩定滲流場孔隙水壓力等值線分布圖

對比分析圖4和圖5，在飽和穩定滲流場的分析中，以浸潤線為分界線分割飽和區，若僅認為浸潤線以下的飽和區發生滲流，忽略浸潤線以上的飽和區域，忽略對部分壓力總水頭和孔隙水壓力的計算分析，這樣的分析是不合理的。對比分析了非飽和區在滲流場中的作用，浸潤線是飽和-非飽和滲流區的分界線。其下孔隙水壓力、總水頭值均是正值，反之均為負值。

通過對比分析飽和-非飽和穩定滲流場的分析得出壓力總水頭、孔隙水壓的等值線分布規律大致相同，在實際工程中對二者的分析研究必須結合起來。為了使更加顯著區別，對非穩定滲流的課題研究尤為重要。

3 結 論

文章主要是在土石壩飽和-非飽和土體的理論研究基礎上，根據有限單元法的基本原理和數學模型的建立，運用有限元程序模擬了正常工況下，飽和-非飽和穩定滲流的浸潤線變化趨勢，對滲流場中壓力總水頭、孔隙水壓力進行了計算分析，結合工程實際認為計算結果正確合宜，通過等值線分布圖驗證了壩坡的安全穩定性。同時認為滲流分析是飽和滲流與非飽和滲流有本質的區別，二者的結合分析會更貼近工程實際。

[1]顧慰慈．滲流計算原理及應用[M]．北京：中國建材工業出版社，2000：26-37.

[2]劉杰．土石壩滲流控制理論基礎及工程經驗教訓[M]．北京：水利水電出版社，2005：4-12.

[3]鄭敏生，錢鏡林，蘇玉杰．考慮非飽和區的土石壩滲流分析[J]．水力發電學報，2010，29(01)：186-191.

[4]毛昶熙．滲流計算分析與控制[M]．北京：中國水利水電出版社，2003：1-63，89-124，306-342.

[5]王華?。\屏二級水電站閘基深厚覆蓋層滲流分析與控制研究[D]．四川：成都理工大學環境與土木工程學院，2005．

[6]柴軍瑞．大壩工程滲流力學[M]．拉薩：西藏人民出版社，2001：6-7.

[7]朱曉源．考慮非飽和土體的土石壩滲流與壩坡穩定分析研究[D]．杭州：浙江大學，2006.

[8]黃俊，蘇向明．土壩飽和-非飽和滲流數值分析方法研究[J]．巖土工程學報，1990(05) .

Impact Analysis for Saturation and Un-saturation Stability Seepage Field of Earth and Rock Fill Dam

HAN Xue

(Tonghe County Tax Technology Service Center, Tonghe 150900, China)

Earth-rock dam is a common dam type at the present time, and the construction of seepage protection system takes very important roles for the stability of dam body, several safety accidents in water conservancy projects both at home and abroad will arouse enough attention. To analyze the seepage of earth and rock dam needs to be defined as a community of saturation and un-saturation areas, the past experiences have proved them to be supplemented each other. For both the movement characteristics, un-saturation seepage is more complex and the affecting factors are more than saturation seepage. In this paper, combination with the project on-site, the valuable simulation for limited element is applied to draw the changing law of seepage fields in both saturation and un-saturation seepage, concluding a correct theory about seepage filed, and the results studied will have a certain of reference values in constructing the seepage protection system of earth and rock fill dam.

earth-rock dam;un-saturated soil;stable seepage field

1007-7596(2017)05-0023-04

2017-04-18

韓雪(1968-)，女，黑龍江通河人，工程師。

TV641

B