[发明专利]一种在空间‑频率域基于格林理论的鬼波压制方法

权利要求书

1.一种在空间-频率域基于格林理论的鬼波压制方法，包括以下步骤：

S10，对地震波场进行包括去噪的预处理，获得预处理后的地震波场；

S20，根据地震波场的采集方式，计算预处理后的地震波场的垂向导数；

S30，将所述预处理后的地震波场和预处理后的地震波场的垂向导数从时间域变换为频率域；

S40，根据地震波场的采集维度，计算背景介质中的格林函数和格林函数的垂向导数；

S50，根据所述预处理后的地震波场、预处理后的地震波场的垂向导数、背景介质中的格林函数和格林函数的垂向导数，在空间-频率域中基于格林理论从预处理后的地震波场中分离出上行波场，所述从预处理后的地震波场中分离出的上行波场为鬼波压制后的频域波场；

S60，将鬼波压制后的频域波场从频率域变换为时间域，获得鬼波压制后的时域波场。

2.根据权利要求1所述的鬼波压制方法，其特征在于，步骤S10包括以下步骤：

S10.1利用低切滤波、去气泡和压制涌浪噪音的方法对地震波场进行去噪处理；

S10.2利用高分辨率抛物拉东变换法进行波场外推，恢复近偏移距数据。

3.根据权利要求1所述的鬼波压制方法，其特征在于，步骤S20中，所述地震波场的采集方式为双检地震采集或上下缆地震采集。

4.根据权利要求3所述的鬼波压制方法，其特征在于，步骤S20中，所述地震波场的采集方式为双检地震采集，通过下式获得预处理后的地震波场的垂向导数

∂∂zP(r,rs,ω)=iωρVz]]>

其中，P(r,r_s,ω)为预处理后的地震波场，V_z为垂直速度分量，r为检波点位置，r_s为震源位置，ω为圆周频率，ρ为背景介质密度，i为虚数单位。