SWAT与Visual MODFLOW的耦合模型模拟应用与研究

2 模型构建

2.1 SWAT模型

SWAT（soil and water assessment tool）流域分布式水文模型（又称数字水文模型），属于物理模型，其运算效率高且能够进行长期的连续时间模拟 [6]。它利用ARCGIS提供的空间数据信息，并通过DEM提取大量的陆地表面形态信息，并根据一定的算法确定出地表水流路径、河流网络和流域的边界，然后在运用DEM划分的子流域单元上建立水文模型，进而模拟流域内的土壤-植被-大气系统中的水运动，最终实现地表水和地下水循环响应机理的综合演算 [7-8]。SWAT模型涉及的水文物理过程主要有降水、植被截留与蒸散发、下渗、地表和地下径流，如图1所示。其具体模拟的水平衡方程：

S_{w,t} = S_{w0} + \sum_{i=1}^{t} {R_{d,i} - Q_{s,i} - E_{a,i}-W_{s,i}-Q_{s,i}} 

式中， SW_t 表示最后 t 天末的土壤水含量；SW_0 是在第 1 天的初始土壤水含量； t 是时间；R_{d,i} \) 表示第 \( i \) 天的总降雨量；Q_{s,i} \) 表示第 \( i \) 天的总径流量；E_{a,i} \) 表示第 \( i \) 天的总蒸发量；W_{s,i} \) 表示第 \( i \) 天从包气带进入地下的总入渗量；Q_{g,i} \) 表示第 \( i \) 天的总回归量。

本次运用SWAT2009进行模拟，将MODFLOW中单元结构与水文响应单元互换，以实现对SWAT模型的改进 [9-10]。

2.2 Visual MODFLOW数学模型

Visual MODFLOW是综合的可视化地下水模拟软件 [11-12]，具有水质点的向前、向后示踪流线模拟、允许直接接受ARCGIS的输出数据及各种图形文件、简化数值模拟数据前、后处理等特点。近年来，Visual MODFLOW在国内应用也较为广泛。此次模拟运用Visual MODFLOW 4.6进行计算。根据研究区的水文地质条件得出地下水流系统的数学模型 [13]。

式中，Ω为渗流区域；μ为含水介质的给水度；n为边界面的法线方向；h为含水层的水位标高；h0为含水层的初始水位；W为含水层的垂直入渗补给量；Kn为边界面法向方向的渗透系数，m/d；Γ1为渗流区域的侧向和下边界；Kx，Ky，Kz分别是x，y，z向的渗透系数，m/d；q(x,y,z,t)Γ1边界的流量，流入为正，流出为负，隔水边界为零，m/d。

2.3 模型耦合

参考已有研究成果，针对能适应人类强烈活动的SWAT-Visual MODFLOW耦合模型，构建SWAT-MODFLOW的交互界面，并运用ARCGIS 9.3软件实现耦合计算。对模拟区进行概化，对单元格进行分类处理，即：将受人类活动影响大的平原区网格区域用Visual MODFLOW模型进行地下水模拟；对于受人类活动影响小的山区网格用SWAT模型的地下水模块进行模拟，这样既保证模拟精度又保证运算速度。