如何利用GMS构建地下水数值模拟模型？

___Y1

已于 2025-05-13 15:07:14 修改

阅读量819

点赞数 19

文章标签：学习算法

于 2025-05-12 20:44:19 首次发布

本文链接：https://blog.csdn.net/weixin_44259522/article/details/147901758

版权

0.建模准备

1.结构（模型网格剖分）

1.1 剖分

第一步：点图中三维网格的图标

然后点Grid→Create Grid

这个界面是用于创建一个有限差分网格的设置窗口。它允许你定义网格在三维空间中的范围、分辨率以及其他相关属性。这个界面提供了一个直观的方式来定义一个三维的有限差分网格，你需要指定每个维度上的范围、划分的单元格数量，以及一些控制单元格大小分布和网格方向的参数。这个网格通常会被用于数值模拟，例如地下水流模拟。在点击 "OK" 之前，请务必仔细检查所有参数是否符合你的建模需求。

1. X-Dimension (X 轴维度)

Origin (原点): X 轴方向的起始坐标，当前设置为 0.0 米。
Length (长度): X 轴方向的总长度，当前设置为 100.0 米。
Number cells (单元数量): X 轴方向上的单元格数量，当前设置为 10 个。
Cell size (单元尺寸): 根据给定的长度和单元数量自动计算出的每个单元格在 X 轴方向上的尺寸，当前为 10.0 米 (100.0 米 / 10 个单元)。
Bias (偏移): 用于控制单元格尺寸的非均匀性。当 Bias 不为 1.0 时，单元格的尺寸会逐渐增大或减小。当前设置为 1.0，表示单元格尺寸均匀。
Limit cell size (限制单元尺寸): 用于限制由于 Bias 造成的过大或过小的单元格尺寸。当前设置为 50.0 米。

2. Y-Dimension (Y 轴维度)

Origin (原点): Y 轴方向的起始坐标，当前设置为 0.0 米。
Length (长度): Y 轴方向的总长度，当前设置为 100.0 米。
Number cells (单元数量): Y 轴方向上的单元格数量，当前设置为 10 个。
Cell size (单元尺寸): 根据给定的长度和单元数量自动计算出的每个单元格在 Y 轴方向上的尺寸，当前为 10.0 米 (100.0 米 / 10 个单元)。
Bias (偏移): 用于控制单元格尺寸的非均匀性，当前设置为 1.0 (均匀)。
Limit cell size (限制单元尺寸): 用于限制单元格尺寸，当前设置为 50.0 米。

3. Z-Dimension (Z 轴维度)

Origin (原点): Z 轴方向的起始坐标，当前设置为 0.0 米。
Length (长度): Z 轴方向的总长度，当前设置为 10.0 米。
Number cells (单元数量): Z 轴方向上的单元格数量，当前设置为 1 个。
Cell size (单元尺寸): 根据给定的长度和单元数量自动计算出的每个单元格在 Z 轴方向上的尺寸，当前为 4.0 米 (注意这里可能是笔误，10.0 米 / 1 个单元应该是 10.0 米，显示为 4.0 米可能需要检查)。
Bias (偏移): 用于控制单元格尺寸的非均匀性，当前设置为 1.0 (均匀)。
Limit cell size (限制单元尺寸): 用于限制单元格尺寸，当前设置为 20.0 米。

4. Orientation / type (方向 / 类型)

MODFLOW: 下拉菜单，指示该网格是为 MODFLOW 模型创建的。MODFLOW 是一种广泛使用的地下水流模型。
Orientation (方向): 下拉菜单，允许你选择网格的方向。当前显示为 "—"，表示未选择或默认方向。
Type (类型): 下拉菜单，指定单元格的类型。当前设置为 "Cell centered" (单元中心)，意味着模型变量（如水头）存储在每个单元格的中心。

5. Rotation about Z-axis (绕 Z 轴旋转)

一个数值输入框，允许你设置网格绕 Z 轴旋转的角度，当前设置为 0.0 度，表示没有旋转。

6. 底部按钮

Help... (帮助...): 点击后可能会打开相关的帮助文档或教程。
OK (确定): 点击后将使用当前设置创建有限差分网格并关闭窗口。
Cancel (取消): 点击后将放弃当前的设置并关闭窗口。

1.2 新建MODFLOW

GMS MODFLOW/Basic Simulation

好的，我们来详细解读一下这个名为 "MODFLOW Global/Basic Package" 的界面。这个窗口通常用于设置 MODFLOW 地下水流模型的一些全局参数和基本属性。

1. 顶部标题栏

MODFLOW Global/Basic Package: 表明这是用于配置 MODFLOW 模型全局设置和基础包的窗口。

2. Heading (标题)

有两个文本输入框，允许你为 MODFLOW 模型输入两行标题信息。
- 第一行当前显示为 "GMS MODFLOW Simulation"。
- 第二行当前显示为 "12 May 2025"。这两行标题通常会出现在 MODFLOW 的输出文件中，用于标识模拟。

3. MODFLOW version (MODFLOW 版本)

一组单选按钮，用于选择使用的 MODFLOW 版本。当前选中了 "2000"。其他选项包括：
- 2005
- NWT (Newton-Raphson Formulation)
- LGR (Local Grid Refinement)
- USG (UnStructured Grid)
- USG Transport (非结构网格水流与溶质运移)
复选框:
- Save native text copy: 如果选中，将保存 MODFLOW 的原始文本输入文件。
- Save MODFLOW 6 copy: 如果选中，将保存 MODFLOW 6 格式的输入文件（如果所选版本支持）。
- Use custom Run dialog: 如果选中，将使用自定义的运行对话框。
- Double precision: 如果选中，将使用双精度浮点数进行计算，提高精度但可能增加计算时间。
- Parallel: 如果选中，将尝试使用并行计算来加速模拟。
Set As Default (设置为默认): 一个按钮，点击后可能会将当前的 MODFLOW 版本和相关设置保存为默认值。

4. Model type (模型类型)

一组单选按钮，用于选择模拟的类型：
- Steady state (稳态): 选中后，模型假设水流条件不随时间变化。
- Transient (瞬态): 选中后，模型将模拟水流条件随时间的变化。当前未选中。
Stress Periods... (应力期...): 如果选择了瞬态模型，点击此按钮将打开一个对话框，用于定义不同的应力期（模拟时段）及其持续时间和时间步长。

5. Run options (运行选项)

一组单选按钮，用于选择模型的运行方式。当前选中了 "Forward" (正演模拟)，即标准的地下水流模拟。其他选项包括：
- Parameter Estimation (参数估计)：用于根据观测数据反演模型参数。
- Stochastic (随机模拟)：用于考虑模型参数的不确定性。
- Stochastic Inverse (随机反演)：结合随机模拟和参数估计。
Sto. method: Parameter randomization (随机方法：参数随机化): 当选择 "Stochastic" 或 "Stochastic Inverse" 时，会显示相关的随机方法。当前显示为 "Parameter randomization" (参数随机化)。
Stochastic Options... (随机选项...): 点击后可能会打开一个对话框，用于设置随机模拟的具体参数。

6. Packages (包)

Packages... (包...): 点击此按钮将打开一个对话框，用于选择和配置 MODFLOW 的各个 "包" (Packages)。包是 MODFLOW 中用于模拟不同水文过程（如井、河流、降雨入渗等）的模块。
Units... (单位...): 点击此按钮将打开一个对话框，用于定义模型中使用的长度、时间和流量单位。
IBOUND... (IBOUND...): 点击此按钮将打开一个对话框，用于定义模型的边界条件类型。IBOUND 数组用于指定模型区域的活动单元、固定水头单元和非活动单元。
Starting Heads... (初始水头...): 点击此按钮将打开一个对话框，用于设置模拟开始时的水头分布。
Top Elevation... (顶部高程...): 点击此按钮将打开一个对话框，用于定义模型顶部的海拔高程。
Bottom Elevation... (底部高程...): 点击此按钮将打开一个对话框，用于定义模型底部的海拔高程。
Confining Beds... (弱透水层...): 点击此按钮将打开一个对话框，用于定义模型中的弱透水层（夹层）的属性。
LGR Options... (局部网格加密选项...): 如果 MODFLOW 版本选择了 LGR，点击此按钮将打开局部网格加密的设置对话框。
Layers... (层...): 点击此按钮将打开一个对话框，用于定义模型的水文地质分层及其属性。
Porosity... (孔隙度...): 点击此按钮将打开一个对话框，用于定义模型的孔隙度分布。

7. 其他复选框和输入框

Starting heads equal grid top elevations (初始水头等于网格顶部高程): 如果选中，模型将自动将每个单元格的初始水头设置为该单元格顶部的海拔高程。当前已选中。
Force convergence (CONVERGE) (强制收敛 (CONVERGE)): 如果选中，模型将尝试强制达到收敛，即使某些收敛标准可能没有完全满足。
No flow head (HNOFLO) (无流量水头 (HNOFLO)): 一个数值输入框，用于指定在非活动单元中分配的水头值。当前设置为 -999.0。

8. 底部按钮

Help... (帮助...): 点击后可能会打开与 MODFLOW 全局和基础包相关的帮助文档。
OK (确定): 点击后将保存当前设置并关闭窗口。
Cancel (取消): 点击后将放弃当前所做的更改并关闭窗口。

总结

这个 "MODFLOW Global/Basic Package" 界面是设置 MODFLOW 模型全局属性的关键部分。你可以在这里选择 MODFLOW 版本、模型类型（稳态或瞬态）、运行选项，并访问其他对话框来定义模型的几何形状、水文地质参数、边界条件、初始条件以及用于模拟不同水文过程的 "包"。仔细配置这些选项对于成功运行 MODFLOW 模型至关重要。

MODFLOW Packages/Processes

好的，我们来详细解读一下这个名为 "MODFLOW Packages / Processes" 的界面。这个窗口用于选择和激活 MODFLOW 模型中需要用到的各种 "包" (Packages) 和求解器 (Solver)。这些包是 MODFLOW 模拟不同水文过程和控制模型行为的关键模块。

1. 顶部标题栏

MODFLOW Packages / Processes: 表明这是用于选择 MODFLOW 模型中使用的包和求解器的窗口。

2. Flow package (水流包)

这是一个单选按钮组，用于选择一个且只能选择一个水流计算的核心包。当前选中了 "LPF - Layer Property Flow"。其他选项包括：
- BCF - Block-Centered Flow: 最早的水流包之一，基于块中心有限差分法。
- HUF - Hydrogeologic Unit Flow: 允许基于水文地质单元而不是层来定义模型属性，更灵活地处理复杂地质结构。
- UPW - Upstream Weighting: 也是一种水流包，通常用于模拟非饱和带或存在显著密度差异的水流。

3. Solver (求解器)

这是一个单选按钮组，用于选择用于求解水流方程组的数值方法。当前选中了 "PCG - Preconditioned Conjugate-Gradient"。其他选项包括：
- DE4 - Direct: 直接法求解器，适用于中小规模问题，但内存需求较高。
- GMG - Geometric Multigrid: 几何多重网格求解器，适用于大型复杂问题，收敛速度快。
- LMG - Link-AMG: 链接代数多重网格求解器，也是一种高效的大型问题求解器。
- NWT - Newton: 基于牛顿迭代法的求解器，通常与 UPW 包一起使用，适用于非线性问题。
- PCGN - PCG with Imp. Nonlin.: 预处理共轭梯度法的变体，适用于某些非线性问题。
- SIP1 - Strongly Impl. Proc.: 强隐式过程迭代求解器，较早期的求解器。
- SMS - Sparse Matrix \ LMG-USG: 稀疏矩阵求解器，可能与非结构网格 (USG) 相关。
- SOR1 - Slice Succ. Overrel.: 逐片超松弛迭代求解器，较早期的求解器。

4. Optional packages / processes (可选包 / 过程)

这是一个复选框列表，允许你选择模型中需要模拟的其他水文过程或添加额外的功能。你可以根据你的建模需求选择一个或多个包。以下是一些常见的包及其功能：
- BCT - Block Centered Transport: 块中心溶质运移包。
- BFH - Boundary Flow and Head: 指定随时间变化的边界流量或水头。
- CHD1 - Time-Variant Specified-Head: 指定随时间变化的固定水头边界。
- CLN - Connected Linear Network Process: 模拟管道网络中的水流。
- DPT - Dual Porosity Transport: 双孔隙介质溶质运移。
- DRN1 - Drain: 模拟排水系统。
- DRT1 - Drain Return: 模拟排水回流。
- ETS1 - Evapotranspiration Segments: 将蒸散划分为多个高程段进行模拟。
- EVT1 - Evapotranspiration: 模拟蒸散发。
- GAGE - Gage: 用于计算模型中的流量。
- GHB1 - General-Head Boundary: 更通用的固定水头边界条件。
- GNC - Ghost Node Correction: 用于处理某些数值问题。
- HFB6 - Horizontal Flow Barrier: 模拟水平方向的低渗透性障碍。
- LAK3 - Lake: 模拟湖泊与地下水的相互作用。
- MDT - Matrix Diffusion Transport: 基质扩散溶质运移。
- MNW1 - Multi-Node Well: 多节点井包的早期版本。
- MNW2 - Multi-Node Well: 多节点井包的中间版本。
- MNW1 - Multi-Node Well Information: 可能与 MNW 包相关的信息。
- PCB - Prescribed Concentration: 指定固定的溶质浓度边界。
- PEST-ASP: 与模型参数估计软件 PEST 相关。
- RCH1 - Recharge: 模拟降雨入渗补给。