分布式水文模型EasyDHM在海河阜平流域的应用.docVIP

0 引言 水文模型以水文系统为研究对象, 根据降雨和径流在自然界的运动规律建立数学模型, 通过计算机快速分析、数值模拟、图像显示和实时预测各种水体的存在、循环和分布及物理和化学特性[1] 。分布式流域水文模型最基本的特征是根据流域各处气候信息和下垫面特性要素信息的不同, 将流域划分为若干小单元; 在每一个单元上用一组参数反映其流域特征,具有从机理上考虑降雨和下垫面条件空间分布不均匀对流域降雨径流形成影响的功能。 分布式水文模型的模拟领域涵盖了地表水与地下水计算、水资源数量与质量的联合评价、洪水预报、土地覆盖与土地利用变化对水文过程的影响、生态需水、农田灌溉与农业节水等多方面的模拟计算, 并且可以通过尺度转换与大气环流模式耦合来预测全球变化对水文水资源的影响, 从而纳入全球变化水文研究前沿。由此可见, 分布式水文模型同时具有理论上的前沿性和应用上的广泛性。分布式水文模型在实践中的应用主要包括: 水资源评价; 洪水预报; 干旱评估; 土壤侵蚀及水沙迁移; 水源污染影响土地利用变化的影响; 水生态环境演变; 气候变化影响; 水利工程影响等等[2] 。 本文利用雷晓辉等最新开发的分布式水文模型EasyDHM 模型[3-5] 对阜平流域进行水文模拟, 为下一步对该流域进行洪水预报及水资源评价工作奠定基础。 1 EasyDHM 模型概况 随着水文科学、遥感技术和计算机技术的发展, 从技术可能性上讲, 分布式水文模型已经逐渐可以满足实际水管理业务中的各类需求。为了推进分布式水文模型技术在实际水资源管理与洪水管理中的应用, 雷晓辉等于2009 年自主开发了一个分布式水文模型(Easy DHM Easy Distributed Hydrological Model) 。该模型具有以下几个主要特点。 在水平尺度上, Easy DHM 模型采用子流域内部单元的水平结构; 采用模块化编程的思想, 既保证了模型代码的可读性,也实现了模型的可扩展性; 同时支持多种产流计算算法; 同时支持多种汇流计算算法; 支持全局参数敏感性分析及参数自动率定算法; 提出了参数分区、计算分区的概念。 提供了一套完整的软件平台, 包括了4 大部分: 自主开发的面向水文分析的通用GIS 软件; 前处理模块; 模型计算模块; 后处理模块。 Wetspa 产流算法是EasyDHM 模型采用的产流算法之一。WetSpa 模型是De Smedt 等[ 3] 开发的基于GIS 的分布式水文模型, 可用于流域尺度上洪水预报及水量平衡计算。 EasyDHM 模型在Wet Spa 模型基础上进行改进, 简化了WetSpa 模型的产流参数推求方法。WetSpa 模型的产流概化方法, 在垂向上划分了4 层: 植被冠层、地表层、土壤层和地下水含水层。降雨先进入植被冠层, 发生冠层截留, 穿过冠层的水分会进入地表层, 超渗的水量发生地表填洼, 其他水分入渗进入土壤, 土壤水渗漏则进入地下水含水层。其中地表填洼形成地表径流, 土壤水侧向流出形成壤中流, 地下水回归形成地下径流, 三部分径流之和即为总径流。而植被冠层、地表填洼、土壤水和地下水都会形成蒸发。其中各参数推求算法参见文献[4] 。本研究采用马斯京干法对河道汇流进行模拟。 2 研究区概况 阜平流域地处河北省西部、太行山东麓,位于东经113b45c 至114b32c ,北纬38b39c至39b8c 之间(见图1) 。流域总面积约2 239 km2,其中森林覆盖率达到31. 3%, 植被覆盖率为70%。 图1 阜平流域位置及数字高程模型 Fig. 1 T he locat ion and DEM of Fu ping Basin 阜平流域为全山区流域, 地势西北高东南低, 境内山高坡陡, 沟深谷长, 素有/ 九山半水半分田0 之称, 从大的地形地貌上可分为亚高山、中山、低山和丘陵。2000 m以上的山峰4座, 最高山峰( 歪头山) 海拔2 286. 2 m; 最低海拔200 m, 地貌差异大。山地相对高差200~ 500 m, 山地坡度一般大于25b。 阜平属暖温带半湿润地区, 气候为北温带大陆性季风气候, 四季分明, 主要特点为冬春干旱少雨, 夏季多暴雨, 光热充足, 雨热同季。年平均绝对湿度9. 7mbar, 相对湿度52%。年平均气温12. 7 e ,无霜期140~ 190 d, 降雨量350~ 790 mm, 年平均日照2 801.3 h, 年平均光辐射量136. 6kcal/ cm,可以满足各种农作物及树木生长发育对光和辐射量的需要。 3 数据准备 3. 1 数字高程模型DEM 阜平流域的数字高程模型(DEM) 从全国数字高程模型中提取, 其栅格大小为1 km@1 km(见图1) , 并对其进行河网修正、填洼、流向计算、累积量计算、坡度计