1、1WEP 模型的开发与分布式流域水循环模拟贾仰文中国水利水电科学研究院水资源研究所摘要: 本文介绍了 WEP 模型的开发和验证情况，论述了分布式流域水循环模拟的意义和作用。在分析分布式流域水循环模拟的研究状况及面临问题的基础上，考虑到我国流域尺度大、人类活动影响深、环境复杂多变的实际情况，作者认为，基于物理概念和变时空步长的、模拟对象包括水循环的各要素过程而不仅是产汇流过程的分布式流域水循环模型将是未来的发展方向。作者主张，应加快开发适应我国自然地理特征与气候特点的各类基于 GIS 的耦合式应用系统，为流域水资源管理、洪水预报调度、水环境评价、水土流失监督治理及水生态环境分析等提供强力支持。关
2、键词: 流域、水循环、水文、分布式、模型、WEP一、分布式流域水循环模拟的意义与作用地球环境变化和人类活动的影响改变了水的自然循环规律, 加剧了我国水资源的供需矛盾，许多地区出现了水环境与水生态恶化的严重局势。地表水、地下水及人工侧支循环水等各类水资源转化频繁，狭义的水资源概念与传统的水资源评价方法已显不适。20 世纪 80 年代中期以来，随着计算机技术、地理信息系统和遥感技术的发展，从水循环过程的物理机制入手并考虑水文变量的空间变异性问题，即分布式流域水文（水循环）模型或称“白箱”模型的研究在国内外受到广泛重视，涌现出许多分布式或半分布式模型,如SHE 模型、IHDM 模型及 TOPMODE
3、L 模型等（参见文献 1） 。另外，全球大循环（GCM）研究对陆地地表过程模拟提出了越来越高的要求，土壤-植物-大气连续体（SPAC）研究受到重视，出现了各类 SVATS（土壤-植物-大气通量交换方案）模型，从另一方面加强了水循环的研究。本文使用“流域水循环模拟”而不是“流域水文模拟” ，意在强调需要将流域水循环系统的所有要素过程联系起来研究而不仅仅是产汇流模拟。分布式流域水循环模拟能够回答水在时空间上如何移动和转化、什么样的工程与管理措施才能减少无效耗水以及人与生态如何分水等问题，而且其模型参数具有物理意义、可根据测量和下垫面条件进行推算。因此，基于物理机制的分布式流域水循环模型的研究与开发
4、具有重要意义，在以下几个方面具有不可替代的作用：（1）预测未来环境变化下的流域水资源演变趋势， （2）分析人类活动的影响与各类对策的效果， （3）借助各类遥测技术在缺乏地面观测资料流域进行水文分析与预测， （4）为流域水资源评价与配置、洪水预报调度、水环境评价、水土流失监督治理及水生态环境分析等各专业应用提供强力支持。二、WEP 模型的开发与验证本文作者从 1995 年起从事分布式流域水循环模拟研究，开发了网格分布式流域水循环模型 WEP (Water and Energy transfer Process) 模型（参见文献 2 至 4） 。该模型以长方形或正方形网格为计算单元，便于使用 GI
5、S 和卫星遥感数据，并具有物理概念强、计算精度高和速度快等特点，已在日本谷田川等多个流域得到验证，正在日本战略性创造研究推进事2业项目(CREST)“都市生态圈、大气圈和水圈中的水量能量交换”课题中使用，并正在我国的几个流域进行验证中。WEP 模型 2002 年 10 月获日本国著作权登录，并可从互联网上下载，详见 http:/www.pwri.go.jp/team/suiri/yata-r/index_e.html。虽然 WEP 模型还包括水质模拟模块，受篇幅所限，这里仅就 WEP 模型的水循环模拟部分的开发与验证情况做简要介绍。11 WEP 模型的开发为提高计算效率，WEP 模型对非饱和土
6、壤水运动的模拟采取了比 SHE 模型简化的算法，但强化了对植物生态耗水与热输送过程的模拟，对水热输送各过程的描述大都是基于物理概念。（1）模型结构。各网格单元的铅直方向结构如图-1（a）所示。从上到下包括植被或建筑物截留层、地表洼地储留层、土壤表层、过渡带层、浅层地下水层和深层地下水层等。状态变量包括植被截留量、洼地储留量、土壤含水率、地表温度、过渡带层储水量、地下水位及河道水位等。主要参数包括植被最大截留深、土壤渗透系数、土壤水分吸力特征曲线参数、地下水透水系数和产水系数、河床的透水系数及坡面和河道的糙率等。为考虑网格内土地利用的不均匀性，采用了“马赛克”法即把网格内的土地归成数类，分别计算
7、各类土地类型的地表面水热通量，取其面积平均值为网格单元的地表面水热通量。土地利用首先分为水域、裸地植被域、不透水域三大类。裸地植被域又分为裸地、草地与耕地、树木类、不透水域分为都市地表面与都市建筑物。另外，为反映表层土壤的含水率随深度的变化和便于描述土壤蒸发、草或作物根系吸水和树木根系吸水，将裸地植被域的表层土壤分割成 3 层。（a） （b）SurfacerunofSoil-VegetationGroupWater BodyGroup ImperviousArea GroupWater useleakagePumpedgroundwaterRechargeDepresionLayerTop S
8、oil Layer2nd Soil Layer3rd Soil LayerTransition LayerSubsurfacerunofInterceptionLayerSuctiondifsiTranspirationInfiltrationEvaporation Anthropogenicenergy sourcePrecipitationheatfluxesConfined Aquifer 1 Unconfined AquiferFlow in FlowoutPercolation 1Aquitard 2Confined Aquifer 2 Percolation 2Exchange b
9、ysepaeoroutflowFlow inFlow in FlowoutFlowoutShortwaveradiation LongwaveradiationheatfluxesAquitard 11D River FlowMain River Tributary1D Overland Flow图-1 WEP 模型的结构：（a）网格单元内的铅直方向结构， （b）平面结构WEP 模型的平面结构如图-1(b)所示。首先，为追迹计算坡面径流，根据流域数字高程（DEM）及数字化实际河道等，设定网格单元的汇流方向（落水线） 。然后，将坡面径流沿着落水线用维运动波法由流域的最上游端追迹计算至最下游端。关
10、于各支流及干流的河道汇流计算，视有无下游边界条件采用维运动波法或动力波法由上游端至下游端追迹计算。地3下水流动采用多层模型进行数值解析，并考虑其与地表水、土壤水及河道水的水量交换。(2) 水循环过程的模拟。蒸发蒸腾包括植被截留蒸发、土壤蒸发、水面蒸发和植被蒸腾等。WEP 模型按照土壤-植被-大气通量交换方法(SVATS)、采用 Penman-Monteith 公式详细计算了蒸发蒸腾。由于蒸发蒸腾过程和能量交换过程客观上融为一体，地表附近的辐射、潜热、显热、热传导及地表温度的计算不可缺少。为减轻计算负担，热传导及地表温度的计算采用了强制复原法(FRM)。GREEN-AMPT 入渗模型物理概念明确
11、，所用参数可由土壤物理特性推出，并已得到大量应用验证，因此，WEP 模型采用 GREEN-AMPT 铅直一维入渗模型模拟降雨入渗及超渗坡面径流。GREEN-AMPT 模型仅适用于降雨入渗过程。而非降雨期的表层土壤(通常是非饱和状态)水分量的再分配将影响到降雨入渗时的初期水分量、土壤和植被的蒸发蒸腾和对浅层地下水的补给等，为减轻计算负担，WEP 模型将表层土壤分成数层，按照非饱和状态的达西定律和连续方程进行计算。 在山地丘陵等地形起伏地区，同时考虑坡向壤中径流及土壤渗透系数的各向变异性。地下水流动采用多层模型进行数值解析。浅层地下水运动按照BOUSINESSQ 方程进行二维数值计算，源项包括表层
12、土壤的降雨补给、地下水取水、深层渗漏及地下水溢出(或来自河流的补给)等。在河流下部及周围，河流水和地下水的相互补给量根据其水位差与河床材料的特性等按达西定律计算。为考虑包气带层过厚可能造成的地下水补给滞后问题，在表层土壤与浅层地下水之间设一过渡层，用储流函数法处理。另外，WEP 还考虑了雨水人工储留渗透设施的模拟、防灾调节池的计算及水田的模型化等。22 WEP 模型的验证WEP 模型先后在日本东京的多摩川中部流域（578 km2） 、千叶县海老川流域（27 km2）及茨城县谷田川流域（166 km2）得到验证和应用。其中，海老川流域是高度都市化的流域，谷田川流域是农地与人工林地为主的自然流域，
13、多摩川中部流域是半都市化半自然的流域。WEP 模型的模拟结果示例见图2 至 4。可以看出，WEP 模型不仅对流量，而且对地下水位及土壤水分等均有良好的模拟结果。验证后的 WEP 模型曾用来分析都市化对东京都水热收支及水热通量的空间分布的影响，评价雨水人工储留渗透设施和防灾调节池对流域水循环的改善作用，研究水田的维持河川枯水流量及滞洪效果等（参见文献 2 至 4） 。WEP 模型具有较高的计算效率。以谷田川流域的计算为例，共有 16661 个计算网格单元，计算时段步长采用小时，在 CPU 为 1.4GHZ 的微机上，一年的计算时间约为小时。05101520253035400/10/1 0/1/3
14、0 01/01/30 01/04/1 01/06/1 01/08/1 01/10/1 01/1/30Year/Month/DayHourlyDiscahrge （m3/s)02040608010120140160HourlyRain（m）Rain (PWRI)Observed at KojirabashiSimulated at Kojirabashi图 WEP 模型的流量模拟结果示例（谷田川流域）40369121995-1-1 1995-5-1 1995-8-29 1995-12-27 1996-4-25 1996-8-23 1996-12-21Daily rain (mm)05010015
15、0200Groundwater level (m)Rain Observed Simulated图 WEP 模型的地下水位模拟结果示例（海老川流域）0204060801001 25 49 73 97 121 145 169 193 217 241Hours from Nov.26 1:00 AM, 1997Saturation ofsoilmoisture（） 0510152025303540 Hourly rain（mm）Simulated(top layer averaged)Observed(10cm below surface)图 WEP 模型的土壤水分模拟结果示例（海老川流域）三、分
16、布式流域水循环模拟面临的难题与对策分布式流域水循环模拟在我国推广应用所面临的主要难题有：（1）水文变量及参数的空间变异性与尺度问题。我国流域尺度大、人类活动影响深。可根据流域不同地区的地形地貌特点，分区选取不同的计算网格步长，然后根据网格内土壤等参数的概率分布规律考虑其空间变异性对产汇流的影响。 （2）水循环的动力学机制的描述和计算量大之间的矛盾。水循环的许多过程如降雨时的入渗和地表径流过程变化快，描述这些过程常需要日以下的时间步长。如果所有过程所有时期均采用很短的时间步长，计算量将很大。因此，采用变时间步长，即针对不同过程及同一过程的不同时期采用不同的时间步长，将是缓解矛盾的对策之一。（3）
17、下包气带过厚滞后了降雨对浅层地下水的补给问题。我国许多地区特别是干旱半干旱地区的浅层地下水位往往很深，和地表之间存在很厚的包气带，滞后了降雨对浅层地下水的补给。可通过典型调查和观测，采取滞后曲线法、储留函数法等方法来解决。 （4）资料收集难与数据不足问题。分布式水循环模拟需要大量的基础数据。虽然我国的水文气象观测、地质调查与资料整编等基础工作开展较早、质量较高，但目前仍存在资料收集难与数据不足问题。四、结束语分布式流域水循环模拟和 GIS、DEM 和各类遥测技术相结合，解决水资源评价、洪水预报调度、水土流失、水污染以及水生态等各种生产实际问题，近年来已成为跨学科的国际研5究前沿。国际水文学会（
18、IAHS）2002 年将“观测资料缺乏流域的预测（PUBs） ”提议为下一个国际水文十年研究计划。欧美国家已开发出分布式流域水循环模拟与流域水资源管理、污染物运移或土壤侵蚀流失计算等耦合的应用系统，如美国 USGS 的 MMS 系统、欧洲的SHETRAN 模型等。因此，加快开发适应我国自然地理特征与气候特点的各类基于 GIS 的耦合式应用系统显得十分重要。此外，考虑到我国流域尺度大、人类活动影响深、环境复杂多变的实际情况，虽然传统的以率定参数为本的集总式水文模型无法客观地描述产汇流机制和预测人类活动带来的影响，但完全按数学物理方程模拟又受计算量的限制和尺度问题的困扰，因此基于物理概念和变时空步
19、长的分布式流域水循环模型将是未来的发展方向。参考文献1 Singh V.P. and D.A Woolhiser, Mathematical modeling of watershed hydrology，Journal of Hydrologic Engineering, ASCE, Aug.，2002，Vol.7, No.4, 270-2922 Jia Yangwen, Guangheng Ni, Yoshihisa Kawahara and Tadashi Suetsugi, Development of WEP Model and Its Application to an Urban
20、 Watershed, Hydrological Processes, May, 2001, Vol.15, 2175-21943 Jia Y., G. Ni, T. Kinouchi, J. Yoshitani, Y. Kawahara and T. Suetsugi, Assessing the Effects of Infiltration and Storm-water Detention Facilities in an Urban Watershed Using a Distributed Hydrological Model, Proc. of 29th Congress of
21、International Association of Hydraulic Research (IAHR), Beijing, China, Sept., 2001, Theme C, pp.379-3844 Jia Yangwen and Nobuyuki Tamai，Integrated Analysis of Water and Heat Balances in Tokyo Metropolis with a Distributed Mode，J. Japan Soc. Hydrol. & Water Resour. , March, 1998, Vol.11 No.2, pp.150-163作者简介 贾仰文，1965 年 9 月出生，1997 年 9 月获日本东京大学工学博士学位，现在中国水利水电科学研究院水资源研究所工作，院聘教授级高工，联系地址：北京市车公庄西路 20 号，邮政编码：100044，联系电话：10-68415522-3505，传真：10-68483367，电子信箱：为黄河 973 项目课题资助。
