(3)、制作逐日面雨量表 3.2 次降雨量的计算
把一次降雨量累加起来就是该段时间的次降雨量,分别把所有时段的次降雨量求出列入一栏表中。 3.3 逐日面雨量、逐日平均流量过程线;
选中时间再选中对应的流量值和面雨量绘画逐日面雨量、逐日平均流量过程线图。
4 产流预报方案评定
具体方法
(1)、在三年中选出退水流量标上对应的洪号,再通过计算求出消退系数Cg=各时段比值的均值和的均值,k=-1/Ln Cg
(2)、选出若干段由小到大到小的流量,标注洪号,求实测径流: Ro=(流量总和-(第一个流量-最后一个流量)*K)*24*3600/1000/F。
(3)、假设k值(0.8-0.99),求前期影响雨量Pa:Pa1=5月下旬流量*5月流量*0.281 Pa2?Pan=IF(P+ Pa1)*k>100,100, (P+ Pa1) *k) ?
(4)、计算该洪号对应的∑次降雨量+前期影响雨量Pa
(5)、制作表头:洪号、次降雨P、实测径流Ro、Pa、Pa+ P、R计 (6)、绘制 实测径流Ro和Pa+ P相关图,从P+Pa曲线中查处R计值。 (7)合格率的分析;方案比较。
①、制作表头:洪号、次降雨P、实测径流Ro、Pa、Pa+ P、R计、误差、允许误差、合格率。 ②、误差=ABS(实测径流Ro- R计),允许误差=IF(误差<=3,3,IF(误差>=20,20,0.2*C3)),合格率=IF(误差<=允许误差,\合格\不合格\。
表3-2 降雨径流预报方案评定表
预报值 序号 洪号 (m/s) 3实测值 (m/s) 3预报误差 (m/s) 3允许误差 合格否 (m/s) 3- 6 -
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 57822 57828 59718 59723 60621 60807 60824 61707 61804 61812 61823 62721 63721 63723 63725 63815 64722 64807 64815 64820 65728 65803 51.5 38 29 32 15.8 67.7 38 67.5 16.9 19.4 20.4 24.5 22.1 41.2 55.5 18.9 20.5 30.5 41.7 67 39 42.5 55.2 43.9 30.5 23.4 18.8 68.4 43.7 66.9 15.6 19.9 20.1 20.6 16 39.2 47 21.9 21.1 32.9 45 79.7 18.1 40.9 3.7 5.9 1.5 8.6 3 0.7 5.7 0.6 1.3 0.5 0.3 3.9 6.1 2 8.5 3 0.6 2.4 3.3 12.7 20.9 1.6 11.04 合格 8.78 合格 6.1 合格 4.68 不合格 3.76 合格 13.68 合格 8.74 合格 13.38 合格 3.12 合格 3.98 合格 4.02 合格 4.12 合格 3.2 不合格 7.84 合格 9.4 合格 4.38 合格 4.22 合格 6.58 合格 9 合格 15.94 合格 3.62 不合格 8.18 合格 - 7 -
23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 65813 65905 66812 66815 66817 67731 68721 69719 69725 70716 70730 71612 71707 71721 74626 75703 77709 80615 80616 80707 80710 81627 52 60 20 23.3 53 26.5 14.5 15.7 22.3 33 21.8 14.4 27.5 25.9 25 30 33 18.1 63.7 17.6 39 43.7 47 63.9 9.2 23.4 51.1 29.4 18.8 32.9 21.9 31.4 20.3 47.2 29 30.7 32.9 34.6 31 19.9 42.2 24 37.7 47.8 5 3.9 10.8 0.1 1.9 2.9 4.3 17.2 0.4 1.6 1.5 32.8 1.5 4.8 7.9 4.6 2 1.8 21.5 6.4 1.3 4.1 9.4 合格 12.78 合格 1.84 不合格 4.68 合格 10.22 合格 5.88 合格 3.76 不合格 6.58 不合格 4.38 合格 6.28 合格 4.06 合格 9.44 不合格 5.8 合格 6.14 合格 6.58 不合格 6.92 合格 6.2 合格 3.98 合格 8.44 不合格 4.8 不合格 7.54 合格 9.56 合格 - 8 -
45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 81709 81718 82829 83717 83723 83804 84802 84813 85807 85815 86715 87704 87720 92730 94704 97831 39 67 55.2 50.5 38.3 24 7.9 37 43 33.9 19.5 45 43 52.5 68 70 47.2 87.7 31.8 49.4 32.3 16 15 38.3 46.1 33.1 26 25.4 37.2 34 52.3 60.8 8.2 20.7 23.4 1.1 6 8 7.1 1.3 3.1 0.8 6.5 19.6 5.8 18.5 15.7 9.2 9.44 合格 17.54 不合格 6.36 不合格 9.88 合格 6.46 合格 3.2 不合格 3 不合格 7.66 合格 9.22 合格 6.62 合格 5.2 不合格 5.08 不合格 7.44 合格 6.8 不合格 10.46 不合格 12.16 合格 允许误差公式 允许误差上限20mm 方案合格率 允=Qm*20% 允许误差下限3mm 70.00%
③计算成果及其分析:
产流预报方案模型建立时,判断条件较多,要反复的试算,才能获得最优的参数组合,确立优等方案。
- 9 -
计算工作量大,而且按传统方法计算,要查图,带来人为误差,违背模型最初的宗旨。本方法完全按原理用公式计算,避免了这一现象的发生。提高计算机精度,减少计算工作量。方案建立后,对于率定合格的参数、采用相应的计算模型,在发生降雨后,同样采用本方法进行产流预报,非常简便易行,在预报过程中,如果发现诸如土壤含水量等有明显的误差,还可以做现实改正。如果方案合格率较低或者想进一步提高合格率,可以在产流模型和蒸发模型计算表中及方案评定表中修改参数直到满意为止。
5 单位线的推求
5.1 基本概念
在单位时段△T内,由均匀分布在流域上的一个单位地表净雨深,在流域出口断面所形成的地表径流过程线,称为单位过程线,简称单位线。在我国,单位地表净雨深常为10mm。时段△T一般选用3、6、12小时等。在分析单位线时。常基于如下假定:(1)流域的下垫面因素不变,单位线形状不变。(2)如果单位时段内的净雨深是nmm,它所形成的出流过程,总历时与单位线总历时相同,流量Qi是单位线流量qi的n/10倍。(3)如果净雨历时是m个时段,则各时段地表净雨深所形成的出流过程互不干扰,出口断面的流量等于m个出流过程错加时段叠加而成,把流域汇流看成是一个线性时不变系统。它反映了流域特性对单位净雨的响应,属于一种典型的黑箱子分析方法,它依据倍比和迭加假定,由已知的输入、输出来确定系统的作业特性,不顾及由净雨转换为洪水过程时的各环节间的物理成因联系。单位线是指单位时段内,均匀分布的单位净雨量在流域出口断面形成的地面径流过程线。利用单位线来推求洪水汇流过程线,称单位线法。 瞬时单位线集中反映了流域的汇流特性 ,便于流域暴雨洪水关系的理论研究和地区综合 ,在流域防洪兴利工程计算领域具有重要意义。目前常用矩法估计瞬时单位线 ,据此还原的洪水过程与实测洪水过程之间 ,特别是洪峰点之间 ,常具有明显的差别。为此 ,提出通过用加速遗传算法求解还原洪水过程与同次实测洪水过程加权残差绝对值和最小化问题 ,来优化估计瞬时单位线的一套数值方法。实例研究的结果表明了该方法的有效性 ,计算精度高于矩法。单位线时段选3h、6h、12h、 24h等单位净雨深一般采用10mm.时段单位线(UH)的推求传统的的时段单位线推求方法有分析法、图解法、试错法、最小二乘法。
单位线--指在单位时段内流域上均匀分布的单位净雨,在流域出口断面所形成的径流过程线。 单位线是一种特定的流量过程线,它反映了流域特性对单位净雨的响应,属于一种典型的黑箱子分析方法,它依据倍比和迭加假定,由已知的输入、输出来确定系统的作业特性,不顾及由净雨转换为洪水过程时的各环节间的物理成因联系。
分析和应用单位线时,迭加和倍比假定具体叙述如下。
(1)倍比假定 同一流域上,若两次净雨的历时相同,但净雨深不同,各为h1,h2,则二者所形成的地表径流
过程线形状相似.即洪水总历时、涨洪和退水历时完全相同,相应时段的流量座标与净雨量成正比,Qa1/Qb1=h1/h2。
(2)迭加假定 同一流域上,两相邻单位时段Δt的净雨深h1、h2各自在出口断面形成的地表径流过程线
- 10 -