[发明专利]一种高保真声波干扰压制方法

说明书

技术领域

本发明属于地震勘探技术领域，更具体地讲，涉及一种用于石油地震资料的高保真声波干扰压制方法。

背景技术

在山地地区，由于激发和接收条件的限制，大量地采用浅井组合及坑炮激发，因而导致声波较为发育。声波会在地震记录数据上形成尖锐的强初至，呈窄带出现却比较稳定。声波的特点是能量强、集中在中高频率段且传播速度稳定，因此声波的分布范围极易划定。如果这种强能量的噪声不能得到有效地压制，将直接影响后续的反褶积、速度分析等处理，对提高地震资料的信噪比、分辨率以及可信度是极为不利的。

目前用于压制声波干扰的手段仍然仅局限于一维滤波、二维滤波、人工编辑等，例如：(1)反褶积方法，其主要通过反褶积技术对声波进行压制，但在声波主频较高时，此压制方法往往不是很理想；(2)切除法：使用内切除法将声波完全剔除，其虽能从根本上消除声波对地震数据的影响并可更好第提高信噪比，但湮没在强噪声干扰中的有效信号也会损失掉，同时还会大大增加工作量，延长处理周期。

以上压制声波干扰的方法尽管有一些效果，但也极容易伤害和损失有效反射信息。因此，寻找一种有效的衰减声波方法，无论在提高处理质量还是在缩短处理周期，减轻劳动强度方面都是非常具有实际意义的。

发明内容

针对现有技术中存在的不足，本发明的目的在于解决上述技术问题中的一个或多个。

本发明基于声波具有非常稳定的传播速度以及能量强的特点，提供了一种数据量小、计算效率高且压制声波干扰效果明显的高保真声波干扰压制方法。

为了实现上述目的，本发明提供了一种高保真声波干扰压制方法，其特征在于，所述方法包括以下步骤：A)在共炮点道集的每道地震数据中划分出含有声波干扰的时窗；B)将每道地震数据中含有声波干扰的时窗内的数据进行小波分解并等分为多个频带；C)统计每道地震数据的每个频带的所有样点的振幅值并求取每道地震数据的每个频带的振幅加权中值；D)将每道地震数据的每个频带的每个样点的噪声分贝值分别与预设阈值进行比较和处理得到每道地震数据的多个频带的衰减后数据，其中，利用所述样点的振幅值与所述样点所在频带的振幅加权中值的比值计算所述样点的噪声分贝值，当样点的噪声分贝值大于或等于预设阈值时则取所述样点之前的3～5个样点的振幅的平均值代替所述样点的振幅值，当样点的噪声分贝值小于预设阈值时则保留所述样点的振幅值；E)将所述每道地震数据的多个频带的衰减后数据进行小波重构得到声波干扰压制后的地震数据。

根据本发明的高保真声波干扰压制方法的一个实施例，在步骤A中，包括以下子步骤：A1)在所述共炮点道集中拾取声波干扰起始点的起始时间t₀以及起始偏移距r₀，再利用式1计算第i道地震数据的声波干扰时间t_i：

t_i＝t₀+abs(r_i-r₀)/v——式1，式中，v为声波的传播速度；

A2)给定时移量d，划分第i道地震数据的含有声波干扰的时窗为(t_i-d，t_i+d)，其中，i为地震数据的道号且1≤i≤N，N为所述共炮点道集的总道数；t_i+d≤t_max，t_max为所述共炮点道集的最大记录时间。