垃圾填埋場防滲屏障的服役壽命評價(jià)

席永慧 楊 帆 蔡策毅

(同濟(jì)大學(xué)土木工程學(xué)院建筑工程系,上海 200092)

0 引 言

隨著經(jīng)濟(jì)的快速發(fā)展,人民生活水平提高,生活垃圾的產(chǎn)生量也不斷增加。填埋是生活垃圾處理的主要方式之一,國家統(tǒng)計(jì)局最新數(shù)據(jù)顯示,2014年全國垃圾填埋量達(dá)10 744.3萬噸[1]。為了避免填埋場中的滲濾液滲出,導(dǎo)致地下水、土壤或混凝土遭受污染,在已存在或者正在建設(shè)中的填埋場地應(yīng)用污染物防滲屏障技術(shù),其目的是確保結(jié)構(gòu)的完整性,并減少因滲濾液遷移造成的危害。

填埋是生活垃圾處理的主要方法之一。垃圾填埋場滲濾液滲漏引起的地下環(huán)境污染問題越來越受到社會(huì)的關(guān)注和重視。垂直防滲屏障是最主要污染控制措施,已獲得廣泛應(yīng)用。防滲屏障在填埋場中有廣泛的應(yīng)用,用于生活垃圾填埋場滲濾液污染控制的垂直防滲屏障的滲透系數(shù)宜在10-7cm/s量級[2],其類型可選用水泥-膨潤土墻、土-膨潤土墻、塑性混凝土墻、HDPE土工膜-膨潤土復(fù)合墻等。

目前國內(nèi)外關(guān)于防滲屏障的研究主要集中于屏障材料特性和配比對防滲性能影響的試驗(yàn)研究[3-6];陳云敏[7]提出了基于污染物濃度擊穿時(shí)間的屏障設(shè)計(jì)方法,屏障擊穿時(shí)間需要大于垃圾穩(wěn)定固化時(shí)間和填埋場運(yùn)行時(shí)間之和;劉偉[8]對懸掛式及嵌入式塑性混凝土防滲帷幕服役壽命進(jìn)行了研究,通過建立填埋場含塑性混凝土防滲帷幕的二維運(yùn)移數(shù)值分析模型,分析了不同因素對塑性混凝土材料類型的防滲帷幕防污性能的影響;詹良通等[9]研究了垃圾填埋場污染物擊穿豎向防滲帷幕時(shí)間的影響因素,以COD作為代表性污染物,在蘇州七子山填埋場地質(zhì)條件勘查基礎(chǔ)上建立了有限差分計(jì)算模型,對防滲帷幕被污染物擊穿時(shí)間的影響因素進(jìn)行了計(jì)算分析;D.G.Ruffing[10]在膨潤土中參入粉煤灰制成的泥漿,使鎘離子在防滲屏障中的擴(kuò)散系數(shù)減小60%,對屏障的擊穿時(shí)間增加了5.3倍;也有對于其他防滲基材的使用壽命進(jìn)行研究的,Abdelaal[11]對不同滲濾液情況下,土工膜的機(jī)械性能的研究,并用阿倫尼烏斯模型對土工膜的破壞時(shí)間進(jìn)行推斷;國內(nèi)外對防滲屏障的服役壽命的研究仍較少。

對于我國填埋場來說,滲濾液水位往往較高,達(dá)到了十幾米甚至更高。在如此高的水位條件下,防滲屏障的服役壽命尤其值得關(guān)注。防滲屏障服役壽命研究對于無控制填埋場的污染狀況及風(fēng)險(xiǎn)的評估有重要的作用,而且可以對防滲屏障的設(shè)計(jì)提出要求和建議。本文基于垃圾填埋場遷移方程的解析解,進(jìn)行了Matlab曲線的模擬。采用本項(xiàng)目支持的實(shí)驗(yàn)室擴(kuò)散試驗(yàn)得到的試驗(yàn)數(shù)據(jù),結(jié)合浦東老港垃圾填埋場對其服役壽命進(jìn)行評價(jià)與分析,并結(jié)合某工業(yè)搬遷場地對其防滲屏障進(jìn)行設(shè)計(jì)。

1 一維對流-彌散遷移方程及邊界條件

人們關(guān)心的是填埋場中污染物離子的遷移對地下水及周圍土壤的危害。描述污染物離子在土壤中的遷移,通常采用一維數(shù)學(xué)表達(dá)式,這類偏微分方程統(tǒng)稱為遷移方程。污染場地的滲濾液中含有多種污染物質(zhì),這些污染物質(zhì)在黏土、水泥土(水泥摻量較小)襯墊或防滲屏障中遷移時(shí),將受到對流運(yùn)移、分子擴(kuò)散、機(jī)械彌散、吸附作用等遷移機(jī)理的綜合作用。綜合考慮這些遷移機(jī)理的作用,滲濾液中污染物在防滲屏障中遷移的計(jì)算模型包含如下假定:

(1) 土壤多孔介質(zhì)系統(tǒng)是均質(zhì)的、半無限的;

(2) 土壤是飽和的,地下水在土壤中的流動(dòng)符合達(dá)西定律;

(3) 溶質(zhì)沿一個(gè)方向遷移;

(4) 防滲屏障中滲透系數(shù)k是常數(shù);

(5) 污染物在防滲屏障中的遷移只發(fā)生在一個(gè)方向。

在上述假定的基礎(chǔ)上,污染物在填埋場防滲屏障中的遷移就可用一維對流-彌散模型來描述平流、彌散與吸附等作用對滲濾液中污染物在防滲屏障中遷移的影響。如果這些假定被接受,且既不考慮污染物的源匯,又不考慮污染物的蛻變,則污染物溶質(zhì)通過飽和土壤的一維遷移模型由著名的平流-彌散方程[12-13]來表達(dá):

(1)

式中:Rf是阻滯因子;D為污染物在土中的分子擴(kuò)散系數(shù);v是滲流速度。

假定:①土壤多孔介質(zhì)系統(tǒng)是半無限的;②土體開始沒有受污染;③污染源濃度是常數(shù),則邊界條件為

c(∞,t)=0

c(x,0)=0

c(0,t)=c0

(2)

式中:c0為常數(shù),土/滲濾液交界面處的溶質(zhì)濃度。

Freez和Cherry(1979)給出了遷移方程在上述邊界條件下的解析解:

(3)

對于反應(yīng)性物質(zhì)(Rf≠1):

;

(4)

式中:v*=v/Rf;De=D/Rf(De為有效擴(kuò)散系數(shù))。擴(kuò)散模型示意圖如圖1所示。

圖1 擴(kuò)散模型Fig.1 Model of diffusion

這個(gè)解被廣泛地應(yīng)用在野外和實(shí)驗(yàn)室數(shù)據(jù)分析中。本文用MATLAB進(jìn)行計(jì)算,繪制了不同De值、v值組合下,對有限厚度(L)的防滲屏障,給定c,c0,vs,v,De值,在屏障出口處污染液的濃度可將L代入方程中的x計(jì)算求得,這樣可得到屏障出口處相對濃度與時(shí)間的關(guān)系圖c/c0～t/Rf曲線圖2、圖3。如果給定t,則可求得屏障在經(jīng)過給定時(shí)間t時(shí),屏障出口處的相對濃度,這樣可畫出c/c0-x曲線圖(圖4和圖5)。

圖2 De=0.001 4 m2/a,v=0.003 m/a情況下c/c0-t/Rf曲線圖Fig.2 Curves of c/c0-t/Rf under the condition of De=0.001 4 m2/a,v=0.003 m/a

2 填埋場服役性能的評價(jià)

《生活垃圾衛(wèi)生填埋場巖土工程技術(shù)規(guī)范》(CJJ 176—2012)[14]中提到防滲屏障的服役壽命應(yīng)該大于填埋場運(yùn)行時(shí)間和固廢穩(wěn)定化時(shí)間的時(shí)間總和,即

圖3 De=0.000 47 m2·a-1,v=0.003 m·a-1情況下c/c0-t/Rf曲線圖Fig.3 Curves of c/c0-t/Rf under the condition of De=0.000 47 m2/a,v=0.003 m/a

圖4 De=0.001 4 m2·a-1,v=0.003 m·a-1情況下c/c0-x曲線圖Fig.4 Curves of c/c0-x under the condition of De=0.001 4 m2/a,v=0.003 m/a

圖5 De=0.000 47 m2·a-1,v=0.003 m·a-1情況下c/c0-x曲線圖Fig.5 Curves of c/c0-x under the condition of De=0.000 47 m2/a,v=0.003 m/a

(5)

填埋場的運(yùn)行時(shí)間,可用防滲屏障的擊穿時(shí)間代表。確定防滲屏障的擊穿時(shí)間可按照下列步驟進(jìn)行:

(1) 選擇設(shè)計(jì)標(biāo)準(zhǔn):首先,選定設(shè)計(jì)厚度;其次,確定在設(shè)施設(shè)計(jì)壽命期中土屏障中最大的滲流速度。

(2) 確定或估計(jì)設(shè)計(jì)中滲濾液的初始濃度c0。

(3) 計(jì)算擊穿水平ce/c0;關(guān)于擊穿濃度標(biāo)準(zhǔn)沒有統(tǒng)一說法,有按照擊穿濃度相對值確定[15],有按照地下水標(biāo)準(zhǔn)確定,本研究用ce/c0=10%來確定擊穿濃度ce,當(dāng)溶液原濃度c0較大時(shí),根據(jù)《國外水和空氣質(zhì)量標(biāo)準(zhǔn)》,對不同污染物離子有不同的擊穿濃度限制。

(4) 估計(jì)有效擴(kuò)散系數(shù),用已報(bào)道的研究成果或通過擴(kuò)散試驗(yàn)確定。

(5) 確定阻滯因子Rf?！渡罾l(wèi)生填埋場巖土工程技術(shù)規(guī)范》[14](CJJ 176—2012)中提到在重金屬離子污染的土體中,阻滯因子可取3～40,計(jì)算時(shí),保守起見取阻滯因子為3。 因?yàn)檠芯勘砻?阻滯因子的影響因素很多,如混和液和單種溶液的有效擴(kuò)散系數(shù)和阻滯因子是有差別的,相伴離子對擴(kuò)散的影響也是明顯的。

(6) 利用相關(guān)曲線圖估計(jì)屏障的擊穿時(shí)間。利用所求得的c/c0-t/Rf曲線,在Matlab程序中求出擊穿水平為10%的點(diǎn),從而求得污染物離子擊穿防滲屏障的時(shí)間。

3 案例分析

3.1 案例一

現(xiàn)有上海老港衛(wèi)生填埋場,填埋深度3.5 m。為防止垃圾滲瀝液擴(kuò)散污染土壤與地下水,在垃圾體底部設(shè)有防滲襯墊,垃圾體周邊設(shè)置豎向防滲墻,墻深12 m,厚300 mm,墻體采用混合材料,膨潤土∶水泥∶砂=2.57∶1∶1.2,水固比取1.4。已知垃圾滲瀝中銅離子濃度為2 mg/L,不考慮防滲屏障兩側(cè)的水頭差。

通過我們最新的室內(nèi)擴(kuò)散試驗(yàn)結(jié)果測得銅離子在上述混合材料的有效擴(kuò)散系數(shù)D/Rf=4.5×10-7cm2/s (0.001 4 m2/a),滲濾液滲流速度為0.003 m2/a,查圖2的局部放大圖(圖6),擊穿水平為10%時(shí),在Matlab程序中找到相對應(yīng)的點(diǎn),可知t/Rf=10 a (年),因?yàn)樽铚蜃訛?,所以銅離子屏障擊穿的時(shí)間t=30 a (年)。

圖6 De=0.001 4 m2/a,v=0.003 m/a情況下c/c0～t/Rf曲線選點(diǎn)圖Fig.6 Curves of c/c0～t/Rf under the condition of De=0.001 4 m2/a,v=0.003 m/a

王羅春等[16]對老港垃圾填埋場穩(wěn)定化的研究表明:老港填埋場達(dá)到三、二、一級穩(wěn)定化狀態(tài)所需的時(shí)間分別為4 a、10 a、32 a。根據(jù)李華和趙由才[17]對上海市老港填埋場的研究,垃圾填埋24 a后,垃圾中可生物降解含量(BDM)為2.55%,填埋場表面沉降量0.009 m2/a,滲濾液COD濃度約為5 mg/L,即可認(rèn)為達(dá)到穩(wěn)定化。據(jù)此本文假定,老港填埋場的穩(wěn)定化時(shí)間為24 a。

根據(jù)公式(5)和銅離子的有效擴(kuò)散系數(shù),老港填埋場的服役壽命為屏障擊穿時(shí)間與老港填埋場的穩(wěn)定化時(shí)間總和,即為54 a (年){30 a (屏障擊穿時(shí)間)+24 a (填埋場穩(wěn)定化時(shí)間)},滿足填埋場的一般設(shè)計(jì)使用年限50年的要求。

3.2 案例二

現(xiàn)有某工業(yè)搬遷場地,經(jīng)調(diào)查場地土壤地下水存在重金屬等污染,其中地下水銅離子濃度500 mg/L,最大污染深度4 m。擬對污染區(qū)域采用水泥土攪拌樁墻進(jìn)行污染阻隔,水泥摻量為12%,需設(shè)計(jì)此污染場地防滲屏障的厚度。

根據(jù)我們最新的室內(nèi)擴(kuò)散試驗(yàn)結(jié)果,銅離子在水泥摻量為12%的防滲屏障中的有效擴(kuò)散系數(shù)Rf=1.5×10-7cm2/s (0.000 47 m2/a)為參考,滲濾液滲流速度取0.003 m2/a。根據(jù)擬合的曲線(圖7),在Matlab中選擇擊穿水平c/c0為10%的點(diǎn),同樣取穩(wěn)定化時(shí)間為24年,根據(jù)曲線圖假定t=35 a。從圖中查得所需的防滲屏障最小厚度x=0.202 m,結(jié)合結(jié)構(gòu)使用的長久性和經(jīng)濟(jì)性,防滲屏障厚度宜設(shè)計(jì)為0.25 m。

圖7 De=0.000 47 m2/a,v=0.003 m/a情況下c/c0-x曲線選點(diǎn)圖Fig.7 Curves of c/c0-x under the condition of De=0.000 47 m2/a,v=0.003 m/a

4 結(jié) 論

從一維對流一彌散方程的數(shù)值模擬結(jié)果可以看出:

(1) 隔離屏障的服役壽命大小主要取決于隔離屏障的擊穿時(shí)間,而對隔離屏障取的De值越大,擊穿時(shí)間越短,De提高3倍可將擊穿時(shí)間增加1.5倍。所以如果要增加隔離屏障的服役壽命,則可以通過改變防滲屏障材料來降低金屬離子在防滲屏障中的De值。

(2) 填埋場隔離屏障服役壽命是填埋場穩(wěn)定化時(shí)間與隔離屏障擊穿時(shí)間之和,通過對上海老港填埋場0.30 m厚防滲屏障的分析,查閱c/c0～t/Rf曲線圖得出上海老港填埋場的服役壽命為54 a。

(3) 填埋場隔離屏障服役壽命是填埋場穩(wěn)定化時(shí)間與隔離屏障擊穿時(shí)間之和,通過對屏障擊穿時(shí)間的假定,對案例二中的某工業(yè)場地的防滲屏障設(shè)計(jì),假定其擊穿時(shí)間為35年,通過查閱c/c0～x的曲線圖,厚度宜設(shè)計(jì)為0.25 m。