频谱分解方法及装置

技术领域

本发明涉及油气勘探技术领域，尤其涉及一种频谱分解方法及装置。

背景技术

频谱分解技术是一种频率域的解释方法，是地震属性分析中的重要组成部分，用频率随时间的变化关系来研究薄层结构、判断沉积环境。

目前，频谱分解技术的相关算法都是按照地震道进行逐道分析的，但是，在对地震道进行逐道分析的过程中，会受到地震数据中噪音的影响，从而导致频谱分解结果出现横向不稳定的问题，影响后续储层厚度预测、储层识别等作业精度。

发明内容

本发明实施例提供一种频谱分解方法，用以解决频谱分解结果出现横向不稳定的问题，进而保证后续储层厚度预测、储层识别等的作业精度，该方法包括：

根据预设工区的地震数据获取地震倾角数据；

对所述地震数据进行频谱分解计算；

根据频谱分解计算的结果，结合所述地震倾角数据获取频谱分解结果。

可选的，根据频谱分解计算的结果，结合所述地震倾角数据获取频谱分解结果，包括：

根据频谱分解计算的结果和所述地震倾角数据配置滤波器；

根据配置后的滤波器获取频谱分解结果。

可选的，所述地震倾角数据包括：地震倾角信息和地震倾向信息。

可选的，对所述地震数据进行频谱分解计算采用的方法为匹配追踪算法。

本发明实施例还提供一种频谱分解装置，用以解决频谱分解结果出现横向不稳定的问题，进而保证后续储层厚度预测、储层识别等的作业精度，该装置包括：

数据获取模块，用于根据预设工区的地震数据获取地震倾角数据；

计算模块，用于对所述地震数据进行频谱分解计算；

频谱分解模块，用于根据频谱分解计算的结果，结合所述地震倾角数据获取频谱分解结果。

可选的，频谱分解模块进一步用于：

根据频谱分解计算的结果和所述地震倾角数据配置滤波器；

根据配置后的滤波器获取频谱分解结果。

可选的，所述地震倾角数据包括：地震倾角信息和地震倾向信息。

可选的，对所述地震数据进行频谱分解计算采用的方法为匹配追踪算法。

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

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。

本发明实施例提供的频谱分解方法，通过根据预设工区的地震数据获取地震倾角数据，并对地震数据进行频谱分解计算，进而根据频谱分解计算的结果，结合地震倾角数据获取频谱分解结果，保证了频谱分解是在地震倾角数据的约束下进行了，解决了频谱分解结果出现横向不稳定的问题，进而保证了后续储层厚度预测、储层识别等的作业精度。

附图说明

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

图1为本发明实施例中频谱分解方法的流程图；

图2为本发明实施例中频谱分解装置的结构示意图；

图3为本发明实施例中预设工区的叠后地震体的地震剖面示意图；

图4为本发明实施例中预设工区的时间切片示意图；

图5为采用现有技术的频谱分解方法得到的预设工区数据剖面示意图；

图6为采用现有技术的频谱分解方法得到的预设工区数据平面示意图；

图7为本发明实施例中预设工区数据剖面示意图；

图8为本发明实施例中预设工区数据平面示意图；

图9为本发明实施例中时窗示例图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

图1为本发明实施例提供的频谱分解方法的流程图，如图1所示，该方法包括：

步骤101、根据预设工区的地震数据获取地震倾角数据。

在本实施例中，地震倾角数据包括：地震倾角信息和地震倾向信息。

具体地，该地震倾角和地震倾向的计算公式如下：

地震倾角s：

s＝(p

地震倾向φ：

其中，p为直角坐标系中x方向的倾角，q为直角坐标系中y方向的倾角。

步骤102、对所述地震数据进行频谱分解计算。

在本实施例中，为了确保频谱分解计算精度较高，可以采用匹配追踪算法对所述地震数据进行频谱分解计算。

步骤103、根据频谱分解计算的结果，结合所述地震倾角数据获取频谱分解结果。

在本实施例中，步骤103包括：

根据频谱分解计算的结果和所述地震倾角数据配置滤波器；

根据配置后的滤波器获取频谱分解结果。

具体实施时，可以用窗函数发定义滤波器，如图9所示，

由图1可知，本发明实施例提供的频谱分解方法，通过根据预设工区的地震数据获取地震倾角数据，并对地震数据进行频谱分解计算，进而根据频谱分解计算的结果，结合地震倾角数据获取频谱分解结果，保证了频谱分解是在地震倾角数据的约束下进行了，解决了频谱分解结果出现横向不稳定的问题，进而保证了后续储层厚度预测、储层识别等的作业精度。

下面以预设工区的具体示例对本发明进行说明：

图3和图4分别表示叠后地震体的地震剖面和时间切片。图3中的箭头指示强反射点，这是由冲积沉积系统的差异压实作用引起的，图4中的箭头所指为特征的主河道。

基于匹配-追踪方法，进行了谱分解之后，我们分别得到了图5和图6，波峰频率的地震剖面和和时间切片。由于谱分析对沉积异常比较敏感，我们发现河道具有强反射，可以在图5中划出河道边界，与周围区域相比，其特征为波峰频率响应较低。同时，由于尖峰异常的存在，我们也勉强可以进行横向构造追踪，当然，尖峰异常也可在图6的时间切片中区分主河道系统。因此，在地震解释中，有必要实现水平约束，以便在地震解释中赋予直接的地质意义。

图7和图8表示基于匹配-追踪方法，使用用户定义的不同滤波参数，使构造约束谱分解得到增强。图7中表示由河道差异压实作用引起强反射得到保留，尖峰异常大幅降低(用户定义滤波参数α＝0.25)，这使在图8的时间切片中识别河道成为可能。如果我们把用户定义滤波参数α设为1.0，用横向变化约束造导向，在图8中初始化了很多地质意义，这与图3中的地震剖面是相对应的。图8中的主河道成为可见河道，为解释人员提供了另一种进行地震解释的可靠证据。

基于同一发明构思，本发明实施例中还提供了一种频谱分解装置，如下面的实施例所述。由于频谱分解装置解决问题的原理与频谱分解方法相似，因此，频谱分解装置的实施可以参见频谱分解方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图2为本发明实施例提供的频谱分解装置的结构示意图，如图2所示，该装置包括：

数据获取模块201，用于根据预设工区的地震数据获取地震倾角数据。

计算模块202，用于对所述地震数据进行频谱分解计算。

频谱分解模块203，用于根据频谱分解计算的结果，结合所述地震倾角数据获取频谱分解结果。

在本发明实施例中，频谱分解模块203进一步用于：

根据频谱分解计算的结果和所述地震倾角数据配置滤波器；

根据配置后的滤波器获取频谱分解结果。

在本发明实施例中，所述地震倾角数据包括：地震倾角信息和地震倾向信息。

在本发明实施例中，对所述地震数据进行频谱分解计算采用的方法为匹配追踪算法。

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

本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

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