SWAT水文模型

SWAT水文模型介绍

1概述

SWAT (Soil and Water Assessment Tool)模型是美国农业部(USDA)农业研究局(ARS)开发的基于流域尺度的一个长时段的分布式流域水文模型。它主要基于SWRRB模型，并吸取了CREAMS、GLEAMS、EPIC和ROTO的主要特征。SWAT 具有很强的物理基础，能够利用GIS和RS提供的空间数据信息模拟地表水和地下水的水量和水质，用来协助水资源管理，即预测和评估流域内水、泥沙和农业化学品管理所产生的影响。该模型主要用于长期预测，对单一洪水事件的演算能力不强，模型主要由8个部分组成：水文、气象、泥沙、土壤温度、作物生长、营养物、农业管理和杀虫剂。SWAT模型拥有参数自动率定模块，其采用的是Q.Y.Duan等在1992年提出的SCE-UA算法。模型采用模块化编程，由各水文计算模块实现各水文过程模拟功能，其源代码公开，方便用户对模型的改进和维护。

2模型原理

SWAT模型在进行模拟时，首先根据DEM把流域划分为一定数目的子流域，子流域划分的大小可以根据定义形成河流所需要的最小集水区面积来调整，还可以通过增减子流域出口数量进行进一步调整。然后在每一个子流域内再划分为水文响应单元HRU。HRU是同一个子流域内有着相同土地利用类型和土壤类型的区域。每一个水文响应单元内的水平衡是基于降水、地表径流、蒸散发、壤中流、渗透、地下水回流和河道运移损失來计算的。地表径流估算一般采用SCS径流曲线法。渗透模块采用存储演算方法，并结合裂隙流模型來预测通过每一个土壤层的流量，一旦水渗透到根区底层以下则成为地下水或产生回流。在土壤剖面中壤中流的计算与渗透同时进行。每一层土壤中的壤中流采用动力蓄水水库来模拟。河道中流量演算采用变动存储系数法或马斯金根演算法。模型中提供了三种估算潜在蒸散发戢的计算方法一HargreavesPriestley-Taylor 和Penman-Monteith。每一个子流域内侵蚀和泥沙量的估算采用改进的USLE方程，河道内泥沙演算采用改进的Bagnold泥沙运移方程。植物吸收的氮采用供需方法计

算，植物的氮口需求量是植物与生物量中氮浓度的函数。土壤中向植物供给氮，当需求超过供给时, 出现营养物压力。地表径流、壤中流和渗透过程运移的硝态氮量由水量和土壤层中的平均硝态氮浓度来估计。泥沙中运移的有机氮采用McElroy et al.开发的负荷方程，后经进一步改进。该负荷方程基于土壤表层的有机氮浓度、泥沙星和富集率来估计径流中的有机氮损失。植物吸引的磷采用与氮相似的供需方法。径流中带走的可溶解磷采用土壤表层中的不稳定磷、径流量和磷土分离系数来计算。泥沙运移的磷采用与有机氮运移相同的方程。河道中营养物的动态模拟采用QUAL2E 模型。

2.1产汇流模型

模型中采用的水量平衡表达式为：

t

SW t = SWo + 但切一Qsurf - Ea - W seep - Qgw )

i = l

式中：Sg为土壤最终含水最，mm； SWo为土壤前期含水量，mm： t为时间步

长，d；只切为第i天降水最，mm； Q^f为第i天的地表径流，mm： E&为第i 天的蒸发屋，mm； Wseep为第i天存在于土壤剖而底层的渗透星和侧流量，mm； Q卵为第i天地下水含帚:，mm。SWAT模型水文循环陆地阶段主要有水文、天气、沉积、土壤温度、作物产量、营养物质和农业管理等部分组成。模型径流量产生计算流程图如下图2-1所示。

图2-1 SWAT 模型产流计算流程图

2.1.1地表径流

当落到地表的降水星多余入渗量时产生地表径流。SWAT 模型采用SCS 径流 曲线法计算。

SCS 曲线方程自上世纪50年代逐渐得到广泛使用，属于经验模型，是对全 美小流域降水与径流关系20多年的研究成果。模型能反应不同土壤类型和土地 利用方式及前期土壤含水最对降雨径流的影响，它是基于流域的实际入渗最（F ） 与实际径流星（Q ）之比等于流域该场降雨前的绘大可能入渗ft （S ）与最大可 能径流量（Q m ）之比的假定基础上建立的

SCS 模型的降雨毎流基本关系表达式如下：

F _ S

Q = Q^

式中：假定潜在径流量（QQ 为降水量（P ）与由径流产生前植物截留、初渗和 填洼蓄水构成的流域初损（I a ）的差值。由此推到上式有：

1

地下水位 I

初损一受土地利用、耕作方式、灌溉条件、冠层截留、下渗、填洼等因素的

影响，它与土壤最大可能入渗量S呈一定的正比关系，美国农业部土壤保持局在

分析了大量长期的实验结果基础上，提出了二者最合适的比例系数为0.2,即：

I a = 0.2 S 由此可得SCS方程为：

=(P-0.2S)2

P > 0.2S

Q _ P+0.8S

P < 0.2S

流域当时最大可能滞留量S在空间上与土地利用方式、土壤类型和坡度等下垫而因素密切相关，模型引入的CN值可较好地确定S,公式如下：

25400

S = ^N—254

CN是一个无量纲参数，CN值是反映降雨前期流域特征的一个综合参数，它

是前期土壤湿度、坡度、土地利用方式和土壤类型状况等因素的综合。

2.1.2蒸散发

模型考虑的蒸散发是指所有地表水转化为水蒸气的过程，包括树冠截留的水分蒸发、蒸腾和升华及土壤水的蒸发。蒸散发是水分转移出流域的主要途径，在许多江河流域，蒸发量都大于径流量。准确地评价蒸散发量是估算水资源量的关键，也是研究气候和土地覆盖变化对河川径流影响的关键。

（1）潜存蒸散发

模型提供了Penman-Monteith> Priestley-Taylor 和Hargreaves 三种计算潜在蒸散发的方法，另外还可以使用实测资料或已经计算好的逐I I潜在蒸散发资料。一般采用Penman-Monteith方法来计算流域的潜在蒸散发。

（2）实际蒸散发

实际蒸散发以潜在蒸散发为计算基础。在计算流域实际蒸散发量的时候，模型首先计算植物冠层截留水分的蒸发，然后计算最大蒸腾量、最大升华量和最大土壤蒸发量，最后计算实际的升华量和土壤水分蒸发量。

（3）冠层截留蒸发量

模型在计算实际蒸发时假定尽可能蒸发冠层截留的水分，如果潜在蒸发量Eo 小于冠层截留的自由水S：E Iwr ,则：

E a = Ecai = Eo

E INT(f) = E INT(i) 一E can