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

本发明专利技术公开了一种压制地震资料涌浪噪声的方法、存储介质及设备，读取原始地震资料，并对其进行高通滤波；在不损失地震资料有效信号的情况下滤除一部分低频涌浪噪声；针对原始地震资料中有效信号的特点选取固定稀疏变换基；针对滤波得到的高频成分，基于半二次优化技术和线性化Bregman方法构造交替迭代优化方法在压制涌浪噪声的同时得到地震资料有效信号在稀疏变换基下的稀疏表示；直到原始地震资料中涌浪噪声压制完成。本发明专利技术的压制地震资料涌浪噪声的相关熵诱导度量鲁棒稀疏表示方法，可以有效压制高振幅涌浪噪声，提高地震资料信噪比。

【技术实现步骤摘要】
一种压制地震资料涌浪噪声的方法、存储介质及设备
本专利技术属于地震勘探数据处理
，具体涉及一种压制地震资料涌浪噪声的相关熵诱导度量鲁棒稀疏表示方法、存储介质及设备。
技术介绍
按照噪声在地震剖面上出现的特征，将其分为随机噪声和相关噪声。地震资料成像过程中，检波器接收到的反射信号含有大量的随机噪声。随机噪声的存在，直接影响了地震资料时间剖面同相轴的连续性，极大增加了地震资料解释的难度。为了满足目前地震资源勘探的要求，压制随机噪声，提高地震资料信噪比是地震数据处理的首要任务。涌浪噪声是海洋地震资料中常见的噪声类型，它是海况条件时进行地震资料采集时所导致的，特别是拖缆沉放深度较浅时更严重。涌浪噪声的特点是振幅强度较高，频率较低，在单炮记录上显示为干扰波振幅几乎不随时间而衰减。它通常影响相邻连续几道，在原始炮集记录上的形态呈条带状，或一串斑点，对地震资料的中深层形成强烈的低频噪声背景。低频滤波可以消除一定量的涌浪噪声，但不能完全消除。现有技术：中值滤波方法。该方法基于噪声特性设定一个窗口长度控制阈值,再根据地震数据本文档来自技高网...

【技术保护点】
1.一种压制地震资料涌浪噪声的方法，其特征在于，包括以下步骤：/nS1、读取原始地震资料，并对其进行高通滤波；在不损失地震资料有效信号的情况下滤除一部分低频涌浪噪声；/nS2、针对原始地震资料中有效信号的特点选取固定稀疏变换基；/nS3、针对滤波得到的高频成分，基于半二次优化技术和线性化Bregman方法构造交替迭代优化方法在压制涌浪噪声的同时得到地震资料有效信号在稀疏变换基下的稀疏表示；/nS4、重复步骤S3直到原始地震资料中涌浪噪声压制完成。/n

【技术特征摘要】
1.一种压制地震资料涌浪噪声的方法，其特征在于，包括以下步骤：
S1、读取原始地震资料，并对其进行高通滤波；在不损失地震资料有效信号的情况下滤除一部分低频涌浪噪声；
S2、针对原始地震资料中有效信号的特点选取固定稀疏变换基；
S3、针对滤波得到的高频成分，基于半二次优化技术和线性化Bregman方法构造交替迭代优化方法在压制涌浪噪声的同时得到地震资料有效信号在稀疏变换基下的稀疏表示；
S4、重复步骤S3直到原始地震资料中涌浪噪声压制完成。


2.根据权利要求1所述的压制地震资料涌浪噪声的相关熵诱导度量鲁棒稀疏表示方法，其特征在于，步骤S1中，将原始地震资料中二维地震剖面y表示为有效信号s与涌浪噪声e的线性叠加，采用高通滤波，在不损失地震资料有效信号的情况下滤除一部分低频涌浪噪声


3.根据权利要求2所述的压制地震资料涌浪噪声的相关熵诱导度量鲁棒稀疏表示方法，其特征在于，低频涌浪噪声具体为：



其中，为压制一部分低频涌浪噪声以后的地震资料。


4.根据权利要求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表示位移。


5.根据权利要求4所述的压制地震资料涌浪噪声的相关熵诱导度量鲁棒稀疏表示方法，其特征在于，Curvelet反变换的实现步骤如下：
S205、将每个尺度和角度方向的CD(j,l,k)序列做二维FFT变换得到如下傅里叶...

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