一种地震记录的方波压制方法

技术领域

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

背景技术

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

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

发明内容

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

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

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

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

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

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

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

设地震记录为s

对于第i个样点，在C

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

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

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

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

将s

若M＞N，令IM＝M-L，则在序列A

由s′

按照本发明的另一方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：

利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；本发明克服了常规方法的边界效应及对有效信号的破坏，保真度较高。

附图说明

图1是本发明实施例提供的一种地震记录的方波压制方法的流程示意图；

图2是本发明实施例提供的一种压制方波前的地震数据图；

图3是本发明实施例提供的一种压制方波后的地震数据图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

如图1所示，本发明实施例提供的一种地震记录方波压制方法，主要包括以下步骤：

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

在本发明实施例中，步骤S1可以通过以下方式实现：

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

具体地，设地震记录为s

设置滑动时窗长度参数为2*L+1，识别方波的门槛因子为C，令：

当C

S2：如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值不变，接着处理下一个样点，以此类推。

在本发明实施例中，步骤S2可以通过以下方式实现：

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

具体地，取步骤S1中的同一个时窗内样点数据，令：

将s

若M＞N，令IM＝M-L，则在序列A

s′i＝s

s′

图2是压制方波前的地震数据图。由于仪器漏码引起的方波能量非常强，几乎将有效信号完全淹没，如果不做去噪处理将严重影响后续子波处理、速度分析以及叠加效果。

图3是压制方波后的地震数据图。与图2相比，方波被完全去除，掩盖在方波下的地震信息也显现出来。去噪效果非常明显。

综上，本发明的一种地震记录方波压制方法主要改进包括三点：(1)利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；(2)在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；(3)该技术克服了常规方法的边界效应及对有效信号的破坏，保真度较高。

本申请还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，程序被处理器执行时实现方法实施例中的地震记录的方波压制方法。

需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。