洪峰流量预测

1、2013年河南科技大学数学建模选拔赛承 诺 书我们仔细阅读了数学建模选拔赛的规则.我们完全明白，在做题期间不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与选拔题有关的问题。我们知道，抄袭别人的成果是违反选拔规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守选拔规则，以保证选拔的公正、公平性。如有违反选拔规则的行为，我们将受到严肃处理。我们选择的题号是（从A/B/C中选择一项填写）： A 队员签名 ： 日期：2013年8月23日2013年河南科技大学数学建模竞赛选拔

2、编 号 专 用 页评阅编号（评阅前进行编号）：评阅记录（评阅时使用）：评阅人评分备注洪峰流量预测摘要本文讨论了关于百年一遇的定义、百年中最大洪峰流量预测、怎样根据现有的年最大洪峰样本来预测未来三年的洪峰值等问题。首先对于问题一，通过查阅资料并客观的分析，可认为“百年一遇”不能简单机械的理解为“一百年只出现一次”。百年一遇应该理解为数学统计中的“概率”，指的是随机事件发生的概率，表现在时间上就是我们平常所说的“时间周期”。百年一遇，就是指超过一定标准的洪峰流量超过一定标准发生的概率是1%。对于问题二，考虑到洪峰流量的分布，在较大值和较小值分布较少，处于一般情况下的分布最多，因此，考虑正态分布建立

3、模型。我们通过spss进行正态分布拟合的探索。然后我们对拟合的正态分布的正态性进行检验，结果发现正态性良好，我们建立的正态分布模型成立。然后我们求得此正态分布概率密度函数，然后对正态分布概率积分函数进行求解，得到概率积分为0.99的临界点，作为“百年一遇”的观测值，即超过此值，我们就称为“百年一遇”。此观测点为百年一遇的洪峰流量值。对于问题三，题中要求对后三年的洪峰量进行预测，再分析了题目后我们选择了用灰色理论对数据进行预测，建立灰色模型求解预测结果三年内结果均无明显变化，对模型检验发现模型并不合理，后来我们又对模型进行了改进，用马尔科夫与灰色理论混合模型能够比较好的通过检验并预测了后三年的洪

4、峰量为：年份201320142015洪峰值583.67516.67513.95关键字:正态分布 马尔科夫预测 灰色预测问题重述外界对三峡工程“万年一遇”“千年一遇”“百年一遇”等防洪标准说法不一提出质疑，并将相关报道整理如下：现假设附件中是某水文站的每年的最大洪峰流量观测值。1、有人说，百年一遇就是一百年内只出现过一次，能否这样理解？给出你对“百年一遇”的定义。2、 能否计算出附件中水文站百年一遇的洪峰流量值；如能，给出具体的计算过程及结果。3、 预测该水文站将来3年内每一年的最大洪峰流量。问题分析对于问题一，如何定义“百年一遇”，首先否定了“百年一遇就是一百年内只出现过一次”的理解。经过查相

5、关资料所得：所谓“百年一遇”这是一个关于频率的概念。它表示在许多次试验中某一事件重复出现的时间间隔的平均数。需要特别指出的是所谓“重现期”并不是说正好多少年中出现一次，它带有统计平均的意义，说得更确切一点是表示某种水文变量大于或等于某一指定值，每出现一次平均所需的时间间隔数。对此，可以通过随机模拟来对此进行解释，通过MATLAB随机产生一百个随机数来代表这一百年的洪峰，用其中的最大洪峰值来代表百年一遇的最大洪峰值，再用随机抽样的方法确定最大洪峰值发生的概率。从而可以确定上述说法的错误在于把平均发生的次数认为了为一定会发生。准确定义应用概率来定义为：在一百年中每天最大洪峰发生的概率为1%。对于问

6、题二求最大洪峰值,可以根据现有的某水文站的数据，将其分成间距相同的组，计算出各个组中的频数，通过计算出这些样本的平均值和方差来估计整体的均值和方差，从而确定洪峰值分布的概率分布曲线。通过确定洪峰值大于等于百年洪峰值得概率为1%所对应于的洪峰值作为百年一遇的最大洪峰值。对于问题三，可以采用灰色理论预测模型，预测该水文站将来3年内每一年的最大洪峰流量。基本假设1、所给数据真实有效，水文站的检测结果可靠无差。2、每年最大洪峰流量只考虑时间的影响，附表没给出影响因素，不加考虑。3、不考虑环境剧烈变化所引起的对预测结果的影响。符号说明带估参数第年份的最大洪峰值第个叠加序列的值模型的建立与求解对于问题一：

7、经查阅相关资料可以得到关于百年一遇的理解: 所谓“百年一遇”这是一个关于频率的概念。往往用“重现期”来替代“频率”，它表示在许多次试验中某一事件重复出现的时间间隔的平均数。需要特别指出的是所谓“重现期”并不是说正好多少年中出现一次，它带有统计平均的意义，说得更确切一点是表示某种水文变量大于或等于某一指定值，每出现一次平均所需的时间间隔数。经分析可以用随机模拟的方法来解释，先用MATLAB随机产生一百个数如附录：用随机抽样的方法对数据进行单次抽样，进行一百次，分别代表一百年，每个数分别代表每一年的洪峰值，最大的数可以作为百年一遇的最大洪峰值。用excel进行一百次的随机抽样并统计不同阶段的的概率

8、得：流量0-100100-200200-300300-400400-500500-600600-700700-800800-900频率0.010.020.150.430.310.070.0100并得出其直方图：可以得出最大流量值得概率为：P=0.01。则发生次数均值为，百年中发生次数为1仅仅是在大量实测统计的情况下，百年发生的均值。但是在实际中每一年都有0.01的概率发生。所以不能说是一百年中一定会发生。我们期望找到一个的洪峰流量值，作为一个观测点，即衡量的标准，超过此洪峰流量，我们就称之为“百年一遇”。此观测点，应能保证有99%每年最大洪峰流量都小于此观测点，即大于此观测点的每年最大洪峰流量

9、出现的概率为1%，引入重现期，即重现期为100年。对于问题二：考虑到每年的洪峰流量的值，分布的范围在，很大或很小的情况很少，而大多数都分布在洪峰流量既不很大，也不很小的范围内，由此我们联系到所学的正态分布，可以考虑建立起一个关于洪峰流量的正态分布的模型。1、 模型的建立建立以下正态分布的模型：2、模型求解洪峰流量638.10658.99577.11502.42384.55455.04405.57555.06510.22578.33470.90洪峰流量834.00703.92597.79521.35535.82348.53463.45502.04597.32472.09433.74最大洪峰流量6

10、10.67512.23465.35555.56485.52672.23通过spss中分析探索描述性统计，可得到正态分布的，。以下是spss处理的结果：可得到，。正态性检验：经查询资料，得知，正态分布的正态性的两个评判标准为偏度和峰度，以及Q-Q图。（1）、偏度，峰度检验方法通过计算和及其标准误及然后作检验。两种检验同时得出，即在95%的置信概率下，可以认为正态性良好。下面给出spss做出的偏度、峰度检验的结果：描述统计量标准误最大洪峰流量均值537.42519.7669均值的 95% 置信区间下限496.867上限577.9835% 修整均值532.936中值516.790方差10940.44

11、0标准差104.5966极小值348.5极大值834.0范围485.5四分位距130.9偏度.733.441峰度1.122.858可得 ； 计算得：偏度检验： 峰度检验：由检验结果可知，在95%的置信概率下，偏度、峰度检验均合格，即在偏度、峰度检验下，正态分布的正态性良好。（2）、Q-Q图检验以样本的分位数作为横坐标，以按照正态分布计算的相应分位点作为纵坐标，把样本表现为指教坐标系的散点。如果资料服从正态分布，则样本点应该呈一条围绕第一象限对角线的直线。我们根据spss作出Q-Q图如下：由图可知，28个观测点中，由27个都是很好的围绕（几乎完全分布）在第一象限对角线上，只有第28个观测点，偏离

12、稍远。因此，我们完全可以认为Q-Q图检验的结果为，其探索的正态分布的正态性良好。（4）、正态分布模型即： 我们此时是由样本来估计总体，由样本求出的，均是对总体，的无偏估计。故我们建立的这一正态分布模型可以推广到总体更一般的情形。 则，总体最大洪峰流量的概率密度函数为：下面是我们用matlab做出的总体正态分布曲线图3、对“百年一遇”的最大洪峰流量值得计算（1）、 引入新的变量，对于一般的正态分布,的值可以利用以下变换公式计算：注：此变换公式为摘录，其中变量与我们建立模型所用变量含义可能不同，仅供参考。故我们做类似变换得到：不同的值，可由正态分布积分表查出。由我们在第一问给出的定义，“百年一遇”

13、即大于此洪峰流量的概率小于1%，即“百年一遇”的临界值，小于此洪峰流量值的概率积分为99%查表得出，得出了“百年一遇”的最大洪峰流量值临界点代入数据： 得：由此，我们得出了总体的“百年一遇”的洪峰流量值的临界点，即当某年的最大洪峰流量大于此临界点，我们就称为百年一遇。对于问题三（1）模型建立通过对问题的分析,得出对洪峰值的预测可以建立灰色理论预测模型。建立某水文站从1985到2012年最大洪峰值的原始时间序列为：将原始数据按如下方式累加：得到新的序列：令为的紧邻均值生成序列：则GM(1,1)的灰微分方程为：其中：称为发展系数；称为内生控制灰数。设为带估参数，可利用最小二乘法求解。解得： 其中

14、用MATLAB解得: 可得:带入参数得:还原值可得：则可通过MATLAB求得2013、2014、2015年最大洪峰值:对灰色模型的检验：年份误差年份误差19850.00%19992.31%198623.25%20000.41%19877.95%200134.69%19886.01%200213.13%198928.05%20035.88%199014.85%200412.00%199124.09%200511.47%19923.9%200618.65%19934.48%200714.55%19948.29%20083.9%199511.82%200912.68%199656.20%20104.

15、26%199731.86%20118.87%199811.99%201226.19%（2）模型反思与改进我们通过灰色理论预测2013,2014,2015的最大洪峰流量：发现预测的结果，每年最大洪峰流量的值几乎不变，这与我们实际每年最大洪峰流量的波动不符，并且误差较大，所以我们决定对这一模型，进行改进。改进思路：考虑到每一年的最大洪峰流量和下一年的最大洪峰流量之间有较大的相关性，所以我们采用马尔科夫状态转移矩阵来改进。首先进行数据处理，把洪峰流量进行一个分级：800-900 一级；700-800 二级；600-700 三级；500-600 四级；400-500 五级；300-400 六级。我们可

16、以由附表得到每一年的最大洪峰流量等级和下一年的最大洪峰流量等级之间存在一定的联系，所以建立由该水文站上一年的最大洪峰值转移到下年最大洪峰值时转移概率矩阵如下：次年当年一级二级三级四级五级六级一级010000二级000100三级001/32/300四级0001/21/31/6五级1/801/43/81/40六级000010因此可得到马尔科夫状态转移矩阵为：由附表我们知道：2012年最大洪峰流量值为672,23，属于三级，故对应的2013年最大洪峰流量值出现在三级的概率为，出现在四级的概率为，出现在其他等级的概率为0。故预测2013年最大洪峰流量值为（650为第三等级范围中值，550为第四范围等级

17、的中值）2013年预测的最大洪峰流量值为583.33，属于四级。由MATLAB可得： 故对应的2014年最大洪峰流量值出现在三级的概率为0.1111，出现在四级的概率为0.5556，出现在五级的概率为0.2222，出现在六级的概率为0.1111。故预测2014年最大洪峰流量值为故对应的2015年最大洪峰流量值出现在一级的概率为0.0278，出现在三级的概率为0.096，出现在四级的概率为0.4352，出现在五级的概率为0.3519，出现在六级的概率为0.0926。故预测2015年最大洪峰流量值为模型评价模型优点：1、本文对百年一遇的解释巧妙地利用了随机模拟的模型，充分的解释了百年一遇是一种关于

18、概率的概念。2、本文对问题二最大洪峰值得预测充分的利用了正态分布曲线的概率模型，巧妙地利用概率统计的方法得出最大洪峰值。模型的不足：1、对于用灰色模型来预测未来三年的洪峰值，经过验证并不理想。3、对洪峰值得等级分配会带来一定的误差。参考文献1、党耀国,灰色预测与决策模型研究,科学出版社,2009-12-01 2、彭放 杨瑞琰 罗文强 肖海军 何水明,数学建模方法,科学出版社 20073、费业泰,误差理论与数据处理,机械工业出版社,2010-6附录随机数：625.33 591.82 313.76 506.70 841.59 398.10 500.41 395.00 539.82 517.01 4

19、44.54 535.33 449.61 448.73 623.74 293.73 451.28 585.02 510.02 514.66 547.89 533.79 679.11 420.26 681.66 385.85 372.47 448.65 354.36 505.72 480.47 453.94 425.28 801.64 426.74 572.31 590.56 502.40 507.55 539.84 569.17 643.98 637.94 710.59 488.41 578.36 566.92 595.24 450.47 542.79 474.63 523.49 550.40

20、569.59 508.93 584.67 540.93 646.11 435.00 623.83 588.67 462.69 687.70 405.93 652.32 523.80 397.93 420.52 416.47 697.14 614.76 678.78 332.32 446.90 508.36 556.64 655.36 669.29 481.62 586.26 716.48 513.91 516.75 518.96 610.80 487.62 574.05 606.47 327.96 515.49 517.12 475.81 411.09 620.20 322.81 627.59

21、 506.14 530.33 638.28 602.82 用excel抽样数表：841.5939638.2805586.2628569.5923539.8355517.1203506.1436451.2803446.8996420.5233710.5851627.5899585.0197569.5923539.8355517.0054481.6165451.2803444.5395405.9342710.5851627.5899585.0197556.6383539.8355514.664481.6165451.2803444.5395332.3207697.1422623.8287584.6

22、693556.6383539.8192514.664480.4689450.4666444.5395332.3207697.1422620.2047578.3594556.6383530.3265508.9255475.8143450.4666435.003332.3207681.6613614.7602578.3594547.8944530.3265508.3604474.6326449.6065426.7425322.8112678.7759610.8041572.3092542.7898518.9601507.5468474.6326448.7345426.7425322.8112678.7759602.818572.3092540.9269518.9601507.5468453.9395448.7345425.2808322.8112646.1108588.6739572.3092540.9269517.12