逆时偏移中的假像去除的制作方法

优先权声明

本申请要求于2017年9月11日递交的美国专利申请no.15/701,038的优先权，其全部内容通过引用并入本文。

本公开涉及地震数据处理，并且更具体地涉及逆时偏移中的假像去除。

背景技术：

地球地下层的高分辨率深度图像对于油藏勘探、圈定和开发很重要。基于波动方程的地震深度偏移技术(例如，逆时偏移(rtm))是石油行业中(尤其是在诸如盐下勘探之类的复杂环境中)用于地震深度成像的合适的技术。然而，由于沿射线路径的非反射点处的源波场和接收机波场的不期望的互相关，rtm中存在成像噪声。

技术实现要素：

本公开描述了用于去除逆时偏移(rtm)中的假象的方法和系统，包括计算机实现的方法、计算机程序产品和计算机系统。一种计算机实现的方法包括：将源波场和接收机波场的体积划分为多个窗口；对于多个窗口中的每个窗口并且在特定时间值处，计算源波场的下行分量和接收机波场的上行分量，源波场的下行分量是源波场的傅里叶变换的正分量与源波场的傅里叶变换的正分量的希尔伯特变换之和的傅里叶逆变换的实分量，接收机波场的上行分量是接收机波场的傅里叶变换的负分量与接收机波场的傅里叶变换的负分量的希尔伯特变换之和的傅里叶逆变换的实分量；通过将多个窗口的每个窗口中的源波场的下行分量相加来计算特定时间值处的源波场的下行分量；以及通过将多个窗口的每个窗口中的接收机波场的上行分量相加来计算特定时间值处的接收机波场的上行分量。

这个方面的其他实现方式包括对应的计算机系统、装置和记录在一个或多个计算机存储设备上的计算机程序，每一个均被配置为执行方法的动作。一个或多个计算机的系统可以被配置为通过在系统上安装的、在操作中使得系统执行动作的软件、固件、硬件或软件、固件或硬件的组合来执行特定的操作或动作。一个或多个计算机程序可以被配置为通过包括指令来执行特定操作或动作，所述指令在被数据处理装置执行时使得该装置执行动作。

前述和其他实现方式可以各自可选地以单独或组合的方式包括以下特征中的一个或多个：

与总的实现方式可组合的第一方面，包括：计算源波场的时间导数的傅里叶变换的正分量，以及计算接收机波场的时间导数的傅里叶变换的负分量。

与前述方面中的任何方面可组合的第二方面，其中，多个窗口是多个重叠的窗口。

与前述方面中的任何方面可组合的第三方面，其中，将每个窗口中的源波场的下行分量相加以及将每个窗口中的接收机波场的上行分量相加包括：以重叠区域消减的方式将每个窗口中的源波场的下行分量相加以及将每个窗口中的接收机波场的上行分量相加。

与前述方面中的任何方面可组合的第四方面，其中，傅里叶变换是1d空间傅里叶变换。

与前述方面中的任何方面可组合的第五方面，包括：基于特定时间值处的源波场的下行分量和接收机波场的上行分量来计算图像项。

本说明书的主题的一个或多个实现方式的细节在附图和以下描述中阐述。通过说明书、附图和权利要求书，所述主题的其他特征、方面和优点将变得显而易见。

附图说明

图1示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法。

图2示出了根据一些实现方式的用于基于源波场的下行分量和接收机波场的上行分量来计算图像的示例方法。

图3示出了根据一些实现方式的在去除逆时偏移(rtm)中的假象之前和之后的示例图像。

图4是示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法的流程图。

图5是示出了根据一些实现方式的用于提供与如本公开中描述的所述算法、方法、功能、处理、流程和过程相关联的计算功能的示例计算机系统的框图。

各附图中相似的附图标号和标记指示相似的元素。

具体实施方式

以下详细描述描述了去除由在具有错误传播方向组合的源波场和接收机波场上施加了逆时偏移(rtm)图像条件引起的假像，并提出以下详细描述以使本领域技术人员能够制造和使用一个或多个特定实现方式的上下文中所公开的主题。可以对所公开的实现方式进行各种修改、变更和置换，这对于本领域技术人员而言将是显而易见的，并且在不脱离本公开的范围的情况下，所定义的总的原理可以应用于其他实现方式和应用。因此，本公开不旨在限于所描述或图示的实现方式，而是应被赋予与所公开的原理和特征一致的最宽泛范围。

地球地下层的高分辨率深度图像对于成功进行油藏勘探、圈定和开发很重要。基于波动方程的地震深度偏移技术(例如rtm)是石油行业中(尤其是在诸如盐下勘探之类的复杂环境中)用于地震深度成像的合适的技术。然而，由于沿射线路径的非反射点处的源波场和接收机波场的不期望的互相关，rtm中存在成像噪声。

在rtm中，通过应用常规图像条件来获得图像

其中，是源波场，是接收机波场，表示位置向量，t表示时间并且范围在0和tmax之间(例如，最大记录时间)。通常，根据liuetal.，2007，图像可以分解为四个分量：

其中，idu是下行源波场和上行接收机波场的互相关，idd是下行源波场和下行接收机波场的互相关，iuu是上行源波场和上行接收机波场的互相关，iud是上行源波场和下行接收机波场的互相关。另外，四个分量表示为：

其中，和分别是源波场和接收机波场的下行分量，和分别是源波场和接收机波场的上行分量。

可以通过从图像中消除和分量来抑制rtm中的低频成像噪声。为了计算剩余的图像可以沿着深度坐标z对波场应用希尔伯特变换，之后可以计算图像融合-组合(imageintegration-combination)。分量也可以形成假像，也需要将其消除。结果是，仅在图像中保留了分量，以便去除rtm中的假像。为了仅保留分量，将源波场或接收机波场中的至少一个波场拆分为下行分量和上行分量。因此，先前描述的希尔伯特变换可以不产生图像

将波场拆分为下行分量和上行分量的一种方法是在频率-波数域内执行拆分。例如，可以(例如，经由2d傅里叶变换)将源波场变换为其中，ω表示频率，表示波数向量。或中的下行分量是ωkz＞0的部分，而域中的上行分量是ωkz＜0的部分。然而，在rtm的实际实现方式中，通常通过有限差分方法来计算波场，并且由于数据量大，因此仅能存储波场的几个时间片。因此，不能执行从时间t到频率ω的波场的傅里叶变换、或用于挑选正/负ω部分的希尔伯特变换。

从高层面来看，所描述的方法提供了一种将波场拆分为下行分量和上行分量的近似方法。所描述的近似方法仅利用例如波场的当前时间片和当前时间导数片的知识将波场拆分为下行分量和上行分量。可以根据rtm中的波场的有限差分正演建模获取波场的当前时间片和当前时间导数片。基于该拆分，可以计算图像并且可以从中消除由错误分量组合(例如，)产生的假像，以增强rtm的最终图像。结果是，rmt的最终图像可以用于例如更好地估计地球地下层的特性，以用于石油和天然气勘探。

图1示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法100。例如，基于rtm正演建模，源波场在t＝ti处的时间片为接收机波场在t＝ti处的时间片为另外，可以获得t＝ti处的源波场的时间导数片和接收机波场的时间导数片，如下：

利用波场的时间片和时间导数片，可以通过图1中所示的工作流将波场拆分为下行分量和上行分量。

如图1所示，示例方法100开始于将波场(例如，源波场和/或接收机波场)的体积划分为多个重叠的窗口。例如，该窗口可以在3d空间中定义，并提供有限的3d子体积。在3d空间中的每个维度中，窗口与相邻的窗口重叠一半。通常，窗口大小可以是例如其维度大小的十分之一。在本公开中，假设窗口内介质的速度是近似均匀的，并且对于所有波数向量波长的频散关系为：

其中，v是窗口内介质的平均速度，ω是圆频率。

对于多个重叠窗口中的每个窗口，对源波场的时间片(即)、源波场的时间导数片(即)、接收机波场的时间片(即)和接收机波场的时间导数片(即)应用局部空间傅里叶变换(例如，一维傅里叶变换)。所得到的傅里叶变换分别为和可以获得和的kz＞0的部分，以及和的kz＜0的部分，如下：

如果仅对应于例如两个圆频率(如等式(9)中所示的一个正的和一个负的)，则等式(10)中的可以分解为如下两项：

基于傅里叶变换的特性，可以获得下式：

其中，j是虚数单位。类似地，等式(12)中的也可以被分解为如下两项：

基于傅里叶变换的特性，可以获得下式：

根据等式(15)和(17)，可以近似导出和的希尔伯特变换如下：

其中，是波数向量的范数，ε为用于避免等式中除以0而选择的用户指定的较小的数。例如，ε的值可以是0.0001或更小。

根据等式(14)和(18)，可以获得下式(即，具有kz＞0且ω＞0的分量)：

通过傅立叶变换的共轭特性，可以获得特定窗口的下行源波场：

其中，ifft()是从到的空间傅里叶逆变换，re()是返回实部的函数。

类似地，根据等式(16)和(19)，可以获得下式(即，具有kz＜0且ω＜0的分量)：

通过傅立叶变换的共轭特性，可以获得特定窗口的上行接收机波场：

为了获得整个下行源波场将的所有窗口结果(例如，针对多个重叠窗口中的每一个窗口经由等式(21)计算的结果)相加。类似地，为了获得整个上行接收机波场将的所有窗口结果(例如，针对多个重叠窗口中的每一个窗口经由等式(23)计算的结果)相加。如上所述，图1提供了计算源波场的下行分量和接收机波场的上行分量的近似方法。如等式(18)至(20)和(22)中所示，在一些实现方式中，仅利用波场的时间片和波场的时间导数片拆分波场的关键是选择ω＞0和/或ω＜0分量，而不在时域中使用傅里叶变换。

图2示出了根据一些实现方式的用于基于源波场的下行分量和接收机波场的上行分量来计算图像的示例方法200。例如，示例方法200开始于框205，其中，i和图像均被初始化为0。在框210处，从源波场和接收机波场(即和)获得时间ti处的源波场和接收机波场的时间片。另外，还(例如，经由等式(7)和(8))获得时间ti处的源波场和接收机波场的时间导数片。在框215处，可以使用图1中所描述的近似方法来计算整个下行源波场(即)和整个上行接收机波场(即)。在框220处，将添加到图像在框225处，i增加1。在框230处，确定ti是否大于tmax。tmax是例如源波场和接收机波场的最大记录时间。如果ti不大于tmax，则方法进行到框210。如果ti大于tmax，则在框235处输出最终的图像尽管本公开参考图2中的方法，但这是出于示例的目的，该文档的主题可以使用其他方法，基于使用图1中描述的近似方法计算的下行源波场和上行接收机波场来计算图像

图3示出了根据一些实现方式的在去除逆时偏移(rtm)中的假象之前和之后的示例图像300。在图3中，在顶部示出了使用liuet.al.，2007提出的方法的rtm脉冲响应结果305。结果305示出了偏移噪声(例如，从和生成的)以及由三个箭头所指的假像(例如，从生成的)。在图3的底部示出了使用所描述的近似方法的结果310。结果310示出了偏移噪声和假像的去除。

图4是示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法400的流程图。为了说明的清楚，下面的描述总体上描述了在本说明书中的其他附图的上下文下的方法400。例如，方法400可以由图5中所描述的计算机系统执行。然而，应理解，方法400可以视情况而定例如由任何合适的系统、环境、软件和硬件执行，或由系统、环境、软件和硬件的组合执行。在一些实现方式中，方法400的各个步骤可以并行、组合、循环或以任意顺序运行。

方法400开始于框405，其中将源波场和接收机波场的体积划分为多个窗口。在一些实现方式中，多个窗口是多个重叠的窗口。对于多个窗口中的每个窗口且在特定的时间值处(例如，0和tmax之间的ti)，执行下面的操作(即框410至框425)。

在框410处，计算源波场的傅里叶变换的正分量(例如，等式(10))和接收机波场的傅里叶变换的负分量(例如，等式(12))。在一些实现方式中，傅里叶变换是在时域中不执行的1d空间傅里叶变换。在一些实现方式中，正分量是kz＞0的部分，负分量是kz＜0的部分。

在框415处，计算源波场的傅里叶变换的正分量的希尔伯特变换和接收机波场的傅里叶变换的负分量的希尔伯特变换。在一些实现方式中，基于源波场的时间导数的傅里叶变换的正分量来计算源波场的傅里叶变换的正分量的希尔伯特变换(例如，等式(18))。在一些实现方式中，基于接收机波场的时间导数的傅里叶变换的负分量来计算接收机波场的傅里叶变换的负分量的希尔伯特变换(例如，等式(19))。

在框420处，计算特定窗口中的源波场的下行分量(例如，)。在一些实现方式中，首先计算源波场的傅里叶变换的正分量和源波场的傅里叶变换的正分量的希尔伯特变换之和(例如，等式(20))。接下来，对该和应用傅里叶逆变换。然后，将特定窗口中的源波场的下行分量计算为和的傅立叶逆变换的实分量(例如，等式(21))。

在框425处，计算特定窗口中的接收机波场的上行分量(例如，)。在一些实现方式中，首先计算接收机波场的傅里叶变换的负分量和接收机波场的傅里叶变换的负分量的希尔伯特变换之和(例如，等式(22))。接下来，对该和应用傅里叶逆变换。然后，将特定窗口中的接收机波场的上行分量计算为和的傅立叶逆变换的实分量(例如，等式(23))。

在框430处，计算特定时间值处的源波场的下行分量和接收机波场的上行分量(例如，和)。通过将多个窗口的每个窗口中的源波场的下行分量(例如，在等式(21)中计算的)相加来计算特定时间值处的源波场的下行分量。通过将多个窗口的每个窗口中的接收机波场的上行分量(例如，在等式(23)中计算的)相加来计算特定时间值处的接收机波场的上行分量。在一些实现方式中，以重叠区域消减(tapering)的方式将多个窗口的每个窗口中的源波场的下行分量相加。在一些实现方式中，以重叠区域消减的方式将多个窗口的每个窗口中的接收机波场的上行分量相加。

可以修改或重新配置图4所示的示例方法400，以包括附加的、更少的或不同的步骤(图4中未示出)，这些步骤可以按所示顺序或以不同顺序执行。例如，在框430之后，可以基于特定时间值处的源波场的下行分量和接收机波场的上行分量来计算图像(例如，)。例如，在一些实现方式中，可以重复或迭代图4中所示的步骤，直到达到终止条件为止。在一些实现方式中，图4中所示的各个步骤中的一个或多个步骤可以作为多个分离的步骤来执行，或者图4中所示的步骤的一个或多个子集可以组合并且作为单个步骤来执行。在一些实现方式中，也可以从示例方法400中省略图4所示的各个步骤中的一个或多个步骤。

图5是示出了根据实现方式的用于提供与如本公开中描述的所述算法、方法、功能、处理、流程和过程相关联的计算功能的示例计算机系统500的框图。所示出的计算机502旨在包括任意计算设备，例如服务器、台式计算机、膝上型/笔记本计算机、无线数据端口、智能电话、个人数字助理(pda)、平板计算设备、这些设备内的一个或多个处理器、或任意其他合适的处理设备(包括计算设备的物理或虚拟实例(或这两者))。附加地，计算机502可以包括以下这样的计算机，该计算机包括可以接受用户信息的输入设备(例如键区、键盘、触摸屏或其他设备)以及输出设备，该输出设备传达与计算机502的操作相关联的信息，包括数字数据、视觉或音频信息(或信息的组合)或图形用户界面(gui)。

计算机502可以用作用于执行本公开中描述的主题的计算机系统的客户端、网络组件、服务器、数据库或其他持久性或任意其他组件(或它们的组合)。所示出的计算机502可通信地与网络530耦接。在一些实现方式中，计算机502的一个或多个组件可以被配置为在包括基于云计算、局部、全局、或其他环境在内的环境(或者环境的组合)中操作。

从高层面来看，计算机502是可操作用于接收、发送、处理、存储或管理与所描述的主题相关联的数据和信息的电子计算设备。根据一些实现方式，计算机502还可以包括或可通信地耦接到应用服务器、电子邮件服务器、web服务器、缓存服务器、流传输数据服务器或其他服务器(或服务器的组合)。

计算机502可以通过网络530从客户端应用(例如，在另一计算机上执行的应用)接收请求，并通过使用适当的软件应用处理所接收的请求来响应所述请求。此外，还可以从内部用户(例如，从命令控制台或通过其他适当的访问方法)、外部或第三方、其他自动化应用以及任何其他适当的实体、个人、系统或计算机向计算机502发送请求。

计算机502的组件中的每个组件可以使用系统总线503通信。在一些实现方式中，计算机502的任意或所有组件(硬件和/或软件(或硬件和软件的组合))可以使用应用编程接口(api)512或服务层513(或api512和服务层513的组合)，通过系统总线503彼此交互或与接口504交互(或两者的组合)。api512可以包括针对例程、数据结构和对象类的规范。api512可以是独立于或依赖于计算机语言的，并且可以指代完整的接口、单个功能或甚至是一组api。服务层513向计算机502或可通信地耦接到计算机502的其他组件(无论是否被示出)提供软件服务。计算机502的功能可以是对于使用该服务层的所有服务消费者是可访问的。软件服务(例如由服务层513提供的软件服务)通过定义的接口提供可重用的、定义的业务功能。例如，接口可以是以java、c++或以可扩展标记语言(xml)格式或其他合适格式提供数据的其他合适语言编写的软件。尽管被示为计算机502的集成组件，但是备选实现方式可以将api512或服务层513示作为相对于计算机502的其他组件或可通信地耦接到计算机502的其他组件(无论是否被示出)的独立组件。此外，在不脱离本公开的范围的情况下，api512或服务层513的任意或所有部分可以被实现为另一软件模块、企业应用或硬件模块的子模块或副模块。

计算机502包括接口504。尽管在图5中被示为单个接口504，但是可以根据计算机502的特定需要、期望或特定实现方式而使用两个或更多个接口504。接口504由计算机502使用以与分布式环境中的连接到网络530的其他系统(无论是否示出)通信。通常，接口504包括以软件或硬件(或软件和硬件的组合)编码的逻辑，并且可操作用于与网络530通信。更具体地，接口504可以包括支持与通信相关联的一个或多个通信协议的软件，使得网络530或接口的硬件可操作用于在所示出的计算机502内部和外部传送物理信号。

计算机502包括处理器505。尽管在图5中被示出为单个处理器505，但是可以根据计算机502的特定需要、期望或特定实现方式而使用两个或更多个处理器。通常，处理器505执行指令并操纵数据，以执行计算机502的操作以及如本公开中所描述的任何算法、方法、功能、处理、流程和过程。

计算机502还包括数据库506，其可以保存用于计算机502或可以连接到网络530的其他组件(无论是否被示出)(或两者的组合)的数据。例如，数据库506可以是存储与本公开一致的数据的内部存储器、常规或其他类型的数据库。在一些实现方式中，根据计算机502的特定需要、期望或特定实现方式和所描述的功能，数据库506可以是两个或更多个不同数据库类型(例如，混合的内部存储器和常规数据库)的组合。尽管在图5中被示出为单个数据库506，但是可以根据计算机502的特定需要、期望或特定实现方式和所描述的功能而使用(相同类型或多个类型的组合的)两个或更多个数据库。尽管数据库506被示出为计算机502的集成组件，但是在备选实现方式中，数据库506可以在计算机502的外部。如图所示，数据库506保存地震数据516、源波场和接收机波场518、源波场和接收机波场的时间导数片520以及下行源波场和上行接收机波场522。

计算机502还包括存储器507，其保存用于计算机502或可以连接到网络530的其他组件(无论是否被示出)(或两者的组合)的数据。例如，存储器507可以是存储与本公开一致的数据的随机存取存储器(ram)、只读存储器(rom)、光、磁等。在一些实现方式中，根据计算机502的特定需要、期望或特定实现方式和所描述的功能，存储器507可以是两个或更多个不同类型的存储器的组合(例如，ram和磁存储器的组合)。尽管在图5中被示出为单个存储器507，但是可以根据计算机502的特定需要、期望或特定实现方式和所描述的功能而使用(具有相同类型或多个类型组合的)两个或更多个存储器507。尽管存储器507被示出为计算机502的集成组件，但是在备选实现方式中，存储器507可以在计算机502的外部。

应用508是根据计算机502的特定需要、期望或特定实现方式提供功能(特别是关于本公开中描述的功能)的算法软件引擎。例如，应用508可以用作一个或多个组件、模块、应用等。此外，尽管被示出为单个应用508，但是应用508可以被实现为计算机502上的多个应用508。另外，尽管被示出为与计算机502集成在一起，但是在备选实现方式中，应用508可以在计算机502的外部。

可以存在与包含计算机502的计算机系统相关联或在其外部的任意数量的计算机502，每个计算机502通过网络530进行通信。此外，在不脱离本公开的范围的情况下，术语“客户端”、“用户”和其他适当的术语可以视情况而定互换使用。此外，本公开包含许多用户可以使用一个计算机502，或者一个用户可以使用多个计算机502。

在本说明书中描述的主题和功能操作的实现可以在数字电子电路中、在有形地被实现的计算机软件或固件中、在计算机硬件(包括在本说明书中公开的结构及其结构等同物)中、或在它们中的一个或多个的组合中实施。在本说明书中描述的主题的实现可以被实现为一个或多个计算机程序，即，在有形的非暂时性计算机可读计算机存储介质上编码的计算机程序指令的一个或多个模块，以用于被数据处理装置执行或控制数据处理装置的操作。备选地或附加地，程序指令可以编码在人工生成的传播信号(例如，机器生成的电、光或电磁信号)中/上，所述信号被生成以对信息进行编码，以传输给合适的接收机装置，以供数据处理装置执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或计算机存储介质的组合。

术语“数据处理装置”、“计算机”或“电子计算机设备”(或本领域普通技术人员所理解的等同物)指代数据处理硬件，并且包括用于处理数据的各种装置、设备和机器，例如，包括可编程处理器、计算机、或多个处理器或计算机。所述装置还可以是或还包括专用逻辑电路，例如，中央处理单元(cpu)、fpga(现场可编程门阵列)或asic(专用集成电路)。在一些实现方式中，数据处理装置或专用逻辑电路(或数据处理装置或专用逻辑电路的组合)可以是基于硬件或基于软件(或基于硬件和基于软件的组合)的。可选地，装置可以包括为计算机程序创建执行环境的代码，例如，构成处理器固件、协议栈、数据库管理系统、操作系统或者执行环境的组合的代码。本公开考虑具有或不具有常规操作系统(例如，linux、unix、windows、macos、android、ios或任意其他合适的常规操作系统)的数据处理装置的使用。

可以以任何形式的编程语言来写计算机程序(也可以称作或描述为程序、软件、软件应用、模块、软件模块、脚本或代码)，所述编程语言包括：编译或解译语言、或者声明或程序语言，并且可以以任何形式来部署计算机程序，包括部署为单独的程序或者部署为适合于用于计算环境的模块、组件、子例程、或者其他单元。计算机程序可以但不是必须与文件系统中的文件相对应。程序可以被存储在保存其他程序或数据(例如，存储在标记语言文档中的一个或多个脚本)的文件的一部分中、被存储在专用于所讨论的程序的单个文件中、或者被存储在多个协同文件(例如，存储一个或多个模块、子程序或代码部分的文件)中。计算机程序可以被部署为在一个计算机上或者在位于一个站点或分布在多个站点并且通过通信网络互联的多个计算机上执行。尽管各图中所出示的程序的部分被示出为通过各种对象、方法或其他过程实施各种特征和功能的各个模块，但是视情况而定程序可以替代地包括多个子模块、第三方服务、组件、库等。相反，各种组件的特征和功能可以视情况而定组合成单个组件。可以统计地、动态地或者统计地且动态地确定用于进行计算确定的阈值。

本说明书中描述的方法、处理或逻辑流可以由一个或多个可编程计算机来执行，所述一个或多个可编程计算机执行一个或多个计算机程序以通过操作输入数据并且生成输出来执行功能。方法、处理或逻辑流也可以由专用逻辑电路(例如cpu、fpga或asic)来执行，并且装置也可以被实现为专用逻辑电路(例如cpu、fpga或asic)。

适合于执行计算机程序的计算机可以基于通用或专用微处理器、这两者或任何其他类型的cpu。通常，cpu将从只读存储器(rom)或随机存取存储器(ram)或者这二者接收指令和数据。计算机的必不可少的元件是用于执行指令的cpu和用于存储指令和数据的一个或更多个存储器设备。通常，计算机还将包括用于存储数据的一个或多个大容量存储设备(例如，磁盘、磁光盘或光盘)，或可操作地耦接以便从所述一个或多个大容量存储设备接收或向其发送数据或两者。然而，计算机不必非要具有这些设备。此外，计算机可以嵌入在另一设备中，例如，移动电话、个人数字助理(pda)、移动音频或视频播放器、游戏机、全球定位系统(gps)接收机或者便携式存储设备(例如，通用串行总线(usb)闪存驱动器)，这仅是举几个例子。

适合于存储计算机程序指令和数据的计算机可读介质(暂时或非暂时的，视情况而定)包括所有形式的非易失性存储器、介质和存储器设备，包括例如半导体存储器设备，例如可擦除可编程只读存储器(eprom)、电可擦除可编程只读存储器(eeprom)和闪存设备；磁盘(例如内部硬盘或可移动盘)；磁光盘；以及cdrom、dvd+/-r、dvd-ram和dvd-rom盘。存储器可以存储各种对象或数据，包括：高速缓存区、类(class)、框架、应用、备份数据、工作、网页、网页模板、数据库表格、存储动态信息的资源库、以及包括任意参数、变量、算法、指令、规则、约束、对其的引用在内的任意其他适当的信息。此外，存储器可以包括任何其他适当的数据，例如，日志、策略、安全或访问数据、报告文件等等。处理器和存储器可以由专用逻辑电路来补充或者并入到专用逻辑电路中。

为了提供与用户的交互，本说明书中描述的主题的实现可以实现在计算机上，该计算机具有用于向用户显示信息的显示设备(例如，crt(阴极射线管)、lcd(液晶显示器)、led(发光二极管)或等离子监视器)和用户可以向计算机提供输入的键盘和指点设备(例如，鼠标、轨迹球或轨迹板)。还可以使用触摸屏(例如，具有压敏性的平板计算机表面、使用电容或电感测的多点触摸屏或其他类型的触摸屏)向计算机提供输入。其他种类的设备也可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈，例如，视觉反馈、听觉反馈或触觉反馈；以及可以以任意形式(包括声音、语音或触觉输入)来接收来自用户的输入。此外，计算机可以通过向用户使用的设备发送文档或者从该设备接收文档来与用户交互；例如，通过响应于从用户客户端设备上的web浏览器接收到的请求而向所述web浏览器发送网页。

术语“图形用户界面”或“gui”可以以单数或复数形式使用，以描述一个或多个图形用户界面和特定图形用户界面的每一次显示。因此，gui可以表示任意图形用户界面，包括但不限于web浏览器、触摸屏或处理信息并且高效地地向用户呈现信息结果的命令行界面(cli)。通常，gui可以包括多个用户界面(ui)元素，其中一些或全部元素与web浏览器相关联，诸如交互式字段、下拉列表和按钮。这些和其他ui元素可以与web浏览器的功能相关或表示web浏览器的功能。

本说明书中描述的主题的实现可以在计算系统中实现，该计算系统包括后端组件(例如，数据服务器)、或包括中间件组件(例如，应用服务器)、或者包括前端组件(例如，具有用户通过其可以与本说明书中描述的主题的实现进行交互的图形用户界面或者web浏览器的客户端计算机)、或者一个或多个此类后端组件、中间件组件或前端组件的任意组合。系统的组件可以通过任何形式的有线或无线数字数据通信(或数据通信的组合)或者有线或无线数字数据通信的任何介质(例如通信网络)互连。通信网络的示例包括局域网(lan)、无线电接入网(ran)、城域网(man)、广域网(wan)、全球微波接入互操作性(wimax)、使用例如802.11a/b/g/n或802.20(或802.11x和802.20的组合或与本公开一致的其他协议)的无线局域网(wlan)、互联网的全部或一部分、或一个或多个位置处的任意其他通信系统(或通信网络的组合)。网络可以在网络地址之间传输例如互联网协议(ip)分组、帧中继帧、异步传输模式(atm)小区、语音、视频、数据或其他合适信息(或通信类型的组合)。

计算系统可以包括客户端和服务器。客户端和服务器一般相互远离并且通常通过通信网络进行交互。客户端和服务器的关系通过在相应计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生。

尽管本说明书包含许多特定实现细节，然而这些细节不应被解释为对可以要求保护的范围或任何发明的范围的限制，而是作为可以专用于特定发明的特定实现方式的特征的描述。在单个实现方式中，还可以组合实现本说明书中在独立实现方式的上下文中描述的某些特征。反之，在单个实现方式的上下文中描述的各种特征也可以在多个实现方式中分开地或以任何合适的子组合来实现。此外，尽管可以将先前描述的特征描述为以某些组合起作用并且甚至最初要求如此保护，但是在一些情况下，来自所要求保护的组合的一个或多个特征可以从组合中删除，并且所要求保护的组合可以涉及子组合或子组合的变型。

已经描述了本主题的特定实现方式。对于本领域技术人员显而易见的是，所描述的实现方式的其他实现方式、改变和置换在下文的权利要求的范围内。尽管在附图或权利要求中以特定顺序描述了操作，但这不应被理解为：为了实现希望的结果，要求按所示出的特定顺序或按相继的顺序来执行这些操作，或者要求执行所有所示出的操作(一些操作可以被认为是可选的)。在一些情况下，多任务或并行处理(或者多任务和并行处理的组合)可以是有利的并且视情况来执行。

此外，在前述的实现中的各种系统模块和组件的分离或集成不应被理解为在所有实现方式中要求这样的分离或集成，并且应该理解的是，所描述的程序组件和系统一般可以一起集成在单个软件产品中或封装为多个软件产品。

因此，先前描述的示例实现方式不限定或限制本公开。在不脱离本公开的精神和范围的情况下，还可以有其他改变、替换和变化。

此外，任何要求保护的实现方式被认为适用于至少一种计算机实现的方法；存储用于执行计算机实现的方法的计算机可读指令的非暂时性计算机可读介质；以及计算机系统，该系统包括与硬件处理器可互操作地耦接的计算机存储器，所述硬件处理器被配置为执行计算机实现的方法或存储在非暂时性计算机可读介质上的指令。

参考文献

faqiliu，guanquanzhang，scotta.mortonandjacquesp.leveille，2011，aneffectiveimagingconditionforreverse-timemigrationusingwavefielddecomposition：geophysics，vol.76，no.1，s29-s39.

faqiliu，guanquanzhang，scotta.morton，jacquesp.leveille，2007，reverse-timemigrationusingone-waywavefieldimagingcondition：77thannualinternationalmeeting，segexpandedabstracts.

tongw.fei，yiluoandfuhaoqin，2014，anendemicprobleminreverse-timemigration：84thannualinternationalmeeting，segexpandedabstracts。