[发明专利]一种偏移距道集的获取方法及装置

说明书

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

技术领域

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

背景技术

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

发明内容

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

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

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

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

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

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

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

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

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

进一步地，基于所述逆时偏移算法，根据每个单炮的单炮数据，获取对应的第一成像结果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时刻的单道检波点波场；

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

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