专利名称：衰减采集观测系统在地震数据中导致的假频的系统和方法
技术领域：
本发明涉及衰减由检波器间隔导致的地震数据中的假频。
背景技术：
采集地震数据时，通常以“放炮”将地震能量发送到受关注的地质体中，然后在靠近震源排列的多个检波器记录地震能量。这多个检波器典型情况下被排列在规则的预定网格的网格点处。对于受关注地质体内的给定同相轴，存在着最高频率，在该最高频率以上由于检波器间隔而呈现假频。衰减这样的假频的常规技术在共炮点域中可能是不足够或不实用的，其中大炮检距数据往往作用在地震数据的叠加道中。
发明内容
本发明的一方面涉及地震数据去假频的方法，其中所述方法在计算机系统中实施，该计算机系统具有被配置为执行一个或多个计算机程序模块的一个或多个处理器，。在一个实施例中，所述方法包括(a)获得表示从位于单个震源位置的震源穿过受关注地质体到位于根据检波器间隔分开的检波器位置的检波器的地震能量传播的第一道地震数据，所述第一道地震数据进一步根据时间表示多个检波器波束波至；(b)确定所述第一道地震数据中表示的对应于至少一个检波器波束波至的中心点倾角，其中所述中心点倾角根据所述至少一个检波器波束波至的到达时间确定；以及(C)确定要应用到地震数据道的滤波器，以便从已经由检波器间隔导致了假频的所述第一道地震数据中去除地震能量的频率，所述滤波器根据所述对应于至少一个波束波至的所述中心点倾角和所述检波器间隔确定。本发明的另一方面涉及被配置为地震数据去假频的系统。在一个实施例中，所述系统具有被配置为执行计算机程序模块的一个或多个处理器。在一个实施例中，所述计算机程序模块包括数据模块、中心点倾角模块和滤波器确定模块。所述数据模块被配置为获得地震数据道，包括表示从位于单个震源位置的震源穿过受关注地质体到位于根据检波器间隔分开的检波器位置的检波器的地震能量传播的第一道地震数据，所述第一道地震数据进一步根据时间表示多个检波器波束波至。所述中心点倾角模块被配置为确定所述第一道地震数据中表示的所述对应于至少一个检波器波束波至的中心点倾角，所述中心点倾角模块被配置为根据到达时间确定中心点倾角。所述滤波器确定模块被配置为确定要应用到已经由所述检波器间隔形成假频的所述第一道地震数据的滤波器，以便去除已经由检波器间隔导致了假频的地震能量的频率，所述滤波器确定模块被配置为根据所述对应于至少一个波束波至的所述中心点倾角和所述检波器间隔来确定所述滤波器。当使用地震数据来模拟地震能量波束经由受关注地质体传播时，由本文介绍的系统和方法所确定的滤波器减少、衰减和/或去除了在地震采集期间由检波器位置之间间隔所导致的假频。此外，这还可以减少和/或去除根据地震能量的模拟波束所形成的在所述受关注地质体图像中出现的假象。参考附图(其全部形成本说明书的一部分，其中相同附图标记指明多幅图中对应的部分)考虑了以下说明和附带的权利要求书后，本发明的这些和其他目的、特点和特征，以及相关构造元素的操作方法和功能以及部件的组合和制造节约措施将变得更显然。不过，应当明确地理解，这些附图仅仅是为了展示和说明，而不试图作为限制本发明的定义。正如在说明书和权利要求书中的使用，单数形式的“某”和“所述”包括了复数对象，除非语境清楚地另外规定。
图I展示了根 据本发明的一个或多个实施例，被配置为处理表示地震能量经由受关注地质体传播的地震数据的系统；图2展示了根据本发明的一个或多个实施例，地震能量经由受关注地质体的传播；图3展示了根据本发明的一个或多个实施例，震源和检波器的观测系统；图4展示了根据本发明的一个或多个实施例，倾角时差椭圆；图5展示了根据本发明的一个或多个实施例，震源和检波器的观测系统；图6展示了根据本发明的一个或多个实施例，中心点倾角的确定；图7展示了根据本发明的一个或多个实施例，处理表示地震能量经由受关注地质体传播的地震数据的方法。
具体实施例方式本发明可以在系统和由计算机执行的计算机方法的一般语境中介绍和实施。这样的计算机可执行的指令可以包括若干程序、例程、对象、组件、数据结构和计算机软件技术，它们能够被用于执行具体的任务和处理抽象的数据类型。本发明的软件实施可以以不同的语言编码，用于各种各样计算平台和环境中的应用。应当认识到，本发明的范围和基本原理不限于任何具体的计算机软件技术。不仅如此，本领域的技术人员将认识到，使用硬件和软件配置的任何一种或其结合都可以实践本发明，包括但是不限于具有单个和/或多个处理器的计算机处理器系统的系统、手持设备、可编程消费者电子设备、小型计算机、大型计算机等。在分布式计算环境中也可以实践本发明，其中任务由若干服务器或其他处理设备执行，它们经由一个或多个数据通信网络链接。在分布式计算环境中，程序模块位于本机和远程计算机存储介质中都可以，包括内存存储设备。同样，用于计算机处理器的制成品比如⑶、预记录磁盘或其他等效设备，可以包括计算机程序存储介质和其上记录的程序装置，用于指导计算机处理器促进本发明的实施和实践。这样的设备和制成品也落入本发明的实质和范围之内。现在将参考附图介绍本发明的若干实施例。本发明能够以无数种方式实施，包括例如作为系统(包括计算机处理系统)、方法(包括计算机实施的方法)、装置、计算机可读介质、计算机程序产品、图形用户界面、网络入口或者在计算机可读存储器中有形地固定的数据结构。以下讨论了本发明的几个实施例。附图仅仅展示了本发明的典型实施例所以不应当认为对其范围和广度的限制。图I展示的系统10被配置为处理表示地震能量穿过受关注地质体传播的地震数据。地震能量穿过受关注地质体传播，从在或接近受关注地质体的一个或多个震源位置到在或接近受关注地质体的一个或多个检波器位置。在处理地震数据时，系统10将地震能量模拟为波束(如高斯波束)。系统10执行的处理包括确定地震数据所用的滤波器，它减少、衰减和/或去除由于检波器位置之间的间隔可能在地震数据中出现的假频。系统10可以对地震数据应用滤波器以衰减假频效应，它们在成像和其他后续处理期间可能导致假象和/或其他不准确。在一个实施例中，系统10包括电子存储器12、用户接口 14、一个或多个信息资源16、一个或多个处理器18以及/或者其他组件。
在一个实施例中，电子存储器12包括电子存储器介质，它以电子方式存储信息。电子存储器12的电子存储器介质可以包括与系统10整合地提供(即本质上非可拆卸)的系统存储器以及/或者经由例如某端口(如USB端口、火线端口等)或驱动器(如磁盘驱动器等)可拆卸地连接到系统10的可拆卸存储器之一或双方。电子存储器12可以包括一种或多种光学地可读的存储器介质(如光盘等)，磁性地可读的存储器介质(如磁带、磁性硬盘驱动器、软盘驱动器等)、基于电荷的存储器介质(如EEPR0M、RAM等)、固态存储器介质(如闪存驱动器等)以及/或者其他以电子方式可读的存储器介质。电子存储器12可以存储软件算法、由处理器18确定的信息、经由用户界面14收到的信息、从信息资源16收到的信息以及/或者使系统10能够正常工作的其他信息。电子存储器12可以是系统10内分开的组件，电子存储器12也可以与系统10的一个或多个组件(如处理器18)整合地提供。用户接口 14被配置为提供系统10与用户之间的接口，用户可以通过它对系统10提供信息和从系统10接收信息。这使得数据、结果和/或指令以及任何其他可通信项——统称为“信息”——能够在用户与系统10之间通信。正如本文使用，术语“用户”可以指单个个人，也可以指可能正在协同工作的一组个人。适合包括在用户接口 14中的接口设备的实例包括小键盘、按钮、开关、键盘、旋钮、控制杆、显示屏、触摸屏、扬声器、话筒、指示灯、听觉报警和打印机。在一个实施例中，用户接口 14实际上包括多个分开的接口。应当理解，其他通信技术，无论是硬连线的还是无线的，同样被本发明设想为用户接口。例如，本发明设想用户接口 14可以与由电子存储器12提供的可拆卸存储接口集成。在这个实例中，可以将信息从可拆卸存储器(如智能卡、闪存驱动器、可拆卸磁盘等)加载到系统10中，这使得用户能够定制系统10的实施。适合系统10作为用户接口 14使用的其他示例性输入设备和技术包括，但是不限于RS-232端口、RF链接、IR链接、调制解调器(电话、电缆或其他)。简而言之，本发明设想了作为用户接口 14与系统10交流信息的所有技术。信息资源16包括与受关注地质体以及/或者产生受关注地质体图像的过程有关的一个或多个信息源。作为非限制实例，信息源16之一可以包括在或接近受关注地质体处采集的地震数据、从其导出的信息以及/或者与该采集有关的信息。地震数据可以包括地震数据的各个道，或者在一个通道上记录的从震源穿过受关注地质体传播的地震能量的数据。从地震数据导出的信息可以包括例如速度模型、与用于模拟地震能量穿过受关注地质体传播的波束相关联的波束特性、与用于模拟地震能量穿过受关注地质体传播的波束相关联的格林函数以及/或者其他信息。与地震数据采集有关的信息可以包括，例如与地震能量源的位置和/或朝向有关的数据、一个或多个地震能量检波器的位置和/或朝向、能量由震源产生并被导入受关注地质体中的时间以及/或者其他信息。
处理器18被配置为提供系统10中的信息处理能力。因此，处理器18可以包括一个或多个数字处理器、模拟处理器、设计为处理信息的数字电路、设计为处理信息的模拟电路、状态机以及/或者以电子方式处理信息的其他机构。尽管图I所示的处理器18为单个实体，但是这仅仅为了展示目的。在某些实施中，处理器18可以包括多个处理单元。这些处理单元可以物理上位于同一设备或计算平台之内，或者处理器18可以表示协同运行的多台设备的处理功能。正如图I所示，处理器18可以被配置为执行一个或多个计算机程序模块。这一个或多个计算机程序模块可以包括一个或多个数据模块20、中心点倾角模块22、滤波器确定模块24、滤波器应用模块26以及/或者其他模块。处理器18可以被配置为由软件；硬件；固件；软件、硬件和/或固件的某种结合；以及/或者配置处理处理器18上功能的其他机构来执行模块20、22、24以及/或者26。应当认识到，尽管图I展示的模块20、22、24和26共同位于单个处理单元之内，但是在处理器18包括多个处理单元的若干实施中，一个或多个模块20、22、24和/或26可以 位于远离其他模块的位置。以下描述的对由不同模块20、22、24和/或26提供的功能的说明是为了展示性目的，而不意在限制，因为模块20、22、24和/或26的任何一个都可以提供比所描述的更多或更少的功能。例如，可以除去一个或多个模块20、22、24和/或26，并且其某些或全部功能可以由其他的模块20、22、24和/或26提供。作为另一个实例，处理器18可以被配置为执行一个或多个附加的模块，它们可以执行以下归结为模块20、22、24和/或26之一的某种或全部功能。数据模块20被配置为获得穿过受关注地质体的地震数据道。由数据模块20所得到给定的地震数据道根据时间表示来自地震能量源的单炮地震能量的地震能量的检测器波束的波至，以特定的检波器倾角在不同的到达时间到达检波器位置。给定道的检波器倾角是关于水平线的角度，包括在给定道中的所述检波器波束波至以该角度到达该检波器位置。作为非限制展示，图2描绘了检波器波束路径28，检波器波束沿该路径到达检波器30的检波器位置。检波器波束路径28具有的检波器倾角为Pd。沿着检波器波束路径28到达检波器30的检波器波束包括对应于来自震源32的地震能量的多个震源波束34 (图2显示为震源波束34a、震源波束34b、震源波束34c以及震源波束34d)的检波器波束。震源波束34从震源32穿过受关注地质体传播。地质体内的同相轴(如表面)使与震源波束34相关联的地震能量在某点移动到检波器波束路径28上，并到达检波器30，作为在检波器30和检波器倾角Pd所对应的道中表示的检波器波束波至。对应于震源波束34的每个检波器波束波至都具有其自己的到达时间。给定检波器波束波至的到达时间是与给定检波器波束波至相关联的地震能量穿过受关注地质体的路径以及与给定检波器波束波至相关联的地震能量穿过受关注地质体的速度的函数。在检波器30记录的地震数据道将根据时间表示在检波器30收到的地震能量。基于根据时间表示地震能量，能够区分由地震数据道所表示的检波器波束波至的到达时间。返回图1，由数据模块20获得的地震数据道包括叠加的数据。这意味着这些地震道的单道包括在受关注地质体处实际记录的地震数据道的聚集。为了降低噪声并改进整体数据质量而聚集(如相加、平均等)这些道。在一个实施例中，由数据模块20获得的地震数据道在共炮点域中叠加。这是指从同一震源位置产生一炮或多炮地震能量期间在多个检波器位置记录的地震数据道叠加形成的道。在一个实施例中，数据模块20被配置为获得先前已经产生的地震数据的叠加道(如通过叠加过程)。数据模块20可以从电子存储器12、从信息资源16之一以及/或者从地震数据的其他来源获得地震数据道。在一个实施例中，数据模块20被配置为通过处理地震数据获得地震数据的叠加道。例如，数据模块20可以叠加以及/或者以其他方式处理地震数据，以获得地震数据的适当叠加道。中心点倾角模块22被配置为确定在由数据模块20所获得的地震数据中表示的对应于检波器波束波至的中心点倾角。给定检波器波束波至的中心点倾角是此给定检波器波束波至的检波器倾角与对应于此给定检波器波束的地震能量的震源波束的震源倾角之和。
通过展示性说明，再次参考图2，震源波束34c被显示为具有震源倾角Ps。沿着检波器波束路径28到达检波器30的震源波束34c所对应的检波器波束波至具有根据以下关系(等式I)确定的中心点倾角Pm ( I) Pm=Pd+Ps返回图1，在一个实施例中，中心点倾角模块22被配置为通过射线追踪确定给定检波器波束波至的中心点倾角。在这个实施例中，基于给定检波器波束波至的到达时间和受关注地质体的模型，中心点倾角模块22被配置为经由射线追踪确定此给定对应于检波器波束波至的震源波束。一旦通过射线追踪已经确定了震源波束，就能够根据地震数据道的检波器倾角和已确定的震源波束的震源倾角，确定此给定检波器波束波至的中心点倾角。返回参考图2，从沿着检波器波束路径28到达检波器30的给定检波器波束波至的到达时间，能够确定震源波束34中对应的一个。根据受关注地质体的模型(如速度模型等)，进行射线追踪以识别会产生检波器波束在此给定检波器波束波至的到达时间到达检波器30的到达的震源波束。例如，震源波束34c可能被确定为对应此给定检波器波束波至。一旦由射线追踪已经确定了震源波束(如震源波束34c)，那么震源波束的震源倾角Ps就能够连同检波器倾角Pd使用，确定此给定检波器波束波至的中心点倾角(见以上等式(I))。对地震数据道内的许多检波器波束波至都能够重复这个过程。在一个实施例中，中心点倾角模块22不是使用射线追踪确定中心点倾角，而是以解析方式确定至少一个检波器波束波至的中心点倾角。在这个实施例中，中心点倾角模块22利用先前已知的受关注地质体处或附近的同相轴的时差速度，以及一个或多个解析射线追踪方程确定中心点倾角。例如，可以使用水底和/或受关注地质体处或附近的其他反射层的先前已知的时差速度。除了同相轴的时差速度，还可能知道该同相轴的到达时间。图3至图6展示了一组解析方程，用于根据受关注地质体处或附近的同相轴的已知参数(如时差速度、到达时间等)，确定中心点倾角。图3显示了震源32和检波器30的鸟瞰图。正如在图3中可见，震源32与检波器30之间的距离等于2c，其中c是震源32和检波器30的位置的偏移距。检波器30的位置被认为在坐标系的原点(0，0，0)。震源32的位置的坐标被表示为(_2hx，2hy，0)，其中hx和hy是X和y方向半炮检距的值。图4显示了震源32和检波器30坐落在表面36上，受关注地质体的侧面高程图。图4进一步包括同相轴(如水底等)的倾角时差椭圆38。对于给定的有限炮检距(如图4中的2c)、倾角时差椭圆38，垂直维度由零炮检距传播时间定义。零炮检距传播时间是从震源直达倾角时差椭圆38所对应的同相轴，然后返回到表面的到达时间。因此，零炮检距传播时间能够被表示为等式(2):
2氺^(2) tQ =-
V0
其中h表示零炮检距传播时间、Vtl表示时差速度而Ztl由图4中b定义。倾角时差椭圆38的水平轴由有限的炮检距(如2c)定义。使用时差速度能够将有限炮检距的传播时间联系到零炮检距的传播时间。因此，通过定义零炮检距时间(如通过已知同相轴的深度)和时差速度，对图4所示的有限炮检距定义了倾角时差椭圆38。图5显示了受关注地质体的鸟瞰图，包括在检波器30处从表面沿着检波器波束路径28 (向倾角时差椭圆)行进到受关注地质体中的向量DE。正如图5展示，向量DE包括沿着X轴的分量(DEX)、沿着y轴的分量(DEy)和沿着z轴的分量(DEZ)。从检波器30向倾角时差椭圆(图5中未显示)延伸的向量DE的三个分量不是独立的。这些参数之间的关系能够表示为等式(3)( 3 ) DEv * DE'. +DE,. *DE1 + DE, * DE7
' ' 3 y图6显示了受关注地质体的侧面高程图，向量DE从检波器30沿着检波器波束路径28向倾角时差椭圆38延伸。应当认识到，如果DEX、DEy和v已知，那么就定义了从检波器30到倾角时差椭圆的向量DE(如等式(3)能够解出DEzX由于检波器30和震源32的坐标已知，表面与检测器波束路径28之间的角度(如检测器波束路径28的检测器波束倾角)已知，并且Vtl能够用作V，所以就定义了向量DE。从已知的倾角时差椭圆38和DE所对应的检波器波束波至的已知(或假设)的参数，就能够通过解析几何分析确定在倾角时差椭圆38上DE所对应的检波器波束波至所接触的点(图6中标注为E)。从这个点，就能够确定从震源32行进到点E的震源波束对应的向量SE。表面与向量SE之间的角度将是向量SE所对应的震源波束的震源波束倾角。一旦向量SE所对应的震源波束的震源波束倾角已知，该震源波束倾角就能够连同检测器波束路径28的检波器倾角一起用于确定对应于SE的震源波束和对应于DE的检波器波束波至的中心点波束倾角(如根据等式(I)中表示的关系)。这个中心点倾角将是沿着检测器波束路径28到达检波器30的、由对应于倾角时差椭圆38的同相轴所产生的检波器波束波至的中心点倾角。返回图1，滤波器确定模块24被配置为确定要被应用到地震数据道的滤波器，以
便衰减在地震数据采集期间检波器位置之间的间隔所导致的假频。对于给定的地震数据
道，由滤波器确定模块24所确定的滤波器取决于此给定地震数据道中表示的一个或多个
检波器波束波至所对应的中心点倾角(如由中心点倾角模块22确定的倾角)、检波器位置之
间的间隔以及/或者其他参数。作为非限制实例，滤波器确定模块24可以根据以下关系(等
式(4))对给定地震数据道确定滤波器
,, k_7] ( 4 ) Ima, =
m
其中Fmax表示通过滤波器的最高频率，k表示常数，而d表示检波器位置的空间采样(或采集观测系统的相关参数)。在一个实施例中，从空间采样或者检波器间隔的先前记录的值自动地确定d。在一个实施例中，根据指定检波器间隔的用户选择的空间采样值确定cl。在中心点倾角模块22对单道内表示的多个检波器波束波至确定中心点(如通过射线追踪)的实施例中，在等式(4)中阐明的滤波器将对Fmax产生多个值。在这样的实施例中，滤波器确定模块24被配置为把滤波器确定为时变滤波器。例如，滤波器确定模块24可以为由中心点倾角模块22所确定的每个中心点确定F_。然后，滤波器确定模块24可以使用确定的中心点所对应的到达时间，作为到达时间的函数，即F_(t)来确定F_。滤波器应用模块26被配置为将滤波器确定模块24所确定的滤波器应用到适当的地震数据道。这引起从已经由检波器间隔形成假频的地震数据道中去除频率(如在Fmax以上的全部频率)。图7展示了衰减地震数据中假频的方法40。以下呈现的方法40的操作是示意性的。在某些实施例中，可以用一种或多种未描述的补充操作，以及/或者没有讨论的一种或多种操作来实行方法40。另外，图7中展示和以下描述的方法40的操作顺序不意味着是限制。在某些实施例中，方法40可以在一个或多个处理器件(如数字处理器、模拟处理器、设计为处理信息的数字电路、设计为处理信息的模拟电路、状态机以及/或者以电子方式处理信息的其他机构)中实施。一个或多个处理器件可以包括响应在电子存储介质上存储的电子指令执行方法40的某些或全部操作的一个或多个器件。一个或多个处理器件可以包括被配置为通过专门设计的硬件、固件和/或软件执行方法40的一种或多种操作的一个或多个器件。在操作42，获得第一道地震数据。这第一道地震数据表示从位于震源位置的震源穿过受关注地质体到达位于根据检波器间隔分开的检波器位置处的检波器的地震能量传播。第一道地震数据表示根据时间的、以第一检波器位置的第一检波器倾角，在不同到达时间到达第一检波器位置的地震能量的检波器波束波至。在一个实施例中，操作42由与(图I所示并且在以上描述的)数据模块20类似或相同的数据模块执行。在操作44，确定在第一道地震数据中表示的至少一个检波器波束波至所对应的中心点倾角。这可以包括对第一道地震数据中表示的多个检波器波束波至分别确定中心点倾角。操作44可以包括根据受关注地质体处或附近的同相轴的先前已知的时差速度，仅仅对第一道地震数据中单个检波器波束波至确定中心点倾角。在一个实施例中，操作44由与(图I所示并且在以上介绍的)中心点倾角模块22类似或与相同的中心点倾角模块执行。在操作46，确定用于第一道地震数据的滤波器。该滤波器可被应用到第一道地震数据以便从已经由检波器间隔形成假频的第一地震数据道中去除频率。确定该滤波器时根据在操作44确定的中心点倾角并根据检波器间隔。如果在操作44对多个检波器波束波至 已经确定了中心点倾角，操作46就可以包括确定时变滤波器。如果对单个检波器波束波至已经确定了中心点倾角，操作46就可以包括确定不时变的滤波器。在一个实施例中，操作46由与(图I所示并且在以上介绍的)滤波器确定模块24类似或相同的滤波器确定模块执行。
在操作48，在操作46确定的滤波器被应用到第一道地震数据。对第一道地震数据应用此滤波器引起从已经由检波器间隔形成假频的第一道地震数据中去除频率。在一个实施例中，操作48由与(图I所示并且在以上描述的)滤波器应用模块26类似或相同的滤波器应用模块执行。在循环50，在第一道的检波器位置处沿着所有可用的检波器波束倾角，对地震数据道执行操作42、44、46和48。在循环52，对所有可用的检波器波束位置都执行循环50。 尽管根据当前认为是最实用的和优选的实施例为了展示目的已经详细描述了本发明，但是应当理解，这样的细节仅仅是为了这种目的，并且本发明不限于所公开的实施例，而是相反，意在涵盖附带的权利要求书的实质和范围内的修改和等效布局。例如，应当理解，本发明设想在可能的范围，任何实施例的一种或多种特征都能够与任何其他实施例的一种或多种特征结合。
权利要求
1.一种地震数据去假频的方法，其中所述方法在计算机系统中实施，该计算机系统具有被配置为执行一个或多个计算机程序模块的一个或多个处理器，所述方法包括 (a)获得表示从位于单个震源位置的震源穿过受关注地质体到位于根据检波器间隔分开的检波器位置的检波器的地震能量传播的第一道地震数据，所述第一道地震数据进一步根据时间表示多个检波器波束波至，以及 (b)确定所述第一道地震数据中表示的对应于至少一个检波器波束波至的中心点倾角，其中所述中心点倾角根据所述至少一个检波器波束波至的到达时间确定；以及 (C)确定要应用到地震数据道的滤波器，以便从已经由检波器间隔导致了假频的所述第一道地震数据中去除地震能量的频率，所述滤波器根据所述对应于至少一个波束波至的所述中心点倾角和所述检波器间隔确定。
2.根据权利要求I的方法，其中所述第一道地震数据对应于第一检波器位置和第一 检波器倾角，以及其中所述方法进一步包括对根据时间表示从所述单个震源以第二检波器倾角到达所述第一检波器位置的地震能量的检波器波束波至的第二道地震数据重复操作(a)-(c)。
3.根据权利要求I的方法，进一步包括对根据时间表示从所述单个震源在第二检波器位置以第二检波器倾角到达所述第二检波器位置的地震能量的检波器波束波至的第二道地震数据重复操作(a)-(C)。
4.根据权利要求I的方法，进一步包括对所述第一道地震数据应用所述滤波器，以便从已经由检波器间隔导致了假频的所述第一道地震数据中去除地震能量的频率。
5.根据权利要求I的方法，其中所述中心点确定步骤包括确定以不同到达时间到达所述第一检波器位置并且在所述第一道地震数据中表示的多个检波器波束波至的每一个对应的所述中心点倾角，以及其中确定所述滤波器包括根据对应于所述多个检波器波束波至的所述中心点倾角确定所述滤波器。
6.根据权利要求5的方法，其中所述滤波器是时变滤波器。
7.根据权利要求I的方法，其中所述中心点确定步骤包括根据受关注地质体处或附近的同相轴的先前已知的时差速度，确定对应于单个检波器波束波至的所述中心点倾角。
8.根据权利要求7的方法，其中，受关注地质体处或附近的所述同相轴是水底。
9.一种被配置为地震数据去假频的系统，所述系统包括 被配置为执行计算机程序模块的一个或多个处理器，所述计算机程序模块包括 数据模块，被配置为获得地震数据道，包括表示从位于单个震源位置的震源穿过受关注地质体到位于根据检波器间隔分开的检波器位置的检波器的地震能量传播的第一道地震数据，所述第一道地震数据进一步根据时间表示多个检波器波束波至； 中心点倾角模块，被配置为确定所述第一道地震数据中表示的所述对应于至少一个检波器波束波至的中心点倾角，所述中心点倾角模块被配置为根据到达时间确定中心点倾角；以及 滤波器确定模块，被配置为确定要应用到所述第一道地震数据的滤波器，以便去除已经由检波器间隔导致了假频的地震能量的频率，所述滤波器确定模块被配置为根据所述对应于至少一个波束波至的所述中心点倾角和所述检波器间隔来确定所述滤波器。
10.根据权利要求9的系统，进一步包括滤波器应用模块，被配置为对对应的第一道地震数据应用所述滤波器，以便从已经由检波器间隔导致了假频的所述第一道地震数据中去除地震能量的频率。
11.根据权利要求9的系统，其中所述中心点倾角模块被配置为确定在所述第一道地震数据中表示的对应于多个检波器波束波至的中心点倾角。
12.根据权利要求11的系统，其中所述滤波器确定模块被配置为确定时变滤波器。
13.根据权利要求9的系统，其中所述中心点倾角模块被配置为根据受关注地质体处或附近的同相轴的先前已知的时差速度，确定在所述第一道地震数据中表示的对应于检波器波束波至的中心点倾角。
14.根据权利要求13的系统，其中，受关注地质体处或附近的所述同相轴是水底。
15.根据权利要求9的系统，其中，在所述地震数据中表示的地震能量的检波器波束已经由从所述震源发射的地震能量的震源波束的反射产生，使得给定检波器波束已经由在所述震源位置处具有某个震源倾角的对应震源波束产生，以及对应于给定检波器波束波至的所述中心点倾角是所述给定检波器波束波至的检波器倾角与对应于给定检波器波束波至的所述震源波束的震源倾角之和。
全文摘要
处理表示穿过受关注地质体的地震能量传播的地震数据。地震能量从受关注地质体处或附近的一个或多个震源位置穿过受关注地质体到达受关注地质体处或附近的一个或多个检波器位置。在处理地震数据时，地震能量被模拟为波束(如高斯波束)。所述处理包括确定地震数据的滤波器，该滤波器衰减由于检波器位置之间的间隔可能在地震数据中出现的假频。
文档编号G01V1/36GK102656481SQ201180004935
公开日2012年9月5日 申请日期2011年2月4日 优先权日2010年2月25日
发明者王躍 申请人:雪佛龙美国公司