成像方法、装置、电子设备和存储介质与流程

1.本发明涉及成像技术领域，尤其涉及一种成像方法、装置、电子设备和存储介质。


背景技术：

2.在海上工程地质下，介质具有较强的吸收衰减属性的勘察区，例如，含气溶洞、浅层气、采空区等，该类勘察区的地震资料往往表现出低频和弱振幅的特性。传统不考虑吸收衰减补偿的成像算法在该类勘察区的成像结果会出现分辨率下降的情况。所以，为了提高该类工程地质勘察的成像精度，成像过程中需要进行吸收衰减效应的补偿。
3.相关技术中，通常利用双程波外推算子的逆时偏移算法进行吸收衰减效应的补偿并成像。但该算法的计算复杂，从而导致成像效率低。


技术实现要素：

4.针对现有技术存在的问题，本发明提供一种成像方法、装置、电子设备和存储介质。
5.本发明提供一种成像方法，包括：
6.基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
7.基于所述中心射线的旅行时确定射线束的旅行时；
8.对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述复旅行时包括成像点的旅行时和成像点的辅助旅行时；
9.基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
10.基于每个所述成像点的旅行时和每个所述补偿后的辅助旅行时进行成像。
11.根据本发明提供的一种成像方法，所述基于炮点坐标确定中心射线的旅行时，包括：
12.基于公式(1)至公式(5)确定所述中心射线的旅行时；
[0013][0014][0015][0016][0017][0018]
其中，x为所述中心射线在笛卡尔坐标系中的水平坐标，z为所述中心射线在笛卡尔坐标系中的垂直坐标，p
x
为中心射线传播方向的水平分量，pz为中心射线传播方向的垂直
分量，v为中心射线路径上的声波速度，s为中心射线路径，τ为中心射线的旅行时。
[0019]
根据本发明提供的一种成像方法，所述基于所述中心射线的旅行时确定射线束的旅行时，包括：
[0020]
基于公式(6)确定所述射线束的旅行时；
[0021]
τ(q,s)＝τ(s)+m(s)q2/2
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0022]
其中，τ(q,s)为所述射线束在射线中心坐标系中的旅行时，τ(s)为所述中心射线的旅行时，τ(s)＝τ(0,s)，m为旅行时关于坐标q的二阶导数，p和q为经典动力学射线追踪的辅助变量，线追踪的辅助变量，v0为所述中心射线上的声波速度。
[0023]
根据本发明提供的一种成像方法，所述对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时，包括：
[0024]
基于公式(7)对所述射线束的旅行时进行振幅和相位补偿，得到所述射线束的复旅行时；
[0025][0026]
其中，x为笛卡尔坐标系下的成像点，τc(x)为射线束在笛卡尔坐标系下的复旅行时，τ(x)为所述射线束在笛卡尔坐标系中的旅行时，τq(x)为辅助旅行时，i为虚数，w0为参考频率，w为频率，q为表征吸收衰减的品质因子，为振幅补偿，为相位补偿，
[0027]
根据本发明提供的一种成像方法，所述基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿，包括：
[0028]
基于公式(8)对每个成像点的辅助旅行时进行补偿；
[0029][0030]
其中，lf-1
为短时傅里叶反变换，lf为短时傅里叶变换，d(τ)为补偿前的地震道，τc为所述射线束的复旅行时，τq为辅助旅行时，d(τ,τq)为不同辅助旅行时补偿后的地震道。
[0031]
根据本发明提供的一种成像方法，所述基于每个所述旅行时和每个所述补偿后的辅助旅行时进行成像，包括：
[0032]
基于公式(9)进行成像；
[0033][0034]
其中，i
shot
(x)为成像结果，ps为炮点出射的射线束的射线参数，pr为检波点出射的射线束的射线参数，w(x)为振幅加权脉冲函数，d为时间域波场，xs(ps)为炮点坐标，x
r0
为检波点的中心坐标，ω为成像孔径，τs为上行波旅行时，τr为下行波旅行时，τ
qs
为上行波辅助旅行时，τ
qr
为下行波辅助旅行时。
[0035]
本发明还提供一种成像装置，包括：
[0036]
第一确定单元，用于基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
[0037]
第二确定单元，用于基于所述中心射线的旅行时确定射线束的旅行时；
[0038]
第一补偿单元，用于对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述复旅行时包括成像点的旅行时和辅助旅行时；
[0039]
第二补偿单元，用于基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
[0040]
成像单元，用于基于每个所述旅行时和每个所述补偿后的辅助旅行时进行成像。
[0041]
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述成像方法的步骤。
[0042]
本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述成像方法的步骤。
[0043]
本发明还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述成像方法的步骤。
[0044]
本发明提供的一种成像方法、装置、电子设备和存储介质，在确定射线束的旅行时时，首先对射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；再基于正反短时傅里叶变换对射线束的复旅行时中的每个成像点的辅助旅行时进行补偿；最终基于每个成像点的旅行时和每个补偿后的辅助旅行时进行成像。可知，本发明是基于正反短时傅里叶变换在频率域对辅助旅行时进行补偿后再成像，从而在提高成像分辨率的同时，提高了成像效率。
附图说明
[0045]
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0046]
图1是本发明提供的成像方法的流程示意图；
[0047]
图2是本发明提供的射线中心坐标系的示意图；
[0048]
图3是本发明提供的成像装置的结构示意图；
[0049]
图4是本发明提供的电子设备的实体结构示意图。
具体实施方式
[0050]
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0051]
下面结合图1和图2描述本发明的成像方法。
[0052]
图1是本发明提供的成像方法的流程示意图，应用于粘声介质的成像场景中，如图1所示，该成像方法包括以下步骤：
[0053]
步骤101、基于炮点坐标确定中心射线的起始位置并通过射线追踪得到中心射线旅行时。
[0054]
示例地，从地震数据中获取每个炮点坐标，基于每个炮点坐标确定炮端的中心射线的起始位置，进而基于中心射线的起始位置采用射线追踪方法计算得到中心射线的旅行时。
[0055]
可选地，基于公式(1)至公式(5)确定所述中心射线的旅行时。
[0056][0057][0058][0059][0060][0061]
其中，x为中心射线在笛卡尔坐标系中的水平坐标，z为中心射线在笛卡尔坐标系中的垂直坐标，p
x
为中心射线传播方向的水平分量，pz为中心射线传播方向的垂直分量，v为中心射线路径上的声波速度，s为中心射线路径，τ为中心射线的旅行时。
[0062]
示例地，地震数据成像可以分成两个部分：正算子的波场传播和成像条件的施加。本发明利用傍轴近似简化的射线束构造正传播算子来计算旅行时。
[0063]
具体地，利用特征线法求解波动方程，可以得到表征地震波场运动学特征的射线追踪系统，利用该射线追踪系统可以获得中心射线的坐标和中心射线的旅行时。具体采用上述公式(1)至公式(5)求解中心射线的旅行时。
[0064]
步骤102、基于所述中心射线的旅行时确定射线束的旅行时。
[0065]
其中，射线束是指中心射线附近一定范围的区域。
[0066]
可选地，基于公式(6)确定所述射线束的旅行时。
[0067]
τ(q,s)＝τ(s)+m(s)q2/2
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0068]
其中，τ(q,s)为所述射线束在射线中心坐标系中的旅行时，τ(s)为所述中心射线的旅行时，τ(s)＝τ(0,s)，m为旅行时关于坐标q的二阶导数，p和q为经典动力学射线追踪的辅助变量，线追踪的辅助变量，v0为所述中心射线上的声波速度，其中，τ(0,s)中的q等于0，则说明τ(s)为中心射线。
[0069]
示例地，图2是本发明提供的射线中心坐标系的示意图，在得到中心射线的旅行时的情况下，在图2所示的射线中心坐标系中，利用傍轴射线理论计算中心射线附件的射线束的旅行时。具体是基于公式(6)计算射线束的旅行时，其中，辅助变量p和q的值可以通过求
解经典动力学射线追踪系统获得，即基于和获得。
[0070]
需要说明的是，图2中的x0表示中心射线的起始位置，s是沿中心射线的弧长，q表示垂直于中心射线的距离。
[0071]
步骤103、对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时。
[0072]
其中，所述复旅行时包括成像点的旅行时和成像点的辅助旅行时。
[0073]
可选地，基于公式(7)对所述射线束的旅行时进行振幅和相位补偿，得到所述射线束的复旅行时。
[0074][0075]
其中，x为笛卡尔坐标系下的成像点，τc(x)为射线束在笛卡尔坐标系下的复旅行时，τ(x)为射线束在笛卡尔坐标系中的实旅行时，τq(x)为辅助旅行时，i为虚数，w0为参考频率，w为频率，q为表征吸收衰减的品质因子，也就是上述所提到的辅助变量，为振幅补偿，为相位补偿，
[0076]
其中，品质因子q是通过其他方法得到的，例如层析和质心频移等方法，参考频率w0为地震数据的主频，频率w的范围为地震数据的频带范围。
[0077]
示例地，在粘声介质中，可以通过引入复速度的方式来计算一个复旅行时，该复旅行时可以采用公式(7)来表达，可以看出，公式(7)中包括沿射线的实旅行时τ(x)和辅助旅行时τq(x)，所以，该复旅行时可以通过求解实旅行时τ(x)和辅助旅行时τq(x)获得，具体求解关于τ(x)和τq(x)的这两个公式即可获得τ(x)和τq(x)，将τ(x)和τq(x)代入公式(7)中，即可计算得到射线束的复旅行时。
[0078]
需要说明的是，τ(x)和τ(q,s)实质上指示一个意思，只是τ(q,s)为射线束在射线中心坐标系中的旅行时，沿着射线看；而τ(x)为射线束在笛卡尔坐标系中的实旅行时，是在笛卡尔坐标系下看。
[0079]
步骤104、基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿。
[0080]
可选地，基于公式(8)对每个成像点的辅助旅行时进行补偿。
[0081][0082]
其中，lf-1
为短时傅里叶反变换，lf为短时傅里叶变换，d(τ)为补偿前的地震道，τc为所述射线束的复旅行时，τq为辅助旅行时，d(τ,τq)为不同辅助旅行时补偿后的地震道。
[0083]
示例地，射线束成像可以看成基尔霍夫(kirchhoff)积分成像的一个推广，经典的基于炮道集的射线束成像方程可以采用以下公式(10)表示：
[0084][0085]
其中，i
shot
(x)为成像结果，ps为炮点出射的射线束的射线参数，pr为检波点出射的射线束的射线参数，w(x)为振幅加权脉冲函数，xs(ps)为炮点坐标，x
r0
为检波点的中心坐标，ω为成像孔径，d表示频域局部平面波场。
[0086]
在声介质条件下，旅行时是实数，公式(10)的频域积分可以由时间域波场d与一个振幅加权脉冲函数褶积得到。换言之，可以通过引入激励成像条件τ＝τs+τr而省略公式(10)的频域积分，具体采用以下公式(11)表示：
[0087][0088]
其中，d为时间域波场。
[0089]
但是，上述射线束的复旅行时是复数，在粘声介质中不能将公式(10)的频域积分通过在时间域波场与振幅加权脉冲函数做褶积来表达，且需要在时间频率域进行补偿，其时频变换和补偿的处理在成像内层循环，因此，粘声介质成像的计算复杂度远大于公式(11)的计算，从而严重影像成像效率。
[0090]
傅里叶分析是频域信号处理的常用方法，但其对非平稳信号的适用性有一定的限制。利用局部高斯函数加窗的短时傅里叶变换，相对于傅里叶变换对非平稳信号有更好的适应性，而且加窗的短时傅里叶变换还能够提高成像效率。具体利用短时傅里叶变换在频率域进行补偿的过程可采用以下公式(12)和公式(13)来表示：
[0091]
d(w)＝lf(d(τ))＝∫
τ
d(τ)e-iwt
dτ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0092]
d(τc)＝lf-1
(d(w))＝∫wd(w)e
iwc
dw
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0093]
其中，d(τ)为补偿前的地震道，d(τc)为补偿后的复地震道，通常，取d(τc)的实部代入上述公式(11)就可以进行基于空变q模型补偿的射线束成像。但是，上述公式(12)和公式(13)的短时傅里叶正反变换是在成像过程的内循环中进行的，这样时间和频率的积分会大大降低成像效率。
[0094]
鉴于此，本发明基于短时傅里叶变换对每个成像点的辅助旅行时进行补偿，具体采用上述公式(8)对辅助旅行时进行补偿，从公式(8)可以看出，补偿后的数据不仅是τ的函数，还是τc的函数，该数据能够体现不同程度的补偿结果，补偿后的数据从一维变为二维，最后将该二维数据全部保存在设备内存中，便于后续成像时使用。
[0095]
本发明提供的成像方法，基于传统基尔霍夫(kirchhoff)积分偏移实现框架，引入高效的短时傅里叶变换算子对数据进行时频变换，在频率域进行吸收衰减效应的补偿，即振幅和相位的补偿。
[0096]
步骤105、基于每个所述成像点的旅行时和每个所述补偿后的辅助旅行时进行成像。
[0097]
可选地，基于公式(9)进行成像。
[0098][0099]
其中，i
shot
(x)为成像结果，ps为炮点出射的射线束的射线参数，pr为检波点出射的射线束的射线参数，w(x)为振幅加权脉冲函数，d为时间域波场，xs(ps)为炮点坐标，x
r0
为检波点的中心坐标，ω为成像孔径，τs为上行波旅行时，τr为下行波旅行时，τ
qs
为上行波辅助旅行时，τ
qr
为下行波辅助旅行时。
[0100]
示例地，在利用正算子进行射线束传播的过程中，基于空变模型计算炮点到成像点准确的旅行时和辅助旅行时，其中，模型数值随着空间x和z变化，基于空变模型计算得到
的准确的旅行时和辅助旅行时用于成像；再对满足条件τq＝τ
qs
+τ
qr
的数据利用激励成像条件τ＝τs+τr进行成像，具体可采用上述公式(9)进行成像。对比公式(9)和公式(10)可以发现二者具有相似的计算复杂度，由于τ和τq都是关于炮检点到地下成像点x的函数，傅里叶变换到频率域进行补偿与τq有关，因此，必须是在成像空间x的内循环进行傅里叶变换；然后补偿再做傅里叶反变换；公式(9)的策略是提前补偿好辅助旅行时，按照计算好的τq(基于空变模型计算得到)去找补偿后的数据就可以了，所以补偿可以在成像空间x循环之外。即公式(9)的成像策略将成像空间循环内的时频变换转移到了成像空间循环外，从而极大的提高了成像效率。
[0101]
本发明提供的一种成像方法，在确定射线束的旅行时时，首先对射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；再基于正反短时傅里叶变换对射线束的复旅行时中的每个成像点的辅助旅行时进行补偿；最终基于每个成像点的旅行时和每个补偿后的辅助旅行时进行成像。可知，本发明是基于正反短时傅里叶变换在频率域对辅助旅行时进行补偿后再成像，从而在提高成像分辨率的同时，提高了成像效率。即本发明在粘声介质偏移过程中采用先补偿后成像的策略，该策略使得成像过程中避免频率积分，计算复杂度与传统不考虑吸收衰减补偿的射线类成像方法的计算复杂度基本一致，从而提高了成像效率。
[0102]
下面对本发明提供的成像装置进行描述，下文描述的成像装置与上文描述的成像方法可相互对应参照。
[0103]
图3是本发明提供的成像装置的结构示意图，如图3所示，该成像装置包括第一确定单元3013、第二确定单元302、第一补偿单元303、第二补偿单元304和成像单元305；其中：
[0104]
第一确定单元301，用于基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
[0105]
第二确定单元302，用于基于所述中心射线的旅行时确定射线束的旅行时；
[0106]
第一补偿单元303，用于对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述复旅行时包括成像点的旅行时和辅助旅行时；
[0107]
第二补偿单元304，用于基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
[0108]
成像单元305，用于基于每个所述旅行时和每个所述补偿后的辅助旅行时进行成像。
[0109]
本发明提供的一种成像装置，在确定射线束的旅行时时，首先对射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；再基于正反短时傅里叶变换对射线束的复旅行时中的每个成像点的辅助旅行时进行补偿；最终基于每个成像点的旅行时和每个补偿后的辅助旅行时进行成像。可知，本发明是基于正反短时傅里叶变换在频率域对辅助旅行时进行补偿后再成像，从而在提高成像分辨率的同时，提高了成像效率。
[0110]
基于上述任一实施例，所述第一确定单元301具体用于：
[0111]
基于公式(1)至公式(5)确定所述中心射线的旅行时；
[0112]
[0113][0114][0115][0116][0117]
其中，x为中心射线在笛卡尔坐标系中的水平坐标，z为中心射线在笛卡尔坐标系中的垂直坐标，p
x
为中心射线传播方向的水平分量，pz为中心射线传播方向的垂直分量，v为中心射线路径上的声波速度，s为中心射线路径，τ为中心射线的旅行时。
[0118]
基于上述任一实施例，所述第二确定单元302具体用于：
[0119]
基于公式(6)确定所述射线束的旅行时；
[0120]
τ(q,s)＝τ(s)+m(s)q2/2
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0121]
其中，τ(q,s)为射线束在射线中心坐标系中的旅行时，τ(s)为所述中心射线的旅行时，τ(s)＝τ(0,s)，m为旅行时关于坐标q的二阶导数，p和q为经典动力学射线追踪的辅助变量，踪的辅助变量，v0为所述中心射线上的声波速度。
[0122]
基于上述任一实施例，所述第一补偿单元303具体用于：
[0123]
基于公式(7)对所述射线束的旅行时进行振幅和相位补偿，得到所述射线束的复旅行时；
[0124][0125]
其中，x为笛卡尔坐标系下的成像点，τc(x)为射线束在笛卡尔坐标系下的复旅行时，τ(x)为射线束在笛卡尔坐标系中的实旅行时，τq(x)为辅助旅行时，i为虚数，w0为参考频率，w为频率，q为表征吸收衰减的品质因子，为振幅补偿，为相位补偿，
[0126]
基于上述任一实施例，所述第二补偿单元304具体用于：
[0127]
基于公式(8)对每个成像点的辅助旅行时进行补偿；
[0128][0129]
其中，lf-1
为短时傅里叶反变换，lf为短时傅里叶变换，d(τ)为补偿前的地震道，τc为所述射线束的复旅行时，τq为辅助旅行时，d(τ,τq)为不同辅助旅行时补偿后的地震道。
[0130]
基于上述任一实施例，所述成像单元305具体用于：
[0131]
基于公式(9)进行成像；
[0132][0133]
其中，i
shot
(x)为成像结果，ps为炮点出射的射线束的射线参数，pr为检波点出射的射线束的射线参数，w(x)为振幅加权脉冲函数，d为时间域波场，xs(ps)为炮点坐标，x
r0
为检波点的中心坐标，ω为成像孔径，τs为上行波旅行时，τr为下行波旅行时，τ
qs
为上行波辅助旅行时，τ
qr
为下行波辅助旅行时。
[0134]
图4是本发明提供的电子设备的实体结构示意图，如图4所示，该电子设备可以包括：处理器(processor)410、通信接口(communications interface)420、存储器(memory)430和通信总线440，其中，处理器410，通信接口420，存储器430通过通信总线440完成相互间的通信。处理器410可以调用存储器430中的逻辑指令，以执行成像方法，该方法包括：基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
[0135]
基于所述中心射线的旅行时确定射线束的旅行时；
[0136]
对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述复旅行时包括成像点的旅行时和成像点的辅助旅行时；
[0137]
基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
[0138]
基于每个所述成像点的旅行时和每个所述补偿后的辅助旅行时进行成像。
[0139]
此外，上述的存储器430中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0140]
另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的成像方法，该方法包括：基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
[0141]
基于所述中心射线的旅行时确定射线束的旅行时；
[0142]
对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述复旅行时包括成像点的旅行时和成像点的辅助旅行时；
[0143]
基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
[0144]
基于每个所述成像点的旅行时和每个所述补偿后的辅助旅行时进行成像。
[0145]
又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的成像方法，该方法包括：基于炮点坐标确定中心射线的起始位置，并通过射线追踪得到所述中心射线的旅行时；
[0146]
基于所述中心射线的旅行时确定射线束的旅行时；
[0147]
对所述射线束的旅行时进行振幅和相位补偿，得到射线束的复旅行时；其中，所述
复旅行时包括成像点的旅行时和成像点的辅助旅行时；
[0148]
基于正反短时傅里叶变换对每个成像点的辅助旅行时进行补偿；
[0149]
基于每个所述成像点的旅行时和每个所述补偿后的辅助旅行时进行成像。
[0150]
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
[0151]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
[0152]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。