【技术实现步骤摘要】
本专利技术涉及海洋地球物理勘探领域，更具体地，涉及一种自由表面多次波预测方法、装置、电子设备及介质。
技术介绍
1、在海洋地震勘探中，由于海平面的存在，地震波在海平面和地下反射界面间多次反射形成自由表面多次波，自由表面多次波是影响海洋地震资料信噪比和分辨率的主要因素之一。因此，衰减或消除自由表面多次波是海洋地震资料处理中一个关键环节。为了消除多次波的干扰，提高资料分辨率，物探界提出了两类多次波压制方法：一类是基于一次波与多次波之间特性差异的滤波法；另一类是基于波动理论的预测相减法。目前，预测相减法无需先验信息，是压制多次波主要方法。主要包括反馈迭代法(srme)和逆散射级数法(iss)，针对自由表面多次波压制，由于没有去除震源鬼波，反馈迭代法(verschuur等，1992)仅可以预测出多次波准确的时间和近似的振幅，而逆散射级数法完全数据驱动，可以预测出多次波准确的时间和振幅，然后通过直接相减压制多次波。逆散射级数法从weglein等(1997，2003)提出以来得到了迅速的发展，matson等(1997)在海洋资料应用中取得了较好的效果；李翔等(2009)分析了逆散射级数法在实际应用中的子波估计和近道波场恢复问题；yang等(2013，2016)分析了震源子波和辐射方向模式对自由表面多次波预测的影响，改进了逆散射级数法多次波预测的精确度。
2、目前尚需开发一种二维震源一维介质的自由表面多次波预测方法。
3、公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认
技术实现思路
1、本专利技术提出了一种自由表面多次波预测方法、装置、电子设备及介质，本专利技术在地下介质一维近似的情况下，通过算法优化，在保证多次波预测精度的前提下，提高计算效率。改进的算法保持了原有方法的优越性，完全数据驱动，无需已知地下速度模型。模型数据的处理结果验证了本方法可以预测出自由表面多次波的准确相位和振幅，有效地压制自由表面多次波。
2、第一方面，本公开实施例提供了一种自由表面多次波预测方法，包括：
3、对预处理后的地震数据从炮集转换cmp道集；
4、将cmp道集数据从空间时间域变换到波数频率域；
5、计算中心点水平波数，进而计算中心点垂向波数；
6、对各阶自由表面多次波进行预测；
7、对预测的各阶自由表面多次波进行二维傅里叶反变换从波数频率域变回到空间时间域；
8、计算自由表面多次波压制后的地震数据。
9、优选地，所述预处理包括低切滤波、去气泡、压制涌浪噪音、去除直达波、压制检波器鬼波和震源鬼波、近道恢复、规则化。
10、优选地，对cmp道集数据进行二维傅里叶变换从空间时间域变换到波数频率域。
11、优选地，通过公式(1)计算所述中心点水平波数：
12、kh＝(kg+ks)/2 (1)
13、其中，kh为中心点水平波数，ks和kg为震源和检波器的水平波数。
14、优选地，通过公式(2)计算所述中心点垂向波数：
15、
16、其中，q为中心点垂向波数，ω为圆周频率，c0为背景介质中的水速。
17、优选地，通过公式(3)对各阶自由表面多次波进行预测：
18、
19、其中，ω为圆周频率，i为虚数单位，c0为背景介质中的水速，a(ω)为震源子波；kh为震源和检波器的中心点水平波数，q为中心点垂向波数；εs和εg为震源和检波器的深度；d1为输入地震数据，dn为预测的第n-1阶自由表面多次波地震数据。
20、优选地，通过公式(4)计算所述自由表面多次波压制后的地震数据：
21、
22、其中，d为自由表面多次波压制后的地震数据，dn为预测的第n-1阶自由表面多次波地震数据。
23、作为本公开实施例的一种具体实现方式，
24、第二方面，本公开实施例还提供了一种自由表面多次波预测装置，包括：
25、预处理模块，用于对预处理后的地震数据从炮集转换cmp道集；
26、变换模块，用于将cmp道集数据从空间时间域变换到波数频率域；
27、第一计算模块，用于计算中心点水平波数，进而计算中心点垂向波数；
28、预测模块，用于对各阶自由表面多次波进行预测；
29、反变换模块，用于对预测的各阶自由表面多次波进行二维傅里叶反变换从波数频率域变回到空间时间域；
30、第二计算模块，用于计算自由表面多次波压制后的地震数据。
31、优选地，所述预处理包括低切滤波、去气泡、压制涌浪噪音、去除直达波、压制检波器鬼波和震源鬼波、近道恢复、规则化。
32、优选地，对cmp道集数据进行二维傅里叶变换从空间时间域变换到波数频率域。
33、优选地，通过公式(1)计算所述中心点水平波数：
34、kh＝(kg+ks)/2 (1)
35、其中，kh为中心点水平波数，ks和kg为震源和检波器的水平波数。
36、优选地，通过公式(2)计算所述中心点垂向波数：
37、
38、其中，q为中心点垂向波数，ω为圆周频率，c0为背景介质中的水速。
39、优选地，通过公式(3)对各阶自由表面多次波进行预测：
40、
41、其中，ω为圆周频率，i为虚数单位，c0为背景介质中的水速，a(ω)为震源子波；kh为震源和检波器的中心点水平波数，q为中心点垂向波数；εs和εg为震源和检波器的深度；d1为输入地震数据，dn为预测的第n-1阶自由表面多次波地震数据。
42、优选地，通过公式(4)计算所述自由表面多次波压制后的地震数据：
43、
44、其中，d为自由表面多次波压制后的地震数据，dn为预测的第n-1阶自由表面多次波地震数据。
45、第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
46、存储器，存储有可执行指令；
47、处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的自由表面多次波预测方法。
48、第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的自由表面多次波预测方法。
49、本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。
本文档来自技高网...【技术保护点】
1.一种自由表面多次波预测方法，其特征在于，包括：
2.根据权利要求1所述的自由表面多次波预测方法，其中，所述预处理包括低切滤波、去气泡、压制涌浪噪音、去除直达波、压制检波器鬼波和震源鬼波、近道恢复、规则化。
3.根据权利要求1所述的自由表面多次波预测方法，其中，对CMP道集数据进行二维傅里叶变换从空间时间域变换到波数频率域。
4.根据权利要求1所述的自由表面多次波预测方法，其中，通过公式(1)计算所述中心点水平波数：
5.根据权利要求4所述的自由表面多次波预测方法，其中，通过公式(2)计算所述中心点垂向波数：
6.根据权利要求5所述的自由表面多次波预测方法，其中，通过公式(3)对各阶自由表面多次波进行预测：
7.根据权利要求6所述的自由表面多次波预测方法，其中，通过公式(4)计算所述自由表面多次波压制后的地震数据：
8.一种自由表面多次波预测装置，其特征在于，包括：
9.一种电子设备，其特征在于，所述电子设备包括：
10.一种计算机可读存储介质，其特征在于，该计算机可读存
...【技术特征摘要】
1.一种自由表面多次波预测方法，其特征在于，包括：
2.根据权利要求1所述的自由表面多次波预测方法，其中，所述预处理包括低切滤波、去气泡、压制涌浪噪音、去除直达波、压制检波器鬼波和震源鬼波、近道恢复、规则化。
3.根据权利要求1所述的自由表面多次波预测方法，其中，对cmp道集数据进行二维傅里叶变换从空间时间域变换到波数频率域。
4.根据权利要求1所述的自由表面多次波预测方法，其中，通过公式(1)计算所述中心点水平波数：
5.根据权利要求4所述的自由表面多次波预测方法，其中，通过公式(2)计算所述中...
【专利技术属性】
技术研发人员：杨金龙，张印堂，裴云龙，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。