[发明专利]一种基于时空域降频模型的二维地震勘探噪声去除方法

说明书

技术领域

本发明涉及的是实际地震勘探环境下数据采集所得二维地震勘探记录的随机噪声消减方法，实际地震勘探记录中反射信号的高频成分极易在传统消减方法中受到损失。本发明建立时空域的二维地震勘探记录降频模型，计算出满足给定降频条件的时空变换方向，并对时空域提取出的降频信号进行非线性滤波处理，能够有效去除地震勘探过程中的随机噪声，同时能够更好的保护反射同相轴信息，提高实际二维地震勘探记录的信噪比。为后续反射振幅、速度及频率信息的准确提取提供有力保障，进而有利于准确估算油气储量及分布范围。

背景技术

实际地震勘探工作中，不论是风吹草动、检波器随机振动、，还是叠前、叠后的处理都不可避免的会在地震勘探记录中留下随机噪声。随机噪声的存在影响对反射信号的辨认，进而影响一系列后续的地质解释工作。由于地震勘探随机噪声是由很多不可预知的因素综合作用而成，因而在地震勘探记录中随机出现，频带较宽、速度不定且没有规则的传播方向。

当前对实际地震勘探记录进行随机噪声消减时，有很多种处理方法。然而，大多数去噪方法在消减随机噪声时可以保护较低频反射信号的能量，却无法保护较高频反射能量，致使处理后信噪比提高有限，处理结果达不到地震勘探的相关要求。

时频峰值滤波技术是当前较常用的一种能够消减实际地震勘探记录随机噪声的方法。它的原理如下：实际二维地震勘探记录中沿时间方向的一道信号包含反射波和随机噪声，表示为下式：

s(t)＝x(t)+n(t)

其中x(t)表示反射波有效信号，n(t)表示随机噪声。随后进入编码过程，将含噪地震勘探记录调制为某解析信号的瞬时频率，得到解析信号：

可以证明，一定条件下，解析信号z(t)的伪Wigner-Ville分布(PWVD)的峰值是有效信号x(t)的无偏估计。因此可将x(t)从含噪信号s(t)中提取出来，即求得达到信号去噪的目的，公式如下：

使用该方法对地震勘探记录进行去噪的效果较好，但是PWVD的参数需要依据反射信号的主频率而做出调节，参数较大时能获得更好的随机噪声消减效果，但对较高主频反射能量会造成严重的损失。因而该方法只适用于反射波主频率25Hz以下的地震勘探记录去噪处理。而当反射波主频较大时，无偏估计条件得不到满足，反射信息与随机噪声一并被去除，影响有效信号同相轴的识别，进而影响地震勘探记录信噪比的进一步提高。

发明内容

发明提供一种基于时空域降频模型的二维地震勘探噪声去除方法，用于解决的技术问题是在地震勘探记录采集后，利用二维地震勘探记录中邻近地震道与道之间的子波相关特性，不在原始的时间-偏移矩域而是变换至时间-空间域，通过构建时空域降频模型寻找到能将实际地震勘探记录高频信号转化为低频信号的时空方向。在计算出的时空方向上提取时空降频信号并设置合理参数进行随机噪声消减处理，并且与当前地震勘探去噪技术结合使用，达到提高实际地震勘探记录信噪比的目的。为后续反射振幅、速度及频率信息的准确提取提供有力保障，进而有利于准确估算油气储量及分布范围。

本发明采取的技术方案是，包括下列步骤：

1)二维地震勘探记录获取：

采用与地质构造走向垂直的方向来安排测线，沿地震测线等间距布置N个检波器组，选择中间放炮排列方式来接收地震波信号，根据观测仪器的不同，检波器组的数量N可少至几十个，多至上千个，激发炮点后，每个检波器接收的信号通过放大器和记录器，得到一道地震勘探波形记录，记录的长度与接收时间成正比，与检波器的采样频率成反比，N道记录组成一幅二维地震勘探记录；

2)定义时空域降频模型：

时空域降频模型的建立来源于径向道变换思想，对一幅二维地震勘探记录进行径向道变换可使某些变换域信号的视时间拉长，周期增加，频率降低，即信号的线性度获得降低，定义时空域降频模型是为了获取低线性度地震信号的时空方向，一幅二维地震勘探记录共N道，设局部范围内各道的地震子波时距曲线近似满足直线方程，该直线l₀与时间方向夹角为单个地震子波的周期为T，相邻道地震子波的时间差为t₀，第一道子波初始位置与第N道子波结束位置的连线记为l_R，与时间方向夹角为l₀与l_R之间的夹角为默认二维地震勘探记录处理时的时间间隔与空间间隔一致，即时域采样间隔及相邻道的间隔均为△t，模型分为和两种情况；

3)二维地震勘探记录降频模型参数选取方法：