第一章基本任务本次课程设计从2013年7月8日至2012年7月12日，主要任务是对广东省东江一级支流西枝江白盆珠水库的上游宝口流域编制预报方案与产汇流计算。

其基本任务为：任务一：根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较；每人计算两年。

任务二：根据已给设计暴雨资料、参数及做过的习题,自己编写程序，将流域作为整体进行次洪产流量、划分水源、直接径流汇流、地下径流汇流计算；绘出直接径流过程、地下径流过程、总的流量过程。

第二章基本资料2.1 流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积8562。

流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。

暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。

降雨年际间变化大，年内分配不均，多年平均降雨量为1800，实测年最大降雨量为3417，汛期4—9月降雨量占年降雨量的81％左右；径流系数0.5-0.7。

流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。

台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。

流域上游有宝口水文站，流域面积为5532，占白盆珠水库坝址以上集雨面积的64.6％。

白盆珠水库有6年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。

雨量站分布较均匀，有6年逐日降水资料和时段降水资料；宝口水文站具有6年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3小时（h），有时更短；一次洪水总历时2~5天（d）。

2.2 基本资料2.2.1 任务一参数及相关实测资料计算流域为白盆珠水库上游的宝口流域，其流域面积为5532。

该流域内有四个雨量站：禾多布、马山、高潭、宝口，其权重系数分别为：0.33、0.14、0.33、0.20。

本次课设提供的资料为1987年至1992年该流域径流量、蒸发量及各雨量站的日实测资料，其中我计算年份为1987年与1988年，所使用的相关参数值如下表所示：其中蒸散发折算系数需编程优选，优选范围为0.90-1.30。

优选原则为：计算的2年内每年的年径流相对误差尽可能不超过5%。

2.2.2 任务二参数及相关实测资料任务二产流参数为任务一中日模优选的值，为10.03h ，产流参数用日模优选的值，其他参数取任务1中各组的参数，初始张力水蓄量取各组的容量值。

任务二汇流参数为0.968，0=55.3m 3，单位线序号1至11对应的Q （m 3）分别为：0、40、80、130、100、80、48、20、10、5、02.3 产流方式论证宝口流域地处我国南方湿润地区，气候暖热，雨量充沛，多年平均降雨量为1800>1000，年径流系数在0.5—0.7之间，大于0.4；流域内土质疏松，植被良好，不易超渗；一次洪水的流量过程陡涨缓落，持续时间2~5d 左右。

从流域的气象条件，下垫面条件和流量过程的分析知，该流域降雨径流关系具有蓄满产流的特点，可以按蓄满产流建立产流量预报方案。

第三章 计算公式() () () BCW()() () () 14020600.2 0.16 0.001 110104060蒸散发月日时(mm)禾多布马山高潭宝口92312 1.3 6.29.921.617.315 1.37.61620.612.618 1.3 6.2 6.414.915.921 1.38.817.229.418.52424 1.22534.835.324.630.929.929.243.937.860.938.624.846.93390.9 6.97.5 6.112.3120.928.329.934.228.5150.925.642.739.875.4180.993.9137.612413.2210.985.390.88575.925240.851.547.449.238.53 1.139.870.342.197.76 1.143.247.361.545.99 1.120.513.315.813.112 1.110.58 1.8 3.315 1.17.48.47.610.918 1.1 1.8 2.8 2.1 4.621 1.10.200.302624 1.200003 2.100006 2.100009 2.1000012 2.1000015 2.1000018 2.10000212000日 期降 雨 量(mm)3.1 降雨量计算由已知资料知，该地区雨量站分布均匀，且宝口以上有四个雨量站：禾多布、马山、高潭、宝口，其权重系数分别为：0.33、0.14、0.33、0.20，所以宝口流域平均降雨量0.33*P1+0.14*P2+0.33*P3+0.20*P4。

3.2 蒸散发量计算蒸散发计算采用三层蒸发计算模式，即： 上层蒸发量： 下层蒸发量：* 深层蒸发量：* 总蒸发量：式中：为流域蒸发能力（）；为下层土壤含水量（）；为下层土壤含水量（）；C 为蒸发扩散系数。

三层蒸发模式按照先上层后下层的次序，具体计算为： 1)当>时，002)当<，>*时，()*03)当<, C*()<<C*时，*()0 4)当<, <C*()时，*()3.3 降雨产流量计算1）a ＋≤时，＋W －＋（1－WMMa PE +)12）a ＋＞时，－（－W ）3.4 水源划分计算通过稳渗率可划分产流中的直接径流和地下径流。

次洪的各水源分量为：∑∑∑>≤>-=+=FCPE i ii c i FCPE i i FCPE i iici i i PE r f PE rPE r f RG )(RS 总直接径流：总地下径流：3.5 汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。

其计算公式：直接径流过程：(i)(i)*U地下径流过程：(i)*(1)+(1)*(i)*U 总的流量过程：Q(i)(i)(i)第四章计算结果4.1 任务一结果经过程序优选，0.98年份实测R（）计算R（）绝对误差（）相对误差（%）1987 1080 1100 19.3 1.791988 824 813 -10.2 -1.23 4.2 任务二结果t(i)P(i)R(i)(i)(i)(i)(i)Q(i) 2004/9/23 12:001412.7 2.710069.969.9 2004/9/23 15:0014.112.9 2.91010.884.194.9 2004/9/23 18:00119.709.733.297.3130 2004/9/23 21:0018.717.47.41058.3111169 2004/9/24 0:0029.728.518.51094.3123218 2004/9/24 3:003635.125.1101841363202004/9/24 6:00 38.3 37.5 27.5 10 381 148 529 2004/9/24 9:00 7.8 6.9 0 6.9 645 154 799 2004/9/24 12:00 30.5 29.6 19.6 10 799 166 965 2004/9/24 15:00 42.6 41.7 31.7 10 874 177 **** ****/9/24 18:00 93.8 92.9 82.9 10 864 188 **** ****/9/24 21:00 84.1 83.3 73.3 10 1220 198 **** ****/9/25 0:00 47.6 46.8 36.8 10 1770 208 1980 2004/9/25 3:00 56.4 55.3 45.3 10 2370 218 2590 2004/9/25 6:00 50.4 49.3 39.3 10 2650 227 2870 2004/9/25 9:00 16.5 15.4 5.4 10 2600 236 2840 2004/9/25 12:00 5.8 4.8 0 4.8 2360 237 2600 2004/9/25 15:00 8.3 7.2 0 7.2 1860 241 2100 2004/9/25 18:00 2.6 1.5 0 1.5 1250 236 1480 2004/9/25 21:00 0.2 0 0 0 774 228 1000 2004/9/26 0:00 0 0 0 0 396 221 617 2004/9/26 3:00 0 0 0 0 168 214 382 2004/9/26 6:00 0 0 0 0 72.8 207 280 2004/9/26 9:00 0 0 0 0 25.0 200 225 2004/9/26 12:00 0 0 0 0 2.70 194 197 2004/9/26 15:00 0 0 0 0 0.00 188 188 2004/9/26 18:00 0 0 0 0 0.00 182 182 2004/9/26 21:000 0 0 0 0.00 176 1765001000150020002500300035002004/9/23 0:002004/9/2312:002004/9/24 0:002004/9/2412:002004/9/25 0:002004/9/2512:002004/9/26 0:002004/9/2612:002004/9/27 0:002004/9/2712:00Rs Rg R第五章 误差统计与分析5.1 精度评定从计算结果可见，年产流量绝对误差均小于100，所产流量的相对误差均小于5%。

精度统计表明，率定的模型参数是基本合理的。

但由于课设时间限制与任务要求，每位同学只用两年的资料来率定，不满足《水文情报预报规范》中规定：洪水预报方案要求使用样本数量不少于10年的水文气象资料，其中应包括大、中、小水各种代表性年份，并保证有足够代表性的场次洪水资料。

显然0.98的结果还是存在一定问题的。

5.2 误差分析影响流域降雨径流过程的因素很多，利用蓄满产流新安江模型的结构与参数能够很好反映湿润地区降雨径流过程的主要规律和特点，因而能获得较好的精度。

但是模型本身以及模型计算中有很多概化，会造成误差。

造成宝口流域产汇流计算方案误差来源主要有以下几个方面：1.资料代表性的影响。

我们每位同学只用两年实测资料来率定，资料必不能满足时期要求及代表性要求。

实际操作中，应要求有12年以上的连续的对未来有代表性的实测资料，其中10年为率定期，2年为检验期，且这种资料应具备丰、平、枯水年代表性，资料系列前后一致，受人类影响较小。