MODFLOW 6 / Flopy 中的瞬态与稳态仿真

Question

如何在 Flopy 中将 MODFLOW6 模型定义为瞬态或稳态。根据我在这里here的发现，瞬态计算是通过使用stress_period_data来强制执行的。这是否意味着，一旦一个包使用stress_period_data 以及相应的 MFTransientList 对象，整个模型就会以瞬态术语进行评估？

Answer 1

不。您需要使用 STO（存储）包中的应力周期类型来定义流动状态（0=稳定，1=瞬态）。您可以在 TDIS 包中定义压力周期，然后在 STO 包中定义其类型。

软盘笔记本可能是阐明这些概念的最佳方式，并且它们最近已移至软盘文档中。

看看这个： https://flopy.readthedocs.io/en/latest/Notebooks/mf6_complex_model_example.html

在单元格 [7] 中，您将找到 STO 包的定义。

sto = flopy.mf6.ModflowGwfsto(
    gwf,
    pname="sto",
    save_flows=True,
    iconvert=1,
    ss=ss,
    sy=sy,
    steady_state={0: True},
    transient={1: True},
)

此外，包文档在这些情况下总是非常有帮助。 https://flopy.readthedocs.io/en/latest/_modules/flopy/mf6/modflow/mfgwfsto.html

请注意，在 MODFLOW6 中，您定义了此类第一个应力周期的瞬态特征，这将应用后续应力周期，直到您定义不同的东西（即另一个稳态 sp）。

Answer 2

我认为答案并不明确。为什么同时存在布尔条目“稳定”和布尔条目“瞬态”。这似乎是多余的。如果steady_state={0: True) 意味着transient={1: False}，反之亦然。函数中需要有两个条目吗？很混乱。