一种偏移距道集的获取方法及装置

摘要

本发明实施例提供一种偏移距道集的获取方法及装置。方法包括：对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；汇总每个单炮对应的地表地表偏移距道集，获取所述待测区的偏移距道集。本发明实施例提供的方法及装置，应用GPU计算增加30％计算量即可得到待测区的偏移距道集，节省了计算资源，从而使得计算资源利用率更高。

权利要求

1.一种偏移距道集的获取方法，其特征在于，包括：

对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；

根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；

汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集；

所述波动方程偏移算法为逆时偏移算法；相应地，

基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果，包括：获取每个单炮在各时刻的检波点波场和在各时刻的炮点波场，并作为每个单炮的单炮数据；

基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果；

基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果I(x,s)的计算公式为：其中，x代表成像位置，s代表单炮，t代表时刻，n为正整数，u(x,s,t)为单炮s在t时刻的检波点波场，s(x,s,t)为单炮s在t时刻的炮点波场；

基于波动方程偏移算法，获取每个单炮的调制单炮数据对应的第二成像结果，包括：将所述第一成像结果I(x,s)的计算公式进行变形，得到变形后的第一成像结果I(x,s)的计算公式如下：其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场；

根据偏移距属性，对每个单炮在各时刻的单道检波点波场进行调制，得到调制后的单道检波点波场；

将每个单炮在各时刻的调制后的单道检波点波场和在各时刻的炮点波场，作为每个单炮的调制单炮数据；

基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果；

基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果Ih(x,s)的计算公式为：其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场，offset为单炮s到检波点r的距离。

2.根据权利要求1所述的偏移距道集的获取方法，其特征在于，根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集I(x,h|s)的计算公式为：I(x,h|s)＝I(x,s)δ(h-hA(x,s))其中，

其中，x代表成像位置，s代表单炮，I(x,s)为第一成像结果，Ih(x,s)为第二成像结果，h代表偏移距，hA(x,s)为单炮s的第二成像结果相对于第一成像结果的各成像点对应的偏移距。

3.根据权利要求2所述的偏移距道集的获取方法，其特征在于，汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集的计算公式为：其中，x代表成像位置，s代表单炮，I(x,h|s)为单炮s对应的地表偏移距道集。

说明书

一种偏移距道集的获取方法及装置

技术领域

[0001] 本发明涉及地震勘探技术领域，尤其涉及一种偏移距道集的获取方法及装置。

背景技术

[0002] 勘探地震数据深度成像方法包括Kirchhoff偏移，单程波偏移和逆时偏移等，单程波偏移和逆时偏移都属于波动方程偏移。在成像速度一定时，Kirchhoff偏移方法计算精度最低，单程波和逆时偏移方法成像效果更好。Kirchhoff偏移计算可以输出偏移距道集，偏移距道集可以用来进行AVO分析，速度更新等，是深度偏移方法必备输出数据。但当构造复杂区成像时，需应用精度更高的波动方程偏移，但因为其输入数据为单炮记录，不能输出偏移距道集，这将限制了该类数据深度的信息发掘和诸如速度更新等。目前计算波动方程偏移距道集的方法主要是将单炮数据均分成不同的偏移距段进行分步偏移，这相当于把一组数据分成若干组再计算，高强度的计算限制了其应用。

发明内容

[0003] 针对现有技术存在的问题，本发明实施例提供一种偏移距道集的获取方法及装置。

[0004] 第一方面，本发明实施例提供一种偏移距道集的获取方法，包括：

[0005] 对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；

[0006] 根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；

[0007] 汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0008] 进一步地，所述波动方程偏移算法为逆时偏移算法；相应地，

[0009] 基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果，包括：

[0010] 获取每个单炮在各时刻的检波点波场和在各时刻的炮点波场，并作为每个单炮的单炮数据；

[0011] 基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果。

[0012] 进一步地，基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果I(x,s)的计算公式为：

[0013]

[0014] 其中，x代表成像位置，s代表单炮，t代表时刻，n为正整数，u(x,s,t)为单炮s在t时刻的检波点波场，s(x,s,t)为单炮s在t时刻的炮点波场。

[0015] 进一步地，基于波动方程偏移算法，获取每个单炮的调制单炮数据对应的第二成像结果，包括：

[0016] 将所述第一成像结果I(x,s)的计算公式进行变形，得到变形后的第一成像结果I(x,s)的计算公式如下：

[0017]

[0018] 其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场；

[0019] 根据偏移距属性，对每个单炮在各时刻的单道检波点波场进行调制，得到调制后的单道检波点波场；

[0020] 将每个单炮在各时刻的调制后的单道检波点波场和在各时刻的炮点波场，作为每个单炮的调制单炮数据；

[0021] 基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果。

[0022] 进一步地，基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果Ih(x,s)的计算公式为：

[0023]

[0024] 其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场，offset为单炮s到检波点r的距离。

[0025] 进一步地，根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集I(x,h|s)的计算公式为：

[0026] I(x,h|s)＝I(x,s)δ(h-hA(x,s))

[0027] 其中，

[0028]

[0029] 其中，x代表成像位置，s代表单炮，I(x,s)为第一成像结果，Ih(x,s)为第二成像结果，h代表偏移距，hA(x,s)为单炮s的第二成像结果相对于第一成像结果的各成像点对应的偏移距。

[0030] 进一步地，汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集的计算公式为：

[0031]

[0032] 其中，x代表成像位置，s代表单炮，I(x,h|s)为单炮s对应的地表偏移距道集。

[0033] 第二方面，本发明实施例提供一种偏移距道集的获取装置，包括：

[0034] 成像结果获取模块，用于对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；

[0035] 地表偏移距道集获取模块，用于根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；

[0036] 偏移距道集获取模块，用于汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集

[0037] 第三方面，本发明实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如第一方面所提供的方法的步骤。

[0038] 第四方面，本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如第一方面所提供的方法的步骤。

[0039] 本发明实施例提供的一种偏移距道集的获取方法及装置，对于待测区内的每个单炮，根据偏移距属性，对单炮数据进行调制，并基于波动方程偏移算法，根据调制前的单炮数据和调制后的单炮数据，得到调制前的单炮数据对应的第一成像结果和调制后的单炮数据对应的第二成像结果，然后基于两个成像结果获取单炮对应的地表偏移距道集，最后根据每个单炮对应的地表偏移距道集，获取待测区的偏移距道集。本发明实施例提供的方法，无需将单炮数据均分成若干组分别利用不同的GPU进行计算，仅使用一个GPU计算即可得到待测区的偏移距道集，节省了计算资源，从而使得计算资源利用率更高。

附图说明

[0040] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0041] 图1为本发明实施例提供的一种偏移距道集的获取方法流程图；

[0042] 图2为本发明实施例提供的一种基于逆时偏移的地表偏移距道集的计算流程图；

[0043] 图3为本发明实施例提供的一种基于不同偏移方法获得的偏移Marmousi道集示意图；

[0044] 图4为本发明实施例提供的一种基于不同偏移方法获得的偏移距道集的示意图；

[0045] 图5为本发明实施例提供的一种偏移距道集的获取装置的结构示意图；

[0046] 图6为本发明实施例提供的一种电子设备的实体结构示意图。

具体实施方式

[0047] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0048] 图1为本发明实施例提供的一种偏移距道集的获取方法流程图，该方法的执行主体为偏移距道集的获取装置，例如服务器，本发明实施例对此不作具体限定，如图1所示，该方法包括：

[0049] 步骤101，对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的。

[0050] 步骤102，根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集。

[0051] 步骤103，汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0052] 具体地，首先需要说明的是，地震勘探指的是人工激发所引起的弹性波利用地下介质弹性和密度的差异，通过观测和分析人工地震产生的地震波在地下的传播规律，推断地下岩层的性质和形态的地球物理勘探方法。

[0053] 将需要进行地震勘探的地表范围称为待测区，为了获取待测区的地震勘探数据，需要在待测区中布置多个单炮(炮点)，并设置多个检波器(检波点)。需要说明的是，本发明实施例对单炮的设置以及检波器的设置不作具体限定。

[0054] 对于待测区中的每个单炮，将单炮激发的弹性波的波场称为炮点波场，将每个检波器检测到的波场称为检波点波场，将单炮在各时刻的检波点波场和在各时刻的炮点波场称为单炮的单炮数据。

[0055] 基于波动方程算法，获取每个单炮的单炮数据对应的第一成像结果。其中，波动方程算法可以为单程波偏移算法或逆时偏移算法等，本发明实施例对此不作限定。

[0056] 对于待测区中的每个单炮，通过偏移距属性对单炮数据进行调制得到调制单炮数据，其中，偏移距属性指的是单炮到各检波器的距离值，通过其对单炮数据中的检波点波场进行调制，将单炮在各时刻的调制后的检波点波场和在各时刻的炮点波场称为单炮的调制单炮数据。

[0057] 基于波动方程算法，获取每个单炮的调制单炮数据对应的第二成像结果。其中，波动方程算法可以为单程波偏移算法或逆时偏移算法等，本发明实施例对此不作限定。需要说明的是，获取第一成像结果和获取第二成像结果的波动方程算法一致。

[0058] 根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集。汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0059] 本发明实施例提供的方法，对于待测区内的每个单炮，根据偏移距属性，对单炮数据进行调制，并基于波动方程偏移算法，根据调制前的单炮数据和调制后的单炮数据，得到调制前的单炮数据对应的第一成像结果和调制后的单炮数据对应的第二成像结果，然后基于两个成像结果获取单炮对应的地表偏移距道集，最后根据每个单炮对应的地表偏移距道集，获取待测区的偏移距道集。本发明实施例提供的方法，使用GPU计算地表偏移距道集，计算量只增加30％，节省了计算资源，从而使得计算资源利用率更高。

[0060] 基于上述任一实施例，本发明实施例中的所述波动方程偏移算法为逆时偏移算法；相应地，

[0061] 基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果，包括：

[0062] 获取每个单炮在各时刻的检波点波场和在各时刻的炮点波场，并作为每个单炮的单炮数据；

[0063] 基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果。

[0064] 其中，基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果I(x,s)的计算公式为：

[0065]

[0066] 其中，x代表成像位置，s代表单炮，t代表时刻，n为正整数，u(x,s,t)为单炮s在t时刻的检波点波场，s(x,s,t)为单炮s在t时刻的炮点波场。

[0067] 具体地，本发明实施例主要说明如何根据调制前的单炮数据(单炮数据)计算单炮对应的第一成像结果。

[0068] 在本发明实施例中，将单炮的检波点波场和炮点波场称为调制前的单炮数据(单炮数据)，将每个单炮的检波点波场和炮点波场相乘，得到相乘后的结果，并将所有单炮的相乘后的结果相加，得到第一成像结果。其中，检波点波场和炮点波场均通过波场传播计算获得。

[0069] 基于上述任一实施例，基于波动方程偏移算法，获取每个单炮的调制单炮数据对应的第二成像结果，包括：

[0070] 将所述第一成像结果I(x,s)的计算公式进行变形，得到变形后的第一成像结果I(x,s)的计算公式如下：

[0071]

[0072] 其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场；

[0073] 根据偏移距属性，对每个单炮在各时刻的单道检波点波场进行调制，得到调制后的单道检波点波场；

[0074] 将每个单炮在各时刻的调制后的单道检波点波场和在各时刻的炮点波场，作为每个单炮的调制单炮数据；

[0075] 基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果。

[0076] 其中，基于所述逆时偏移算法，根据每个单炮的调制单炮数据，获取对应的第二成像结果Ih(x,s)的计算公式为：

[0077]

[0078] 其中，x代表成像位置，s代表单炮，t代表时刻，r代表检波点，n为正整数，s(x,s,t)为单炮s在t时刻的炮点波场，u(x,s,r,t)为单炮s在t时刻的单道检波点波场，offset为单炮s到检波点r的距离。

[0079] 具体地，本发明实施例主要说明如何根据调制后的单炮数据(调制单炮数据)计算单炮对应的第二成像结果。

[0080] 在本发明实施例中，首先将上述实施例中第一成像结果的计算公式进行变形，即，将上述公式改写成单道计算形式的公式，如下：

[0081]

[0082] 然后根据单炮到各检波点的距离对检波点的波场进行单道偏移距的调制，得到第二成像结果Ih(x,s)的计算公式为：

[0083]

[0084] 基于上述任一实施例，根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集I(x,h|s)的计算公式为：

[0085] I(x,h|s)＝I(x,s)δ(h-hA(x,s))

[0086] 其中，

[0087]

[0088] 其中，x代表成像位置，s代表单炮，I(x,s)为第一成像结果，Ih(x,s)为第二成像结果，h代表偏移距，hA(x,s)为单炮s的第二成像结果相对于第一成像结果的各成像点对应的偏移距。

[0089] 其中，汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集的计算公式为：

[0090]

[0091] 其中，x代表成像位置，s代表单炮，I(x,h|s)为单炮s对应的地表偏移距道集。

[0092] 具体地，本发明实施例主要说明如何根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集，以及如何根据每个单炮对应的地表偏移距道集获取待测区的偏移距道集。

[0093] 作为一个优选的实施例，本发明实施例结合图2，对上述实施例中的方法进行进一步具体说明，图2为本发明实施例提供的一种基于逆时偏移的地表偏移距道集的计算流程图，如图2所示，计算流程如下：

[0094] 1、CPU中，偏移参数准备；

[0095] 2、CPU中，偏移速度读取；

[0096] 3、CPU中，随机边界条件设置；

[0097] 4、将CPU中的上述数据拷贝至GPU中；

[0098] 5、GPU中，令it＝0，根据接收到的数据，进行炮点正传；

[0099] 6、GPU中，插入炮点波场数据；

[0100] 7、GPU中，令it＝it+1；

[0101] 8、GPU中，判断：it是否小于itmax；

[0102] 9、GPU中，若是，则进行炮点波场反传；

[0103] 10、GPU中，原检波波场传播，插入检波点波场；调制检波场传播，插入检波点波场；

[0104] 11、GPU中，提取偏移距；

[0105] 12、令it＝it-1；

[0106] 13、判断it是否大于或等于0；

[0107] 14，若否，则重新执行步骤8；若否，则拷贝数据到CPU；

[0108] 15、CPU中，根据拷贝的数据，提取偏移距道集。

[0109] 需要说明的是，图2中，如果去掉调制检波波场传播部分，就是逆时偏移计算流程，因此从图2中可以看出，主要包含炮点波场正、反传，检波点反传，调制检波点反传，因为只增加了调制检波点反传，因此在炮点和检波点孔径一致的情况下，计算量只增加了1/3，因此该方法是一种高效的计算方法，节省了计算资源，从而使得计算资源利用率更高。

[0110] 图3为本发明实施例提供的一种基于不同偏移方法获得的偏移Marmousi道集示意图，其中，(a)部分为利用本发明实施例中提供的方法获得的偏移Marmousi道集示意，(b)部分为利用Kirchhoff偏移方法获得的偏移Marmousi道集示意，显然，本发明实施例中提供的方法相对于Kirchhoff偏移成像效果更好。图4为本发明实施例提供的一种基于不同偏移方法获得的偏移距道集的示意图，其中，(a)部分为利用本发明实施例中提供的方法获得的偏移距道集的示意，(b)部分为利用Kirchhoff偏移方法获得的偏移距道集的示意，显然，本发明实施例中提供的方法相对于Kirchhoff偏移产生的道集更平滑，分辨率更高。

[0111] 基于上述任一实施例，图5为本发明实施例提供的一种偏移距道集的获取装置的结构示意图，如图5所示，该装置包括：

[0112] 成像结果获取模块501，用于对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；地表偏移距道集获取模块502，用于根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；偏移距道集获取模块503，用于汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0113] 具体地，本发明实施例提供的装置具体用于执行上述方法实施例，本发明实施例对此不再进行赘述。本发明实施例提供的装置，对于待测区内的每个单炮，根据偏移距属性，对单炮数据进行调制，并基于波动方程偏移算法，根据调制前的单炮数据和调制后的单炮数据，得到调制前的单炮数据对应的第一成像结果和调制后的单炮数据对应的第二成像结果，然后基于两个成像结果获取单炮对应的地表偏移距道集，最后根据每个单炮对应的地表偏移距道集，获取待测区的偏移距道集。本发明实施例提供的方法，无需将单炮数据均分成若干组分别利用不同的GPU进行计算，仅使用一个GPU计算即可得到待测区的偏移距道集，节省了计算资源，从而使得计算资源利用率更高。

[0114] 图6为本发明实施例提供的一种电子设备的实体结构示意图，如图6所示，该电子设备可以包括：处理器(processor)601、通信接口(Communications Interface)602、存储器(memory)603和通信总线604，其中，处理器601，通信接口602，存储器603通过通信总线604完成相互间的通信。处理器601可以调用存储在存储器603上并可在处理器601上运行的计算机程序，以执行上述各实施例提供的方法，例如包括：对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0115] 此外，上述的存储器603中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

[0116] 本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的传输方法，例如包括：对于待测区内多个单炮中的每个单炮，基于波动方程偏移算法，获取每个单炮的单炮数据对应的第一成像结果以及每个单炮的调制单炮数据对应的第二成像结果；其中，所述调制单炮数据是根据偏移距属性对所述单炮数据进行调制后得到的；根据每个单炮对应的第一成像结果和第二成像结果，获取每个单炮对应的地表偏移距道集；汇总每个单炮对应的地表偏移距道集，获取所述待测区的偏移距道集。

[0117] 以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

[0118] 通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。

[0119] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。