2.4 Numerical Modeling

该研究使用模拟代码MODFLOW-SURFACT和FEMWATER来模拟ETA排水田。接下来将对每个模拟代码的历史、公式和应用进行回顾。

2.4.1 MODFLOW-SURFACT

2.4.1.1 MODFLOW-SURFACT: Introduction and history

MODFLOW-SURFACT（Hydrogeologic，Inc.，1996年）是一个完全集成的流动和输运代码，基于美国地质调查局（USGS）地下水建模软件MODFLOW（McDonald和Harbaugh，1988年）。MODFLOW-SURFACT在原始MODFLOW代码的基础上提供了额外的模块，增强了MODFLOW的物理模拟能力，可以执行复杂的饱和-非饱和、地下流分析和污染物运输运算。MODFLOW-SURFACT版本1.2于1996年首次公开发布。版本1.2包含了增强原始MODFLOW代码处理复杂场问题的流动包。版本1.2中的附加包提供了MODFLOW执行以下任务的能力。

•轻松适应含水层系统的脱饱和/重新饱和条件。

•执行轴对称分析。

•改进多层井的模拟。

•使用额外的健壮数值和矩阵解法技术处理复杂问题。

•模拟规定的积水补给和渗漏面边界条件。

•使用自适应时间步进方案高效执行瞬态模拟。

1998年，版本2.1公开发布。版本2.1具有健壮的数值方案和额外的物理流动建模能力，以下简要总结。

•水流的严格饱和/非饱和建模。

•气流的严格饱和/非饱和建模。

•用于稳健性的牛顿-拉夫逊线性化。

最新版本，MODFLOW-SURFACT版本2.2，对前几个版本进行了以下流动建模增强。

•在每个模型层内的可变水平各向异性。

•当输入网格块的垂直导水性时，内部计算层间饱和漏水量，从而在计算和手动输入漏水值时节省时间。

•在面积平面上生成非矩形网格。

•在模型域中的预选节点处绘制水头或污染物浓度的突破曲线。

•在瞬态模拟中以用户定义的间隔执行质量平衡输出，而不考虑时间步长或涉及的时间步数。

MODFLOW-SURFACT的所有先前版本，包括版本2.2，都允许在原始MODFLOW和MODFLOW-SURFACT之间混搭模块。

MODFLOW-SURFACT模块与原始MODFLOW中的模块完全兼容。

原始MODFLOW的模块结构得以保持，增强功能的模块仍然从主程序中调用。

在附加包中实施新的公式和计算方案时，通过使用分析解和其他数值模型进行验证。

MODFLOW-SURFACT的输入格式结构简单直观，类似于原始MODFLOW的格式结构（Hydrogeologic，Inc.，1996）。

2.4.1.2 MODFLOW-SURFACT Version 2.2 

版本2.2包括几个新的流动包，这些包对原始MODFLOW版本进行了补充或增强，被称为SURF包。

以下是对SURF包功能的简要描述。

BCF4包是原始MODFLOW版本中基于块的流动包的替代方案。BCF4包允许在MODFLOW中模拟可变饱和流和轴对称模拟。版本2.2的BCF4包包含了使用Richards方程模拟空气和水的非饱和区域湿度移动等增强功能。

FWL4包增强了原始MODFLOW对建模裂缝井和抽水井的能力。FWL4包将井孔模拟为高导热性的裂缝管。高导热性裂缝管表示将与井相关的含水层单元与1-D有限直径圆柱形井单元连接起来，从而克服了原始井包的几个问题，并真正表示了多层井。<