一种大气污染物溯源方法及系统与流程

1.本技术涉及空气质量监测技术领域，尤其涉及一种大气污染物溯源方法及系统。


背景技术：

2.空气环境质量和人们的健康息息相关,及时发现影响空气质量的污染源有助于环保部门进行及时决策，防止环境质量持续恶化。目前,实现大气污染溯源的方式主要分为两种，一种是基于站点监测的方式，在区域内布设精细化的监控网格，实时监测区域内的环质量变化趋势，当发生污染时，可根据数据的变化趋势路径进行污染溯源。另一种是基于空气质量模型，在模拟区域气象的基础上，进行污染物扩散以及各个污染物之间的生化反应的计算对空气质量进行模拟,当某个站点的空气质量出现异常时，可通过模型网格的输出结果和传输过程定位到污染源，从而进行污染溯源。
3.基于站点监测的方法依赖于园区及其周边建立精密的监测网格，网格越精密，溯源效果越好，但同时运维网格站点需要大量专业技术人员，维护成本较高难以实施。
4.基于空气质量模型进行溯源的方式，其局限性主要是区域三维精细化网格气象场的模拟需要占用较多的计算资源，导致空气质量模型的计算速度受到制约，往往不能满足污染溯源的时效性。而且模型适用于较大尺度区域，对于较小区域内的溯源效果精准性较差，而且模型移植性难度较大，需要重新率定模型，对于数据缺失的区域模型难以调优。


技术实现要素：

5.本技术的目的在于提供一种大气污染物溯源方法及系统，该方法基于污染指标浓度数据进行污染物的溯源，溯源效果越好、计算速度较快、对于较小区域内的溯源精准性高。
6.为达到上述目的，本技术提供一种大气污染物溯源方法，该方法包括：对监测站点的污染指标浓度数据进行趋势异常判断，获取数据趋势异常的监测站点作为污染站点；采集污染站点、监测站点及企业的污染指标浓度数据；根据污染站点的污染指标浓度数据，获取污染站点的污染起始时刻；根据监测站点及企业的污染指标浓度数据，获取污染站点的污染起始时刻前一段时间内，监测站点及企业的污染指标浓度数据满足高值点的监测站点集合和企业集合；将污染站点、监测站点集合中的监测站点和企业集合中的企业构建成传输轨迹，形成传输轨迹集合；根据传输轨迹集合，获取污染站点的污染源。
7.如上的，其中，获取污染站点的污染起始时刻的方法为：获取污染站点a的连续时间段内的污染指标浓度序列；对污染指标浓度序列进行检验，若该污染指标浓度序列呈上升趋势，并且污染指标浓度序列的趋势统计量值大于或等于预设阈值，则定义连续时间段的初始时刻为污染站点的污染起始时刻。
8.如上的，其中，获取监测站点集合的方法包括：获取从t0开始到之前的t时刻内，与污染站点距离为d范围内的所有监测站点的第一污染指标浓度数据；根据第一污染指标浓度数据，获取所有监测站点的第一污染指标浓度序列；根据第一污染指标浓度序列，获取满
足连续时刻数为m的第一污染指标浓度序列趋势统计量值大于或等于预设阈值的监测站点，形成监测站点集合。
9.如上的，其中，获取企业集合的方法包括：获取从t0开始到之前的t时刻内，与污染站点距离为d范围内的所有企业的第二污染指标浓度数据；根据第二污染指标浓度数据，获取所有企业的第二污染指标浓度序列；根据第二污染指标浓度序列，获取满足连续时刻数为m的第二污染指标浓度序列趋势统计量值大于或等于预设阈值的监测站点，形成企业集合。
10.如上的，其中，形成传输轨迹集合的方法包括：以污染站点为起始点、企业为结束点，监测站点为中间连接点，构建完整轨迹，形成轨迹集合；将轨迹集合中的所有轨迹转换到同一直角坐标系中，形成基础传输轨迹集合；计算基础传输轨迹集合中所有轨迹的稳定概率；根据轨迹的稳定概率，消除基础传输轨迹集合中满足预设第三条件的轨迹，形成最终的传输轨迹集合。
11.如上的，其中，形成轨迹集合的方法包括：污染站点的经纬度和企业集合中的每一个企业的经纬度均构成一条完整轨迹；污染站点的经纬度和监测站点集合中的每一个监测站点的经纬度均构成一条缺失轨迹；将满足预设第一条件的站点和满足预设第二条件的企业添加到缺失轨迹中，构建成完整轨迹。
12.如上的，其中，形成基础传输轨迹集合的方法包括：以污染站点的经纬度为坐标系原点，构建直角坐标系；将监测站点集合中的监测站点和企业集合中的企业的经纬度坐标转换为直角坐标系下的坐标点；在直角坐标系中，对于轨迹集合中的点利用贝塞尔曲线形成基础传输轨迹。
13.如上的，其中，计算基础传输轨迹集合中所有轨迹的稳定概率的方法包括：计算轨迹集合的轨迹中除污染站点外的监测站点和企业的传输概率；根据监测站点和企业的传输概率，计算轨迹集合中所有轨迹的稳定概率。
14.如上的，其中，根据传输轨迹集合，获取污染站点的污染源的方法为：对于传输轨迹集合k中的所有轨迹，以轨迹结束点的企业为起始点，以污染站点a为结束点，中间的监测站点为控制点，绘制贝塞尔曲线，得到的曲线即为传输路径；其中，传输路径起始点的企业即为造成污染站点a的污染源。
15.一种大气污染物溯源系统，该系统包括：污染站点获取模块，用于对监测站点的污染指标浓度数据进行趋势异常判断，获取数据趋势异常的监测站点作为污染站点；采集装置，用于采集污染站点、监测站点及企业的污染指标浓度数据；第一获取模块，用于根据污染站点的污染指标浓度数据，获取污染站点的污染起始时刻；第二获取模块，用于根据监测站点及企业的污染指标浓度数据，获取污染站点的污染起始时刻前一段时间内，监测站点及企业的污染指标浓度数据满足高值点的监测站点集合和企业集合；传输轨迹形成模块，用于将污染站点、监测站点集合中的监测站点和企业集合中的企业构建成传输轨迹，形成传输轨迹集合；污染源获取模块，用于根据传输轨迹集合，获取污染站点的污染源。
16.本技术实现的有益效果如下：
17.(1)本技术污染溯源方法简单可行，且容易推广，适用于尺度较小的区域，溯源的准确度较高，溯源效果越好，提高了溯源效率。
18.(2)本技术在找到污染源的同时，也绘制出了污染源的传输路径，为污染的传输过
程提供数据支持。
附图说明
19.为了更清楚地说明本技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术中记载的一些实施例，对于本领域技术人员来讲，还可以根据这些附图获得其他的附图。
20.图1为本技术实施例的一种大气污染物溯源方法的流程图。
21.图2为本技术实施例的获取污染站点的污染起始时刻的方法流程图。
22.图3为本技术实施例的获取监测站点集合和企业集合的方法流程图。
23.图4为本技术实施例的形成传输轨迹集合的方法流程图。
24.图5为本技术实施例的形成轨迹集合的方法流程图。
25.图6为本技术实施例的形成基础传输轨迹集合的方法流程图。
26.图7为本技术实施例的计算基础传输轨迹集合中所有轨迹的稳定概率的方法流程图。
27.图8为本技术实施例的一种大气污染物溯源系统的结构示意图。
28.附图标记：10
‑
污染站点获取模块；20
‑
采集装置；30
‑
第一获取模块；40
‑
第二获取模块；50
‑
传输轨迹形成模块；60
‑
污染源获取模块；100
‑
大气污染物溯源系统。
具体实施方式
29.下面结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
30.实施例一
31.如图1所示，本技术提供一种大气污染物溯源方法，该方法包括如下步骤：
32.步骤s0，对监测站点的污染指标浓度数据进行趋势异常判断，获取数据趋势异常的监测站点作为污染站点。
33.具体的，采集监测站点的污染指标浓度数据，并利用mann—kenddall(曼
‑
肯德尔法)检验法，对监测站点的污染指标浓度数据进行检验，若监测站点的污染指标浓度数据呈上升趋势，并且污染指标浓度数据的趋势z值(趋势的统计量)不小于(大于或等于)预设阈值zn，则定义该监测站点为污染站点。
34.步骤s1，采集污染站点、监测站点及企业的污染指标浓度数据。
35.具体的，采集污染站点、与污染站点距离为d范围内的所有监测站点及与污染站点距离为d范围内的所有企业的污染指标浓度数据。
36.步骤s2，根据污染站点的污染指标浓度数据，获取污染站点的污染起始时刻。
37.如图2所示，步骤s2包括如下子步骤：
38.步骤s210，获取污染站点a的连续时间段内的污染指标浓度序列。
39.具体的，获取污染站点a的前n时刻内的污染指标浓度序列为：c＝[c0,c1,...,c
n
]。
[0040]
步骤s220，利用mann—kenddall(曼
‑
肯德尔法)检验法，对污染指标浓度序列进行检验，若该污染指标浓度序列呈上升趋势，并且污染指标浓度序列的趋势z值(趋势的统计量)不小于(大于或等于)预设阈值zn，则定义连续时间段的初始时刻t0为污染站点的污染起始时刻。
[0041]
步骤s3，根据监测站点及企业的污染指标浓度数据，获取污染站点的污染起始时刻前一段时间内，监测站点及企业的污染指标浓度数据满足高值点的监测站点集合和企业集合。
[0042]
具体的，从t0开始到之前的t时刻内，计算污染指标p的连续时刻数为m的浓度序列趋势z值不小于zn的所有满足条件的监测站点集合s和企业集合f，并记录每个监测站点和企业的污染指标浓度序列的污染起始时刻。
[0043]
如图3所示，步骤s3包括如下子步骤：
[0044]
步骤s310，获取从t0开始到之前的t时刻内，与污染站点距离为d范围内的所有监测站点的污染指标浓度序列，及与污染站点距离为d范围内的所有企业的污染指标浓度序列。
[0045]
步骤s320，根据监测站点的污染指标浓度序列，利用mann—kenddall(曼
‑
肯德尔法)检验法，获取满足连续时刻数为m的监测站点的污染指标浓度序列趋势z值不小于预设阈值zn的监测站点，形成监测站点集合s。
[0046]
步骤s330，根据企业的污染指标浓度序列，利用mann—kenddall(曼
‑
肯德尔法)检验法，获取满足连续时刻数为m的企业的污染指标浓度序列趋势z值不小于zn的企业，形成企业集合f。
[0047]
步骤s4,将污染站点、监测站点集合中的监测站点和企业集合中的企业构建成传输轨迹，形成传输轨迹集合。
[0048]
如图4所示，步骤s4包括如下子步骤：
[0049]
步骤s410，以污染站点为起始点、企业为结束点，监测站点为中间连接点，构建完整轨迹，形成轨迹集合。
[0050]
如图5所示，步骤s410包括如下子步骤：
[0051]
步骤s411，污染站点a的经纬度和企业集合f中的每一个企业的经纬度均构成一条完整轨迹，类型为h。
[0052]
具体的，以污染站点a的经纬度为起始点，以企业集合f中的每一个企业的经纬度为结束点，构成多条完整轨迹，类型为h。
[0053]
步骤s412，污染站点a的经纬度和监测站点集合s中的每一个监测站点的经纬度均构成一条缺失轨迹。
[0054]
具体的，污染站点a的经纬度监测站点集合s中的所有监测站点的经纬度形成多条缺失轨迹。
[0055]
步骤s413，将满足预设第一条件的站点和满足预设第二条件的企业添加到缺失轨迹中，构建成完整轨迹。
[0056]
步骤s413包括如下子步骤：
[0057]
步骤s4131，将满足预设第一条件的监测站点添加到缺失轨迹中。
[0058]
具体的，预设第一条件为:对于监测站点集合s中的监测站点g，只要监测站点集合
s中有污染起始时刻小于监测站点g的监测站点f，并且该监测站点f和监测站点g在污染站点a的同向半平面内，将该监测站点f添加到监测站点g的缺失轨迹中。
[0059]
步骤s4132，若缺失轨迹无可加入的监测站点，则将满足预设第二条件的企业添加到缺失轨迹的结束点，形成一条完整轨迹，类型为w。
[0060]
作为本发明的实施例，第二条件包括：第一：企业的污染起始时刻小于最后添加到缺失轨迹中监测站点的污染起始时刻；第二：企业和缺失轨迹的最后一个监测站点g在站点a的同向半平面内。
[0061]
步骤s420，将轨迹集合中的所有轨迹转换到同一直角坐标系中，形成基础传输轨迹集合。
[0062]
如图6所示，步骤s420包括如下子步骤：
[0063]
步骤s421，以污染站点a的经纬度为坐标系原点，构建直角坐标系。
[0064]
步骤s422，将监测站点集合s中的监测站点和企业集合f中的企业的经纬度坐标转换为直角坐标系下的坐标点。
[0065]
步骤s423，在直角坐标系中，对于轨迹集合中的点利用贝塞尔曲线形成基础传输轨迹。
[0066]
步骤s430，计算基础传输轨迹集合中所有轨迹的稳定概率。
[0067]
如图7所示，步骤s430包括如下子步骤：
[0068]
步骤s431，计算轨迹集合的轨迹中除污染站点a外的点(包括监测站点和企业)的传输概率。
[0069]
具体的，点的传输概率计算公式如下：
[0070][0071][0072]
其中，p(γ
ir
)表示轨迹γ
ir
的第r个点的传输概率；w
v
表示点污染起始时刻的风向；α表示轨迹曲线中点的切线与正北方向之间的夹角；0
°
≤α＜360
°
；该切线位于点的传输方向一侧；w
s
表示风速；ε为一个趋近于0的正数，sgn()为符号函数；abs()表示求绝对值运算；σ
p
为参数。
[0073]
步骤s432，根据点的传输概率，计算轨迹集合中所有轨迹的稳定概率。
[0074]
具体的，轨迹的稳定概率的计算公式为：
[0075][0076]
其中，γ
i
表示第i条完整轨迹；k表示轨迹γ
i
中除污染站点a外所有点(包括监测站点和企业)的总个数；m表示所有完整轨迹的条数；p(γ
ir
)表示轨迹γ
ir
的第r个点(包括监测站点和企业)的传输概率；z
r
为第i条完整轨迹γ
i
中第r个点的趋势检验z值。
[0077]
步骤s440，根据轨迹的稳定概率，消除基础传输轨迹集合中满足预设第三条件的轨迹，形成最终的传输轨迹集合k。
[0078]
具体的，预设第三条件为：对于轨迹h，如果存在一条比轨迹h只多一个监测站点的轨迹w，并且p(γ
w
)≥p(γ
h
)，则删除轨迹h，直到轨迹集合中所有的轨迹点不同。
[0079]
步骤s5，根据传输轨迹集合，获取污染站点的污染源。
[0080]
具体的，对于传输轨迹集合k中的所有轨迹，以轨迹结束点的企业为起始点，以污染站点a为结束点，中间的监测站点为控制点，绘制贝塞尔曲线，得到的曲线即为传输路径；其中，传输路径起始点的企业即为造成污染站点a的污染指标呈上升趋势的污染源。
[0081]
实施例二
[0082]
如图8所示，本技术提供一种大气污染物溯源系统100，该系统包括：
[0083]
污染站点获取模块10，用于对监测站点的污染指标浓度数据进行趋势异常判断，获取数据趋势异常的监测站点作为污染站点；
[0084]
采集装置20，用于采集污染站点、监测站点及企业的污染指标浓度数据；
[0085]
第一获取模块30，用于根据污染站点的污染指标浓度数据，获取污染站点的污染起始时刻；
[0086]
第二获取模块40，用于根据监测站点及企业的污染指标浓度数据，获取污染站点的污染起始时刻前一段时间内，监测站点及企业的污染指标浓度数据满足高值点的监测站点集合和企业集合；
[0087]
传输轨迹形成模块50，用于将污染站点、监测站点集合中的监测站点和企业集合中的企业构建成传输轨迹，形成传输轨迹集合；
[0088]
污染源获取模块60，用于根据传输轨迹集合，获取污染站点的污染源。
[0089]
本技术实现的有益效果如下：
[0090]
(1)本技术污染溯源方法简单可行，且容易推广，适用于尺度较小的区域，溯源的准确度较高，溯源效果越好，提高了溯源效率。
[0091]
(2)本技术在找到污染源的同时，也绘制出了污染源的传输路径，为污染的传输过程提供数据支持。
[0092]
上所述仅为本发明的实施方式而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理的内所做的任何修改、等同替换、改进等，均应包括在本发明的权利要求范围之内。