[发明专利]一种地震数据随机噪声压制方法及装置

权利要求书

1.一种地震数据随机噪声压制方法，其特征在于，包括以下步骤：

对获取的地震数据进行互补经验态分解CEEMD，获得本征模函数IMF分量序列；

确定所述IMF分量序列中含有随机噪声的IMF分量和不含有随机噪声的IMF分量；

对所述含有随机噪声的IMF分量进行曲波最优迭代阈值去噪处理，获得处理后的IMF 分量；

利用所述处理后的IMF分量以及所述不含有随机噪声的IMF分量进行重构获得去除 随机噪声的地震数据。

2.如权利要求1所述的地震数据随机噪声压制方法，其特征在于，所述对含有随机 噪声的IMF分量进行曲波最优迭代阈值去噪处理，获得处理后的IMF分量，包括：

根据所述含有随机噪声的IMF分量中所含随机噪声的不同，对所述含有随机噪声的 IMF分量选取不同的阈值进行曲波最优迭代阈值去噪处理，获得处理后的IMF分量。

3.如权利要求1所述的地震数据随机噪声压制方法，其特征在于，所述确定IMF分 量序列中含有噪声的IMF分量和不含有噪声的IMF分量，包括以下步骤：

计算所述地震数据和所述IMF分量序列之间的互相关系数，获得互相关系数序列；

确定所述互相关系数序列中的极小值；

根据所述极小值对应的IMF分量位置确定含有随机噪声的IMF分量和不含有随机噪 声的IMF分量。

4.如权利要求3所述的地震数据随机噪声压制方法，其特征在于，所述计算地震数 据和所述IMF分量之间的互相关系数，具体公式为：

R(j)=Rdata,IMFj(0)Rdata(0)RIMFj(0)]]>

其中，R(j)为互相关系数，R_data(0)为地震数据自相关后零时刻所对应的最大峰值点， 为第j个IMF分离的自相关后零时刻所对应的最大峰值点，为原始数 据与第j个IMF分量互相关后零时刻所对应的值。