专利名称：用于三维表面相关的多次波消除的动态孔径确定的方法
技术领域：
本发明一般涉及地球物理勘探领域。更具体而言，本发明涉及海洋地震(seismic)调查中的表面相关多次波消除(multiple elimination)的领域。
背景技术：
在石油和天然气工业中，一般使用地球物理勘探来帮助搜索和评估地下地层。地 球物理勘探技术提供地球的地下结构的知识，这对寻找和提取有价值的矿产资源、特别是 诸如石油和天然气等油气藏是有用的。地球物理勘探的一种众所周知的技术是地震调查。处理在执行地震调查时获得的所得到的地震数据以提供关于调查区域中的地下 地层的地质结构和性质的信息。处理所处理地震数据以便显示并分析这些地下地层的潜在 油气含量。地震数据处理的目的是从地震数据中提取尽可能多的关于地下地层的信息以便 充分地形成地下地质的图像。为了识别存在找到石油聚集的概率的地球地下的位置，在收 集、处理、和解释地震数据方面花费了大量的金钱。根据记录的地震数据来构造限定感兴趣 的地下地层的反射表面的过程以深度或时间为单位提供地球的图像。产生地球地下结构的 图像以便使得解释者能够选择具有最大的具有石油聚集的概率的位置。在海洋地震调查中，使用地震能源来生成地震信号，该地震信号在传播到地表中 之后至少部分地被地下地震反射物(seismic reflector)反射。此类地震反射物通常是具 有不同弹性特性、尤其是声波速度和岩石密度的地下地层之间的界面，所述不同的弹性特 性引起界面处的声阻抗的差异。由地震传感器(也称为地震接收器)来检测并记录被反射 的地震能量。在海洋地震调查中用于生成地震信号的适当震源通常包括用船牵引并被周期性 地激活以产生声波场的水中震源。产生波场的震源通常是气枪或气枪的空间分布阵列。适当类型的地震传感器通常包括安装在牵引地震等浮电缆(也称为地震电缆)内 的粒子速度传感器(在本领域中称为地震检波器)和水压力传感器(在本领域中称为水听 器)。可以自行部署地震传感器，但是更常见的是将其在等浮电缆内部署成传感器阵列。在反射波到达地震传感器之后，该波继续传播到水表面处的水/空气界面，在那 里波被向下反射，并再次被传感器检测到。反射波继续传播且可以再次被水底或其它地下 地层界面向上反射。不止一次被反射的反射波称为“多次波(multiple)”且通常被视为噪 声。特定类别的噪声包括从水表面反射至少一次并称为表面相关多次波的多次波。三维表面相关多次波消除(3D SRME)力图通常通过预测和减法过程来使表面相关 多次波衰减。首先根据地震数据来估计表面相关多次波并随后从地震数据减去预测的多次 波而留下衰减过噪声的信号。此过程中的第一步骤包括构造用于源-接收器轨迹(trace) 的多次波影响集合(multiple contribution gather)，其涉及计算成对轨迹在称为孔径 (aperture)的空间区域内的卷积。第二步骤包括构造主要包含多次波反射的预测的多次波 轨迹，其涉及将所有多次波影响轨迹堆叠在用于源-接收器轨迹的多次波影响集合中。第 三步骤包括从原始地震数据减去许多此类预测的多次波轨迹。
最容易实现的多次波影响集合的孔径的选择是具有预选纵向(inline)和横向 (crossline)维度的矩形空间区域，其以将对其进行多次波轨迹预测的源-接收器轨迹的 中点位置为中心。当至少多次波影响集合中的所有有影响(contributing)事件的顶点落 在孔径内时，孔径对于特定多次波影响集合而言是足够的。这意味着对应于顶点的表面反 射点落在孔径内且Fresnel堆叠用于在顶点时使有影响事件瓦解成相应的多次波事件。因 此，多次波影响集合的孔径的简单选择将不可避免地是成本与精确度之间的权衡。对于某 些多次波影响集合而言，孔径将太大且对于其它的而言将太小。尺寸和形状最佳的孔径将 是尽可能小以使计算工作量最小化且足够大以捕捉所有表面反射点的所有顶点，并且对于 相长干涉而言足够大以用于那些顶点的Fresnel堆叠。因此，需要一种用于高效地确定多次波影响集合的孔径的适当形状和尺寸的方 法。优选地，可以根据地震数据本身来动态地确定孔径的最佳形状和最佳尺寸
发明内容

本发明是一种用于使用多次波影响集合的动态孔径确定来将地震数据变换成衰 减过噪声的数据的方法。通过以下内容来预测表示地下地层的地震数据中的多次波轨迹。 在多次波影响集合的孔径边界附近，根据通过在要评估的地球地下区域附近部署多个地震 传感器获取的地震数据，计算一系列的相邻多次波影响轨迹组的倾角(dip)，所述地震传感 器响应于地震能量而生成电学和光学信号中的至少一个。基于所计算的倾角，递归地扩展 多次波影响集合的孔径的边界。将具有扩展孔径边界的多次波影响集合中的多次波影响轨 迹堆叠以生成预测的多次波轨迹。从地震数据减去预测的多次波轨迹以生成可用于形成地 下地层的图像的表面相关多次波衰减数据。


通过参照以下详细说明和附图，可以更容易地理解本发明及其优点，在附图中图1是本发明的方法中的用于源-接收器轨迹的多次波影响集合的构造的示意 图；图2是示出本发明的方法中的多次波影响集合的扩展孔径的示意图；图3是示出用于使用多次波影响集合的动态孔径确定来将地震数据变换成衰减 过噪声的数据的本发明的实施例的流程图；图4是示出多次波影响集合的动态孔径确定的本发明的另一实施例的流程图；图5是示出关于多次波影响集合的构造的本发明的要素的流程图；以及图6是示出关于递归孔径扩展的本发明的要素的流程图。虽然将结合本发明的优选实施例来描述本发明，但应理解的是本发明不限于此。 相反，本发明意图涵盖可以被包括在随附权利要求所定义的本发明范围内的所有替换、修 改、和等价物。
具体实施例方式本发明是一种用于使用多次波影响集合的动态孔径确定来将地震数据变换成衰 减过噪声的数据的方法。本发明可以在三维预测和减法过程中用来使地震数据中的表面相关多次波衰减。由源-接收器轨迹来定义将预测并减去多次波的位置。图1是本发明的方 法中的用于源-接收器轨迹的多次波影响集合的构造的示意图。由特定源位置11和特定 接收器位置12的组合来确定源_接收器轨迹10。根据用于此源-接收器轨迹10的地震数 据预测表面相关多次波的过程包括对成对轨迹求卷积并随后累计(accumulate)表面区域 13内的卷积结果。每对卷积轨迹中的第一轨迹14具有与源-接收器轨迹10和所选联接 (1 inkage)位置15处的接收器位置相同的特定源位置11，且该对卷积轨迹中的第二轨迹16 具有在与源_接收器轨迹10相同的联接位置15和相同的特定接收器位置12处的其源位 置。用于所有成对的第一轨迹14和第二轨迹16的联接位置15分布在表面区域13上。每 个联接位置15对应于与此源-接收器轨迹10相关的表面相关多次波的潜在表面反射点。 由成对的第一轨迹14和第二轨迹16的卷积产生的轨迹称为多次波影响轨迹并将这些多次 波影响轨迹(在累计之前)的系综(ensemble)称为多次波影响集合。在预测中结合(在 累计之后)的表面反射点的表面区域13的范围或覆盖范围统称为预测的孔径。图2是示出本发明的方法中的多次波影响集合的扩展孔径的示意图。多次波组成 集合20被示为在顶部处具有任意的空间维度21和在左侧向下进行的时间维度22。在某 一时间发生的属于源-接收器组合的源-接收器轨迹中的多次波事件将在多次波影响集合 20表现为在同一时间具有多次波影响集合20内的顶点24的表面23和稍后时间以下的表 面23其余部分。此外，顶点24的空间位置将对应于该特定多次波事件的表面反射点的空 间位置。多次波影响轨迹到预测的表面相关的多次波中的累计包括一般称为Fresnel堆叠 的过程。每个顶点24处和周围的相长干涉增强期望的预测的多次波且远离每个顶点24的 破坏性干涉减小信号的其余部分(即，不是预测的多次波)。多次波预测的质量取决于所有 相关表面反射点被多好地结合到多次波影响集合20的孔径中和Fresnel堆叠多好地对多 次波影响轨迹起作用。两者均取决于选择正确的孔径和孔径内的联接位置的间距。返回图1，可以将孔径选择为最简单的情况，矩形表面区域13围绕用于特定源位 置11和特定接收器位置12的源-接收器轨迹10的中点17。可以用纵向维度18和横向维 度19或替换地用纵向孔径和横向孔径来对矩形表面区域13进行参数化。而且，可以改变 矩形表面区域13的取向。然而，孔径的表面区域13不必是矩形区域。孔径的形状可以是 例如椭圆形或不规则形状。在任何情况下，孔径的最佳表面区域13应尽可能地小以使计算 工作量最小化并足够大到以捕捉所有顶点并足以使那些顶点的相长干涉起作用。联接位置 15之间的间距或距离的情况也是如此。返回图2，示出了多次波影响集合20的初始孔径25。本发明的方法对多次波影响 集合20的边界附近的轨迹的(重叠)部分采用倾角确定方法以产生扩展孔径26。如果发 现沿着跨越轨迹的事件计算的倾角27是正的，亦即事件本身指向上，意味着朝向较早的时 间，并向内，然后，相应顶点24将落在多次波影响集合20内部。如果所计算的倾角28是负 的，则相应顶点29将落在多次波影响集合20外面并应沿着边界的那些部分增大初始孔径 25以产生扩展孔径26。本发明的这种方法允许动态地确定扩展孔径26。因此，可以使适当 孔径区域的选择是完全或部分数据驱动的。从小初始孔径25开始，基于所计算的倾角的某 些度量，可以动态地扩展边界。作为一个示例，倾角 的初始度量可以是在应用倾角滤波器之后剩余的能量。更复 杂的度量还可以将倾角的相对振幅、时间、和斜率考虑在内。可以根据水底深度、其局部或全局变化或更具体而言地质复杂性的观察或先验知识来进行度量。对于典型的矩形孔径区域而言，通常将产生四个相关度量，每个边上一个。对于形状更复杂的区域，类似于例如椭 圆形或卵形区域，可以根据边界的曲率进行许多的度量且其沿着边界展开。图3 6示出图解多次波影响集合的动态孔径确定的本发明的实施例的流程图。 图3是示出本发明的一般实施例的流程图。图4是示出本发明的更详细实施例的流程图。 图5 6是进一步示出如图3 4所描述的本发明的要素的流程图。图3是示出用于使用多次波影响集合的动态孔径确定来将地震数据变换成衰减 过噪声的数据的本发明的实施例的流程图。在方框30处，根据表示在多次波影响集合的孔 径的边界附近的地下地层的地震数据，计算一系列的相邻多次波影响轨迹组的倾角。通过 在要评估的地球的地下区域附近部署多个地震传感器来获取地震能量，所述地震传感器响 应于地震能量而生成电学或光学信号中的至少一个。在方框31处，基于在方框30中计算 的倾角来递归地扩展多次波影响集合的孔径的边界。在方框32处，将具有来自方框31的 扩展孔径边界的多次波影响集合中的多次波影响轨迹堆叠以生成预测的多次波轨迹。方框 30 32预测地震数据中的多次波轨迹。在方框33处，从地震数据减去来自方框32的预测 的多次波轨迹以生成对形成地下地层的图像有用的表面相关多次波衰减数据。图4是示出多次波影响集合的动态孔径确定的本发明的另一实施例的流程图。图 4所示的实施例包括比图3所示的实施例更详细的说明。在方框40处，构造多次波影响集合。在下文的讨论中相对于图5中的流程图更详 细地描述了多次波影响集合的构造。在方框41处，选择用于在方框40中确定的多次波影响集合的初始孔径。可以以 适合于来自方框40的多次波影响集合的任何形状和尺寸来选择初始孔径。形状一般包括 但不限于矩形、椭圆形和不规则形状。在一个实施例中，初始孔径的选择可以将来自先前确 定的多次波影响集合的最终孔径考虑在内。在另一实施例中，可以做出对用于轨迹的初始孔径的选择以作为轨迹相关属性的 函数而变，所述轨迹相关属性包括但不限于偏移，意指轨迹的源和接收器坐标之间的距离； 方位角，意指通过源和接收器坐标的线与预定义参考方向之间的角度；以及纵向和横向偏 移，意指分别沿着参考方向和垂直于参考方向的偏移距离。在方框42处，为在方框41中选择的初始孔径选择最大孔径。可以基于为其进行 多次波预测的源-接收器轨迹的最大速度、长度、和该源-接收器轨迹的偏移(源与接收器 之间的距离)来定义最大孔径。在方框43处，根据在方框41中所选的初始孔径和在方框42中所选的最大孔径来 选择最大孔径扩展。最大孔径扩展可以沿边界而变化，特别是当多次波影响集合的孔径的 初始边界具有不规则形状或由于先前仅沿着边界的一部分应用的扩展而变为不规则形状 时。可以基于包括但不限于被度量轨迹中的时间位置、倾角的振幅和角度、以及来自 边界附近(并包括边界)的多次波影响集合的振幅、能量或其它导出值的因素来构建度量。 可以使用此度量来将递归约束至孔径的可行更新。在另一实施例中，还可以在选择用于未 来多次波影响集合的初始孔径时将此度量考虑在内。在方框44处，递归地扩展孔径。在下文的讨论中将相对于图6中的流程图来更详细地描述此递归孔径扩展。在方框45处，确定当前孔径是否足够。如果当前孔径不足，则过程继续至方框46。 如果当前孔径足够，则过程继续至方框47。足够的孔径意指应由所估计的多次波轨迹来适当地估计存在于多次波影响集合 的源-接收器轨迹中的所有多次波反射。因此，被在方框40中的多次波影响集合的构造中 使用的孔径所覆盖的表面反射点应对应于所有这些多次波反射。每个表面反射点具有多次 波影响集合中的相应顶点。顶点是对应于空间变化的表面的最小时间的空间位置 ，由此当 远离顶点移动时，表面的时间增加。在替换实施例中，可以用维度深度来代替维度时间。另外，当Fresnel堆叠对于那些顶点适当地作用时，孔径是足够的。Fresnel堆叠 是数据的堆叠以获得特定反射事件的更高的信噪比，由此，要堆叠的数据量与相应反射物 的一部分相关，来自该部分的反射能可以在第一反射能的二分之一波长内到达检测器。在方框46处，根据在方框44中应用的空间扩展的量来减小在方框43中所选的最 大孔径扩展。在方框47处，确定更多的多次波影响集合是否仍具有所确定的孔径。如果仍存在 更多的多次波影响集合，则过程返回到方框40以选择另一多次波影响集合。如果没有更多 的多次波影响集合需要检验，则过程结束。图5是示出关于多次波影响集合的构造的本发明的要素的流程图。图5更详细地 示出上文在图4的方框40中所讨论的本发明的部分。在方框50处，选择源位置。在方框51处，选择接收器位置。来自方框50的源位置和来自方框51的接收器位 置定义源-接收器轨迹。在方框52处，选择联接位置。联接位置对应于与源_接收器轨迹相关的多次波的 潜在表面反射点。联接位置取决于如参照图4中的流程图的讨论所构造的孔径和联接位置 之间的适当间距。在方框53处，从具有与在方框50中所选的源位置基本一致的源位置和与在方框 52中所选的联接位置基本一致的接收器位置的可用轨迹中选择第一轨迹。因此，第一轨迹 是用于在方框50中所选的源位置的公共源集合的成员。在方框54处，从具有与在方框52中所选的联接位置基本一致的源位置和与在方 框51中所选的接收器位置基本一致的接收器位置的可用轨迹中选择第二轨迹。因此，第二 轨迹是用于在方框51中所选的接收器位置的公共接收器集合的成员。公共源集合中的第一轨迹的接收器坐标与公共接收器集合中的第二轨迹的源坐 标基本上一致。由于记录数据中的源和接收器位置很少是精确一致的，所以在另一实施例 中可能需要使数据组规则化为标称几何结构以便实现源和接收器位置的此一致性。在另一 实施例中，随后在使用规则化数据组的同时针对原始几何结构来预测的多次波。在方框55处，相互地对来自方框53的第一轨迹和来自方框54的第二轨迹求卷 积。可以在时域进行此卷积，或者，可以作为频域中的乘法来执行卷积。从卷积得到的轨迹 是多次波影响轨迹。在这里生成的所有多次波影响轨迹组形成多次波影响集合。在方框56处，确定是否仍要对任何其它成对的第一和第二轨迹求卷积。如果仍有 更多成对的第一和第二轨迹，则过程返回到方框52以选择另一对。如果不存在其它成对的第一和第二轨迹，则过程结束。图6是示出关于递归孔径扩展的本发明的要素的流程图。图6更详细地示出上述 在图4的方框44中所讨论的本发明的部分。在方框60处，在多次 波影响集合的孔径的边界附近选择来自多次波影响集合的 相邻多次波影响轨迹组。如将在以下方框61中进行的那样，由需要多少轨迹以计算用于倾 角的可靠估计来确定在相邻多次波影响轨迹组中选择多少轨迹。在一个示例性实施例中，相邻多次波影响轨迹组包括两个相邻的多次波影响轨 迹。在另一示例性实施例中，相邻多次波影响轨迹组包括多于两个的相邻多次波影响轨迹 序列。在另一示例性实施例中，相邻多次波影响轨迹组包括相邻多次波影响轨迹块。例如， 该块可以在纵向5X横向5的轨迹块中包括25个相邻多次波影响轨迹。本发明的方法并 不意图局限于相邻多次波影响轨迹组的这些示例，而是可适用于任何此类相邻多次波影响 轨迹组。在方块61处，计算在方框60中所选的相邻多次波影响轨迹组的倾角。可以使用 集合中的倾角来获悉顶点在集合中的什么位置。倾角是同一地震事件的集合中的两个相邻 多次波影响轨迹之间的斜率。作为时间差与空间差之间的比(类似于dt/dx)来测量该倾 角。虽然倾斜事件通常意指非水平的，但不排除水平事件。可以通过检验集合的边界处的倾角来检测过小的多次波影响集合的孔径并自动 地将其放大。可以递归地进行检测和放大，由此，所使用的最后一个孔径对用于确定新孔径 的检测和放大的过程进行初始化。如果在多次波影响集合的边界的某一部分附近存在负的 倾角，则应扩展该部分边界，因为在该集合外面存在顶点。在方框62处，确定在方框61中所计算的倾角是否是负的。负的倾角意味着从多 次波影响集合的外部到内部的事件的斜率在时间和深度方面也增大。这意味着该事件的顶 点在集合的外面，因为顶点将始终位于较早的时间。如果倾角是负的，则过程继续至方框63 以扩展孔径的边界。如果倾角不是负的，则过程继续至方框65。在方框63，确定边界扩展的量。可以使得孔径的边界扩展量根据以下因素而变 该因素包括但不限于时间位置、与源-接收器轨迹相关的水深、和在边界附近检测到的最 大倾角。在方框64处，将在方框60中所选的相邻多次波影响轨迹组附近将孔径扩展在方 框63中确定的边界扩展量。在方框65处，确定是否仍有任何更多的相邻多次波影响轨迹组。如果仍有更多的 相邻多次波影响轨迹组，则过程返回到方框60以选择孔径边界附近的另一相邻多次波影 响轨迹组。如果不再有相邻多次波影响轨迹组，则过程结束。在替换实施例中，将在边界的扩展之后会被添加到多次波影响集合的多次波影响 轨迹堆叠以生成修正轨迹。然后可以将此修正轨迹应用于通过累计(通过Fresnel堆叠) 边界扩展之前的多次波影响集合中的多次波影响轨迹所确定的预测的多次波轨迹。在另一实施例中，可以用层间多次波消除(Interbed MultipleElimination, IME) 来代替消多次波(demultiple)方法SRME，层间多次波消除指的是预测层间或内部多次波 的扩展，由此使用两个卷积和轨迹的三元组(triplet)作为该方法的一部分。上文已仅仅出于说明的目的将本发明作为一种方法进行讨论，但还可以将其实现为系统。优选地借助于计算机、特别是数字计算机、以及其它传统数据处理设备来实现本发 明的系统。在本领域中众所周知的此类数据处理设备将包括计算机处理设备的任何适当组 合或网络，包括但不限于硬件(处理器、临时和永久性存储设备、以及任何其它适当的计算 机处理设备)、软件(操作系统、应用程序、数学程序库、以及任何其它适当软件)、连接(电 学、光学、无线、或其它方式)、以及外围设备(诸如键盘、定位设备、以及扫描仪的输入和输 出设备；诸如监视器和打印机的显示设备；诸如磁带、磁盘、和硬盘驱动器的计算机可读存 储介质、以及任何其它适当设备)。在另一实施例中，可以将本发明实现为上述方法，尤其是通过使用可编程计算机 来执行该方法而执行的。在另一实施例中，可以将本发明实现为存储在计算机可读介质中 的计算机程序，该程序具有可用于促使可编程计算机执行上述方法的逻辑。在另一实施例 中，可以将本发明实现为具有存储在介质上的计算机程序的计算机可读介质，使得该程序 具有可用于促使可编程计算机执行上述方法的逻辑。应理解的是前述内容仅仅是本发明的特定实施例的详细说明且在不脱离本发明 的范围的情况下可以依照这里的公开内容对所公开的实施例进行许多变更、修改、和替换。 因此，前述说明并不意图限制本发明 的范围。相反，仅仅由随附权利要求及其等价物来确定 本发明的范围。
权利要求
一种用于将地震数据变换成衰减过噪声的数据的方法，包括通过以下各项来预测表示地下地层的地震数据中的多次波轨迹根据通过在要评估的地球地下区域附近部署多个地震传感器获取的地震数据，计算一系列的相邻多次波影响轨迹组的倾角，所述地震传感器响应于在多次波影响集合的孔径边界附近的地震能量而生成电学和光学信号中的至少一个；基于所计算的倾角，递归地扩展多次波影响集合的孔径的边界；以及将具有扩展孔径边界的多次波影响集合中的多次波影响轨迹堆叠以生成预测的多次波轨迹；以及从地震数据减去所预测的多次波轨迹以生成对形成地下地层的图像有用的表面相关的多次波衰减数据。
2.权利要求1的方法，还包括 构造多次波影响集合；选择多次波影响集合的初始孔径；选择初始孔径的最大孔径；选择初始孔径和最大孔径的最大孔径扩展；递归地扩展孔径；以及基于孔径的扩展来减小最大孔径扩展。
3.权利要求2的方法，其中，构造多次波影响集合包括 选择源位置；选择接收器位置； 选择联接位置组；对联接位置组中的每个联接位置重复以下各项选择具有与所选源位置基本一致的源位置和与联接位置基本一致的接收器位置的第 一轨迹；选择具有与联接位置基本一致的源位置和与所选接收器位置基本一致的接收器位置 的第二轨迹；以及对第一轨迹与第二轨迹求卷积，生成多次波影响轨迹；以及 由多次波影响轨迹生成多次波影响集合。
4.权利要求2的方法，其中，递归地扩展孔径包括从多次波影响集合的孔径边界附近的多次波影响集合中选择一系列的相邻多次波影 响轨迹组；以及对所述一系列的相邻多次波影响轨迹组中的每个相邻多次波影响轨迹组重复以下各项计算所述相邻多次波影响轨迹组的倾角；以及 对于负的倾角重复以下各项 确定孔径边界扩展量；以及将所述相邻多次波影响轨迹组附近的边界扩展所述边界扩展量。
5.权利要求4的方法，还包括将在边界扩展之后添加到多次波影响集合的多次波影响轨迹堆叠，生成修正轨迹；以及应用所述修正轨迹来修正在边界扩展之前对多次波影响集合确定的预测的多次波轨迹。
6.一种多次波影响集合的动态孔径确定的方法，包括 使用可编程计算机来执行以下各项计算多次波影响集合的孔径边界附近的一系列的相邻多次波影响轨迹组的倾角；以及 基于所计算的倾角，递归地扩展多次波影响集合的孔径边界。
7.权利要求6的方法，还包括 使用可编程计算机来执行以下各项 构造多次波影响集合；选择多次波影响集合的初始孔径；选择初始孔径的最大孔径；选择初始孔径和最大孔径的最大孔径扩展；递归地扩展孔径；以及基于孔径的扩展来减小最大孔径扩展。
8.权利要求7的方法，其中，构造多次波影响集合包括 使用可编程计算机来执行以下各项选择源位置； 选择接收器位置； 选择联接位置组；对联接位置组中的每个联接位置重复以下各项选择具有与所选源位置基本一致的源位置和与联接位置基本一致的接收器位置的第 一轨迹；选择具有与联接位置基本一致的源位置和与所选接收器位置基本一致的接收器位置 的第二轨迹；以及对第一轨迹与第二轨迹求卷积，生成多次波影响轨迹；以及 由多次波影响轨迹生成多次波影响集合。
9.权利要求7的方法，其中，递归地扩展孔径包括 使用可编程计算机来执行以下各项从多次波影响集合的孔径边界附近的多次波影响集合中选择一系列的相邻多次波影 响轨迹组；以及对所述一系列的相邻多次波影响轨迹组中的每个相邻多次波影响轨迹组重复以下各项计算所述相邻多次波影响轨迹组的倾角；以及 对于负的倾角重复以下各项 确定孔径边界扩展量；以及将所述相邻多次波影响轨迹组附近的边界扩展所述边界扩展量。
10.权利要求9的方法，还包括 使用可编程计算机来执行以下各项将在边界扩展之后添加到多次波影响集合的多次波影响轨迹堆叠，生成修正轨迹；以及应用所述修正轨迹来修正在边界扩展之前对多次波影响集合确定的预测的多次波轨迹。
11.一种具有存储在上面的计算机程序的计算机可读介质，所述程序具有可用于促使 可编程计算机执行包括以下的步骤的逻辑计算多次波影响集合的孔径边界附近的一系列的相邻多次波影响轨迹组的倾角；以及 基于所计算的倾角，递归地扩展多次波影响集合的孔径边界。
12.权利要求11的介质，还包括 构造多次波影响集合；选择多次波影响集合的初始孔径；选择初始孔径的最大孔径；选择初始孔径和最大孔径的最大孔径扩展；递归地扩展孔径；以及基于孔径的扩展来减小最大孔径扩展。
13.权利要求12的介质，其中，构造多次波影响集合包括 选择源位置；选择接收器位置； 选择联接位置组；对联接位置组中的每个联接位置重复以下各项选择具有与所选源位置基本一致的源位置和与联接位置基本一致的接收器位置的第 一轨迹；选择具有与联接位置基本一致的源位置和与所选接收器位置基本一致的接收器位置 的第二轨迹；以及对第一轨迹与第二轨迹求卷积，生成多次波影响轨迹；以及 由多次波影响轨迹生成多次波影响集合。
14.权利要求12的介质，其中，递归地扩展孔径包括从多次波影响集合的孔径边界附近的多次波影响集合中选择一系列的相邻多次波影 响轨迹组；以及对所述一系列的相邻多次波影响轨迹组中的每个相邻多次波影响轨迹组重复以下各项计算所述相邻多次波影响轨迹组的倾角；以及 对于负的倾角重复以下各项 确定孔径边界扩展量；以及将所述相邻多次波影响轨迹组附近的边界扩展所述边界扩展量。
15.权利要求14的介质，还包括将在边界扩展之后添加到多次波影响集合的多次波影响轨迹堆叠，生成修正轨迹；以及应用所述修正轨迹来修正在边界扩展之前对多次波影响集合确定的预测的多次波轨迹。
全文摘要
本发明涉及用于三维表面相关的多次波消除的动态孔径确定的方。根据表示在多次波影响集合的孔径边界附近的地下地层的地震数据，计算一系列的相邻多次波影响轨迹组的倾角，所述地震数据是通过在要评估的地球地下区域附近部署多个地震传感器获取的，所述地震传感器响应于地震能量生成电学和光学信号中的至少一个。基于所计算的倾角，递归地扩展多次波影响集合的孔径的边界。将具有扩展孔径边界的多次波影响集合中的多次波影响轨迹堆叠以生成预测的多次波轨迹。从地震数据减去预测的多次波轨迹以生成可用于形成地下地层的图像的表面相关多次波衰减数据。
文档编号G01V1/32GK101937099SQ201010222060
公开日2011年1月5日 申请日期2010年6月30日 优先权日2009年6月30日
发明者R·F·赫格, R·G·范博尔塞伦 申请人:Pgs地球物理公司