01 MODFLOW-SURFACT概述

1.1 简介和版本历史

HydroGeoLogic, Inc. 很高兴呈现我们的 MODFLOW-SURFACT Version 3.0，这是一个完全集成的流动和输运代码，基于美国地质调查局（USGS）地下水建模软件 MODFLOW（McDonald 和 Harbaugh，1988）。MODFLOW-SURFACT 包含了对 MODFLOW 的额外模块，以提高其稳健性，并增加其物理模拟能力，包括复杂的饱和-不饱和地下流分析，以及在一个全面的、已经建立良好的仿真框架中进行的污染物命运和输运计算。MODFLOW-SURFACT 的地下水流动和溶质输运模块在代码文档的第 1 和第 2 卷中有详细介绍。本文档提供了代码安装、注册和操作的简要概述，适用于单独获取 MODFLOW-SURFACT 的用户。
MODFLOW-SURFACT的首次公开版本是1996年的1.2版。在MODFLOW-SURFACT 1.2版中，流动包（一般称为SURF包）扩展了原始MODFLOW代码的适用范围，以处理复杂的现场问题。

•    MODFLOW能够轻松适应含水层系统的脱水/再湿润化化条件。
•    MODFLOW具有执行轴对称分析的能力。
•    改进了MODFLOW模拟多层井的能力。
•    复杂问题的额外稳健数值和矩阵解决技术。
•    MODFLOW能够模拟指定的ponding补给和渗透面边界条件。
•    MODFLOW具有使用自适应时间步长方案高效执行瞬态模拟的能力。

版本2.1随后在1998年发布，提供了额外的物理流模拟功能和强大的数值方案，其中包括:

•    水流的严格饱和-非饱和建模。
•    空气流的严格饱和-非饱和建模。
•    Newton-Raphson线性化用于鲁棒性。

识别到与MODFLOW兼容的污染物传输代码很少，HydroGeoLogic已开发并将我们的新污染物传输分析（ACT）模块集成到MODFLOW-SURFACT中。版本1.2的传输建模功能包括：

•    水中高达五种物种的对流-扩散运输。
•    线性/非线性阻滞
•    土壤和/或水中的一级放射性衰变和生物化学降解，包括转化产物的生成。

Version 2.1 提供了额外的物理传输建模功能，包括：

•    在饱和-不饱和系统的所有相中进行质量分配，其中在活动相中进行对流+弥散输运，在非活动相中进行扩散输运。
•    从固定的（单一或多种）非极性液体（NAPL）相向另外两个相（空气和水）的质量分配，同时通过主动相中的对流+扩散和非活动相中的扩散进行污染物的传输，随后减少NAPL相。

这份文件介绍了MODFLOW-SURFAT的3.0版本。与2.2版本相比，版本3.0的增强功能包括：

•    可变各向异性选项 - 此选项允许在逐单元基础上包含水平各向异性。请注意，MODFLOW要求各模型层内的各向异性是均匀的。
•    垂直导电性输入选项 - 使用此选项，输入网格块的垂直导电性，并在内部计算层间饱和渗透系数。请注意，MODFLOW要求输入层间的饱和渗透系数。
•    MODFLOW-SURFACT版本还吸收了美国地质调查局（USGS）为MODFLOW-96和MODFLOW-2000的后续版本开发的其他模块。从MODFLOW-96中吸收的模块包括瞬时泄漏模块（TLK）、直接求解器模块（DE4）、河流路由模块（STR1）、水力流隔断模块（HFB）、互层储存模块（IBS）、恒定（时变）水头模块（CHD）、一般有限差分模块
（GFD）和流量与水头边界（时变）模块（FHB）。从MODFLOW 2000中吸收的模块包括水库模块（RES1）、河流流量路由模块（SFR1）、湖泊模块（LAK 3）、测量模块（GAGE）和蒸发蒸腾段模块（EST1）。
•    附加的模块
LAK3、SFR1、GAGE、TLK、STR1、HFB、ETS1、CHD、FHB 和 RES1已在各自的 USGS 文档中有所记录。当这些模块的输入在 MODFLOW-SURFACT 框架中被改变时（ETS1、FHB），这些模块的输入文档将在本文档中提供。
•    在此处实现的 LAK3 包中，当使用 LAYCON 的值时，应将其设置为无约束选项（LAYCON = 3、13、23、33、43），以使 LAK 3 包能够与此处的BCF 包正确配合使用。如果未满足此条件，输出列表文件中将反映出此警告。
•    在此处对 LAK3 包的再润湿过程进行了修改，使得湖泊水位在再润湿时略高于湖底。原始算法将湖水位设置为邻近地下水水头的平均值，可能导致再润湿时湖水位远高于湖底，从而在再润湿步骤中人为地增加了大量水
量。原始的 LAK3 包由 Merritt 和 Konikow（2000年）记录。

•   LAK2 软件包由GeoTrans 的Greg Council 于1996 年开发，已在作者许可下实施。 该包的输入说明和实现细节可从 GeoTrans 获得。
•  已包含观测节点包，用于列出在任何规定的观测位置随时间变化的水头或污染物浓度突破情况。
• 通过创建NAME 文件来允许MODFLOW96 样式输入，该文件打开模拟的所有相关包的输入文件。
• 对将MODPATH 与FWL4、FWL5 和ATO4 软件包一起使用进行了修改。 此增强功能将时间步信息提取到二进制 *.mph 文件（使用 ATO4 包时），以便与 MODPATH 的修改版本一起使用。 MODPATH 的原始版本仅利用 OC 包中的时间步长信息，并且不识别 MODFLOW-SURFACT 逐个单元输出的 FWL4 和 FWL5。
• MODPATH 的文档可在原始USGS 报告中找到（Pollock，1994）。 当前文档讨论了增强功能并为修改后的 MODPATH 提供了必要的输入指令。
• 支链衰减选项- 此选项允许使用适当的化学计量进行支链衰减计算。 最初的选项仅允许直链衰减，其中一个单位的母体衰减为一个单位的子产品。
• 增强了MODFLOW 蒸发蒸腾(EVT) 包的传输模拟。 以前的版本允许水由于 ET 而离开系统，并留下污染物质量。 增强功能包括提供根吸收因子，该因子确定用 ET 通量去除多少污染物质量，从而能够根据用户指定的不同程度去除由于 ET 造成的污染物。 此功能也已扩展到 MODFLOW 2000 中包含的蒸发蒸腾段包。
• BTN 文件中提供的用于输运模拟的IUNCCAD 选项考虑了仅在地下水流模拟情况下无约束网格块的饱和部分中发生的吸附。
• 当自适应时间步进包用于瞬态模拟时，现在可以按照用户定义的时间间隔执行质量平衡打印输出，而与头和浓度信息的打印无关。 （分别参见 AT04 包中的流量和传输质量平衡标记 IBUDFL 和 ICONBD）。
• 提供逐个电池质量平衡打印输出选项（参见AT04 包中的标志IUNITMB），以获取每个电池的绝对和归一化质量平衡。
• 包含可选的双域传输模块，可为聚合或断裂设置提供双孔隙度模拟功能。
• 选项NOMATRIX 已添加到传输方程的解中并通过BTN 文件输入。 关闭时（NOMATRIX = 0），系数矩阵将写入磁盘以供后续需要时使用。 当 NOMATRIX =1 时，该矩阵不会写入磁盘，而是在需要时进行计算。 根据计算环境的 CPU 和磁盘速度，其中一种可能比另一种更快。
• 密度相关传输模块可作为独家选项，用于沿海系统或流场受污染物密度影响的其他环境中的盐水入侵分析。
• 自适应时间步进包(ATO) 已修改为包含从数据集的某个中间点重新启动模拟的选项。

•    空气流量方程的可压缩存储项已得到修正。 该术语表达为：
 

更正为： 

仅在瞬态气流模拟中才能感受到这种变化的影响。 此外，这种效应仅对于没有显着强制空气运动的模拟才显着，其中该术语的不正确版本允许比正确形式更快地平衡到稳态。