CN104850734B - 一种基于时空分布特征的空气质量指数预测方法 - Google Patents
一种基于时空分布特征的空气质量指数预测方法 Download PDFInfo
- Publication number
- CN104850734B CN104850734B CN201510188119.0A CN201510188119A CN104850734B CN 104850734 B CN104850734 B CN 104850734B CN 201510188119 A CN201510188119 A CN 201510188119A CN 104850734 B CN104850734 B CN 104850734B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- air quality
- mtd
- quality index
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种空气质量指数预测方法,尤其涉及一种基于时空分布特征的空气质量指数预测方法。本方法采用时间序列预测方法预测时间方向,即未来某个时间的空气质量指数;再采用克里格插值法,利用已知的监测站点经纬度坐标,结合时间序列预测结果,插值估计整个区域任意地点的空气质量指数。因此,本发明具有如下优点:1.降低了模型的复杂度,并且缩短了整体的计算时间,同时又保证了模型的准确性。2.可以更准确地预测区域内每个地方、一天内多个时段的空气质量指数。
Description
技术领域
本发明涉及一种空气质量指数预测方法,尤其涉及一种基于时空分布特征的空气质量指数预测方法。
背景技术
空气是地球上的生物赖以生存的物质,是必不可少的一种物质。环境空气质量与人们的日常生活息息相关,同时也在城市环境综合评价中占有重要地位。但是随着人类文明和经济的发展,空气污染越来越严重,如何改善空气质量、合理进行大气环境质量预测预警变得越来越重要,根据当地的空气质量预测结果,人们可以采取相应措施如带口罩,尽量避免外出等,保护自己免受空气污染物的侵害。
空气质量指数用AQI表示,是一个用来定量描述空气质量水平的数值。它由空气中二氧化硫(SO2)、二氧化氮(NO2)、一氧化碳(CO)、臭氧(O3)、悬浮颗粒物PM10和PM2.5,共6中污染物的含量计算出来的。AQI的取值范围位于0–500之间。由于空气质量情况与气候因素、地理位置有关,在数据角度上对应的就是空气质量指数与时间、空间分布都相关。因此需要有一定的技术方法,利用区域内多个监测站观测并计算得到的空气质量指数数据,预测整个地区未来的空气质量指数。
但是,传统的空气质量预测方法存在以下问题:
(1)不能准确的对没有观测站的区域进行空气指数预测:一般利用各地的多个空气质量监测站历史的检测数据,来预测该站点(或附近)未来一段时间的空气质量指数,而不能预测没有观测站的地方未来的空气质量情况;也有一种时空克里格插值方法可以仅根据离散分布的观测点的历史数据资料,预测整个区域范围内某个时间的空气质量指数,但是时空克里格插值法是内插法,即当要预测的未来时间点超出统计资料样本数据的最大时间时,预测误差将会变得很大。因此研究出新的方法可以比较准确的预测整个区域内每个地方在未来某个时间的空气质量指数是有必要的。
(2)不能进行细粒度时间段的空气指数预测:传统的空气质量预测通常是对接下来一整天的空气质量进行预测,这种粗粒度时间段的空气质量预测存在一个弊端,假设明天的空气质量预测结果是优,表明对空气污染物较敏感的人群可以进行户外活动,而实际上某些时间段(如上午8:00-9:00,下午5:00-6:00)的空气质量有可能是中度污染或者重度污染,对空气污染较敏感人群在这些时间段不宜进行户外活动。因此,对接下来一天的空气质量进行细粒度时间段的预测更为合理。
发明内容
本发明主要是解决现有技术所存在的不能准确的对没有观测站的区域进行空气指数预测以及不能进行细粒度时间段的空气指数预测的技术问题;提供了一种基于时空分布特征的空气质量指数预测方法。该方法将空间克里格插值法与时间序列预测方法相结合,把空间插值估计过程和时间序列预测过程分开进行,降低了模型的复杂度,并且缩短了整体的计算时间,同时又保证了模型的准确性;并且本方法可以更准确地预测区域内每个地方、一天内多个时段的空气质量指数,满足人们挑选空气质量较好的时间段进行外出活动的需要。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种基于时空分布特征的空气质量指数预测方法,包括以下步骤:
样本获取步骤,用于获取待预测区域中的空气质量监测站的经纬度坐标,每隔一定时间段记录区域中的各个空气质量监测站在同一时刻监测得到的空气质量指数,得到空气质量样本数据;
时间序列预测步骤,用于根据样本获取步骤中得到空气质量样本数据建立各个空气质量监测站的时间序列函数,选择置信度较高的时间序列函数预测该函数所对应的空气质量检测站在未来某个时间采样点的第一监测点空气指数预测集;
克里格插值步骤,用于根据第一监测点空气指数据预测集拟合得到待预测区域中的空间变异函数,根据所述空间变异函数以及第一监测点空气指数预测集进行克里格插值估计。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述时间序列预测步骤进一步包括以下子步骤:
时间序列划分子步骤,用于将空气质量样本数据划分为m个时间序列,每个时间序列包括来自同一监测站点的n个空气质量样本,其中:m是空气质量监测站的个数,n是空气质量监测站的数据记录次数;
模型训练子步骤,对于每个时间序列,根据已有样本中n个时间,选择最后的ι个时间对应的ι×m个样本作为验证样本,即实际值realij(i=1,2,…,m;j=1,2,…,ι);选择前n-ι个数据作为时间序列训练样本,利用训练样本训练得到各个站点的移动平均自回归模型,利用移动平均自回归模型预测验证样本中的ι个时间对应的空气质量指数,记为zij(i=1,2,…,m;j=1,2,…,ι),分别计算各个zij值的80%置信区间和90%的置信区间;
模型验证子步骤,对m站点的时间序列mi(i=1,2,…,m),若其ι个实际值分别落在zij的80%置信区间和90%的置信区间内部的数量比例分别达到预设的第一阈值和预设的第二阈值,则将该时间序列mi对应的移动平均自回归模型保留,共得到s个观测点的移动平均自回归模型;
时间序列预测子步骤,利用保留的s个观测点的移动平均自回归模型预测未来t时刻该s个站点的空气质量指数Z(xi)(i=1,2,…,s),将Z(xi)作为第一监测点空气指数预测集。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述第一阈值为85%,所述第二阈值为75%。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述克里格插值步骤进一步包括以下子步骤:
区域划分步骤,用于将待预测区域划分成不相交的正方形网格;
变异函数拟合子步骤,利用第一监测点空气指数预测集以及相应的监测点坐标,并基于以下公式计算得到待预测区域内的变异函数离散值:
式中,h为两监测点空间分隔距离,N(h)表示距离为h的监测点的对数,Z(xi)和Z(xi+h)分别监测点xi和监测点xi+h处第一监测点空气指数预测值,其中:i∈[1,2,...,N(h)];选择适当的圆滑曲线来对变异函数离散值进行拟合,将样本均方差值最小的函数作为待预测区域内的变异函数;
范围搜索子步骤,利用拟合得到的变异函数的变程a,搜索对网格中心处x0有影响的监测点,x0处经纬度坐标为(x0,y0),计算第一监测点空气指数预测集中的s个站点与x0的距离di(i=1,2,…,s),取di小于等于a的所有的有影响的观测点作为有效观测点构成有效观测点集合s0;
加权系数计算子步骤,利用变异函数、观测点集合s0并基于以下公式计算得到有效范围内的有效观测点集s0的空气质量指数在对某个网格X0进行克里格估值计算时的权值λi:
式中:,
γ(xi,xj)为有效观测点xi与xj之间的变异函数值,γ(xi,x0)为待估点x0与有效观测点xi之间的变异函数值;μ为拉格朗日乘数;
空间克里格估值子步骤,利用有效观测集合s0对应的第一监测点空气质量指数预测值Z(xi)(i=1,2,…,s0),并基于以下公式对网格X0进行预测:
式中,Z(xi)是网格X0的有效观测点集s0对应的第一监测点空气质量指数预测值,λi是s0个有效观测点对网格x0权值;对区域内其余的每个网格,重复范围搜索子步骤、加权系数计算子步骤和空间克里格估值子步骤,最后获得区域内各个网格的空气质量指数预测值。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
式中,C0为块金常数,(C0+C)为基台值,C为拱高,a为变程。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
式中,C0为块金常数,C为拱高,但a为三分之一的变程。
优化的,上述的一种基于时空分布特征的空气质量指数预测方法,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
式中,C0为块金常数,C为拱高,a为三分之一的变程。
因此,本发明具有如下优点:
1.本方法将空间克里格插值法与时间序列预测方法相结合,把空间插值估计过程和时间序列预测过程分开进行,降低了模型的复杂度,并且缩短了整体的计算时间,同时又保证了模型的准确性。
2.本方法可以更准确地预测区域内每个地方、一天内多个时段的空气质量指数,满足人们挑选空气质量较好的时间段进行外出活动的需要。
附图说明
附图1是本发明的流程图;
附图2是本发明的时间序列的移动平均自回归模型预测图;
附图3是本发明拟合的球形变异函数曲线。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:
图1是本发明的基于时空分布特征的空气质量指数预测方法。该方法的主要包括三个大的步骤,分别是样本获取步骤、时间序列预测步骤、克里格插值步骤,下面分别进行介绍。
1、样本获取
获取待预测区域中的空气质量监测站的经纬度坐标,每隔一定时间段记录区域中的各个空气质量监测站在同一时刻监测得到的空气质量指数,得到空气质量样本数据;本实施例中,各采样点之间的时间间隔优选为6小时,为保证训练样本完整性,监测站点有数据历史数据有缺失情况,填补缺失记录,获取n×m个样本。
2、时间序列预测
根据样本获取步骤中得到空气质量样本数据建立各个空气质量监测站的时间序列函数,选择置信度较高的时间序列函数预测该函数所对应的空气质量检测站在未来某个时间采样点的第一监测点空气指数预测集;该步骤进一步包括以下子步骤:
2.1时间序列划分
用于将空气质量样本数据划分为m个时间序列,每个时间序列包括来自同一监测站点的n个空气质量样本,其中:m是空气质量监测站的个数,n是空气质量监测站的数据记录次数;
2.2模型训练
对于每个时间序列,根据已有样本中n个时间,选择最后的ι个时间对应的ι×m个样本作为验证样本,即实际值realij(i=1,2,…,m;j=1,2,…,ι);选择前n-ι个数据作为时间序列训练样本,利用训练样本训练得到各个站点的移动平均自回归模型作为该站点的时间序列函数。
移动平均自回归模型,简称为ARIMA(p,d,q)模型,是目前时间序列的众多模型中比用通用,且预测效果较佳的模型。该模型含有p,d,q三个参数,通过调整3个参数可保证时间序列的平稳性。因此对以上m个时间序列,采用该模型时p,d,q三个参数各不相同。
利用训练得到的移动平均自回归模型预测验证样本中的ι个时间对应的空气质量指数,记为zij(i=1,2,…,m;j=1,2,…,ι),其中,zij的80%置信区间上下限分别为80_lowerij和80_upperij,95%置信区间上下限分别为95_lowerij和95_upperij(i=1,2,…,m;j=1,2,…,ι)。
2.3模型验证
对站点的时间序列mi(i=1,2,…,m),计算ι个实际值分别落在80%置信区间内部的数量比例80_ratioi和落在95%置信区间内部的数量比例95_ratioi:
当80_ratioi>=80%,并且95_ratioi>=75%,就认为时间序列mi(i=1,2,…,m)的移动平均自回归模型ARIMA的预测精度较高,并且对应的观测点予以保留,m个观测点中保留的有s个,并记录这s个观测点的经纬度坐标为(sxi,syi)i=1,2,…,s)。
如图2所示,是一个时间序列的ARIMA模型的预测值、置信区间和实际值的比较图,此ARIMA模型的3个参数分别是p=1,d=q=0,右边阴影部分中的线条是预测结果,阴影部分分别表示80%置信区间和95%置信区间,圆圈代表实际值,此图比较直观的展示了ARIMA模型的预测效果。
2.4时间序列预测
利用保留的s个观测点的移动平均自回归模型预测未来t时刻该s个站点的空气质量指数Z(xi)(i=1,2,…,s),将Z(xi)作为第一监测点空气指数预测集。
3、克里格插值
根据第一监测点空气指数据预测集拟合得到待预测区域中的空间变异函数,根据所述空间变异函数以及第一监测点空气指数预测集进行克里格插值估计。
在知道区域内离散分布的观测点的同一时刻的空气质量指数时,要估计整个区域任意地方的空气质量指数。由于空气质量情况是有空间分布特征的,因此空气质量指数也有空间分布特征,对于有空间分布特征的数据,采用克里格插值法可以无偏估计出区域内任意点的空气质量指数。主要步骤如下:
3.1区域划分
将待预测区域划分成不相交的正方形网格,在同一时刻每个网格内的各处的空气质量指数相等。克里格插值过程中网格单元越小,插值估计经度也越高,因此在划分网格时,网格的面积很小,使每个网格内部的空气质量指数近似相等的。
区域G的经纬度的最大最小值记为:经度最小值xmin,经度最大值xmax,纬度最小值ymin,纬度最大值ymax。每个网格g的经纬跨度记为cell,那么按照cell为单元来划分,区域G的经度可以划分为XX=(xmax-xmin)/cell份,其中XX取整数,经度方向可得到xmin,xmin+cell,xmin+2*cell,…,xmin+XX*cell序列;区域G的纬度可以划分为YY=(ymax-ymin)/cell份,其中YY取整数。
纬度方向可划分为ymin,ymin+cell,ymin+2*cell,…,ymin+YY*cell序列;将这些经纬度序列进行组合,总共有N=XX*YY个经纬度坐标,即空间坐标。这些坐标依次记为(xi,yi)(i=1,2,…,N),每个坐标点是网格g的中心,最终目标就是要预测每个格网g的中心点处的空气质量指数。
3.2变异函数拟合
空间变异的意思是,变量在空间中随着空间位置的变化而不同;变异函数是描述这种变异性质的表达式。
利用第一监测点空气指数预测集以及相应的监测点坐标,并基于以下公式计算得到待预测区域内的变异函数离散值:
式中,h为两监测点空间分隔距离,N(h)表示距离为h的监测点的对数,Z(xi)和Z(xi+h)分别监测点xi和监测点xi+h处第一监测点空气指数预测值,其中:i∈[1,2,...,N(h)];
为了考察空间位置x处的值Z(x)受多大范围内的空间点值的影响,所以要对上面离散的变异函数γ*h进行拟合,采用加权多项式回归拟合法进行拟合,拟合的光滑的变异函数曲线,与空间点位置xi无关,只与空间点距离h有关。
本实施例中的变异函数的拟合可采用以下几种方法:
(1)球形变异函数模型:
式中:C0为块金常数,(C0+C)为基台值,C为拱高,a为变程。当空间点的距离h>=变程a时,变异函数γ(h)值不再变化,即以点x处为中心,半径为a范围内的点xi处的值Z(xi)对x处的值Z(x)有影响,超出这个范围的无影响。也就是说克里格差值过程中,只选择与所预测的网格距离小于等于变程的观测点。
采用球形变异函数模型进行拟合的结果如图3所示。
(2)指数模型:
式中:C0,C意义同前,但a不是变程,由于1-e-3=1-0.05=0.95≈1,则变程为3a
(3)高斯模型:
式中:C0,C意义同前,但a不是变程,由于1-e-3=1-0.05=0.95≈1,则变程为3a。
3.3范围搜索
利用拟合得到的变异函数的变程a,搜索对网格中心处x0有影响的监测点,x0处经纬度坐标为(x0,y0),计算第一监测点空气指数预测集中的s个站点(经纬度坐标为(sxi,syi))与x0的距离di(i=1,2,…,s),取di小于等于a的所有的有影响的观测点作为有效观测点构成有效观测点集合s0。
3.4加权系数计算
利用变异函数、观测点集合s0并基于以下公式计算得到有效范围内的有效观测点集s0的空气质量指数在对某个网格X0进行克里格估值计算时的权值λi:
式中:γ(xi,xj)为有效观测点xi与xj之间的变异函数值,γ(xi,x0)为待估点x0与有效观测点xi之间的变异函数值;μ为拉格朗日乘数;此公式仅为置信度高且在网格的有效范围内的监测站对一个网格的权值计算,对每个网格的权值需重复计算该步骤,其中拉格朗日常数,也称为朗格朗日乘子,是一种寻找变量受一个或多个条件所限制的多元函数的极值的方法。这种方法将一个有n个变量与k个约束条件的最优化问题转换为一个有n+k个变量的方程组的极值问题,因而可用求解极值方法求解条件方程组。
3.5空间克里格估值
利用有效观测集合s0对应的第一监测点空气质量指数预测值Z(xi)(i=1,2,…,s0),并基于以下公式对网格X0进行预测:
式中,Z(xi)是网格X0的有效观测点集s0对应的第一监测点空气质量指数第一监测点空气质量指数预测值,λi是s0个有效观测点对网格x0权值;对区域内其余的每个网格,重复范围搜索子步骤、加权系数计算子步骤和空间克里格估值子步骤,最后获得区域内各个网格的空气质量指数预测值。
通过以上方案,本方法采用空间克里格插值法与时间序列预测方法相结合,把空间插值估计过程和时间序列预测过程分开进行,降低了模型的复杂度,并且缩短了整体的计算时间。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了等术语,但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质;把它们解释成任何一种附加的限制都是与本发明精神相违背的。
Claims (7)
1.一种基于时空分布特征的空气质量指数预测方法,其特征在于,包括以下步骤:
样本获取步骤,用于获取待预测区域中的空气质量监测站的经纬度坐标,每隔一定时间段记录区域中的各个空气质量监测站在同一时刻监测得到的空气质量指数,得到空气质量样本数据;
时间序列预测步骤,用于根据样本获取步骤中得到空气质量样本数据建立各个空气质量监测站的时间序列函数,选择置信度较高的时间序列函数预测该函数所对应的空气质量检测站在未来某个时间采样点的第一监测点空气指数预测集;
克里格插值步骤,用于根据第一监测点空气指数据预测集拟合得到待预测区域中的空间变异函数,根据所述空间变异函数以及第一监测点空气指数预测集进行克里格插值估计。
2.根据权利要求1所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述时间序列预测步骤进一步包括以下子步骤:
时间序列划分子步骤,用于将空气质量样本数据划分为m个时间序列,每个时间序列包括来自同一监测站点的n个空气质量样本,其中:m是空气质量监测站的个数,n是空气质量监测站的数据记录次数;
模型训练子步骤,对于每个时间序列,根据已有样本中n个时间,选择最后的ι个时间对应的ι×m个样本作为验证样本,即实际值realij,i=1,2,…,m;j=1,2,…,ι;选择前n-ι个数据作为时间序列训练样本,利用训练样本训练得到各个站点的移动平均自回归模型,利用该移动平均自回归模型预测验证样本中的ι个时间对应的空气质量指数,记为zij,i=1,2,…,m;j=1,2,…,ι,分别计算各个zij值的80%置信区间和90%的置信区间;
模型验证子步骤,对m站点的时间序列ti′,i=1,2,…,m,若其ι个实际值分别落在zij的80%置信区间和90%的置信区间内部的数量比例分别达到预设的第一阈值和预设的第二阈值,则将该时间序列ti′对应的移动平均自回归模型保留,共得到s个观测点的移动平均自回归模型;
时间序列预测子步骤,利用保留的s个观测点的移动平均自回归模型预测未来t时刻该s个站点的空气质量指数Z(xi),i=1,2,…,s,将Z(xi)作为第一监测点空气指数预测集。
3.根据权利要求2所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述第一阈值为85%,所述第二阈值为75%。
4.根据权利要求1所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述克里格插值步骤进一步包括以下子步骤:
区域划分步骤,用于将待预测区域划分成不相交的正方形网格;
变异函数拟合子步骤,利用第一监测点空气指数预测集以及相应的监测点坐标,并基于以下公式计算得到待预测区域内的变异函数离散值:
<mrow>
<msup>
<mi>&gamma;</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
式中,h为两监测点空间分隔距离,N(h)表示距离为h的监测点的对数,Z(xi)和Z(xi+h)分别监测点xi和监测点xi+h处第一监测点空气指数预测值,其中:i∈[1,2,...,N(h)];选择适当的圆滑曲线来对变异函数离散值进行拟合,将样本均方差值最小的函数作为待预测区域内的变异函数;
范围搜索子步骤,利用拟合得到的变异函数的变程a,搜索对网格中心处x0有影响的监测点,x0处经纬度坐标为(x0,y0),计算第一监测点空气指数预测集中的s个站点与x0的距离di,i=1,2,…,s,取di小于等于a的所有的有影响的观测点作为有效观测点构成有效观测点集合s0;
加权系数计算子步骤,利用变异函数、观测点集合s0并基于以下公式计算得到有效范围内的有效观测点集s0的空气质量指数在对某个网格X0进行克里格估值计算时的权值λi:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
<mo>=</mo>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,γ(xi,xj)为有效观测点xi与xj之间的变异函数值,γ(xi,x0)为待估点x0与有效观测点xi之间的变异函数值;μ为拉格朗日乘数;
空间克里格估值子步骤,利用有效观测集合s0对应的第一监测点空气质量指数预测值Z(xi),i=1,2,…,s0,并基于以下公式对网格X0进行预测:
<mrow>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
</munderover>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>i</mi>
<mo>&CenterDot;</mo>
</mrow>
</msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中,Z(xi)是网格X0的有效观测点集s0对应的第一监测点空气质量指数预测值,λi是s0个有效观测点对网格X0权值;对区域内其余的每个网格,重复范围搜索子步骤、加权系数计算子步骤和空间克里格估值子步骤,最后获得区域内各个网格的空气质量指数预测值。
5.根据权利要求4所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>h</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>a</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mi>a</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>h</mi>
<mo>&le;</mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>C</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>></mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,C0为块金常数,C0+C为基台值,C为拱高,a为变程。
6.根据权利要求4所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mi>h</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,C0为块金常数,C为拱高,但a1为三分之一的变程。
7.根据权利要求4所述的一种基于时空分布特征的空气质量指数预测方法,其特征在于,所述变异函数拟合子步骤中变异函数的拟合基于以下公式:
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mrow>
<mi>a</mi>
<msup>
<mn>1</mn>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,C0为块金常数,C为拱高,a1为三分之一的变程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510188119.0A CN104850734B (zh) | 2015-04-21 | 2015-04-21 | 一种基于时空分布特征的空气质量指数预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510188119.0A CN104850734B (zh) | 2015-04-21 | 2015-04-21 | 一种基于时空分布特征的空气质量指数预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104850734A CN104850734A (zh) | 2015-08-19 |
CN104850734B true CN104850734B (zh) | 2017-09-15 |
Family
ID=53850375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510188119.0A Active CN104850734B (zh) | 2015-04-21 | 2015-04-21 | 一种基于时空分布特征的空气质量指数预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104850734B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550784B (zh) * | 2016-01-20 | 2020-01-14 | 中科宇图科技股份有限公司 | 一种空气质量监测站优化布点方法 |
CN106228010A (zh) * | 2016-07-22 | 2016-12-14 | 上海海洋大学 | 一种北太平洋鱿鱼资源补充量预测方法 |
CN106770532B (zh) * | 2016-12-22 | 2019-01-08 | 广州市真量纯科技有限公司 | 一种气体浓度预测装置、及气体浓度实时检测和预测*** |
CN107492130B (zh) * | 2017-01-09 | 2020-06-23 | 辽宁工程技术大学 | 一种基于层面渲染的非连续等值云图生成方法 |
CN107239592B (zh) * | 2017-04-24 | 2019-11-29 | 清华大学 | 一种空气污染监测数据的时间序列增长分析方法及*** |
CN107798425A (zh) * | 2017-10-16 | 2018-03-13 | 中国科学院地理科学与资源研究所 | 一种基于大数据的时空混淆暴露度评估***及方法 |
CN108120661B (zh) * | 2017-12-19 | 2021-02-05 | 北京理工大学 | 一种城市空气中颗粒物含量时空分布测定方法 |
CN108053071A (zh) * | 2017-12-21 | 2018-05-18 | 宇星科技发展(深圳)有限公司 | 区域空气污染物浓度预测方法、终端及可读存储介质 |
CN107992454B (zh) * | 2017-12-25 | 2020-11-17 | 安徽大学 | 一种基于在线顺序回归的空气质量等级预测方法 |
CN108062454A (zh) * | 2018-01-19 | 2018-05-22 | 宁波市镇海规划勘测设计研究院 | 污染物时空分布不确定性特征分析方法、***及存储介质 |
CN108564263B (zh) * | 2018-04-02 | 2021-10-15 | 国网安徽省电力有限公司电力科学研究院 | 一种用于飑线环境下输电线路杆塔受灾预测方法 |
WO2020043029A1 (zh) * | 2018-08-25 | 2020-03-05 | 司书春 | 一种通过向量化分析预测大气污染的方法 |
CN113988351A (zh) * | 2020-07-27 | 2022-01-28 | 海南省气象科学研究所 | 一种基于时空局部模型的数值模式气温预报产品释用方法 |
CN112541302B (zh) * | 2020-12-23 | 2024-02-06 | 北京百度网讯科技有限公司 | 空气质量预测模型训练方法、空气质量预测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799772A (zh) * | 2012-07-03 | 2012-11-28 | 中山大学 | 面向空气质量预报的样本优化方法 |
CN103514366A (zh) * | 2013-09-13 | 2014-01-15 | 中南大学 | 一种城市空气质量浓度监测缺失数据的修复方法 |
KR101382507B1 (ko) * | 2012-10-19 | 2014-04-10 | 사단법인대기환경모델링센터 | 대기질 예측 및 관리 시스템 |
CN104200103A (zh) * | 2014-09-04 | 2014-12-10 | 浙江鸿程计算机***有限公司 | 一种基于多领域特征的城市空气质量等级预测方法 |
CN104504475A (zh) * | 2014-12-24 | 2015-04-08 | 上海大学 | 基于ar*-svm混合建模的雾霾时间序列预测方法 |
-
2015
- 2015-04-21 CN CN201510188119.0A patent/CN104850734B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799772A (zh) * | 2012-07-03 | 2012-11-28 | 中山大学 | 面向空气质量预报的样本优化方法 |
KR101382507B1 (ko) * | 2012-10-19 | 2014-04-10 | 사단법인대기환경모델링센터 | 대기질 예측 및 관리 시스템 |
CN103514366A (zh) * | 2013-09-13 | 2014-01-15 | 中南大学 | 一种城市空气质量浓度监测缺失数据的修复方法 |
CN104200103A (zh) * | 2014-09-04 | 2014-12-10 | 浙江鸿程计算机***有限公司 | 一种基于多领域特征的城市空气质量等级预测方法 |
CN104504475A (zh) * | 2014-12-24 | 2015-04-08 | 上海大学 | 基于ar*-svm混合建模的雾霾时间序列预测方法 |
Non-Patent Citations (1)
Title |
---|
北京地区空气质量指数时空分布特征及其与气象条件的关系;张建忠 等;《气象与环境科学》;20140331;第37卷(第1期);第33-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104850734A (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104850734B (zh) | 一种基于时空分布特征的空气质量指数预测方法 | |
US10877012B2 (en) | Automatic siting for air quality monitoring stations | |
Delle Monache et al. | Kalman filter and analog schemes to postprocess numerical weather predictions | |
CN110264709A (zh) | 基于图卷积网络的道路的交通流量的预测方法 | |
CN113919231B (zh) | 基于时空图神经网络的pm2.5浓度时空变化预测方法及*** | |
CN109308543A (zh) | 基于ls-svm和实时大数据的地铁短期客流预测方法 | |
CN107103754A (zh) | 一种道路交通状况预测方法及*** | |
CN111665575B (zh) | 一种基于统计动力的中长期降雨分级耦合预报方法及*** | |
CN107563574A (zh) | 一种基于长短时记忆模型循环神经网络的地磁场预测方法 | |
CN105760970A (zh) | 空气质量指数的预测方法 | |
CN107610021A (zh) | 环境变量时空分布的综合分析方法 | |
CN113177657B (zh) | 轨道交通客流预测方法及装置 | |
CN106407633A (zh) | 基于时空回归克里金模型估算地面pm2.5的方法及*** | |
CN109376935A (zh) | 一种基于神经网络的公交客流分时段组合预测方法 | |
CN102253426A (zh) | 采用空间密度聚类进行雷电运动预测的方法 | |
CN109272169A (zh) | 交通流量预测方法、装置、计算机设备及存储介质 | |
CN105974495A (zh) | 利用分类拟合法预判目标区域未来平均云量的方法 | |
Si et al. | Effects of urbanization on daily temperature extremes in North China | |
CN114494865A (zh) | 一种山地生态***未来情景的cset方法 | |
CN115544889A (zh) | 基于深度学习的数值模式降水偏差校正方法 | |
CN108564228A (zh) | 一种基于时序特征预测轨道交通od客流量的方法 | |
CN109190797A (zh) | 一种基于改进k最近邻的大规模路网状态短时预测方法 | |
CN109615414A (zh) | 房产预估方法、装置及存储介质 | |
Lin et al. | Building autocorrelation-aware representations for fine-scale spatiotemporal prediction | |
CN107194507A (zh) | 一种基于组合支持向量机的风电场短期风速预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 430000 Wuda science and Technology Park, Jiangxia Avenue, Miaoshan District, Donghu Development Zone, Wuhan City, Hubei Province Patentee after: Geospace Information Technology Co., Ltd. Address before: 430223 Wuda science and Technology Park, Jiangxia Avenue, Miaoshan community, Donghu Development Zone, Wuhan City, Hubei Province Patentee before: WUDA GEOINFORMATICS Co.,Ltd. |
|
CP03 | Change of name, title or address |