网格层析反演反射点确定方法及装置与流程

本发明涉及地震资料处理领域，尤其涉及一种网格层析反演反射点确定方法及装置。

背景技术：

基于网格层析反演的深度域速度建模过程中，涉及多个步骤。其中步骤之一就是需要确定网格层析反演射线追踪的反射点，也就是射线追踪的起点。射线的起点常选择在叠前深度偏移得到的叠加剖面强同相轴的地方。

在早期自动拾取技术成熟之前，常常采用手工拾取的方式进行反射点拾取。二维情况下，在叠加剖面上进行人工层位拾取，拾取的多套层位需要进行封闭处理，层位封闭的好处是：一方面它利于速度填充，不会使速度体出现空洞现象；另一方面能形成层位模型，层位模型上的层位值就是层析反演需要的反射点。三维情况下，也需要在叠加剖面上进行人工层位拾取，间隔几十条测线拾取一条测线，然后把拾取的同套层位在不同测线间进行插值成层位面，再对有相交的层位进行封闭成块，所有的交线封闭成层位体，俗称实体模型(solidmodel)。该方法有三个不足：1)人工拾取费时费力，消耗大量的人工，使深度域速度建模效率低下；2)拾取的层位是大套地层，对两个地层之间的强同相轴没有考虑，缺少了层析反演的数据，使得层析反演的效果不佳，因此拾取的反射点准确度低；3)在层位封闭过程中，会进行层位外推，外推出来的层位，可信度不高，同时，外推的层位常常是在工区的边缘，该处的地震数据信噪比低，不利于层析反演，反而会降低层析反演的效果，因此拾取的反射点准确度低。

技术实现要素：

本发明实施例提出一种网格层析反演反射点确定方法，用以确定网格层析反演射线追踪的反射点，效率高，准确度高，该方法包括：

从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；

基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；

根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；

基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；

确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。

本发明实施例提出一种网格层析反演反射点确定装置，用以确定网格层析反演射线追踪的反射点，效率高，准确度高，该装置包括：

波峰获得模块，用于从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；

第一过滤模块，用于基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；

相干能量谱获得模块，用于根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；

第二过滤模块，用于基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；

反射点确定模块，用于确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。

本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述网格层析反演反射点确定方法。

本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述网格层析反演反射点确定方法的计算机程序。

在本发明实施例中，从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。在上述过程中，通过第一次过滤处理，过滤掉一些不符合要求的波峰，然后获得过滤剩下的波峰对应的相关能量谱数据值，并进行第二次过滤处理，过滤掉一些不符合要求的相干能量谱数据值，因此，最后用获得的相干能量谱数据值对应的多个波峰，来确定反射点使得确定的反射点的准确度高；且上述过程自动完成，不需要人工拾取反射点，效率高。

附图说明

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

图1为本发明实施例中网格层析反演反射点确定方法的流程图；

图2为本发明实施例提出的网格层析反演反射点确定方法的详细流程图；

图3为本发明实施例中网格层析反演反射点确定装置的示意图。

具体实施方式

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

在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。

图1为本发明实施例中网格层析反演反射点确定方法的流程图，如图1所示，该方法包括：

步骤101，从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；

步骤102，基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；

步骤103，根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；

步骤104，基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；

步骤105，确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。

在本发明实施例中，在上述过程中，通过第一次过滤处理，过滤掉一些不符合要求的波峰，然后获得过滤剩下的波峰对应的相关能量谱数据值，并进行第二次过滤处理，过滤掉一些不符合要求的相干能量谱数据值，因此，最后用获得的相干能量谱数据值对应的多个波峰，来确定反射点使得确定的反射点的准确度高；且上述过程自动完成，不需要人工拾取反射点，效率高。

具体实施时，在步骤101中，波峰的获得方法是如果当前深度点上的地震数据振幅值同时大于上一个深度点的地震数据振幅值和下一个深度点的地震数据振幅值，则认为该深度点上对应着一个波峰，即

ai＞ai-1且ai＞ai+1(1)

其中，ai为第i个深度点对应的地震数据振幅值。其中2≤i＜n-1，n为最大深度点数。

在一实施例中，在基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰之前，还包括：

对所有波峰的振幅值进行归一化处理；

基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰，包括：

基于第一阈值，对归一化处理后的所有波峰进行过滤处理，获得过滤后的多个波峰。

在上述实施例中，对所有波峰的振幅值进行归一化处理的目的是便于对波峰进行过滤处理，对所有波峰进行归一化处理的过程为：

采用如下公式，寻找最大波峰的振幅值：

bmax＝max{b1,b2,…bi…,bn}(2)

bmax为最大波峰的振幅值，bi为第i个波峰对应的地震数据振幅值，n为波峰的个数。

用每个波峰的振幅值除以该最大波峰振幅值，即为归一化后的每个波峰的振幅值，公式如下：

其中，ci为归一化后的第i个波峰的振幅值。

然后基于第一阈值，对归一化处理后的所有波峰进行过滤处理，获得过滤后的多个波峰，其中，第一阈值m1可以是一个百分比值，对归一化后的波峰振幅值小于该第一阈值m1的波峰进行剔除，对归一化后的波峰振幅值大于该第一阈值m1的波峰振幅进行保留，从而获得过滤后的多个波峰。

发明人发现，叠前深度偏移后的叠加剖面是对道集进行叠加的结果，若直接从上述过滤后的多个波峰拾取反射点，而不考虑道集中存在不合理波峰振幅的影响，则确定的反射点的准确度还是不够高，因此，本发明实施例考虑了上述道集。

在一实施例中，根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值，包括：

对过滤后的每一波峰，提取该波峰对应的深度位置的满足预设要求的多个alpha值；

对该波峰对应的深度位置的满足预设要求的多个alpha值对应的道集进行扫描，获得该波峰的多个相干能量谱数据值；

对该波峰的多个相干能量谱数据值进行加权平均处理，获得该波峰对应的相干能量谱数据值。

在上述实施例中，alpha值可以是多个，为了获得更准确的反射点，这里的满足预设要求的多个alpha值可以是最优alpha值以及该最优alpha值临近的alpha值；且临近的alpha值可以有多个，alpha值在相干能量谱数据值中自动拾取得到，它是网格层析反演流程中的一步；另外，多个相干能量谱数据值在网格层析反演流程中已经扫描得到，它是网格层析反演流程中的一步。之后，对该波峰的多个相干能量谱数据值进行加权平均处理，获得该波峰的相干能量谱数据值。这里，在提取该波峰对应的深度位置的满足预设要求的多个alpha值时，还可以考虑该波峰对应的cmp点编号，即提取该波峰对应的cmp点编号和深度位置的满足预设要求的多个alpha值，cmp点编号的作用是标识，方便记录波峰、深度位置和alpha的关系。

假设最优的alpha为α，其临近的前面一个alpha为α前，其临近的后面一个alpha为α后，对三个alpha对应的道集进行扫描，得到三个相干能量谱数据值为s，s前，s后，对其进行加权平均：

其中ω1，ω2和ω3为加权系数，且ω1+ω2+ω3＝1；

sα为加权平均后的相干能量谱数据值。

此处选取的是最优alpha前后临近的alpha，本发明也可以选取临近更多的alpha值对应的相干能量谱数据值进行加权计算。

当然可以理解的是，还可以有多种对该波峰的多个相干能量谱数据值进行处理，从而得到该波峰对应的相干能量谱数据值的办法，相关变化例均应落入本发明的保护范围。

在一实施例中，在基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值之前，还包括：

对过滤后的多个波峰对应的相干能量谱进行归一化处理；

基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值，包括：

基于第二阈值，对归一化处理后的多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值。

在上述实施例中，对过滤后的多个波峰对应的相干能量谱数据值进行归一化处理的目的是便于对相干能量谱数据值进行过滤处理，过程包括：

采用如下公式，寻找最大的相干能量谱数据值：

sα，max＝max{sα，1,sα，2,…sα，i…,sα，n}(5)

sα，max为最大相干能量谱数据值，sα，i为第i个相干能量谱数据值，n为相干能量谱数据值的总个数。

用每一个过滤后的波峰位置对应的相干能量谱数据值除以该最大相干能量谱数据值，公式如下：

其中dα，i为归一化后的第i个相干能量谱数据值。

在步骤104中，第二阈值m2是一个百分比值，对归一化后的相干能量谱数据值小于第二阈值m2的相干能量谱数据值进行剔除，对归一化后的相干能量谱数据值大于该第二阈值m2的相干能量谱数据值进行保留，获得过滤后的多个相干能量谱数据值。

基于上述实施例，本发明提出如下一个实施例来说明网格层析反演反射点确定方法的详细流程，图2为本发明实施例提出的网格层析反演反射点确定方法的详细流程图，如图2所示，在一实施例中，网格层析反演反射点确定方法的详细流程包括：

步骤201，从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；

步骤202，对所有波峰的振幅值进行归一化处理；

步骤203，基于第一阈值，对归一化处理后的所有波峰进行过滤处理，获得过滤后的多个波峰；

步骤204，对过滤后的每一波峰，提取该波峰对应的深度位置的满足预设要求的多个alpha值；

步骤205，对该波峰对应的深度位置的满足预设要求的多个alpha值对应的道集进行扫描，获得该波峰的多个相干能量谱数据值；

步骤206，对该波峰的多个相干能量谱数据值进行加权平均处理，获得该波峰对应的相干能量谱数据值；

步骤207，对过滤后的多个波峰对应的相干能量谱数据值进行归一化处理；

步骤208，基于第二阈值，对归一化处理后的多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；

步骤209，确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。

当然，可以理解的是，上述网格层析反演反射点确定方法的详细流程还可以有其他变化例，相关变化例均应落入本发明的保护范围。

综上所述，在本发明实施例提出的方法中，从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。在上述过程中，通过第一次过滤处理，过滤掉一些不符合要求的波峰，然后获得过滤剩下的波峰对应的相关能量谱数据值，并进行第二次过滤处理，过滤掉一些不符合要求的相干能量谱数据值，因此，最后用获得的相干能量谱数据值对应的多个波峰，来确定反射点使得确定的反射点的准确度高；且上述过程自动完成，不需要人工拾取反射点，效率高。

基于同样的发明构思，本发明实施例还提供了一种网格层析反演反射点确定装置，如下面的实施例所述。由于这些解决问题的原理与网格层析反演反射点确定方法相似，因此装置的实施可以参见方法的实施，重复之处不在赘述。

图3为本发明实施例中网格层析反演反射点确定装置的示意图，如图3所示，该装置包括：

波峰获得模块301，用于从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；

第一过滤模块302，用于基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；

相干能量谱获得模块303，用于根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；

第二过滤模块304，用于基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；

反射点确定模块305，用于确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。

在一实施例中，所述装置还包括第一预处理模块306，用于：

对所有波峰的振幅值进行归一化处理；

第一过滤模块302具体用于：

基于第一阈值，对归一化处理后的所有波峰进行过滤处理，获得过滤后的多个波峰。

在一实施例中，相干能量谱获得模块303具体用于：

对过滤后的每一波峰，提取该波峰对应的深度位置的满足预设要求的多个alpha值；

对该波峰对应的深度位置的满足预设要求的多个alpha值对应的道集进行扫描，获得该波峰的多个相干能量谱数据值；

对该波峰的多个相干能量谱数据值进行加权平均处理，获得该波峰对应的相干能量谱数据值。

在一实施例中，所述装置还包括第二预处理模块307，用于：

对过滤后的多个波峰对应的相干能量谱数据值进行归一化处理；

第二过滤模块304具体用于：

基于第二阈值，对归一化处理后的多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值。

综上所述，在本发明实施例提出的装置中，从叠前深度偏移后的叠加剖面上，获得地震数据的所有波峰；基于第一阈值，对所有波峰进行过滤处理，获得过滤后的多个波峰；根据过滤后的每个波峰对应的深度位置，获得过滤后的每个波峰对应的相干能量谱数据值；基于第二阈值，对多个波峰对应的相干能量谱数据值进行过滤处理，获得过滤后的多个相干能量谱数据值；确定过滤后的多个相干能量谱数据值对应的多个波峰的位置为网格层析反演反射点。在上述过程中，通过第一次过滤处理，过滤掉一些不符合要求的波峰，然后获得过滤剩下的波峰对应的相关能量谱数据值，并进行第二次过滤处理，过滤掉一些不符合要求的相干能量谱数据值，因此，最后用获得的相干能量谱数据值对应的多个波峰，来确定反射点使得确定的反射点的准确度高；且上述过程自动完成，不需要人工拾取反射点，效率高。

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

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

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

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

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