[发明专利]一种压制叠前地震资料线性干扰方法、存储介质及设备

权利要求书

1.一种压制叠前地震资料线性干扰方法，其特征在于，包括以下步骤：

S1、读取原始地震炮集资料；

S2、扫描分析原始地震炮集数据X中任一需要压制的线性干扰波视倾角信息，沿线性干扰波视倾角方向拉平干扰波并沿空间方向对拉平后数据进行随机抽样或者随机重排得到X′；

S3、对X′每道数据沿着时间方向使用一维离散傅里叶变换，将空间方向的每个频点排成Hankel矩阵，对该矩阵进行奇异值分解，保留前三个范围内的奇异值重构数据，对重构的数据做一维逆离散傅里叶变换，得到线性干扰波的估计；将估计得到的线性干扰波数据反拉平后从原始炮集数据中减去，完成对应视倾角方向线性干扰波的压制；

S4、重复步骤S2和S3，直到原始炮集数据中所有需要压制的不同倾角方向的线性干扰波均被处理，压制线性干扰后的炮集数据用于后期速度分析、叠加处理、成像或地质结构分析。

2.根据权利要求1所述的方法，其特征在于，步骤S2中，读取的原始地震炮集数据X包括两部分，X＝X_S+(X_N1+X_N2+…)，X_N1,X_N2…表示不同视倾角的线性干扰波，X_S表示其它信号；

首先，对炮集数据进行视倾角扫描分析，估计X_N1,X_N2…中某一线性干扰波的视倾角方向，沿此方向将X中的线性干扰波同相轴拉平再沿空间方向随机抽样或者随机重排得到新数据X′表达式。

3.根据权利要求2所述的方法，其特征在于，数据X′表达式如下：

其中，x′_n,m为X′中第n行第m列的数据。

4.根据权利要求1所述的方法，其特征在于，步骤S3具体为：

S301、对X′每道数据沿时间方向做一维离散傅里叶变换，得到对应的傅里叶系数矩阵F；

S302、用第i个频率切片数据堆砌Hankel矩阵H_i；

S303、使用奇异值分解对矩阵H_i进行分解；

S304、取前三个范围内的奇异值重构，得到H′_i，对H′_i中数据进行反对角线平均，得到[f′_i,1 … f′_i,m]，重复步骤S302～S304，直到每一个频率切片均被处理，得到F′，对F′的每个空间位置使用一维逆离散傅里叶变换后得到X″，将X″反拉平后得到X_N1,X_N2…中某一线性干扰数据，从原始炮集数据X中减去，重复上述步骤，直到X_N1,X_N2…均被减去，得到压制线性干扰波的结果X_S。

5.根据权利要求4所述的方法，其特征在于，步骤S301中，一维傅里叶矩阵F的表达式如下：

其中，f_n,m为F中第n行第m列的数据。

6.根据权利要求4所述的方法，其特征在于，步骤S302中，矩阵H_i如下：

其中，q为大于0小于m的任意正整数，f_i,m为F中第i行第m列的数据。

7.根据权利要求4所述的方法，其特征在于，步骤S303中，使用奇异值分解如下：

其中，[u₁ … u_q]，λ_γ，γ＝1,…,m-q+1分别为H_i的左右奇异矩阵和奇异值，q为大于0小于m的任意正整数。