CN111971587B - 逆时偏移中的假像去除 - Google Patents

逆时偏移中的假像去除 Download PDF

Info

Publication number
CN111971587B
CN111971587B CN201880072723.9A CN201880072723A CN111971587B CN 111971587 B CN111971587 B CN 111971587B CN 201880072723 A CN201880072723 A CN 201880072723A CN 111971587 B CN111971587 B CN 111971587B
Authority
CN
China
Prior art keywords
wavefield
fourier transform
component
receiver
source
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201880072723.9A
Other languages
English (en)
Other versions
CN111971587A (zh
Inventor
杨佳瑞
费同
骆毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Saudi Arabian Oil Co
Original Assignee
Saudi Arabian Oil Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Saudi Arabian Oil Co filed Critical Saudi Arabian Oil Co
Publication of CN111971587A publication Critical patent/CN111971587A/zh
Application granted granted Critical
Publication of CN111971587B publication Critical patent/CN111971587B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32—Transforming one recording into another or one representation into another
    • G01V1/325—Transforming one representation into another
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32—Transforming one recording into another or one representation into another
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/44—F-k domain
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/48—Other transforms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling
    • G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics (AREA)
  • Image Processing (AREA)
  • Complex Calculations (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • User Interface Of Digital Computer (AREA)
  • Image Analysis (AREA)

Abstract

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

Description

逆时偏移中的假像去除
优先权声明
本申请要求于2017年9月11日递交的美国专利申请No.15/701,038的优先权，其全部内容通过引用并入本文。
技术领域
本公开涉及地震数据处理，并且更具体地涉及逆时偏移中的假像去除。
背景技术
地球地下层的高分辨率深度图像对于油藏勘探、圈定和开发很重要。基于波动方程的地震深度偏移技术(例如，逆时偏移(RTM))是石油行业中(尤其是在诸如盐下勘探之类的复杂环境中)用于地震深度成像的合适的技术。然而，由于沿射线路径的非反射点处的源波场和接收机波场的不期望的互相关，RTM中存在成像噪声。
发明内容
本公开描述了用于去除逆时偏移(RTM)中的假象的方法和系统，包括计算机实现的方法、计算机程序产品和计算机系统。一种计算机实现的方法包括：将源波场和接收机波场的体积划分为多个窗口；对于多个窗口中的每个窗口并且在特定时间值处，计算源波场的下行分量和接收机波场的上行分量，源波场的下行分量是源波场的傅里叶变换的正分量与源波场的傅里叶变换的正分量的希尔伯特变换之和的傅里叶逆变换的实分量，接收机波场的上行分量是接收机波场的傅里叶变换的负分量与接收机波场的傅里叶变换的负分量的希尔伯特变换之和的傅里叶逆变换的实分量；通过将多个窗口的每个窗口中的源波场的下行分量相加来计算特定时间值处的源波场的下行分量；以及通过将多个窗口的每个窗口中的接收机波场的上行分量相加来计算特定时间值处的接收机波场的上行分量。
这个方面的其他实现方式包括对应的计算机系统、装置和记录在一个或多个计算机存储设备上的计算机程序，每一个均被配置为执行方法的动作。一个或多个计算机的系统可以被配置为通过在系统上安装的、在操作中使得系统执行动作的软件、固件、硬件或软件、固件或硬件的组合来执行特定的操作或动作。一个或多个计算机程序可以被配置为通过包括指令来执行特定操作或动作，所述指令在被数据处理装置执行时使得该装置执行动作。
前述和其他实现方式可以各自可选地以单独或组合的方式包括以下特征中的一个或多个：
与总的实现方式可组合的第一方面，包括：计算源波场的时间导数的傅里叶变换的正分量，以及计算接收机波场的时间导数的傅里叶变换的负分量。
与前述方面中的任何方面可组合的第二方面，其中，多个窗口是多个重叠的窗口。
与前述方面中的任何方面可组合的第三方面，其中，将每个窗口中的源波场的下行分量相加以及将每个窗口中的接收机波场的上行分量相加包括：以重叠区域消减的方式将每个窗口中的源波场的下行分量相加以及将每个窗口中的接收机波场的上行分量相加。
与前述方面中的任何方面可组合的第四方面，其中，傅里叶变换是1D空间傅里叶变换。
与前述方面中的任何方面可组合的第五方面，包括：基于特定时间值处的源波场的下行分量和接收机波场的上行分量来计算图像项。
本说明书的主题的一个或多个实现方式的细节在附图和以下描述中阐述。通过说明书、附图和权利要求书，所述主题的其他特征、方面和优点将变得显而易见。
附图说明
图1示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法。
图2示出了根据一些实现方式的用于基于源波场的下行分量和接收机波场的上行分量来计算图像的示例方法。
图3示出了根据一些实现方式的在去除逆时偏移(RTM)中的假象之前和之后的示例图像。
图4是示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法的流程图。
图5是示出了根据一些实现方式的用于提供与如本公开中描述的所述算法、方法、功能、处理、流程和过程相关联的计算功能的示例计算机系统的框图。
各附图中相似的附图标号和标记指示相似的元素。
具体实施方式
以下详细描述描述了去除由在具有错误传播方向组合的源波场和接收机波场上施加了逆时偏移(RTM)图像条件引起的假像，并提出以下详细描述以使本领域技术人员能够制造和使用一个或多个特定实现方式的上下文中所公开的主题。可以对所公开的实现方式进行各种修改、变更和置换，这对于本领域技术人员而言将是显而易见的，并且在不脱离本公开的范围的情况下，所定义的总的原理可以应用于其他实现方式和应用。因此，本公开不旨在限于所描述或图示的实现方式，而是应被赋予与所公开的原理和特征一致的最宽泛范围。
地球地下层的高分辨率深度图像对于成功进行油藏勘探、圈定和开发很重要。基于波动方程的地震深度偏移技术(例如RTM)是石油行业中(尤其是在诸如盐下勘探之类的复杂环境中)用于地震深度成像的合适的技术。然而，由于沿射线路径的非反射点处的源波场和接收机波场的不期望的互相关，RTM中存在成像噪声。
在RTM中，通过应用常规图像条件来获得图像
Figure BDA0002484231620000031
Figure BDA0002484231620000032
其中，
Figure BDA0002484231620000033
是源波场，
Figure BDA0002484231620000034
是接收机波场，
Figure BDA0002484231620000035
表示位置向量，t表示时间并且范围在0和Tmax之间(例如，最大记录时间)。通常，根据Liu et al.，2007，图像
Figure BDA0002484231620000041
可以分解为四个分量：
Figure BDA0002484231620000042
其中，Idu是下行源波场和上行接收机波场的互相关，Idd是下行源波场和下行接收机波场的互相关，Iuu是上行源波场和上行接收机波场的互相关，Iud是上行源波场和下行接收机波场的互相关。另外，四个分量表示为：
Figure BDA0002484231620000043
Figure BDA0002484231620000044
Figure BDA0002484231620000045
Figure BDA0002484231620000046
其中，
Figure BDA0002484231620000047
和
Figure BDA0002484231620000048
分别是源波场和接收机波场的下行分量，
Figure BDA0002484231620000049
和
Figure BDA00024842316200000410
分别是源波场和接收机波场的上行分量。
可以通过从图像
Figure BDA00024842316200000411
中消除
Figure BDA00024842316200000412
和
Figure BDA00024842316200000413
分量来抑制RTM中的低频成像噪声。为了计算剩余的图像
Figure BDA00024842316200000414
可以沿着深度坐标z对波场应用希尔伯特变换，之后可以计算图像融合-组合(image integration-combination)。
Figure BDA00024842316200000415
分量也可以形成假像，也需要将其消除。结果是，仅在图像
Figure BDA00024842316200000416
中保留了
Figure BDA00024842316200000417
分量，以便去除RTM中的假像。为了仅保留
Figure BDA00024842316200000418
分量，将源波场或接收机波场中的至少一个波场拆分为下行分量和上行分量。因此，先前描述的希尔伯特变换可以不产生图像
Figure BDA00024842316200000419
将波场拆分为下行分量和上行分量的一种方法是在频率-波数域内执行拆分。例如，可以(例如，经由2D傅里叶变换)将源波场
Figure BDA00024842316200000420
变换为
Figure BDA00024842316200000421
其中，ω表示频率，
Figure BDA00024842316200000422
表示波数向量。
Figure BDA00024842316200000423
或中的下行分量是ωkz＞0的部分，而
Figure BDA00024842316200000424
域中的上行分量是ωkz＜0的部分。然而，在RTM的实际实现方式中，通常通过有限差分方法来计算波场，并且由于数据量大，因此仅能存储波场的几个时间片。因此，不能执行从时间t到频率ω的波场的傅里叶变换、或用于挑选正/负ω部分的希尔伯特变换。
从高层面来看，所描述的方法提供了一种将波场拆分为下行分量和上行分量的近似方法。所描述的近似方法仅利用例如波场的当前时间片和当前时间导数片的知识将波场拆分为下行分量和上行分量。可以根据RTM中的波场的有限差分正演建模获取波场的当前时间片和当前时间导数片。基于该拆分，可以计算图像
Figure BDA0002484231620000051
并且可以从
Figure BDA0002484231620000052
中消除由错误分量组合(例如，
Figure BDA0002484231620000053
)产生的假像，以增强RTM的最终图像。结果是，RMT的最终图像可以用于例如更好地估计地球地下层的特性，以用于石油和天然气勘探。
图1示出了根据一些实现方式的用于计算源波场的下行分量和接收机波场的上行分量的示例方法100。例如，基于RTM正演建模，源波场
Figure BDA0002484231620000054
在t＝ti处的时间片为
Figure BDA0002484231620000055
接收机波场
Figure BDA0002484231620000056
在t＝ti处的时间片为
Figure BDA0002484231620000057
另外，可以获得t＝ti处的源波场
Figure BDA0002484231620000058
的时间导数片和接收机波场
Figure BDA0002484231620000059
的时间导数片，如下：
Figure BDA00024842316200000510
Figure BDA00024842316200000511
利用波场的时间片和时间导数片，可以通过图1中所示的工作流将波场拆分为下行分量和上行分量。
如图1所示，示例方法100开始于将波场(例如，源波场和/或接收机波场)的体积划分为多个重叠的窗口。例如，该窗口可以在3D空间中定义，并提供有限的3D子体积。在3D空间中的每个维度中，窗口与相邻的窗口重叠一半。通常，窗口大小可以是例如其维度大小的十分之一。在本公开中，假设窗口内介质的速度是近似均匀的，并且对于所有波数向量
Figure BDA00024842316200000512
波长的频散关系为：
Figure BDA00024842316200000513
其中，V是窗口内介质的平均速度，ω是圆频率。
对于多个重叠窗口中的每个窗口，对源波场的时间片(即
Figure BDA0002484231620000061
)、源波场的时间导数片(即
Figure BDA0002484231620000062
)、接收机波场的时间片(即
Figure BDA0002484231620000063
)和接收机波场的时间导数片(即
Figure BDA0002484231620000064
)应用局部空间傅里叶变换(例如，一维傅里叶变换)。所得到的傅里叶变换分别为
Figure BDA0002484231620000065
Figure BDA0002484231620000066
和
Figure BDA0002484231620000067
可以获得
Figure BDA0002484231620000068
和
Figure BDA0002484231620000069
的kz＞0的部分，以及
Figure BDA00024842316200000610
和
Figure BDA00024842316200000611
的kz＜0的部分，如下：
Figure BDA00024842316200000612
Figure BDA00024842316200000613
Figure BDA00024842316200000614
Figure BDA00024842316200000615
如果
Figure BDA00024842316200000616
仅对应于例如两个圆频率(如等式(9)中所示的一个正的和一个负的)，则等式(10)中的
Figure BDA00024842316200000617
可以分解为如下两项：
Figure BDA00024842316200000618
基于傅里叶变换的特性，可以获得下式：
Figure BDA00024842316200000619
其中，j是虚数单位。类似地，等式(12)中的
Figure BDA00024842316200000620
也可以被分解为如下两项：
Figure BDA0002484231620000071
基于傅里叶变换的特性，可以获得下式：
Figure BDA0002484231620000072
根据等式(15)和(17)，可以近似导出
Figure BDA0002484231620000073
和
Figure BDA0002484231620000074
的希尔伯特变换如下：
Figure BDA0002484231620000075
Figure BDA0002484231620000076
其中，
Figure BDA0002484231620000077
是波数向量
Figure BDA0002484231620000078
的范数，ε为用于避免等式中除以0而选择的用户指定的较小的数。例如，ε的值可以是0.0001或更小。
根据等式(14)和(18)，可以获得下式(即，具有kz＞0且ω＞0的分量)：
Figure BDA0002484231620000079
通过傅立叶变换的共轭特性，可以获得特定窗口的下行源波场：
Figure BDA00024842316200000710
其中，IFFT()是从
Figure BDA00024842316200000711
到
Figure BDA00024842316200000712
的空间傅里叶逆变换，Re()是返回实部的函数。
类似地，根据等式(16)和(19)，可以获得下式(即，具有kz＜0且ω＜0的分量)：
Figure BDA0002484231620000081
通过傅立叶变换的共轭特性，可以获得特定窗口的上行接收机波场：
Figure BDA0002484231620000082
为了获得整个下行源波场
Figure BDA0002484231620000083
将
Figure BDA0002484231620000084
的所有窗口结果(例如，针对多个重叠窗口中的每一个窗口经由等式(21)计算的结果)相加。类似地，为了获得整个上行接收机波场
Figure BDA0002484231620000085
将
Figure BDA0002484231620000086
的所有窗口结果(例如，针对多个重叠窗口中的每一个窗口经由等式(23)计算的结果)相加。如上所述，图1提供了计算源波场的下行分量和接收机波场的上行分量的近似方法。如等式(18)至(20)和(22)中所示，在一些实现方式中，仅利用波场的时间片和波场的时间导数片拆分波场的关键是选择ω＞0和/或ω＜0分量，而不在时域中使用傅里叶变换。
图2示出了根据一些实现方式的用于基于源波场的下行分量和接收机波场的上行分量来计算图像的示例方法200。例如，示例方法200开始于框205，其中，i和图像
Figure BDA0002484231620000087
均被初始化为0。在框210处，从源波场和接收机波场(即
Figure BDA0002484231620000088
和
Figure BDA0002484231620000089
)获得时间ti处的源波场和接收机波场的时间片。另外，还(例如，经由等式(7)和(8))获得时间ti处的源波场和接收机波场的时间导数片。在框215处，可以使用图1中所描述的近似方法来计算整个下行源波场(即
Figure BDA00024842316200000810
)和整个上行接收机波场(即
Figure BDA00024842316200000811
)。在框220处，将
Figure BDA00024842316200000814
添加到图像
Figure BDA00024842316200000813
在框225处，i增加1。在框230处，确定ti是否大于Tmax。Tmax是例如源波场和接收机波场的最大记录时间。如果ti不大于Tmax，则方法进行到框210。如果ti大于Tmax，则在框235处输出最终的图像
Figure BDA0002484231620000091
尽管本公开参考图2中的方法，但这是出于示例的目的，该文档的主题可以使用其他方法，基于使用图1中描述的近似方法计算的下行源波场和上行接收机波场来计算图像
Figure BDA0002484231620000092
图3示出了根据一些实现方式的在去除逆时偏移(RTM)中的假象之前和之后的示例图像300。在图3中，在顶部示出了使用Liu et.al.，2007提出的方法的RTM脉冲响应结果305。结果305示出了偏移噪声(例如，从
Figure BDA0002484231620000093
和
Figure BDA0002484231620000094
生成的)以及由三个箭头所指的假像(例如，从
Figure BDA0002484231620000095
生成的)。在图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处，计算特定窗口中的源波场的下行分量(例如，
Figure BDA0002484231620000101
)。在一些实现方式中，首先计算源波场的傅里叶变换的正分量和源波场的傅里叶变换的正分量的希尔伯特变换之和(例如，等式(20))。接下来，对该和应用傅里叶逆变换。然后，将特定窗口中的源波场的下行分量计算为和的傅立叶逆变换的实分量(例如，等式(21))。
在框425处，计算特定窗口中的接收机波场的上行分量(例如，
Figure BDA0002484231620000102
)。在一些实现方式中，首先计算接收机波场的傅里叶变换的负分量和接收机波场的傅里叶变换的负分量的希尔伯特变换之和(例如，等式(22))。接下来，对该和应用傅里叶逆变换。然后，将特定窗口中的接收机波场的上行分量计算为和的傅立叶逆变换的实分量(例如，等式(23))。
在框430处，计算特定时间值处的源波场的下行分量和接收机波场的上行分量(例如，
Figure BDA0002484231620000103
和
Figure BDA0002484231620000104
)。通过将多个窗口的每个窗口中的源波场的下行分量(例如，在等式(21)中计算的)相加来计算特定时间值处的源波场的下行分量。通过将多个窗口的每个窗口中的接收机波场的上行分量(例如，在等式(23)中计算的)相加来计算特定时间值处的接收机波场的上行分量。在一些实现方式中，以重叠区域消减(tapering)的方式将多个窗口的每个窗口中的源波场的下行分量相加。在一些实现方式中，以重叠区域消减的方式将多个窗口的每个窗口中的接收机波场的上行分量相加。
可以修改或重新配置图4所示的示例方法400，以包括附加的、更少的或不同的步骤(图4中未示出)，这些步骤可以按所示顺序或以不同顺序执行。例如，在框430之后，可以基于特定时间值处的源波场的下行分量和接收机波场的上行分量来计算图像(例如，
Figure BDA0002484231620000105
)。例如，在一些实现方式中，可以重复或迭代图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(或API 512和服务层513的组合)，通过系统总线503彼此交互或与接口504交互(或两者的组合)。API 512可以包括针对例程、数据结构和对象类的规范。API 512可以是独立于或依赖于计算机语言的，并且可以指代完整的接口、单个功能或甚至是一组API。服务层513向计算机502或可通信地耦接到计算机502的其他组件(无论是否被示出)提供软件服务。计算机502的功能可以是对于使用该服务层的所有服务消费者是可访问的。软件服务(例如由服务层513提供的软件服务)通过定义的接口提供可重用的、定义的业务功能。例如，接口可以是以JAVA、C++或以可扩展标记语言(XML)格式或其他合适格式提供数据的其他合适语言编写的软件。尽管被示为计算机502的集成组件，但是备选实现方式可以将API 512或服务层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、MAC OS、ANDROID、IOS或任意其他合适的常规操作系统)的数据处理装置的使用。
可以以任何形式的编程语言来写计算机程序(也可以称作或描述为程序、软件、软件应用、模块、软件模块、脚本或代码)，所述编程语言包括：编译或解译语言、或者声明或程序语言，并且可以以任何形式来部署计算机程序，包括部署为单独的程序或者部署为适合于用于计算环境的模块、组件、子例程、或者其他单元。计算机程序可以但不是必须与文件系统中的文件相对应。程序可以被存储在保存其他程序或数据(例如，存储在标记语言文档中的一个或多个脚本)的文件的一部分中、被存储在专用于所讨论的程序的单个文件中、或者被存储在多个协同文件(例如，存储一个或多个模块、子程序或代码部分的文件)中。计算机程序可以被部署为在一个计算机上或者在位于一个站点或分布在多个站点并且通过通信网络互联的多个计算机上执行。尽管各图中所出示的程序的部分被示出为通过各种对象、方法或其他过程实施各种特征和功能的各个模块，但是视情况而定程序可以替代地包括多个子模块、第三方服务、组件、库等。相反，各种组件的特征和功能可以视情况而定组合成单个组件。可以统计地、动态地或者统计地且动态地确定用于进行计算确定的阈值。
本说明书中描述的方法、处理或逻辑流可以由一个或多个可编程计算机来执行，所述一个或多个可编程计算机执行一个或多个计算机程序以通过操作输入数据并且生成输出来执行功能。方法、处理或逻辑流也可以由专用逻辑电路(例如CPU、FPGA或ASIC)来执行，并且装置也可以被实现为专用逻辑电路(例如CPU、FPGA或ASIC)。
适合于执行计算机程序的计算机可以基于通用或专用微处理器、这两者或任何其他类型的CPU。通常，CPU将从只读存储器(ROM)或随机存取存储器(RAM)或者这二者接收指令和数据。计算机的必不可少的元件是用于执行指令的CPU和用于存储指令和数据的一个或更多个存储器设备。通常，计算机还将包括用于存储数据的一个或多个大容量存储设备(例如，磁盘、磁光盘或光盘)，或可操作地耦接以便从所述一个或多个大容量存储设备接收或向其发送数据或两者。然而，计算机不必非要具有这些设备。此外，计算机可以嵌入在另一设备中，例如，移动电话、个人数字助理(PDA)、移动音频或视频播放器、游戏机、全球定位系统(GPS)接收机或者便携式存储设备(例如，通用串行总线(USB)闪存驱动器)，这仅是举几个例子。
适合于存储计算机程序指令和数据的计算机可读介质(暂时或非暂时的，视情况而定)包括所有形式的非易失性存储器、介质和存储器设备，包括例如半导体存储器设备，例如可擦除可编程只读存储器(EPROM)、电可擦除可编程只读存储器(EEPROM)和闪存设备；磁盘(例如内部硬盘或可移动盘)；磁光盘；以及CD ROM、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)小区、语音、视频、数据或其他合适信息(或通信类型的组合)。
计算系统可以包括客户端和服务器。客户端和服务器一般相互远离并且通常通过通信网络进行交互。客户端和服务器的关系通过在相应计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生。
尽管本说明书包含许多特定实现细节，然而这些细节不应被解释为对可以要求保护的范围或任何发明的范围的限制，而是作为可以专用于特定发明的特定实现方式的特征的描述。在单个实现方式中，还可以组合实现本说明书中在独立实现方式的上下文中描述的某些特征。反之，在单个实现方式的上下文中描述的各种特征也可以在多个实现方式中分开地或以任何合适的子组合来实现。此外，尽管可以将先前描述的特征描述为以某些组合起作用并且甚至最初要求如此保护，但是在一些情况下，来自所要求保护的组合的一个或多个特征可以从组合中删除，并且所要求保护的组合可以涉及子组合或子组合的变型。
已经描述了本主题的特定实现方式。对于本领域技术人员显而易见的是，所描述的实现方式的其他实现方式、改变和置换在下文的权利要求的范围内。尽管在附图或权利要求中以特定顺序描述了操作，但这不应被理解为：为了实现希望的结果，要求按所示出的特定顺序或按相继的顺序来执行这些操作，或者要求执行所有所示出的操作(一些操作可以被认为是可选的)。在一些情况下，多任务或并行处理(或者多任务和并行处理的组合)可以是有利的并且视情况来执行。
此外，在前述的实现中的各种系统模块和组件的分离或集成不应被理解为在所有实现方式中要求这样的分离或集成，并且应该理解的是，所描述的程序组件和系统一般可以一起集成在单个软件产品中或封装为多个软件产品。
因此，先前描述的示例实现方式不限定或限制本公开。在不脱离本公开的精神和范围的情况下，还可以有其他改变、替换和变化。
此外，任何要求保护的实现方式被认为适用于至少一种计算机实现的方法；存储用于执行计算机实现的方法的计算机可读指令的非暂时性计算机可读介质；以及计算机系统，该系统包括与硬件处理器可互操作地耦接的计算机存储器，所述硬件处理器被配置为执行计算机实现的方法或存储在非暂时性计算机可读介质上的指令。
参考文献
Faqi Liu，Guanquan Zhang，Scott A.Morton and Jacques P.Leveille，2011，Aneffective imaging condition for reverse-time migration using wavefielddecomposition：Geophysics，vol.76，no.1，S29-S39.
Faqi Liu，Guanquan Zhang，Scott A.Morton，Jacques P.Leveille，2007，Reverse-time migration using one-way wavefield imaging condition：77th AnnualInternational Meeting，SEG Expanded Abstracts.
Tong W.Fei，Yi Luo and Fuhao Qin，2014，An endemic problem in reverse-time migration：84th Annual International Meeting，SEG Expanded Abstracts。

Claims (16)

1.一种用于去除逆时偏移RTM中的假象的方法，包括：
由一个或多个硬件处理器，将源波场
Figure FDA0003337259920000011
和接收机波场
Figure FDA0003337259920000012
的体积划分为多个窗口，其中，
Figure FDA0003337259920000013
表示位置向量，并且t表示时间；
对于所述多个窗口中的每个窗口且在特定时间值t＝ti处：
由所述一个或多个硬件处理器，计算所述源波场的傅里叶变换的正分量
Figure FDA0003337259920000014
其中，
Figure FDA0003337259920000015
表示波数向量，其中，
Figure FDA0003337259920000016
并且其中
Figure FDA0003337259920000017
是所述源波场的所述傅里叶变换；
由所述一个或多个硬件处理器，计算所述源波场的时间导数的傅里叶变换的正分量
Figure FDA0003337259920000018
其中，
Figure FDA0003337259920000019
是所述源波场的时间导数的傅里叶变换；
由所述一个或多个硬件处理器，计算所述源波场的傅里叶变换的正分量的希尔伯特变换
Figure FDA00033372599200000110
由所述一个或多个硬件处理器，将特定窗口中所述源波场的下行分量
Figure FDA00033372599200000111
计算为所述
Figure FDA00033372599200000112
与所述
Figure FDA00033372599200000113
之和的傅里叶逆变换的实分量；
由所述一个或多个硬件处理器，计算所述接收机波场的傅里叶变换的负分量
Figure FDA00033372599200000114
其中
Figure FDA00033372599200000115
并且其中
Figure FDA00033372599200000116
是所述接收机波场的所述傅里叶变换；
由所述一个或多个硬件处理器，计算所述接收机波场的时间导数的傅里叶变换的负分量
Figure FDA00033372599200000117
其中，
Figure FDA00033372599200000118
是所述接收机波场的时间导数的傅里叶变换；
由所述一个或多个硬件处理器，计算所述接收机波场的傅里叶变换的负分量的希尔伯特变换
Figure FDA00033372599200000119
以及
由所述一个或多个硬件处理器，将所述特定窗口中所述接收机波场的上行分量
Figure FDA0003337259920000021
计算为所述
Figure FDA0003337259920000022
与所述
Figure FDA0003337259920000023
之和的傅里叶逆变换的实分量；
由所述一个或多个硬件处理器，通过将所述多个窗口的每个窗口中的所述源波场的下行分量相加来计算所述特定时间值处的所述源波场的下行分量
Figure FDA0003337259920000024
由所述一个或多个硬件处理器，通过将所述多个窗口的每个窗口中的所述接收机波场的上行分量相加来计算所述特定时间值处的所述接收机波场的上行分量
Figure FDA0003337259920000025
以及
由所述一个或多个硬件处理器，基于所计算的
Figure FDA0003337259920000026
和
Figure FDA0003337259920000027
生成图像，以去除逆时偏移RTM中的假象。
2.根据权利要求1所述的方法，其中，
Figure FDA0003337259920000028
并且
Figure FDA0003337259920000029
其中，V是所述特定窗口中的平均速度，
Figure FDA00033372599200000210
是所述波数向量的范数，ε是用户指定的数，并且j是虚数单位。
3.根据权利要求2所述的方法，其中，
Figure FDA00033372599200000211
并且
Figure FDA00033372599200000212
4.根据权利要求1所述的方法，其中，所述多个窗口是多个重叠的窗口。
5.根据权利要求4所述的方法，其中，将每个窗口中的所述源波场的下行分量相加以及将每个窗口中的所述接收机波场的上行分量相加包括：以重叠区域消减的方式将每个窗口w中的
Figure FDA00033372599200000213
相加以及将每个窗口w中的
Figure FDA00033372599200000214
相加。
6.根据权利要求1所述的方法，其中，所述傅里叶变换是1D空间傅里叶变换。
7.根据权利要求1所述的方法，还包括：
基于所述
Figure FDA00033372599200000215
和所述
Figure FDA00033372599200000216
来计算图像项
Figure FDA00033372599200000217
8.一种用于去除逆时偏移RTM中的假象的设备，包括：
存储器；以及
一个或多个硬件处理器，被布置为执行包括如下步骤的操作：
将源波场
Figure FDA0003337259920000031
和接收机波场
Figure FDA0003337259920000032
的体积划分为多个窗口，其中，
Figure FDA0003337259920000033
表示位置向量，并且t表示时间；
对于所述多个窗口中的每个窗口且在特定时间值t＝ti处：
计算所述源波场的傅里叶变换的正分量
Figure FDA0003337259920000034
其中，
Figure FDA0003337259920000035
表示波数向量，其中，
Figure FDA0003337259920000036
并且其中
Figure FDA0003337259920000037
是所述源波场的所述傅里叶变换；
计算所述源波场的时间导数的傅里叶变换的正分量
Figure FDA0003337259920000038
其中，
Figure FDA0003337259920000039
是所述源波场的时间导数的傅里叶变换；
计算所述源波场的傅里叶变换的正分量的希尔伯特变换
Figure FDA00033372599200000310
将特定窗口中所述源波场的下行分量
Figure FDA00033372599200000311
计算为所述
Figure FDA00033372599200000312
与所述
Figure FDA00033372599200000313
之和的傅里叶逆变换的实分量；
计算所述接收机波场的傅里叶变换的负分量
Figure FDA00033372599200000314
其中
Figure FDA00033372599200000315
并且其中
Figure FDA00033372599200000316
是所述接收机波场的所述傅里叶变换；
计算所述接收机波场的时间导数的傅里叶变换的负分量
Figure FDA00033372599200000317
其中，
Figure FDA00033372599200000318
是所述接收机波场的时间导数的傅里叶变换；
计算所述接收机波场的傅里叶变换的负分量的希尔伯特变换
Figure FDA00033372599200000319
以及
将所述特定窗口中所述接收机波场的上行分量
Figure FDA00033372599200000320
计算为所述
Figure FDA00033372599200000321
与所述
Figure FDA00033372599200000322
之和的傅里叶逆变换的实分量；
通过将所述多个窗口的每个窗口中的所述源波场的下行分量相加来计算所述特定时间值处的所述源波场的下行分量
Figure FDA00033372599200000323
通过将所述多个窗口的每个窗口中的所述接收机波场的上行分量相加来计算所述特定时间值处的所述接收机波场的上行分量
Figure FDA0003337259920000041
以及
基于所计算的
Figure FDA0003337259920000042
和
Figure FDA0003337259920000043
生成图像，以去除逆时偏移RTM中的假象。
9.根据权利要求8所述的设备，其中，
Figure FDA0003337259920000044
并且
Figure FDA0003337259920000045
其中，V是所述特定窗口中的平均速度，
Figure FDA0003337259920000046
是所述波数向量的范数，ε是用户指定的数，并且j是虚数单位。
10.根据权利要求9所述的设备，其中，
Figure FDA0003337259920000047
并且
Figure FDA0003337259920000048
11.根据权利要求8所述的设备，其中，所述多个窗口是多个重叠的窗口。
12.根据权利要求11所述的设备，其中，将每个窗口中的所述源波场的下行分量相加以及将每个窗口中的所述接收机波场的上行分量相加包括：以重叠区域消减的方式将每个窗口w中的
Figure FDA0003337259920000049
相加以及将每个窗口w中的
Figure FDA00033372599200000410
相加。
13.根据权利要求8所述的设备，其中，所述傅里叶变换是1D空间傅里叶变换。
14.根据权利要求8所述的设备，所述操作还包括：
基于所述
Figure FDA00033372599200000411
和所述
Figure FDA00033372599200000412
来计算图像项
Figure FDA00033372599200000413
15.一种存储指令的非暂时性计算机可读介质，所述指令能够由计算机系统执行以执行操作，所述操作包括：
将源波场
Figure FDA00033372599200000414
和接收机波场
Figure FDA00033372599200000415
的体积划分为多个窗口，其中，
Figure FDA00033372599200000416
表示位置向量，并且t表示时间；
对于所述多个窗口中的每个窗口且在特定时间值t＝ti处：
计算所述源波场的傅里叶变换的正分量
Figure FDA00033372599200000417
其中，
Figure FDA0003337259920000051
表示波数向量，其中，
Figure FDA0003337259920000052
并且其中
Figure FDA0003337259920000053
是所述源波场的所述傅里叶变换；
计算所述源波场的时间导数的傅里叶变换的正分量
Figure FDA0003337259920000054
其中，
Figure FDA0003337259920000055
是所述源波场的时间导数的傅里叶变换；
计算所述源波场的傅里叶变换的正分量的希尔伯特变换
Figure FDA0003337259920000056
将特定窗口中所述源波场的下行分量
Figure FDA0003337259920000057
计算为所述
Figure FDA0003337259920000058
与所述
Figure FDA0003337259920000059
之和的傅里叶逆变换的实分量；
计算所述接收机波场的傅里叶变换的负分量
Figure FDA00033372599200000510
其中
Figure FDA00033372599200000511
并且其中
Figure FDA00033372599200000512
是所述接收机波场的所述傅里叶变换；
计算所述接收机波场的时间导数的傅里叶变换的负分量
Figure FDA00033372599200000513
其中，
Figure FDA00033372599200000514
是所述接收机波场的时间导数的傅里叶变换；
计算所述接收机波场的傅里叶变换的负分量的希尔伯特变换
Figure FDA00033372599200000515
以及
将所述特定窗口中所述接收机波场的上行分量
Figure FDA00033372599200000516
计算为所述
Figure FDA00033372599200000517
与所述
Figure FDA00033372599200000518
之和的傅里叶逆变换的实分量；
通过将所述多个窗口的每个窗口中的所述源波场的下行分量相加来计算所述特定时间值处的所述源波场的下行分量
Figure FDA00033372599200000519
通过将所述多个窗口的每个窗口中的所述接收机波场的上行分量相加来计算所述特定时间值处的所述接收机波场的上行分量
Figure FDA00033372599200000520
以及
基于所计算的
Figure FDA00033372599200000521
和
Figure FDA00033372599200000522
生成图像，以去除逆时偏移RTM中的假象。
16.根据权利要求15所述的介质，其中，所述多个窗口是多个重叠的窗口。
CN201880072723.9A 2017-09-11 2018-09-07 逆时偏移中的假像去除 Active CN111971587B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US15/701,038 US10571586B2 (en) 2017-09-11 2017-09-11 False image removal in reverse time migration
US15/701,038 2017-09-11
PCT/US2018/049834 WO2019051151A1 (en) 2017-09-11 2018-09-07 ELIMINATION OF A MISSING IMAGE IN REPEATED MIGRATION

Publications (2)

Publication Number Publication Date
CN111971587A CN111971587A (zh) 2020-11-20
CN111971587B true CN111971587B (zh) 2022-04-01

Family

ID=63684581

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880072723.9A Active CN111971587B (zh) 2017-09-11 2018-09-07 逆时偏移中的假像去除

Country Status (7)

Country Link
US (1) US10571586B2 (zh)
EP (1) EP3682270A1 (zh)
JP (1) JP7254778B2 (zh)
CN (1) CN111971587B (zh)
CA (1) CA3075260A1 (zh)
SA (1) SA520411480B1 (zh)
WO (1) WO2019051151A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115453621B (zh) * 2022-09-14 2024-03-22 成都理工大学 基于一阶速度-应力方程的纵横波解耦分离假象去除方法
CN115755175B (zh) * 2022-11-18 2024-06-11 成都理工大学 基于Hilbert变换的复合波场弹性波逆时偏移成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103018777A (zh) * 2012-12-03 2013-04-03 中国石油天然气股份有限公司 一种能够消除逆时偏移噪声的成像方法及装置
CN103091710A (zh) * 2013-01-15 2013-05-08 中国石油天然气股份有限公司 一种逆时偏移成像方法及装置
CN105445786A (zh) * 2014-08-04 2016-03-30 中国石油化工股份有限公司 用于基于gpu获取叠前逆时偏移的方法及装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8116168B1 (en) 2008-06-18 2012-02-14 Halliburton Energy Services, Inc. Hybrid one-way and full-way wave equation migration
US20120051176A1 (en) 2010-08-31 2012-03-01 Chevron U.S.A. Inc. Reverse time migration back-scattering noise removal using decomposed wavefield directivity
MX346092B (es) 2010-11-12 2017-03-07 Chevron Usa Inc Sistema y método para investigar características subsuperficiales de una formación de roca.
US10371841B2 (en) * 2011-04-12 2019-08-06 Cgg Services Sas Device and method for calculating 3D reverse time migration in tilted orthorhombic media
ES2553741T3 (es) 2011-05-23 2015-12-11 Total Sa Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie
US9575194B2 (en) 2013-05-01 2017-02-21 Cgg Services Sas Method apparatus and system for migration noise attenuation and image enhancement
US9766357B2 (en) 2013-11-05 2017-09-19 Exxonmobil Upstream Research Company Seismic image dip decomposition estimation and recomposition
US10281604B2 (en) 2015-05-20 2019-05-07 Exxonmobil Upstream Research Company Method of removing noise in seismic reverse-time migration

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103018777A (zh) * 2012-12-03 2013-04-03 中国石油天然气股份有限公司 一种能够消除逆时偏移噪声的成像方法及装置
CN103091710A (zh) * 2013-01-15 2013-05-08 中国石油天然气股份有限公司 一种逆时偏移成像方法及装置
CN105445786A (zh) * 2014-08-04 2016-03-30 中国石油化工股份有限公司 用于基于gpu获取叠前逆时偏移的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An effective imaging condition for reverse-time migration using wavefield decomposition;FAQI LIU ET AL;《Geophysics》;20110101;第76卷(第1期);S29-S39 *
Angle gathers from reverse time migration using analytic wavefield propagation and decomposition in the time domain;JIANGTAO HU ET AL;《Geophysics》;20160101;第81卷(第1期);S1-S9 *

Also Published As

Publication number Publication date
CN111971587A (zh) 2020-11-20
JP7254778B2 (ja) 2023-04-10
US10571586B2 (en) 2020-02-25
JP2020533591A (ja) 2020-11-19
WO2019051151A1 (en) 2019-03-14
US20190079206A1 (en) 2019-03-14
SA520411480B1 (ar) 2022-09-15
EP3682270A1 (en) 2020-07-22
CA3075260A1 (en) 2019-03-14

Similar Documents

Publication Publication Date Title
CN110945386B (zh) 基于逆时偏移建模角度域共成像道集
CN110914718B (zh) 计算频域中地震速度反演的与振幅无关的梯度
US11275190B2 (en) Generating diffraction images based on wave equations
US10288755B2 (en) Seismic processing workflow for broadband single-sensor single-source land seismic data
US10775524B2 (en) Simultaneous wavefield reconstruction and receiver deghosting of seismic streamer data using an L1 inversion
CN111971587B (zh) 逆时偏移中的假像去除
US20210199827A1 (en) Full waveform inversion using time delayed seismic data
CN110612461B (zh) 降低表面散射噪声
CN110799856B (zh) 使用与地下方位角和反射角相关的全波形反演产生速度模型
JP7286667B2 (ja) 波動場分離への、直交化フィルタリングの適用
JP7358385B2 (ja) 地震画像の向上
US10310112B2 (en) Processing geophysical data using 3D norm-zero optimization for smoothing geophysical inversion data
US10996358B2 (en) Image-guided velocity interpolation using a mask cube
US12123994B2 (en) Enhancement of seismic data
US20220268957A1 (en) Enhancement of seismic data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant