大气污染多尺度融合溯源方法、装置和可读存储介质与流程

1.本发明涉及大气污染防治与溯源技术领域，具体而言，涉及一种大气污染 多尺度融合溯源方法、装置和可读存储介质。


背景技术：

2.随着我国经济快速发展和工业化、城镇化进程不断加速，能源消耗增加， 出现了一系列的大气环境问题。大气污染防治最有效的方法就是控制内因，核 心是找准大气污染源头，实现靶向治理，提高控制效率。
3.大气污染溯源通过动力学模式追溯排放源及贡献率的时空分布，结合排放 源成分解析结果和排放源清单可追溯到具体的排放源，并实现对预报结果的同 步溯源。源模型法是以不同尺度数值模式方法定量描述大气污染物从源到受体 所经历的物理化学过程，定量估算不同地区和不同类别污染源排放对环境空气 中颗粒物的贡献。当前基于源模型的大气污染溯源技术主要包括空气质量模型 cmaq(community multiscale air quality，第三代空气质量预报和评估装置 的模式)、camx(comprehensive air quality model with extensions，第三代 三维空气质量模型)；颗粒物污染精准溯源hysplit(hybrid single particle lagrangian integrated trajectory model，拉格朗日混合单粒子轨道模型)等。 但这些模型都只能做到大尺度区域污染溯源，需要气象模型wrf(the weather research and forecasting model，天气预报模式)和排放源模型smoke(大气 排放源清单处理模型)的支持，存在操作难度大、计算量大等问题，且无法实 现多尺度，特别是小尺度精准溯源。


技术实现要素：

4.本发明旨在解决上述技术问题的至少之一。
5.为此，本发明的第一目的在于提供一种大气污染多尺度融合溯源方法。
6.本发明的第二目的在于提供一种大气污染多尺度融合溯源装置。
7.本发明的第三目的在于提供一种大气污染多尺度融合溯源装置。
8.本发明的第四目的在于提供一种可读存储介质。
9.为实现本发明的第一目的，本发明的技术方案提供了一种大气污染多尺度 融合溯源方法，包括：基于监控站点的大气污染浓度值超标，得到大气污 染浓度值超标信息，获取大气污染浓度值超标前第一时间段内的气象数据； 根据气象数据，模拟计算得到污染来源轨迹和主要轨迹贡献比率；根据污 染来源轨迹，计算得到污染潜在源区和污染潜在源区贡献率；根据污染潜在 源区，模拟计算得到源区内企业对监控站点的污染浓度贡献值；根据主要 轨迹贡献比率对污染潜在源区贡献率进行融合修正，得到修正后的潜在源 区贡献率；基于修正后的污染潜在源区贡献率对企业对监控站点的污染浓 度贡献值进行修正，得到修正后的企业污染浓度贡献值。
10.本技术方案的一种多尺度逐层递进的溯源方法，从大方向到小区域再 到企业点
源逐层精细化溯源，可以满足不同环境治理决策需要。
11.另外，本发明提供的技术方案还可以具有如下附加技术特征：
12.上述技术方案中，根据气象数据，模拟计算得到污染来源轨迹和主要轨 迹贡献比率，具体包括：根据气象数据，以监控站点位置为起点，采用拉 格朗日大气扩散模型进行后向轨迹分析，计算污染来源轨迹；对污染来源 轨迹进行聚类，得到主要轨迹；基于主要轨迹，得到每个主要轨迹贡献比 率。
13.本技术方案中，通过第一重大尺度溯源可以看出污染来源主要方向、输送 轨迹及主要轨迹贡献比率，可以对污染来源进行大方向上的溯源。
14.上述任一技术方案中，采用拉格朗日大气扩散模型进行后向轨迹分析， 计算污染来源轨迹，具体包括：假设污染来源为质点，质点的轨迹即为污染 来源轨迹，计算公式为：
15.p'(t+δt)＝p(t)+v(p,t)
·
δt
16.p(t+δt)＝p(t)+0.5[v(p,t)+v(p',t+δt)]δt
[0017]
其中，p(t)表示t时刻质点位置，p'(t+δt)表示t+δt时刻质点位置，p(t+δt) 表示t+δt时刻质点位置通过矢量速度在时间和空间上线性插值得出，v(p,t) 表示质点在p位置和t时刻的速度，v(p',t+δt)表示质点在p
′
位置和t+δt时刻 的速度，δt表示时间差，t表示时间。
[0018]
本技术方案给出了污染来源轨迹的具体计算公式，计算过程简单，容易 实现。
[0019]
上述任一技术方案中，根据污染来源轨迹，计算得到污染潜在源区和污 染潜在源区贡献率，具体包括：根据监控站点确定污染溯源范围；对污染 溯源范围区域进行网格划分；获取监控站点过去第二时间内的污染监测浓 度值；计算得到污染潜在源区贡献率，具体包括：统计每个网格内经过的 污染来源轨迹总数，得到第一轨迹个数；每条污染来源轨迹对应一个监控 站点的污染物浓度，对污染物浓度设定阈值，统计每个网格内污染轨迹浓 度值超过阈值的轨迹个数，得到第二轨迹个数；计算第二轨迹个数与第一 轨迹个数比值：
[0020][0021]
其中，m
ij
表示超过设定阈值的污染来源轨迹在网格(i,j)内被截获的第 二轨迹个数，n
ij
表示所有污染来源轨迹被网格(i,j)截获的第一轨迹个数；
[0022]
获取权重因子：
[0023][0024]
其中，n
ave
表示溯源范围内所有网格的平均轨迹个数，当3n
ave
＜n
ij
时，w
ij
＝1， 当1.5n
ave
＜n
ij
＜3n
ave
时，w
ij
＝0.7，当n
ave
＜n
ij
＜3n
ave
时，w
ij
＝0.42，当n
ij
＜n
ave
时，w
ij
＝0.05；
[0025]
根据权重因子和比值，得到污染潜在源区的加权贡献率：
[0026][0027]
其中，表示网格(i,j)对监控站点的污染贡献率，即污染潜在源区贡 献率。
[0028]
本技术方案通过对污染溯源范围进行网格划分，完成中尺度潜在区域溯 源，满足不同环境治理决策需要。
[0029]
上述任一技术方案中，大气污染多尺度融合溯源方法，还包括：设定不 同污染潜在源区贡献率对应的颜色；根据污染潜在源区贡献率，对每个网 格进行颜色填充，进行可视化展示。
[0030]
本实施中，通过对网格进行颜色填充，可以更加直观的显示不同污染 潜在源区贡献率，满足不同环境治理决策需要。
[0031]
上述任一技术方案中，大气污染多尺度融合溯源方法，还包括：
[0032]
获取每个网格的污染物对监控站点的平均权重浓度值：
[0033][0034]
其中，c
ij
表示网格(i,j)的平均权重浓度值，c
l
表示污染来源轨迹l经过 网格(i,j)时的大气污染浓度值，t
ijl
表示污染来源轨迹l经过网格(i,j)的停留时 间，m表示污染来源轨迹l的数量；
[0035]
根据平均权重浓度值和权重因子，得到相对权重浓度值：
[0036][0037]
其中，表示网格(i,j)的相对权重浓度值；
[0038]
根据相对权重浓度值和比值，更新污染潜在源区贡献率。
[0039]
本技术方案弥补潜在源区贡献率分析无法区分相同贡献值的网格对报 警点浓度贡献大小的缺陷。
[0040]
上述任一技术方案中，根据污染潜在源区，模拟计算得到源区内企业对 监控站点的污染浓度贡献值，具体包括：根据污染潜在源区内的企业排放 数据、气象数据和地形数据，利用非稳态拉格朗日烟团模型装置，计算企 业对监控站点的污染浓度贡献值。
[0041]
本技术方案为小尺度溯源，通过企业排放数据、气象数据和地形数据， 得到企业对监控站点的污染浓度贡献值，采用非稳态拉格朗日烟团模型系 统，计算过程简单，准确率高。
[0042]
上述任一技术方案中，根据主要轨迹贡献比率对污染潜在源区贡献率进 行融合修正，得到修正后的潜在源区贡献率；基于修正后的污染潜在源区 贡献率对企业对监控站点的污染浓度贡献值进行修正，得到修正后的企业 污染浓度贡献值，具体包括：
[0043]
确定网格(i,j)落入的污染来源轨迹l，采用污染来源轨迹l的主要轨迹 贡献比率，对污染来源轨迹l的污染潜在源区贡献率进行修正：
[0044][0045]
其中，表示修正后的污染潜在源区贡献率，表示第l条主 要轨迹上的原始污染潜在源区贡献率，表示第l条主要轨迹的贡献比率；
[0046]
确定企业落入的网格，采用修正后的污染潜在源区贡献率，对企业对 监控站点的污染浓度贡献值进行修正，得到修正后的企业污染浓度贡献值：
[0047][0048]
其中，表示修正后的第k个企业的贡献浓度值，表示网格(i,j) 区域内第k个企业的原始浓度贡献值。
[0049]
本技术方案给出了进行多尺度融合的具体计算方法，可以弥补由于单 个模型计算误差、奇异点和数据不完整等带来的计算结果不准确问题，提 升潜在区域贡献率和企业污染浓度贡献度的准确性。
[0050]
上述任一技术方案中，大气污染多尺度融合溯源方法，还包括：对企业 的修正后的企业污染浓度贡献值进行排名，进行可视化展示。
[0051]
本实施对得到的修正后的企业污染浓度贡献值进行排名，通过排名可 以简单直观的找出具体的企业污染源头，对企业污染源头进行直观显示。
[0052]
为实现本发明的第二目的，本发明的技术方案提供了一种大气污染多尺 度融合溯源装置，包括：数据获取模块、第一获取模块、第二获取模块、第 三获取模块和第四获取模块；数据获取模块基于监控站点的大气污染浓度 值超标，得到大气污染浓度值超标信息，获取大气污染浓度值超标前第一 时间内的气象数据；第一获取模块根据气象数据，模拟计算得到污染来源 轨迹和主要轨迹贡献比率；第二获取模块根据污染来源轨迹，计算得到污 染潜在源区和污染潜在源区贡献率；第三获取模块根据污染潜在源区，模 拟计算得到企业对监控站点的污染浓度贡献值；第四获取模块根据主要轨 迹贡献比率对潜在源区贡献率进行融合修正，得到修正后的潜在源区贡献 率；基于修正后的潜在源区贡献率对企业对监控站点的污染浓度贡献值进 行融合修正，得到修正后的企业污染浓度贡献值。
[0053]
本技术方案充分利用不同尺度溯源结果进行融合分析，使用大尺度溯 源结果对下一级尺度溯源结果进行修正，可以弥补由于单个模型计算误差、 奇异点和数据不完整等带来的计算结果不准确问题，提升潜在源区贡献率 和企业污染浓度贡献度的准确性。
[0054]
为实现本发明的第三目的，本发明的技术方案提供了一种大气污染多尺度 融合溯源装置，包括：存储器和处理器，存储器存储有程序或指令，处理 器执行程序或指令；其中，处理器在执行程序或指令时，实现如本发明任 一技术方案的大气污染多尺度融合溯源方法的步骤。
[0055]
本技术方案提供的大气污染多尺度融合溯源装置实现如本发明任一技术 方案的大气污染多尺度融合溯源方法的步骤，因而其具有如本发明任一技术 方案的大气污染多尺度融合溯源方法的全部有益效果，在此不再赘述。
‘
[0056]
为实现本发明的第四目的，本发明的技术方案提供了一种可读存储介质， 可读存储介质存储有程序或指令，程序或指令被执行时，实现上述任一技 术方案的大气污染多尺度融合溯源方法的步骤。
[0057]
本技术方案提供的可读存储介质实现如本发明任一技术方案的大气污 染多尺度融合溯源方法的步骤，因而其具有如本发明任一技术方案的大气污 染多尺度融合溯源方法的全部有益效果，在此不再赘述。
[0058]
本发明的附加方面和优点将在下面的描述部分中变得明显，或通过本发明 的实
践了解到。
附图说明
[0059]
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中 将变得明显和容易理解，其中：
[0060]
图1为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之一；
[0061]
图2为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之二；
[0062]
图3为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之三；
[0063]
图4为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之四；
[0064]
图5为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之五；
[0065]
图6为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之六；
[0066]
图7为根据本发明一个实施例的大气污染多尺度融合溯源装置组成示意 图；
[0067]
图8为根据本发明一个实施例的大气污染多尺度融合溯源装置组成示意 图；
[0068]
图9为根据本发明一个实施例的大气污染多尺度融合溯源方法原理示意 框图；
[0069]
图10为根据本发明一个实施例的污染来源轨迹图示意图；
[0070]
图11为根据本发明一个实施例的污染来源轨迹聚类结果展示图；
[0071]
图12为根据本发明一个实施例的溯源区域确定及网格划分示意图；
[0072]
图13为根据本发明一个实施例的潜在源区分析可视化展示示意图；
[0073]
图14为根据本发明一个实施例的小尺度点源溯源结果展示；
[0074]
图15为根据本发明一个实施例的大气污染多尺度融合溯源方法流程示意 图之七。
[0075]
其中，图7至图9中附图标记与部件名称之间的对应关系为：
[0076]
100：大气污染多尺度融合溯源装置，110：数据获取模块，120：第一 获取模块，130：第二获取模块，140：第三获取模块，150：第四获取模块， 200：大气污染多尺度融合溯源装置，210：存储器，220：处理器，300：监控 站点报警时间、地点、浓度值和气象数据，310：大尺度溯源，320：中尺度溯 源，330：小尺度溯源，340：多尺度融合。
具体实施方式
[0077]
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和 具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情 况下，本技术的实施例及实施例中的特征可以相互组合。
[0078]
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发 明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范 围并不受下面公开的具体实施例的限制。
[0079]
下面参照图1至图15描述本发明一些实施例的大气污染多尺度融合溯源 方法、装置和可读存储介质。
[0080]
实施例1：
[0081]
如图1所示，本实施例提供了一种大气污染多尺度融合溯源方法，包括 以下步骤：
[0082]
步骤s102，基于监控站点的大气污染浓度值超标，得到大气污染浓度 值超标信
息，获取大气污染浓度值超标前第一时间段内的气象数据；
[0083]
步骤s104，根据气象数据，模拟计算得到污染来源轨迹和主要轨迹贡 献比率；
[0084]
步骤s106，根据污染来源轨迹，计算得到污染潜在源区和污染潜在源 区贡献率；
[0085]
步骤s108，根据污染潜在源区，模拟计算得到源区内企业对监控站点 的污染浓度贡献值；
[0086]
步骤s110，根据主要轨迹贡献比率对污染潜在源区贡献率进行融合修 正，得到修正后的潜在源区贡献率；基于修正后的污染潜在源区贡献率对 企业对监控站点的污染浓度贡献值进行修正，得到修正后的企业污染浓度 贡献值。
[0087]
相关技术的大气污染多尺度融合溯源方法，一般采用大方向大区域溯 源或小尺度企业关联分析溯源。无法满足不同的治理环境决策需求。
[0088]
本实施例中，首先，获取大气污染浓度值超标的监测站点的气象数据， 根据气象数据进行大尺度溯源，即获取污染来源轨迹和主要轨迹贡献比率， 然后，进行中尺度溯源，得到污染潜在源区和污染潜在源区贡献率，第三， 进行小尺度溯源，得到企业对监控站点的污染浓度贡献值，最后，进行多 尺度融合，采用主要轨迹贡献比率，对污染潜在源区贡献率进行修正，采 用修正后的污染潜在源区贡献率，对企业对监控站点的污染浓度贡献值进 行修正，实现多尺度融合修正，得到修正后的企业污染浓度贡献值。
[0089]
本实施例中，第一时间段根据实际需求进行设定，举例而言，可以为 一周。
[0090]
本实施例的一种多尺度逐层递进的溯源方法，从大方向到小区域再到 企业点源逐层精细化溯源，可以满足不同环境治理决策需要。
[0091]
本实施例充分利用不同尺度溯源结果进行融合分析，使用大尺度溯源 结果对下一级尺度溯源结果进行修正，可以弥补由于单个模型计算误差、 奇异点和数据不完整等带来的计算结果不准确问题，提升潜在源区贡献率 企业污染浓度贡献度的准确性。
[0092]
实施例2：
[0093]
如图2所示，本实施例提供了一种大气污染多尺度融合溯源方法，除上 述实施例的技术特征以外，本实施例进一步地包括了以下技术特征：
[0094]
根据气象数据，模拟计算得到污染来源轨迹和主要轨迹贡献比率，具 体包括以下步骤：
[0095]
步骤s202，根据气象数据，以监控站点位置为起点，采用拉格朗日大 气扩散模型进行后向轨迹分析，计算污染来源轨迹；
[0096]
步骤s204，对污染来源轨迹进行聚类，得到主要轨迹；
[0097]
步骤s206，基于主要轨迹，得到每个主要轨迹贡献比率。
[0098]
本实施例中，根据对气象数据进行大尺度溯源，采用拉格朗日大气扩散 模型，可以获取准确度更高的污染来源轨迹，计算方法简单，易于实现。
[0099]
本实施例中，通过第一重大尺度溯源可以看出污染来源主要方向、输送轨 迹及主要轨迹贡献比率，可以对污染来源进行大方向上的溯源，满足不同环 境治理决策需要。
[0100]
实施例3：
[0101]
本实施例提供了一种大气污染多尺度融合溯源方法，除上述实施例的技 术特征以外，本实施例进一步地包括了以下技术特征：
[0102]
采用拉格朗日大气扩散模型进行后向轨迹分析，计算污染来源轨迹， 具体包括：
假设污染来源为质点，质点的轨迹即为污染来源轨迹，计算公式 为：
[0103]
p'(t+δt)＝p(t)+v(p,t)
·
δt
[0104]
p(t+δt)＝p(t)+0.5[v(p,t)+v(p',t+δt)]δt
[0105]
其中，p(t)表示t时刻质点位置，p'(t+δt)表示t+δt时刻质点位置，p(t+δt) 表示t+δt时刻质点位置通过矢量速度在时间和空间上线性插值得出，v(p,t) 表示质点在p位置和t时刻的速度，v(p',t+δt)表示质点在p
′
位置和t+δt时刻 的速度，δt表示时间差，t表示时间。
[0106]
本实施例给出了污染来源轨迹的具体计算公式，计算过程简单，容易实 现。
[0107]
实施例4：
[0108]
如图3所示，本实施例提供了一种大气污染多尺度融合溯源方法，除上 述实施例的技术特征以外，本实施例进一步地包括了以下技术特征：
[0109]
根据污染来源轨迹，得到污染潜在源区和污染潜在源区贡献率，具体包 括以下步骤：
[0110]
步骤s302，根据监控站点确定污染溯源范围；
[0111]
步骤s304，对污染溯源范围区域进行网格划分；
[0112]
步骤s306，获取监控站点过去第二时间内的污染监测浓度值；
[0113]
步骤s308，获取污染潜在源区贡献率。
[0114]
获取污染潜在源区贡献率，具体包括：：统计每个网格内经过的污染 来源轨迹总数，得到第一轨迹个数；每条污染来源轨迹对应一个监控站点 的污染物浓度，对污染物浓度设定阈值，统计每个网格内污染轨迹浓度值 超过阈值的轨迹个数，得到第二轨迹个数；计算第二轨迹个数与第一轨迹 个数比值：
[0115][0116]
其中，m
ij
表示超过设定阈值的污染来源轨迹在网格(i,j)内被截获的第 二轨迹个数，n
ij
表示所有污染来源轨迹被网格(i,j)截获的第一轨迹个数；
[0117]
获取权重因子：
[0118][0119]
其中，n
ave
表示溯源范围内所有网格的平均轨迹个数，当3n
ave
＜n
ij
时，w
ij
＝1， 当1.5n
ave
＜n
ij
＜3n
ave
时，w
ij
＝0.7，当n
ave
＜n
ij
＜3n
ave
时，w
ij
＝0.42，当n
ij
＜n
ave
时，w
ij
＝0.05；
[0120]
根据权重因子和比值，得到加权污染潜在源区贡献率：
[0121][0122]
其中，表示网格(i,j)对监控站点的污染贡献率，即污染潜在源区贡 献率。
[0123]
本实施例中，首先，根据监控站点划定污染溯源范围，范围的大小根据 实际情况选定，举例而言，可以是100公里
×
100公里等。然后，将选定 的污染溯源范围划分为网格，
网格大小根据实际需要动态设置，举例而言， 可以是5公里
×
5公里等。第三，获取监控站点过去第二时间内的污染监 测浓度值，第二时间可以根据实际情况进行设定，举例而言，可以为过去 一周以小时为分辨率的污染监测浓度值。最后，得到污染潜在源区贡献率。
[0124]
本实施例通过对污染溯源范围进行网格划分，完成中尺度潜在区域溯 源，满足不同环境治理决策需要。
[0125]
实施例5：
[0126]
如图4所示，本实施例提供了一种大气污染多尺度融合溯源方法，除上 述实施例的技术特征以外，本实施例进一步地包括了以下技术特征：
[0127]
大气污染多尺度融合溯源方法，还包括以下步骤：
[0128]
步骤s402，设定不同污染潜在源区贡献率对应的颜色；
[0129]
步骤s404，根据污染潜在源区贡献率，对每个网格进行颜色填充，进 行可视化展示。
[0130]
本实施中，污染溯源范围被划分为网格，根据每个网格中污染潜在源区 贡献率不同，对网格进行颜色填充，填充颜色后的网格能够进行可视化展 示，通过对网格进行颜色填充，可以更加直观的显示不同污染潜在源区贡 献率，满足不同环境治理决策需要。
[0131]
实施例6：
[0132]
本实施例提供了一种大气污染多尺度融合溯源方法，除上述实施例的技 术特征以外，本实施例进一步地包括了以下技术特征：
[0133]
大气污染多尺度融合溯源方法，还包括：
[0134]
获取每个网格的污染物对监控站点的平均权重浓度值：
[0135][0136]
其中，c
ij
表示网格(i,j)的平均权重浓度值，c
l
表示污染来源轨迹l经过 网格(i,j)时的大气污染浓度值，t
ijl
表示污染来源轨迹l经过网格(i,j)的停留时 间，m表示污染来源轨迹l的数量；
[0137]
根据平均权重浓度值和权重因子，得到相对权重浓度值：
[0138][0139]
其中，表示网格(i,j)的相对权重浓度值；
[0140]
根据相对权重浓度值和比值，更新污染潜在源区贡献率。
[0141]
本实施例中，获取每个网格的平均权重浓度值，再根据平均权重浓度值得 到相对权重浓度值，通过上述计算，弥补潜在源区贡献率分析无法区分相 同贡献值的网格对报警点浓度贡献大小的缺陷。
[0142]
实施例7：
[0143]
如图5所示，本实施例提供了一种大气污染多尺度融合溯源方法，除上 述实施例的技术特征以外，本实施例进一步地包括了以下技术特征：
[0144]
根据污染潜在源区，模拟计算得到源区内企业对监控站点的污染浓度 贡献值，具
体包括以下步骤：
[0145]
步骤s502，根据污染潜在源区内的企业排放数据、气象数据和地形数 据，利用非稳态拉格朗日烟团模型系统，计算企业对监控站点的污染浓度 贡献值。
[0146]
本实施例为小尺度溯源，通过企业排放数据、气象数据和地形数据，得 到企业对监控站点的污染浓度贡献值，采用非稳态拉格朗日烟团模型系统， 计算过程简单，准确率高。
[0147]
实施例8：
[0148]
本实施例提供了一种大气污染多尺度融合溯源方法，除上述实施例的技 术特征以外，本实施例进一步地包括了以下技术特征：
[0149]
根据主要轨迹贡献比率对污染潜在源区贡献率进行融合修正，得到修 正后的潜在源区贡献率；基于修正后的污染潜在源区贡献率对企业对监控 站点的污染浓度贡献值进行修正，得到修正后的企业污染浓度贡献值，具 体包括：
[0150]
确定网格(i,j)落入的污染来源轨迹l，采用污染来源轨迹l的主要轨迹 贡献比率，对污染来源轨迹l的污染潜在源区贡献率进行修正：
[0151][0152]
其中，表示修正后的污染潜在源区贡献率，表示第l条主 要轨迹上的原始污染潜在源区贡献率，表示第l条主要轨迹的贡献比率；
[0153]
确定企业落入的网格，采用修正后的污染潜在源区贡献率，对企业对 监控站点的污染浓度贡献值进行修正，得到修正后的企业污染浓度贡献值：
[0154][0155]
其中，表示修正后的第k个企业的贡献浓度值，表示网格(i,j) 区域内第k个企业的原始浓度贡献值。
[0156]
本实施例给出了进行多尺度融合的具体计算方法，采用主要轨迹贡献 比率，对污染潜在源区贡献率进行修正，采用修正后的污染潜在源区贡献 率，对企业对监控站点的污染浓度贡献值进行修正，实现多尺度融合修正， 得到修正后的企业污染浓度贡献值，可以弥补由于单个模型计算误差、奇 异点和数据不完整等带来的计算结果不准确问题，提升企业污染浓度的准 确性。
[0157]
实施例9：
[0158]
如图6所示，本实施例提供了一种大气污染多尺度融合溯源方法，除上 述实施例的技术特征以外，本实施例进一步地包括了以下技术特征：
[0159]
大气污染多尺度融合溯源方法，还包括以下步骤：
[0160]
步骤s602，对企业的修正后的企业污染浓度贡献值进行排名，进行可 视化展示。
[0161]
本实施对得到的修正后的企业污染浓度贡献值进行排名，通过排名可 以简单直观的找出具体的企业污染源头，对企业污染源头进行直观显示。
[0162]
实施例10：
[0163]
如图7所示，本实施例提供了一种大气污染多尺度融合溯源装置100， 包括：数据
获取模块110、第一获取模块120、第二获取模块130、第三获 取模块140和第四获取模块150；数据获取模块110基于监控站点的大气 污染浓度值超标，得到大气污染浓度值超标信息，获取大气污染浓度值超 标前第一时间内的气象数据；第一获取模块120根据气象数据，模拟计算 得到污染来源轨迹和主要轨迹贡献比率；第二获取模块130根据污染来源 轨迹，计算得到污染潜在源区和污染潜在源区贡献率；第三获取模块140 根据污染潜在源区，模拟计算得到企业对监控站点的污染浓度贡献值；第 四获取模块150根据主要轨迹贡献比率对潜在源区贡献率进行融合修正， 得到修正后的潜在源区贡献率；基于修正后的潜在源区贡献率对企业对监 控站点的污染浓度贡献值进行融合修正，得到修正后的企业污染浓度贡献 值。
[0164]
本实施例的大气污染多尺度融合溯源装置100，从大方向到小区域再到 企业点源逐层精细化溯源，可以满足不同环境治理决策需要。
[0165]
本实施例充分利用不同尺度溯源结果进行融合分析，使用大尺度溯源 结果对下一级尺度溯源结果进行修正，可以弥补由于单个模型计算误差、 奇异点和数据不完整等带来的计算结果不准确问题，提升潜在源区贡献率 企业污染浓度贡献度的准确性。
[0166]
实施例11：
[0167]
如图8所示，本实施例提供了一种大气污染多尺度融合溯源装置200， 包括：存储器210和处理器220，存储器210存储有程序或指令，处理器 220执行程序或指令；其中，处理器220在执行程序或指令时，实现如本 发明任一实施例的大气污染多尺度融合溯源方法的步骤。
[0168]
实施例12：
[0169]
本实施例提供了一种可读存储介质，可读存储介质存储有程序或指令， 程序或指令被处理器220执行时，实现上述任一实施例的大气污染多尺度融 合溯源方法的步骤。
[0170]
实施例13：
[0171]
本实施例提供了一种大气污染多尺度融合溯源方法，针对某监控站污染 报警进行3重尺度溯源并进行融合分析：
[0172]
第1重尺度溯源(即大尺度溯源)，根据气象数据，以监控站点位置为起 点，利用hysplit(hybrid single
‑
particle lagrangian integrated trajectory model，拉格朗日大气扩散模型)进行后向轨迹分析，计算出污染来源轨迹， 并聚类出主要轨迹，得出各主要轨迹的贡献比率(第l条轨迹的贡献率)。通 过大尺度溯源可以看出污染来源主要方向、输送轨迹及其贡献率。
[0173]
第2重尺度溯源(即中尺度溯源)，基于第1重尺度溯源得到的污染来源 轨迹，根据概率论和统计原理，计算出污染潜在源区、潜在源区的贡献率 (网格(i,j)的贡献率)和权重浓度值，并进行潜在源区排名。通过中尺度溯源 可以找出对报警点污染超标造成影响的潜在源区，并计算出潜在源区对报警点 的贡献率及重权浓度值。
[0174]
第3重尺度溯源(即小尺度溯源)，基于第2重尺度溯源得到的潜在源区， 根据潜在源区内的企业排放数据、气象数据和地形数据，利用calpuff(非 稳态拉格朗日烟团模型系统)计算出各企业对报警点的污染浓度贡献(第k 个企业的贡献浓度)，并进行企业贡献度排名。通过小尺度溯源可以找出具体 的企业污染源头，并计算出各企业对报警点的浓度
贡献。
[0175]
多尺度融合修正分析，利用大尺度溯源计算出的各主要轨迹占比对中尺度 潜在源区贡献率进行修正：
[0176][0177]
其中，表示修正后的污染潜在源区贡献率，表示第l条主 要轨迹上的原始污染潜在源区贡献率，表示第l条主要轨迹的贡献比率。
[0178]
利用修正后的潜在源区贡献率对小尺度各企业浓度贡献进 行修正：
[0179][0180]
其中，表示修正后的第k个企业的贡献浓度值，表示网格(i,j) 区域内第k个企业的原始浓度贡献值。
[0181]
实施例14：
[0182]
本实施例的一种大气污染多尺度融合溯源方法(即多尺度融合修正的污 染精准溯源方法)，技术路线如图9所示，首先，得到监控站点报警时间、地 点、浓度值和气象数据300，根据上述数据，进行大尺度溯源310(即污染来 源轨迹溯源)，中尺度溯源320(即污染潜在区域溯源)和小尺度溯源330(即 污染来源企业溯源)，然后，进行多尺度融合340(即多尺度融合修正分析)。
[0183]
多尺度融合修正的污染精准溯源方法主要包括第一重大尺度方向轨迹溯 源、第二重中尺度潜在区域溯源、第三重小尺度企业溯源和多尺度融合修正分 析四大步骤，流程如图15所示，具体为：
[0184]
步骤s702：第一重大尺度方向轨迹溯源；
[0185]
当某监测站点监测到污染超标报警时，记录污染报警时间、地点和污染浓 度值。获取报警点区域以报警时间点为终点的前一段时间的气象数据，如获取 报警点区域近一周的气象数据。
[0186]
根据获取的气象数据，以监控站点位置为起点，利用hysplit进行后向 轨迹模拟，计算出所有污染来源轨迹。hysplit中假定质点的轨迹是随着风 场而运动的，轨迹是质点在空间和时间上的积分。质点所在位置的矢量速度在 时间和空间上都是线性插值得出的，其具体计算公式如下：
[0187]
p'(t+δt)＝p(t)+v(p,t)
·
δt
[0188]
p(t+δt)＝p(t)+0.5[v(p,t)+v(p',t+δt)]δt
[0189]
其中，p(t)表示t时刻质点位置，p'(t+δt)表示t+δt时刻质点位置，p(t+δt) 表示t+δt时刻质点位置通过矢量速度在时间和空间上线性插值得出，v(p,t)表 示质点在p位置和t时刻的速度，v(p',t+δt)表示质点在p
′
位置和t+δt时刻的速 度，δt表示时间差，t表示时间。p(
·
)表示质点的位置，v(
·
)表示质点的速度， 下时刻指定点的位置是由上一时刻的平均速度和猜值所在点的速度平均后与 时间步长的乘积而得到。本实施例以
污染报警时间前一周气象数据为例进行污 染来源轨迹模拟分析，模拟时间间隔为4小时，即每天计算6次，模拟高度分 别取100m、500m和1000m，因此一共可得到126条轨迹(126＝6次
×
3高度
×
7天)， 如图10所示。
[0190]
对计算出的126条污染来源轨迹，利用聚类算法进行轨迹聚类，得出主要 的4条污染来源及其贡献比率其中，代表某一种污染来源的概率，是 大尺度溯源结果，l＝1,2,3,4，图11给出了大尺度溯源结果示例。
[0191]
通过第一重大尺度溯源可以看出污染来源主要方向、输送轨迹及其贡献 率。
[0192]
步骤s704：第二重中尺度潜在区域溯源；
[0193]
首先，根据污染超标报警点确定污染溯源范围，举例而言，以报警点所在 城市为范围确定溯源区域为100公里
×
100公里。然后，将确定的溯源范围区 域进行网格划分，网格大小根据实际需要动态设置，本实施例以5公里
×
5公 里为例，如图12所示，其中黑点为报警点，大矩形区域为某市(其大小为100 公里
×
100公里)，每个小网格大小为5公里
×
5公里。每个5公里
×
5公里网 格便是本实施例定义的中尺度溯源区域。
[0194]
导入监控站点过去一周以小时为分辨率的污染监测浓度值，基于步骤 s702得到污染来源轨迹，计算轨迹经过网格(i,j)内超出某一设定阈值的网格 端点数m
ij
与网格总端点数n
ij
之比：
[0195][0196]
对于那些较偏远的网格，由于轨迹经过的端点数n
ij
太少，可能会给上式的 计算带来较大的不确定性。因此引入权重w
ij
来降低这些偏远网格所带来的影 响：
[0197][0198]
权重因子w
ij
的赋值主要取决于轨迹经过的某一网格端点数n
ij
与该网格平 均端点数n
ave
。因此，每个网格对报警点的潜在贡献率定义为：
[0199][0200]
其中，为网格(i,j)对报警点的污染贡献率，即潜在源区贡献率。
[0201]
进一步的，根据的值对每个网格进行颜色填充和可视化展示，网格的 贡献率值越大，其颜色越深，越有可能是潜在源区，如图13所示，图中右侧 标记的方形表示贡献率范围，从上至下，逐渐增加，从a至e可以对应不同 的颜色，颜色越深其潜在源区贡献率越大。
[0202]
优选的，进一步计算每个网格(i,j)的污染物对报警点相对权重浓度值，弥 补潜在源区贡献率分析无法区分相同贡献值的网格对报警点浓度贡献大小的 缺陷，其计算原理如下：
[0203][0204]
其中，c
ij
是网格(i,j)的平均权重浓度，c
l
是轨迹l经过网格(i,j)时的污染 物浓度，t
ijl
是轨迹l经过网格(i,j)的停留时间。采用与潜在源区贡献率分析相 同的权重公式和方法，可获得
[0205]
步骤s706：第三重小尺度企业溯源；
[0206]
针对每一个潜在源区(网格点(i,j)),获取网格点内企业排放信息，包括但 不限于企业经纬度、烟囱源高、烟囱出口直径、烟囱出口流速、烟囱出口烟气 温度等。获取网格点到报警点之间区域的地形数据和气象数据，其中气象数据 包括地面气象数据和高空模拟数据。将上面得到的企业排放信息、地形数据和 气象数据输入calpuff模型计算出网格内各企业排放污染到达报警点处时 的浓度值(k表示第k个企业)，其中，便是小尺度企业点源溯源结果。
[0207]
进一步的，根据每个企业对报警点的浓度值，计算区域内企业的贡献浓度 比例，并进行排名可视化展示，如图14所示。
[0208]
步骤s708：多尺度融合修正分析；
[0209]
在进行第二重中尺度溯源时，针对各个潜在源区(网格(i,j))：首先，确 定网格(i,j)落在哪个主要轨迹l上，其中，轨迹l由第一重大尺度溯源得到；然 后，利用大尺度溯源计算出的各主要轨迹占比对中尺度潜在源区贡献率进行修 正：
[0210][0211]
其中，为修正后的潜在源区贡献率，为第l条主要轨迹上的 原始潜在源区贡献率，为第l条主要轨迹的贡献比率。
[0212]
在进行第三重小尺度溯源时，针对各个企业：首先，确定企业k落在哪个 网格(i,j)内；然后，利用修正后的潜在源区贡献率对小尺度各企业贡献 浓度进行修正：
[0213][0214]
其中，为修正后第k个企业的贡献浓度值，为第(i,j)区域内第 个企业的原始浓度贡献值，为修正后的第(i,j)源区贡献率。
[0215]
目前现有的大气污染多尺度融合溯源方法要么只能进行大方向大区域溯 源，要么只能进行小尺度企业关联分析溯源。本实施例提出一种多尺度逐层递 进的溯源方法，从大方向到小区域再到企业点源逐层精细化溯源，满足不同环 境治理决策需要。
[0216]
综上，本发明实施例的有益效果为：
[0217]
1.本实施例充分利用不同尺度溯源结果进行融合分析，使用大尺度溯 源结果对下一级尺度溯源结果进行修正，可以弥补由于单个模型计算误差、 奇异点和数据不完整等带来的计算结果不准确问题，提升潜在源区贡献率 企业污染浓度贡献度的准确性。
[0218]
2.本实施例中，通过第一重大尺度溯源可以看出污染来源主要方向、输 送轨迹及主要轨迹贡献比率，可以对污染来源进行大方向上的溯源，通过对 污染溯源范围进行网格划分，完成中尺度潜在区域溯源，通过计算企业对 监控站点的污染浓度贡献值，完成小尺度溯源，可以满足不同环境治理决策 需要。
[0219]
在本发明中，术语“第一”、“第二”、“第三”仅用于描述的目的，而不能理 解为指示或暗示相对重要性；术语“多个”则指两个或两个以上，除非另有明确 的限定。术语“安装”、“相连”、“连接”、“固定”等术语均应做广义理解，例如，
ꢀ“
连接”可以是固定连接，也可以是可拆卸连接，或一体地连接；“相连”可以是 直接相连，也可以通过中间媒介间接相连。对于本领域的普通技术人员而言， 可以根据具体情况理解上述术语在本发明中的具体含义。
[0220]
本发明的描述中，需要理解的是，术语“上”、“下”、“左”、“右”、“前”、“后
”ꢀ
等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描 述本发明和简化描述，而不是指示或暗示所指的装置或单元必须具有特定的方 向、以特定的方位构造和操作，因此，不能理解为对本发明的限制。
[0221]
在本说明书的描述中，术语“一个实施例”、“一些实施例”、“具体实施例
”ꢀ
等的描述意指结合该实施例或示例描述的具体特征、结构、材料或特点包含于 本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述 不一定指的是相同的实施例或实例。而且，描述的具体特征、结构、材料或特 点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0222]
以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的 技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内， 所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。