西藏地区由于近几年经济的发展,生活水平的逐渐提高以及旅游业的兴起,产生大量的生活垃圾和旅游垃圾[1]。对于城市生活垃圾的处理,在西藏地区主要还是以卫生填埋为主[2]。而拉萨市垃圾填埋场使用接近20年,随着使用年限的增加,垃圾填埋场的防渗层可能出现一些安全隐患[3]。垃圾填埋场在使用中会产生大量的气体和渗滤液,如果防渗措施功能不佳,则会对周围的大气、地下水、地表水和土壤等自然环境有所污染[4]。则防渗层的损坏会造成污染物的泄漏,这些污染物中包含大量的重金属。如果地下水中出现铅的泄漏,会造成铅的超标和引起多种疾病,例如肾损伤、不孕、流产、铅脑病、腹绞痛、溶血性贫血等[5,6]。而目前主要处理含铅的地下水方法有物理屏蔽法、抽出处理法和原位修复法[7-11]。
文中采用AHP(层次分析法)和TOPSIS法(逼近理想解排序法/优劣解距离法)相结合分析。其中,AHP决策分析法是由美国运筹学家Saaty提出的一种定性与定量相结合的决策分析方法[12]。TOPSIS法是由Hwang等[13]提出,根据有限个评价对象与理想化目标的接近程度进行排序的方法,是多目标决策中经常使用的方法。例如,张婧等[14]运用AHP-TOPSIS方法选取填埋场区域氨氮污染地下水的修复方案,结果表明,空气注入+原位修复技术更适合地下水修复。王纪洋等[15]用AHP-TOPSIS模型评价了危险化工工艺的风险等级。张士宽等[16]运用AHP和MCDA进行地下水污染修复技术优选排序,表明抽出处理技术、原位微生物修复技术较为理想。霍攀等[17]用AHP和模糊评判法对垃圾填埋场选址的应用进行对比,表明2种方法结合可提高评判结果的准确度。吴军年等[18]也运用AHP法选取废渣库的选址确定环境经济最适宜场址。陈海滨等[19]用AHP确定权重和节约法对收运路线进行优化,从而选取村镇生活垃圾收运路线。李丹丹等[20]以瓮安县农产品基地农用地土壤作为研究对象,用GIS、地统计分析法和组合赋权TOPSIS模型法进行了综合评价。赵国存等[21]用AHP-TOPSIS方法评价装备保障信息,提到该模型的优点在于不损失指标初始信息,使评价结果与实际非常接近。
本文分析拉萨市垃圾填埋场地下水铅的修复技术,是基于拉萨市的自然环境和社会环境来选取建设条件、技术条件、经济条件、环境条件、污染特征及水化学特征6个方面作为评价指标体系。先利用AHP确定权重,然后用TOPSIS法对指标体系的优劣进行排序分析。以AHP和TOPSIS相结合的方法完成综合分析,最后通过评价结果确定了最适宜的修复方案。
AHP主要用于分析多目标、多要素和多层次的复杂决策问题,将问题分解为若干具有条理性的层次和因素,建立层次分析的模型,然后比较每一层次的要素,得到相对重要的标度建立矩阵,从而计算出最大特征值和特征向量,即可得出不同方案的权重。AHP决策分析的主要步骤如下:
1)分析问题。了解问题所涉及的范围和因素关系,掌握充分的实际信息。
2)建立层次结构模型。对问题的要素进行分组分层次分析:目标层—准则层—措施层排列,建立层次结构模型。
3)建立判断矩阵。该步骤中评定层次中元素间的相对重要程度判断,进行两两比较。元素间的相对重要性的判断值(表1)采用1~9标度法。
表1 相对重要性判断值
Table 1 The judgement value criteria of relative significance
同等重要略微重要重要更重要极重要相邻中值135792,4,6,8
4)层次单排序。对上层次中某元素而言,确定本层次与之间有联系的各元素重要性次序的权重值,是层次总排序的基础[12]。这一层次计算矩阵的特征根和特征向量,计算出最大特征根,再计算权重值。为了使矩阵合理化,需要进行一致性检验。计算出的一致性指标要与表2中RI平均一致性指标比较。其中的比值CR若<0.1,则矩阵具有一致性;若>0.1,则矩阵需要进行调整满足一致性。
表2 RI平均随机一致性指标
Table 1 RI average random consistency index
阶数123456789RI000.580.901.121.241.321.411.45
5)层次总排序。计算所有元素的权重值,对每个层次进行逐层排序。判断层次总排序是否具有一致性,也需要进行一致性检验。
而AHP分析方法主要是计算最大特征值和特征向量,有2种方法(方根法、和积法),这里主要采用和积法[12,14,15,19]。
TOPSIS是在归一化后的原始数据矩阵中,找出方案中最优方案和最劣方案,然后计算各评价对象与各方案距离,获得评价对象与最优方案的相对接近程度[14,15, 21],其步骤如下:
层次总排序与目标层矩阵X进行规范化,得到标准综合矩阵R:
(1)
加权规范化综合矩阵V:属性的权重向量ω和规范化矩阵R结合:
V=Rω
(2)
根据V,选择指标得分最多的为正理想解V+,得分少的为负理想解V-,然后计算方案和正负理想解的距离D。
(3)
(4)
确定相对接近度C,进行排序:
(5)
Ci按大小排序,值最大者最优的方法。
根据拉萨市垃圾填埋场项目特点和该项目所在地的地质环境状况,在该填埋场周围6个监测井中选取4个没有被异物堵住的监测井进行监测。按照HJ/T164—2004《地下水环境监测技术规范》和GB/T 14848—2017《地下水质量标准》,确定主要监测的指标有pH值、电导率、溶解氧、高锰酸盐指数、氨氮、化学需氧量、石油类、挥发酚、铬(六价)、氰化物、阴离子表面活性剂、砷、汞、硒、镉、铅、溶解性总固体、总硬度、硝酸盐、亚硝酸盐、氟化物、铜、锌、氯化物,水质检测结果如表3所示。依照GB/T 14848—2017 Ⅲ类标准值计算污染物超标率,可以看出重金属中仅有铅超标。
表3 地下水质量监测数据统计
Table 3 Statistics of groundwater quality monitoring data mg/L(除pH外)
监测项目样本量检出率/%2号3号5号6号MAXMINAVERAGEMAXMINAVERAGEMAXMINAVERAGEMAXMINAVERAGE超标率/%标准值Ⅲ类pH值12100.008.107.747.969.488.158.638.787.838.228.47.988.2416.676.5~8.5电导率12100.008358733781372320.6782672326.3385570333.330无溶解氧12100.006.926.076.367.045.116.057.385.976.487.225.976.410无高锰酸盐指数12100.001.700.501.033.201.302.102.401.001.533.200.901.7716.67≤3.0氨氮120.00————————————0≤0.5化学需氧量120.00————————————0≤3.0石油类120.00————————————0无挥发酚120.00————————————0≤0.002铬(六价)1241.670.0070.0070.0070.0180.0150.017———0.0060.0060.0060≤0.05氰化物120.00————————————0≤0.05阴离子表面活性剂1233.330.170.170.170.170.170.170.160.160.160.150.150.150≤0.3砷1216.670.00150.00150.00150.00580.00580.0058——————0≤0.01汞120.00————————————0≤0.001硒120.00————————————0≤0.01镉128.330.00020.00020.0002—————————0≤0.005铅1233.330.0160.0160.0160.0060.0060.0060.0040.0040.0040.0040.0040.0048.33≤0.01溶解性总固体12100.00246104164.67276183242.6728888171.3326668158.670≤1000总硬度12100.00172127146.3319332.10108.2318834.40106.631661261400≤450硝酸盐(以N计)1291.677.081.203.341.400.280.841.000.170.592.271.301.790≤20.0亚硝酸盐(以N计)120.00————————————0≤1.0氟化物1258.330.480.420.451.560.310.940.960.260.610.250.250.258.33≤1.0铜120.00————————————0≤1.0锌120.00————————————0≤1.0氯化物1233.336.996.996.9981.5081.5081.504.994.994.9935.5035.5035.500≤250
注:“—”表示监测结果低于方法检出限和无监测值。
地下水铅污染修复技术主要由物理屏蔽、抽出处理和原位修复技术。其中,物理屏蔽法主要通过隔离受污染的水体,来减少周围水体污染,缺点是处理范围小和地下水污染初期治理[6]。抽出处理是目前最为普遍的方法,将受污染地下水抽出,然后对其净化处理,该方法可以防治污染水体向周围迁移,但是存在浪费大量水资源和能量的缺点[6, 22]。而原位修复技术是目前主要的发展方向,包括渗透反应格栅、原位生物修复、动电修复技术,这3种技术分别具有无需外加动力、节省空间、经济便捷,成本低、不破坏生态平衡,人工需求少、经济高效等优点[6]。水化学特征中,pH会影响化学和生物技术修复过程,会对污染物产生影响[23]。因此,综合考虑主要从场地建设条件、技术条件、经济条件、环境条件、污染特征及水化学特征6个方面选取了17 指标,建立指标体系,如图1所示。自上而下建立了目标层(A)、准则层(B)、子准则层(S)、措施层(P)。
图1 修复技术决策分析层次模型
Figure 1 Hierarchical model for decision analysis of the repair technologies
根据表1中1~9标度法的判断值,对同一层次中的因素进行两两比较,构造矩阵。因此,构造目标层A对准则层B的判断矩阵A-B,见表4。计算出向量AW=[1.473, 1.496, 0.846, 0.578, 1.625, 0.196]T,CI=0.055,而RI参选表2选取1.24,得出CR=0.044<0.1,说明此矩阵具有可接受的一致性。由表4可看出:准则层中污染特征、技术条件、建设条件是修复过程中主要考虑的因素,3个条件权重分别为0.266、0.249、0.235。但在修复过程中经济条件也要考虑,尤其是修复中的设备投资、成本运行等因素,其权重为0.129,其重要性处于第4位。在修复过程中也要考虑对周边环境的影响,其权重为0.09,其重要性位于第5位。最后在修复过程中要特别考虑到水化学特征,因为水中的pH对修复过程中的化学和生物技术会产生一定的影响。张婧等[14]在通过AHP-TOPSIS方法对填埋场地下水污染氨氮的空气注入+原位生物修复和P&T+生物法两种修复方法的比选中,也认为技术指标和场地条件要比经济指标和环境指标更为重要。张伯强等[24]在采用AHP法和MCDA法对填埋场地下水修复技术优选中表明技术指标权重要比经济指标和社会环境指标权重大[16]。张伯强等[24]用MCDA选取沙漠地区污染地下水修复技术中也得出技术指标权重要比经济指标和社会环境指标权重大。
表4 目标层对准则层判断矩阵
Table 4 Judgment matrix of target layer to criteria layer
AB1B2B3B4B5B6WB11132150.235B21123170.249B31/31/2131/350.129B41/21/31/311/350.090B51133170.266B61/51/71/51/51/710.031
通过公式计算出:准则层B对子准则层S的判断矩阵B1-S1-2、B2-S3-7、B3-S8-11、B4-S12-14、B5-S15-16、B6-S17,如表5所示。可知:对修复影响最大的是污染物类型,其权重为0.228。其次,地质条件在修复过程中也有重要的影响,其权重为0.157。修复过程中的技术成熟度也相当重要,其权重为0.088;修复过程中技术可实施性也存在影响,其权重为0.079;修复中的水文条件也是一个因素,其权重均为0.078;在修复过程中设备投资也有一定影响作用,其权重为0.054;修复中也要注意水资源的保护,其权重为0.045;在修复中,修复周期也有影响,其权重为0.041。污染物浓度、后期管理费用、pH、生态保护也起到一定作用,其权重分别为0.038、0.035、0.031、0.030。当然,在修复过程也要注意运行成本、修复效率、工人健康安全、监测费用、副产品危害等问题,其权重分别为0.025、0.024、0.017、0.016、0.015。
构造次准则层对措施方案层的判断矩阵,计算出目标层对措施层的优化矩阵,见表6。方案层对次准则层的矩阵为3阶,RI=0.58。通过AHP的计算,总权重0.426(原位修复法)>0.293(抽出处理法)>0.282(物理屏蔽法),因此原位修复法适合该地地下水铅的修复技术。有学者研究填埋场地下水污染修复方案的比选中,也显示原位修复类的方法优于其他方法[14,16]。
通过式(1)得到R,再通过式(2),将R与权重ω相结合得到V。在V中选择各措施方案的得分最多的正理想解V+和得分最少的负理想解V-,如下所示:
表5 层次总排序
Table 5 Hierarchical total sorting for APH method
AB1B2B3B4B5B60.2350.2490.1290.0900.2660.031W排序S10.3330.0785S20.6670.1572S30.3520.0883S40.1650.0418S50.0980.02414S60.3180.0794S70.0670.01715S80.1210.01616S90.2690.03510S100.4200.0546S110.1900.02513S120.1670.01517S130.3330.03012S140.5000.0457S150.8570.2281S160.1430.0389S1710.03111合计1CI00.0530.0380000.018RI01.120.90.58000.447CR00.0470.0420000.041
表6 目标层对方案层的优化矩阵
Table 6 Optimal matrix of target layer to scheme layer
S1S2S3S4S5S6S7S8S9S10S11S12S13S14S15S16S17W总P10.5940.5400.1570.1050.5280.2490.1740.1630.1630.5280.1400.4930.1250.1090.1400.1400.1140.282P20.1570.2970.2490.2580.3330.1570.1920.2970.2970.3330.3330.1960.3750.3090.3330.5280.4050.293P30.2490.1630.5940.6370.1400.5940.6340.5400.5400.1400.5280.3110.5000.5820.5280.3330.4810.426CI0.0270.0050.0270.0190.0270.0270.0050.0050.0050.0270.0270.0270.0000.0020.0270.0270.015RI0.58,(n=3)CR0.0460.0080.0460.0330.0460.0460.0080.0080.0080.0460.0460.0460.0000.0030.0460.0460.025
所示:
V+=[0.070,0.133,0.078,0.037,0.020,0.071,
0.015,0.013,0.029,0.045,0.020,0.012,0.023,
0.039,0.188,0.031,0.023]
V-=[0.019,0.040,0.021,0.006,0.005,0.019,
0.004,0.004,0.009,0.012,0.005,0.005,0.006,
0.007,0.050,0.008,0.006]
再由式(3)—(5)计算各措施方案的D+、D-和相对接近度Ci,如表7所示。
表7 各方案措施的距离计算和相对接近度
Table 7 Distance calculation results and relative proximity of repairing schemes
方案D+D-Ci物理屏蔽法0.0280.0130.311抽出处理法0.0180.0080.301原位修复法0.0120.0290.713
根据相对接近度确定修复方案优劣排序,由表7可看出:相对接近度Ci(原位修复法)>Ci(物理屏蔽法)>Ci(抽出处理法),所以原位修复法在该地区更适合地下水铅的修复。目标层对方案层的分析优化矩阵也表明,在技术条件中的技术成熟度、修复周期、技术可实施性、工人健康安全方面,原位修复技术比抽出处理法和物理屏蔽法更适宜。通过最后指标权重的优劣排序,也可以看出原位修复技术更具有优势,比较适合该地区地下水铅的修复。
1)对拉萨垃圾填埋场地下水污染物的调查,垃圾填埋场监测井中出现污染物铅,然后选取铅污染修复技术,包括物理屏蔽法、抽出处理法和原位修复法3种方法。根据拉萨市的生态环境、社会环境等方面,从建设条件、技术条件、经济条件、环境条件、污染特征、水化学特征6个方面确定了影响填埋场地下水修复的17个指标,构建了层次分析评价指标体系。
2)通过AHP层次分析法进行权重值确定,认为在地下水污染修复过程中主要考虑的因素有污染特征、技术条件、建设条件。其中,选取的指标中污染物类型、地质条件、技术成熟度等作为权重值较大,对污染物修复技术的选取也有一定影响。最后通过优化矩阵,得出原位修复技术优势较明显,适合该地区地下水修复。
3)在采用TOPSIS方法进行优劣排序对各修复技术进行综合分析后,原位修复技术相对接近度值最大,最终确定原位修复法适合该地垃圾填埋场铅污染的修复。
[1] 周文武,陈冠益,周鹏,等.西藏拉萨市区餐厨垃圾现状分析及建议[J].环境与可持续发展,2019,44(1):60-63.
[2] 拉萨市人民政府.2016年拉萨市固体废物管理及拉萨市固体废物污染防治信息公告[EB/OL].[2017-06/20],http://www.lasa.gov.cn.
[3] 中国膜工业协会.拉萨市生活垃圾填埋场建成膜法渗滤液处理站效果良好[EB/OL].[2015-06-16],http://www.membranes.com.cn.
[4] 薛红琴, 速宝玉, 盛金昌,等. 垃圾填埋场渗滤液的防渗措施和地下水的污染防护[J]. 安全与环境学报, 2002, 2(4):18-22.
[5] 袁宝珊,吴宜群.环境铅污染与儿童健康[J].国外医学(卫生学分册),1998(4):3-8.
[6] 姜楠, 王鹤立, 廉新颖. 地下水铅污染修复技术应用与研究进展[J]. 环境科学与技术, 2008, 31(2):56-60.
[7] COULTHARD T L, SAMIA FADL. Adsorption of Water Pollutants by a Coal Sorption Process [R]. For presentation at the CSAE Annual Meeting, 1973,21(8).
[8] VIDIC R D, POHLAND F G. Treatment Walls Technology evaluation report TE-96-01 [R]. Pittsburgh P A: Groundwater Remediation Technologies Analysis Center, 1996.
[9] LEUSCH A, HOLAN Z R, VOLESKY B. Solution and particle effects on the biosorption of heavy metals by seaweed biomass[J]. Applied Biochemistry and Biotechnology, 1997, 61(3):231-249.
[10] RAHMUNI G N, H STEMBERG S P K. Bioremoval of lead from water using Lemna minor[J]. Bioresour Technol, 1999, 70(3): 225-230.
[11] 胡宏韬, 侯玲娟. 铅污染地下水电动修复研究[J]. 水处理技术, 2009, 35(9):56-59.
[12] 徐建华. 计量地理学[M]. 北京:等教育出版社, 2006.
[13] HWANG C L,YOON K S. Multiple Attribute Decision Making [M]. Berlin: Spring-Verlag,1981.
[14] 张婧, 皮鎏, 崔佳鑫, 等. 垃圾填埋场区域氨氮污染地下水的修复方案比选[J]. 环境保护科学, 2017,43(3):125-131.
[15] 王纪洋,张明广,王雪栋.基于AHP-TOPSIS模型的危险化工工艺风险等级评价研究[J].安全与环境工程,2016,23(6):100-105.
[16] 张士宽,王月,安达,等. 垃圾填埋场地下水污染修复技术优选研究[J].环境工程技术学报,2017,7(4):463-469.
[17] 霍攀,曹丽文,田艳凤. AHP与模糊评判法在垃圾填埋场选址中的应用对比[J].环境工程,2015,33(3):131-135.
[18] 吴军年,张后辉,邢景敏. 基于层次分析法的危废渣库选址的环境经济适宜性分析[J].环境工程,2011,29(2):99-103.
[19] 陈海滨,李文杰,杨龑,等. 基于节约法的村镇生活垃圾收运路线优化研究[J].环境工程,2016,34(12):122-125.
[20] 李丹丹,周忠发,但雨生,等. 基于组合赋权TOPSIS模型的土壤养分空间分析及综合评价:以瓮安县为例[J].环境工程,2018,36(8):183-188.
[21] 赵国存,刘占岭. 基于AHP-TOPSIS的装备保障信息定量评价研究[J].装备环境工程,2012,9(4):118-122.
[22] 李君, 常莉. 我国城市地下水污染状况与治理对策[J]. 开封大学学报, 2006, 20(4):89-91.
[23] 鄂佳楠, 周睿, 郑龙日,等. 基于蒙特卡罗法和层次分析法的污染场地地下水修复技术筛选方法研究[J]. 环境污染与防治, 2017, 39(5):499-503.
[24] 张伯强,安达,王月,等.基于MCDA的沙漠地区污染场地地下水修复技术优化方法[J].环境工程学报,2016,10(10):5521-5527.