[发明专利]一种海上变深度缆地震数据鬼波压制的方法

权利要求书

1.一种海上变深度缆地震数据鬼波压制的方法，其包括以下步骤：

1)对变深度缆接收的地震数据进行抽道集，形成共中心点道集；

2)将共中心点地震数据进行动校正，并将动校正后的共中心点地震数据进行傅里叶变换，得到含有鬼波的原始地震数据的频谱；

3)根据检波点的不同深度，将得到的含有鬼波的原始地震数据的频谱进行简化，得到H；

4)采用最小二乘谱拟合方法计算得到压制鬼波后的地震数据的频谱；

5)将压制鬼波后的地震数据的频谱进行傅里叶逆变换，得到压制鬼波后的地震数据。

2.如权利要求1所述的一种海上变深度缆地震数据鬼波压制的方法，其特征在于：所述步骤3)中的H(f)＝1-e^-i2πfτ，其中，τ为地震波从检波器到水表的旅行时。

3.如权利要求1或2所述的一种海上变深度缆地震数据鬼波压制的方法，其特征在于：所述步骤4)采用最小二乘谱拟合方法计算得到压制鬼波后的地震数据的频谱具体过程为：

鬼波压制相当于求解如下公式的极小值；

||HS^r-S^*||  (1)

其中，S^*是已知的变深度缆实际采集到的地震数据的矩阵表示，S^r是地下反射回来的地震信号的矩阵表示；

在公式(1)基础上加入体现不同道数据相互补充的约束项，则公式(1)变为：

||HS^r-S^*||+α||MS^r||  (2)

其中，||MS^r||为横向约束条件，α为约束因子；

在公式(2)的基础上加入稳定泛函得：

||HS^r-S^*||+α||MS^r||+β||IS^r||  (3)

其中，β表示稳定因子；

求解公式(3)的极小值，对公式(3)关于S^r求导数并令导数等于零，得：

H^T(HS^r-S^*)+αM^TMS^r+βIS^r＝0  (4)

整理得：

(H^TH+αM^TM+βI)S^r＝H^TS^*  (5)

根据地震数据的特点给定α和β参数，将α和β代入公式(5)得到一个方程组，解这个方程组可以求解得到S^r，得到压制鬼波后的地震数据的频谱。

4.如权利要求3所述的一种海上变深度缆地震数据鬼波压制的方法，其特征在于：α和β两个参数的取值至关重要，α和β根据地震资料的信噪比来确定，假如地震资料的信噪比为b，则α取为β取为