绕射波分离方法、装置和电子设备制造方法及图纸

本发明专利技术提供了一种绕射波分离方法、装置和电子设备，其中，该方法包括：获取零偏移距道集；对零偏移距道集进行处理，得到地震波的瞬时相位；瞬时相位包括多个采样点；对每个采样点执行以下步骤：对采样点进行转换，得到采样点对应的多个能量值；根据多个能量值，确定采样点的最大能量值；根据最大能量值确定绕射点的位置信息；根据位置信息对绕射点的绕射波进行提取，得到提取结果；将多个采样点的提取结果进行叠加，得到分离后的绕射波。本发明专利技术通过利用绕射波走时曲线的双曲特征，采用改进的移动顶点拉东变换法，准确得到绕射点的位置信息，并对绕射波进行有效提取，有效压制噪音和反射波，无需速度信息便可较好的实现绕射波的保极性分离。

【技术实现步骤摘要】
绕射波分离方法、装置和电子设备
本专利技术涉及地震勘探
，尤其是涉及一种绕射波分离方法、装置和电子设备。
技术介绍
地震波在传播过程中，遇到断点、地层尖灭点、裂缝、陷落柱等地质异常体时，会以球面波的形式继续向前传播，形成绕射波，绕射波信息可用来实现地质异常体的高分辨率探测，在煤炭、石油等矿产资源的开采中具有重要的研究意义。但是绕射波能量较弱，常被反射波压制，难以获得理想的成像结果，采用平面波分解技术，处理精度难以准确把握，存在噪音干扰，难以对绕射点进行准确定位，利用绕射波的运动学和动力学特征进行处理，难以区分边缘绕射和点绕射。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种绕射波分离方法、装置和电子设备，以准确得到绕射点的位置信息，并对绕射波进行有效提取，有效压制噪音和反射波，较好的实现绕射波的保极性分离。第一方面，本专利技术实施例提供了一种绕射波分离方法，其中，包括：获取零偏移距道集；对零偏移距道集进行处理，得到地震波的瞬时相位；瞬时相位包括多个采样点；对每个采样点执行以下步骤：对采样点进行转换，得到采样点对应的多个能量值；根据多个能量值，确定采样点的最大能量值；根据最大能量值确定绕射点的位置信息；根据位置信息对绕射点的绕射波进行提取，得到提取结果；将多个采样点的提取结果进行叠加，得到分离后的绕射波。结合第一方面，本专利技术实施例提供了第一方面的第一种可能的实施方式，其中，对所述零偏移距道集进行处理，得到地震波的瞬时相位，包括：采用中值滤波法对零偏移距道集进行处理，得到处理后的零偏移距道集；对处理后的零偏移距道集进行希尔伯特变换，得到处理后的零偏移距道集的正交道集；根据处理后的零偏移距道集和正交道集，得到地震波的瞬时相位。结合第一方面，本专利技术实施例提供了第一方面的第二种可能的实施方式，其中，对采样点进行转换，得到采样点对应的多个能量值，包括：获取多个扫描曲率；根据多个扫描曲率，分别利用改进的移动顶点拉东变换法对采样点进行转换，得到每个扫描曲率对应的能量值。结合第一方面的第二种可能的实施方式，本专利技术实施例提供了第一方面的第三种可能的实施方式，其中，改进的移动顶点拉东变换法为：其中，x0表示采样点的坐标；t表示采样点的时间；q表示扫描曲率；τ表示转换后的采样点的时间；xi表示除采样点以外的剩余采样点的坐标；M表示采样点坐标之前的道数；N表示采样点坐标之后的道数。结合第一方面，本专利技术实施例提供了第一方面的第四种可能的实施方式，其中，位置信息包括绕射点的道号信息、时间信息和扫描曲率信息。结合第一方面的第四种可能的实施方式，本专利技术实施例提供了第一方面的第五种可能的实施方式，其中，根据位置信息对绕射点的绕射波进行提取，得到提取结果，包括：根据道号信息，利用第一计算公式进行双支拉东变换；第一计算公式为：其中，x0表示位置信息对应的采样点的坐标；t表示位置信息对应的采样点的时间；xi表示除采样点以外的剩余采样点的坐标；M表示采样点坐标之前的道数；N表示采样点坐标之后的道数；根据时间信息和扫描曲率信息，利用第二计算公式对(T，Q)处的区域进行拉东反变换，得到提取结果；第二计算公式为：其中，τ表示时间信息T；q表示扫描曲率信息Q；x0表示位置信息对应的采样点的坐标；M表示采样点坐标之前的道数；N表示采样点坐标之后的道数。第二方面，本专利技术实施例还提供一种绕射波分离装置，其中，包括：获取模块，用于获取零偏移距道集；处理模块，用于对零偏移距道集进行处理，得到地震波的瞬时相位；瞬时相位包括多个采样点；提取模块，用于对每个采样点执行以下步骤：对采样点进行转换，得到采样点对应的多个能量值；根据多个能量值，确定采样点的最大能量值；根据最大能量值确定绕射点的位置信息；根据位置信息对绕射点的绕射波进行提取，得到提取结果；叠加模块，用于将多个采样点的提取结果进行叠加，得到分离后的绕射波。结合第二方面，本专利技术实施例提供了第二方面的第一种可能的实施方式，其中，处理模块，包括：滤波模块，用于采用中值滤波法对零偏移距道集进行处理，得到处理后的零偏移距道集；变换模块，用于对处理后的零偏移距道集进行希尔伯特变换，得到处理后的零偏移距道集的正交道集；计算模块，用于根据处理后的零偏移距道集和正交道集，得到地震波的瞬时相位。第三方面，本专利技术实施例还提供一种电子设备，包括存储器、处理器，存储器中存储有可在处理器上运行的计算机程序，其中，处理器执行计算机程序时实现上述第一方面所述的方法的步骤。第四方面，本专利技术实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，其中，程序代码使处理器执行上述第一方面所述方法。本专利技术实施例带来了以下有益效果：本专利技术提供了一种绕射波分离方法、装置和电子设备，其中，该方法包括：获取零偏移距道集；对零偏移距道集进行处理，得到地震波的瞬时相位；瞬时相位包括多个采样点；对每个采样点执行以下步骤：对采样点进行转换，得到采样点对应的多个能量值；根据多个能量值，确定采样点的最大能量值；根据最大能量值确定绕射点的位置信息；根据位置信息对绕射点的绕射波进行提取，得到提取结果；将多个采样点的提取结果进行叠加，得到分离后的绕射波。本专利技术通过利用绕射波走时曲线的双曲特征，采用改进的移动顶点拉东变换法，准确得到绕射点的位置信息，并对绕射波进行有效提取，有效压制噪音和反射波，无需速度信息便可较好的实现绕射波的保极性分离。本专利技术的其他特征和优点将在随后的说明书中阐述，或者，部分特征和优点可以从说明书推知或毫无疑义地确定，或者通过实施本专利技术的上述技术即可得知。为使本专利技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施方式，并配合所附附图，作详细说明如下。附图说明为了更清楚地说明本专利技术具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例提供的一种绕射波分离方法的流程图；图2为本专利技术实施例提供的一种对每个采样点进行处理的流程图；图3为本专利技术实施例提供的另一种绕射波分离方法的流程图；图4为本专利技术实施例提供的一种零偏移距道集的结构示意图；图5为本专利技术实施例提供的一种瞬时相位的多个采样点对应的最大能量值的结构示意图；图6为本专利技术实施例提供的一种分离后的绕射波的结构示意图；图7为本专利技术实施例提供的一种绕射波分离装置的结构示意图；图8为本专利技术实施例提供的一种电子设备的结构示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。目前，利用平面波分解技术，分离得到的绕射波存在噪音干扰，而且并未准确指出绕射点的位置，利用绕射波的运动学和动力学特征，难以实现极性反转情况下的绕射波提取，且无法区分边缘绕射和点绕射，基于此，本专利技术实施例提供的一种绕射波分离方法、装置和电子设备，可以应用于极性反转情况下的绕射波提取。为便于本文档来自技高网...

【技术保护点】
1.一种绕射波分离方法，其特征在于，包括：获取零偏移距道集；对所述零偏移距道集进行处理，得到地震波的瞬时相位；所述瞬时相位包括多个采样点；对每个所述采样点执行以下步骤：对所述采样点进行转换，得到所述采样点对应的多个能量值；根据所述多个能量值，确定所述采样点的最大能量值；根据所述最大能量值确定绕射点的位置信息；根据所述位置信息对所述绕射点的绕射波进行提取，得到提取结果；将所述多个采样点的提取结果进行叠加，得到分离后的绕射波。

【技术特征摘要】
1.一种绕射波分离方法，其特征在于，包括：获取零偏移距道集；对所述零偏移距道集进行处理，得到地震波的瞬时相位；所述瞬时相位包括多个采样点；对每个所述采样点执行以下步骤：对所述采样点进行转换，得到所述采样点对应的多个能量值；根据所述多个能量值，确定所述采样点的最大能量值；根据所述最大能量值确定绕射点的位置信息；根据所述位置信息对所述绕射点的绕射波进行提取，得到提取结果；将所述多个采样点的提取结果进行叠加，得到分离后的绕射波。2.根据权利要求1所述的方法，其特征在于，所述对所述零偏移距道集进行处理，得到地震波的瞬时相位，包括：采用中值滤波法对所述零偏移距道集进行处理，得到处理后的零偏移距道集；对所述处理后的零偏移距道集进行希尔伯特变换，得到所述处理后的零偏移距道集的正交道集；根据所述处理后的零偏移距道集和所述正交道集，得到所述地震波的瞬时相位。3.根据权利要求1所述的方法，其特征在于，所述对所述采样点进行转换，得到所述采样点对应的多个能量值，包括：获取多个扫描曲率；根据多个所述扫描曲率，分别利用改进的移动顶点拉东变换法对所述采样点进行转换，得到每个所述扫描曲率对应的所述能量值。4.根据权利要求3所述的方法，其特征在于，所述改进的移动顶点拉东变换法为：其中，x0表示所述采样点的坐标；t表示所述采样点的时间；q表示所述扫描曲率；τ表示转换后的所述采样点的时间；xi表示除所述采样点以外的剩余采样点的坐标；M表示所述采样点坐标之前的道数；N表示所述采样点坐标之后的道数。5.根据权利要求1所述的方法，其特征在于，所述位置信息包括所述绕射点的道号信息、时间信息和扫描曲率信息。6.根据权利要求5所述的方法，其特征在于，所述根据所述位置信息对所述绕射点的绕射波进行提取，得到提取结果，包括：根据所述道号信息，利用第一计算公式进行双支拉东变换；所述第一计算公式为：其中，x0表示所述位...

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