CN102576088A - 集群地震能波至以增强地下成像的系统和方法 - Google Patents
集群地震能波至以增强地下成像的系统和方法 Download PDFInfo
- Publication number
- CN102576088A CN102576088A CN2010800474012A CN201080047401A CN102576088A CN 102576088 A CN102576088 A CN 102576088A CN 2010800474012 A CN2010800474012 A CN 2010800474012A CN 201080047401 A CN201080047401 A CN 201080047401A CN 102576088 A CN102576088 A CN 102576088A
- Authority
- CN
- China
- Prior art keywords
- ripple
- cluster
- seismic energy
- extremely
- point place
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000003384 imaging method Methods 0.000 title abstract description 20
- 210000002706 plastid Anatomy 0.000 claims description 98
- 238000006116 polymerization reaction Methods 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 10
- 230000002889 sympathetic effect Effects 0.000 claims description 7
- 230000000644 propagated effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 23
- 230000004931 aggregating effect Effects 0.000 abstract 1
- 230000002776 aggregation Effects 0.000 abstract 1
- 238000004220 aggregation Methods 0.000 abstract 1
- 238000007792 addition Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 230000003278 mimic effect Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 3
- 238000007616 round robin method Methods 0.000 description 3
- 230000010365 information processing Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000000379 polymerizing effect Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Navigation (AREA)
- Geophysics And Detection Of Objects (AREA)
- User Interface Of Digital Computer (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Image Generation (AREA)
Abstract
处理在感兴趣地质体处或附近采集的地震数据。这种处理可以包括从地震数据中形成感兴趣地质体的图像。地震数据可以通过聚合能量波至来处理，以减少确定感兴趣地质体的图像所必须执行的成像过程的数目。这种聚合可以基于在本文中被称为集群的能量波至的分组。
Description
技术领域
本发明涉及处理在感兴趣地质体处或附近采集的地震数据以形成感兴趣地质体的图像。
背景技术
根据在感兴趣地质体处或附近采集的地震数据来成像感兴趣地质体的技术是众所周知的。在一些传统技术中，将地震能波至模拟成波束，然后将各个波束用于将来自感兴趣地质体内的粗网格点的图像数据分配给周围的细网格点。在这些技术中，给定粗网格点处的每个模拟波至需要用于扩展来自给定粗网格点的图像数据的分立成像过程。这可能导致从信息处理的观点来看昂贵的成像技术。
在一些传统技术中，不是对给定粗网格点处的多个波至执行分立成像过程，而是可能选择和处理给定粗网格点处的单个波至。虽然这降低了与成像相关联的处理成本，但成像的准确性和/或精确性可能受损。
发明内容
本发明的一个方面涉及一种被配置成处理与感兴趣地质体相关联的地震数据的系统。在一个实施例中，该系统包含电子存储器和一个或多个处理器。该电子存储器被配置成存储代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息。该一个或多个处理器被配置成执行多个计算机程序模块。所述计算机程序模块包含波至模块、集群模块、聚合模块、以及图像模块。所述波至模块被配置成获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数。感兴趣地质体内的粗网格点包含第一网格点，以及所述波至模块被配置成获得第一网格点处的地震能波至的一个或多个参数。所述集群模块被配置成将粗网格点处的地震能波至分组成粗网格点的波至集群。所述波至集群包括第一波至集群和第二波至集群，第一波至集群包括第一网格点处的一个或多个地震能波至，第二波至集群包括第一网格点处的一个或多个地震能波至。所述聚合模块被配置成确定各个波至集群的聚合数据，使得第一波至集群的聚合数据反映包括在第一波至集群中的第一网格点处的每个地震能波至的参数，以及第二波至集群的聚合数据反映包括在第二波至集群中的第一网格点处的每个地震能波至的参数。所述图像模块被配置成实现波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体。所述图像模块被配置成实现第一波至集群的聚合数据和第二波至集群的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体。
本发明的另一个方面涉及一种处理与感兴趣地质体相关联的地震数据的计算机实现方法，其中，该方法是在计算机系统中实现的，所述计算机系统包含配置成执行一个或多个计算机程序模块的一个或多个处理器。在一个实施例中，该方法包含：将代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息存储到一个或多个处理器可访问的电子存储器中；在一个或多个处理器上获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数，其中，感兴趣地质体内的粗网格点包含第一网格点，以及其中，获得第一网格点处的地震能波至的一个或多个参数；在一个或多个处理器上将粗网格点处的地震能波至分组成粗网格点的波至集群，其中，所述波至集群包括第一波至集群和第二波至集群，所述第一波至集群包括第一网格点处的一个或多个地震能波至，所述第二波至集群包括第一网格点处的一个或多个地震能波至；在一个或多个处理器上确定各个波至集群的聚合数据，使得第一波至集群的聚合数据反映包括在第一波至集群中的第一网格点处的每个地震能波至的参数，以及第二波至集群的聚合数据反映包括在第二波至集群中的第一网格点处的每个地震能波至的参数；以及在一个或多个处理器上实现波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体，其中，第一波至集群的聚合数据和第二波至集群的聚合数据被实现，以在围绕第一网格点的细网格点处成像感兴趣地质体。
本发明的又一个方面涉及一种处理与感兴趣地质体相关联的地震数据的计算机实现方法，其中，该方法是在计算机系统中实现的，所述计算机系统包含被配置成执行一个或多个计算机程序模块的一个或多个处理器。在一个实施例中，所述方法包含：将代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息存储到一个或多个处理器可访问的电子存储器中；在一个或多个处理器上获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数，其中，感兴趣地质体内的粗网格点包含第一网格点和第二网格点，其中，获得第一网格点处的地震能波至的一个或多个参数，以及获得第二网格点处的地震能波至的一个或多个参数；在一个或多个处理器上确定粗网格点处的波至的聚合数据，使得第一网格点处的地震能波至的聚合数据反映第一网格点处的每个地震能波至的参数，以及第二网格点处的地震能波至的聚合数据反映第二网格点处的每个地震能波至的参数；以及在一个或多个处理器上实现波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体，其中，实现第一网格点处的地震能波至的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体，以及其中，实现第二网格点处的地震能波至的聚合数据，以便在围绕第二网格点的细网格点处成像感兴趣地质体。
通过参考相同标号表示各个图形中的相应部件、形成本说明书一部分的附图对如下描述和所附权利要求书加以考虑，可以更加清楚地了解本发明的这些和其它目的、特征以及特性，以及操作方法、结构的相关元件的功能、部件的组合和制造的经费。但是，显而易见，这些附图只是为了例示和描述的目的，并未打算作为限制本发明的定义。正如用在说明书和权利要求书中的那样，“一个”、“一种”、和“该”等的单数形式也包括复数指示物，除非上下文另有明确指出。
附图说明
图1例示了依照本发明的一个或多个实施例，被配置成处理地震数据的系统；
图2例示了依照本发明的一个或多个实施例，感兴趣地质体内的网格点处的地震能波至；
图3例示了依照本发明的一个或多个实施例，感兴趣地质体内的网格点处的地震能波至；
图4例示了依照本发明的一个或多个实施例的处理地震数据的方法；
图5例示了依照本发明的一个或多个实施例的处理地震数据的方法；以及
图6例示了依照本发明的一个或多个实施例的处理地震数据的方法。
具体实施方式
图1例示了被配置成处理在感兴趣地质体处或附近采集的地震数据的系统10。这种处理可包括从地震数据中形成感兴趣地质体的图像。系统10通过聚合能量波至来处理地震数据，以便减少确定感兴趣地质体的图像所必须执行的成像过程的数目。这种聚合可以基于本文被称为集群的能量波至分组。在一个实施例中，系统10包含电子存储器12、用户接口14、一个或多个信息资源16、一个或多个处理器18、和/或其它部件。
在一个实施例中，电子存储器12包含电子地存储信息的电子存储介质。电子存储器12的电子存储介质可以包括与系统10集成配备(即，基本上不可拆卸)的系统存储器和/或经由，例如，端口(例如，USB端口、火线端口等)或驱动器(例如，盘驱动器等)可拆卸地与系统10连接的可拆卸存储器之一或两者。电子存储器12可以包括光可读存储介质(例如，光盘等)、磁可读存储介质(例如，磁带、磁硬盘驱动器、软盘驱动器等)、基于电荷的存储介质(例如，EEPROM、RAM等)、固态存储介质(例如，闪速驱动器等)和/或其它电子可读存储介质的一种或多种。电子存储器12可以存储软件算法、处理器18所确定的信息、经由用户接口14接收的信息、从信息资源16接收的信息和/或使系统10能够适当起作用的其它信息。电子存储器12可以是系统10内的分立部件，或电子存储器12可以与系统10的一个或多个其它部件(例如，处理器18)集成配备。
用户接口14被配置成提供用户可以向系统10提供信息和从系统10接收信息的系统10与用户之间的接口。这使得可在用户与系统10之间传送被统称为“信息”的数据、结果和/或指令以及任何其它可传送项目。如本文所使用，术语“用户”可以指单个个体或可能协同工作的一群个体。适合包括在用户接口14中的接口设备的例子包括键板、按钮、开关、键盘、旋钮、手柄、显示屏、触摸屏、扬声器、麦克风、指示灯、声音报警器、和打印机。在一个实施例中，用户接口14实际上包括多个分立接口。
应该明白，本发明也可以将其它通信技术，无论是硬连线的还是无线的，设想成用户接口14。例如，本发明设想可以将用户接口14与电子存储器12所提供的可拆卸存储接口集成。在本例中，可以将信息从使用户能够定制系统10的实现的可拆卸存储器(例如，智能卡、闪速驱动器、移动盘等)装入系统10中。适合作为用户接口14与系统10一起使用的其它示范性输入设备和技术包括但不限于RS-232端口、RF链路、IR链路、调制解调器(电话、电缆等)。简而言之，本发明将与系统10传送信息的任何技术都设想成用户接口14。
信息资源16包括与感兴趣地质体和/或生成感兴趣地质体的图像的过程有关的一个或多个信息源。举非限制性的例子来说，信息资源16之一可以包括在感兴趣地质体处或附近采集的地震数据、从中导出的信息、和/或与采集有关的信息。地震数据可以包括地震数据的各条迹线、或在从震源传播通过感兴趣地质体的地震能的一个通道上记录的数据。从地震数据中导出的信息可以包括，例如，速度模型、与用于模拟地震能通过感兴趣地质体的传播的波束相关联的波束参数、与用于模拟地震能通过感兴趣地质体的传播的波束相关联的格林(Green)函数、和/或其它信息。与地震数据的采集有关的信息可以包括，例如，与地震能的震源的位置和/或取向、地震能的一个或多个检测器的位置和/或取向、震源产生能量并将能量引向感兴趣地质体的时间、和/或其它信息有关的数据。
处理器18被配置成在系统10中提供信息处理能力。这样，处理器18可以包括数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电子地处理信息的其它机构中的一个或多个。尽管处理器18在图1中被显示成单个实体，但这只是为了例示的目的。在一些实现中，处理器18可以包括多个处理单元。这些处理单元可以物理地处在同一个设备或计算平台内，或处理器18可以代表协同运行的多个设备的处理功能。
如图1所示，处理器18可被配置成执行一个或多个计算机程序模块。所述一个或多个计算机程序模块可以包括网格模块20、数据模块22、波至模块24、集群模块26、特征波至模块27、聚合模块28、图像模块30、和/或其它模块中的一个或多个。处理器18可被配置成通过软件；硬件；固件；软件、硬件和/或固件的某种组合；和/或用于配置处理器18上的处理能力的其它机构，来执行模块20，22，24，26，27，28和/或30。
应当懂得，尽管模块20，22，24，26，27，28和30在图1中被例示成共处在单个处理单元内，但在处理器18包括多个处理单元的实现中，模块20，22，24，26，27，28和/或30中的一个或多个可以处在远离其它模块的地方。对下面所述的不同模块20，22，24，26，27，28和/或30所提供的功能的描述是为了例示的目的，而不是为了限制的目的，因为模块20，22，24，26，27，28和/或30中的任何一个都可以提供比所述多或少的功能。例如，可以省略模块20，22，24，26，27，28和/或30中的一个或多个，并且可由模块20，22，24，26，27，28和/或30中的其它几个来提供其一些或所有功能。作为另一个例子，处理器18可被配置成执行可以执行在下面分配给模块20，22，24，26，27，28和/或30之一的一些或所有功能的一个或多个附加模块。
网格模块20被配置成获得通过感兴趣地质体的网格的多个网格点的位置。网格和/或网格点的位置可由网格模块20存储到电子存储器12中。网格点的位置可以通过坐标(例如，三维坐标)来指定。在一个实施例中，网格模块20被配置成从处理器18外部的源(例如，从信息资源16之一，经由用户接口14从用户那里等)获得带有网格点的位置的网格。
网格模块20所获得的网格包括粗网格点和细网格点。粗网格点比细网格点较不密集地分布在整个感兴趣地质体上。在一个实施例中，细网格点以规则间隔分布在粗网格点之间。
数据模块22被配置成获得地震数据和与之有关的信息。数据模块22从，例如，信息资源16之一、经由用户接口14从用户那里、和/或从其它源获得这样的数据和信息。地震数据是在感兴趣地质体处或附近采集的地震数据。在一个实施例中，获得的地震数据包括在感兴趣地质体处或附近采集地震数据期间记录的地震数据的各条迹线。地震数据的迹线可以是“原始的”，或者这些迹线可能在以前已被处理过。例如，以前可能对地震数据的迹线进行了加权(例如，高斯波束加权)和/或堆叠(例如，局部倾斜堆叠)。
波至模块24被配置成获得通过感兴趣地质体的网格中的粗网格点处的地震能波至。获得地震能波至包括获得描述在地震数据采集期间地震能主体通过感兴趣地质体中的网格点的传播的参数。对于各个粗网格点，波至模块24获得相应地震能波至组。例如，对于第一粗网格点，波至模块24获得描述第一粗网格点处的第一组地震能主体波至的参数。获得的参数分别描述第一粗网格点处的第一组波至中的每个波至。对于第二网格点，波至模块24获得分别描述第二粗网格点处的第二组地震能波至的参数。
在一个实施例中，波至模块24被配置成通过确定描述地震能主体到粗网格点的传播的参数来确定地震能波至。波至模块24可从数据模块22所获得的地震数据中确定所述参数。波至模块24可从描述通过感兴趣地质体的地震能主体的参数的函数中确定所述参数。例如，所述函数可以包括描述地震能主体通过感兴趣地质体的传播的格林函数。所述函数可由波至模块24确定，或可由波至模块24从外部源(例如，从信息资源16，从用户接口14等)获得。
在一个实施例中，波至模块24被配置成从存储或存取以前确定的描述地震能主体到粗网格点的传播的参数的外部源(例如，从信息资源16，从用户接口14等)，获得描述地震能主体的传播的参数。
在一个实施例中，将地震能主体模拟成诸如高斯波束的波束。在这个实施例中，由波至模块24获得的描述给定粗网格点处的给定波至的参数可以包括中心射线路径、行进时间(实的和/或虚的)、振幅、相位、围绕中心射线路径的波束形成和/或其它波束参数中的一个或多个。
集群模块26被配置成对于各个粗网格点，形成地震能波至的一个或多个集群的组，其中，地震能波至的一个集群是一群具有相似性质的地震能波至。例如，给定粗网格点处的地震能波至的一个集群包括具有指示包括在波至集群中的地震能主体具有相似传播历史(例如，中心射线路径)、运动和/或动态性质、和/或其它类似性质的参数的地震能波至。
举例来说，图2例示了生成传播通过感兴趣地质体34的地震能的震源32。由于该地震能传播通过感兴趣地质体34，图2示出了到达粗网格点42的三个能量主体36、38和40。能量主体36、38和40中的每一个具有相似的传播历史(例如，基本上直的中心射线路径)。这样，与集群模块26(如图1所示和如本文所示)相似或相同的集群模块可将能量主体36、38和40中的每一个的波至一起分组到一个集群中。
另一方面，图3例示了从震源32传播到粗网格点42，但具有几乎不同于能量主体36、38和40(如图2所示和如上所示)的传播历史的传播历史的地震能主体44。取代沿着相对直的中心射线路径传播，能量主体44的中心射线路径通过感兴趣地质体34内具有与感兴趣地质体34的其余部分不同的成分的区域46。取代直接通过区域46，能量主体44被区域46折射，并以有点迂回的路径行进到粗网格点42。
由于这条相对曲折路径，能量主体44将具有稍微不同于能量主体36、38和40(如图2所示和如上所示)的性质。这样，在一个实施例中，粗网格点42处的能量主体44的波至将被包括在粗网格点42处的与包括能量主体36、38和40(如图2所示和如上所示)的波至集群分开的波至集群中。粗网格点42处的包括区域46的波至的波至集群可以包括一个或多个其它波至(未示出)，或该波至集群可能只包括粗网格点42处的区域46的波至。
返回到图1，在一个实施例中，集群模块26根据行进时间和/或其空间导数的分析，使各个粗网格点处的波至形成各个粗网格点处的波至集群。通过将具有相似行进时间和/或行进时间的空间导数的波至分组在一起，将给定粗网格点处的地震能波至分组成给定粗网格点的波至集群。这无意成为限制性的，因为可以实现将地震能波至划分成集群的其它基础，而不偏离本公开的范围。
在一个实施例中，将地震能波至分组成集群的一个或多个方面可由用户配置(例如，经由用户接口14的输入)。举非限制性的例子来说，根据用户选择，集群模块26可以设置每个粗网格点的最大集群数量、每个粗网格点的最小集群数量、和/或每个网格点的绝对集群数量。举另一个例子来说，根据用户选择，集群模块26可以设置将波至分组成集群所要分析的地震能波至的参数。将地震能波至分组成集群的其它方面可由用户配置。
要懂得的是，尽管已经在在各个网格点处创建多个集群的上下文中描述了将粗网格点处的波至分组成集群，但这无意成为限制性的。在一个实施例中，集群模块26在每个粗网格点处创建单个波至集群(或不包括集群模块26，并且在随后处理中，简单地认为每个粗网格点处的所有波至被包括在单个组或集群中)。
特征波至模块27被配置成确定各个集群内的波至的特征地震能主体。在将地震能主体模拟成波束的实施例中，给定粗网格点处的给定波至集群的特征地震能主体是到达给定粗网格点的特征地震能波束。在一个实施例中，特征波束是通过在给定集群中对给定粗网格点处的所有波束波至的波束参数求平均来确定的。这种平均可被加权或不加权。作为加权平均的非限制性例子，可以使用各个波束波至的行进时间、振幅和/或传播路径来加权各个波束波至的参数以求平均。例如，在一个实施例中，特征波至模块27识别虚行进时间最小的波束波至。然后由特征波至模块27将用于确定波至集群的波束参数的加权平均的应用于给定波束波至的波束参数的权重确定成具有最小虚行进时间的波束波至与给定波束波至之间的传播角度差的余弦函数。
在一个实施例中，取代对给定集群内的波束波至求平均来确定特征波束，特征波至模块27简单地选择波束波至之一作为特征波束。举非限制性的例子来说，特征波至模块27可以选择具有最小虚行进时间的波束波至、具有最高振幅的波束波至、和/或满足一些其它准则的波束。
聚合模块28被配置成聚合集群内的波至。这种聚合导致各个集群的聚合数据。给定集群的聚合数据可被用在后续成像中，使得可以根据聚合数据对整个集群进行成像处理，而不是对给定集群中的每个波至分别进行成像处理。给定集群的聚合数据反映给定集群中的每个单独地震能波至，而不仅仅是从集群中选择单个波至。
在一个实施例中，为了聚合集群内的波至，聚合模块28通过时移、缩放和相加与集群中的波束波至相关联的地震数据迹线来确定集群的聚合数据。在这个实施例中，为了确定给定粗网格点处的给定波至集群的聚合数据，聚合模块28获得与已被分组成给定集群的给定粗网格点处的地震能波至相关联的地震数据迹线。聚合模块28可以从数据模块22获得这些地震数据迹线。
为了时移与获得的地震能波至相关联的地震数据迹线，聚合模块28比较给定网格点处的各个地震能波至的行进时间与给定集群的特征波束的行进时间。具体地，将与包括在给定集群中的第一波至相关联的地震数据迹线时移第一时移量，以及将与包括在给定集群中的第二波至相关联的地震数据迹线时移第二时移量。第一时移量是根据给定粗网格点处的第一波至的行进时间与给定粗网格点处的特征波束波至的行进时间之间的时间差来确定的。在一个实施例中，第一时移量是给定粗网格点处的第一波至的行进时间与特征波束波至的行进时间之间的时间差。第二时移量是根据给定粗网格点处的第二波至的行进时间与给定粗网格点处的特征波束波至的行进时间之间的时间差来确定的。
为了缩放与获得的地震能波至相关联的地震数据迹线，聚合模块28使用地震能波至的振幅和/或虚行进时间。虚行进时间被用在频域中的指数函数中以获得缩放值。然后，这些缩放值和振幅被用于乘以由粗网格间隔、波至行进时间和波至行进时间的空间导数定义的有限时间窗内的相关数据迹线。
一旦时移和/或缩放了与给定粗网格点处的给定集群的波至相关联的地震数据迹线，聚合模块28相加被时移和/或缩放的迹线，使得给定集群与代表被时移、缩放和/或相加的地震数据迹线的地震数据的单个通道相关联。然后在图像处理(例如，如下面针对图像模块30所述)中实现这个与给定波至集群相关联的代表被时移、缩放和/或相加的地震数据迹线的地震数据的单个通道。
在一个实施例中，聚合模块28并不通过相加各条地震数据迹线来生成集群的聚合数据。在这个实施例中，聚合模块28确定与各个集群相对应的聚合子波，然后可将聚合子波实现成由子波描述的波场。
举非限制性的例子来说，在将地震能模拟成地震能波束的实施例中，给定集群内的给定粗网格点处的每个波束波至具有由波至模块24确定的波束参数。这些波束参数可包括行进时间、振幅和传播方向。对于带限震源子波(与采集时地震能震源的实际位置和/或取向相对应)和给定粗网格点，可将给定集群中的每个波束波至表达成以波束行进时间为中心的短子波。短子波具有其自身的中心到达时间、振幅、相位、和/或其它参数。然后，将表达给定集群中的波束波至的短子波聚合成给定粗网格点处的集群子波。子波的聚合可包括相加各个子波。所述相加可被加权或不被加权。作为一个非限制性例子，可以作为每个单独波至与特征波至之间的传播角度差的余弦函数来计算权重。将这个权重值用于通过相乘来缩放子波。然后，可以在后续图像处理(例如，如下面相对于图像模块30所述的)中实现所述子波。
图像模块30被配置成实现图像模块30为粗网格点处的波至集群生成的聚合数据以成像感兴趣地质体。成像感兴趣地质体包括扩展粗网格点处的聚合数据以成像感兴趣地质体内的细网格点。
如上面所讨论的，图像模块30可以通过时移、缩放和/或相加与集群中的地震能波至相关联的地震数据迹线来生成集群的聚合数据。这为各个集群生成了代表与相应各个集群相关联的被时移、缩放和/或相加的地震数据迹线的地震数据的相应单个通道。在一个实施例中，图像模块30通过应用与给定粗网格点处的给定波至集群相对应的被时移、缩放和/或相加的地震数据以在围绕给定粗网格点的细网格点上成像，来实现这种聚合地震数据。为了应用被时移、缩放和/或相加的地震数据以在细网格点上成像，图像模块30使用特征波至模块27为给定粗网格点处的给定集群确定的特征波束波至的波束性质(和/或其空间和/或时间导数)。应当懂得，这使得地震数据的所有相加迹线能够在单个成像过程中被应用于围绕给定粗网格点的成像过程，而不是每条迹线的各个成像过程。
如上面所讨论的，聚合模块28可以通过确定与粗网格点处的各个集群相对应的子波来生成粗网格点处的波至集群的聚合数据。在一个实施例中，图像模块30实现与粗网格点处的各个波至集群相对应的子波以成像感兴趣地质体。为了达到这个目的，图像模块30堆叠与各个集群相对应的子波，以便导出粗网格点处的子波，并且从从数据模块22获得的地震数据迹线和波场中，形成可被扩展到围绕粗网格点的细网格点的图像迹线。
例如，在给定粗网格点处，图像模块30将堆叠的集群子波与地震数据迹线互相关，以导出给定粗网格点处的图像迹线。然后，图像模块30获得通过给定粗网格点的地震数据迹线，并实现互相关的集群子波，以从地震数据迹线中导出图像迹线。可对通过给定粗网格点的多条获得的地震数据迹线重复这个过程。
为了成像围绕给定粗网格点的细网格点，扩展通过给定粗网格点的图像迹线。可以通过特征波至模块27为给定粗网格点确定的特征波至的空间导数行进时间，将图像迹线扩展到细网格点。
图4例示了依照本发明的一个或多个实施例，处理地震数据以便获得感兴趣地质体的图像的方法48。下面展示的方法48的操作旨在例示。在一些实施例中，方法48可以通过未描述的一个或多个附加操作，和/或没有所讨论的一个或多个操作来完成。另外，在图4中例示和下面描述方法48的操作的次序无意成为限制性的。
在一些实施例中，方法48可以在一个或多个处理设备(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电子地处理信息的其它机构)中实现。一个或多个处理设备可以包括响应于电子地存储在电子存储介质上的指令来执行方法48的一些或所有操作的一个或多个设备。一个或多个处理设备可以包括通过专门被设计成执行方法48的一个或多个操作的硬件、固件和/或软件而配置的一个或多个设备。
在操作50中，获得通过感兴趣地质体的网格。获得通过感兴趣地质体的网格包括获得感兴趣地质体内的多个粗网格点和多个细网格点的位置。在一个实施例中，操作50由与网格模块20(如图1所示和如上所述)相同或相似的网格模块来执行。
在操作52中，获得给定粗网格点处的地震能波至。获得给定粗网格点处的地震能波至包括获得描述给定网格点处的波至的性质的一个或多个参数。例如，可将地震能波至模拟成诸如高斯波束的波束，以及这些参数可包括波束的波束参数。在一个实施例中，操作52由与波至模块24(如图1所示和如上所述)相同或相似的波至模块来执行。
在操作54中，将给定粗网格点处的地震能波至分组成一个或多个集群。根据在操作52中获得的参数的相似性和/或获得的参数所描述的性质，将地震能波至分组成一个或多个集群。将地震能波至分组成一个或多个集群可部分基于一个或多个用户输入(例如，最大集群、最小集群、限定的集群数目等)。在一个实施例中，操作54由与集群模块26(如图1所示和如上所述)相同或相似的集群模块来执行。
在操作56中，为给定粗网格点处的给定波至集群确定特征波至。特征波至是根据为给定波至集群内的地震能波至获得的参数来确定的。举非限制性的例子来说，可对给定波至集群内的地震能波至的参数求平均来确定特征波至的参数。这种平均可被加权或不被加权。作为另一个非限制性例子，可以选择波至集群内的波至之一作为特征波至。在一个实施例中，操作56由与特征波至模块27(如图1所示和如上所述)相同或相似的特征波至模块来执行。
在操作58中，确定给定粗网格点处的给定波至集群的聚合数据。给定集群的聚合数据使得能够进行计及给定波至集群内的所有地震能波至的给定粗网格点处的统一成像处理。换句话说，根据给定集群的聚合数据，进行不计及给定波至集群内的各个波至的图像处理。取而代之，在后续图像处理中，使用聚合数据来代替给定波至集群内的各个波至。在一个实施例中，操作58由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作60中，在给定网格点处以及在围绕给定网格点的细网格点处成像感兴趣地质体。这种成像不是通过将与给定粗网格点处的各个波至相关联的图像信息扩展到细网格点的分立成像过程来进行的。取而代之，操作60所进行的成像充分利用为给定粗网格点处的波至集群获得的聚合数据，以减少在细网格点处导出图像信息所需的处理量。
在一个实施例中，方法48对于操作54所创建的给定网格点处的所有波至集群在操作56、58和60上循环。在一个实施例中，一旦完成了对于给定网格点处的所有波至集群的操作56、58和60，方法48对于在操作50中获得的每个粗网格点在操作52，54，56，58和60上循环。结果是形成感兴趣地质体的图像。
图5例示了依照本发明的一个或多个实施例，处理地震数据以便获得感兴趣地质体的图像的方法62。下面展示的方法62的操作旨在例示。在一些实施例中，方法62可以通过未描述的一个或多个附加操作和/或没有所讨论的一个或多个操作来完成。另外，在图5中例示和下面描述方法62的操作的次序无意成为限制性的。
在一些实施例中，方法62可以在一个或多个处理设备(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电子地处理信息的其它机构)中实现。一个或多个处理设备可以包括响应于电子地存储在电子存储介质上的指令来执行方法62的一些或所有操作的一个或多个设备。一个或多个处理设备可以包括通过专门被设计成执行方法62的一个或多个操作的硬件、固件和/或软件而配置的一个或多个设备。
在一个实施例中，方法62在与方法48(如图4所示和如上所述)相同或相似的整体方法内实现。具体地说，方法62可被实现成图4中的方法48内的操作54和56。
具体参照图5，在操作64中，对于给定粗网格点以及对于给定粗网格点处的给定波至集群，获得地震数据迹线。该地震数据迹线通过给定粗网格点。在一个实施例中，操作64由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作66中，时移在操作66中获得的地震数据迹线。地震数据迹线是根据给定波至集群中的波至的行进时间与给定波至集群的特征波至的行进时间之间的时间差(例如，如上面关于显示在图4中的方法48的操作42所述的被先前确定的)而被时移的。在一个实施例中，操作66由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作68中，缩放地震数据迹线。地震数据迹线可根据给定网格点处的给定波至集群当中的波至的振幅而被缩放。这种缩放可能是绝对的，或基于给定波至集群中的波至和/或给定波至集群的特征波至的相对比较。在一个实施例中，操作68由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作70中，相加被时移和/或被缩放的地震数据迹线(例如，与以前处理的地震数据迹线)。这种相加可被加权或不被加权。在一个实施例中，操作70由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
然后，方法62在操作64、66、68和70上循环。一旦方法62对所有适当地震数据迹线都完成了操作64、66、68和70上的循环，则操作70中对地震数据迹线的相加导致计及给定粗网格点处的给定波至集群的所有波至的地震数据的单个通道。这种地震数据的单个通道是给定网格点处的给定集群的聚合数据。
在操作72中，将给定粗网格点处的给定波至集群的聚合数据扩展到围绕给定粗网格点的细网格点以成像细网格点。为了扩展给定粗网格点处的给定波至集群的聚合数据，使用给定波至集群的特征波至、特征波至的波束参数、和/或特征波至的行进时间的空间导数。
然后，方法62在给定粗网格点处的所有集群上循环。结果是给定粗网格点和周围的细网格点的图像。作为较大的整体方法(例如，图4的方法48)的一部分，可以对于每个给定粗网格点再次循环方法62，以便成像整个感兴趣地质体内的细网格点。
图6例示了依照本发明的一个或多个实施例，处理地震数据以便获得感兴趣地质体的图像的方法74。下面展示的方法74的操作旨在例示。在一些实施例中，方法74可以通过未描述的一个或多个附加操作和/或没有所讨论的一个或多个操作来完成。另外，在图6中例示和下面描述方法74的操作的次序无意成为限制性的。
在一些实施例中，方法74可以在一个或多个处理设备(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电处理信息的其它机构)中实现。一个或多个处理设备可以包括响应于电子地存储在电子存储介质上的指令来执行方法62的一些或所有操作的一个或多个设备。一个或多个处理设备可以包括通过专门被设计成执行方法62的一个或多个操作的硬件、固件和/或软件而配置的一个或多个设备。
在一个实施例中，方法74在与方法48(如图4所示和如上所述)相同或相似的整体方法内实现。具体地说，方法74可以实现成图4中的方法48内的操作54和54。
具体参照图6，在操作76中，获得感兴趣地质体中的给定粗网格点处的地震能波至的集群。获得波至集群可以包括获得波至集群中的各个波至的波束参数、和/或先前为集群波至确定的特征波至。在一个实施例中，操作76由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作78中，对于在操作76中获得的波至集群内的给定波至，确定中心在给定粗网格点的子波。所述子波是从给定波至的参数中确定的震源子波。在一个实施例中，操作78由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
在操作80中，将为给定粗网格点处的给定波至确定的子波与从在操作76中获得的波至集群中的波至中确定的其它子波堆叠。在一个实施例中，操作80由与聚合模块28(如图1所示和如上所述)相同或相似的聚合模块来执行。
方法74对于在操作76中获得的给定粗网格点处的波至集群中的所有波至，在操作78和80上循环。结果是对于给定粗网格点处的波至集群，给定粗网格点处的堆叠子波。这个堆叠子波构成给定粗网格点处的波至集群的聚合数据。也就是说，后续图像处理实现堆叠子波，而不用分别回头参照波至集群中的地震能波至。
在操作82中，获得给定粗网格点的地震数据迹线。该地震数据迹线可以是地震数据的“原始”迹线，或已经在之前被处理过(例如，被倾斜堆叠和/或被波束加权)。在一个实施例中，操作82由与数据模块22(如图1所示和如上所述)相同或相似的数据模块来执行。
在操作84中，实现在操作80中确定的堆叠子波和在操作82中获得的地震数据迹线，以确定给定粗网格点的图像迹线。例如，可以使堆叠子波与地震数据迹线互相关以确定图像迹线。在一个实施例中，操作82由与图像模块30(如图1所示和如上所述)相同或相似的图像模块来执行。
在一个实施例中，取代上面在操作80中堆叠子波，可将每个子波与一条相关数据迹线互相关以形成一条图像迹线。对于每个子波，相关数据迹线可以是同一个或可以不同。对于关于操作78所述的若干子波，可以创建若干图像迹线。然后，可将这些图像迹线堆叠在一起，以形成可在后续步骤中使用的单条图像迹线。
在操作86中，将在操作84中确定的图像迹线扩展到围绕给定粗网格点的细网格点。为了将图像迹线扩展到细网格点，实现波至集群的特征波至的参数和/或特征波至的行进时间的空间导数。在一个实施例中，操作86由与图像模块30(如图1所示和如上所述)相同或相似的图像模块来执行。
一旦将图像迹线扩展到围绕给定粗网格点的细网格点，方法74对于可用于给定粗网格点的所有地震数据迹线在操作82、84和86上循环。这导致确定被扩展以成像围绕给定粗网格点的细网格点的多条图像迹线(可获得多个地震数据迹线)。
在为所有可用地震数据迹线循环了操作82、84和86之后，为存在用于给定粗网格点的任何其它波至集群循环操作76、78、80、82、84和86。在一个实施例中，在为给定粗网格点处的波至集群循环了操作76、78、80、82、84和86之后，可以对于感兴趣地质体内的多个粗网格点再次循环方法74。多个粗网格点的循环方法可以是整体方法(例如，如图4所示和如上所述的方法48)的一部分。
尽管为了例示的目的，根据当前认为最实用的优选实施例对本发明作了详细描述，但应当明白，这样的细节仅仅为了那个目的，本发明不局限于所公开的实施例，而是相反，旨在涵盖在所附权利要求书的精神和范围之内的所有修改和等效安排。例如，要明白的是，本发明设想，可以尽可能地将任何实施例的一个或多个特征与任何其它实施例的一个或多个特征组合。
Claims (15)
1.一种被配置成处理与感兴趣地质体相关联的地震数据的系统，该系统包含：
电子存储器，被配置成存储代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息；以及
一个或多个处理器，被配置成执行多个计算机程序模块，该计算机程序模块包含：
波至模块，被配置成获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数，
其中，感兴趣地质体内的粗网格点包含第一网格点，以及所述波至模块被配置成获得第一网格点处的地震能波至的一个或多个参数；
集群模块，被配置成将粗网格点处的地震能波至分组成粗网格点的波至集群，其中，所述波至集群包括第一波至集群和第二波至集群，所述第一波至集群包括第一网格点处的一个或多个地震能波至，所述第二波至集群包括第一网格点处的一个或多个地震能波至；
聚合模块，被配置成确定各个波至集群的聚合数据，使得第一波至集群的聚合数据反映包括在第一波至集群中的第一网格点处的每个地震能波至的参数，以及第二波至集群的聚合数据反映包括在第二波至集群中的第一网格点处的每个地震能波至的参数；以及
图像模块，被配置成实现波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体，
其中，所述图像模块被配置成实现第一波至集群的聚合数据和第二波至集群的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体。
2.如权利要求1所述的系统，其中，所述地震能波至被模拟成波束，以及其中，所述波至模块被配置成获得感兴趣地质体内的粗网格点处的各个地震能波至的波束参数。
3.如权利要求1所述的系统，进一步包含特征波至模块，被配置成确定与各个地震能波至集群相对应的特征波至，其中，确定特征波至包含确定特征波至的参数。
4.如权利要求3所述的系统，其中，所述聚合模块被配置成：
实现所述特征波至模块为第一地震能波至集群确定的特征波至的参数，以确定第一地震能波至集群的聚合数据；以及
实现所述特征波至模块为第二地震能波至集群确定的特征波至的参数，以确定第二地震能波至集群的聚合数据。
5.如权利要求3所述的系统，其中，所述图像模块被配置成：
实现(i)所述特征波至模块为第一地震能波至集群确定的特征波至的参数，和(ii)第一地震能波至集群的聚合数据，以便在围绕第一网格点的细网格点处成像感兴趣地质体；以及
实现(i)所述特征波至模块为第二地震能波至集群确定的特征波至的参数，和(ii)第二地震能波至集群的聚合数据，以便在围绕第一网格点的细网格点处成像感兴趣地质体。
6.一种处理与感兴趣地质体相关联的地震数据的计算机实现方法，其中，该方法是在计算机系统中实现的，该计算机系统包含被配置成执行一个或多个计算机程序模块的一个或多个处理器，该方法包含：
将代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息存储到一个或多个处理器可访问的电子存储器中；
在一个或多个处理器上获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数，
其中，感兴趣地质体内的粗网格点包含第一网格点，以及其中，获得第一网格点处的地震能波至的一个或多个参数；
在一个或多个处理器上将粗网格点处的地震能波至分组成粗网格点的波至集群，其中，所述波至集群包括第一波至集群和第二波至集群，所述第一波至集群包括第一网格点处的一个或多个地震能波至，所述第二波至集群包括第一网格点处的一个或多个地震能波至；
在一个或多个处理器上确定各个波至集群的聚合数据，使得第一波至集群的聚合数据反映包括在第一波至集群中的第一网格点处的每个地震能波至的参数，以及第二波至集群的聚合数据反映包括在第二波至集群中的第一网格点处的每个地震能波至的参数；以及
在一个或多个处理器上实现波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体，
其中，实现第一波至集群的聚合数据和第二波至集群的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体。
7.如权利要求6所述的方法，其中，地震能波至被模拟成波束，以及其中，获得地震能波至包含获得感兴趣地质体内的粗网格点处的各个地震能波至的波束参数。
8.如权利要求6所述的方法，进一步包含：在一个或多个处理器上确定与各个地震能波至集群相对应的特征波至，其中，确定特征波至包含确定特征波至的参数。
9.如权利要求8所述的方法，其中，确定第一地震能波至集群的聚合数据基于为第一地震能波至集群确定的特征波至的参数，以及其中，确定第二地震能波至集群的聚合数据基于为第二地震能波至集群确定的特征波至的参数。
10.如权利要求8所述的方法，其中，在围绕第一网格点的细网格点处成像感兴趣地质体包含：
在一个或多个处理器上实现(i)为第一地震能波至集群确定的特征波至的参数，和(ii)第一地震能波至集群的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体；以及
在一个或多个处理器上实现(i)为第二地震能波至集群确定的特征波至的参数，和(ii)第二地震能波至集群的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体。
11.一种处理与感兴趣地质体相关联的地震数据的计算机实现方法，其中，该方法是在计算机系统中实现的，该计算机系统包含被配置成执行一个或多个计算机程序模块的一个或多个处理器，该方法包含：
将代表通过感兴趣地质体从一个或多个能量源传播到感兴趣地质体处或附近的一个或多个能量接收器的地震能的信息存储到一个或多个处理器可访问的电子存储器中；
在一个或多个处理器上获得位于感兴趣地质体内的粗网格点处的多个地震能波至的一个或多个参数，使得对于各个粗网格点，获得相应地震能波至组的参数，
其中，感兴趣地质体内的粗网格点包含第一网格点和第二网格点，其中，获得第一网格点处的地震能波至的一个或多个参数，以及其中，获得第二网格点处的地震能波至的一个或多个参数；
在一个或多个处理器上确定粗网格点处的波至的聚合数据，使得第一网格点处的地震能波至的聚合数据反映第一网格点处的每个地震能波至的参数，以及第二网格点处的地震能波至的聚合数据反映第二网格点处的每个地震能波至的参数；以及
在一个或多个处理器上实现所述波至集群的聚合数据，以在围绕粗网格点的细网格点处成像感兴趣地质体，
其中，实现第一网格点处的地震能波至的聚合数据以在围绕第一网格点的细网格点处成像感兴趣地质体，以及其中，实现第二网格点处的地震能波至的聚合数据以在围绕第二网格点的细网格点处成像感兴趣地质体。
12.如权利要求11所述的方法，其中，地震能波至被模拟成波束，以及其中，获得地震能波至包含获得感兴趣地质体内的粗网格点处的各个地震能波至的波束参数。
13.如权利要求11所述的方法，进一步包含：在一个或多个处理器上确定与粗网格点处的地震能波至相对应的特征波至，其中，确定特征波至包含确定特征波至的参数。
14.如权利要求13所述的方法，其中，确定第一网格点处的地震能波至的聚合数据基于为第一网格点处的地震能波至确定的特征波至的参数，以及其中，确定第二网格点处的地震能波至的聚合数据基于为第二网格点处的地震能波至确定的特征波至的参数。
15.如权利要求8所述的方法，其中，在围绕粗网格点的细网格点处成像感兴趣地质体包含：
在一个或多个处理器上实现(i)为第一网格点处的地震能波至确定的特征波至的参数，和(ii)第一网格点处的地震能波至的聚合数据，以在围绕第一网格点的细网格点处成像感兴趣地质体；以及
在一个或多个处理器上实现(i)为第二网格点处的地震能波至确定的特征波至的参数，和(ii)第二网格点处的地震能波至的聚合数据，以在围绕第二网格点的细网格点处成像感兴趣地质体。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/582,902 | 2009-10-21 | ||
| US12/582,902 US20110093203A1 (en) | 2009-10-21 | 2009-10-21 | System and method for clustering arrivals of seismic energy to enhance subsurface imaging |
| PCT/US2010/043860 WO2011049656A1 (en) | 2009-10-21 | 2010-07-30 | System and method for clustering arrivals of seismic energy to enhance subsurface imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102576088A true CN102576088A (zh) | 2012-07-11 |
Family
ID=43879963
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010800474012A Pending CN102576088A (zh) | 2009-10-21 | 2010-07-30 | 集群地震能波至以增强地下成像的系统和方法 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US20110093203A1 (zh) |
| EP (1) | EP2491427A1 (zh) |
| CN (1) | CN102576088A (zh) |
| AU (1) | AU2010308503B2 (zh) |
| BR (1) | BR112012008689A2 (zh) |
| CA (1) | CA2777319A1 (zh) |
| EA (1) | EA201270579A1 (zh) |
| WO (1) | WO2011049656A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459769A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震图像增强方法 |
| CN112859166A (zh) * | 2019-11-27 | 2021-05-28 | 中国石油天然气集团有限公司 | 地震数据处理方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5274605A (en) * | 1992-06-26 | 1993-12-28 | Chevron Research And Technology Company | Depth migration method using Gaussian beams |
| US6317695B1 (en) * | 2000-03-30 | 2001-11-13 | Nutec Sciences, Inc. | Seismic data processing method |
| US20050143921A1 (en) * | 2003-12-12 | 2005-06-30 | Exxonmobil Upstream Research Company | Method for seismic imaging in geologically complex formations |
| US7095677B2 (en) * | 2003-05-27 | 2006-08-22 | Paradigm Geophysical | Crossplot analysis of A.V.O. anomolies in seismic surveying |
| US7151857B2 (en) * | 2000-11-28 | 2006-12-19 | Monolith Co., Ltd. | Image interpolating method and apparatus |
| US7383133B1 (en) * | 2007-04-11 | 2008-06-03 | Pgs Onshore, Inc. | Diffuse seismic imaging systems and methods |
| US20090240438A1 (en) * | 2008-03-24 | 2009-09-24 | Chevron U.S.A. Inc. | System and method for migrating seismic data |
Family Cites Families (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3662325A (en) * | 1969-08-27 | 1972-05-09 | Western Geophysical Co | Method of displaying seismic data |
| US4703426A (en) * | 1984-11-08 | 1987-10-27 | Texas Instruments Incorporated | Method and apparatus for reviewing and modifying horizons of a three-dimensional seismic volume |
| USH482H (en) * | 1986-10-20 | 1988-06-07 | Seismic migration method | |
| US5012675A (en) * | 1989-07-25 | 1991-05-07 | Amoco Corporation | Integrating multiple mappable variables for oil and gas exploration |
| US5530679A (en) * | 1993-05-10 | 1996-06-25 | Western Atlas International, Inc. | Method for migrating seismic data |
| US6018497A (en) * | 1997-02-27 | 2000-01-25 | Geoquest | Method and apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore |
| US6574566B2 (en) * | 1999-12-27 | 2003-06-03 | Conocophillips Company | Automated feature identification in data displays |
| WO2002003099A2 (en) * | 2000-06-30 | 2002-01-10 | Exxonmobil Upstream Research Company | Method for imaging discontinuities in seismic data using dip-steering |
| FR2818387B1 (fr) * | 2000-12-18 | 2003-02-14 | Inst Francais Du Petrole | Methode pour obtenir des temps de trajet en reflexion a partir d'une interpretation de donnees sismiques en ondes cylindriques migrees |
| US7069149B2 (en) * | 2001-12-14 | 2006-06-27 | Chevron U.S.A. Inc. | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
| GB2398124B (en) * | 2003-02-08 | 2006-10-25 | Abb Offshore Systems Ltd | Estimating the time of arrival of a seismic wave |
| US7542373B2 (en) * | 2003-10-28 | 2009-06-02 | Baker Hughes Incorporated | Vector 3-component 3-dimensional kirchhoff prestack migration |
| US7480206B2 (en) * | 2004-09-13 | 2009-01-20 | Chevron U.S.A. Inc. | Methods for earth modeling and seismic imaging using interactive and selective updating |
| US7783060B2 (en) * | 2005-05-10 | 2010-08-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Deconvolution methods and systems for the mapping of acoustic sources from phased microphone arrays |
| WO2007056278A2 (en) * | 2005-11-03 | 2007-05-18 | Saudi Arabian Oil Company | Continuous reservoir monitoring for fluid pathways using 3d microseismic data |
| US7224642B1 (en) * | 2006-01-26 | 2007-05-29 | Tran Bao Q | Wireless sensor data processing systems |
| US7411863B2 (en) * | 2006-03-10 | 2008-08-12 | Westerngeco L.L.C. | Marine seismic data acquisition systems and methods |
| US8120991B2 (en) * | 2006-11-03 | 2012-02-21 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
| FR2909775A1 (fr) * | 2006-12-11 | 2008-06-13 | Inst Francais Du Petrole | Methode de construction d'un modele geologique d'une formation du sous-sol contraint par des donnees sismiques |
| US20080137480A1 (en) * | 2006-12-07 | 2008-06-12 | Woodside Energy Limited | Method of Building a Subsurface Velocity Model |
| US7773455B2 (en) * | 2006-12-13 | 2010-08-10 | Westerngeco L.L.C. | Time-lapse seismic acquisition |
| WO2008087505A2 (en) * | 2007-01-20 | 2008-07-24 | Spectraseis Ag | Time reverse reservoir localization |
| US20080225642A1 (en) * | 2007-03-16 | 2008-09-18 | Ian Moore | Interpolation of Irregular Data |
| US7646671B2 (en) * | 2007-12-11 | 2010-01-12 | Pgs Geophysical As | Method for processing marine towed streamer seismic data from regular multi-azimuth surveys |
| US8077540B2 (en) * | 2008-06-13 | 2011-12-13 | The United States Of America As Represented By The Secretary Of The Navy | System and method for determining vector acoustic intensity external to a spherical array of transducers and an acoustically reflective spherical surface |
| US20100054082A1 (en) * | 2008-08-29 | 2010-03-04 | Acceleware Corp. | Reverse-time depth migration with reduced memory requirements |
-
2009
- 2009-10-21 US US12/582,902 patent/US20110093203A1/en not_active Abandoned
-
2010
- 2010-07-30 EP EP10825354A patent/EP2491427A1/en not_active Withdrawn
- 2010-07-30 CA CA2777319A patent/CA2777319A1/en not_active Abandoned
- 2010-07-30 WO PCT/US2010/043860 patent/WO2011049656A1/en active Application Filing
- 2010-07-30 BR BR112012008689A patent/BR112012008689A2/pt not_active IP Right Cessation
- 2010-07-30 AU AU2010308503A patent/AU2010308503B2/en not_active Ceased
- 2010-07-30 EA EA201270579A patent/EA201270579A1/ru unknown
- 2010-07-30 CN CN2010800474012A patent/CN102576088A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5274605A (en) * | 1992-06-26 | 1993-12-28 | Chevron Research And Technology Company | Depth migration method using Gaussian beams |
| US6317695B1 (en) * | 2000-03-30 | 2001-11-13 | Nutec Sciences, Inc. | Seismic data processing method |
| US7151857B2 (en) * | 2000-11-28 | 2006-12-19 | Monolith Co., Ltd. | Image interpolating method and apparatus |
| US7095677B2 (en) * | 2003-05-27 | 2006-08-22 | Paradigm Geophysical | Crossplot analysis of A.V.O. anomolies in seismic surveying |
| US20050143921A1 (en) * | 2003-12-12 | 2005-06-30 | Exxonmobil Upstream Research Company | Method for seismic imaging in geologically complex formations |
| US7383133B1 (en) * | 2007-04-11 | 2008-06-03 | Pgs Onshore, Inc. | Diffuse seismic imaging systems and methods |
| US20090240438A1 (en) * | 2008-03-24 | 2009-09-24 | Chevron U.S.A. Inc. | System and method for migrating seismic data |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459769A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震图像增强方法 |
| CN104459769B (zh) * | 2013-09-22 | 2017-06-20 | 中国石油化工股份有限公司 | 一种地震图像增强方法 |
| CN112859166A (zh) * | 2019-11-27 | 2021-05-28 | 中国石油天然气集团有限公司 | 地震数据处理方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20110093203A1 (en) | 2011-04-21 |
| EA201270579A1 (ru) | 2012-09-28 |
| AU2010308503A1 (en) | 2012-03-29 |
| AU2010308503B2 (en) | 2013-12-05 |
| CA2777319A1 (en) | 2011-04-28 |
| EP2491427A1 (en) | 2012-08-29 |
| WO2011049656A1 (en) | 2011-04-28 |
| WO2011049656A8 (en) | 2012-04-05 |
| BR112012008689A2 (pt) | 2016-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105205441A (zh) | 用于从点云提取特征区域的方法和设备 | |
| US8279707B2 (en) | Fusing geophysical data representing a geophysical space | |
| CN102066978B (zh) | 用于偏移地震数据的系统和方法 | |
| CN102792186B (zh) | 具有吸收边界和随机边界的逆时迁移 | |
| WO2012030468A2 (en) | Reverse time migration back-scattering noise removal using decomposed wavefield directivity | |
| CN102395902A (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
| US8849574B2 (en) | System and method for perturbing an initial horizon-picking solution to follow local features of a volume | |
| CN102576088A (zh) | 集群地震能波至以增强地下成像的系统和方法 | |
| EP2496968B1 (en) | System and method for seismic beam formation that accounts for equipment misalignment | |
| CN102216810B (zh) | 利用基于射线的原理和有限元原理两者来推导地震波场的系统和方法 | |
| AU2011277060B2 (en) | Shot gather data beamer and debeamer | |
| CN102656481B (zh) | 衰减采集观测系统在地震数据中导致的假频的系统和方法 | |
| CN104662447B (zh) | 由亚最佳质量的地震图像来推测地层的系统和方法 | |
| CN104704393A (zh) | 用于产生地下目标的局部图像的系统和方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120711 |