水文预报课设报告

基本任务

任务一：根据已给数据资料及参数（本报告采用89-90年的历史数据），将流域作为整体（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc。

任务二：根据已给的设计暴雨资料（宝口流域23日至26日暴雨）和任务一率定的Kc，将流域作为整体进行（1）次洪产流量，划分水源.直接径流汇流，地下径流汇流计算；（2）绘出直接径流过程，地下径流过程，总的流量过程。

基本资料

白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期4—9月降雨量占年降雨量的81％左右：径流系数0.5一0.7。

流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。

流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。雨量站分布较均匀．有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一

般2—3h，有时更短；一次洪水总历时2~5d。

计算公式

（1）蒸散发计算

根据流域特点，蒸散发计算采用的是三层蒸散发计算模式。三层蒸发模式的具体计算如下：

1)当WU+P>=E P,

EU=E p,EL=0,ED=0;

2)当WU+P=C.WLM,

EU=WU+P,EL=(E P-EU)*WL/WLM,ED=0;

3)当WU+P

EU=WU+P,EL=C*(E P-EU),ED=0;

4)当WU+P

EU=WU+P,EL=WL,ED=C*(E P-EU)-EL.

（2）产流量计算：

根据流域特点，产流量计算系根据蓄满产流理论得出的。蓄满产流，即任一地点上，土壤含水量达田间持水量前，降雨量全部补充土壤含水量，不产流；当土壤蓄满后，其后续降雨量全部产生径流。

a=WMM(1-(1-W0/WM)^(1/(b+1))

土壤三层含水量WU WL WD 计算

1.PE>0，则产流，即R>0；此时土

壤中上中下层含水量会发生变化，如果

R ΔW

产流，则PE-R后的值就给三层土壤来

分，上层的土壤先补充水，如果上层

WU>WUM，中层开始补充水，如果WL>WLM，

下层开始补充水，如果三层都达到最大

值，则土壤水分保持最大值不变。

2.否则不产流，R=0。即P-E<=0 此时土壤中水分含量随着蒸发减少，上层土壤先开始蒸发，然后中层，最后下层。

产流时，

当PE+a

R=PE+W-WM+WM*(1-(PE+a)/WMM)^(b+1)

当PE+a

R=PE+W-WM

（3）水源划分：

流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄

作用不同，各径流成分在流量过程线上的反应是不一样的。

各水源分量如下：

1）当PE〈=FC时：RD=0.0

RG=R

2）当PE〉FC时：RG=FC*R/PE

RD=R-RG

（4）汇流计算：

根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。

地上产生的径流用单位线推流，地下部分用线性水库。

计算公式：

QS(I)=RS(I)*U

QG(I)=CG*RG(I)*U+(1-CG)*QG(I-1)

Q(I)=QS(I)+QG(I)

上式中RS即为地表径流Rd.

基本数据

任务一：

1.计算流域为白盆珠水库上游的宝口流域，流域面积553 km2。

2.流域内有四个雨量站：禾多布、马山、高潭、宝口，权中重分别为0.55、0.14、0.33、0.20。

3. 资料文件名：87-92data csv，数据格式为：T(i) Q(i) E(i) P1(i) P2(i) P3(i)

P4(i)

4.计算参数：

任务二：

1.计算参数：

WM=140,UM=20,LM=60,DM=60

IM=0.002,FC=11.0

W=110,WU=10,WL=40,WD=60

b=0.3,C=0.16 Cg=0.978 Qg=55.3 2.计算基本数据如下表：

计算结果任务一：kc优选计算表

1991年KC 优选

1992年KC 优选

1991，1992两年KC 优选

由于个人认为KC优选应该以单独一年为单位，因此本程序只是对单独一年的KC 优选做了程序，如果以两年为单位使用范围会很小。因此，两年共同的KC优选，在本程序中的解决方法是，将单独一年KC从0.9到1.3间隔0.01的值全部带入计算并得出所有R与实测R计算相对误差的值输出并可以导入EXCEL，因此可以把两年所有KC和相对误差导入EXCEL。通过绘制图表来反映共同KC时，相对误差之间的关系。以选出尽可能好的KC值。

两年KC优选的图表：