目录一、课程设计的目的和任务..............................2二、课程设计预报流域概况及使用的资料..................2(一)流域概况..............................................2(二)基本资料..............................................2三、课程设计的原理和内容..............................3(一)流域产流方式的初步论证................................3(二)径流预报方案..........................................31、用分析法计算径流量..................................32、由实测资料计算次洪径流量R实........................103、成果分析...........................................12(三)直接径流汇流过程预报方案............................131、选用资料...........................................132、从实测流量过程线上分割地下径流.....................143、确定单位线时段长∆t.................................184、从实测降雨过程求净雨过程...........................185、单位线分析.........................................276、分析由两次洪水所得出的单位线,确定单位线预报方案...29(四)地下径流汇流过程预报方案。..........................29(五)降雨径流过程预报方案的校核预报......................30四、设计的体会与建议..........................................34《水文预报》课程设计报告—密塞流域降雨径流预报方案的制作一、课程设计的目的和任务课程设计的任务是制作湿润地区中小流域的降雨径流预报方案。通过本次课程设计,掌握制作预报方案的全部基本操作及计算工作,并对设计成果能运用课堂所学理论加以分析论证,将课堂理论知识应用于实践,从而对所学知识有更好的理解。二、课程设计预报流域概况及使用的资料(一)流域概况本次预报流域为浙江省衢江密塞流域,流域面积797km2,流域内有六个雨量站:齐溪、马金、燕溪、大溪边、皇林坑和密塞,一个蒸发站:衢县蒸发站(80cm套盆式)。流域系温湿山丘区,雨量丰沛,植被良好。年雨量在1500~2000mm左右,年径流系数在0.7~0.8左右,每年4~6月的梅雨期为本地区的主要雨季,降雨量可达全年的60%以上。梅雨的特点是雨日多,雨强不大,地区分布较不均匀,7月上旬梅雨结束,进入盛夏,降雨以热雷雨为主,形成一个相对干旱的季节,但遭强台风侵袭,也可形成大暴雨,热雷雨和台风雨往往雨强大,雨时短,地区分布也不均匀。11月至次年2月,通常降雨量稀少,为河流的枯水季节。由于流域植被良好,土壤发育,表土下渗率较高,地下径流较丰富,约占年径流的40%左右。(二)基本资料本次课程设计编制方案使用的资料为1986年:(1)流域出口断面逐日平均流量;(2)流域内雨量站逐日降水量;(3)实测蒸发皿逐日蒸发量;(4)实测流域出口断面逐时流量;(5)流域内雨量站逐时降水量;(6)实测蒸发皿逐时蒸发量。三、课程设计的原理和内容(一)流域产流方式的初步论证由于密塞流域系温湿山丘区,每年4~6月的梅雨期为本地区的主要雨季,降雨量可达全年的60%以上,这段时间雨日多,雨强不大,雨量丰沛,年雨量大于1000mm,在1500~2000mm左右,年径流系数大于0.4,在0.7~0.8左右。而该流域植被良好,土壤发育,表土下渗率较高,年蒸发量不大,且地下径流较丰富,约占年径流的40%左右。属典型的蓄满产流区。(二)径流预报方案1、用分析法计算径流量(1)产流系数的分析在蓄满产流情况下,流域产流的过程及径流量的计算公式可由流域蓄水容量分布曲线来表示。据大量经验分析,蓄水容量曲线可由如下指数方程近似描述:α=1−(1−WM′WMM⁄)b(3−1)式中b是常数,反映了流域包气带蓄水容量分布的不均匀性,b值越小表示越均匀,当b=0时表示流域内包气带蓄水容量均匀不变,而b值越大表示越不均匀。据上式,流域平均蓄水容量WM为WM=WMM(1+b)⁄(3−2)经计算,流域平均的初始土壤含水量W的最大值a为a=WMM[1−(1−WWM⁄)1(1+b)⁄(3−3)在初始土湿为W条件下,降雨量PE的产流量可由下式计算:在全流域蓄满前为R=PE−W−WM+WM(1−(PE+a)WMM⁄)b+1(a+PE≤WMM)(3−4)在全流域蓄满后为R=PE−(WM−W)(a+PE𝑊𝑀𝑀)(3−5)为此,需要确定的参数有流域平均最大缺水量WM、流域蓄水容量分布曲线指数b、流域蒸发能力Ep及深层扩散系数C。a、WM的确定WM相当于历史上一次降雨的最大损失量。由一次降雨停止时,能使土壤含水量达到田间持水量时的包气带的水量平衡原理可知:WM=P−E−W0−R(3−6)当W0=0时,WM=P-E-R。因此,选用几次W0很小,P很大,估计雨止时土壤含水量达到最大值的实测洪水资料,及可分析出WM。在南方湿润地区,若用三层蒸发模式计算蒸发,WM可采用140mm(WUM=20mm,WLM=60mm,WDM=60mm)左右。本次课程设计选用WM=140mm。b、b的确定b是个地形地质因素,表示流域蓄水量分布的不均匀性。当全流域蓄满后,这个因素不再起作用,b可取0.3~0.4左右。本次课程设计选用b=0.37。c、蒸发能力Ep的推求流域内若有蒸发皿资料,则可以直接用蒸发皿实测蒸发值乘折减系数K,作为流域的Ep值。K值完全靠优选,当计算的年径流误差最小时,K值就是最优。如原始资料用E80,K值大约为0.85~0.9左右,原始资料用E601,K值大约为0.9~1.0左右。由于密塞流域衢县蒸发站采用80cm套盆式蒸发皿,本次课程设计K选0.86。d、C的确定C可取1/6。(2)计算各年逐日的流域平均降水量,绘制各年逐日流量过程线及流域平均日降水过程线。结果如图1所示。(要求计算1986年3月~10月资料)图11986年逐日流域流量过程和平均日降水过程线(3)计算W0根据水量平衡原理,W0的计算公式为:W0,t+1=W0,t+Pt−Et−Rt(3−7)Rt可由产流计算公式计算;Et采用三层蒸发计算模型计算;W0的起算值,可以以久旱以后(W0=0)或大雨后(W0=WM)开始。03-0103-3104-3005-3006-2907-2908-2809-2710-2701002003004005006007000.020.040.060.080.0100.0120.0140.0降水P(mm)流量Q(m3/s)日期(日-月)表1:密塞流域降水径流计算表(1986年3月至10月)单位:(mm)(WM=140K=0.86b=0.37)日期(月日)各雨量站雨量¯pEwEpEuELEDEP-EaR△WWUWLWdWOP1P2P3P4P5P618/May0.443.140.083.619/May172.8119.4106.5120.5115.577.0118.63.32.82.80.002.8115.893.059.356.420.060.060.0140.020/May34.840.439.537.147.544.840.743.43.40.003.437.2191.837.20.020.060.060.0140.021/May0.300000.10.10.60.50.50.000.5-0.4191.80.0-0.419.660.060.0139.622/May0000000.04.13.53.50.003.5-3.5188.90.0-3.516.060.060.0136.023/May0000000.054.34.30.004.3-4.3177.50.0-4.311.760.060.0131.724/May0000000.04.43.83.80.003.8-3.8167.50.0-3.87.960.060.0127.925/May0000000.03.22.82.80.002.8-2.8159.80.0-2.85.260.060.0125.226/May0000000.04.13.53.50.003.5-3.5154.60.0-3.51.760.060.0121.727/May0000000.05.54.71.73.104.7-4.7148.30.0-4.70.056.960.0116.928/May0.81.12.12.104.41.83.83.31.81.403.2-1.4140.40.0-1.40.055.560.0115.529/May0000000.03.42.90.02.702.7-2.7138.00.0-2.70.052.860.0112.830/May0000000.03.42.90.02.602.6-2.6133.80.0-2.60.050.260.0110.231/May1.82.75.40.18.43.23.64.84.13.60.404.0-0.4129.80.0-0.40.049.860.0109.81/Jun27.127.142.321.162.947.338.03.12.72.70.002.735.3129.214.820.510.360.060.0130.32/Jun0.22.50.49.93.418.85.900.00.00.000.05.9164.53.12.713.060.060.0133.03/Jun0000000.00.20.20.20.000.2-0.2170.30.0-0.212.960.060.0132.94/Jun0000000.04.53.93.90.003.9-3.9169.90.0-3.99.060.060.0129.05/Jun0.20.40.40.20.300.36.55.65.60.005.6-5.3161.80.0-5.33.660.060.0123.66/Jun0000000.02.52.22.20.002.2-2.2151.80.0-2.21.560.060.0121.5日期(月日)各雨量站雨量¯pEwEpEuELEDEP-EaR△WWUWLWdWOP1P2P3P4P5P67/Jun0000000.03.22.81.51.302.8-2.8148.00.0-2.80.058.760.0118.78/Jun0000000.04.13.50.03.503.5-3.5143.30.0-3.50.055.360.0115.39/Jun0.5000.100.80.22.72.30.21.902.2-1.9137.70.0-1.90.053.460.0113.410/Jun7.33.96.24.342.64.72.21.91.90.001.92.8134.71.01.80.055.260.0115.211/Jun17.76.812.55.562.28.521.71.70.001.76.7137.52.64.10.059.360.0119.312/Jun9.49.886.99.87.18.50.80.70.70.000.77.8144.23.34.53.860.060.0123.813/Jun0000000.00.90.80.80.000.8-0.8152.00.0-0.83.060.060.0123.014/Jun0000000.02.62.22.20.002.2-2.2150.70.0-2.20.860.