0

我正在创建一个 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'), ],
                        )

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

4

1 回答 1

0

您需要初始化 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)
于 2021-03-15T07:26:14.013 回答