基于统计自适应方法计算突发污染源和外来污染输送贡献与流程
本发明涉及一种空气质量监测中的突发污染源贡献和外来污染输送贡献计算方法，特别是涉及一种基于统计的自适应方法计算突发污染源和外来污染输送贡献，属于环境监测技术领域。
背景技术：
突发污染源和外来污染输送贡献，传统的计算方法主要采用数值模型模拟和污染物来源解析，但该方法存在以下两个方面的缺点：一方面，计算突发污染源和外来污染输送贡献，均需要污染源信息；另一方面，由于污染源谱的不确定性和模式物理化学过程与实际大气过程有差异，导致计算时间长、成本和技术要求较高，计算结果也存在较大不确定性。
技术实现要素：
本发明要解决的技术问题是提供一种计算方法简单，成本低，可靠性高，而且不需要污染源信息的基于统计自适应方法计算突发污染源和外来污染输送贡献。
为解决上述技术问题，本发明采用这样一种基于统计自适应方法计算突发污染源和外来污染输送贡献，包括以下步骤：
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，该突然升高的数据段a1仅为单个站点突变，其中，9点的浓度55μg/m3为突变点也即突变期间污染物浓度，突变前时刻污染物浓度即8点的浓度19μg/m3，突变后时刻污染物浓度即10点的浓度24.5μg/m3，突变前后时刻污染物浓度线性插值即8点与10点之间的平均浓度21.8μg/m3，则突发污染源贡献由突变期间污染物浓度55μg/m3减去突变前后时刻污染物浓度线性插值21.8μg/m3计算得出33.2μg/m3。
参见图2所示的外来污染输送贡献一种优选实施例计算示意图，首先，通过设置在每个空气污染检测站点内的空气污染检测仪器，监测得到每个空气污染检测站点外来污染输送过程的每隔1小时的污染物浓度数据序列，例如图2中所示单个空气污染检测站点外来污染输送过程的1点的浓度为68μg/m3，2点的浓度为78μg/m3，3点的浓度为87μg/m3，4点的浓度为100μg/m3，5点的浓度为108μg/m3，6点的浓度为118μg/m3，7点的浓度为132μg/m3，8点的浓度为139μg/m3，9点的浓度为147μg/m3，10点的浓度为155μg/m3，11点的浓度为158μg/m3，12点的浓度为169μg/m3，13点的浓度为174μg/m3，14点的浓度为172μg/m3，15点的浓度为169μg/m3，16点的浓度为167μg/m3，17点的浓度为146μg/m3，18点的浓度为161μg/m3，19点的浓度为185μg/m3，20点的浓度为182μg/m3，21点的浓度为168μg/m3，22点的浓度为158μg/m3，23点的浓度为147μg/m3，24点的浓度为124μg/m3，选取该数据序列初始时段往前推24小时内的最低浓度例如64μg/m3为参考浓度，通过上述污染物浓度—小时数据序列依次减去参考浓度得到1点～24点的浓度依次为4μg/m3、14μg/m3、23μg/m3、36μg/m3、44μg/m3、54μg/m3、68μg/m3、75μg/m3、83μg/m3、91μg/m3、94μg/m3、105μg/m3、110μg/m3、108μg/m3、105μg/m3、103μg/m3、82μg/m3、97μg/m3、121μg/m3、118μg/m3、104μg/m3、94μg/m3、83μg/m3和60μg/m3，从而得到单个空气污染检测站点外来污染输送过程的污染物浓度增长曲线b；通过设置在每个空气污染检测站点内的空气污染检测仪器，监测得到每个空气污染检测站点外来污染传输前1日的每隔1小时的污染物浓度数据序列，例如图2中所示上述单个空气污染检测站点外来污染传输前1日的1点的浓度为48.5μg/m3，2点的浓度为47.5μg/m3，3点的浓度为49.5μg/m3，4点的浓度为50.5μg/m3，5点的浓度为53μg/m3，6点的浓度为56.5μg/m3，7点的浓度为61.5μg/m3，8点的浓度为67μg/m3，9点的浓度为72.5μg/m3，10点的浓度为72.5μg/m3，11点的浓度为67μg/m3，12点的浓度为58.5μg/m3，13点的浓度为56μg/m3，14点的浓度为59μg/m3，15点的浓度为65.5μg/m3，16点的浓度为73μg/m3，17点的浓度为81.5μg/m3，18点的浓度为84μg/m3，19点的浓度为82μg/m3，20点的浓度为77.5μg/m3，21点的浓度为67.5μg/m3，22点的浓度为61μg/m3，23点的浓度为58.5μg/m3，24点的浓度为64μg/m3，选取该数据序列初始时段往前推24小时内的最低浓度例如48μg/m3为参考浓度，通过上述污染物浓度—小时数据序列依次减去参考浓度得到1点～24点的浓度依次为0.5μg/m3、-0.5μg/m3、1.5μg/m3、2.5μg/m3、5μg/m3、8.5μg/m3、13.5μg/m3、19μg/m3、24.5μg/m3、24.5μg/m3、19μg/m3、10.5μg/m3、8μg/m3、11μg/m3、17.5μg/m3、25μg/m3、3.5μg/m3、36μg/m3、34μg/m3、29.5μg/m3、19.5μg/m3、13μg/m3、10.5μg/m3和16μg/m3，从而得到上述单个空气污染检测站点外来污染传输前1日的污染物浓度增长曲线c；则上述单个空气污染检测站点的外来污染输送贡献是由单个空气污染检测站点外来污染输送过程的污染物浓度增长曲线b减去外来污染传输前的污染物浓度增长曲线c计算得出，即图2中所示阴影部分d；整个空气质量检测区域的外来污染输送贡献是由同一区域内所有空气污染检测站点的外来污染输送贡献之和然后算术平均得出，图中未示。
经过试用，本发明在不需要污染源信息的情况下，通过人工手动或电脑自动快速计算出突发污染源贡献和外来污染输送贡献，计算方法简单，计算结果准确，成本低，取得了良好的效果。
- 还没有人留言评论。精彩留言会获得点赞！