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

灰色马尔可夫模型在遂宁地区降水量中的应用陆鸿彬 游泳（遂宁气象局，四川遂宁，629000） 摘要：本文根据遂宁地区1961-2000 年的降水量资料，通过年降水量与灰色预测结论的比较进行分级建立降水序列的分级标准，采用规范化的各阶自相关系数为权重，通过状态转移概率矩阵预测未来时段的降水状态，并应用级别特征值计算具体的降水量。通过建立的马尔科夫模型对遂宁地区2001-2011 年降水量进行预测和检验，准确率较高，能够较好地满足中长期降水预报的误差要求。该方法客观、准确、可靠、简便，为遂宁区中长期降水的预测提供了新的解决途径。关键词：灰色预测；加权马尔科夫链；年降水量；预测1 引言 降水量是影响区域气候的重要气象因子，与蒸发量、地表返照率、积雪量等各种气象要素密切相关，同时也对人类生产生活活动有重要影响。在气象学、水文学中，降水量预测都是一项重要的研究工作。由于气象系统的复杂性、多样性，使得降水过程具有不确定性、较难精确预测的特点，尤其是对于较长时间段内的降水量的预测，更难以通过物理成因的分析来准确预测该时段内的降水量值。鉴于此，本文引入了灰色-马尔科夫链模型，以遂宁地区1961-2010 年的全年降水量资料为例，探讨了灰色-马尔科夫链模型在遂宁地区年降水量预测中的应用。2 、物理机制2.1 GM( 1 ,1) 模型建模机理 灰色系统是既含有已知的、又含有未知的或非确定的信息系统。对灰色系统进行预测的方法即为灰色预测法。灰色系统理论是从信息的非完备性出发研究和处理复杂系统的理论，它不是从系统内部特殊的规律出发去研究系统，而是通过对系统某一层次的观测资料加以数学处理，达到在更高层次上了解系统内部变化趋势、相互关系等机制的。灰色预测通过鉴别系统因素之间发展趋势的相依程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程模型，从而预测事物未来发展的趋势。2.2 马尔科夫链 马尔科夫过程一种典型的随机过程。该过程是研究一个系统的状况及其转移的理论。它是通过对不同状态的初始概率以及状态之间的转移概率的研究，来确定状态的变化趋势，从而达到对未来进行预测的目的。马尔科夫过程有两个基本特征：一是“无后效性”，即事物将来的状态及其出现的概率的大小，只取决于该事物现在所处的状态，而与以前时间的状态无关；二是“遍历性”，是指不管事物现在出于什么状态，在较长时间内，马尔科夫过程逐渐趋于稳定状况，而且与初始状况无关。用数学语言描述马尔科夫过程就是：设为随机过程，若在时刻对观测得到相应的观测值满足条件或 则称此类过程为具有马尔科夫性质的过程或马尔科夫过程，简称马氏过程。其中代表在的条件下时刻取值得条件分布函数。若把时刻看成“现在”，因为 则就可以看成“将来”，就当做“过去”。因此上述定义可表述为现在的状态取值为的条件下，将来状态的取值于过去状态的取值是无关的。2.3 加权马尔科夫链预测模型的方法和步骤计算年降水量误差序列的各阶自相关系数rk表示第k 阶自相关系数， 表示第t 年的年降水量， 表示近51 年北京市年降水量的平均值，n 为年降水量序列的长度。规范化各阶自相关系数，即：wk为各种滞时的马尔科夫链的权。以前面的m个时段的状态为初始状态，转移步数为1，2，m。结合相应的状态概率矩阵所对应的行向量，组成一个新的概率矩阵。将同一个状态的m项预测概率与所对应的加权相乘之后求和，即：max( pi，i E) 所对应的状态即为该时段指标值的预测状态。3、实例分析遂宁市地处四川盆地中部丘陵区，位于涪江中游，是典型的农业区，属四川盆地亚热带湿润季风气候区，气候温和，季风气候显著，具有冬暖春旱、夏热秋凉的特征。雨量较为充沛，但分布不均，年季变幅大； 410月为雨季，降雨量占全年雨量的89%左右；盛夏（78月）易产生局地对流性天气，降雨最为集中，常有暴雨洪涝灾害；因其地形地势、气候变化和区域环流等的共同作用，也常造成春旱、夏旱、伏旱和冬干发生，其干旱具有种类齐全、发生频繁、危害面积广、损失程度重、季节性强等特点。由于降水时间分布极为不均，极端事件频发，给社会经济发展和粮食农作物的生产造成严重损失，因此深入分析遂宁地区近几十年来的降水特征、变化规律，对该地区短期气候预测和保障农业持续稳定增长具有重要意义。3.1 灰色预测选取1961-2011年遂宁市年降水量的数据作为预测的原始时间序列，记为，其中表示第年的遂宁年降水量，。表1 19612010年遂宁市年降水量 （单位：毫米）年份原始值年份原始值年份原始值年份原始值年份原始值1961891.51971753.719811191.7199110062001877.819621092.91972771.419821146.71992975.120021097.51963935.319731186.71983889.419931215.120031072.11964925.619741096.419841064.91994722.12004837.419651313.61975964.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为了消除序列的随机性，对原始序列进行一次累加，得到累加生成列，如图1所示。从图中能明显看出序列的平滑的指数增长形式。 图1 年降水量累加为此，利用上述灰色预测法，作序列的预测模型。计算过程通过在MATLAB软件中进行。 图2年降水量实际值与灰色预测值对比从图2中可以看出，灰色预测的结果恰当的反映了遂宁年降水量的长期趋势。在进行模型的残差检验时，模型的相对残差序列中有较大的数，导致相对误差较大的原因是由于指标的上下波动，而不是预测值与实际值曲线的严重偏离，基于此，不采用灰色预测法进行模型修正，因为灰色预测得到的残差已经不是一组符号相同的数，根据灰色预测的原理，这样的修正已经没有什么意义了。因此本文根据年降水量波动的特点，采用马尔科夫法进行模型修正。 3.2利用马尔科夫状态转移矩阵修正遂宁年降水量预测值从图2中可以获得，灰色模型预测得到的预测结果呈现递减的指数分布，恰当的反映了实际值的长期趋势，但没有反映出实际值本应出现的波动情况。此问题可以通过对灰色预测模型的结果进行马尔科夫链改进来提高其预测的准确性。马尔科夫链过程是以其发现者俄国著名的数学家马尔科夫命名的，其目的是根据某些变量的现在状态及其变化趋向，来预测它在未来某一特定期间可能出现的状态，从而提供某种决策的依据。 在本文中，将利用长度为40年的降水时间序列对年降水量进行预测，并与实际观测值进行对比，评估加权马尔科夫链模型对遂宁地区降水量的预测效果。首先，本文以19612000年遂宁年降水量序列来预测2001年的年降水量，并与2001年的年降水量实际观测值进行对比。经计算，遂宁19612000 年年降水量的均值为932.6mm，均方差为166.3mm，为了获得状态转移矩阵，首先需要确定遂宁年降水量的状态。根据马尔科夫链分析方法的应用经验和实际情况，按照年降水量与灰色预测结论的比较，可以划分为5种状态，详细情况列表如表2：表2 遂宁年降水量状态划分状态名称状态特点状态区间包括年份年份总数旱1966 1971 1976 1994 1995 1996 1997 20068年偏旱1961 1967 1969 1972 1978 19831986 1990 20049年正常1963 1964 1968 1970 1975 1977 1987 1988 1989 1999 2000 2001 2007 2008 201115年偏涝1962 1974 1979 1980 1984 1991 1992 1998 2002 2003 2005 200912年涝1965 1973 1981 1982 1985 1993 20107年注：表中表示实际年降水量，表示灰色预测值，表示实际年降水量均值。 在进行降水量的预测时，主要考虑利用步长为15 的概率转移矩阵进行计算。转移矩阵P(k) 的第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。据1996-2000年降水量及其相应的状态转移概率矩阵对2001年的降水量状态及其降水量进预测。计算结果见表3。 表3 2001年加权马尔科夫预测结果初始年滞时权重12345200010.53110.09660.14480.193100.0483199920.10450.00950.01900.02850.01900.0095199830.03920.00980.00490.009800.0098199740.2871000.16400.08200199650.03810.010900.010900.0054 pi 加权求和0.12680.16870.40640.10100.0730由表3可知：当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年遂宁年降水量预测和检验 表4 20012011 年降水量预测值与实际观测值对比表（单位：mm）年份马尔科夫预测值实况值误差(%)2001年877.5877.8-0.042002年910.11097.5-17.12003年966.91072.1-9.82004年971.6837.4162005年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 年是遂宁出现罕见高温伏旱天气，但仍预报出偏旱状况。预测起始年份状态以及各类状态间转移的概率，若待预测年份降水量为大旱或大涝，则模型预测将可能会有较大误差。但同时由于马尔科夫链具有很强的“自我调整能力”，故其中某年的预测误差偏大并不会影响到后续年份的预测。4 结语 降水量的预测在气候预测中占有重要地位，但是中长期的降水量预测仍然是当前气候预测中的难题。本文利用灰色-马尔科夫链模型和遂宁地区50 年的年降水量资料，对20012011 十一年的年降水量进行了预测并与实际观测值进行对比，可得到以下结论：1）灰色模型预测得到的预测结果呈现递减的指数分布，恰当的反映了实际值的长期趋势，但没有反映出实际值本应出现的波动情况。2） 利用年降水量与灰色预测结论的比较分级法对年降水量进行分级，并确定各年降水量的状态，较好把遂宁年降水量分为5级；3）马尔科夫链模型预测以各种步长的概率转移矩阵为依据，用各种步长的马尔可夫链加权来预测降水量状态； 4）马尔科夫链模型预测能较好的预测年降水量，以遂宁11年时间长度的降水量预测为例，预测误差在15%以内的年份占81.8%；5） 对于大旱或大涝年份降水量的预测，马尔科夫链模型存在着预测误差偏大的缺陷，需要在后续研究中进一步调整模型预测方法，以期得到更好的预测效果。参考文献： 1 张曙红等.灰色马尔科夫SCGM(1,1)预测模型J.佛山科学技术学院学报：自然科学版，2004，3（1）:16-19. 2 孙才志,张