一种PS波地震道数据的OVT道集处理方法、装置及电子设备与流程

一种ps波地震道数据的ovt道集处理方法、装置及电子设备
技术领域
1.本文件涉及地震勘探技术领域，尤其涉及一种ps波地震道数据的ovt道集处理方法、装置及电子设备。


背景技术：

2.宽方位高密度的地震勘探已成为地震勘探技术重要的发展。ovt技术是一种面向宽方位地震数据的新型面元抽取叠加技术，它是地下均匀照明成像的最小数据体，保留了宽方位采集中的偏移距和方位角信息，在pp波成像上已经取得到较好的应用。但在ps波成像尚缺少相应的成熟技术，原因主要是在于以下两个方面：1)首先ps波的下行波为p波、上行为s波，对应射线路径不对称，无法直接套用pp波ovt技术；2)其次ps波的转换点的位置随深度变化而变化，该点也无法套用ovt技术。
3.为此，有必要针对ps波开发ovt道集的应用，以在全新的数据角度的对地下构造特征进行成像。


技术实现要素：

4.本说明书实施例目的是提供一种ps波地震道数据的道集处理方法、装置及电子设备，能够对ps波地震道数据进行ovt道集抽取，从而得到全新角度的地震道集数据以对地下地质体进行成像。
5.为了实现上述目的，本说明书实施例是这样实现的：
6.第一方面，提供一种ps波地震道数据的道集处理方法，包括：
7.在ccp的面元网格上映射ps波地震道数据，并以所述ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距；
8.基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线由非双曲向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上；
9.在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
10.第二方面，提供一种ps波地震道数据的道集处理装置，包括：
11.ccp道集抽取模块，用于在ccp的面元网格上映射ps波地震道数据，并以所述ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距；
12.时序矫正模块，用于基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线向双曲时距关系进行校正，使得所述ps波地震道
数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上；
13.ovt道集抽取模块，用于在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
14.第三方面，提供一种电子设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行：
15.在ccp的面元网格上映射ps波地震道数据，并以所述ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距；
16.基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上；
17.在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
18.第四方面，提供一种算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如下步骤：
19.在ccp的面元网格上映射ps波地震道数据，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距，并以所述ps波地震道数据的成像点相对地表的投影点为中心；
20.基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上；
21.在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
22.本发明对ps波先进行ccp道集抽取，在ccp道集的基础上对ps波的时距曲线向双曲时距关系进行校正，使ps波的转换点的位置与虚拟炮检点的中心点相一致，能够以虚拟炮检点的中心点表示转换点的位置，从而在ovt的炮检线网格上，利用pp波的方式映射时距曲线校正后的ps波地震道数据，以适用于对ps波地震道数据的ovt道集抽取，可得到均匀覆盖的ps波ovt道集。本发明不仅以全新的数据角度对地下构造特征进行成像、刻画，此外，若基于抽取到的ps波ovt道集进行叠加成像，在维度信息上也要多于传统的共偏移距偏移成像方式，因此具有更高的成像精度。
附图说明
23.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明实施例中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
24.图1为本发明实施例提供的ps波地震道数据的道集处理方法的流程示意图。
25.图2为ccp道集抽取示意图。
26.图3为等效c波转换示意图。
27.图4为ps波地震数据的ovt道集示意图。
28.图5为ovt面元偏移距和方位角分布图。
29.图6为psv波偏移剖面。
30.图7为ps波地震道数据的ovt面元覆盖次数图。
31.图8为ps波地震道数据的ovt方位角道集示意图。
32.图9为各ps波地震道数据向异性校正后偏移剖面。
33.图10为ps波地震道数据模型ovt成像点道集示意图。
34.图11为ps波地震道数据的实际数据ovt成像道集示意图。
35.图12为ps波地震道数据的成像剖面对比示意图。
36.图13为ps波地震道数据的快慢横波剖面对比视图。
37.图14为本发明实施例提供的ps波地震道数据的道集处理装置的结构示意图。
38.图15为本发明实施例提供的电子设备的结构示意图。
具体实施方式
39.为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。
40.如前文所述，目前ps波在ovt道集成像缺少相应的成熟技术的原因主要是在于以下两个方面：1)首先ps波的下行波为p波、上行为s波，对应射线路径不对称，无法直接套用pp波ovt技术；2)其次ps波的转换点位置随深度变化而变化，该点也无法套用pp波ovt技术。
41.针对上述问题，本发明旨在针对ps波提出一种ovt道集的应用方案，具体能够抽取ps波均匀覆盖的ovt道集以用于叠加成像。不仅以从全新的数据角度对地下构造进行成像，相对于ps波传统的共偏移距偏移成像方式，由于维度信息得到了增加，因此具有更高的成像精度。
42.一方面，本发明实施例提供一种ps波地震道数据的道集处理方法，图1是该道集处理方法的流程图，包括如下步骤：
43.s102，在ccp的面元网格上映射ps波地震道数据，并以ps波地震道数据的成像点相对地表的投影点为中心，确定ps波地震道数据相对面元网格的虚拟炮检点，虚拟炮检点与ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距。
44.s104，基于ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对ps波地震道数据的时距曲线由非双曲向双曲时距关系进行校正，使得ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上。
45.应理解，校正后的地震道数据的共中心点即等效c波的共中心点。
46.s106，在ovt的炮检线网格上映射时距曲线校正后的ps波地震道数据，对ps波地震
道数据进行ovt道集抽取，得到ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
47.本发明实施例的方法对ps波先进行ccp道集抽取，在ccp道集的基础上对ps波的时距曲线向双曲时距关系进行校正，使ps波的转换点的位置与虚拟炮检点的中心点相一致，能够以虚拟炮检点的中心点表示转换点的位置，从而在ovt的炮检线网格上，利用pp波的方式映射时距曲线校正后的ps波地震道数据，以适用于对ps波地震道数据的ovt道集抽取，可得到均匀覆盖的ps波ovt道集。本发明不仅以全新的数据角度对地下构造特征进行成像、刻画，此外，若基于抽取到的ps波ovt道集进行叠加成像，在维度信息上也要多于传统的共偏移距偏移成像方式，因此具有更高的成像精度。
48.下面对本发明实施例方法的原理进行介绍。
49.本发明要解决的问题在于提供一种基于精确转点的ps波(本文ps波地震数据简称ps波)ovt道集抽取与成像的新方法，在成像精度得以提升的同时，ps波还能使用类似pp波的相关技术。其主要内容包括：1)精确转换点计算与等效c波转换；以及，ps波ovt道集抽取与成像两大内容。
50.(1)ps波精确转换点计算与等效c波转换
51.ps波的转换点位置随着深度发生变化，转换点在平面上的投影也会随着深度发生变化。在ps波ccp道集抽取过程中，一般的做法是利用面元分割的方法寻找转换点，快速抽取ccp道集。
52.如图2示，在将ps波地震数据映射ccp面元网格后，渐近转换点acp和检波点receiver之间的连线，与成像网格可确定出一系列交点，每个交点可以作为转换点，每两个交点所属的面元即为地震波反射面元(如图中阴影面元)。根据深度域上的转换点轨迹方程(如图中绿色曲线所示)可以计算得到每个交点对应的ps波反射时间t，作为该交点的旅行时。
53.旅行时t计算公式如下：
54.其中，为p波均方根速度，z0为转换点的深度，x
p
为转换点到原始炮检点中震源点之间的距离，x
′
为转换点到原始炮检点中检波点之间的距离。不同反射时间的地震数据便分割到不同的面元上。
55.由于转换点不是炮检连线的中点，并且不止一个，因此不能直接用炮检点所在的炮线和检波线来确定十字排列。利用了等效c波的方法，构建虚拟炮检点来解决这个问题。如图3所示，震源s激发的地震波在地下成像点o反射后传至地表检波器r点，o点在地表的投影点为c，深度为z0，炮检距为2h。保持炮检距大小不变，将炮检点向检波器一侧平移至s’、r’位置，s
′
、r
′
分别为s、r的虚拟点，使c为s
′
与r
′
的中点。假设原始炮检点旅行时为t，则从虚拟震源s
′
处激发，地震波传播至成像点o，再反射回地表虚拟检波器r
′
的旅行时为t
′
为的计算公式为：
56.57.其中，t0＝2z0/vc，v
p
为p波速度，vs为s波速度，z0为转换点的深度。
58.通过构建等效c波，将时间t对应点的值移动到时间t
′
的对应点，实现了ps波时距曲线由非双曲向双曲的转换。
59.2)ps波地震数据ovt道集抽取与成像
60.炮检坐标平移后的中心点即为反射点的位置，该反射点仍位于炮检移动前ccp成像面元内。利用平移后炮检坐标信息，便可采用pp波地震数据的方式进行ovt道集的抽取。与pp波不同的是，反射点位置在平面上的投影会随着深度发生横向偏移，因此来自同一组炮检对的地震数据会分配到不同的十字排列上。
61.图4为ps波地震数据ovt道集抽取示意图，假设反射点位于面元1内，炮检坐标平移后，便可找到虚拟炮点、检波点相应的炮线、检波线，并确定相应的十字排列中心1’处。同理，当反射点位于面元2、3、4等时，所确定的一系列十字排列中心位于蓝色圆圈处。通过上述方式，来自同一组炮检对的地震数据便分配到不同的十字排列上。
62.通过ps波ovt道集抽取，三维地震数据(纵向cmpx、横向cmpy、时间t)便分选到了五维(偏移距offset、方位角azimuth、cmpx、cmpy、t)。ps波地震数据基于ovt十字排列网格的叠加成像公式为：
[0063][0064]
其中，i(τ,x)为自激自收旅行时τ-空间坐标x＝(x,y,z)之间的关系图像，k为炮线号，l为检波线号，i为转换点面元的编号，y
x
表示cmp点的横坐标，w
k,l
(x,yy,y
x
)为震源点s＝(s
x
,sy)、检波点r＝(r
x
,ry)在空间x＝(x,y,z)处的振幅加权函数，u为地震数据，δ()为脉冲函数，t
#
(x,sy,r
x
,i；vc)为空间坐标x处以等效c波的速度vc，从震源点到成像点再到检波点的时长。
[0065]
基于该ps波地震道数据对应ovt道集的叠加成像结果，在保留方位角和偏移距信息的同时，还可对对地下构造特征进行识别。
[0066]
3)ps波方位各向异性校正
[0067]
ps波传播到各向异性介质中时会发生横波分裂现象，产生快横波ps1和慢横波ps2；对于含有单组垂直裂缝的hti介质，横波分裂还会与方位有关。为了进一步提高成像质量，需要校正各向异性的影响。假设裂缝方位为φ，在炮检方位为θ时，r分量、t分量与分裂后的ps1波、ps2波之间的关系为：
[0068]
其中，ps1(t)表示分裂后的快横波，ps2(t-δt)表示分裂后的慢横波，δt为ps1波与ps2波之间的时间延迟，裂缝方位φ是对u
t
(t,θ)进行极值求解所得到的,t为地震波的旅行时，r(t)和t(t)分别为沿着径向和切向方向接收到的ps波地震数据。
[0069]
基于上述公式，可以对ps波地震道数据中的快横波和慢横波进行分离，实现ps波各向异性的校正；
[0070]
其中，φ与δt的求取可通过参数扫描，求解目标函数u
t
(t,θ)的极值e(φ,δt)的
来获得：
[0071][0072]
为保持ps1波与ps2波之间的时差，将裂缝方位φ带入ps1波、ps2波之间的关系公式，得到分裂后的ps1波和ps2波。
[0073]
ps波地震数据的ovt道集抽取与叠加技术流程如图5所示。本发明在计算精确转换点的基础上，运用等效c波假设来进行ps波地震数据ovt道集抽取。该方法在ps波ovt处理与pp波ovt处理间建立起联系，为ps波地震数据ovt处理提供了便捷、高效的新实现方法与技术路线。
[0074]
在方法原理和技术流程的指导下，本发明运用理论模拟数据和实际数据来进行测试，通过和传统处理进行对比分析，来进一步证明本方法的可行性以及优越性。
[0075]
(1)数值模型测试
[0076]
为验证上述方法的正确性，设计了一个三维倾斜断层模型(图6)，模型参数如表1所示，层位1的p波速度、s波速度以及密度均为线性增大，层厚间隔为80m。采集时采用三维正交宽方位观测系统，炮线距、检波线距为50m，炮间距、道间距为50m；观测系统的一个patch示意图如图7所示；模拟采用12hz的雷克子波，采样间隔2ms，共1251个采样点。
[0077]
表1模型参数表
[0078][0079][0080]
在ovt道集叠加成像前，对数据进行了rt旋转、波场分离(lu et al，2012)预处理工作，获得了转换psv波。常规处理只包含偏移距域信息，而ovt的每个面元都有特定的偏移距和方位角范围，因此除偏移距外还包含了方位角信息。以该模型为例，当确定好ovt面元大小后，便可根据面元序号计算出如图5所示的面元中心位置对应的偏移距offset和方位角azimuth。
[0081]
本次ovt道集抽取时选取的面元大小为100*100m。原始psv波数据炮检分布规则、空间采样均匀，直接偏移便可获得较好的成像结果如图6的左边部分所示。对psv波每道地震数据进行转换点计算和等效c波转换后，便可进行ovt道集的分选。在完成ovt道集分选后运用公式(τ,x)进行成像，所得成像剖面图6的中间部分所示。经过本方法处理后的成像剖面与常规处理的剖面基本相同，但得到了更小的成像噪音。
[0082]
本方法与bale等传统基于acp点的ps波ovt处理方法相比，不仅在空间上实现了均匀照明，还提高了成像质量。图7的左边部分为pp波ovt面元的覆盖次数图，其在整个工区内都达到了均匀、完整的覆盖。图7的中间部分和右边部分为传统acp方法(即，基于acp点的ps波ovt处理方法)与本方法ovt面元覆盖次数结果。在acp方法中，选定v
p
/vs＝2，即ovt面元大
小为150m*75m。通过对比可发现，本方法与pp波ovt在空间上均到达了均匀覆盖，而acp法的覆盖次数在不同位置处差异较大。此外，由于acp点并不能准确的描述转换点的位置，在浅层位置ovt面元内反射点一致性变差，使得不能准确的进行同相叠加，进而影响成像效果。在以上因素的共同影响下，以acp点位置为基准的ovt道集抽取与本方法的相比，从成像剖面图上可以明显观察到在第一个反射界面处(对应图6箭头处)成像存在明显的不足。
[0083]
此外，由于包含有方位角信息，受第二层hti介质的影响，ovt成像道集上(图11)还可以观察到明显的各向异性特征。由于ps波具有低频、信噪比的特点，加之成像道集信噪比相对较低，为使计算的裂缝方位更加准确，各向异性校正时选择了在ovt方位角道集上进行。我们选取了30
°
的间隔，对偏移后的道集进行叠加。在ovt方位角道集上，同相轴(1800ms)表现为sv波呈“三角函数”状抖动如图8左上部分所示，且周期为180
°
，sh波每隔90
°
发生一次相位反转如图8右上部分所示。我们利用公式发生一次相位反转如图8右上部分所示。我们利用公式进行快慢波分离，经过快慢横波时差校正，方位各向异性影响得以消除，同相轴重新拉平如图8左下部分和右下部分所示。将消除周期性抖动后的方位角道集进行叠加，便最终得到了如图9所示的成像剖面。各向异性校正后，相对于各向异性校正前(图6中间部分)，在最后一个反射界面位置处(图6中间部分的箭头)纵向分辨率得以明显提升。
[0084]
(2)实际数据测试
[0085]
将方法运用到于安徽淮南煤田全方位三维三分量地震数据ps波处理。工区采集时设计炮间距、道间距为40m，炮线距、检波线距为90m，炮检坐标分布如图10所示。在进行ovt道集抽取前，已经对三分量地震数据进行了rt旋转等预处理。为与pp波抽取ovt道集方式相同，设置ps波的ovt面元大小也为两倍的炮线距，即180m*180m。
[0086]
r分量每道地震数据进行转换点计算和等效c波转换后，便可进行ovt道集的分选。在经过本方法ovt分选后，空间采样变得均匀，运用公式i(τ,x)进行成像，成像道集如图11所示，远近偏移距的能量差异变小。ovt处理叠加成像剖面如图12右边部分所示，相比于图12左边部分所示的直接偏移的叠加剖面，偏移噪声减弱，剖面更加清晰，分辨率也有所提高。
[0087]
实际地下介质情况复杂，煤系地层从浅到深往往发育着不同尺度和方位的裂缝，深部地层分裂出来的ps1波和ps2波经过了上覆地层的多次分裂。因此我们采用了滑动时窗的层剥离方法，利用公式通过60ms的滑动时窗逐层剥离和校正，来减小各向异性的影响。位于相邻测线位置处，ovt成像道集快慢横波分离前后的叠加剖面如图13所示。与分离前的r、t分量成像剖面对比，剖面质量得以提高，尤其浅层成像效果得到了明显改善。
[0088]
此外，对于图1所示的方法，本说明书实施例还提供一种ps波地震道数据的道集处理装置。图14是本说明书实施例的道集处理装置1400的结构示意图，包括：
[0089]
ccp道集抽取模块1410，用于在ccp的面元网格上映射ps波地震道数据，并以所述
ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距。
[0090]
时序矫正模块1420，用于基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上。
[0091]
ovt道集抽取模块1430，用于在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
[0092]
显然，本说明书实施例的装置可以作为上述图1所示方法的执行主体，因此能够实现该方法在图1所实现的步骤和功能。由于原理相同，本文不再赘述
[0093]
图5是本说明书的一个实施例电子设备的结构示意图。请参考图5，在硬件层面，该电子设备包括处理器，可选地还包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(random-access memory，ram)，也可能还包括非易失性存储器(non-volatile memory)，例如至少1个磁盘存储器等。当然，该电子设备还可能包括其他业务所需要的硬件。
[0094]
处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是isa(industry standard architecture，工业标准体系结构)总线、pci(peripheral component interconnect，外设部件互连标准)总线或eisa(extended industry standard architecture，扩展工业标准结构)总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图5中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
[0095]
存储器，用于存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。
[0096]
其中，处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行，在逻辑层面上形成上述ps波地震道数据的道集处理装置。对应地，处理器，执行存储器所存放的程序，并具体用于执行以下操作：
[0097]
在ccp的面元网格上映射ps波地震道数据，并以所述ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距。
[0098]
基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上。
[0099]
在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
[0100]
上述如本说明书图1所示实施例揭示的方法可以应用于处理器中，由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤
可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器，包括中央处理器(central processing unit，cpu)、网络处理器(network processor，np)等；还可以是数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field－programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
[0101]
应理解，本发明实施例的电子设备可以使业务处理装置实现对应于图1所示方法中的步骤和功能。由于原理相同，本文不再赘述。
[0102]
当然，除了软件实现方式之外，本说明书的电子设备并不排除其他实现方式，比如逻辑器件抑或软硬件结合的方式等等，也就是说以下处理流程的执行主体并不限定于各个逻辑单元，也可以是硬件或逻辑器件。
[0103]
此外，本发明实施例还提出了一种计算机可读存储介质，该计算机可读存储介质存储一个或多个程序，该一个或多个程序包括指令。
[0104]
其中，上述指令当被包括多个应用程序的便携式电子设备执行时，能够使该便携式电子设备执行图1所示方法的步骤，包括：
[0105]
在ccp的面元网格上映射ps波地震道数据，并以所述ps波地震道数据的成像点相对地表的投影点为中心，确定所述ps波地震道数据相对面元网格的虚拟炮检点，所述虚拟炮检点与所述ps波地震道数据相对面元网格的原始炮检点具有相同的炮检距。
[0106]
基于所述ps波地震道数据在面元网格反映的虚拟炮检点的旅行时，对所述ps波地震道数据的时距曲线由非双曲向双曲时距关系进行校正，使得所述ps波地震道数据相对面元网格反映的转换点位于校正后的ps波地震道数据的共中心点上。
[0107]
在ovt的炮检线网格上映射时距曲线校正后的所述ps波地震道数据，对所述ps波地震道数据进行ovt道集抽取，得到所述ps波地震道数据的ovt道集，其中，ovt道集抽取的信息包括：偏移距、方位角、cmp点的横坐标、cmp点的纵坐标和旅行时。
[0108]
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0109]
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
[0110]
以上仅为本说明书的实施例而已，并不用于限制本说明书。对于本领域技术人员来说，本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书的权利要求范围之内。此外，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本文件的保护范围。