江雨濛,曹思远,陈思远,马敏瑶,郑 铎,黄 芳,曹国明

(1.中国石油大学(北京)油气资源与探测国家重点实验室,北京102249;2.中国石油大学(北京)地球物理学院,北京102249;3.中国电子科技集团第二十九研究所,四川成都610036;4.中国石油大港油田公司,天津300280)

地震波在地下介质中传播时,由于地层吸收衰减和噪声干扰等因素的影响,子波高频成分能量衰减更快,因此实际采集的地震数据往往是非稳态的,且地震子波的频谱形态会随传播时间变化。由于地震子波是反射地震勘探的重要参数[1-2],也是地震资料处理和解释的基础,在很多环节中都起着关键作用,因此地震子波提取方法的精度直接影响后续地震资料处理及解释的精度和准确性。传统的地震子波提取方法大多基于时不变子波的假设条件,忽略了地震数据的非稳态特征,使得子波估计结果存在误差,进而限制了后续反演结果的分辨率。为了满足高精度地震勘探的需求,如何从非稳态地震数据中准确提取时变子波成为近年来研究的热点之一[3-4]。

针对上述问题,前人提出了基于分段处理的时变子波估计方法[5-10],即首先将地震记录划分为若干段,并假设每一段内地震子波是稳态的,进而结合不同方法提取每一段的时不变子波。此类方法依赖于分段方法和分段长度的选择,且每段提取的平均意义上的子波与实际子波存在一定的误差,使得后续反演结果不精确[11]。近年来,为了消除时窗长度对提取子波的影响,提高子波估计精度,王蓉蓉等[12-13]和张漫漫等[14]相继提出基于时频谱模拟估计时变混合相位子波,实现逐点提取子波。……