[发明专利]一种起伏海面的虚反射压制方法和装置

权利要求书

1.一种起伏海面的虚反射压制方法，其特征在于，所述方法包括：

基于起伏海面反射系数计算公式，对获取的频率波数域的地震数据对应的角频率和波数、与所述频率波数域的地震数据相对应的海浪高度、以及相对应的海水速度进行计算得到反射系数，所述起伏海面反射系数计算公式以预设参数拟合获得；

利用所述频率波数域的地震数据对应的角频率计算得到白噪声系数；

利用所述反射系数、所述白噪声系数分别计算得到预设数量的拖缆扫描深度所对应的频率波数域的虚反射压制算子；

将所述频率波数域的虚反射压制算子分别与频率波数域的地震数据相乘得到预设数量的频率波数域压制数据；

将所述预设数量的频率波数域压制数据进行时间空间参量的二维傅里叶反变换，得到预设数量的时间空间域压制数据；

以预设的时窗间隔划分所述时间空间域压制数据，分别计算同一时窗内所述预设数量的拖缆扫描深度所对应的频域数据的绝对值的和，将绝对值的和最小所对应的拖缆扫描深度的时间空间域压制数据作为相应时窗的压制虚反射后地震数据。

2.根据权利要求1所述的方法，其特征在于，所述起伏海面反射系数计算公式包括如下计算公式：

上式中，r表示反射系数；k_x表示频率波数域地震数据在x方向对应的波数；ω表示频率波数域的地震数据对应的角频率；h表示海浪高度；c表示海水速度；e表示自然常数；π表示圆周率。

3.根据权利要求1或2任意一项所述的方法，其特征在于，所述预设数量的拖缆扫描深度包括：

根据预设扫描范围、扫描间隔和获取的检波器深度确定出的拖缆扫描深度。

4.根据权利要求1或2任意一项所述的方法，其特征在于，所述利用所述频率波数域的地震数据对应的角频率计算得到白噪声系数包括采用下述公式计算：

上式中，white表示白噪声系数；ω表示频率波数域的地震数据对应的角频率。

5.根据权利要求1或2任意一项所述的方法，其特征在于，所述利用所述反射系数、所述白噪声系数分别计算得到预设数量的拖缆扫描深度所对应的频率波数域的虚反射压制算子包括采用下述公式计算：

上式中，A(ω,k_x)_z表示拖缆扫描深度为z时频率波数域的虚反射压制算子，z＝(i-1)*d_z，i表示扫描序号，i取值范围为[0，n]；n表示扫描次数，z₁表示获取的拖缆扫描的起始深度；z₂表示获取的拖缆扫描的终止深度；d_z表示预设扫描间隔；white表示拖缆扫描深度为z时的白噪声系数；G(ω,k_x,r)_z表示拖缆扫描深度为z时包括虚反射的地震波传播算子，表示G(ω,k_x,r)_z的共轭；r表示反射系数；k_x表示频率波数域地震数据在x方向对应的波数；ω表示频率波数域的地震数据对应的角频率。

6.根据权利要求1或2任意一项所述的方法，其特征在于，所述获取的频率波数域的地震数据对应的角频率和波数包括：

对采集的地震数据中单炮地震数据分别进行时间空间参量的二维傅里叶变换后得到的频率波数域的地震数据对应的角频率和波数。

7.根据权利要求1或2任意一项所述的方法，其特征在于，所述预设数量的拖缆扫描深度所对应的频域数据包括：

对所述预设数量的时间空间域压制数据进行一维傅里叶变换后得到的数据。