水利工程论文-月径流序列的多层递阶预报研究.doc

水利工程论文-月径流序列的多层递阶预报研究摘要：月径流序列是一类具有周期变化的非平稳时间序列.本文根据其特点，建立了多层递阶预报模型，文中对此类非平稳时间序列的建模及预报方法进行了深入研究.关键词：月径流序列多层递阶预报时变参数非平稳序列ResearchontheMulti-LayerHierarchicalForecastofMonthlyRunoffSeriesAbstractMonthlyrunoffseriesisasortofnon-stationarytimesequencewithperiodicity.Inthispaper,accordingtotheircharacteristics,amulti-layerhierarchicalforecastmodelisestablished,andmodellingandforecastingmethodsarethoroughlyresearched.Keywordsmonthlyrunoffseries；multi-layerhierarchicalforecast；time-varyingparametric；non-stationaryseries1多层递阶预报的基本原理1.1后验残差一致小准则令Yk-1=y(0),y(1),y(k-1)，Uk=u(0),u(1),u(k)为系统输入；l(k)为随机噪声；(k)是m维的参数向量，y(k)是系统输出，则(1)系统输出向前一步的预报公式是(2)一般的自适应算法仅能根据Yk-1,Uk得到(k-1)的估计值(k-1)，若系统是快时变的，(k)与(k-1)将有较大差别.为了减少预报误差，应设法求出(k)的估计值*(k).此时，(3)上式称为含有参数预报值的预报公式.相应的预报残差为(4)式中(k*(k)称为后验残差.可以证明1，对于任何的0，在适当条件下必有N0和0存在，使得当kN时，由算法所确定的(k)(5)恒满足称此为后验残差一致小准则.若则(6)多层递阶预报方法，就是采用式(5)或式(6)的算法，对时变参数进行预测，使估值序列(k)满足后验残差一致小准则.1.2多层递阶建模与预报多层递阶预报，是将动态系统预报分为时变参数预报与系统状态预报两部分.其中参数预报包括时变参数跟踪估计以及对参数未来值进行预测.1)时变参数跟踪估计时变系统不仅状态在变化，而且参数也在变化，多层递阶预报充分考虑了系统的这种时变特性.由掌握的历史资料，按式(5)或式(6)进行参数跟踪估计，可得参数(k)的一系列估计值：(1),(2)，,(N)，N为资料长度.2)时变参数预测对参数估计序列(k),k=1,2,N，进行分析，寻找其变化规律，采用适当的数学方法建立时变参数的预报公式，进而可求得时变参数一系列的预测值*(N+1)，*(N+2)，,*(N+h)，h为预报步长.常采用的时变参数预测方法有：定常增量法、常数因子法、分段周期法、周期增量法、多层自回归法等2.3)系统状态预报系统状态向前多步的预报公式为(7)式(7)中的*(N+h)是用各种方法求得的时变参数预测值.这种方法由于充分考虑了模型参数的时变性，使预报精度有所提高.2月径流序列的多层递阶预报月径流序列是季节性的非平稳时间序列，具有更强的随机性和非平稳性，对此时间序列进行中长期预报有一定难度.2.1数据预处理与建模对月径流序列的季节性变化所引起的非平稳性，可进行中心化、标准化、季节性差分处理.本文在对淮河某大型水库月径流序列预报时，进行了季节性差分：式中：Xt为月径流序列；Xt=Xt-Xt-1；12为月度差分算子，12Xt=Xt-Xt-12.现以新序列ut作为多层递阶预报的建模序列.经模型识别、参数估计、模型检验，得ut序列的二阶自回归模型AR(2)为：ut=a1(t)ut-1+a2(t)ut-2+t(8)式中a1(t),a2(t)为模型参数，已采用最小二乘法求得.2.2时变参数估计对于ut的AR(2)模型可写成y(k)=(k)T(k)+e(k)(9)(k)=u(k-1),u(k-2)T(k)=a1(k),a2(k)根据后验残差一致小准则，时变参数(k)的估值算法可按式(6)写成如下递推形式，称为(k)的跟踪公式.(10)按分量展开，即i=1,2；k=2,3,，k,N(11)式中：(1)和2(1)可采用原参数估计值.按上述跟踪公式计算，可得1(2),1(3),1(N)；2(2),2(3),2(N)的估计值.2.3时变参数预测根据月径流序列的特点，现采用如下方法建立时变参数预测公式.1)分段周期法.月径流序列存在以年为周期的季节性变化，可按下式推求各月时变参数的预测值，即(12)式中：为第i个预报参数进行h步预报的参数预测值；h为预报步长，h=1,2,，12；N为样本长度，M为样本资料年限.2)多层自回归方法.对参数估值序列i(k)建立自回归模型AR(2)：若采用跟踪公式求得的变化不大，可用其平均值值，进而可求得参数估计值*i(k)，(13)(14)若变化较大，可建立多层自回归模型，然后预报出3)周期增量法.对于具有周期增量的序列，其参数估计公式为其中:p满足h=pT+q,qT.若p=0,T=1,h=q时，(15)式中：T为正整数；N为预报的起始时刻；h为预报步长；i(j)为时变参数的估值增量.3应用实例淮河某大型水库的月径流预报，采用了多层递阶预报方法.在实际应用中进行了多种方法的对比研究，并针对月径流序列的特点，对模型及时变参数估计和预测进行了改进.多层递阶预报方法的关键，在于时变参数的估计和预测.为了探讨月径流序列预报模型参数的时变规律，对时变参数估计采用了跟踪公式法和直接参数估计法；时变参数预测采用了多层自回归方法、周期增量法、分段周期法.鉴于月径流序列的模型参数时变性较强，不可直接采用时段平均值作为参数的预报值.多层自回归方法中的高层自回归模型参数改用了参数估计方法；周期增量法考虑了相邻时段时变参数增量的变化，对外推预报时各时段增量i(j)的计算进行了改进.根据水库运行调度的需要，编制了两种预报调度方式的计算机程序：年度预报和逐月预报.年度预报是在年底预报出下一年12个月的月径流量，为编制下年度的水库调度计划提供依据.逐月预报是在年度预报的基础上，根据前期月径流的实测值，逐月递推预报出下一个月的月径流量，是实施水库优化调度的基础.4模型评定与检验4.1模型评定模型评定是用建模序列进行模拟预报，根据预报误差评定模型的有效性.设是第i年第j月的h步预报值，xij为相应时段的实测值；MSE(j,h)为第j月预报值的均方误差；S2j为实测序列xij第j月的方差，则确定性系数R2(j,h)为(16)显然，确定性系数R2愈接近于1，模型精度愈高.表1为淮河某大型水库月径流序列多层递阶预报的模型评定表.表1淮河某水库月径流序列多层递阶预报评定表月份123456789101112R20.960.970.960.900.940.820.770.760.920.930.980.96等级甲甲甲甲甲乙乙乙甲甲甲甲由表1可见，汛期68月的R2在0.700.89间，属乙等，其余月份的R2值均在0.90以上，属甲等.说明模型评定精度尚好，所建立的预报模型是有效的.4.2模型检验模型检验是用未参加建模的预留序列进行外推预报，以检验模型的预报精度和可靠性.表2是淮河某大型水库月径流序列多层递阶预报模型的预报值与实测值对照表.表2淮河某水库多层递阶预报模型检验表(106m3)月份项目年份1234567891011121989实测值9.6311.340.433.379.9160.079.7141.125.710.110.611.0预报值10.821.438.656.472.4106.4120.783.617.414.715.010.61990实测值11.0389.061.438.353.047.764.633.511.19.029.92.96预报值10.334.544.142.469.868.4138.963.629.120.519.88.3由于径流变化的随机性，月径流预报作到完全准确是困难的，特别是水文序列出现丰枯交替或年内分配异常时.尽管该模型对月径流序列进行了季节性差分处理，利用了以年为采样周期和以月为采样周期的相关信息，并对时变参数进行了预测，使月径流预报精度有所提高，但仍然未能完全克服统计模型的弱点.就该水库的月径流序列而言，