网站首页  词典首页

请输入您要查询的论文:

 

标题 时间序列分析方法在郑州市降水量预报中的应用
范文 基金项目: 河南平原地区地下水污染调查评价(淮河流域):国土资源大调查项目(1212010634502)
作者简介:吕志涛(1968),男,河南禹州人,高级工程师,主要从事气候变化与水文地质工作。 Email:lztzmd@sina.comDOI:10.13476/j.cnki.nsbdqk.
摘要:根据1971年-2013年郑州市的降水量资料,采用二次多项式拟合提取降水量的趋势分量,采用谐波分析法提取降水量的周期成分,利用自回归模型求解随机成分,最后将三者叠加,构建了郑州市降水量的预报模型。模型计算结果与实测数据对比可知,应用预报模型对降水量进行预报精度较高。因此利用建立的预报模型对2014年-2016三年的降水量进行了预测,为该区水资源的管理提供依据。
关键词:降水量;时间序列分析;预报;郑州市
中图分类号:TV121.1文献标志码:A文章编号:16721683(2014)04003503
Application of Time Series Analysis on the Annual Precipitation of Zhengzhou city
LVU Zhitao1,2
(1.College of Resources and Environment,North China University of Water Resources and
Electric Power,Zhengzhou 450011,China;2.Henan Institute of Meteorological Sciences,Zhengzhou 450003,China)
Abstract:Time series analysis method has comprehensive importance and prospects in precipitation forecast.After analyzing precipitation data from 1971 to 2009 in Zhengzhou city,precipitation forecasting model was constructed by time series analysis method in the paper.When establishing forcasting model,First,the trend component of precipitation is picked up by quadratic polynomial calibration,the periodic component is extracted by spectrum analysis and the stochastic component is simulated by using autoregressionmodel.Finally,the forecasting model is established through superposition of these components,and the method has been used to forecast the next three years precipitation in the Zhengzhou city.In comparison with model calculation results and measured data,the precipitation forecast model accuracy meets the requirements,we can use the model to forecast precipitation in the future,this paper gives the precipitation prediction results of three years of 2014~2016.The results would be helpful for water management in this region.
Key words:precipitation;time series analysis;forecast;Zhengzhou city
降水量预报方法很多,如多年平均值法、皮尔逊Ⅲ型概率分布曲线等。上述降水量预报方法,虽然在某程度上反映了地区降水量分布特征,但还不能反映降水分布的本质特征。实际上,一个地区降水量的年内分布具有随机性、不确定性,然而其多年分布却存在统计规律性[1]。人们把某一变量或指标的数值或观测值,按照其出现时间的先后次序,以相同的间隔时间排列一组数值,称为时间序列。目前,分析时间序列的演变规律是掌握系统动力学特性的重要手段。该方法易于掌握,计算工作量小,易于应用推广[2],因此在国民经济宏观控制、区域综合发展规划、企业经营管理、市场潜量预测、气象预报、水文预报、地震前兆预报、农作物病虫灾害预报、环境污染控制、生态平衡、天文学和海洋学等领域已有广泛应用。
近年来,时间序列分析法在降水量预报方面取得了十分显著的效果。董晓萌等应用时间序列模型过对杨凌张家岗气象站1935年-1955年的四季降水量建立了自回归模型,利用所求得的模型对1956年的四季降水量进行了预报[3];许晓艳根据浑河流域长系列的降水量资料,运用时间序列分析原理建立了浑河流域降水量预测模型,预测该流域的各站降水量,经检验其预测的中长期结果具有较高的精度[4]。本文在分析1971年-2013年郑州市降水量(面雨量)的年际变化特征的基础上,应用时间序列分析法对该区的降水量进行了预报,力求为该区水资源的管理提供依据。
1时间序列分析方法的基本原理
1.1预报模型
时间序列分析法也称历史资料延伸预测法,通过对已观测的一系列历史数据进行曲线拟合和参数估计,掌握该序列的发展过程、方向以及趋势,从而对其进行类推或延伸,并来预测下一段时间可能达到的水平。时间序列可分解成趋势成分、周期成分和随机性成分,每一种成分都可以通过一定的数学方法来确定;将三者线性叠加,即可得到总的预报模型;经检验合格后,可用于预报[5]。预报模型可用下式表示:
H(t)=X(t)+P(t)+R(t)(1)
式中:X(t)为趋势项,反映H(t)随时间的变化趋势;P(t)为周期项,反映H(t)的周期性变化;R(t)为随机项,反映随机要素的影响。
1.2趋势项分析
计算趋势成分的方法主要有两种:平滑法(移动平均法、指数平均法)和多项式拟合法。本文采用多项式拟合法进行趋势项分析,表达式为:
t=C0+C1t+C2t2+…+Cktk=∑kk=0Cktk(2)
式中:t为数据的趋势项,C0、C1、C2、…、Ck为常系数,t为实测数据。
可采用多元回归方法确定待定系数。一般采用二次多项式Xt=b0+b1t+b2t2拟合趋势成分,应用最小二乘法确定多项式系数,即可确定Xt的表达式。为检验拟合结果,需在一定的显著性水平下计算趋势曲线拟合的相关系数R,若无最佳拟合函数则认为无趋势项或趋势项不显著。
R=1-∑nt=1(xt-t)2∑nt=1(xt-x)2(3)式中:n为实测数据序列的总个数,Xt为实测序列,t为计算序列,X为实测序列的平均值。R愈接近1时,表明y与xi(i=1,2,…,k)的线性关系愈密切,反之愈差。
1.3周期项分析
本文采用谐波分析法,将剩余序列Pt(t=1,2 ,…,n)看成是由不同周期的规则波叠加而成。因而在分离周期时,逐步分解出一些比较明显的周期波,然后叠加起来作为该时间序列的周期项,把这个周期项进行外推可用于预报。对序列Pt可用L个波叠加的形式表示其估计值为:
t=a02+∑Lk=1[akcos2πktn+bksin2πktn](4)
式中:t为序列Pt的估计值;L为谐波个数,取n/2的整数部分;K为谐波序号,(k=1,2,…,L);ak,bk为傅里叶系数,其计算式为:
Vn+2=Vn+1=0
Vt=x(t)+Vt+1·2cos(ωj)-Vt+2
(t=n,n-1,…,2,1)(5)
a0=1n∑nt=1x(t)ak=2n[V1cos(ωj)-V2]
bkj=2n[V1sin(ωj)](j=1,2,…,k) (6)
通常在L个波中选取波动比较显著的几个谐波相加来估计Pt,在实际应用中只需选取前6个显著谐波即可满足精度要求。若满足式(7),则认为第k个波显著,否则不显著。
sk2=ak2+bk2>4s2lnkan(7)
式中:α为检验的显著性水平(一般取5%) ;s2为系列的方差。
1.4随机项分析
消除趋势项和近似周期项后的剩余序列为平稳随机系列项R(t), R(t)=H(t)-X(t)-P(t),可直接对其用自回归模型求解。考虑到R(t)与t时刻的取值和它前1个到p个时间间隔t-1, t-2,…,t-p的取值有关,其自回归模型为:
t=Φ0+Φ1rt-1+Φ2rt-2+…+Φprt-p(8)
式中:t为估计值;Φi为模型自回归系数;p为模型阶数,用AIC(p)准则确定。
当AIC(p)值最小时,所对应的p值为最佳阶数。
AIC(p)=nlnp2+2p(9)
式中:n为序列数据总个数; p2为AR(p)残差的方差。
1.5预报模型的建立和检验
将上述趋势分量、周期分量、随机分量线性叠加,即可得到降水量的总预测模型[68]:
t=t+t+t(10)
建立的模型要经过检验才能判定其是否合理,是否合格。本文采用后验差方法来进行检验模型的精度,即通过后验差比值c和小误差频率P来检验。设动态样本中前n个数据的标准差为S2,后验数据残差为ej,残差均值为,则后验差比值c及小误差频率p的计算公式为:
c=s2s1(11)
p={|ej-e|<0.6745s1}(12)
表1后验预测法预测效果表[9]
Table 1Effects of posterior prediction
预测效果好合格勉强不合格p>0.95>0.80>0.70≤0.70c<0.35<0.50<0.65≥0.65如果p、c值都在允许范围内见表1,则模型可用于计算预报值,否则需要对模型进行检查、分析和重新调参。
2时间序列分析在降水预报中的应用
本文利用的郑州市1971年-2013年降水量统计资料来源于郑州国家气候观象台。其中1971年-2005年35个降水量数据用于建立数学模型, 2006年-2013年8个降水量数据用于检验数学模型的精度。
2.1趋势项的提取
采用二次多项式拟合法,根据前文的趋势项提取方法,可得趋势项函数为:
t=726.6591-13.2598t+0.3567t2(13)
为检验趋势拟合效果,计算趋势曲线拟合的相关系数R:
图1趋势项计算值与实测值比较
Fig.1Comparison of Trend term calculated and measured values
R=S回S总=1-∑nt=1(xt-t)2∑nt=1(xt-)2=0.1975(14)
显然,趋势项并不明显,可以忽略。
2.2周期项的提取
取L=(35-1)/2=17。计算系数a0,ak,bk(k=1,2,…,17),其中a0=-0.014705927。根据周期成分系数计算结果,当k=5、7、16时,波显著。则模型为:
t=a02+a5cos(2π35·5t)+b5sin(2π35·5t)+a7cos(2π35·7t)+b7sin(2π35·7t)+a16cos(2π35·16t)+b16sin(2π35·16t)(15)
2.3平稳随机项的拟合
平稳随机项可以通过自回归模型AR(P)来拟合。用AIC(p)准则确定阶数,最终确定 p=6。平稳随机项为:
r(t)=h(t)-p(t)(16)
t=Φ0+Φ1rt-1+Φ2rt-2+Φ3rt-3+Φ4rt-4+Φ5rt-5+Φ6rt-6
=2536.2305-0.0868rt-1-0.167rt-2-1.5880rt-3-0.1890rt-4-0.463rt-5-0.98282rt-6(17)
因此,把周期项和平稳随机项项叠加起来,得到降水量的预报方程为:
ht=a0+a5cos(2πn·5t)+b5sin(2πn·5t)+a7cos(2πn·7t)+b7sin(2πn·7t)+a16cos(2πn·16t)+b16sin(2πn·16t)+Φ0+Φ1rt-1+Φ2rt-2+Φ3rt-3+Φ4rt-4+Φ5rt-5+Φ6rt-6 (18)
从图2中模型计算值和实测值的比较可知,拟合效果较好,说明建立的模型较符合实际。
2.4模型检验及预报
令表示后验预测的残差,即:
ej=hn+1-hn+j(j=1,2,3,4,5)
设动态样本中1971年-2005年35个数据的标准差为,2006年-2013年8个数据的后验预测残差的标准差为s2。
根据公式计算得:
s1=134∑34t=1(ht-h)2=61.0572
s2=15∑5j=1(ej-e)2=18.5718
c=s2s1=0.3042
p={|ej-e|<0.6745s1}=1
图2降水量实测值与计算值比较
Fig.2Comparison of precipitation calculated and measured values
表2后验预测误差表
Table 2Error of posterior prediction
n实测值/mm计算值/mmej|ej-e|36692692.23 -0.232.2937596607.75 -11.7513.8138658643.5114.4912.4339764738.33 25.6723.6140600.3589.1511.159.0941706.5691.3715.1313.0742498.7524.31-25.6127.6743390.5417.29-26.7928.85查表1可知,后验预测效果等级为“好”,可以进行降水量的预报。2014年-2016年的降水量预测值见表3。
表3降水量预报结果
Table 3Prediction of precipitation
时间/年201420152016降水量/mm687.264705.905547.8803结语
本文运用时间序列分析方法,通过对降水量趋势、周期项和随机项的提取,建立了郑州市降水量和时间之间的数学模型。在模型检验表明,建立的数学模型的预报精度较高,但随着预测周期的增加,预报误差也在逐渐增大。因此为使预测更准确,预报时间长度不宜太长。参考文献(References):
[1]汤成友,官学文,张世明.现代中长期水文预报方法及其应用[M].北京:中国水利水电出版社,2008.(TANG Chengyou,GUANG Xuewen,ZHANG Shiming.Modern Midtolongterm Hydrological Forecasting Method and its Application[M].China Water Power Press.2008.(in Chinese))
[2]Martinez J C ,Maleski J J,Miller F M.Trends in Precipitation and Temperature in Florida[J].USA J Hydrol,2012,452453,259281.
[3]董晓萌,罗凤娟,郭满才,等.时间序列的自回归模型在杨凌地区降水量预报中的应用[J].中国农学通报,2007,23(11) :403407.(DONG Xiaomeng,LUO Fengjuan ,GUO Mancai et al.The Autoregressive Model of Time Sequence Apply to Rainfall Forecast in Yangling Area[J].Chinese Agricultural Science Bulletin,2007,23(11) :403407.(in Chinese))
[4]许晓艳.基于时间序列的浑河流域降水量预测模型[J].水土保持应用技术,2011(12):2627.(XU Xiaoyan.Prediction of Precipitation Model of Hunhe River Based on Time Series[J].Technology of Soil and Water Conservation,2011,(12):2627.(in Chinese))
[5]刘杨,徐淑琴,董丽丽,等.时间序列模型在查哈阳农场降雨量预测中应用[J].黑龙江水利科技,2012,40(2):57.(LIU Yang,XU Shuqin,DONG Lili,et al.Application of Time Series Model in Chahayang Farm in Rainfall Prediction[J].Heilongjiang Science and Technology of Water Conservancy,2012,40(2):57.(in Chinese))
[6]Wang S,Feng J,Liu G.Application of Seasonal Time Series Model in the Precipitation Forecast[J].Mathematical and Computer Modelling,2013,58(34):677683.
[7]何书元.应用时间序列分析[M].北京:北京大学出版社,2004.(HE Shuyuan.Application of time series analysis[M].Beijing.Peking University press,2004.(in Chinese))
[8]SayemuzzamanM,Jha M.Seasonal and Annual Precipitation Time Series Trend Analysis in North Carolina[J].United States.Atmospheric Research,2014(137):183194.
[9]Vilibic I,Leder N.Longterm Variations in the Mediterranean Sea level Calculated by Spectral Analysis[J].Oceanographic Literature Review,1997,44(5):422.参考文献(References):
[1]汤成友,官学文,张世明.现代中长期水文预报方法及其应用[M].北京:中国水利水电出版社,2008.(TANG Chengyou,GUANG Xuewen,ZHANG Shiming.Modern Midtolongterm Hydrological Forecasting Method and its Application[M].China Water Power Press.2008.(in Chinese))
[2]Martinez J C ,Maleski J J,Miller F M.Trends in Precipitation and Temperature in Florida[J].USA J Hydrol,2012,452453,259281.
[3]董晓萌,罗凤娟,郭满才,等.时间序列的自回归模型在杨凌地区降水量预报中的应用[J].中国农学通报,2007,23(11) :403407.(DONG Xiaomeng,LUO Fengjuan ,GUO Mancai et al.The Autoregressive Model of Time Sequence Apply to Rainfall Forecast in Yangling Area[J].Chinese Agricultural Science Bulletin,2007,23(11) :403407.(in Chinese))
[4]许晓艳.基于时间序列的浑河流域降水量预测模型[J].水土保持应用技术,2011(12):2627.(XU Xiaoyan.Prediction of Precipitation Model of Hunhe River Based on Time Series[J].Technology of Soil and Water Conservation,2011,(12):2627.(in Chinese))
[5]刘杨,徐淑琴,董丽丽,等.时间序列模型在查哈阳农场降雨量预测中应用[J].黑龙江水利科技,2012,40(2):57.(LIU Yang,XU Shuqin,DONG Lili,et al.Application of Time Series Model in Chahayang Farm in Rainfall Prediction[J].Heilongjiang Science and Technology of Water Conservancy,2012,40(2):57.(in Chinese))
[6]Wang S,Feng J,Liu G.Application of Seasonal Time Series Model in the Precipitation Forecast[J].Mathematical and Computer Modelling,2013,58(34):677683.
[7]何书元.应用时间序列分析[M].北京:北京大学出版社,2004.(HE Shuyuan.Application of time series analysis[M].Beijing.Peking University press,2004.(in Chinese))
[8]SayemuzzamanM,Jha M.Seasonal and Annual Precipitation Time Series Trend Analysis in North Carolina[J].United States.Atmospheric Research,2014(137):183194.
[9]Vilibic I,Leder N.Longterm Variations in the Mediterranean Sea level Calculated by Spectral Analysis[J].Oceanographic Literature Review,1997,44(5):422.参考文献(References):
[1]汤成友,官学文,张世明.现代中长期水文预报方法及其应用[M].北京:中国水利水电出版社,2008.(TANG Chengyou,GUANG Xuewen,ZHANG Shiming.Modern Midtolongterm Hydrological Forecasting Method and its Application[M].China Water Power Press.2008.(in Chinese))
[2]Martinez J C ,Maleski J J,Miller F M.Trends in Precipitation and Temperature in Florida[J].USA J Hydrol,2012,452453,259281.
[3]董晓萌,罗凤娟,郭满才,等.时间序列的自回归模型在杨凌地区降水量预报中的应用[J].中国农学通报,2007,23(11) :403407.(DONG Xiaomeng,LUO Fengjuan ,GUO Mancai et al.The Autoregressive Model of Time Sequence Apply to Rainfall Forecast in Yangling Area[J].Chinese Agricultural Science Bulletin,2007,23(11) :403407.(in Chinese))
[4]许晓艳.基于时间序列的浑河流域降水量预测模型[J].水土保持应用技术,2011(12):2627.(XU Xiaoyan.Prediction of Precipitation Model of Hunhe River Based on Time Series[J].Technology of Soil and Water Conservation,2011,(12):2627.(in Chinese))
[5]刘杨,徐淑琴,董丽丽,等.时间序列模型在查哈阳农场降雨量预测中应用[J].黑龙江水利科技,2012,40(2):57.(LIU Yang,XU Shuqin,DONG Lili,et al.Application of Time Series Model in Chahayang Farm in Rainfall Prediction[J].Heilongjiang Science and Technology of Water Conservancy,2012,40(2):57.(in Chinese))
[6]Wang S,Feng J,Liu G.Application of Seasonal Time Series Model in the Precipitation Forecast[J].Mathematical and Computer Modelling,2013,58(34):677683.
[7]何书元.应用时间序列分析[M].北京:北京大学出版社,2004.(HE Shuyuan.Application of time series analysis[M].Beijing.Peking University press,2004.(in Chinese))
[8]SayemuzzamanM,Jha M.Seasonal and Annual Precipitation Time Series Trend Analysis in North Carolina[J].United States.Atmospheric Research,2014(137):183194.
[9]Vilibic I,Leder N.Longterm Variations in the Mediterranean Sea level Calculated by Spectral Analysis[J].Oceanographic Literature Review,1997,44(5):422.
随便看

 

科学优质学术资源、百科知识分享平台,免费提供知识科普、生活经验分享、中外学术论文、各类范文、学术文献、教学资料、学术期刊、会议、报纸、杂志、工具书等各类资源检索、在线阅读和软件app下载服务。

 

Copyright © 2004-2023 puapp.net All Rights Reserved
更新时间:2025/2/11 4:05:37