1.本发明涉及大气治理技术领域。特别是涉及一种工业园区大气 污染时空解析方法。
背景技术：2.大气污染物排放和pm
2.5
浓度之间呈现显著的非线性变化。目前, 重点区域、城市开展大气污染时空解析工作多利用空气质量模型,空 气质量模型考虑了二次生成的颗粒物,而且可以反映污染物经过传输、 沉降、物理化学作用后的浓度分布特征,对不同排放特征的来源进行 更可靠的解析.但空气质量模型法适用条件复杂、数据种类和数量多, 硬件要求高.近年来统计学模型如机器学习和深度学习应用广泛,相 比空气质量模型,统计模型更容易且高效.但是传统机器学习、深度学 习模型绝大多数只采用空气质量监测数据和气象条件数据进行训练, 缺少污染排放数据从而不能鉴别空气污染形成因素,无法有效管控。
3.
技术实现要素：4.为了解决现有大气污染时空解析方法中存在的上述缺陷，本发 明提出了一种大气污染时空解析方法。
5.本发明利用卷积神经网络和双向长短期记忆网络模拟工业园区 气象、排放源强和污染物浓度之间的变化规律,针对具有代表性的污 染时段,设定区域内不同大气污染排放量,通过敏感性分析研究大气污 染物在不同时刻不同排放情景下的污染削峰效果，获得其污染时空 特征，以削减大气重污染峰值、减缓重污染的发生和发展。
6.本发明通过以下技术方案实现：
7.一种工业园区大气污染时空解析方法，包括如下步骤：
8.1.预处理数据。数据包括气象历史数据、大气污染物浓度历史 数据和排放历史数据，预处理包括异常值处理和缺失值填补；异常 值处理采用3σ法则去除离群值，缺失值填补采用随林森林的链式 方程式多重插补(miceforest，multiple imputation by chainedequations for random forest)，通过一系列迭代的预测模型来插补数 据集中的缺失数据。在每次迭代中，将使用数据集中的其他变量来 估算数据集中的每个指定的变量，这些迭代持续运行，直到满足收 敛为止。
9.2.构建cnn-bilstm模拟模型。模拟模型由卷积神经网络 (cnn)和双向长短期记忆神经网络(bilstm)组成，cnn包括卷积 层、池化层和全连接层；bilstm包括前向lstm层和后向lstm 层。模型输入为气象数据、排放和大气污染物时序数据，输出为未 来1h平均pm
2.5
浓度。卷积层接收输入矩阵
[0010][0011]
t＝{t1,t2,t3,
…
,tn}，代表时间步数量，p＝{p1,p2,p3,
…
,ps}，代 表每一条时间序列的特征数量，通过卷积运算和池化提取其空间特 征，lstm层用于捕获时间序列上的长期依赖，考虑到lstm只能 利用过去的信息进行预测，采用bilstm对过去和未来的信息同时 处理以提高预测精度。cnn提取后的特征集输入到前向lstm层和 后向lstm层中，得到当前时刻其前向层和后向层的隐藏状态和 然后隐藏状态拼接得到神经元输出值o
t
，再经过全连接层后得 到预测值。当前时刻其前向层和后向层的隐藏状态和的计算方 式为：
[0012][0013][0014]
式中，f(
·
)为前向信息提取函数；b(
·
)为后向信息提取函数； wf是前向lstm层的权重矩阵；bf是前向lstm层的偏置；wb是 后向lstm层的权重矩阵；bb是后向lstm层的偏置；
[0015]
当前时刻神经元输出值的计算公式为：
[0016][0017]
式中，为前向lstm层的隐藏状态，为反向lstm层的 隐藏状态，wy是bilstm预测输出的权重矩阵；by是bilstm预测 输出的偏置；
[0018]
3.设计敏感性实验方案。针对代表性的污染时段(pm
2.5
超过中 国3级污染浓度限值的时刻),基于cnn-bilstm模型,对n种主要大 气污染物分别设置s种减排比例,设计ns组对比实验,通过未来 1hpm
2.5
浓度残差作为评价指标衡量减排效果,进而探索工业园区主要 大气污染排放物与空气质量之间的影响关系。
[0019]
4.分析大气污染物时空贡献度，利用ward聚类对敏感性实验 结果在时间和排放量两个维度上进行聚类分析，得到不同时间下工 业园区污染物排放和外部排放对工业园区空气质量的贡献度。
[0020]
本发明利用深度学习模型拟合工业园区内的气象、排放和污染 物浓度关系和设置不同减排指数来模拟不同时段污染物不同排放量 对工业园区空气质量的影响，实现大气污染时空解析。
[0021]
本发明的优点是：利用基于空气污染监测站点的历史数据构建 深度学习模型，通过设置敏感性实验,模拟工业园区气象、排放源强 和污染物浓度之间的变化规律及复合型污染特点，结合ward聚类分 析在时间和排放量两个维度上工业园区污染物排放和外部排
放对工 业园区空气质量的贡献度，实现大气污染时空解析，方法简单，通 用性好。
附图说明
[0022]
图1是cnn-bilstm框架图。
[0023]
图2(a)～图2(f)是不同模型预测值与真实值拟合对比图其 中，图2(a)是svr模型，图2(b)是cnn模型，图2(c)是 lstm模型，图2(d)是bilstm模型，图2(e)是训练数据中不 含污染源排放数据的cnn-bilstm模型，图2(f)是cnn
‑ꢀ
bilstm模型。
[0024]
图3(a)是二氧化硫在不同减排指数下各时段减排效果对比图。
[0025]
图3(b)是氮氧化物在不同减排指数下各时段减排效果对比图。
[0026]
图3(c)是烟尘在不同减排指数下各时段减排效果对比图。
具体实施方式
[0027]
下面结合附图对本发明进行进一步描述。
[0028]
大气污染物浓度数据为2020年12月1日-2020年12月31日濮 阳市工业园区空气自动监测站点逐小时监测数据,包括pm
2.5
、 pm
10
、co、so2及o3质量浓度。其中气象数据来源于濮阳市地面 气象站,主要包括气温、相对湿度、气压、降水、风速风向及露点温 度等。排放数据来源于河南省重点污染源自动监控系统,主要排放物 包括二氧化硫、氮氧化物及烟尘等。
[0029]
一种工业园区大气污染物协同减排评估方法，包括以下步骤：
[0030]
1.预处理数据。数据包括气象历史数据、大气污染物浓度历史 数据和排放历史数据，预处理包括异常值处理和缺失值填补；异常 值处理主要利用3σ法则去除离群值，缺失值填补采用随机森林的链 式方程式多重插补，拟合pm
2.5
浓度和气象因素及其他污染物之间的 关系模型从而进行缺失值填补；气象数据如气温、相对湿度、气 压、降水、风速风向及露点温度等。排放数据如二氧化硫、氮氧化 物及烟尘等大气排放物，大气污染物浓度数据如pm
2.5
、pm
10
、 co、so2及o3等质量浓度；
[0031]
2.构建cnn-bilstm模拟模型。预测模型由cnn和bilstm 神经网络组成，cnn包括卷积层、池化层和全连接层；bilstm包 括前向lstm层和后向lstm层。模型输入为asi、排放和大气污 染物时序数据，输出为未来24h平均pm
2。5
浓度。卷积层通过反向 传播学习到合理的权值,池化层减少空间信息参数,提高运算效率,防 止过拟合。分别在bilstm层和全连接层的末尾加上dropout层随 机丢弃节点,防止过拟合。
[0032]
如图1所示，数据建模包括对原始数据划分时间步和归一化处 理，数据归一化采用min-max标准化方法。预测建模包括cnn卷 积层将时序数据映射为n*s矩阵，其中，n代表时序数据样例数，s 代表数据维度，经过滤波器处理后进入池化层，提取局部特征，减 少冗余信息,将提取后的多组特征向量输入bilstm层、dropout层 和全连接层中，得到输出结果。此外,为了增强cnn-bilstm网络 模型的泛化性和鲁棒性,选择adam作为优化算法。
[0033]
不同模型的拟合效果和评价指标如图2所示,图2(a)～(f)分别代 表svr模型、cnn模型、lstm模型、bilstm模型、训练数据中 不含污染源排放数据的cnn-bilstm模型、cnn-bilstm模型,从 可决系数来看,训练数据包含污染源排放数据的cnn-bilstm模 型拟合效果最好,评价指标最优.其r2值为0.93,预测值与真实值的相 关性极高,相比不包含污染源
排放数据的传统模型,精度提高4％.证明 融合污染源排放数据后cnn-bilstm模型有更高的预测精度,表明 cnn-bilstm模型利用工业园区气象、排放量和空气质量数据成功 模拟了pm
2.5
浓度的变化趋势，因此基于该模型的分析结果具有可靠 性。
[0034]
3.设计敏感性实验方案。针对代表性的污染时段(pm
2.5
超过中 国3级污染浓度限值的时刻),基于cnn-bilstm模型,对n种主要大 气污染物分别设置s种减排比例,设计ns组对比实验,通过未来 1hpm
2.5
浓度残差作为评价指标衡量减排效果,进而探索工业园区主要 大气污染排放物之间的相互影响关系。
[0035]
4.分析大气污染物时空贡献度，利用ward聚类对敏感性实验 结果在时间和排放量两个维度上进行聚类分析，得到不同时间下工 业园区污染物排放和外部排放对工业园区空气质量的贡献度。
[0036]
如图3所示,图3a到图3c分别代表二氧化硫、氮氧化物和烟尘 在不同时刻不同减排指数下未来1hpm
2.5
浓度预测值与真实值残差对 比,利用ward聚类在时间和减排指数两个维度上对结果进行聚类分 析.结果表明,在减排指数维度上,二氧化硫最优减排指数为10,其次为 4和7；氮氧化物最优减排指数为6,其次为3；烟尘最优减排指数为 1和10.在时间维度上,二氧化硫和氮氧化物在6-7时和17-18时的减 排效果最好,其次是15-17时；烟尘在15-16时减排效果最好,其次为 6-7和17-18时.这表明大气污染物排放量削减会降低工业园区pm
2.5
浓度,不同大气污染物在空气中的响应程度与排放量呈非线性关 系,pm
2.5
浓度在6-7时和15-18时减排下降幅度最大,这是由于逆温现 象一般出现在夜间和早晨,冬季傍晚因为地面温度急剧降低,近地面的 下层大气温度降低,上层空气降温相对较慢,这种下冷上热的逆温层导 致空气无法上下对流,污染物不易扩散；随着太阳升高,地面温度逐渐 上升,静稳层结被打破,污染物开始扩散；因此早晨和傍晚大气自净能 力最强,此时减排效果显著.
[0037]
空间上工业园区内二氧化硫、氮氧化物和烟尘在减排指数为 10、10、6时pm
2.5
浓度下降幅度最大,对pm
2.5
浓度的平均最大贡献 度为22％、20％、23％,可见外来排放占35％左右；时间上由于逆温 效应，在6-7时和15-18时减排污染削峰效果最优，pm
2.5
下降比例 分别为24％和23％。
[0038]
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列 举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式， 本发明的保护范围也及于本领域技术人员根据本发明构思所能够想 到的等同技术手段。
技术特征：1.一种工业园区大气污染时空解析方法，包括如下步骤：s1.预处理数据；数据包括气象历史数据、大气污染物浓度历史数据和排放历史数据，预处理包括异常值处理和缺失值填补；异常值处理采用3σ法则去除离群值，缺失值填补采用随林森林的链式方程式多重插补(miceforest，multiple imputation by chained equations for random forest)，通过一系列迭代的预测模型来插补数据集中的缺失数据；在每次迭代中，将使用数据集中的其他变量来估算数据集中的每个指定的变量，这些迭代持续运行，直到满足收敛为止；s2.构建cnn-bilstm模拟模型；模拟模型由卷积神经网络(cnn)和双向长短期记忆神经网络(bilstm)组成，cnn包括卷积层、池化层和全连接层；bilstm包括前向lstm层和后向lstm层；模型输入为气象数据、排放和大气污染物时序数据，输出为未来1h平均pm
2.5
浓度；卷积层接收输入矩阵卷积层接收输入矩阵t＝{t1,t2,t3,
…
,t
n
}，代表时间步数量，p＝{p1,p2,p3,
…
,p
s
}，代表每一条时间序列的特征数量，通过卷积运算和池化提取其空间特征，lstm层用于捕获时间序列上的长期依赖，考虑到lstm只能利用过去的信息进行预测，采用bilstm对过去和未来的信息同时处理以提高预测精度；cnn提取后的特征集输入到前向lstm层和后向lstm层中，得到当前时刻其前向层和后向层的隐藏状态和然后隐藏状态拼接得到神经元输出值o
t
，再经过全连接层后得到预测值；当前时刻其前向层和后向层的隐藏状态和的计算方式为：的计算方式为：式中，f(
·
)为前向信息提取函数；b(
·
)为后向信息提取函数；w
f
是前向lstm层的权重矩阵；b
f
是前向lstm层的偏置；w
b
是后向lstm层的权重矩阵；b
b
是后向lstm层的偏置；当前时刻神经元输出值的计算公式为：式中，为前向lstm层的隐藏状态，为反向lstm层的隐藏状态，w
y
是bilstm预测输出的权重矩阵；b
y
是bilstm预测输出的偏置；s3.设计敏感性实验方案；针对代表性的污染时段,基于cnn-bilstm模型,对n种主要大气污染物分别设置s种减排比例,设计n
s
组对比实验,通过未来1hpm
2.5
浓度残差作为评价指标衡量减排效果,进而探索工业园区主要大气污染排放物与空气质量之间的影响关系；s4.分析大气污染物时空贡献度，利用ward聚类对敏感性实验结果在时间和排放量两
个维度上进行聚类分析，得到不同时间下工业园区污染物排放和外部排放对工业园区空气质量的贡献度。2.如权利要求1所述的一种工业园区大气污染时空解析方法，其特征在于：步骤s3所述的代表性的污染时段是pm
2.5
超过中国3级污染浓度限值的时刻。
技术总结一种工业园区大气污染时空解析方法，包括：S1.预处理气象历史数据、大气污染物浓度历史数据和排放历史数据，包括异常值处理和缺失值填补；S2.构建CNN-BILSTM模拟模型；S3.设计敏感性实验方案；针对代表性的污染时段,基于CNN-BILSTM模型,对N种主要大气污染物分别设置s种减排比例,设计N
技术研发人员：董红召 廖世凯 金灿 方浩杰 王桢 余滔滔
受保护的技术使用者：浙江工业大学
技术研发日：2022.06.15
技术公布日：2022/11/1
