传统新安江模型在流域分布式水文模拟中的解决方案.docxVIP

1

传统新安江模型在流域分布式水文模拟中的解决方案

赵海伟1，夏达忠1，张行南1，唐雷彬2，夏希1

1、河海大学水文资源与水利工程科学国家重点实验室，南京，210098

2、辽宁省水文水资源勘测局，沈阳，110003

摘要：大流域特征十分复杂，如何描述和体现流域特征的时空变化一直是水文工作者致力解决的问题。本文结合传统新安江模型的特点，基于DEM数据，借助GIS空间分析功能，提出了采用“流域－产汇流分区－单元面积”的流域空间离散方法，并同时给出了模型参数空间化的方法，以期使该模型能更好地适应流域分布式水文模拟。

关键词：分布式水文模拟，空间离散，参数空间化

1．引言

近几十年以来，伴随空间信息科学、计算机技术、地理信息系统、数字高程模型和遥感技术的发展，为流域水文模拟研究提供了崭新的技术手段，为研制和建立分布式水文物理模型提供了强大的技术支撑，使得分布式水文物理模型成为水文学家研究的前沿热点之一[1]，越来越多的分布式水文物理模型被开发，并广泛应用于水利相关领域，给传统新安江模型带来了挑战。

要使传统新安江模型能适应流域分布式水文模拟，主要需解决两个问题：①如何对流域进行空间离散，以考虑降雨和流域下垫面特性的空间变异性；②如何使得模型参数空间化，以解决无资料或缺资料地区的水文模拟。

2．传统新安江流域水文模型

水文模型是对流域上发生的水文过程进行模拟所建立的数学模型，其实质是对流域上发生的水文过程进行模拟计算[2]。水文模型可分为确定性模型和随机性模型。流域水文模型一般是指确定性模型。从反映水流运动物理规律的科学性和复杂程度而言，水文模型可分为数学物理模型、概念性模型和系统理论模型。从反映水流运动空间变化能力而言，水文模型又可分为集总式模型和分布式模型[1]。目前，具有物理基础的紧密耦合型分布式水文模型是流域水文模拟的重要研究方向，国内外很多学者正致力于这方面的研究，并已成功在水利相关领域得以成功实践。

传统的新安江模型是具有物理概念的分散性模型，其属于概念性模型范畴，又是介于集总模型与分布式模型之间，可称之为“准分布式模型”。其在单元面积上采用集总式，但相对于流域又划分了单元面积。

3．流域空间离散

对于大流域水文模拟来讲，流域内部的地理要素和地理过程存在着极大的空间变异性，气候因子和下垫面因子的空间分布极不均匀，所以需对流域进行空间离散[3]。流域的空间离散中，离散单元可以是规则的，也可以是不规则的。离散单元之间可以是有严格的空间拓扑关系，也可以是相当概化的空间关系[4]。

传统新安江模型一般采用泰森多边形分块的方法，其能较好地考虑降雨的空间变异，但也存在诸多不合理之处。随着数字高程模型数据的普及和GIS技术的发展，使得按自然流域分块的方法更为方便。

2

本文提出了“流域－产汇流分区－单元面积”的空间离散化形式。要实现这种形式，需基于数字高程模型，进行流域流水网的推求，再进行流域单元面积的划分，最终建立流域产汇流分区，从而完成对流域的空间离散。

3.1基于DEM的流域流水网推求

地形对流域汇流过程起着决定性作用。根据DEM数据，对汇流过程进行描述的最好方法是流域流水网[5]。根据水往低处流的公理，每个网格上的径流，将流向相邻网格中相对较低的网格，据此可确定各网格的方向。将各网格方向首尾相连，即为流域流水网。

目前，利用DEM生成流域流水网的方法很多，一般可采用的软件有：由Martz和Garbrecht[6]研制的数字高程流域水系模型(DEDNM)；美国RSI公司开发的RiverTools；美国环境系统研究所(ESRI)开发的ArcGIS软件等等。

本文采用ArcGIS软件作为推求流域流水网的主要工具。图1是基于流域DEM数据,使用ArcGIS软件生成的流域流水网。

当然，任何一种软件和方法导出的流水网不尽相同，可能存在错误或误差。因此，需对流水网进行合理性检验。只有合理可靠的流水网，才能以此为依据，进行空间分析。

3.2流域单元面积划分

在流域水文模拟中，为了考虑降雨及流域下垫面水文特性的空间变异，需将流域按一定形式分成若干单元面积。

(a)原始的DEM

(b)流域流水网

图1西江中下游流域

3.2.1泰森多边形法

传统新安江模型分单元面积的方法是以雨量站为中心的泰森多边形法。其主要是从考虑雨量代表性出发，即尽可能采用直线距离最近的雨量站作为代表站，而没有考虑地形等其它影响因素。由此带来的最主要问题是单元面积很可能跨分水岭。

如图2圆圈处所示，其主要存在两处不合理：①分水岭一边的雨量站来作为分水岭另一边部分面积的代表站，造成面雨