裂隙承压含水层中延伸井椭圆流的精确解析解：在页岩气开发区地下水水头解释中的应用

《Processes》：Exact Analytical Solutions for Elliptical Flow Toward Extended Wells in Fractured Confined Aquifers: Application to Groundwater-Head Interpretation in Shale-Gas Development Areas

　　

本研究开发了瞬态椭圆地下水流向各向异性裂隙承压含水层中延伸井的精确解析解，并讨论了由此产生的水力响应如何支持页岩气开发区的地下水水头解释。环境联系是在含水层保护尺度上建立的：该模型不是页岩气储层生产模型，也不直接求解污染物传输。相反，它提供了一个水力解释框架，用于估算各向异性、等效裂隙长度、井筒存储效应以及可能泄漏点、旧井、裂隙或监测井周围水头传播的优先方向。基于Mathieu函数理论（Mathieu-function theory）和分离变量法（separation-of-variables method），研究人员在拉普拉斯空间（Laplace space）中推导了恒定速率和恒定水头解，并通过Stehfest算法（Stehfest algorithm）反演到时域。解析结果与COMSOL5.2有限元模拟进行了验证，并通过降深（drawdown）和流量典型曲线（type curves）分析了各向异性系数和井筒存储的影响。包含一个合成但具有现场风格的水头示例，以展示如何将监测记录转换为降深，拟合到椭圆流解，并用于划定初步水力响应区。结果表明，各向异性主要控制早中期响应，而井筒存储可能掩盖早期水头变化并延迟裂隙连通性的识别。因此，该解最好被视为一种快速水力筛选和监测设计工具，可以在页岩气地下水保护研究中先于但不会取代特定地点的污染物传输建模。相关的技术问题包括可能与地表泄漏、返排水处理、旧井、断层和裂隙连通含水层相关的水头扰动和优先地下水流路径。由于研究人员无法获得经过验证的当地现场监测记录，因此应用示例明确描述为合成的现场风格演示；它用于展示工作流程及其局限性，而不是声称对污染物浓度进行特定地点预测。

页岩气开发中水平钻井和水力压裂技术提升了低渗透页岩层的产气效率，但同时也引发了地下水环境保护的公众与科学关注。潜在的水相关问题包括压裂液或返排液泄漏、井筒完整性失效、废水储存以及优先裂隙或断层作为通往含水层的路径等。尽管这些问题并不意味着每个场地都会发生污染，但它们表明地下水监测和水力解释是页岩气区环境管理的必要组成部分。现有研究多集中于经典径向流或线性流模型，但裂隙承压含水层中有限长度延伸井（如与主裂隙连通的井）会产生椭圆形水力响应，而非圆形对称。传统径向流解无法表征这种非径向响应，而线性流模型假设无限长裂隙，数值模型则需大量场地数据与计算资源。因此，研究人员提出一种基于Mathieu函数理论（Mathieu-function theory）和分离变量法（separation-of-variables method）的解析模型，专门针对有限延伸井的椭圆流问题，旨在为页岩气开发区的地下水水头解释和快速水力筛选提供理论工具。该研究发表于《Processes》。

研究人员采用了以下主要关键技术方法：（1）基于椭圆坐标系，建立裂隙承压含水层中延伸井的瞬态流控制方程，通过坐标变换将各向异性问题转化为等效各向同性形式；（2）利用Mathieu函数理论在拉普拉斯域（Laplace domain）分离变量，推导恒定速率（CRT）和恒定水头（CHT）条件下的解析解，并通过Stehfest算法（Stehfest algorithm）反演到时域；（3）采用COMSOL5.2有限元模拟进行验证，评估各向异性系数和井筒存储效应的影响；（4）通过合成但具有现场风格的水头数据示例，展示如何将监测记录转换为降深并拟合解析解，用于划定水力响应区并指导监测井布设。研究未涉及具体场地样本队列，所有验证基于数值模拟和合成数据。

研究结果如下：

3.1 Validation with Numerical Simulation（数值模拟验证）
通过COMSOL5.2有限元模拟对比，验证了椭圆流解析解的准确性。在无井筒存储（C_D=0）和有井筒存储（C_D=20）条件下，无量纲井筒降深随时间的变化曲线吻合良好，定量误差指标（如均方根误差约10^-2，最大相对偏差低于5%）表明解析解能正确描述椭圆流行为。验证还确认了无存储解在C_D=0时自动简化为标准形式。降深分布图显示等势线呈椭圆形，不同于圆形对称，说明水力扰动具有方向性。

3.2 Effect of Anisotropy Coefficient ξ_w on Drawdown Curves（各向异性系数ξ_w对降深曲线的影响）
各向异性系数ξ_w与水力传导率比K_x/K_y相关，ξ_w越小表示各向异性越强。恒定速率试验中，不考虑井筒存储时，ξ_w较小（即强各向异性）导致更大的降深，尤其在早中期（t_D < 100）。考虑井筒存储（C_D=5）时，早期曲线几乎重合，但中期差异显现，ξ_w仍控制降深幅度。因此，各向异性主要影响早中期响应，井筒存储则会掩盖早期信号。

3.3 Wellbore-Storage Control on Drawdown Behavior（井筒存储对降深行为的控制）
井筒存储系数C_D显著影响降深曲线：在t_D < 1000时，降深受井筒存储控制，早期呈现近似单位斜率的直线；随着时间增大，存储效应消失，各曲线趋同。这意味着监测中若井筒存储较大，早期弱响应不应误判为水力隔离，需进行存储校正或延长测试时间。

3.4 Effect of Anisotropy Coefficient ξ_w on Flow Rate Curves（各向异性系数ξ_w对流量曲线的影响）
恒定水头试验中，无量纲流量随时间递减，ξ_w越小流量越低，该效应在t_D < 1000时显著，后期减弱。恒定水头条件可理想化表示固定水头源或汇（如泄漏源），该解可用于评估裂隙连通源与含水层之间的水交换速率，适用于返排水储存或废弃井的泄漏潜力筛选。

3.5 Application to Water-Head Data in a Shale-Gas Groundwater Monitoring Scenario（页岩气地下水监测场景中水头数据的应用）
通过一个合成案例展示了工作流程：假设一口恒定速率抽水井（EW-1）和三个观测井（MW-A、MW-B、MW-C），其中MW-A沿裂隙方向。生成48小时水头记录并转换为降深，拟合椭圆流解后得到参数：K_x/K_y=4.5，等效裂隙半长约65 m，井筒存储系数不可忽略。拟合残差多在±0.03 m以内，表明解能有效捕捉主要响应。基于拟合参数绘制24小时水力响应图，显示降深沿裂隙方向拉长，且优先响应区由s/s_max≥0.05阈值定义，指导额外监测井布设于该区下游方向。

在讨论部分，研究人员强调该模型不直接模拟污染物传输，其水力响应区不可等同于污染羽；它仅作为快速筛选工具，优先用于监测设计，而非替代场地特异性数值模拟。实际应用需结合水样采集、地层资料、示踪试验等。

研究结论翻译如下：
（1）通过COMSOL数值模拟验证了解的有效性，等势线呈椭圆形，表明裂隙连通井或泄漏点可产生椭圆流；误差分析和网格敏感性测试增强了数学可信度。
（2）恒定速率试验中，无井筒存储时各向异性显著影响早中期流态，有存储时主要影响中期；短时监测可能因存储掩盖而低估裂隙连通性。
（3）当t_D < 1000时，降深曲线受井筒存储影响；泵送时间足够长时存储效应消失，中期响应更适合识别主要水力路径。
（4）恒定水头试验中，无量纲流量随时间递减，各向异性在早中期显著；这可用作固定水头源/汇的理想化表示。
（5）水头应用示例表明，水头记录可转化为降深并拟合解析解，用于估算水力各向异性、等效裂隙长度和监测井布设；该模型不能替代污染物传输模拟，但提供快速筛选工具。因控制方程线性，多井响应可叠加，但复杂裂隙网络或非线性条件需扩展模型。