《Petroleum Science》：Numerical investigation of wave responses in poroelastic borehole acoustic reflection surveys using a mixed-precision method
编辑推荐：
复杂碳酸盐岩储层对油气聚集和地质CO2封存均具有重要潜力，但其井外异常体会产生复杂的波场响应。井孔声反射勘探利用低频偶极源（dipole source）激发地层并接收反射波场，为刻画远离井孔的结构提供了关键途径。Biot理论（Biot t
复杂碳酸盐岩储层对油气聚集和地质CO2封存均具有重要潜力，但其井外异常体会产生复杂的波场响应。井孔声反射勘探利用低频偶极源（dipole source）激发地层并接收反射波场，为刻画远离井孔的结构提供了关键途径。Biot理论（Biot theory）考虑固-流耦合，因此较纯弹性模型更真实地描述此类含流体碳酸盐岩地层；然而，传统孔隙弹性（poroelastic）模拟在长历时、大规模远探测问题中受限于低计算效率和大内存需求。为此，研究人员提出BEHP-CRDS（Biot’s Equation Half Precision-Compressed Residual Dynamic Scaling），一种结合指数缩放（exponential scaling）、压缩残差补偿（compressed residual compensation）和时间相关动态缩放（time-dependent dynamic scaling）的混合精度（mixed-precision）方法，在保持数值精度的同时显著降低内存消耗。利用该方法，在时间域和时频域（time-frequency domain）分析了代表性模型的孔隙弹性响应。对于不同孔隙度的白云岩和石灰岩，高频成分和晚至波包被削弱后，时频能量更集中于较低频率的主响应周围；在两种岩性中，孔隙度增大均削弱弯曲波（flexural wave）振幅并延迟其到达，其中石灰岩更敏感，尤其在振幅衰减方面。以更复杂的白云岩为背景，SH波（SH wave）对井外异常体显示出更强响应，更适合分析。对于平面薄互层（planar thin interbed），时域到达时间和振幅随距离呈现清晰且单调的变化，而时频属性相对稳定，因此其响应在时域更敏感；相比之下，对于球形孔洞（spherical vug），更复杂的散射使时频分析在定位孔洞方面具有明显优势。总体而言，本研究为井孔声反射勘探提供了一种更高效的数值模拟方法，并为复杂碳酸盐岩地层检测提供了有益见解和理论指导。
研究背景方面，油气勘探开发日益转向深层和复杂储层，常规声波测井主要刻画紧邻井孔的地层，其局限日益明显；井孔声反射勘探（borehole acoustic reflection survey）可探测井周数十米范围内的裂缝、孔洞和薄互层等异常，显著扩展声波测井的探测范围和应用场景。已有进展多集中于单极P波（P-wave）和偶极S波（S-wave）井孔声反射勘探、反射波成像处理以及水平井、斜井和套管井环境中的应用，并在弹性波数值模拟框架下形成了关于单极和偶极源辐射接收机制、典型异常响应、不同井配置适用性和反射波成像处理策略的相对系统认识。然而，实际油气储层通常具有明显孔隙结构和流体赋存，固-流耦合会显著影响波传播、衰减和反射，纯弹性波模型不足以完整揭示复杂储层中井孔声反射勘探的响应机制；已有研究虽开始关注孔隙介质中的成像系统特征及波传播衰减，但孔隙弹性井孔声反射响应的系统研究仍有限。特别是，孔隙弹性响应敏感性随井孔至目标距离增大而衰减的定量分析仍缺乏，目标距井超过10 m时尤其如此；异常体几何属性与围岩如何共同控制不同偏移距记录响应仍不清楚；近井异常响应受直达波及其尾波污染程度尚未系统刻画。这些问题直接关系到可探测极限、有效探测半径、反射波提取窗口选择、近井盲区刻画以及复杂储层异常解释可靠性。除物理机制未决外，孔隙弹性井孔声反射勘探数值模拟计算成本高：波场需在井孔附近、围岩和数十米外异常体之间长距离传播，三维模拟通常需要大计算域和细空间离散以分辨弱反射和复杂散射，计算与存储需求大；与弹性波方程相比，孔隙弹性公式涉及更多耦合变量及孔隙度、渗透率、流体模量等参数，进一步增加复杂性和内存需求。因此，有必要发展面向孔隙弹性井孔声反射勘探的高效计算方法。
研究人员开展了BEHP-CRDS混合精度方法研究，以Biot孔隙弹性波方程为控制方程，构建可在FP16（16位浮点格式）存储并保持FP32（32位浮点格式）有限差分核运算的混合精度格式，并用于偶极源井孔声反射勘探中复杂碳酸盐岩地层孔隙弹性响应分析。研究结论表明，BEHP-CRDS在保持数值精度的同时显著降低内存消耗，并可为长历时、大规模孔隙弹性井孔声反射模拟提供较高效手段；白云岩和石灰岩中孔隙度变化主要引起波形振幅减弱、晚波包抑制、主频向低频移动及弯曲波到达延迟和振幅衰减，石灰岩对孔隙度变化更敏感；以白云岩为背景时SH波更适合分析井外异常，平面薄互层在时域更敏感，球形孔洞则更适合用时频属性定位。其意义在于为复杂碳酸盐岩地层检测提供高效数值模拟方法和理论指导。
关键技术方法概括不超过250字：研究人员采用低频Biot孔隙弹性波方程描述固-流耦合，在交错网格有限差分（staggered-grid finite-difference）框架下求解；提出BEHP-CRDS，将指数缩放、压缩残差补偿和时间相关全局动态缩放结合，以FP16存储缩放波场并保留FP32运算，并以8位有符号整数残差补偿局部舍入信息；采用PML（perfectly matched layer）吸收边界和CFL稳定性条件；设计6个验证算例，与FP32及RAIM（real-axis integration method）对比，并在四块NVIDIA A100 40 GiB平台上测试内存和运行时间；对偶极源四分量记录采用Choi-Williams（CW）时频分布和归一化时频熵H
norm分析，并通过两个正交水平位移分量相减提取SH波响应。本研究无样本队列来源。
研究结果如下。2.2 BEHP-CRDS理论验证：通过6个覆盖不同模型尺度、源主频、源类型、空间离散、时间步、总迭代长度和激发配置的算例，将BEHP-CRDS与常规FP32及RAIM对比，得出BEHP-CRDS与FP32结果高度一致，与RAIM的微小偏差主要来自有限差分离散而非混合精度存储；在较粗空间网格、较高频率复杂全波形、较大模型长历时以及交叉偶极剪切源长历时条件下均保持数值稳定和精度。内存测试中，在1344×192×592网格、0.01 m空间步长、对应13.4×1.92×5.92 m模型上，BEHP-CRDS将总内存由FP32的52.93 GiB降至30.16 GiB；相对原BEHP方案内存略增，但以FP32为参考的最大绝对误差由6.31×10
?3降至5.00×10
?4，降低一个数量级以上。性能测试中，328×128×128网格、10个时间步下，串行方法需934.94 s，并行FP32为7.29 s，BEHP为11.69 s，BEHP-CRDS为7.65 s，说明其在显著节省内存的同时保持与FP32相近的运行时间。
3.1 数值模型与模拟参数：采用三维垂直井模型，井轴平行z方向，水充填井孔半径0.10 m，网格间距0.01 m，PML厚度0.40 m，时间步0.9×10
?6 s；使用主频3 kHz偶极声源，源距3.2 m，对应XMAC-F1仪器配置，沿井孔提取四分量记录用于分析。背景考虑白云岩和石灰岩两类代表性碳酸盐岩，薄互层和孔洞模拟为水充填异常；水充填异常可用相应流体性质参数的弹性波方程近似。由此建立后续孔隙度、孔洞距离和薄互层距离响应分析的基础模型。
3.2 白云岩和石灰岩孔隙度效应：通过比较孔隙度0.05至0.30的XX分量记录及时频属性，得出两种岩性中孔隙度增大均逐渐削弱波形振幅并抑制晚波包，低孔隙度时振幅较大且晚波包更明显，高孔隙度时主事件后快速衰减；这与Biot公式中孔隙度改变孔隙弹性系数、控制应力—应变关系、孔隙压力响应和固-流惯性耦合相一致。CW时频属性显示主响应随孔隙度增大向低频移动，峰值频率由约5 kHz降至约4.2–4.4 kHz，归一化时频熵总体随孔隙度降低，表明高频成分和晚波包削弱后剩余时频能量更集中于较低频主响应。弯曲波结果表明，两种岩性中孔隙度增大均削弱弯曲波振幅并延迟到达；线性拟合的归一化变化率显示，石灰岩到达时间约0.862%/0.05孔隙度，高于白云岩约0.678%/0.05孔隙度；振幅变化率石灰岩约?16.92%/0.05孔隙度，白云岩约?9.75%/0.05孔隙度，负值表示随孔隙度增加而降低，说明石灰岩弯曲波响应对孔隙度变化更敏感，尤其振幅衰减；岩性差异主要体现于弯曲波到达和振幅特征，并更多受剪切波性质差异控制而非较小P波速度差。
3.3 孔洞距井距离效应：在孔隙度10%的白云岩中设置半径0.10 m球形孔洞，孔洞中心沿井孔方向位于源上方1 m，并使其径向远离井轴，井孔至孔洞距离为1–8 m；通过四分量记录提取SH波响应并分析到达时间和振幅。硬碳酸盐岩中，5%–30%孔隙度范围内同一孔洞距离下最大SH波振幅平均约为最大SV波振幅的2.5倍，且P波和SV波到达在时间和振幅上重叠、难以可靠分离，因此采用SH波分析。结果表明，随孔洞远离井壁，SH反射波到达时间逐渐延迟，振幅因几何扩散、散射损失和较长传播路径衰减而减小，较大距离时振幅变化趋缓；CW属性显示SH反射能量随距离减弱并向晚时移动，散射响应峰值频率主要分布在约5 kHz而非源主频3 kHz；2–6 m距离内CW峰值频率和归一化时频熵变化相对缓慢，较大距离时反射能量弱、可检测性受低信号水平限制。由此得出SH波响应可保留清晰的距离相关到达延迟和振幅衰减，适合估计硬地层中小型井外孔洞的径向位置。
3.4 薄互层距井距离效应：在相同孔隙度10%地层中设置厚0.1 m、平行井轴的垂直薄互层，使其径向远离井轴，井孔至薄互层距离为1–8 m，并同样提取SH波响应。波场显示穿过薄互层的波被强烈衰减，同时在薄互层处发生散射。SH波结果表明，随薄互层远离井孔，反射到达时间几乎线性增加，振幅在近距离较高后逐渐下降；CW属性显示最近距离时峰值频率较高且归一化时频熵较低，随距离增大峰值频率下降后稳定在较低水平，归一化时频熵略升后近于常数，说明平面薄互层反射SH波能量相对集中，不随距离显著分散，区别于球形孔洞的复杂散射。相同径向距离下，薄互层SH反射振幅为数倍于球形孔洞，原因是平面薄互层在扩展界面上相干反射入射波场，而小孔洞近似点状散射体；因此平面薄互层更易检测、有效探测半径更大，而小孔洞可检测性更快受振幅衰减限制。平面薄互层响应在时域更敏感，球形孔洞则更适合用CW时频分析研究几何位置，为区分扩展平面反射体与局部散射体提供线索。
讨论与结论部分，原文结论指出：为实现井孔声反射勘探长历时、大规模数值模拟并研究孔隙流体对复杂碳酸盐岩地层响应的影响，研究人员提出BEHP-CRDS，这是一种在交错网格有限差分上求解Biot孔隙弹性波方程的高效混合精度方法；通过结合指数缩放、压缩残差补偿和时间相关动态缩放，该方法扩展了基于Biot的公式，并在不同源频率、模型尺寸、模拟时长和源类型的四个算例中产生准确稳定结果。在1344×192×592网格上，BEHP-CRDS将常规FP32方案的总内存消耗由52.93 GiB降至30.16 GiB；与原BEHP方案相比，其内存略增但更稳定，长历时模拟中最大绝对误差降至5.00×10
?4，比BEHP低一个数量级以上。随后考察白云岩和石灰岩孔隙弹性响应：随孔隙度增加，波形振幅减弱且晚波包更明显，高孔隙度时主波后响应快速衰减；时频域峰值频率降低、能量向低频集中，与Biot系数物理意义一致。XX分量提取的弯曲波显示两种岩性中到达时间随孔隙度增加，石灰岩总体增幅更大（约0.862%/0.05孔隙度），振幅随孔隙度降低，石灰岩变化率约?16.92%/0.05孔隙度，表明石灰岩弯曲波对孔隙度更敏感。最后固定岩性为白云岩，分析本研究模拟条件下井外孔洞和薄互层响应：在本研究考虑的硬碳酸盐岩、3 kHz偶极源和5%–30%孔隙度范围内，SH波比P波和SV波更适合研究远探测井外异常。对单个球形孔洞，SH波反射到达时间随距离逐渐延迟、振幅随距离减小，其较复杂散射产生更分散的时频分布，使时频分析在定位孔洞方面具有明显优势；对单个垂直薄互层，SH波反射到达时间随距离几乎线性增加，与传播距离增加一致，振幅近距离较高后逐渐下降，反射SH波能量相对集中且不随距离显著分散，与孔洞更分散的时频分布形成对比。在本研究测试的井孔至异常1–8 m距离范围内，两类异常均保持近线性到达时间增加和稳定远场特征，因此SH波分量是估计其径向位置的可靠依据；平面薄互层产生更强、更相干且衰减更慢的反射，更易检测并具有更大有效探测半径，而小孔洞可检测性更快受振幅衰减限制。对于所考虑的具体异常几何，薄互层在时域更敏感，孔洞更适合用时频属性刻画，为区分扩展平面反射体与局部散射体提供实用线索。总体而言，该研究为井孔声反射勘探提供更高效数值模拟方法，并为复杂碳酸盐岩地层检测提供有益见解和理论指导。
