CN112785206A - 一种金属行业涉污企业影响范围模拟与风险评估方法 - Google Patents
一种金属行业涉污企业影响范围模拟与风险评估方法 Download PDFInfo
- Publication number
- CN112785206A CN112785206A CN202110233167.2A CN202110233167A CN112785206A CN 112785206 A CN112785206 A CN 112785206A CN 202110233167 A CN202110233167 A CN 202110233167A CN 112785206 A CN112785206 A CN 112785206A
- Authority
- CN
- China
- Prior art keywords
- pollution
- risk
- adjustment distance
- distance
- related enterprise
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 229910052751 metal Inorganic materials 0.000 title claims abstract description 46
- 239000002184 metal Substances 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000011156 evaluation Methods 0.000 claims abstract description 18
- 238000012502 risk assessment Methods 0.000 claims description 27
- 238000001556 precipitation Methods 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 19
- 229910001385 heavy metal Inorganic materials 0.000 claims description 15
- 238000003723 Smelting Methods 0.000 claims description 12
- 238000005065 mining Methods 0.000 claims description 12
- 239000002689 soil Substances 0.000 claims description 12
- 230000002265 prevention Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000002776 aggregation Effects 0.000 claims description 5
- 238000004220 aggregation Methods 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 231100000331 toxic Toxicity 0.000 claims description 4
- 230000002588 toxic effect Effects 0.000 claims description 4
- 239000003344 environmental pollutant Substances 0.000 abstract description 9
- 231100000719 pollutant Toxicity 0.000 abstract description 9
- 238000009792 diffusion process Methods 0.000 abstract description 6
- 238000012795 verification Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000010606 normalization Methods 0.000 description 5
- 238000011109 contamination Methods 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000012550 audit Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 239000009519 fu-yuan Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000012950 reanalysis Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002912 waste gas Substances 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例公开了一种金属行业涉污企业影响范围模拟与风险评估方法，针对金属类涉污企业，基于金属污染物扩散的特征，通过GIS空间分析方法，模拟涉污企业的影响距离并生成影响范围，开展涉污企业的风险评估，识别出高风险区，并进行精度验证。本发明实施例提供一种金属行业涉污企业影响范围模拟与风险评估方法，以解决区域风险评估现有技术中污染半径设定单一和无法表征真实污染范围的问题。
Description
技术领域
本发明实施例涉及污染模拟评估技术领域，具体涉及一种金属行业涉污企业影响范围模拟与风险评估方法。
背景技术
金属相关行业类涉污企业，如金属尾矿库、金属冶炼厂等，在从事生产活动过程中不合理的污染物堆放、排放等行为，导致有毒有害物质随雨水淋溶、大气沉降等作用向外扩散，造成涉污企业周边一定范围内的土壤、耕地、人口等受到不同程度的污染。该区域范围定义为涉污企业的影响范围。
从宏观层面精准、详细地掌握涉污企业对环境可能造成的风险对工业企业空间布局优化、风险管控、环境保护管理决策制定等具有重要意义。然而，现有的污染风险评估严重依赖于实地采样和人工上报的方式进行，不仅数据获取时间周期长，且人力和财力消耗成本巨大，极大程度上限制了大尺度范围风险评估以及宏观决策的制定。
目前的污染风险评估主要有两种途径，一种是基于单个污染场地的微观评估，另一种是基于行政单元的区域评估。对于第一种，多是针对单个、典型的涉污企业开展的，这样的微观评估主要通过结合实地污染物采样与实验室测量的方式，评估个体涉污企业周围的污染风险情况。然而，这样的小范围风险评估无法体现区域整体的风险情况，不利于宏观上指导产业布局的部署和优化，且如果仅依赖单个涉污企业微观污染风险评估方法获得区域尺度上风险状态，将无疑消耗大量的人力、物力和财力，给政府的财政带来巨大压力。同时，这样的评估往往采用统一的标准半径作为涉污企业对其周边受体产生影响的作用范围，并没有结合涉污企业局地的地形、气候等条件，无法反映真实的污染物影响范围。另一种是基于行政单元的区域尺度风险评估，多以行政单元为空间尺度，如采用地级市尺度统计数据(重工业企业数量、废气排放量、人口等)评估各地级市的风险水平。这种风险评估将行政单元作为均一单元，不能反映风险在行政单元内部的空间差异，更无法表征高风险区的精准空间位置，导致对政府的决策支撑力不足。
发明内容
为此，本发明实施例提供一种金属行业涉污企业影响范围模拟与风险评估方法，以解决区域风险评估现有技术中污染半径设定单一和无法表征真实污染范围的问题。
为了实现上述目的，本发明实施例提供如下技术方案：
根据本发明实施例公开的一种金属行业涉污企业影响范围模拟与风险评估方法，该方法包括以下步骤：
S1、根据涉污企业的行业类型确定所述涉污企业的基本影响距离；
S2、根据涉污企业的经营时长和所处地理位置的多年平均降水量对涉污企业的基本影响距离进行调整，设置涉污企业的污染半径，并以涉污企业为中心绘制圆形缓冲区；
S3、根据涉污企业近十年每月的平均风场环境绘制风向椭圆，共计绘制120个风向椭圆；
S4、将120个风向椭圆与以涉污企业为中心绘制的圆形缓冲区进行合并，合并后再取包络面为涉污企业的影响范围；
S5、采用基于敏感受体的栅格尺度风险评估方法在1km栅格尺度上开展，且认为涉污企业仅对影响范围内的敏感受体有影响，评估只在模拟的影响范围内进行；
S6、建立评估指标体系，所述风险指数为R，影响所述风险指数的第一层指标包括涉污企业的危险性H、路径扩散性P和敏感受体的脆弱性V，第i个像元的风险指数Ri＝Hi×Pi×Vi；
S7、影响危险性H的第二层指标包括核密度K和经营时长DT，核密度代表了涉污企业的聚集程度，第n个涉污企业的危险性Hn＝DTn+Kn；影响路径扩散性P的第二层指标包括土壤质地ST、多年平均降雨量p、多年平均风速WS和地势起伏度RA，第i个像元的路径扩散性Pi＝1/4(STi+pi+WSi+RAi)；影响敏感受体的脆弱性V的第二层指标包括耕地产量；
S8、根据估算出的风险指数R确定高风险区域，并将识别出的高风险区域与国家划定的重金属重点防控区进行对比，验证估算精度。
进一步地，在步骤S1中，所述涉污企业的行业类型包括金属矿采选型和金属冶炼型，所述金属矿采选型企业的基本影响距离为1km，所述金属冶炼型企业的基本影响距离为1.5km。
进一步地，在步骤S2中，金属矿采选型企业基本影响距离受影响因素为：经营时长<15年调整距离0km，经营时长≥15调整距离+0.5km，多年平均降水量<400mm调整距离+1km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km；
金属冶炼型基本影响距离受影响因素为：经营时长<5年调整距离0km，经营时长5-15年调整距离+1km，经营时长>15年调整距离+2km，多年平均降水量<400mm调整距离+0.5km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km。
进一步地，在步骤S3中，每个风向椭圆以涉污企业的空间位置为顶点，月平均风向的下风方向为长轴方向；金属矿采选型企业的风向椭圆长短轴在基本影响距离的基础上调整为：月平均风速<3m/s，短轴长度调整距离0km，月平均风速3-5m/s调整距离+0.5km，月平均风速5-7m/s调整距离+1km，月平均风速>7m/s调整距离+2km，长轴长度在短轴长度基础上+0.5km；金属冶炼型企业的风向椭圆长短轴在基本影响距离的基础上调整为：月平均风速<2m/s调整距离+0.5km，月平均风速2-4m/s调整距离0km，月平均风速>4m/s调整距离-0.5km，长轴长度在短轴长度基础上+1km。
进一步地，敏感受体是指受有毒有害物质影响的受害主体，敏感受体的脆弱性是指受体面对污染的易损性情况，主要的敏感受体包括耕地、土壤、水源或人口。
进一步地，为将点位尺度的危害性指标转换为栅格尺度，以影响范围为插值空间，将危害性H插值为1km的空间数据。
进一步地，在步骤S8中，识别出的高风险区域与国家划定的重金属重点防控区吻合度越高，则识别与评估的估算精度越高。
本发明实施例具有如下优点：
1、基于风场的影响范围空间模拟模型考虑了各个涉污企业所处环境的风场特殊性，模拟了具有场地针对性的影响范围，较以往单一给定影响范围的操作更具有科学性。
2、基于敏感受体的栅格尺度风险评估方法，实现了高风险区空间位置精准识别，可以反映风险在行政单元内部的空间差异。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。
本说明书所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本发明可实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容能涵盖的范围内。
图1为本发明实施例提供的一种属行业涉污企业影响范围模拟与风险评估方法流程图；
图2为本发明实施例1的影响范围模拟示意图；
图3为本发明实施例1的影响范围模拟结果；
图4为本发明实施例1的涉污企业核密度分布示意图；
图5为本发明实施例1的各涉污企业核密度归一化值示意图；
图6为本发明实施例1的涉污企业经营时长归一化结果示意图；
图7为本发明实施例1的危害性评估结果示意图；
图8为本发明实施例1的土壤质地归一化结果示意图；
图9为本发明实施例1的多年平均降水归一化结果示意图；
图10为本发明实施例1的多年平均风速归一化结果示意图；
图11为本发明实施例1的地势起伏度归一化结果示意图；
图12为本发明实施例1的扩散性评估结果示意图；
图13为本发明实施例1的影响范围内的扩散性示意图；
图14为本发明实施例1的耕地产量归一化结果示意图；
图15为本发明实施例1的1km栅格尺度评估结果示意图；
图16为本发明实施例1的县级尺度评估结果示意图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
结合图1所示，本发明具体实施例公开了一种金属行业涉污企业影响范围模拟与风险评估方法，针对金属类涉污企业，基于金属污染物扩散的特征，通过GIS空间分析方法，模拟涉污企业的影响距离并生成影响范围模型，开展涉污企业的风险评估，识别出高风险区，并进行精度验证。该方法具体包括以下步骤：S1、根据涉污企业的行业类型确定所述涉污企业的基本影响距离；S2、根据涉污企业的经营时长和所处地理位置的多年平均降水量对涉污企业的基本影响距离进行调整，设置涉污企业的污染半径，并以涉污企业为中心绘制圆形缓冲区；S3、根据涉污企业近十年每月的平均风场环境绘制风向椭圆，共计绘制120个风向椭圆；S4、将120个风向椭圆与以涉污企业为中心绘制的圆形缓冲区进行合并，合并后再取包络面为涉污企业的影响范围；S5、采用基于敏感受体的栅格尺度风险评估方法在1km栅格尺度上开展，且认为涉污企业仅对影响范围内的敏感受体有影响，评估只在模拟的影响范围内进行；S6、建立评估指标体系，所述风险指数为R，影响所述风险指数的第一层指标包括涉污企业的危险性H、路径扩散性P和敏感受体的脆弱性V，第i个像元的风险指数Ri＝Hi×Pi×Vi；S7、影响危险性H的第二层指标包括核密度K和经营时长DT，核密度代表了涉污企业的聚集程度，第n个涉污企业的危险性Hn＝DTn+Kn；影响路径扩散性P的第二层指标包括土壤质地ST、多年平均降雨量p、多年平均风速WS和地势起伏度RA，第i个像元的路径扩散性Pi＝1/4(STi+pi+WSi+RAi)；影响敏感受体的脆弱性V的第二层指标包括耕地产量；S8、根据估算出的风险指数R确定高风险区域，并将识别出的高风险区域与国家划定的重金属重点防控区进行对比，验证估算精度。本发明基于风场的影响范围空间模拟模型考虑了各个涉污企业所处环境的风场特殊性，模拟了具有场地针对性的影响范围，较以往单一给定影响范围的操作更具有科学性。基于敏感受体的栅格尺度风险评估方法，实现了高风险区空间位置精准识别，可以反映风险在行政单元内部的空间差异。
表1大气沉降影响范围：
涉污企业的行业类型包括金属矿采选型和金属冶炼型，所述金属矿采选型企业的基本影响距离为1km，所述金属冶炼型企业的基本影响距离为1.5km。金属矿采选型企业基本影响距离受影响因素为：经营时长<15年调整距离0km，经营时长≥15调整距离+0.5km，多年平均降水量<400mm调整距离+1km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km；金属冶炼型基本影响距离受影响因素为：经营时长<5年调整距离0km，经营时长5-15年调整距离+1km，经营时长>15年调整距离+2km，多年平均降水量<400mm调整距离+0.5km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km。每个风向椭圆以涉污企业的空间位置为顶点，月平均风向的下风方向为长轴方向。风向椭圆短轴和长轴的长度依据当月平均风速，在基本影响距离的基础上进行调整，长轴的长度则根据行业类型的不同比短轴长0.5km或1.0km；最后，将120个风向椭圆与初始的圆形缓冲区进行合并，并取包络面作为最后的影响范围。以某一金属矿采选业涉污企业为例，假设其经营时长为18年，所处地理位置的多年降水量为200mm，则其基本影响距离为1.0km，经营时长影响因素的调整距离为+0.5km，多年降水量影响因素的调整距离为+1.0km，圆形缓冲区的半径为1.0+0.5+1.0＝2.5km；若该涉污企业某一月的平均风速为6m/s，则当月的风向椭圆的短轴距离在基本影响距离1.0km的基础上，增加1.0km，最终为1.0+1.0＝2.0km，长轴距离比短轴长0.5km，为2.5km。
以往开展的大多数区域尺度风险评估，多是以行政单元为评估尺度。这样的评估由于没有考虑污染源和敏感受体之间的空间关系而容易引发误差。如当涉污企业空间分布远离人口或农田时，涉污企业对人口或耕地可能造成的风险应是较小的，但基于行政尺度的评估以整个区域的统计数据为评判标准，忽略了污染源和敏感受体之间的空间距离，从而导致得到的评估结果可能偏大于实际值。此外，基于行政单元的评估结果无法反映行政单元内部的风险差异，更无法确切反映高风险区的精准位置。
本发明提出的基于敏感受体的栅格尺度风险评估方法在1km栅格尺度上开展，且认为涉污企业仅对影响范围内的敏感受体有影响，因此评估只在模拟的影响范围内进行。评估主要分为以下两个步骤：
第一步构建评估指标体系：根据风险三要素理论，风险主要由污染源危害性、扩散路径扩散性和敏感受体脆弱性决定。涉污企业作为污染物的释放主体，是主要的污染源，其经营时长、空间聚集程度影响了危害性的高低。区域的水、土、地环境介质影响着污染物扩散的难易程度，如降水量，周围土壤质地，地表硬化覆盖程度等，都影响着污染物迁移扩散的难易程度。敏感受体是指受有毒有害物质影响的受害主体，敏感受体的脆弱性是指受体面对污染的易损性情况。主要的敏感受体有耕地、土壤、水源、人口等，此处主要以耕地敏感受体为例。当某片耕地受污染时，耕地产量越高，造成的经济损失越大，可能受害的人群越多，因此从该角度，认为耕地的生产力越高，其面对污染的脆弱性越高。
表2风险评估指标体系
表3指标分级与赋值
第二步构建风险评估模型，根据风险三要素理论，风险的形成必须具备污染源、扩散路径和敏感受体三个要素，因此风险R的估算遵循以下公式：
Ri＝Hi×Pi×Vi (1)
Hn＝DTn+Kn (2)
Pi＝1/4(STi+pi+WSi+RAi) (3)
式(1)表示，第i个像元的风险R取决于该像元上涉污企业的危险性H、由区域环境决定的污染物扩散性P和敏感受体的脆弱性V。单个涉污企业的危害性由公式(2)表示，n代表第n个涉污企业，核密度K代表了涉污企业聚集的程度。为将点位尺度的危害性指标转换为栅格尺度，以影响范围为插值空间，将危害性H插值为1km的空间数据。根据公式(3)，第i个像元的路径扩散性P由土壤质地ST、多年平均降雨量p、多年平均风速WS，以及地势起伏度RA决定。最终，具有相对高风险值的区域为高风险区。
本发明提出的金属行业影响范围模拟与风险评估方法，不需要人工实地去采集数据，因此完成速度快，能有效克服传统人工采样中消耗巨大人力、财力与时间等缺点。提出的影响范围模拟模型能够解决以往区域风险评估中影响范围设定单一、无法表征真实污染范围的问题，此外，相较于行政尺度的风险评估，本发明聚焦的栅格尺度评估不仅考虑了涉污企业和敏感受体的空间相互作用，还能获得精细的高风险区空间分布信息。模拟影响范围时，历史风场数据可以使用近十年，或者近三十年等不同历史长度的风场数据；风场的数据源可以是不同类型的栅格数据，如ERA再分析数值模拟数据，或风云系列卫星风场遥感反演产品等。同理，多年平均降水数据可以是近十年，或者近三十年等不同时间长度，数据源可以是由气象站点数据插值的数据，或数值模拟数据，或卫星遥感反演数据。开展风险评估时，可根据关注问题的不同，将耕地更替为水体、老年人口、儿童等其他敏感受体；涉污企业危害性指标可根据数据的可获得性，更替为涉污企业废气废水排放量等其他指标。
实施例1：
依据图2的模拟示意图，根据模拟方法及调整距离确定影响范围，并得到具体的影响范围的模拟结果。
计算危害性H：计算涉污企业核密度，归一化后，获取各企业的核密度归一化值Kn。并将各企业的经营时长归一化，得到DTn。根据Hn＝DTn+Kn计算各涉污企业的危害性Hn。为将点位尺度的危害性指标转换为栅格尺度，以影响范围为插值空间，将危害性H插值为1km的空间数据。
计算扩散性P：归一化土壤质地ST、多年平均降雨量p、多年平均风速WS、地势起伏度RA，根据Pi＝1/4(STi+pi+WSi+RAi)计算第i个像元的路径扩散性，且只保留影响范围内的像元。
计算脆弱性V：将影响范围内的耕地，根据产量高、中、低三个等级进行归一化。
计算风险指数R：计算1km栅格尺度的风险指数时，影响范围内的第i个像元的风险指数Ri＝Hi×Pi×Vi。计算县级尺度的风险指数时，基于1km栅格尺度的风险指数，统计各县内风险指数的总和，作为县的风险水平值。
风险评估结果：包括1km栅格尺度和县级尺度两种结果。1km栅格尺度评估结果显示，云南省东部地区比西部地区面临着更大的重金属污染威胁，其中世界级“锡都”——个旧市面临着最大的威胁。将县的风险等级由大到小排序后，根据自然分段法，划分出高风险区。本研究识别出的高风险区共12个，包括个旧区、腾冲县、安宁区、易门县、东川区、会泽县、马关县、文山市、建水县、蒙自市、砚山县及宣威市。
精度验证：县级评估结果与国家划定的重金属重点防控区对比。
本方法识别出的高风险区共12个，包括个旧区、腾冲县、安宁区、易门县、东川区、会泽县、马关县、文山市、建水县、蒙自市、砚山县及宣威市。根据国家生态环境部(参考资料[1]中华人民共和国生态环境部.云南部署重金属污染防治盯紧重点区域和重点行业。http://www.mee.gov.cn/ywdt/hjnews/201105/t20110520_210913.shtml)、云南省生态环境厅(参考资料[2]云南省生态环境厅.红河州公开省委、省政府环境保护督察反馈意见问题整改情况。http://sthjt.yn.gov.cn/zwxx/xxyw/xxywrdjj/201907/t20190730_191420.html.)、曲靖市人民政府(参考资料[3]曲靖市人民政府《陆良县西桥工业片区重金属污染防治方案》通过2015年重金属污染防治重点区域竞争性评审。http://www.qj.gov.cn/html/2015/hbj_0625/27153.html)及最新的研究结果(参考资料[4]Tan，Z.，Lu，S.，Zhao，H.et al，2018.Magnetic，geochemical characterization and healthrisk assessment of road dust in Xuanwei and Fuyuan，China.Environ.Geochem.Health 40(4)，1541–1555。参考资料[5]王红英.陆良县重金属污染防治主要措施及成效[J].环境科学导刊，2017，36(S2):139-142.)明确的高风险区共计14个，本方法识别精度为85.7％。
国家于2011年，在云南省划定了11个重金属重点防控区，包括：昆明市东川区、安宁区、红河州个旧区、金平区、曲靖市会泽县、陆良县、文山州马关县、文山市、保山市腾冲县、玉溪市易门县、怒江州兰坪县。
对比结果显示：成功识别个旧、腾冲、安宁、易门、东川、会泽、马关、文山市8个地区，补充识别建水、蒙自、砚山、宣威4个地区，未能识别兰坪、金平、陆良3个地区。补充识别的4个地区中，建水和蒙自是红河州重金属污染重点调查区；砚山于2016年被列为重金属污染防治试点；宣威受重金属污染影响，是云南肺癌高发区。未能识别的3个地区中，兰坪、金平危害性较高，但脆弱性低；陆良县于2015年因治理得当而取消了“重金属污染防控区”的帽子。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。
Claims (7)
1.一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于，该方法包括以下步骤：
S1、根据涉污企业的行业类型确定所述涉污企业的基本影响距离；
S2、根据涉污企业的经营时长和所处地理位置的多年平均降水量对涉污企业的基本影响距离进行调整，设置涉污企业的污染半径，并以涉污企业为中心绘制圆形缓冲区；
S3、根据涉污企业近十年每月的平均风场环境绘制风向椭圆，共计绘制120个风向椭圆；
S4、将120个风向椭圆与以涉污企业为中心绘制的圆形缓冲区进行合并，合并后再取包络面为涉污企业的影响范围；
S5、采用基于敏感受体的栅格尺度风险评估方法在1km栅格尺度上开展，且认为涉污企业仅对影响范围内的敏感受体有影响，评估只在模拟的影响范围内进行；
S6、建立评估指标体系，所述风险指数为R，影响所述风险指数的第一层指标包括涉污企业的危险性H、路径扩散性P和敏感受体的脆弱性V，第i个像元的风险指数Ri＝Hi×Pi×Vi；
S7、影响危险性H的第二层指标包括核密度K和经营时长DT，核密度代表了涉污企业的聚集程度，第n个涉污企业的危险性Hn＝DTn+Kn；影响路径扩散性P的第二层指标包括土壤质地ST、多年平均降雨量p、多年平均风速WS和地势起伏度RA，第i个像元的路径扩散性Pi＝1/4(STi+pi+WSi+RAi)；影响敏感受体的脆弱性V的第二层指标包括耕地产量；
S8、根据估算出的风险指数R确定高风险区域，并将识别出的高风险区域与国家划定的重金属重点防控区进行对比，验证估算精度。
2.根据权利要求1所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：在步骤S1中，所述涉污企业的行业类型包括金属矿采选型和金属冶炼型，所述金属矿采选型企业的基本影响距离为1km，所述金属冶炼型企业的基本影响距离为1.5km。
3.根据权利要求2所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：在步骤S2中，金属矿采选型企业基本影响距离受影响因素为：经营时长<15年调整距离0km，经营时长≥15调整距离+0.5km，多年平均降水量<400mm调整距离+1km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km；
金属冶炼型基本影响距离受影响因素为：经营时长<5年调整距离0km，经营时长5-15年调整距离+1km，经营时长>15年调整距离+2km，多年平均降水量<400mm调整距离+0.5km，多年平均降水量400–800mm调整距离0km，多年平均降水量>800mm调整距离-0.5km。
4.根据权利要求2所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：在步骤S3中，每个风向椭圆以涉污企业的空间位置为顶点，月平均风向的下风方向为长轴方向；金属矿采选型企业的风向椭圆长短轴在基本影响距离的基础上调整为：月平均风速<3m/s，短轴长度调整距离0km，月平均风速3-5m/s调整距离+0.5km，月平均风速5-7m/s调整距离+1km，月平均风速>7m/s调整距离+2km，长轴长度在短轴长度基础上+0.5km；金属冶炼型企业的风向椭圆长短轴在基本影响距离的基础上调整为：月平均风速<2m/s调整距离+0.5km，月平均风速2-4m/s调整距离0km，月平均风速>4m/s调整距离-0.5km，长轴长度在短轴长度基础上+1km。
5.根据权利要求1所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：敏感受体是指受有毒有害物质影响的受害主体，敏感受体的脆弱性是指受体面对污染的易损性情况，主要的敏感受体包括耕地、土壤、水源或人口。
6.根据权利要求1所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：为将点位尺度的危害性指标转换为栅格尺度，以影响范围为插值空间，将危害性H插值为1km的空间数据。
7.根据权利要求1所述的一种金属行业涉污企业影响范围模拟与风险评估方法，其特征在于：在步骤S8中，识别出的高风险区域与国家划定的重金属重点防控区吻合度越高，则识别与评估的估算精度越高。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110233167.2A CN112785206B (zh) | 2021-03-03 | 2021-03-03 | 一种金属行业涉污企业影响范围模拟与风险评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110233167.2A CN112785206B (zh) | 2021-03-03 | 2021-03-03 | 一种金属行业涉污企业影响范围模拟与风险评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112785206A true CN112785206A (zh) | 2021-05-11 |
| CN112785206B CN112785206B (zh) | 2023-06-30 |
Family
ID=75762243
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110233167.2A Active CN112785206B (zh) | 2021-03-03 | 2021-03-03 | 一种金属行业涉污企业影响范围模拟与风险评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112785206B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976061A (zh) * | 2010-08-06 | 2011-02-16 | 中国环境科学研究院 | 重大环境风险源监控系统构建方法 |
| CN102062769A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院地理科学与资源研究所 | 一种陆地石油开采生态风险识别及风险值的计算方法 |
| CN105321122A (zh) * | 2015-09-23 | 2016-02-10 | 国家电网公司 | 一种污区分布图的现场验证方法 |
| CN106269837A (zh) * | 2016-08-26 | 2017-01-04 | 北京航天环境工程有限公司 | 一种挥发性有机污染土壤的修复系统及方法 |
| CN107469563A (zh) * | 2017-09-28 | 2017-12-15 | 广州京海科技有限公司 | 一种城市道路大气污染净化装置 |
| CN109154017A (zh) * | 2016-03-17 | 2019-01-04 | 不列颠哥伦比亚大学 | 细胞分泌分析的装置和方法 |
| CN110009226A (zh) * | 2019-04-03 | 2019-07-12 | 中国科学院地理科学与资源研究所 | 一种场地污染风险等级评测及敏感受体影响识别系统 |
| CN110414153A (zh) * | 2019-07-31 | 2019-11-05 | 辽宁工程技术大学 | 一种确定露天煤矿粉尘扩散范围的方法 |
| CN110849837A (zh) * | 2019-11-26 | 2020-02-28 | 南昌工程学院 | 一种大气污染叠加风险区识别及评价方法 |
| CN111507610A (zh) * | 2020-04-14 | 2020-08-07 | 江苏省环科院环境科技有限责任公司 | 一种土壤重金属污染的累积性环境风险预警方法 |
-
2021
- 2021-03-03 CN CN202110233167.2A patent/CN112785206B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976061A (zh) * | 2010-08-06 | 2011-02-16 | 中国环境科学研究院 | 重大环境风险源监控系统构建方法 |
| CN102062769A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院地理科学与资源研究所 | 一种陆地石油开采生态风险识别及风险值的计算方法 |
| CN105321122A (zh) * | 2015-09-23 | 2016-02-10 | 国家电网公司 | 一种污区分布图的现场验证方法 |
| CN109154017A (zh) * | 2016-03-17 | 2019-01-04 | 不列颠哥伦比亚大学 | 细胞分泌分析的装置和方法 |
| CN106269837A (zh) * | 2016-08-26 | 2017-01-04 | 北京航天环境工程有限公司 | 一种挥发性有机污染土壤的修复系统及方法 |
| CN107469563A (zh) * | 2017-09-28 | 2017-12-15 | 广州京海科技有限公司 | 一种城市道路大气污染净化装置 |
| CN110009226A (zh) * | 2019-04-03 | 2019-07-12 | 中国科学院地理科学与资源研究所 | 一种场地污染风险等级评测及敏感受体影响识别系统 |
| CN110414153A (zh) * | 2019-07-31 | 2019-11-05 | 辽宁工程技术大学 | 一种确定露天煤矿粉尘扩散范围的方法 |
| CN110849837A (zh) * | 2019-11-26 | 2020-02-28 | 南昌工程学院 | 一种大气污染叠加风险区识别及评价方法 |
| CN111507610A (zh) * | 2020-04-14 | 2020-08-07 | 江苏省环科院环境科技有限责任公司 | 一种土壤重金属污染的累积性环境风险预警方法 |
Non-Patent Citations (1)
| Title |
|---|
| 兰冬东，等: "区域环境污染事件风险分区技术及其应用", 《应用基础与工程科学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112785206B (zh) | 2023-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Improved enrichment factor model for correcting and predicting the evaluation of heavy metals in sediments | |
| CN110009226B (zh) | 一种场地污染风险等级评测及敏感受体影响识别系统 | |
| Ban et al. | Effect of urbanization on aerosol optical depth over Beijing: Land use and surface temperature analysis | |
| Khuzestani et al. | Quantification of the sources of long-range transport of PM2. 5 pollution in the Ordos region, Inner Mongolia, China | |
| Sun et al. | Development of emissions inventory and identification of sources for priority control in the middle reaches of Yangtze River Urban Agglomerations | |
| CN105631537A (zh) | 基于气象服务平台的空气质量预报业务系统 | |
| CN102628852A (zh) | 基于污染物来源识别技术的大气污染源分级方法 | |
| Yang et al. | Human activities and the natural environment have induced changes in the PM2. 5 concentrations in Yunnan Province, China, over the past 19 years | |
| CN110658307A (zh) | 评估污染源对环境空气质量影响的方法 | |
| CN111882243A (zh) | 一种喀斯特地区国土空间开发适宜性评价方法及系统 | |
| Lin et al. | Environmental risk mapping of potential abandoned uranium mine contamination on the Navajo Nation, USA, using a GIS-based multi-criteria decision analysis approach | |
| Cao et al. | Spatially resolved risk assessment of environmental incidents in China | |
| CN115902114B (zh) | 基于半定量法的小尺度大气污染溯源方法 | |
| Qiu et al. | Effect of terrain gradient on cadmium accumulation in soils | |
| Wu et al. | Application of an integrated Weather Research and Forecasting (WRF)/CALPUFF modeling tool for source apportionment of atmospheric pollutants for air quality management: A case study in the urban area of Benxi, China | |
| Li et al. | Application of MM5‐CAMx‐PSAT Modeling Approach for Investigating Emission Source Contribution to Atmospheric SO2 Pollution in Tangshan, Northern China | |
| Cheng et al. | Application of trajectory clustering and source apportionment methods for investigating trans-boundary atmospheric PM10 pollution | |
| Huang et al. | Ecological risk assessment from the perspective of soil heavy metal accumulations in Xiamen city, China | |
| Kang et al. | Identifying the potential soil pollution areas derived from the metal mining industry in China using MaxEnt with mine reserve scales (MaxEnt_MRS) | |
| Shan et al. | The pattern and mechanism of air pollution in developed coastal areas of China: From the perspective of urban agglomeration | |
| CN114219345B (zh) | 一种基于数据挖掘的二次空气质量预报优化方法 | |
| Wang et al. | Risk assessment of typhoon disaster chains in the Guangdong–Hong Kong–Macau greater bay area, China | |
| Zhao et al. | Potential source regions and transportation pathways of reactive gases at a regional background site in northwestern China | |
| Garg et al. | Determining the contribution of long-range transport, regional and local source areas, to PM10 mass loading in Hessen, Germany using a novel multi-receptor based statistical approach | |
| Li et al. | Multipath diffusion process and spatial accumulation simulation of Cd in lead-zinc mining areas |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |