1、科研与管理水利规划与设计 年第 期：.基于网格新安江模型的汇流方法研究罗玉玲，张甜甜，章玉霞，杜若愚（.河海大学水文水资源学院，江苏 南京；.江苏省徐州市水务局，江苏 徐州；.南京河海大学出版社有限公司，江苏 南京）摘要：为使网格新安江模型更好应用于无资料地区洪水模拟及预报，文章采用网格水滴汇流方法构建网格新安江模型的汇流部分，以湿润地区屯溪流域为研究流域，将网格水滴模型与扩散波汇流模型及马斯京根汇流模型对比验证其模拟精度，并将网格水滴模型应用于半湿润半干旱地区东湾流域验证其适用性。结果表明，采用网格水滴汇流方法的洪水模拟过程在屯溪流域的模拟精度与其他 种方法基本相当，均可以达到甲级预报精度；
2、同时在半湿润半干旱地区东湾流域的应用达到乙级预报精度。网格水滴汇流方法具有一定适用性。关键词：网格新安江模型；网格水滴；汇流演算；分布式水文模型中图分类号：文献标识码：文章编号：（）收稿日期：基金项目：国家自然科学基金（）作者简介：罗玉玲（年），女，硕士研究生。：.无资料地区的洪水预报是目前水文工作者面临的挑战之一。年国际水文学会提出了水文资料匮乏流域的洪水预报 年计划，期间涌现了大量相关研究成果。其中流域水文模型作为研究水文现象的有效工具是解决无资料流域水文预报的代表性方法之一，是无资料地区水文洪水预报的重要解决手段。分布式水文模型因其对降雨和下垫面条件的空间变异性的考虑，可以更好地利用现代
3、化技术条件，是目前水文模型的重要研究方向之一。汇流过程是洪水预报的重要环节，影响最终预报的精度，是水文模型研究的重点之一。目前的分布式水文模型研究，在汇流方法上，主要是基于圣维南方程组及其简化算法的水力学方法如运动波、扩散波等，还有基于水量平衡和槽蓄方程的水文学方法，如线性水库法、马斯京根法。水力学方法具有较强的物理性，虽然可以有效反映汇流过程中的诸如河道断面变化等的影响，但需要较多的河道实测资料，且因其计算复杂，目前尚没有方法求得其精确解析解，在无资料地区运用受限。水文学方法是水力学的简化，计算更加简单，资料需求也相对较少，但在无资料地区的运用仍然受到一定限制。为更好地利用下垫面条件，对无资
4、料地区进行快速洪水预报，本文引入网格水滴汇流方法，进行分布式水文模型的汇流演算，以更综合地考虑下垫面条件和水文特征。网格新安江模型是分布式水文模型，其自创建以来就有着较为广泛的应用。故本文在网格新安江基础上，同时构建采用网格水滴汇流网格新安江模型，将其应用于湿润地区屯溪流域，与原有基于扩散波方法的网格新安江模型及基于马斯京根法的网格新安江模型进行模拟精度对比，验证其可行性，分析网格新安江模型应用网格水滴汇流方法的合理性；并将基于网格水滴的网格新安江模型应用于半湿润半干旱地区东湾流域，验证其适用性。模型与方法.网格水滴汇流方法“网格水滴”是流域上将有降雨并发生产流的网格看作水滴，水滴受到重力和阻
5、力的影响，沿坡度方向运动，其到达出口断面的时间主要取决于汇流路径及汇流速度。.坡面汇流根据 汇流方法确定网格水滴汇流的时间和路径，坡面汇流速度公式采用由美国水土保持局工程手册提出的流速与地表坡度的函数关系式：（）式中，坡面流水流速度，；反应地表糙 年第 期水利规划与设计科研与管理率的经验系数；沿水流方向的地表坡度。.河道汇流河道汇流采用根据曼宁公式演变的河道汇流速度公式：（）式中，渠道流水流速度，；反应糙率的参数；经验指数；网格汇水面积，；河道网格坡度。.汇流时间矩阵根据汇流路径和汇流速度，将各个网格单元所需要的时间累加得到其到达流域出口所需要的时间，即为该网格的汇流时间。（）所有网格的汇流时
6、间组合一起，即为流域的汇流时间矩阵。.模型构建.基于扩散波的网格新安江模型基于扩散波的网格新安江模型（模型）是一个分布式水文模型，以 网格为计算单元，主要分为蒸散发、蓄满产流、分水源、汇流 个计算模块，其中汇流模块采用一维扩散波汇流方法，依次将各种水源演算至流域出口，该方法物理基础较好，可以利用下垫面信息估算参数。.基于逐网格马斯京根法的网格新安江模型基于逐网格马斯京根法的网格新安江模型（模型）其模块区分与 模型一致，其汇流方法采用的是逐网格马斯京根法，依次将地表径流、壤中流、地下径流演算至流域出口。该方法具有计算简便、实用性强的特点。.基于网格水滴的网格新安江模型基于网格水滴的网格新安江模型
7、（模型），产流以蓄满产流原理为基础，计算每个网格单元的产流量，采用自由水蓄水库结构，将其划分为地表径流、壤中流及地下径流 种水源，再根据汇流演算规则演算至出口断面。采用的网格水滴模型进行汇流演算时假设栅格单元是由坡面及河道组成，坡面栅格分为地面径流、地下径流及壤中流，采用网格水滴汇流方法，根据网格水滴汇流速度、汇流路径及时间矩阵，将径流演算至流域出口。研究流域和数据处理.研究流域屯溪流域如图 所示，位于安徽省，该流域靠近中国东南沿海地区，属于亚热带季风气候区、流域内四季分明，气候温和，多年平均气温约为。屯溪流域面积为，地势西高东低，相对高差较大，水流方向自西向东。流域内植被良好，雨量充沛，多年
8、平均降雨为，降水在年内、年际分配不均匀，汛期内的降雨量占全年总降雨量；为典型的湿润流域。图 屯溪流域水系及水文站点图东湾流域（如图 所示）伊河河源地区，流域面积。属于大陆性季风气候，是半湿润地区。降雨量分布不均匀，流域内多年平均降水量为，日降水量一般在 以上，最大值可达。暴雨主要集中在 月，占总降水量的，多年平均蒸发量为，流域地势西高东低，上游地区植被良好，主要为林地，下游河道附近裸地较多。图 东湾流域水系及站点图.数据本文模型所需下垫面数据有 数字高程数据、土壤类型和土地利用数据。数据通过地理空间数据云网站下载，采用 分辨率数据。科研与管理水利规划与设计 年第 期土壤类型数据为世界土壤数据库
9、 中 万土壤栅格数据。土地利用数据为美国地质勘探局网站下载，是分辨率为 的矢量文件。将研究流域数据包括降雨、蒸发及流量枢机们进行 插值。选取屯溪流域内 年 场洪水资料进行模拟，其中 场洪水进行率定，场洪水进行验证；选取东湾流域 年 场洪水资料进行模拟，其中 场洪水进行率定，场洪水进行验证。屯溪流域流量数据来自屯溪水文站，降雨资料来自月谭等 个雨量站（如图 所示）；东湾流域流量数据来自东湾水文站，降雨资料来自谭头等 个雨量站（如图 所示）。结果与讨论.参数率定在应用模型进行屯溪流域洪水模拟时，采用的栅格单元分辨率为 （），网格水滴汇流参数主要参数为 及。其中 个流域的 率定为.，屯溪流域的 值率
10、定为.，东湾流域的 值率定为.。坡面汇流公式及河道汇流公式与流域下垫面输入极其相关，由于在处理 填洼时，对于洼地的处理为将其处理为与附近网格 值相当，如此会造成单元网格的坡度过小导致汇流速度过慢影响最终汇流，为解决此问题，采取措施为坡度小于.时，网格采用流域坡面网格的平均坡度进行速度计算。.时间矩阵分析根据 模型，可以得到屯溪流域时间矩阵。如图 所示，从屯溪流域的时间分布图来看，汇流时间的分布为河道网格的汇流时间小于坡面汇流时间，且相对河道距离近的地方，汇流时间更短（如图 所示）。造成如此汇流时间分布的主要原因是，图 屯溪流域汇流时间分布图图 汇流时间与水系分布图河道的汇流速度比于坡面汇流速度
11、大，水滴在河道汇流所花费时间更短，坡面水滴在汇流过程受到坡度及汇流路径的影响大，部分坡度较小网格，汇流速度较慢，且因其远离河道，坡面汇流路径长，导致汇流时间相应变长。随着坡面水滴流入河道，汇流速度也相应增大，故距离河道近的地方汇流时间相对较短。.模拟结果分析对比依据 水 文 情 报 预 报 规范，选取洪峰相对误差、径流深相对误差、峰现时间误差和确定性系数 个指标对模拟效果进行综合评价。.屯溪流域结果对比根据屯溪水文站实测流量数据，采用 模型、模型及 模型 种不同汇流模型模拟屯溪流域洪水过程的合格率、确定性系数均值的统计结果见表，图 是 种模型洪水模拟结果与实测流量的对比图。表 屯溪流域各模型模
12、拟结果统计洪水期水文模型径流深合格率 洪峰合格率 峰现时间合格率 平均确定性系数率定期.验证期.率定期共有 场洪水，率定期的径流深模拟，模型、模型和 模型模拟相对误差分别介于 .、.及.之间。模型、模型和 模型的洪峰相对误差分别介于 .、.和 .年第 期水利规划与设计科研与管理.；模型的峰现时间误差分别在 ，及 之间。模型、模型和 模型的确定性系数分别介于.、.和.之间。验证期共有 场洪水，验证期的的径流深模拟，模型、模型及 模型的相对误差介于.、.和.；模型、模型和 模型的洪峰相对误差分别介于 .、.和 .之间；种汇流方法的时间合格率都为.，误差介于 、及 之间。模型、模型和 模型的确定性系
13、数范围在.、.及.之间图 屯溪流域 号和 号洪水实测与模拟洪水过程线对比根据 种模型的模拟结果及图 三种模型模拟洪水过程与实测流量对比图，表明 种汇流方法在屯溪流域洪水过程的模拟精度相当，从径流深合格率、洪峰合格率和平均确定性系数上看均达到甲级预报精度，采用这 种方法进行网格新安江的汇流计算，均可以取得较好的应用效果，但 号洪水，种汇流方法的峰现时间都相对提前，方法的峰现时间较为接近实际峰现时间，方法峰现时间较为提前一些，方法居于中间，种汇流方法精度没有较大区别，都有较好的模拟效果。种汇流方法都能很好的反应洪水的走势，对于降雨的反应也较好，能够很好的应用于湿润地区。在物理意义上，模型的参数具有
14、较为明确的物理意义，对于下垫面特征信息有着良好的考虑；模型是水文学方法，对下垫面情况的考虑较少，模型物理意义较为明确，且对下垫面特征信息考虑较多；在计算效率方面，模型计算一次所花费时间比 模型及 模型长，计算效率较低，模型的计算效率与 模型的计算效率相当，均具有较高的计算效率。随着数据信息的增多，如分辨率的提高，模型与 模型及 模型计算时间的差异增大。模型物理意义较为明确，适用性较为广泛，在平原地区也可以较好应用，对下垫面条件对汇流影响考虑较为充分，在无资料地区应用具有一定的适用性，但是难以对无资料地区洪水进行快速预报；模型计算简便，运算效率较高，实用性强，但较为依赖降雨径流的观测资料，且在无
15、资料地区的适用性较弱；模型的参数物理意义较为明确，计算简便，具有较高的运算效率，且其输入项可以根据下垫面特征提取，参数估算较为简单，故对于无资料地区的快速洪水预报具有一定的适用性。.东湾流域结果分析根据东湾流域流量数据，采用 模型对东湾流域洪水过程进行模拟，验证网格水滴汇流模型在半湿润半干旱地区的适用性，模型在东湾流域的模拟结果误差统计见表。从表 可以看出，东湾流域用于率定的 场洪水种，场洪水的径流深误差小于，率定期平均径流深合格率为；模拟的平均洪峰合格率为；平均峰现时间合格率为；平均确定性系数为.；达到了乙级预报精度。率定期的平均径流深合格率为，平均洪峰合格率为，平均峰现时间合格率为，平均确
16、定性系数为.；图 为东湾流域 号洪水的过程，率定期的 场洪水对汇流参数进行了很好的率定，使得率定期和验证期次洪模拟达到了较好的精度，结果表明，网格水滴汇流模型在半湿润半干旱地区的平均模拟精度可以达到乙级标准，即网格水滴模型的汇流方法在半湿润半干旱地区具有一定的适用性。科研与管理水利规划与设计 年第 期表 东湾流域 模型模拟结果误差统计洪水期洪号径流深相对误差 洪峰相对误差 峰现时间误差 确定性系数率定期.验证期.图 东湾流域 号洪水 结语引入网格水滴汇流（模型）的网格新安江模型，将其应用于湿润地区屯溪流域及半湿润半干旱流域东湾流域。研究结果表明，与采用扩散波方法（模型）的网格新安江模型及采用逐网格马斯京根方法（模型）的网格新安江模型相比，模型在计算精度相当的情况下，可以兼具计算效率、实用性及较为明确的物理意义的特点；模型在湿润半湿润地区东湾流域的适用性良好，为解决无资料地区洪水的快速预报问题提供了一个较好的解决思路。参考文献 ，（），：，（）姚成，章玉霞，李致家，等 无资料地区水文模拟及相似性分析 河海大学学报（自然科学版），（）：李治军，王华凡，侯岳，等 我国无资料地区水文预报研究进
