[发明专利]基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法

摘要

本发明涉及一种基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法，首先结合矩阵的克罗内克积性质，改写3D鬼波压制算子形式，将二维大型矩阵反演问题简化为两次小型矩阵反演；然后，将预动校正条件下的三维海上拖缆资料变换至3D Radon域，估计鬼波延迟时间的搜索范围，逐道计算不同延迟时间对应峰值振幅因数，确定压制鬼波的最佳延迟时间；最后，利用构建算子进行鬼波压制。显著提高了3D Radon域鬼波压制方法的运算效率及效果；结合Radon域波场会聚性及峰值振幅因数优化策略，实现了运算效率与压制效果的平衡，同时对起伏海面的资料性有良好的适应性；最终的运算时间是常规算法的20余倍，这有利于3D鬼波压制算法在工业界的广泛应用。

主权项

1.一种基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法，其特征在于，包括以下步骤：a、输入地震资料、简单观测系统信息以及Radon变换参数；b、对数据沿inline、crossline方向进行预动校正处理；c、利用三维快速抛物线Radon变换将地震资料变换至抛物线Radon域，使用频率域三维Radon变换：式中，M和D代表时空域和Radon域数据的频率切片矩阵，和代表inline、crossline方向Radon变换算子，和为对应算子的共轭转置；d、确定鬼波延迟时间的搜索范围，计算不同延迟时间求得的伪分离记录对应的峰值振幅因数，式中，m_i(q_x,q_y,t)表示第i个鬼波延迟时间求得的一次波数据，rme表示m_i(q_x,q_y,t)的均方差,q_x、q_y、t表示Radon域数据的坐标；e、利用峰值振幅因数最大值对应的鬼波延迟时间Δ_op构建鬼波压制算子并将其变换回时空域，使用式中，表示时空域一次波记录的频率切片，表示角频率；f、对时空域一次波记录进行反动校正，得到最终拓频结果。