基于優(yōu)勢(shì)流效應(yīng)的填埋場(chǎng)滲瀝液回灌過程預(yù)測(cè)

梁 冰， 張 柴， 劉 磊

(1.遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院， 阜新 123000； 2.中國科學(xué)院武漢巖土力學(xué)研究所， 巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室， 武漢 430071；3.中國科學(xué)院武漢巖土力學(xué)研究所-香港理工大學(xué)固體廢棄物科學(xué)聯(lián)合實(shí)驗(yàn)室， 武漢 430071)

滲濾液是經(jīng)降雨入滲和地下水補(bǔ)給在垃圾填埋場(chǎng)內(nèi)部流動(dòng)形成的高濃度、有毒有害特質(zhì)液體[1]，生物反應(yīng)器填埋場(chǎng)將經(jīng)導(dǎo)排系統(tǒng)收集后的滲濾液循環(huán)回灌到填埋場(chǎng)中，加強(qiáng)了垃圾、微生物和水分間的相互接觸和作用，有助于促進(jìn)微生物降解速率、減少滲濾液的處理成本和對(duì)環(huán)境污染的危害，但回灌效果受垃圾自身的水力特性、回灌井布水方式、回灌速率、回灌頻率和回灌量等多方面因素影響[2]。部分學(xué)者將垃圾視為均勻、各向同性的多孔介質(zhì)，用單孔隙度模型描述水分在垃圾堆體中的遷移規(guī)律[3-5]，但廣泛原位試驗(yàn)和工程應(yīng)用中表明模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)具有較大誤差，這歸因于各地區(qū)填埋垃圾壓實(shí)方式、組成成分各不相同，復(fù)雜的顆粒形狀、粒徑大小、排列結(jié)構(gòu)引起滲流通道的不均勻，以大孔隙優(yōu)勢(shì)流為主導(dǎo)[6]。因此，考慮垃圾的優(yōu)勢(shì)流效應(yīng)去設(shè)計(jì)生物反應(yīng)器填埋場(chǎng)滲濾液回灌工藝參數(shù)更具合理性、準(zhǔn)確性和科學(xué)性。

中外學(xué)者對(duì)垃圾中的優(yōu)勢(shì)流和滲濾液回灌運(yùn)移規(guī)律展開研究。優(yōu)勢(shì)流效應(yīng)方面，Gerke 等[7]基于對(duì)流-彌散方程建立了描述裂隙流和基質(zhì)流的雙滲透率模型，并考慮兩域間水分質(zhì)量交換；Kode?ová等[8]設(shè)計(jì)了優(yōu)先流路徑的控制入滲和排水實(shí)驗(yàn)，預(yù)測(cè)了非均質(zhì)土體的水流分布和溶質(zhì)運(yùn)移；Han等[9]、Audebert等[10]采用中尺度實(shí)驗(yàn)裝置研究了滲濾液入滲或排出的遷移規(guī)律，表明水分分布演化表現(xiàn)出多域流動(dòng)行為。滲濾液回灌方面，White等[11]假設(shè)垃圾為均質(zhì)和各向同性，數(shù)值模擬和現(xiàn)場(chǎng)滲濾液回灌試驗(yàn)了水平溝或垂直井的響應(yīng)性能；Feng等[12]考慮含水率、揮發(fā)性脂肪酸(volatile fatty acid,VFA)濃度和溫度對(duì)垃圾降解影響的條件下填埋場(chǎng)的穩(wěn)定化進(jìn)程的最佳回灌參數(shù)。以上研究為本文模擬的構(gòu)建提供了借鑒。

為了進(jìn)一步分析典型回灌工況對(duì)堆體水分分布及遷移過程的影響，基于垃圾堆體內(nèi)部的優(yōu)勢(shì)流效應(yīng)，運(yùn)用Hydrus-1D軟件模擬了不同回灌工況下滲濾液的遷移與分布，對(duì)比分析了不同回灌速率、回灌頻率、回灌量和不同初始含水率對(duì)填埋垃圾水分影響深度、流速、域間質(zhì)量交換、貯水率的影響規(guī)律。

1 數(shù)學(xué)模型及求解

1.1 數(shù)學(xué)模型

1.1.1 單孔隙度模型

由連續(xù)性原理和達(dá)西定律推導(dǎo)的一維狀態(tài)非飽和土壤Richards方程[13]可表述為

(1)

式(1)中：θ為含水率；t為時(shí)間；K(h)為滲透系數(shù)；z為位置；h為水壓頭；S為源匯項(xiàng)。

土水特征曲線(水壓頭h與含水率θ)、非飽和滲透系數(shù)K(θ)的本構(gòu)關(guān)系是求解方程(1)的關(guān)鍵，相關(guān)學(xué)者提出了一系列描述土水特征關(guān)系、非飽和導(dǎo)水率的經(jīng)驗(yàn)數(shù)學(xué)模型，其中Van Genuchten (VGM)模型[14]在土壤和垃圾中應(yīng)用最為廣泛，可寫為

(2)

(3)

式中:Se為有效飽和度；α為進(jìn)氣壓力值的倒數(shù)；n為孔徑分布參數(shù)；m為土水特征曲線對(duì)稱性參數(shù)；Ks為飽和導(dǎo)水率；l為扭曲因子。

1.1.2 雙滲透率模型

雙滲透率模型(double permeability model，DPeM)假定多孔介質(zhì)域由兩個(gè)子域組成，包括一個(gè)孔隙尺寸較大的裂隙域，以及一個(gè)孔隙尺寸小的基質(zhì)域，在大孔裂隙域中水分流動(dòng)被假定為垂直的、層流的，只由重力驅(qū)動(dòng)，沒有毛細(xì)效應(yīng)，兩域內(nèi)的水分均可流動(dòng)且允許兩域之間存在質(zhì)量交換，由一維非線性Richards方程定義裂隙域和基質(zhì)域控制方程分別為[7]

(4)

式(4)中：θf、θm分別為裂隙域、基質(zhì)域水分體積除以總體積；Kf、Km分別為裂隙域、基質(zhì)域?qū)?；hf、hm分別為裂隙域、基質(zhì)域水頭；Sf、Sm為源匯項(xiàng)；wf為裂隙域體積除以總流動(dòng)域體積，0

(5)

式(5)中：β為材料形狀決定的幾何因子；a為有效擴(kuò)散路徑長度，即從基質(zhì)中心到裂隙邊界的距離；γw為經(jīng)驗(yàn)尺度因子；Ka為裂隙域與基質(zhì)域交界面有效滲透系數(shù)。

1.2 參數(shù)設(shè)置

雙滲透模型對(duì)于描述垃圾中的滲濾液優(yōu)勢(shì)流動(dòng)具有很強(qiáng)的適用性，但參數(shù)復(fù)雜，難以獲取，研究相對(duì)較少。Audebert等[10]假設(shè)垃圾是由裂縫網(wǎng)絡(luò)組成的非均質(zhì)介質(zhì)，由電阻率層析成像(electrical resistivity tomography，ERT)反演法系統(tǒng)性探討了雙滲透率模型的各個(gè)參數(shù)取值，具有較高的準(zhǔn)確性，因此選用Audebert具有一定工程背景下的雙滲透率模型相關(guān)參數(shù)，如表1所示。表1中，θr為殘余含水率，θs為飽和含水率。

1.3 數(shù)值模擬

模擬填埋深度為30 m的生物反應(yīng)器填埋場(chǎng)，以點(diǎn)源形式進(jìn)行滲濾液回灌，上、下邊界分別設(shè)定為通量邊界、滲透面邊界，初始條件為含水率。根據(jù)國內(nèi)填埋場(chǎng)基本物理特性的調(diào)查研究顯示，中國垃圾初始含水率在21%～48%，設(shè)定初始含水率為30%和35%，并考慮垃圾成層性分布的影響[15]。馮世進(jìn)等[15]、李守升等[16]研究了中外生物反應(yīng)器填埋場(chǎng)典型回灌工況分析表面噴灑回灌的滲流情況，建議回灌強(qiáng)度在0.1～0.2 m/d；劉丹[17]模擬填埋場(chǎng)導(dǎo)排層堵塞時(shí)分析滲濾液回灌強(qiáng)度為0.28～0.57 m/d適宜；陳馨等[18]、楊強(qiáng)[19]指出厭氧-準(zhǔn)好氧聯(lián)合型生物反應(yīng)器填埋場(chǎng)不同產(chǎn)酸階段回灌頻率在1～4 d/次?；诖耍欣诶盥駡?chǎng)穩(wěn)定和滲濾液的利用率的最佳回灌速率在0.1～0.5 m/d，配水次數(shù)為1～4 d/次，模擬方案如表2所示。

表1 數(shù)值模擬的參數(shù)取值[10]Table 1 Parameters of numerical simulation[10]

表2 滲濾液回灌模擬方案

2 優(yōu)勢(shì)流效應(yīng)下滲濾液回灌遷移規(guī)律

2.1 回灌速率對(duì)滲濾液遷移的影響

回灌速率可控制滲濾液入滲過程與填埋垃圾既有非飽和孔隙接觸的時(shí)間，以低速率回灌時(shí)水分流動(dòng)慢且影響深度低，不能滿足各位置垃圾層微生物降解所需的最佳水分；而以高速率回灌時(shí)會(huì)減少滲濾液在填埋場(chǎng)內(nèi)部的滯留時(shí)間，還可能引起垃圾堆體中微生物被沖刷。

圖1所示為不同回灌速率下雙滲透率模型(DPeM)和單孔隙度模型(single porosity model，SPM)預(yù)測(cè)的含水率與埋深曲線，增大回灌速率促進(jìn)了滲濾液向底部延伸，以第10天為例，回灌速率從0.15 m/d增大0.35 m/d時(shí)，DPeM預(yù)測(cè)滲濾液影響深度從埋深6 m延伸至14 m，這是由于垃圾堆體的非飽和滲透性與含水率正相關(guān)，回灌速率的增大使?jié)B透性增強(qiáng)，進(jìn)而加速了干燥區(qū)的濕潤，但回灌結(jié)束時(shí)速率0.35 m/d卻比0.15 m/d的影響深度小5 m，這是由于控制同等回灌量下，低回灌速率比高回灌速率要長17 d，時(shí)間的加劇效應(yīng)較回灌速率占主導(dǎo)，與馮世進(jìn)等[15]的研究相符合。單域SPM相同條件下影響深度、流速僅占DPeM的64.3%、53.5%，要小于Feng等[12]SPM預(yù)測(cè)影響深度、流速為DPeM的80%，滲濾液達(dá)到底部DPeM比SPM平均快124 d。

圖1 不同回灌速率影響深度曲線Fig.1 Influence depth curve with different injection rate

圖2所示為回灌速率為0.25 m/d的水分質(zhì)量交換，初期滲濾液先從滲透性較高的裂隙域進(jìn)入形成優(yōu)勢(shì)流，在第15天、埋深15 m處達(dá)到0.017 d-1，由于持水性差，流向持水性高的基質(zhì)域，裂隙域含水量下降，基質(zhì)域含水量升高，基質(zhì)域水分達(dá)到一定限度且裂隙域水分向深部流動(dòng)時(shí)，基質(zhì)域的水分再次回流到裂隙域，但此時(shí)因?yàn)榛|(zhì)域吸持了一定的水分，只有一部分質(zhì)量交換項(xiàng)回流到裂隙域，質(zhì)量轉(zhuǎn)移項(xiàng)即隨深度而減小。

圖3所示為不同回灌速率單位面積上、下邊界累積流量隨時(shí)間變化曲線，入滲曲線先線性增大到回灌量后趨于不變，回灌速率越大，斜線愈陡。DPeM預(yù)測(cè)滲濾液幾乎全是由裂隙域流入，而下邊界近96%滲濾液從基質(zhì)域中流出，大部分滲濾液在垃圾堆體內(nèi)部發(fā)生了質(zhì)量交換。當(dāng)回灌速率增大時(shí)，不同的回灌速率單位面積累積流出量幾乎保持一致，意味著貯水率相等。將貯水率定義滲濾液流入與流出的差值占回灌總量的比值，如圖4所示，SPM預(yù)測(cè)一年模擬結(jié)束時(shí)平均流出1.35 m3/m2，貯水率均達(dá)70%以上，比DPeM預(yù)測(cè)貯水率要高出56%，極大高估了垃圾堆體對(duì)滲濾液的吸持性。

圖2 回灌速率為0.25 m/d質(zhì)量交換曲線Fig.2 Mass transfer curve with injection rate 0.25 m/d

圖3 不同回灌速率單位面積流量曲線Fig.3 Flow per unit area curve with different injection rate

圖4 不同回灌速率貯水率曲線Fig.4 Water storage rate curve with different injection rate

2.2 回灌頻率對(duì)滲濾液遷移的影響

回灌頻率即滲濾液的配水次數(shù)，當(dāng)配水次數(shù)過低時(shí)，滲濾液的影響深度只局限于填埋場(chǎng)上部垃圾的降解，不能緩慢釋流延伸到底部；而當(dāng)頻率過高時(shí)，反復(fù)收集、導(dǎo)排和回灌的運(yùn)行成本會(huì)明顯增多。

圖5所示為不同回灌頻率下DPeM、SPM預(yù)測(cè)的含水率與埋深曲線，滲濾液的影響深度隨著回灌頻率的減小而逐漸縮短，仍以第10天為例，回灌頻率從1 d/次減小至4 d/次時(shí)，DPeM預(yù)測(cè)滲濾液的影響深度由埋深7 m減小至2 m，配水次數(shù)n=1、2、4回灌結(jié)束時(shí)滲濾液的影響深度為埋深20、25、30 m，與回灌頻率呈負(fù)相關(guān)，配水頻次的降低雖減少了入滲量、濕潤面積，流速變慢，但增加了滲濾液在垃圾堆體內(nèi)部滲流的時(shí)間，使得影響深度更大。同等回灌頻率下滲濾液到達(dá)底部DPeM預(yù)測(cè)的時(shí)間要比SPM快127 d，因?yàn)樵诨毓嗲捌趦?yōu)勢(shì)流模型流速是均質(zhì)性模型流速的1.67倍，SPM流速已漸漸趨近與0，但DPeM預(yù)測(cè)底部仍以0.1 m/d的最大流速進(jìn)行釋水。

圖5 不同回灌頻率影響深度曲線Fig.5 Influence depth curve with different injection frequency

圖6所示為回灌頻率為2 d/次的水分質(zhì)量交換曲線，在第4天、埋深1 m處達(dá)到0.045 d-1，同時(shí)發(fā)生從基質(zhì)域到裂隙域的質(zhì)量交換，這是由于當(dāng)回灌具有一定時(shí)間間隔后，裂隙域快速釋水，優(yōu)勢(shì)入滲，含水量降低，儲(chǔ)存在基質(zhì)域中的水分向裂隙域中回流，在第61天具有最大的質(zhì)量交換項(xiàng)0.08 d-1，而到第16天～第61天，埋深10～20 m處，發(fā)現(xiàn)質(zhì)量交換幾乎為0，推測(cè)堆體內(nèi)部形成穩(wěn)定滲流。

圖7所示為不同回灌頻率單位面積上、下邊界累積流量隨時(shí)間變化曲線，滲濾液入滲隨著回灌頻率的不同呈階梯化上升后不變。減小回灌頻率后單位面積的累積流出量也逐漸減小，合適的回灌頻率可以延長回灌過程中滲濾液的滯留時(shí)間，進(jìn)而為增加垃圾堆體既有的非飽和孔隙接觸吸持、微生物降解消耗提供可能。當(dāng)回灌頻率減小，貯水率明顯減小的更慢，如圖8所示，以1、2、4 d/次頻率回灌時(shí)的DPeM貯水率在60%以上的時(shí)間分別是前74、103、161 d，模擬結(jié)束時(shí)仍有16%的貯水率，而均質(zhì)化的SPM預(yù)測(cè)幾乎每種回灌頻率下滲濾液貯水率都能達(dá)到60%。

圖6 回灌頻率為2 d/次質(zhì)量交換曲線Fig.6 Mass transfer curve with injection frequency 2 d/次

圖7 不同回灌頻率單位面積流量曲線Fig.7 Flow per unit area curve with different injection frequency

圖8 不同回灌頻率貯水率曲線Fig.8 Water storage rate curve with different injection frequency

2.3 回灌量對(duì)滲濾液遷移的影響

回灌量即從導(dǎo)排系統(tǒng)中收集的滲濾液重新回灌到垃圾堆體中的總量，回灌量過低時(shí)不能滿足微生物降解所需要的水分和養(yǎng)分，而高量的再循環(huán)滲濾液可能會(huì)引發(fā)環(huán)境巖土災(zāi)害及給滲濾液收集系統(tǒng)帶來沉重的工作負(fù)荷。

圖9所示為不同回灌量DPeM預(yù)測(cè)含水率與埋深曲線，當(dāng)單位面積回灌量3 m3/m2增大1倍時(shí)，DPeM預(yù)測(cè)影響深度由埋深13 m延伸至27 m，3種回灌量下滲濾液到達(dá)底部的時(shí)間分別為60、41、41 d，回灌量的增加對(duì)于滲濾液向底部延伸具有一定促進(jìn)作用，但單位面積回灌量為4.5 m3/m2和6 m3/m2時(shí)水分到達(dá)底部的時(shí)間相同，這是由于垃圾堆體的最大滲流能力是恒定的，高量回灌滲濾液并不會(huì)產(chǎn)生良好的回灌效果。SPM預(yù)測(cè)回灌量3、4.5、6 m3/m2水分到達(dá)底部時(shí)間較DPeM結(jié)果分別晚305、124、54 d，孔隙均質(zhì)化后流速要小得多，優(yōu)勢(shì)流模型流速是均質(zhì)化模型的1.2倍，靠近下邊界DPeM有突然釋水、流速產(chǎn)生突變，而SPM此時(shí)水分流速已減小至0～0.3 m/d，與埋深擬線性分布。

圖9 不同回灌量影響深度曲線Fig.9 Influence depth curve with differenttotal injection

圖10所示為回灌量為6 m3/m2的水分質(zhì)量交換曲線，回灌初期基質(zhì)域孔隙開始吸持大量滲濾液，在第8天、第17天時(shí)水分的質(zhì)量交換最大，埋深5～15 m處以0.015 d-1質(zhì)量交換項(xiàng)從基質(zhì)域向裂隙域回流?；毓噙^程在45th、埋深25 m處再次產(chǎn)生較大的質(zhì)量交換，這是由于下邊界出水量小于受重力作用向底部流動(dòng)的滲濾液聚集量，裂隙域中水分被吸持，貯存到基質(zhì)域中。

圖10 回灌量6 m3/m2質(zhì)量交換曲線Fig.10 Mass transfer curve with total injection 6 m3/m2

圖11 不同回灌量單位面積流量曲線Fig.11 Flow per unit area curve with different total injection

圖12 不同回灌量貯水率曲線Fig.12 Water storage rate curve with different total injection

圖11所示為不同回灌量單位面積上、下邊界累積流量隨時(shí)間的變化，單位面積的入滲曲線均保持相等斜率線性增大后保持不變。隨著回灌量的增大，DPeM中的裂隙域和基質(zhì)域中的流出量也隨之增大，當(dāng)回灌量從3 m3/m2增大1倍，DPeM從基質(zhì)域中流出的滲濾液從2.2 m3/m2增大到5.14 m3/m2，垃圾凈吸持量僅增加5%，回灌量4.5、6 m3/m2模擬結(jié)束時(shí)最終貯水率均為14%。SPM回灌量為3、4.5 m3/m2時(shí)貯水率均達(dá)到60%以上，6 m3/m2時(shí)達(dá)到60%貯水率為前300 d，如圖12所示，表明提高回灌量可以有效縮短回灌周期，但回灌量過大只會(huì)增加滲濾液收集導(dǎo)排系統(tǒng)的負(fù)擔(dān)。

2.4 初始含水率對(duì)滲濾液遷移的影響

垃圾的初始含水高，滲濾液回灌時(shí)水分遷移速度快，吸持、貯存的滲濾液少；含水率低，非飽和區(qū)域大，滲流速度慢，所能吸持的回灌滲濾液量也就較大，在設(shè)計(jì)滲濾液回灌工藝時(shí)必須考慮垃圾自身初始含水率。

圖13所示為不同初始含水率下DPeM和SPM預(yù)測(cè)含水率與埋深變化曲線，仍以第10天為例，初始含水率從0.3增加至0.35時(shí)，DPeM預(yù)測(cè)滲濾液的影響深度分別由埋深4 m延伸至6 m，SPM流速與之相比要小0.02 m/d，在滲濾液到達(dá)深部時(shí)，DPeM的流速與垃圾層深度擬線性分布，而SPM流速已趨近于0。不同深度、不同時(shí)間SPM預(yù)測(cè)的垃圾最大含水率始終比DPeM要大，由于DPeM中預(yù)測(cè)含水率為裂隙域和基質(zhì)域含水率的加權(quán)平均值，而裂隙域沒有持水性，只有基質(zhì)域持水，含水率經(jīng)過加權(quán)平均后會(huì)減小，這也符合垃圾堆體內(nèi)部滲濾液優(yōu)先大孔隙通道流動(dòng)但不貯存的遷移特征。

圖13 不同初始含水率影響深度曲線Fig.13 Influence depth curve withdifferent initial moisture content

圖14 初始含水率0.3～0.35質(zhì)量交換曲線Fig.14 Mass transfer curve with initial moisture content 0.3～0.35

圖14所示為考慮垃圾分層初始含水率0.3～0.35的水分質(zhì)量交換?？梢钥闯觯诘?0天、第31天時(shí)裂隙域和基質(zhì)域之間水分的質(zhì)量交換最大，達(dá)到0.24 d-1，由于埋深0～15 m的含水率低(θ0=0.3)，滲流速度慢，回灌的滲濾液從裂隙域向基質(zhì)域中流動(dòng)。在第31天～第50天期間天埋深5～15 m質(zhì)量交換項(xiàng)幾乎為0，第50天后埋深20～30 m有頻繁的水分質(zhì)量交換，這是由于下邊界的不斷流出，上部裂隙域和基質(zhì)域水分隨流出量的多少而動(dòng)態(tài)轉(zhuǎn)化。

圖15所示為不同初始含水率單位面積上、下邊界累積流量隨時(shí)間的變化曲線，單位面積累積入滲流量時(shí)間曲線重合。隨著初始含水量的增大，單位面積的累積流出量也隨之增大，初始含水率增大0.05時(shí)DPeM預(yù)測(cè)下邊界多流出0.678 m3/m2，這是由于垃圾的持水性是恒定的，所能吸持滲濾液也有限。在考慮垃圾成層性分布特征、初始含水率θ0介于0.3～0.35時(shí)，DPeM和SPM單位面積流出曲線趨勢(shì)介于θ0為0.3、0.35的均值，初始含水率愈小，貯水率隨時(shí)間減小變慢，如圖16所示。

圖15 不同初始含水率單位面積流量曲線Fig.15 Flow per unit area curve with different initial moisture content

圖16 不同初始含水率貯水率曲線Fig.16 Water storage rate curve with different initial moisture content

3 結(jié)論

(1)θ0=0.3時(shí)最佳回灌工藝為v=0.25 m/d、n=2 d/次、q=4.5 m3/m2，一定限度內(nèi)增大回灌速率、回灌頻率、回灌量、初始含水量等條件可提高影響深度、入深量、貯水率。

(2)回灌時(shí)幾乎所有滲濾液從裂隙域入滲到垃圾堆體內(nèi)部，與均質(zhì)化的VGM模型入滲曲線具有高度重合性，但在下邊界疏水過程近95%滲濾液均是從基質(zhì)域中流出，只有5%的滲濾液從裂隙域流出。

(3)水分質(zhì)量交換項(xiàng)在回灌初期由裂隙域流向基質(zhì)域，且從裂隙域流向基質(zhì)域的水分始終要比從基質(zhì)域回流至裂隙域的水分大，隨深度而減小。

(4)雙滲透率模型具有更高的水力滲透性和流速，滲濾液可以更快地通過裂隙入滲、流動(dòng)和排出，SPM模型低估了垃圾填埋場(chǎng)的滲流能力，不能準(zhǔn)確地描述液體在生物反應(yīng)器填埋場(chǎng)中的運(yùn)動(dòng)和滯留。

(5)優(yōu)勢(shì)流效應(yīng)是垃圾填埋場(chǎng)滲濾液遷移的重要運(yùn)動(dòng)形式，在填埋場(chǎng)的實(shí)際工程中須基于不同垃圾堆體的水力學(xué)特性，再綜合考慮回灌工藝因素才可實(shí)現(xiàn)最佳回灌效果。