基于等效Q值与重心频率单调关系构建正演算子的非平稳地震反射率反演方法

《Applied Sciences》：Centroid-Frequency-Guided Nonstationary Reflectivity Inversion and Its Application for Top and Bottom Interface Identification of Coalbed Methane Reservoirs

　　非平稳卷积模型被广泛用作非平稳地震反射率反演（NSRI）的正演公式。该正演公式由时变子波与反射率的卷积确定。以往的非平稳反射率反演算法仅关注根据不同的稀疏正则化策略进行反射率估计，正演算子——时变子波矩阵通常使用用户定义的常数Q值来生成。然而，Q值随时间和空间变化，因此传统NSRI方法将产生不准确的反射率反演结果。为了提高反演反射率的精度，我们提出了一种基于等效Q值（Q_e）与重心频率（CF）之间单调关系的新正演算子构建方法。我们首先从地震信号的Gabor时频谱中提取重心频率，然后通过CF-Q_e模板获得Q_e。因此，避免了选择Q值的盲目性，可以构建出更精确的正演算子。合成数据和实际数据实例表明，所提出的方法通过构建更可靠的正演算子，提供了更准确的反射率估计。此外，增强后的反射率剖面能够更清晰地划分煤层气储层的上、下界面，为煤层解释和后续开采规划提供高分辨率地震支持。