第 5 章：自适应时间步进和输出控制 (ATO4) 包

5.1 概述

在 MODFLOW 代码中，时域（用于瞬态流模拟）使用后向差分公式进行离散化，时间步值通过指定以下参数预先确定：时间步数 (NSTP)、应力周期的长度 （PERLEN），以及使用几何级数增加时间步长的乘法因子（TSMULT）。 下面列出了使用这种方案的限制。

• 如果矩阵求解器或非线性迭代方案（对于无约束情况）未能在给定时间步内收敛，则计算将中止。
• 即使获得了收敛的数值结果，预选的时间步长对于整体解决方案来说可能效率低下。
• 必须保持时间步长的几何级数。 结果，用户在分配特定时间值来检查流动系统方面缺乏灵活性。 因此，只能在时间级数上的时间值请求输出，这需要由用户提前确定。
为了克服上述困难并提高解决方案的效率，开发了自适应时间步进和输出控制包（ATO4）。 自适应时间步长方案根据给定计算的系统的预期非线性来选择时间步长大小。 如果预期的非线性不显着，则选择较大的时间步长以积极地推进模拟。 如果预期的非线性很严重，则选择较小的时间步长以确保该时间步长的收敛。 如果解在给定的时间步长内未能收敛，则进一步减小时间步长，并重复求解。
对于自动时间步进方案，仿真结果的时间值是未知的； 因此，ATO4 软件包还包括用于输出控制的新模块，该模块确定所需输出的时间值。 调用 ATO4 包后，不再需要使用 MODFLOW 之前的输出控制 (OC) 包。 请注意，原始 MODFLOW 的输出控制实现起来很繁琐且麻烦。 另一方面，新 (ATO4) 方案的输出控制则简单明了。 用户只需输入任何给定压力期间所需输出的时间值（以先前使用的时间单位），以及所请求输出的详细信息（即，头、回撤、质量预算和预算中的一项或全部）。 细胞间的流动项是否被输出）。 ATO4 包自动识别输出的选定时间值，并调整时间步长以在这些打印时间执行计算。 因此，例如，用户可以在 15 个月的压力期模拟的 10、11 和 15 个月后检查结果。 ATO4 包不仅在整个模拟过程中提供高效的时间步进，而且还确保在 10、11 和 15 个月（请求输出的时间值）执行计算。
最后，ATO 包提供了从非零起始时间值重新启动仿真的选项。 MODFLOW 的 OC 包不允许重新启动，并且开发新的输入文件来执行重新启动可能会变得极其繁琐，特别是在涉及多个压力周期时。

5.2 ATO4 包中使用的自适应时间步长的制定

ATO4 包中实现的自适应时间步长选择方案简单而有效。 用户提供初始时间步长 (DELT)、最大时间步长 (TMAX)、最小时间步长 (TMIN)、时间步长放大因子 (TSMULT) 和时间步长缩小因子 (TSDIV) 对于长度为 PERLEN 的每个应力周期。 初始时间步长 DELT 用于应力周期的第一个时间步长。 如果在最大迭代次数 (MXITER) 的 35% 以内实现收敛，则时间步长将增加乘数 TSMULT。 对于 MXITER 35% 到 65% 之间的收敛率，时间步长与之前的 DELT 相比没有变化，并且系统的收敛行为被认为对于该 DELT 值是最佳的。 然而，如果实现收敛所需的迭代次数超过 MXITER 的 65%，则时间步长将减小除法因子 TSDIV。 如果特定时间步长未实现收敛，则时间步长将减小 5.0 倍，并对较小的时间步长重新进行计算。 另外，如果达到打印时间(即，需要输出打印/保存的目标时间值)，则选择时间步长以包括打印时间值。 如果 DELT 的值减小到小于 TMIN，则模拟中止。 如果遇到严重困难，该技术提供了退出程序模拟的方法。 最后，DELT的值被限制为小于或等于TMAX。 即使非线性轻微或不显着，此约束对于时间离散化的准确性也很重要。 用户可以灵活地使用 TMVEC 阵列在指定时间进行打印，或在指定数量的时间步长 (NPSTP) 后进行打印，或两者兼而有之。

5.3 验证示例

对无承压含水层中的抽水进行瞬态模拟，以验证 ATO4 包并演示其实施。 BCF4 包文章中描述了所选问题，其中使用原始 MODFLOW 的时间步进和输出控制来执行模拟。 在计算效率和鲁棒性方面对这两次模拟运行进行了比较。
本例中的瞬态分析考虑了 300 英尺厚的无承压含水层，如图 5.1 所示。 建模域对应于尺寸为 75,000 英尺 x 75,000 英尺的正方形。含水层的顶部和底部分别位于 50 英尺和 -250 英尺的海拔处。 整个区域含水层的初始地下水位为零。 该域受到 3.0 x 10–9 ft/s 的均匀连续补给。 图中的左（西）边界是恒定（零）水头边界，其余边界是无流边界。 在含水层底部 100 英尺范围内筛选了 15 口井，其位置如图 5.1 所示。 每口井的抽水速度高达 3.85 ft 3/s。 四年后，井中的抽水就会停止。 （见图 5.2。）然后允许含水层恢复 12 年。 含水层参数为：

Horizontal hydraulic conductivity, K xx. Kyy =10–4 ft/s

Vertical hydraulic conductivity, K zz=10–5 ft/s

Specific Yield, S y=0.03
Specific Storage, S s=10–6 ft–1