波的走时获得方法及装置、成像方法及装置制造方法及图纸

本发明专利技术实施例公开了一种波的走时获得方法，包括：获得第一点至第二点的水平距离；将水平距离输入第一预置模型，得到波在第一点和第二点之间的射线参数；将该射线参数输入第二预置模型，得到波从第一点传播至第二点的走时；其中，第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，和/或，第二预置模型是波的传播时间多项式经至少一次Shanks变换得到。由于第一预置模型和第二预置模型中任意一个是根据上述对应的多项式经至少一次Shanks变换得到的，经Shanks变换得到的第一预置模型和/或第二预置模型具有更高的精度和收敛速度，结果较根据原始计算公式得到的结果精度更高、计算耗时更短，能够满足地质勘探的要求。

【技术实现步骤摘要】
波的走时获得方法及装置、成像方法及装置本申请要求于2017年10月31日提交中国专利局、申请号为201711047776.9、专利技术名称为“波的走时获得方法及装置、成像方法及装置”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。
本专利技术涉及地质勘探
，尤其涉及一种波的走时获得方法及装置、成像方法及装置。
技术介绍
波的走时，指的是波从震源传播至观测点所需的时间。例如地质探勘领域所涉及的地震波测绘，地震波的走时指的是地震波从炮点(即震源点)传播至成像点(即观测点)所需的时间。在实际生产中，一般利用地面传感器记录地震波的走时、波形及能量，并根据上述记录对地下结构进行反演，推测出地下地质体的真实位置及地下岩石的物理性质。根据反演出的地下结构，可以对地下的含油气性进行推测，为油气资源的勘探提供数据支持。但是，现有的走时计算方法一般是对走时计算多项式的结果进行预估，存在较大的误差，在一些对走时精度要求高的场景中，如采用绕射进行隔夹层描述以及深层大偏移观测时，走时的误差会严重影响成像及速度分析的结果，使得结果的准确度及精度不能满足地质勘探的需求。
技术实现思路
有鉴于上述背景，本专利技术提供了一种地震波走时获得方法及装置、地下结果测绘方法及装置，能够解决现有技术中走时计算误差而造成的地质观测结果不能满足勘探需求的问题。本专利技术实施例提供的一种波的走时获得方法，包括：获得第一点至第二点的水平距离；将所述水平距离输入第一预置模型，得到波在所述第一点和所述第二点之间的射线参数；将所述射线参数输入第二预置模型，得到波从所述第一点传播至所述第二点的走时；其中，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型为波的传播时间多项式；或者，所述第一预置模型为波的传播距离多项式，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到；或者，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到。可选的，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，具体为：当所述波的传播距离多项式包括至少四项时，所述波的传播距离多项式经两次Shanks变换得到所述第一预置模型；所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到，具体为：当所述波的传播时间多项式包括至少四项时，所述波的传播时间多项式经两次Shanks变换得到所述第二预置模型。可选的，所述波的传播距离多项式具体是根据波的传播距离公式经多项式展开后得到；所述波的传播时间多项式具体是根据波的传播时间公式经多项式展开后得到；其中，x为波传播的水平距离，t为波的走时，ai和bj均为系数，t0为时间校正参数，m和n均为大于2的整数，p为波的射线参数，v为波的速度。可选的，当m＝3时，所述波的传播距离多项式经一次Shanks变换得到的第一阈值模型，具体为：当n＝4时，所述波的传播时间多项式经一次Shanks变换得到的第二阈值模型，具体为：本专利技术实施例提供的一种成像方法，包括：遍历成像区域中所有成像点，利用上述上述实施例所述的波的走时获得方法，获得每个成像点对应的第一走时和第二走时；其中，第一走时为波从震源点传播至该成像点的走时，第二走时为所述波的反射波从该成像点传播至检波点的走时；基于所有成像点的对应的第一走时和第二走时，反演所述成像区域的地下结构，得到所述成像区域的地下结构图。本专利技术实施例提供的一种波的走时获得装置，包括：获取模块、第一计算模块和第二计算模块；所述获取模块，用于获得第一点至第二点的水平距离；所述第一计算模块，用于将所述获取模块获得的水平距离输入第一预置模型，得到波在所述第一点和所述第二点之间的射线参数；所述第二计算模块，用于将所述第一计算模块获得的射线参数输入第二预置模型，得到波从所述第一点传播至所述第二点的走时；其中，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型为波的传播时间多项式；或者，所述第一预置模型为波的传播距离多项式，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到；或者，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到。本专利技术实施例提供的一种成像装置，包括：走时计算模块和成像模块；所述走时计算模块，用于遍历成像区域中所有成像点，利用上述实施例所述的波的走时获得方法，获得每个成像点对应的第一走时和第二走时；其中，第一走时为波从震源点传播至该成像点的走时，第二走时为所述波的反射波从该成像点传播至检波点的走时；所述成像模块，用于基于所述走时计算模块获得的成像点对应的第一走时和第二走时，反演所述成像区域的地下结构，得到所述成像区域的地下结构图。本专利技术实施例还提供了一种计算机可读存储介质，所述机算机可读存储介质中存储有指令，当所述指令在终端设备上运行时，使得所述终端设备执行上述实施例所述的波的走时获得方法。本专利技术实施例还提供了一种终端设备，包括：存储器和处理器；所述存储器，用于存储程序代码；所述处理器，用于获得所述程序代码执行上述实施例所述的波的走时获得方法。与现有技术相比，本专利技术至少具有以下优点：在本专利技术实施例中，在获得第一点和第二点的水平距离后，将该水平距离输入第一预置模型得到波在第一点和第二点之间的射线参数，再将该射线参数输入第二预置模型，从而得到波从第一点传播至第二点的走时。由于第一预置模型和第二预置模型中任意一个是根据原始计算公式经至少一次Shanks变换得到的，经Shanks变换得到的第一预置模型和/或第二预置模型较原始计算公式阶数低，具有更高的精度和收敛速度，使得得到的射线参数和/或走时较根据原始计算公式得到的结果精度更高、计算耗时更短，使得走时计算能够满足地质勘探所需的高精度、快速的要求。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。图1为本专利技术实施例提供的一种波的走时获得方法的流程示意图；图2为本专利技术实施例提供的一种成像方法的流程示意图；图3a为现有走时计算方法得到的地下结构示意图；图3b为本专利技术具体实施例提供的一种成像方法得到的地下结构示意图；图4为本专利技术实施例提供的一种波的走时获得装置的结构示意图；图5为本专利技术实施例提供的一种成像装置的结构示意图。具体实施方式为了使本
的人员更好地理解本专利技术方案，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。首先需要说明的是，本专利技术实施例提供的波的走时获得方法及装置，不仅可以用于获得地震波的走时，还可以用于获得其他形式的波的走时，例如光波、雷达波、声波等。为了便于理解本文档来自技高网...

【技术保护点】
1.一种波的走时获得方法，其特征在于，所述方法，包括：获得第一点至第二点的水平距离；将所述水平距离输入第一预置模型，得到波在所述第一点和所述第二点之间的射线参数；将所述射线参数输入第二预置模型，得到波从所述第一点传播至所述第二点的走时；其中，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型为波的传播时间多项式；或者，所述第一预置模型为波的传播距离多项式，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到；或者，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到。

【技术特征摘要】
2017.10.31 CN 20171104777691.一种波的走时获得方法，其特征在于，所述方法，包括：获得第一点至第二点的水平距离；将所述水平距离输入第一预置模型，得到波在所述第一点和所述第二点之间的射线参数；将所述射线参数输入第二预置模型，得到波从所述第一点传播至所述第二点的走时；其中，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型为波的传播时间多项式；或者，所述第一预置模型为波的传播距离多项式，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到；或者，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到。2.根据权利要求1所述的方法，其特征在于，所述第一预置模型是波的传播距离多项式经至少一次Shanks变换得到，具体为：当所述波的传播距离多项式包括至少四项时，所述波的传播距离多项式经两次Shanks变换得到所述第一预置模型；所述第二预置模型是波的传播时间多项式经至少一次Shanks变换得到，具体为：当所述波的传播时间多项式包括至少四项时，所述波的传播时间多项式经两次Shanks变换得到所述第二预置模型。3.根据权利要求1所述的方法，其特征在于，所述波的传播距离多项式具体是根据波的传播距离公式经多项式展开后得到；所述波的传播时间多项式具体是根据波的传播时间公式经多项式展开后得到；其中，x为波传播的水平距离，t为波的走时，ai和bj均为系数，t0为时间校正参数，m和n均为大于2的整数，p为波的射线参数，v为波的速度。4.根据权利要求3所述的方法，其特征在于，当m＝3时，所述波的传播距离多项式经一次Shanks变换得到的第一阈值模型，具体为：当n＝4时，所述波的传播时间多项式经一次Shanks变换得到的第二阈值模型，具体为：5.一种成像方法，其特征在于，所述方法，包括：遍历成像区域中所有成像点，利用上述权利要求1-4任意一项所述的波的走时获得方...

【专利技术属性】
技术研发人员：王真理，
申请(专利权)人：中国科学院地质与地球物理研究所，
类型：发明
国别省市：北京,11