如何在Flopy中实现跨多层的井（Multi-layer well screen）？

我需要在模型域中实现覆盖多个层的井，但Flopy的WEL包设置时，每条应力期数据只能指定单个层的井，示例代码如下：

# 创建井包
pumping_rate = -500.0
wel_sp = [[0, nrow / 2 - 1, ncol / 2 - 1, pumping_rate]]
stress_period_data = {0: wel_sp}
wel = flopy.modflow.ModflowWel(mf, stress_period_data=stress_period_data)

我试过在下层再设一口井并将抽水量减半：

# 创建井包
pumping_rate = -500.0
wel_sp1 = [[0, nrow / 2 - 1, ncol / 2 - 1, pumping_rate/2]]
wel_sp2 = [[1, nrow / 2 - 1, ncol / 2 - 1, pumping_rate/2]]
stress_period_data = {0: wel_sp1, 1: wel_sp2}
wel = flopy.modflow.ModflowWel(mf, stress_period_data=stress_period_data)

但这种均分抽水量的方式在各层传导系数不同时并不合理，想知道Flopy有没有内置方法实现跨多层的井？

Flopy本身没有内置的跨层井定义功能，因为底层MODFLOW的WEL包就是基于单层层位设计的。不过可以通过按层导水系数比例分配抽水量的方式，来更合理地模拟跨层井，替代简单均分的做法：

• 先计算目标层的导水系数（渗透系数×层厚，即Transmissivity）；
• 按各层导水系数占总导水系数的比例，分配总抽水量；
• 将多个层的井记录放在同一个应力期数据中，而非不同应力期。

示例代码如下：

# 总抽水量
total_pumping = -500.0

# 计算目标层的导水系数
# 假设已通过模型的LPF和DIS包获取渗透系数与层厚
t_layer0 = mf.lpf.hk[0] * mf.dis.thickness[0]
t_layer1 = mf.lpf.hk[1] * mf.dis.thickness[1]
total_transmissivity = t_layer0 + t_layer1

# 按导水系数比例分配抽水量
rate_layer0 = total_pumping * (t_layer0 / total_transmissivity)
rate_layer1 = total_pumping * (t_layer1 / total_transmissivity)

# 同一应力期内添加两个层的井记录
wel_stress_data = [
    [0, nrow//2 - 1, ncol//2 - 1, rate_layer0],
    [1, nrow//2 - 1, ncol//2 - 1, rate_layer1]
]
stress_period_data = {0: wel_stress_data}
wel = flopy.modflow.ModflowWel(mf, stress_period_data=stress_period_data)

如果井的过滤器在不同层的长度有差异，还可以结合过滤器长度进一步调整分配比例（比如用「导水系数×过滤器长度」的占比），让抽水量分配更贴合实际工况。

内容的提问来源于stack exchange，提问作者Aeonian