趙 義，施建勇

（河海大學教育部巖土力學與堤壩工程重點實驗室，江蘇南京210098）

群井回灌滲濾液條件下填埋場水分遷移規(guī)律研究

趙 義，施建勇

（河海大學教育部巖土力學與堤壩工程重點實驗室，江蘇南京210098）

為研究群井回灌滲濾液條件下填埋場內部水分遷移規(guī)律，基于已有的考慮生化降解的城市固體廢棄物（MSW）力－液－氣耦合一維沉降模型，并考慮到垃圾的分層填埋造成滲濾液的水平向滲透系數(shù)大于豎向滲透系數(shù)，忽略滲濾液的豎向流動，建立了填埋場群井回灌滲濾液時的液相運移方程，推求出常水頭和常流量補水條件下的解析解，并以常水頭回灌滲濾液為例分析了填埋場內水壓變化規(guī)律。結果表明：填埋場內部生化反應造成水壓初期降低，之后垃圾降解速率減慢和滲濾液回灌充分造成水壓增大，且回灌水頭越大、井距越小，水壓增大越快。

豎井；滲濾液回灌；解析解；水壓

1 引 言

回灌滲濾液是現(xiàn)代衛(wèi)生填埋場加速填埋穩(wěn)定的常見措施，回灌方式分為豎井回灌，水平溝回灌，表面噴灌等方式。國內外學者針對豎井回灌滲濾液條件下填埋場中水分遷移規(guī)律進行了一些有意義的研究。

Jain［1］假定垃圾體是均質且不可壓縮的介質，不考慮氣相、生化反應等因素影響，利用SEEP／W程序計算分析填埋場用豎井回灌滲濾液時井的尺寸、垃圾體的水力特性、回灌水壓等對回灌影響區(qū)域、回灌穩(wěn)定所需時間的影響，并通過參數(shù)的無量綱化繪制計算圖表以方便進行回灌設計。Khire［2］在不考慮溫度、生化反應等條件影響并假定垃圾體均質和各向同性的情況下，利用HYDRUS－2D分析了單井回灌時，回灌速率、垃圾體滲透系數(shù)、井徑等因素對浸濕寬度和襯墊上部壓力水頭的影響。Al－Thani［3］學者在不考慮溫度變化和填埋場生化反應等條件并假定滲透系數(shù)隨深度變化、豎井貫穿填埋體的情況下，利用MODELFOLW－SUREACT計算分析了群井條件下浸潤面對抽排滲濾液效果的影響。王洪濤［4］針對填埋場滲濾液回灌建立了三維飽和－非飽和非穩(wěn)定水分遷移數(shù)值模型，該模型考慮到垃圾體的非均質和各向異性和填埋介質在降解和壓力作用下的沉降，通過該模型作者給出不同回灌條件下的計算方法，并比較了豎井回灌和水平溝回灌效果?？梢钥吹?，群井下滲濾液回灌的研究較少，且較多研究未能考慮填埋場內部生化反應的影響。本文將著重研究群井回灌滲濾液下填埋場內部水分的遷移。

2 計算模型的建立和求解

填埋場中垃圾的分層填埋會造成滲濾液水平向滲透性大于豎向滲透性［5］，同時，考慮到填埋場中氣體消散相對較快，因此本文忽略了填埋場中滲濾液的豎向流動和氣相影響。建立控制方程前引入假設條件：①豎井貫穿填埋體，滲濾液只發(fā)生水平向遷移；②垃圾體的液相滲透系數(shù)為常量；③填埋場上覆壓力不變，無空氣和雨水入滲。劉曉東［6］曾提出產氣峰值模型，根據(jù)非飽和土力學理論，建立了填埋場力－氣－液耦合一維沉降模型，這里，本文只關心滲濾液的運移規(guī)律，因此只需將其沉降模型中的液相運移方程變換成平面軸對稱形式，即得控制方程：。各符號意義如下：uw為液壓（Pa）；kw為液相滲透系數(shù)（m／d）；V0為垃圾土初始總體積（m3）；t為時間（d）；k是MSW中纖維素和半纖維素的比例；Mb0為初始可降解有機物干重（kg）；ρw為水的密度；A和B是和溫度T有關的降解參數(shù)。為相應于基質吸力變化的液體體積變化系數(shù)（Pa－1）。

填埋場進行滲濾液回灌時有常水頭和常流量兩種補水方式，分別計算如下：

2.1 常水頭回灌

填埋場進行群井下常水頭回灌，邊界條件和初始條件如下：

式（2）中rw為井徑，re為影響半徑。將邊界條件齊次化，令

uw（r，t）＝u（r，t）＋P，則 u（r，t）滿足如下條件：

為求解上述定解問題，采用特征函數(shù)法先求解對應的齊次方程，設u＝T（t）R（r），代入齊次方程得到特征函數(shù)系

R（r）＝Y0（βnrw）J0（βnr）－J0（βnrw）Y0（βnr），β由Y0（βnrw）J1（βnre）－J0（βnrw）Y2（βnre）確定。

將 u＝T（t）［Y0（βnrw）J0（βnr）－J0（βnrw）Y0（βnr）］

代入方程并結合初始條件，解得

2.2 常流量回灌

填埋場進行群井下常流量回灌，邊界條件和初始條件如下：

式中：q2為回灌流量（m3／d），L為井深（m），其余參數(shù)意義同上。

變換得：

考慮到虛宗量貝塞爾函數(shù)在自變量為實數(shù)時發(fā)散，因此需確定ε為0和虛數(shù)時方程解答，據(jù)文獻［9］，得ˉh的Laplace逆變換為：

則

其中βn由

J1（βnrw）Y1（βnre）－Y1（βnrw）J1（βnre）＝0確定。

3 計算分析

3.1 對比驗證

考慮到文中公式較復雜，為保證計算結果的準確性，現(xiàn)對比驗證如下：

（1）對于常水頭回灌，這里將方程（1）中的耗水項簡化成常數(shù)，得到和文獻［9］類似的方程，將兩者進行對比，以驗證本文結果的可靠性。文獻［9］的定解問題如下：

易解得

其中An和βn意義同式（4）。取一潛水井實例計算［9］，參數(shù)如下：初始水位 H0＝20 m；進井水位 h0＝14 m；rw＝1 m；re＝200 m；水位傳導系數(shù) B1＝850 m2／d；入滲強度M ＝0.04m／d。

將計算結果和文獻［9］進行對比，如圖1所示。

圖1 t＝10 d時水頭曲線對比圖

圖1為t＝10 d時的水頭對比圖，可以看出，計算結果很接近，存在差異的原因是計算所取項數(shù)不同，表明公式（4）具有可靠性。

（2）對于常流量回灌，本文相比文獻［9］只是多了一個耗水項，不考慮耗水項，定解問題如下：

令式（6）中D＝0，即得式（9）的定解：

取計算參數(shù)：rw＝0.305 m，re＝2000 m，H0＝85m，q0＝1.145457，B1＝6×106m2／d。

這里仍取t＝10 d的計算結果進行對比，如圖2所示。

圖2 水頭曲線對比圖

可以看到，結果接近，存在差異的原因在于文獻［9］計算未計入級數(shù)項，表明公式（6）是可靠的。

3.2 滲濾液遷移規(guī)律

以常水頭回灌為例，各計算參數(shù)取值為［6－7］：溫度T＝293.6 K；降解物質干量Mb0＝28.6 kg；體積V0＝0.1m3；纖維素和半纖維素比例K＝4；降解參數(shù)A＝0.03，B＝11.13；體積變化參數(shù)＝2× 10－6Pa－1；初始水壓q＝7989 Pa；液相滲透系數(shù)Kw＝0.01m／d；rw＝0.5m；re＝25m；回灌水壓P＝20000 Pa。計算結果見圖3和圖4。

圖3和圖4表示的是超靜孔隙水壓隨時間和半徑的變化關系?？梢钥闯鏊畨弘S時間先降低再增大，且隨半徑增大而減小，這種變化規(guī)律是由填埋場的生化反應決定的，垃圾體的降解是耗水產氣的過程，初始時段降解耗水多，回灌的滲濾液尚未影響到計算半徑處，造成水壓降低。隨著時間的增大，垃圾土降解速率降低，而回灌的滲濾液補給充分，水壓逐漸增大。下面，再分別分析回灌水頭和井距對回灌水壓的影響。

圖3 孔隙水壓隨半徑和時間的變化曲線圖

圖4 不同半徑處孔隙水壓隨時間變化圖

3.2.1 回灌水頭對水壓的影響

取水頭分別為2m、4 m、6 m，其他參數(shù)同上，取r＝15 m的計算結果對比，見圖5。

圖5 水頭高度對水壓的影響圖

由圖5可以看出，回灌水頭高度越大，同一時刻水壓越大，且增大得更快。

3.2.2 井距對水壓的影響

取井距為30m、50 m、80 m，其他參數(shù)同上，取r＝15m的計算結果對比，見圖6。

由圖6可以看出，井距越大，同一時刻水壓越小，且水壓增大越慢。通過以上分析，填埋回灌滲濾液時，應選擇合理的井距和水頭，水頭越大，井距越小，回灌效果越顯著，但填埋場滲濾液水壓的增加會影響其穩(wěn)定性，因此生化填埋場回灌滲濾液時應選擇合適的井距和水頭，以免造成填埋場失穩(wěn)。

圖6 井距對水壓的影響圖

4 結論及展望

基于劉曉東提出的城市固體廢棄物（MSW）力－液－氣耦合一維沉降模型，在忽略液相豎向遷移的基礎上，得到填埋場回灌滲濾液條件下液壓控制方程，得到常水頭和常流量補水條件下的解答，計算結果表明：

（1）填埋場內部的生化反應造成水壓初期先減小，之后回灌滲濾液補給充足且垃圾降解速率減慢，水壓逐漸增大，同時水壓隨半徑增大而減??；

（2）回灌常水頭越大，井距越小，水壓增大越快，但水壓增大過快不利于填埋場的穩(wěn)定性，因此填埋場回灌滲濾液要選擇合適的水頭和井距。

（3）本文在建立群井回灌滲濾液條件下的液相控制方程時，做了較多簡化，實際上填埋場回灌滲濾液是三維飽和－非飽和滲流問題，水分沿豎向也會有遷移，外部降雨也會滲入填埋場，同時滲透系數(shù)會隨著滲濾液的回灌發(fā)生變化，而填埋場的降解及沉降也會影響回灌效果，因此群井回灌滲濾液時水壓的變化規(guī)律還有待深入研究。

［1］ Jain P，Townsend TG，Tolaymat TM.Steady-state design of verticalwells for liquids addition at bioreactor landfills［J］.Waste Management，2010，30（11）：2022-2029.

［2］ Khire M V，Mukherjee M.Leachate injection using vertical wells in bioreactor landfills［J］.Waste Management，2007，27（9）：1233-1247.

［3］ Al-Thani A A，Beaven R P，White JK.Modeling flow to leachatewells in landfills［J］.Waste Management，2004，24（3）：271-276.

［4］ 王洪濤，殷 勇.滲濾液回灌條件下生化反應器填埋場水分運移數(shù)值模擬［J］.環(huán)境科學，2003，24（2）：66-72.

［5］ 錢學德，施建勇，劉曉東，等.現(xiàn)代衛(wèi)生填埋場的設計與施工（第二版）［M］.北京：中國建筑工業(yè)出版社，2011.

［6］ 劉曉東.城市固體廢棄物（MSW）考慮生化降解的力－液－氣耦合沉降模型及試驗研究［D］.南京：河海大學，2012.

［7］ Fredlund D G，Rahardjo H.Soilmechanics for unsaturated soils［M］.New York：JohnWiley and Sons，1993.

［8］ Rowe K K，Nadarajah P.Estimating leachate drawdown due to pumping wells in landfills［J］.Journal of Canadian Geotechnical，1996，33（1）：1-10.

［9］ 李佩成.地下水非穩(wěn)定滲流解析法［M］.北京：科學出版社，1990.

Analysis of Moisture Movement in Landfill Site under Leachate Recirculation w ith VerticalWells

ZHAO Yi，SHIJian-yong
（Key Laboratory ofMinistry of Education for Geomechanics and Embankment Engineering，HohaiUniversity，Nanjing，Jiangsu 210098，China）

To study themoisturemigration of the landfill site under the condition of leachate recirculation with vertical wells，the flow equation of leachate is established based on one-dimensional settlementequation by considering hydraulicgaseous coupling effect and degradation ofmunicipal solid waste（MSW）aswell as ignoring the vertical flow of leachate when the leachate permeability in horizontal direction is greater than that in vertical direction for layered landfills.Then，the analytical solution under constanthead and constant flow recirculation is got，and the transport law of leachate is studied with the constanthead recirculation as an example.The results show that the liquid pressure is decreased firstby biochemical reactions，and then increased due to the leachate recirculation and the decrease ofwaste degradation rate.Simultaneously，the liquid pressure is increased faster with the increase of recirculation head and decrease ofwell space.

verticalwell；leachate recirculation；analytical solution；liquid pressure

X705

A

1672—1144（2014）01—0188—05

10.3969／j.issn.1672－1144.2014.01.039

2013-08-18

2013-09-16

趙 義（1990—），男，江蘇南京人，碩士研究生，研究方向為環(huán)境巖土工程。