[发明专利]一种高效的地震数据随机噪声衰减方法

说明书

技术领域

本发明涉及地球物理勘探地震数据处理技术范畴，是一种高效的地震数据随机噪声衰减方法。

背景技术

地震数据在采集过程中，不可避免的受到各种噪声的干扰，致使有效信号失真，提高地震数据的信噪比是地震勘探中一项重要而基础的工作。频率-空间域预测随机噪声衰减技术是压制三维叠前地震数据随机噪声的一种行之有效的处理方法，对于提高地震数据的信噪比有着显著的效果。但是，现有的方法计算量巨大，处理效率低下，处理周期很长，严重限制了该技术在石油物探实际生产中的推广应用。

现有频率-空间域预测随机噪声衰减技术的一大弊端是对地震资料各方向做一致处理。众所周知，三维叠前地震资料是各向异性的，空间各方向反射波复杂程度不同，各方向间可能存在巨大差异。现有技术忽视了三维叠前地震资料的这一重要特征，没有对不同方向区分处理，在某些空间方向计算量浪费严重，导致这一随机噪声衰减技术难以推广。与此同时，由于各方向一致处理的限制，现有技术难以处理空间分布不均的地震数据，例如宽线数据等。

发明内容

本发明目地在于提供一种提高地震数据信噪比、提升处理效率、节约处理成本的一种高效的地震数据随机噪声衰减方法。

本发明采用以下步骤实现：

1）采集三维叠前地震数据；

2）根据地震数据空间各方向反射波复杂程度的不同，分别给定各方向的预测算子长度，根据预测算子长度给定各空间方向计算窗长度，并给定时间窗大小；

所述的预测算子长度如果某方向反射波较复杂,则可适当增加这个方向预测算子长度;如果某方向反射波较平缓,则可以适当减少该方向预测算子长度。

所述的预测算子长度为5－8。

所述的给定各空间方向计算窗长度选取各空间预测算子长度的四倍加一。

时间窗为500毫秒一个计算时间窗。

3）根据各方向计算窗大小将数据分成小计算窗；

4）在当前计算窗内沿时间方向做Fourier变换得到频率空间域原始数据

{s_i,j,k;f}，其中i,j,k分别表示X、Y、O三个空间方向，f表示频率；

5）用下式计算三维预测算子序列{x_l,m,n}：

Σl,m,n(Σi,j,ksi-l,j-m,k-n;f·s‾i-p,j-q,k-r;f)xl,m,n=Σi,j,ksi,j,k;f·s‾i-p,j-q,k-r;f]]>