[发明专利]一种适用于海上斜缆数据的f-p域鬼波压制方法

摘要

本发明公开了一种适用于海上斜缆数据的f‑p域鬼波压制方法，属于海洋勘探技术领域，作为一种端到端的处理方法，输入为斜缆地震数据，输出即为鬼波压制后的地震数据。通过对传统的线性Radon变换改进得到适用于斜缆地震数据的线性Radon变换算子，随后基于改进的Radon变换算子将地震波分解为不同入射方向的平面波，从平面波角度给出了一次波与鬼波的时差，以该时差公式为基础在f‑p域建立了一次波、虚反射与总波场的关系式，鬼波压制问题也就转换为求解线性方程组的问题。本发明不再需要对鬼波进行预测和匹配相减，简化了鬼波压制的处理流程，提高了处理效率；由于不受中间环节的影响，得到处理结果的信噪比更高。

主权项

1.一种用于斜缆数据鬼波压制的方法，其特征在于：包括如下步骤：步骤1：分别对线性Radon变换算子进行时间校正和炮检距校正，包括如下步骤：步骤1.1：将一次波、虚反射沿着传播路径分别延拓至海平面，等效于将检波器置于海平面来接收地震波场，延拓的时差表示如下：其中：Zi为第i个检波点的深度；vw为海水速度；Δhi,j为海平面延拓接收点到实际接收点的水平距离，并表示如下：Δhi,j＝Zi tanθ   (2)；sinθ＝vwpj    (3)；其中：θ为地震波入射方向或反射方向与垂直方向的夹角；pj为射线参数；步骤1.2：计算海面延拓点的炮检距xi+Δhi,j，并替换线性Radon变换算子中的炮检距xi，校正之后表示如下：τpr＝(xi+Δhi,j)pj‑Δτpi      (4)；将τpr代入到Radon变换算子中，得到适用于斜缆数据的线性Radon变换算子，表示如下：Lu＝exp(‑iωτpr)    (5)；步骤2：构建海水中一次波和相应的海平面鬼波之间关系的平面波模型；步骤3：基于步骤2建立的平面波模型，建立斜缆观测系统下一次波与鬼波的时差公式，其是与入射角度、海水速度以及检波器深度相关的，表示如下：其中，z为接收点相对海平面的深度，Δτ为鬼波比一次波晚到达接收点的时间；步骤4：基于步骤1构建的线性Radon变换算子和步骤3构建的时差关系，在f‑p域建立一次波、虚反射与全波场的关系式，进而构建待求解的线性方程组：P+G＝S      (7)；其中：P为一次波波场；G为鬼波波场；S为全波场；通过步骤3构建的一次波与鬼波的时差关系，得到P与G的关系如下：G＝Re‑iωΔτP    (8)；由于Δτ是与角度相关的，因此必须在f‑p域才能对其进行表征，通过步骤1构建的线性Radon算子将(7)式映射到f‑p域，进而得到最终的求解形式：LuP+RLue‑iωΔτP＝LuS      (9)；其中：R为海面反射系数；Lu为改进之后的线性Radon变换算子；通过求解式(9)得到最终的P，即得到只包含一次波数据的地震数据。