基于波场延拓的反Q滤波方法比较.docVIP

基于波场延拓的反Ｑ滤波方法比较 俞岱，何志军，孙渊，王颖 （1.长安大学 地质工程与测绘学院，陕西 西安710054；.广东华路交通科技有限公司，广东 广州10420） 摘 要： 波场延拓法是最常用的反Ｑ滤波方法，但其振幅因子的不稳定性致使深层的高频信息放大，从而影响反Ｑ滤波的效果。本文通过推导证明了完全稳定振幅补偿只和 Glim（增益限制）有关难以提高深层地震波分辨率，系统分析了波场延拓反Ｑ滤波算法中不同振幅补偿算子的滤波特性，增益控制法对某一频率之后既不放大也不压制，正则化滤波和迭代滤波法不仅具有较强的抗噪声干扰能力，而且对深层地震波具有较好的振幅补偿能力。通过对理论模型和井间地震实际资料应用不同振幅补偿算子，发现正则化振幅补偿因子和迭代法振幅补偿因子不仅对高频干扰有较强的压制效果，而且对于地震记录的中深层部分和低Ｑ能量衰减严重的地震记录，都有较好的补偿效果。 关键词： 波场延拓 反Ｑ滤波 增益控制法 完全稳定法 正则化法 迭代法振幅补偿算子 0 引言 由于地下介质的非均匀性和粘弹性，地震波在地下传播时，地层对振幅能量进行衰减，使得地震波能量降低，同时对地震子波相位进行改变，使得子波延续时间变长，从而导致地震资料的主频向低频方向移动，吸收衰减的程度可用Ｑ来表示，即地层的品质因子。反Ｑ滤波是大地吸收衰减的逆过程，通过反Ｑ滤波可有效地补偿因地层吸收衰减的能量。反Ｑ滤波在补偿地下非均匀粘弹介质的吸收衰减时，在一定程度上能提高资料的分辨率，改善同相轴的连续性，同时可以有效压制一些噪声。1982年Hale提出反Ｑ滤波方法，它基于Futterman 数学模型,即衰减因子是有关频率的指数项，和时间成正比和Ｑ值的大小成反比，其相位项是由振幅项经希尔伯特变换得到的,之后，Varela对 Hale算法进行改进，提出了变型 算法，该算法分为频率域和时间域两种计算方式；上述两种方法都是使用级数展开作高频近似补偿的算法，存在运算量大的缺陷。2002年，Yanghua Wang提出了一种基于波场延拓加增益控制的反Ｑ滤波方法，对于层状Ｑ模型来说，采用变量代换，并且将振幅补偿因子分解，利用快速傅里叶加快运算速率；2003年，姚振兴利用粘弹性介质的理论基础，推导出深度域的反Ｑ滤波算法，该算法在一定程度上可以提高资料的分辨率和信噪比，而且可以保留原始资料能量的分布情况；2005年，刘财基于 stokes 方程重新建立吸收衰减方程，从而提出一种新的基于频率域的补偿方法；2006年，Wang将波场延拓反Ｑ滤波算法推广到连续Ｑ变换的模型中，并且利用Gabor谱变换来加快运算速率；2008年，王珺通过求解时间域的Ｑ模型方程，提出一种新的带状矩形解算器反Ｑ滤波算法，运算速率比较快；2009年， 严红勇等对Wang的波场延拓算法进行改进，提出了一种按射线路径进行反Ｑ滤波的算法，可以对叠前数据进行反Ｑ滤波，但运算速率很低；2013年，张瑾基于梁文锦提出的正则化滤波方法，提出了采用正则化滤波算子来代替完全稳定的振幅补偿因子，同时提出了用迭代法来求取反Ｑ滤波中的振幅补偿算子， 并给出迭代法的收敛性。 文中通过推导证明 Wang 的完全稳定振幅补偿算子难以达到提高深层地震波分辨率的需求，系统分析了波场延拓反Ｑ滤波算法中不同振幅补偿算子的滤波特性。通过对理论模型和井间地震实际资料应用不同振幅补偿算子，对比发现正则化振幅补偿因子和迭代法振幅补偿因子不仅对高频干扰有较强的压制效果，而且对于地震记录的中深层部分和低Ｑ能量衰减严重的地震记录，都有较好的补偿效果。 １ 波场延拓反Ｑ滤波基本原理 反Ｑ滤波方程为： ：传播时间；时间间隔；频率；主频；地层的品质，为频域波场，两个指数分别是相位因子 振幅因子 因此Ｑ后的波场为： 。 利用方程（１）和方程（２）就可以得到波场延拓算法的反Ｑ滤波方程。方程（２）即为成像条件，类似在相偏移中的做法 把所有频率值相加得到 ｔ等于零时刻的波场值，最后输出结果。但注意到振幅因子是不稳定的在延拓过程中会产生高频的数据噪声。针对波场延拓反Ｑ滤波振幅补偿算子的高通滤波特性，Wang提出了完全稳定法和增益控制法来压制高频干扰。 张瑾在前人的基础上，将迭代滤波法及正则化方法引入反 Ｑ 滤波算法中，与Wang提出的反Ｑ滤波相比，这两种方法能更好地解决早期反Ｑ滤波振幅补偿不稳定和深层地震记录无法有效补偿的问题。 1.1完全稳定法 Ｗang提出了一种稳定的反Ｑ滤波方法，其振幅补偿算子为： , (3) 。 （） 是稳定因子，为一个固定的值；β（τ,ω）为振幅衰