[发明专利]气枪阵列子波方向性反褶积方法、装置、介质和设备

权利要求书

1.一种基于波场分解的气枪阵列子波方向性反褶积方法，其特征在于，该方法包括：

将共接收点地震数据变换至频率波数域，得到频率波数域地震数据；

在频率波数域计算确定每一频率和每一波数对应的气枪阵列远场子波的出射角，基于所述气枪阵列远场子波的出射角，以及所述气枪阵列的空间位置参数，合成不同方向的远场子波，基于所述不同方向的远场子波，生成频率波数域方向性反褶积滤波器；

基于所述频率波数域方向性反褶积滤波器，对所述频率波数域地震数据进行滤波处理，得到方向性反褶积后的地震数据；

对所述方向性反褶积后的地震数据进行傅里叶逆变换，得到最终的地震数据。

2.根据权利要求1所述的气枪阵列子波方向性反褶积方法，其特征在于，所述基于所述不同方向的远场子波，生成频率波数域方向性反褶积滤波器，包括：

基于所述不同方向的远场子波，构建方向性反褶积滤算子；

将所述方向性反褶积滤算子映射到频率波数域，得到频率波数域方向性反褶积滤波器。

3.根据权利要求2所述的气枪阵列子波方向性反褶积方法，其特征在于，所述基于所述不同方向的远场子波，构建方向性反褶积滤算子，包括：

基于以下公式构建方向性反褶积滤波算子q(ω，θ)：

其中，ε为稳定性因子，取值范围为0.001-0.005；W(ω，θ)为远场子波。

4.根据权利要求1～3任一项所述的气枪阵列子波方向性反褶积方法，其特征在于，所述在频率波数域计算确定每一频率和每一波数对应的气枪阵列远场子波的出射角，包括：

根据以下公式计算每一频率和每一波数对应的气枪阵列远场子波的出射角θ：

其中，v为地震波在海水中的传播速度，f为频率，k为水平视波数。

5.根据权利要求4所述的气枪阵列子波方向性反褶积方法，其特征在于，所述基于所述气枪阵列远场子波的出射角，以及所述气枪阵列的空间位置参数，合成不同方向的远场子波，包括：

根据以下公式，在频率域合成不同方向的远场子波W(ω，θ)：

其中，j为虚数单位；ω＝2πf，为圆频率；θ为远场位置相对于气枪阵列中心的角度即远场子波的出射角；n为气枪个数；p_i(ω)为气枪阵列中第i支气枪的假想子波；为第i支气枪到远场位置上的矢量距离；为第i支气枪其虚震源到远场位置上的矢量距离；v为地震波在海水中的传播速度；σ为海水面虚反射系数。

6.根据权利要求4所述的气枪阵列子波方向性反褶积方法，其特征在于，所述基于所述频率波数域方向性反褶积滤波器，对所述频率波数域地震数据进行滤波处理，得到方向性反褶积后的地震数据，包括：

根据以下公式计算频率波数域方向性反褶积滤波器的反向滤波器H(f，k_x)：

其中，δ为白噪稳定因子，取值范围为0.001-0.005；G(f，k_x)为频率波数域方向性反褶积滤波器；

基于所述反向滤波器H(f，k_x)对所述频率波数域地震数据进行滤波处理，得到方向性反褶积后的地震数据。

7.根据权利要求1所述的气枪阵列子波方向性反褶积方法，其特征在于，对所述方向性反褶积后的地震数据进行傅里叶逆变换，得到最终的地震数据，包括：

对所述方向性反褶积后的地震数据进行二维傅里叶逆变换，将方向性反褶积后的地震数据反变换回时间空间域，以得到最终的地震数据。