[发明专利]一种基于MVMD和MSSA的三维地震数据随机噪音压制方法

权利要求书

1.一种基于MVMD和MSSA的三维地震数据随机噪音压制方法，其特征在于，包括以下步骤：

S1，计算地震数据的傅里叶频谱；

S2，计算地震数据的尺度空间表示结果SSR并确定模态分解的个数；

S3，确定模态分解参数及模态分解；

S4，采用多通道奇异谱分析MSSA方法处理各模态数据IMF并得到各模态数据IMF的去噪结果，转换模态数据IMF数据体到频率域，为每个地震剖面分别构造汉克尔矩阵，使用奇异值分解和秩缩方法保留前a个奇异值并进行矩阵重构，重构后的矩阵从频域转换到时间域即可得到去噪后的模态数据IMF；

S5，对去噪后IMF的求和，完成噪音压制。

2.根据权利要求1所述的一种基于MVMD和MSSA的三维地震数据随机噪音压制方法，其特征在于，S1中，首先读取三维地震数据，选取地震数据的一个二维剖面，并计算此剖面的傅里叶频谱。

3.根据权利要求1所述的一种基于MVMD和MSSA的三维地震数据随机噪音压制方法，其特征在于，S2中，通过尺度空间表示结果SSR确定待分解地震数据的模态个数，对于给定的地震数据的傅里叶频谱，计算得到全局变化趋势，通过计算全局变化趋势的峰值个数得到模态分解的个数。

4.根据权利要求1所述的一种基于MVMD和MSSA的三维地震数据随机噪音压制方法，其特征在于，S3中，对于待分析的地震道数量为C的地震数据，计算输入地震数据的第c道地震信号，以及输入地震数据的第m个IMF的第c道信号，使用拉格朗日乘子法和交替方向乘数法求解约束优化问题得到模态分解后的IMF。

5.根据权利要求1所述的一种基于MVMD和MSSA的三维地震数据随机噪音压制方法，其特征在于，S5中，为每个使用MSSA去噪处理后的模态数据IMF求和，求和后的结果即为地震数据体压制随机噪声后的滤波结果。