中长期水文预报方法

1、水文预报大作业水文预报大作业中长期水文预报方法作 者 姓 名： * 学科、 专业： 水文预报理论与方法 学 号： * 指 导 教 师： * 完 成 日 期： * 摘 要本文运用历史演变法、单要素模糊推理法、平稳时间序列法三种中长期预报方法，分别对某水库19672011年7月份平均流量进行处理，进而预报2012年7月份平均流量。三种方法的预报结果分别为：333.15m/s、180 m/s、 258.6 m/s。文中具体介绍了各方法的计算过程，并在文末对各方法的优缺点进行了评价。关键词：历史演变法；单要素模糊推理法；平稳时间序列法II目 录摘 要I前言- 1 -1预报方法原理- 1 -1.1历史演

2、变法- 1 -1.2单要素模糊推理法- 2 -1.3 平稳时间序列法- 2 -2三种预报方法的实例应用52.1 历史演变法52.2 单要素模糊推理法72.3 平稳时间序列法123.方法结果对比与分析14- III -前言根据前期水文气象要素，用成因分析与数理统计的方法，对未来长时期(一旬至一年)的水文要素进行科学的预测称为中长期水文预报。中长期水文预报通常泛指预见期超过流域最大汇流时间，且在3天以上，1年以内的水文预报。其预报内容主要有：径流量、河口的水位、冰清、河道及冲於变化和旱涝趋势等。本文主要运用历史演变法、单要素模糊推理法及平稳时间序列法对某水库1967年2011年7月份平均流量进行分

3、析，并预测2012年7月份的平均流量。1预报方法原理1.1 历史演变法（1） 基本思想任一水文气象要素的长期历史实测值，全面地反映了这一要素的时序变化规律。尽管目前还不能辨认影响要素变化的内因与外因或影响程度，但是诸因素的综合影响却已全面地融汇于要素的时序变化之中，若能找出该要素历史变化规律。便可用它进行预报（2）基本思路任一水文气象要素的长期历史实测值，全面地反映了这一要素的时序变化规律。尽管目前还不能辨认影响要素变化的内因与外因或影响程度，但是诸因素的综合影响却已全面地融汇于要素的时序变化之中，若能找出该要素历史变化规律，便可用它进行预报。（3）基本出发点任何一个水文要素的长期记录反映了这

4、一要素全面的变化过程。只要我们能够找出这一要素的演变规律，就可利用这些规律来进行预报。（4）基本规律持续性：连续上升或下降，升降超出平均值某一数值 后将改变原趋势的特性。 相似性：两个水文时期水文要素变化趋势的相似性。 周期性：丰、平、枯周期变化。 最大最小可能性：指水文气象要素经常出现的范围或年上升、下降的极限幅度。 转折点：两个时期之间的转折期为转折点。1.2 单要素模糊推理法单要素模糊推理建立的基本思想为：（1）统计待预报的水文要素x（年平均流量或年平均降雨量等）的历史资料，取其最大、最小值域（或适当外延）作为论域U。（2）将论域U上的x值离散化，并分成l级（如l=5级，则表示丰、偏丰、

5、正常、偏枯、枯），每级便对应一个模糊子集，Yk U，k=1，2，3，l。模糊子集的分布可取用多种形式，而取梯形分布较简单，隶属函数为Yk（x）。（3）如果有n年历史演变资料，则每个值都是论域U上的一个模糊子集，即Yi，k U，i=1，2，3，l。（4）设想利用前期连续m年模糊信息（亦称预报因子）作输入，以模糊推理法输出（亦称预报量）第m+1年模糊量，便可根据n年资料归纳出（n-m）条m次复合模糊条件语句，即:（1.1）式（1）就是具有（n-m）重m次复合模糊蕴涵命题特性的“单要素模糊推理模式”，式中，Yi，j是代表第i个模式，第j个因子的模糊子集，具有如下明显特点，即当i1，则Yi，j=Yi+

6、1，j-1。（5）根据确定的 m 值，将所有年份径流资料全部统计在内建立模糊控制器，将未来年份前 m年年径流量资料做输入，预报未来年年径流量。1.3 平稳时间序列法（1） 基本概念平稳时间序列是指某种随机变量的时间序列， 它的前期演变过程的统计相关规律在未来的一段时间内是不变的。也就是说它的数学期望与方差是不变的， 它的相关函数只与时间间隔有关而与时间无关。符合上述条件的时间序列是围绕着一条水平线作平稳摆动的。平稳时间序列法预报是利用上述特点， 根据已出现的时间序列推求未来的预报值。（2） 平稳时间序列的自相关函数今有某水文要素观测值的距平值时间序列， 如:X 1 = X 1 - ， X 2

7、= X2 - ，X n =X n - 式中：X1、X2Xn为观测值，X为平均值。将距平值按下列方式排列:X1， X2， X 3，Xn - 2， Xn- 1， X nX2， X3， X 4X n- 1， X nX3， X4， X 5 X n可以看出， 第1行与第2行是时间间隔为1 个单位的两个时间序列； 第1行与第3行是时间间隔为2个单位的两个时间序列。第1行与第2行之间， 有n - 1项是对应的。它们之间的相关函数为R (1) :第1行与第3行两个时间序列之间的相关函数是时间间隔为2个单位的自相关函数， 即二阶自相关R (2)：R（2）=以此类推， 得到时间间隔为个单位的自相关函数R（）=（3

8、） 预报方程的建立利用平稳时间序列作预报的基本依据是， 认为未来的预报值X t+ 1与前期已出现的观测值有关， 即:X t+ 1 = b1X t + b2X t- 1 + + bmX t- m + 1式中: b1、b2、 bm 为预报系数。预报系数需通过解线性方程组求得， 利用实测资料确定预报系数后， 可得预报方程。当建立的预报方程取5项前期量， 即: X t+ 1 = b1X t + b2X t- 1 + b3X t- 2 + b4 X t- 3 + b5 X t- 4时， 线性方程组的形式为:随着所取前期量的项数不同， 方程的形式按上述类推。在实际工作中一般采用小于资料项数的1 /4， 根

9、据实际最佳拟合情况来确定。142三种预报方法的实例应用本文主要通过恒仁水库19672011年7月份的径流资料（表2.1），运用历史演变法、单因素模糊推理法、平稳时间序列法预报2012年7月份的径流量。表2.1 恒仁水库19672011年7月径流资料年份7月年份7月年份7月1967494.5 198285.7199795.51968272.1 1983469.31998317.71969314.2 1984124.31999107.31970144.4 1985896.3200039.41971616.31986570.62001357.2197286.61987137.82002142.519

10、73334.61988222.22003332.01974185.21989612.52004669.41975526.31990396.72005605.8197692.81991665.82006254.21977197.01992223.22007132.1197892.91993179.32008434.71979528.01994247.32009241.21980329.219951052.12010433.01981184.31996619.3201179.4 2.1 历史演变法（1）五种基本规律分析某水库各年7月份的流量变化图如图2.1所示：图2.1 某水库1967到2011年

11、7月份的流量变化图持续性当7月平均流量500m3/s时，第二年的7月平均流量很可能下降。持续时间为1到2年。出现8次。当7月平均流量100m3/s时，第二年的7月平均流量很可能上升至150m3/s以上。出现6次。 当7月份流量连续下降2年或3年，且降到200m3/s以下时，流量很可能上升。出现5次。 当7月份流量连续上涨两年，流量很可能下降。出现3次。 相似性1984年到1988年、1994年到1998年线型为升-降-降-升，且后面的升上去的流量要比前面的的高。 1988年到1992年、2007年到2011年线型为升-降-升-降，且后面的峰值和谷值都比前面的低。 周期性1968年到1972、1

12、973到1976、1977到1980年，三短时间为升-降-升-降型，周期为4年。如果结合1967到2011年的年均流量（图2.2），可以大概确定出枯-平-丰水年周期即：1967年到1975年为枯-平-丰水年第一周期；1976年到1986年为枯-平-丰水年第二周期；1987年到1995年为枯-平-丰水年第三周期；1996年到2005年为枯-平-丰水年第四周期；2006年到2010年为枯-平-丰水年第五周期；2011年为新一周期的开始，即枯水期。图2.2 1967到2011年的年均流量最大最小可能需要定出预测年份7月份可能出现在什么样的水期中，然后再定出相应水期的最大可能性上下限，然后再将上下限平

13、均即为预测值。特征转折点由历年7月降雨量图（图2.1）可观察特征转折点为:1987、1997、2006年与年均降雨量的周期转折点基本类似。（2）2012年7月份降雨量预报分析 根据分析，2011年的7月份平均流量小于100m3/s，因此2012年7月份的平均流量很可能上升，且很可能上升至150m3/s以上。 根据分析，2010年为枯-平-丰水年第五周期的结束，2011年为新一周期的开始，也即枯水期的开始，恒仁水库经常出现连续的枯水年组2，故2012年很可能仍为枯水期，因此可以推断其7月份降雨不会太高突破500m3/s。与2011年7月份流量相似的流量年份为1972年、1976年、1978年、1

14、982年、1997年、1999年、2000年，因此从各年份的第二年的流量中求最大可能性的上限和下限。 去掉超过500m3/s和低于100m3/s的两个值（这两个值代表最小可能性），在剩下的值中取469.3m3/s为最大可能性上限，197m3/s为最大可能性下限。 用上限和下限的算术平均值作为预报值，即333.15m3/s。 2.2 单要素模糊推理法（1）论域的选定 根据水文情报预报规范（GBT 224822008） 延展流量最大值、最小值的20%作为论域范围，最大流量为1052m/s，最小流量为39m/s，则论域为：U=30，1300 （2）模糊子集的确定将论域U离散化，并分5级，即枯、偏枯、

15、平、偏丰、丰。各级间的界限具有模糊性，每级对应一个模糊集Yk，k=1，2，5（k为量级值）。通过各级的P-频率，求出各级的等级值，并根据等级值划分各子集范围，结果见表2.1。表2.1 各模糊子集及范围模糊子集 P-频率 等级值 子集范围 枯 Y1 90%80(30，130 偏枯 Y2 70%180(130，230 平 Y3 50%285(230，340 偏丰 Y4 30%415(340，490 丰 Y5 10%660(490，1300 注：单位 m/s（3）模糊子集隶属度函数及确定模糊分布函数常用的线性分布有梯形分布、右半梯形分布(降半梯形分布)、左半梯形分布(升半梯形分布)和三角分布等，各模

16、糊集的分布形式及隶属度函数如表2.2。表2.2 各模糊集的分布形式及隶属度函数等级模糊集分布形式隶属度函数枯水Y1右半梯形偏枯水Y2三角形续表2.2等级模糊集分布形式隶属度函数平水Y3三角形偏丰水Y4三角形丰水Y5左半梯形（4） 各模糊集代表值代表值隶属度由各模糊集的隶属度函数可求得各级内代表值的隶属度，见表2.3。表2.3论域U离散化分级模糊集代表值隶属度论域U离散化分级模糊集代表值隶属度模糊集123456891011121314区段=80120140180220240285340360415520570660=Y110.60.40000000000(30，130Y200.40.610.60

17、.40000000(130，230Y300000.40.610.60.40000(230，340Y400000000.40.610.60.40(340，490Y500000000000.40.61(490，1300（5） 小前提的确定选择周期m=6。即20062011年为小前提。 根据模糊子集区段确定各年的量级值K分别为： （3、2、4、3、4、1） ，其对应的模糊子集： （Y1*、Y2*、Y3*、Y4*、Y5*、Y6*） （6） 大前提：复合条件语句在m=6前提下，归纳出19672005年39重6次复合模糊条件语句，建立其模糊推理模式，见表2.4（部分）。 表2.4模糊推理模式39重6次复合

18、模糊条件语句（部分）年份iif Yi，1 and Yi，2 and Yi，3 andand Yi，m then Yi，m+1Yi，m+1m+i年12345619671533251319731968233251321974196933251325197519704251325119761971551325121977197261325121197819737325121519791974825121531980197595121532198119761012153211982（7）特征系数的计算根据特征系数计算公式，计算aij和 aij ，结果见表2.5（部分）。表2.5特征系数成果表大前提序号

19、i小前提量级值k3 2 4 3 4 1 最小值100.40.40.40.410210.4000003110.400.40.4040.4001000500.40.40.40.410600.40000.407110.4001080.4000.4000900.4000.4001001000.40.40110.40.40.4101012000.40.400013010.40.40.400141100.4000150.40.4100.4001600000.40.40170.40.40.4000.4018000.40.4000190000.40.400200100100210.410.40.40.40.4

20、0.4220.401000.4023000.40.4000240.4000.40.4002501010.400（8）推理法则大前提描述的前m年要素（预报因子）模糊集与第m+1年待预报要素模糊集之间的模糊关系为： 预报第m+1年要素，把前m年已知的要素模糊集（Y1*、Y2*、Y3*、Y4*、Y5*、Y6*）作为小前提，则模糊推理结论为： 其中： aij为Yj*关于Yij的特征系数（9）结果计算通过特征系数表分析，有信息输出年的年份只有一年即1993年，相应的量级值k=2，即Y2012*的预测量级值为2 。则预测年份2012年属模糊子集（130,230，取平均值，预测流量为：180m/s。 2.3

21、 平稳时间序列法（1）预报方程的建立取用资料系列20062011年6项前期量的距平值预报2012年的流量，预报方程为：Qt+ 1 = b1Qt + b2Qt- 1 + b3Qt- 2 + b4Qt- 3 + b5Qt- 4+b6Qt-5 Qt+ 1 ：预报年距平值 Qt+ 1 ：预报年距平值 即： Q2012= b1Q2011 + b2Q2010+ b3Q2009+ b4Q2008+ b5Q2007+b6Q2006 （2）相关函数r()的计算 另=0、1、2、3、4、5、6，用公式：结合“平稳时间序列法相关函数计算表”（表2.6）计算相关函数，结果见表2.7。表2.6平稳时间序列法相关函数计算

22、表（部分）年份7月流量Qt距平值Qt=0Q2t=1=5=6Qt+1Qt*Qt+1Qt+5Qt*Qt+5Qt+6Qt*Qt+61967494.5 158 24979 -64 -10179 -250 -39493 -2 -303 1968272.1 -64 4148 -22 1434 -2 124 -151 9746 1969314.2 -22 495 -192 4275 -151 3368 190 -4224 1970144.4 -192 36896 280 -53736 190 -36453 -244 46811 1971616.3 280 78263 -250 -69905 -244 -68

23、177 -139 -39025 197286.6 -250 62440 -2 480 -139 34858 -244 60864 2004669.4 333 110835 269 89649 -95 -31725 96 32121 2005605.8 269 72512 -82 -22175 96 25981 -257 -69220 2006254.2 -82 6781 -204 16835 -257 21168 2007132.1 -204 41797 98 -20072 2008434.7 98 9639 -95 -9356 2009241.2 -95 9081 96 -9194 2010433.0 96 9309 -257 -24801 201179.4 -257 66076 均值341.3 38090 -8211 -10039 4596 表2.7 相关函数计算结果0123456r()1-0.0