weixin_44259522的博客

1.1 通用1.1 通用。

MODFLOW代码采用了基于区块的有限差分方法来解决地下水流方程。流动域被离散为行、列和层，使得每个节点表示一个多孔材料的矩形块，被称为单元。结果有限差分网格中的节点表示为不流、可变头或恒定头单元，与单元相关的任何水力特性都相对于相应的节点指定。对于非围限流模拟，原始MODFLOW代码的BCF1包（McDonald和Harbaugh，1988）将脱水的可变头单元转换为不流单元。这导致模型模拟永久性排除了流动域的部分部分。

ThePlot AxesAxestab of thePlot Axestab of thePlot AxesAxesPlot Axesand the。

MODFLOW代码通过使用Recharge（RCH1）包来容纳地下水补给边界条件。然而，RCH1包在表示不受限系统时有一定的局限性。如果不受限的含水层饱和到地表，含水层吸收供应的补给水的能力会减小，剩余的水会以地表径流的形式流失。RCH1包无法处理这种情况。相反，它继续为含水层提供补给水（就像受限系统一样），使水头持续上升，远远超过地表。MODFLOW-SURFACT的新RSF4包旨在消除这种物理上不现实的情况。RSF4包允许将补给水引入地下水系统，前提是水位低于用户规定的池（积水）高程。

在MODFLOW代码中，时间域（用于瞬态流动模拟）使用反向差分公式进行离散化，通过指定以下参数预先确定时间步长值：时间步数（NSTP），应力周期长度（PERLEN）以及一个乘法因子（TSMULT），该因子使用几何级数递增时间步长大小。使用这种方案的局限性如下。• 如果矩阵求解器或非线性迭代方案（用于非均质情况）在给定的时间步长内无法收敛，则计算将中止。• 预选的时间步长大小可能对整体解决方案效率低效，即使获得了收敛的数值结果。• 必须保持时间步长大小的几何级数。

MODFLOW代码以有限差分形式生成描述地下水流系统的方程组。通常通过使用Slice Successive Over-Relaxation (SSOR)或Strongly Implicit Procedure (SIP) Package来求解这组代数方程。对于复杂的场地模拟，这些求解程序可能缺乏效率和稳健性。Hill（1994）提供了一个Preconditioned Conjugate Gradient (PCG2) Package，为MODFLOW提供了另一种求解器。

对于水或空气作为活动相的严格非饱和带流模拟通常是高度非线性的，需要小的时间步长和多次迭代来解决实际情况。在一些情况下，对于非饱和情况，使用MODFLOW的Picard迭代方案，或者对于非饱和流（空气或水）情况，使用MODFLOW-SURFACT/MODHMS（第2章）的修改Picard方案，或者在不需要进行严格非饱和流模拟的不受限制的含水层中，一步稳态模拟可能无法收敛。提供二次收敛的牛顿-拉夫逊方案可以极大缓解收敛困难。

在MODFLOWSURFACT/MODHMS中包含了一个观测节点包，以简化瞬态模拟期间在域内任意节点处的突破曲线的后处理。流动模拟输出显示了在模拟的请求节点处的压力-时间关系，而传输输出包括每种污染物模拟的浓度-时间关系。会创建一个单独的输出文件，其中包含突破输出，可以轻松进行后处理以提供突破曲线。

MODFLOW-SURFACT/MODHMS is a fully integrated groundwater flow and solute transport code developed by HydroGeoLogic, Inc., based on the widely popular U.S. Geological Survey modular three-dimensional (3-D) groundwater flow modeling code, MODFLOW. This volum

ACT模块被整合到MODFLOW-SURFACT/MODHMS中，用于在瞬态或稳态流场中对最多五种溶质物种进行对流-扩散传输方程的数值求解。在传输分析中使用的数值方案严格保守质量。允许每种物种的线性或非线性（Freundlich）平衡吸附。允许在流体相和土壤中进行一阶生化降解。还实现了直链降解（或放射性衰变）并生成转化产物。提供了选项，允许只在活动相中进行传输；在空气和水相中进行传输（其中一个对流方程无效，其中扩散是传输机制）；以及在空气和水相中进行传输，其中污染物从不动的耗尽的NAPL相中溶解和挥发。