第3章 基于 SWAT的海伦市地表水径流模拟

        SWAT (Soil and Water assessment Tool) 模型是反映流域下垫面空间分布特征和地表径流模拟的半分式水文模型。

通过将 ArcSWAT、AvSWAT 模块嵌入 ArcGIS 软件，实现了地表径流模型的可视化效果。相比传统水文模型，其空间数据管理、数据处理和展示功能大大增强 [75]。模型建立所需的基础资料包括土壤特征、流域的地形、气象资料、土地利用、水文资料等 [76]。

SWAT 模型与 ArcGIS 紧密结合，在建模前需对研究的基础数据进行整理与分析，在此基础上构建模型需要的用户数据库，例如土壤、土地利用类型、天气发生器等。

本文以通肯河流域上游为研究区，联合利用 ArcGIS 和 SWAT 软件，建立海伦市基础数据库和 SWAT 数据库。

3.1模拟理论与运算过程

SWAT 模型由许多变量（蒸散发量、河道汇流量、地表径流量、土壤水分运移、地下水补给）组成，实现对流域内的水文循环的计算与模拟 [78]。

受流域下垫面性质和降雨过程的影响，将 SWAT 模型分为子流域计算模块和河道产汇流模块 [79]。

子流域计算模块对流域河道内水的输入量进行模拟，其受气候、地表产流和地表覆盖等因素的影响；河道汇流模块计算水、泥沙、污染物等物质从个子流域向流域总出口的输送情况[80]。

水量均衡是模型计算和统计流域水分随时间变化的根本依据，SWAT 模型内置的水量平衡均衡模块，运算式如下[81-84]：

其中：SW0为初始土壤水量，单位mm；SWt为t时间节点土壤含水量，单位mm；t为时间，单位d；Rday为第i天的降水量，单位mm；Qsurf为地表汇水量，单位mm；E为蒸发量，单位mm；Wseep为土壤层的旁侧流量，单位mm；Qgw为地下水的回归流量，单位mm。

3.1.1地表径流

在SWAT模型中，通过将所有HRU的计算结果汇聚到子流域，得到河道出水口的径流量。地表径流的多少受土壤水分的影响较大。如果土壤处于干燥状态，此时下渗率最大，地表径流不会形成；随着土壤逐渐湿润且水分的增多，下渗率逐渐减小，逐渐达到填洼临界点，即地表径流产生[85]。对于SWAT模型，模型自带地表径流估算方法：径流曲线法(SCS)和格林-安姆普特(GREEN-AMPT)。本次研究获取到日尺度数据，因此使用SCS径流曲线法来进行地表径流计算。计算公式如下[86-90]：

 Q_{surf} = \frac{(R_{day}-I_a)^2}{(R_{day}-I_a+S)}

S = 25.4 \times (1000 / CN - 10)

其中: Qsurf—地表汇水量,单位mm;            Rday—当天降水量,单位mm;

Ia—降雨初损,单位mm; S—潜在最大滞留量,mm;            CN—曲线系数,流域降雨前期特征的参数,无量纲。 

3.1.2蒸散发量

        SWAT模型中的蒸散发包括:植物蒸腾与升华、树冠截留的水分蒸发、土壤水的蒸发。流域水量输出的主要途径是蒸散发。蒸散发计算:Priestley-Taylor、Hargreaves、Penman-Monteith方法。本次研究选用Penman-Monteith方法[91-93],详细计算公式如下: 

\lambda E_t = \frac{\Delta(H_{net} - G) + \gamma k_1 (0.622 \frac{\lambda \rho_{air}}{P}) e^0_s - e_s/\gamma_a}{\Delta +\lambda(1+\gamma_c/\gamma_a)}

其中:SW0—热通量密度,MJ/m2·d;

E—潜发深率,单位mm/d;

Δ—曲线斜率,kpa/℃;

Hnet—净辐射,单位MJ/m2·d;

G—热通量密度,MJ/m2·d;

ρair—空气密度,单位kg/m3;

Cp—热量,MJ/kg·℃;

P—大气压,单位KPa;

k1—8.64×104;

es—大气水汽压,单位KPa;

es0—饱和水蒸汽压;

γc—植物叶孔阻抗;

γa—空气动力学阻抗。 

3.1.3 土壤水分运移

        入渗到土壤中的水存在多种运动形式，有一部分通过蒸散发损失掉，其中一部分通过入渗补给进入地下水含水层，余下为壤中流。SWAT利用动力贮水方法计算 [94-96]：