图案引导倾角估计制造技术

本说明书中描述的主题的创新方面可以体现在包括获取地震数据图像的方法中。将第一平面波破坏滤波器倾角估计应用于地震数据图像以生成初始倾角模型。使用初始倾角模型将第二平面波破坏滤波器倾角估计应用于地震数据图像以生成图案引导倾角估计。图案引导倾角估计存储在数据存储装置中。存储在数据存储装置中。存储在数据存储装置中。

【技术实现步骤摘要】
【国外来华专利技术】图案引导倾角估计
[0001]在先申请的优先权的引用
[0002]本申请要求于2019年2月25日递交的美国专利申请No.16/284,767的优先权，其全部内容通过引用并入本文。


[0003]本公开涉及地震数据和图像处理。

技术介绍

[0004]结构信息往往是地震图像最重要的内容，可以提取结构信息进行面向结构的数据处理，如平滑、插值和选取。沿结构平滑可以增强结构特征，同时保留重要的不连续性，例如断层或通道。这种平滑也被用作基于逆(inversion)的地震成像中的模型约束。沿结构插值可以重建具有有意义的地质意义的地震图像/模型，并且它已被用于测井插值。沿结构选取可以选取用于结构解释的层位，或选取共成像集(CIG)中的残余时差用于叠前成像。
[0005]结构信息可以用局部倾角属性来表征，可通过几种方法估计，包括相似扫描方法和局部结构张量。倾角估计方法可以从不包括强烈冲突和陡峭结构的图像中估计准确的倾角。

技术实现思路

[0006]本公开讨论并入结构取向的图案信息以指导倾角估计和减轻走样问题。当图像具有非常陡峭的结构时，这种倾角估计方法可能会遇到走样问题，这会使倾角估计具有挑战性。来自走样部分的估计主导倾角是不希望的，因为它遵循虚假走样结构而不是真实结构。在一些实现方式中，图案信息被合并以提供平面波破坏(PWD)滤波器的初始倾角。PWD滤波器求解非线性逆问题，并且解法基于提供于逆的初始倾角。如果初始倾角与图像有关，那么估计更接近真实结构的倾角的可能性很大，而不是走样(aliased)结构的倾角。因此，PWD会导致包含图案信息的新倾角，并称为图案引导倾角。在一些示例中，当合成数据走样时，可以通过提供不同的初始倾角来估计倾角属于真实事件或走样事件。
[0007]本说明书中描述的主题的创新方面可以体现在包括获取地震数据图像的方法中。将第一平面波破坏滤波器倾角估计应用于地震数据图像以生成初始倾角模型。使用初始倾角模型将第二平面波破坏滤波器倾角估计应用于地震数据图像以生成图案引导倾角估计。图案引导倾角估计存储在数据存储装置中。
[0008]这些方面的其他实现方式包括：对应的系统、装置和计算机程序，被配置为执行在计算机存储设备上编码的方法的动作。
[0009]这些实现方式和其他实现方式可以各自可选地包括以下特征中的一个或多个。例如，该方法可以包括：基于地震数据图像计算相干图；用预定阈值裁剪相干图；以及响应于裁剪生成掩膜算子。应用第一平面波破坏滤波器倾角估计可以包括：将第一加权因子应用于地震数据图像的走样影响区域，并且将第二加权因子应用于地震数据图像的无走样区
域。第一加权因子的值为0，并且第二加权因子的值为1。初始倾角模型仅包括地震数据图像的无走样区域。图案引导倾角估计是非线性逆估计。使用图案引导倾角估计，估计面向结构的插值目标图像。
[0010]可以实施本说明书中描述的主题的特定实现方式，以实现以下优点中的一个或多个。例如，主题的实现方式减轻了在面向结构的数据处理中可能存在的走样问题。
[0011]在本公开中描述的主题的一个或多个实现方式的细节在附图和描述中阐述。通过说明书、附图和权利要求书，所述主题的其他特征、方面和优点将变得明显。
附图说明
[0012]图1是用于图案引导倾角估计的系统的示意图。
[0013]图2A是原始地震数据图像。
[0014]图2B是作为面向结构的插值的输入图像的子采样图像。
[0015]图3A是包括估计倾角的图像。
[0016]图3B是插值结果的图像。
[0017]图4A是使用掩膜算子的初始倾角估计的图像。
[0018]图4B是使用初始倾角估计的新倾角估计的图像。
[0019]图5A是用于测量倾角估计的可靠性的相干图。
[0020]图5B是通过裁剪相干图的掩膜算子的图像
[0021]图6A是插值结果的图像。
[0022]图6B是插值结果与原始图像的差异的图像。
[0023]图7表示图案引导倾角估计的流程图。
[0024]图8表示用于实现本文描述的技术的示例计算环境。
具体实施方式
[0025]本公开描述了用于减轻在面向结构的数据处理中存在的走样问题的图案引导倾角估计。具体而言，图案引导倾角估计可以包括三个平面波破坏(PWD)滤波器。第一倾角估计可以用于生成掩膜算子(mask operator)以区分无走样数据和受走样影响的数据。第二倾角估计仅使用无走样数据进行，并输出用于第三(最终)倾角估计逆的初始模型。在一些示例中，例如对表面偏移共成像集(CIG)识别残余时差(residual moveout)，掩膜算子的生成(第一倾角估计)可以是可选的。
[0026]本公开描述了用于图1所示的图案引导倾角估计的计算系统100。计算系统100包括可以通过一个或多个网络(未示出)与一个或多个其他计算系统(未示出)通信的计算设备102。系统100还包括数据存储装置106，其中计算设备102与数据存储装置106通信连接。
[0027]倾角估计可以表述为正则化非线性逆问题，并定义为等式1中的最小二乘法：
[0028]MD(σ)≈0
ꢀꢀꢀꢀꢀ
[1][0029]在等式1中，在对倾角σ平滑正则化目标的情况下，D(σ)是PWD滤波器，D表示已知数据，M是掩膜算子，并且近似相等表示MD(σ)的幂的最小化。如果估计的倾角准确，则应用PWD滤波器后的数据应该没有能量(或具有最小能量)。为此，可以使用解析线性化来求解PWD滤波器D中的倾角σ。
[0030]图2A表示深度迁移的三维图像200的一部分，包括结构折叠和角度不整合。如图2B所示，通过对原始图像(图像200)进行子采样来生成输入图像202。图2B表示包括如区域204所指出的走样问题在内的图像202。图3A表示包括使用PWD滤波器的主测线(inline)倾角估计的图像300。图3B表示在图像300的面向结构的插值之后的图像302。为此，如图3B所示，走样会导致插值结果中事件的错误方向；例如，2700米深度和10，500米主测线距离附近的事件的倾角与相邻事件的倾角相反。
[0031]返回参考图1，计算设备102可以接收地震数据120。地震数据120可以包括地震数据图像122。计算设备102可以将第一平面波破坏(PWD)滤波器倾角估计132应用于地震数据图像122以生成初始倾角模型134。图4A表示图像400，其图示了已经通过使用图案信息平滑正则化来调整(填充)的走样影响区域中的倾角。在一些示例中，在第一PWD滤波器倾角估计132中仅使用无走样数据，其中通过在PWD的逆中平滑正则化，通过“图案倾角”来自动扩展受走样影响的区域中的倾角。换言之，在第一PWD滤波器倾角估计132中应用掩膜算子130。即，掩膜算子130在第一PWD滤波器倾角估计中将值1应用于地震数据图像122的无走样区域，并且将值0应用于地震数据图像122的走样影响区域。在一些示例中，初本文档来自技高网...

【技术保护点】

【技术特征摘要】
【国外来华专利技术】1.一种方法，包括：获取地震数据图像；将第一平面波破坏滤波器倾角估计应用于所述地震数据图像以生成初始倾角模型；使用所述初始倾角模型将第二平面波破坏滤波器倾角估计应用于所述地震数据图像以生成图案引导倾角估计；和将所述图案引导倾角估计存储在数据存储装置中。2.如权利要求1所述的方法，还包括：根据所述地震数据图像计算相干图；用预定的阈值裁剪所述相干图；和响应于所述裁剪，生成掩膜算子。3.如权利要求2所述的方法，还包括：在所述第一平面波破坏滤波器倾角估计中应用所述掩膜算子，包括：将第一加权因子应用于所述地震数据图像的走样影响区域，以及将第二加权因子应用于所述地震数据图像的无走样区域。4.如权利要求3所述的方法，其中，所述第一加权因子的值为0，并且所述第二加权因子的值为1。5.如权利要求3所述的方法，其中，所述初始倾角模型仅包括所述地震数据图像的无走样区域。6.如权利要求1所述的方法，其中，所述图案引导倾角估计是非线性逆估计。7.如权利要求1所述的方法，使用所述图案引导倾角估计来估计面向结构的插值目标图像。8.一种系统，包括：一个或多个处理器；以及非暂时性计算机可读存储介质，耦合到所述一个或多个处理器并存储由所述一个或多个处理器执行的编程指令，所述编程指令指示所述一个或多个处理器执行以下操作：获取地震数据图像；将第一平面波破坏滤波器倾角估计应用于所述地震数据图像以生成初始倾角模型；使用所述初始倾角模型将第二平面波破坏滤波器倾角估计应用于所述地震数据图像以生成图案引导倾角估计；和将所述图案引导倾角估计存储在数据存储装置中。9.如权利要求8所述的系统，所述操作还包括：根据所述地震数据图像计算相干图；用预定的阈值裁剪所述相干图；和响应于所述裁剪，生成掩膜算子。10.如权利要求9所述的系统，所述操作还包括：在所述第一平...

【专利技术属性】
技术研发人员：赵杨，张厚竹，
申请(专利权)人：沙特阿拉伯石油公司，
类型：发明
国别省市：