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径流曲线法。渗透模块采用存储演算方法,并结合裂隙流模型来预测通过每一个土壤层的流量,一旦水渗透到根区底层以下则成为地下水或产生回流。在土壤剖面中壤中流的计算与渗透同时进行。每一层土壤中的壤中流采用动力蓄水水库来模拟。河道中流量演算采用变动存储系数法或马斯金根演算法。模型中提供了三种估算潜在蒸散发量的计算方法—Hargreaves、Priestley-Taylor与Penman-Monteith。每一个子流域内侵蚀与泥沙量的估算采用改进的USLE方程,河道内泥沙演算采用

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

2、1产汇流模型

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

式中:为土壤最终含水量,mm;为土壤前期含水量,mm;t为时间步长,d;为第i天降水量,mm;为第i天的地表径流,mm;为第i天的蒸发量,mm;为第i天存在于土壤剖面底层的渗透量与侧流量,mm;为第i天地下水含量,mm。SWAT模型水文循环陆地阶段主要有水文、天气、沉积、土壤温度、作物产量、营养物质与农业管理等部分组成。模型径流量产生计算流程图如下图2-1所示。

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

2、1、1地表径流

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

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

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

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

初损受土地利用、耕作方式、灌溉条件、冠层截留、下渗、填洼等因素的影响,它与土壤最大可能入渗量S呈一定的正比关系,美国农业部土壤保持局在分析了大量长期的实验结果基础上,提出了二者最合适的比例系数为0、2,即:

由此可得SCS方程为:

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

CN就是一个无量纲参数,CN值就是反映降雨前期流域特征的一个综合参数,它就是前期土壤湿度、坡度、土地利用方式与土壤类型状况等因素的综合。

2、1、2蒸散发

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

(1)潜在蒸散发

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

(2)实际蒸散发

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

土壤蒸发量,最后计算实际的升华量与土壤水分蒸发量。

(3)冠层截留蒸发量

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

式中:为某日流域的实际蒸发量,mm;为某日冠层自由水蒸发量,mm;为某日的潜在蒸发量,mm;为某日植被冠层自由水初始含量,mm;为某日植被冠层自由水终止含量,mm。

如果潜在蒸发量大于冠层截留的自由水含量则:

当植被冠层截留的自由水被全部蒸发掉,继续蒸发所需的水分就会从植被与土壤中得到。

(4)植物蒸腾

假设植物生长在一个理想的条件下,植物蒸腾可用以下表达式计算:

当0≤LAI≤3、0时

当LAI＞3、0时

式中:为某日最大蒸腾量,mm;为植被冠层自由水蒸发调整后的潜在蒸发,mm;LAI为叶面积指数。因为没有考虑到植物下面图层的含水量问题,由此公式计算处的蒸腾量可能比实际蒸腾量要大一些。

(5)土壤水分蒸发