用于估计和减弱地震数据中的噪声的系统和方法

摘要

一种用于减弱代表感兴趣的地下区域的地震数据中的噪声的系统和方法包括：接收地震数据；将所述地震数据变换到其中所述地震数据具有稀疏的或可压缩的表示的域中，以生成变换的地震数据；将所述域分成窗口，其中所述窗口代表所述地震数据中已知的空间-时间位置；确定每个窗口中所述变换的地震数据的统计数据；基于所述变换的地震数据的统计数据确定针对每个窗口的滤波器；对每个窗口中所述变换的地震数据应用针对每个窗口的滤波器以生成滤波的变换的地震数据；以及执行滤波的地震数据的逆变换以生成滤波的地震数据。

说明书

技术领域

本发明一般涉及用于估计和减弱地震数据中的噪声的系统和方 法，并且具体涉及允许空间和时间变化适应可变噪声水平的用于估计 和减弱地震数据中的噪声的系统和方法。

背景技术

含烃储层的探测和开采可以借助于地震数据有效地进行，所述地 震数据必须被恰当处理以便允许对地下特征的解释。实际上，地震数 据经常被本质上可以相干或非相干(例如，随机)的噪声所污染。此 外，噪声水平可以在空间上和时间上这两者变化。

传统噪声抑制方法在估计和去除空间上和时间上变化的噪声方面 经常面临困难。传统方法在经常使用类似自动增益控制(AGC)的算 法的减弱步骤以前，可以试着跨地震数据归一化幅度。这可以在具有 强信号和弱噪声的区域中导致对信号的错误抑制。

需要用于估计和减弱地震数据中空间上和时间上变化的噪声的高 效和有效的方法，以改善最终地震图像并且允许对地下特征的恰当解 释。

发明内容

在此描述的是用于地震数据中噪声估计和减弱的计算机实施的方 法的各种手段的实施方式。

公开了用于减弱代表感兴趣的地下区域的地震数据中的噪声的计 算机实施的方法。所述方法包括：接收地震数据；将所述地震数据变 换到其中所述地震数据具有稀疏的或可压缩的表示的域中以生成变换 的地震数据；将所述域分成窗口，其中所述窗口代表所述地震数据中 已知的空间-时间位置；确定每个窗口中所述变换的地震数据的统计数 据；基于所述变换的地震数据的统计数据为每个窗口确定滤波器；对 每个窗口中所述变换的地震数据应用针对每个窗口的滤波器以生成滤 波的变换的地震数据；以及执行滤波的变换的地震数据的逆变换以生 成滤波的地震数据。

在一种实施例中，所述域可以是曲波(curvelet)域。在另一种 实施例中，所述域可以是小波域。

在一种实施例中，所述滤波器可以是阈值。

在一种实施例中，所述滤波的地震数据是噪声减弱的地震数据。 在另一种实施例中，所述滤波的地震数据是噪声模型，其可以接着被 从所述地震数据中减去以生成噪声减弱的地震数据。

在另一种实施例中，公开了一种计算机系统，包括数据源或存储 设备、用于实施减弱所述地震数据中的噪声的方法的至少一个计算机 处理器和用户接口。

在又一种实施例中，公开了一种制品，包括其上具有计算机可读 代码的计算机可读介质，所述计算机可读代码配置成实施减弱代表感 兴趣的地下区域的地震数据中的噪声的方法。

以上发明内容部分被提供以简单的形式介绍在以下的具体实施方 式部分中进一步描述的构思的选择。本发明内容不旨在指出所要求保 护的主题的关键特征或本质特征，也不旨在被用于限制所要求保护的 主题的范围。而且，所要求保护的主题不限于解决本公开的任何部分 所指出的任何或所有缺点的实施方式。

附图说明

参考以下说明书、权利要求和附图，本发明的这些或其他特征会 变得更好理解，其中

图1是说明根据本发明的一种实施例的方法的流程图；

图2示出了噪声合成地震数据集以及用于所述数据集的曲波系数 的直方图；

图3是在稀疏的或可压缩的域中的窗口的表示；

图4A使用合成数据以将本发明的一种实施例的结果与传统方法 比较；

图4B将本发明的一种实施例的结果中的误差与传统方法的结果 中的误差比较；

图5使用真实地震数据以将本发明的一种实施例的结果与传统方 法比较；以及

图6示意性说明用于执行根据本发明的一种实施例的方法的系 统。

具体实施方式

本发明可以被描述或实施在由计算机执行的系统和计算机方法的 一般上下文中。此类计算机可执行指令可以包括可用于执行特定任务 并且处理抽象数据类型的程序、例程、对象、组件、数据结构、和计 算机软件技术。本发明的软件实施方式可以针对多种计算平台和环境 中的应用被编码为不同语言。将要理解，本发明的范围和基础原理不 限于任何特定的计算机软件技术。

而且，本领域技术人员会理解本发明可以使用硬件和软件配置的 任意一种或组合来实践，包括但不限于具有单个和/或多个处理器计算 机、手持设备、可编程消费电子、迷你计算机、大型计算机等的系 统。本发明还可以在其中任务由通过一个或多个数据通信网络连接的 服务器或其他处理设备执行的分布式计算环境中。在分布式计算环境 中，程序模块可位于包括存储器存储设备的本地和远程计算机存储介 质这两者中。本发明还可被实践为井下传感器或测量设备的一部分或 实验室测量设备的一部分。

另外，与计算机处理器一同使用的制品，诸如CD、预记录的盘 或其他等价设备，可以包括其上记录的用于使计算机处理器促进本发 明的实践和实施的计算机程序存储介质和程序装置。此类设备和制品 也落入本发明的精神和范围内。

现在参考附图，本发明的实施例将被描述。本发明可以各种方式 被实施，包括例如作为系统(包括计算机处理系统)、方法(包括计 算机实施的方法)、装置、计算机可读介质、计算机程序产品、图形 用户界面、网页入口、或者有形地固定在计算机可读存储器中的数据 结构。本发明的若干实施例在以下被讨论。附图仅说明了本发明的典 型实施例并且因此不应被认为限制它的范围和广度。

本发明涉及估计和减弱地震数据中的噪声。本发明的一种实施例 被示为图1中的方法100。在此实施例中，在10处接收地震数据。仅 作为例子并且不是限制，地震数据可以来自陆上地震调查、海底地震 调查或合成地震数据。所述地震数据可以是1D数据、2D数据、2.5D 数据、3D数据和/或时延或4D数据。

所述地震数据在11处被变换到一个域中，在该域中变换的数据 具有稀疏的或可压缩的表示。在此域中，信号由相对小数量的显著系 数代表，而随机噪声由大量的小值的或零值的系数代表。此域的例子 包括小波域和曲波域。

图2包括噪声合成数据20的例子。所述数据包括5个线性事件 20A并且充满随机噪声。在变换为曲波域以后，创建代表所述变换的 数据的曲波系数的直方图21。曲波系数的近似中值被指示为22。在此 情况中，本领域技术人员会能够分辨所述中值代表具有大量小系数的 噪声23，而直方图的其余部分代表具有勉强可观察的数量的大系数的 信号24。所述噪声和信号存在某些交叠25。

所述数据被变换到的域被划分成窗口12。诸如所述划分的部分的 尺寸和形状之类的此划分的性质依赖于所述变换的数据的分布和域。 例如，如果所述域是曲波域，则所述划分可以是划分为每个尺度上在 每个有角楔形中的窗口。每个尺度代表了所述数据的带通滤波的版本 以及尺度内的每个有角楔形对应于有限范围的倾斜(dip)。说明了针 对单一尺度的8个有角楔形30A-30H的图3展示了这些。代表相同或 大体相同的空间-时间位置的窗口32A-32H被显示在每个楔形中。

一旦所述域已被划分，就可以在13处确定针对特定尺度的所有 有角楔形中的每个窗口的所述变换的地震数据的统计数据。所述统计 数据可以包括例如所述变换的数据的中值或其他分位点。针对所述窗 口的所述统计数据将与针对全体数据集合所见到的那些类似，如图2 中所示，意味着诸如中值之类的值可以表示噪声。针对每个窗口的所 述统计数据可以针对每个尺度来计算。

一旦每个窗口中的所述变换的数据的所述统计数据在方法100的 步骤13处被确定，就可以在方法100的步骤14处基于那些统计数据 确定滤波器。所述滤波器可以是阈值。通过在方法100的步骤15处应 用所述滤波器，所述滤波器可被用以从每个窗口中的所述变换的地震 数据中的所述信号中分离噪声。所述滤波器对所述域的每个窗口可以 不同，意味着所述滤波器可以在空间和时间这两者都变化。在每个窗 口都应用了它的滤波器后，在步骤16处所述滤波的数据被逆变换回到 原始域。所述滤波器可被设计为抑制所述噪声或者所述信号，这意味 着所述逆变换可分别产生信号模型或噪声模型。

图4A中可以看到方法100的结果。此处纯信号被显示为单个正 弦信号40。在左上角所述信号与相对弱的噪声结合并且到右下角增加 强度以产生合成输入数据集合41。使用跨越空间和时间恒定的阈值的 传统噪声减弱方法被应用于输入数据集合41以产生传统输出42。使 用在空间和时间变化的阈值实施的方法100已经被用于隔离并且抑制 所述噪声以产生新输出44。线43和45分别指示用于42和44的每个 轨迹的RMS幅度。通常，传统噪声减弱方法和方法100都被参数化以 从整个轨迹集合中移除同样量的能量。图4B显示了结果中的误差；在 46处可见传统方法的结果42与纯信号40的差别以及在47处可见本发 明的实施例的结果44与纯信号40的差别。因为传统方法在时间和空 间中不改变它的抑制，所以它比产生44的方法100的实施例移除更多 的信号46A。所述传统方法还显示出比方法100的实施例略微更高的 平均误差，如线48和49所示，线48示出了传统误差46的平均RMS 水平，并且线49示出了本发明得出的结果中的误差的平均RMS水 平。用于所述误差的每个轨迹的RMS幅度被示出为用于传统方法误差 46的线48A以及用于本发明的结果中的误差47的线49A。

图5在真实地震数据上将方法100的结果与传统噪声抑制方法比 较。地震数据50包括空间上和时间上变化的噪声并且具有若干事件和 包括区域50A的强信号与弱信号的区域。地震数据50被输入到图1的 方法100中。使用方法100移除的噪声被示出在图5的51中。地震数 据50还被输入到使用自动增益控制(AGC)的传统噪声抑制方法以在 尝试噪声抑制之前跨所述地震数据归一化幅度。由使用跨所有轨迹恒 定的阈值的传统方法移除的噪声可以在图5的53中可见。强信号区域 50A在51和53中被分别指示为区域50B和50C。在53中有显著更多 “噪声”能量，它实际上是已经被传统方法错误认定为噪声的强信号 的一部分。在传统方法中决定噪声减弱有多强的AGC归一化标量从信 号与噪声的和导出。由于在50A中所述信号强，因此这已经导致了在 53的区域50C中过度侵略性的噪声减弱。在另一方面，如在51的 50B中可见，本发明显示了对信号的存在几乎不敏感(强反射器)、 估计的信号是弱的、有效阈值是小的、以及强反射器是未损坏的。

图5中噪声的能量还可如用于本发明的线52以及用于传统方法 的线54所见。这也示出由传统方法移除的“噪声”在整个区域50B较 高。但是，如线55所示，由两种方法移除的噪声的平均能量是相同 的。

图6示意性地说明了用于执行图1的方法100的系统600。所述 系统包括数据源/存储设备60，其可以包括数据存储设备或计算机存储 器以及其他。设备60可以包含记录的地震数据和/或合成地震数据。 来自设备60的数据可以对诸如可编程通用目的计算机之类的处理器 61可用。处理器61被配置以执行实施方法100的计算机模块62。这 些计算机模块可以包括变换模块62A以执行方法100的步骤11，包括 划分窗口模块62B以执行步骤12，包括统计模块62C以执行步骤 13，包括确定模块62D以执行步骤14，包括应用模块62E以执行步骤 15以及逆变换模块62F以执行步骤16。所述系统可以包括诸如用户接 口63之类的接口组件。用户接口63可以被用以显示数据和处理的数 据产品以及允许用户从用于实施所述方法的方面的选项中作选择。仅 作为例子并且非限制，在处理器61上计算的移除的噪声和噪声减弱的 地震数据可以显示在用户接口63上，存储在数据存储设备或存储器 60上，或者既被显示又被存储。

尽管在上述说明中本发明已经根据它的某些优选实施例被描述， 并且许多细节已经被公开用于说明目的，但是对本领域技术人员将会 清晰的是，本发明是容易改变的并且此处描述的某些其他细节在不背 离本发明的基本原理下可以较大地变化。此外，应当理解，在此处的 任意一个实施例中描述或显示的结构特征或方法步骤也可以被其他实 施例使用。