专利名称：针对源分离的地震勘测优化的制作方法
技术领域：
本发明总体上涉及针对源分离的地震勘测优化。
背景技术：
地震勘探涉及为寻找碳氢化合物储藏而勘测地下的地质构造。勘测典型地涉及在 预定的位置布置地震源和地震传感器。所述源产生地震波，该地震波传播到地质构造中，沿 其路线产生压力变化和振动。地质构造的弹性属性的变化使地震波散开，改变了它们的传 播方向和其它属性。由所述源发射出的部分能量到达地震传感器。一些地震传感器对压力 变化敏感(水听器)，其它地震传感器对粒子运动敏感(例如，地震检波器)，并且工业勘测 可以只布置一种类型的传感器或者两种都布置。响应检测到的地震事件，传感器产生电信 号，以产生地震数据。然后，对地震数据的分析可以指示存在或者不存在可能的碳氢化合物 储藏位置。一些勘测称为“海洋”勘测，因为它们是在海洋环境中进行的。但是，“海洋”勘测 不仅能在咸水环境下进行，而且也可以在淡水或者半咸水环境下进行。在称为“拖曳阵列” 勘测的一种海洋勘测类型中，包含地震传感器的拖缆和源的阵列被拖曳在勘测船后面。

发明内容
在本发明的一种实施方式中，一种技术包括确定表征地震勘测的至少一个参数， 在所述地震勘测中使多个发生干扰的地震源发射，而地震传感器感测由地震源产生的能 量。参数的确定包括优化该地震勘测，以便分离所感测到的对应于地震源的能量。在本发明的一种实施方式中，一种系统包括地震源和地震传感器，所述地震传感 器适于感测由地震源的发射所产生的能量。对该系统进行优化，以便分离所感测到的对应 于地震源的能量。在本发明的又一种实施方式中，一种系统包括存储器和处理器。所述存储器存储 指令，当所述指令被所述处理器执行时，使所述处理器处理表征地震勘测的至少一个参数， 以便优化地震勘测，从而分离所感测到的由多个发生干扰的地震源产生的能量。从以下附图、描述和权利要求，本发明的优点和其它特征将变得显而易见。


图1是根据本发明的实施方式基于海洋的地震获取系统的示意图。图2、3和4是示出根据本发明的实施方式设计地震勘测的技术的流程图。图5是根据本发明的实施方式的数据处理系统的示意图。
具体实施例方式图1示出了根据本发明的一些实施方式、基于海洋的地震数据获取系统的实施方 式10。在系统10中，勘测船20在船20的后面拖曳一条或多条地震拖缆30(图1中示出了一条示例拖缆30)。应当指出，拖缆30可以展开布置，其中多条拖缆30在相同的深度在大 致相同的平面中拖曳。作为另一个非限制性例子，拖缆可以在多个深度(例如以上/下分 布的形式)拖曳。地震拖缆30可以有几千米长而且可以包含各种支持电缆(未示出)及可以用于 支持沿拖缆30的通信的配线和/或电路(未示出)。总的来说，每条拖缆30包括主线缆， 将记录地震信号的地震传感器安装到该主线缆中。拖缆30包含地震传感器58，依赖于本 发明的特定实施方式，地震传感器58可以是获取压力数据的水听器(作为一种非限制性例 子)或者多部件传感器。对于其中传感器58是多部件传感器(作为另一种非限制性例子) 的本发明的实施方式，每个传感器都能够检测压力波场及与接近传感器的声学信号关联的 粒子运动的至少一个分量。粒子运动的例子包括粒子移位的一个或多个分量、粒子速度的 一个或多个分量(纵测线(inline) (χ)、横测线(crossline) (y)和垂直(ζ)分量(例如，见 轴59))和粒子加速度的一个或多个分量。依赖于本发明的特定实施方式，多部件地震传感器可以包括一个或多个水听器、 地震检波器、粒子移位传感器、粒子速度传感器、加速计、压力梯度传感器或者其组合。例如，根据本发明的一些实施方式，特定的多部件地震传感器可以包括用于测量 压力的水听器和用于测量传感器附近的粒子速度和/或加速度的三个对应正交分量的三 个正交对准的加速计。应当指出，依赖于本发明的特定实施方式，多部件地震传感器可以实 现为单个设备(如图1所示)或者可以实现为多个设备。特定的多部件地震传感器还可以 包括压力梯度传感器，所述压力梯度传感器构成了另一种类型的粒子运动传感器。每个压 力梯度传感器测量在特定点处、关于特定方向的压力波场的变化。例如，一个压力梯度传感 器可以获得指示在特定点的关于横测线方向的压力波场的偏导数的地震数据，而另一个压 力梯度传感器可以获得在特定点的指示关于纵测线方向的压力数据的地震数据。海洋地震数据获取系统10包括地震源40(图1中示出了两个示例地震源40)，例 如气枪，等等。在本发明的一些实施方式中，地震源40可以耦合到勘测船20，或者被其拖 曳。另选地，在本发明的其它实施方式中，地震源40可以独立于勘测船20运行，其中地震 源40可以耦合到其它的船或者浮标，这只是举了一些例子。当地震拖缆30被拖曳在勘测船20后面时，通常被称为“射击”的声学信号42 (图 1中示出了示例声学信号42)由地震源40产生，并向下通过水柱44引导到水底面M之下 的地层62和68中。声学信号42从各种地下的地质构造(例如图1中所示的示例地质构 造65)反射。由源40产生的入射声学信号42产生对应的反射声学信号，或者压力波60，该信号 被地震传感器58感测。应当指出，被地震传感器58接收和感测的压力波包括无反射地传 播到传感器58的“上行”压力波及由压力波60从空气-水边界31反射产生的“下行”压 力波。地震传感器58产生称为“轨迹(trace)，，的信号(例如，数字信号)，该信号指示 所获得的压力波场和粒子运动的测量结果。根据本发明的一些实施方式，轨迹被记录并可 以至少部分地被布置在勘测船20上的信号处理单元23处理。例如，特定的地震传感器58 可以提供轨迹，该轨迹对应于由其水听器的对压力波场进行的测量结果；而且(依赖于本 发明的特定实施方式)传感器58可以提供对应于粒子运动的一个或多个分量的一个或多个轨迹。地震获取的目的是建立勘测区域的图像，以便识别地下的地质构造，例如示例地 质构造65。对图像的后续分析可以揭示地下地质构造中碳氢化合物储藏的可能位置。依 赖于本发明的特定实施方式，对图像的部分分析可以在地震勘测船20上执行，例如通过信 号处理单元23执行。根据本发明的其它实施方式，所述图像可以由地震数据处理系统来处 理，其中地震数据处理系统可以例如位于陆地上或者船20上。因此，许多变化都是可能的， 而且都在所附权利要求的范围之内。特定的地震源40可以由地震源元件(例如气枪)的阵列形成，这些元件可以布置 在阵列的各行(例如，枪的行)中。另选地，特定的地震源40可以由阵列中的一个或者预 定数量的气枪形成，也可以由多个阵列形成，等等。不管地震源的特定组成如何，源都可以 在勘测过程中以特定的时间序列发射。在获取地震数据时有许多的物理约束，例如记录长度与获取效率的关系。更具 体而言，由地震传感器感测到的、由于地震源的给定发射或者“射击”所造成的能量一般 是在以下的时间间隔记录的，所述时间间隔从射击发生的时间开始，并在下次射击发生的 时间稍前一点结束。为了提高效率，例如在1999年7月13日发布且标题为“METH0DSF0R ACQUIRING AND PROCESSING SEISMIC DATA” 的美国专利 No. 5，9 , 049 中所公开的技术允 许多个地震源同时发射，并且基本上允许同时记录多于一个的记录。因此，地震源40可以按一定顺序发射，使得多个地震源40可以同时或者在短时间 间隔内几乎同时发射，从而由地震传感器58感测到的合成能量包含来自多于一个的地震 源40的显著量的能量。换句话说，地震源彼此干扰，使得合成能量信号不容易分离成归属 于特定地震源的信号。但是，源分离技术可以用来处理所获得的地震数据，以便形成每个 都与一个地震源40关联的数据集，使得每个数据集都理想地指示合成地震能量信号中归 属于关联地震源40的分量。源分离一般并不完美，使得一些感测到的能量(称为“残余能 量”)最终不能归属到任何一个地震源，并且一些感测到的能量(称为“泄漏能量”)可能被 分离处理归属到错误的地震源。地震勘测具有多个特征参数，这些参数影响源分离的质量。更具体而言，例如支配 地震源发射的定时序列、源的几何结构(例如，气枪的横测线和纵测线间距)和接收器的几 何结构(作为例子，分布的类型，以及地震传感器58的横测线和纵测线间距)的这类参数， 可能影响源能量多有效地分开。根据在此描述的本发明的实施方式，在图2中所示的技术 150用于设计地震勘测的目的。参考图2，根据技术150，设计地震勘测(块154)，其中使多 个发生干扰的地震源发射。为了分离感测到的与地震源对应的源能量，对获取的数据进行 处理。根据技术150的块158，优化地震勘测，以便分离源能量。换句话说，技术150包括为 了优化源分离而确定勘测的一个或多个参数。如在此进一步描述的，优化的勘测参数可以是与地震源的几何结构相关的参数， 例如地震源的个数、横测线源间隔、纵测线源间隔、具体的源位置，等等；接收器几何结构； 地震源相对于彼此发射的时间(即，源发射的定时序列)；源发射的定时序列和频率之间的
关系；等等。尽管作为例子在此描述了拖曳的海洋地震勘测，但是应当理解，在此描述的技术 和系统还可以同样地应用到任何其它类型的具有发生干扰的地震源的勘测，例如非拖曳的海洋勘测、基于陆地的勘测、基于海底线缆的勘测、震动源勘测，等等。例如，根据本发明的 其它实施方式，例如震动源勘测的源频率、源振幅和源发射定时序列的这类参数可以针对 源分离而进行优化。作为另一个例子，根据本发明的实施方式，为了优化其中地震源和/或地震传感 器可以设置于井孔中的钻孔勘测系统，可以应用在此描述的系统和技术。因此，许多变化都 是预期的，并且在所附权利要求的范围之内。图3示出了根据本发明的一些实施方式用于设计地震勘测的系统200。依赖于本 发明的特定实施方式，系统200可以由在基于一个或多个处理器的系统上执行的软件单独 实施，或者可以作为软件和硬件的组合来实施。当勘测的地况已知(例如，已经从先前勘测 的结果获知的地况)或者至少在可以确定对勘测的地况的合理估计时，可以使用系统200。 例如，根据本发明的一些实施方式，系统200可以基于ID估计来使用，其中ID估计假设地 况仅仅在垂直方向变化而不在水平方向变化。为了针对源分离而优化勘测参数，系统200包括应用数字处理技术的基于计算机 的模拟器220。更具体而言，根据本发明的一些实施方式，基于计算机的模拟器220执行蒙 特卡罗模拟，该模拟基于随机或者伪随机产生的输入来对勘测系统建模。但是，模拟和输入 可能受到各种约束。因此，为了确定用于源分离的最优勘测参数，例如发射时间、源几何结 构、接收器几何结构(即，地震传感器的几何结构)等这类参数都可以在预定的范围内随机 地或者伪随机地变化。更具体而言，如作为例子在图3中所描述的，基于计算机的模拟器220可以接收指 示勘测的地况的数据216及分别指示源几何结构、接收器几何结构和发射时间的数据214、 215和218。应当指出，依赖于本发明的特定实施方式，源几何结构数据214、接收器几何结 构数据215和发射时间数据218可以分别由随机、伪随机或者不随机的源几何结构206、接 收器几何结构205和源发射时间204发生器产生。因此，根据本发明的一些实施方式，发生 器204、205和206可以产生预定范围之内随机的源几何结构、接收器几何结构和发射时间。依赖于本发明的特定实施方式，产生源几何结构、接收器几何结构和发射时间的 方式可以变化。作为例子，在确定最优的发射时间时源和接收器几何结构可以保持恒定；源 几何结构、接收器几何结构和发射时间可以同时随机变化；在确定最优的源几何结构参数 时发射时间和接收器几何结构可以保持恒定；等等。因此，许多变化都是预期的，并在所附 权利要求的范围之内。基于计算机的模拟器220产生合成数据集224，该数据集是预测为实际勘测中要 由地震传感器获取的地震数据集，其中勘测是由当前的勘测地况和作为输入被基于计算机 的模拟器220接收的勘测参数定义的。合成数据集2M指示预测的所感测的由发生干扰的 地震源40产生的合成能量信号。基于合成数据集224，系统200的源能量分离器230产生 多个数据集，每个数据集都归属于一个特定的地震源。应当指出，源能量分离器230可以采用数值反演(inversion)算法(通过在计算 机上执行的指令来执行)，该算法涉及线性系统的反演(作为非限制性例子)。源能量分离 可能不能够将所有的能量都归属到地震源中的一个，这意味着由源能量分离器230进行的 处理会产生残余能量。而且，存在泄漏能量，这是由一个地震源产生但被分离处理归属到不 同地震源的能量。从数学上讲，如果所记录的数据是d = C^d2，且所估计的数据是d/和d2’，则残余是Cl1，+d2，-d ；对于地震源S1泄漏是Cl1，-Cl1 ；对于地震源S2泄漏是d2，-d2。因此，残余和泄漏能量指示源分离中的误差程度，理想地，残余和泄漏能量为零， 因为，理想地，所有感测到的能量都在地震源之间划分并正确地划分到各地震源。因此，根 据本发明的一些实施方式，当地震勘测由最小化残余和泄漏能量或者至少产生低于所选阈 值的残余和泄漏能量的一组参数表征时，发生地震勘测的优化。根据本发明的实施方式，系统200的控制器240从源能量分离器230接收指示残 余能量的数据234。控制器240 (作为非限制性例子，通过改变源发射时间或者源几何结构) 改变一个或多个勘测参数，直到基于残余能量数据234和泄漏能量数据235，控制器240确 定这些能量已经被充分地最小化了。例如，控制器240可以继续运行接收器几何结构205、 源几何结构206和发射时间204的发生器，控制器240可以改变输入范围，控制器240可以 保持一些输入恒定而改变其它输入，等等。一旦控制器240确定残余和泄漏能量已经充分 地最小化了，控制器240就提供指示最优勘测参数的数据244。在一些情况下，勘测的地况可能是未知的或者可能无法获得对地况的可靠估计。 对于这种情况，基于表征勘测系统的线性系统，可以对勘测参数进行优化以用于源能量分 离。作为更具体的例子，参考图4，根据本发明的一些实施方式，优化勘测参数的技术250包 括将地震数据集建模(块254)成可以从由勘测的地况与勘测的参数表征的线性系统确定。 每个数据集与感测到的唯一地归属于不同地震源的能量关联。该线性系统可以针对地震数 据集进行反演。根据块258，通过最大化用于数据集恢复的反演的精度，对勘测参数进行优 化以用于源分离。作为更具体的例子，根据本发明的一些实施方式，源颤振(即，利用其中源发射时 间彼此之间有轻微偏移的源发射定时序列)可以用于使发生干扰的地震源的分离，如在 2007 年 12 月洸提交且标题为 “SEPARATING SEISMIC SIGNALS PRODUCED BYINTERFERING SEISMIC SOURCES”的美国专利申请序列No. 11/964，402(律师卷号57. 0820)中所描述的， 该申请的全部内容通过引用包含于此。为了将所感测到的合成能量信号分解成每个都唯一 地与特定地震源关联的信号，线性算子变换可以在勘测系统的模型中使用。为了从源分离 的角度确定最优的勘测设计，该系统矩阵的可逆性(irwertablity)被最大化。更具体而言，地震数据(在此被称为“地震数据向量d”)被认为是由于N个(S卩，多 个)地震源的发射而由地震传感器获取的。因此，地震源的同时或者几乎同时发射造成来 自所有这些发射的重要能量都存在于地震数据向量d中。描述影响源能量的地况的模型与 描述源机制的物理性质的线性算子、波传播和勘测几何结构关联。于是，地震数据向量d可 以表征为该模型与线性算子的函数。因此，理论上，该函数可以对模型联合反演(invert)，
这允许将地震数据向量d分离成N个地震数据集屯、d2、d3.....dN，使得每个数据集都唯一
地归属于一个地震源。换句话说，每个数据集都代表所感测到的合成能量信号的一个分量， 该分量唯一地归属于一个地震源。作为更具体的例子，假设地震数据向量d是由于称为“S/’和"S2”的两个地震源几 乎同时发射而获取的。对于这个例子，地震源S1和&根据一个定时序列发射，该定时序列 可以基于预定的定时模式或者可以基于随机或伪随机的时间。不管具体的定时方案如何， 对这个例子假设对于所有轨迹地震源S1都是在地震源&之前发射的，而且还假设轨迹的零 时刻对应于S1的发射时间。因此，轨迹的零时刻在"S1的时间”中。对地震源S1和&的偏移或者向量，分别称为“X1”和“X2”。地震源&的定时延迟，标记为“t”，对于每条轨迹都是 已知的。对这个例子，假设轨迹的集合使得t的值为随机的。在实践中，对CMP、接收器或者 共同的偏移集合来说，情况就是这样。为了简化这种讨论，假设每个集合中的轨迹可以分别 利用称为"X1i”和“X2/’的标量关于地震源S1和地震源&定位。在这种标记法中，下标“i” 指示集合中的轨迹编号。作为更具体的例子，对于CMP集合，“χ1/’可以是相对于地震源S1 的标量偏移，而且这些量在下面都称为偏移。类似地，“、”指示对于第i条轨迹的定时延 迟。对地震源S1所记录的能量可以通过向称为“m/’的未知模型施加称为“L/’的线性 算子(该算子代表地震源S1的物理性质，与源S1关联的波传播和与地震源S1关联的勘测 几何结构)来建模，其中“m/’描述了影响从地震源S1传播的能量的地况。模型HI1对模型 空间中的每个参数包含一个元素。一般来说，模型空间可以通过慢度(slowness)或者其平 方来参数化，这分别对应于线性或者双曲线/抛物线Radon变换。线性算子L1是对于源S1 的偏移、表征模型空间的参数和时间或频率的函数。地震数据向量Cl1对每条轨迹都包含一 个元素(在每个时间或者频率)，而且是地震数据d中与地震源S1关联的分量。换句话说， 地震数据向量Cl1表示归属于地震源S1的数据集。地震数据向量Cl1可以如下描述Cl1 = L1Hi1公式 1与地震源&关联的能量在地震数据向量d中不连贯地出现。但是，该能量与连贯 的数据集相关，其中，通过对轨迹应用时间偏移ti;地震源&的发射时间是在时间零(即， 地震源&时间)。称为“D2”的对角线线性算子可以用于描述这些时间偏移，使得地震数据 向量d中与地震源&关联且称为“d2”的分量可以如下描述d2 = D2L2Hi2 公式 2在公式2中，称为"L2”的线性算子表示地震源&的物理性质，与地震源&关联的 波传播和与地震源&关联的勘测几何结构。而且，在公式2中，称为“m2”的模型描述了影 响从地震源&传播的能量的地况。由地震传感器记录的合成地震能量信号归属于地震源S1和&。因此，地震数据向 量d(即，所记录的数据)是地震数据向量Cl1和d2的组合，如下所述d = Cl^d2公式 3由于公式1、2和3中的关系，地震数据向量d可以表示为如下的线性系统
mI
公式4 d = [L1 D2L2.
W1因此，公式4可以利用如最小二乘方算法的标准技术对模型向量m(即，Oii1 ;m2)) 求解；并且，在知道模型向量m之后，公式1和2可以应用模型Hi1和m2，以便将地震数量向 量d分离成地震数据向量Cl1和d2，即，分离成指示归属于每个地震源的测量结果的数据集。公式4可以在频率(ω)域中反演。在那种情况下，(D2)jk = MpGicotj) δ jk且 (Ls) Jk = exp (- ω tsJk)，其中t、是在与&关联的模型空间中与偏移力和用于第k条轨道 的参数关联的时间偏移。对于由慢度参数化的线性fcidon变换psk，tsjk = XsjPV对于由曲
率参数化的抛物线fcidon变换qsk，ts1k = (Xs1)V
k，L jk νΛ j/ 4 k°
8
以上描述的源分离技术的成功依赖于用于分离与两个源关联的能量的变换的能 力。不像Radon变换的大多数应用那样，成功不依赖于将对于正确模型参数的能量聚集到 Hi1或者m2中的能力。当随机或者伪随机的时间延迟用在源发射之间时，用于两个模型域的 基本函数(〖\和tj+tsjk)是非常不同的，而且这实现了源的极其有效的分离。根据本发明的实施方式，为了优化勘测以便分离源，确定最大化矩阵L= [L1 D2L2] (见公式4)的可逆性的勘测参数。应当指出，如本领域技术人员可以认识到的，矩阵可逆 性的程度可以由例如评估赫赛函数(Hessian)、矩阵的条件数、矩阵的稀疏程度和矩阵的特 征值分布这样的技术来评估。因此，最优的勘测参数(源发射的定时序列、源几何结构，等 等)是最大化矩阵可逆性的参数。应当指出，使L可逆对于优化而言是足够的，但不是必需的。如果L是可逆的，则
d = Lm可以对m进行求解，其中
m=
mI ^m 2
但是Hl1和Hl2不需要明确地重新获得。如由公式
1和2给出的，重新获得Cl1和d2就足够了。可以重新获得“错误的”Hi1和m2，但是重新获得 正确的Cl1和d2。假如HI1建模了与地震源S1关联的所有能量，而m2建模了与地震源&关联 的所有能量，等等，即使它在模型中处于错误的位置，也会发生这种情况。总体上说，最大化 L的可逆性是重要的。 从数学上讲，如果M定义为
Mj
M
2 J
即L的通用逆，则&
，Ui,
的估计如下
L1 0 YMj
O D2L2
J\
M2
(L1 m} +D2 L2 公式 5 数据Cl1被认为具有L1Hi1的形式等，因此M1和M2满足以下关系
Li m}
L1MjL1
L1M1D2L2
D2L2 m2 J [D2L2M2Lj D2L2M2D2L2 j
或者
fLl(MlLl-I) D2L2M2L1
Jn 2 j
公式6
LjMlD2L2 D2L2(M2D2L2-I)
V
mI m2j
公式7公式7 对于所有的 Hi1 和 m2 都满足，因此 L1 (M1L1-I) = 0、L1M1D2L2 = 0、D2L2M2L1 = 0且D2L2 (M2D2L2-I) =0，其中“I”代表单位矩阵。这比使M成为L的逆更通用且更容易，其 中使M成为L的逆需要M1L1 = KM1D2L2 = (KM2L1 = 0和M2D2L2 = I。例如，M1D2L2需要在L1 的零空间中，但不需要为零。尽管通过例子在此描述了线性系统，但是应当指出，根据本发明的其它实施方式， 勘测可以利用非线性系统进行优化以用于源分离。参考图5，根据本发明的一些实施方式，为了优化地震勘测以便分离源，数据处理 系统320可以执行在此所公开的一种或多种技术的至少一些部分。根据本发明的一些实施 方式，系统320可以包括处理器350，例如一个或多个微处理器和/或微控制器。依赖本发 明的特定实施方式，处理器350可以位于拖缆30上(图1)、位于船20上或者位于基于陆地 的处理设施上(作为例子)。为了接收如所获得的地震数据、地况测量数据、地况估计数据、指示对勘测参数的
9约束的勘测数据、勘测参数的范围等的这类数据，处理器350可以耦合到通信接口 360。根 据在此描述的本发明的一些实施方式，当执行存储在地震数据处理系统320的存储器中的 指令时，处理器350可以实施在图3中所示的一个或多个处理块或者可以帮助执行图2和 4中所示的一个或多个处理步骤。作为例子，通信接口 360可以是通用串行总线(USB)接口、网络接口、可移除介质 (例如闪存卡、CD-ROM，等等)接口或者磁存储器接口(IDE或者SCSI接口，作为例子)。因 此，依赖于本发明的特定实施方式，通信接口 360可以采取多种形式。根据本发明的一些实施方式，通信接口 360可以耦合到系统320的存储器340并 且可以存储例如在确定最优勘测数据中所涉及的各种输入和/或输出数据集。根据本发明 的一些实施方式，存储器340可以存储程序指令344，当程序指令344被处理器350执行时， 程序指令344可以使处理器350执行在此所公开的一种或多种技术和系统的各种任务，其 中的技术和系统例如技术150或250和系统200，并在系统320的显示器(在图5中未示 出)上显示通过该技术/系统所获得的结果。其它实施方式也在所附权利要求的范围之内。例如，尽管以上已经描述了拖曳的 基于海洋的地震获取系统，但是，在此描述的优化勘测以进行源分离的技术和系统也可以 同样地应用到其它类型的地震获取系统。作为非限制性的例子，在此描述的技术和系统可 以应用到基于海床、钻孔和陆地的地震获取系统。因此，依赖于本发明的特定实施方式，地 震传感器和源可以是固定的或者可以被拖曳。作为本发明其它实施方式的附加例子，除以上描述那些之外的技术也可以用于优 化勘测以进行源分离。作为例子，如本领域技术人员可以认识到的，在本发明的其它实施方 式中，如遗传算法、模拟退火算法、最速下降算法或共轭梯度算法这样的技术可以用于确定 最优的勘测参数。尽管本发明已经关于有限数量的实施方式进行了描述，但是受益于本公开内容的 本领域技术人员将认识到根据其进行的多种修改和变化。希望所附权利要求覆盖属于本发 明真正主旨和范围类的所有这类修改和变化。
权利要求
1.一种方法，包括确定地震勘测的至少一个特征参数，其中使多个地震源发射并且地震传感器感测由所 述地震源产生的能量，所述确定包括优化所述地震勘测，以便分离所感测到的与地震源对 应的能量。
2.如权利要求1所述的方法，其中所述优化的动作包括 优化支配所述地震源发射的定时序列；或者优化所述地震源的几何结构；或者 优化所述地震传感器的几何结构。
3.如权利要求2所述的方法，其中所述优化的动作包括优化所述地震源或者所述地震 传感器的纵测线间隔或者横测线间隔。
4.如权利要求1所述的方法，其中所述地震勘测包括震动源勘测、海洋勘测、拖曳勘测 或者钻孔勘测。
5.如权利要求1所述的方法，其中所述优化包括 确定与勘测关联的地况；基于所确定的地况，预测勘测中要由地震传感器获得的数据集；基于预测的数据集，确定在能量分离之后不与地震源中任何一个关联的残余能量；基于预测的数据集，确定与错误地震源关联的泄漏能量；以及确定所述至少一个参数，以最小化所述残余能量和所述泄漏能量。
6.如权利要求5所述的方法，其中所确定的地况包括估计的地况。
7.如权利要求5所述的方法，其中所确定的地况包括从勘测确定的地况。
8.如权利要求1所述的方法，其中所感测到的能量是通过针对地震数据集反演线性系统公式来分离的，每个数据集归属 于一个源，以及所述优化的动作包括确定所述至少一个参数，以最大化数据集反演的准确度。
9.如权利要求8所述的方法，其中所述确定所述至少一个参数的动作包括确定矩阵的 条件数、稀疏程度或者特征值分布。
10.如权利要求1所述的方法，其中所述优化的动作包括将地震勘测建模为非线性系统。
11.一种系统，包括 地震源；及地震传感器，适于感测由所述地震源发射而产生的能量， 其中，使用如权利要求1-10中所述的方法对该系统进行优化以执行勘测。
12.如权利要求11所述的系统，还包括 拖曳所述地震传感器和所述地震源的船。
13.一种系统，包括 存储指令的存储器；以及执行所述指令的处理器，所述指令执行如权利要求1-10中所述的方法，以便优化地震 勘测，从而分离由发生干扰的地震源产生的能量。
全文摘要
一种技术，包括确定表征地震勘测的至少一个参数，在地震勘测中使多个发生干扰地震源发射且地震传感器感测由所述地震源产生的能量。参数的确定包括优化该地震勘测，以便分离所感测到的与地震源相对应的能量。
文档编号G01V1/16GK102138082SQ200980133743
公开日2011年7月27日 申请日期2009年7月14日 优先权日2008年7月16日
发明者C·J·贝斯利, I·穆尔, 小W·H·德拉戈塞特 申请人:格库技术有限公司