基于正则化非平稳自回归技术的f-x域随机噪声衰减方法.pdfVIP

507 SPG/SEG 深圳2011 国际地球物理会议论文集 基于正则化非平稳自回归技术的f-x 域随机噪声衰减方法 1* 1 2 1 刘国昌 , 陈小宏, 杜婧 , 吴开龙 1. 中国石油大学CNPC 物探重点实验室 2. 中国石化胜利油田物探研究院 摘要在正则化非平稳自回归技术的基础上, 提出了一种新颖的 f-x 域随机噪声衰减方法. 该方法可以自适应的预测 振幅或同相轴倾角随空间变化的地震信号, 从而压制随机噪声. 传统的 f-x 域预测滤波方法采用滑动窗口技术克服 地震数据非平稳的特征, 在本文提出的方法中, 空间方向不需要滑动窗口，而是通过正则化技术完成自适应预测滤 波. 我们约束预测滤波器系数随空间变化且光滑, 采用迭代方法求解最小平方反演问题, 提高了计算速度. 最后用 模拟和实际数据验证了提出方法的有效性. 方法原理 地震随机噪声通过时间-空间域或频率-空间域预测滤波压制 (Abma and Clarebout, 1995), 常规 的预测滤波方法假设地震同相轴倾角恒定不变. Fomel (2009) 提出了时间域正则化非平稳自回归技 术, 本文将其扩展到频率域并应用于随机噪声的衰减. 考虑因果的f-x 域预测滤波器, 我们有 M M  (f ) S ( f )  a S ( f )  a S ( f ) , (1) n n  i n i  i n i i 1 i 1 其中S ( f ) 为f-x 域地震数据, a 为预测滤波器系数. 为修改公式(1)为非平稳滤波器，我们将滤波器 n i 系数修改为an ,i ( f ) , 可以看到修改后的系数随频率和空间变化，显然修改后最小平方预测问题是不适 定的反演问题，因此需要强加约束来求解该问题。加约束后的最小平方问题为 M min || S ( f )   a ( f )S ( f ) ||2 R[a ( f )] . (2) ,  2 n ,i n n i n i an ,i ( f ) i M , i 0 我们采用整形正则化技术迭代求解上述最小平方问题, 针对复数的情况, 在求解过程中强加的约束 为实数型光滑算子, 即复数的实部和虚部分别加光滑约束. 通过改变滤波器系数为空间有关的函数 可以适应同相轴振幅和倾角随空间变化的情况。 我们分别用一个模拟的叠前道集和实际的叠后剖面验证了该方法. 与传统的t-x 和f-x 预测技术 相比，该方法更能提高地震信号的信噪比. 本文方法的计算时间为O