周传争水文预报课程设计报告剖析.docx

周 传 争 水文预报课程设计报告 水文预报课程设计 （文天学院 11 级） 班级 学号 姓名 周 成绩 河海大院 2014年 1月 4日 1 周 传 争 水文预报课程设计报告 一 基本任务 任务一： 根据已给数据资料（本文选取 1995、1996 年的历史数据）、参数，编写程序将流域作为整体： （1）进行日模型产流量计算； （2）将计算年径流与实测年径流进行比较，优选蒸发折算系数 Kc； 任务二 根据已给设计暴雨资料、参数 , 编写程序，将流域作为整体： 进行次洪产流量、划分水源 直接径流汇流、地下径流汇流计算； （3）绘出直接径流过程、地下径流过程、总的流量过程 ( 绘制成图 ) 。 二 基本资料 1）研究流域概况 呈村流域控制面积为 290 km2，地势南高北低，相对高差较大，平均海拔高 程为 583 m，流域河道平均坡度为 0.95%，最大汇流路径长度为 36 km。流域内 植被良好，雨量充沛，多年平均降雨量约为 2100 mm，流域降水在年内年际分配 极不均匀，为典型的湿润流域。 该流域植被类型主要包括常绿针叶林、 落叶阔叶 林、混合林、森林地、林地草原、牧草地与作物地，土壤类型主要为黏壤土。流 域内地势平缓， 土壤主要有黄壤和砂壤， 具有明显的腐殖层， 淀积层和母质土等 层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植 农作物和经济作物为主，植被良好。 2 周 传 争 水文预报课程设计报告 三 基本数据 (1) 计算的呈村流域面积 290 km2 。 (2) 流域内有 10 个雨量站：呈村、董坑呜、棣甸、樟源口 、汪村、用功城、左 龙、冯村、田里、大连 , 各站面积比例分别为： 0.1, 0.1, 0.10, 0.1 ，0.1, 0.1, 0.1, 0.1, 0.1, 0.1 资料放在 / 周传争 Data/ 文件夹下，以年份命名的文件， 数据格式为： T（ i ）、 Q(i) 、E(i) 、P1(i) 、 P2(i) 、P3(i) 、P4(i) 、P5(i) 、P6(i) 、P7(i) 、P8(i) 、 P9(i) 、 P10(i) 。 暴雨资料放在 / 周传争 Data/ 文件夹下 pdata.csv 单位线数据在 / 周传争 Data/ 文件夹下 unitgraph.csv 计算参数如下表 3 周 传 争 水文预报课程设计报告 表一 任务一计算参数 计算年份 参数 初始张力水容量 WM WUM WLM WDM W WU WL WD 120 20 40 60 1995-1996b C IM Fc 90 10 20 60 0.4 0.18 0.01 2 表二 任务二参数表 计算时间 参数 初始张力水容量 WM WUM WLM WDM W WU WL WD 120 20 40 60 1989052020-1989052621 Fc CG QG 98 18 20 30 2 0.989 4 表三 各参数物理意义 层次 参数符号 参数意义 Kc 流域蒸散发折算系数 WUM 上层张力水容量（ mm） 第一层次 蒸散发计算 WLM 下层张力水容量（ mm） WDM 深层张力水容量（ mm） C 深层蒸散发折算系数 WM 流域平均张力水容量（ mm ） 第二层次 产流计算 IM 不透水面积比例 b 流域蓄水容量分布曲线指数 第三层次 水源划分 Fc 流域时段稳定入渗水量（ mm ） 第四层次 汇流计算 CG 地下水消退系数 QG 地下水出流系数 四 计算公式 蒸散发计算 根据流域特点，蒸散发计算采用的是三层蒸散发计算模式。 三层蒸发模式的具体计算如下： （1）当 WU+P≥ Ep 时 EU= Ep EL=0 ED=0 4 周 传 争 水文预报课程设计报告 （2）当 WU+PEp，WL/WLM ≥C 时 EU=WU+P EL=(Ep-EU)* WL/WLM ED=0 （3）当 WU+PEp，C*(EP-EU) ≤WLC*WLM 时 EU=WU+P EL=C*(Ep-EU) ED=0 （4）当 WU+PEp，WLC*(Ep-EU) 时 EU=WU+P EL=WL ED=C*(Ep-EU)-EL 土壤蒸发量： E=EU+EL+ED ( 同时刻相加 ) 产流量计算 根据流域特点，产流量计算系根据蓄满产流理论得出的。蓄满产流，即任 一地点上，土壤含水量达田间持水量前，降雨量全部补充土壤含水量，不产流； 当土壤蓄满后，其后续降雨量全部产生径流。 PE0，产流；否则不产流。 1 (1 WM ) b WMM WMM )dWM WM (1 0 WMM WM b 1 5 周 传 争 水文预报课程设计报告 初始初始土湿分布与计算 a W (1 )dWM 0 W WM[1 (1 a )b 1 ] WMM