[发明专利]一种压制地震噪声的方法、装置及系统

说明书

本申请实施例公开了一种压制地震噪声的方法、装置及系统，所述方法包括将原始地震数据进行曲波变换处理，获得曲波系数；对所述曲波系数求取梯度获得局部波向量；对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数；对所述阈值处理后的曲波系数进行逆曲波变换，获得压制噪声后的地震数据。利用本申请各个实施例，可以更加有效的压制地震噪声，提高地震数据的信噪比。

技术领域

本发明涉及地震资料处理技术领域，特别地，涉及一种压制地震噪声的方法、装置及系统。

背景技术

地震资料数据处理作为地震勘探的关键环节之一，在其“高信噪比、高分辨率和高保真度”的三大要求中，信噪比是基础。由于种种原因，有效信号被各种各样的噪声污染，造成地震资料信噪比低，给速度分析、偏移成像带来极大困难,严重影响成像的真实性和可靠性，难以达到勘探需求。因此压制噪声、提高资料的信噪比成为地震数据处理的首要任务。

曲波变换因其能够很好地处理高维信号的线奇异性，而被引入到地震资料数据处理中。现有技术中利用曲波变换进行去噪处理时，一般都是基于噪声与有效信号曲波系数数值大小不同，通过设定阈值来压制噪声。但是地震数据中有些噪声(如相干噪声)与有效反射波的曲波系数很难直接通过数值进行有效区分，从而影响最终获得的地震数据的信噪比。

发明内容

本申请实施例的目的在于提供一种压制地震噪声的方法、装置及系统，可以更加有效的压制地震噪声，提高地震数据的信噪比。

本申请提供的一种压制地震噪声的方法、装置及系统是通过包括以下方式实现的：

一种压制地震噪声的方法，包括：

将原始地震数据进行曲波变换处理，获得曲波系数；

对所述曲波系数求取梯度获得局部波向量；

对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数；

对所述阈值处理后的曲波系数进行逆曲波变换，获得压制噪声后的地震数据。

本申请实施例的压制地震噪声的方法，所述方法还包括：

获取绝对值大于等于第一预设阈值的曲波系数；

相应的，所述对所述曲波系数求取梯度获得局部波向量包括：对绝对值大于等于第一预设阈值的曲波系数求取梯度，获得局部波向量。

本申请实施例的压制地震噪声的方法，所述对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数，包括：

对所述局部波向量的大小进行第一阈值处理，获得第一阈值处理后的局部波向量；

根据预先存储的局部波向量与曲波系数的索引关系，获得第一阈值处理后的局部波向量对应的曲波系数作为所述阈值处理后的曲波系数。

本申请实施例的压制地震噪声的方法，所述对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数，包括：

根据对地震噪声与有效地震波的局部波向量数据特征差异性的分析结果，对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数。

本申请实施例的压制地震噪声的方法，所述对所述曲波系数求取梯度获得局部波向量，包括：

根据如下公式获得不同位置处的局部波向量：

其中，W_d(a,θ,b)表示地震数据d的曲波系数，a表示曲波系数的尺度参数，θ表示曲波系数的方向参数，b表示曲波系数的位置参数，表示在位置b上对曲波系数做梯度，i表示虚数单位，v_d(a,θ,b)表示位置b上的局部波向量。