水文预报课程设计报告书.docx

1、水文预报课程设计报告书水文预报课程设计 姓 名： 学 号： 学 院：水利与环境学院 班 级：指导老师： 时 间：2013.1.23水文预报课程设计说明书1、设计目的1、流域水文模型的用途：洪水预报方案是现代实时洪水预报调度系统的核心部分，是提高预报精度和增长预见期的关键技术。对水资源可持续利用：流域水文模型是水资源评价、开发、利用和管理的理论基础。对水环境和生态系统保护：流域水文模型是构建面污染模型和生态评价模型的主要平台。流域水文模型还是分析研究气候变化和人类活动对洪水、水资源和水环境影响的有效工具。本次课程设计的目的是通过一个具体的降雨径流预报方案的制作，使学生了解生产单位对预报任务的要求

2、。2、通过课程设计，要求掌握如下容：1) 流域综合退水曲线、地下水退水曲线的制作以及次洪分割方法；2) 熟悉降雨径流相关图编制的完整过程；3) 新安江两水源模型结构及产流参数率定方法；4) 流域经验单位线的推求方法；5) 洪水预报方案精度评定方法；6) 利用预报方案进行实行洪水预报方法；7) 利用马斯京根分段连续演算法进行长河段洪水演进预报。2、设计基本资料该流域集水面积1884.6km2，干流河长约273km。流域气候温湿，年降雨量在1700毫左右。地下水位较高，且随季节变幅小，因次，一般情况下，土壤含水量较大。根据流域自然地理条件情况和气候条件，以及洪水流量过程线分析，可知流域产流规律符合

3、湿润地区的蓄满产流模型特征。采用降雨径流相关图制作流域蓄满产流方案，用二层蒸发模型计算蒸发，水源划分考虑两水源划分。3、课程设计资料1. 一场历史洪水的流量过程，相应的各雨量站（3个雨量站）时段雨量与权重（时段长为3h）2. 洪水的前期日降雨量，日最大蒸发量资料（历史洪水与实时预报洪水）3. XX场洪水的退税过程4. 部分场次洪水降雨、径流特征值成果表5. 干流河段的马斯京根参数及分段数（3段）4、产流计算4.1 面平均雨量计算4.11 计算方法泰森多边形法 (1)其中为i站的面积权重，为i站的实测降雨4.12 各测站基本资料如表1、表2所示表1 各测站降雨量 单位:mmA站降雨B站降雨C站降

4、雨18.181818188.51063829812.9032258190.9090909142.5531914964.51612903168.181818278.72340426119.354838745.4545454521.2765957432.2580645236.3636363617.021276625.80645161表2 各测站权重A站权重B站权重C站权重0.220.470.314.13 计算步骤将表1、表2中数据代入公式(1)中并使其结果除以3即得面平均降雨量如表3所示表3 面平均雨量 单位:mm流域面雨量4.0020.0037.0010.008.00再由实测洪水过程可画出流量过

5、程线如图1图1 实测洪水流量过程线4.2 次洪分割4.21 计算流域平均消退系数Cgg由最小二乘法公式 (2)算得Cg=0.837(从流量开始消退时算起)同理可算得地下水消退系数Cgg=0.746(从退水开始算起)4.22 划分水源由地下水退水过程线不变的规律，结合实测洪水过程(无后续洪水)和实测洪水过程可画出图2图2 分割复式洪水再由洪峰(如此题中的683.91)与Cg的乘积的递推，算出一系列的Q值，用径流深计算公式可累加出一系列的R值，即得蓄泻关系曲线，如图3图3 蓄泻关系曲线用图3可查出此场洪水对应的值，则由公式 (2)算得再算得 另附以下几图表表4 地面与地下流量过程 单位:mmQQg

6、Qs1.75711.757101.95184.23320.857217.35716.709410.6477117.00039.1856107.8146361.847311.6612350.1856620.159514.1379606.0216683.912416.6140667.2984666.279219.0902647.1887558.234321.5663536.6681430.399124.0425406.3565315.250426.5187288.7317235.767328.9948206.7724168.251531.4710136.7805112.756833.947178.

7、809662.117136.423325.693838.899438.89940图5 流域平均退水曲线图6 地下水退水曲线表5 蓄泻关系曲线表 单位:m3/sQR0.000688910.000002360.000722620.000007210.000757980.000012290.000795080.000017620.000833990.000023210.00087480.000029070.000917620.000035220.000962520.000041680.001009630.00004844图7 相邻时段流量关系图4.3 流域初始土壤蓄水量的计算4.31 前期影响雨量P

8、a的计算 (3) (4)4.32 计算实例已知资料见表6表6 Pa基本数据 单位:mm日降雨量Pa蒸发能力WMPa,t+11.4609.78053.95525047.413185.946.848951.942.838140.738.259144.4 降雨径流相关图制作及误差检验4.41 计算原理：（1）流域蓄满产流模型计算公式 (5) (6)（2）蒸发计算（两层蒸发模式） (7) (8)表7 降雨径流相关表 单位:mmW0010203038.2591440506063.17080P a010.163820.69567731.674741.1599843.2165655.50625868.879

9、81773.342584.11941106.400000-7.1E-15-7.1E-1500-7.1E-157.11E-1502.78750.012120.103810.20606710.322230.432390.4577020.62213490.8359370.9190621.1589142.7875100.158460.495780.8730911.303471.7136951.8082882.42648653.24322833.5664124.5281311015.750.398260.940891.5495262.246352.9136843.0680824.08308145.446

10、48165.9967617.70498715.75200.648621.349052.13641353.040643.9100634.1117975.4449117.2638028.01305210.482462021.5750.757641.517962.37338513.356974.3042244.5242755.98158077.98347738.81600811.6763521.575301.49582.590713.82905445.263536.6590536.9856889.180772212.3585313.812522030402.730814.259366.0016437

11、8.0433710.0633510.5425213.861032023.13040504.392116.404898.723036711.485214.2934614.9761820.1387943033.1405052.54.879547.023189.499904412.46715.5139616.2622122.532.535.642.552.557.36255.916068.3277811.13416914.539718.1303819.0401427.362537.362540.462547.362557.3625606.529359.0949312.09388515.764819.

12、7178620.75181304043.1506066.158.1081811.060314.55640918.958624.4091426.1536.1546.1549.2556.1566.15709.2097812.427916.28005121.271328.2591430405053.1607072.937510.113413.549817.70562723.286531.1966432.937542.937552.937556.037562.937572.937581.2512.998217.15122.42526931.2539.5091441.2551.2561.2564.357

13、1.2581.258012.531116.564521.6330863038.2591440506063.170809016.652721.8387304048.2591450607073.1809010021.903230405058.2591460708083.190100110.000 30.000 40.000 50.000 60.000 68.259 70.000 80.000 90.000 93.100 100.000 110.000 120.000 40.000 50.000 60.000 70.000 78.259 80.000 90.000 100.000 103.100 1

14、10.000 120130.000 50.000 60.000 70.000 80.000 88.259 90.000 100.000 110.000 113.100 120.000 130.000 140.000 60.000 70.000 80.000 90.000 98.259 100.000 110.000 120.000 123.100 130.000 140.000 图8 降雨径流相关图4.4.2 误差检验用其他场次PER关系查降雨径流相关图，校核是否合格并算出合格率与等级，见表8表8 径流深校核 单位：mmPa有效降雨PE次洪径流深演算径流深许可误差实际误差合格否合格率等级15.

15、245.710.369738246.45286893433.916869311否0.88甲级48.858.229.95190859275.9903817192.951908594是43.925.22.4062514166.1604104853-3.75415907否4681.546.9607386347.59.392147726-0.53926137是40.1113.371.1697806671.5591429514.23395613-0.3893623是45.36631.9997257731.36.3999451530.699725765是35.764.727.5650327120.7011

16、18335.5130065426.863914379否1071.910.5153988413.145397783-2.62999893是29.876.231.84636564266.3692731295.846365644是44.4114.277.2041834578.615.44083669-1.39581655是44.981.944.844248446.88.96884968-1.9557516是5165.938.7288917736.97.7457783531.828891765是80101.7103.827372101.7202.127372048是10.737.64.76295699

17、73.92356886230.839388135是52.4100.568.4976318472.913.69952637-4.40236816是42.774.735.3260480737.47.065209614-2.07395193是39.6150106.9960346109.620-2.60396537是27.860.513.5689830515.106506262.71379661-1.53752321是3894.950.5537405352.910.11074811-2.34625947是52.485.158.7773674857.511.75547351.277367482是47.9

18、46.719.7176971316.683075653.9435394273.034621487是11.5109.440.147997940.98.029599579-0.7520021是56.4104.57958512.91687690631.662708194是44.718.96.6571116344.38046924632.276642387是5、汇流计算5.1、fc的试算5.1.1试算方法 (9) (10)5.1.2 试算实例(1) 基本数据如表9表9 单位:mm 流域面雨(mm)EpPE累计PE累积净雨量R时段净雨量riri/PE4.001.21252.792.790.43240.4

19、3240.1551174120.001.212518.7921.584.30423.87180.2060856437.001.212535.7957.3618.130413.82620.3863404710.001.21258.7966.1524.40916.27880.714510338.001.21256.7972.9431.19666.78751再假设，则有故划分时段净雨如表10所示表10 划分地面及地下净流深 单位:mmRgRs0.30 0.13 0.40 3.47 0.76 13.07 1.40 4.88 1.96 4.83 5.2 计算单位线(试错法)5.2.1 演算实例如表11所

20、示表11 科林法试算表 单位:m3/s时段时段净雨ri(mm)实测地面径流量Qs(m3/s)0.133.474.884.83各部分之和均值13.07试算的UH(mm)假定的UH(mm)平均的UH(mm)调整的UH(mm)(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)00.1300.000.000.000.0000.00013.470.8570.450.000.450.220.630.4735.0017.7335.158213.0710.6482.1712.1414.317.383.272.41169.0085.70169.86734.88107.

21、8154.0258.610.0062.6438.4769.3451.10313.00182.05270.7644.83350.3.45108.5517.080.00129.0895.86254.33187.42268.00227.71240.2975606.2.9292.9582.4616.90195.22162.15443.87327.10227.00277.05227.1666667.2982.4478.73152.7281.61315.50255.36411.94303.57190.00246.78185.9117647.1891.7065.90130.76151.14349.50332

22、.50314.69231.90132.00181.95132.8238536.6681.0545.78110.76129.41287.00318.25218.42160.9682.00121.4882.9406.3570.8628.4492.71109.61231.62259.31147.05108.3667.0087.6867.72710288.7320.0023.2464.4191.75179.39205.5083.2361.330.0030.6744.74911206.7720.000.0040.0163.74103.75141.5765.2048.050.0024.02012136.7

23、810.0032.6939.6072.2988.0248.7635.9317.971378.8100.0032.3532.3552.3226.4919.529.761425.6940.000.000.0016.189.527.013.511500.000.000.000.000.000.00再用调整的单位线计算地面流量，比较预测与实际值间的误差大小。5.2.2 用计算出的单位线推出地面流量如表12所示表12 用计算出的单位线推出地面流量 净雨量ri(mm)调整后单位线UH(mm)0.1285543.46816513.069424.8792174.828756Qs(m3/s)总径流Q(m3/s)

24、0.12855400003.46816535.1580.451972035.6099735.6099713.06942169.8672.18371612.193370184.2441184.24414.879217270.763.48074158.9126845.949460379.1029Qg10379.10294.828756240.2973.08912693.90404222.006317.154350576.4508Qg211.23017587.681227.1662.92032183.33896353.867582.8817916.97694767.1516Qg323.30145790.453185.9112.38996978.78492314.0542132.109782.02463795.2744Qg445.36054840.6349132.8231.70749964.477296.8927117.2461130.7434743.8897Qg585.58246829.472282.1.05509846.0