一种压制地震资料异常噪声的方法、存储介质及设备技术

本发明专利技术公开了一种压制地震资料异常噪声的方法、存储介质及设备，读取含有高振幅异常噪声的原始地震资料，根据地震资料中噪声的复杂程度选取混合高斯分布的高斯成分个数和迭代次数；针对原始地震资料中有效信号的特点选取固定稀疏变换基，利用线性化Bregman方法得到地震资料中有效信号的稀疏表示，根据稀疏表示与固定稀疏变换基的乘积和原始地震资料得到关于有效信号的估计和异常噪声的估计，将异常噪声建模为混合高斯分布，利用EM算法更新每个采样点上异常噪声的权重矩阵；直到原始地震资料中异常噪声压制完成或者达到迭代次数。本发明专利技术的压制地震资料异常噪声的混合高斯鲁棒稀疏表示方法，可以有效压制异常噪声，提高地震资料信噪比。

【技术实现步骤摘要】
一种压制地震资料异常噪声的方法、存储介质及设备
本专利技术属于地震勘探数据处理
，具体涉及一种压制地震资料异常噪声的混合高斯鲁棒稀疏表示方法、存储介质及设备。
技术介绍
按照噪声在地震剖面上出现的特征，将其分为随机噪声和相关噪声。地震资料成像过程中，检波器接收到的反射信号含有大量的随机噪声。随机噪声的存在，直接影响了地震资料时间剖面同相轴的连续性，极大增加了地震资料解释的难度。为了满足目前地震资源勘探的要求，压制随机噪声，提高地震资料信噪比是地震数据处理的首要任务。异常噪声属于随机噪声。陆上和海上地震资料中通常包含异常噪声。气流、电力线和其他人文噪声、记录和奇偶性误差、未修正的极性反转、孤立噪声或耦合不良的检波器、风、雨和其他影响都可能导致空间方向上的异常噪声。异常噪声比一般随机噪声幅值大，通常为非高斯分布，从而对高斯分布假设下的经典去噪方法产生了无限制的失真影响。现有技术：中值滤波方法。该方法基于噪声特性设定一个窗口长度控制阈值,再根据地震数据与阈值的关系选取合适的窗口进行中值滤波,以取得有效信号与噪声之本文档来自技高网...

【技术保护点】
1.一种压制地震资料异常噪声的方法，其特征在于，包括以下步骤：/nS1、读取含有高振幅异常噪声的原始地震资料，根据地震资料中噪声的复杂程度选取混合高斯分布的高斯成分个数和迭代次数；/nS2、针对原始地震资料中有效信号的特点选取固定稀疏变换基；/nS3、利用线性化Bregman方法得到地震资料中有效信号的稀疏表示，根据稀疏表示x

【技术特征摘要】
1.一种压制地震资料异常噪声的方法，其特征在于，包括以下步骤：
S1、读取含有高振幅异常噪声的原始地震资料，根据地震资料中噪声的复杂程度选取混合高斯分布的高斯成分个数和迭代次数；
S2、针对原始地震资料中有效信号的特点选取固定稀疏变换基；
S3、利用线性化Bregman方法得到地震资料中有效信号的稀疏表示，根据稀疏表示x*与固定稀疏变换基的乘积和原始地震资料得到关于有效信号的估计和异常噪声估计，将异常噪声建模为混合高斯分布，利用EM算法更新每个采样点上异常噪声的权重矩阵w；
S4、重复步骤S3直到原始地震资料中异常噪声压制完成，。


2.根据权利要求1所述的方法，其特征在于，步骤S1中，将原始地震资料中二维地震剖面y表示为有效信号s与异常噪声e的线性叠加；利用混合高斯分布拟合异常噪声e，根据异常噪声e的复杂程度确定迭代次数。


3.根据权利要求1所述的方法，其特征在于，步骤S2中，基于wrapping方法的离散Curvelet，其正变换步骤如下：
S201、对二维信号做二维FFT，得到信号的二维傅里叶矩阵如下：



其中，(n1,n2)表示数据中的采样点；
S202、对每个尺度和角度方向的频率矩阵做加窗处理如下：



其中，为窗函数，j表示第j个尺度，l表示第l个角度；
S203、对步骤S202中得到的每个矩阵做环绕处理，标准遵循Candes准则得到如下：



其中，W表示对加窗后的频率矩阵做环绕处理；
S204、对步骤S203中得到的做二维FFT处理获得Curvelet变换系数矩阵CD(j,l,k)，其中，k表示位移。


4.根据权利要求3所述的方法，其特征在于，Curvelet反变换的实现步骤如下：
S205、将每个尺度和角度方向的CD(j,l,k)序列做二维FFT变换得到如下傅里叶序列如下：



S206、将...

【专利技术属性】
技术研发人员：陈文超，徐威威，周艳辉，王晓凯，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61