一种用于抑制杂波的快速滤波方法及存储介质

摘要

本发明公开了用于抑制杂波的快速滤波的方法及存储介质。用于抑制杂波的快速滤波方法，包括：确定地震数据滤波方向，可以对地震数据的数据体进行多个数据剖面的划分，并对每个数据剖面进行滤波；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制地震数据的杂波。本发明实施例的用于抑制杂波的快速滤波的方法及存储介质，可以有效去除地震数据杂波，进而为矿产、常规油气、页岩油气及煤层气的勘探和开发服务。

说明书

技术领域

本发明属于矿产、油气及煤层气勘探与开发中的地震数据处理技术领域，特别涉及一种用于抑制杂波的快速滤波方法及存储介质。

背景技术

勘探是为了寻找和查明地下资源，而利用各种勘探手段了解地下的地质状况，确定地下资源聚集的有利地区。为了进行勘探与开发，可以采用地震勘探的方法。地震勘探是根据地质学和物理学的原理，采用人工方法引起地壳振动，再用精密仪器记录下爆炸后地面上各点的震动情况，把记录下来的地震数据经过处理、解释，推断地下地质构造的特点，寻找可能的地下资源构造。

在处理地震数据时，会存在大量的杂波，这会干扰地震数据处理的效果。尤其是，针对中国石化自主研发的地震综合解释软件-NEWS地震综合解释系统来说，需要开发一种快速滤波技术为地震数据信号提频服务。

发明内容

本发明所要解决的技术问题之一是如何提供一种用于抑制杂波的快速滤波方法及存储介质。

为了解决上述技术问题，本申请的实施例首先提供了一种用于抑制杂波的快速滤波方法，包括：确定地震数据滤波方向；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制所述地震数据的杂波。

优选地，所述确定地震数据滤波方向，包括：基于地震数据的数据体确定数据剖面；根据所述数据剖面确定滤波方向。

优选地，所述数据剖面包括线方向数据剖面、道方向数据剖面和切片方向数据剖面。

优选地，所述针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数，包括：针对所述地震数据滤波方向，确定滤波本位点；以所述滤波本位点为中心，定义n×n的网格，所述网格为所述时窗范围，其中，n为奇数；定义所述时窗滑动步数为大于等于1。

优选地，n为3、5、7或9。

优选地，所述时窗滑动步数等于1，则不做标识；所述时窗滑动步数大于1，则标识结果数据为稀疏重采样。

优选地，对所述时窗范围内的地震数据进行滤波，包括：利用算数平均对所述时窗范围内的地震数据进行滤波。

优选地，对所述时窗范围内的地震数据进行滤波，包括：利用均方根平均对所述时窗范围内的地震数据进行滤波。

优选地，对所述时窗范围内的地震数据进行滤波，包括：利用加权平均对所述时窗范围内的地震数据进行滤波。

本发明实施例还公开了一种存储介质，其存储有计算机程序，所述计算机程序被处理器执行时实现上述任意一项所述方法。

与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：

用于抑制杂波的快速滤波方法，包括：确定地震数据滤波方向，可以对地震数据的数据体进行多个数据剖面的划分，并对每个数据剖面进行滤波；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制所述地震数据的杂波。通过该方法可以有效去除地震数据杂波，进而为矿产、常规油气、页岩油气及煤层气的勘探和开发服务。

本发明的其他优点、目标，和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书，权利要求书，以及附图中所特别指出的结构来实现和获得。

附图说明

附图用来提供对本申请的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本申请实施例的附图与本申请的实施例一起用于解释本申请的技术方案，但并不构成对本申请技术方案的限制。

图1是根据本发明实施例一的用于抑制杂波的快速滤波的方法的流程示意图。

具体实施方式

以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本发明的保护范围之内。

用于抑制杂波的快速滤波方法，包括：确定地震数据滤波方向，可以对地震数据的数据体进行多个数据剖面的划分，并对每个数据剖面进行滤波；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制地震数据的杂波。通过该方法可以有效去除地震数据杂波，进而为矿产、常规油气、页岩油气及煤层气的勘探和开发服务。

下面以具体的实施例对本发明作进一步说明。

如图1所示，为本发明实施例的一种用于抑制杂波的快速滤波方法，其包括步骤S1至S3。

S1，确定地震数据滤波方向。

确定地震数据滤波方向，包括：基于地震数据的数据体确定数据剖面；根据数据剖面确定滤波方向。其中，数据剖面包括线方向数据剖面、道方向数据剖面和切片方向数据剖面。

当对地震数据进行快速滤波时，可以将地震数据根据数据体划分为线方向、道方向或切片方向的数据剖面，其中，线方向、道方向和切片方向为数据体的三个相互正交的方向，并分别对各数据剖面进行滤波处理，从各数据剖面方向对地震数据进行滤波来抑制地震数据的杂波。当处理完全部数据剖面的数据后，就完成了整个数据体的处理。

S2，针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数。

在一个实施例中，针对地震数据滤波方向，定义时窗范围和时窗滑动步数，包括：针对地震数据滤波方向，确定滤波本位点；以滤波本位点为中心，定义n×n的网格，网格为时窗范围，其中，n为奇数；定义时窗滑动步数为大于等于1。

具体地，针对地震数据滤波方向，确定滤波本位点，以滤波本位点定义3×3或者5×5或者7×7或者9×9的网格，首先在这个时窗范围内进行滤波处理，即基于网格的能量进行处理，并将滤波处理的结果填到滤波本位点，然后根据时窗滑动步数滑动时窗范围，并在这个滑动后的时窗范围内进行滤波处理并将滤波处理的结果填到滤波本位点，持续这个过程直到这个滤波方向的地震数据全部处理完毕。处理完此滤波方向后可以继续处理其他滤波方向的地震数据，直到完成整个数据体的处理。

优选地，为了快速滤波处理，n小于等于11。例如，n取1、3、5、7、9、11等奇数，分别形成1×1、3×3、5×5、7×7、9×9、11×11的网格。利用网格作为时窗范围，将时窗范围内的能量进行处理，并将处理结果填写到滤波本位点。

其中，时窗滑动步数等于1，则不做标识，即处理结果为全采样，不需要进行标识；时窗滑动步数大于1，则标识结果数据为稀疏重采样，即处理结果为稀疏采样，对处理结果进行标识。

S3，对时窗范围内的地震数据进行滤波以抑制地震数据的杂波。

根据上述方法对时窗范围内的地震数据进行滤波，可以去除地震数据的杂波，滤波处理后的地震数据可以进行用于后续地震处理，用于分析矿产、常规油气、页岩油气及煤层气的情况。

在一个实施例中，对时窗范围内的地震数据进行滤波，包括：

利用算数平均(如下式所示)对时窗范围内的地震数据进行滤波。

在一个实施例中，对时窗范围内的地震数据进行滤波，包括：

利用均方根平均(如下式所示)对时窗范围内的地震数据进行滤波。

在一个实施例中，对时窗范围内的地震数据进行滤波，包括：

利用加权平均(如下式所示)对时窗范围内的地震数据进行滤波。

上述三个对时窗范围内的地震数据进行滤波实施例的方案，可以任选其中一种对网格中的能量进行处理并回填到滤波本位点。进而，时窗范围滑动从而对全部地震数据进行滤波，滤波处理后的地震数据中杂波明显减少。利用本发明申请的实施例，可以有效去除地震数据的杂波，滤波处理后的地震数据可以为矿产、常规油气、页岩油气及煤层气的勘探、开发服务。

本发明实施例公开了一种存储介质，其存储有计算机程序，计算机程序被处理器执行时实现上述任一项所述方法。其中，计算机程序包括：确定地震数据滤波方向；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制所述地震数据的杂波的步骤。

在一个实施例中，所述确定地震数据滤波方向，包括：

基于地震数据的数据体确定数据剖面；

根据所述数据剖面确定滤波方向。

在一个实施例中，所述数据剖面包括线方向数据剖面、道方向数据剖面和切片方向数据剖面。

在一个实施例中，所述针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数，包括：

针对所述地震数据滤波方向，确定滤波本位点；

以所述滤波本位点为中心，定义n×n的网格，所述网格为所述时窗范围，其中，n为奇数；

定义所述时窗滑动步数为大于等于1。

在一个实施例中，n为3、5、7或9。

在一个实施例中，所述时窗滑动步数等于1，则不做标识；

所述时窗滑动步数大于1，则标识结果数据为稀疏重采样。

在一个实施例中，对所述时窗范围内的地震数据进行滤波，包括：

利用算数平均对所述时窗范围内的地震数据进行滤波。

在一个实施例中，对所述时窗范围内的地震数据进行滤波，包括：

利用均方根平均对所述时窗范围内的地震数据进行滤波。

在一个实施例中，对所述时窗范围内的地震数据进行滤波，包括：

利用加权平均对所述时窗范围内的地震数据进行滤波。

通过确定地震数据滤波方向，可以对地震数据的数据体进行多个数据剖面的划分，并对每个数据剖面进行滤波；针对所述地震数据滤波方向，定义时窗范围和时窗滑动步数；对所述时窗范围内的地震数据进行滤波以抑制地震数据的杂波，可以有效去除地震数据杂波，进而为矿产、常规油气、页岩油气及煤层气的勘探和开发服务。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉该技术的人员在本发明所揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。