灰色-马尔可夫模型在遂宁地区降水量中的应用1

1、灰色马尔可夫模型在遂宁地区降水量中的应用陆鸿彬游泳（遂宁气象局，四川遂宁，629000 ）摘要：本文根据遂宁地区1961-2000 年的降水量资料，通过年降水量与灰色预测结论的比较进行分级建立降水序列的分级标准， 采用规范化的各阶自相关系数为权重， 通过状态转移概率矩阵预测未来时段的降水状态， 并应用级别特征值计算具体的降水量。 通过建立的马尔科夫模型对遂宁地区 2001-2011 年降水量进行预测和检验， 准确率较高， 能够较好地满足中长期降水预报的误差要求。该方法客观、准确、可靠、简便，为遂宁区中长期降水的预测提供了新的解决途径。关键词：灰色预测；加权马尔科夫链；年降水量；预测引言降水量是

2、影响区域气候的重要气象因子，与蒸发量、地表返照率、积雪量等各种气象要素密切相关，同时也对人类生产生活活动有重要影响。在气象学、水文学中， 降水量预测都是一项重要的研究工作。由于气象系统的复杂性、多样性，使得降水过程具有不确定性、较难精确预测的特点，尤其是对于较长时间段内的降水量的预测，更难以通过物理成因的分析来准确预测该时段内的降水量值。鉴于此，本文引入了灰色- 马尔科夫链模型，以遂宁地区1961-2010年的全年降水量资料为例，探讨了灰色- 马尔科夫链模型在遂宁地区年降水量预测中的应用。2 、物理机制2.1 GM( 1 ,1)模型建模机理灰色系统是既含有已知的、又含有未知的或非确定的信息系统

3、。对灰色系统进行预测的方法即为灰色预测法。灰色系统理论是从信息的非完备性出发研究和处理复杂系统的理论，它不是从系统内部特殊的规律出发去研究系统，而是通过对系统某一层次的观测资料加以数学处理， 达到在更高层次上了解系统内部变化趋势、相互关系等机制的。灰色预测通过鉴别系统因素之间发展趋势的相依程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律， 生成有较强规律性的数据序列，然后建立相应的微分方程模型，从而预测事物未来发展的趋势。2.2马尔科夫链马尔科夫过程一种典型的随机过程。该过程是研究一个系统的状况及其转移的理论。它是通过对不同状态的初始概率以及状态之间的转移概率的研究，来确定状

4、态的变化趋势，从而达到对未来进行预测的目的。马尔科夫过程有两个基本特征：一是 “无后效性”， 即事物将来的状态及其出现的概率的大小，只取决于该事物现在所处的状态，而与以前时间的状态无关；二是“遍历性”，是指不管事物现在出于什么状态，在较长时间内， 马尔科夫过程逐渐趋于稳定状况，而且与初始状况无关。用数学语言描述马尔科夫过程就是：设 X ( t), tT为随机过程，若在t1, t2 , tn 1, tn (t1t2tn 1tnT )时刻对X (t)观测得到相应的观测值x1 , x2 , xn1 , xn满足条件P X ( t n )x n X ( t n 1 )x n 1 , X ( t n 2

5、 )x n 2 , X ( t 1 )x 1P X ( t n )x n X ( t n 1 )x n 1或F Xxn ; t n x n 1 , x n 2 , x 2 , x1 ; t n 1 , t n 2 , t 2 , t1F Xx n ; t n x n 1 ; t n 1则称此类过程为具有马尔科夫性质的过程或马尔科夫过程，简称马氏过程。其中FX xn ;tn xn 1, xn 2 , x2, x1 ;tn 1 ,t n 2, t2, t1代表在 X (tn 1 )xn1, X (tn 2 ) xn 2 , X ( t1 )x1 的条件下时刻X (tn ) 取 xn 值得条件分布函

6、数。若把 tn 1时刻看成“现在” ，因为 t1t2tn 1 tn则 tn 就可以看成“将来”，t1, t2 , , tn 2 就当做“过去” 。因此上述定义可表述为现在的状态X (tn 1 ) 取值为 xn 1 的条件下，将来状态X (t) 的取值于过去状态X (t ), X (t), X ( tn2) 的取值是无关的。n122.3加权马尔科夫链预测模型的方法和步骤计算年降水量误差序列的各阶自相关系数n kn k( xt x)2rk(xt x)( xt k x) /t 1t 1rk 表示第 k 阶自相关系数，表示第 t年的年降水量，表示近 51年北京市年降水量的平均值， n 为年降水量序列的

7、长度。规范化各阶自相关系数，即：wk 为各种滞时的马尔科夫链的权。mwkrk/rk(m5)k 1以前面的 m个时段的状态为初始状态，转移步数为 1，2，?m。结合相应的状态概率矩阵所对应的行向量，组成一个新的概率矩阵。将同一个状态的m项预测概率与所对应的加权相乘之m后求和，即： piwk pi (k )k 1max( pi ， i E)所对应的状态即为该时段指标值的预测状态。3、实例分析遂宁市地处四川盆地中部丘陵区，位于涪江中游，是典型的农业区，属四川盆地亚热带湿润季风气候区，气候温和，季风气候显著，具有冬暖春旱、夏热秋凉的特征。雨量较为充沛，但分布不均，年季变幅大； 410 月为雨季，降雨量

8、占全年雨量的 89%左右；盛夏（ 78 月）易产生局地对流性天气，降雨最为集中，常有暴雨洪涝灾害；因其地形地势、气候变化和区域环流等的共同作用，也常造成春旱、夏旱、伏旱和冬干发生，其干旱具有种类齐全、发生频繁、危害面积广、损失程度重、季节性强等特点。由于降水时间分布极为不均，极端事件频发， 给社会经济发展和粮食农作物的生产造成严重损失，因此深入分析遂宁地区近几十年来的降水特征、变化规律， 对该地区短期气候预测和保障农业持续稳定增长具有重要意义。3.1灰色预测选 取1961-2011 年遂宁市年降水量的数据作为预测的原始时间序列，记为X (0)( x(0) (1),x(2)(2), , x(0)

9、 (51) ，其中x(0) (k) 表示第 1960 k 年的遂宁年降水量，k 1,2,51 。表 1 19612010 年遂宁市年降水量（单位：毫米）年份原始值年份原始值年份原始值年份原始值年份原始值1961891.51971753.719811191.7199110062001877.819621092.91972771.419821146.71992975.120021097.51963935.319731186.71983889.419931215.120031072.11964925.619741096.419841064.91994722.12004837.419651313.61

10、975964.519851195.91995602.920051058.91966741.51976735.31986833.81996648.52006549.81967790.91977907.11987929.81997610.42007933.71968909.51978832.21988954.919981037.820088591969857.919791013.61989913.21999892.320099901970959.91980997.61980866.72000871.920101146.9为了消除序列的随机性，对原始序列进行一次累加，得到累加生成列X (1) ，如图

11、 1所示。从图中能明显看出序列X (1) 的平滑的指数增长形式。图 1年降水量累加为此，利用上述灰色预测法，作序列X (1) 的预测模型。计算过程通过在MATLAB软件中进行。图 2 年降水量实际值与灰色预测值对比从图2 中可以看出，灰色预测的结果恰当的反映了遂宁年降水量的长期趋势。在进行模型的残差检验时， 模型的相对残差序列中有较大的数，导致相对误差较大的原因是由于指标的上下波动， 而不是预测值与实际值曲线的严重偏离，基于此， 不采用灰色预测法进行模型修正， 因为灰色预测得到的残差已经不是一组符号相同的数，根据灰色预测的原理， 这样的修正已经没有什么意义了。因此本文根据年降水量波动的特点，采

12、用马尔科夫法进行模型修正。3.2 利用马尔科夫状态转移矩阵修正遂宁年降水量预测值从图 2中可以获得，灰色模型预测得到的预测结果呈现递减的指数分布，恰当的反映了实际值的长期趋势， 但没有反映出实际值本应出现的波动情况。此问题可以通过对灰色预测模型的结果进行马尔科夫链改进来提高其预测的准确性。马尔科夫链过程是以其发现者俄国著名的数学家马尔科夫命名的，其目的是根据某些变量的现在状态及其变化趋向，来预测它在未来某一特定期间可能出现的状态，从而提供某种决策的依据。在本文中，将利用长度为40年的降水时间序列对年降水量进行预测，并与实际观测值进行对比，评估加权马尔科夫链模型对遂宁地区降水量的预测效果。首先，

13、本文以19612000年遂宁年降水量序列来预测2001年的年降水量，并与 2001年的年降水量实际观测值进行对比。经计算，遂宁 19612000年年降水量的均值为932.6mm，均方差为 166.3mm，为了获得状态转移矩阵， 首先需要确定遂宁年降水量的状态。根据马尔科夫链分析方法的应用经验和实际情况，按照年降水量与灰色预测结论的比较，可以划分为5种状态，详细情况列表如表 2：表2遂宁年降水量状态划分状态状态特点状态区间包括年份年份总数名称X?196619711976F 1旱X1994199519968年20%X1997 2006?196119671969F 2偏旱20%X X5%1972 1

14、978 19839年X1986 1990 2004F3正常F4偏涝?XX5%5%X?XX5%20% X1963 196419681970 197519771987 1988198915年1999 200020012007 2008201119621974197919801984199112年19921998200220032005 2009196519731981F 5涝?XXX20%19821985 19937年2010注：表中?X 表示实际年降水量均值。X 表示实际年降水量， X 表示灰色预测值，在进行降水量的预测时，主要考虑利用步长为15 的概率转移矩阵进行计算。转移矩阵 P(k)的第

15、i 行第 j 列元素表示由状态 i 经 k 步转移至状态 j 的概率。统计得到步长为15 的各概率转移矩阵。3.3对 2001 年遂宁年降水量预测和检验19962000 年分别经过 51 步（年） 状态转移即到达 2001 年，为了求得 2001 年的年降水量预测值，可根据 19962000 年降水量所处状态。根据公式，可以得到自相关系数， r1=0.1962 ，r2=-0.0386 ，r3=-0.0145. ，r4=-0.106. ，r5=0.0141 ，从而得到权重， w1=0.5311 ，w2=0.1045， w3=0.0392 ， w4=0.2871 ，w5=0.0381 。据 199

16、6-2000 年降水量及其相应的状态转移概率矩阵对 2001 年的降水量状态及其降水量进预测。计算结果见表3。表 32001年加权马尔科夫预测结果初始年滞时权重12345200010.53110.09660.14480.193100.0483199920.10450.00950.01900.02850.01900.0095199830.03920.00980.00490.009800.0098199740.2871000.16400.08200199650.03810.010900.010900.0054pi加权求和0.12680.16870.40640.10100.0730由表 3可知：当

17、i=3 时， pi=0.4064 为最大值，这说明2001年的降水量状态为 3（正常），即降水量区间为 868.4,959.8，而 2001年实况降水量为877.8mm。当 =1.15 时，级别特征值为2.901 ，可以求出预测降水量为877.5mm，误差为 0.04%。3.3 2001年到 2011年遂宁年降水量预测和检验表420012011年降水量预测值与实际观测值对比表（单位：mm）年份马尔科夫预测值实况值误差 (%)2001年877.5877.8-0.042002 年910.11097.5-17.12003 年966.91072.1-9.82004 年971.6837.4162005

18、年971.91058.9-8.22006 年851.2549.854.82007 年844.3933.7-9.62008 年888.98593.52009 年889.1990-102010 年983.01146.9-14.32011 年833.9892.4-6.6分析表 4中的数据可知，若将预测值与观测值的误差控制在15%以内的预测认为是“有效预测” ，则对于遂宁地区20012011 年降水量的预测的有效率为81.8%，其中预测误差在10%以内的年份占 63.6%。预测误差较大的3 年分别为 2002、 2004 和 2006 年，尤其是误差较大的2006 年是遂宁出现罕见高温伏旱天气，但仍预

19、报出偏旱状况。预测起始年份状态以及各类状态间转移的概率， 若待预测年份降水量为大旱或大涝， 则模型预测将可能会有较大误差。但同时由于马尔科夫链具有很强的 “自我调整能力” ，故其中某年的预测误差偏大并不会影响到后续年份的预测。4 结语降水量的预测在气候预测中占有重要地位，但是中长期的降水量预测仍然是当前气候预测中的难题。 本文利用灰色 - 马尔科夫链模型和遂宁地区50 年的年降水量资料， 对 20012011十一年的年降水量进行了预测并与实际观测值进行对比，可得到以下结论：1）灰色模型预测得到的预测结果呈现递减的指数分布，恰当的反映了实际值的长期趋势，但没有反映出实际值本应出现的波动情况。2）

20、 利用年降水量与灰色预测结论的比较分级法对年降水量进行分级， 并确定各年降水量的状态，较好把遂宁年降水量分为 5级；3）马尔科夫链模型预测以各种步长的概率转移矩阵为依据，用各种步长的马尔可夫链加权来预测降水量状态；4 ）马尔科夫链模型预测能较好的预测年降水量，以遂宁11年时间长度的降水量预测为例，预测误差在 15%以内的年份占 81.8%；5） 对于大旱或大涝年份降水量的预测，马尔科夫链模型存在着预测误差偏大的缺陷，需要在后续研究中进一步调整模型预测方法，以期得到更好的预测效果。参考文献： 张曙红等 . 灰色马尔科夫 SCGM(1,1)预测模型 J. 佛山科学技术学院学报：自然科学版， 200

21、4，3（ 1） :16-19.孙才志 , 张戈 , 林学钰 . 加权马尔可夫模型在降水丰枯状况预测中的应用J. 系统工程理论与实践 , 2003, (4): 100-105.马占青 , 徐明仙 , 俞卫阳等 . 年降水量统计马尔科夫预测模型及其应用 J. 自然资源学报, 2010, 25(6): 1033-1041.5权马尔科夫链在北京市年降水量预科协论坛 2010年第 1期（下）作者简介陆鸿彬出生1980年3月 29日现在遂宁市气象局从事天气预报与服务工作，气象台副台长、 工程师，地址遂宁市明月路557号，邮编 629000，电话， E-mail ：出师表两汉：诸葛亮先帝创业未半而中道崩殂，今天下三分， 益州疲弊， 此诚危急存亡之秋也。然侍卫之臣不懈于内，忠志之士忘身于外者