《水文预报课程设计.doc》由会员分享，可在线阅读，更多相关《水文预报课程设计.doc（25页珍藏版）》请在三一文库上搜索。
1、水 文 预 报课 设 设 计指导老师：王*系别：水资源工程系班级：水文0801学号：20090803*姓名：* 目 录目 录- 0 -一、设计任务- 1 -二、设计资料- 1 -三、流域自然地理概况- 1 -四、设计步骤及技术要求- 3 -1绘制汛期栾川站流量过程线和相应的降雨量过程线- 3 -2 计算流域平均次降雨量P- 3 -3 分析栾川站流量过程的退水规律，制作退水方案- 4 -4 划分洪水，计算各次洪水的实测径流深R- 4 -5.初定蓄满产流模型参数- 4 -6 应用蓄满产流模型计算各次洪水的径流深R- 7 -7.对方案进行精度评定- 7 -8确定蓄满产流模型参数- 11 -附 表-
2、12 -表1 各站累积雨量摘录表- 12 -表2 退水资料摘录表- 13 -表3 退水流量相应径流深计算表- 14 -表4 实测次洪径流深计算表- 14 -表5 产流计算表- 16 -表6 次洪径流深精度统计表- 18 -表7 t=2h时段流量过程摘录表- 19 -表8 由蓄满产流模型计算时段径流量- 20 -表9 用蓄满产流模型计算时段径流深R修正- 20 -一、设计任务本次设计的任务主要包括以下两个方面：（一）编制伊河栾川以上流域降雨径流量预报方案（二）编制伊河栾川以上流域降雨径流过程预报方案二、设计资料 编制方案采用的资料有：1流域出口断面逐时流量；2流域内雨量站逐日降水量；3流域内雨量
3、站逐时降水量；4实测蒸发皿日蒸发量。三、流域自然地理概况伊河栾川水文站建于东径110036，北纬33047，是伊河上游的第一个水文站，集水面积340平方公里，河长36.9公里，流域内属石山、林区、纵坡较陡、森林茂密，植被率在50%左右，土层覆盖较薄，表层土为壤土，下层为沙卵石，下渗能力较大，河道比降陡，河网密度大，河系呈扇形。伊河流域的气候变化，完全受季风的支配，每年的9月份以后季风逐渐南退，地方冷空气不断南下，则为北方冷高压所控制造成季风寒冷，雨雪稀少的干燥气候，6月份以后，夏季风开始活跃，南方暖湿空气不断袭来，形成了夏季炎热，而潮湿的多雨气候。 流域内各站多年平均降水量700800毫米，6
4、10月为汛期，汛期雨量占全年的60%以上，78月为暴雨期，暴雨历时短，强度大，一般降雨历时一至两天，超过三天以上很少出现，暴雨中心多在栾川、陶湾。伊河洪水为暴雨所形成，受降雨特性影响，洪水过程陡涨陡落，洪水多发生在 710月份，特大洪水发生在78月份。流域内水系和站网分布见下图。流域内水系和站网分布图四、设计步骤及技术要求（一）编制产流量预报方案1绘制汛期栾川站流量过程线和相应的降雨量过程线根据1964年和1978年栾川站洪水水文要素摘录表可以绘制出汛期栾川站流量过程线，见图如下。 图1 流量过程线2 计算流域平均次降雨量P由各站累积雨量摘录表（表1各站累积雨量摘录表如下）能够算出1964年和
5、1978年流域平均次降雨量P1964年=74.8mm，P1978年=56.2mm。由各站累积雨量摘录表知道各个时刻的降雨量，然后把各个时刻的雨量相加的该时间段的降雨总量，用降雨总量除以该段时间的天数得到该流域的平均次降雨量。3 分析栾川站流量过程的退水规律，制作退水方案根据已知资料，算出栾川站流量的退水流量过程，其详细计算成果表见表2所示。4 划分洪水，计算各次洪水的实测径流深R水源的划分：蓄满产流条件下，产流量是产流面积上的降雨量，能形成地下径流的也只是产流面积上由稳定下渗率fc形成的水量。因此水源公式为当hfc时，rg= 当hfc时, 一次洪水地下水和直接径流总产量分别为 式中 、就是划分
6、时段和次洪水源的公式。由以上公式可以计算出表3和表4的结果，其成果见附表中的表3及表4所示。5.初定蓄满产流模型参数蓄满产流模型参数主要有蓄水容量曲线参数b，Wm，蒸发模型参数Kc，C分水原参数fc共五个。其中对产流量预报影响最大的是Kc。模型参数在预报中起重要作用，故在建立产流模型前要确定其初值，供模型的调试分析作用。（1）流域稳定下渗率fc的推求Fc值可用实测降雨径流资料反推而得，方法较多，常用试算法。对一次洪水，首先用上一步所述方法求出实测径流深，然后根据、，假设不同的fc，算出，直至由降雨资料计算得与实测值相等，则fc值求定。为了使fc值有一定的精度，选择洪水时应选用次洪产流量计算误差
7、较小和能达全流域产流的洪水，并分析多次，如差别不大，可取平均值。有时，需对fc值作古算，精度要求不高时，可选择雨初基本能达全流域产流而雨强又较大的洪水，算出，并统计出降雨历时T，则由多次洪水分析的fc值可能会差别较大，应该具体分析，不能一概简单地取平均值。如因未计时段内下渗历时，会使fc值受时段均化而变小，时段越长，这种现象越明显，但fc有自身规律，尚需研究。（2） 蒸发模型参数Kc的推求通过分析流域的降雨径流资料，可得到较为准确的初值，应用水量平衡原理，进行流域蒸发接近于Ep时期的水平衡计算，令湿润期为T，在T时段内水量平衡方程为式中为T时段内流域蒸发量，,R为总降雨量和由其形成的产流量，均
8、为实测，为时段始末土壤蓄水量之差。如在汛期中某湿润期时段始末都能蓄满，则=0，湿润期的蒸发也接近于Ep，并假定其相等。另外，统计出T时期蒸发器观测的水面蒸发量,按Kc的意义，可用下式估算出 估算的Kc值中综合了多种误差，故必须多分析几次洪水，求其平均值才比较可靠。另外，也可用多年平均值，据式算出Kc值，但只能作为Kc的下限值，因为实际蒸发量E总是小于Ep的。（3） 的推求一般天然流域，值的确定，通常由实测降雨径流资料分析而得。是流域缺水容量，故应尽量选取前期很干旱（0，缺水量最大），本次降雨又较大，能达到全流域产流（缺水量全部能不足）的洪水，其降雨损失量即值，可多分析几次洪水取其较大者作为初值
9、。有时不易找到满足上述条件的洪水资料，则可找出长期干旱以后的多次洪水，取至达全流域蓄满时作为计算时段，进行水量平衡计算，用下式定出。但此时由于时间较长，故需估算出蒸发期值，可参考同期确定。（4） b的推求B值在部分流域面积产流时起作用，全流域产流时，由于为1.故b值大小与产流量无关。在湿润地区，b只影响小洪水或大洪水的前期产流量（不影响总产流量），故对R影响不大；但干旱地区则局部产流较多而全流域蓄满情况较少，故b值影响较大，根据经验，一般取值为0.20.4，山区大于平坦地区，b值一般随流域面积增大而增大。（5） C的推求C值的确定，由于研究不够，常取经验数据，约0.10.2.一般作物多的地区旱
10、季深根植物散发量较大，C值也较大。6 应用蓄满产流模型计算各次洪水的径流深R（1）确定雨量站，计算流域平均雨量P；（2）分析流域退水曲线，划分Qt，计算次洪实测产流量R0；（3）初定参数b，Kc，C,fc；（4）计算次洪Rc，根据蒸散发计算模型，计算逐日E，并计算出，由结合PR关系（或查图）算出逐日Rc，据此统计出次洪水径流深Rc值，即R值。由以上四个步骤即可计算出各次洪水的径流深R，其计算成果见下表4所示。7.对方案进行精度评定模型的误差有三个方向。一是原始资料的观测误差，二是模型计算方法和参数确定方面的误差，三是模型的理论误差。原始资料的误差，除了较明显的外，一般难以理解，但是具有这样的概
11、念，有助于分析判断。计算方法的误差涉及面比较广，主要是降雨径流关系各要素P，R0，W0的计算及其中的取定，常见一下几个方面。（1） 流域平均降雨量的计算很多流域表明，面雨量的计算误差常成为点据误差的主要因素。首先是测站的控制性问题，在雨量分布不均匀的地区及季节，现有站网往往难以控制面雨量，这个问题实际上很难完全解决。其次，面雨量的各种计算方法基本依据是点雨量对其附近的面积往往有代表性，但是，由于雨型空间多变，加上测站控制性问题，采用何种方法最为合适或如何结合雨型采用不同的方法，须对多年计算结果进行对比分析，了解计算误差后，才能做出判断和改进。一般流域上游降雨较大，而站网稀疏，易漏掉暴雨而使算得
12、的流域平均雨量偏小，这种情况可通过洪水资料的分析比较来确定。分块计算的流域，其子流域的平均雨量也存在上述类似的情况。（2） 次洪R0的计算误差来自采用统一的地下水退水曲线，没有考虑不同退水规律的差异，地下水退水本身也有季节性变化。此外，有些小洪水，虽流量较小，却并非全是地下水，还包括直接径流，往往出现小洪水比大洪水退水快，孤立洪水比连续洪水退水快的现象。因此，在退水曲线变化较大或退水分割误差较大时，要考虑实际情况，取不同的退水曲线或改进计算方法。（3） 蒸散发量计算影响蒸散发计算的参数主要有两个即Kc和Wm。 KcKc决定了Ep，常对蒸散发能力计算起主要作用。如取得不准，将使点据出现系统偏差。
13、如Kc取得过小，Ep也小，则W0偏大，算的R0偏大，点据将偏离理论线上方，反之则偏于下方。若取得合适，则无系统偏离。有些地区Kc还成季节性变化，与植物散发有关，一般是夏季大于冬季。Wm的影响不如Wm明显，它通过蒸发计算来影响蒸发量。若Wm取得过小，又采用一层蒸发模型计算时，会使损失量（Wm-W0）偏小，而使计算的产流量偏大，主要影响久旱后的洪水。一般采用分层模型时，Wm较大，上述情况不易出现。 bB值只影响局部流域产流时的产流量，b值越大，R越大，对于湿润地区，大部分洪水为全流域蓄满，但对半湿润地区，则常有局部产流情况，需要注意的是b值对R的影响。（4） 模型简化生产的误差 将流域作为整体计算
14、如将降雨径流各项要素都取流域平均值，不考虑流域上的不均性，则流域越大，不均匀性程度越剧烈，将影响计算成果。首先是p在面上分布的不均匀性。虽不均匀，但在最小面积上雨末也能蓄满，因全流域都蓄满，损失达最大值，为一常数，则对R计算无影响。因此，只有在p不均匀而形成部分流域产流时，用平均 p计算R与实际才会不同，其点据将偏于理论线下方，即计算径流深偏小。当p在面上分布不均但最终达到全流域产流时，其产流量与p分布不均匀时相同，但产流量时程分配是不同的。其次，在降雨分布明显不均匀的流域，往往引起流域缺水量也不均匀，其影响与p分布不均是相似的。为消除误差，宜于各雨量站分别计算R，再求流域平均R值，这种分站计
15、算法可以大致考虑到p分布不均和下垫面W不均的影响，提高精度。 模型主要考虑下渗损失量模型着重考虑下渗损失量，当其他损失项目突出时，会产生误差，常见的是中小水库蓄水，稻田蓄水等人类活动。这些因素往往在短暂时间内较大地影响流域蓄水量。如久旱之后，由于上列原因产生径流量比正常情况下小。如衢县流域，根据调查，其影响量可达2030mm。显然，也会因后期水库泄水、农田放水而造成相反的误差。这在半湿润地区影响尤为明显。模型的理论误差，即发生超渗产流情况。湿润地区在久旱后如遇短历时暴雨，可能发生超渗产流现象，其中R0大于用蓄满产流计算出的Rc值，点据偏于理论线下方。但根据经验，这种情况并不多见，南方地区并非每
16、年都能发生这种情况。以上分析了蓄满产流模型误差来源及主要影响因素，但对于有误差点据，首先应找出可能产生这种误差的原因，然后逐个分析，排除无关因素，找出主因。如Rc偏大，则可能是R0偏小，p偏小，Kc偏小，b偏大等。然后逐个分析其可能性和论据。其次，分析时要注意是否有系统性，要对所有资料作分析，因为有些因素是有系统影响的，如Kc值，对所有蒸发计算都有影响，b值则对局部产流的洪水有影响。明确了误差的性质和找到原因后，需修改参数或计算方法，对全部资料重新计算，一般需反复多次才能使误差最小，最后确定方案和模型参数。综合以上分析，根据所得的各次洪水径流深对方案进行精度评定，见下表6。 8确定蓄满产流模型
17、参数（1） 确定计算时段T由已知图表知：T=2h（2） 按T摘录流量过程经上一步的计算结果知T=2h，所以按T=2h时段摘录流量过程见下表7。 (3）由蓄满产流模型计算时段径流量由蓄满产流模型计算时段径流量表见表8。（4） 对时段径流量进行修正对时段径流量修正表见表9。】附 表表1 各站累积雨量摘录表日期时间站名时间站名年月日起迄PP起迄PP19647269:0016:003.13.1 11:4018:009.5 9.5 16:0018:0011.014.1 18:0020:006.5 16.0 18:0020:006.420.5 20:004:007.7 23.7 20:0022:005.7
18、26.2 4:006:0015.5 39.2 22:004:003.429.6 6:008:0018.0 57.2 274:006:007.036.6 8:0010:0010.6 67.8 6:008:0014.651.2 10:0012:004.5 72.3 8:0016:008.960.1 12:0014:005.5 77.8 14:0016:006.7 84.5 16:0018:002.5 87.0 18:0018:402.5 89.5 平均次降雨量74.819787420:0022:0019.519.520:0022:0016.4 16.4 22:0023:009.629.122:00
19、0:0024.0 40.4 23:000:0017.546.60:002:009.0 49.4 0:001:007.854.42:004:006.5 55.9 1:002:000.6552:005:001567:008:000.556.5平均次降雨量56.2表2 退水资料摘录表时段t=6小时640727750808780805790712810715830908079.572.693.28076.951.8158.36048.8434838.6244.24836.527.23628.8338.63827.4183024423.434.51313.424.518.4529.531.510.210
20、.32116.1624.826.59.59.61815.9723.722.58.79.11715.7822.3207.38.215.815.5920.118.56.87.314.515.51019186.36.313.215.41118.4165.55.212.515.21216.815511.9151315.314.54.511.214.91413.814410.613.21512.5133.59.611.61610.81238.510.7171010.52.57.59.6188.1102.56.59.4197.39.76.2206.88.55.8216.35.6225.4234.8244.
21、4253.42627282930表3 退水流量相应径流深计算表时间Qt（m3/s）Qt（m3/s）Re(mm)050.0 309.919.7 139.0 259.916.5 231.5 220.914.0 326.0 189.412.0 421.8 163.410.4 518.7 141.69.0 616.1 122.97.8 714.0 106.86.8 812.4 92.85.9 911.2 80.45.1 1010.0 69.24.4 119.0 59.23.8 128.1 50.23.2 137.2 42.12.7 146.3 34.92.2 155.8 28.61.8 165.0 2
22、2.81.4 174.5 17.81.1 184.0 13.30.8 193.5 9.30.6 203.0 5.80.4 212.8 2.80.2 表4 实测次洪径流深计算表1964年时间流量时段长平均流量ww67.216:0032.518:0041.0236.7573.50 73.50 20:0053.6247.394.60 168.10 0:0064.2458.9235.60 403.70 2:0075.5269.85139.70 543.40 4:0083.0279.25158.50 701.90 6:00106.0294.5189.00 890.90 7:00146.01126126.
23、00 1016.90 8:00198.01172172.00 1188.90 9:00251.01224.5224.50 1413.40 10:00322.01286.5286.50 1699.90 10:30322.00.5322161.00 1860.90 10:36323.00.1322.532.25 1893.15 11:00320.00.4321.5128.60 2021.75 11:15320.00.2532080.00 2101.75 11:30308.00.2531478.50 2180.25 12:00302.00.5305152.50 2332.75 13:00285.01
24、293.5293.50 2626.25 14:00258.01271.5271.50 2897.75 15:00235.01246.5246.50 3144.25 16:00207.01221221.00 3365.25 17:00190.01198.5198.50 3563.75 18:00174.01182182.00 3745.75 20:00145.02159.5319.00 4064.75 0:00104.04124.5498.00 4562.75 2:0094.0299198.00 4760.75 4:0079.5286.75173.50 4934.25 6:0070.0274.7
25、5149.50 5083.75 8:0066.8268.4136.80 5220.55 10:0058.3262.55125.10 5345.65 12:0052.8255.55111.10 5456.75 16:0044.2448.5194.00 5650.75 18:0043.0243.687.20 5737.95 20:0039.7241.3582.70 5820.65 22:0038.6239.1578.30 5898.95 0:0036.4237.575.00 5973.95 4:0033.4434.9139.60 6113.55 8:0030.5431.95127.80 6241.
26、35 12:0028.5429.5118.00 6359.35 16:0025.7427.1108.40 6467.75 0:0024.3825200.00 6667.75 8:0023.0823.65189.20 6856.95 20:0018.71220.85250.20 7107.15 1978年时间流量时段长平均流量ww27.820:0011.522:0018.9215.2 30.4 30.4 23:0041.6130.3 30.3 60.7 23:1671.50.27 56.6 15.1 75.7 23:2485.30.13 78.4 10.5 86.2 23:36930.289.2
27、 17.8 104.0 0:00930.493.0 37.2 141.2 0:301270.5110.0 55.0 196.2 1:001960.5161.5 80.8 277.0 1:302430.5219.5 109.8 386.7 2:002430.5243.0 121.5 508.2 3:002351239.0 239.0 747.2 5:001782206.5 413.0 1160.2 6:241231.4150.5 210.7 1370.9 7:001070.6115.0 69.0 1439.9 8:0093.21100.1 100.1 1540.0 9:0086.1189.7 8
28、9.7 1629.7 10:0072.2179.2 79.2 1708.8 12:0056.8264.5 129.0 1837.8 14:0048.8252.8 105.6 1943.4 17:0039343.9 131.7 2075.1 20:0036.5337.8 113.3 2188.4 0:0030.7433.6 134.4 2322.8 4:0024427.4 109.4 2432.2 8:0015.8419.9 79.6 2511.8 20:0010.21213.0 156.0 2667.8 8:008.7129.5 113.4 2781.2 表5 产流计算表时间pE水EmEuEl
29、EP-EWuWLWOR2010.4 0.40.4 0.4 00.4 10.0 000.1 212.82.7 2.7 2.7 -2.7 9.9 9.9 2221.9 1.9 1.9 -1.9 7.3 7.3 233.63.5 3.5 3.5 -3.5 5.4 5.4 245.95.7 1.9 1.9 -1.9 1.9 1.9 2554.8 264.24.0 273.6 3.13.0 3.0 3.0 0.6 0.0 2854.8 0.6 0.6 -0.6 0.6 0.6 298.4 5.35.1 5.1 5.1 3.3 0.01302.9 3.73.6 3.6 3.6 -0.7 3.3 3.3 1
30、54.8 2.6 2.6 -2.6 2.6 2.6 23.53.4 35.55.3 41.2 2.92.8 1.2 1.2 55.55.3 65.95.7 725.0 43.8 3.8 3.8 21.1 0.5 83.23.1 3.1 3.1 -3.1 20.7 20.7 93.33.2 3.2 3.2 -3.2 17.6 17.6 105.45.2 5.2 5.2 -5.2 14.4 14.4 1154.8 4.8 4.8 -4.8 9.2 9.2 1265.8 4.9 4.9 -4.9 4.4 4.4 1376.7 0.0 0.0 0.0 147.26.9 1576.7 160.2 4.8
31、4.6 0.2 0.2 170.4 3.33.2 0.4 0.4 186.96.6 194.64.4 205.35.1 217.77.4 227.26.9 2387.7 248.58.2 256.46.1 264.84.6 270.7 4.34.1 0.7 0.7 280.1 32.9 0.1 0.1 2924.5 2.32.2 2.2 2.2 22.3 0.5 3017.3 0.30.3 0.3 0.3 17.0 21.8 21.8 1.2 3111.0 1.0 1.0 -1.0 37.6 37.6 0.0 13.43.3 3.3 3.3 -3.3 36.6 36.6 0.0 22.6 54
32、.8 4.8 4.8 -2.3 33.4 33.4 36.86.5 6.5 6.5 -6.5 31.1 31.1 45.85.6 5.6 5.6 -5.6 24.6 24.6 543.8 3.8 3.8 -3.8 19.0 19.0 65.45.2 5.2 5.2 -5.2 15.2 15.2 76.56.2 6.2 6.2 -6.2 10.0 10.0 85.14.9 4.9 4.9 -4.9 3.7 3.7 99.59.1 0.6 0.6 -0.6 109.38.9 0.0 119.89.4 0.0 121.3 6.36.0 1.3 1.3 1311.8 11.8 0.1 142.82.7
33、 2.7 2.7 -2.7 11.6 11.6 153.8 8.27.9 7.9 7.9 -4.1 8.9 8.9 166.66.3 4.9 4.9 -4.9 4.8 4.8 1743.8 186.56.2 191.0 6.36.0 1.0 1.0 2032.9 216.56.2 22109.6 237.97.6 2419.5 0.10.1 0.1 0.1 19.4 0.4 2514.5 0.50.5 0.5 0.5 14.0 19.0 19.0 0.8 262.82.7 2.7 2.7 -2.7 32.2 32.2 274.24.0 4.0 4.0 -4.0 29.5 29.5 2834.6
34、 7.16.8 6.8 6.8 27.7 25.5 25.5 2.6 29 5.14.9 4.9 4.9 -4.9 50.6 50.6 309.3 87.7 7.7 7.7 1.6 45.7 45.7 0.2 130.9 0.40.4 0.4 0.4 30.5 47.1 47.1 5.2 218.4 1.91.8 1.8 1.8 16.5 72.4 72.4 4.4 313.0 11.0 1.0 1.0 12.0 84.5 84.5 4.2 459.9 3.93.7 3.7 3.7 56.2 85.0 7.3 92.3 48.5 5 4.0 3.8 3.7 3.7 -3.7 85.0 15.0 100.0 表6 次洪径流深精度统计表产流方案误差分析表洪号实测径流深计算径流深允许误差实际误差是否合格74100254.652.310.9 2.378071014.315.82.9 -1.581071864.8