不连续属性体的提取方法、装置及电子设备制造方法及图纸

本发明专利技术提供了一种不连续属性体的提取方法、装置及电子设备，涉及地震成像技术领域，该方法包括：获取待处理的地震偏移数据；计算地震偏移数据的倾角场信息；基于地震偏移数据、倾角场信息和反偏移核函数，生成反偏移结果；其中，反偏移核函数包括反极性滤波公式；对反偏移结果进行偏移处理，以提取出地震偏移数据对应的不连续属性体。本发明专利技术可以有效提高不连续属性体的提取精度。

Extraction Method, Device and Electronic Equipment of Discontinuous Attribute Body

【技术实现步骤摘要】
不连续属性体的提取方法、装置及电子设备
本专利技术涉及地震成像
，尤其是涉及一种不连续属性体的提取方法、装置及电子设备。
技术介绍
根据惠更斯原理，当将地下地质体视为震源向地表激发地震波时，不连续地质体和连续地质体会以不同形式进行扩散。不连续地质体对油气运移和煤炭开采的安全问题紧密相关，通过对不连续地质体进行准确识别，可以有助于提高油气开采效率，减少开采成本，以及降低煤炭开采过程中的地质危害。其中，通过对地震偏移结果进行反偏移-偏移操作，从而实现不连续属性体的提取，但是由于在反偏移过程中连续属性体和不连续属性体的波场将同时被接收，导致得到的不连续属性体中仍然存在连续属性体，从而导致不连续属性体的提取精度较低。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种不连续属性体的提取方法、装置及电子设备，可以有效提高不连续属性体的提取精度。第一方面，本专利技术实施例提供了一种不连续属性体的提取方法，包括：获取待处理的地震偏移数据；计算所述地震偏移数据的倾角场信息；基于所述地震偏移数据、所述倾角场信息和反偏移核函数，生成反偏移结果；其中，所述反偏移核函数包括反极性滤波公式；对所述反偏移结果进行偏移处理，以提取出所述地震偏移数据对应的不连续属性体。结合第一方面，本专利技术实施例提供了第一方面的第一种可能的实施方式，其中，所述计算所述地震偏移数据的倾角场信息的步骤，包括：利用倾角估计法计算所述地震偏移数据的倾角场信息。结合第一方面，本专利技术实施例提供了第一方面的第二种可能的实施方式，其中，所述反偏移核函数表示为：U(t,s,r)＝∫w(s,x,r)I(x)δ(t-td(s,x,r))dx其中，U(t,s,r)表示所述反偏移结果，t表示第一双程旅行时间，s表示炮点信息，r表示检波点信息，x表示地质成像点信息，w(s,x,r)表示所述反极性滤波公式，I(x)表示地震偏移数据，δ(t-td(s,x,r))表示狄拉克函数，td(s,x,r)表示第二双程旅行时间。结合第一方面的第二种可能的实施方式，本专利技术实施例提供了第一方面的第三种可能的实施方式，其中，所述反极性滤波公式w(s,x,r)表示为：w(s,x,r)＝sgn(θ(s,x,r)-θ0(x))sin(θ(s,x,r)-θ0(x))其中，sgn(θ(s,x,r)-θ0(x))表示符号函数，θ(s,x,r)表示观测倾角，θ0(x)表示所述倾角场信息。结合第一方面的第三种可能的实施方式，本专利技术实施例提供了第一方面的第四种可能的实施方式，其中，所述观测倾角θ(s,x,r)表示为：其中，Ps表示所述炮点到所述地质成像点的矢量，Pr表示所述检波点到所述地质成像点的矢量，z表示单位矢量。结合第一方面的第二种可能的实施方式，本专利技术实施例提供了第一方面的第五种可能的实施方式，其中，所述第二双程旅行时间td(s,x,r)表示为：其中，xs表示所述所述炮点的横坐标，ys表示所述炮点的纵坐标，x0表示所述地质成像点的横坐标，y0表示所述地质成像点的纵坐标，yr表示所述检波点的纵坐标，h表示所述地质成像点的深度值，vrms表示所述地质成像点的均方根速度。结合第一方面的第二种可能的实施方式，本专利技术实施例提供了第一方面的第六种可能的实施方式，其中，所述对所述反偏移结果进行偏移处理，以提取出所述地震偏移数据对应的不连续属性体的步骤，包括：按照以下公式提取所述地震偏移数据对应的不连续属性体D(x)：D(x)＝∫U(t,s,r)δ(t-td(s,x,r))dtdrds。第二方面，本专利技术实施例还提供一种不连续属性体的提取装置，包括：偏移数据获取模块，用于获取待处理的地震偏移数据；倾角场计算模块，用于计算所述地震偏移数据的倾角场信息；反偏移结果生成模块，用于基于所述地震偏移数据、所述倾角场信息和反偏移核函数，生成反偏移结果；其中，所述反偏移核函数包括反极性滤波公式；不连续属性体提取模块，用于对所述反偏移结果进行偏移处理，以提取出所述地震偏移数据对应的不连续属性体。第三方面，本专利技术实施例还提供一种电子设备，所述电子设备包括存储器以及处理器，所述存储器用于存储支持处理器执行第一方面至第一方面的第六种可能的实施方式任一项所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。第四方面，本专利技术实施例还提供一种计算机存储介质，用于储存为第一方面至第一方面的第六种可能的实施方式任一项所述方法所用的计算机软件指令。本专利技术实施例带来了以下有益效果：本专利技术实施例提供的一种不连续属性体的提取方法、装置及电子设备，首先获取待处理的地震偏移数据，并计算该地震偏移数据的倾角场信息，然后基于地震偏移数据、倾角场信息以及包含有反极性滤波公式的反偏移核函数生成反偏移结果，通过对反偏移结果进行偏移处理，以提取出地震数据对应的不连续属性体。本专利技术实施例利用不连续属性体不受菲涅尔带限制的性质，对反偏移核函数进行修改，得到包含有反极性滤波公式的反偏移核函数，并通过反极性滤波公式压制连续属性体，从而使得到的不连续属性体中不会掺杂有连续属性体，进而有效提高了不连续属性体的提取精度。本专利技术的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本专利技术而了解。本专利技术的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。为使本专利技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。附图说明为了更清楚地说明本专利技术具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例提供的一种不连续属性体的提取方法的流程示意图；图2为本专利技术实施例提供的另一种不连续属性体的提取方法的流程示意图；图3为本专利技术实施例提供的一种不连续属性体的提取装置的结构示意图；图4为本专利技术实施例提供的一种电子设备的结构示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合实施例对本专利技术的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。目前，由于在反偏移过程中连续属性体和不连续属性体的波场将同时被接收，导致得到的不连续属性体中仍然存在连续属性体，从而导致不连续属性体的提取精度较低，基于此，本专利技术实施提供一种不连续属性体的提取方法、装置及电子设备，可以有效提高不连续属性体的提取精度。为便于对本实施例进行理解，首先对本专利技术实施例所公开的一种不连续属性体的提取方法进行详细介绍，参见图1所示的一种不连续属性体的提取方法的流程示意图，该方法可以包括以下步骤：步骤S102，获取待处理的地震偏移数据。其中，地震偏移数据是通过对地震数据进行处理后得到的，在具体实施时，可以采用均值滤波方法对地震数据进行噪音压制，并基于噪音压制后的地震数据得到地震偏移数据。步骤S104，计算地震偏移数据的倾角场信息。因为地震偏移数据的倾角场信息对反偏移本文档来自技高网...

【技术保护点】
1.一种不连续属性体的提取方法，其特征在于，包括：获取待处理的地震偏移数据；计算所述地震偏移数据的倾角场信息；基于所述地震偏移数据、所述倾角场信息和反偏移核函数，生成反偏移结果；其中，所述反偏移核函数包括反极性滤波公式；对所述反偏移结果进行偏移处理，以提取出所述地震偏移数据对应的不连续属性体。

【技术特征摘要】
1.一种不连续属性体的提取方法，其特征在于，包括：获取待处理的地震偏移数据；计算所述地震偏移数据的倾角场信息；基于所述地震偏移数据、所述倾角场信息和反偏移核函数，生成反偏移结果；其中，所述反偏移核函数包括反极性滤波公式；对所述反偏移结果进行偏移处理，以提取出所述地震偏移数据对应的不连续属性体。2.根据权利要求1所述的方法，其特征在于，所述计算所述地震偏移数据的倾角场信息的步骤，包括：利用倾角估计法计算所述地震偏移数据的倾角场信息。3.根据权利要求1所述的方法，其特征在于，所述反偏移核函数表示为：U(t,s,r)＝∫w(s,x,r)I(x)δ(t-td(s,x,r))dx其中，U(t,s,r)表示所述反偏移结果，t表示第一双程旅行时间，s表示炮点信息，r表示检波点信息，x表示地质成像点信息，w(s,x,r)表示所述反极性滤波公式，I(x)表示地震偏移数据，δ(t-td(s,x,r))表示狄拉克函数，td(s,x,r)表示第二双程旅行时间。4.根据权利要求3所述的方法，其特征在于，所述反极性滤波公式w(s,x,r)表示为：w(s,x,r)＝sgn(θ(s,x,r)-θ0(x))sin(θ(s,x,r)-θ0(x))其中，sgn(θ(s,x,r)-θ0(x))表示符号函数，θ(s,x,r)表示观测倾角，θ0(x)表示所述倾角场信息。5.根据权利要求4所述的方法，其特征在于，所述观测倾角θ(s,x,r)表示为：其中，Ps表示所述炮点到所述地质成像点的矢量，Pr表示所述检...

【专利技术属性】
技术研发人员：李闯建，彭苏萍，赵惊涛，崔晓芹，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京,11