最佳混合相位反滤滤器的估计

作     者：Bjorn Ursin 冯诗齐 

作者机构：不详 上海海洋局计算中心，上海浦东200120 

出 版 物：《石油物探译丛》 

年 卷 期：2001年第6期

页      面：68-76页

学科分类：081801[工学-矿产普查与勘探] 081802[工学-地球探测与信息技术] 08[工学] 0818[工学-地质资源与地质工程] 

主　　题：混合相位反滤波器 子波估计 反演 反褶积 YW方程组 地震勘探 

摘      要：反滤波用于地震数据以消除子波效应和获得反射序列的估计。在许多情况下，子波是未知的，而仅能估算出子波的自相关函数（ACF）。求解Yule-Walker方程能获得符合最小延迟子波的反滤波器。当子波是混合延迟时，这个反滤波器产生的结果较差。通过求解滤波器方程组的主对角线上ACF的延迟量为α的扩展Yule-Walker方程组，可将反滤波器分解为一个有限长度的滤波器与一个无限长度滤波器的褶积。我们在先前的一篇论文中曾提出，在一个混合延迟的反滤波器中，有限长度滤波器是最大延迟的，而无限长度滤波器则是最小延迟的。本文通过分析有限长度滤波器的Z变换多项式的根改进了上述方法。通过改变混合延尺反滤波器单位圆内根的数目，可获得最多2^a个不同的滤波器。将每个滤波器应用于小的数据集（如一个CMP道集），并选择使得输出结果具有最大L^p(p=5)范数的滤波器为最佳滤波器。这是通过提高a的值来获得最终最佳滤波器。从这一最佳波滤器出发，可以很容易构造出反子波，它可作为地震子波的一个估计来使用。新的算法已用于合成子波和气枪子波以检验其性能和验证重构的子波与初始子波的接近程度。本算法还用于海上叠前地震数据，与用量小延迟滤波器的处理结果比较，本方法改善了叠加剖面。