安徽气象新安江流域旱涝趋势预测.doc

基于谱分析和ARIMA时间序列的新安江流域旱涝趋势预测 刘裕禄 凌来寿 （黄山市气象台，245021 歙县气象局,245200) 提 要 本文用1941～2000年资料详细研究了新安江流域旱、涝趋势周期变化特征，然后建立相应的预报模型，并以2001年～2005年作预报检验。主要结果：用谐波分析新安江流域年降水表现出2.6年和20年两个振动周期,小波分析经过滤波降噪再进行重建后的降水只呈现出20年振动周期，原降水序列经过差分标准化处理后的ARIMA时间序列模型，对2001年～2005年预报结果表明，ARIMA模型效果比前两者更好，更接近实际，三种模型都预报出01年～05年新安江流域平均降水处偏少期。 关键词：降水预测、谐波 小波分析 ARIMA、时间序列 引言 新安江流域主要位于黄山市境内，中低纬度的过渡地带，洪涝与干旱是主要的自然灾害。在全球变暖的大背景下，气候变化异常，洪涝与干旱频率周期发生了一定的变化；旱涝趋势特点又如何？本文利用谱分析和时间序列模型对新安江的不同尺度旱涝周期作分析研究并得出数学模型，做出5年趋势定量预报以便检验。资料选取1941年～2000年新安江流域60年降水资料，由新安江电厂提供，新安江流域雨量站包括祁门、黟县、休宁、屯溪、黄山、歙县、千岛湖、淳安等组成，新安江年降水量是求其各站平均值所得。 1 谐波分析预测 某一气象要素的时间序列，都包含有各种复杂的波动。但是不管什么样的复杂波动曲线，总可以分解为由不同的正弦波所构成。通常这些正弦波有各种周期波动，且相互正交，最长周期等于序列长度为基波，其余正弦波为谐波，时间序列谐波分析是把不同频率的振动，按其方差贡献大小进行分解，从而研究在序列中哪一种频率振动占优势。根据数学方法，一般谐波分析模型为： X(t)= (1) 各谐波方差贡献： （2) 系数ak,bk表达式为： （3） 其中：t=0,1,2,3,……n 本文选取新安江流域1941年～2000年连续60年降水资料，时间序列样本n=60。K=60/2=30，首先利用3式计算各谐波系数ak,和bk，然后计算1到30波数各谐波方差，结果如下： 3005 62.042 11138 7549.3 472.45 1868 4534.5 443.25 4323.4 2846.6 2803.2 1829.3 301.68 167.95 1816.3 3969.6 6553.1 3385.2 8581.2 4604.6 82.753 163.77 12078 1755.3 2116.5 486.78 379.23 233.33 5140.6 5164.7 从以上可看出当k=23,3时，即周期T=n/k=2.6,20年，方差贡献都较大，因此可以说新安江流域近60年降水有20,2.6年两个主要周期。作显著性检验，结果k=3时F=3.663.5，k=23时F=4.0133.5都通过信度=0.05显著性检验，证明新安江流域近60年降水呈20年和2.6年的周期振动。将此时ak,bk 代入（1）式得降水预报回归方程： 2 小波分析滤波预测 以上利用谐波周期分析是以2.6年和20年周期振动拟合得出的降水回归方程，为了更精细说明新安江年降水振动周期，小波分析比传统的谐波周期分析方法更加优越，小波技术可以显示出传统频域分析方法一般难以实现的特性，因为小波系数能同时反映时域参数ｔ和频域参数ａ的特性。小波变换后得到的小波系数过零点可能是气候时间序列上的突变点利用1941年～200年0 年降水资料对它作连续小波变换，取Morlet小波作为母小波函数 a=17 10年 11年 a=20 10年 10年 5个过0点的突变点，可能就是旱涝趋势转换点，以a=14年为例，分别是42年、50年、59年、69年、80年、92年，平均每隔10年为一旱涝期，也就是说新安江流域降水有20年旱涝振动周期，42—50年平均降水偏少，51—59年平均降水偏多，60—69年平均降水偏少，70—80年平均降水偏多，81—92年平均降水偏少，92—2000年平均降水偏多，与实际也相符合。2000年以后进入降水偏少期。 为了进一步说明降水旱涝周期振动，以下再进行具体分析。如果对降水资料进行小波变换技术的信号分解和重建两次处理，也就是滤波降噪和反演重新得出降水资料，再进行谐波周期分析，预报未来的旱涝趋势，可