真振幅的偏移成像方法、装置、设备及存储介质与流程

1.本公开属于油气勘探技术领域，特别涉及一种真振幅的偏移成像方法、装置、设备及存储介质。


背景技术：

2.地震勘探方法是石油、天然气及其他能源勘探的重要方法，而地震数据真振幅偏移成像方法作为现代地震勘探数据处理的关键环节，为地震解释和岩性反演提供必要的技术支撑。
3.相关技术中，地震数据真振幅偏移成像方法中，首先获取地震数据，该地震数据包括地表的波场值，然后依据波动方程，利用地表的波场值成像地下构造。
4.然而，由于波动方程含有深度方向的二阶偏导数，数学上要使该波动方程能够求解，必须确定地表的波场值和地表波场在深度方向的偏导数。而地震数据在获取时，得不到地表波场在深度方向的偏导数。因此，只能对波动方程进行各种近似处理，以便通过地表波场值对地下构造进行成像。而对波动方程的近似处理会导致成像的失真，以致真振幅偏移图像不能真实的反映地下构造的岩性变化，不利于岩性勘探。


技术实现要素：

5.本公开实施例提供了一种真振幅的偏移成像方法、装置、设备及存储介质，能够大幅提高地震数据的真振幅偏移成像效率。所述技术方案如下：
6.本公开实施例提供了一种真振幅的偏移成像方法，所述偏移成像方法包括：
7.确定发射点在三维坐标系中的位置，所述发射点为地震波的震源所在的位置；
8.将模拟多次波以所述发射点为衍射点进行衍射，所述模拟多次波是基于实际多次波生成的，所述实际多次波是指所述地震波经过地层多次衍射所形成的波；获取所述模拟多次波在衍射过程中的衍射波形图；基于所述衍射波形图，得到地震数据的真振幅偏移图像。
9.在本公开的又一种实现方式中，所述获取所述模拟多次波在衍射过程中的衍射波形图，包括：
10.在所述模拟多次波以所述发射点为衍射点进行衍射的过程中，以相同时间间隔捕捉所述模拟多次波，得到多个衍射捕捉点；确定所述衍射捕捉点在所述三维坐标系中的位置；根据所述衍射捕捉点在所述三维坐标系中的位置，得到所述衍射捕捉点所形成的衍射捕捉图；根据所述衍射捕捉图，得到所述模拟多次波在所述三维坐标系中不同时间段所对应的分波形图，每个时间段对应不同视角的多个分波形图；根据所述分波形图，得到所述多次波在所述三维坐标系中的衍射波形图。
11.在本公开的又一种实现方式中，所述确定所述衍射捕捉点在所述三维坐标系中的位置，包括：
12.根据以下公式得到所述衍射捕捉点在所述三维坐标系中的z轴方向坐标值；
[0013][0014]
其中，fa为衍射捕捉点在所述三维坐标系中的z轴方向的坐标值，也为衍射捕捉点的振幅值；in为所述发射点在所述三维坐标系中的z轴方向的坐标值；e
xw
为扰动波场值；cav为衍射捕捉点捕捉时在所述三维坐标系中对应的z轴方向的坐标值，所述z轴方向用于表示垂直于地面的方向。
[0015]
在本公开的又一种实现方式中，所述分波形图包括x视角分波形图、y视角分波形图和z视角分波形图，所述x视角分波形图为沿着所述三维坐标系中的x轴方向观察到的衍射捕捉图，所述y视角分波形图为沿着所述三维坐标系中的y轴方向观察到的衍射捕捉图，所述z视角分波形图为沿着所述三维坐标系中的z轴方向观察到的衍射捕捉图；所述根据所述分波形图，得到所述多次波在所述三维坐标系中的衍射波形图，包括：
[0016]
将不同时间段所对应的所述分波形图以所述分波形图的起始位置所对应的坐标排列至所述三维坐标系的xoy平面中，得到所述x视角分波形图、所述y视角分波形图和所述z视角分波形图之间的交叉点的信息；根据所述交叉点的信息，得到不同时间段的合波形图的起始位置和终止位置，所述合波形图为所述x视角分波形图、所述y视角分波形图和所述z视角分波形图合在一起形成的合波形图；根据不同时间段的所述合波形图的起始位置和终止位置，得到所述多次波在所述三维坐标系中的衍射波形图。
[0017]
在本公开的又一种实现方式中，所述根据所述交叉点的信息，得到不同时间段的合波形图的起始位置和终止位置，包括：
[0018]
若所述交叉点的个数为零，则将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的起始位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图起始位置的坐标值，将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的终止位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图终止位置的坐标值；
[0019]
若所述交叉点的个数为一，则将所述交叉点对应的位置作为所述合波形图的起始位置和终止位置；
[0020]
若所述交叉点的个数为二，则将两个所述交叉点的位置分别作为所述合波形图的起始位置和终止位置；
[0021]
若所述交叉点的个数大于二，则选取所述交叉点的位置中对应的y轴方向坐标差最大的两个交叉点的位置，分别作为所述合波形图的起始位置和终止位置。
[0022]
在本公开的又一种实现方式中，所述基于所述衍射波形图，得到地震数据的真振幅偏移图像，包括：
[0023]
对所述衍射波形图进行保幅成像处理，得到所述真振幅偏移图像。
[0024]
在另一方面，还提供一种真振幅的偏移成像装置，所述偏移成像装置包括：发射点确定模块，用于确定发射点在三维坐标系中的位置，所述发射点为地震波的震源所在的位置；衍射模块，用于将模拟多次波以所述发射点为衍射点进行衍射，所述模拟多次波是基于实际多次波生成的，所述实际多次波是指所述地震波经过地层多次衍射所形成的波；衍射波形图获取模块，用于获取所述模拟多次波在衍射过程中的衍射波形图；成像模块，用于基于所述衍射波形图，得到地震数据的真振幅偏移图像。
[0025]
在本公开的又一种实现方式中，还提供一种计算机设备，所述计算机设备包括处
理器和被配置为存储所述处理器可执行指令的存储器；所述处理器被配置为执行以上方法。
[0026]
在本公开的又一种实现方式中，还提供一种计算机存储介质，其上存储有计算机指令，所述计算机指令被处理器执行时实现以上方法。
[0027]
本公开实施例提供的技术方案带来的有益效果是：
[0028]
通过本公开实施例提供的真振幅的偏移成像方法在对地震资料进行成像时，由于该成像模拟方法是通过对多次波以所述发射点为衍射点进行衍射得到多次波在三维坐标系中的衍射波形图，并根据衍射波形图形成地震数据的真振幅偏移图像，所以无需对波动方程进行近似求解，便可得到地震数据的真振幅偏移成像，大大提高了成像的准确性，使得真振幅偏移图像能够更真实的反映地下构造的岩性变化。
附图说明
[0029]
为了更清楚地说明本公开实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0030]
图1是本公开实施例提供的一种真振幅的偏移成像方法流程图；
[0031]
图2是本公开实施例提供的另一种真振幅的偏移成像方法流程图；
[0032]
图3是本公开实施例提供的一种真振幅的偏移成像装置的结构示意图；
[0033]
图4是本公开实施例提供的一种计算机设备的结构示意图。
具体实施方式
[0034]
为使本公开的目的、技术方案和优点更加清楚，下面将结合附图对本公开实施方式作进一步地详细描述。
[0035]
本公开实施例提供了一种真振幅的偏移成像模拟方法。该方法可以由计算机设备执行。如图1所示，该方法包括：
[0036]
s101：确定发射点在三维坐标系中的位置，发射点为地震波的震源所在的位置。
[0037]
地震波的震源为波点发射基站。
[0038]
s102：将模拟多次波以发射点为衍射点进行衍射，模拟多次波是基于实际多次波生成的，实际多次波是指地震波经过地层多次衍射所形成的波。
[0039]
在本公开实施例中，衍射也可以称为绕射。模拟多次波经过位于发射点的障碍物(障碍物的尺寸不大于模拟多次波的波长)后进行衍射。
[0040]
障碍物的尺寸可以自由设置，只要满足衍射条件即可。
[0041]
本实施例中，实际多次波是在地面建立波点发射基站，然后通过波点发射基站向地层发射传播地震波(例如夯锤模拟，夯锤锤击地面)，地震波在地层中多次反射而形成的。
[0042]
实际多次波的获取是通过对地震数据进行分离得到的。模拟多次波是根据实际多次波的传播过程，对实际多次波的传播特性进行模拟得到的。模拟多次波具有与实际多次波相同的物理特性，比如具有相同的振幅、相位、传播速度。本公开实施例对基于实际多次波生成模拟多次波的方式不做限定，可以采用任意计算机仿真算法和数学算法等实现。
[0043]
s103：获取模拟多次波在衍射过程中的衍射波形图。
[0044]
s104：基于衍射波形图，得到地震数据的真振幅偏移图像。
[0045]
通过本公开实施例提供的真振幅的偏移成像方法在对地震资料进行成像时，由于该成像模拟方法是通过对模拟多次波以发射点为衍射点进行衍射得到模拟多次波在三维坐标系中的衍射波形图，并根据衍射波形图形成地震数据的真振幅偏移图像，这样就可以用模拟多次波作为输入数据，通过形成衍射波形图，将地表实际多次波以衍射波形图的形式来逆向反映各个地层界面，即将实际多次波归位到地层的真实构造上，以此便可得到地震数据的真振幅偏移图像。以上方法无需对波动方程进行近似求解，便可得到地震数据的真振幅偏移成像，大大提高了成像的准确性，使得真振幅偏移图像能够更真实的反映地下构造的岩性变化。
[0046]
图2是本公开实施例提供的另一种真振幅的偏移成像方法流程图。该方法可以由计算机设备执行。结合图2，该方法包括：
[0047]
s201：确定发射点在三维坐标系中的位置。
[0048]
示例性地，三维坐标系是笛卡尔坐标系，包括x、y、z轴构成。其中z轴所在的方向为垂直于地面的方向，简称为z轴方向。
[0049]
本实施例中，发射点可以记为γ点，γ点在三维坐标系中的位置可以定位在三维坐标系中z轴方向坐标值为非零的位置。
[0050]
s202：将模拟多次波以发射点为衍射点进行衍射。
[0051]
其中，模拟多次波的相关内容参见s102，在此省略详细描述。
[0052]
s203：在模拟多次波以发射点为衍射点进行衍射的过程中，以相同时间间隔捕捉模拟多次波，得到多个衍射捕捉点在三维坐标系中的位置。
[0053]
示例性地，时间间隔可以根据实际需要设置，通常为毫秒级。例如，为1毫秒、2毫秒等。本实施例中，每间隔1毫秒对模拟多次波进行捕捉，可以得到多个衍射捕捉点在三维坐标系中的位置。
[0054]
s204：根据衍射捕捉点在三维坐标系中的位置，得到衍射捕捉点所形成的衍射捕捉图。
[0055]
在一些示例中，该步骤s204包括：
[0056]
2041：根据扰动波场对所述衍射捕捉点在所述三维坐标系中的位置进行修正。
[0057]
2041包括：
[0058]
(1)基于扰动波场值确定修正因子。
[0059]
该修正因子用于表示扰动波场值对模拟多次波在衍射过程中的影响程度。
[0060]
扰动波场值是指扰动波场中某个位置在某个时刻的振幅值。在本公开实施例中，该扰动波场值是指发射点在采样时刻的扰动波场值。该扰动波场值可以根据采样时刻该发射点对应环境中的波场值、声波慢度场和震源子波，计算得到。具体计算过程可以参见相关技术中的射线追踪法等手段，这里不再赘述。
[0061]
其中，采样时刻该发射点对应环境中的波场值、声波慢度场和震源子波可以通过测量得到。本实施例中，每隔设定时长(例如15-20min)重新确定环境中的波场值、声波慢度场和震源子波，从而计算得到不同时段对应的扰动波场值。
[0062]
本实施例中，该修正因子可以等于发射点的位置与扰动波场值的比例值。比例值
越大，扰动波场值对模拟多次波在衍射过程中的影响越大。
[0063]
修正因子可以按照以下公式计算：
[0064]
sc＝in÷exw
；
ꢀꢀꢀ
(1)
[0065]
公式(1)中：sc为比例值，in为发射点在三维坐标系中的z轴方向坐标值，e
xw
为扰动波场值。
[0066]
(2)根据修正因子，修正多个衍射捕捉点在三维坐标系中的位置。
[0067]
多个衍射捕捉点在三维坐标系中的位置通过以下公式确定：
[0068]
fa＝sc×
cav；
ꢀꢀꢀ
(2)
[0069]
公式(2)中，fa为衍射捕捉点修正后的位置在三维坐标系中的z轴方向的坐标值，也为衍射捕捉点的振幅值；cav为衍射捕捉点捕捉时的位置在三维坐标系中对应的z轴方向的坐标值。
[0070]
也即是，该步骤2041可以包括：
[0071]
根据以下公式(3)得到所述衍射捕捉点在所述三维坐标系中的z轴方向坐标值；
[0072][0073]
公式(3)中各参数的含义参见前述公式(1)和(2)。
[0074]
衍射捕捉点经过修正后的位置在三维坐标系中的x轴方向的坐标值以及y轴方向的坐标值认为修正之前对应的坐标值。
[0075]
2042：根据修正后的衍射捕捉点在三维坐标系中的位置，得到衍射捕捉点所形成的衍射捕捉图。
[0076]
可选地，2042包括：
[0077]
将单位时间长度内的多个衍射捕捉点的修正后的位置连成一条线，作为这多个衍射捕捉点对应的衍射捕捉图。
[0078]
示例性地，单位时间长度为前述时间间隔的整数倍，可以根据实际需要设置。例如，一个单位时间长度包括1000个时间间隔。当时间间隔为1毫秒时，单位时间长度为1秒。
[0079]
在另一些示例中，直接根据衍射捕捉点在三维坐标系中的位置，得到衍射捕捉点所形成的衍射捕捉图。也即是，将单位时间长度内的多个衍射捕捉点的位置连成一条线，作为这多个衍射捕捉点对应的衍射捕捉图。
[0080]
s205：根据衍射捕捉图，得到模拟多次波在三维坐标系中不同时间段所对应的分波形图。
[0081]
每个时间段对应的分波形图均包括x视角分波形图、y视角分波形图和z视角分波形图，x视角分波形图为沿着三维坐标系中的x轴方向观察到的衍射捕捉图，y视角分波形图为沿着三维坐标系中的y轴方向观察到的衍射捕捉图，z视角分波形图为沿着三维坐标系中的z轴方向观察到的衍射捕捉图。
[0082]
衍射捕捉图在三维坐标系中为三维立体图，所以，沿着三维坐标系中的不同坐标轴可以得到对应的x视角分波形图、y视角分波形图和z视角分波形图。
[0083]
这里，一个时间段包括连续的多个单位时间长度，不同的时间段互不重合。在本公开实施例中，一个时间段所包含的单位时间长度的数量可以根据实际需要设置，例如，10个或者20个单位时间长度等。例如，当单位时间长度为1秒时，一个时间段包括10个单位时间
长度，为10秒。
[0084]
本实施例中，分波形图的获取为每间隔10秒进行。即首先通过获取连续10秒内的衍射捕捉图，再得到与该连续10秒内的衍射捕捉图对应的三个分波形图。
[0085]
s206：根据分波形图，得到模拟多次波在三维坐标系中的衍射波形图。
[0086]
可选地，s206包括：
[0087]
2061：将不同时间段所对应的分波形图以分波形图的起始位置所对应的坐标排列至三维坐标系的xoy平面中，得到x视角分波形图、y视角分波形图和z视角分波形图之间的交叉点的信息。这里，交叉点的信息至少包括交叉点的数量和位置。
[0088]
2062：根据交叉点的信息，确定不同时间段的合波形图的起始位置和终止位置，合波形图为x视角分波形图、y视角分波形图和z视角分波形图合在一起形成的合波形图。
[0089]
比如，当得到第1个十秒内记录的衍射捕捉图，将该衍射捕捉图中的x视角分波形图、y视角分波形图和z视角分波形图分别按照对应的起始位置的坐标排列至三维坐标系的xoy平面中，这样便可得到该三个分波形图之间的交叉点的个数。
[0090]
比如，获取到第1个十秒内记录的衍射捕捉图后，分别得到对应的x视角分波形图、y视角分波形图和z视角分波形图，然后将x视角分波形图、y视角分波形图和z视角分波形图对应到xoy平面内，此时，x视角分波形图的各点位置的坐标为由(0、y、z)变为xoy平面内(x、y)，即y值变为x，z值变为y。y视角分波形图的各点位置的坐标为由(x、0、z)变为xoy平面内(x、y)，即x值变为x，z值变为y。z视角分波形图的各点位置的坐标(x、y、0)直接取为(x、y)。这样便可将不同时间段所对应的分波形图沿x轴方向进行排列。本实施例中，通过以下方式获得合波形图的起始位置和终止位置：
[0091]
(1)若交叉点的个数为零，则将x视角分波形图、y视角分波形图和z视角分波形图的起始位置对应的不同坐标轴的坐标值的总和的平均值作为合波形图起始位置的坐标值，将x视角分波形图、y视角分波形图和z视角分波形图的终止位置对应的不同坐标轴的坐标值的总和的平均值作为合波形图终止位置的坐标。
[0092]
比如，当交叉点个数为零时，x视角分波形图的起始位置对应的坐标为(x1，y1)，y视角分波形图的起始位置对应的坐标为(x2，y2)，和z视角分波形图的起始位置对应的坐标为(x3，y3)，此时，合波形图的起始位置的坐标值为((x1+x2+x3)/3，(y1+y2+y3)/3)。合波形图的终止位置的坐标值与上述类似。
[0093]
(2)若交叉点的个数为一，则将交叉点对应的位置作为合波形图的起始位置和终止位置。
[0094]
(3)若交叉点的个数为二，则将两个交叉点的位置分别作为合波形图的起始位置和终止位置。
[0095]
(4)若交叉点的个数大于二，则选取交叉点的位置对应的y轴方向坐标值的差值最大的两个交叉点的位置，分别作为合波形图的起始位置和终止位置。
[0096]
2063：根据不同时间段的合波形图起始位置和终止位置，得到多次波在三维坐标系中的衍射波形图。
[0097]
本实施例中，通过每十秒记录的模拟多次波的分波形图，得到对应的交叉点，然后通过交叉点，得到合波形图。再将每分钟内的合波形图的起始位置和终止位置对应的点，顺着x轴方向连为一体，形成每分钟的合波形图的整图，再将各分钟的合波形图的整图连为一
体，便可对应得到衍射波形图。
[0098]
以上构建波形图的方法，能够简单且快速的确定出衍射波形图，使得衍射波形图的获取既简单又方便。
[0099]
s207：基于衍射波形图，得到地震数据的真振幅偏移图像。
[0100]
对得到的衍射波形图进行保幅成像剖面处理，得到真振幅偏移图像。这样能够使得整体波组特征更好，资料归位更加准确、断面形态更加清晰、内幕成像也得到改善，有效地提高了偏移成像的质量，同时也可以提高地震资料的垂向分辨率和横向分辨率。
[0101]
保幅成像剖面处理可以参见相关技术，在此不做具体限制。
[0102]
保幅成像剖面除了能够有效地补偿了中深层的振幅外，也可以提高部分浅层的资料质量，这样可以使得地震数据偏移成像的准确性高，确保移地震数据偏移成像的质量。
[0103]
图3是本公开实施例提供的一种真振幅的偏移成像装置的结构示意图，结合图3，该装置具有实现上述方法示例的功能，功能可以由硬件实现，也可以由硬件执行相应的软件实现。偏移成像装置包括：发射点确定模块31、衍射模块32、衍射波形图获取模块33和成像模块34。
[0104]
发射点确定模块31用于确定发射点在三维坐标系中的位置，发射点为地震波的震源所在的位置。衍射模块32用于将模拟多次波以发射点为衍射点进行衍射，模拟多次波是基于实际多次波生成的，实际多次波是指地震波经过地层多次衍射所形成的波。衍射波形图获取模块33用于获取模拟多次波在衍射过程中的衍射波形图。成像模块34用于基于衍射波形图，得到地震数据的真振幅偏移图像。
[0105]
可选地，衍射模块32用于：在模拟多次波以发射点为衍射点进行衍射的过程中，以相同时间间隔捕捉模拟多次波，得到多个衍射捕捉点；确定衍射捕捉点在三维坐标系中的位置；根据衍射捕捉点在三维坐标系中的位置，得到衍射捕捉点所形成的衍射捕捉图；根据衍射捕捉图，得到模拟多次波在三维坐标系中不同时间段所对应的分波形图，每个时间段对应不同视角的多个分波形图；根据分波形图，得到多次波在三维坐标系中的衍射波形图。
[0106]
可选地，衍射模块32还用于根据公式(3)确定衍射捕捉点在三维坐标系中的位置。
[0107]
可选地，衍射模块32还用于：将不同时间段所对应的分波形图以分波形图的起始位置所对应的坐标排列至三维坐标系的xoy平面中，得到x视角分波形图、y视角分波形图和z视角分波形图之间的交叉点的个数；根据交叉点的个数，得到不同时间段的合波形图的起始位置和终止位置，合波形图为x视角分波形图、y视角分波形图和z视角分波形图合在一起形成的合波形图；根据不同时间段的合波形图的起始位置和终止位置，得到多次波在三维坐标系中的衍射波形图。
[0108]
可选地，衍射模块32还用于：若所述交叉点的个数为零，则将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的起始位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图起始位置的坐标值，将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的终止位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图终止位置的坐标值；若所述交叉点的个数为一，则将所述交叉点对应的位置作为所述合波形图的起始位置和终止位置；若所述交叉点的个数为二，则将两个所述交叉点的位置分别作为所述合波形图的起始位置和终止位置；若所述交叉点的个数大于二，则选取所述交叉点的位置中对应的y轴方向坐标差最大的两个交叉点的位置，分别作为所述合波
形图的起始位置和终止位置。
[0109]
可选地，成像模块34用于对衍射波形图进行保幅成像剖面处理，得到真振幅偏移图像。
[0110]
需要说明的是：上述实施例提供的真振幅的偏移成像在进行偏移成像时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的真振幅的偏移成像装置与真振幅的偏移成像方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
[0111]
图4是本公开实施例提供的一种计算机设备的结构示意图，结合图4，计算机设备400可以包括以下一个或多个组件：处理器401、存储器402、通信接口403和总线404。
[0112]
处理器401包括一个或者一个以上处理核心，处理器401通过运行软件程序以及模块，从而执行各种功能应用以及信息处理。
[0113]
存储器402和通信接口403通过总线404与处理器401相连。存储器402可用于存储至少一个指令，处理器401用于执行该至少一个指令，以实现上述方法中的各个步骤。
[0114]
此外，存储器402可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，易失性或非易失性存储设备包括但不限于：磁盘或光盘，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，静态随时存取存储器(sram)，只读存储器(rom)，磁存储器，快闪存储器，可编程只读存储器(prom)。
[0115]
本公开实施例还提供了一种非临时性计算机可读存储介质，当所述存储介质中的指令由计算机设备的处理器执行时，使得计算机设备能够执行本公开实施例中提供的真振幅的偏移成像方法。
[0116]
本公开实施例还提供了一种计算机程序产品，包括计算机程序/指令，所述计算机程序/指令被处理器执行时实现本公开实施例中提供的真振幅的偏移成像方法。
[0117]
以上所述仅为本公开的可选实施例，并不用以限制本公开，凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。

技术特征：
1.一种真振幅的偏移成像方法，其特征在于，所述偏移成像方法包括：确定发射点在三维坐标系中的位置，所述发射点为地震波的震源所在的位置；将模拟多次波以所述发射点为衍射点进行衍射，所述模拟多次波是基于实际多次波生成的，所述实际多次波是指所述地震波经过地层多次衍射所形成的波；获取所述模拟多次波在衍射过程中的衍射波形图；基于所述衍射波形图，得到地震数据的真振幅偏移图像。2.根据权利要求1所述的偏移成像方法，其特征在于，所述获取所述模拟多次波在衍射过程中的衍射波形图，包括：在所述模拟多次波以所述发射点为衍射点进行衍射的过程中，以相同时间间隔捕捉所述模拟多次波，得到多个衍射捕捉点在所述三维坐标系中的位置；根据所述衍射捕捉点在所述三维坐标系中的位置，得到所述衍射捕捉点所形成的衍射捕捉图；根据所述衍射捕捉图，得到所述模拟多次波在所述三维坐标系中不同时间段所对应的分波形图，每个时间段对应不同视角的多个分波形图；根据所述分波形图，得到所述多次波在所述三维坐标系中的衍射波形图。3.根据权利要求2所述的偏移成像方法，其特征在于，所述根据所述衍射捕捉点在所述三维坐标系中的位置，得到所述衍射捕捉点所形成的衍射捕捉图，包括：根据扰动波场对所述衍射捕捉点在所述三维坐标系中的位置进行修正；根据修正后的所述衍射捕捉点在所述三维坐标系中的位置，得到所述衍射捕捉点所形成的衍射捕捉图。4.根据权利要求3所述的偏移成像方法，其特征在于，所述根据扰动波场对所述衍射捕捉点在所述三维坐标系中的位置进行修正，包括：根据以下公式得到所述衍射捕捉点在所述三维坐标系中的z轴方向坐标值；其中，f
a
为衍射捕捉点在所述三维坐标系中的z轴方向的坐标值，也为衍射捕捉点的振幅值；i
n
为所述发射点在所述三维坐标系中的z轴方向的坐标值e
xw
为扰动波场值；c
a
v为衍射捕捉点捕捉时在所述三维坐标系中对应的z轴方向的坐标值，所述z轴方向用于表示垂直于地面的方向。5.根据权利要求2所述的偏移成像方法，其特征在于，所述分波形图包括x视角分波形图、y视角分波形图和z视角分波形图，所述x视角分波形图为沿着所述三维坐标系中的x轴方向观察到的衍射捕捉图，所述y视角分波形图为沿着所述三维坐标系中的y轴方向观察到的衍射捕捉图，所述z视角分波形图为沿着所述三维坐标系中的z轴方向观察到的衍射捕捉图；所述根据所述分波形图，得到所述多次波在所述三维坐标系中的衍射波形图，包括：将不同时间段所对应的所述分波形图以所述分波形图的起始位置所对应的坐标排列至所述三维坐标系的xoy平面中，得到所述x视角分波形图、所述y视角分波形图和所述z视角分波形图之间的交叉点的信息；
根据所述交叉点的信息，得到不同时间段的合波形图的起始位置和终止位置，所述合波形图为所述x视角分波形图、所述y视角分波形图和所述z视角分波形图合在一起形成的合波形图；根据不同时间段的所述合波形图的起始位置和终止位置，得到所述多次波在所述三维坐标系中的衍射波形图。6.根据权利要求5所述的偏移成像方法，其特征在于，所述根据所述交叉点的信息，得到不同时间段的合波形图的起始位置和终止位置，包括：若所述交叉点的个数为零，则将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的起始位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图起始位置的坐标值，将所述x视角分波形图、所述y视角分波形图和所述z视角分波形图的终止位置对应的不同坐标轴的坐标值的总和的平均值作为所述合波形图终止位置的坐标值；若所述交叉点的个数为一，则将所述交叉点对应的位置作为所述合波形图的起始位置和终止位置；若所述交叉点的个数为二，则将两个所述交叉点的位置分别作为所述合波形图的起始位置和终止位置；若所述交叉点的个数大于二，则选取所述交叉点的位置中对应的y轴方向坐标差最大的两个交叉点的位置，分别作为所述合波形图的起始位置和终止位置。7.根据权利要求1所述的偏移成像方法，其特征在于，所述基于所述衍射波形图，得到地震数据的真振幅偏移图像，包括：对所述衍射波形图进行保幅成像剖面处理，得到所述真振幅偏移图像。8.一种真振幅的偏移成像装置，其特征在于，所述偏移成像装置包括：发射点确定模块，用于确定发射点在三维坐标系中的位置，所述发射点为地震波的震源所在的位置；衍射模块，用于将模拟多次波以所述发射点为衍射点进行衍射，所述模拟多次波是基于实际多次波生成的，所述实际多次波是指所述地震波经过地层多次衍射所形成的波；衍射波形图获取模块，用于获取所述模拟多次波在衍射过程中的衍射波形图；成像模块，用于基于所述衍射波形图，得到地震数据的真振幅偏移图像。9.一种计算机设备，其特征在于，所述计算机设备包括处理器和被配置为存储所述处理器可执行指令的存储器；所述处理器被配置为执行权利要求1至7任一项所述的真振幅的偏移成像方法。10.一种计算机存储介质，其上存储有计算机指令，其特征在于，所述计算机指令被处理器执行时实现权利要求1至7任一项所述的真振幅的偏移成像方法。

技术总结
本公开提供了一种真振幅的偏移成像方法、装置、设备及存储介质，属于油气勘探技术领域。所述偏移成像方法包括：确定发射点在三维坐标系中的位置，所述发射点为地震波的震源所在的位置；将模拟多次波以所述发射点为衍射点进行衍射，所述模拟多次波是基于实际多次波生成的，所述实际多次波是指所述地震波经过地层多次衍射所形成的波；获取所述模拟多次波在衍射过程中的衍射波形图；基于所述衍射波形图，得到地震数据的真振幅偏移图像。本公开通过该成像方法，能够大幅提高地震数据的真振幅偏移成像效率。像效率。像效率。


技术研发人员：韩嵩 陈康 屠志慧 杨广广 曾鸣 孔令霞 彭浩天 汤聪
受保护的技术使用者：中国石油天然气股份有限公司
技术研发日：2022.03.03
技术公布日：2023/9/13