CN1488078A - 对源编码和谐波消除采用级联扫描的方法 - Google Patents

对源编码和谐波消除采用级联扫描的方法 Download PDF

Info

Publication number
CN1488078A
CN1488078A CNA018190588A CN01819058A CN1488078A CN 1488078 A CN1488078 A CN 1488078A CN A018190588 A CNA018190588 A CN A018190588A CN 01819058 A CN01819058 A CN 01819058A CN 1488078 A CN1488078 A CN 1488078A
Authority
CN
China
Prior art keywords
signal
seismic
scanning
scanning section
sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CNA018190588A
Other languages
English (en)
Other versions
CN1254693C (zh
Inventor
瑞纳・莫瑞格
瑞纳·莫瑞格
里克・J・巴尔
弗雷德里克·J·巴尔
李・耐兰德
戴维·李·耐兰德
亚瑟・斯通
加里·亚瑟·斯通
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Westerngeco Seismic Holdings Ltd
Original Assignee
Westerngeco Seismic Holdings Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Westerngeco Seismic Holdings Ltd filed Critical Westerngeco Seismic Holdings Ltd
Publication of CN1488078A publication Critical patent/CN1488078A/zh
Application granted granted Critical
Publication of CN1254693C publication Critical patent/CN1254693C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design
    • G01V1/005Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/37Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy specially adapted for seismic systems using continuous agitation of the ground, e.g. using pulse compression of frequency swept signals for enhancement of received signals
    • G01V1/375Correlating received seismic signals with the emitted source signal
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering
    • G01V2210/3246Coherent noise, e.g. spatially coherent or predictable

Abstract

一种利用一个或多个由扫描信号驱动的震动地震能源的地震勘测方法。确定强度足以导致一扫描段的明显谐波畸变的最高阶谐波。选择超过源的数量的扫描段数量。为每个地震能源的每个扫描段选择初始相位角,从而抑制掉直到并且包括该最高阶谐波的所有谐波。利用级联的扫描,获得并处理基本衰减掉直到该选定阶数的谐波的地震数据。

Description

对源编码和谐波消除采用级联扫描的方法
技术领域
本发明涉及地球的物理勘探,尤其涉及采集地震震动器数据。
背景技术
地震勘探业利用声脉冲把声振动传递到地球以便为矿物勘探和开采描绘地下结构。这些声脉冲可以来自爆炸、内破裂、扫频(线性调频脉冲)或随机源。利用对从该源传播到接收器的声反射和折射波前的记录产生地震场记录。这些场记录中的反射和折射事件的传播时间变化指示地球内反射表面的位置。在地震数据处理中对一个或多个场记录里的事件的分析及相关产生表明地下结构的声像。利用这些声像寻找有价值的矿床。
扫频或线性调频脉冲型地震源可以采用相对长的例如2到15秒的导频信号以确保对地球施加足够的能量。扫频或线性调频脉冲型源方法依赖于信号压缩以便压缩信号并且确保足够的垂直分辨率,从而解出地下反射物的位置。信号压缩通常被称为去卷积,它的许多技术在地震数据处理中是周知的。扫描或线性调频脉冲信号的去卷积把源信号压缩成代表地下反射边界的短得多的信号。任何去卷积技术的精度和有效性直接依赖于对源信号的知道或了解程度。大多数去卷积运算子是从实际源波形的统计估计导出的。
利用扫频型源,能量是在地震频率范围内按有规律地增加(向上扫描)或降低(向下扫描)频率的扫描形式发射的。可以通过能够控制地震信号的频率和相位的控制信号控制振动。
当只利用一个地震源时,地震勘测可能是非常耗时的。利用现代信号处理方法,如果可以同时使用多于一个的地震源,则可以缩短该周期。如果能提供某种区分不同源发出的信号的手段,则可以使用多个源。另外,振动源产生谐波,在某些情况下这些谐波可能具有接近甚至超过基频的能量,并且它们可以和来自其它源的信号串馈,从而当处理这些信号以便分开来自每个源的信号时会给出错误结果。此外,这些谐波是噪音源并且可能掩盖来自更深层的弱反射信号。
授予Anstey等的美国3885225号专利提出一种区分多个源的方法。Anstey等的专利涉及一种利用数个同时发出信号的震动器进行宽线地震剖面测定的方法和设备。把额定发射频率带宽划分成数个部分,把这些部分分配给各个按顺序独立发射的震动器,并且使这些震动器在任何时候发射互相排它的频率。根据代表来自每个震动器的相应信号的频率分离检测到的信号。对每个震动器的频率限制降低勘测的灵敏性。各个震动器中的谐波畸变或者它们和大地的耦合会抑制对来自不同发生器的信号的分离。该225号专利未试图去掉或减小该畸变。
在授予Landrum的美国4715020号专利中公开另一种利用不同扫描上的信号相位偏移分离来自多个震动源的信号的方法。但是,该专利不解决非线性畸变和串馈问题。
如美国4042910号专利中公开那样,Reitsch得出一种衰减因各地震震动器的谐波能量输出产生的谐波相关噪音的方法。该方法包括顺序地生成多个扫描信号的步骤,其中每个相继的扫描信号的相位相对于前一个扫描信号偏移预定的按2B的几分之一的相位角。在按常规方式相加或堆叠之前,通过反向相位偏移分别记录和变换这些生成的信号。利用该方法,通过多个震动器(每个记录一个)输出一组只是相位不同的N个扫描段。直至第N个谐波的相关噪音被衰减。该方法利用相位偏移提供一种抑制谐波的方法,但是只对单个震动源,而且不解决串馈问题。
授予Ward的美国4823326号专利请求保护一种具有减小的谐波畸变的对从多个同时运行的震动地震源导出的地震记录进行分离的方法。
授予Edington和Khan的美国4982374号专利公开了一种方法,该方法为任何数量的同时运行的震动地震源减小任何选定价数的谐波的畸变和串馈,同时提供对来自不同源的信号的分离,而且用于改进信噪比。在确定可能引起畸变的最高阶谐波后,选择每个位置处的每个源的扫描数量。该数量取决于源的数量以及要抑制的最高阶谐波。接着为每个源的每个扫描选择初始相位角以允许信号分离并且同时抑制直到并且包括该最高价谐波的各阶谐波。
授予Andersen的美国5410517号专利公开一种级联扫描段以便抑制不希望的谐波能量的方法。该方法使用具有不同相位角的扫描段从而衰减相关弱波中的谐波能量。依据该方法,生成包含N个首尾相连的扫描段的第一级联扫描序列。除了该扫描序列内的每个扫描段的初始相位角按约为360/N度的固定相位增量累进转动外,这N个扫描段基本相同。第二级联扫描序列按如下方式生成:(i)N个首尾相连的与所述第一级联扫描序列对应的连续扫描段,以及(ii)一个和这N个连续扫描段连接的附加扫描段,其被定位并且定相成在相关期间大致抑制掉谐波叠影(ghost)。这些级联扫描序列之一用作为震动器扫描序列而另一个用作为相关基准序列。
发明内容
一种利用通过可级联在一起的扫描信号激励的震动地震能源的地震勘测方法。选择若干地震源,选择能造成明显畸变的最高阶谐波能,确定并且选择若干要组合在一起的扫描或扫描段以便抑制不希望的谐波并且分离源信号。通过在相关前测到的地面力信号能调整相关运算子的幅值,以便平衡和衰减谐波能。该相关运算子可以是导频扫描信号、地面力信号或者导频或地面力的派生分量。
附图说明
通过参阅下面的详细说明和各附图会更好地理解本发明和它的优点,附图中:
图1示出本发明的总结构以及它的各实体。
图2a示出一个带有相角转动的级联扫描。
图2b示出含有谐波能的场记录。
图3a示出一个带有可用于抑制直到五阶谐波的相角转动的级联扫描序列。
图3b是本发明的一实施例的示意说明,其中相关信号包括由一个扫描段扩充的扫描信号。
图3c是本发明的一实施例的示意说明,其中相关信号包括不扩充的扫描信号。
图4示出应用本发明后的场记录。
图5示出可用于二个实现并行数据采集的源的二个扫描序列。
图6示出应用本发明之前的通过二个源获得的场记录。
图7a示出本发明的用来从图6中显示的并行采集的一个源中分离信号的应用。
图7b示出本发明的用来从图6中显示的采集的第二个源分离信号的应用。
图8示出借助导频信号的震动器信号的相关。
图9示出借助地面力信号的震动器信号的相关。
图10示出对谐波抑制的扫描相位转动后与地面力的相关。
图11示出在添加噪音情况下在相位转动扫描段组上与导频扫描的相关。
图12示出添加噪音情况下在相位转动扫描段组上与地面力信号的相关。
图13示出当使用地面力信号时的相关处理的示意图。
尽管通过它的各优选实施例说明本发明,应理解本发明不受它们的限制。相反,本发明意图复盖所有包含在本发明的由附后的权利要求书定义的精神和范围之内的替代、修改和等同物。
具体实施方式
通过参照图1可最佳地了解本发明的总结构。这会有助于理解后面的详细说明。如图1中所示的本发明的方法利用级联扫描获取地震数据1。本发明的目的是能够从多个源中分离出数据以及衰减谐波3。在本发明的一实施例中,这是通过使反射的地震信号和一理论的预定信号相关实现的,该理论的预定信号和用于驱动震动器的导频扫描有关5。这可在任何数量的源和任何所希望的谐波清除级别下完成9。概括地说,当相关信号和导频扫描有关时,每个震动器用N个扫描段扫描并且相关信号由N+1个段构成13。
在本发明的一替代实施例中,在每个震动器处测量地面力信号,并且相关信号和地面力信号有关7。如第一方法一样,可能使用任何数量的源以及任何程度的谐波消除11。和第一实施例一样,用N个扫描段扫描震动器;但是,相关信号包括N+2个段的记录的地面力数据。这将在后面说明。
本发明是一种地震数据采集方法,其采用可以对多个同时运行的源使用的相位转动扫描和级联相位转动扫描段。本发明在同时使用多个源的情况下提供对谐波叠影以及串馈能量的抑制,以便减小勘测采集时间并改进数据质量。从而,下面的说明的范围特别涉及本发明的特定实施例或特定使用,但它只是说明性的而不是对本发明范围构成限制。
当场效应涉及把通常在震动器点(VP)输出的扫描组合成例如图2a中所示的扫描段级联序列时,震动器工作人员的工作会更加有效。当使用一个源时,该组带有相位序列0、72、144、216和288度的五个级联扫描段可用于抑制直到第五阶的谐波(二阶、三阶、四阶和五阶)。记录这五个扫描段的数据所需的场时间远短于分别记录它们所需的时间,因为只需要单个收听时间。但是,即使利用用来消除谐波的相位转动法建立级联的扫描段并且利用相同的级联扫描导频信号进行交叉相关。在相关数据中的每个小弱波之前和之后都产生谐波相关噪音。在图2b的相关记录中由级联扫描段组造成的相关噪音是明显的。
1994年授予Andersen的美国5410517号专利(该专利全文收录为本文的参考文献)公开了一种衰减由图2a的级联扫描段组造成的谐波相关噪音的方法。发现通过向图2a中示出的序列添加另一个扫描段(即在图2a的第一扫描段的前面放上最后一个扫描段的复制品)形成相关导频信号,从而衰减谐波。Andersen的方法和权利要求规定组成级联扫描的各个扫描段的相位转动序列必须为累进的+-2B/N,其中N是扫描段的数量。
本发明将采用多个源时的源编码和谐波消除方法与用于级联相位转动震动器的各扫描段的技术相结合。在本发明的一实施例中,在基准导频或相关信号包括N或N+1个级联扫描段方面,级联方式和美国5410517号专利类似。但是‘517’专利中为级联导频扫描段确定相位排序的方法是严格的,即相位转动是累进的并且只由2B除以扫描段数量来确定。遵循‘517’专利提出的方式的教导,将不具有为实现多个同时运行的源所必须的灵活性。由于不能实现串馈谐波消除要求,‘517’专利的方法实际上排除同时使用多个源的可能性。将级联相位转动扫描段的方法和用于产生源信号分离的源编码技术相结合允许同时使用多个源并且抑制谐波叠影和串馈量。
在本发明的一实施例中,一组或多组震动器产生要比包括N个级联的长度额定的扫描段的常规输出长。若使用多于一组的震动器,它们位于不同的VP。记录时间等于N乘以额定长度的扫描时间加上单个收听时间(相关记录长度)。从而该记录时间要比记录N个扫描段中每一个的收听时间短得多。对N个级联扫描段中的每个进行相位转动以用于源编码和抑制谐波。由于缩短了总勘测时间,利用该方法可使现场工作更加有效。
例如,通过如图3a中示出的相位转动为0、144、288、72和216度的五个扫描段可以利用一组震动器源并且抑制到直到第五阶谐波。利用相位转动为0、102.86、205.71、308.57、51.43、154.29和257.14度的7个扫描段可以实现对直至第七阶谐波的抑制。为了用8个扫描段抑制直到第8阶的谐波,这些扫描段的相位转动序列可以为0、135、270、45、180、315、90和225度。为达到现场工作的效率,震动器扫描将包括N个段并且相关运算子将包括N+1个段。
图3a示出用于本发明的采用一组震动器并消除直至第五阶的谐波的一优选实施例的级联扫描序列。该扫描序列不是累进的。利用图3中示出的扫描段组的交叉相关数据(带有一个附加的扫描)的结果在图4示出,其中这些谐波已经消除。和不对相关运算子级联该附加扫描段的图2b的交叉相关数据相比,该结果显示出谐波已被消除。在图3b中示出用于一个震动器的该实施例的方法。数据段A′、B′、C′和LISTEN对应3个扫描段和一个跟在该扫描的结束处的收听时间。相关信号包括扫描段C、A、B、C。
在本发明的另一优选实施例中,可使用例如图3a的扫描序列获取数据。在数据处理期间,在第一扫描段记录的数据部分可以和扫描序列之后的记录阶段“收听时间”组合。替代地,可以在采集期间使数据记录的该第一扫描部分和现场数据记录的“收听时间”部分组合。将此应用于多个同时运行的源是直接的。在图3c中示出的该方法中,如图3b中那样,存在三个后面跟着一个收听时间LISTEN的数据段A′、B′、C′。在扩大的数据集中,对收听时间期间内的数据添加来自阶段A′的数据并且利用相关运算子A、B、C执行相关。
如果要同时使用多个震动器源,就需要一种分离各个源信号的方法。在美国4982374号专利中公了的Edington和Khan得出的一种相位转动方法,该方法不设想级联的扫描段而是提供同时的源编码并且分离多个源而且消除交叉相关的谐波噪音。美国4982374号专利收录在此作为参考文献。
在本发明中,对每组震动器用级联扫描导频信号来交叉相关记录数据。相关运算子基准信号或震动器扫描信号可包含N或N+1个级联扫描段,但这些扫描段的相位排序被选择成能分离来自同时运行的源的信号。不论是把附加的扫描段放在相关运算子上还是放在震动器扫描序列中,包含该附加段的扫描段序列的第一个和最后一个扫描段将具有相同的初始和最后扫描段。接着交叉相关产生从每组震动器得到的压缩数据并且从而可把各个源信号做成是分离的数据集。源谐波以及级联扫描段二者造成的相关噪音被大大衰减,这使得采用级联的扫描段是可接受的并且可应用于通过多于一个源的同时数据采集。
Edington和Khan的专利教导了一种确定用于信号分离的相位转动的方法。图5示出用于二组震动器的允许进行信号分离并且能消除直到第五阶谐波(即,消除二阶、三阶、四阶和五阶谐波)的级联扫描序列。对于该例,第一震动器或震动器组Vib 1的相位序列是0、102.86、205.71、308.57、51.43、154.29和257.14度。注意,如前面在单个震动器例子中指出的那样,Vib 1序列独自可抑制直到第七阶的谐波。同时的第二扫描震动器或震动器组Vib 2将具有0、257.14、154.29、51.43、308.57、205.71和102.86度的扫描序列相位转动。在同时使用Vib 1和Vib 2级联扫描段的情况下,和单个震动器例子对比只抑制直到第五阶的谐波。图6示出利用扫描序列Vib 1和Vib 2获得的未相关的场记录的第一部分,其指出这二个源位于记录跨距的相对二端处。图7a和7b显示在消除谐波噪音情况下图6的每个震动器源的分离的记录。
扫描段序列已示出是顺序相邻的,但是本发明还可以用附带着重叠的扫描序列段实现。例如,可以在一个10秒长的扫描段尚存2秒结束时开始下一个扫描段。对于包含七个10秒长的扫描的扫描序列,这会在采集工作上显示出明显的时间节约。
构成级联扫描段序列的每个扫描段除了它们各自的相位转动外是相同的。但是,它们可以为任何类型,包括线性或非线性向上扫描、向下扫描或者甚至伪随机扫描。伪随机扫描的例子例如在Goder,A,的“Minimum phase correlation of pseudo-random vibrator signals”(64thAnn.Internat.Mtg:Soc.of Expl.Geophys.,884-885,1994)和Cunningham,A.B.的“Some alternate vibrator signals”(Geophysics,V.44,pp1901-1921,1979)中说明。
允许信号分离的相位转动序列等同于添加着常量的相位转动序列。换言之,可以通过简单地对相位转动序列的每个相位项加上一个常量使该序列更加普通。例如,如果S1=0、144、288、72和216并且S2=0、216、72、288和144,则具有相同谐波消除和信号分离的等同序列是S3=10、154、298、82和226以及S4=10、226、82、298和154。应注意,当使用单个震动器源时,不应采用Edington规定的表中的第一列和最后一列。
上面的讨论和本发明的其中相关信号与用来驱动震动器的导频扫描信号关联的各实施例相关。在本发明的另一实施例中,相位转动扫描段序列和地面力(GF)信号而不是和导频或带有一重复段的导频交叉相关。GF信号代表震动器实际发送到地球中的能量,并且可以和导频基准扫描不同。通常,由于和地表面的非线性耦合,GF信号包含谐波和副谐波。可以通过授于Barr等的美国6161076号专利(收录为参考文献)中公开的方法或者任何其它方法测量地面力信号。GF信号测量是一种公用的实践,并且存在若干技术上周知的获取这些测量的方法。
例如,可以通过一种或多种方法从震动器上或者其附近的不同位置测量GF信号。在获取勘测数据期间可以把各种测量放置在辅助信道上。把GF测量保存到磁带或存储器允许例如在处理中心的进一步处理。可以通过技术上已知的各种方法调整GF信号,例如实质上均衡化各震动器或各震动器组之间的测量,或者组合各信号以建立地面力的加权的和表达。然后可以使用这些GF信号或它们的进一步处理的形式,供与现场记录相关或者供建立其它滤波器。
早就知道将地震记录和测到的地面力信号相关产生不希望有的基波和高阶谐波之间的交叉相关项。可以从把导频用作为相关运算子的图8和其中利用GF信号执行相关的图9中看出该效应。来自谐波的交叉相关项的谐波能从每个小弱波的中央向外扩散并且可能掩盖较弱的反射。
和导频信号相关的主要问题是,震动器输出的许多能量未呈现在得到的记录中。换言之,从震动器输出的谐波能未被积极利用。本发明利用来自震动器的谐波能。
图10表示利用GF信号的结果,由于把相关运算子应用到其中已转动各扫描从而抑制掉不希望的谐波能的无噪音模型上,该结果包含较高阶的谐波能。如果谐波能的振幅在各次扫描之间不变,会消除掉该谐波能;如果振幅不同,会衰减谐波能。
可以利用相关运算子(不论运算子是震动器扫描导频或是测到的GF)的振幅定标来补偿不同扫描期间随时间变化的震动器输出之间的差异。震动器输出是利用地面力信号测量的,从而可以建立定标运算子。例如,可以为一个震动器源建立用于平衡不同扫描间出现的振幅变化的定标运算子。对输出扫描的振幅的平衡或均衡将趋于均衡谐波能;不同扫描之间的谐波能均衡加大不希望的谐波能的衰减。
可以形成另一种振幅定标运算子,以便加强各源之间的信号分离而不只是对某特定源的扫描序列消除谐波。有效的和有用的源信号分离至少部分地取决于从同时运行的源发出的能量的振幅关系以及它们的相位。可以形成和用来在特定源的各扫描之间进行振幅补偿相类似的振幅补偿运算子,这可确保在同时运行的源之间恰当地出现所需的消除(并且从而实现信号分离)。
利用带有适当转动扫描的GF信号产生的结果优于利用导频扫描相关运算子产生的结果。当和图12的使GF信号与相位转动扫描的结果相关相比,如可从图11中看出,使用导频的信噪比较低。
把GF信号用作为相关运算子还可以应用于类似于上面概述的技术的相位转动级联扫描段方法。但是,一般情况下要求把二个扫描段添加到相关运算子上,一个在该运算子的开头处,另一个在相反的结束处。由N个段构成的级联扫描会变成一个带有N+2个段的级联GF相关运算子。例如,对于包含相位为0、144、288、72和216的震动器扫描段序列,应从测到的GF信号形成相关运算子。应复制震动器扫描GF信号的初始段和结束段并且添加到相反的二端上。从该示例的五段震动器扫描的输出记录的测到的GF信号应变成包括下述相位序列的七扫描段GF相关运算子的基础:216、0、144、288、72、216和0。图13中示出用于三个数据段A′、B′、C′(LISTEN)的相关处理,其中相关信号包括C、A、B、C、A。
应理解,不应把本发明严格地限制在为说明目的上述描述。对于本领域技术人员来说,在不背离附后权利要求中所规定的本发明的真正范围的情况下,各种修改和更替是明显的。

Claims (38)

1.一种利用通过扫描信号激励的震动地震能源的地震勘测方法,该方法包括:
(a)选择多个要用于勘测的地震能源,并且选择超过所述第一多个的第二数量的扫描段;
(b)选择要从利用所述地震源获取的处理数据中衰减的最高阶谐波的第三较高的阶数;
(c)定义和每个所述源关联的扫描信号,每个扫描信号包括多个扫描段,每个地震源的每个扫描段具有确定的相位;
(d)用其关联的扫描信号驱动每个所述源以产生传播到地球中的地震波并且记录响应的地震信号;以及
(e)利用多个与每个地震源相关的扩充信号处理所述记录的地震信号以产生多个处理过的地震信号,所述多个处理过的地震信号中的每一个和所述地震源中的一个不同的源关联并且使直到所述选定的第三较高阶数的各次谐波衰减。
2.如权利要求1的方法,还包括为每个所述扫描信号的每个所述段确定所述相位,其中所述相位按如下确定:
(i)构建一个N行和N列的表,其中N为该多个扫描段的数量,按以下公式计算每个项eH,M
eH,M=hm模N
其中h是行号,m是列号;
(ii)选择和源的数量相等的列数,以使每个选定列的第一行中的数在行H=F+1之前不在任何选定列中再次出现,其中F是确定为产生明显串馈的最高阶谐波的阶数;
(iii)对每个地震能源分配一个选定的列,该地震能源的第h个扫描的初始相位角用该列的第h行中的数乘以2B/N得到的数表示。
3.如权利要求1的方法,其中每个所述关联的扫描信号包括若干扫描段,所述扫描段顺序地相连或者重叠。
4.如权利要求3的方法,其中每个所述扩充的扫描信号还包括所述关联扫描信号以及(i)所述关联序列的第一扫描段和(ii)所述关联序列的最后段中的至少一个。
5.如权利要求1的方法,还包括测量生成的地震信号的电平,并且利用所述电平在相关所述记录地震数据之前调整所述其它的所述扩充信号中的每个的振幅。
6.如权利要求1的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波列。
7.如权利要求6的方法,其中所述频率的所述增大或减小随时间是线性变化的。
8.如权利要求1的方法,其中所述扫描段还包括伪随机扫描序列。
9.如权利要求8的方法,其中所述伪随机扫描序列按选定增量相位转动。
10.一种用于记录和处理震动源地震数据的方法,该方法包括:
(a)生成一个包括N个顺序地相连或重叠的扫描段的级联扫描序列,其中N大于或等于2,除了所述N个扫描段的初始相位角按大约为2mB/N弧度的固定相位增量累进转动外,所述N个扫描段基本相同,其中m为整数且|m|>1;
(b)利用所述级联扫描序列驱动一个震动源从而在一选定位置处把一个地震波传播到地球中;
(c)记录与所述地震波关联的地面力信号;
(d)记录至少一个来自地球内某位置处的对所述地震波做出反应的反射信号;以及
(e)利用所述记录的地面力信号处理记录的信号以产生处理信号。
11.如权利要求10的方法,其中处理记录信号包括使该记录信号与一个包括所述地面力信号以及所述地面力信号的第一个和最后一个扫描段的信号交叉相关。
12.如权利要求10的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波列。
13.如权利要求12的方法,其中所述频率的所述增大或减小随时间是线性变化的。
14.如权利要求10的方法,其中所述扫描段还包括伪随机的扫描序列。
15.一种利用n个由扫描信号驱动的震动地震能源的地震勘测方法,该方法包括:
(a)选择用于勘测的地震能源的数量;
(b)选择强度足以导致作为扫描序列的一部分的一个扫描段的明显谐波畸变的最高阶谐波的阶数;
(c)定义和每个所述源关联的基准信号,每个基准信号包括数量为N的扫描段,N大于n,所述扫描段是相连的或是重叠的;
(d)为每个基准信号的每个扫描段选择初始相位角,从而基本抑制掉直到并包括所述最高阶谐波的所有谐波,并且用各自的扫描序列驱动每个震动源和用各个关联基准信号驱动每个所述源,从而产生传播到地球中的地震波;
(e)记录和每个基准信号关联的地面力信号;
(f)记录至少一个来自地球内某位置处的对所述地震波做出反应的反射信号;
(g)利用所述n个地面力信号处理所述至少一个记录的反射信号。
16.如权利要求15的方法,其中处理所述至少一个记录信号包括使该至少一个记录信号和多个处理信号相关,每个所述处理信号包括一个基准信号组合该基准信号的至少二个段。
17.如权利要求15的方法,其中选择初始相位角还包括:
(I)构建一个N行和N列的表,从公式“ehm=hm模N”来计算每个项ehm,其中h为行号而m为列号;
(ii)选择和源的数量相等的列数,以使在行H=F+1之前使每个选定的列中的数都不在任何其它选定的列中再次出现;以及
(iii)对每个地震能源分配一个选定的列,该地震能源的第h个扫描段的初始相位角用该列的第h行中的数乘以2B/N得到的数表示。
18.如权利要求15的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波列。
19.如权利要求18的方法,其中所述频率的所述增大或减小随时间是线性变化的。
20.如权利要求15的方法,其中所述扫描段还包括伪随机扫描序列。
21.如权利要求20的方法,其中所述伪随机扫描序列按选定增量相位转动。
22.一种在大致相同的位置处同时进行n次地震勘测的方法,其中n为整数,该方法包括:
(a)同时把与地面力信号关联的n个震动源发射到地球中,每个所述地面力信号包括基波信号及其谐波,所述基波信号具有一个初始相位;
(b)记录所述n个地面力信号;
(c)记录至少一个来自地球内某位置对所述多个地面力信号作出反应的反射信号以给出第一记录信号;
(d)重复步骤(a)-(c)m-1次,其中按2B/m弧度偏移所述初始相位以给出总共m个记录信号和n×m个地面力信号;
(e)利用该n×m个地面力信号处理该m个记录信号以给出和每个震动源关联的处理信号。
23.如权利要求22的方法,其中处理该m个记录信号还包括:
(i)使第一记录信号和每个关联的地面力信号相关以给出n个中间信号;
(ii)重复步骤(i)m-1次以给出总共n×m个中间处理信号;
(iii)相加与每个震动源关联的m个中间处理信号以给出和每个震动源关联的处理信号。
24.如权利要求23的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波列。
25.如权利要求24的方法,其中所述频率的所述增大或减小随时间是线性变化的。
26.如权利要求22的方法,其中所述扫描段还包括伪随机扫描序列。
27.一种利用由扫描信号驱动的震动地震能源进行地震勘测的方法,该方法包括:
(a)选择用于勘测的地震能源的数量;
(b)选择强度足以导致作为扫描序列的一部分的一个扫描段的明显谐波畸变的最高阶谐波的阶数;
(c)选择超过所使用的地震能源的数量的扫描段的数量;
(d)为每个地震能源的每个扫描段选择初始相位角从而基本抑制掉直到并包括所述最高阶谐波的所有谐波,该选择初始相位角的步骤还包括:
(i)构建一个N行和N列的表,其中N为扫描段的数量,按公式“eH,M=hm模N”来计算每个项eH,M,其中h为行号而m为列号;
(ii)选择和源的数量相等的列数,以便在行H=F+1之前使每个选定列的第一行中的数不在任何选定列中再次出现,其中F是确定为产生明显串馈的最高阶谐波的阶数;
(iii)对每个地震能源分配一个选定的列,该地震能源的第h个扫描的初始相位角用该列的第h行中的数乘以2B/N得到的数表示;
(e)生成由具有在步骤(d)中得到的初始相位角的扫描段构成的并和每个能源关联的第一扫描序列,所述扫描段是顺序相连或重叠的;
(f)利用其关联的扫描序列驱动每个所述能源从而产生传播到地球中的地震信号;
(g)记录所述地震信号从地下地质结构层反射产生的未相关的地震数据;
(h)组合第一扫描段期间记录的所述未相关地震数据部分和“收听时间”期间记录的所述记录的未相关地震数据,以便形成组合的未相关地震数据记录;
(i)使所述组合的未相关地震数据记录和每个所述第一序列相关以给出和每个所述能源关联的相关数据。
28.如权利要求27的方法,还包括测量生成的地震信号的电平,并且利用所述电平在相关所述组合的记录地震数据之前调整所述其它的所述级联扫描序列中的每个序列的振幅。
29.如权利要求27的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波列。
30.如权利要求29的方法,其中所述频率的所述增大或减小随时间是线性变化的。
31.如权利要求27的方法,其中所述扫描段还包括伪随机扫描序列。
32.如权利要求31的方法,其中所述伪随机扫描序列按用权利要求27(d)的方法选择的增量进行相位转动。
33.一种利用由扫描信号驱动的震动地震能源进行地震勘测的方法,该方法包括:
(a)选择强度足以导致作为扫描序列的一部分的一个扫描段的明显谐波畸变的最高阶谐波;
(b)选择多个扫描段和每个扫描段的初始相位角,以便基本抑制掉直到并且包括所述最高阶谐波的所有谐波;
(c)产生和每个能源关联的第一扫描序列,所述扫描段为顺序相连或重叠的;
(d)利用该扫描序列驱动所述能源,从而产生传播到地球中的地震信号;
(e)记录所述地震信号从地下地质结构层反射产生的未相关地震数据;
(f)组合第一扫描段期间记录的所述未相关地震数据部分和“收听时间”期间记录的所述记录的未相关地震数据,以便形成组合的未相关地震数据记录;以及
(g)使所述组合的未相关地震数据记录和所述序列相关以给出相关数据。
34.如权利要求33的方法,其中选择所述初始相位角还包括:
(i)构建一个N行和N列的表,其中N为扫描段的数量,按公式
“eh,m=hm模N”
计算每个项eH,M,其中h为行号而m为列号;
(ii)选择和源的数量相等的列数,以便在行H=F+1之前使每个选定列的第一行中的数不在任何选定列中再次出现,其中F是确定为产生明显串馈的最高阶谐波的阶数;
(iii)对每个地震能源分配一个选定的列,该地震能源的第h个扫描的初始相位角用该列的第h行中的数乘以2B/N得到的数表示。
35.如权利要求33的方法,还包括测量生成的地震信号的电平,并且利用所述电平在相关所述组合的记录地震数据之前调整所述其它的所述级联扫描序列的每个序列的振幅。
36.如权利要求33的方法,其中所述扫描段还包括一个频率随时间单调增大或者随时间单调减小的正弦波群。
37.如权利要求36的方法,其中所述频率的所述增大或减小随时间是线性变化的。
38.如权利要求33的方法,其中所述扫描段还包括伪随机扫描序列。
CNB018190588A 2000-10-17 2001-10-17 对源编码和谐波消除采用级联扫描的方法 Expired - Fee Related CN1254693C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US24121100P 2000-10-17 2000-10-17
US60/241,211 2000-10-17

Publications (2)

Publication Number Publication Date
CN1488078A true CN1488078A (zh) 2004-04-07
CN1254693C CN1254693C (zh) 2006-05-03

Family

ID=22909719

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB018190588A Expired - Fee Related CN1254693C (zh) 2000-10-17 2001-10-17 对源编码和谐波消除采用级联扫描的方法

Country Status (6)

Country Link
US (2) US6687619B2 (zh)
EP (1) EP1327162A2 (zh)
CN (1) CN1254693C (zh)
AU (2) AU1176802A (zh)
CA (1) CA2426160A1 (zh)
WO (1) WO2002033442A2 (zh)

Families Citing this family (74)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2387226C (en) * 2002-04-06 2008-05-12 Westerngeco Ltd A method of seismic surveying
US6891776B2 (en) * 2002-09-04 2005-05-10 Westerngeco, L.L.C. Vibrator sweep shaping method
CA2520640C (en) * 2003-04-01 2012-10-23 Exxonmobil Upstream Research Company Shaped high frequency vibratory source
MXPA06001607A (es) * 2003-08-11 2006-05-19 Exxonmobil Upstream Res Co Metodo para exploracion y separacion continuas de multiples vibradores sismicos.
US20050128874A1 (en) * 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
GB2416033B (en) * 2004-07-10 2006-11-01 Westerngeco Ltd Seismic vibratory acquisition method and apparatus
US20070274155A1 (en) * 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
CA2664352C (en) * 2006-09-28 2011-09-27 Exxonmobil Upstream Research Company Iterative inversion of data from simultaneous geophysical sources
US8000168B2 (en) * 2006-12-08 2011-08-16 Conocophillips Company Dynamic source parameter selection for seismic vibrator data acquisition
DE102007013413A1 (de) * 2007-03-20 2008-10-09 GeoForschungsZentrum Potsdam Stiftung des öffentlichen Rechts Seismische Quelle mit adaptiver Regelung und entsprechendes Verfahren
WO2008123920A1 (en) 2007-04-10 2008-10-16 Exxonmobil Upstream Research Company Separation and noise removal for multiple vibratory source seismic data
US7859945B2 (en) * 2007-07-06 2010-12-28 Cggveritas Services Inc. Efficient seismic data acquisition with source separation
GB2451630B (en) * 2007-08-04 2009-12-09 Westerngeco Seismic Holdings Composite sweeps of high and low frequency part
US7869304B2 (en) * 2007-09-14 2011-01-11 Conocophillips Company Method and apparatus for pre-inversion noise attenuation of seismic data
US7639567B2 (en) * 2007-09-17 2009-12-29 Ion Geophysical Corporation Generating seismic vibrator signals
WO2009075961A1 (en) * 2007-12-12 2009-06-18 Exxonmobil Upstream Research Company Method and apparatus for evaluating submarine formations
BRPI0908156A2 (pt) 2008-03-21 2015-08-11 Exxonmobil Upstream Company Métodos para invesão de dados geofísicos medidos e para produzir hidrocarbonetos de uma região de subsuperfície
US8892410B2 (en) * 2008-08-11 2014-11-18 Exxonmobil Upstream Research Company Estimation of soil properties using waveforms of seismic surface waves
US8947976B2 (en) * 2008-10-03 2015-02-03 Westerngeco L.L.C. Harmonic attenuation using multiple sweep rates
FR2941055B1 (fr) 2009-01-12 2011-02-25 Cggveritas Services Sa Procede d'acquisition de donnees vibrosismiques concernant une zone du sous-sol, et procede d'exploration sismique incluant un tel procede
US8553496B2 (en) * 2010-02-09 2013-10-08 Ion Geophysical Corporation Seismic source separation
US8537638B2 (en) * 2010-02-10 2013-09-17 Exxonmobil Upstream Research Company Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration
US8400873B2 (en) 2010-02-22 2013-03-19 Westerngeco L.L.C. Vibroseis seismic acquisition technique
US8223587B2 (en) * 2010-03-29 2012-07-17 Exxonmobil Upstream Research Company Full wavefield inversion using time varying filters
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8756042B2 (en) 2010-05-19 2014-06-17 Exxonmobile Upstream Research Company Method and system for checkpointing during simulations
WO2011144215A2 (fr) * 2010-05-20 2011-11-24 Entreprise Nationale De Geophysique - Enageo- Filiale Du Groupe Sonatrach Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence
WO2012021290A1 (en) * 2010-08-10 2012-02-16 Conocophillips Company Method for creating an improved sweep for a seismic source
WO2012021293A1 (en) * 2010-08-11 2012-02-16 Conocophillips Company Unique seismic source encoding
US8767508B2 (en) 2010-08-18 2014-07-01 Exxonmobil Upstream Research Company Using seismic P and S arrivals to determine shallow velocity structure
EP2622457A4 (en) 2010-09-27 2018-02-21 Exxonmobil Upstream Research Company Simultaneous source encoding and source separation as a practical solution for full wavefield inversion
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
CN103238158B (zh) 2010-12-01 2016-08-17 埃克森美孚上游研究公司 利用互相关目标函数进行的海洋拖缆数据同时源反演
CN103703391B (zh) 2011-03-30 2017-05-17 埃克森美孚上游研究公司 使用频谱整形的全波场反演的系统和计算机实施的方法
CN103460074B (zh) 2011-03-31 2016-09-28 埃克森美孚上游研究公司 全波场反演中小波估计和多次波预测的方法
ES2640824T3 (es) 2011-09-02 2017-11-06 Exxonmobil Upstream Research Company Utilización de la proyección sobre conjuntos convexos para limitar la inversión del campo de onda completa
FR2981746B1 (fr) 2011-10-19 2014-11-21 Cggveritas Services Sa Source et procede d'acquisition sismique marine
FR2981758B1 (fr) 2011-10-19 2013-12-06 Cggveritas Services Sa .
FR2981759B1 (fr) 2011-10-19 2014-07-18 Cggveritas Services Sa Procede et dispositif pour determiner un signal de commande pour des sources marines vibrosismiques
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
US9348041B2 (en) 2012-02-15 2016-05-24 Westerngeco L.L.C. Phase modulation and noise minimization for simultaneous vibroseis acquisition
US9453928B2 (en) 2012-03-06 2016-09-27 Westerngeco L.L.C. Methods and computing systems for processing data
CA2861863A1 (en) 2012-03-08 2013-09-12 Exxonmobil Upstream Research Company Orthogonal source and receiver encoding
WO2013152174A1 (en) 2012-04-04 2013-10-10 Inova Ltd. Devices and systems for controlling harmonic distortion in seismic sources
US8619497B1 (en) 2012-11-15 2013-12-31 Cggveritas Services Sa Device and method for continuous data acquisition
US8724428B1 (en) 2012-11-15 2014-05-13 Cggveritas Services Sa Process for separating data recorded during a continuous data acquisition seismic survey
WO2014084945A1 (en) 2012-11-28 2014-06-05 Exxonmobil Upstream Resarch Company Reflection seismic data q tomography
CA2909105C (en) 2013-05-24 2018-08-28 Ke Wang Multi-parameter inversion through offset dependent elastic fwi
US10459117B2 (en) 2013-06-03 2019-10-29 Exxonmobil Upstream Research Company Extended subspace method for cross-talk mitigation in multi-parameter inversion
US9702998B2 (en) 2013-07-08 2017-07-11 Exxonmobil Upstream Research Company Full-wavefield inversion of primaries and multiples in marine environment
EP3033638B1 (en) 2013-08-12 2021-10-06 The University of Houston Low frequency seismic acquisition using a counter rotating eccentric mass vibrator
DK3036566T3 (en) 2013-08-23 2018-07-23 Exxonmobil Upstream Res Co SIMILAR SOURCE APPLICATION DURING BOTH SEISMIC COLLECTION AND SEISMIC INVERSION
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
US9910189B2 (en) 2014-04-09 2018-03-06 Exxonmobil Upstream Research Company Method for fast line search in frequency domain FWI
AU2015256626B2 (en) 2014-05-09 2017-10-12 Exxonmobil Upstream Research Company Efficient line search methods for multi-parameter full wavefield inversion
US10185046B2 (en) 2014-06-09 2019-01-22 Exxonmobil Upstream Research Company Method for temporal dispersion correction for seismic simulation, RTM and FWI
RU2016150545A (ru) 2014-06-17 2018-07-17 Эксонмобил Апстрим Рисерч Компани Быстрая вязкоакустическая и вязкоупругая инверсия полного волнового поля
US10838092B2 (en) 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10422899B2 (en) 2014-07-30 2019-09-24 Exxonmobil Upstream Research Company Harmonic encoding for FWI
US10386511B2 (en) 2014-10-03 2019-08-20 Exxonmobil Upstream Research Company Seismic survey design using full wavefield inversion
AU2015337108B2 (en) 2014-10-20 2018-03-01 Exxonmobil Upstream Research Company Velocity tomography using property scans
EP3234659A1 (en) 2014-12-18 2017-10-25 Exxonmobil Upstream Research Company Scalable scheduling of parallel iterative seismic jobs
US10520618B2 (en) 2015-02-04 2019-12-31 ExxohnMobil Upstream Research Company Poynting vector minimal reflection boundary conditions
CA2972028C (en) 2015-02-13 2019-08-13 Exxonmobil Upstream Research Company Efficient and stable absorbing boundary condition in finite-difference calculations
RU2017132164A (ru) 2015-02-17 2019-03-18 Эксонмобил Апстрим Рисерч Компани Многоэтапный процесс инверсии полного волнового поля, при выполнении которого образуют массив свободных от многократных волн данных
SG11201708665VA (en) 2015-06-04 2017-12-28 Exxonmobil Upstream Res Co Method for generating multiple free seismic images
US10838093B2 (en) 2015-07-02 2020-11-17 Exxonmobil Upstream Research Company Krylov-space-based quasi-newton preconditioner for full-wavefield inversion
CN108139499B (zh) 2015-10-02 2020-02-14 埃克森美孚上游研究公司 Q-补偿的全波场反演
WO2017065889A1 (en) 2015-10-15 2017-04-20 Exxonmobil Upstream Research Company Fwi model domain angle stacks with amplitude preservation
US10768324B2 (en) 2016-05-19 2020-09-08 Exxonmobil Upstream Research Company Method to predict pore pressure and seal integrity using full wavefield inversion
CN109991590B (zh) * 2019-02-21 2021-02-02 中国船舶重工集团公司第七一五研究所 一种在有限空间压力罐内测试换能器低频发射特性的系统与方法
US11681064B2 (en) * 2019-10-08 2023-06-20 Saudi Arabian Oil Company Look-ahead VSP workflow that uses a time and depth variant Q to reduce uncertainties in depth estimation ahead of a drilling bit
CN113127799B (zh) * 2020-01-15 2023-12-01 中国石油天然气集团有限公司 一种扫描信号校正方法及装置
US11573346B2 (en) 2021-04-15 2023-02-07 Saudi Arabian Oil Company Determining a seismic quality factor for subsurface formations for marine vertical seismic profiles

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1423366A (en) 1972-07-21 1976-02-04 Seiscom Ltd Broad line seismic profiling using simultaneously radiating sources
US4168485A (en) * 1974-08-12 1979-09-18 Continental Oil Company Simultaneous use of pseudo-random control signals in vibrational exploration methods
DE2448007C3 (de) * 1974-10-09 1978-11-16 Deutsche Texaco Ag, 2000 Hamburg Verfahren zum Erzeugen und Aufzeichnen von Seismogrammen, bei dem eine Serie von Vibratorsignalen je einiger Sekunden Dauer und monoton veränderlicher Frequenz in den zu untersuchenden Untergrund abgestrahlt wird und zwecks Erzeugung von oberwellenarmen Seismogrammen
US4295213A (en) * 1979-10-09 1981-10-13 Exxon Production Research Company Composite seismic signal
US4715020A (en) * 1986-10-29 1987-12-22 Western Atlas International, Inc. Simultaneous performance of multiple seismic vibratory surveys
US4751687A (en) * 1986-07-10 1988-06-14 Industrial Vehicles International, Inc. Methods of reducing harmonic distortion in continuous wave seismic exploration
US4823326A (en) * 1986-07-21 1989-04-18 The Standard Oil Company Seismic data acquisition technique having superposed signals
US4864546A (en) * 1987-03-23 1989-09-05 Terra Marine Engineering, Inc. Vibrator seismic data correlation system
US4982374A (en) * 1989-10-23 1991-01-01 Halliburton Geophysical Services, Inc. Method of source coding and harmonic cancellation for vibrational geophysical survey sources
US5901112A (en) * 1994-04-11 1999-05-04 Walker; David A. Signal energy enhancement for seismic exploration
US5410517A (en) * 1994-05-13 1995-04-25 Exxon Production Research Company Method for cascading sweeps for a seismic vibrator
US5721710A (en) * 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation

Also Published As

Publication number Publication date
US6865488B2 (en) 2005-03-08
US20040148104A1 (en) 2004-07-29
US6687619B2 (en) 2004-02-03
CA2426160A1 (en) 2002-04-25
EP1327162A2 (en) 2003-07-16
WO2002033442A3 (en) 2003-01-16
CN1254693C (zh) 2006-05-03
US20020091487A1 (en) 2002-07-11
AU1176802A (en) 2002-04-29
AU2002211768B2 (en) 2005-08-04
WO2002033442A2 (en) 2002-04-25

Similar Documents

Publication Publication Date Title
CN1254693C (zh) 对源编码和谐波消除采用级联扫描的方法
US6842701B2 (en) Method of noise removal for cascaded sweep data
RU2126983C1 (ru) Способ формирования сейсмических данных с использованием сейсмического вибратора (варианты)
CN1166963C (zh) 多个震动器单步反演分离方法
CN1133882C (zh) 源分离的高保真度震动源地震探矿法
EP0861448B1 (en) Method for separation of a plurality of vibratory seismic energy source signals
Bagaini Overview of simultaneous vibroseis acquisition methods
CA2534519C (en) Method for continuous sweeping and separation of multiple seismic vibrators
US20070195644A1 (en) Methods and Systems for Efficient Compaction Sweep
CN1073236C (zh) 处理具有多次反射噪声的地震数据的方法
AU2010213635A1 (en) Multiple seismic signal inversion
EA007911B1 (ru) Профилированный высокочастотный вибрационный источник
AU1397899A (en) Seismic data acquisition and processing using non-linear distortion in a groundforce signal
EA003494B1 (ru) Усовершенствованный способ сейсмической разведки
KR100412097B1 (ko) 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치
US4189704A (en) Method and apparatus for determination of optimum velocity in processing seismic signals from low energy sources
US5153857A (en) Method for selecting seismic traces for higher efficiency of pre-stack two dimensional or three dimensional depth migration
AU2005202928B2 (en) Method of using cascaded sweeps for source coding and harmonic cancellation
CN105319583A (zh) 基于分频动态编码的可控震源混叠数据成像方法
SU940096A1 (ru) Способ сейсмической разведки
Bagaini An Overview of Efficient Vibroseis Acquisition Methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20060503

Termination date: 20111017