线性干扰衰减方法及系统与流程

本发明涉及地震数据处理技术领域，具体地，涉及一种线性干扰衰减方法及系统。

背景技术：

线性干扰是地震原始道集中的常见干扰，对线性干扰进行噪声衰减是地震数据处理过程中的常规步骤，线性干扰衰减需要提前获取线性干扰的频带和速度范围参数。

为了有效地衰减线性干扰，处理人员需要在工区范围内在不同的位置人工统计线性干扰的特征，统计过程效率低、误差大，受制于时间和人力限制也无法对全工区每个道集进行统计。

技术实现要素：

本发明实施例的主要目的在于提供一种线性干扰衰减方法及系统，以提高线性干扰去噪的效率和精度。

为了实现上述目的，本发明实施例提供一种线性干扰衰减方法，包括：

对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集；

对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅；

确定每个时间样点的多个振幅中的最大值对应的射线参数；

根据每个道集中的每个频带的射线参数形成第一射线参数数据体；

获取其中一个道集的线性干扰速度和该道集对应的射线参数；

根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体；

根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

本发明实施例还提供一种线性干扰衰减系统，包括：

分频单元，用于对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集；

线性拉东变换单元，用于对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅；

射线参数确定单元，用于确定每个时间样点的多个振幅中的最大值对应的射线参数；

第一射线参数数据体单元，用于根据每个道集中的每个频带的射线参数形成第一射线参数数据体；

射线参数获取单元，用于获取其中一个道集的线性干扰速度和该道集对应的射线参数；

第二射线参数数据体单元，用于根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体；

线性干扰衰减单元，用于根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现所述的线性干扰衰减方法的步骤。

本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现所述的线性干扰衰减方法的步骤。

本发明实施例的线性干扰衰减方法及系统先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例中线性干扰衰减方法的流程图；

图2是本发明实施例中异常能量干扰衰减前的不同时窗及属性平面图；

图3是本发明实施例中异常能量干扰衰减后的不同时窗及属性平面图；

图4是本发明实施例中线性干扰衰减系统的结构框图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。

鉴于现有技术的统计过程效率低、误差大，受制于时间和人力限制也无法对全工区每个道集进行统计，本发明实施例提供一种线性干扰衰减方法，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。以下结合附图对本发明进行详细说明。

图1是本发明实施例中线性干扰衰减方法的流程图。如图1所示，线性干扰衰减方法包括：

s101：对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集。

其中，道集可以为炮集、检波点道集、共中心点道集或十字排列道集，道集可以分为左右排列分别处理。频带宽度可以为10hz-20hz。

s102：对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅。

其中，线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数，或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。线性拉东变换结果的纵坐标为时间，横坐标为射线参数。

s103：确定每个时间样点的多个振幅中的最大值对应的射线参数。

s104：根据每个道集中的每个频带的射线参数形成第一射线参数数据体。

s105：获取其中一个道集的线性干扰速度和该道集对应的射线参数。

s106：根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体。

其中，在执行s106之前，还可以包括：对第一射线参数数据体沿空间方向进行平滑处理以消除部分异常结果。平滑半径可以为道集位置横向间隔的3-5倍。

此时s106包括：根据线性干扰速度、该道集对应的射线参数和经过平滑处理的第一射线参数数据体，获得第二射线参数数据体。

具体实施时，第二射线参数数据体为经过平滑处理的第一射线参数数据体与比例系数的积；比例系数为线性干扰速度与该道集对应的射线参数的商。即，线性干扰速度除以该道集对应的射线参数，乘以经过平滑处理的第一射线参数数据体，即可得到第二射线参数数据体。

s107：根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

图1所示的线性干扰衰减方法的执行主体可以为计算机。由图1所示的流程可知，本发明实施例的线性干扰衰减方法先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

本发明实施例的具体流程如下：

1、对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集。

2、对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅。线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数，或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。

3、确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体。

4、对第一射线参数数据体沿空间方向进行平滑处理以消除部分异常结果。平滑半径可以为道集位置横向间隔的3-5倍。

5、获取其中一个道集的线性干扰速度和该道集对应的射线参数。

6、令线性干扰速度除以该道集对应的射线参数，并乘以经过平滑处理的第一射线参数数据体，得到第二射线参数数据体。

7、根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

图2是本发明实施例中线性干扰衰减前的不同时窗及射线参数属性平面图；图3是本发明实施例中线性干扰衰减后的不同时窗及射线参数属性平面图。由图2可以看出在线性干扰衰减前，线性干扰的倍数分布不均匀且横向差异较大，如果通过人工调查来获得线性干扰视速度数据，不但工作量大且精度低。由图3可以看出经过线性干扰衰减后，地震数据射线参数属性横向趋于一致，说明地震数据中在线性干扰视速度范围内已经没有明显的线性信号，达到了线性干扰衰减的目的。

综上，本发明实施例的线性干扰衰减方法先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

基于同一发明构思，本发明实施例还提供了一种线性干扰衰减系统，由于该系统解决问题的原理与线性干扰衰减方法相似，因此该系统的实施可以参见方法的实施，重复之处不再赘述。

图4是本发明实施例中线性干扰衰减系统的结构框图。如图4所示，线性干扰衰减系统包括：

分频单元，用于对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集；

线性拉东变换单元，用于对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅；

射线参数确定单元，用于确定每个时间样点的多个振幅中的最大值对应的射线参数；

第一射线参数数据体单元，用于根据每个道集中的每个频带的射线参数形成第一射线参数数据体；

射线参数获取单元，用于获取其中一个道集的线性干扰速度和该道集对应的射线参数；

第二射线参数数据体单元，用于根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体；

线性干扰衰减单元，用于根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

在其中一种实施例中，线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数，或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。

在其中一种实施例中，还包括：

平滑单元，用于对第一射线参数数据体进行平滑处理；

第二射线参数数据体单元具体用于：

根据线性干扰速度、该道集对应的射线参数和经过平滑处理的第一射线参数数据体，获得第二射线参数数据体。

在其中一种实施例中，道集为炮集、检波点道集、共中心点道集或十字排列道集。

综上，本发明实施例的线性干扰衰减系统先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时可以实现基于线性干扰衰减方法的全部或部分内容，例如，处理器执行计算机程序时可以实现如下内容：

对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集；

对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅；

确定每个时间样点的多个振幅中的最大值对应的射线参数；

根据每个道集中的每个频带的射线参数形成第一射线参数数据体；

获取其中一个道集的线性干扰速度和该道集对应的射线参数；

根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体；

根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

综上，本发明实施例的计算机设备先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时可以实现线性干扰衰减方法的全部或部分内容，例如，处理器执行计算机程序时可以实现如下内容：

对原始地震数据中的每个道集进行分频处理，每个道集均包括多个频带；原始地震数据包括多个道集；

对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果；线性拉东变换结果包括多个时间样点，每个时间样点对应多个振幅；

确定每个时间样点的多个振幅中的最大值对应的射线参数；

根据每个道集中的每个频带的射线参数形成第一射线参数数据体；

获取其中一个道集的线性干扰速度和该道集对应的射线参数；

根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体；

根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。

综上，本发明实施例的计算机可读存储介质先对原始地震数据中的每个道集进行分频处理，再对每个频带的原始地震数据进行线性拉东变换，得到每个频带的线性拉东变换结果，接着确定每个时间样点的多个振幅中的最大值对应的射线参数，根据每个道集中的每个频带的射线参数形成第一射线参数数据体；然后获取其中一个道集的线性干扰速度和该道集对应的射线参数，并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体，获得第二射线参数数据体，最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据，无需参数实验且易于实施，可以提高线性干扰去噪的效率和精度。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrativelogicalblock)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability)，上述的各种说明性部件(illustrativecomponents)，单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。

本发明实施例中所描述的各种说明性的逻辑块，或单元，或装置都可以通过通用处理器，数字信号处理器，专用集成电路(asic)，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。

本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于ram存储器、闪存、rom存储器、eprom存储器、eeprom存储器、寄存器、硬盘、可移动磁盘、cd-rom或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于asic中，asic可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。

在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于ram、rom、eeprom、cd-rom或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(dsl)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、dvd、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。