弯曲射线叠前时间偏移速度求取方法及装置与流程

1.本发明涉及勘探地球物理技术领域，尤其涉及弯曲射线叠前时间偏移速度求取方法及装置。


背景技术：

2.本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
3.现阶段，叠前时间偏移技术仍然是地震资料成像处理中最主要的一种成像手段，其偏移结果被广泛地应用于地震解释及储层预测。时间偏移处理主要有两种方法，一种是传统的直射线叠前时间偏移，一种是弯曲射线叠前时间偏移。其中：直射线法使用rms速度，偏移算子为对称形，不能反映速度横向变化，因此在构造复杂、地层倾角大的工区往往不能得到理想的效果。实践中发展起来的弯曲射线偏移使用时间域层速度场，偏移算子为非对称形，算法在一定程度上考虑了速度场横向变速问题，实际应用中比传统的直射线法成像更合理。
4.现阶段的弯曲射线法速度更新方法通过拾取剩余时差，借助于dix公式实现速度更新，但是由于拾取误差的存在，加上垂向拾取没有考虑地层速度的横向变化，这样的速度更新会导致速度纵横向变化不合理，例如产生的横向速度突变等不合理现象。这种不合理现象目前的解决方案是对更新后的速度场进行平滑处理或使用约束速度反演的方式来消除这种速度场的病态变化。虽然经过处理后的速度场更符合地层的速度变化规律，但因为没有有效的质控手段，这样的更新方案也带来了一个问题：它们同时也会抹去合理的速度变化细节，导致迭代次数增多，大大的增加了工作量。
5.针对上述问题，目前尚未提出有效的解决方案。


技术实现要素：

6.本发明实施例提供一种弯曲射线叠前时间偏移速度求取方法，用以高效精确地求取弯曲射线叠前时间偏移速度，该方法包括：
7.获取给定弯曲射线叠前时间偏移速度；
8.以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：
9.在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；
10.将时间域数据转换到深度域，得到深度域数据；
11.根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；
12.将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；
13.在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度。
14.本发明实施例还提供一种弯曲射线叠前时间偏移速度求取装置，用以高效精确地求取弯曲射线叠前时间偏移速度，该装置包括：
15.获取单元，用于获取给定弯曲射线叠前时间偏移速度；
16.确定单元，用于以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：
17.在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；
18.将时间域数据转换到深度域，得到深度域数据；
19.根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；
20.将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；
21.在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度。
22.本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述弯曲射线叠前时间偏移速度求取方法。
23.本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述弯曲射线叠前时间偏移速度求取方法的计算机程序。
24.本发明实施例中，弯曲射线叠前时间偏移速度求取方案，与现有技术中通过拾取剩余时差，借助于dix公式实现速度更新，不准确且效率低的技术方案相比，通过：获取给定弯曲射线叠前时间偏移速度；以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；将时间域数据转换到深度域，得到深度域数据；根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度，实现了通过时深转换方法，将时间偏移的速度更新问题转换为深度域偏移速度更新问题，利用深度域速度网格层析技术迭代更新速度，速度在深度域优化更新后再转换到时间域，避免了直接在时间域优化弯曲射线叠前时间偏移速度时产生的横向速度突变等不合理现象，进而实现了高效准确地求取弯曲射线叠前时间偏移速度，对地震解释及储层预测具有重要的指导意义。
附图说明
25.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
26.图1为本发明实施例中弯曲射线叠前时间偏移速度求取方法的流程示意图；
27.图2为时间偏移域rms速度场；
28.图3是对图2用dix公式生成的层速度场；
29.图4是对图3层速度场的平滑结果；
30.图5是对图2用约束速度反演生成的层速度场；
31.图6是对图5层速度场的平滑结果；
32.图7是时间域初始层速度场；
33.图8是时间域初始层速度比例到深度域；
34.图9是时间偏移域层位显示；
35.图10是时间偏移域层位比例到深度域的结果；
36.图11是时间偏移域道集；
37.图12是时间偏移域道集比例到深度域的结果；
38.图13是深度域剩余谱横向拾取示例；
39.图14是深度域剩余谱纵向拾取示例；
40.图15是网格层析得到更新后的速度场(深度域)；
41.图16是图17的速度场比例到时间域的结果；
42.图17是速度优化前pstm剖面；
43.图18是速度优化后的stm剖面；
44.图19a为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之一(时间偏移域拾取剩余延迟，用dix公式更新速度，并对更新的速度场进行小平滑)建立的速度模型；
45.图19b为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之二(时间偏移域拾取剩余延迟，在dix公式基础上用约束速度反演方法更新速度)建立的速度模型；
46.图19c为采用本发明的稳定、高效的弯曲射线叠前时间偏移精确速度求取方法建立的速度模型；
47.图20是采用本发明稳定、高效的弯曲射线叠前时间偏移精确速度求取方法建立的速度模型；
48.图21为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之二(时间偏移域拾取剩余延迟，在dix公式基础上用约束速度反演方法更新速度)建立速度模型的偏移结果；
49.图22为采用本发明稳定、高效的弯曲射线叠前时间偏移精确速度求取方法建立的速度模型的偏移结果；
50.图23为本发明施例中弯曲射线叠前时间偏移速度求取的原理示意图；
51.图24为本发明施例中弯曲射线叠前时间偏移速度求取装置的结构示意图。
具体实施方式
52.为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
53.在介绍本发明实施例之前，首先对发明人发现技术问题从而提出本发明实施例方
案的思路进行介绍。
54.发明人发现：近几年，叠前深度域层速度网格层析建模技术发展迅速，并在大量的实际应用中见到了良好的效果。剩余延迟的自动拾取和高效的网格层析速度反演使得深度偏移的速度更新更加灵活和高效，纵横向剩余谱的拾取和质量控制使得深度偏移网格层析的质量得到很好的控制。因此，深度域层速度的网格层析迭代方法则比时间域的dix公式加约束速度反演更加有效，更新后的速度模型更合理。
55.本发明实施例主要为了克服目前在时间域通过拾取剩余时差，借助dix公式实现弯曲射线叠前时间偏移速度更新中存在的问题：1、垂向拾取无法考虑速度横向变化，导致速度更新后会在纵横向存在明显不合理；2、通过平滑或者约束反演消除速度更新产生不合理现象的同时会抹去合理的速度细节。因此发明人提出一种弯曲射线叠前时间偏移速度求取方案，该方案通过借助于目前比较成熟的深度域速度网格层析反演迭代技术，有效地完成弯曲射线叠前时间偏移速度的精确更新。即发明人将深度域速度更新的优点发明应用于时间域速度场更新，克服在时间域直接更新弯曲射线偏移速度模型时存在难点，借用深度域速度网格层析，完成时间域层速度迭代优化。该方案为一种弯曲射线法叠前时间偏移的精确速度求取方法，属于勘探地球物理领域中对地震资料进行处理解释的方法，适用于地震资料处理阶段的叠前时间偏移速度建模处理，与时间域弯曲射线叠前偏移相配套。下面对该弯曲射线叠前时间偏移速度求取方案进行详细介绍。
56.图1为本发明实施例中弯曲射线叠前时间偏移速度求取方法的流程示意图，如图1所示，该方法包括如下步骤：
57.步骤101：获取给定弯曲射线叠前时间偏移速度；
58.步骤102：以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：
59.步骤1021：在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；
60.步骤1022：将时间域数据转换到深度域，得到深度域数据；
61.步骤1023：根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；
62.步骤1024：将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；
63.步骤1025：在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度。
64.本发明实施例提供了一种高效、稳定的弯曲射线叠前时间偏移精确速度求取方法。通过时深转换，借助于目前比较高效、稳定的深度域网格层析速度优化迭代技术，在深度域完成速度优化更新后再转换到时间域，从而克服在时间域直接更新速度无法适应速度横向快速变化的缺点。经过多次时深转换和速度迭代优化后，获得精确的弯曲射线叠前偏移速度。
65.在一个实施例中，在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据，可以包括：在时间域应用当前周期速度进行目标弯曲射线叠前时间偏移处理，得到当前周期弯曲射线叠前时间域偏移速度、时间域偏移层位和时间域偏移道集；
66.将时间域数据转换到深度域，得到深度域数据，可以包括：将时间域数据转换到深度域，得到当前周期弯曲射线叠前深度域偏移速度、深度域偏移层位和深度域偏移道集；
67.根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度，可以包括：
68.根据当前周期弯曲射线叠前深度域偏移层位和深度域偏移道集，拾取目标弯曲射线剩余延迟数据；
69.基于当前周期弯曲射线叠前深度域偏移速度和剩余延迟数据，进行深度域速度网格层析反演，得到深度域速度网格层析反演结果；
70.根据深度域速度网格层析反演结果，进行深度域速度优化更新处理，得到优化更新后的深度域速度。
71.具体实施时，上述得到时间域数据，得到深度域数据，以及应用深度域速度网格层析反演技术优化更新后深度域速度的实施方式进一步提高了弯曲射线叠前时间偏移速度求取的精度和效率。
72.在一个实施例中，获取给定弯曲射线叠前时间偏移速度，可以包括：
73.利用rms速度来转换得到时间域偏移rms速度；
74.对所述时间域偏移rms速度进行约束速度反演，得到约束速度反演后的速度；
75.对约束速度反演后的速度进行野值编辑处理，得到所述给定弯曲射线叠前时间偏移速度。
76.具体实施时，上述得到所述给定弯曲射线叠前时间偏移速度的实施方式进一步提高了弯曲射线叠前时间偏移速度求取的精度。
77.在一个实施例中，将时间域数据转换到深度域，得到深度域数据，可以包括：对时间域数据进行构造约束层析反演到深度域，在深度域保留时间域构造形态。
78.具体实施时，对时间域偏移结果进行构造约束层析反演到深度域，在深度域仍然保留时间域构造形态，上述得到所述给定弯曲射线叠前时间偏移速度的实施方式进一步提高了弯曲射线叠前时间偏移速度求取的精度。
79.下面再结合附图2至图23举一例子，以便于理解本发明如何实施。
80.本发明提供的一种弯曲射线叠前时间偏移精确速度求取方法，借助于目前比较成熟的深度域速度网格层析优化技术，通过时深转换在深度域自动拾取网格点道集的剩余延迟，建立大型层析方程组并约束反演速度差，在深度域进行速度优化后再转换到时间域。相对于直接在时间域进行速度优化，本发明的速度求取方法能够适应速度在纵横向的变化，优化后的速度更加合理，速度变化细节更加丰富。
81.具体地，该弯曲射线叠前时间偏移精确速度求取方法借鉴深度域网格层析速度建模思路，方法包括应用初始叠前时间偏移速度进行目标线弯曲射线叠前时间偏移，生成目标线偏移道集及剖面；然后将初始速度、层位及偏移结果转换到深度域，再在深度域拾取目标线剩余延迟并进行网格层析反演，对初始速度进行优化更新；最后将更新后速度转换到时间域进行下一轮的速度更新迭代。本次发明充分考虑到了时间域速度场更新的难点，借助深度域速度网格层析优化更新的优势，很好的完成了时间域弯曲射线叠前时间偏移速度的优化迭代。下面进行详细介绍。
82.一、首先，介绍上述步骤101。
83.在上述步骤101中，为保证叠前时间偏移结果的合理性，需建立合理的时间域层速度场。一般情况下，用rms速度(图1)来转换得到时间域层速度；对比图2、图3、图4和图5，我们发现用约束速度反演(cvi)方法得到的速度更加符合地质规律，我们用对图5平滑处理后的速度(图6)做去野值编辑后得到初始的时间偏移层速度场(图7)。
84.二、接着，介绍上述步骤102。
85.1.首先，介绍上述步骤1021，即图23中的步骤201。
86.具体实施时，在时间域用图7的速度进行目标线弯曲射线叠前时间偏移，生成时间偏移道集(图9)和偏移剖面。在时间偏移剖面上拾取层位(图13)。
87.2.接着，介绍上述步骤1022，即图23中的步骤202。
88.具体实施时，将初始速度、层位、偏移结果转换到深度域，形成深度域道集、速度、层位：将时间偏移域数据道集比例(转换)到深度域(图12)，将时间偏移域速度比例到深度域(图8)，将时间偏移域层位比例到深度域(图10)。
89.3.接着，介绍上述步骤1023，该步骤在深度域拾取网格点剩余延迟并建立层析方程组，进行约束反演求解，应用反演求得速度差在深度域进行速度优化更新，即图23中的步骤203-205。
90.具体实施时，在步骤203中，进行深度域拾取目标线剩余延迟或剩余速度。可以在深度域横向沿层拾取剩余谱(图13)，也可以在深度域纵向拾取剩余谱(图14)。得到剩余延迟场或剩余速度场。
91.具体实施时，在步骤204中，进行深度域速度网格层析反演。基于图8的速度场和步骤203的剩余场，试验多组参数，进行网格层析运算。得到一组备选的速度场。
92.具体实施时，在步骤205中，进行深度域速度更新。根据网格层析试验结果，对比选取出合适的速度场，并根据需要对速度场进行平滑和去野值处理，得到更新后的速度场。
93.具体实施时，在该步骤1023中，对于转换到深度域的道集进行剩余延迟拾取并建立地震层析方程，通过约束层析反演求取速度差(剩余速度)，得到深度域更新后的速度，相较于在时间域的剩余能量谱拾取，剩余延迟拾取精度高，约束层析反演稳定性好，通过多次迭代可以获取精确速度。
94.具体实施时，地震层析指的是用从许多不同方向入射的地震波经地下物体反射和透射后的数据，对地下物体断面进行成像，跟医学ct或核磁共振成像原理相同。
95.为了便于理解，下面举一例子说明本步骤。
96.地震走时一定的情况下，一个小的速度扰动，在深度域将会产生一个剩余量与之对应，在网格层析中，用下式求得使远偏道集拉平的剩余速度δvi：
[0097][0098]
式中：为地层的真实深度，zh为炮检距h处深度，θ是反射的方位角，是倾角，v是初始速度。令z0为零偏处深度，为了得到真实速度，必修使z
h-z0最小，于是：
[0099][0100]
式中，th指地震波到达炮检距h处旅行时间，v0指的是炮检距为0处速度，vh指的是
炮检距h处的速度，同理θh指的是炮检距h处的反射方位角，θ0指的是炮检距0处的反射方位角。
[0101]
具体实施时，为了得到δvi，利用上式在每个地震网格点逐点求解，得到每个网格点的速度变量后再插值，得到剩余速度体，将这个剩余速度体与初始速度相加得到网格层析更新后速度模型，通过多次迭代，最终得到一个高精度的速度模型，实现了深度域的速度优化更新。
[0102]
4.接着，介绍上述步骤1024，即图23中的步骤206。
[0103]
具体实施时，将优化更新后深度域速度转换到时间域。用步骤206得到的速度进行弯曲射线叠前时间偏移，进入下一轮迭代。
[0104]
5.接着，介绍上述步骤1025，即图23中的步骤207。
[0105]
在满足预设迭代终止条件后，得到迭代优化后的速度，即最优弯曲射线叠前时间偏移速度，应用迭代优化后的速度进行最终时间域弯曲射线叠前偏移处理，即应用最优弯曲射线叠前时间偏移速度进行地震解释及储层预测，从而科学指导进行石油天然气钻井。
[0106]
下面对比弯曲射线叠前时间偏移现有速度更新技术和本发明的结果，可以了解本发明的应用效果。
[0107]
对比数据一：时间域弯曲射线叠前偏移应用现有速度更新技术之一(时间偏移域拾取剩余延迟，用dix公式更新速度，并对更新的速度场进行小平滑)的建模结果
[0108]
图19a-图19c是应用现有速度更新技术(时间偏移域拾取剩余延迟，用dix公式更新速度，并对更新的速度场进行小平滑)的建立的速度模型
[0109]
图19a为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之一(时间偏移域拾取剩余延迟，用dix公式更新速度，并对更新的速度场进行小平滑)建立的速度模型。
[0110]
图20为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之一(时间偏移域拾取剩余延迟，用dix公式更新速度，并对更新的速度场进行小平滑建立的速度模型的偏移结果。
[0111]
对比数据二：时间域弯曲射线叠前偏移应用现有速度更新技术之二(时间偏移域拾取剩余延迟，在dix公式基础上用约束速度反演方法更新速度)的建模结果
[0112]
图19b为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之二(时间偏移域拾取剩余延迟，在dix公式基础上用约束速度反演方法更新速度)建立的速度模型。
[0113]
图21为采用现有的时间域弯曲射线叠前偏移应用现有速度更新技术之二(时间偏移域拾取剩余延迟，在dix公式基础上用约束速度反演方法更新速度)建立速度模型的偏移结果。
[0114]
实施例数据：稳定、高效的弯曲射线叠前时间偏移精确速度求取方法
[0115]
图19c为采用本发明的稳定、高效的弯曲射线叠前时间偏移精确速度求取方法建立的速度模型。
[0116]
图22为采用本发明稳定、高效的弯曲射线叠前时间偏移精确速度求取方法建立的速度模型的偏移结果。
[0117]
对比图19a、图19b、图19c的速度场的结果和图20、图21和图22的偏移结果可见，传统方法一的建立的速度模型的野值控制能力差，偏移结果存在层位畸变；传统方法二的建
立的速度模型的野值控制能力好，符合地质规律，但缺少细节变化，偏移结果好但迭代次数多；本发明通过时深转换，借用深度域网格层析建模方法，得到的速度场野值控制好，符合地质规律，偏移结果且迭代次数少。图20、图21和图22中横纵坐标代表主测线和时间。
[0118]
综上，本发明实施例提供的弯曲射线叠前时间偏移速度求取方法的优点是：在时间域直接进行弯曲射线叠前时间偏移精确速度求取时，对剩余速度谱垂向拾取误差比较大，而且不能考虑速度横向变化，会导致更新后速度存在明显的异常，更新后速度质量并不理想。本发明中的一种弯曲射线叠前时间偏移精确速度更新方法，克服了常规方法的缺陷，在深度域进行速度网格层析反演后在转换回时间域，提高了弯曲射线叠前时间偏移速度精确度。
[0119]
本发明实施例中还提供了一种弯曲射线叠前时间偏移速度求取装置，如下面的实施例所述。由于该装置解决问题的原理与弯曲射线叠前时间偏移速度求取方法相似，因此该装置的实施可以参见弯曲射线叠前时间偏移速度求取方法的实施，重复之处不再赘述。
[0120]
图24为本发明施例中弯曲射线叠前时间偏移速度求取装置的结构示意图，如图24所示，该装置包括：
[0121]
获取单元01，用于获取给定弯曲射线叠前时间偏移速度；
[0122]
确定单元02，用于以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：
[0123]
在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；
[0124]
将时间域数据转换到深度域，得到深度域数据；
[0125]
根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；
[0126]
将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；
[0127]
在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度。
[0128]
在一个实施例中，在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据，包括：在时间域应用当前周期速度进行目标弯曲射线叠前时间偏移处理，得到当前周期弯曲射线叠前时间域偏移速度、时间域偏移层位和时间域偏移道集；
[0129]
将时间域数据转换到深度域，得到深度域数据，包括：将时间域数据转换到深度域，得到当前周期弯曲射线叠前深度域偏移速度、深度域偏移层位和深度域偏移道集；
[0130]
根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度，包括：
[0131]
根据当前周期弯曲射线叠前深度域偏移层位和深度域偏移道集，拾取目标弯曲射线剩余延迟数据；
[0132]
基于当前周期弯曲射线叠前深度域偏移速度和剩余延迟数据，进行深度域速度网格层析反演，得到深度域速度网格层析反演结果；
[0133]
根据深度域速度网格层析反演结果，进行深度域速度优化更新处理，得到优化更新后的深度域速度。
[0134]
在一个实施例中，所述获取单元具体用于：
[0135]
利用rms速度来转换得到时间域偏移rms速度；
[0136]
对所述时间域偏移rms速度进行约束速度反演，得到约束速度反演后的速度；
[0137]
对约束速度反演后的速度进行野值编辑处理，得到所述给定弯曲射线叠前时间偏移速度。
[0138]
在一个实施例中，将时间域数据转换到深度域，得到深度域数据，包括：对时间域数据进行构造约束层析反演到深度域，在深度域保留时间域构造形态。
[0139]
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述弯曲射线叠前时间偏移速度求取方法。
[0140]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述弯曲射线叠前时间偏移速度求取方法的计算机程序。
[0141]
本发明实施例中，弯曲射线叠前时间偏移速度求取方案，与现有技术中通过拾取剩余时差，借助于dix公式实现速度更新，不准确且效率低的技术方案相比，通过：获取给定弯曲射线叠前时间偏移速度；以给定弯曲射线叠前时间偏移速度作为初始弯曲射线叠前时间偏移速度迭代确定最优弯曲射线叠前时间偏移速度，每个迭代周期均执行以下操作：在时间域应用当前周期弯曲射线叠前时间偏移速度进行目标弯曲射线叠前时间偏移处理，得到时间域数据；将时间域数据转换到深度域，得到深度域数据；根据深度域数据，应用深度域速度网格层析反演技术，得到优化更新后的深度域速度；将优化更新后的深度域速度转换到时间域，得到优化更新后的时间域速度作为下一周期弯曲射线叠前时间偏移速度；在满足预设迭代终止条件后，得到最优弯曲射线叠前时间偏移速度，实现了通过时深转换方法，将时间偏移的速度更新问题转换为深度域偏移速度更新问题，利用深度域速度网格层析技术迭代更新速度，速度在深度域优化更新后再转换到时间域，避免了直接在时间域优化弯曲射线叠前时间偏移速度时产生的横向速度突变等不合理现象，进而实现了高效准确地求取弯曲射线叠前时间偏移速度，对地震解释及储层预测具有重要的指导意义。
[0142]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0143]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0144]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或
多个方框中指定的功能。
[0145]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0146]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。