CN105425290A - 一种叠前时间偏移的方法及装置 - Google Patents

Abstract

本发明提供一种叠前时间偏移的方法及装置。所述方法包括：从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集；根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时；提取所述记录道在所述旅行时处对应的振幅值；根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。利用所述本发明实施例提供的技术方案，通过从成像点的实际位置处计算偏移振幅，避免了将地震数据校正到统一基准面上而带来的计算旅行路径的误差，提高了叠前时间偏移的成像效率和成像质量。

Description

一种叠前时间偏移的方法及装置

技术领域

本发明涉及地震成像处理技术领域，尤其涉及一种叠前时间偏移的的方法及装置。

背景技术

叠前时间偏移可以作为复杂构造成像和速度分析的重要手段，已经广泛应用于石油勘探的叠前偏移技术中。叠前时间偏移具有较高的成像精度和较低的计算成本，从石油公司角度来说，叠前时间偏移可以提高石油开发效益，减少勘探开发风险。

现有技术中通常采用积分法计算叠前时间偏移，主要包括两个步骤：第一步是按照记录道中地震波的旅行路径求走时，第二步是按照所述走时求取所有记录道的地震波能量之和，再将所述地震波能量之和转移到成像点上。从上述两个步骤考虑，影响叠前时间偏移精度的计算步骤主要是第一步，即计算地震波的走时。地震波的走时主要取决于所引用成像点处的均方根速度的精度及地震波的旅行路径与地震波真实的传播路径的接近程度。现有技术中通常采用直射线走时计算法求取地震波的走时，所述直射线走时计算法基于统一水平基准面，也就是可以将地震数据和速度场校正到统一基准面上。但是，由地表到统一水平基准面的校正往往会给速度分析的精度以及地震波的旅行路径造成误差，进而造成走时及速度误差，影响了成像的精度。尤其对于起伏地表的校正，往往需要的校正量较大，影响成像精度的程度也越大。

现有技术中利用积分法计算叠前时间偏移需要将地表校正到统一水平基准面上，造成走时及速度误差，影响成像的精度，因此，现有技术亟需一种成像精度较高的叠前时间偏移方法。

发明内容

本申请的目的在于提供一种叠前时间偏移的方法及装置，可以提高叠前时间偏移的成像精度。

本申请提供的一种叠前时间偏移的方法及装置具体是这样实现的：

一种叠前时间偏移的方法，所述方法包括：

从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集；

根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时；

提取所述记录道在所述旅行时处对应的振幅值；

根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

可选的，在本发明的一个实施例中，所述记录道的地震数据包括：

炮点、检波点、成像点的坐标、地表高程；

成像点到地表的单程垂直时间；

成像点处的速度。

可选的，在本发明的一个实施例中，所述计算所述记录道在所述成像点处的旅行时的计算公式包括下式：

式中，T_ij为所述成像点第i个共中心点道集中第j个记录道在成像点处的旅行时；t为成像点到地表的单程垂直时间，v为成像点处的速度；Eleve_cmpijs为第i个共中心点道集中第j个记录道炮点的地面高程；Eleve_cmpijr为第i个共中心点道集中第j个记录道检波点的地面高程；Eleve_crp可以为成像点处的地表高程；X_ijs为第i个共中心点道集中第j个记录道炮点到成像点的水平距离；X_ijr为第i个共中心点道集中第j个记录道检波点到成像点的水平距离。

可选的，在本发明的一个实施例中，在所述从地震数据中获取成像点的共中心点以及共中心点道集之前，还包括：

建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线间距、最小炮检距、最大炮检距；

根据所述观测系统采集地震数据。

可选的，在本发明的一个实施例中，所述偏移孔径的大小的范围包括：所述偏移孔径的大小为所述最大炮检距的两倍。

一种叠前时间偏移的装置，所述装置包括：

道集获取单元，用于从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集；

旅行时计算单元，用于根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时；

振幅值提取单元，用于提取所述记录道在所述旅行时处对应的振幅值；

偏移振幅获取单元，用于根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

可选的，在本发明的一个实施例中，所述记录道的地震数据包括：

炮点、检波点、成像点的坐标、地表高程；

成像点到地表的单程垂直时间；

成像点处的速度。

可选的，在本发明的一个实施例中，所述计算所述记录道在所述成像点处的旅行时的计算公式包括下式：

式中，T_ij为所述成像点第i个共中心点道集中第j个记录道在成像点处的旅行时；t为成像点到地表的单程垂直时间，v为成像点处的速度；Eleve_cmpijs为第i个共中心点道集中第j个记录道炮点的地面高程；Eleve_cmpijr为第i个共中心点道集中第j个记录道检波点的地面高程；Eleve_crp为成像点处的地表高程；X_ijs为第i个共中心点道集中第j个记录道炮点到成像点的水平距离；X_ijr为第i个共中心点道集中第j个记录道检波点到成像点的水平距离。

可选的，在本发明的一个实施例中，所述装置还包括：

观测系统建立单元，用于建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线间距、最小炮检距、最大炮检距；

地震数据采集单元，用于根据所述观测系统采集地震数据。

可选的，在本发明的一个实施例中，所述偏移孔径的大小的范围包括：所述偏移孔径的大小为所述最大炮检距的两倍。

由此可见，本申请实施例的技术方案通过通过从成像点的实际位置处计算偏移振幅，避免了将地震数据校正到统一基准面上而带来的计算旅行路径的误差，提高了叠前时间偏移的成像效率和成像质量。所述方法及装置还为后续的例如时间域的地质解释、地质构造成图、深度域地质模型、初始深度域速度场、深度偏移等资料处理提供准确度可靠的成像资料。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本发明提供的叠前时间偏移的方法的一种实施例的方法流程图；

图2是本发明提供的从炮检点位置进行叠前时间偏移的走时图；

图3是本发明提供的叠前时间偏移的理论模型；

图4是本发明应用场景中对理论模型采用叠前深度偏移方法偏移得到的结果；

图5是本发明应用场景中采用实施例方法得到的偏移结果；

图6是本发明提供的叠前时间偏移装置的一种实施例的模块结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

由上述分析可知，由地表到统一水平基准面的校正往往会给速度分析的精度以及地震波的旅行路径造成误差，如果从炮检点的真实位置进行叠前时间偏移将省去校正的步骤，理论上可以减小由于校正而带来的误差。图1是本发明提供的叠前时间偏移方法的一种实施例的方法流程图，所述方法包括以下步骤：

步骤S1：从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集。

通过激发人工地震可以采集到大量的地震数据，所述地震数据蕴含着丰富的地层信息。在实际的采集过程中，由于经济和地理因素的制约，采集到的地震数据往往是有限的，通过有限的地震数据直观地显示大范围内的地层信息通常需要对地震样本数据作空间插值得到网格地震数据，例如生成的地质剖面图、等值线图等。可以在所述地质剖面图的任意位置选取一个成像点，可以根据最大炮检距、所述成像点所在地层的倾角或者深度确定所述成像点的偏移孔径范围。在本发明的一个实施例中，所述偏移孔径范围可以包括以所述成像点为中心、以偏移孔径为半边长的正方形区域，或者以所述成像点为中心、以偏移孔径为半径的圆形区域。那么，在所述偏移孔径范围内包含至少一个的网格地震数据，可以从所述网格地震数据中选取共中心点，所述偏移孔径范围可以包含一个或者多个共中心点。将不同炮集中拥有所述共中心点的记录道(也称为地震道)提取出来，可以形成所述共中心点的共中心点道集。

本发明的一个实施例中，在所述从地震数据中获取成像点的共中心点以及共中心点道集之前，还包括以下步骤：

步骤SS1：建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线距离、最小炮检距、最大炮检距。

所述观测系统的正确与否可以直接影响地震数据的采集质量、资料处理的精度。定义所述观测系统的参数可以建立合适的观测系统，所述参数包括：网格间距、炮线间距、最小炮检距。在本发明的一个实施例中，选取的网格间距大小适中，在本发明的观测系统中，所述网格间距可以包括10m×10m；所述炮线间距可以表示炮点线向前滚动的距离，在本发明的观测系统中，所述炮线间距可以包括80米到100米；较大的炮检距可以压制干扰波，但是，过大的炮检距可以不必要地压制反射波，本发明的观测系统中，所述最小炮检距和最大炮检距可以根据经验所得。

所述偏移孔径与所述最大炮检距的关系较为密切，在本发明的一个实施例中，所述偏移孔径的大小的范围可以包括：所述偏移孔径的大小可以为所述最大炮检距的两倍。

步骤SS2：根据所述观测系统采集地震数据。

在网格化后的地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集，为所述成像点的叠前时间偏移提供必要的数据基础。在获取地震数据之前定义准确的观测系统，有利于提高地震数据的采集质量和数据资料的处理精度。

步骤S2：根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时。

本发明中的共中心点道集至少包括一个记录道，所述记录道包括从炮点到共中心点的路径以及从共中心点到检波点的路径，可以根据记录道的地震数据计算所述记录道在所述成像点处的旅行时。本发明中的记录道的地震数据可以包括：炮点、检波点、成像点的坐标以及地表高程；成像点到地表的单程垂直时间；成像点处的速度。所述计算所述记录道在所述成像点处的旅行时的计算公式包括下式(1)：

式中，T_ij可以为所述成像点第i个共中心点道集中第j个记录道在成像点处的旅行时；t可以为成像点到地表的单程垂直时间，v可以为成像点处的速度；Eleve_cmpijs可以为第i个共中心点道集中第j个记录道炮点的地面高程；Eleve_cmpijr可以为第i个共中心点道集中第j个记录道检波点的地面高程；Eleve_crp可以为成像点处的地表高程；X_ijs可以为第i个共中心点道集中第j个记录道炮点到成像点的水平距离；X_ijr可以为第i个共中心点道集中第j个记录道检波点到成像点的水平距离。

下式(2)是现有技术中将炮检点从真实位置校正到统一水平基准面后，利用直线法计算记录道到成像点旅行时的计算公式：

式中，t_i可以为地震波从炮点到成像点以及从成像点到接收点的走时，t可以为成像点的单程垂直时间，v(t)可以为成像点的叠前时间偏移速度，x₁可以为炮点至成像点的水平距离，x₂可以为检波点至成像点的水平距离。对比式(1)和式(2)，可以发现式(1)相对于式(2)在成像点的单程垂直时间上增加了一项校正量，所述校正量随炮点、检波点、成像点到底面的高程的不同而不同。因此，式(1)相对式(2)多了三个参变量Elevecmp_ijs、Elevecmp_ijr、Eleve_crp，所述三个参变量可以体现成像点与参与成像的炮点、检波点真实位置的关系，可以实现从以炮检点真实位置为基准面的叠前时间偏移的分析。

举例说明，图2是本发明提供的从炮检点位置进行叠前时间偏移的走时图。如图2所示，CMP1、CMP2可以是以P为成像点的偏移孔径范围内的两个共中心点道集，S1、R1可以为CMP1中的一个记录道，S2、R2可以为CMP2中的一个记录道，其中，S1、S2可以为炮点位置，R1、R2可以为检波点位置。Crp可以是成像点P在以炮检点真实位置为基准面上的投影，t、v可以分别是从以炮检点真实位置为基准面到成像点处的单程时间和速度，S1、R1、S2、R2、Crp的地表高程分别为Eleve_cmp11s、Eleve_cmp11r、Eleve_cmp22s、Eleve_cmp22r、Eleve_crp。X_1s、X_1r、X_2s、X2r分别为成像点P与S1、R1、S2、R2的水平距离，过S1作水平线与线段P-Crp延长线交于P1，过R1、S2、R2作水平线与线段P-P1分别交于P2、P3、P4，S1、R1到成像点P的走时之和为T1，S2、R2到成像点P的走时之和为T2，那么T1、T2的计算公式分别如下式(3)、(4)所示：

步骤S3：提取所述记录道在所述旅行时处对应的振幅值。

对于式(1)而言，所述记录道在所述成像点处的旅行时T_ij后，可以提取在成像点偏移孔径范围内的第i个共中心点道集的第j个记录道在旅行时T_ij所对应的样点值，获取所述样点值的振幅值，至此可以完成第i个共中心点的第j个记录道的偏移。

步骤S4：根据所述振幅值计算所述共中心点道集的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

假设A_ij为第i个共中心点道集的第j个记录道的偏移振幅值，那么所述共中心点道集的振幅平均值可以表达为下式(5)：

式(5)中，A_p可以表示共中心点道集的振幅平均值，M可以表示共中心点或者共中心点道集的个数，N可以表示第i个共中心点道集中记录道的个数。可以将所述共中心点道集的振幅平均值A_p作为所述成像点的偏移振幅。

下面通过一个具体场景应用所述叠前时间偏移方法，图3是本发明提供的叠前时间偏移的理论模型。如图3所示，基准面位于X轴上，炮检点的位置曲线为z＝x+800，所述位置曲线的地层倾角为α＝45度，cmp1和cmp2为偏移孔径内的两个共中心点，成像点为Crp，成像点在地面的投影为P点，Crp与P点的距离为3000米，t0为成像点Crp与地面的单程吹响时间，v0为成像点Crp处的速度。从共中心点cmp1中抽取一个记录道S1R1，从共中心点cmp2中抽取一个记录道S2R2，这两个记录道的参数数据如下表1所示。

表1共中心点和成像点的参数数据表

根据表1的参数数据，应用本发明实施例中的公式(1)，计算成像点处的旅行时，并和理论计算结果进行对比，对比结果如表2所示。可以发现，采用公式(1)计算得到的成像点处的旅行时和理论计算的结果相同，因此，本发明实施例中计算成像点的旅行时精度较高。

表2本发明实施例方法计算成像点处旅行时和理论旅行时的结果对比表

下面通过另一个具体场景应用所述叠前时间偏移方法，某工区理论模型的炮检点的位置曲线为z＝-0.28x+2500，根据偏移孔径范围获取共中心点的范围可以包括120米到678米，共中心点的间距可以为25米。该工区的地层可以是均匀介质，介质速度为3000m/s。图4是本发明应用场景中对理论模型采用叠前深度偏移方法偏移得到的结果，图5是本发明应用场景中采用实施例方法得到的偏移结果。由于叠前深度偏移方法得到的成像结果可以表示真实的地质体，通过对图4和图5的对比发现，叠前深度偏移方法偏移得到的结果和本发明实施例得到的结果相匹配，进一步证实本申请成像精度的准确性。

本发明另一方面还提供一种叠前时间偏移的装置，图6是本发明提供的叠前时间偏移装置的一种实施例的模块结构示意图，结合附图6，装置60可以包括：道集获取单元61、旅行时计算单元62、振幅值提取单元63、偏移振幅获取单元64，其中，

道集获取单元61，用于从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集。

旅行时计算单元62，用于根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时。

振幅值提取单元63，用于提取所述记录道在所述旅行时处对应的振幅值。

偏移振幅获取单元64，用于根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

如图6所示，所述装置还可以包括：观测系统建立单元65和地震数据采集单元66，其中，

观测系统建立单元65，用于建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线间距、最小炮检距、最大炮检距。

地震数据采集单元66，用于根据所述观测系统采集地震数据。

由此可见，本发明一种叠前时间偏移的方法及装置的技术方案通过从成像点的实际位置处计算偏移振幅，避免了将地震数据校正到统一基准面上而带来的计算旅行路径的误差，提高了叠前时间偏移的成像效率和成像质量。所述方法及装置还为后续的例如时间域的地质解释、地质构造成图、深度域地质模型、初始深度域速度场、深度偏移等资料处理提供准确度可靠的成像资料。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

上述实施例阐明的装置或模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然，也可以将实现某功能的模块由多个子模块或子单元组合实现。

本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(ApplicationSpecificIntegratedCircuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC625D、AtmelAT91SAM、MicrochipPIC18F26K20以及SiliconeLabsC8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的硬件的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。

本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种叠前时间偏移的方法，其特征在于，所述方法包括：

从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集；

根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时；

提取所述记录道在所述旅行时处对应的振幅值；

根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

2.根据权利要求1所述的一种叠前时间偏移的方法，其特征在于，所述记录道的地震数据包括：

炮点、检波点、成像点的坐标、地表高程；

成像点到地表的单程垂直时间；

成像点处的速度。

3.根据权利要求2所述的一种叠前时间偏移的方法，其特征在于，所述计算所述记录道在所述成像点处的旅行时的计算公式包括下式：

$\begin{array}{c}{T}_{ij}=\sqrt{{(t+({\mathrm{Eleve}}_{cmpijs}-{\mathrm{Eleve}}_{crp})/v)}^2+\frac{{X_{ijs}}^2}{v^2}}+\\ \sqrt{{(t+({\mathrm{Eleve}}_{cmpijr}-{\mathrm{Eleve}}_{crp})/v)}^2+\frac{{X_{ijr}}^2}{v^2}}\end{array}$

式中，T_ij为所述成像点第i个共中心点道集中第j个记录道在成像点处的旅行时；t为成像点到地表的单程垂直时间，v为成像点处的速度；Eleve_cmpijs为第i个共中心点道集中第j个记录道炮点的地面高程；Eleve_cmpijr为第i个共中心点道集中第j个记录道检波点的地面高程；Eleve_crp可以为成像点处的地表高程；X_ijs为第i个共中心点道集中第j个记录道炮点到成像点的水平距离；X_ijr为第i个共中心点道集中第j个记录道检波点到成像点的水平距离。

4.根据权利要求1所述的一种叠前时间偏移的方法，其特征在于，在所述从地震数据中获取成像点的共中心点以及共中心点道集之前，还包括：

建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线间距、最小炮检距、最大炮检距；

根据所述观测系统采集地震数据。

5.根据权利要求4所述的一种叠前时间偏移的方法，其特征在于，所述偏移孔径的大小的范围包括：所述偏移孔径的大小为所述最大炮检距的两倍。

6.一种叠前时间偏移的装置，其特征在于，所述装置包括：

道集获取单元，用于从地震数据中获取成像点在偏移孔径范围内的共中心点以及共中心点道集；

旅行时计算单元，用于根据所述共中心点道集中记录道的地震数据计算所述记录道在所述成像点处的旅行时；

振幅值提取单元，用于提取所述记录道在所述旅行时处对应的振幅值；

偏移振幅获取单元，用于根据所述振幅值计算所述共中心点的振幅平均值，将所述振幅平均值作为所述成像点的偏移振幅。

7.根据权利要求6所述的一种叠前时间偏移的装置，其特征在于，所述记录道的地震数据包括：

炮点、检波点、成像点的坐标、地表高程；

成像点到地表的单程垂直时间；

成像点处的速度。

8.根据权利要求7所述的一种叠前时间偏移的装置，其特征在于，所述计算所述记录道在所述成像点处的旅行时的计算公式包括下式：

$\begin{array}{c}{T}_{ij}=\sqrt{{(t+({\mathrm{Eleve}}_{cmpijs}-{\mathrm{Eleve}}_{crp})/v)}^2+\frac{{X_{ijs}}^2}{v^2}}+\\ \sqrt{{(t+({\mathrm{Eleve}}_{cmpijr}-{\mathrm{Eleve}}_{crp})/v)}^2+\frac{{X_{ijr}}^2}{v^2}}\end{array}$

式中，T_ij为所述成像点第i个共中心点道集中第j个记录道在成像点处的旅行时；t为成像点到地表的单程垂直时间，v为成像点处的速度；Eleve_cmpijs为第i个共中心点道集中第j个记录道炮点的地面高程；Eleve_cmpijr为第i个共中心点道集中第j个记录道检波点的地面高程；Eleve_crp为成像点处的地表高程；X_ijs为第i个共中心点道集中第j个记录道炮点到成像点的水平距离；X_ijr为第i个共中心点道集中第j个记录道检波点到成像点的水平距离。

9.根据权利要求6所述的一种叠前时间偏移的装置，其特征在于，所述装置还包括：

观测系统建立单元，用于建立观测系统，定义所述观测系统的参数，所述参数包括：网格间距、炮线间距、最小炮检距、最大炮检距；

地震数据采集单元，用于根据所述观测系统采集地震数据。

10.根据权利要求9所述的一种叠前时间偏移的装置，其特征在于，所述偏移孔径的大小的范围包括：所述偏移孔径的大小为所述最大炮检距的两倍。