向 銳， 徐 亞， 劉景財(cái)， 劉玉強(qiáng)*， 董 路， 雷國元， 黃啟飛

1.中國環(huán)境科學(xué)研究院， 環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室， 北京 100012 2.中國環(huán)境科學(xué)研究院固體廢物污染控制技術(shù)研究所， 北京 100012 3.武漢科技大學(xué)， 湖北 武漢 430081

工業(yè)固體廢物尤其是危險(xiǎn)廢物的環(huán)境管理和無害化處置是全球性環(huán)境問題[1]. 世界范圍內(nèi)，尤其在發(fā)展中國家填埋依然是危險(xiǎn)廢物無害化處置的主要手段[2-4]. 為了降低危險(xiǎn)廢物的泄露風(fēng)險(xiǎn)和不利影響，使得危險(xiǎn)廢物填埋場(hazardous waste landfill, HWL)的設(shè)計(jì)過于繁雜. 然而，機(jī)械損壞及化學(xué)老化等因素均可能引起土工膜性能退化，最終導(dǎo)致工程襯層系統(tǒng)失效[5-6]，所以填埋場泄漏及地下水環(huán)境污染事件仍時有發(fā)生[7]. 因此，填埋場工程屏障的優(yōu)化設(shè)計(jì)[8-9]和地質(zhì)屏障的選擇至關(guān)重要[10].

合理的工程屏障，如導(dǎo)排系統(tǒng)和防滲系統(tǒng)，可以減少滲濾液滲漏量[11-12]. 有利的地質(zhì)屏障，如厚且滲透系數(shù)小的包氣帶可以增加污染物的穿透時間，使其充分降解；更新能力強(qiáng)的含水層可以快速稀釋污染物[13]，降低污染物濃度. 因此，提高工程或地質(zhì)屏障要求可以有效地降低地下水污染風(fēng)險(xiǎn)，但同時也需要考慮工程費(fèi)用[14]. 因?yàn)楫?dāng)工程和地質(zhì)屏障升至一定程度時，工程費(fèi)用與污染控制效益比將顯著降低. 因此，識別工程和地質(zhì)屏障的關(guān)鍵因子，并確定其合理取值顯得尤為重要.

為控制填埋場滲漏風(fēng)險(xiǎn)，許多國家都對填埋場的工程和地質(zhì)屏障制定了強(qiáng)制性標(biāo)準(zhǔn). 如美國[15]和中國[16-17]均規(guī)定導(dǎo)排層坡度≥0.02，導(dǎo)排層滲透系數(shù)≥1.0×10-1cms；德國[18]和中國[16-17]均規(guī)定庫底與地下水位之間的距離≥1 m. 但是這些取值背后的科學(xué)背景并不明確，同時是否適用于我國實(shí)際情況，也未經(jīng)系統(tǒng)研究和論證. 一些學(xué)者通過基于數(shù)學(xué)模型的過程模擬方法[19]和指標(biāo)體系方法[20]研究了工程和地質(zhì)屏障對其環(huán)境污染風(fēng)險(xiǎn)的影響，如Kotowska等[21]通過內(nèi)梅羅指數(shù)和化學(xué)計(jì)量學(xué)方法識別了關(guān)鍵的參數(shù)，并給出了各參數(shù)的權(quán)重. 但是大部分權(quán)重都是基于主觀方法給出，背后的物理含義并不明確，因此合理性有待商榷. 另外，基于指標(biāo)體系的方法也無法體現(xiàn)同一參數(shù)在不同取值區(qū)間對填埋場滲漏風(fēng)險(xiǎn)的影響.

與指標(biāo)體系方法相比，過程模擬方法考慮了整體影響因素，可以給出綜合評估結(jié)果，并在其他領(lǐng)域被廣泛應(yīng)用于關(guān)鍵參數(shù)及敏感區(qū)間的識別[22-23]. 然而，基于該方法開展的危險(xiǎn)廢物填埋場關(guān)鍵參數(shù)和敏感區(qū)間識別，以及優(yōu)化設(shè)計(jì)參數(shù)的研究鮮見報(bào)道. 鑒于此，該研究擬采用過程模擬方法，結(jié)合參數(shù)敏感性分析，開展工程和地質(zhì)屏障參數(shù)優(yōu)化研究，識別不同參數(shù)的敏感性和敏感區(qū)間，推薦最適宜取值，評估各參數(shù)的相對重要性排名，以期為國內(nèi)危險(xiǎn)廢物填埋場工程和地質(zhì)屏障參數(shù)的選擇及相關(guān)標(biāo)準(zhǔn)和規(guī)范的制定提供理論依據(jù).

1 模型與方法

1.1 代表性場景構(gòu)建

為從理論上分析不同工程和地質(zhì)屏障參數(shù)對地下水污染風(fēng)險(xiǎn)的影響，識別其重要程度和敏感性程度，需要建立代表性的填埋場，模擬滲濾液滲漏與地下水污染的場景. 填埋場結(jié)構(gòu)根據(jù)GB 18598—2001《危險(xiǎn)廢物填埋場污染控制標(biāo)準(zhǔn)》的要求確定，從上到下劃分為四部分，即封場覆蓋系統(tǒng)、危險(xiǎn)廢物層、主滲濾液防滲和導(dǎo)排系統(tǒng)、次滲濾液防滲和導(dǎo)排系統(tǒng)(見表1).

表1 代表性的填埋場剖面結(jié)構(gòu)及尺寸

圖1 危險(xiǎn)廢物填埋場滲濾液產(chǎn)生、滲漏及地下水污染過程Fig.1 Generation, leakage and groundwater pollution of leachate in hazardous waste landfill

危險(xiǎn)廢物填埋場滲濾液的產(chǎn)生、滲漏及地下水污染過程如圖1所示. 由于雨水淋溶及廢物自身重力釋水所產(chǎn)生的滲濾液，在重力作用下垂直下滲進(jìn)入到填埋場導(dǎo)排層內(nèi)，一部分經(jīng)過導(dǎo)排系統(tǒng)收集并處理達(dá)標(biāo)后排放；另一部分以滲漏(通過防滲膜漏洞)或滲透(通過土工膜完整部分以蒸汽擴(kuò)散形式泄露)形式進(jìn)入到包氣帶中. 經(jīng)過包氣帶的截污和含水層的稀釋擴(kuò)散后進(jìn)入到周邊水井中，通過飲水、洗澡等途徑暴露于人體，對人體健康構(gòu)成危害. 對填埋場的地下水污染風(fēng)險(xiǎn)進(jìn)行評估實(shí)際上就是對防滲層破損條件下的滲濾液滲漏以及滲漏導(dǎo)致的地下水污染程度進(jìn)行評估.

1.2 基于過程模擬的滲濾液滲漏及滲漏污染預(yù)測

1.2.1滲濾液滲漏預(yù)測模型

通過設(shè)計(jì)雨水防滲系統(tǒng)、滲濾液導(dǎo)排和防滲系統(tǒng)，盡量減少危險(xiǎn)廢物填埋場滲濾液的產(chǎn)生和滲漏. 然而，滲漏事故仍然頻繁發(fā)生. 主要原因是：作為防滲材料核心的HDPE膜(high-density polyethylene geomembrane, 高密度聚乙烯膜)在安裝和施工過程中不可避免地會造成一些破損；滲濾液通過HDPE膜破損的滲漏受襯層結(jié)構(gòu)(雙襯層、單襯層等)、漏洞大小、漏洞上方的滲濾液飽和液位、HDPE膜下方介質(zhì)的滲透系數(shù)及其與HDPE膜接觸情況等諸多因素影響，過程極為復(fù)雜. 為簡化滲濾液滲漏計(jì)算，Giroud等[24]通過室內(nèi)滲漏模擬實(shí)驗(yàn)和仿真模型分析，建立了不同條件下的滲漏速率計(jì)算公式. 針對我國危險(xiǎn)廢物填埋場的主要防滲結(jié)構(gòu)〔兩層HDPE膜和CCL (compacted clay liner, 壓實(shí)黏土襯層)〕，根據(jù)Giroud等[25]開發(fā)的經(jīng)驗(yàn)?zāi)Ｐ汀惨娛?1)〕來估計(jì)滲漏速率：

(1)

式中：Q為滲漏速率，m3/s；a為HDPE中的破損面積，m2；ks為CCL的水力傳導(dǎo)率，m/s；Ls為CCL的厚度，m；hw為HDPE上的滲濾液高度，m；N為HDPE中的破洞密度，holes/(104m2)(見表2)；S為垃圾填埋場的底部面積，104m2；βc為HDPE和CCL之間的接觸系數(shù)，良好和不良時分別為0.21和1.15.

表2 HDPE膜常規(guī)破損情況

1.2.2滲濾液濃度的衰變模擬

危險(xiǎn)廢物填埋場中的滲濾液是危險(xiǎn)廢物堆放和填埋過程中，由于壓實(shí)和雨水滲流作用產(chǎn)生的一種高濃度的有機(jī)或無機(jī)成份的液體. 顯然，隨著雨水不斷淋溶，危險(xiǎn)廢物中的有害組分產(chǎn)生不同耗損，滲濾液中的污染物濃度也逐漸降低，這一過程通常用源項(xiàng)衰減模型〔見式(2)(3)〕[24-25]描述：

Ct=C0-λt

(2)

(3)

式中：Ct為t時滲濾液中有害物質(zhì)的濃度，mg/L；C0為滲濾液中有害物質(zhì)的初始濃度，mg/L；λ為滲濾液中有害物質(zhì)濃度變化系數(shù)，10-7m-3a-1;t為時間，a；i為滲透速率，mm/a；Wd為最終廢物的深度，m；Wfc為廢物的現(xiàn)場容量，104m3.

1.2.3滲濾液及其危害組分在土壤-地下水中遷移轉(zhuǎn)化預(yù)測

污染物在包氣帶和地下水介質(zhì)中的遷移和分布主要受水流的稀釋效應(yīng)和衰減效應(yīng)(即降解和吸附)的影響. 在均勻和各向同性的土壤-水系統(tǒng)中，上述效應(yīng)可以通過一維對流-擴(kuò)散方程〔見式(4)(5)〕[29]描述.

(4)

DL=av+Dm

(5)

式中：x為沿流動方向的路徑距離，m；c為距離為x、時間為t時污染物的濃度，mg/L；v為地下水流速，m/s；n為有效孔隙度；R為衰減因子；γ為一階衰減率，s-1；DL為流體動力學(xué)縱向彌散系數(shù)，m2/s (見表3)；a為介質(zhì)的分散度，m；Dm為分子擴(kuò)散系數(shù)，m2/s.

表3 含水層參數(shù)取值[6]

顯然，對于不同的邊界條件，上述對流-擴(kuò)散方程的解析解形式不同. 對于滲漏這一邊界條件，假設(shè)污染物濃度按式(2)所示的規(guī)律衰減，則可通過拉普拉斯變換獲得其解析解[30]：

(6)

(7)

式中：c(x,t)為模擬包氣帶或含水層中污染物遷移的初始濃度，mg/L；cd為模擬包氣帶或含水層中污染物遷移后的濃度，mg/L；μ為特征根.

式(6)(7)適用于滲漏條件下，防滲層下方包氣帶和含水層中的污染物遷移分析. 利用式(6)(7)計(jì)算得到包氣帶底部的污染物濃度(將該濃度作為含水層的污染物初始濃度)和地下水中污染物濃度. 式(6)(7)中基本參數(shù)取值參考表3、4.

表4 計(jì)算滲濾液濃度和滲漏的模型參數(shù)取值

注： 1) 為表2中數(shù)據(jù)平均值.

1.3 模擬參數(shù)的選擇

為討論參數(shù)設(shè)置的必要性與取值的合理性，同時為參數(shù)的最適宜取值提供參考，從填埋場選址與設(shè)計(jì)的角度對現(xiàn)有填埋場選址與設(shè)計(jì)標(biāo)準(zhǔn)中的參數(shù)進(jìn)行敏感性分析[31]. 由表5可見：從設(shè)計(jì)的角度，由于導(dǎo)排層具有快速導(dǎo)流滲濾液、降低滲漏風(fēng)險(xiǎn)的功能，對滲漏量與污染物濃度影響相對較大，建議選取導(dǎo)排層參數(shù)中的導(dǎo)排支管間距、導(dǎo)排層坡度和導(dǎo)排層滲透系數(shù)3個參數(shù)；從選址的角度，主要考慮滲濾液的滲漏途徑，選取天然襯層滲透系數(shù)、含水層厚度和地下水流速3個參數(shù).

表5 在地下介質(zhì)中模擬污染物運(yùn)輸?shù)哪Ｐ蛥?shù)取值

圖2 不同導(dǎo)排支管間距下地下水中污染物濃度的空間分布Fig.2 Spatial distribution of pollutant concentration in groundwater under different spacing of guide and exhaust branch pipes

2 各參數(shù)對污染物濃度的影響

2.1 導(dǎo)排支管間距對污染物濃度的影響

圖3 支管間距對滲濾液飽和水位的影響[32]Fig.3 Influence of branch pipe spacing on saturated leachate water level[32]

圖2為不同導(dǎo)排支管間距下地下水中污染物濃度的空間分布. 從圖2可以看出，隨著導(dǎo)排管間距增大，監(jiān)測井中污染物的濃度也逐漸增加，以240 m處為例，導(dǎo)排支管間距從5 m增至10、25、50、75、100 m過程中，污染物濃度分別增加了4.00×10-6、8.24×10-6、1.26×10-5、1.55×10-5、1.65×10-5mgL. 這是因?yàn)閷?dǎo)排支管間距增大，使得滲濾液在導(dǎo)排層中的滲流途徑增長(見圖3)[32]，防滲膜上滲濾液飽和水位上升，水頭壓也上升. 根據(jù)達(dá)西滲流定律，滲透流量與水頭壓差成正比，所以隨著滲漏量的增加，地下水污染程度將加劇.

2.2 導(dǎo)排層坡度對污染物濃度的影響

圖4為不同導(dǎo)排層坡度下地下水中污染物濃度的空間分布. 從圖4可以看出，隨著導(dǎo)排層坡度增大，監(jiān)測井處污染物的濃度減小[6]. 以240 m處為例，導(dǎo)排層坡度從0.01增至0.02、0.03、0.04、0.06、0.10和0.16過程中，污染物濃度分別降低了5.67×10-6、8.23×10-6、9.67×10-6、1.13×10-5、1.36×10-5mgL. 主要原因是，隨著導(dǎo)排層坡度增大，導(dǎo)排介質(zhì)中的水力坡度會增大，使得導(dǎo)排層的側(cè)向?qū)拍芰υ鰪?qiáng)、防滲膜上滲濾液飽和水位降低、滲漏量減少.

圖4 不同導(dǎo)排層坡度下地下水中污染物濃度的空間分布Fig.4 Spatial distribution of pollutant concentration in groundwater under different slope of drainage layer

圖5 不同導(dǎo)排層滲透系數(shù)下地下水中污染物濃度的空間分布Fig.5 Spatial distribution of pollutant concentration in groundwater under different permeability coefficients of guide layers

2.3 導(dǎo)排層滲透系數(shù)對污染物濃度的影響

根據(jù)美國EPA標(biāo)準(zhǔn)[33]規(guī)定，導(dǎo)排層滲透系數(shù)應(yīng)不小于0.1 cms，但根據(jù)文獻(xiàn)[34]，滲濾液中的鈣鎂離子、懸浮顆粒、有機(jī)質(zhì)等均可導(dǎo)致導(dǎo)排層發(fā)生淤堵，導(dǎo)排層滲透系數(shù)減小1～3個數(shù)量級. 圖5為不同導(dǎo)排層滲透系數(shù)下地下水中污染物濃度的空間分布. 從圖5可以看出，隨著導(dǎo)排層滲透系數(shù)減小，監(jiān)測井中污染物的濃度增加，含水層中其他位置的污染物濃度也增加. 以240 m處為例，導(dǎo)排層滲透系數(shù)從0.1 cms 減至0.05、0.01、0.005、0.001、0.000 5、0.000 1 cms的過程中，污染物濃度分別降低了4.78×10-4、5.91×10-4、7.25×10-4、7.79×10-4、8.15×10-4、8.22×10-4mgL. 主要原因是，導(dǎo)排層滲透系數(shù)減小、側(cè)向排水能力減弱，導(dǎo)致防滲膜上滲濾液飽和水位增加，滲漏量增加.

2.4 天然襯層滲透系數(shù)對污染物濃度的影響

圖6 不同天然襯層滲透系數(shù)下地下水中污染物濃度的空間分布Fig.6 Spatial distribution of pollutant concentration in groundwater under different permeability coefficient of natural linings

圖6為不同天然襯層滲透系數(shù)下地下水中污染物濃度的空間分布. 從圖6可以看出，隨著天然襯層滲透系數(shù)增加，監(jiān)測井中污染物的濃度增大，含水層中其他位置的污染物濃度也增大. 以240 m處為例，天然襯層滲透系數(shù)從5.0×10-5cms減至1.0×10-5、5.0×10-6、1.0×10-6、5.0×10-7、1.0×10-7和5.0×10-8cms過程中，污染物濃度分別降低了6.50×10-4、8.86×10-4、9.06×10-4、9.36×10-4、9.36×10-4、9.36×10-4mgL. 根據(jù)達(dá)西滲透定律，水分在多孔介質(zhì)中的滲透速率和流量與多孔介質(zhì)的天然襯層滲透系數(shù)成正比，因此天然襯層的滲透系數(shù)越大，從防滲膜上漏洞滲漏的滲濾液能以更大的流量和流速到達(dá)包氣帶和地下水中[35].

2.5 含水層厚度對污染物濃度的影響

圖7為不同含水層厚度下地下水中污染物濃度的空間分布. 從圖7可以看出，隨著含水層厚度增加，含水層中其他位置處的污染物濃度逐漸降低. 以240 m處為例，含水層厚度從2.7 m增至5.4、10.8、21.6和43.2 m過程中，污染物濃度分別降低了2.73×10-5、2.87×10-5、2.96×10-5、2.96×10-5mgL. 推測原因可能是：隨著橫向彌散系數(shù)增大，污染物在垂直方向上的擴(kuò)散效應(yīng)增強(qiáng)；含水層越厚其稀釋和自凈能力越強(qiáng)，污染物濃度越低[36].

圖7 不同含水層厚度下地下水中污染物濃度的空間分布Fig.7 Spatial distribution of pollutant concentration in groundwater under different aquifer thickness

圖8 不同地下水流速下地下水中污染物濃度的空間分布Fig.8 Spatial distribution of concentration in groundwater atunder different underground water flow velocities

2.6 地下水流速對污染物濃度的影響

圖8為不同地下水流速下地下水中污染物濃度的空間分布. 從圖8可以看出，地下水流速越大，含水層中的其他位置處的污染物濃度越低，污染程度越小. 以240 m處為例，地下水流速從10 ma增至20、50、100和150 ma過程中，污染物濃度分別降低了1.96×10-5、3.71×10-5、4.59×10-5、4.98×10-5mgL. 主要是因?yàn)榈叵滤魉僭酱?，含水層的自凈和稀釋能力越?qiáng)[37].

3 討論

3.1 各參數(shù)的敏感性分析

對各參數(shù)進(jìn)行敏感性分析的結(jié)果如圖9所示. 由圖9(a)可見，當(dāng)導(dǎo)排支管間距從5 m增至100 m時，污染物濃度增至1.17×10-5mgL，增加速率為1.23×10-7mg(L·m). 導(dǎo)排支管間距主要影響填埋場內(nèi)的滲濾液導(dǎo)排途徑長度，在同等面積的填埋場中，導(dǎo)排支管間距越大，滲濾液導(dǎo)排途徑越長，使得滲濾液最高水位增加，水壓增大，導(dǎo)致滲濾液滲漏風(fēng)險(xiǎn)增大，地下水中污染物濃度增大[38].

由圖9(b)可見：當(dāng)導(dǎo)排層坡度從0.01增至0.04時，污染物濃度迅速減至6.31×10-6mgL，單位坡度減小速率為2.10×10-4mgL；而當(dāng)導(dǎo)排層坡度從0.04增至0.16時，單位坡度減小速率降低，僅為2.08×10-5mgL；而當(dāng)導(dǎo)排層坡度繼續(xù)增大時，其變化對監(jiān)測井中污染物濃度的影響已經(jīng)可以忽略. 由圖9(c)可見：當(dāng)導(dǎo)排層滲透系數(shù)從0.000 1 cms增至0.01 cms時，污染物濃度迅速減至1.96×10-4mgL，減小速率為0.20 (mg·s)(L·cm)；而當(dāng)導(dǎo)排層滲透系數(shù)從0.01 cms增至0.1 cms時，減小速率降低，僅為3.15×10-4(mg·s)(L·cm)；而當(dāng)導(dǎo)排層滲透系數(shù)繼續(xù)增大時，其變化對監(jiān)測井中污染物濃度的影響已經(jīng)可以忽略. 研究[39-40]表明，導(dǎo)排層坡度和導(dǎo)排層滲透系數(shù)是滲濾液導(dǎo)排能力的主要影響因素，導(dǎo)排層坡度和導(dǎo)排層滲透系數(shù)增大會提高滲濾液的導(dǎo)排效率，降低滲濾液飽和液位，進(jìn)而降低滲濾液滲漏源強(qiáng)，使得地下水中污染物濃度減小.

由圖9(d)可見：當(dāng)天然襯層滲透系數(shù)從5.0×10-8cms增至1.0×10-6cms時，污染物濃度迅速增至3.30×10-5mgL，增加速率為37.74 (mg·s)(L·cm)；而當(dāng)天然襯層滲透系數(shù)從1.0×10-6cms增至5.0×10-5cms 時，增加速率降低，僅為11.75 (mg·s)(L·cm). 可見，天然襯層滲透系數(shù)對于滲漏污染風(fēng)險(xiǎn)控制具有重要作用，會直接影響滲濾液進(jìn)入地下水的穿透時間，以保證含水層具有更充分的時間對污染物進(jìn)行稀釋緩沖[41-42].

由圖9(e)可見：當(dāng)含水層厚度從2.7 m增至5.4 m時，污染物濃度迅速減至1.55×10-5mgL，減小速率為5.74×10-6mg(L·m)；而當(dāng)含水層厚度從5.4 m增至43.2 m時，減小速率降低，僅為5.58×10-7mg(L·m)；而當(dāng)厚度繼續(xù)增大時，其變化對監(jiān)測井中污染物濃度的影響已經(jīng)可以忽略. 由圖9(f)可見：當(dāng)?shù)叵滤魉購?0 ma增至50 ma時，污染物濃度迅速減至1.61×10-5mgL，減小速率為4.03×10-7(mg·a)(L·m)；而當(dāng)?shù)叵滤魉購?0 ma增至150 ma時，減小速率降低，僅為7.7×10-8(mg·a)(L·m)；而當(dāng)?shù)叵滤魉倮^續(xù)增大時，其變化對監(jiān)測井中污染物濃度的影響已經(jīng)可以忽略. 含水層厚度與地下水流速對污染物濃度的主要影響在于稀釋擴(kuò)散作用，隨著稀釋擴(kuò)散作用形成的污染羽面積逐漸增大，污染物濃度逐漸減小，且減小速率會降低，當(dāng)污染羽面積達(dá)到最大時，污染物濃度降至最低[43-44].

為進(jìn)一步細(xì)化分析污染物濃度對各參數(shù)的敏感性，對各參數(shù)進(jìn)行函數(shù)擬合. 污染物濃度與導(dǎo)排支管間距、導(dǎo)排層坡度、導(dǎo)排層滲透系數(shù)、天然襯層滲透系數(shù)、含水層厚度、地下水流速均存在冪函數(shù)關(guān)系，相關(guān)系數(shù)(R2)分別為 1.000 0、0.997 1、0.991 1、1.000 0、0.883 5、0.988 6(見圖9)，說明擬合效果好，可信度高. 然后，對變化參數(shù)取值，根據(jù)污染物濃度變化確定敏感性指數(shù). 對參數(shù)敏感性排序(見表6)可知，導(dǎo)排支管間距>含水層厚度>導(dǎo)排層坡度>天然襯層滲透系數(shù)>導(dǎo)排層滲透系數(shù)>地下水流速，其中導(dǎo)排層支管間距敏感性指數(shù)為正值，污染物濃度會隨著導(dǎo)排支管間距的增加而升高，其余為負(fù)值，污染物濃度會隨著參數(shù)值的增加而降低.

圖9 各參數(shù)敏感性分析Fig.9 Sensitivity analysis of parameters

表6 各參數(shù)敏感性指數(shù)及排序

3.2 各參數(shù)的敏感區(qū)間與推薦取值

根據(jù)3.1節(jié)的討論結(jié)果，確定了6個因素的敏感區(qū)間(見表7). 在敏感區(qū)間中選取合適的參數(shù)值是控制填埋場污染風(fēng)險(xiǎn)的關(guān)鍵. 對于導(dǎo)排層參數(shù)，在經(jīng)濟(jì)情況允許的范圍內(nèi)，建議在敏感區(qū)間中選擇敏感度較高的取值，導(dǎo)排支管間距取10～25 m、導(dǎo)排層坡度取0.02～0.03、導(dǎo)排層滲透系數(shù)取 0.000 5～0.001 cm/s.

對于水文地質(zhì)參數(shù)，需要通過實(shí)地考察，參考相關(guān)地質(zhì)材料確定所選廠址的相關(guān)參數(shù)數(shù)值，綜合各方面因素考慮，選擇敏感區(qū)間中相對合適的取值，建議天然襯層滲透系數(shù)取5.0×10-7～1.0×10-6cm/s、含水層厚度取2.7～5.4 m、地下水流速取20～50 m/a.

與國內(nèi)外相關(guān)標(biāo)準(zhǔn)相比，對于導(dǎo)排支管間距、含水層厚度和地下水流速3個參數(shù)，并沒有作出明確的指示，該研究給出了比較合理的推薦取值. 對于其他3個參數(shù)，導(dǎo)排層坡度的推薦取值與國內(nèi)標(biāo)準(zhǔn)相吻合，美國與德國要求更加嚴(yán)格；導(dǎo)排層滲透系數(shù)與德國標(biāo)準(zhǔn)相吻合，中國和美國要求更加嚴(yán)格；天然襯層滲透系數(shù)與中國和美國標(biāo)準(zhǔn)相吻合，德國要求更加嚴(yán)格(見表7).

表7 各參數(shù)敏感區(qū)間與推薦取值

4 結(jié)論

a) 當(dāng)導(dǎo)排支管間距>25 m、導(dǎo)排層坡度<2%、導(dǎo)排層滲透系數(shù)<0.000 1 cms、天然襯層滲透系數(shù)>1×10-6cms時，地下水污染風(fēng)險(xiǎn)急劇增大；反之，當(dāng)導(dǎo)排支管間距<10 m、導(dǎo)排層坡度>3%、導(dǎo)排層滲透系數(shù)>0.01 cms、天然襯層滲透系數(shù)<5×10-7cms時，對地下水污染風(fēng)險(xiǎn)的管控并無明顯影響.

b) 建議將標(biāo)準(zhǔn)中設(shè)計(jì)參數(shù)的最低技術(shù)要求為導(dǎo)排支管間距≥25 m、導(dǎo)排層坡度≥2%、導(dǎo)排層滲透系數(shù)≥0.01 cms、天然襯層滲透系數(shù)取≥1.0×10-6cms.

c) 各參數(shù)對地下水污染風(fēng)險(xiǎn)的敏感性不同，由大到小依次為導(dǎo)排支管間距>含水層厚度>導(dǎo)排層坡度>天然襯層滲透系數(shù)>導(dǎo)排層滲透系數(shù)>地下水流速，利用指標(biāo)體系方法進(jìn)行填埋場地下水風(fēng)險(xiǎn)評估時，指標(biāo)權(quán)重可參考此排序.

參考文獻(xiàn)(References)：

[1] OMAR H,ROHANI S.Treatment of landfill waste,leachate and landfill gas:a review[J].Frontiers of Chemical Science and Engineering,2015,9(1):15-32.

[2] 徐亞,劉玉強(qiáng),能昌信,等.基于概率風(fēng)險(xiǎn)分析的填埋場安全防護(hù)距離研究[J].中國環(huán)境科學(xué),2017,37(7):2786-2792.

XU Ya,LIU Yuqiang,NAI Changxin,etal.A probabilistic risk-based method to determine the safety protection distance for solid waste landfill[J].China Environmental Science,2017,37(7):2786-2792.

[3] WANG Xiaojun,JIA Mingsheng,LIN Xiangyu,etal.A comparison of CH4,N2O and CO2emissions from three different cover types in a municipal solid waste landfill[J].Air Repair,2017,67(4):507-515.

[4] BRANCHINI L.Municipal waste overview[M].Berlin:Springer International,2015.

[5] GRUGNALETTI M,PANTINI S,VERGINELLI I,etal.An easy-to-use tool for the evaluation of leachate production at landfill sites[J].Waste Management,2016,55:204-219.

[6] XU Ya,XUE Xiangshan,DONG Lu,etal.Long-term dynamics of leachate production,leakage from hazardous waste landfill sites and the impact on groundwater quality and human health[J].Waste Management,2018,82:156-166.

[7] LU Xinwei,WANG Lijiang,LORETTA Y,etal.Multivariate statistical analysis of heavy metals in street dust of Baoji,NW China[J].Journal of Hazardous Materials,2010,173(123):744-749.

[8] LEWIS T W,PIVONKA P,FITYUS S G,etal.Parametric sensitivity analysis of coupled mechanical consolidation and contaminant transport through clay barriers[J].Computers and Geotechnics,2009,36(12):31-40.

[9] 徐亞,劉玉強(qiáng),胡立堂,等.填埋場井筒效應(yīng)及其對污染監(jiān)測井監(jiān)測效果的影響[J].中國環(huán)境科學(xué),2018,38(8):315-322.

XU Ya,LIU Yuqiang,HU Litang,etal.Wellbore effect and its influence on contamination monitor well-a case study on groundwater monitor well around landfill[J].China Environmental Science,2018,38(8):315-322.

[10] XU Ya,DONG Lu,NAI Changxin,etal.Buffering distance between hazardous waste landfill and water supply wells in a shallow aquifer[J].Journal of Cleaner Production,2019,211:1180-1189.

[11] FLEMING I R,ROWE R K.Laboratory studies of clogging of landfill leachate collection and drainage systems[J].Canadian Geotechnical Journal,2004,41(1):134-153.

[12] BEAVEN R P,HUDSON A P,KNOX K,etal.Clogging of landfill tyre and aggregate drainage layers by methanogenic leachate and implications for practice[J].Waste Management,2013,33(2):431-444.

[13] YANG Y S,LI P,ZHANG X,etal.Lab-based investigation of enhanced BTEX attenuation driven by groundwater table fluctuation[J].Chemosphere,2017,169:678-684.

[14] YUE W,JUAN L,DA A,etal.Site selection for municipal solid waste landfill considering environmental health risks[J].Resources,Conservation and Recycling,2018,138:40-46.

[15] US Environmental Protection Agency.EPA106152 hazardous and solid waste amendment[S].Washington DC:Government Printing Office,1984.

[16] 國家環(huán)境保護(hù)總局.GB 18598—2001 危險(xiǎn)廢物填埋污染控制標(biāo)準(zhǔn)[S].北京:中國標(biāo)準(zhǔn)出版社,2001.

[17] 國家環(huán)境保護(hù)總局.GB 16889—2008 生活垃圾填埋場污染控制標(biāo)準(zhǔn)[S].北京:中國標(biāo)準(zhǔn)出版社,2001.

[18] European Environment Agency.9931EC council directive[S].Brussels:Official Journal of the European Communities,1999.

[19] CAPELLI L,SIRONI S,ROSSO R D,etal.Measuring odours in the environment vs.dispersion modelling:a review[J].Atmospheric Environment,2013,79:731-743.

[20] SHU S,ZHU W,WANG S,etal.Leachate breakthrough mechanism and key pollutant indicator of municipal solid waste landfill barrier systems:centrifuge and numerical modeling approach[J].Science of the Total Environment,2018,612:1123-1131.

[21] KOTOWSKA U,KAPELEWSKA J,STURGULEWSKA J.Determination of phenols and pharmaceuticals in municipal wastewaters from polish treatment plants by ultrasound-assisted emulsification-microextraction followed by GC-MS[J].Environment Science and Polluttion Research International,2014,21(1):660-673.

[22] ZHAN L T,CHEN C,WANG Y,etal.Failure probability assessment and parameter sensitivity analysis of a contaminant′s transit time through a compacted clay liner[J].Computers and Geotechnics,2017,86:230-242.

[23] DJOGO M,RADONIC′J,MIHAJLOVIC′I,etal.Selection of optimal parameters for future research monitoring programmes on MSW landfill in Novi Sad,Serbia[J].Fresenius Environmental Bulletin,2017,26:1329-1332.

[24] GIROUD J P,BONAPARTE R.Leakage through liners constructed with geomembranes:Part Ⅰ.geomembrane liners [J].Geotextiles & Geomembranes,1989,8(1):27-67.

[25] GIROUD J P,BONAPARTE R.Leakage through liners constructed with geomembranes:Part Ⅱ.composite liners[J].Geotextiles & Geomembranes,1989,8(2):71-111.

[26] DARILEK G T,PARRA J O.The electrical leak location method for geomembrane liners[J].Journal of Hazardous Materials,1988,21(2):177-187.

[27] KASTMAN K H.Hazardous waste landfill membrane:design,installation and monitoring[J].Proceedings of the International Conference on Geomembranes,1984,1:215-220.

[28] MCQUADE S J,NEEDHAM A D.Geomembrane liner defects-causes,frequency and avoidance[J].Proceedings of the Institution of Civil Engineers-Geotechnical Engineering,1999,137:203-213.

[29] ROSENBERG N D,SPERA F J.Convection in porous media with thermal and chemical buoyancy:a comparison of two models for solute dispersion[M].New York:Chaotic Processes in the Geological Sciences,1992.

[30] AGENCY E.The development of LandSim 2.5[M].New York:The Development of American Industries:Prentice-Hall,1941.

[31] CHANG H,ZHAO Y,TAN H,etal.Parameter sensitivity to concentrations and transport distance of odorous compounds from solid waste facilities[J].Science of the Total Environment,2019,651:2158-2165.

[32] 徐亞,顏湘華,董路,等.基于Landsim的填埋場長期滲漏的污染風(fēng)險(xiǎn)評價(jià)[J].中國環(huán)境科學(xué),2014,34(5):1355-1360.

XU Ya,YAN Xiang Hua,DONG Lu,etal.Pollution risk assessment of long-term leaking in landfill-based on the Landsim model[J].China Environmental Science,2014,34(5):1355-1360.

[33] Center for Enviromental Research Information Office of Research.EPA625489022 requirement for hazardous waste landfill design,construction,and closure[S].Washington DC:US Environmental Protection Agency,1989.

[34] 徐亞,董路,能昌信,等.危廢填埋場導(dǎo)排層淤堵的時空分布特征[J].中國環(huán)境科學(xué),2016,36(3):849-855.

XU Ya,DONG Lu,NAI Changxin,etal.Spatial and temporal characterization of drainage clogging in hazardous waste landfill[J].China Environmental Science,2016,36(3):849-855.

[35] LIAN X,JIA Y F,YANG Y,etal.Identification of groundwater redox process induced by landfill leachate based on sensitive factor method[J].Environmental Science & Pollution Research,2017,24 (34):1-8.

[36] LIU J G,NIE X Q,ZENG X W,etal.Long-term leaching behavior of phenol in cementactivated-carbon solidifiedstabilized hazardous waste[J].Journal of Environmental Management,2013,115:265-269.

[37] AVISHAI L,SIEBNER H,DAHAN O,etal.Using the natural biodegradation potential of shallow soils for in-situ remediation of deep vadose zone and groundwater[J].Journal of Hazardous Material,2017,324:398-405.

[38] BHOWMICK S,PRAMANIK S,SINGH P,etal.Arsenic in groundwater of West Bengal,India:a review of human health risks and assessment of possible intervention options[J].Science of the Total Environment,2018,612:148-169.

[39] DAN A,MASAO O,YUTA F,etal,Removal of heavy metals from synthetic landfill leachate in lab-scale vertical flow constructed wetlands[J].Science of the Total Environment,2017,584585:742-750.

[40] HAN Z Y,MA H N,SHI G Z,etal.A review of groundwater contamination near municipal solid waste landfill sites in China[J].Science of the Total Environment,2016,569570:1255-1264.

[41] LIOR A,HAGAR S,OFER D,etal.Using the natural biodegradation potential of shallow soils for in-situ remediation of deep vadose zone and groundwater[J].Journal of Hazardous Materials,2017,324:398-405.

[42] SAEEDREZA H,AMITY G,ZIMMER-FAUS T,etal.Remediation of groundwater contaminated with arsenic through enhanced natural attenuation:batch and column studies[J].Water Research,2017,122:545-556.

[43] YING S C,SCHAEFER M V,COCKESTEB A,etal.Depth stratification leads to distinct zones of manganese and arsenic contaminated groundwater[J].Environmental Science & Technology,2017,51:8926-8932.

[44] YANG Y S,LI P,ZHANG X,etal.Lab-based investigation of enhanced BTEX attenuation driven by groundwater table fluctuation[J].Chemosphere,2017,169:678-684.