一种高效的频率域正反τ—ρ变换算法与τ—ρ双曲线速度滤波

摘要

正确的点源记录数据平面波分解(τ-p 变换)的数值计算过程包括 Fourier 变换和 Hallkel 变换。文中对 Hankel 变换中的 Bessel 函数作一逼近后,给出一套在频率域买现正反 r-p 变换的公式。此过程可借助于二维快速傅氏变换(FFT)算法来提高变换的运算速度,其效率比时间域的普通倾斜叠加算法高达二十倍以上,且可大幅度地减少对计算机内存的需要量。我们将文中的高效正反τ-p 变换算法与τ-p 双曲线速度滤波(HVF)相结合,编制了一个在频率域实现的带 HVF 滤波的处理模块。理论合成记录与实际数据的处理结果表明,该方法对压制多种干扰及随机噪声十分有效。