数值模型和线性规划模型相结合的大气污染协同控制方法技术

本发明专利技术公开了一种数值模型和线性规划模型相结合的大气污染协同控制方法，基于源追踪数值模型建立PM

【技术实现步骤摘要】
数值模型和线性规划模型相结合的大气污染协同控制方法


[0001]本专利技术属于气象科学
，尤其涉及一种源追踪数值模型和数学规划模型用以估算大气环境容量的技术。

技术介绍

[0002]空气质量数值模型是当前研究大气污染问题的主要手段。数值模型通过模拟大气中不同空间尺度的物理、化学过程，在技术上能够对历史污染事件进行还原，另外还能对未来污染物浓度变化进行预报预测。当前我国区域复合污染现象突出，为实现对大气污染的控制，已有一些基于数值模型的方法能够根据各污染源的排放强度资料和气象资料，估算污染源对受体的贡献，通过人为调整数值模型内排放源参数，对某种大气污染物进行预报预测，得到该大气污染物的减排控制方案。
[0003]由于近年来我国大气污染格局的变化，颗粒物和臭氧成为影响我国城市和区域空气质量的主要大气污染物，对不同大气污染物（尤其是PM
2.5
和O3）的协同控制成为我国空气质量改善的焦点和打赢蓝天保卫战的关健。另外，线性规划模型虽然在大气环境容量资源问题中的应用具有更强客观性，但其往往作为独立规划模型，在与区域空气质量的污染物协同控制方面未见有研究。
[0004]因此，有必要发展一种数值模型和数学规划模型相结合的区域空气质量达标规划的新方法，。

技术实现思路

[0005]针对现有技术存在的问题，本专利技术提供了一种数值模型和线性规划模型相结合的大气污染协同控制方法，该方法能够在满足城市PM
2.5
和O3同时达标的前提下，给出区域大气环境容量和分城市、分行业区的不同大气污染物减排量，为我国重点城市群长期和短期空气质量管理（PM
2.5
和O3协同控制）提供一定的参考依据。
[0006]为解决上述技术问题，本专利技术采用了以下技术方案：一种数值模型和线性规划模型相结合的大气污染协同控制方法，包括以下步骤：步骤S1，基于空气质量模型CAMx中对颗粒物和臭氧来源的解析，模拟建立研究区域内PM
2.5
、O3及其组分的排放源与受体的非线性响应关系；步骤S2，利用数学规划模型，分别以污染物排放量和区域排放源与受体的响应关系作为输入值，以一次污染源区域允许排放总量最大为目标函数，以PM
2.5
和MDA1 O3浓度同时达标为约束条件，计算得到不同城市、不同行业源中一次污染物优化减排方案的所有可行解；步骤S3，依据步骤S2得到的数学规划模型计算结果，确定污染源减排清单；并利用空气质量模型CAMx重新模拟，判断PM
2.5
和MDA1 O3浓度是否达标，以评估减排方案的有效性和可靠性；步骤S4，筛选满足PM
2.5
和O3同时达标的可行减排方案，给出研究区域不同污染物
最大允许排放量的建议值，以及分区域、分行业排放源的减排比例。
[0007]进一步的，所述一次污染源包括SO2、NO
x
、NH3、VOCs和一次PM
2.5
。
[0008]进一步的，步骤S2中进行区域允许排放总量的计算，采用以下数学规划模型：对于不同物种排放量间的量级差异，使用归一法按对不同物种的排放量加以权重计算不同城市、不同行业源的一次污染源（SO2、NO
x
、NH3、VOCs和一次PM
2.5
）排放总量的最大值，（1）式中，为允许排放总量的最大值，代表不同城市、不同源类、不同物种的控制比例，下角标代表源区城市序号，下角标代表源人为排放源类别序号，下角标代表PM
2.5
或O3的组成物种或对应的一次污染源的序号，为不同物种排放量的权重系数；代表一次污染源的总数，代表受体城市数量，代表考虑的所有源区城市的数量。
[0009]进一步的，步骤S2中，研究区域空气质量目标的约束条件，考虑研究区域的PM
2.5
和O3浓度由背景和非背景贡献的浓度组成，并通过式（2）实现对研究区域空气质量目标的约束，（2）式中，代表单位排放量对受体城市贡献的浓度，代表污染物排放量，代表所有源区自然源对受体城市贡献的浓度，代表目标城市以外区域的人为源对受体城市贡献的浓度，代表初始和边界条件对受体城市贡献的浓度，代表受体城市中PM
2.5
和O3空气质量目标；进一步的，所述不同城市、不同源类、不同物种的控制比例的上下界约束条件符合以下关系式（3）：（3）式中，表示排放源控制比例下界，表示表示排放源控制上界，取值0~100%。。
[0010]有益效果：相对于现有技术，本专利技术能应用于重点城市群PM
2.5
和O3污染的协同控制问题，实现大气环境容量的估算以优化减排建议；本专利技术方法基于源追踪数值模型建立PM
2.5
和O3浓度的区域排放源—受体响应关系，并结合数学规划模型能获得PM
2.5
和O3浓度同时达标情况下的最大允许排放量方案，实现对不同城市、不同行业排放源减排方案的有效性评估，具有更强的客观性。
附图说明
[0011]图1为本专利技术实施例所述某研究地区和源区划分；图2为应用本专利技术方法得到的PM
2.5
及其组分区域排放源—受体响应关系图；图3为应用本专利技术方法得到的MDA1 O3及其组分区域排放源—受体响应关系图；图4为应用本专利技术方法得到不同实验目标函数值的可行解个数以及比较图；图5为应用本专利技术方法得到不同实验一次污染物分行业源排放量和区域排放总量的减排比例比较图；图6为应用本专利技术方法得到Case1实验145种减排方案平均情况下不同城市PM
2.5
和MDA1O3的下降浓度、减排后浓度和下降比例的对比图。
具体实施方式
[0012]下面结合附图并以具体实施例，进一步阐明本专利技术。应理解这些实施例仅用于说明本专利技术而不用于限制本专利技术的范围，在阅读了本专利技术之后，本领域技术人员对本专利技术的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
[0013]本专利技术提供了一种数值模型和线性规划模型相结合的大气污染协同控制方法，包括有以下步骤：第一步：来源解析模块中源追踪区域设置成M个区域，源解析受体为区域内N个城市，源追踪类型设置为5类人为源和自然源，基于空气质量模型CAMx中颗粒物和臭氧来源解析技术（PSAT和OSAT），模拟建立研究区域PM
2.5
、O3及其组分的排放源—受体非线性响应关系；第二步：利用数学规划模型，将第一步模拟建立的污染物排放量和区域排放源—受体响应关系作为输入，以一次污染物（SO2、NO
x
、NH3、VOCs和一次PM
2.5
）区域允许排放总量最大为目标函数，以PM
2.5
和MDA1 O3浓度同时达标为约束条件，计算得到不同城市、不同行业源一次污染物优化减排方案的所有可行解；其中，区域允许排放总量的数学规划模型如下：
①
对于不同物种排放量间的量级差异，使用归一法对不同物种的排放量加以权重计算不同城市、不同行业源一次前本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种数值模型和线性规划模型相结合的大气污染协同控制方法，包括以下步骤：步骤S1，基于空气质量模型CAMx中对颗粒物和臭氧来源的解析，模拟建立研究区域内PM
2.5
、O3及其组分的排放源与受体的非线性响应关系；步骤S2，利用数学规划模型，分别以污染物排放量和区域排放源与受体的响应关系作为输入值，以一次污染源区域允许排放总量最大为目标函数，以PM
2.5
和MDA1 O3浓度同时达标为约束条件，计算得到不同城市、不同行业源中一次污染物优化减排方案的所有可行解；步骤S3，依据步骤S2得到的数学规划模型计算结果，确定污染源减排清单；并利用空气质量模型CAMx重新模拟，判断PM
2.5
和MDA1 O3浓度是否达标，以评估减排方案的有效性和可靠性；步骤S4，筛选满足PM
2.5
和O3同时达标的可行减排方案，给出研究区域不同污染物最大允许排放量的建议值，以及分区域、分行业排放源的减排比例。2.根据权利要求1所述数值模型和线性规划模型相结合的大气污染协同控制方法，其特征在于：所述一次污染源包括SO2、NO
x
、NH3、VOCs和一次PM
2.5
。3.根据权利要求2所述数值模型和线性规划模型相结合的大气污染协同控制方法，其特征在于：步骤S2中进行区域允许排放总量的计算，采用以下数学规划模型：对于...

【专利技术属性】
技术研发人员：王体健，束蕾，
申请(专利权)人：南京智汇环境气象产业研究院有限公司，
类型：发明
国别省市：