危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)評(píng)價(jià)中指標(biāo)權(quán)重計(jì)算方法優(yōu)化比選

韓 旭, 生 賀, 夏 甫, 徐祥健, 尚長(zhǎng)健, 楊 昱, 姜永海

中國(guó)環(huán)境科學(xué)研究院， 國(guó)家環(huán)境保護(hù)地下水污染模擬與控制重點(diǎn)實(shí)驗(yàn)室， 北京 100012

隨著我國(guó)生態(tài)環(huán)保理念的推廣，相關(guān)部門對(duì)場(chǎng)地風(fēng)險(xiǎn)管理也投入了更高的關(guān)注，編制了多項(xiàng)管理辦法及標(biāo)準(zhǔn)作為政策支持. 對(duì)場(chǎng)地進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià)是場(chǎng)地管理的一個(gè)重要方向和環(huán)節(jié)，通過(guò)風(fēng)險(xiǎn)評(píng)價(jià)不僅可以確定場(chǎng)地中的風(fēng)險(xiǎn)點(diǎn)也可以判斷各場(chǎng)地的風(fēng)險(xiǎn)大小[1]，以此確定風(fēng)險(xiǎn)管理的場(chǎng)地優(yōu)先級(jí)排序和重點(diǎn)關(guān)注指標(biāo)[2-4].

由于工業(yè)場(chǎng)地建設(shè)水平參差不齊、場(chǎng)地特性各有差異，并且各個(gè)場(chǎng)地的風(fēng)險(xiǎn)特征差異較大，要想針對(duì)已建成的工業(yè)場(chǎng)地實(shí)施風(fēng)險(xiǎn)評(píng)價(jià)，就需要選擇適用于同一類場(chǎng)地的風(fēng)險(xiǎn)評(píng)價(jià)方法. Visual MODFLOW[5-6]、HYDRUS-1D[7-9]等數(shù)值模擬方法可以考慮場(chǎng)地中裝置、源-匯項(xiàng)等多種因素，并能描述多維、多相態(tài)、多組分、多過(guò)程的污染物遷移轉(zhuǎn)化過(guò)程，形成定量化的描述，為地下水污染風(fēng)險(xiǎn)評(píng)價(jià)提供基礎(chǔ). 數(shù)值模擬方法涵蓋全面、因素多樣、模擬結(jié)果貼合實(shí)際，但過(guò)程復(fù)雜、需要專業(yè)人員操作、工作量大，因此，建立一種簡(jiǎn)便、可操作性強(qiáng)、結(jié)果準(zhǔn)確的污染風(fēng)險(xiǎn)評(píng)價(jià)方法對(duì)于部門管理人員更具有可行性. 場(chǎng)地地下水污染風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)方法是通過(guò)構(gòu)建評(píng)價(jià)指標(biāo)體系、分配指標(biāo)權(quán)重、計(jì)算風(fēng)險(xiǎn)綜合指數(shù)分析地下水污染風(fēng)險(xiǎn)[10]，指標(biāo)數(shù)量少、計(jì)算方法簡(jiǎn)單，且結(jié)果可以直接進(jìn)行排序從而確定風(fēng)險(xiǎn)管控的關(guān)鍵點(diǎn)，最重要的是該方法可以體現(xiàn)污染源釋放特征等生產(chǎn)復(fù)雜性. 權(quán)重計(jì)算是風(fēng)險(xiǎn)指數(shù)計(jì)算中一個(gè)至關(guān)重要的環(huán)節(jié). 常用的權(quán)重計(jì)算方法歸類為主觀賦權(quán)法、客觀賦權(quán)法和綜合賦權(quán)法3種[11]. 主觀賦權(quán)法根據(jù)決策者(專家)主觀上對(duì)各屬性的重視程度來(lái)確定屬性權(quán)重，包含專家調(diào)查法(Delphi法)[12]、層次分析法(AHP)[13-15]等. 客觀賦權(quán)法可以完全采用已有數(shù)據(jù)，不加入人為判斷進(jìn)行權(quán)重計(jì)算，但是缺乏人類經(jīng)驗(yàn)的總結(jié)，常用的包含熵權(quán)法[16-17]、模糊集理論[18]、主成分分析法[19]、離差及均方差法等[20]. 主客觀綜合賦權(quán)法，即分別計(jì)算主觀、客觀權(quán)重，利用加法或乘法將2項(xiàng)權(quán)重進(jìn)行融合[21].

鐘秀等[10]運(yùn)用層次分析法確定了地下水飲用水源地污染源風(fēng)險(xiǎn)等級(jí)評(píng)價(jià)指標(biāo)體系中21個(gè)指標(biāo)的權(quán)重，通過(guò)加權(quán)求和模型，構(gòu)建了污染源風(fēng)險(xiǎn)等級(jí)評(píng)價(jià)指數(shù)的計(jì)算方法，并針對(duì)不同風(fēng)險(xiǎn)區(qū)域提出了修復(fù)技術(shù)要求. 屈飛行等[11]利用層次分析法的4種不同計(jì)算方法進(jìn)行了地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的指標(biāo)權(quán)重計(jì)算，并采用邊界確定的實(shí)際滑坡及崩塌面對(duì)評(píng)價(jià)結(jié)果進(jìn)行驗(yàn)證，確定了地質(zhì)災(zāi)害危害性評(píng)價(jià)時(shí)優(yōu)先選擇特征向量法計(jì)算權(quán)重. 李艷萍等[22]利用主觀賦權(quán)法中的層次分析法對(duì)工業(yè)園區(qū)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)19項(xiàng)指標(biāo)的權(quán)重進(jìn)行了計(jì)算，得出基于不同人群打出的分值而得到的權(quán)重結(jié)果各有不同，體現(xiàn)了人群認(rèn)識(shí)不同，則層次分析法標(biāo)度不同，對(duì)計(jì)算結(jié)果的影響程度較大的觀點(diǎn). 李帥等[23]利用AHP和熵權(quán)法結(jié)合的方法確定指標(biāo)權(quán)重，并對(duì)城市人居環(huán)境質(zhì)量大小進(jìn)行綜合評(píng)價(jià)，確定了5個(gè)城市的人居環(huán)境質(zhì)量大小順序，并相應(yīng)提出了建議.

為了探討層次分析法(主觀賦權(quán)法)、熵權(quán)法(客觀賦權(quán)法)、層次分析-熵權(quán)法(綜合賦權(quán)法)對(duì)于風(fēng)險(xiǎn)評(píng)價(jià)的適用性，以我國(guó)37家危險(xiǎn)廢物填埋場(chǎng)為案例，進(jìn)行地下水污染風(fēng)險(xiǎn)評(píng)價(jià)，通過(guò)風(fēng)險(xiǎn)指數(shù)與地下水污染指數(shù)擬合結(jié)果進(jìn)行分析比較，確定風(fēng)險(xiǎn)指數(shù)計(jì)算中推薦選擇的賦權(quán)方法.

1 數(shù)據(jù)來(lái)源及研究方法

1.1 研究方法

場(chǎng)地地下水污染風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)計(jì)算的基本步驟：①構(gòu)建指標(biāo)體系；②指標(biāo)數(shù)值量化；③分配指標(biāo)權(quán)重；④計(jì)算風(fēng)險(xiǎn)指數(shù)；⑤結(jié)果驗(yàn)證. 具體如下：首先采用層次分析法構(gòu)建指標(biāo)體系，之后利用聚類分析法和賦分法將各指標(biāo)進(jìn)行統(tǒng)一量化，然后分別采用層次分析法、熵權(quán)法、層次分析-熵權(quán)法計(jì)算指標(biāo)權(quán)重，最后采用積和法計(jì)算地下水污染風(fēng)險(xiǎn)綜合指數(shù)，并利用場(chǎng)地調(diào)查獲取的地下水苯濃度進(jìn)行風(fēng)險(xiǎn)指數(shù)與地下水污染指數(shù)的線性擬合，確定最優(yōu)的權(quán)重計(jì)算方法. 其中權(quán)重計(jì)算、風(fēng)險(xiǎn)指數(shù)計(jì)算及結(jié)果驗(yàn)證方法介紹如下.

1.1.1權(quán)重計(jì)算

a) 層次分析法權(quán)重λj的計(jì)算步驟[24-25]：①建立遞階層次結(jié)構(gòu)模型，繪制層次結(jié)構(gòu)圖；②依據(jù)表1所示的標(biāo)度方法構(gòu)造各層次中的成對(duì)比較矩陣；③針對(duì)每個(gè)成對(duì)比較矩陣計(jì)算層次單排序及進(jìn)行一致性檢驗(yàn)；④層次總排序及一致性檢驗(yàn). 一致性檢驗(yàn)中的指標(biāo)CI計(jì)算如式(1)所示，一致性比率(CR)計(jì)算如式(2)所示.

表1 判斷矩陣元素aij的標(biāo)度方法

(1)

式中，λ為最大特征根.

b) 熵權(quán)法[26]的計(jì)算步驟：①按照式(3)將各指標(biāo)的數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，構(gòu)建標(biāo)準(zhǔn)化矩陣xij′；②利用式(4)計(jì)算各指標(biāo)的信息熵ej；③通過(guò)式(5)計(jì)算各指標(biāo)權(quán)重μj.

(3)

(4)

(5)

c) 層次分析-熵權(quán)法采用“加法”集成法[27-28]，具體步驟是：①構(gòu)建2種權(quán)重的距離函數(shù)D(λj,μj)；②計(jì)算分布系數(shù)α和β；③根據(jù)組合積分權(quán)重計(jì)算公式得到綜合權(quán)重ωj.

(6)

D(λj,μj)2=(α-β)2

(7)

α+β=1

(8)

ωj=αλj+βμj

(9)

式中，λj和μj分別為層次分析法權(quán)重和熵權(quán)法權(quán)重計(jì)算結(jié)果，α和β分別為λj和μj的分布系數(shù).

1.1.2風(fēng)險(xiǎn)指數(shù)計(jì)算

地下水污染風(fēng)險(xiǎn)綜合指數(shù)(R)利用綜合指數(shù)法進(jìn)行計(jì)算.

R=R1W1+R2W2+…+RnWn

(10)

式中，Rn為各指標(biāo)分值，Wn為各指標(biāo)權(quán)重，n為指標(biāo)序號(hào).

1.1.3污染指數(shù)與風(fēng)險(xiǎn)指數(shù)擬合

利用距平指數(shù)法構(gòu)造場(chǎng)地地下水污染指數(shù)C，計(jì)算公式如式(11)所示，并利用軟件將C與R進(jìn)行線性擬合.

(11)

1.2 數(shù)據(jù)來(lái)源

該研究以我國(guó)37家危險(xiǎn)廢物填埋場(chǎng)為案例，采用的數(shù)據(jù)均來(lái)自于收集的危險(xiǎn)廢物填埋場(chǎng)可行性研究報(bào)告、地質(zhì)勘探報(bào)告及現(xiàn)場(chǎng)調(diào)研獲得的信息，具體見(jiàn)表2. 37家危險(xiǎn)廢物填埋場(chǎng)分別位于我國(guó)37個(gè)城市，占地面積普遍為1.0×105～1.5×105m2，填埋場(chǎng)填埋年限10～30年，涉及丘陵、平原、臺(tái)地、溝谷4種地形地貌，地下水類型為孔隙水或裂隙水或巖溶水，地下水埋深普遍小于10 m.

2 案例分析

2.1 構(gòu)建指標(biāo)體系

構(gòu)建指標(biāo)體系首先需要明確風(fēng)險(xiǎn)原因，危險(xiǎn)廢物填埋場(chǎng)地下水污染的主要風(fēng)險(xiǎn)來(lái)源于滲濾液的泄露，造成滲濾液泄漏的主要因素一般為場(chǎng)內(nèi)原因，也就是由于填埋場(chǎng)自身構(gòu)造造成的[29-31]. 而從滲濾液泄漏到產(chǎn)生地下水污染風(fēng)險(xiǎn)，涉及填埋場(chǎng)建設(shè)特性、滲濾液性質(zhì)、場(chǎng)地水文地質(zhì)情況、受體距離及管理水平等[32-33]. 借鑒國(guó)內(nèi)外針對(duì)危險(xiǎn)廢物填埋場(chǎng)的風(fēng)險(xiǎn)研究，結(jié)合所收集的我國(guó)37家危險(xiǎn)廢物填埋場(chǎng)資料，利用層次分析法確定危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)評(píng)價(jià)45個(gè)指標(biāo)的初步指標(biāo)體系[34]. 基于危險(xiǎn)廢物滲濾液泄露污染地下水的3個(gè)層面進(jìn)行主要影響因素篩選：一是危險(xiǎn)廢物填埋場(chǎng)本身的風(fēng)險(xiǎn)指標(biāo)，表征由于填埋場(chǎng)自身?xiàng)l件和性質(zhì)而導(dǎo)致地下水受到污染的風(fēng)險(xiǎn)，包括填埋場(chǎng)建設(shè)規(guī)模和滲濾液及防滲系統(tǒng)特性；二是體現(xiàn)危險(xiǎn)廢物填埋場(chǎng)地下水污染物垂向遷移風(fēng)險(xiǎn)指標(biāo)，表征污染物通過(guò)包氣帶介質(zhì)進(jìn)入含水層的風(fēng)險(xiǎn)，包括水文條件和包氣帶介質(zhì)特性；三是體現(xiàn)危險(xiǎn)廢物填埋場(chǎng)地下水污染物水平遷移風(fēng)險(xiǎn)指標(biāo)，表征污染物進(jìn)入含水層后，隨著地下水流動(dòng)遷移擴(kuò)散至飲用水源地的風(fēng)險(xiǎn)，包含含水層介質(zhì)特性和飲用水源地距離填埋場(chǎng)的距離. 最終建立了包含14項(xiàng)指標(biāo)的指標(biāo)體系，見(jiàn)表3.

表3 危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)

2.2 指標(biāo)數(shù)值量化

為了減少指標(biāo)量化賦分差異對(duì)綜合指數(shù)結(jié)果的影響，選擇通過(guò)聚類分析法的方法，將收集到的所有指標(biāo)數(shù)據(jù)分為3個(gè)等級(jí)，部分指標(biāo)分級(jí)標(biāo)準(zhǔn)見(jiàn)表4. 按照單個(gè)指標(biāo)可能造成地下水污染風(fēng)險(xiǎn)等級(jí)高、中、低進(jìn)行賦分：若指標(biāo)屬于高風(fēng)險(xiǎn)則為一級(jí)，分值為50分； 若指標(biāo)屬于中風(fēng)險(xiǎn)則為二級(jí)，分值為30分； 若指標(biāo)為低風(fēng)險(xiǎn)則為三級(jí)，分值為10分.

表4 危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)部分指標(biāo)分級(jí)標(biāo)準(zhǔn)

2.3 分配指標(biāo)權(quán)重

根據(jù)式(1)～(9)計(jì)算14項(xiàng)指標(biāo)的權(quán)重，3種權(quán)重方法對(duì)應(yīng)的權(quán)重結(jié)果分別用λj、μj及ωj表示，結(jié)果見(jiàn)表5.

表5 指標(biāo)權(quán)重計(jì)算結(jié)果

2.4 計(jì)算風(fēng)險(xiǎn)指數(shù)

利用式(10)計(jì)算37家危險(xiǎn)廢物填埋場(chǎng)的地下水污染風(fēng)險(xiǎn)綜合指數(shù)R，3種權(quán)重方法對(duì)應(yīng)的地下水污染風(fēng)險(xiǎn)綜合指數(shù)分別用Rλ、Rμ和Rω表示，計(jì)算結(jié)果取2位有效數(shù)字，見(jiàn)表6.

表6 37個(gè)危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)綜合指數(shù)

3 結(jié)果與討論

3.1 權(quán)重計(jì)算結(jié)果對(duì)比

按照權(quán)重大小排序形成3組權(quán)重對(duì)比圖. 由圖1可知，層次分析法和熵權(quán)法計(jì)算得出的評(píng)價(jià)指標(biāo)權(quán)重有比較明顯的差異. 層次分析法中權(quán)重前5名由大到小依次為含水層滲透系數(shù)、包氣帶滲透系數(shù)、地下水埋深、含水層厚度、防滲膜HDPE厚度；熵權(quán)法中權(quán)重前5名由大到小依次為含水層滲透系數(shù)、包氣帶滲透系數(shù)、地形坡度、年填埋量、滲濾液產(chǎn)生量. 2種方法差異的主要原因在于2種計(jì)算方法的內(nèi)涵，層次分析法反映的是基于專家認(rèn)知得出的指標(biāo)重要性排序，熵權(quán)法是單純通過(guò)指標(biāo)數(shù)值計(jì)算得到的權(quán)重，反映的是指標(biāo)數(shù)據(jù)差異的大小.

圖1 3組權(quán)重對(duì)比Fig.1 Comparison chart of three sets of weights

層次分析-熵權(quán)法是基于層次分析法和熵權(quán)法二者的分配系數(shù)綜合得出的權(quán)重，層次分析-熵權(quán)法中權(quán)重前5名由大到小依次為含水層滲透系數(shù)、包氣帶滲透系數(shù)、含水層厚度、防滲膜HDPE厚度、地下水埋深. 可以看出，層次分析-熵權(quán)法與層次法得出的權(quán)重前5名指標(biāo)完全相同，只是排序不同，這是由于層次分析-熵權(quán)法中層次分析法分配系數(shù)為0.58，熵權(quán)法分配系數(shù)為0.42，在計(jì)算綜合權(quán)重時(shí)，層次分析法的權(quán)重具有更大的比重，而以上5個(gè)指標(biāo)的λj都要大于μj，因此只有指標(biāo)排序發(fā)生了變化.

3種方法的結(jié)果中前5名同時(shí)包含含水層滲透系數(shù)和包氣帶滲透系數(shù)，說(shuō)明這2項(xiàng)指標(biāo)在不同的危險(xiǎn)廢物填埋場(chǎng)的數(shù)值有著明顯差異，也是可能造成危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)的兩大主要因素. 3種方法均得出含水層滲透系數(shù)和包氣帶滲透系數(shù)是可能造成地下水污染風(fēng)險(xiǎn)的兩大指標(biāo)，說(shuō)明在建設(shè)或運(yùn)行管理危險(xiǎn)廢物填埋場(chǎng)時(shí)，為保護(hù)地下水選址需要優(yōu)先考慮含水層及包氣帶介質(zhì)類型，必要時(shí)采取更高性能的防滲技術(shù)手段.

3.2 污染風(fēng)險(xiǎn)綜合指數(shù)結(jié)果對(duì)比

依據(jù)3組權(quán)重得到的37個(gè)危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)綜合指數(shù)如圖2所示，按照風(fēng)險(xiǎn)指數(shù)由大到小排序前15家危險(xiǎn)廢物填埋場(chǎng)如表7所示. 由表7可知，3組地下水污染風(fēng)險(xiǎn)綜合指數(shù)排序前15名危險(xiǎn)廢物填埋場(chǎng)中共有10家同時(shí)出現(xiàn)在3組結(jié)果中，此外，層次分析-熵權(quán)法中的28號(hào)和14號(hào)場(chǎng)地、以及21和11號(hào)場(chǎng)地分別與層次分析法和熵權(quán)法的結(jié)果對(duì)應(yīng)，僅有3號(hào)場(chǎng)地沒(méi)有出現(xiàn)在另外2組結(jié)果中，由此可見(jiàn)層次分析-熵權(quán)法將2種權(quán)重進(jìn)行了融合，而層次分析法和熵權(quán)法雖然有13家重合，但是排序完全不同.

表7 危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)綜合指數(shù)前15名及其計(jì)算結(jié)果

圖2 3組地下水污染風(fēng)險(xiǎn)綜合指數(shù)對(duì)比Fig.2 Comparison chart of three groups of comprehensive index of groundwater pollution risk

進(jìn)一步分析地下水污染風(fēng)險(xiǎn)前5名的場(chǎng)地，按照層次分析法的結(jié)果排序?yàn)?號(hào)>36號(hào)>19號(hào)>18號(hào)>29號(hào)；按照熵權(quán)法的結(jié)果排序?yàn)?2號(hào)>19號(hào)>16號(hào)>7號(hào)>23號(hào)；按照層次分析-熵權(quán)法的結(jié)果排序?yàn)?號(hào)>18號(hào)>9號(hào)>10號(hào)>14號(hào). 在危險(xiǎn)廢物填埋場(chǎng)指標(biāo)分值固定的情況下，地下水污染風(fēng)險(xiǎn)排序因?yàn)?種權(quán)重計(jì)算方法得出的權(quán)重結(jié)果而表現(xiàn)出顯著差異，僅有7號(hào)危險(xiǎn)廢物填埋場(chǎng)在3組結(jié)果中同時(shí)出現(xiàn). 可見(jiàn)，按照不同的權(quán)重計(jì)算方法判斷的場(chǎng)地風(fēng)險(xiǎn)管理優(yōu)先級(jí)是不同的，只有選用能更準(zhǔn)確匹配實(shí)際情況的權(quán)重計(jì)算方法才可以保證風(fēng)險(xiǎn)管理實(shí)施的有效性和準(zhǔn)確性.

3.3 權(quán)重計(jì)算的優(yōu)化比選

3組污染風(fēng)險(xiǎn)指數(shù)(R)與苯污染指數(shù)(C)之間的擬合結(jié)果見(jiàn)圖3. 結(jié)果顯示，3組污染風(fēng)險(xiǎn)綜合指數(shù)均可以和苯濃度形成線性關(guān)系，層次分析法得出的R2為0.75，熵權(quán)法得出的R2為0.51，層次分析-熵權(quán)法得出的R2為0.84，說(shuō)明層次分析-熵權(quán)法優(yōu)于層次分析法，更優(yōu)于熵權(quán)法，層次分析-熵權(quán)法得出的污染風(fēng)險(xiǎn)綜合指數(shù)更能準(zhǔn)確的描述場(chǎng)地的污染情況，即相對(duì)于另外2種權(quán)重計(jì)算方法，綜合權(quán)重法更適用于危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)評(píng)價(jià).

圖3 3種權(quán)重計(jì)算方法下風(fēng)險(xiǎn)指數(shù)與污染指數(shù)的擬合結(jié)果Fig.3 Relationship between risk index and contamination indexunder three weight calculation methods

將表5中3組權(quán)重分別按由大到小排序后，可以看到3組排序完全不同，依次加和得到的總權(quán)重如圖4所示. 當(dāng)總權(quán)重達(dá)到0.96時(shí)，層次分析-熵權(quán)法包含了12項(xiàng)指標(biāo)，層次分析法包含了11項(xiàng)指標(biāo)，熵權(quán)法包含了9項(xiàng)指標(biāo)，層次分析-熵權(quán)法相較于其他2種方法可以在一定權(quán)重范圍內(nèi)體現(xiàn)更多的指標(biāo)，且權(quán)重分配更加均衡. 層次分析-熵權(quán)法既體現(xiàn)了對(duì)地下水污染較大的2個(gè)指標(biāo)的影響，同時(shí)平衡了數(shù)據(jù)差異與專家認(rèn)知之間不同，由此建立的綜合指數(shù)計(jì)算方法更加穩(wěn)定，不易受到單個(gè)指標(biāo)缺失的影響，在一定程度上可以為危險(xiǎn)廢物填埋場(chǎng)地下水污染的風(fēng)險(xiǎn)評(píng)價(jià)及運(yùn)行管理提供支持.

圖4 指標(biāo)權(quán)重加和結(jié)果Fig.4 Sum results of index weights

4 結(jié)論

a) 從不同權(quán)重計(jì)算方法得出的權(quán)重可知，含水層滲透系數(shù)和包氣帶滲透系數(shù)是可能影響地下水污染風(fēng)險(xiǎn)最顯著的2個(gè)指標(biāo)，說(shuō)明在建設(shè)或運(yùn)行管理危險(xiǎn)廢物填埋場(chǎng)時(shí)，為保護(hù)地下水選址時(shí)需要優(yōu)先考慮含水層及包氣帶介質(zhì)類型或在建設(shè)中采取更高性能的防滲方法.

b) 層次分析-熵權(quán)法得出污染風(fēng)險(xiǎn)綜合指數(shù)與場(chǎng)地的苯污染程度擬合度更高，即相對(duì)于層次分析法和熵權(quán)法，層次分析-熵權(quán)法更適用于危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)評(píng)價(jià).

c) 層次分析-熵權(quán)法相對(duì)于其他2種方法權(quán)重分配更為均衡，建立的綜合指數(shù)計(jì)算方法更加穩(wěn)定，不易受到單個(gè)指標(biāo)缺失的影響，在一定程度上可以為危險(xiǎn)廢物填埋場(chǎng)地下水污染的風(fēng)險(xiǎn)評(píng)價(jià)及運(yùn)行管理提供支持.

d) 該研究建立的污染風(fēng)險(xiǎn)評(píng)價(jià)方法適用于危險(xiǎn)廢物填埋場(chǎng)地下水污染風(fēng)險(xiǎn)的批量評(píng)估與分級(jí)管理，該方法操作簡(jiǎn)便、所需數(shù)據(jù)量小且容易獲取，與常規(guī)場(chǎng)地風(fēng)險(xiǎn)評(píng)價(jià)方法相比，大大節(jié)約人力、物力與時(shí)間成本. 在實(shí)施風(fēng)險(xiǎn)管理時(shí)，可配合詳細(xì)調(diào)查分析場(chǎng)地風(fēng)險(xiǎn)特征，進(jìn)而提出有針對(duì)性的風(fēng)險(xiǎn)防控方案.

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

[1] KAJENTHIRA A,HOLMES J,MCDONNELL R.The role of qualitative risk assessment in environmental management:a Kazakhstani case study[J].Scienceof the Total Environment,2012,420:24-32.

[2] 尹榮堯,楊瀟,孫翔,等.江蘇沿海化工區(qū)環(huán)境風(fēng)險(xiǎn)分級(jí)及優(yōu)先管理策略研究[J].中國(guó)環(huán)境科學(xué),2011,31(7):1225-1232.

YIN Rongyao,YANG Xiao,SUN Xiang,etal.Environmental risk rankings and prior management strategies on chemical industry parks in the coastal area of Jiangsu Province[J].China Environmental Science,2011,31(7):1225-1232.

[3] 呂培辰,李舒,馬宗偉,等.中國(guó)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)體系的完善:來(lái)自美國(guó)的經(jīng)驗(yàn)和啟示[J].環(huán)境監(jiān)控與預(yù)警,2018,10(2):1-5.

LU Peichen,LI Shu,MA Zongwei,etal.Improvement of China′s environmental risk assessment system:experience and enlightenment from the United States [J].Environmental Monitoring and Forewarning,2018,10(2):1-5.

[4] 馮霞.環(huán)境風(fēng)險(xiǎn)管理中比較風(fēng)險(xiǎn)評(píng)價(jià)法的應(yīng)用[J].中國(guó)資源綜合利用,2020,38(2):24-26.

FENG Xia.Application of comparative risk evaluation in environmental risk management[J].China Resources Comprehensive Utilization,2020,38(2):24-26.

[5] LYU Sidan,CHEN Weiping,WEN Xuefa,etal.Integration of HYDRUS-1D and MODFLOW for evaluating the dynamics of salts and nitrogen in groundwater under long-term reclaimed water irrigation[J].Irrigation Science,2019,37(1):35-47.

[6] HUSSEIN I A,LIU Junmin,ZHANG Yinqin,etal.Numerical simulation of groundwater of large-scale irrigation in northern China:case of Baoyang Irrigation Area,Shaanxi,China[J].International Journal of Water Resources and Environmental Engineering,2013,5(9):533-545.

[7] BEHERA S K,PANDA R K.Assessing soil and groundwater contamination with HYDRUS-1D:a study from West Bengal[J].Environmental Quality Management,2011,20(3):59-75.

[8] PONTEDERIO E M,HEILBRON P F,PEREZ-GUERRERO J,etal.Reassessment of the Goiania radioactive waste repository in Brazil using HYDRUS-1D[J].Nephron Clinical Practice,2018,66(2):202-210.

[9] ADHIKARI K,PAL S,CHAKRABORTY B,etal.Assessment of phenol infiltration resilience in soil media by HYDRUS-1D transport model for a waste discharge site[J].Environmental Monitoring and Assessment,2014,186(10):6417-6432.

[10] 鐘秀,馬騰,劉林,等.地下水飲用水源地污染源風(fēng)險(xiǎn)等級(jí)評(píng)價(jià)方法研究[J].安全與環(huán)境工程,2014,21(2):104-108.

ZHONG Xiu,MA Teng,LIU Lin,etal.Method research on the risk level evaluation of contamination sources for drinking groundwater wellhead[J].Safety and Environmental Engineering,2014,21(2):104-108.

[11] 屈飛行,湯明高,王自高,等.地質(zhì)災(zāi)害危險(xiǎn)性區(qū)劃中不同權(quán)重計(jì)算方法的對(duì)比分析[J].水利水電技術(shù),2016,47(11):136-140.

QU Feihang,TANG Minggao,WANG Zigao,etal.Comparative analysis on different methods for weight calculation in geohazard risk zoning[J].Water Resources and Hydropower Engineering,2016,47(11):136-140.

[12] LU Peng.Structural effects of participation propensity in online collective actions:based on big data and Delphi methods[J].Journal of Computational and Applied Mathematics,2018,344:288-300.

[13] XUE Rui,WANG Chong,LIU Mengli,etal.A new method for soil health assessment based on analytic hierarchy process and meta-analysis[J].Science of The Total Environment,2018,650:2771-2777.

[14] ABESSI O,SAEEDI M.Hazardous waste landfill siting using GIS technique and analytical hierarchy process[J].Environment Asia,2010,3(2):69-78.

[15] AHYAPRATAMA A,SARNO R.Application of Analytic Hierarchy Process (AHP) and Simple Additive Weighting (SAW) cethods in singer selection process[J].International Conference on Information and Communications Technology,2018,9:234-239.

[16] YUAN Zhong,ZHANG Xianyong,FENG Shan.Hybrid data-driven outlier detection based on neighborhood information entropy and its developmental measures[J].Expert Systems with Applications,2018,112:243-257.

[17] PAN Guangbo,XU Youpeng,YU Zhihui,etal.Analysis of river health variation under the background of urbanization based on entropy weight and matter-element model:a case study in Huzhou City in the Yangtze River Delta,China[J].Environmental Research,2015,139:31-35.

[18] XU Linyu,LIU Guiyou.The study of a method of regional environmental risk assessment[J].Journal of Environmental Management,2009,90:3290-3296.

[19] SAIKIA J,ROY S,BORDOLOI M,etal.Atmospheric aerosols around three different types of coal-based industries:emission parameters,cytotoxicity assay,and principal component analysis[J].Journal of Aerosol Science,2018,126:21-32.

[20] SANIKHANIA H,DEO C R,YASEEN M Z,etal.Non-tuned data intelligent model for soil temperature estimation:a new approach[J].Geoderma,2018,330:52-64.

[21] BAI Libiao,WANG Hailing,HUANG Ning,etal.An environmental management maturity model of construction programs using the AHP-Entropy approach[J].International Journal of Environmental Research & Public Health,2018,15(7):1317.

[22] 李艷萍,喬琦,柴發(fā)合,等.基于層次分析法的工業(yè)園區(qū)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)權(quán)重分析[J].環(huán)境科學(xué)研究,2014,27(3):334-340.

LI Yanping,QIAO Qi,CHAI Fahe,etal.Study on environmental risk assessment index weight of industrial park based on the analytic hierarchy process[J].Research of Environmental Sciences,2014,27(3):334-340.

[23] 李帥,魏虹,倪細(xì)爐,等.基于層次分析法和熵權(quán)法的寧夏城市人居環(huán)境質(zhì)量評(píng)價(jià)[J].應(yīng)用生態(tài)學(xué)報(bào),2014,25(9):2700-2708.

LI Shuai,WEI Hong,NI Xilu,etal.Evaluation of urban human settlement quality in Ningxia based on AHP and the entropy method[J].Chinese Journal of Applied Ecology,2014,25(9):2700-2708.

[24] GHIMIRE P L,KIM Y.An analysis on barriers to renewable energy development in the context of Nepal using AHP[J].Renewable Energy,2018,129:446-456.

[25] 楊昱,廉新穎,馬志飛,等.再生水回灌地下水環(huán)境安全風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)方法研究[J].生態(tài)環(huán)境學(xué)報(bào),2014,23(11):1806-1813.

YANG Yu,LIAN Xinying,MA Zhifei,etal.Risk assessment technology method on groundwater environment safety of reclaimed water injection[J].Ecology and Environmental Sciences,2014,23(11):1806-1813.

[26] BAHREINI F,SOLTANIAN A R.Identification of a gene set associated with colorectal cancer in microarray data using the entropy method[J].Cell Journal,2019,20(4):569-575.

[27] XU Hongshi,MAChao,LIAN Jijian,etal.Urban flooding risk assessment based on an integrated k-means cluster algorithm and improved entropy weight method in the region of Haikou,China[J].Journal of Hydrology,2018,563:975-986.

[28] YANG Yilin,LEI Xiaohui,LONG Yan,etal.A novel comprehensive risk assessment method for sudden water accidents in the Middle Route of the South-North Water Transfer Project (China)[J].Science of the Total Environment,2020,698:134167.

[29] 于可利,劉華峰,李金惠,等.危險(xiǎn)廢物填埋設(shè)施的環(huán)境風(fēng)險(xiǎn)分析[J].環(huán)境科學(xué)研究,2005,18(S1):43-47.

YU Keli,LIU Huafeng,LI Jinhui,etal.The environmental risk assessment of hazardous waste landfill[J].Research of Environmental Sciences,2005,18(S1):43-47.

[30] VRIJHEID M.Health effects of residence near hazardous waste landfill sites:a review of epidemiologic literature [J].Environmental Health Perspectives,2000,108:101-112.

[32] LI Ying,LI Jinhui,CHEN Shusheng,etal.Establishing indices for groundwater contamination risk assessment in the vicinity of hazardous waste landfills in China[J].Environmental Pollution,2012,165:77-90.

[33] 楊昱,廉新穎,馬志飛,等.污染場(chǎng)地地下水污染風(fēng)險(xiǎn)分級(jí)技術(shù)方法研究[J].環(huán)境工程技術(shù)學(xué)報(bào),2017,7(3):323-331.

YANG Yu,LIAN Xinying,MA Zhifei,etal.Risk ranking technology method on groundwater pollution of contaminated sites[J].Journal of Environmental Engineering Technology,2017,7(3):323-331.

[34] 李娟,呂寧磬,楊洋,等.典型場(chǎng)地地下水污染源強(qiáng)分級(jí)評(píng)價(jià)方法研究[J].環(huán)境工程學(xué)報(bào),2014,8(11):4726-4736.

LI Juan,Lü Ningqing,YANG Yang,etal.Study on groundwater pollution source intensity rating assessmentmethod of typical contaminated sites[J].Chinese Journal of Environmental Engineering,2014,8(11):4726-4736.