飞来峡水利枢纽基于DEM的分布式流域水文预报模型.docxVIP

第48卷第3期2012年3月Vol.48,No.3

第48卷第3期

2012年3月

Vol.48,No.3Mar.,2012

GANSUWATERRESOURCESANDHYDROPOWERTECHNOLOGY

·水文水资源·

飞来峡水利枢纽基于DEM的分布式流域水文预报模型

廖燕芬

(广东省飞来峡水利枢纽管理处，广东清远511825)

摘要：利用飞来峡水利枢纽建库运行后的水文资料，建立基于DEM的分布式水文预报模型，对参数进行率定，并进行预报模拟。从模拟结果看，模型确定性系数、洪峰相对误差、径流总量相对误差等指标基本在可接受范围内，对大部分场次洪水的预报精度较好，洪峰流量合格率也较好，但整体上对洪峰的模拟值偏低，预报精度有待进一步提高。

今后随着资料的不断延长和积累，对参数进行重新率定，其预报结果可应用于预报作业中。关键词：分布式水文模型；数字高程模型(DEM);洪水预报；飞来峡水利枢纽

中图分类号：P338*.9文献标志码：A文章编号：2095-0144(2012)03-0001-04

飞来峡水利枢纽位于北江干流中下游的清远辖区，上距英德市50km,下距清远市33km,坝址控制流域34097km2,占北江流域面积的73%。枢纽是以防洪为主，兼航运、发电等综合效益，与北江大堤、酒江滞洪区共同组成北江中下游防洪工程体系，通过堤库联合调度，提高广州、佛山、清远等珠三角地区的防洪减灾能力，将北江大堤防洪能力由100年一遇提高到300年一遇。但飞来峡水库纳入连江高道、干流白石窑、渝江长湖等支流来水，库区汇水面积仅70km2,并且水库库区属于河道型水库，动库容的影响相当大，导致洪水预报和水库调度的难度大。如何做好水库的洪水预报，提高预报的精度是关键问题。本文利用实测降水径流资料，基于DEM建立了分布式水文预报模型进行洪水预报，预报结果基本可用。

1水文物理过程概化

基于DEM的网格产汇流计算是流域分布式水文预报方案的核心，它将应用于流域表面的每一个基于DEM的网格单元上。在DEM的每一个栅格上，假设有3种不同的蓄水单元，即河道、地下土壤和地表。

1.1河道

位于栅格单元上的每段河道，其上下游断面处的流量分别记为Q和Qot。每个栅格周围都有，8个栅格，也许只有一个栅格的高程比当前栅格的高程高，也可能有几个，上游流量Qm实际上是比当前栅格i高程更高的各个相邻栅格j上的河道出流量之

1.2地下土壤

地下土壤的蓄水能力用Sm,实际蓄水量用S。假定流域产流机制为蓄满产流，降水P落在地表后会直接进入地下土壤。与地下土壤蓄水量有关的物理过程还包括蒸散发ET,地下水入流QSm,地下水出流QSo与Qm相似，QSm也是各个相邻栅格的地

下水出流量之和，即

1.3地表

对某些位于低洼处的栅格，在某一时刻汇入土壤的水量会超过其缺水量(Sm-S)。假设来水量使土壤蓄满后，剩余部分冒出地面形成浅层地表水。浅层地表水在重力的作用下产生坡面流，逐渐汇入河道。

浅层地表水体积S,其产生的流量为Qp。土壤未蓄满、蓄满后如图1、2所示。

(1)Sp=0

图1(未蓄满)模型物理过程示意

收稿日期：2012-03-29

作者简介：廖燕芬(1978-),女，广东人，工程师，主要从事水利工程运行调度和水文情报预报。

2012年第3期甘肃水利水电技术第48卷

·2

(2)Sp0

图2(蓄满)模型物理过程示意

2栅格单元产流模型

2.1蒸散发

每栅格蒸散发的计算公式为：

(1)

式中：ET——栅格点的蒸散发量；i——栅格的空间位置；

S——栅格点的实际蓄水量；S——地下土壤的蓄水能力；

PE——该点的蒸散发能力，由当地气象资料来估计，计算方法参见相关文献[1。

2.2地下水出流

地下水出流速率不仅与土壤蓄水量有关，也与地下水力学坡度有关。假定水力学坡度可用地表坡度近似，则计算每个栅格的地下水出流的近似公式为：

(2)

式中：QS——栅格点的地下水流出量；

ST地下水出流门限值，可表示为当地土壤蓄水能力S的百分比，一般为20%;

T——时间常数，反映地下水水流特性；p——栅格的平均向下坡度；

b——经验参数，反映坡度对地下水出流的影响，取值在0和1.0之间，建议取值为0.5。

2.3土壤水量平衡方程

对某一栅格下的地下土壤来说，其输入有降水P和地下水入流QSm,而输出包