VM飞机场例子演示教程 - 图文(12)

个别地增加一些浓度观测井，通过使用左手边菜单中的菜单栏选项，或者输入以 文件的格式输入井的值，这些文件的格式可以是 ASCII (.TXT) 文件, Excel (.XLS)文件, ESRI (.SHP)文件, 或者 MS Access (.MDB)文件。在本节中，是从一个 Excel 表格中输 入浓度观测井的值，如下表所示：

点击顶部菜单栏中的 Wells/Conc. Observation Wells 保存数据，然后点击 [Yes]。

点击左边菜单栏中的[Import Obs.] Open 输入窗口如下所示，在文件 Files of type 中输入 Excelfiles *.XLS,浏览数

地下水数值模拟技术与应用培训

据文件 Concentration wells.XLS.

点击[Open]打开以下窗口：

在此窗口中，在 Use column names in line 内输入 1，修改 Match to Column number 框架中的个数使得：Column# 1 和 Well Name 栏相匹配；Column# 2 与 X-coordinate 匹配；Column# 3 与 Y-coordinate 匹配；Column# 4 与 Screen ID 匹配; Column# 5 与 Screen Elevation 相匹配 field; Column# 6 与 Observation Time 相匹配; Column# 7 与

JP-4-Conc. 相匹配, 如上表所示： 点击[Next],出现以下窗口：

35

中国地质大学（北京）

点击[Finish], 完成了井的输入。 完成浓度观测井的输入后，点击底部菜单栏中的[F6-Zoom Out]按钮刷新屏幕。 选择底部菜单栏中的[F10 - Main Menu]，选择[Yes]保存数据，开始下一个模型 的建立。

地下水数值模拟技术与应用培训

模块III: Visual MODFLOW的运行

在模块 III 中，主要介绍了一些在 Visual MODFLOW 使用到的关于 MODFLOW, MODPATH 和 MT3D 的运行选项。

如果你跳过了 Module I & II, 你可以从打开在 Visual MODFLOW 安装目录的 Tutorial 文件夹下的 Airport3.vmf 工程开始。

继续运行选项，

点击主菜单选项中的[Run]，进入到运行选项的屏幕下，就会出现一个窗口让你去 选择进入 Steady State 还是 Transient 水流模型。

选择 Steady State 点击 [OK]

第十一部分：MODFLOW的运行选项

运行模型之前，对 MODFLOW, MODPATH and MT3D 的默认运行时间进行设置。 运行部分要求建立参数估计(PEST)模型。

MODFLOW 运行—时间选项的默认运行时间的设置包括以下几项： ?在各应力期（只对非稳定流模型有效）中的时间段 times—steps ?初始水头 initial heads 估计 ?结果 solver 选项和设置 ? 补给项 recharge 设置 ? 各层 layers 类型

? 疏干单元 Rewetting 的设置 ? Anisotropy 各项异性的设置 ? Output Control 设置

MODFLOW 计算方法选项 首先对 MODFLOW.的计算方法选项进行检查。 点击顶部菜单栏中 MODFLOW-2000/Solver 默认计算方法是 WHS Solver, 由 Waterloo Hydrogeologic Inc 所有开发的。单击 next 按钮到所需计算方法的缩写中相关联的计算方法参数就显示出来，选择不同的计算方 法选项。

Visual MODFLOW v.4.0 能过使用新的 SAMG 计算方法。在中等模型中该计 算方法要比传统的 PCG 方法要快 3 到 11 倍, 在包含一百万甚至更多单元的大模型中会 更快。 在本例中，使用的是 WHS Solver 来计算水流结果。 点击[OK] 选择 WHS Solver 计算方法。 疏干单元回灌选项

在水流模型中当水位低于单元底板时，单元疏干的现象非常普遍。疏干单元就变成 不活动单元，并且不能够被回灌，除非该回灌单元是活动的。在一个新的模型中，回 灌单元默认是不活动的，因此必须将单元变成活动的才能回灌，如下所示：

点击 MODFLOW/Rewetting，该 Dry Cell Wetting Options 窗口出现。. 选择 Activate cell wetting

37

中国地质大学（北京）

这个模型的在 Dry Cell WettingOptions 中的其他默认 值(如右图所示)的设置是非常好的。窗口底部增加了两个非标 准选项。实例中使用的是 WinPEST 的参数估计程序，它在迭 带过程中进行的自动调参会使有观测资料的一些单元疏干， 尤其会使参数估计变得混乱，因为 MODFLOW 对疏干单元赋 值的计算水头值是 1e-30。但是，Visual MODFLOW 的疏干单 元的水头赋值是该单元的底板高程。

Visual MODFLOW 在模型底层中提供了避免单元疏干的 选项。

点击 [OK]。

在水流模型中有其它一些可以使用的 MODFLOW 运行选 项，如果有时间可以对这里没有提到的选项进行研究。

第十二部分：MT3DMS 的运行选项

本节主要内容是关于水平对流方法的选择，计算方法的设置和溶质运移模型结果 的输出。

MT3DMS 计算方法的设置 点击顶部菜单栏中的 MT3DMS/Solution Method 一个 Solution Method 窗口就会出现，如下图所示：

在 Advection term 图框中, 使用 the Generalized

Conjugate Gradient (GCG) 计算方法中的 UpstreamFinite Difference 计算方法 Upstream FiniteDifference 方法为溶质运移模型在一个相对较短应力期中提供了相对稳 定可靠的计算方法。GCG 计算方法是使用一种内迭代的方法来计算有限差分方程，该 方法比外迭代方法更快。

点击 Advection term

，选择 Upstream Finite Difference。 选

择 Implicit GCG Solver 项

点击 Options]

Implicit GCG Solver Settings 窗口就会出现如下图所示：