中长期水文预报方法

水文预报大作业I水文预报大作业中长期水文预报方法作者姓名：***学科、专业：水文预报理论与方法学号：***指导教师：***完成日期：***水文预报大作业II摘要本文运用历史演变法、单要素模糊推理法、平稳时间序列法三种中长期预报方法，分别对某水库1967—2011年7月份平均流量进行处理，进而预报2012年7月份平均流量。三种方法的预报结果分别为：333.15m³/s、180m³/s、258.6m³/s。文中具体介绍了各方法的计算过程，并在文末对各方法的优缺点进行了评价。关键词：历史演变法；单要素模糊推理法；平稳时间序列法水文预报大作业-III-目录摘要...........................................................................................................................II前言...............................................................................................................................-1-1预报方法原理............................................................................................................-1-1.1历史演变法.....................................................................................................-1-1.2单要素模糊推理法.........................................................................................-2-1.3平稳时间序列法............................................................................................-2-2三种预报方法的实例应用.............................................................................................52.1历史演变法.........................................................................................................52.2单要素模糊推理法.............................................................................................72.3平稳时间序列法...............................................................................................123.方法结果对比与分析...................................................................................................14水文预报大作业-1-前言根据前期水文气象要素，用成因分析与数理统计的方法，对未来长时期(一旬至一年)的水文要素进行科学的预测称为中长期水文预报。中长期水文预报通常泛指预见期超过流域最大汇流时间，且在3天以上，1年以内的水文预报。其预报内容主要有：径流量、河口的水位、冰清、河道及冲於变化和旱涝趋势等。本文主要运用历史演变法、单要素模糊推理法及平稳时间序列法对某水库1967年—2011年7月份平均流量进行分析，并预测2012年7月份的平均流量。1预报方法原理1.1历史演变法（1）基本思想任一水文气象要素的长期历史实测值，全面地反映了这一要素的时序变化规律。尽管目前还不能辨认影响要素变化的内因与外因或影响程度，但是诸因素的综合影响却已全面地融汇于要素的时序变化之中，若能找出该要素历史变化规律。便可用它进行预报（2）基本思路任一水文气象要素的长期历史实测值，全面地反映了这一要素的时序变化规律。尽管目前还不能辨认影响要素变化的内因与外因或影响程度，但是诸因素的综合影响却已全面地融汇于要素的时序变化之中，若能找出该要素历史变化规律，便可用它进行预报。（3）基本出发点任何一个水文要素的长期记录反映了这一要素全面的变化过程。只要我们能够找出这一要素的演变规律，就可利用这些规律来进行预报。（4）基本规律①持续性：连续上升或下降，升降超出平均值某一数值后将改变原趋势的特性。②相似性：两个水文时期水文要素变化趋势的相似性。③周期性：丰、平、枯周期变化。④最大最小可能性：指水文气象要素经常出现的范围或年上升、下降的极限幅度。⑤转折点：两个时期之间的转折期为转折点。水文预报大作业-2-1.2单要素模糊推理法单要素模糊推理建立的基本思想为：（1）统计待预报的水文要素x（年平均流量或年平均降雨量等）的历史资料，取其最大、最小值域（或适当外延）作为论域U。（2）将论域U上的x值离散化，并分成l级（如l=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）就是具有（n-m）重m次复合模糊蕴涵命题特性的“单要素模糊推理模式”，式中，Yi，j是代表第i个模式，第j个因子的模糊子集，具有如下明显特点，即当in-m，j1，则Yi，j=Yi+1，j-1。（5）根据确定的m值，将所有年份径流资料全部统计在内建立模糊控制器，将未来年份前m年年径流量资料做输入，预报未来年年径流量。1.3平稳时间序列法（1）基本概念平稳时间序列是指某种随机变量的时间序列，它的前期演变过程的统计相关规律在未来的一段时间内是不变的。也就是说它的数学期望与方差是不变的，它的相关函数只与时间间隔有关而与时间无关。符合上述条件的时间序列是围绕着一条水平线作平稳摆动的。平稳时间序列法预报是利用上述特点，根据已出现的时间序列推求未来的预报值。（2）平稳时间序列的自相关函数（1.1）水文预报大作业-3-今有某水文要素观测值的距平值时间序列，如:△X1=X1-X——，△X2=X2-X——，……△Xn=Xn-X——式中：X1、X2……Xn为观测值，X为平均值。将距平值按下列方式排列:△X1，△X2，△X3，……△Xn-2，△Xn-1，△Xn△X2，△X3，△X4……△Xn-1，△Xn△X3，△X4，△X5……△Xn……可以看出，第1行与第2行是时间间隔为1个单位的两个时间序列；第1行与第3行是时间间隔为2个单位的两个时间序列。第1行与第2行之间，有n-1项是对应的。它们之间的相关函数为R(1):第1行与第3行两个时间序列之间的相关函数是时间间隔为2个单位的自相关函数，即二阶自相关R(2)：R（2）=以此类推，得到时间间隔为τ个单位的自相关函数R（τ）=（3）预报方程的建立利用平稳时间序列作预报的基本依据是，认为未来的预报值Xt+1与前期已出现的观测值有关，即:Xt+1=b1Xt+b2Xt-1+……+bmXt-m+1式中:b1、b2、……bm为预报系数。预报系数需通过解线性方程组求得，利用实测资料确定预报系数后，可得预报方程。当建立的预报方程取5项前期量，即:水文预报大作业-4-△Xt+1=b1△Xt+b2△Xt-1+b3△Xt-2+b4△Xt-3+b5△Xt-4时，线性方程组的形式为:随着所取前期量的项数不同，方程的形式按上述类推。在实际工作中一般采用小于资料项数的1/4，根据实际最佳拟合情况来确定。水文预报大作业52三种预报方法的实例应用本文主要通过恒仁水库1967—2011年7月份的径流资料（表2.1），运用历史演变法、单因素模糊推理法、平稳时间序列法预报2012年7月份的径流量。表2.1恒仁水库1967—2011年7月径流资料年份7月年份7月年份7月1967494.5198285.7199795.51968272.11983469.31998317.71969314.21984124.31999107.31970144.41985896.3200039.41971616.31986570.62001357.2197286.61987137.82002142.51973334.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.42.1历史演变法（1）五种基本规律分析某水库各年7月份的流量变化图如图2.1所示：图2.1某水库1967到2011年7月份的流量变化图水文预报大作业6①持续性当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、1973到1976、1977到1980年，三短时间为升-降-升-降型，周期为4年。如果结合1967到2011年的年均流量（图2.2），可以大概确定出枯-平-丰水年周期即：1967年到1975年为枯--平--丰水年第一周期；1976年到1986年为枯--平--丰水年第二周期；1987年到1995年为枯--平--丰水年第三周期；1996年到2005年为枯--平--丰水年第四周期；2006年到2010年为枯--平--丰水年第五周期；2011年为新一周期的开始，即枯水期。图2.21967到2011年的年均流量水文预报大作业7④最大最小可能需要定出预测年份7月份可能出现在什么样的水期中，然后再定出相应水期的最大可能性上下限，然后再将上下限平均即为预测值。⑤特征转折点由历年7月降雨量图（图2.1）可观察特征转折点为:1987、1997、2006年与年均降雨量的周期转折点基本类似。（2）2012年7月份降雨量预报分析①根据分析，2011年的7月份平均流量小于100m3/s，因此2012年7月份的平均流量很可能上升，且很可能上升至150m3/s以上。②根据分析，2010年为枯-平-丰水年第五周期的结束，2011年为新一周期的开始，也即枯水期的开始，恒仁水库经常出现连续的枯水年组[2]，故2012年很可能仍为枯水期，因此可以推断其7月份降雨不会太高突破500m3/s。③与