将美国地质调查局模块有限差分地下水流模型中的不流单元转换为可变头单元的方法

A METHOD OF CONVERTING NO-FLOW CELLS TO VARIABLE-HEAD CELLS FOR THE U.S. GEOLOGICAL SURVEY MODULAR FINITE-DIFFERENCE GROUND-WATER FLOW MODEL By Michael G. McDonald, Arlen W. Harbaugh, Brennon R. Orr, and Daniel J. Ackerman

**将美国地质调查局模块有限差分地下水流模型中的不流单元转换为可变头单元的方法**

**作者：Michael G. McDonald, Arlen W. Harbaugh, Brennon R. Orr, and Daniel J. Ackerman**

**摘要：**
美国地质调查局模块地下水流模型（通常称为MODFLOW）采用有限差分法模拟多孔介质中的地下水流。被建模的区域被划分为一个单元网格，每个单元被定义为无流、可变头或定头。模型计算所有可变头单元的头值，而定头单元的头值由用户指定。如果单元包含不透水材料或是非饱和的，就将其指定为无流单元，因此在这些单元中水流不会被模拟。

在最初发布时，MODFLOW可以模拟可变头模型单元的脱饱和，从而将其转换为无流单元，但无法模拟单元的重新饱和。也就是说，无流单元无法转换为可变头单元。然而，在许多情况下，这样的转换是可取的。例如，可能希望模拟一些单元在抽水后变为非饱和状态，然后在停止抽水后水位的恢复。模型已添加了一个允许单元从无流转换为可变头的选项。在此选项中，单元基于相邻单元的头值被转换为可变头。该选项采用FORTRAN 77编写，并与现有模型完全兼容。本报告记录了新选项，包括概念描述、详细输入说明和代码清单。示例问题说明了该选项的实际应用。尽管修改后的流动方程的求解对于模型求解器可能较为困难，但示例问题表明可以解决各种复杂的问题。

**引言：**
美国地质调查局开发了一种三维有限差分地下水流模型（McDonald和Harbaugh，1988），通常称为MODFLOW。在该模型中，地下水流方程中的连续导数在称为节点的点上被有限差分逼近。每个节点周围都有一个单元，其中定义了水力特性，如水力导数和储存系数。结果是一组包含N个未知头值的N个方程，其中N是节点数。时间导数采用向后差分法逼近（McDonald和Harbaugh，1988，第2-16页）。

该程序在每个时间间隔结束时解决每个节点的未知头。

模型单元（i，j，k）的有限差分流动方程，其中i，j和k分别是行、列和层网格索引（McDonald和Harbaugh，1988，第2-26页），为：