分离由干扰震源产生的地震信号的制作方法

专利名称：分离由干扰震源产生的地震信号的制作方法
技术领域：
本发明总体上涉及分离由干扰震源产生的地震信号。
背景技术：
地震勘探包括对地下地质结构进行碳水沉积物勘测。勘测通常包括在预定位置处 部署震源和地震传感器。震源产生地震波，地震波沿着它们的路径传播到产生压力变化和 振动的地质结构中。地质结构的弹性特性的变化散射地震波，改变其传播方向及其它性质。 由震源发射的那部分能量到达地震传感器。某些地震传感器对压力变化灵敏(水听器)，其 它的对质点运动灵敏(例如，地震检波器)，并且工业勘测可以仅部署一种类型的传感器或 两种。响应于所检测的地震事件，传感器生成电信号以产生地震数据。然后，地震数据的分 析可以指示碳水沉积物的可能位置的存在或不存在。某些勘测被称为“海洋”勘测，因为其在海洋环境中进行。然而，“海洋”勘测不仅 可以在咸水环境中进行，而且可以在淡水和微咸水中进行。在称为“牵引阵列”勘测的一类 海洋勘测中，在勘测船后面牵引地震传感器阵列-包含浮缆(streamer)和震源。

发明内容
在本发明的实施例中，一种技术包括获得地震数据，该地震数据指示多个震源的 引爆(fire)产生的复合地震信号的由地震传感器采集的测量结果。该技术包括使描述与 复合地震信号相关的地质情况的模型与线性算子(operator)相关联并将地震数据表征为 该模型和相关联的线性算子的函数。该技术包括基于函数来同时地确定模型并基于所确定 的模式来生成数据集。每个数据集指示复合地震信号的分量并可归因于震源中的不同的一 个。在本发明的另一实施例中，一种系统包括接口和处理器。接口接收地震数据，该地 震数据指示多个震源的引爆所产生的复合地震信号的由地震传感器采集的测量结果。处理 器处理地震数据以使线性算子与描述与复合地震信号相关的地质情况的模型相关联；地震 数据被表征为该模型和相关线性算子的函数；基于该函数来同时地确定模型；并基于所确 定的模型来生成数据集。每个数据集指示复合地震信号的分量并可归因于震源中的不同的 一个。在本发明的另一实施例中，一种产品包括包含指令的计算机可访问存储介质，该 指令在被基于处理器的系统执行时引起基于处理器的系统接收地震数据，该地震数据指示 多个震源的引爆所产生的复合地震信号的由地震传感器采集的测量结果。该指令在被执行 时引起基于处理器的系统处理地震数据以使线性算子与描述与复合地震信号相关的地质 情况的模型相关联；地震数据被表征为该模型和相关线性算子的函数；基于该函数来同时 地确定所述模型；并基于所确定的模型来生成数据集。每个数据集指示复合地震信号的分 量并可归因于震源中的不同的一个。通过以下附图、说明和权利要求书，本发明的优点及其它特征将变得明显。


图1是根据本发明的实施例的基于海洋的地震采集系统的示意图。图2、3和11是描绘根据本发明的实施例的分离由干扰震源产生的地震信号的技 术的流程图。图4、5、6、7、8、9和10是模拟震源和接收机信号，其图解根据本发明的实施例的复 合地震信号到用原始震源可识别的信号的分离。图12是根据本发明的实施例的数据处理系统的示意图。
具体实施例方式图1描绘依照本发明的某些实施例的基于海洋的地震数据采集系统的实施例10。 在系统10中，勘测船20在船20后面牵引一个或多个地震浮缆30 (在图1中描绘一个示例 性浮缆30)。应注意的是可以在其中在同一深度处的大约相同的平面中牵引多个浮缆30的 散布范围内布置浮缆30。作为另一非限制性示例，例如，可以在多个深度处(诸如在范围散 布之上/之下内)牵引浮缆。地震浮缆30可以有几千米长且可以包含各种支撑电缆(未示出)，以及可以用来 支持浮缆30之间的通信的布线和/或电路(未示出)。通常，每个浮缆30包括向其中安装 了记录地震信号的地震传感器的主电缆。浮缆30包含地震传感器58，根据本发明的特定实 施例，其可以是采集压力数据的水听器(水听器是非限制性示例)或多分量传感器。对于 其中传感器58是多分量传感器(作为另一非限制性示例)的本发明的实施例，每个传感器 能够检测与接近于传感器的声学信号相关的质点运动的至少一个分量和压力波场。质点运 动的示例包括质点位移的一个或多个分量、质点速度的一个或多个分量(线内方向(X)、交 叉线方向(y)和竖直方向(ζ)分量)(例如，参见轴59)和质点加速度的一个或多个分量。根据本发明的特定实施例，所述多分量地震传感器可以包括一个或多个水听器、 地震检波器、质点位移传感器、质点速度传感器、加速计、压力梯度传感器、或其组合。例如，依照本发明的某些实施例，特定的多分量地震传感器可以包括用于测量压 力的水听器和测量传感器附近的质点速度和/或加速度的三个相应正交分量的三个正交 对准的加速计。应注意的是可以将多分量地震传感器实现为单个设备(如图1所描绘的) 或者可以实现为多个设备，取决于本发明的特定实施例。特定多分量地震传感器还可以包 括压力梯度传感器，其构成另一种质点运动传感器。每个压力梯度传感器测量特定点处的 压力波场相对于特定方向的变化。例如，压力梯度传感器之一可以采集地震数据，该地震数 据指示特定点处压力波场相对于交叉线方向的偏导数，并且另一个压力梯度传感器可以采 集特定点处的指示相对于线内方向的压力数据的地震数据。海洋地震数据采集系统10包括一个或多个震源40 (在图1中描绘的两个示例性 震源40)，诸如气枪等。在本发明的某些实施例中，可以将震源40耦合到勘测船20或由勘 测船20牵引。或者，在本发明的其它实施例中，仅仅作为几个示例，震源40可以独立于勘 测船20进行操作，因为震源40可以耦合到其它船或浮体。由于在勘测船20后面牵引地震浮缆30，所以由震源40产生常常称为“射击”的声 学信号42 (在图1中描绘的示例性声学信号42)并被向下指引通过水柱44到达水底面24
5下面的地层62和68。声学信号42从各种地下地质结构反射，诸如图1所描绘的示例性结 构65。被震源40采集的入射声学信号42产生被地震传感器58感测的相应反射声学信 号或压力波60。应注意的是由地震传感器58接收和感测的压力波包括在不反射的情况下 传播到传感器58的“上行”压力波以及由压力波60从空气_水界面31的反射产生的“下 行”压力波。地震传感器58产生称为“轨迹”的信号(例如，数字信号)，其指示所采集的压力 波场和质点运动的测量结果。依照本发明的某些实施例，该轨迹被记录且可以至少部分地 由部署在勘测船20上的信号处理单元23进行处理。例如，特定地震传感器58可以提供轨 迹，其对应于由其水听器55进行的压力波场的测量；并且传感器58可以(根据本发明的特 定实施例)提供对应于质点运动的一个或多个分量的一个或多个轨迹。地震采集的目的是出于识别诸如示例性地质结构65的地下地质结构的目的构建 勘测区域的图像。该表示的后续分析可以揭示地下地质结构中的碳水沉积物的可能位置。 根据本发明的特定实施例，可以诸如由信号处理单元23对地震勘测船20执行该表示的部 分分析。依照本发明的其它实施例，可以由可以例如位于陆地上或船20上的地震数据处理 系统(诸如在图12中描绘且下文进一步描述的示例性地震数据处理系统320)来处理该表 示。因此，可以进行许多修改且其在随附权利要求书的范围内。特定震源40可以由可以布置成阵列串(例如，枪串)的震源元件(例如，诸如气 枪)的阵列形成。或者，特定震源40可以由一个或预定数目的气枪的阵列形成，可以由多 个阵列形成，等等。无论震源的特定组成如何，可以在勘测期间按照特定时序对震源进行引
ο如下文更详细地描述的，可以依次对震源40进行引爆，使得可以同时或以短的时 间间隔接近同时地引爆多个震源40，以便由地震传感器58感测的复合能量信号包含来自 多于一个的震源40的大量能量。换言之，震源相互干扰，使得复合能量信号不能容易地分 离成归因于特定震源的信号。如下所述，由地震传感器58采集的数据被分离成每个与震源 40之一相关联的数据集，以便每个数据集指示可归因于相关震源40的复合地震能量信号 的分量。在传统牵引海洋勘测中，在一个震源的引爆与下一个震源的引爆之间引入延迟， 并且该延迟足以允许由一个震源的引爆产生的能量在与下一次震源引爆相关的能量到达 之前衰减至可接受水平。然而，此类延迟的使用对可以采集地震数据的速率施加约束。对 于牵引海洋勘测，这些延迟还暗示最小的线内射击间隔，因为勘测船的最小速度受到限制。因此，将其中来自震源的信号干扰的同时引爆或接近同时引爆震源用于每个记录 的至少一部分在采集效率和线内震源采样方面具有优点。然而，对于这种有用的技术而言， 必须将所采集的地震数据分离成每个唯一地与震源之一相关联的数据集。用于使得能够进行干扰震源的分离的一种传统技术利用震源的引爆之间的相对 小的延迟(例如，随机延迟)(即涉及震源抖动的使用)。结果得到的地震轨迹被收集到包 括每个源的许多引爆的域中。轨迹被对准使得时间零点对应于特定震源的引爆时间，以便 由于特定震源而采集的信号看起来相干，而由于其它震源采集的信号看起来非相干。基于 相干性来分离所采集的信号。
已经注意到表面上非相干的信号在数学上可能不是非相干的，因为使该信号看起 来非相干的震源引爆之间的时间延迟是已知的。因此，依照本文所述的本发明的实施例，将 由于干扰震源引爆而采集的所有能量视为单个复合能量信号；并且出于将复合能量信号分 解成每个唯一地与特定震源相关的信号的目的使用线性算子变换。更具体而言，图2描绘了通常可以用于分离由于干扰震源的引爆而采集的地震传 感器数据的目的的技术110。参照图2，技术110包括获得地震数据(称为“地震数据矢量 d”)，其是由于N个(即多个)震源的引爆而被地震传感器采集的。因此，同时地或以接近 同时的方式引爆震源，使得在地震数据矢量d中存在来自所有这些引爆的大量能量。根据 方框118，使描述影响震源能量的地质情况的模型与线性算子相关联，线性算子描述地震机 制的物理性质、波传播和勘测几何结构。地震数据矢量d被表征为模型和线性算子的函数 (方框122)。然后，对于模型对此函数进行共同反演，这允许将地震数据矢量d分离(方框 130)成N个地震数据集(V-dN，使得每个数据集唯一地可归因于震源中的一个。换言之，每 个数据集表示所感测的复合能量信号的分量，该分量唯一地可归因于震源中的一个。作为更具体的示例，假设由于称为“S/’和“S2”的两个震源的接近同时引爆而采集 地震数据矢量d。对于本示例，按照时序来引爆震源S1和s2，该时序可以基于预定定时模式 或者可以基于随机或伪随机时间。无论特定的定时方案如何，对于本示例假设对于所有轨 迹而言在震源S2之前引爆震源S1，并且还假设轨迹的零点时间对应于S1的引爆时间。因 此，轨迹的零点时间处于“Si时间”。分别将对震源S1和S2的偏移或矢量称为“X1”和“X2”。 对于每个轨迹而言，用于震源S2的用“t”表示的定时延迟是已知的。对于本示例而言假设轨迹的收集使得t值是随机的。实际上，这是用于CMP、接收 机或共偏移集合(gather)的情况。出于简化本讨论的目的，假设可以分别使用称为“χ1/’ 和“χ2/’的标量相对于震源S1和震源S2对每个集合中的轨迹进行定位。在这种符号表示法 中，下标“i”表示集合中的轨迹数目。作为更具体的示例，对于CMP集合而言，“χ1/’可以是 到震源Sl的标量偏移，并且在下文中将这些量称为偏移。同样地，“t/’表示用于第i轨迹 的定时延迟。震源S1的记录能量可以通过对称为“m/’的未知模型应用称为“L/’的线性算子 (其表示震源S1的物理性质、与震源S1相关的波传播和与震源S1相关的勘测几何结构)来 建模，所述未知模型描述影响从震源S1传播的能量的地质情况。模型HI1包含用于模型空间 中的每个参数的一个元素。通常，分别通过对应于线性或双曲线/抛物线拉冬(Radon)变 换的慢度(slowness)或其平方来将模型空间参数化。线性算子L1是到震源S1的偏移、表 征模型空间的参数、和时间或频率的函数。地震数据矢量Cl1包含用于每个轨迹的一个元素 (在每个时间或频率)且是地震数据d的分量，其与震源S1相关。换言之，地震数据矢量Cl1 表示可归因于震源S1的数据集。地震数据矢量Cl1可以如下描述Cl1 = L1Hi1等式 1在地震数据矢量d中，与震源S2相关的能量看起来是非相干的。然而，该能量通 过对轨迹应用时移ti与其中震源S2的引爆时间在时间零点(即，震源S2时间)的相干数 据集相关。可以出于描述这些时移的目的使用称为“D2”的对角线性算子，使得可以如下描 述与震源S2相关并称为“d2”的地震数据矢量d的分量d2 = D2L2Hi2等式 2
在等式2中，称为“L2”的线性算子表示震源S2的物理性质、与震源S2相关的波传 播和与震源S2相关的勘测几何结构。同样在等式2中，称为“m2”的模型描述影响从震源S2 传播的能量的地质情况。由地震传感器记录的复合地震能量信号可归因于震源S1和S2。因此，如下所述， 地震数据矢量d(即，记录的数据)是地震数据矢量dl和d2的组合d = Cl^d2等式 3由于等式1、2和3中的关系，可以将地震数据矢量d表示为以下线性系统

d = [Lx D2L2,
等式4 因此，可以使用诸如最小二乘算法的标准技术针对模型矢量m(即Oii1 ;m2))对等式 4进行求解(即，共同反演)；并且在知道模型矢量m之后，可以出于将地震数据矢量d分离 成地震数据矢量Cl1和d2、即分离成指示可归因于每个震源的测量结果的数据集的目的对模 型叫和m2应用等式1和2。因此，参照图3，依照本发明的某些实施例，可以将技术150用于分离由干扰震源 (对于本示例而言为两个震源)产生的地震数据。按照技术150，获得地震数据矢量d，按 照方框154，其是由于震源的接近同时引爆而采集的。按照方框158，使模型Hi1和m2与描述 震源机制的物理性质、波传播和勘测结构(1^和1^)及震源引爆之间的定时(D2)的线性算 子LpL2和D2相关联。然后，地震数据矢量d被表征为模型Hi1和m2及线性算子LpL2和D2 的函数(方框162)。然后按照方框166，对于模型Hi1和m2而言，对该函数进行共同反演 ’然 后，按照方框170，可以将地震数据矢量d分离成地震数据矢量Cl1和d2。可以在频(ω)域中对等式4进行反演。在这种情况下，(D2)jk = exp(-i ω、)δ Jk 且(Ls) Jk = exp (-i ω tsJk)，其中，tsJk是与偏移Xsj相关的时移和与Ss相关的模型空间中的 第k个轨迹的参数。对于用慢度参数化的线性拉冬变换而言，psk，tsJk = XsjPsko对于用曲 率参数化的抛物线拉冬变换而言，q\，tsJk = (Xsj)2QskO上文所述的震源分离技术的成功取决于所述变换分离与两个震源相关的能量的 能力。不同于拉冬变换的大部分应用，成功不取决于将能量集中于Hl1或!112内的正确模型 参数处的能力。当在震源引爆之间使用随机或伪随机时间延迟时，用于两个模型域(〖\和 tj+t2jk)的基本函数相差悬殊，并且这使得能够极其有效地进行震源分离。模型域的参数化的细节不是重要的，只要有可能使用该域对记录的数据进行建 模。例如，对于线性拉冬变换而言，慢度范围必须覆盖在数据中观察到的范围，并且采样必 须足以避免混淆(aliasing)。通常不期望需要使用高分辨率变换来改善集中。然而，如果 需要，则可以使用高分辨率变换，例如由于由偏移窗口化或采集几何结构问题产生的偏移 采样不良。图5、6、7、8、9和10描绘应用于简单、合成数据集时的技术150的示例。通过分别 添加对应于震源Sl和S2的合成信号206 (参见图5)和210来形成到分离过程的输入信号 200 (参见图4)(即由地震传感器记录的模拟信号)。输入信号200还包含随机噪声，并且 信号200处于S1时间。信号206包含具有随机零偏移时间、振幅和速度及30Hz里克子波 (Ricker wavelet)的10个双曲线事件。对于S2时间而言，输入信号200对应于图7中的 输入信号214。如从图7可以看到的那样，时间延迟的去除使得S2相关信号214相干。
8
分离过程是针对从采集的输入信号200 (图4)恢复Sl输入信号206 (图5)和S2 输入信号210 (图5)。分别在图8 (分离的S1信号218)和9 (分离的S2信号222)中描绘 了结果得到的估计。输入信号200中的几乎所有能量都出现在信号218或信号222中。如 图10的信号224所示，可以通过到S2时间的时移来使得S2相关数据相干。然后，可以分别 使用到S1和S2的偏移，在传统地震数据处理流程中处理输出数据(即，信号218和224)。虽然上文所述的示例使用震源抖动、或震源的非同时引爆，但依照本发明的其它 实施例，可以同时对震源进行引爆。关于此方面，如果使线性算子成为地震数据的更加唯一 的预报器，则对震源引爆的抖动的要求变得次要。换言之，如果存在用于震源位置的基本函 数的较少重叠，则震源抖动可能是次要的。作为更具体的示例，可以出于使线性算子成为地震数据的更唯一预报器的目的将 本文所述的技术与用于震源分离的其它技术组合。例如，可以确定性地估计波场的某些部 分(例如，诸如直达波)并将其减去作为预处理步骤。另外，可以与本文所述的技术组合地 使用诸如倾角滤波(dip-filtering)的方法。作为更具体的示例，可以将从震源S1记录的能量视为由直达波产生的能量和由反 射产生的能量的组合。由此，可以有效地如下表示地震数据矢量Cl1 Cl1 = dn+cV = L1Hi^H1Iiih,等式 5其中，“屯”表示可归因于来自震源S1的直达波(direct arrival)的地震数据； “dlh”表示可归因于由于震源S1而产生的反射的地震数据；“L/’表示与来自震源S1的直达 波相关的线性拉冬算子；"Hi1 ”表示描述影响直达波的地质情况的模型；“H/’表示与由于来 自震源S1的能量而产生的反射相关的双曲线拉冬变换算子；且“mh”表示描述影响由震源产 生的反射的地质情况的模型。类似地，可以如下描述地震数据矢量，其为可归因于从震源S2记录的能量的d2 d2 = d21+d2h = L2Hi^H2Hih,等式 6其中，“d21 ”表示可归因于直达波的地震数据矢量d2的分量；“d2h”表示可归因于 反射的地震数据d2 ；"L2”表示与来自震源S2的直达波相关的线性拉冬变换算子；且"H2”表 示与由于来自震源S2的能量而产生的反射相关的双曲线拉冬变换。由于在等式5和6中阐述的关系，可以如下表示地震数据矢量d，其表示由地震传 感器记录的实际数据d = dn+cU+c^+cU，等式 7因此，可以用以下函数来表示地震数据矢量d，对于模型Hi1和％而言可以对其进 行反演^ = KL1+ L2) (H1^H2)] m',等式 8然后可以应用等式5和6来导出数据矢量Cl1和d2。虽然上文已描述了线性和双曲线拉冬变换，但应注意的是依照本发明的其它实施 例，可以使用其它线性算子。例如，仅仅作为几个其它非限制性示例，依照本发明的其它实 施例，可以使用抛物线算子或迁移算子(migration operator) 0因此，参照图11，处于分离由于由干扰震源产生的能量而采集的地震数据的目的， 依照本发明的某些实施例，可以使用技术200，对于本示例而言所述干扰震源是两个震源S1
9和S2。按照技术200，获得地震数据矢量d(方框204)，其是由于震源的引爆而采集的。使 描述与直达波(Hi1)和反射(mh)相关的地质情况的模型与线性算子L1和L2 (对于直达波) 及H1和H2 (对于反射)相关联(方框208)。按照方框212，地震数据矢量d被表征为模型 Hi1和mh及线性算子Lp L2, H1和H2的函数。然后，对于模型Hi1和mh而言，按照方框216，对 函数进行共同反演。随后，按照方框220，可以将地震数据矢量d分离成数据子集矢量Cl1和 d2。虽然本文所阐述的示例是用于两个震源S1和S2，但可以将该技术扩展至多于两个 的震源。参照图12，依照本发明的某些实施例，地震数据处理系统320可以出于分离由于 干扰震源产生的能量而采集的地震数据的目的执行本文所公开的至少某些技术。依照本发 明的某些实施例，系统320可以包括处理器350，诸如一个或多个微处理器和/或微控制器。 处理器350可以位于浮缆30(图1)上、位于船20上或位于陆地处理设施处(作为示例)， 这取决于本发明的特定实施例。可以出于从地震传感器58接收对应于压力和/或质点运动测量结果的地震数据 的目的将处理器350耦合到通信接口 360。因此，依照本文所述的本发明的实施例，处理器 350在执行存储在地震数据处理系统320的存储器中的指令时可以接收由地震传感器在牵 引中所采集的多分量数据和/或压力传感器数据。应注意的是根据本发明的特定实施例， 数据可以是在正在采集数据时直接从传感器接收的数据(对于其中处理器350是勘测系统 的一部分的情况，诸如是船或浮缆的一部分)，或者可以是先前由地震传感器在牵引中所采 集且被存储并传送到处理器350的传感器数据，例如，处理器350可以在陆地设施中。作为示例，接口 360可以是USB串行总线接口、网络接口、可移动介质(诸如闪存 卡、⑶-ROM等)接口或磁性存储接口(例如，IDE或SCSI接口)。因此，根据本发明的特定 实施例，接口 360可以采取许多形式。依照本发明的某些实施例，接口 360可以耦合到地震数据处理系统320的存储器 340且可以存储例如与结合技术110、150和/或200来处理地震数据相关的各种输入和/ 或输出数据集，如附图标记348所指示的。依照本发明的某些实施例，存储器340可以存储 程序指令344，该程序指令344在被处理器350执行时可以引起处理器350执行本文所公 开的各种任务或更多技术，诸如技术110、150和/或200，并在系统的320的显示器(图12 未示出)上显示经由该技术获得的结果。其它实施例在随附权利要求书的范围内。例如，依照本发明的其它实施例，可以使 用“振幅抖动”来辅助分离。虽然被牵引震源的振幅的控制通常可能具有挑战性，但依照本 发明的实施例，可以通过根据随机或规则模式故意地不引爆选择的震源来控制震源。作为 另一示例，振幅抖动可以包括选在不引爆特定震源的其它元件的同时选择性地引爆该震源 的某些震源元件(例如，诸如枪)以改变振幅。可以将关于振幅抖动的信息并入上述线性算子。实际上，偶尔可能有震源之一引爆失败。当发生这种情况时，可以通过迫使对于相 应轨迹而言算子具有零输出来将关于失败震源的信息包括在相关线性算子中。这些不引爆 又可以使得不同的震源更容易分离。其它实施例在随附权利要求书的范围内。例如，虽然上文已描述了牵引式海洋地震采集系统，但本文所述的用于分离由干扰震源产生的地震信号的技术和系统可以同样地 应用于其它类型的地震采集系统。作为非限制性示例，可以将本文所述的技术和系统应用 于海床、钻孔和陆基地震采集系统。因此，根据本发明的特定实施例，地震传感器和震源可 以是固定的，或者可以被牵引。作为本发明的其它实施例的其它示例，地震传感器可以是采 集质点运动和压力的测量结果的多组分传感器，或者，地震传感器可以仅仅是采集压力测 量结果的水听器。因此，可以预期许多变更且其在随附权利要求书的范围内。
虽然已相对于有限的实施例描述的本发明，但受益于本公开的本领域的技术人员 应认识到其许多修改和变更。意图在于随附权利要求书涵盖在本发明的真实精神和范围内 的所有此类修改和变更。
权利要求
一种方法，包括获得地震数据，该地震数据指示多个震源的引爆所产生的复合地震信号的由地震传感器采集的测量结果；使模型与线性算子相关联，所述模型描述与所述复合地震信号相关联的地质情况；将地震数据表征为所述模型和相关联的线性算子的函数；基于所述函数同时地确定所述模型；以及基于所确定的模型，生成数据集，每个数据集可指示复合地震信号的分量且可归因于所述震源中的不同的一个。
2.根据权利要求1的方法，其中，同时地确定的动作包括针对模型对函数进行共同反演。
3.根据权利要求1的方法，其中，所述震源被同时引爆。
4.根据权利要求1的方法，其中所述震源包括第一震源和与第一震源在不同时间引爆的第二震源，以及相关联的动作包括使所述第二震源与描述第一和第二震源之间的引爆时间差的线性 算子相关联。
5.根据权利要求4的方法，其中，按照预定时间间隔的定时模式，在相对于第一震源的 时间引爆第二震源。
6.根据权利要求4的方法，其中，按照由随机数发生器控制的定时模式，在相对于第一 震源的时间引爆第二震源。
7.根据权利要求1的方法，其中，相关联的动作包括使描述地质情况的模型与由震源 的引爆产生的直达波和反射相关联。
8.根据权利要求1的方法，其中，所述线性算子包括线性拉冬算子和双曲线拉冬算子。
9.根据权利要求1的方法，其中，所述线性算子包括选自以下各项中的至少一个算子 线性拉冬算子、双曲线拉冬算子、抛物线算子和迁移算子。
10.根据权利要求1的方法，其中，震源的振幅以受控方式相对于彼此而变化。
11.根据权利要求10的方法，其中，震源的振幅根据随机或伪随机方式或按照振幅变 化的预定模式而变化。
12.—种系统，包括接口，其接收地震数据，该地震数据指示多个震源的引爆所产生的复合地震信号的由 地震传感器采集的测量结果；处理器，其处理地震数据以使线性算子与描述与复合地震信号相关的地质情况的模型 相关联，地震数据被表征为所述模型和相关联的线性算子的函数，基于所述函数同时地确 定所述模型并基于所确定的模型来生成数据集；其中，每个数据集指示复合地震信号的分量并可归因于所述震源中的不同的一个。
13.根据权利要求12的系统，其中，所述处理器适合于处理地震数据以针对所述模型 对所述函数进行共同反演。
14.根据权利要求12的系统，其中，所述震源被同时引爆。
15.根据权利要求12的系统，其中所述震源包括第一震源和与第一震源在不同时间引爆的第二震源，以及所述处理器被配置为使所述第二震源与描述第一和第二震源之间的引爆时间差的线 性算子相关联。
16.根据权利要求12的系统，其中，所述线性算子包括选自以下各项中的至少一个算 子线性拉冬算子、双曲线拉冬算子、抛物线算子和迁移算子。
17.根据权利要求12的系统，还包括包含地震传感器的至少一个被牵引浮缆，其中，所述处理器位于所述至少一个被牵引 浮缆上。
18.根据权利要求17的系统，还包括牵引所述至少一个被牵引浮缆的船。
19.根据权利要求10的系统，其中，震源的振幅以受控方式相对于彼此而改变。
20.根据权利要求19的系统，其中，震源的振幅根据随机或伪随机方式或按照振幅变 化的预定模式而变化。
21.一种包括计算机可访问存储介质的产品，所述计算机可访问存储介质包含的指令 在被基于处理器的系统执行时引起基于处理器的系统接收地震数据，该地震数据指示多个震源的引爆所产生的复合地震信号的由地震传感 器采集的测量结果；以及处理地震数据以使线性算子与描述与复合地震信号相关的地质情况的模型相关联，将 地震数据表征为所述模型和相关联的线性算子的函数，基于所述函数同时地确定所述模型 并基于所确定的模型来生成数据集；其中，每个数据集指示复合地震信号的分量并可归因于所述震源中的不同的一个。
22.根据权利要求21的产品，所述存储介质包含的指令在被基于处理器的系统执行时 引起基于处理器的系统处理地震数据以针对所述模型对所述函数进行共同反演。
23.根据权利要求21的产品，其中，所述震源被同时引爆。
24.根据权利要求21的产品，其中所述震源包括第一震源和与第一震源在不同时间引爆的第二震源，以及所述存储介质包含的指令在被基于处理器的系统执行时引起基于处理器的系统使所 述第二震源与描述第一和第二震源之间的引爆时间差的线性算子相关联。
25.根据权利要求21的产品，其中，所述线性算子包括选自以下各项中的至少一个算 子线性拉冬算子、双曲线拉冬算子、抛物线算子和迁移算子。
全文摘要
一种技术包括获得地震数据，该地震数据指示多个震源的引爆产生的复合地震信号的由地震传感器获得的测量结果(114)。该技术包括使描述与复合地震信号相关的地质情况的模型与线性算子相关联(118)并根据模型和相关线性算子来表征地震数据(122)。该技术包括基于函数来同时确定模型(126)并基于所确定的模型来生成数据集。每个数据集指示复合地震信号的分量并可归因于不同的震源之一(130)。
文档编号G01V1/36GK101925835SQ200880125359
公开日2010年12月22日 申请日期2008年11月22日 优先权日2007年12月26日
发明者C·考斯托夫, D·E·尼可斯, I·穆尔 申请人:格库技术有限公司