考虑波束方位的用于衰减地震数据中的多次波的系统和方法

摘要

本发明合并模型驱动的和数据驱动方法的使用，以使用包括多次反射的表面相关的地震波的预测模型，减弱地震数据中的多次波。本发明包括考虑波束方位，并卷积预测的多次波波束与建模的短程层间波束的段以获取卷积的多次波波束的波束技术。然后，可以对卷积的多次波波束去卷积，以减弱存在于原始输入波束中的多次波。

说明书

技术领域

本发明涉及地震探测和处理，更具体而言，涉及利用波束的成像， 还涉及通过组合模型驱动的和数据驱动的方法来基于一次波预测多次 波的方法。

背景技术

在石油行业，通常使用地震勘探技术来帮助搜索和评估地下烃类 矿床。在地震勘探中，一个或多个地震能量源向诸如地质层之类的所 关注的地表下区域发射波。这些波进入地层中并可以分散，例如，通 过由地表下地震反射层面（即，具有不同的弹性属性的地下地层之间 界面）的反射或折射。通过一个或多个检波器采样或测量反射信号， 并记录所产生的数据。被记录的样本可以被称为地震数据或一组“地震 道”。可以分析地震数据，以提取正在被探测的地球的区域的结构的细 节和属性。

地震勘探包括三个单独的阶段：数据获取、数据处理和数据解释。 地震勘探操作的成功取决于所有三个阶段的令人满意的完成。

一般而言，地震探测的目的是通过向地面发射能量并记录从下面 的岩石层返回的“反射”或“回波”，形成地球的地表下的一部分（地层） 的图或图像。向地层发射的能量通常是声音，以及横波能量。向下传 播的声能可以来源于各种源，诸如陆地上的爆炸或地震振动器，或在 海洋环境下的气枪。地震探测通常使用一个或多个能量源以及通常还 使用大量的传感器或检测器。可以被用来探测返回的地震能量的传感 器通常是地震检波器（陆地勘测）或水听器（海上勘测）。

在表面地震勘探期间，能量源可以位于所关注的地质构造或地层 上面的地球表面附近的一个或多个位置，被称为爆炸点。每当激活源 时，源生成穿过地球向下传播的地震信号，并至少部分地从地表下的 各种类型的不连续面反射，包括从“岩石层”边界反射。一般而言，每 当地表下材料的弹性属性有变化时，都可以发生地震信号的部分反射。 反射的地震信号被发射回地球的表面，在此，它们作为若干个位置处 的传播时间的函数而被记录。返回信号被数字化，并被记录为时间的 函数（振幅与时间）。

在数据获取阶段由检波器记录的地震能量所存在的一个流行问 题是，地震道常常包含所希望的地震反射（“一次”反射）和会遮蔽或 淹没主要地震反射的不希望有的多次反射。一次反射是从源传递到检 波器的声波，带有从地表下地震反射层面的单次反射。多次反射是在 被检波器检测到之前反射了至少三次（向上，向下和再次向上）的波。 取决于与它们相关联的主要事件的时间延迟，多次波通常表征为暗示 与它们自己的一次反射干扰的短路径，或表现为单独的事件的长路径。

还有各种当前技术已知的多次波事件。存在有在两个强反射体 （自由表面和水底层）之间的水层中被“捕获”的信号。存在有“短程层 间”多次波事件，它们是以紧随发射之后或就在检测之前通过水层的额 外的往返行程为特征的反射。存在有“剩余的”表面相关的多次波事件， 其中，第一和最后一个向上的反射低于第一（水）层，并且，在两者 之间的自由表面上有至少一个反射。还有“间层”多次波，它们具有从 地表下反射体发生的向下反射。

在大多数情况下，多次波不包含不会更轻松地从一次波中提取的 任何有用信息。此外，水底层多次波已经被认为是许多近海地区中的 地震数据处理中的最严重的噪声问题。多次波会严重地掩盖用于结构 成像的一次反射事件，并污染振幅与偏移（“AVO”）信息。因此，去 除多次波，或至少衰减多次波是许多环境中的地震数据处理阶段的必 需的部分，特别是在其中多次波相对于一次波特别强的海的环境中。

在深水数据的情况下，对一阶和下面几阶海底多次波和短程层间 （peg-leg）反射的抑制是非常重要的。这些相当强的多次波可以与目 标反射体的一次反射具有相同传播时间。

取决于所使用的多次波的属性，有多个现有技术的减弱多次波的 方法。一类多次波衰减方法是预测方法，其中，从它们相应的一次波 预测多次波。现有技术的预测多次波衰减技术一般可以被分成两个类 别；模型驱动的方法和数据驱动的方法。模型驱动的方法一般使用地 球模型和已记录的数据，使用估计的海底反射率函数和计算出的振幅 函数，来预测或模拟多次波，以建模水层多次反射，然后，从原始数 据中减去这些预测的多次波。其他模型驱动的技术使用地球模型或反 射率模型以预测固定的多次波。数据驱动的方法利用这样的事实：一 次波和多次波通过卷积关系在物理上相关联，并通过交叉卷积被视为 包含对于多次波的固定的贡献的相关一次波束，来预测多次波。数据 驱动的方法可以一般性地处理复杂的几何形状，并需要很少或不需要 有关地表下的属性的信息。与数据驱动的技术相比，基于模型的技术 通常是成本效益好的，而后者通常更灵活。

某些模型驱动的方法一开始就要求结构信息，即，有关地表下结 构的信息，确定该结构信息是执行地震探测的理由。其他模型驱动的 方法要求震源子波的形状，由于反射和频率带宽限制，该形状将不会 是纯的δ函数。某些模型驱动的方法要求结构和震源子波信息，而其 他方法使用匹配过滤器来考虑扭曲的震源子波。

数据驱动的方法依赖于来自主要成份的多次波的可预测性。事实 上，该方法使用现有的地震数据来生成多次波，然后，从现有数据中 减去那些所生成的多次波。一个这样的数据驱动的现有技术方法被称 为“表面相关的多次波消除”或“SRME”。简而言之，此方法通过利用 现有的数据来操作，以创建只包含存在于数据中的多次波的预测的数 据集。具体而言，该方法试图预测多次波的地震表达，并在适应数据 中的现有多次波之后，从原始数据中减去预测的多次波，只留下初始 能量（至少理论上）。

数据驱动的SRME技术是用于在复杂地质条件下预测多次波的 有吸引力的解决方案，它们不要求地表下的任何先验知识（反射率、 结构和速度）。然而，这些方法要求对于每一个检波器位置的一个炮 点位置，对于大多数三维（“3D”）获取几何形状，情况不是这样的。 由于比较大的炮点间距，窄的传播长度和/或宽线缆间距，SRME方法 一般受到复杂3D多次波的挑战。可以从现有数据内插或外推缺失的 数据，但是，内插或外插在由大的炮点间距和/或检波器间距所导致的 假的地震数据方面有麻烦。高级的内插或外推法也是难以实现和昂贵 的。挑战3D SRME方法的这些复杂的3D数据的常见的原因是含盐顶 部的皱纹。但是，任何类型的复杂的超负载都会导致难以插入的复杂 的三维地震数据。

另一种数据驱动的方法使用预测反卷积，这是一种假设多次波是 周期性的而一次波不是周期性的过滤方法。对于来自小于500毫秒的 水体深度（大致1,200英尺）和大致分层的地下地质的数据，通常满 足此假设。在水体深度大于500毫秒的区域（其中，一次波和多次波 之间的速度差是显著的），可以使用诸如tau-p和f-k过滤之类的速度 -过滤方法（而不是预测方法），其中，变量f表示频率，k表示波数， p表示射线参数，而tau表示零偏移截距时间。

然而，过滤方法一般要求对地表下介质中的表观波传播速度的确 定，或至少按经验的推测，反射的地震波通过地表下介质从震源传递 到检波器。由于地表下结构的变化和岩石属性的组合，这些速度会千 差万别。另外，由于难以分离多次波和一次波束，预测反卷积常常会 导致对一次波束的非故意的损坏。此外，预测反卷积常常不能考虑反 射率中一般由短程层间多次波所导致的非线性因素。

发明内容

本发明的一个方面涉及用于衰减多次反射的、表面相关的地震波 的计算机实现的方法，其中，该方法在包括被配置成执行一个或多个 计算机程序模块的一个或多个处理器的计算机系统中实现。在一个实 施例中，该方法包括获取从涉及地质块体的宽方位地震数据导出的波 束数据集，其中，所述波束数据集包括通过所述地质块体来描述地震 数据的多个波束的属性的信息，并且其中，由所述波束数据集所描述 的所述多个波束的属性包括方位；从由所述波束数据集为其描述属性 的所述多个波束中选定输入波束；确定对应于所述输入波束的短程层 间波至的属性，其中，所述短程层间波至是部分地基于描述所述输入 波束的属性的所述波束数据集中的信息而确定的，所述输入波束的属 性包括所述输入波束的方位，并且其中，所确定的所述短程层间波至 的属性包括所述短程层间波至的方位；标识对应于所述短程层间波至 的一次波束，其中，标识所述一次波束包括标识由所述波束数据集所 描述的具有对应于所确定的所述短程层间波至的属性的属性（包括对 应于所述短程层间波至的波束的波束方位）的波束；以及，调整描述 所述输入波束的属性的所述波束数据集中的信息，以从所述输入波束 中去除与所述一次波束相关联的一个或多个多次波波至。

本发明的另一方面涉及被配置成减弱多次反射的、表面相关的地 震波的计算机系统。在一个实施例中，该系统包括被配置成执行一个 或多个计算机程序模块的一个或多个处理器。一个或多个计算机程序 模块可以包括数据集模块、输入波束模块、短程层间模块、一次波束 模块，以及调整模块。数据集模块被配置成获取从涉及地质块体的宽 方位地震数据导出的波束数据集，其中，所述波束数据集包括通过所 述地质块体来描述地震数据的多个波束的属性的信息，并且其中，由 所述波束数据集所描述的所述多个波束的属性包括方位。输入波束模 块被配置成从波束数据集描述了其属性的多个波束中选定输入波束。 短程层间模块被配置成确定对应于输入波束的短程层间波至的属性， 其中，所述短程层间波至是部分地基于描述所述输入波束的属性的所 述波束数据集中的信息而确定的，所述输入波束的属性包括所述输入 波束的方位，并且其中，所述短程层间模块进一步被配置，使得所确 定的所述短程层间波至的属性包括所述短程层间波至的方位。一次波 束模块被配置成标识对应于短程层间波至的一次波束，其中，所述一 次波束模块被配置成通过标识由所述波束数据集所描述的具有对应于 所确定的所述短程层间波至的属性的属性来标识所述一次波束，所确 定的所述短程层间波至的属性包括对应于所述短程层间波至的波束方 位的波束方位。调整模块被配置成调整描述输入波束的属性的波束数 据集中的信息，以从输入波束中去除与一次波束相关联的一个或多个 多次波波至。从输入波束去除多次波波至提供了基于输入波束（例如， 和从其中去除了多次波波至的其他输入波束一起）对所关注的地质块 体进行成像时的增强的准确性和/或精确性。

本发明的再一个方面涉及用于衰减多次反射的表面相关的地震 波的计算机实现的方法，其中，该方法在包括被配置成执行一个或多 个计算机程序模块的一个或多个处理器的计算机系统中实现。在一个 实施例中，该方法包括获取从涉及地质块体的宽方位地震数据导出的 波束数据集，其中，所述波束数据集包括通过所述地质块体来描述地 震数据的多个波束的属性的信息，并且其中，由所述波束数据集所描 述的所述多个波束的属性包括方位；从由所述波束数据集为其描述属 性的所述多个波束中选定输入波束；确定对应于所述输入波束的第一 短程层间波至的属性，其中，所述第一短程层间波至是部分地基于描 述所述输入波束的属性的所述波束数据集中的信息而确定的，所述输 入波束的属性包括所述输入波束的方位，并且其中，所确定的所述第 一短程层间波至的属性包括所述第一短程层间波至的方位；标识对应 于第一短程层间波至的第一一次波束，其中，标识所述第一一次波束 包括标识由所述波束数据集所描述的具有对应于所确定的所述第一短 程层间波至的属性的属性的波束，所确定的所述第一短程层间波至的 属性包括对应于所述第一短程层间波至的波束方位的波束方位；至少 部分地基于由所述波束数据集所描述的所述第一一次波束的所述属 性，确定关于所述第一一次波束的第一模型短程层间波至的属性，其 中，所述第一模型短程层间波至的属性包括所述第一模型短程层间波 至的方位；将所述第一模型短程层间波至与所述第一一次波束卷积， 以确定第一卷积的多次波波束的属性，其中，所确定的第一卷积的多 次波波束的属性包括方位；确定对应于所述输入波束的第二短程层间 波至的属性，其中，所述第二短程层间波至是部分地基于描述所述输 入波束的属性的所述波束数据集中的信息而确定的，所述输入波束的 属性包括所述输入波束的方位，并且其中，所确定的所述第二短程层 间波至的属性包括所述第二短程层间波至的方位；标识对应于第二短 程层间波至的第二一次波束，其中，标识所述第二一次波束包括标识 由所述波束数据集所描述的具有对应于所确定的所述第二短程层间波 至的属性的属性的波束，所确定的所述第二短程层间波至的属性包括 对应于所述第二短程层间波至的波束方位的波束方位；至少部分地基 于由所述波束数据集所描述的所述第二一次波束的所述属性，确定关 于所述第二一次波束的第二模型短程层间波至的属性，其中，所确定 的所述第二短程层间波至的属性包括所述第二模型短程层间波至的方 位；将所述第二模型短程层间波至与所述第二一次波束卷积，以确定 第二卷积的多次波波束的属性，其中，所确定的第二卷积的多次波波 束的属性包括方位；累积多个卷积的多次波波束，其中，所述多个卷 积的多次波波束包括所述第一卷积的多次波波束和所述第二卷积的多 次波波束；调整描述所述输入波束的属性的所述波束数据集中的信息， 以考虑所述累积的多个卷积的多次波波束。

本发明的这些及其他目标、特点，以及特征，以及结构的相关元 件以及部件的组合的操作方法和功能，以及制造成本的节约，通过下 面的参考附图对所附的权利要求进行的详细描述，将变得更加显而易 见，所有的这些附图构成了此说明书的一部分，其中，相同的引用编 号表示各个图中的对应的部分。然而，应该明确地理解，附图都只用 于说明和描述，不作为对本发明的限制。如在说明书和权利要求中所 使用的，单数形式也包括多个指称对象的对象，除非上下文明确地特 别指出。

附图说明

图1示出了本发明的用于衰减地震数据中的多次波的一个实施例 的流程图。

图2示出了包括将一次波转换为预测的多次波波束的本发明的一 个实施例。

图3示出了本发明的一个实施例，其中，从检测器位置发射射线， 并在海洋环境的水体底部表面反射。

图4示出了本发明的一个实施例，其中，通过射线跟踪，确定对 应于给定短程层间的一次波束。

图5示出了由本发明的一个实施例用来确定固定的短程层间的方 法。

图6示出了本发明的一个实施例所使用的上面的确定固定的短程 层间的方法的示意图。

图7示出了本发明的一个实施例，其中，一次波正在与短程层间 波束卷积。

图8示出了被本发明的一个实施例使用的单个波束解卷积的示意 图。

图9示出了本发明的用于衰减地震数据中的多次波的另一实施例 的流程图。

图10示出了本发明的用于衰减地震数据中的多次波的再一个实 施例的流程图。

图11示出了本发明的用于衰减地震数据中的多次波的一个实施 例的流程图。

图12示出了本发明的用于衰减地震数据中的多次波的另一实施 例的流程图。

图13示出了本发明的一个实施例，其中，在没有射线跟踪的情 况下，确定对应于短程层间的一次波束。

图14示出了根据本发明的一个实施例的被配置成减弱地震数据 中的多次波的系统。

具体实施方式

图1示出了从对所关注的地质块体的成像中去除多次波的方法 10。下面所呈现的方法10的操作只是说明性的。在某些实施例中，方 法10可以利用没有描述的一个或多个附加的操作完成，和/或没有所 讨论的操作中的一个或多个地完成。另外，图1中所示出并且在下面 描述方法10的操作的顺序并不是限制性的。

在某些实施例中，方法10可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法10的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法10 的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多个设 备。

如图1所示，方法10包括操作12，在此，初始化对应于所关注 的地质块体的地球模型。在操作14中，选择从所关注的地质块体中收 集的地震数据导出的波束数据集。地震数据是宽方位数据或其中方位 信息是属性的其他数据类型。如此处所使用的，“方位”或“方位信息” 是指标识地震能量在所关注的地质块体内的传播的水平方向的方向性 朝向的信息。例如，“方位”可以是指源和检测器之间的线和水平平面 中的某些基准方向之间的角度。波束数据集包括通过所关注的地质块 体来描述地震数据的多个波束的属性的信息。由波束数据集所描述的 多个波束的属性包括方位。由波束数据集所描述的其他属性可包括中 央射线路径的中点和偏移坐标，中央射线位置处的倾角、与中央射线 相关联的地震道序列值中的一项或多项，和/或其他属性。

分别在操作16、18以及20中选择来自在操作14中导出的波束 数据集的输入波束，来自在操作12中初始化的地球模型的多次波生成 表面，以及时间门。多次波生成表面可以以三维来描述。

通过使用来自操作16的输入波束，来自操作18的多次波生成表 面，以及来自操作20的时间门，在操作22中确定固定的短程层间。 如下面进一步讨论的，确定固定的短程层间包括沿着所选输入波束的 路径执行射线跟踪，以确定部分地由多次波生成表面处的反射而导致 的固定的短程层间的路径。确定固定的短程层间包括确定固定的短程 层间的属性，包括固定的短程层间的方位。

然后，在操作24中获取对应于固定的短程层间的一次波束。获 取一次波束包括从波束数据集中标识带有对应于固定的短程层间的属 性的属性的波束。在获取一次波束时为对应关系检查的属性之一是方 位。

在操作26中，获取涉及一次波束的建模的短程层间波束。建模 短程层间波束包括确定将在在18操作中选定的多次波生成表面中通 过一次波束（在操作24中获取的）的反射所生成的短程层间波束的属 性。

在操作28中，建模的短程层间波束（在操作26中获取的）与一 次波束（在操作24中获取的）卷积，以确定卷积的多次波波束。然后， 在操作30中，卷积的多次波波束与输入波束去卷积，以去除输入波束 中的多次波。如此处所使用的，对波束的卷积以“确定”或“获取”波束 将被理解为包括执行对已卷积的波束的卷积，导致“已确定的”或“已获 取的”波束的属性的确定。

在一个实施例中，不是在操作24中从波束数据集获取一次波束， 而是可以基于在操作22中确定的固定的短程层间的属性，来预测一次 波束的属性，并在操作28中，可以直接将预测的一次波束与建模的短 程层间波束进行卷积。

在一个实施例中，从检测器位置发射射线，并基于射线之一，确 定固定的短程层间。这样的实施例中的检测器位置基于在操作16中选 定的输入波束。

在一个实施例中，在操作28中，建模的短程层间波束的段与一 次波束卷积，以获取卷积的多次波波束。建模的短程层间波束的段可 以在多次波生成表面中开始，并可以涉及时间间隔，该时间间隔要么 是显式地指定的，要么是由地球模型中的另一个层位或者多次波生成 表面确定的。然后，在操作30中，可以对卷积的多次波波束去卷积， 以去除输入波束中的多次波。

用于从波束中去除多次波的传统技术对于地震道使用了本地倾 斜叠加或其他倾角区别方法，以将记录的波场分离为波束分量。本发 明使用本地倾斜叠加来将记录的波场分离为在位置和倾角方面局部化 的分量。这些分量是将被记录在如图2所示的位置A 36和B 38处的 波束波至32和34的中心处的东西。到达位置B 38处的波束能量32 在被反射面包围的层40中混响。例如，层40可包括水层。作为这种 混响的结果，假设波束能量32作为位置A 36处记录的波束34内的多 次波42而到达。将波束32移动射线跟踪传播时间TAB 44，从位置B 38 到位置A 36，将会把波束32中的事件与波束34中的多次波46和48 对齐。一旦对齐了事件，并标识了多次波46和48，就可以去除多次 波46和48。

当上文所描述的局部倾斜叠加用于3D获取时，记录的能量不能 完全转向到波束，因为波场沿着所有记录方向没有被密集地采样，如 此，有缺失数据的问题或需要准确地确定射线跟踪传播时间T_AB的严 重的别名。例如，如本领域相关技术人员所理解的，局部倾斜叠加是 在共同的偏移域，而不是在共中心点域里实现的。在本发明的一个实 施例中，作出了层叠速度描述了共中心点（“CMP”）集合中主要事件 的倾角的假设。一般而言，层叠速度是来自地表下盐层构造上方的地 质构造的一次波束的合理描述，可以在海洋环境中混响，以变为最强 的多次波。

图4、5和6示出了计算用于确定固定的短程层间的射线跟踪传 播时间所采用的方式。选择带有源位置“S”50和检测器位置“D”52的 输入波束，如图3所示。确定多次波生成表面，在此实施例中，它是 水体底部（“wb”）54。在其中方位数据将被实现以确定固定的短程层 间的情况下，确定多次波生成表面包括确定由多次波生成表面所构成 的三维平面。射线56从检测器D 52发射，它们被从多次波生成表面 wb 54反射，并到达自由表面58。射线56的角间隔被预先确定为在x 方向（从自由表面58向下垂直）为dp_x，在y方向（与板垂直）为dp_y。 选择带有给定射线参数或检测器射线路径倾角p_d的单个射线60，射 线60的波至位置Q 62被确定为在自由表面58上，如图4所示。射线 60的波至方向（射线参数p_q′64)被确定为在波至位置Q 62。确定位置 Q 62处的反射，并计算传出射线参数p_q 66。位置S 50以及Q 62，以 及它们的射线参数p_s 68和p_q 66，以及多次反射面wb 54的形状，确 定对应于带有位置D 52和Q 62和射线参数p_d 60和p_q′64的一次波束。 通过考虑所选输入波束的方位，增强了用来确定位置Q 62的对点D 52 和Q 62之间的波束的选择。

虽然图4和5是以两个维度示出的，但是，可以理解，所描述的 过程也考虑了对应于带有位置D 52和Q 62和射线参数p_d 60和p_q′64 的短程层间的一次波束的三维确定，以短程层间的方位的确定完成。

求出描述位置Q 62处的固定的射线路径要求射线路径搜索。执 行搜索，以求出在各种位置S 50、Q 62和D 52处发生的单独的反射。 在一个实施例中，搜索比较ray traced p_h与计算出的p_h（p_h是偏移倾 角）。射线跟踪的p_h被计算为：

(1)p_h=p_q-p_s

(2)p_m=p_q+p_s

其中，给定的位置处的p_m（中点倾角）对应于特定波束。从正 常moveout公式中获取计算出的p_h：

$\left(3\right)p_h=\frac{\partial t}{\partial h}=\frac{\partial }{\partial h}\sqrt{T_0^2+{\left(\frac{{2h\cos \theta }_m}{V}\right)}^2}$

其中，t是时间，h是半偏移，V是NMO速度，T₀是零偏移传播 时间，并且，$\cos \theta =\sqrt{1-{(\frac{1}{2}V·p_m)}^2}.$

V和T₀是从从地球模型中获取的层叠速度导出的。图5示出了对 p_h的确定，其中，p_h是h-t图74中的曲线72的局部斜率70。带有射 线跟踪的p_h和计算出的p_h之间的最近的匹配的波束被选为固定的短 程层间。从地球模型中获取的层叠速度可以是在方位方面依赖的。

图6示出了其中确定固定的短程层间的方法76。下面所呈现的方 法76的操作只是说明性的。在某些实施例中，方法76可以利用没有 描述的一个或多个附加的操作完成，和/或没有所讨论的操作中的一个 或多个地完成。另外，图6中所示出并且在下面描述方法76的操作的 顺序并不是限制性的。

在某些实施例中，方法76可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法76的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法76 的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多个设 备。

在操作78中，选择对应于给定源位置S和给定检测器位置D的 中点倾角p_m和方位。在操作80中，从检测器位置D选择试验p_d，该 试验p_d和在操作78中选择的p_m一起，根据公式（2）中给出的关系， 固定p_s的值。然后，此值可以被用来根据公式（1）中给出的关系， 计算p_h。在操作80中选择的参数p_d被用来执行从检测器位置D到在 操作82中在源位置S和检测器位置D之间选择的点Q的射线跟踪。 射线跟踪导致射线跟踪的偏移倾角p_h′和方位在操作84中，将它 们与计算出的偏移倾角p_h和方位进行比较。如果射线跟踪的p_h′紧密 地匹配计算出的p_h并且方位匹配计算出的方位，那么，此过程完成， 在操作86中，对应于p_m和方位的波束被选为固定的短程层间。如 果射线跟踪的p_h′不匹配计算出的p_h和/或射线跟踪的方位不匹配计 算出的方位，那么，方法76返回到操作80，选择另一个试验方向p_d和方位直到获得了射线跟踪的p′_h，方位和计算出的p_h，方位之 间的令人满意的匹配。此过程将确定检测器一侧的短程层间。可以理 解，源和检测器的角色可以互换，以确定源一侧的短程层间。

一个实施例包括卷积预测的多次波波束与建模的波束的段，以获 得卷积的多次波波束。如图7所示，预测的多次波波束B 86与在多次 波生成表面90中开始的建模的短程层间波束Q 88的一段卷积。卷积 的结果是波束E 92，它是卷积的多次波波束，该卷积的多次波波束与 输入波束A 94去卷积，以去除存在于输入波束A 94中的多次波。

在图8中提供了此实施例的示意图，其中，波束B 86被用作源 一侧预测96，而波束C 88被用于检测器一侧预测98。此实施例使用 维纳滤波器100和来自波束A 94、B 86和C 88的输入，来生成对存 在于波束A 94中的多次波的估计。维纳滤波器在当前技术中是已知 的。一般而言，它是因果滤波器，取决于某些约束，它将输入尽可能 地紧密地转换为所期望的输出。如本领域相关技术人员所理解的，有 其他可以执行此特定功能的滤波器或装置，它们旨在落入本发明的范 围内。一旦确定了波束A 94中的多次波，就可以从波束A 94中去除 102多次波，以产生输出波束104。

图9示出了用于从输入波束中去除多次波波束的方法106。下面 所呈现的方法106的操作只是说明性的。在某些实施例中，方法106 可以利用没有描述的一个或多个附加的操作完成，和/或没有所讨论的 操作中的一个或多个地完成。另外，图9中所示出并且在下面描述方 法106的操作的顺序并不是限制性的。

在某些实施例中，方法106可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法106的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法 106的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多 个设备。

在操作108中，初始化相互关联到所关注的特定地质块体的地球 模型。在操作110中，还启动根据所关注的地质块体的地震数据确定 的波束数据集。地震数据是宽方位数据或其中方位信息是属性的其他 数据类型。波束数据集包括通过所关注的地质块体来描述地震数据的 多个波束的属性的信息。由波束数据集所描述的多个波束的属性包括 方位。由波束数据集所描述的其他属性可包括中央射线路径的中点和 偏移坐标，中央射线位置处的倾角、与中央射线相关联的地震道序列 值中的一项或多项，和/或其他属性。

从图9可以看出，方法106包括一系列循环，其中，分别在操作 112、114以及116中选择来自波束数据集的输入波束，地球模型中的 多次波生成表面，以及时间门。在操作118中，从基于在操作112中 选择的输入波束的检测器或源位置发射射线。在操作120中，通过将 包括方位的射线的属性与所选输入波束的属性进行比较，从在操作 118中发射的射线中的一个中是固定的短程层间。在操作122中获取 对应于固定的短程层间的一次波束。获取一次波束包括确定一次波束 的属性，包括一次波束的方位，一次波束的中央射线路径的中点和偏 移坐标，中央射线位置处的一次波束的倾角，与一次波束的中央射线 相关联的地震道序列值，和/或基于短程层间的属性的一次波束的其他 参数。在操作124中，通过对应于在操作120中选择的固定的短程层 间的射线的传播时间的漂移，将一次波束变换成预测的多次波波束。 然后，在操作126中生成涉及预测多次波波束的建模的短程层间波束。 生成预测的多次波波束包括确定预测的多次波波束的属性，包括预测 的多次波波束的方位，预测的多次波波束的中央射线路径的中点和偏 移坐标，预测的多次波波束在中央射线位置处的倾角、与预测的多次 波波束的中央射线相关联的地震道序列值，和/或预测的多次波波束的 其他参数。在操作128中，在多次波生成表面开始的建模的短程层间 波束的段与预测的多次波波束卷积，以获取卷积的多次波波束。然后， 在操作130中，卷积的多次波波束与输入波束累积或者去卷积，以去 除输入波束中的多次波。累积的波束可以被用来将多次波预测重构为 地震道，或被用来以后利用输入波束去卷积。

在操作130中，卷积之后的若干个不同的点中任何一个都可以被 用来累积卷积的多次波波束或者去卷积已卷积的多次波波束，以去除 输入波束中的多次波。取决于正在被处理的数据，这些步骤可以在用 于选定时间门的For loop 116-132的结束之前或在该循环116-132之 后立即发生。这些步骤也可以在用于选定多次波生成表面的For loop 114-134之后或在用于选定输入波束的For loop 112-136之后发生。

图10示出了用于从输入波束中去除多次波波束的方法138。下面 所呈现的方法138的操作只是说明性的。在某些实施例中，方法138 可以利用没有描述的一个或多个附加的操作完成，和/或没有所讨论的 操作中的一个或多个地完成。另外，图10中所示出并且在下面描述方 法138的操作的顺序并不是限制性的。

在某些实施例中，方法138可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法138的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法 138的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多 个设备。

在操作140中，初始化相互关联到所关注的特定地质块体的地球 模型。在操作142中，获取根据所关注的地质块体的地震数据确定的 波束数据集。地震数据是宽方位数据或其中方位信息是属性的其他数 据类型。波束数据集包括通过所关注的地质块体来描述地震数据的多 个波束的属性的信息。由波束数据集所描述的多个波束的属性包括方 位。由波束数据集所描述的其他属性可包括中央射线路径的中点和偏 移坐标，中央射线位置处的倾角、与中央射线相关联的地震道序列值 中的一项或多项，和/或其他属性。

从图10可以看出，方法138包括一系列循环，其中，分别在操 作144、146以及148中选择来自波束数据集的输入波束，地球模型中 的多次波生成表面，以及时间门。在操作150中，从基于在操作144 中选择的输入波束的检测器或源位置发射射线。

在操作152中，从选定固定的短程层间周围的最窄的范围中选择 短程层间。在操作154中获取对应于选择的短程层间的一次波束。获 取一次波束包括确定一次波束的属性，包括一次波束的方位，一次波 束的中央射线路径的中点和偏移坐标，中央射线位置处的一次波束的 倾角，与一次波束的中央射线相关联的地震道序列值，和/或基于选择 的短程层间的属性的一次波束的其他参数。在操作156中，将一次波 束变换成预测的多次波波束。在操作158中获取涉及预测的多次波波 束的建模的短程层间波束。然后，在操作160中，预测的多次波波束 与建模的短程层间波束的段卷积，以获取卷积的多次波波束。操作160 所使用的建模的短程层间波束的段在多次波生成表面开始，并在检测 器或源位置处结束。在用于选定短程层间的For loop 152-162内，围 绕固定的短程层间的最窄的范围的预测的或卷积的多次波波束是层叠 的。在操作164，卷积的多次波波束累积或去卷积，以去除输入波束 中的多次波。操作164可以在用于选定短程层间的For loop 152-162， 用于选定时间门的For loop 148-166，用于选定多次波生成表面的For  loop 146-168，或用于选定输入波束144-170的For loop之后发生。

图11示出了用于从输入波束中去除多次波波束的方法172。下面 所呈现的方法172的操作只是说明性的。在某些实施例中，方法172 可以利用没有描述的一个或多个附加的操作完成，和/或没有所讨论的 操作中的一个或多个地完成。另外，图11中所示出并且在下面描述方 法172的操作的顺序并不是限制性的。

在某些实施例中，方法172可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法172的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法 172的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多 个设备。

在操作174中，初始化相互关联到所关注的特定地质块体的地球 模型。在操作176中，获取根据所关注的地质块体的地震数据确定的 波束数据集。地震数据是宽方位数据或其中方位信息是属性的其他数 据类型。波束数据集包括通过所关注的地质块体来描述地震数据的多 个波束的属性的信息。由波束数据集所描述的多个波束的属性包括方 位。由波束数据集所描述的其他属性可包括中央射线路径的中点和偏 移坐标，中央射线位置处的倾角、与中央射线相关联的地震道序列值 中的一项或多项，和/或其他属性。

从图11可以看出，方法172包括一系列循环，其中，分别在操 作178、180，以及182中选择来自波束数据集的输入波束，地球模型 中的多次波生成表面，以及时间门。在操作184中，从基于在操作178 中选择的输入波束的检测器或源位置发射射线。

在操作186中，选择短程层间。在操作186中选定的短程层间不 与固定的短程层间关联。在操作188中，短程层间被用来获取相对应 的一次波束。获取一次波束包括确定一次波束的属性，包括一次波束 的方位，一次波束的中央射线路径的中点和偏移坐标，中央射线位置 处的一次波束的倾角，与一次波束的中央射线相关联的地震道序列值， 和/或基于所选短程层间的属性的一次波束的其他参数。在操作190中， 将一次波束变换成预测的多次波波束，并在操作192中，从预测飞多 次波波束获取建模的短程层间波束。在操作194中，在多次波生成表 面开始的建模的短程层间波束的段与预测的多次波波束卷积，以获取 卷积的多次波波束。在操作198中，在For loop短程层间186-196 内所生成的卷积的多次波波束被累积和/或去卷积，以去除输入波束中 的多次波。操作198可以在用于选定短程层间的For loop 186-196， 用于选定时间门的For loop 182-200，用于选定多次波生成表面的For  loop 180-202，或用于选定输入波束的For loop 178-204之后发生。

可以为各种多次波，确定固定的短程层间和短程层间。可以为源 一侧多次波，确定固定的短程层间或短程层间。可以为检测器一侧多 次波，确定固定的短程层间或短程层间。可以为源一侧和检测器一侧 多次波两者，确定固定的短程层间或短程层间。

如本领域相关技术人员所理解的，可以有其中地球模型不能轻松 地可用的情况，在这样的情况下，本发明仍能预测和减弱多次波。作 为非限制性示例，图12示出了从不包括地球模型的使用的输入波束中 去除多次波的方法206。下面所呈现的方法206的操作只是说明性的。 在某些实施例中，方法206可以利用没有描述的一个或多个附加的操 作完成，和/或没有所讨论的操作中的一个或多个地完成。另外，图12 中所示出并且在下面描述方法206的操作的顺序并不是限制性的。

在某些实施例中，方法206可以以一个或多个处理设备来实现（例 如，数字处理器、模拟处理器、被设计成能处理信息的数字电路、被 设计成能处理信息的模拟电路、状态机，和/或用于以电子方式处理信 息的其他机制）。一个或多个处理设备可包括响应于以电子方式存储 在电子存储介质上的指令来执行方法206的某些或全部操作的一个或 多个设备。一个或多个处理设备可包括通过专门设计用于执行方法 206的操作中的一个或多个的硬件、固件、和/或软件配置的一个或多 个设备。

在操作208中，获取根据所关注的地质块体的地震数据确定的波 束数据集。地震数据是宽方位数据或其中方位信息是属性的其他数据 类型。波束数据集包括通过所关注的地质块体来描述地震数据的多个 波束的属性的信息。由波束数据集所描述的多个波束的属性包括方位。 由波束数据集所描述的其他属性可包括中央射线路径的中点和偏移坐 标，中央射线位置处的倾角、与中央射线相关联的地震道序列值中的 一项或多项，和/或其他属性。

从图12可以看出，方法206包括一对循环，其中，分别在操作 210和212中选择来自波束数据集的输入波束和时间门。在操作214 中，对固定的预测进行无先验的确定。例如，可以通过假设层状地球 模型或确定源和检测器或源位置的覆盖区，来在操作214中作出对固 定的预测的确定。在操作216中利用对固定的预测的这种确定。在操 作218中，获取对应于短程层间的一次波束。在操作220中，计算对 应于一次波束的建模的短程层间波束。在操作222中，利用建模的短 程层间波束，卷积一次波束。在操作226中，在For loop短程层间 216-224内所生成的卷积的多次波波束被累积和/或去卷积，以去除输 入波束中的多次波。操作226可以在用于选定短程层间的For loop 216-224、用于选定时间门的For loop 212-228，或用于选定输入波束 的For loop 210-230之后发生。

现在参考图13，上文所描述的实施例对于位置Q 232的范围以 及射线参数值P_q 234、P_s 236的范围以及方位，卷积多个一次波束和 短程层间波束。通过输入数据的分析，预定Q 232的范围和范围P_q234，P_s 236（或Pd 238），以及方位。范围是如此确定的，以便包括 足够的波束，以包含在固定的贡献与其他非固定贡献相加之后放大的 建模的多次波的固定的贡献。在图13中，选择带有给定方位的波束， 射线参数P_m=P_s236+P_d238，位置Q 232，以及射线参数P_q 234和P_ss 236（或P_d 238）。确定具有一次波束方位并且P_m1=P_s236+P_q234的 一次波束，以及具有短程层间波束方位并且P_m2=P_d238+P_q234的短程 层间波束。在此实施例中，作出了Q 232处的表面反射从P_q 234到–P_q240的假设。然后，一次波束可以与短程层间波束卷积。在此实施例 中，对于从Q 232=(S 242+D 244)/2或某种其他预先确定的位置开始 的位置Q 232值的范围，卷积一次波和短程层间波束。另外，对于射 线参数P_q 234值范围，卷积一次波和短程层间波束。

如上文所描述的，图10-13中所描绘的本发明的各实施例合并了 循环，这些循环示出了取决于正在被处理的数据，可以重复那些实施 例的某些步骤。

上文所描述的本发明的某些实施例包括从检测器或源位置发射 射线，以选择固定的短程层间或短程层间，应该理解，还有选择固定 的短程层间或短程层间的其他手段，例如，随机选择，这些手段也被 视为在本发明的范围内。

图14示出了被配置成通过所关注的地质块体来处理表示地震能 量的传播的地震数据的系统246。地震能量通过所关注的地质块体从 所关注的地质块体上或附近的一个或多个源位置传播到所关注的地质 块体在上或附近的一个或多个检测器或源位置。在对地震数据进行处 理时，系统246将地震能量建模为波束（例如，高斯波束）。系统246 可以被配置成执行图1，6以及9-12中所示出的，并如上文所描述的 方法10、76、106、138、172和/或206中的一个或多个。由系统246 所执行的处理从所选输出波束中去除多次波。在一个实施例中，系统 246包括电子存储器248、用户界面250、一个或多个信息资源252、 一个或多个处理器254，和/或其他组件。

在一个实施例中，电子存储器248包括以电子方式存储信息的电 子存储介质。电子存储器248的电子存储介质可包括与系统246整体 地提供的（即，基本上不可移动）的系统存储器和/或通过，例如，端 口（例如，USB端口、火线端口等等）或驱动器（例如，磁盘驱动器 等等）可移动地可连接到系统246的可移动存储器中的一种或两种。 电子存储器248包括光学可读存储介质（例如，光盘、等等），以磁 性方式可读的存储介质（例如，磁带、磁性硬盘驱动器、软盘驱动器 等等）、基于电荷的存储介质（例如，EEPROM、RAM等等）、固 态存储介质（例如，闪存驱动器等等）和/或其他以电子方式可读的存 储介质中的一种或多种。电子存储器248可以存储软件算法、由处理 器254所确定的信息、通过用户界面250接收到的信息、从信息资源 252接收到的信息，和/或可使系统246正确地运转的其他信息。电子 存储器248可以是系统246内单独的组件，或者电子存储器248可以 与系统246的一个或多个其他组件（例如，处理器254）整体地提供。

用户界面250被配置成在系统246和用户之间提供接口，通过该 接口，用户可以向系统246提供信息并从系统246接收信息。这可使 数据、结果，和/或指令以及任何其他能传递的项目（统称为“信息”） 在用户和系统246之间传递。如此处所使用的，术语“用户”可以表示 单个个人或一组可以互相配合的个人。适于包括在用户界面250中的 接口设备的示例包括小键盘、按钮、开关、键盘、旋钮、手柄、显示 屏幕、触摸屏、扬声器、麦克风、指示灯、可听警报，以及打印机。 在一个实施例中，用户界面250实际包括多个单独的接口。

可以理解，其他通信技术，无论是硬线连接的还是无线的，也被 本发明预期为用户界面250。例如，本发明预期，用户界面250可以 与由电子存储器248所提供的可移动存储器接口集成。在此示例中， 信息可以被从可使用户自定义系统246的实现的可移动存储器（例如， 智能卡、闪存驱动器、可移动磁盘等等）加载到系统246中。适于作 为用户界面250与系统246一起使用的其他示例性输入设备以及技术 包括，但不仅限于，RS-232端口、RF链路、IR链路、调制解调器（电 话、电缆或其他）。简而言之，用于与系统246传递信息的任何技术 也都可以被本发明预期为用户界面250。

信息资源252包括涉及所关注的地质块体和/或生成所关注的地 质块体的图像的过程的一个或多个信息源。作为非限制性示例，信息 资源252中的一个可包括在所关注的地质块体上或其附近获取的地震 数据，从其中导出的信息，和/或涉及获取过程的信息。地震数据可包 括地震数据的单个的道，或在从源通过所关注的地质块体传播的地震 能量的一个沟道中记录的数据。从地震数据导出的信息可包括，例如， 速度模型、波束数据集，包括与用来建模地震能量通过所关注的地质 块体的传播的波束相关联的波束属性，和/或其他信息。涉及地震数据 的获取的信息可包括，例如，涉及地震能量的源的位置和/或朝向的数 据，地震能量的一个或多个检测器的位置和/或朝向，能量被源生成并 被引向所关注的地质块体的时间，和/或其他信息。

处理器254被配置成在系统246中提供信息处理能力。如此，处 理器254可包括数字处理器、模拟处理器、被设计成能处理信息的数 字电路、被设计成能处理信息的模拟电路、状态机中的一种或多种， 和/或用于以电子方式处理信息的其他机制。虽然处理器254在图14 中被示为单个实体，但是，这只用于说明性目的。在一些实现中，处 理器254可包括多个处理单元。这些处理单元可以在物理上位于同一 个设备或计算平台内，或者，处理器254可以表示协调地操作的多个 设备的处理功能。

处理器254可以被配置成执行一个或多个计算机程序模块。一个 或多个计算机程序模块可包括地球模型模块256、波束数据集模块 258、输入波束模块260、表面模块262、时间门模块264、射线跟踪 模块266、短程层间模块268、一次波束模块270、调整模块272中的 一个或多个，和/或其他模块。处理器254可以被配置成通过软件；硬 件；固件；软件、硬件，和/或固件的某种组合；和/或用于配置处理 器254上的处理能力的其他机制，来执行模块256、258、260、262、 264、266、268、270和/或272。

应该理解，虽然在图14中模块256、258、260、262、264、266、 268、270和272被示为共同位于单个处理单元内，但是，在其中处理 器254包括多个处理单元的各实现中，模块256、258、260、262、264、 266、268、270和/或272中的一个或多个可以远离其他模块。对由下 面所描述的不同的模块256、258、260、262、264、266、268、270 和/或272所提供的功能的描述只是说明性的，并不是限制性的，因为 模块256、258、260、262、264、266、268、270和/或272中的任何 一个模块可以提供比所描述的功能多一些或少一些功能。例如，可以 省略模块256、258、260、262、264、266、268、270和/或272中的 一个或多个，某些或全部其功能可以由模块256、258、260、262、264、 266、268、270和/或272中的其他模块来提供。作为另一个示例，处 理器254可以被配置成执行一个或多个额外的模块，这些额外的模块 可以执行归属于模块256、258、260、262、264、266、268、270和/ 或272中的一个的某些或全部功能。

地球模型模块256被配置成获取涉及所关注的地质块体的初始化 地球模型。可以通过在过程254中初始化地球模型，来获取初始化地 球模型，或者，可以从电子存储器248中和/或外部信息资源252中的 一个获取初始化地球模型。在一个实施例中，地球模型模块256被配 置成执行操作12、108、140和/或174中的一个或多个（图1和9-11 中所示出的，并如上文所描述的）。

波束数据集模块258被配置成获取波束数据集。波束数据集可能 是从关于所关注的地质块体的宽方位地震数据或其中方位信息是属性 的其他数据类型导出的。波束数据集包括通过所关注的地质块体来描 述地震数据的多个波束的属性的信息。由波束数据集所描述的多个波 束的属性包括方位。由波束数据集所描述的其他属性可包括中央射线 路径的中点和偏移坐标，中央射线位置处的倾角、与中央射线相关联 的跟踪序列值中的一项或多项，和/或其他属性。波束数据集所描述的 波束可以是高斯波束。获取波束数据集可包括基于宽方位地震数据或 其中方位信息是属性的其他数据类型来确定波束数据集，和/或获取波 束数据集可包括从电子存储器248和/或外部信息资源252中的一个获 取波束数据集。在一个实施例中，波束数据集模块258被配置成执行 操作14、110、142、176和/或208中的一个或多个（图1和9-12所 示出的，并如上文所描述的）。

输入波束模块260被配置成从由波束数据集模块258所获取的波 束数据集描述了其属性的多个波束中选定输入波束。在一个实施例中， 输入波束模块260被配置成执行操作16、112、144、178和/或210中 的一个或多个（图1和9-12所示出的，并如上文所描述的）。

表面模块262被配置成选定由地球模型模块256所获取的地球模 型内的多次波生成表面。在一个实施例中，表面模块262被配置成执 行操作18、114、146和/或180中的一个或多个（图1和9-11所示出 的，并如上文所描述的）。

时间门模块264被配置成为所选输出波束和/或多次波生成表面 选定时间门。在一个实施例中，时间门模块264被配置成执行操作20、 116、148、182和/或212中的一个或多个（图1和9-12所示出的，并 如上文所描述的）。

射线跟踪模块266被配置成跟踪从给定检测器或源位置到另一表 面位置的射线。由射线跟踪模块266所跟踪的射线可包括在由表面模 块262所选定的表面和另一表面之间反射的射线。由射线跟踪模块266 从其中跟踪射线的给定检测器或源位置是对应于由输入波束模块260 所选定的输入波束的检测器或源位置。在一个实施例中，射线跟踪模 块266被配置成执行操作82、118、150和/或184中的一个或多个（图 7和9-11所示出的，并如上文所描述的）。

短程层间模块268被配置成确定对应于输入波束的一个或多个短 程层间波至的属性。短程层间波至是基于描述输入波束的属性（包括 输入波束的方位）、由射线跟踪模块266所跟踪的射线、由表面模块 262所选定的表面，固定的预测的波束数据集中的信息和/或其他因素 确定或选定的。短程层间波至的属性包括短程层间波至的方位。由短 程层间模块268所确定的一个或多个短程层间波至可包括固定的短程 层间和/或不固定的一个或多个短程层间s。在一个实施例中，短程层 间模块268可以执行操作22、120、152、186和/或216中的一个或多 个（图1和9-12所示出的，并如上文所描述的）。

一次波束模块270被配置成从由波束数据集模块258所获取的波 束数据集标识对应于由短程层间模块268所确定和/或选定的短程层间 的一次波束。此标识可以通过将波束数据集中所描述的波束的属性与 短程层间波至的属性进行比较来作出。可以比较某些属性，来查看是 否有匹配，而可以比较其他属性，以确定它们是否彼此对应。例如， 基于反射短程层间波至的表面的形状，一次波束模块270可以确定给 定波束的方位是否对应于短程层间波至的方位。

可以理解，在标识对应于短程层间波至的一次波束时，一次波束 模块270可能不会发现对短程层间波至的“完美的”匹配。对一次波束 的标识可以基于最佳拟合和/或满足一个或多个预定的准则的拟合来 作出。在一个实施例中，一次波束模块270被配置成执行操作24、122、 154、188和/或218中的一个或多个（图1和9-12所示出的，并如上 文所描述的）。

调整模块272被配置成调整描述输入波束的属性的波束数据集中 的信息，以从输入波束中去除与一次波束相关联的一个或多个多次波 波至。在一个实施例中，为实现此，调整模块272(i)从一次波束的属 性启动涉及一次波束的模型短程层间波至的属性，其中，所启动的属 性包括方位，(ii)将模型短程层间波至与一次波束卷积，以确定包括 卷积的多次波波束的方位的卷积的多次波波束的属性，以及，(iii)利 用输入波束来去卷积已卷积的多次波波束，以确定输入波束的属性， 包括方位，没有来自多次波波束的影响。在一个实施例中，调整模块 272被配置成执行操作26、28、30、126、128、130、158、160、164、 192、194和/或198中的一个或多个（图1和9-12所示出的，并如上 文所描述的）。

虽然为了公开基于了当前被认为是最实用和优选的实施例来详 细描述本发明的，但是，应所附权利要求书理解，这样的细节仅仅用 于所附权利要求书目的，本发明不仅限于所公开的实施例，相反，可 以涵盖在所附的权利要求的精神和范围内的修改和等效的方案。例如， 可以理解，本发明预期，在可能的程度上，任何实施例的一个或多个 特点可以与任何其他实施例的一个或多个特点相结合起来。