利用钻孔热恢复量化地下渗透性的方法

《Geochemistry, Geophysics, Geosystems》：Borehole Thermal Recovery as a Method for Quantifying Subsurface Permeability

　　

地下水利学性质的表征是水文地质学与地质资源应用领域面临的挑战，尤其在裂隙或非均质储层中。本研究探索了开放式钻孔在流体循环后的热恢复能否作为一种原位（in situ土地的用途）方法，用于估算渗透性及可渗透带的空间展布。通过耦合流体流动与热量传递的数值模拟，研究人员分析了水力和热扰动在循环过程中如何驱动流体入渗和对流换热，当水平可渗透带贯穿低渗透性岩石时，这些过程会以依赖于可渗透带渗透性和厚度的方式改变 thermal recovery 行为。研究发现，可渗透带的存在使热恢复时间显著延长，且渗透性大于约10^-14 m²时，可通过恢复时间和表观热导率进行定量估算。该技术在钻井或注入试验相关的循环后，于多个深度采集温度时间序列数据时表现尤为有效。研究结果支持了先前解释的可信度，即恢复时间随深度的变化反映了受流体入渗影响的高渗透带。综上，开放式井中的深度分辨温度监测为估算地下渗透性和可渗透带厚度提供了一种前景广阔且高效的工具。

本研究发表于《Geochemistry, Geophysics, Geosystems》，聚焦于一种创新性的原位渗透性表征方法——水文热响应测试（Hydrologic Thermal Response Test）。地下渗透性作为控制流体流动的关键参数，决定着压力分布、有效应力、热传递、溶质运移以及流体注入或抽采的可行性。然而，传统方法存在明显局限：间接方法如核磁共振（Nuclear Magnetic Resonance, NMR）和电阻率法依赖孔隙度-渗透性关系进行推断，而直接方法如试井测试虽能隔离特定深度段，但操作复杂繁琐。在此背景下，研究人员提出利用开放式钻孔中流体循环后的热恢复信号来定量表征渗透性，为地下渗透性测量提供了一种替代方案。

该研究的科学动机源于Fulton等人在日本海沟快速钻探项目（Japan Trench Fast Drilling Project, JFAST）中的现场观测：在IODP第343/343T航次的海底钻孔中，特定深度出现的热恢复延迟现象被推断与高渗透性断层和裂隙带相关，但这一假设尚未得到定量验证。基于此，研究人员旨在建立一套定量框架，通过分析热恢复模式来识别可渗透带并约束其渗透性和厚度。

研究人员采用二维径向轴对称数值模型，模拟开放式钻孔贯穿低渗透性岩石中的水平可渗透带的热-流耦合过程。模型基于COMSOL Multiphysics v5.4平台构建，涵盖流体循环和后续热恢复两个阶段。循环阶段，钻孔内流体流动产生热扰动，若存在可渗透带，冷却的循环流体在压差驱动下入渗，增强对流换热；恢复阶段，地层逐渐回归背景地温场，通过分布式温度传感器监测不同深度的温度演化。模型关键参数包括：可渗透带渗透性（10^-20至10^-12 m²）、深度（50-750 m）、厚度（0.5-7 m）、循环井压（250-3000 kPa）及循环持续时间（2-12 hr）。流体流动在钻孔内采用纳维-斯托克斯方程（Navier-Stokes equations）描述，在多孔介质中采用Brinkman方程描述；热传递通过包含对流项的对流-扩散方程刻画。表观热导率（apparent thermal conductivity, λ_a）和特征恢复时间（characteristic recovery time, τ_rec）被定义为定量表征热恢复的两个核心指标，前者通过拟合修正的Carslaw-Jaeger导热解获取，后者通过指数函数拟合温度时间曲线获取。

研究结果部分首先验证了数值方法的可靠性。通过与纯导热解析解对比，数值模拟在纯导热条件下展现了良好的一致性，确认了方法的准确性。

可渗透带对钻孔热响应的影响：对比有无可渗透带的模拟结果，发现可渗透带的存在显著扩大了热扰动的空间范围（影响半径约达3 m），并导致恢复时间明显延长。这一结果支持了第一个假设——延长的热恢复时间可用于识别可渗透带。

渗透性对热响应的敏感性：系统测试了12个不同渗透性值（10^-20至10^-12 m²）的影响。结果显示，随渗透性增加，最大冷却幅度增大，扰动区域扩展；表观热导率从1.14降至0.13 W/(m·K)，特征恢复时间从12.5天增至31.4天。两个指标均呈现近阈值行为，在约10^-14 m²处发生显著转变，该阈值对应于入渗驱动对流与导热相当的临界条件。佩克莱数（Péclet number, β），即对流与导热热传输之比，被引入以解释这一行为，高渗透性增加达西流（Darcy flux, q_D），从而强化对流换热。

深度、厚度、井压和循环持续时间的影响：额外24组敏感性分析揭示了各参数的作用机制。深度增加导致更大的流体-地层温差，延长恢复时间并降低表观热导率；厚度增加产生类似效应，且在薄层时尤为敏感，因邻近低渗透性岩石的抑制作用显著；井压增加增大压差驱动更强入渗，延迟恢复；循环持续时间延长则增大热扰动幅度，对短持续时间尤为敏感。

讨论部分阐述了水文热响应测试的实际应用要点。实施步骤包括：在开放式钻孔中进行流体循环，利用光纤分布式温度传感（Distributed Temperature Sensing, DTS）、热敏电阻串或离散微型温度记录仪进行多深度温度监测，结合数值模型生成深度归一化的查找曲线（lookup curves），通过匹配实测恢复时间与表观热导率来估算渗透性和厚度。该方法的优势在于能生成详细的深度分辨渗透性剖面，监测过程被动进行，可同时开展数值模拟与解释。然而，需注意热异常局限于近井区域，而压力扩散距离可能更大；实际操作中存在循环速率、入口温度和井压变化等影响因素；数值模型假设了层流、可渗透带半无限宽等简化条件；地质条件和热扩散率的空间变异也会带来不确定性。

研究结论指出，水文热响应测试作为一种原位方法，利用热作为示踪剂，通过流体循环后的钻孔热恢复来估算渗透性和可渗透带厚度。数值模拟证实了延长的热恢复时间可识别可渗透带，且温度时间序列的反演能定量约束渗透性和厚度，特别确定了10^-14 m²的敏感性阈值。尽管基于特定条件建立的查找曲线需视具体场地调整，但该方法为地下水管理、地热开发等 subsurface 应用中的渗透性表征提供了实用高效的替代工具。