[发明专利]一种地震数据鬼波压制方法、装置、设备及介质

权利要求书

1.一种地震数据鬼波压制方法，其特征在于包括：

将含有鬼波的地震数据变换至频率波数域，得到频率波数域地震数据；

根据地震数据采集的炮点深度和检波点深度计算得到鬼波算子；

基于鬼波算子计算零相位化后的频率波数域地震数据；

计算鬼波算子的频谱振幅校正项，并对其进行阈值截断处理得到阈值截断后的鬼波算子频谱振幅校正项；

利用阈值截断后的鬼波算子频谱振幅校正项对零相位化后的频率波数域地震数据进行处理，得到鬼波压制后的频率波数域地震数据；

对鬼波压制后的频率波数域地震数据进行反变换得到鬼波压制后的时空域地震数据。

2.根据权利要求1所述的地震数据鬼波压制方法，其特征在于，将含有鬼波的地震数据通过二维傅里叶FK变换至频率波数域。

3.根据权利要求1所述的地震数据鬼波压制方法，其特征在于，鬼波算子的计算公式为：

式中，G(k_x,ω)为鬼波算子，k_x为水平波数，ω为角频率，R为海表面的反射系数，z_s为地震数据炮点深度，z_r为检波点深度，k_z为垂直波数。

4.根据权利要求3所述的地震数据鬼波压制方法，其特征在于，零相位化后的频率波数域地震数据的计算公式为：

D_zero(k_x,ω)＝D(k_x,ω)*G(k_x,ω)^H

式中，D_zero(k_x,ω)为零相位化后的频率波数域地震数据，D(k_x,ω)为频率波数域的地震数据，H代表共轭转置。

5.根据权利要求4所述的地震数据鬼波压制方法，其特征在于，阈值截断后的鬼波算子频谱振幅校正项的计算公式为：

式中，α为设定阈值。

6.根据权利要求5所述的地震数据鬼波压制方法，其特征在于，鬼波压制后的频率波数域地震数据的计算公式为：

7.根据权利要求1～6任一项所述的地震数据鬼波压制方法，其特征在于，对鬼波压制后的频率波数域地震数据进行反二维傅里叶变换得到鬼波压制后的时空域地震数据。

8.一种地震数据鬼波压制装置，其特征在于，该装置包括：

数据变换单元，被配置为将含有鬼波的地震数据变换至频率波数域，得到频率波数域地震数据；

鬼波算子计算单元，被配置为根据地震数据采集的炮点深度和检波点深度计算得到鬼波算子；

频率波数域地震数据计算单元，被配置为基于鬼波算子计算零相位化后的频率波数域地震数据；

鬼波算子频谱振幅校正单元，被配置为计算鬼波算子的频谱振幅校正项，并对其进行阈值截断处理得到阈值截断后的鬼波算子频谱振幅校正项；

阈值截断单元，被配置为利用阈值截断后的鬼波算子频谱振幅校正项对零相位化后的频率波数域地震数据进行处理，得到鬼波压制后的频率波数域地震数据；

数据反变换单元，被配置为对鬼波压制后的频率波数域地震数据进行反变换得到鬼波压制后的时空域地震数据。

9.一种电子设备，其特征在于，包括：一个或多个处理器、存储器及一个或多个程序，其中一个或多个程序存储在所述存储器中并被配置为所述一个或多个处理器执行，所述一个或多个程序包括用于执行如权利要求1至7所述方法中的任一方法的指令。

10.一种存储一个或多个程序的计算机可读存储介质，其特征在于，所述一个或多个程序包括指令，所述指令当由计算设备执行时，使得所述计算设备执行如权利要求1至7所述方法中的任一方法。