[发明专利]一种海洋地震资料震源鬼波压制方法和系统

权利要求书

1.一种海洋地震资料震源鬼波压制方法，包括以下步骤：

步骤一：根据炮检互换原理交换炮检位置，将检波器鬼波压制后的地震资料从共炮点道集转换到共检波点道集，得到转换后的地震资料；

步骤二：判断地震资料的震源激发方式：

若是上下震源激发，计算转换后的地震资料的垂向导数，将转换后的地震资料及其垂向导数从时间域变换到频率域；

若非上下震源激发，对转换后的地震资料进行波场延拓，得到延拓后的地震资料，计算延拓后的地震资料的垂向导数，将延拓后的地震资料及其垂向导数从时间域变换到频率域；

步骤三：对频率域地震资料及其垂向导数进行震源鬼波压制；

步骤四：将震源鬼波压制后的地震资料从频率域变换到时间域，得到最终的震源鬼波压制地震资料。

2.根据权利要求1所述的海洋地震资料震源鬼波压制方法，其特征在于，

所述步骤二中，若非上下震源激发，通过下式对转换后的地震资料进行波场延拓，得到另一个不同深度震源位置的地震资料，作为延拓后的地震资料：

其中，P_R(r'_g,r',ω)为经波场延拓得到的另一个不同深度震源位置的检波器鬼波压制后的地震资料，∫∫_s.s.为在震源激发面上的面积分，P_R(r'_g,r,ω)为检波器鬼波压制后的地震资料，为满足二重狄利克雷边界条件的格林函数，其在海平面和检波器测量面上的值为零，为检波器垂向导数算子，为震源垂向导数算子，dS为面元，r'_g为预测点位置矢量，r为检波点位置矢量，r'为另一个不同深度震源位置矢量，ω为圆周频率。

3.根据权利要求1所述的海洋地震资料震源鬼波压制方法，其特征在于，

所述步骤三中，通过下式对频率域地震资料及其垂向导数进行震源鬼波压制：

其中，P_SR(r'_g,r′_s,ω)为震源鬼波压制后的地震资料，∫∫_s.s.为在震源激发面上的面积分，P_R(r'_g,r,ω)为检波器鬼波压制后的地震资料，为检波器垂向导数算子，G₀(r,r′_s,ω)为震源到检波点的格林函数，dS为面元，r'_g为预测点位置矢量，r为检波点位置矢量，r′_s为预测震源位置矢量，ω为圆周频率。

4.根据权利要求1至3中任意一项所述的海洋地震资料震源鬼波压制方法，其特征在于，采用傅里叶变换将延拓后的地震资料及其垂向导数从时间域变换到频率域。

5.根据权利要求1至3中任意一项所述的海洋地震资料震源鬼波压制方法，其特征在于，采用反傅里叶变换将震源鬼波压制后的地震资料从频率域变换到时间域。

6.一种海洋地震资料震源鬼波压制系统，其特征在于，包括：

第一转换模块，其用于根据炮检互换原理交换炮检位置，将检波器鬼波压制后的地震资料从共炮点道集转换到共检波点道集，得到转换后的地震资料；

震源判断模块，其用于判断地震资料的震源激发方式，其中：

若是上下震源激发，所述震源判断模块计算转换后的地震资料的垂向导数，将转换后的地震资料及其垂向导数从时间域变换到频率域；

若非上下震源激发，所述震源判断模块对转换后的地震资料进行波场延拓，得到延拓后的地震资料，计算延拓后的地震资料的垂向导数，将延拓后的地震资料及其垂向导数从时间域变换到频率域；

鬼波压制模块，其用于对频率域地震资料及其垂向导数进行震源鬼波压制；

第二转换模块，其用于将震源鬼波压制后的地震资料从频率域变换到时间域，得到最终的震源鬼波压制地震资料。