降雨作用下填埋場防滲系統水分分布的數值模擬研究*

劉 勇，田 宇，周小軍

（1．襄陽市固體廢棄物處理公司，湖北 襄陽 441057；2．武漢市環境衛生科學研究院，湖北 武漢 430015；3．中國科學院武漢巖土力學研究所，湖北 武漢 430071）

降雨作用下填埋場防滲系統水分分布的數值模擬研究*

劉 勇1，田 宇2，周小軍3

（1．襄陽市固體廢棄物處理公司，湖北 襄陽 441057；2．武漢市環境衛生科學研究院，湖北 武漢 430015；3．中國科學院武漢巖土力學研究所，湖北 武漢 430071）

基于填埋場防滲系統的工程特性，分析了填埋場的降雨入滲機理，建立了填埋場防滲系統水分動態數學模型，通過數值分析研究了降雨作用下填埋場防滲系統孔隙水壓力和含水率的變化情況。結果表明：降雨入滲使填埋場結構層滲流條件發生變化，在坡頂、坡面、坡腳處均產生滲透力，且隨著降雨時間的持續，負孔隙水壓力逐漸減小；在水平分布上，孔隙水壓力曲線呈較大的波動性，且在坡腳處出現峰值。受降雨持續作用，填埋場各結構層含水率逐漸增加，導致材料的密度逐漸增大；結構層材料因雨水的入滲逐漸達到其儲水能力而飽和，土體由非飽和向飽和過程轉化。

填埋場；防滲系統；降雨；水分分布；數值模擬

填埋場襯墊系統失穩將引起填埋場滲瀝液潛在性泄漏，污染周圍水土環境，造成嚴重的環境地質災害，已經引起了國際環境巖土工程界的高度關注［1］。填埋場復合襯墊防滲系統是填埋場工程結構中的關鍵部分，是填埋場底部、邊坡或封場覆蓋中用于隔離滲瀝液、填埋氣體等物質的一種重要屏障，同時也是填埋場長期、有效運營的重要保證。CJJ 17—2004生活垃圾衛生填埋技術規范與GB 16889—2008生活垃圾填埋場污染控制標準推薦采用復合襯墊防滲系統。填埋場的襯墊系統應具有良好的防滲功能，其滲透系數必須小于1×10-7cm/s［2］。采用數值模擬分析方法分析襯墊系統的安全穩定特性及其變形趨勢預測，可定量分析和評價在不同環境條件下其長期有效性，為合理設計填埋場襯墊系統與正確評價和預測填埋場襯墊防滲系統的長期穩定性提供理論依據。

1 填埋場降雨入滲機理分析

填埋場產生不均勻沉降變形，很大一部分是雨水下滲對垃圾土層產生作用，逐步改變垃圾土層的力學性質，使填埋場產生沉降變形［3］。填埋場降雨入滲經歷了飽和-非飽和的滲流過程。

地表水在下滲過程中對垃圾土層的作用機理有2種［4］：①地表水下滲對垃圾土層的軟化壓縮作用。當地表水滲入到垃圾土層，垃圾土受到水的浸泡，其小顆粒骨架受到破壞，產生一種軟化壓縮變形作用，在上覆不同重力作用下，其壓縮變形大小不一從而使得填埋場出現不均勻沉降變形。②地表水下滲對土層產生機鑄化學潛蝕作用。地表水滲入垃圾土層后，將土層內一部分細微土粒及可溶解于水的物質發生混濁作用和化學溶解作用，然后順滲漏途徑流失［5-6］。當這種作用反復多次后，逐步使土層內流失一部分物質，土層勢必產生收縮下沉。

2 填埋場防滲系統水分分布的數學模型

2.1 滲流控制方程

Buckingham［7］于1907年在Darcy方程的基礎上經過改進提出非飽和多孔介質滲流方程：

式中：ρ0為純水密度；S為儲水系數；為滲透系數張量；ρ*為流體密度；q為多孔介質內部的源流或匯流；t為時間。

由公式（2）可知造成水流流動的原因有：位置水頭坡度（▽Z），壓力水頭坡度（▽ψ），其他流體的影響 （ρ0/ρ）。

若不考慮多孔介質與流體的壓縮性與流體密度變化，且滲流場中無任何源流或匯流，則可將公式 （2） 簡化成公式 （3）：

2.2 防滲系統幾何模型

取垃圾填埋場封場覆蓋坡體的典型剖面，有限元計算模型及網格剖分如圖1所示，自下而上，分別為填埋層、防滲層、覆蓋層、植被層，填埋場邊坡坡比約為3∶1，剖分單元918。

圖1 填埋場封場覆蓋幾何模型及網格剖分

模型中假設左邊界和右邊界分別為第1類邊界條件，定義初始水頭分別為1、3 m，上邊界為第2類邊界條件，假設降雨強度為暴雨100 mm/d，下邊界為隔水邊界。基本參數取值見表1。持續降雨作用下填埋場變形見圖2，降雨作用下填埋場總水頭變化云圖見圖3，降雨作用下填埋場在x向和y向的變化云圖如圖4、5所示。

表1 基本參數值

圖2 持續降雨作用下填埋場變形

圖3 降雨作用下填埋場總水頭變化云圖

3 降雨作用下孔隙水壓力及含水率的變化規律

圖6、7分別為無降雨和持續降雨作用下填埋場表層孔隙水壓力隨時間和空間的變化。在無降雨入滲時，孔隙水壓力隨時間增加而增加，在水平距離變化時，孔隙水壓力呈波動變化，在坡角和坡頂均呈現最大峰值和最小峰值。在降雨入滲時，雨水通過填埋場邊坡的坡頂、坡面及坡腳向各結構層入滲。隨降雨時間增加，填埋場表層孔隙水壓力呈上升趨勢，第1天降雨為負的孔隙水壓力，隨時間增加，各結構層土體的逐漸飽和而轉變成正的孔隙水壓力；隨水平距離增加，地下水滲流過程對各結構層土顆粒施加壓力，孔隙水壓力在坡角處形成拐點，達到峰值；孔隙水壓力從-5 kPa上升到坡角處的0 kPa，隨水平距離的繼續增加，峰值后的孔隙水壓力呈下降趨勢。

圖4 降雨作用下填埋場x向應力云圖

圖5 降雨作用下填埋場y向應力云圖

圖6 無降雨入滲植被土層的孔隙水壓力

圖7 降雨條件下植被土層的孔隙水壓力

圖8 、9分別顯示的是降雨入滲條件下不同結構層孔隙水壓力隨時間和空間的變化曲線。降雨入滲導致滲流條件的變化，隨著時間的增加孔隙水壓力上升，隨著水平距離的增加孔隙水壓力呈波動變化，對于防滲層表層孔隙水壓力基本均在0 kPa以上，這是因為防滲層滲透率較低，滲透水逐漸進入，因膨脹而增大材料中孔隙，結構層材料慢慢達到飽和狀態，孔隙水壓力迅速增大。孔隙水壓力曲線的變化也反映了土體結構的非飽和-飽和的轉化過程。

圖8 降雨作用覆蓋層的孔隙水壓力曲線

圖9 降雨作用防滲層的孔隙水壓力曲線

圖10 為填埋場底部孔隙水壓力變化曲線。由于防滲系統的作用、填埋深度的影響，垃圾填埋體產生的滲瀝液均聚集在填埋場底部，故均為正的孔隙水壓力。圖11顯示了填埋場表層含水率分布，受降雨入滲的影響，隨時間的增加含水率逐漸增加，含水率的增加而使結構層材料的密度增加，負孔隙水壓力減小，從曲線看出降雨過程中填埋場表層含水率較之無降雨過程有較大的上升；在降雨持續6 d作用后，結構層材料的密度增加直至達到飽和密度，土體因達到其儲水能力而飽和，含水率保持定值。

圖10 降雨作用填埋層的孔隙水壓力曲線

圖11 填埋場植被土層含水率分布

4 結論

1）降雨入滲使填埋場結構層滲流條件發生變化，坡頂、坡面、坡角均產生滲透力。隨降雨的持續，負孔隙水壓力逐漸減小；隨水平距離的增加，孔隙水壓力曲線呈較大的波動性，這是由于坡角處滲透量的聚集，導致孔隙水壓力的驟增而使坡角處出現峰值，因此坡角處要加強排水防護措施。

2）受降雨持續作用，填埋場各結構層含水率逐漸增加，導致材料的密度逐漸增大；結構層材料因雨水的入滲逐漸達到其儲水能力而飽和，土體由非飽和向飽和過程轉化。

3）填埋場封場覆蓋防滲系統變形規律的研究表明，隨著水頭的降低，流體壓力逐漸降低，垃圾層、黏土層和土壤層的孔隙壓力下降，固體骨架有效應力增大，從而導致了填埋場防滲系統位移的增加。

因此，防滲系統防滲層必須適應不均勻沉降，同時，在填埋場設計時有必要考慮滲流作用和地震影響，以提高填埋場穩定性和安全性。

［1］胡敏云，陳云敏．城市生活垃圾填埋場沉降分析與計算［J］．土木工程學報，2001，34（6）：89-92．

［2］張振營，陳云敏．城市垃圾填埋場沉降模型的研究［J］．浙江大學學報，2004，38（9）：1162-1165．

［3］戚國慶．降雨誘發滑坡機理及其評價方法研究［D］．成都：成都理工大學，2004．

［4］李廣信．高等土力學［M］．北京：清華大學出版社，2004．

［5］ Xue Q，Feng X T，Liang B．Slippage Solution of Gas Pressure Distribution in the Process of Landfill Gas Seepage［J］．Appl Math Mech，2005，26（12）：1623-1633．

［6］薛強，趙穎，劉磊，等．垃圾填埋場災變過程的溫度-滲流-應力-化學耦合效應研究［J］．巖石力學與工程學報，2011，30（10）：1971-1988．

［7］薛強，劉磊，梁冰，等．垃圾填埋場沉降變形條件下氣-水-固耦合動力學模型研究［J］．巖石力學與工程學報，2007，26（S1）：3473-3476．

［8］ Van Genuchten M Th．A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils［J］．Soil Sci Soc Am J，1980，44 （5）：892-898．

Numerical Simulation of Moisture Distribution for Anti-seepage System of Waste Landfill Site by Rainfall

Liu Yong1，Tian Yu2，Zhou Xiaojun3
（1．Xiangyang Solid Waste Treatment Company，Xiangyang Hubei441057；2．Wuhan Environmental Sanitation Science Research Institute，Wuhan Hubei430015；3．Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，Wuhan Hubei430071）

Based on the engineering characteristics of anti-seepage system in waste landfill sites，the rainfall infiltration mechanism of the landfill site was analyzed，and the water dynamic mathematical model of the anti-seepage system was set up．Through the numerical analysis，the changes of pore water pressure and moisture content of the anti-seepage system were researched under the condition of rainfall．The results showed that the rainfall infiltration could change the seepage condition of the structure layer，and there is the seepage force on the top，the slope surface and the slope toe．The negative pore water pressure gradually decreased with the rainfall continuing．The pore water pressure curve waved largely in the horizontal distribution，and appeared the peak in the slope toe．The moisture content in each structure layer gradually increased with the rainfall continuing，and it could increase the density of the material．The material of the structure layer gradually reached its storage capacity to the saturation because of the rainfall infiltration，and the soil changed from unsaturation to saturation．

landfill site；anti-seepage system；rainfall；moisture distribution；numerical simulation

X703；TU993．3

A

1005-8206（2012）04-0021-04

國家自然科學基金（51079143）；武漢市科技攻關項目（201060723312）

2012-07-02

劉勇（1972—），高級工程師，主要從事環境治理。

（責任編輯：鄭雯）