[发明专利]一种地震记录的方波压制方法

说明书

本发明公开了一种地震记录的方波压制方法，属于石油勘探开发领域，该方法包括以下步骤：在单炮地震数据体上逐道、逐样点判断方波是否存在；如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值保持不变，接着处理下一个样点。本发明利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；该技术克服了常规方法的边界效应及对有效信号的破坏，保真度较高。

技术领域

本发明属于石油勘探开发领域，更具体地，涉及一种地震记录中方波的压制方法。

背景技术

在地震记录中，方波通常是由于仪器或检波器漏码故障引起，在地震记录上表现为一个非常大的常数干扰值叠加在有效信号上，形同废道废炮。不同道的干扰值相差较大，可能只是影响局部数据，也可能对整道数据都有影响。

常规处理方法往往有去道去炮，内切除，高通滤波，中值滤波等等。去道去炮和内切除都造成了有效信号的丧失和经济上的损失。而高通滤波因为大值的影响会损害有效信号，且会造成严重的边界效应。中值滤波往往耗时较大，当方波连续多道出现时，去噪结果往往存在较大误差，且在没有方波的地方对有效信号有较大影响。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提出了一种地震记录的方波压制方法，用于解决去噪对有效信号损害较大的问题。

为实现上述目的，按照本发明的一个方面，提供了一种地震记录的方波压制方法，包括：

(1)在单炮地震数据体上逐道、逐样点进行处理，以各样点为中心，通过预设长度时窗内的数据进行方波存在与否判断；

(2)如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值保持不变。

在一些可选的实施方案中，所述通过预设长度时窗内的数据进行方波存在与否判断，包括：

通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断。

在一些可选的实施方案中，所述通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断，包括：

设地震记录为s_i，其中，下标i为样点序号，设置滑动时窗长度参数为2*L+1，识别方波的门槛因子为C，令：L为自然数；

对于第i个样点，在C_i＞C时，则第i个样点存在方波，否则第i个样点不存在方波。

在一些可选的实施方案中，所述进行方波预测与压制，包括：

通过均值结合中值技术进行方波预测与压制。

在一些可选的实施方案中，所述通过均值结合中值技术进行方波预测与压制，包括：

取步骤(1)中的同一个时窗内样点数据，令

将s_j，(j＝i-L，...，i+L)分为两个序列A_i，(i＝1，...，M)和B_i，(i＝1，...，N)，M+N＝2*L+1，其中，A_i为序列s_j中大于的样点值，B_i为序列s_j中小于的样点值，M，N为自然数；

若M＞N，令IM＝M-L，则在序列A_i中找到第IM小的值为样点i处方波值，记为V_i；若M≤N，令IN＝N-L，在序列B_i中找到第IN大的值为样点i处方波值，记为V_i；