抑制地震数据中的噪声
阅读说明：本技术 抑制地震数据中的噪声 (Suppressing noise in seismic data ) 是由 秦福豪 康斯坦丁诺斯·钦戈斯 于 2018-05-15 设计创作，主要内容包括：本公开描述了用于抑制地震数据中的噪声的方法和系统,包括计算机实现的方法、计算机程序产品和计算机系统。一种计算机实现的方法包括：在数据处理装置处接收与地下区域相关联的地震数据集合；由数据处理装置根据识别的地震事件对地震数据集合进行平坦化；由数据处理装置将地震数据集合划分为多个空间窗口；由数据处理装置根据随机次序顺序对地震数据集合进行随机化；由数据处理装置对随机化后的地震数据进行滤波；以及由数据处理装置根据随机化前的顺序对滤波后的地震数据进行重新组织。(The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for suppressing noise in seismic data. A computer-implemented method comprising: receiving, at a data processing device, a set of seismic data associated with a subsurface region; flattening, by the data processing device, the seismic data set according to the identified seismic event; dividing, by a data processing device, a seismic data set into a plurality of spatial windows; randomizing, by the data processing device, the seismic data set according to a random order sequence; filtering the randomized seismic data by a data processing device; and reorganizing, by the data processing device, the filtered seismic data according to the order prior to randomization.)
优先权声明
本申请要求于2017年5月17日递交的美国专利申请No.15/597,861的优先权，其全部内容通过引用并入本文。
技术领域
本公开涉及抑制地震数据中的噪声。
背景技术
在地球物理分析中，地震数据被收集并用于分析地理区域的地下地质结构和岩石特性。这些数据以及基于这些数据的分析对石油和天然气行业的勘探、生产和钻井作业至关重要。
发明内容
本公开描述了用于抑制地震数据中的噪声的方法和系统，包括计算机实现的方法、计算机程序产品和计算机系统。一种用于抑制地震数据中的噪声的计算机实现的方法，包括：在数据处理装置处接收与地下区域相关联的地震数据集合；由数据处理装置根据识别的地震事件来对地震数据集合进行平坦化；由数据处理装置将地震数据集合划分为多个空间窗口；由数据处理装置根据随机次序顺序将地震数据集合随机化；由数据处理装置对随机化后的地震数据进行滤波；以及由数据处理装置根据随机化前的顺序对滤波后的地震数据进行重新组织。
该方面的其他实现方式包括对应计算机系统、装置和记录在一个或多个计算机存储设备上的计算机程序，每个计算机程序被配置为执行所述方法的动作。一个或多个计算机的系统可以被配置为通过在系统上安装的、在操作中使得系统执行动作的软件、固件、硬件或软件、固件或硬件的组合来执行特定的操作或动作。一个或多个计算机程序可以被配置为通过包括指令来执行特定操作或动作，所述指令在被数据处理装置执行时使得该装置执行动作。
本说明书的主题的一个或多个实现方式的细节在附图和下面的描述中阐述。通过说明书、附图和权利要求书，所述主题的其它特征、方面和优点将变得显而易见。
附图说明
图1示出了根据实现方式的示例地震数据噪声抑制处理。
图2示出了根据实现方式的地震数据噪声抑制处理的示例效果。
图3示出了根据实现方式的示例叠前(pre-stack)时间偏移(migrated)的地震断面。
图4示出了根据实现方式的滤波后的地震数据的示例结果。
图5示出了根据实现方式的原始地震输入与提取的滤波后的地震断面之间的示例比较。
图6是根据实现方式的地球物理成像系统的高级架构框图。
各附图中相似的附图标记和标号指示相似的要素。
具体实施方式
给出以下描述以使得本领域任何技术人员能够实现和使用所公开的主题，并且在一个或多个特定实现方式的上下文中提供所述描述。对公开的实现方式的各种修改对本领域技术人员是显而易见的，并且在不脱离本公开的范围的情况下，本文所定义的一般原理可以适用于其他实现方式和应用。因此，本公开并非意不旨在限于所描述的和/或示出的实现方式，而应赋予与本文公开的原理和特征一致的最宽范围。
本公开总体上描述了用于抑制地震数据中的噪声的方法和系统，包括计算机实现的方法、计算机程序产品和计算机系统。在某些情况下，地震信号可以通过源设备发送到源位置处的地球的地下。地震信号的示例包括声学信号。地震信号行进通过地下，并且可以被置于接收器位置处的接收器设备接收。在某些情况下，源设备、接收器设备或其组合可以置于地表上。信号可以向下传播，直到到达反射结构并向上朝向地表反射。因为信号通过地下结构被折射和反射，因此所接收信号的特性包含地下结构的信息。所接收的信号可以被收集为地震数据。地震数据可以包括反射波的道(trace)。
在海上勘测中，***和水听器可以分别用作源设备和接收器设备。在采集期间，地震源从***阵列引爆。通过水听器的拖缆采集反射和折射的信号。在陆地采集中，可以使用***作为***源，并且采用地震检波器作为接收器设备。在另一示例中，可以使用振动卡车作为源设备。也可以使用生成和接收地震信号的其他设备。
在一些实现方式中，对地震数据应用偏移操作以对地下结构成像。地震偏移(Seismic migration)是这样一种处理：在空间或时间上以几何方式将地震事件重新定位到事件在地下发生的位置，而不是地表处记录事件的位置，从而创建更准确的地下图像。地震数据的偏移是通过对地震数据进行数值的基于网格的空间卷积来计算倾斜(dipping)事件(其中地质层不平坦)而对平坦地质层假设进行的校正。在偏移处理中，地震数据中包含的回波信息从记录时间转换为地下深度方面的特征。偏移处理使地震特征在横向和垂直方向上都处于其在空间中的位置。
可以使用不同的偏移算法。例如，在深度域中，基于射线和基于波束的基尔霍夫方法由于其计算效率和面向目标的灵活性而在实践中受到欢迎。之后，广泛使用了单向和双向波动方程偏移。在某些情况下，可以在分析偏移中使用反向时间偏移(RTM)。RTM算法包括对源波场的前向传播的计算、接收器波场的后向传播的计算以及两个计算出的波场之间的关联成像条件的计算。在某些情况下，源波场和接收器波场可以分别称为前向波场和后向波场。在时域中，成像通常通过共深度点(CDP)堆叠来实现，CDP堆叠假设平坦的地下结构，然后进行叠后(post-stack)偏移，或者进行类似于深度域中的基尔霍夫偏移的叠前时间偏移。然而，能量从源和接收器传递到成像点的行进时间不能像深度域那样通过射线追踪来计算。它基于正常时差(Normal Move Out)(NMO)速度和Dix公式，该公式假设速度模型由相对平坦和较厚的层形成。
在某些情况下，由于地震数据中的噪声，可能难以获得基于地震数据生成的可靠地震图像。例如，地震数据可能受到与地下结构相关联的层间多次波(internalmultiples)的影响。这些层间多次波可以是由被地理区域的地下中的多个反射界面反弹的地震能量所导致的。这些层间多次波可以被称为地震数据中的噪声。这些多重噪声中的一些噪声可以具有与主要地震事件的视倾斜(apparent dip)类似的视倾斜，因此可能不能通过堆叠或偏移算子而被抑制或区分。在某些情况下，例如在陆地地震勘测中，地震噪声还可以包括来自不完善的静态或近地表模型、不规则的炮点(shot)和接收器分布以及近地表衍射和散射的影响。
在一些实现方式中，可以基于对多次波波至(multiple arrivals)的预测，然后从获取的数据中将其减去，来抑制地震多重噪声。由于良好的数据质量和多重事件与主要事件之间相对较大的视倾斜区分力，这些预测和相减技术可以适用于深海海洋数据中地表相关的多次波。另一方面，由于一些地理区域中的近地表地质复杂性，预测和相减方案可能对陆地地震数据的效果不佳。结果是，来自堆叠或偏移的最终地震图像通常包括大量相干噪声。在某些情况下，噪声可能非常强，以至于可能损害实际的地震反射，从而使地震解释难以理解。
当实际地震事件与噪声的倾斜差很小时，将出现处理困难。在这种情况下，频率-波(F-K)滤波或其他噪声去除方法(例如，中值滤波器)可能无法将噪声与主要事件分开。
可以对叠前数据或叠后数据执行噪声抑制处理。在叠前数据应用中，可以将炮点道集或CDP道集中的地震事件排成一行并展平。可以执行数据随机化和噪声抑制处理以抑制噪声。然而，将叠前数据随机化的效果取决于该处理中使用的正常时差(NMO)速度的精确度。如果NMO速度不精确，则随机化可能无法对齐相干地震事件，且无法有效地抑制噪声。
在叠后数据应用中，对地震数据的最终图像施加操作。例如，数据可以是NMO堆叠、叠前时间偏移或叠前深度偏移结果。解释器将关于数据的事件识别为主要反射或噪声(例如，多重反射)。该识别可以基于一般的地质知识、已知的区域地质特性或现有的井信息。当解释器由于数据质量或缺少其他信息而无法对事件进行明确识别时，这些事件也可以是解释器的试解释场景(trial interpretation scenarios)。
一旦地震事件被识别并被映射，垂直压缩和拉伸地震数据，以便对所识别的事件进行平坦化，然后对平坦化后的数据进行顺序随机化。保留与所识别的事件平行的地震事件的相干性，同时降低其他事件的相干性。然后可以提取相干事件，同时可以使用滤波来抑制其他事件。这些其他事件也可以通过在原始数据和在随机化以及噪声抑制后的滤波后的数据之间进行相减运算来获得，以用于进一步分析。
图1示出了根据实现方式的示例地震数据噪声抑制处理100。为了清楚地呈现，下文的描述在图2至图6的上下文中一般性地描述处理100。然而，应当理解，处理100可以例如由任何其他合适的系统、环境、软件和硬件，或者系统、环境、软件和硬件的组合来执行，视情况而定。在某些情况下，处理100可以在大规模计算机集群、超级计算机或任何其他计算设备或计算设备的集合上执行。在一些实现方式中，处理500的各个步骤可以并行、组合、循环和/或以任意顺序运行。
在102处，在数据处理装置处接收与地下区域相关联的地震数据集合。该地震数据集合可以是针对地下区域获取的接收器信号数据集合。在某些情况下，在地震数据获取期间，源(***、振动卡车、***阵列等)被激活，来自地下地质边界的反射/折射/透射由位于地球表面上的接收器设备记录。顺序或同时针对每个炮点重复进行这种类型的获取，直到针对该勘测区域获取了所有地震数据。这些获取到的地震数据包括在接收器信号数据中。在某些情况下，所获取的地震数据在现场收集，被传送到办公室(经由计算机网络、物理网络或其组合来存储和传输)，并且被用作执行处理100的计算设备的输入。
在某些情况下，地震数据集合可以是叠后数据。例如，地震数据可以是地震数据的最终图像。叠后数据的示例包括NMO堆叠、叠前时间偏移结果和叠前深度偏移结果。
在104处，根据一个或多个识别的地震事件对地震数据集合进行平坦化。地震平坦化是一种用于去除诸如褶皱或断层之类的结构的解释技术，以在地质层沉积时揭示原始的地下结构。平坦化有助于解释器基于一个层位(horizon)来识别地质特征。在解释中，可以通过解释器、解释工具或其组合来识别一个或多个地震事件。可以基于所识别的事件对地震数据集合进行平坦化。所识别的事件的示例可以包括主要断裂，其中可以明确地匹配断裂两侧的事件。为了基于事件对数据进行平坦化，选择道位置(trace location)作为锚固点。所有其他地震道将根据锚固道进行调整。可以对地震道进行压缩或拉伸，以便所有识别的事件都水平排成一行。压缩和拉伸是通过使用原始未平坦化的数据在相邻的识别事件之间***数据来实现的。
在106处，将平坦化的数据划分为多个重叠的空间窗口。如果每个窗口的大小很小，则每个窗口中的相干事件可能是线性的或接近线性的，并且事件的幅度不会显著变化。最佳窗口大小因数据而异。在某些情况下，可以执行测试以确定窗口大小。
在108处，根据随机次序顺序重新布置地震数据集合。在一个实现方式中，可以基于与空间窗口相对应的空间距离来索引每个空间窗口中的数据。可以生成随机顺序。可以根据随机化顺序重新分配每个空间窗口中的数据。在一个示例中，随机化顺序为[4,7,3,1,2,6,5]，并且基于水平方向上的空间距离对地震数据集合进行随机化。在该示例中，基于空间窗口的水平距离为每个空间窗口指派从1至7的索引。在顺序随机化期间，空间距离被索引为4的空间窗口中的数据将重新分配到第一位置，随后是空间距离被索引为7的空间窗口中的数据，再之后是空间距离被索引为3的空间窗口中的数据，以此类推。
在110处，对重新布置的地震数据进行滤波。在某些情况下，可以应用中值滤波器来提取相干层位事件。中值滤波窗口大小通常比数据随机化窗口大小小得多。在某些情况下，可以执行测试以确定中值滤波窗口大小。通过仅在零K分量附近进行保持，FK滤波也可以用于提取层位事件。在某些情况下，可以执行测试以确定FK滤波窗口大小。也可以采用以已知的倾斜提取线性事件的其他技术。
在112处，基于滤波后的数据的随机化前的顺序对其进行重新组织。在步骤108处给出的示例中，在滤波之后，空间距离被索引为4的空间窗口中的数据从第一位置返回到第四位置。
在某些情况下，可以基于重新组织的数据生成地下区域的图像。图像可以由数据处理装置输出，发送给不同的数据处理装置或其组合。
图2示出了根据实现方式的地震数据噪声抑制处理的示例效果。图2包括合成图像202、204、206和208。图像202示出了已经对基于解释而识别的层位事件进行了平坦化的合成地震数据的示例集合。如图所示，图像202示出了平坦化后的层位事件210以及两个其他事件212和214。事件212具有向右的较小倾斜，事件214具有向左的较小倾斜。由于倾斜之间的差很小，因此直接滤波可能无法有效区分这些事件。
图像204示出了对地震数据进行随机化的影响。这里，图像202中的地震数据是在水平方向上随机重新布置的小窗口。基于解释识别的层位事件不受随机化影响，但会破坏倾斜事件的相干性。图像206示出了在使用中值滤波器对图像204中的地震数据进行滤波并根据其原始水平位置进行重新组织之后的噪声抑制效果。这里，所识别的层位事件不受影响，但倾斜事件在很大程度上被抑制。图像208示出了在滤波后的道返回到其原始位置之后对图像204进行F-K滤波的结果。图像208中的噪声残差与图像206中的噪声残差不同。
在某些情况下，如图像206和208所示，在滤波和重新组织之后，可能会因倾斜事件而存在残留能量。然而，残留能量不连续(如图像206所示)或更弱(如图像208所示)。可以使用附加的滤波处理来进一步抑制残留能量。如图像206和208所示，噪声抑制处理成功抑制了倾斜相干噪声，同时保持了解释中识别的层位事件不受影响。在某些情况下，可以执行多于一个的解释，每种解释都可以识别不同的事件。可以在迭代处理中基于不同的解释来应用本文所讨论的噪声抑制处理。在每次迭代中，根据所识别的事件，根据特定解释将地震数据平坦化。这些处理可以根据解释器的选择而生成更清晰的地震图像，并且帮助解释器测试不同的解释场景。
图3示出了根据实现方式的示例叠前(pre-stack)时间偏移(migrated)的地震断面。如图3所示，相干噪声干扰主要反射，这使解释变得困难。
在某些情况下，可以基于地质和井信息进行假设，假设该区域中可能不存在如此剧烈的构造变化和断层活动。因此，图3中断面的粗糙外观可能通常是由于存在相干噪声。在这些情况下，可以基于解释中识别的主要地震反射事件来对地震数据进行平坦化。与主要地震反射事件平行的其他反射事件也可以被平坦化。在所示示例中，1.5秒以下的地震事件用于平坦化断面，因为它们连续且相干。然后，将随机化和滤波过程应用于平坦化后的数据，并去除平坦化。图4示出了根据实现方式的滤波后的地震数据的示例结果。如图4所示，大多数层位是连续的，并符合当地的地质情况。解释器可以使用新的地震图像和该区域的井控技术，容易地追踪地震事件并导出整个断面上的地质结构。
在某些情况下，可以通过从未被抑制的数据(例如，图3所示的数据)中减去抑制后的数据(例如，图4所示的数据)来获得抑制后的事件。图5示出了根据实现方式的原始地震输入与提取的滤波后的地震断面之间的示例比较。如图所示，数据之间的差异主要包括与原始偏移图像上的噪声相对应的逐段事件。可以例如通过解释器对该差异进行分析，以确认这些事件是否反映地震数据中的噪声。
图6是根据实现方式的地球物理成像系统的高级架构框图。从高层面来看，所示的系统600包括与网络630耦接的地球物理图像处理计算机602。所描述的图示仅是所描述的主题的一种可能的实现方式，并且不旨在将本公开限制于该单个描述的实现方式。本领域普通技术人员应理解这样的事实，根据本公开，所描述的组件可以以备选的方式连接、组合或使用。
网络630促进计算机602与其他组件(例如，获得位置的观测数据并将观测数据发送给计算机602的组件)之间的通信。网络630可以是无线或有线网络。网络630还可以是存储器管道、硬件连接或组件之间的任何内部或外部通信路径。
计算机602包括被配置为执行如本文所述的方法的计算系统。在某些情况下，方法的算法可以用可执行计算代码(例如，C/C++可执行代码)来实现。在某些情况下，计算机602可以包括运行批量应用的独立的Linux系统。在某些情况下，计算机602可以包括具有足够的存储器大小以处理每个地球物理数据块的移动或个人计算机。计算机602可以用于实现本公开中讨论的噪声抑制处理。
计算机602可以包括包含输入设备(例如键区、键盘、触摸屏、麦克风、语音识别设备、能够接受用户信息的其他设备)和/或传达与计算机602的操作相关联的信息(包括数字数据、视觉和/或音频信息或GUI)的输出设备在内的计算机。
计算机602可以用作客户端、网络组件、服务器、数据库或其他持久性设备，和/或系统600的任何其他组件。在一些实现方式中，计算机602的一个或多个组件可以被配置为在基于云计算的环境中操作。
从高层面来看，计算机602是可操作用于接收、发送、处理、存储或管理与系统600相关联的数据和信息的电子计算设备。根据一些实现方式，计算机602还可以包括或通信地耦接到应用服务器、电子邮件服务器、Web服务器、缓存服务器、流传输数据服务器、商业智能(BI)服务器和/或其他服务器。
计算机602可以通过网络630从客户端应用(例如，在另一计算机602上执行的应用)接收请求，并通过在适当的软件应用中处理所述请求来响应所接收的请求。此外，还可以将请求从内部用户(例如，从命令控制台或通过另一适当的访问方法)、外部或第三方、其他自动化应用以及任何其他适当的实体、个人、系统或计算机发送给计算机602。
计算机602的每个组件可以使用系统总线603进行通信。在一些实现方式中，计算机602的任意和/或所有组件(硬件和/或软件两者)可以使用应用编程接口(API)612和/或服务层613通过系统总线603彼此对接和/或与接口604对接。API 612可以包括针对例程、数据结构和对象类的规范。API 612可以是独立于或依赖于计算机语言的，并且涉及完整的接口、单个功能或甚至是API集合。服务层613向计算机602和/或系统600提供软件服务。计算机602的功能可以对于使用该服务层的所有服务消费者是可访问的。软件服务(例如，由服务层613提供的软件服务)通过定义的接口提供可重用的、定义的业务功能。例如，接口可以是以JAVA、C++或以可扩展标记语言(XML)格式或其他合适格式提供数据的其他合适语言所编写的软件。虽然被示出为计算机602的集成组件，但是备选实现方式可以将API 612和/或服务层613示出为与计算机602和/或系统600的其他组件有关的独立组件。此外，在不脱离本公开的范围的情况下，API 612和/或服务层613的任意或所有部分可以被实现为另一软件模块、企业应用或硬件模块的子模块或副模块。
计算机602包括接口604。尽管在图6中示出为单个网络接口604，但是可以根据计算机602和/或系统600的特定需求、需要或特定实现方式，使用两个或更多个接口604。计算机602将接口604用于与连接到网络630(无论是否示出)的分布式环境中的其他系统(包括系统600内)通信。通常，接口604包括以合适的组合以软件和/或硬件编码的逻辑，并且可操作用于与网络630通信。更具体地，接口604可以包括支持与通信相关联的一个或多个通信协议的软件，使得网络630或接口的硬件可操作用于在所示出的系统600内部和外部传送物理信号。
计算机602包括处理器605。尽管在图6中示出为单个处理器605，但是可以根据计算机602和/或系统600的特定需求、需要或特定实现方式，使用两个或更多个处理器。通常，处理器605执行指令并操纵数据以执行计算机602的操作。具体地，处理器605执行处理地球物理数据所需的功能。
计算机602还包括存储器606，存储器606保存用于计算机602和/或系统600的其他组件的数据。尽管在图6中示出为单个存储器606，但是可以根据计算机602和/或系统600的特定需求、需要或特定实现方式，使用两个或更多个存储器。虽然存储器606被示出为计算机602的集成组件，但是在备选实现方式中，存储器606可以在计算机602和/或系统600的外部。
应用607是算法软件引擎，其根据计算机602和/或系统600的特定需求、需要或特定实现方式提供功能，特别是关于处理地球物理数据所需的功能。例如，应用607可以用作图1至图5中描述的一个或多个组件/应用。此外，尽管被示出为单个应用607，但是应用607可以被实现为计算机602上的多个应用607。此外，尽管被示出为与计算机602集成在一起，但是在备选实现方式中，应用607可以在计算机602和/或系统600的外部。
可以存在任意数量的与系统600相关联的或在其外部并通过网络630进行通信的计算机602。此外，在不脱离本公开的范围的情况下，术语“客户端”、“用户”和其他适当的术语可以适当地互换使用。此外，本公开考虑到许多用户可以使用一个计算机602，或者一个用户可以使用多个计算机602。
所描述的主题的实现方式可以单独或组合地包括一个或多个特征。例如，在第一实现方式中，一种用于抑制地震数据中的噪声的计算机实现的方法包括：在数据处理装置处接收与地下区域相关联的地震数据集合；由数据处理装置根据识别的地震事件对地震数据集合进行平坦化；由数据处理装置将地震数据集合划分为多个空间窗口；由数据处理装置根据随机次序顺序将地震数据集合进行随机化；由数据处理装置对随机化后的地震数据进行滤波；以及由数据处理装置根据随机化前的顺序对滤波后的地震数据进行重新组织。
先前和其他所述实现方式可以各自可选地包括以下特征中的一个或多个：
可与任何以下特征组合的第一特征，其中地震数据是叠后数据。
可与任何先前或以下特征组合的第二特征，其中地震数据包括叠前时间偏移结果或叠前深度偏移结果中的至少一个。
可与任何先前或以下特征组合的第三特征，其中根据对地下区域的解释来识别所识别的地震事件。
可与任何先前或以下特征组合的第四特征，还包括由数据处理装置基于重新组织的数据生成地下区域的图像。
可与任何先前或以下特征组合的第五特征，还包括由数据处理装置从地震数据中减去重新组织的地震数据的集合。
可与任何先前或以下特征组合的第六特征，还包括由数据处理装置对重新组织的地震数据进行滤波。
在第二实现方式中，一种存储指令的非暂时性计算机可读介质，当指令执行时，使计算设备执行如下操作，包括：在数据处理装置处接收与地下区域相关联的地震数据集合；由数据处理装置根据识别的地震事件对地震数据集合进行平坦化；由数据处理装置将地震数据集合划分为多个空间窗口；由数据处理装置根据随机次序顺序将地震数据集合进行随机化；由数据处理装置对随机化后的地震数据进行滤波；以及由数据处理装置根据随机化前的顺序对滤波后的地震数据进行重新组织。
先前和其他所述实现方式可以各自可选地包括以下特征中的一个或多个：
可与任何以下特征组合的第一特征，其中地震数据是叠后数据。
可与任何先前或以下特征组合的第二特征，其中地震数据包括叠前时间偏移结果或叠前深度偏移结果中的至少一个。
可与任何先前或以下特征组合的第三特征，其中根据对地下区域的解释来识别所识别的地震事件。
可与任何先前或以下特征组合的第四特征，操作还包括由数据处理装置基于重新组织的数据生成地下区域的图像。
可与任何先前或以下特征组合的第五特征，操作还包括由数据处理装置从地震数据中减去重新组织的地震数据的集合。
可与任何先前或以下特征组合的第六特征，操作还包括：由数据处理装置对重新组织的地震数据进行滤波。
在第三实现方式中，一种设备包括：至少一个硬件处理器；以及非暂时性计算机可读存储介质，耦接到至少一个硬件处理器并存储用于由至少一个硬件处理器执行的编程指令，其中编程指令在被执行时使至少一个硬件处理器执行包括如下各项的操作：在至少一个硬件处理器处接收与地下区域相关联的地震数据集合；由至少一个硬件处理器根据所识别的地震事件对地震数据集合进行平坦化；由至少一个硬件处理器将地震数据集合划分为多个空间窗口；由至少一个硬件处理器根据随机次序顺序将地震数据集合进行随机化；由至少一个硬件处理器对随机化后的地震数据进行滤波；以及由至少一个硬件处理器根据随机化前的顺序对滤波后的地震数据进行重新组织。
先前和其他所述实现方式可以各自可选地包括以下特征中的一个或多个：
可与任何以下特征组合的第一特征，其中地震数据是叠后数据。
可与任何先前或以下特征组合的第二特征，其中地震数据包括叠前时间偏移结果或叠前深度偏移结果中的至少一个。
可与任何先前或以下特征组合的第三特征，其中根据对地下区域的解释来识别所识别的地震事件。
可与任何先前或以下特征组合的第四特征，操作还包括由数据处理装置基于重新组织的数据生成地下区域的图像。
可与任何先前或以下特征组合的第五特征，操作还包括由数据处理装置从地震数据中减去重新组织的地震数据的集合。
可与任何先前或以下特征组合的第六特征，操作还包括由数据处理装置对重新组织的地震数据进行滤波。
本说明书中所描述的主题和功能操作的实现方式可以实现在下述形式中：数字电子电路、有形体现的计算机软件或固件、计算机硬件，包括在本说明书中公开的结构及其结构等同物，或它们中的一个或多个的组合。本说明书中所描述的主题的实现方式可以实现为一个或多个计算机程序，即，在有形的非暂时性计算机可读计算机存储介质上编码的计算机程序指令的一个或多个模块，以用于由数据处理装置执行或控制数据处理装置的操作。备选地或附加地，程序指令可以编码在人工生成的传播信号(例如，机器生成的电、光或电磁信号)中/上，生成所述信号以对信息进行编码以传输到合适的接收器装置，以供数据处理装置执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备或计算机存储介质的组合。
术语“实时”、“实时(快速)(RFT)”、“接近实时(NRT)”、“准实时”或类似术语(如本领域的普通技术人员所理解的)意味着动作和响应在时间上接近，使得个人感知动作和响应实质上同时发生。例如，在个人做出了访问数据的动作之后对数据显示的响应的时间差(或用于启动显示)可以小于1ms、小于1秒或小于5秒。虽然所请求的数据不需要被即时显示(或启动以显示)，但是考虑到所描述的计算系统的处理限制和例如收集、准确测量、分析、处理、存储或发送数据所需的时间，在没有任何有意的延迟的情况下对该数据进行显示(或启动以显示)。
术语“数据处理装置”、“计算机”或“电子计算机设备”(或本领域普通技术人员所理解的等同物)指代数据处理硬件，并且包括用于处理数据的各种装置、设备和机器，例如包括可编程处理器、计算机、或多个处理器或计算机。装置还可以是或还可以包括专用逻辑电路，例如，中央处理单元(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)小区、语音、视频、数据或其它合适信息(或通信类型的组合)。
计算系统可以包括客户端和服务器。客户端和服务器一般彼此远离并且通常通过通信网络进行交互。客户端和服务器的关系通过在相应计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生。
尽管本说明书包含许多特定的实现细节，然而这些细节不应被解释为对任何发明的范围或可以要求保护的事项的范围构成限制，而是用于描述特定于具体发明的具体实现的特征。在单个实现中，还可以组合实现本说明书中在独立实现的上下文中描述的某些特征。相反的，在单个实现的上下文中描述的各种特征也可以在多个实现中单独实现，或以任何合适的子组合来实现。此外，尽管前述特征可以被描述为在某些组合中起作用并且甚至最初如此要求保护，但是来自所要求保护的组合的一个或多个特征在一些情况下可以从组合中删除，并且所要求保护的组合可以针对子组合或子组合的变体。
已经描述了本主题的特定实现。对于本领域技术人员显而易见的是，所描述的实现的其他实现、改变和置换在随后的权利要求的范围内。尽管在附图和权利要求中以特定顺序描述了操作，但这不应被理解为：为了实现期望的结果，要求按所示的特定顺序或按相继的顺序来执行这些操作，或者要求执行所有所示的操作(一些操作可以看作是可选的)。在一些情况下，多任务或并行处理(或者多任务和并行处理的组合)可以是有利的并且视情况来执行。
此外，在前述的实现中的各种系统模块和组件的分离或集成不应被理解为在所有实施方式中要求这样的分离或集成，并且应该理解的是，所描述的程序组件和系统一般可以一起集成在单个软件产品中或封装为多个软件产品。
因此，前述示例实现不限定或限制本公开。在不脱离本公开的精神和范围的情况下，还可以存在其他改变、替换和变化。
此外，任何要求保护的实现被认为适用于至少一种计算机实现的方法；存储用于执行计算机实现的方法的计算机可读指令的非暂时性计算机可读介质；以及计算机系统，该计算机系统包括与硬件处理器可互操作地耦接的计算机存储器，所述硬件处理器被配置为执行计算机实现的方法或存储在非暂时性计算机可读介质上的指令。