基于统计自适应方法计算突发污染源和外来污染输送贡献技术

本发明专利技术公开了一种基于统计自适应方法计算突发污染源和外来污染输送贡献，包括：a.在一个空气质量检测区域内设置若干个检测站点，获得每个站点的污染物浓度—小时或分钟数据序列，并获得污染物浓度增长曲线；b.检测突然升高数据，并进行判识；c.若为单个站点突变则判定为突发污染源，若为所有站点同时或依次突变则判定为外来污染输送；d.突发污染源贡献是由突变期间浓度减去突变前后时刻浓度线性插值结果计算得出；单个站点的外来污染输送贡献是由外来污染输送过程的增长曲线减去外来污染传输前的增长曲线计算得出，整个区域的贡献是由所有站点的外来污染输送贡献之和然后算术平均得出。本发明专利技术计算结果准确，且不需要污染源信息。

Calculation of Sudden Pollution Sources and Foreign Pollution Transport Contributions Based on Statistical Adaptive Method

The invention discloses a statistical adaptive method for calculating the contribution of sudden pollution sources and external pollution transportation, which includes: A. setting up several detection stations in an air quality detection area, obtaining the pollutant concentration-hour or minute data sequence of each station, and obtaining the pollutant concentration growth curve; B. detecting the sudden increase data, and identifying if it is a single station; Station mutation is judged as sudden pollution source, and if all stations mutate simultaneously or sequentially, it is judged as external pollution transport; D. The contribution of sudden pollution source is calculated by linear interpolation of concentration before and after mutation subtracted from concentration during mutation; the contribution of external pollution transport of a single station is calculated by the growth curve of external pollution transport process subtracted from the growth curve before transmission of external pollution. Line calculation shows that the contribution of the whole region is the sum of the contribution of external pollution transport from all stations and then arithmetically averaged. The calculation result of the invention is accurate and does not need the information of the pollution source.

【技术实现步骤摘要】
基于统计自适应方法计算突发污染源和外来污染输送贡献
本专利技术涉及一种空气质量监测中的突发污染源贡献和外来污染输送贡献计算方法，特别是涉及一种基于统计的自适应方法计算突发污染源和外来污染输送贡献，属于环境监测

技术介绍
突发污染源和外来污染输送贡献，传统的计算方法主要采用数值模型模拟和污染物来源解析，但该方法存在以下两个方面的缺点：一方面，计算突发污染源和外来污染输送贡献，均需要污染源信息；另一方面，由于污染源谱的不确定性和模式物理化学过程与实际大气过程有差异，导致计算时间长、成本和技术要求较高，计算结果也存在较大不确定性。
技术实现思路
本专利技术要解决的技术问题是提供一种计算方法简单，成本低，可靠性高，而且不需要污染源信息的基于统计自适应方法计算突发污染源和外来污染输送贡献。为解决上述技术问题，本专利技术采用这样一种基于统计自适应方法计算突发污染源和外来污染输送贡献，包括以下步骤：a.在一个空气质量检测区域内设置若干个间隔分散分布的空气污染检测站点，通过设置在每个空气污染检测站点内的空气污染检测仪器，获得每个空气污染检测站点的污染物浓度—小时或分钟数据序列，并由所述污染物浓度—小时或分钟数据序列，获得每个空气污染检测站点的污染物浓度增长曲线；b.对每个空气污染检测站点内的所述污染物浓度—小时或分钟数据序列变化趋势进行突变检测，并相互比对，根据所述污染物浓度—小时或分钟数据序列变化趋势，检测突然升高的数据，并对该突然升高的数据段进行判识；c.通过对所述污染物浓度—小时或分钟数据序列变化趋势的突变检测，判断所述突然升高的数据段是单个空气污染检测站点突变，还是所有空气污染检测站点同时或依次突变，若为单个站点突变则将所述突然升高的数据段判定为突发污染源影响，若为所有站点同时或依次突变则将所述突然升高的数据段判定为外来污染输送影响；d.在步骤c中，若判定为突发污染源影响，则突发污染源贡献是由突变期间污染物浓度减去突变前后时刻污染物浓度线性插值结果计算得出；若判定为外来污染输送影响，则单个空气污染检测站点的外来污染输送贡献是由单个空气污染检测站点外来污染输送过程的污染物浓度增长曲线减去外来污染传输前的污染物浓度增长曲线计算得出；整个空气质量检测区域的外来污染输送贡献是由同一区域内所有空气污染检测站点的外来污染输送贡献之和然后算术平均得出。在本专利技术中，在步骤a中，由所述污染物浓度—小时或分钟数据序列减去参考浓度得到每个空气污染检测站点的污染物浓度增长曲线，所述污染物浓度—小时或分钟数据序列初始时段往前推24小时内的最低浓度或者前2日的平均最低浓度为所述参考浓度；在步骤d中，所述外来污染传输前的污染物浓度增长曲线为选取外来污染传输前1～2日的污染物浓度增长曲线，或者外来污染传输过程前一段时间的污染物浓度增长曲线。采用上述计算方法后，本专利技术具有以下有益效果：本专利技术采用统计方法，根据污染物浓度—小时或分钟数据序列及其波动增长曲线特征，对空气质量监测数据直接进行统计分析，计算出突发污染源贡献和外来污染输送贡献，不需要污染源信息，计算方法简单，成本低，可靠性高。本专利技术大幅缩短了计算时间，并且计算结果准确，同时对操作人员的技术要求较低。附图说明以下结合附图对本专利技术的具体实施方式作进一步的详细说明。图1为本专利技术中突发污染源贡献的一种计算示意图。图2为本专利技术中外来污染输送贡献的一种计算示意图。具体实施方式参见图1和图2，本专利技术提供了一种基于统计自适应方法计算突发污染源和外来污染输送贡献，包括以下步骤：a.在一个空气质量检测区域内设置若干个间隔分散分布的空气污染检测站点，所述的一个空气质量检测区域可以是某个城市，也可以是某个城市的一个行政区域等，所述的间隔分散分布可以是相互间隔数公里，例如5～10公里等，相互间隔的距离可根据城市大小实际情况确定，通过设置在每个空气污染检测站点内的空气污染检测仪器，所述仪器可采用市售的ThermoFisher5030、MetOne1020等各类型颗粒物监测设备，包含β射线原理、震荡天平原理等各类型监测设备，获得每个空气污染检测站点的污染物浓度—小时或分钟数据序列，并由所述污染物浓度—小时或分钟数据序列，获得每个空气污染检测站点的污染物浓度增长曲线；在本专利技术中，所述污染物浓度—小时或分钟数据序列可以是每隔1小时或15分钟或30分钟等时间段得到的污染物浓度数据序列，本专利技术优选为每隔1小时；b.对每个空气污染检测站点内的所述污染物浓度—小时或分钟数据序列变化趋势进行突变检测，并相互比对，根据所述污染物浓度—小时或分钟数据序列变化趋势，检测突然升高的数据，如图1中的A1段，并对该突然升高的数据段进行判识；c.通过对所述污染物浓度—小时或分钟数据序列变化趋势的突变检测，判断所述突然升高的数据段是单个空气污染检测站点突变，还是所有空气污染检测站点同时或依次突变，若为单个站点突变则将所述突然升高的数据段判定为突发污染源影响，若为所有站点同时或依次突变则将所述突然升高的数据段判定为外来污染输送影响；d.在步骤c中，若判定为突发污染源影响，则突发污染源贡献是由突变期间污染物浓度减去突变前后时刻污染物浓度线性插值结果计算得出；若判定为外来污染输送影响，则单个空气污染检测站点的外来污染输送贡献是由单个空气污染检测站点外来污染输送过程的污染物浓度增长曲线减去外来污染传输前的污染物浓度增长曲线计算得出；整个空气质量检测区域的外来污染输送贡献是由同一区域内所有空气污染检测站点的外来污染输送贡献之和然后算术平均得出。在本专利技术中，在步骤a中，由所述污染物浓度—小时或分钟数据序列减去参考浓度得到每个空气污染检测站点的污染物浓度增长曲线，所述污染物浓度—小时或分钟数据序列初始时段往前推24小时内的最低浓度或者前2日的平均最低浓度为所述参考浓度；在步骤d中，所述外来污染传输前的污染物浓度增长曲线优选为选取外来污染传输前1～2日的污染物浓度增长曲线，或者外来污染传输过程前一段时间的污染物浓度增长曲线，所述前一段时间可以是前数小时，例如前6小时或12小时或24小时等。参见图1所示的突发污染源贡献一种优选实施例计算示意图，首先，根据空气污染检测仪器监测得到的每隔1小时的污染物浓度数据序列，例如图1中所示的1点的浓度为6μg/m3，2点的浓度为5μg/m3，3点的浓度为7.2μg/m3，4点的浓度为7.9μg/m3，5点的浓度为8.8μg/m3，6点的浓度为10.2μg/m3，7点的浓度为13.5μg/m3，8点的浓度为19μg/m3，9点的浓度为55μg/m3，10点的浓度为24.5μg/m3，11点的浓度为19μg/m3，12点的浓度为10.5μg/m3，13点的浓度为8μg/m3，14点的浓度为11μg/m3，15点的浓度为17.5μg/m3，16点的浓度为25μg/m3，17点的浓度为33.5μg/m3，18点的浓度为36μg/m3，19点的浓度为34μg/m3，20点的浓度为29.5μg/m3，21点的浓度为19.5μg/m3，22点的浓度为13μg/m3，23点的浓度为10.5μg/m3，24点的浓度为9μg/m3，由上述污染物浓度—小时数据序列得到污染物浓度曲线A，在图中，8点～10点之间的数据段是突然升高的数据段A1，本文档来自技高网...

【技术保护点】
1.一种基于统计自适应方法计算突发污染源和外来污染输送贡献，其特征在于，包括以下步骤：a. 在一个空气质量检测区域内设置若干个间隔分散分布的空气污染检测站点，通过设置在每个空气污染检测站点内的空气污染检测仪器，获得每个空气污染检测站点的污染物浓度—小时或分钟数据序列，并由所述污染物浓度—小时或分钟数据序列，获得每个空气污染检测站点的污染物浓度增长曲线；b. 对每个空气污染检测站点内的所述污染物浓度—小时或分钟数据序列变化趋势进行突变检测，并相互比对，根据所述污染物浓度—小时或分钟数据序列变化趋势，检测突然升高的数据，并对该突然升高的数据段进行判识；c. 通过对所述污染物浓度—小时或分钟数据序列变化趋势的突变检测，判断所述突然升高的数据段是单个空气污染检测站点突变，还是所有空气污染检测站点同时或依次突变，若为单个站点突变则将所述突然升高的数据段判定为突发污染源影响，若为所有站点同时或依次突变则将所述突然升高的数据段判定为外来污染输送影响；d. 在步骤c中，若判定为突发污染源影响，则突发污染源贡献是由突变期间污染物浓度减去突变前后时刻污染物浓度线性插值结果计算得出；若判定为外来污染输送影响，则单个空气污染检测站点的外来污染输送贡献是由单个空气污染检测站点外来污染输送过程的污染物浓度增长曲线减去外来污染传输前的污染物浓度增长曲线计算得出；整个空气质量检测区域的外来污染输送贡献是由同一区域内所有空气污染检测站点的外来污染输送贡献之和然后算术平均得出。...

【技术特征摘要】
1.一种基于统计自适应方法计算突发污染源和外来污染输送贡献，其特征在于，包括以下步骤：a.在一个空气质量检测区域内设置若干个间隔分散分布的空气污染检测站点，通过设置在每个空气污染检测站点内的空气污染检测仪器，获得每个空气污染检测站点的污染物浓度—小时或分钟数据序列，并由所述污染物浓度—小时或分钟数据序列，获得每个空气污染检测站点的污染物浓度增长曲线；b.对每个空气污染检测站点内的所述污染物浓度—小时或分钟数据序列变化趋势进行突变检测，并相互比对，根据所述污染物浓度—小时或分钟数据序列变化趋势，检测突然升高的数据，并对该突然升高的数据段进行判识；c.通过对所述污染物浓度—小时或分钟数据序列变化趋势的突变检测，判断所述突然升高的数据段是单个空气污染检测站点突变，还是所有空气污染检测站点同时或依次突变，若为单个站点突变则将所述突然升高的数据段判定为突发污染源影响，若为所有站点同时或依次突变则将所述突然升高的数据段判定为外来污染输送影响；d.在步骤c...

【专利技术属性】
技术研发人员：何涛，李艳萍，夏京，彭燕，杨卫芬，
申请(专利权)人：常州市环境监测中心，
类型：发明
国别省市：江苏,32