[发明专利]一种基于时空特征的实时估算细颗粒物浓度的方法

说明书

技术领域

本发明涉及实时地细颗粒物浓度估算，更具体地涉及基于时空特征的实时细颗粒物浓度估算方法。

背景技术

近年来由于雾霾天气不断增多，人们对空气质量健康越来越关注。细颗粒物(PM2.5)浓度是评估空气质量的重要指标之一。同时细颗粒物具有粒径小，活性强，可吸入性，易附带有毒有害物质的特征，因而对人体健康的影响更大。如何快速实时获取细颗粒物浓度成了人们日常生活的迫切需要。虽然政府部门在城市内设置了空气质量监测站点，但是由于土地成本，仪器成本，人力成本的各方面因素的综合考虑，无法在城市内部布设过多的空气质量监测站点。因此在远离监测站点的区域，获取实时的可靠的空气污染物浓度变得十分困难。

细颗粒污染物浓度是一种时空现象，但是传统的基于空间特征的估计细颗粒物浓度的方法只考虑到污染物在空间上的相关性，比如反距离法和普通克里金法。反距离法仅仅只以相互之间距离的倒数值为已知浓度站点对未知浓度站点的权重系数从而估算未知浓度站点的浓度值。普通克里金法只考虑到细颗粒物浓度在空间的分布特征，在此基础上求算已知站点对未知站点影响的权重值，最终估算未知点细颗粒物浓度值。本发明在考虑空间分布的基础上同时考虑到细颗粒物浓度在时间维度的变异特征，并在时空变异的基础上引入与细颗粒物浓度相关的三大气象因子：湿度、风速和大气压强，进而估算未知站点处的细颗粒物浓度。

发明内容

本该发明弥补了传统空间估算方法的不足，在考虑时空相关性的基础上，同时加入风速、湿度和大气压强三个气象因子加权计算细颗粒物浓度，在监测站点无法覆盖的区域可以提供实时可靠的污染物浓度值，为人们出行提供空气质量参考信息。

为实现上述目的，本发明提供一种基于时空特征的实时估算细颗粒物浓度的方法，其特征在于：包含3个步骤，分别为数据获取与预处理步骤、拟合时空变异曲线步骤，估算细颗粒物浓度步骤；

(1)所述数据获取与预处理步骤包括：

获取以当前小时为起点向前72小时n个空气质量监测站点的细颗粒物浓度，以及各个站点的地理位置信息，共72*n组数据；

获取当前小时m个气象监测站点的湿度、风速和大气压强值，以及气象站点地理位置信息；

将m个气象监测站点构建泰森多边形，将每个空气质量监测站点与气象监测站点构成的泰森多边形一一匹配，为每个空气监测站点赋予湿度、风速以及大气压强值；用x，y值表示空气质量监测站点的空间位置，M₁、M₂和M₃分别表示该空气质量监测站点的湿度、风速和大气压强值；

由于细颗粒物浓度、湿度、风速和大气压强值具有不同的量纲，为消除数据不同量纲的影响性，使数据具有可比性，使用极差法对数据进行标准化处理，极差法公式如下：

————公式(1)

式中x^*为标准化的新数值，x_max为样本数据最大值，x_min为样本数据最小值。

(2)所述拟合时空变异曲线步骤包括：

对当前小时n个空气监测站点细颗粒物浓度进行空间变异函数拟合，以空气监测站点间空间距离为自变量，具体为：

●计算各个空气质量监测站点间的距离值以及变异函数值，计算公式分别如下：

————公式(2)

式中(x_i，y_i)为监测站点i的空间位置，(x_j，y_j)为监测站点j的空间位置；

————公式(3)

式中Z(i)和Z(j)分别为站点i和站点j细颗粒物浓度值；

●以距离h为横坐标，空间变异函数γ_s(h_s)值为纵坐标，绘制空间变异函数散点图；

●对上述空间变异函数散点图进行球状模型拟合，拟合公式如下：

————公式(4)

式中C₀，a，C均为常量

对72小时n个空气监测站点细颗粒物浓度进行时间变异函数拟合，以时间间隔为自变量，具体为：

●由于72小时细颗粒物浓度是一组具有周期的时间序列值，可分解为以下公式：Z(t)＝μ(t)+s(t)+ε(t)————公式(5)

其中，μ(t)为趋势项，s(t)为周期项，ε(t)为随机误差项；

●对72小时细颗粒物浓度数据进行时序分解操作，去除周期项s(t)；