本发明专利技术提供了一种不连续地质体检测方法、装置及电子设备，获取待处理区域的地震波场共偏移距数据；基于预先确定的局部倾角场的估算值，确定地震波场共偏移距数据中的绕射波数据；其中，局部倾角场的估算值基于局部平面波方程以及预先确定的滤波原理确定；基于绕射波数据，确定待处理区域的成像结果。本发明专利技术中，从理论出发确定局部倾角场的估算值，降低了实际数据对局部倾角场的干扰，提高了对绕射波的分离效果，从而提高了不连续地质体检测的质量。从而提高了不连续地质体检测的质量。从而提高了不连续地质体检测的质量。
【技术实现步骤摘要】
不连续地质体检测方法、装置及电子设备
[0001]本专利技术涉及地震数据处理
，尤其是涉及一种射波成像方法、装置及电子设备。
技术介绍
[0002]相关技术中，通常依赖于地震数据来估计局部倾角，然后基于局部倾角对地震波场共偏移距数据中的绕射波进行分离。然而，该方式得到的局部倾角受数据质量影响较大，算法稳定性较差，对绕射波分离的效果较差，导致不连续地质体检测质量较差。
技术实现思路
[0003]有鉴于此，本专利技术的目的在于提供一种射波成像方法、装置及电子设备，以提高对绕射波的分离效果，提高不连续地质体检测的质量。
[0004]第一方面，本专利技术实施例提供了一种不连续地质体检测方法，包括：获取待处理区域的地震波场共偏移距数据；基于预先确定的局部倾角场的估算值，确定地震波场共偏移距数据中的绕射波数据；局部倾角场的估算值基于局部平面波方程以及预先确定的滤波原理确定；基于绕射波数据，确定待处理区域的成像结果。
[0005]进一步地，上述局部倾角场的估算值通过以下方式确定：基于局部平面波方程，确定地震波场数据的偏导算子；基于预先确定的滤波原理，确定偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数；基于偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数，确定局部倾角场的估算值。
[0006]进一步地，基于预先确定的滤波原理，确定偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数的步骤，包括：基于有限脉冲响应滤波器，分别计算偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数。
[0007]进一步地，上述偏导算子的频率响应函数包括x方向偏导算子的频响函数以及y方向偏导算子的频响函数；希尔伯特的频率响应函数包括y方向的希尔伯特变换的频响函数以及x方向的希尔伯特变换的响应函数；局部倾角场包括多个设定频率的局部倾角；基于偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数，确定局部倾角场的估算值的步骤，包括：将x方向偏导算子的频响函数的傅里叶逆变换参数以及y方向偏导算子的频响函数的傅里叶逆变换参数之商，确定为局部倾角场在设定频率的局部倾角；基于偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数之间的关系，将x方向的希尔伯特变换的响应函数与y方向的希尔伯特变换的频响函数之商确定为局部倾角的估算值。
[0008]进一步地，基于预先确定的局部倾角场的估算值，确定地震波场共偏移距数据中的绕射波数据的步骤，包括：基于局部倾角场的估算值，生成平面波分解算子；基于平面波分解算子及地震波场共偏移距数据，得到绕射波数据。
[0009]第二方面，本专利技术实施例还提供一种不连续地质体检测装置，包括：数据获取模块，用于获取待处理区域的地震波场共偏移距数据；绕射波确定模块，用于基于预先确定的
局部倾角场的估算值，确定地震波场共偏移距数据中的绕射波数据；局部倾角场的估算值基于局部平面波方程以及预先确定的滤波原理确定；成像模块，用于基于绕射波数据，确定待处理区域的成像结果。
[0010]进一步地，上述装置还包括局部倾角场确定模块；局部倾角场确定模块用于：基于局部平面波方程，确定地震波场数据的偏导算子；基于预先确定的滤波原理，确定偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数；基于偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数，确定局部倾角场的估算值。
[0011]进一步地，上述绕射波确定模块还用于：基于局部倾角场的估算值，生成平面波分解算子；基于平面波分解算子及地震波场共偏移距数据，得到绕射波数据。
[0012]第三方面，本专利技术实施例还提供一种电子设备，包括处理器和存储器，存储器存储有能够被处理器执行的机器可执行指令，处理器执行机器可执行指令以实现上述不连续地质体检测方法。
[0013]第四方面，本专利技术实施例还提供机器可读存储介质，机器可读存储介质存储有机器可执行指令，机器可执行指令在被处理器调用和执行时，机器可执行指令促使处理器实现上述不连续地质体检测方法。
[0014]本专利技术实施例带来了以下有益效果：
[0015]本专利技术实施例提供了一种不连续地质体检测方法、装置及电子设备，获取待处理区域的地震波场共偏移距数据；基于预先确定的局部倾角场的估算值，确定地震波场共偏移距数据中的绕射波数据；其中，局部倾角场的估算值基于局部平面波方程以及预先确定的滤波原理确定；基于绕射波数据，确定待处理区域的成像结果。该方法中，从理论出发确定局部倾角场的估算值，降低了实际数据对局部倾角场的干扰，提高了对绕射波的分离效果，从而提高了不连续地质体检测的质量。
[0016]本专利技术的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本专利技术而了解。本专利技术的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
[0017]为使本专利技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
[0018]为了更清楚地说明本专利技术具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施方式，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0019]图1为本专利技术实施例提供的一种不连续地质体检测方法的流程图；
[0020]图2为本专利技术实施例提供的另一种不连续地质体检测方法的流程图；
[0021]图3为本专利技术实施例提供的一种不连续地质体检测装置的结构示意图；
[0022]图4为本专利技术实施例提供的一种电子设备的结构示意图。
具体实施方式
[0023]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0024]地下不连续地质体，如空洞、断层、裂缝等，往往和矿产资源分布具有紧密的关系，对非均匀不连续地质体的精确定位，可有效提高勘探成功率，降低成本，避免可能发生的地质灾害，规避风险。绕射波是小尺度地质体的地震响应，包含了小尺度地质体的构造信息，可以用来精确定位非均匀不连续地质体，提供更强的地下空间照明度。但相对反射波而言，绕射波在传播过程中衰减较快，能量较弱，容易被强能量的反射波掩盖。因此，需要将绕射波从反射波中分离，以增强绕射响应，从而进行绕射波高精度成像，对小尺度地质体精确定位。
[0025]现有方法中，传统的局部倾角估计方法依赖于数据本身，受数据质量影响较大，算法稳定性较差，特别是对于含噪音的低信噪比数据，绕射波分离结果不理想。
[0026]基于此，本专利技术实施例提供的一种不连续地质体检测方法、装置以及电子设备，可以应用于各种地震数据处理场景。本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种不连续地质体检测方法，其特征在于，包括：获取待处理区域的地震波场共偏移距数据；基于预先确定的局部倾角场的估算值，确定所述地震波场共偏移距数据中的绕射波数据；所述局部倾角场的估算值基于局部平面波方程以及预先确定的滤波原理确定；基于所述绕射波数据，确定所述待处理区域的成像结果。2.根据权利要求1所述的方法，其特征在于，所述局部倾角场的估算值通过以下方式确定：基于局部平面波方程，确定地震波场数据的偏导算子；基于预先确定的滤波原理，确定所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数；基于所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数，确定局部倾角场的估算值。3.根据权利要求2所述的方法，其特征在于，基于预先确定的滤波原理，确定所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数的步骤，包括：基于有限脉冲响应滤波器，分别计算所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数。4.根据权利要求2所述的方法，其特征在于，所述偏导算子的频率响应函数包括x方向偏导算子的频响函数以及y方向偏导算子的频响函数；所述希尔伯特的频率响应函数包括y方向的希尔伯特变换的频响函数以及x方向的希尔伯特变换的响应函数；所述局部倾角场包括多个设定频率的局部倾角；基于所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数，确定局部倾角场的估算值的步骤，包括：将所述x方向偏导算子的频响函数的傅里叶逆变换参数以及所述y方向偏导算子的频响函数的傅里叶逆变换参数之商，确定为所述局部倾角场在设定频率的局部倾角；基于所述偏导算子的频率响应函数，以及希尔伯特变换的频率响应函数之间的关系，将所述x方向的希尔伯特变换的响应函数与所述y方向的希尔伯特变换的频响函数之商确定为所述局部倾角的...
【专利技术属性】
技术研发人员：林朋，彭苏萍，向阳，崔晓芹，李闯建，柳倩男，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。