我如何在 FloPy Modflow6 中输出所有时间步长的 MAW 水头值？

Question

我正在创建一口 MAW 井，想将其用作观察井以便稍后将其与现场数据进行比较，它应该在多层上进行筛选。但是，我只在输出文件的最后一个时间步长的井中获取水头值。关于如何在输出中获取所有时间步长的任何想法？

FloPy 手册说了一些关于它需要在输出控制中的内容，但我不知道该怎么做：

print_head (boolean) –   print_head (boolean) keyword to indicate that the list of multi-aquifer well heads will be printed to the listing file for every stress period in which “HEAD PRINT” is specified in Output Control. If there is no Output Control option and PRINT_HEAD is specified, then heads are printed for the last time step of each stress period.

在MODFLOW6手册中我看到可以进行连续输出： modflow6

我的 MAW 定义如下所示：

maw = flopy.mf6.ModflowGwfmaw(gwf, 
                      nmawwells=1, 
                      packagedata=[0, Rwell, minbot, wellhead,'MEAN',OBS1welllayers], 
                      connectiondata=OBS1connectiondata, 
                      perioddata=[(0,'STATUS','ACTIVE')], 
                      flowing_wells=False, 
                      save_flows=True, 
                      mover=True, 
                      flow_correction=True, 
                      budget_filerecord='OBS1wellbudget', 
                      print_flows=True, 
                      print_head=True,
                      head_filerecord='OBS1wellhead',
                      )

我的输出控件是这样的：

oc = flopy.mf6.ModflowGwfoc(gwf,
                        budget_filerecord=budget_file,
                        head_filerecord=head_file,
                        saverecord=[('HEAD', 'ALL'), ('BUDGET', 'ALL'), ],
                        )

希望一切都清楚，有人可以帮助我，谢谢！

Answer 1

您需要初始化 MAW 观测文件...它在 OC 包中没有完成。

您可以在此处的 MF6 文档中找到三个 MAW 示例的脚本： https://github.com/MODFLOW-USGS/modflow6-examples/tree/master/notebooks

看起来像这样：

obs_file = "{}.maw.obs".format(name)
csv_file = obs_file + ".csv"
obs_dict = {csv_file: [
                ("head", "head", (0,)),
                ("Q1", "maw", (0,), (0,)),
                ("Q2", "maw", (0,), (1,)),
                ("Q3", "maw", (0,), (2,)),
                ]}
maw.obs.initialize(filename=obs_file, digits=10, print_input=True, continuous=obs_dict)

我如何在 FloPy Modflow6 中输出所有时间步长的 MAW 水头值？

How do I in FloPy Modflow6 output MAW head values for all timesteps?

flopy