CN117556984A - 一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 - Google Patents
一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 Download PDFInfo
- Publication number
- CN117556984A CN117556984A CN202311322998.2A CN202311322998A CN117556984A CN 117556984 A CN117556984 A CN 117556984A CN 202311322998 A CN202311322998 A CN 202311322998A CN 117556984 A CN117556984 A CN 117556984A
- Authority
- CN
- China
- Prior art keywords
- risk
- classification
- site
- index
- pollution risk
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000003673 groundwater Substances 0.000 title claims abstract description 33
- 238000003900 soil pollution Methods 0.000 title claims abstract description 26
- 238000003745 diagnosis Methods 0.000 title claims abstract description 19
- 238000003895 groundwater pollution Methods 0.000 title claims abstract description 19
- 238000010801 machine learning Methods 0.000 title claims abstract description 19
- 238000013528 artificial neural network Methods 0.000 claims abstract description 40
- 238000013145 classification model Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- 238000007726 management method Methods 0.000 claims description 44
- 239000003344 environmental pollutant Substances 0.000 claims description 29
- 231100000719 pollutant Toxicity 0.000 claims description 29
- 239000002689 soil Substances 0.000 claims description 25
- 238000004519 manufacturing process Methods 0.000 claims description 22
- 238000012544 monitoring process Methods 0.000 claims description 21
- 230000007774 longterm Effects 0.000 claims description 18
- 238000012549 training Methods 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 12
- 229910052751 metal Inorganic materials 0.000 claims description 12
- 239000002184 metal Substances 0.000 claims description 12
- 230000008439 repair process Effects 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 9
- 238000003723 Smelting Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 7
- 238000013508 migration Methods 0.000 claims description 6
- 230000005012 migration Effects 0.000 claims description 6
- 210000002569 neuron Anatomy 0.000 claims description 6
- 230000006872 improvement Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 claims description 4
- 239000003208 petroleum Substances 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 239000013064 chemical raw material Substances 0.000 claims description 3
- 239000003245 coal Substances 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 239000003802 soil pollutant Substances 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims description 2
- 239000003209 petroleum derivative Substances 0.000 claims description 2
- 238000003490 calendering Methods 0.000 claims 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims 1
- 239000003651 drinking water Substances 0.000 claims 1
- 235000020188 drinking water Nutrition 0.000 claims 1
- 239000003345 natural gas Substances 0.000 claims 1
- 238000012216 screening Methods 0.000 description 13
- 238000012706 support-vector machine Methods 0.000 description 12
- 238000010276 construction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 6
- 238000011835 investigation Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000033228 biological regulation Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 5
- 239000002994 raw material Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000011478 gradient descent method Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000035699 permeability Effects 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000009189 diving Effects 0.000 description 3
- 229910001385 heavy metal Inorganic materials 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000008020 evaporation Effects 0.000 description 2
- 238000001704 evaporation Methods 0.000 description 2
- 239000002957 persistent organic pollutant Substances 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 239000002910 solid waste Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000010754 BS 2869 Class F Substances 0.000 description 1
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009713 electroplating Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 229910052755 nonmetal Inorganic materials 0.000 description 1
- 230000033116 oxidation-reduction process Effects 0.000 description 1
- 235000013550 pizza Nutrition 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002407 reforming Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000002076 thermal analysis method Methods 0.000 description 1
- 230000005619 thermoelectricity Effects 0.000 description 1
- 238000007514 turning Methods 0.000 description 1
- 239000012855 volatile organic compound Substances 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Molecular Biology (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Evolutionary Biology (AREA)
- Game Theory and Decision Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Primary Health Care (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明公开了一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法,包括:构建场地污染风险分类分级案例库,建立场地污染风险分类和分级指标体系,确定场地污染风险分类和分级指标规则;对场地污染风险分类和分级的各指标进行赋分,并运用层次分析法中1~9标度法对各指标的分值进行一致性检验;通过调整BP神经网络和SVM的各项参数,优化最佳的分类、分级模型,基于BP神经网络进行风险分类,基于SVM进行风险分级,并在MATLAB中完成目标案例的风险分类分级预测,根据风险分类分级的预测结果进行风险管控。本发明通过BP神经网络和SVM算法,构建分类分级模型,实现快速、智能、准确地场地污染风险分类分级,为风险管控提供可靠参考。
Description
技术领域
本发明涉及土壤重金属污染分析技术领域,特别是涉及一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法。
背景技术
随着城市布局调整和产业结构升级,大量化工、石油、塑料、炼焦、热电、钢铁、电镀、有色治炼等重污染企业搬迁或关停,遗留大量的污染场地。污染场地中的重金属或有机污染物对土壤、地下水形成长期污染,极大威胁人类健康。污染场地修复成为环保行业的热点问题。但是,由于场地污染深度大、地理环境特殊等原因,部分污染场地无法在短期内实施达标修复,需采用风险管控方法阻断场地内的污染物迁移或挥发途径,以避免污染物扩散。现有的污染场地风险管控技术无法对中长期场地污染风险进行准确的分类分级,不利于风险管控。
发明内容
本发明的目的是提供一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法,通过建立场地污染风险分类、分级指标体系以及指标规则,明确分类分级标准,并通过BP神经网络和SVM算法,构建分类分级模型,实现快速、准确地场地污染风险分类分级,为风险管控提供可靠参考。
为实现上述目的,本发明提供了如下方案:
一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法,该方法包括以下步骤:
构建场地污染风险分类分级案例库,包括有风险案例和无风险案例;
根据场地污染风险判定需求及场地污染风险分类分级案例,分别建立场地污染风险分类和分级指标体系;
分别确定场地污染风险分类和分级指标规则;
根据场地污染风险分类和分级指标规则,对场地污染风险分类和分级的各指标进行赋分,并运用层次分析法中1~9标度法对各指标的分值进行一致性检验;
通过调整BP神经网络和SVM的各项参数,优化最佳的分类、分级模型,基于BP神经网络进行风险分类,基于SVM进行风险分级,并在MATLAB中完成目标案例的风险分类分级预测,根据风险分类分级的预测结果进行风险管控。
进一步地,所述根据风险分类分级的预测结果进行风险管控,具体包括:
针对在产+高风险场地,从源头污染物减排、生产技术升级改造、污染物治理能力提升三方面开展风险管控;针对在产+低风险场地,进行长期监测,在厂区不同功能区布设长期监测点位并进行样品采集及分析测试;针对具备修复条件的关闭搬迁+高风险场地,从治理修复和被动修复-衰减两方面开展风险管控;针对不具备修复条件的关闭搬迁+高风险场地,采取工程控制、长期监测和制度控制相结合的方法;针对关闭搬迁+低风险场地,采取长期监测和准入管理相结合的方法。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,针对污染场地风险管控的技术要求,结合有风险案例和无风险案例,构建场地污染风险分类、分级指标体系以及指标规则,明确分类分级标准,并综合BP神经网络和SVM算法,构建分类、分级模型,大幅减少风险评价指标数量,缩短传统风险评价工作流程(如省去初步调查、详细调查环节),缩减传统风险评价工作量,节省工作成本投入,缩短决策时间,实现快速、智能、准确地场地污染风险分类分级,为风险管控提供可靠参考,本发明操作简单、简单实用、技术可靠、易于工程化实施,可用于有机污染场地、重金属污染场地和复合污染场地的中长期风险管控,具有重要的环保和社会意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于机器学习的场地地下水和土壤污染风险诊断与管控方法的流程图;
图2为训练函数对BP神经网络预测准确性的影响示意图;
图3为学习率对BP神经网络预测准确性的影响示意图;
图4为神经元个数对BP神经网络预测准确性的影响示意图;
图5为本发明实施例BP神经网络支持的风险分类结果示意图;
图6为本发明实施例支持向量机支持的风险分类结果示意图;
图7为本发明实施例BP神经网络支持的风险分级结果示意图;
图8为本发明实施例支持向量机支持的风险分级结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法,通过建立场地污染风险分类、分级指标体系以及指标规则,明确分类分级标准,并通过BP神经网络和SVM算法,构建分类分级模型,实现快速、准确地场地污染风险分类分级,为风险管控提供可靠参考。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,包括以下步骤:
收集和整理场地污染风险分类分级案例,构建场地污染风险分类分级案例库,包括有风险案例和无风险案例;
根据场地污染风险判定需求及场地污染风险分类分级案例,分别建立场地污染风险分类和分级指标体系;
分别确定场地污染风险分类和分级指标规则;
根据场地污染风险分类和分级指标规则,对场地污染风险分类和分级的各指标进行赋分,并运用层次分析法中1~9标度法对各指标的分值进行一致性检验;
通过调整BP神经网络和SVM的各项参数,优化最佳的分类分级模型,基于BP神经网络进行风险分类,基于SVM进行风险分级,并在MATLAB中完成目标案例的风险分类分级预测,根据风险分类分级的预测结果进行风险管控。
示例地,所述收集和整理场地污染风险分类分级案例,具体包括:
借助全国土壤环境信息平台、200个试点示范项目、文献资料等,收集和整理场地土壤污染风险分类分级案例(包括初步调查报告、详细调查报告、风险评估报告),获得274个有风险案例和55个无风险案例,涉及化学原料和化学制品制造业,有色金属冶炼和压延加工业,黑色金属冶炼和压延加工业,金属制品业,医药制造业,石油、煤炭及其他燃料加工业等。
参考《重点行业企业用地调查信息采集技术规定(试行)》附录3中关闭搬迁企业地块信息调查表,确定场地污染风险分类分级案例框架;将收集的案例按照该案例框架进行整理;最后,将所有案例进行汇总形成案例库。
本发明所应用的大数据算法的基本原理如下:
1)BP神经网络
BP是误差反向后传(BackPropagation)的英文缩写,原理是利用输出后的误差估计输出层的直接前导层,再利用这个误差估计更前一层的误差,如此一层层地反传下去,就可以获得所有其他各层的误差估计。
2)支持向量机
SVM是一类按监督学***面。SVM采用不同的内积函数(核函数)将导致不同的算法。在使用SVM对数据进行训练时,在准确性符合一定要求条件下模型的泛化能力应尽可能更好,以防止在训练过程中发生过拟合现象。因此,在训练时,需要不断的调整gamma和c值,并对数据不断地进行交叉验证,以找到合适的gamma和c值。gamma值越大,映射的维度越高,训练的结果越好,但是越容易引起过拟合,即泛化能力低;gamma值越小,训练的结果准确率低,拟合程度差。另一方面,c值越大,分类越严格,正确率越高,但相应的泛化能力降低;c值越小,意味着有更大的错误容忍度,但正确率越低。
基于案例库,所述基于机器学习的场地地下水和土壤污染风险诊断与管控方法,具体包括:
1)确定场地污染风险管控源案例
从案例库中选取65个有风险案例和55个无风险案例,作为场地污染风险分类分级源案例,涉及化学原料和化学制品制造业53个,非金属矿物制品业4个,有色金属冶炼和压延加工业8个,黑色金属冶炼和压延加工业12个,金属制品业20个,石油、煤炭及其他燃料加工业3个和其他行业20个。
2)构建场地污染风险分类、分级指标体系
参考《重点行业企业用地调查信息采集技术规定(试行)》(环办土壤〔2017〕67号)中附录3的关闭搬迁企业地块信息调查表,分别建立场地污染风险分类和分级指标体系。
3)确定场地污染风险分类和分级指标规则
基于《土壤环境质量建设用地土壤污染风险管控标准(试行)》(GB 36600-2018),参照《土地利用现状分类》(GB/T 21010-2017)和《地下水水质标准》(DZ/T 0290-2015)中Ⅲ类水标准,分别确定场地污染风险分类和分级指标规则。
4)指标赋分一致性检验
运用层次分析法中1~9标度法(如表1所示)对各指标分值进行一致性检验,通过检验则认为分值可接受。
表1因子相对重要性标度及其对应含义
在MATLAB软件中,根据各特征指标信息的特征值中最大值和特征指标信息的数目计算一致性指标,如公式(1)所示。
式中,IC为一致性指标,λmax为判断矩阵的最大特征值,n为特征指标信息的个数。
根据一致性指标和预设随机一致性指标,计算一致性比率,如公式(2)所示:
式中,RC为一致性比率,IR为预设随机一致性指标(如表2所示)。
表2随机一致性指标分值
n | 1 | 2 | 3 | 4 | 5 | 6 | ... | 22 | 23 | 24 | 25 |
IR | 0 | 0 | 0.5303 | 0.8845 | 1.1143 | 1.2525 | ... | 1.6401 | 1.6482 | 1.6511 | 1.6550 |
5)大数据算法优化
通过调整BP神经网络和SVM的各项参数,优化最佳的分类和分级模型,并在MATLAB中完成目标案例的分类分级预测。选用相关系数(R2)、平均相对误差(MRE)及均方根误差(RMSE)作为评价标准,评价分类和分级模型性能,如公式(3)和公式(4)所示。
式中:n为测试样本数,A(i)为样本的实测值,P(i)为样本的预测值。
其中,BP神经网络和SVM参数优化的具体过程如下:
(1)BP神经网络优化
选取合适的BP神经网络参数对风险判定结果准确性有重要影响,影响BP神经网络参数有训练函数、学习率、隐含层神经元个数。目前,常用的训练函数有梯度下降法(traingd)、有动量的梯度下降法(traingdm)、自适应学习率梯度下降法(traingda)、自适应学习率动量梯度下降法(traingdx)、共轭梯度法(traincg)、拟牛顿算法(trainbfg)及Levenberg-Marquardt算法(trainlm)等7种。本发明中,训练函数对BP神经网络预测准确性的影响见图2。trainlm具有最新小的MSE(0.00043),且不存在最大错误次数的影响(如图2)。因此,选取trainlm为BP神经网络的训练函数。
学习率对BP神经网络预测准确性的影响见图3。在训练次数接近一致情况下,学习率为0.1时,模型的MSE为0.000289(图3)。因此,选取模型的学习率为0.1。
神经元个数对BP神经网络预测准确性的影响见图4。当隐含层神经元个数为6时,模型均方误差最小(0.000314)(图4)。
综上所述,通过调整模型超参数,在trainlm为训练函数、隐含层神经元取6、学习率取0.1时,模型的性能最好,拟合效果好,收敛速度快。
(2)SVM参数优化
如表3所示,模型2与其他两个模型相比,正确率最高,达73.3%。可见,核函数为linear、cost为1、gamma为0.08时模型效果最好。
表3支持向量机不同参数对比
模型 | 核函数 | 最优参数 | 正确率(%) |
模型1 | radial | cost为5gamma为0.1 | 66.7 |
模型2 | linear | cost为1gamma为0.08 | 73.3 |
模型3 | polynomial | cost为0gamma为1 | 66.7 |
本发明涉及场地污染风险分类以及场地污染风险分级,采用的指标体系以及模型不同,具体包括:
一、场地污染风险分类
(1)指标体系构建
根据场地污染风险判定需求及收集的有风险案例和无风险案例信息,构建场地污染风险分类指标体系,对场地污染是否存在风险进行判定(表4)。表格中,Es表示土壤污染超标总倍数。
表4风险分类指标体系
(2)计算规则建立
1)干湿指数规则
干湿指数规则见公式(5):
式中,DWI表示干湿指数(无量纲);P表示多年平均降水量(mm);E0表示多年平均蒸发量(mm),计算降水量和蒸发量时所取的年份数量可根据实际需求进行调整。
2)特征污染物规则
土壤特征污染物和地下水特征污染物中SVOCs、VOCs、重金属和无机物、石油烃主要关注《土壤环境质量建设用地土壤污染风险管控标准(试行)》(GB 36600-2018)中85项污染物。土壤污染物倍数计算见公式(6)。
式中,n为土壤中浓度超过筛选值的污染物种类数量;Csi为浓度超过筛选值的第i种污染物浓度,取场地所有土壤样品测试数据中最高值(mg/kg);RSsi为第i种污染物的筛选值(mg/kg),参照建设用地土壤污染风险筛选指导值。
地下水污染物倍数计算见公式(7):
式中,n为地下水中浓度超过筛选值的污染物种类数量;Cgwi为浓度超过筛选值的第i种污染物浓度,取场地所有地下水样品测试数据中最高值(mg/L);RSgwi为第i种污染物的限值(mg/L),参照《地下水水质标准》(DZ/T 0290-2015)中Ⅲ类水质标准。
3)土地利用规划
基于《土地利用现状分类》(GB/T 21010-2017)和《土壤环境质量建设用地土壤污染风险管控标准(试行)》(GB 36600-2018),根据场地未来规划确定的场地土地利用规划等级建立土地利用规划规则(表5)。
表5土地利用规划规则
4)城市等级
根据《关于调整城市规模划分标准的通知》(国发〔2014〕51号)确定城市等级,其中A类等级城市包括北京、上海、广州和深圳4个;B类等级城市包括重庆、天津、杭州和南京等21个城市;C类等级城市包括太原、呼和浩特、贵阳和兰州等41个城市;D类等级城市包括银川、西宁、拉萨和三亚等62个城市;E类等级城市包括肇庆、清远、汕尾和梅州等190个城市;F类为其他城市(表6)。
表6城市等级赋分规则
城市等级 | 赋分 |
A | 6 |
B、C | 4 |
D、E、F | 1 |
5)包气带的影响和潜水含水层的影响
结合源案例和目标案例场地内包气带和潜水含水层实际,确定包气带渗透系数最大的岩性赋分规则和潜水含水层最主要岩性和赋分规则,包气带渗透系数最大的岩性及其赋值的对应关系见表7和表8。
表7包气带渗透系数赋分规则
表8含水层介质规则
岩性类型 | 渗透系数 | 赋分 |
粘土或不确定 | 1×10-11-4.7×10-9 | 2 |
粉土 | 6×10-7-6×10-6 | 2 |
粉砂 | 6×10-6-1.2×10-5 | 2 |
细砂 | 1.2×10-5-6×10-5 | 4 |
中砂 | 6×10-3-2.4×10-2 | 4 |
粗砂 | 2.4×10-2-6×10-2 | 4 |
砾石 | 6×10-2-1.8×10-1 | 6 |
(3)指标赋分
根据指标规则,给场地污染风险分类指标进行赋分(表9)。
表9场地污染风险分类指标分值
/>
/>
(4)分值一致性检验
根据风险分类指标赋分一致性检验矩阵,如公式(8),进行判断矩阵的一致性检验,确定一致性指标(IC)为11.4517、随机一致性指标(IR)为1.5079、一致性比率(RC)为0.0299<0.1,认为风险分类指标赋分可接受。
(5)分类模型筛选
BP神经网络和支持向量机均能较好地进行分类预测,表现为R2分别为0.9974和0.9103;另外,BP神经网络的MRE(0.12%)和RMSE(0.0376)均小于支持向量机的MRE(49.1%)和RMSE(13.7232),如表10。显然,BP神经网络的预测性能最好,所以在风险分类时选择BP神经网络进行预测(图5和图6)。
表10风险判定误差结果对比
评价指标 | BP神经网络 | 支持向量机 |
R2 | 0.9974 | 0.9103 |
MRE | 0.12% | 49.1% |
RMSE | 0.0376 | 13.7232 |
二、场地污染风险分级
(1)指标体系构建
参考《企业地块风险筛查与风险分级技术规定》,分别建立风险分级有关土壤指标体系和地下水指标体系(表11和表12)。
表11风险分级有关土壤指标体系
/>
/>
表12风险分级有关地下水指标体系
/>
(2)指标赋分
参考《企业地块风险筛查与风险分级技术规定》,进行场地污染风险分级有关土壤指标和地下水指标赋分(表13和表14)。
表13风险分级土壤指标赋分
/>
表14风险分级地下水指标赋分
/>
/>
(3)分值一致性检验
根据场地污染风险分级土壤指标赋分检验矩阵,如公式(9),在进行风险分级矩阵一致性检验时,IC为11.5279、IR为1.5079、RC为0.035<0.1,认为风险分级阶段指标分值可以接受。
(4)场地风险分值计算与分级
根据《企业地块风险筛查与风险分级技术规定》和收集到的污染地块基础信息,计算场地风险分值(公式10),并根据分值将场地分为低风险、中风险和高风险(表15)。
式中:S为地块风险筛查总分,Ss为地块土壤得分,Sgw为地块地下水得分。
表15风险等级划分规则
地块风险筛查总分 | 地块风险等级划分 |
S≥70分 | 高风险地块 |
50≤S<70分 | 中风险地块 |
S<50分 | 低风险地块 |
(5)分级模型筛选
如图7和图8所示,BP神经网络和支持向量机均能较好地进行分类预测,表现为R2分别为0.9700和0.9977;另外,虽然支持向量机的MRE(0.50%)和MSE(0.4341)大于BP神经网络的MRE(0.05%)和RMSE(0.0358),但是误差值亦属于允许范围内(表16)。显然,综合考虑R2、MRE和RMSE,支持向量机的预测性能较好,所以在风险分级时选择支持向量机进行预测。
表16风险分级误差结果对比
评价指标 | BP神经网络 | 支持向量机 |
R2 | 0.970 | 0.997 |
MRE | 0.05% | 0.50% |
RMSE | 0.0358 | 0.4341 |
进一步地,所述根据风险分类分级的预测结果进行风险管控,具体包括:
(1)场景1:在产+高风险场地。主要从源头污染物减排、生产技术升级改造、污染物治理能力提升等三方面开展风险管控。其中,源头污染物减排主要采取环保生产原料和辅料替换,从生产线源头减少污染物产生;生产技术改造主要包括生产技术升级、落后产能淘汰、生产设备改造等措施,提高生产效率、减少源头原料用量,从而减少后端污染物产生和排放;污染物治理能力提升包括企业治污设施升级改造、固废堆存场所整治等,降低污染物排放浓度、减少环境扰动。
(2)场景2:在产+低风险场地。主要采取长期监测措施,防止风险发生。在厂区不同功能区布设长期监测点位并进行样品采集及分析测试,实现控制潜在暴露风险的目的。长期监测点位布设时,需在生产车间、固废和原料堆存场所、排污管线布设区域进行加密布点,从而得以及时评估和发现环境风险。
(3)场景3:具备修复条件的关闭搬迁+高风险场地。主要从治理修复和被动修复-衰减两方面开展风险管控。其中,治理修复是利用一系列工程技术开展场地污染治理,具体工程技术包括客土、换土、翻土、热解析、固化稳定化、无害化填埋、土壤淋洗、氧化还原、植物修复等;被动修复-衰减通过利用工程技术及场地内自然条件,依靠自然发生的生物、化学或物理过程对污染物进行固定、削减或去除,具体技术包括渗透反应格栅、强化监测自然衰减等,主要适用于有机污染物迁移较慢、治理修复预算相对较低的场地。
(4)场景4:不具备修复条件的关闭搬迁+高风险场地。主要采取工程控制、长期监测和制度控制相结合的方法,防止风险扩散。其中,工程控制主要利用各种阻隔工程,限制污染物迁移,切断污染源与受体之间的暴露途径,达到降低污染暴露风险和保护受体安全的目的,具体工程控制技术包括泥浆墙、灌浆墙、板桩墙、土壤原位搅拌、土工膜、衬层等;长期监测主要环节与在产低+风险场地长期监测一致,主要区别为关闭搬迁+高风险场地监测点位设置时,需与场地调查结果相结合,在污染物浓度较高的地区加密布点;在采取长期监测措施情况下,同步开展的制度控制包括政府控制、所有权控制、强制与许可制度、信息公开(如设立标识牌)等,确保场地周边土壤环境和人居环境安全。
(5)场景5:关闭搬迁+低风险场地。主要采取长期监测和准入管理相结合的方法,防范环境风险。其中,长期监测点位布设时,需在曾经的生产区、原料堆存区及场地周边敏感目标区域设置监测点位;在用途变更时,需当地生态环境管理部门出具场地无风险证明后,方可进入用地审批程序。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,包括以下步骤:
构建场地污染风险分类分级案例库,包括有风险案例和无风险案例;
根据场地污染风险判定需求及场地污染风险分类分级案例,分别建立场地污染风险分类和分级指标体系;
分别确定场地污染风险分类和分级指标规则;
根据场地污染风险分类和分级指标规则,对场地污染风险分类和分级的各指标进行赋分,并运用层次分析法中1~9标度法对各指标的分值进行一致性检验;
通过调整BP神经网络和SVM的各项参数,优化最佳的分类、分级模型,基于BP神经网络进行风险分类,基于SVM进行风险分级,并在MATLAB中完成目标案例的风险分类分级预测,根据风险分类分级的预测结果进行风险管控。
2.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述场地污染风险分类分级案例库中的场地污染风险分类分级案例涉及化学原料和化学制品制造业,有色金属冶炼和压延加工业,黑色金属冶炼和压延加工业,金属制品业,医药制造业,石油、煤炭及其他燃料加工业。
3.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述场地污染风险分类指标体系包括一、二、三级指标,其中,一级指标包括场地概况、污染源、迁移途径、敏感受体;
场地概况下属的二级指标包括所属行业、土地利用规划、场地是否位于工业园区或集聚区、城市等级,其中,所属行业下属的三级指标包括金属制品业、化学原料和化学制品制造业、有色金属冶炼和延压工业、黑色金属冶炼和延压工业、石油和天然气开采业;土地利用规划下属的三级指标包括一类用地、二类用地;场地是否位于工业园区或集聚区下属的三级指标包括是、否;城市等级下属的三级指标包括≤0.5、0.5-1、≥1;
污染源下属的二级指标包括土壤污染物、地下水污染物,其中,土壤污染物下属的三级指标包括Es≥100、50≤Es<100、10≤Es<50、0≤Es<10;地下水污染物下属的三级指标包括Es≥100、50≤Es<100、10≤Es<50,其中,Es表示土壤污染超标总倍数;
迁移途径下属的二级指标包括干湿指数,干湿指数下属的三级指标包括≥1、0.5-1、≤0.5;
敏感受体下属的二级指标包括人群进入可能性、地块周边500m范围内人口、地块距离最近饮用水源,人群进入可能性下属的三级指标包括可能、不可能,地块周边500m范围内人口下属的三级指标包括>500、<500,地块距离最近饮用水源下属的三级指标包括>500m、<500m。
4.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述场地污染风险分级指标体系包括风险分级有关土壤指标体系和风险分级有关地下水指标体系,所述风险分级有关土壤指标体系和风险分级有关地下水指标体系均包括一、二、三级指标,其中,一级指标均包括污染源、污染物迁移途径、敏感受体。
5.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述场地污染风险分类指标规则包括:干湿指数规则、特征污染物规则、土地利用规划、城市等级、包气带的影响和潜水含水层的影响。
6.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述通过调整BP神经网络和SVM的各项参数,优化最佳的分类分级模型,包括:
优化BP神经网络的参数,得到最佳的分类模型,其中,优化BP神经网络的参数包括:
选取Levenberg-Marquardt算法为BP神经网络的训练函数;
选取模型的学习率为0.1;
选取隐含层神经元个数为6。
7.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述通过调整BP神经网络和SVM的各项参数,优化最佳的分类分级模型,包括:
优化SVM参数,得到最佳的分级模型,其中,优化SVM参数包括:核函数为linear、cost为1、gamma为0.08。
8.根据权利要求1所述的基于机器学习的场地地下水和土壤污染风险诊断与管控方法,其特征在于,所述根据风险分类分级的预测结果进行风险管控,具体包括:
针对在产+高风险场地,从源头污染物减排、生产技术升级改造、污染物治理能力提升三方面开展风险管控;
针对在产+低风险场地,进行长期监测,在厂区不同功能区布设长期监测点位并进行样品采集及分析测试;
针对具备修复条件的关闭搬迁+高风险场地,从治理修复和被动修复-衰减两方面开展风险管控;
针对不具备修复条件的关闭搬迁+高风险场地,采取工程控制、长期监测和制度控制相结合的方法;
针对关闭搬迁+低风险场地,采取长期监测和准入管理相结合的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311322998.2A CN117556984A (zh) | 2023-10-13 | 2023-10-13 | 一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311322998.2A CN117556984A (zh) | 2023-10-13 | 2023-10-13 | 一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117556984A true CN117556984A (zh) | 2024-02-13 |
Family
ID=89815533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311322998.2A Pending CN117556984A (zh) | 2023-10-13 | 2023-10-13 | 一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117556984A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10004065A1 (de) * | 1999-12-23 | 2001-07-12 | Jena Geos Ingenieurbuero Gmbh | Verfahren zur Bestimmung der Verhältnismäßigkeit von Maßnahmen der Altlastensanierung beim Boden- und/oder Gewässerschutz |
CN110969345A (zh) * | 2019-11-20 | 2020-04-07 | 北京淖尔科技有限公司 | 一种基于土壤重金属污染途径分析的风险评估方法 |
CN111080097A (zh) * | 2019-12-03 | 2020-04-28 | 中国环境科学研究院 | 农业面源和重金属污染风险综合评估方法 |
CN112150009A (zh) * | 2020-09-25 | 2020-12-29 | 生态环境部环境规划院 | 一种场地污染风险管控与修复模式推荐的方法及装置 |
CN115660419A (zh) * | 2022-10-24 | 2023-01-31 | 北京市生态环境保护科学研究院 | 一套工业园区土壤与地下水污染监测、评估及预警方法 |
CN115907481A (zh) * | 2022-12-26 | 2023-04-04 | 吉林大学 | 一种固化/稳定化重金属场地地下水污染风险评估方法 |
CN116385689A (zh) * | 2023-06-02 | 2023-07-04 | 北京建工环境修复股份有限公司 | 一种场地污染数据的可视化信息管理方法、***及介质 |
-
2023
- 2023-10-13 CN CN202311322998.2A patent/CN117556984A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10004065A1 (de) * | 1999-12-23 | 2001-07-12 | Jena Geos Ingenieurbuero Gmbh | Verfahren zur Bestimmung der Verhältnismäßigkeit von Maßnahmen der Altlastensanierung beim Boden- und/oder Gewässerschutz |
CN110969345A (zh) * | 2019-11-20 | 2020-04-07 | 北京淖尔科技有限公司 | 一种基于土壤重金属污染途径分析的风险评估方法 |
CN111080097A (zh) * | 2019-12-03 | 2020-04-28 | 中国环境科学研究院 | 农业面源和重金属污染风险综合评估方法 |
CN112150009A (zh) * | 2020-09-25 | 2020-12-29 | 生态环境部环境规划院 | 一种场地污染风险管控与修复模式推荐的方法及装置 |
CN115660419A (zh) * | 2022-10-24 | 2023-01-31 | 北京市生态环境保护科学研究院 | 一套工业园区土壤与地下水污染监测、评估及预警方法 |
CN115907481A (zh) * | 2022-12-26 | 2023-04-04 | 吉林大学 | 一种固化/稳定化重金属场地地下水污染风险评估方法 |
CN116385689A (zh) * | 2023-06-02 | 2023-07-04 | 北京建工环境修复股份有限公司 | 一种场地污染数据的可视化信息管理方法、***及介质 |
Non-Patent Citations (2)
Title |
---|
黄国鑫等: "基于机器学习的区域土壤重金属污染风险诊断研究", 《环境工程》, 13 June 2023 (2023-06-13) * |
黄沈发;杨洁;吴健;喻恺;: "城市再开发场地污染风险管控研究及实践", 《环境保护》, no. 1, 29 February 2028 (2028-02-29) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rezaeisabzevar et al. | Landfill site selection using multi criteria decision making: Influential factors for comparing locations | |
Nadiri et al. | Groundwater vulnerability indices conditioned by supervised intelligence committee machine (SICM) | |
Yan et al. | Risk assessment for construction of urban rail transit projects | |
CN109886568B (zh) | 一种危化品道路运输流域水环境风险评估方法及*** | |
Ward et al. | A multi-objective optimisation model for sewer rehabilitation considering critical risk of failure | |
Zhou et al. | Source mapping and determining of soil contamination by heavy metals using statistical analysis, artificial neural network, and adaptive genetic algorithm | |
CN112990739B (zh) | 基于多介质污染特征识别的磷矿区水环境风险评估方法 | |
Malinowska | A fuzzy inference-based approach for building damage risk assessment on mining terrains | |
CN113762784B (zh) | 污水处理***布局规划方案的定量化综合评价方法及*** | |
Bayatzadeh Fard et al. | Use of artificial intelligence techniques to predict distribution of heavy metals in groundwater of Lakan lead-zinc mine in Iran | |
CN111859779A (zh) | 一种燃气管网防第三方施工破坏风险预警的方法及装置 | |
Zhang et al. | Water pipe failure prediction using AutoML | |
CN111612360A (zh) | 一种地下水潜在污染风险源识别方法 | |
Al-Waeli et al. | ANN-based model to predict groundwater salinity: A case study of West Najaf–Kerbala region | |
CN110796375B (zh) | 城市搬迁地园林绿化中生态技术筛选布局方法 | |
CN117556984A (zh) | 一种基于机器学习的场地地下水和土壤污染风险诊断与管控方法 | |
CN115577956A (zh) | 一种工业园区重点污染源地下水污染风险综合评估方法 | |
Arani et al. | Environmental impact assessment of a steel industry development plan using combined method involving Leopold matrix and RIAM | |
Wu et al. | Hse risk assessment of major sewage transport tunnel projects at the construction stage based on the structuralentropy weight method and the cloud model | |
Marceta et al. | Multi-criteria decision-making model to support landfill prioritization: methane risk assessment | |
CN116823011A (zh) | 一种废弃硫铁矿区环境风险分级评价方法 | |
Zhao et al. | Regional ecological risk assessment of chemical industry stress under China's coastal development strategy | |
Cheng et al. | Safety grade evaluation of aqueduct structure based on fuzzy cloud theory analysis | |
Tian et al. | Vulnerability analysis for reusing an old industrial area: a case study of the Winter Olympic Park | |
Pacheco et al. | A framework model to determine groundwater contamination risk based on a L-Matrix of aquifer vulnerability and hazardous activity indices |
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 |