基于迭代收缩高分辨率Radon变换的地震道插值

摘要

在地震数据采集过程中,由于受经济预算的限制以及施工条件的影响,沿空间方向常会出现数据缺失、不规则采样或稀疏采样现象.地震数据的缺失或稀疏采样会不同程度地影响地震资料多道处理的结果,尤其对偏移成像、表面多次波压制以及谱估计的影响最为严重.因此,地震数据的插值重建在地震资料处理中显得尤为关键.常见的地震道插值方法有基于预测滤波的插值重建方法,基于数据变换的地震数据插值,如Fourier变换,曲波变换,Radon变换等,基于波动方程的插值重建方法.其中基于高分辨Radon变换(RT)的地震道插值方法因无需精确的NMO速度，并且具有较强的抗噪能力，因此在地震数据插值中得到了广泛的应用(Sacchi等，1995)。本文提出了一种基于迭代收缩高分辨率Radon变换的地震道插值方法，它将2D迭代模型收缩算法代替传统的混合域高分辨率RT中的IRLS算法，在获取高分辨率Radon模型解的基础上，极大地提高了计算效率，并将其用于地震数据的插值重建中获得了更准确的插值结果。rn 分别利用合成模型和实际数据对提出的插值算法进行了测试。通过测试结果可以发现，基于迭代收缩高分辨率radon变换的地震道插值方法，相对于基于LSRT的地震道插值、基于FSRT的地震道插值结果具有更高的准确性，插值的地震同相轴连续性更好;同时它可以用于非规则采样数据的重建，具有明显的抗假频能力，实验结果表明了其相对于传统的f-x预测滤波插值算法，以及快速广义Fourier变换(fgft)插值算法的优越性。