基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法技术

本发明专利技术公开了一种基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法。考虑到鬼波参数随时间和空间变化，影响鬼波去除效果。故首先用滑动时间窗口对海洋地震数据进行划分，降低鬼波参数时变的影响；交替枚举鬼波与一次波的延迟时间和振幅差异系数；通过频率慢度域波场延拓算子对时间窗口内的数据进行延拓，再用原始海洋地震数据分别加上和减去延拓一次的记录；以运算所得新记录乘积绝对值之和最小为目标，迭代寻找最优鬼波延迟时间和振幅差异系数，进而压制鬼波；在每个窗口内的数据压制鬼波后，将所有窗口数据整合成一个完整的记录。该方法能有效地压制鬼波，拓宽频带，提高地震资料的分辨率。震资料的分辨率。震资料的分辨率。

【技术实现步骤摘要】
基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法


[0001]本专利技术涉及鬼波压制
，尤其涉及一种基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法。

技术介绍

[0002]海上勘探中，拖缆通常沉放于海水面以下，受海水面强反射的影响，海洋地震资料中地层反射波后会伴随反相位的鬼波。鬼波的陷波效应使得地震记录有效频带变窄，压制鬼波是拓宽地震资料频带的重要途径。
[0003]在频率慢度域鬼波压制方法中，在海面水平的条件下，检波器的波场在同一慢度时，鬼波与一次波的延迟时间是检波器深度的函数，且因为海水与空气巨大的波阻抗差异，常将海面反射系数假定为
‑
1。通过同一慢度下鬼波与一次波的关系建立两者的关系式，再以最小二乘法对关系式进行求解得到上行波场。但该方法中鬼波与一次波的延迟时间和振幅差异系数在实际采集中有较大限制，实际资料中延迟时间和振幅差异系数是浮动的。
[0004]因此，如何精确获取鬼波与一次波的延迟时间和振幅差异系数是亟待解决的技术问题。

技术实现思路

[0005]为了解决上述技术问题，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法，其特征在于，包括以下步骤：S1：确定滑动时间窗口大小，对海洋地震数据进行划分；S2：将窗口内海洋地震数据中鬼波和一次波的延迟时间和振幅差异系数设为鬼波参数，确定延迟时间和振幅差异系数变化范围，并对延迟时间和振幅差异系数进行等间隔划分；S3：选定延迟时间和振幅差异系数，通过频率慢度域波场延拓算子对窗口内的海洋地震数据进行延拓，再用原始海洋地震数据分别加上和减去延拓一次记录，以运算所得新记录乘积绝对值之和为目标函数，并计算目标函数值；S4：保持振幅差异系数不变，对划分的所有延迟时间重复步骤S3，取目标函数最小值对应的延迟时间为当前最优延迟时间，而后更新延迟时间；S5：保持延迟时间不变，对划分的所有振幅差异系数重复步骤S3，取目标函数最小值对应的振幅差异系数为当前最优振幅差异系数，而后更新振幅差异系数；S6：重复步骤S4至S5，更新最优延迟时间和最优振幅差异系数，直至达到最大迭代次数，以最终得到的最优延迟时间和最优振幅差异系数去除鬼波；S7：对每个窗口内的海洋地震数据重复步骤S3至S6，最后将所有窗口内的数据整合成一个完整的记录，即得到鬼波压制后的海洋地震数据。2.根据权利要求1所述的基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法，其特征在于，步骤S1具体包括：确定滑动时间窗口大小，每个矩形包含部分地震数据，窗口沿着时间轴将整个海洋地震数据离散化，尽量将鬼波和上行波分在同一个时窗，前后滑动窗口之间部分重叠，保证海洋地震数据的全覆盖。3.根据权利要求1所述的基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法，其特征在于，步骤S2具体包括：将窗口内海洋地震数据中鬼波和一次波的延迟时间和振幅差异系数设为鬼波参数，确定延迟时间间隔为α，延迟时间变化范围为Δτ
min
‑
Δτ
max
，振幅差异系数间隔为β，振幅差异系数变化范围为R
min
‑
R
max
，对变化范围内的延迟时间和振幅差异系数以给定间隔进行等距划分。4.根据权利要求1所述的基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法，其特征在于，步骤S3具体包括：S31：将窗口内的海洋地震数据d(x,t)进行一维傅里叶变换到D(x,ω)，再经频率域Radon正变换到S32：根据频率慢度域鬼波压制方法获取频率慢度域波场延拓算子，具体包括：在频率慢度域鬼波压制方法中，对于每一个频率ω有D(x,ω)＝G(x,p,ω)U(p,ω)式中，x为检波器偏移距离，p为水平慢度参数，D(x,ω)为总波场d(x,t)的频率空间域记录，是有N个元素的列向量；U(p,ω)为一次波u(x,t)的频率慢度域记录，是有M个元素的列向量；G(x,p,ω)为鬼波滤波算子在频率慢度域的表示形式，是N
×
M的矩阵，其元素为
式中，G
U
(x,p,ω)与U(p,ω)相乘得到一次波，U(p,ω)和
‑
G
G
(x,p,ω)相乘得到鬼波，θ
j
是第j个水平慢度p
j
对应平面波的出射角度，j＝1,2,...,M，M为慢度个数；Δτ为鬼波与一次波的延迟时间；x
i
、z
i
分别为第i个检波器的炮检距离和深度，i＝1,2,...,N，N是记录的道数；R是假定的鬼波与一次波振幅差异系数；V
w
是海水波场传播速度；由D(x,ω)的表达式得D(x,ω)＝(G
U
(x,p,ω)
‑
G
G
(x,p,ω))U(p,ω)将D(x,ω)简写为D，U(p,ω)简写为U，G
G
(x,p,ω)简写为G
G
，G
U
(x,p,ω)简写为G
U
...

【专利技术属性】
技术研发人员：顾汉明，青杰，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：