第 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）。

 图 3.1 压裂井的物理情况 (a) 和概念化 (b)。

其中 Δx 和 Δy 是包含孔的单元格的水平网格块尺寸。 然后，FWL5 程序包会同时求解井中的水头和多孔基质块节点，这些节点彼此不同； 因此，多孔基质网格块不需要在井周围进行细化，以获得更好的井流量、水头和提取条件的准确性。
FWL4 和 FWL5 软件包要求用户提供的信息包括：井位、特定井筛选的顶层和最底层的标识、井底高程、井半径、抽取 (-ve) 或注入 (+ve) ）来自井的速率（L3 T–1 ），以及井筒特定储存值（L–1）。 此外，FWL5 软件包还需要输入表面摩擦系数和在任何应力周期激活任何井时的井内初始水头。 FWL4 和 FWL5 模块
组织几何信息，使用方程（25）计算Kl，并将离散裂缝井方程（24）应用于各个井节点。 然后，在每个应力周期，将抽取或注入应用于井中最底部的筛选层。 该公式将应力适当地分布在井节点之间。 最后，当需要时，井节点处的边界条件（即，从规定的通量到规定的水头）的切换自动执行，如下所述。
每口井按规定的速率继续抽取，直至井内水位达到井底。 井的底部标高实际上位于被筛选的最低节点层中的任何位置。 当最低井节点处的测压水头对应于井的底部高程时，水头值保持固定，并且允许井产量低于需求并达到零。 如果在后期，含水层系统恢复并且井产量增加到大于需求，或者如果需求下降到产量以下，则井恢复以规定的抽水速率抽水。 输出文件中报告了总井产量和每层对井产量的贡献。 在以下部分中，提供说明性示例来验证和演示 FWL4 和 FWL5 软件包的应用。

3.3 验证示例

通过将数值结果与 Theis (1935) 以及 Papadopulos 和 Cooper (1967) 获得的解析解进行比较，验证了压裂井包 (FWL4) 选项，该解析解针对具有恒定流量井的均质含水层中的水位下降。 与 Theis 解不同，Papadopulos 和 Cooper 的解考虑了井筒储存。 因此，进行了两种模拟，一种有井眼存储，另一种没有井眼存储，以测试 FWL4 中使用的数值方案。
Sudicky 等人描述了一个验证问题。 （1995）。 在此问题中，半径为 0.1 m 的全穿透井位于尺寸为 1,000 m x 1,000 m x 5 m 的均质各向同性含水层的中心。 该井的抽水速度为 10 –3 m3/s。 横向边界规定了 20 m 的恒定水头，并且假定含水层的顶部和底部边界是不渗透的。 含水层的水力传导率为 10–4 m/s，比容为 10–4 m–1。 使用 FWL4 软件包执行瞬态仿真。

 建模方法和结果

该域被离散为 5 层，每层 1 英尺厚，具有 35 行和 35 列可变网格间距（从井附近的 111.1 英尺到 0.1 英尺）。 含水层的顶部和底部高程分别设置为 0 和 –5 m。 MODFLOW 要求的 VCONT 值（垂直导水率除以任意两层之间的层间距离）计算为 10–4 s–1。 含水层的主要储存系数（特定储存时间块厚度）指定为10 –4。 含水层侧面规定恒定水头为 20 m，并假设初始水头均匀为 20 m。 0.1 m 半径井位于区域中心（第 18 行、第 18 列），允许以 10-3 m3/s 的速率泵送 1000 秒。 由于井已完全贯通，井顶和井底分别位于1层和5层，井底标高指定为–5 m。 为了将井压降的时间变化与 Theis 和 Papadopulos 的解析解进行比较，在具体井筒储水量分别为零和 0.1 m –1 的情况下进行了两次瞬态模拟。 对于这两种模拟，PCG4 软件包用于求解矩阵方程。 图 3.2 绘制了井面的压降与时间曲线以及 Theis 和 Papadopulos 解析解。 如图所示，使用FWL4 Package获得的MODFLOW-SURFACT/MODHMS结果与各自的解析解非常吻合，从而验证了FWL4 Package中采用的压裂井公式。