SWAT水文模型
SWAT (Soil and Water Assessment Tool)是美国农业部(USDA)农业研究局(ARS)开发的流域尺度模型,用于模拟地表水和地下水的水质和水量,长期预测土地管理措施对具有多种土壤、土地利用和管理条件的大面积复杂流域的水文、泥沙和农业化学物质产量的影响。
水平衡在SWAT流域模拟中十分重要,流域的水文模拟分为两个部分:一是 水循环的陆地阶段,控制进入亚流域的水、沉积物、富营养物质和杀虫剂的数量;二是 水循环的演算阶段,定义通过流域水网到流域出口的水、沙等物质的运动。
1 SWAT水文模型原理
SWAT模型中主要含有水文过程子模型、土壤侵蚀子模型和污染负荷子模型,以下分别介绍三个子模型的原理。
1.1 产汇流模型
模型产流部分模拟的径流成分包括:地表径流、壤中流、地下径流。
1.1.1 水文循环的陆地阶段
1.1.2 水文循环的演算阶段
SWAT模型水文循环的演算阶段分为主河道和水库两个部分。
- 主河道的演算包括河道洪水演算、河道沉积演算以及河道营养物质和农药演算
- 水库演算主要包括水库水平衡和演算、水库泥沙演算、水库营养物质和农药演算。
SWAT 中的水流运动潜在路径如下图所示:
1.1.2.1 主河道演算
- 洪水演算:
- 河道沉积演算:沉积演算模型包含同时运行的两个部分(沉积和降解),沉积部分依靠沉降速度,降解部分依靠Bagnold的河流功率概念。从亚流域到流域出口的渠道和泛滥平原的沉积依靠沉积颗粒的沉降速度,沉降速度用Stokes Law粒径平方方程来计算。河道的沉降深度是沉降速度和河段行程时间的乘积,每一一个粒径的输送速率是沉降速度、行程时间和水流深度的线性函数,河流功率用来预测演算河段的降解。
- 河道营养物质和农药演算:目前,模型中还没有模拟在河道中营养物质和农药输移与降解的部分,并假设可溶性化学物质是保守物质,吸附到沉积物上的化学物质同沉积物一起沉降。
1.1.2.2 水库演算
水库水平衡包括入流、出流、表面的降水、蒸发、从库底渗漏、引水等,水库演算包括水库出流、沉积演算、水库营养物质和有毒物质演算。
- 水库水平衡和演算:水库水平衡包括入流、出流、地表降水、蒸发、库底渗漏、引水和回归流。目前有两种方法来评估出流:第一个方法简单的读入测量的出流,让模型模拟水平衡的其他部分;第二种方法用于小的不受控制的水库,当水量超过基本库容时,以特定的释放速率发生出流,超过紧急溢洪道的水量在一天内被释放,对于加大控制的水库,采用月目标水量方法。
- 水库沉积演算:对于水库和池塘的入流沉积量用MUSLE方程来计算。出流量用出流水量和沉积物浓度的乘积来计算,出流浓度根据入流量和浓度以及池塘储量的简单连续方程来估算。
- 水库营养物质和农药演算:使用Thomann和Mueller 的简单磷物质平衡模型,模型假定湖泊或水库内物质完全混合,可以用总磷来衡量营养状态。
1.1.3 地表径流
降水径流计算是计算土壤侵蚀的基础,当无降水过程资料而只有降水总量资料时,可以采用SCS曲线数法来进行降水径流模拟。SCS曲线方程自20世纪50年代逐渐得到广泛使用,属于经验模型,是对全美小流域降水与径流关系20多年的研究成果,模型用于计算不同土壤类型和土地利用条件下连续下垫面的径流量。
1.1.4 蒸散发量
模型考虑的蒸散发是指所有地表水转化为水蒸气的过程,包括树冠截留的水分蒸发、蒸腾和升华及土壤水的蒸发。蒸散发是水分转移出流域的主要途径,在许多江河流域及除南极洲以外的大陆,蒸发量都大于径流量。
准确地评价蒸散发量是估算水资源量的关键,也是研究气候和土地覆被变化对河川径流影响的关键。
1.1.4.1 潜在蒸散发
模型提供了Penman-Monteith、 Priestley -Taylor和Hargreaves三种计算潜在蒸散发能力的方法,另外还可以使用实测资料或已经计算好的逐日潜在蒸散发资料。一般可采用Penman-Monteith 方法来计算流域的潜在蒸发。
- Penman-Monteith方法
- Priestley -Taylor方法
- Hargreaves方法
1.1.4.2 实际蒸散发
在潜在蒸散发的基础上计算实际蒸散发。SWAT模型中,首先从植被冠层截留的蒸发开始计算,然后计算最大蒸腾量、最大升华量和最大土壤水分蒸发量,最后计算实际的升华量和土壤水分蒸发量。
1.1.4.3 层截留蒸发量
1.1.4.4 植物蒸腾
假设植被生长在一个理想的条件下,植物蒸腾可用以下表达式计算:
1.1.4.5 土壤水分蒸发
1.1.5 土壤水
下渗到土壤中的水以不同的方式运动着。土壤水可以被植物吸收或蒸腾而损耗,可以渗漏到土壤底层最终补给地下水,也可以在地表形成径流,即壤中流。
1.1.6 地下水
1.2 土壤侵蚀模型
1.3 污染负荷模型
SWAT模型可以模拟不同形态氮的迁移转化过程,地表径流流失、入渗淋失、化肥输入等物理过程,有机氮矿化、反硝化等化学过程以及作物吸收等生物过程,氮可以分为有机氮、作物氮和硝酸盐氮三种化学状态,氮的生物固定、有机氮向无机氮的转化以及溶解性氮随侧向壤中流的迁移等过程,有机氮又被划分为活泼有机氮和惰性有机氮两种状态,以及氨态氮挥发过程的模拟。
2 SWAT模型数据库
2.1 土壤属性数据
SWAT模型中用到的土壤数据主要包括两大类:物理属性数据和化学属性数据。
- 土壤的物理属性决定了土壤剖面中水和气的运动情况,并且对HRU中的水循环起着重要的作用。物理属性数据主要包括土层厚度、砂石、黏土、容积密度、有机碳、有效含水量、饱和水力传导率等。
- 土壤的化学属性主要用来给模型赋初始值。
- 其中物理属性是必需的,化学属性是可选的。
土壤输入文件.sol为土壤的各个层定义了模型模拟过程所需要的物理属性,.chm 文件为土壤的各个层定义了所需要的化学属性。
下表列出了SWAT模型土壤属性.sol、.chm 文件的各变量及含义,这些变量根据需要可定义10个层。
表1 SWAT模型土壤物理属性输入文件
表2 SWAT模型土壤化学属性输入文件
2.2 土地利用与植被覆盖
SWAT模型中有关土地利用和植被覆盖的数据通过文件crop.dat进行存储和计算,具体变量及定义见下表:
表3 SWAT 模型土地利用和植被覆盖属性输入文件
下载地址如下:
2.3 气象数据库
SWAT模型需要的气象数据主要包括流域的日降水量、最高/最低气温、太阳辐射、风速和相对湿度。这些数据可以是统计数据,也可以通过SWAT模型的天气模拟程序(Weather Generator)生成,或者是统计和模拟生产数据的结合。
理想情况下SWAT模型需要至少20年的天气资料,如果缺少必要的资料则可用天气生成程序进行补充。
在时间尺度上,模型的模拟时间步长可以为年、月、日。
2.4 DEM数据
DEM (Digital Elevation Model) 是描述区域地貌形态空间分布的GIS数据。DEM用一组有序数值阵列形式表示地面高程,是一定范围内规则格网点的平面坐标(X,Y)及其高程(Z)的数据集。在利用SWAT模型建模的过程中,将利用DEM数据进行流域 水系生成 以及 子流域划分 等工作。
2.5 水系图
2.6 径流数据(率定)
径流校准采用数字滤波技术将实测径流分割为直接径流(地表径流与壤中流)和基 流两部分, 然后分别对直接径流和基流进行校准。
参考
1.书籍-非点源污染模型——理论方法与应用
主要参考其中关于SWAT水文模型基本原理的讲解。