MODFLOW-SURFACT/MODHMS 第二卷：传输模块

摘要

MODFLOW-SURFACT/MODHMS是由HydroGeoLogic，股份有限公司开发的一种完全集成的地下水流动和溶质输送代码，基于广受欢迎的美国地质调查局模块化三维（3-D）地下水流动建模代码MODFLOW。本卷（代码文件第二卷）描述了MODFLOW-SURFACT/MODHMS中包含的溶质传输模块，以模拟单物种和多组分污染物传输。详细介绍了传输计算模块的物理和数学概念，并讨论了这些概念在MODFLOW-SURFACT/MODHMS模块结构中的实现。还提供了运输模拟输入文件的准备说明。
为了验证数值技术并演示MODFLOW-SURFACT/MODHMS的现场应用，提出了三个级别的代码测试。对于第一级，考虑了五个具有可用分析解决方案的问题。这些问题的范围从半无限柱中的一维（1-D）输运的简单问题到由一对补给井和抽油井产生的非均匀流场中的相对复杂的输运问题。使用了各种案例来评估数值解的准确性，并证明了在MODFLOW-SURFACT/MODHMS中使用的总变差递减（TVD）方案的某些优势。对于二级测试，考虑了三个复杂的实际问题（例如，材料不均匀性、瞬态地下水位和泵送条件），这些问题排除了分析解决方案。这些示例问题旨在测试代码鈥檚 实用性和灵活性，并强调相关的模型功能（包括代码鈥檚 产生质量守恒和数值稳定且精确的解的能力），这对于真实的场模拟是必要的。另外四个测试问题需要进行二级测试，以验证MODFLOW-SURFACT/MODHMS第2版中包含的多相传输配方。最后，第三级代码测试侧重于实际的现场应用程序。提供了两个现场仿真实例。第一个例子涉及加拿大部队基地
（CFB）Borden垃圾填埋场下方浅层无侧限含水层系统中氯化物输送的横截面模拟。第二个例子涉及美国废物处理场放射性核素运输的三维模拟。本研究获得的结果与现场观测和先前的建模研究进行了很好的比较。

第一章：引言

1.1概述

MODFLOW-SURFACT/MODHMS是一个基于美国地质调查局模块化地下水流动模型MODFLOW的完全集成的地下水流动和溶质传输代码。自开发以
来，MODFLOW已被广泛使用，并通过加入额外的模块来增强其功能。最近的增强功能包括径流路由包、水平流动屏障包、夹层存储包和时变指定水头包。还对基本的有限差分程序包进行了修改，以避免细胞去饱和/再饱和带来的困难。修改中包括用于块间透射率平均的附加选项。所有这些在MODFLOW的早期版本中出现的增强功能都已被纳入MODFLOW-SURFACT/MODHMS
中。HydroGeoLogic，股份有限公司进一步扩展了MODFLOW的模拟能力，包括对无侧限系统进行严格的可变饱和处理、稳健而准确的解决方案、规定的积水补给和渗流面条件，以及对多层井的严格处理。我们的新流程包通常被称为“SURF包”，并在代码文档的第一卷中介绍。
认识到与MODFLOW一致的污染物运输代码很少，HydroGeoLogic，股份有限公司开发了兼容的运输模块，并将其与原始MODFLOW模块和SURF包相结合，以产生完全集成的MODFLOW-SURFACT/MODHMS代码。本报告构成了规范文件的第二卷，描述了新的污染物传输分析
（ACT）模块的配方、能力和使用。ACT模块提供的主要模拟功能包括以下内容：

•    稳态或瞬态流场中多种化学物质的对流分散传输
•    每个物种的线性或非线性延迟
•    水和土壤中污染物的一级衰变和生物化学降解
•    转型产品的生成。

版本2具有其他建模功能，包括

•    在非活性相中具有扩散质量运动的单一或多组分污染物的质量分配
•    单一或多组分污染物从固定的NAPL相中的质量分配，在活性相中具有平流和分散传输，在非活性相中具有扩散传输。
MODFLOW-SURFACT/MODHMS使用最先进的数值方案求解传输方程。代码中包括总变差递减（TVD）通量限制方案，旨在提供准确、物理正确和严格质量守恒的数值解。使用自适应隐式方案来最小化时间离散化误差。用户还可以控制解决方案的隐含程度，从完全隐含到显式。上游、中心和组合的上游-中心差分方案也被提供作为更快解决方案的选项（但通常不如TVD解决方案准确）。使用有效的Orthomin矩阵求解器来求解由有限差分近似得到的矩阵方程。传输方程的原始质量守恒形式用于所有方案，提供了严格质量守恒的数值解。所提到的各种方案旨在为新手建模师提供高质量的结果，并为经验丰富的建模师提供对解决方案选项的完全控制。
MODFLOW-SURFACT/MODHMS的输入/输出结构是MODFLOW格式的扩
展。MODFLOW-SURFACT/MODHMS在数值程序、数据结构、空间和时间离散化以及边界条件方面保持了流量和输运分析之间的完全兼容性。MODFLOW-SURFACT/MODHMS使用所有MODFLOW数据集进行运输模拟，并消除了流量和运输模拟的重复数据输入。关于传输参数和第一类型传输边界条件的信息通过附加的输入文件提供给代码。通过对MODFLOW边界数据文件的简单修改来提供第二类和第三类传输边界条件。

1.2 MODFLOW-SURFACT/MODHMS中传输模块的属性

MODFLOW-SURFACT/MODHMS为用户提供了以下地下水流动和污染物传输建模的关键优势：

•    利用计算机内存的动态分配进行地下水流动和污染物传输模拟的综合能力。此外，使用MODFLOW作为基础为需要污染物传输模拟能力的MODFLOW用户提供了快速学习。
•    严格解决传输方程的原始形式，同时保持与MODFLOW及其所有模拟选项的完全兼容性。其他基于MODFLOW的传输代码忽略了由压缩性效应、瞬态地下水位和泵送条件等因素引起的瞬态流量存储效应，这些因素会导致网格单元的排水和重新湿润。这种假设可能导致不切实际的运输模拟。
•    利用转化（子）产物处理污染物的运输，考虑一级生化降解或具有线性或非线性平衡吸附的放射性链衰变。

•    处理污染物的多相传输，包括不动的、消耗的NAPL阶段。
•    轴对称流和传输分析选项的可用性。这些选项允许用户获得含水层测试、单井示踪剂测试以及涉及井注入和泵送的污染物传输场景的高效分析。

•    准确、稳健、高效的大规模保守运输解决方案。这些方案无缝集成到MODFLOW-SURFACT/MODHMS中，适用于在微型计算机上模拟特定场地的现场特征和修复问题。
•    一致的运输数字解决方案选项。一种选择使用标准的上游加权方案，该方案是高效的，但仅一阶精确。第二种方案使用了计算要求更高、精度更高的TVD方案和van Leer磁通限制器。给定对数值精度和计算速度的某些要求，用户可以容易地选择合适的选项。

1.3 MODFLOW-SURFACT/MODHMS操作和输入选项

MODFLOW-SURFACT/MODHMS可以在独立模式下操作，也可以作为包括图形用户界面的完整软件包的一部分操作。如果用户仅获取MODFLOW-SURFACT/MODHMS或希望在用户界面系统之外执行模拟，则第一个选项是相关的。
本节重点介绍第一个选项（即在独立模式下运行MODFLOW-SURFACT/MODHMS）。在运行代码之前，用户需要为给定的模拟问题创建输入数据文件。对于每次模拟运行，MODFLOW-SURFACT/MODHMS执行地下水流动和溶质运移分析。鉴于此，用户在设置输入文件时应采取以下步骤：
1. 参考原始MODFLOW用户手册（McDonald和Harbaugh，1988），并按照手册中的说明设置相关流量包的输入文件。
2. 如果需要使用新的SURF软件包，请参阅MODFLOW-SURFACT/MODHMS代码文件第一卷的相关章节，并遵循所选流量软件包的输入说明。
3. 请参阅本文件（第二卷）第3章，并遵循运行ACT模块所需的相关软件包的输入说明。
注意，已经具有用于地下水流动模拟的现有MODFLOW输入文件的用户可以跳过步骤1和/或2。

1.4文件组织和使用指南

本文件分为八章。下文概述了其余章节及其目的。
第2章介绍了ACT的基本模块。对MODFLOW-SURFACT/MODHMS使用的传输模型公式、数值技术和质量平衡计算方案相关信息感兴趣的用户应阅读本章。
第3章介绍了基本ACT模块的输入说明。本章简要讨论了MODFLOW-SURFACT/MODHMS代码的功能操作。提供了设置输入文件的基本过程以及传输分析需要调用的单个包的特定输入指令。用户需要阅读第3章的相关章节，以便深入了解代码操作和传输输入准备。
第4章介绍了模型验证及其在地下水污染问题中的应用。用户应阅读本章，以获得与MODFLOW-SURFACT/MODHMS的分析解决方案测试、与其他可用数值模型（如MT3D和FTWORK）的比较以及现场应用相关的信息。
第5章提供了本文件前几章引用的参考文献列表。
第6章是附录部分，提供了与MODFLOW-SURFACT/MODHMS的双域模拟能力相关的公式、相关输入指令、示例应用程序和参考文献。有兴趣将断裂条件的复杂性纳入输送模拟的用户应阅读本章。
第7章是附录部分，提供了与MODFLOW-SURFACT/MODHMS的密度相关流量和传输模拟能力相关的公式、相关输入说明、示例应用和参考资料。对盐水入侵或污染物密度影响流场的相关问题感兴趣的用户应阅读本章。
第8章是附录部分，提供了与MODFLOW-SURFACT/MODHMS通用反应包相关的配方、相关输入说明、示例应用程序和参考文献。该包包括对工业污染物的生物降解感兴趣的几个特定反应，以及在几个可移动和不可移动的化学成分物种之间结合用户定义的反应的能力。它已被开发为RT3D（Clement，1977）反应模块的扩展，该模块无缝集成到MODFLOW-SURFACT/MODHMS的污染物传输程序中。