基于相对地质时间与多道约束克里金优化的地震波形相似性反演方法

《Petroleum Science》：A Seismic Waveform Similarity Inversion Method Based on Relative Geological Time with Multi-Trace Constrained Kriging Optimization

　　

地质统计学反演是石油勘探开发中的关键技术之一，能够实现有效的高分辨率储层表征。传统的地质统计学反演方法依赖变差函数（variogram）来描述空间相关性，限制了地震信息的充分利用，常导致反演结果横向连续性差。近年来，基于地震波形相似性的反演方法通过用波形相似性替代变差函数，在一定程度上增强了对地震信息的利用。然而，这些方法通常采用逐道独立匹配策略，缺乏有效的空间结构约束，导致地震数据信噪比（SNR）较低时反演结果稳定性不足，且在复杂地质条件下横向不连续。为解决这些问题，研究人员提出了一种基于相对地质时间（RGT）约束的多道联合克里金优化地震波形相似性反演方法。该方法利用RGT数据体表征等时地层结构，并将空间结构约束引入波形匹配过程；同时，将传统的逐道独立匹配方案扩展为多道联合反演框架，提高了波形匹配的稳定性并增强了地质一致性。此外，为增强相似性度量的鲁棒性，研究人员引入了多种波形属性构建混合相似性度量，实现对地震数据及其包络属性的综合表征。最后，将多道联合约束嵌入克里金插值框架，实现了多源信息的协同整合，用于地震波形相似性反演。合成数据和实际数据实例均表明，该方法改善了波形匹配性能、横向连续性和反演稳定性，为高分辨率储层反演提供了有效解决方案。

论文解读文章

**研究背景与问题**

在石油勘探开发中，地震反演是储层预测和流体识别的核心手段，通过从地震数据中提取弹性参数来表征地下储层的空间分布与物理性质。随着勘探精度要求的提高，对地震反演的分辨率、可靠性和稳定性提出了更高要求。传统地震反演方法分为确定性反演和统计反演两大类。确定性反演将模型参数视为确定未知量，依赖正则化技术提升稳定性，但受限于地震分辨率和频带宽度，难以满足薄层勘探需求。统计反演基于贝叶斯理论，通过结合先验信息和似然函数得到后验概率分布，能够更好地表征不确定性，但严格遵循经典贝叶斯框架并采用采样方法遍历解空间在计算上往往不可行。地质统计学反演作为统计反演的重要分支，将地质统计学方法与地震反演结合，有效融合弹性参数的相关性和概率分布等先验地质信息，但严重依赖准确的变差函数建模，计算量大，且存在随机性强、横向连续性差、建模伪影明显等问题，同时未能充分利用井间地震数据中的丰富信息。

基于地震波形相似性的反演方法（WFSI）通过用波形相似性替代传统变差函数，在一定程度上克服了上述局限，减少了建模伪影并增强了横向连续性。然而，现有波形相似性反演方法大多基于逐道独立匹配，在低信噪比或复杂地质条件下缺乏有效的空间结构约束，导致匹配过程不稳定，且未充分利用多道联合优化，反演结果的横向连续性和地质一致性难以保证。

**研究内容与结论**

研究人员提出了一种基于相对地质时间（RGT）约束的多道联合克里金优化地震波形相似性反演方法。该方法的核心思想是利用RGT数据体表征等时地层结构，并将其转化为空间结构约束：一方面，RGT用于定义波形匹配窗口，确保匹配过程沿等时地层界面进行；另一方面，基于RGT等时面信息，使相邻道的垂直影响范围与目标道保持一致，从而在多道联合克里金优化反演中保持垂直一致性，增强道间耦合，抑制跨层干扰。在相似性度量方面，研究人员结合地震振幅和包络属性构建了混合波形相似性度量，提高了波形匹配的稳定性和鲁棒性。在克里金插值框架内，引入压缩感知方法和多道联合约束，基于混合波形相似性选择用于插值的井段并计算相应权重，实现弹性参数的空间估计。反演输出作为随机反演的初始模型，在Metropolis-Hastings（M-H）采样框架内迭代更新，优化弹性参数估计。

合成数据和实际数据测试表明，该方法在波形匹配性能、横向连续性和反演稳定性方面均优于传统地质统计学反演（GSI）和波形相似性反演（WFSI）方法，为高分辨率储层反演提供了有效解决方案。论文发表在《Petroleum Science》。

**主要技术方法**

研究人员采用的关键技术方法包括：基于优化有限差分（OFD）的局部倾角估计与预测绘画法构建RGT数据体；结合皮尔逊相关系数与曼哈顿距离的混合相似性度量，并融合地震包络属性构建双源混合相似性度量；基于RGT约束的多道联合克里金优化框架，引入压缩感知正则化与多道横向一致性约束，采用交替方向乘子法（ADMM）求解权重；将波形相似性反演结果作为先验信息，在贝叶斯框架下通过Metropolis-Hastings MCMC（MH-MCMC）随机反演进行弹性参数最终估计。合成数据测试采用Marmousi2模型局部区域，包含970个时间采样点和1400道地震道，添加20%随机噪声，均匀选取6口伪井；实际数据采用某油田工区地震数据，时间范围2200–2850 ms，采样间隔2 ms，共326个时间采样点，520条Inline线×1001条Crossline线，共8口井作为约束和验证（其中4井和8井作为盲井）。

**研究结果**

**3. 合成数据测试（Synthetic data test）**
通过Marmousi2模型测试，将所提方法与GSI和WFSI方法对比。结果表明，所提方法在纵波速度、横波速度和密度反演中均能更准确地刻画异常体的位置和几何形态，层边界更清晰，横向连续性更好，尤其在薄层和复杂构造区域优势明显。误差分布显示所提方法误差幅度更小且空间分布更均匀，有效避免了传统方法中的条带状伪影。单道反演结果（第820道）显示所提方法曲线更贴近真实模型，在速度和密度突变界面位置准确，并有效抑制了振荡。全局统计指标（平均绝对百分比误差MAPE和相关系数）表明，所提方法在三个参数上均取得最低MAPE和最高相关系数，且误差曲线更平滑。RGT误差敏感性分析显示，对原始RGT施加不同尺度高斯平滑后，反演精度随平滑尺度增大而逐渐下降，但中等误差下仍能保持可接受的结构连续性和参数精度，表明所提方法对RGT适度不确定性具有一定容忍度。噪声鲁棒性测试表明，在信噪比从20 dB降至5 dB的过程中，所提方法相关系数逐渐下降，但MAPE增长平缓，表现出良好的噪声鲁棒性。

**4. 实际数据测试（Field data test）**
在某油田工区实际地震数据上进行反演，井4和井8作为盲井验证。第一阶段波形相似性反演结果中，GSI方法在反射杂乱区域受噪声影响明显，横向不连续且存在不稳定跳变；WFSI方法有所改善但仍有局部不连续和伪影；所提方法在三个参数的反演中均表现出更高的稳定性和更好的横向连续性。第二阶段MH-MCMC随机反演结果中，基于所提方法初始模型的反演结果在盲井4和盲井8处均表现出更好的结构一致性、横向连续性和整体准确性，相关系数分别为0.84（纵波速度）、0.88（横波速度）和0.84（密度），高于GSI（0.76、0.75、0.78）和WFSI（0.81、0.82、0.82）。盲井8处所提方法相关系数为0.88（纵波速度）、0.78（横波速度）和0.88（密度），同样高于GSI（0.77、0.71、0.78）和WFSI（0.79、0.77、0.79）。计算效率方面，所提方法耗时6533.1秒，高于GSI（1396.6秒）和WFSI（3427.5秒），但综合考虑精度提升，计算成本在可接受范围内。

**讨论与结论**

讨论部分指出，该研究的主要贡献在于将RGT结构约束和多道协同克里金优化策略引入波形相似性反演框架，有效增强了横向连续性和整体稳定性。RGT约束建立了统一的地层时间框架，将参数更新限制在同一沉积层内，使反演更符合沉积和构造模式。多道约束克里金优化加强了道间耦合，使反演参数空间变化更平滑、更符合地质规律。超参数敏感性分析表明，所提方法在合理参数范围内保持相对稳定，最终选取参数为ω₁=0.02，ω₂=0.3，σ₀=3。尽管有改进，仍存在局限性：RGT精度在复杂构造（如断层、地层不连续）中可能产生系统误差；波形相似性反演更适用于相对稳定的沉积环境，多道约束可能降低对局部异常的敏感性，导致强异质性区域过度平滑。未来工作可聚焦于多尺度融合或数据驱动方法提升RGT精度，以及引入多尺度或自适应约束机制以平衡横向连续性与局部异常捕捉能力。

**研究结论**
在相对稳定的沉积环境中，地下构造及其相关弹性参数通常表现出较强的横向连续性。基于这一前提，所提方法将RGT结构约束和多道协同克里金优化策略引入传统波形相似性反演框架，通过整合地质结构信息与空间统计约束，将参数更新过程从简单空间距离驱动转变为地层引导传播，实现了有效的多道信息融合，提高了地质一致性。合成模型和实际数据实验均表明，所提方法在三个参数的反演中取得了良好性能，在横向连续性和整体稳定性方面优于传统方法，且井点反演精度高于GSI和WFSI方法，表明该方法在不牺牲局部精度的前提下增强了空间连续性。总体而言，所提方法为地震波形相似性反演提供了新视角，为储层表征和勘探提供了可靠的技术支持。