第 3 章：压裂井（FWL4 和 FWL5）套件

3.1 概述

HydroGeoLogic, Inc. 在美国地质调查局模块化地下水流模型 MODFLOW 中添加了新的井包，可将井筒模拟为高电导率压裂管。 这些包是 MODFLOW 中现有井包（WEL1 包）的替代方案，被称为 Fracture-Well（FWL4 和 FWL5）包。 Fracture-Well 封装旨在克服与原始 WEL1 封装相关的多个问题。
首先，FWL4 和 FWL5 封装模拟通过多层孔的流动。 这是因为裂缝管表示允许将与井相关的含水层单元（节点）与一维有限直径圆柱形井单元连接。 此外，还严格考虑了该井规定的总提取率。 来自与压裂井相关的每个单独节点的体积通量由代码自动计算，以求和到从井中提取的总流
量。 另一方面，原始的 WEL1 包仅在代表井的含水层节点处实现源/汇项。 因此，它要求用户确定每个节点对穿透两个或多个节点的物理井的贡献。 通常估计这些节点通量贡献，以便根据贡献井节点的透射率来分配规定的提取率。 这种估计本身可能是一个严重的错误，从而歪曲了物理系统。 此
外，对于无约束系统，根据 WEL1 包规定节点通量可能会导致计算困难。 然而，FWL4 和 FWL5 软件包会自动将净井提取量（由用户输入）分配给井节点。 对于每个时间步长，计算出的节点通量可能会有很大变化，就像各个层去饱和和再饱和时的情况一样。
新FWL4和FWL5套件的另一个特点是，对于超抽无侧限系统，当井中水位达到井底时，总井抽水量会自动调整。 在此阶段，底部井节点的水头固定在井底条件下，并计算总抽取量以对应于该井在这些条件下可以供应的可行数量。 如果发生抽水需求小于井在井底压力条件下可供应的情况（即，在特定时期减少抽水或增加补给），井节点处的水头将反弹 为了反映这些新的情况，该井将提供全部的开采需求。 当水位低于井底时，以前的 WEL1 套件不会物理调整抽水速率。 相反，它会继续计算到非物理领域，井节点处的水头值不断下降到井底标高以下。 在其他情况下，WEL1 包会在迭代之间循环干燥和润湿泵单元，从而导致系统在计算上变得不稳定。
FWL4 和 FWL5 软件包还实现了井眼存储，这在 MODFLOW 的早期版本中被忽略。 FWL4 和 FWL5 软件包之间的差异在于，FWL4 软件包假设叠加，将井流方程添加到多孔基质的方程中，而 FWL5 软件包使用对数井函数来描述从基质块到井节点的流动。 因此，如果用户有兴趣使用 FWL4 软件包准确预测局部水位下降，则仍然需要在每口井附近进行精细的面积离散化。
最后，FWL4 和 FWL5 封装仅适用于 BCF4 封装的严格可变饱和模拟选项（即输入 LAYCON = 40 或 43）。 本章的以下部分进一步阐述了 FWL4 和 FWL5 包的功能，并包括对公式和演示示例问题的简要讨论。 压裂井公式的数值实现也根据文献中可用的分析解决方案进行了验证，然后是各个包的输入说明。

3.2 裂缝井方程的建立

通过将代表有限半径井眼筛选部分的 1-D 线元素叠加到代表多孔介质的 3-D 单元上，为 FWL4 包制定压裂井，如图 3.1 所示。 苏迪基等人。 (1995) 详细介绍了压裂井单元的推导和使用，并提供了示范性示例。 出于一般性的考虑，我们将他们的技术应用于 BCF4 包中实现的可变饱和流量公式。 该公式进一步扩展到 FWL5 包，其中使用井函数在多孔基质块和井节点之间发送流量，而不是使用叠加假设。 这将在本节稍后讨论。 FWL4 和 FWL5 套件中通过管状裂缝的流量均由一维可变饱和流量方程控制，如下所示:

式中，l为沿井筛的距离（L），Kl'为井的导水率（LT-1），Krw'为井内的相对渗透
率，Sw'为井单元的无量纲饱和厚度（ 饱和厚度与单元厚度之比），Ss' 为井筒比储存量（L–1），h' 为井内水头（L）。
伪土壤函数用于定义压裂井中的不饱和流动行为，如 BCF4 包文档（第 2 章）中所述，以便井单元内的地下水位条件可以通过垂直平衡条件准确表示。 井筒的饱和水力传导率，根据层流理论，应用 Hagen-Poiseuille 方程估算，可表示为 (Sudicky et al., 1995)：

其中ρ是水的密度(ML–3)，g是重力加速度(LT-2)，µ是水的粘度(ML–1 T–1)，rw是井的半径(L) 。 通常，对于大直径井，井水力传导率值比多孔基质的块间水力传导率大几个数量级。 因此，仿真结果对使用等式(25)可能产生的误差不敏感。 如前所述，FWL4 包使用井节点与多孔基体节点的叠加。 另一方面，FWL5 软件包将它们保持为单独的节点，并使用 Theim 方程井函数提供矩阵块与其关联井节点之间的流动。 将 Bennett 等人 (1982) 提出的概念应用到任意层 k 的裂缝井公式中，基质块与其井节点之间的通量由下式给出：

其中 Qmwk 是 k 层井通量矩阵； krk 为相对渗透率，Kk 和 Thickk 为基质块的水平电导和厚度，hmk 为 k 层基质块中的水头，hwk 为 k 层井节点的水头，Sk 为表皮系数 re 是矩阵块的有效半径，定义为（Halford 和 Hanson，2002；Neville 和 Tonkin，2004）。