SWAT水文模型

1 概述

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

2 模型原理

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

发量的计算方法— Hargreaves 、 Priestley-Taylor 和 Penman-Monteith。每一个子流域内侵蚀和泥沙量的估算采用改进的USLE方程，

河道内泥沙演算采用改进的Bagnold泥沙运移方程。植物吸收的氮采用供需方法计算，植物

的氮日需求量是植物与生物量中氮浓度的函数。土壤中向植物供给氮，

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

和磷土分离系数来计算。泥沙运移的磷采用与有机氮运移相同的方程。河道中营养物的动态模拟采用QUAL2庄莫型。

产汇流模型

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

式中：SWi.为土壤最终含水量，mm 为土壤前期含水量，mm t为时间步长, d; Rda,为第i天降水量，mm h$urr为第i天的地表径流，mm 为第i天的蒸发量，mm叫卿为第i天存在于土壤剖面底层的渗透量和侧流量， mm。回为第i 天地下水含量，mm SWAT莫型水文循环陆地阶段主要有水文、天气、沉积、土

壤温度、作物产量、营养物质和农业管理等部分组成。模型径流量产生计算流程图如下图2-1所示

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

地表径流

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

SCS 曲线方程自上世纪50年代逐渐得到广泛使用，属于经验模型，是对全 美小流域降

水与径流关系20多年的研究成果。模型能反应不同土壤类型和土地 利用方式及前期土壤含水量对降雨径流的影响， 它是基于流域的实际入渗量（F ）

与实际径流量（Q ）之比等于流域该场降雨前的最大可能入渗量（S ）与最大可能 径流量

（i 」）之比的假定基础上建立的

SCS 模型的降雨-径流基本关系表达式如下:

I F s

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

（p - J

Q 二 s - i a

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

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

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

I a = 0. 2S

由此可得SCS 方程为：

r (P - 0, 2S)-

Q ~ F + 0-8S

9=0

流域当时最大可能滞留量£在空间上与土地利用方式、土壤类型和坡度等下

垫面因素密切相关，模型引入的卅值可较好地确定•，公式如下:

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

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

蒸散发

模型考虑的蒸散发是指所有地表水转化为水蒸气的过程， 包括树冠截留的水

分蒸发、蒸腾和升华及土壤水的蒸发。蒸散发是水分转移出流域的主要途径， 在 许多江河s 二 25400

CN 254

流域，蒸发量都大于径流量。准确地评价蒸散发量是估算水资源量的关键，也是研究气候和土地覆盖变化对河川径流影响的关键。

(1)潜在蒸散发

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

(2)实际蒸散发

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

(3)冠层截留蒸发量

模型在计算实际蒸发时假定尽可能蒸发冠层截留的水分，如果潜在蒸发量卜』

小于冠层截留的自由水量两，贝

ElNT(f)二ElNT(i) - Ecan

式中：E』为某日流域的实际蒸发量，mm比餉为某日冠层自由水蒸发量，mm Eo为某日的潜在蒸发量，mm E1KI⑴为某日植被冠层自由水初始含量，mm ⑴为某日植被冠层自由水终止含量，mm

如果潜在蒸发量I珂大于冠层截留的自由水含量二DT则：

Ecan 二ElNT(i)

二0

当植被冠层截留的自由水被全部蒸发掉，继续蒸发所需的水分就会从植被和