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模型的地下水模块进行模拟,这样既保证模拟精度又保证运算速度。