基于层次分析法-熵权法的污...法——以农药类污染地块为例_胡星路.pdf
《基于层次分析法-熵权法的污...法——以农药类污染地块为例_胡星路.pdf》由会员分享,可在线阅读,更多相关《基于层次分析法-熵权法的污...法——以农药类污染地块为例_胡星路.pdf(14页珍藏版)》请在咨信网上搜索。
1、文章栏目:土壤污染防治DOI10.12030/j.cjee.202302002中图分类号X53文献标识码A胡星路,熊杰,张超艳,等.基于层次分析法-熵权法的污染地块异味污染物优先度评价方法以农药类污染地块为例J.环境工程学报,2023,17(6):1923-1936.HUXinglu,XIONGJie,ZHANGChaoyan,etal.Priorityevaluationmethodofodorpollutantsincontaminatedsitesbasedonanalytichierarchyprocessandentropyweightmethod:Acasestudyofpesti
2、cidecontaminatedsitesJ.ChineseJournalofEnvironmentalEngineering,2023,17(6):1923-1936.基于层次分析法-熵权法的污染地块异味污染物优先度评价方法以农药类污染地块为例胡星路1,2,熊杰2,张超艳2,郭晓欣2,周友亚2,李海明1,1.天津科技大学海洋与环境学院,天津300457;2.生态环境部土壤与农业农村生态环境监管技术中心,北京100012摘要土壤中异味污染物的风险管控已成为我国农药化工等行业污染地块环境管理的迫切需求。为筛查土壤优先控制异味污染物,研究基于土壤中挥发性异味污染物的迁移、暴露过程和危害效应的系统量
3、化与分析,建立了污染地块异味污染物优先度评价指标体系;采用层次分析法-熵权法确定指标综合评价权重,分析了异味污染物优先度的影响因素,并选取 3 个农药类案例地块进行验证性评价。结果表明,构成地块异味污染物优先度排序指标体系第二层次的 4 个准则中,影响程度从高到低为:危害性、迁移性、防护性与暴露性;第三层次15 个要素主客观综合权重范围为 0.03360.1496,其中嗅阈值权重最高。案例地块异味污染物基于层次分析法-熵权法的优先度评价等级与基于实测数据通过异味活度值和健康风险影响评价的等级具有一致性。本研究结果可为异味污染地块优先控制异味污染物的筛查提供参考。关键词污染地块;异味污染物;层次
4、分析法;熵权法;优先度评价农药化工类行业生产原辅料及产物中含有大量异味污染物,因跑冒滴漏和突发事故可能对土壤与地下水造成污染1。现阶段我国污染地块主要采用基于健康风险的管理思路,暂未对异味污染物提出明确的风险管控要求2。实际地块修复与再开发过程中,施工扰动可能造成土壤中异味污染物挥发扩散,并对周边环境和人群造成污染与影响3-5。因此,开展异味污染物风险管控已成为污染地块环境管理的迫切需求。异味污染物种类众多、且对人体嗅觉刺激不一6,提出地块优先控制异味污染物筛查和评价方法,制定土壤异味污染物优先控制名录,是指导异味污染物调查与评估工作、有效避免和消除异味污染的重要前提。国内外提出的大气和水环境
5、优先污染物名录确定过程一般为筛选待评估污染物清单、对清单中的污染物开展环境风险评估并进行优先度排序、确定应优先管理的污染物7。优先度排序的具体做法是基于污染物的危害性、暴露性、风险性等因素选择评价指标8,建立优先污染物的评价指标体系,进而采用主观赋权法(如层次分析法、专家打分法等)9-11、客观赋权法(如熵权法、模糊集理论、主成分分析法等)12-13和综合赋权法14-17等确定评价指标的权重并进行打分收稿日期:2023-02-01;录用日期:2023-04-04基金项目:国家自然科学基金资助项目(41977377);国家重点研发计划土壤专项资助项目(2018YFC1803001)第一作者:胡星
6、路(1997),男,硕士研究生,;通信作者:李海明(1971),男,博士,教授,环境工程学报Chinese Journal ofEnvironmental Engineering第 17 卷 第 6 期 2023 年 6 月Vol.17,No.6Jun.2023http:/E-mail:(010)62941074排序。其具体指标的选取主要考虑区域或行业污染物产生、排放与暴露的总体特征,而较少考虑地块水文地质条件、企业环境管理水平等方面的个体差异。孟洁等18和马妍等19采用异味活度值与异味贡献率筛查了土壤中的主要异味污染物,具体做法是借助土壤气样品的全扫描结果,确定异味污染物清单;利用近土壤环境
7、空气的各污染物质量浓度,计算其异味活度值并进行排序;计算总异味活度值和各污染物的异味贡献率,最终以累计贡献率超过 90%为依据确定主要异味污染物。但这一方法需借助大量现场监测数据,所需时间长、成本高,且异味污染物筛查滞后于调查采样工作,无法有效指导现场调查工作。本研究旨在综合考虑地块土壤特性、水文地质条件及企业环境管理水平等地块个性化要素,借助层次分析法与熵权法,考虑主客观因素,对场地污染土壤中异味污染物的危害程度、暴露水平、在各介质中的迁移分配以及地层防护能力等相关要素进行系统量化、分析和评价,拟建立地块尺度的异味污染物优先级排序评价方法,并在 3 个农药异味地块开展实证研究,以期为地块优先
8、控制异味污染物的筛查和名录建立提供参考。1材料与方法1.1指标体系构建地块异味污染物优先度评价指标体系的构建可通过 3 个层次实现:第一层次为目标层(A),即预定的目标或理想结果;第二层次为准则层(B),包含了为实现目标所涉及的中间环节,包括所需考虑的准则和子准则;第三层次为要素层(C),即为实现目标所考虑准则下的具体要素。农药行业地块土壤中异味污染物对暴露人群的健康和嗅觉器官造成影响需经历以下过程:1)生产过程泄漏等造成污染物下渗并污染土壤20-21;2)土壤中污染物解吸挥发至气相,在非饱和带土壤中向上迁移;3)气态异味污染物从地表进入大气,在呼吸区混合迁移,被人群吸入22。基于上述污染物产
9、生、迁移和暴露过程,设置危害性、暴露性、迁移性和防护性 4 个准则评价异味污染物的优先度,分别表征污染物的异味潜在影响与健康毒性、污染物潜在暴露量、污染物迁移能力和地块防护土壤污染的能力。每一准则可细分为 35 项主要评价指标,评价指标体系如图 1 所示。1.2指标分级与赋值将上述每个评价指标分为 5 级,按影响程度分别赋值 15 分,其中 1 为无影响、2 为低度影响、3 为中度影响、4 为高度影响、5 为极高影响23-30,对于数据不足的指标赋分 3 分。指标的分级与赋值方法如表 1 所示。1.3评价指标权重采用层次分析法与熵权法相耦合的综合赋权法(层次分析法-熵权法)确定评价指标权重值,
10、其优点是既考虑了基于实际经验的人为主观判断排序,又考虑了基于真实数据的客观权重,减少了主客观因素带来的误差,使评价结果更符合客观实际31-32。WjWj1)指标主观权重。采用层次分析法计算指标的主观权重,可基于人为主观判断各评价指标间的相对重要性排序,步骤包括构建层次结构模型、建立判断矩阵、一致性检验和权重计算33-36。根据层次分析法的应用准则,以及目标层、准则层与要素层的相对关系,建立污染地块异味污染物优先度评价的递阶层次结构模型(图 1)。采用 19 标度理论,根据评价指标之间的重要程度进行评价指标两两比较,构造各层级间的判断矩阵。根据 SAATY37提出的 RI 值表,计算判断矩阵的一
11、致性比例 CR,包括层次单排序一致性检验与层次总排序一致性检验,当 CR0.1 则认为判断矩阵满足一致性条件,说明权重分配合理。将准则层权重值与要素层权重值相乘,得到指标的主观权重值。CR、CI 与的计算如式(1)式(3)所示。CR=CIRI(1)1924环境工程学报第17卷图1污染地块异味污染物优先度评价指标体系Fig.1Priorityevaluationindexsystemofodorpollutantsincontaminatedsites表1评价指标的分级与赋值Table1Evaluationindexclassificationandvalueassignment评价指标赋值分级
12、依据参考文献准则层要素层单位12345危害性嗅阈值mgm31021102102110410210423致癌斜率因子(mgkg1d1)1524非致癌参考浓度mgm3124暴露性年使用量(产量)t10425生产时间a5024企业环境管理水平无量纲好较好中等较差差27迁移性饱和蒸汽压kPa1023分配系数cm3g11035001031025001102124亨利常数无量纲124降解性(土壤半衰期)d20026防护性土壤渗透系数cms110228土壤含水率%4030402530202520152010155.010529土壤孔隙度%5029地下水水位m4020401020210230第6期胡星路等:基
13、于层次分析法-熵权法的污染地块异味污染物优先度评价方法以农药类污染地块为例1925式中:CR 为判断矩阵的一致性比例,无量纲;RI 为平均随机一致性指标,无量纲;CI 为判断矩阵一致性指标,无量纲。CI=maxnn1(2)式中:max为判断矩阵最大特征根,无量纲;n 为判断矩阵阶数,无量纲。Wj=wBzwCj(3)WjwBzwCj式中:为要素层 j 指标对于目标层的层次总排序权重值,无量纲;为准则层中 z 准则对于目标层的层次单排序权重值,无量纲;为要素层中 j 指标对于所对应准则层的层次单排序权重值,无量纲。xijXijEjWjXijEjWj2)指标客观权重。采用熵权法计算指标客观权重,可反
14、应评价指标真实数据的贡献,具体步骤包括构建原始数据矩阵、原始数据矩阵标准化、信息熵计算及客观权重值确定38-39。选取土壤污染信息,对 15 个评价指标进行赋值,构成原始数据矩阵。对原始数据矩阵进行标准化处理,构建标准化矩阵。计算各指标的信息熵。计算各指标权重值。、和的计算公式如式(4)式(6)所示。Xij=xijmin(xj)max(xj)min(xj)(4)xijXijmax(xj)min(xj)式中:为第 i 个地块的第 j 项指标的实际取值构成的原始数据矩阵,无量纲;为标准化矩阵,无量纲;为第 j 项指标在 i 个地块中的上限值,无量纲;为第 j 项指标在 i 个地块中的下限值,无量纲
15、。Ej=1ln(n)ni=1(Xijni=1XijlnXijni=1Xij)(5)Ej式中:为信息熵,无量纲;n 为污染地块的数量,无量纲。Wj=1Ejmj=11Ej(6)Wj式中:为基于熵权法计算的客观权重,无量纲;m 为要素层中评价指标个数,无量纲。WjWjRjRjRj3)指标综合权重。综合层次分析法确定的主观权重值与熵权法确定的客观权重值,计算各指标的综合权重值,构建综合权重向量,反应各指标对异味污染物优先度的影响程度,越大表示该指标对异味污染物优先度的影响程度越高。的计算公式如式(7)所示。Rj=Wj+(1)Wj(7)Rj式中:为基于层次分析法-熵权法计算的综合权重值,无量纲;为分布系
16、数,无量纲。在偏差的平方和最小的前提下,主观权重和客观熵权重分别占 50%,主客观权重值对评价指标重要性是一致的40-41,即=0.5。1.4优先度等级将异味污染物评价指标赋值与该指标权重的乘积进行加和,得到污染地块异味污染物优先度评价分值,如式(8)所示。G=RjCj(8)Cj式中:G 为异味污染物优先度评价分值;为异味污染物评价指标赋值。优先度等级根据异味污染物对人体健康与嗅觉产生影响的程度采用 3 级评分制,各等级中具体评价分值区间的对应关系如表 2 所示。根据计算的异味污染物优先度评价分值,按照表 2 对应1926环境工程学报第17卷关系进行综合评价。1.5优先度验证分析方法研究选取山
17、东某农药厂(地块 a)、广州某农药化工厂(地块 b)与苏州某农药化工厂(地块 c)3 个异味污染地块作为验证案例,通过查询案例地块土壤污染状况调查与风险评估报告等资料,获取土壤中异味污染物质量分数实测值。采用异味活度值和健康风险影响度表征土壤中异味污染物对暴露人群的嗅觉感官和健康风险的影响程度,划分并确定案例地块异味污染物的影响等级,并与基于优先度评价结果进行对比验证。1)异味活度值计算方法。采用异味活度值表征土壤中异味污染物对人体嗅觉感官的影响程度,异味活度值表征污染物在空气中的暴露质量浓度与其嗅阈值之间的比值,异味活度值越大说明该污染物异味影响越明显2。污染物异味活度值计算如式(9)所示。
18、OVA=CairCthr(9)式中:OVA 为异味污染物异味活度值,无量纲;Cair为异味污染物空气暴露质量浓度,mgm3;Cthr为异味污染物嗅阈值,mgm3。将污染物室内空气暴露质量浓度与土壤污染物质量分数的比值定义为挥发因子,故 Cair计算如式(10)式(14)所示。Cair=CsVFsubia(10)VFsubia1=1KswH|1+Def fsLBER186400Ls+DeffsLcrackDeffsLs|LBER186400DeffsLs103(11)VFsubia2=dsubbLBER18640031536000103(12)Deffs=Da3.33air2+DwH3.33wa
19、ter2(13)Ksw=water+Koc focb+Hairb(14)式中:Cs为土壤中异味污染物质量分数,mgkg1;VFsubia1为下层土壤中污染物扩散进入室内空气的挥发因子(算法 1),kgm3;VFsubia2为下层土壤中污染物扩散进入室内空气的挥发因子(算法 2),kgm3;VFsubia为土壤中污染物扩散至室内空气的挥发因子(取算法 1 与算法 2 的较小值),kgm3;Dseff为土壤中气态污染物的有效扩散系数 cm2s1;Ksw为土壤-水中污染物分配系数 cm3g1;H 为污染物亨利常数,无量纲;LS为下层污染土壤上表面到地表距离,cm;dsub为下层污染土壤厚度,cm;D
20、a为污染物在空气中的扩散系数,cm2s1;Dw为污染物在水中的扩散系数,cm2s1;Koc为土壤-有机碳分配系数,cm3g1,其余参数定义及取值如表 3 所示。2)健康风险评价方法。健康风险影响度是土壤中污染物质量分数与风险筛选值之间的比值。健康风险影响度越大,说明该污染物产生的健康风险对暴露人群的影响程度越大42。健康风险影表2优先度等级与分值对应关系Table2Correspondencebetweenprioritygradeandscore异味污染物优先度等级评价分值优先度等级对应的污染物描述高优先级3.55.0对人体健康具有高风险且异味强烈难以接受中优先级2.53.5对人体健康具有较
21、高风险且能感觉到明显异味低优先级02.5对人体健康风险较低(或风险可接受)且异味较轻(或无异味)第6期胡星路等:基于层次分析法-熵权法的污染地块异味污染物优先度评价方法以农药类污染地块为例1927响度的计算如式(15)所示。HRV=CsoilCsoil,screeing(15)式中:HRV 为污染物健康风险影响度,无量纲;Csoil为土壤中异味污染物质量分数,mgkg1;Csoil,screeing为异味污染物的风险筛选值,mgkg1,参考土壤环境质量建设用地土壤污染风险管控标准(试行)(GB36600-2018)43中土壤污染风险筛选值和依据建设用地土壤污染状况调查技术导则(HJ25.3-2
22、019)42技术方法推导计算的土壤污染风险筛选值。2结果与讨论2.1指标权重值计算结果1)主观权重计算结果。以准则层 B 为例,通过专家评分结果,构造目标层 A 与准则层 B 中4 种指标的判断矩阵 D1,进行一致性检验与准则层权重值计算,如表 4 所示。经验证,判断矩阵随机一致性比率 CR0.1,因此所构建判断矩阵一致性可接受,主观权重分配合理,主观权重值如表 5 所示。表3计算模型参数取值Table3Valuesofmodelparameters参数定义单位取值取值来源地块a地块b地块cwater非饱和土层土壤中孔隙水体积比无量纲0.350.360.30实测值air非饱和土层土壤中孔隙空气
23、体积比无量纲0.070.080.12实测值foc土壤有机碳质量分数无量纲0.00250.004260.0029实测值b土壤容重kgdm-31.511.521.58实测值非饱和土层土壤中孔隙空气体积比无量纲0.450.080.42实测值LB室内空间体积与气态污染物入渗面积比cm22022030041ER室内空气交换速率次d-112122041地基裂隙表面积占室内地表面积比例无量纲0.0050.0050.00541Lcrack室内地基或墙体厚度cm35353541气态污染物入侵持续时间a30302541表4目标层 A 与准则层 Bi 的判断矩阵 D1Table4JudgmentmatrixD1fo
24、rlayerAandlevelBi目标层危害性迁移性暴露性防护性主观权重值一致性检验一致性是否可接受危害性12330.4548max=4.0098可接受迁移性1/21220.2634CI=0.0038暴露性1/31/2110.1409RI=0.89防护性1/31/2110.1409CR=0.00370.1表5评价指标权重分配Table5Theweightvaluedistributionofevaluationindicators目标层准则层要素层Wj主观权重Wj客观权重Rj综合权重污染地块异味污染物优先度评价危害性嗅阈值0.23640.06290.1496致癌斜率因子0.12580.0792
25、0.1025非致癌参考浓度0.09250.05850.07551928环境工程学报第17卷xij2)客观权重计算结果。选取地块 ac 的土壤污染信息,取其平均值,对 15 个评价指标进行赋值,构成原始数据矩阵。xij=|1.300.140.7120 00049212.708370.301.720.000 328.984.930.231.502.940.075.2120 00026382.71627.620.231.450.000 0225.166.000.445.672.171.110.8830 00050366.96593.660.271.860.000 0424.004.550.464.7
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 层次 分析 熵权法 农药 污染 地块 胡星路
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。