分布式水文模型研究概况由于传统的流域水文模型本身所具有的局限性，同时随着水文循环中各个组成要素的深入研究，以及计算机、地理信息系统（GIS）和遥感技术的迅速发展，使构造具有一定物理基础的流域分布式水文模型成为可能。流域分布式水文模型充分考虑流域下垫面空间分布不均对水文循环的影响。在水平方向上将流域划分成许多单元网格和子流域（一般基于DEM，在垂直方向上将土壤分层，并依据流域产汇流的特性，使用一些物理的、水力学的微分方程 （如连续方程与动量方程）求解径流的时空变化。与传统的流域概念性集总水文模型相比具 有以下显著的优点：具有物理机理，能描述流域内水文循环的时空变化过程；其分布式结构，容易与GCM嵌套，研究自然和气候变化对水文循环的影响；由于建立在 DEM基础之上，所以能及时地模拟人类活动和下垫面因素变化对流域水文循环过程的影响。下面我简单介绍一下国内外的著名的分布式水文模型。主要从模型名称，模型结构， 输入输出变量，网格还是子流域为计算单元，适用性和范围等方面来描述。、分布式水文模型研究的发展现状-国际在国外，分布式水文模型的研究可以认为始于Freeze和Harlan 于1969年发表的一个具有物理基础数值模拟的水文响应模型的蓝图这篇文章。该文章提出了分布式水文物理模型的基本概念和框架。随后，Hewlett和Troenale 在1975年提出了森林流域的变源面积模拟模型 （简称VSAS）。1979年Bevenh和Kirbby提出了以变源产流为基础的TOPMODE模 型(TOPgraphy based hydrological MODEL )。1 TOPMODE模 型：1 (TOPgraphy based hydrological MODEL )该模型基于DEMt求地形指数，并利用地形指数来反映下垫面的空间变化对流域水文循 环过程的影响，模型的参数具有物理意义，能用于无资料流域的产汇流计算。TOPMOD模型的基础是变动产流面积的概念。流域降水满足冠层截留和填洼等初损以后，下渗进入土壤包气带，包气带分为土壤水带、 中间带和毛细水带。 只有包气带的含水量达到田间持水量后， 多余的水分中才有一部分以重力水的形式，通过大空隙直接进入饱和地下水层，所以入渗没有马上引起地下水位抬升至地表面。TOPMODE模型把全流域按DEM网格分块，每一个网格称为一个水文单元。大的流域又可被分成若干个子流域（或单元流域）。对每一个单元流域进行产汇流计算。产流计算包括不饱和层水分运动、饱和层水分运动及地表径流。地表径流和地下径流均视为在空间上相等，可通过等流时线方法进行汇流演算，求出单元流域出口处的流量过程。通过河道汇流得出流域总出口断面流量过程。河道演算多采用近似运动波的常波速洪水演算方法。TOPMODEL型可以根据地形指数和流域平均缺水量计算出各点的缺水量， 直观地反映源面积的大小和分布。地形指数的空间分布即反映了流域蓄水容量的分布。蓄水容量曲线与地形指数的功能相同，都是为了计算径流。在应用时，仅需DEMO和基本的水文资料（降水、蒸发、流量），甚至可用于无资料地区。填洼鹽理IJ时段内爭英地抵的F世单元网格产说计算:可网而输七结果 （出口址流虽）图1 TOPMODEI模型结构2 MIKE SHE 模型Beven 等，1980年；Abbott 等，1986年；Bathurst 等，1995 年；Refsgaard and Storm, 1995年；Chapters等等对SHE模型进行了改进，生成了很多不同版本的SHE模型匚2,4如：20世纪90年代初由丹麦水利学研究院 （DHI）在SHE模型上进一步发展研制出 MIKESHE模型。MIKE SHE模型功能上体现三维空间特性，包括了陆地区全部的水循环过程，同时对地下水资源和地下水环境问题分析、规划和管理是它的一大特色；采用多模块耦合的方式来模拟水循环中 几乎所有主要的水文过程，包括了大气循环、水流运动、溶质和泥沙输移等；通过连续计算四个不同且相互影响的储水层的含水量来模拟产汇流过程，这几个储水层代表了流域内不同的物理单元。这些储水层是：积雪储水层、地表储水层、土壤或植物根区储水层、地下水储 水层。因此它被广泛应用于流域管理、土地利用变化影响评价、地下水模拟、水质污染模拟、灌溉及农作物生长对水分和污染物质在非饱和带运移的影响等众多研究领域。模型中一个流域被沿水平方向划分成一系列的相互联系单元（grid ），各自具有不同的物理参数，而在垂直方向又被划分成若干层（zon e）,包括冠层、不饱和层和饱和层。它所反映的流域水文过 程主要包括降水、蒸散发、含植物冠层截留、地表汇流、河道汇流、非饱和壤中流和饱和地 下径流等过程，每个子过程分别进行计算建模。图2 MIKE SHE模型结构3 ANSWERS(AreaNonpoint Source Watershed Environment ResponseSimulation)和AGNPS(Agricultural Nonpoint Source)ANSWER模型5包括水文模型、泥沙分散-输送模型和几个描述坡面、亚表面、渠中的 水流路径的组件，采用概念模型模拟水文，用泥沙连续性方程模拟侵蚀 ，用方形网格划分研究 区域，可供水质规划者或其他用户模拟土地利用方式对水文和侵蚀响应的影响，对控制非点源污染进行规划。ANSWER模型中，首先项目管理器从用户那儿收集信息，然后利用GIS提取数据,产生一个输入文件，并将ANSWER模型的输出读入新的 GIS层,项目管理器的应用使输 入数据赋值时间减少了 710倍。AGNPS模型6采用方形网格划分单元，模型包含水文、侵 蚀和泥沙输送、氮磷和 COD的输移等内容，其中径流量用CREAM濮型中使用的公式计算，侵 蚀用RUSLE预测,化学物质的输移采用 CREAM濮型和一个饲育场评价模型中的方法，并对土壤结构的影响方面做了改进 ，化学物质的输移计算分为溶解相和泥沙吸附相的计算 ，溶解相 的计算与径流量有关，而泥沙吸附相的计算则与产沙量有关;WEPP ( Water Erosion Prediction Project )模型估计了陆地和水渠的径流和侵蚀、保护措施的影响，是一种计算山坡和集水区土壤侵蚀和泥沙输送的新技术，模型包括气候、表面和亚表面的水文、冬季冻融过程、灌溉、残余物的降解、沟渠和蓄水坑中泥沙的分散、输送和沉积等部分。最早的非点源污染模型是集总模型如CREAMS (Chemicals Ru nof and Erosion fromAgricultural Management Systems)模型，主要用于研究土地管理对水、泥沙、营养物和杀虫剂的影响，其中预测径流使用的是 SCS法(美国农业部土壤保持局曲线 )，产沙子模型 采用经验公式USLE(通用水土流失方程),预测污染物负荷采用的是概念模型。限于集总模型不考虑时空变异性，适用流域面积小，因而提出了用网格划分流域、可以模拟时空变异性的分布式模型。按模拟事件的时间连续性可分为单个事件模型和时间连续性模型8。单个事件模型出现较早，在模拟暴雨事件时，不考虑亚表面流、蒸散发、植物生长等水文因素，ANSWERS(Areal Nonpoint Source Watershed EnvironmentResponse Simulation) 和AGNPS(Agri2cultural Nonpoint Source)模型属于这一类；时间连续性模型则考虑了亚表面流、蒸散发、植物生长等水文因素，HSPF(Hydrological Simulation Program2Fortra n)、SWRRB(Simulator for Water Resources in Rural Basins)、SWAT(Soil Water AssessmentTool) 、WEPP (Water Erosion Prediction Project)属于这一类。图3 ANSWERS 口 AGNP水文模型结构专幣DEM1地剧用惦|土崛彊里留4K L. SPffi1rrL.K也莎布圈G FONI图4 ANSWERS口 AGNP水文模泥沙分散-输送模型结构4 VIC ( Variable In filtration Capacity)模型VIC( Variable In filtration Cap acity)模型9,10是一个大尺度陆面水文模型，它同新安江模型和通用产流模型一样，利用一个空间分布函数表示次网格内的土壤蓄水能力的变 化。VIC模型的参数分为四类：土壤参数、植被参数、水文参数和背景参数。背景参数指流域地理特性（纬度、经度和高程）和气候参数（近地表的多年平均气温和多年平均降水量）土壤参数、植被参数和背景参数具有比较明确的物理基础，可以直接给定，对于每一个网格的土壤和植被参数，分别基于全球10km的土壤数据库和全球 1km的陆地覆盖类型数据库来确定；模型中有7个待率定的水文参数，分别为B（入渗能力形状参数）、Ds（当基流非线性增长发生时，所占 Dm的比例）、Dm （这是底层土壤一天内产生基流的最大值）、Ws （当基流非线性增长发生时底层土壤含水量与最大土壤含水量的比值）、D1 （表层土壤厚度）、D2（第二层土壤厚度）、D3（第三层土壤厚度）。模型率定采用流域出口断面的实测流量过程线，因为它是流域水文过程的总体响应。VIC模型在中国的东部、东南部、中部、南部和中南部的湿润和半湿润地区的模拟效果较好，在西部和西北干旱区模拟结果不太理想。一方面因为干旱条件下多以超渗产流为主，而VIC模型是基于蓄满产流的；另一方面由于这些地区条件恶劣、人口稀少和经济欠发达，比湿润和半湿润地区的气象站网密度较稀疏，增加了模型输入的不确定性。基于 VIC （ Variable Infiltration Capacity）大尺度水文模型，吴志勇，陆桂华，张建云11等采用实测的日降水和日最高、最低气温数据，模拟了近35年（1971-2005）全国范围30km 30km分辨率的逐日土壤含水量。文章利用43个流域的实测流量资料，率定模型水文参数，建立和验证水文参数移用公式，确定无资料的水文参数，同时通过28个站的实测土壤含水量资料验证结果。variaklo Infiltration U 已(V ICIVImu 厂o u 曰 lero logic M ocle I1片/7 Vfefzruf冃(门“旧严冃口与tiy di-o_waslinfon_elu/l jftewimaiei/IVIortelsXVTC/VXrhoiT：Layor oLiLayor 2frif/ftration Ct/rvoAWuWu7 EW SiZuirf WOf、K ra ctional Ar& aWu图5 VIC模型结构图5 SWA水文模型由于SWAT模型是分布式水文模型中的一个杰出代表。所以下面我重点介绍一下SWAT水文模型理论12。5.1SWAT水文模型简介1985 年，美国农业部农业研究中心的 Alonso和De C