02 MODFLOW-SURFACT卷一：流动模块

摘要

MODFLOW-SURFACT是由HydroGeoLogic, Inc.开发的完全集成的地下水流和溶质运移代码，基于广受欢迎的美国地质调查局模块化三维（3-D）地下水流模拟代码MODFLOW。本卷（代码文档的第I卷）描述了添加到MODFLOW中的新流动模块，以增强其地下水流模拟能力和计算鲁棒性。详细说明了新增增强功能的物理和数学概念，并讨论了这些概念如何实现到MODFLOW的模块化结构中。提供示例问题以验证代码，并使用户熟悉其应用。预期读者熟悉原始MODFLOW文档
（McDonald和Harbaugh，1988）。
MODFLOW使用块心有限差分方法来模拟地下水流。可以执行受限和无限制层的全或准三维模拟。通过其新的流动包，MODFLOW-SURFACT增强了执行无限制模拟的方案，以严格模拟含水层的除湿/再湿润，并克服了先前版本的MODFLOW遇到的数值困难。MODFLOW-SURFACT提供了使用轴对称几何离散化域的选项，以便有效地模拟抽水试验、水库抽水/恢复试验等。通常允许的外部应力包括恒定头、恒定通量、面积补给、蒸发蒸腾、排水和河流。此
外，MODFLOW-SURFACT提供了一个严格的井提取包、无限制的补给边界条件和渗透面边界条件。最后，MODFLOW-SURFACT包括自适应时间步长和输出控制程序的选项，以及额外的预处理共轭梯度（PCG）解包。MODFLOW-SURFACT的第二版包含额外的功能，包括严格的饱和-非饱和湿度运动模拟能力、空气流动模拟能力和用于改进鲁棒性的牛顿-拉弗逊线性化包。MODFLOW-SURFACT的新流动包是用FORTRAN 77编写的，并使用CompaqVisual Fortran Professional Edition 6.6®编译器编译。建议最低配置为带有64Mb RAM的任何基于Pentium的个人电脑上无需修改即可运行。

第一章：介绍

1.1 概述

美国地质调查局模块化流模型（MODFLOW）可能是全球最流行的地下水代码，被用作支持地下水调查的仿真工具。 MODFLOW已经通过HydroGeoLogic, Inc.大大增强，以允许用户以更强大和高效的方式处理复杂的现场问题。 增强的代码被命名为MODFLOW-SURFACT（即具有新的流动和污染物传输包）。 给新流包分配的SURF首字母缩写代表以下内容：

•    S -健壮而高效（矩阵解法、非线性迭代与自适应时间步长、输出控制以及轴对称）方案
•    U -Unconfined flow 的严格处理
•    R -" 非积水补给" 和" 渗漏面边界条件"
•    F -将裂缝井表示法提升到可以严格处理井条件。

HydroGeoLogic, Inc.开发的七个新模块包已经添加到原始的MODFLOW代码中（McDonald和Harbaugh，1988年）。这些包括：

(1)    BCF4- 这个包括轴对称分析选项的块中心流动包，采用可变饱和度表述和伪土
壤函数，对无限制流动进行严格处理。该包还包括 Version 2 的严格水流和气流模拟功能。
(2)    FWL4/FWL5-Fracture-Well Packages（FWL4和FWL5）使用一维（1-D）断裂管元素对井采条件进行严格处理，以模拟井。
(3)    RSF4-"可容纳非塘水充注和渗流边界条件的补给-渗流面包"
(4)    ATO4 -自适应时间步长和输出控制包是一个使用积极的时间步进方案进行瞬态解决方案的软件包，具有自动生成和控制时间步长以及更好的输出控制功能。
(5)    PCG4 -Preconditioned Conjugate Gradient Package包含了一种新的PCG矩阵求解器，作为MODFLOW中现有的Strongly Implicit Procedure（SIP）和Slice Successive Over-Relaxation（SSOR）求解器的稳健替代方案。
(6)    NRB1 -Newton-Raphson 线性化包结合回溯策略，提高了对非受限或非饱和流（空气或水）模拟的鲁棒性。
(7)    OBS1- 观测节点软件包，用于记录选定节点处的水文曲线。

这些附加的流动包在MODFLOW-SURFACT中得到了记录。以下章节提供了它们的简要描述、公式、验证和应用示例、输入说明以及样本数据文件。本文档的用户应当具有对原始MODFLOW代码及其文档的合理了解。对于不熟悉MODFLOW的用户，应参考以下参考文献：
McDonald, M.G., and A.W. Harbaugh (1988), Techniques of Water-Resources Investigations of the United States Geological Survey, Chapter A1: A Modular Three-Dimensional Finite-Difference Ground-Water Flow Model, Book 6, Modeling Techniques.

1.2 MODFLOW-SURFACT中新流包的属性

新的流程套餐和增强功能提供以下关键好处，如下所述。
网格块的完全脱水和重新饱和处理-原始的美国地质调查局MODFLOW的BCF包在瞬态模拟中存在网格块干燥/重新湿润的困难。
准确地勾勒和跟踪地下水位位置，考虑到非饱和区域的流动、延迟产流以及垂直流组成部分 - 原始的美国地质调查局 MODFLOW 的 BCF 包忽略了这些特征。
轴对称流动模拟选项 - 原始的美国地质调查局（USGS）MODFLOW的BCF包忽略了这些特性。
多层井的总流量自动且正确地分配到井节点——原始的美国地质调查局
（USGS）MODFLOW软件中的WEL模块要求对来自多个含水层或多个模型网格层的井抽水进行先验分配。
当井筒储存和超抽水井的需求大于井的供水时，原始的USGS MODFLOW中的WEL包会失败。
处理非集水或规定的集水重新充足条件 - 原始的美国地质调查局MODFLOW软件的RCH模块忽略了这些特征。
处理渗流边界条件 - 美国地质调查局原始的MODFLOW RCH模块忽略了这些特征。
自动生成和控制时间步长的自适应时间步长方案——原始的美国地质调查局
（USGS）MODFLOW的OC包忽略了这些功能。
B更好地组织和控制模拟输出 - 原始的USGS MODFLOW的OC程序忽略了这些功能。
强大高效的PCG矩阵解决方案选项。
对非饱和状态下的水分运移进行建模的能力。
使用双峰水力导度函数实现非平衡优先不饱和流的能力。
模拟非饱和空气流动的能力。
稳健高效的牛顿-拉弗逊线性化选项。
可变各向异性选项-此选项允许按单元格基础包括水平各向异性。请注意，MODFLOW要求各向异性在每个模型层内保持均匀。
垂直导水性输入选项-使用此选项，输入网格块的垂直导水性，并在内部计算层间饱和渗透系数。请注意，MODFLOW要求输入层间饱和渗透系数。
美国地质调查局（USGS）为MODFLOW的后续版本（包括MODFLOW-96和MODFLOW-2000）开发的其他模块也已纳入此版本的MODHMS中。从MODFLOW-96中纳入的模块包括瞬态渗漏包（TLK）、直接求解器包（DE4）、河流路由包（STR1）、水力流动障碍包（HFB）、层间储存包（IBS）、恒定（时变）水头包（CHD）、一般有限差分包（GFD）以及流量和水头边界（时变）包（FHB）。

  从MODFLOW 2000中纳入的模块包括水库包（RES1）、河流流量路由包（SFR1）、湖泊包（LAK3）、测量点包（GAGE）和蒸发蒸腾片段包
（ETS1）。请注意，为了使LAK3包与本文中的BCF包正确配合工
作，LAYCON应设置为非围限选项（LAYCON=3, 13, 23, 33, 43）。如果未满足此条件，则输出清单文件中将反映此警告。

•    用于将MODPATH与FWL4、FWL5和ATO4模块配合使用的修改。此改进将时间步信息提取到二进制*.mph文件中（当使用ATO4模块时），以便与MODPATH的修改版本一起使用。原始版本的MODPATH仅利用OC模块的时间步信息，并不识别MODHMS中FWL4和FWL5模块的单元格输出：请注意，由于MODPATH不包括表面水模块，因此只能将其用于MODHMS的地下部分。
•    曲线网格选项 - 此选项允许在平面上使用非矩形网格。因此，对于每个层中的每个节点，需要输入DELR和DELC，并且在第一层下方垂直堆叠多个层。请注意，原始的MODFLOW网格在平面上是矩形的。
•    LAK2模块，由GeoTrans的Greg Council在1996年开发，已获得作者许可并已实施。有关该模块的输入说明和实施细节可从GeoTrans获取。
•    已包含观测节点模块，可列出在任何规定的观测位置随时间的水头或污染物浓度的突破情况。
MODFLOW-SURFACT的所有模块均完全兼容，并且可以在单个模拟中混合使用来自原始MODFLOW和MODFLOW-SURFACT的模块。 MODFLOW的模块化结构得到了保留，并且增强功能的模块从主程序中调用。 新的公式和计算方案已经通过解析解和其他数值模型进行验证，并且
在附加包中实施。 MODFLOW-SURFACT的附加包与原始MODFLOW代码之间保持完全兼容性。 MODFLOW-SURFACT的输入准备简单明了，并遵循MODFLOW格式结构。 流动建模的输入数据基本上是原始MODFLOW所需的数据，增强功能需要的额外数据很少。 传输模拟利用所有MODFLOW数据集，并且消除了提供的输入信息的重复。 传输参数通过附加数据集提供给代码，并且边界条件通过MODFLOW边界数据集实施。

1.3 MODFLOW-SURFACT操作和输入选项

MODFLOW-SURFACT可以以独立模式运行，也可以作为一个包含图形用户界面的完整套件的一部分。如果用户仅获取MODFLOW-SURFACT本身，或希望在用户界面系统之外执行模拟，则第一种选项是相关的。
本节重点介绍第一种选项（即以独立模式运行MODFLOW-SURFACT）。在运行代码之前，用户需要为给定的模拟问题创建输入数据文件。为了构建地下水流模拟的输入数据文件，应采取两个关键步骤：

1.    参考原始的MODFLOW用户手册（McDonald和Harbaugh，1988年），并按照手册中的说明为相关的流动模块创建输入文件。
2.    如果希望使用新的SURF模块，则继续阅读本文档的相关章节，并根据所选模块的附加输入说明编辑相应的数据文件。
如果用户已经有现有的MODFLOW输入文件，则只需进行第2步。

1.4 文档组织和使用指南

该文档分为九个章节。以下是其余章节及其目的的概述。
第2章介绍了BCF4模块，其中包含新的可变饱和和轴对称流动公式。该章节包括公式描述、验证和应用示例，以及BCF4模块的输入说明。由于BCF4替代了以前的BCF模块，第2章为用户提供了在IUNIT(1)中设置块中心流输入文件的最一般的指导。
第3章介绍了FWL4和FWL5模块，这些模块包含了处理多层井和特殊井条件（例如，井筒储存和超泵水）的新方案。对于希望严格处理多层井并考虑特殊井条件的用户，应查阅本章。
第4章介绍了RSF4模块，该模块包含了处理充水边界条件（带有规定或零蓄水高程）
和渗流边界条件的新方案。对于处理具有复杂边界条件的非围限流问题的用户，应查阅本章。
第5章介绍了ATO4模块，该模块包含自适应时间步长方案和时间步长尺寸的自动控制，以及模拟输出控制。这些方案使得模拟能够更有效地进行。因此，对于对此类选项感兴趣的用户，应阅读本章。
第6章介绍了PCG4模块，该模块包含预处理共轭梯度（PCG）和Orthomin矩阵解法方案。对于有兴趣尝试PCG求解器的用户，应查阅本章。
第7章介绍了NRB1模块，用于牛顿-拉夫逊线性化。请注意，PCG4模块需要
与NRB1模块一起使用，因为它是MODFLOW中唯一能处理不对称矩阵的求解器。
第8章介绍了OBS1模块，用于记录指定观测节点的时间变化水头。第9章提供了前面章节引用的参考文献列表。
文档还包括五个附录（A至E），其中包含了描述第2至5章中所述示例问题的附加信息和输入数据文件。用户可以选择使用这些文件进行测试运行，并深入了解MODFLOW-SURFACT提供的新功能。另外四个附录（F至I）提供了MODFLOW-SURFACT/MODHMS扩展模块的输入说明。