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

摘要

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

主权项

一种高保真声波干扰压制方法，其特征在于，所述方法包括以下步骤：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为所述共炮点道集的最大记录时间。