Nothing Special   »   [go: up one dir, main page]

CN101620276A - 用于在地震数据中衰减多次反射的方法 - Google Patents

用于在地震数据中衰减多次反射的方法 Download PDF

Info

Publication number
CN101620276A
CN101620276A CN200910139598A CN200910139598A CN101620276A CN 101620276 A CN101620276 A CN 101620276A CN 200910139598 A CN200910139598 A CN 200910139598A CN 200910139598 A CN200910139598 A CN 200910139598A CN 101620276 A CN101620276 A CN 101620276A
Authority
CN
China
Prior art keywords
filtering
frequency
signal
seismic
reflection
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
CN200910139598A
Other languages
English (en)
Other versions
CN101620276B (zh
Inventor
E·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.)
PGS Geophysical AS
Original Assignee
PGS Geophysical AS
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 PGS Geophysical AS filed Critical PGS Geophysical AS
Publication of CN101620276A publication Critical patent/CN101620276A/zh
Application granted granted Critical
Publication of CN101620276B publication Critical patent/CN101620276B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808Seismic data acquisition, e.g. survey design
    • 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. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/20Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/21Frequency-domain filtering, e.g. band pass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/56De-ghosting; Reverberation compensation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

用于衰减地震信号中多次反射的一种方法包括相干滤波该地震信号。在选定的频率低通滤波和高通滤波该地震信号,使得多次反射实质上仅仅存在于该频率以上的地震信号中。产生具有单位值的掩码,其中该高通滤波的信号幅度超出所选定的阈值。该掩码对于所有其他的滤过的信号具有单位零值。将所述低通滤波信号与该掩码值相乘。将所述相乘的信号加入高通滤波信号并且应用相干滤波器以除去剩余的初始能量以产生所述多次反射的模型。从该相干性滤波的地震信号中减去该多次反射模型。

Description

用于在地震数据中衰减多次反射的方法
技术领域
本发明通常涉及地震勘探领域。更具体地说,本发明涉及用于处理地震数据以衰减多次反射效应的方法。
背景技术
在地震勘探中,为获得地震数据,通过将声能注入到地球表面附近并检测从表面下的地层结构的不同层之间边界反射的声能。当相邻层间对边界的声阻抗不同时,声能被反射。解析表示检测到的声能的信号来推断表面下地球结构的构造和组成。
在海洋地震勘探中,一般采用诸如气枪或者气枪阵列这样的地震能量源将声能注入地球。当气枪或阵列通过船拖曳时,气枪或阵列一般在水中选定的深度启动。相同的或者不同的船在水中拖拽一个或多个被称作“横幅”的地震传感器缆索。通常,横幅在船后沿着其被拖曳的方向延伸。典型地,横幅包括多个设置在缆索上的水听器(位于沿着缆索的相隔一定距离的已知位置)。本领域已知的水听器是一种检测器,该检测器产生对应于水压或者水中压力的时间梯度的光信号或者电信号。拖拽一个或多个横幅的船通常包括记录装置以记录(以时间为索引)水听器产生的信号,从而响应于检测到的声能。如上所述,处理信号的记录以推断位于进行地震测量的位置以下的地层构造和组成。
海洋地震数据常常包括两种特殊的假象,它们需要技术来解析以便能更精确地推断表面下地层的构造和组成。这二种被称为虚反射(ghosting)和水层多次反射(water layer multiplereflection)的假象的出现是由于水与水面上的空气相比具有基本上不同的声阻抗以及水与水底(或者海底)的地层结构相比通常具有基本上不同的声阻抗。
虚反射和水层多次反射可以理解如下。当启动气枪或者气枪阵列时,声能通常向下辐射,它穿过海底并进入表面下地层。如前所述,一些声能在地层之间的表面下声阻抗边界被反射。反射的声能通常向上地传播最终被一个或多个横幅上的地震传感器(水听器)所检测。然而在该反射能量到达横幅之后,它继续向上传播直到到达水面。水面相对于向上传播的声能具有几乎全反射率(反射系数等于1)。因此,几乎所有向上传播的声能将会从水面反射并再一次向下传播。从水面反射的声能相对于向上传播的入射的声能还将相移大约180度。表面反射的向下传播的声能通常被称为“虚(ghost)”信号。在水听器检测到的声能中,虚信号在受限频率内造成能量的明显的“切口”或者衰减。本领域众所周知,在检测到的声信号中切口的频率与设置横幅的选定深度有关。
从水面反射的向下传播的声能以及直接从地震能量源散发的声能可能从水底反射并向上传播,在那里被水听器检测。该同一向上传播的声能也将会从水面反射并再一次向下传播。因此,声能在被衰减以前可能在水面和水底两者之间多次反射,导致所谓的水层回响(reverberations)。这种回响在整个检测到的声能内可具有相当的幅度,从而掩盖了从表面下地层边界反射的声能,并因此使得利用地震数据推断地下结构和组成变得更困难。
有许多“去虚反射”和水层效应衰减的技术。授予Ruehle的美国专利号4,486,865描述了一种这样的技术。多对检测器的每对包含地震检波器(geophone)和水听器。滤波器被用于各对中的地震检波器或者水听器的至少一个的输出端以便调整该滤波的信号的频率组成。调整该频率组成使得当该滤波的信号与来自其它检测器的信号结合时,消除虚反射。
授予Moldovenu的美国专利号5621700也披露了在衰减虚反射和水层回响的方法中使用至少一对检测器。
授予Sanders等人的美国专利号4935903披露一种降低水层回响效应的方法,该方法包括在垂直间隔距离的深度上测量压力,或者通过使用检测器对来测量压力和质点运动。该方法包括通过加入虚反射数据来增强用于叠前处理的原始反射数据。
美国专利号4979150披露一种用于海洋地震勘探的方法,其中基本上并列的水听器和地震检波器的输出经过比例系数(scalefacor)的处理。该并列的水听器和地震检波器可以被置于海底或者海底之上。
水体以下的许多地表下部分使用水底缆索来勘测是不切实际的,此外,目前可行的海洋地震采集技术使用了水听器检测器。更进一步,存在大量这样的水听器海洋地震数据,它们有益于改进用于分离多次反射的技术。因此,仍然需要技术来衰减水层多次反射对地震数据的效应。
发明内容
根据本发明的一方面,一种用于衰减地震信号中的多次反射的方法包括相干滤波(coherency filtering)该地震信号。地震信号为一定频率处的低通滤波和高通滤波地震信号,该频率的选择使得在该频率之上的地震信号中基本上只存在多次反射。产生具有单位值(unity value)的掩码(mask),这里高通滤波的信号幅度超出所选定的阈值。该掩码对于所有其他的滤波信号具有单位值。将低通滤波信号与掩码值相乘。将相乘的信号加入高通滤波信号以产生多次反射的模型。从所述相干滤波的地震信号中减去该多次反射模型。
本发明的另一个方面是存储在计算机可读介质中的计算机程序。该程序包括可操作的逻辑,使得可编程计算机执行特定的动作,包括在选定的频率上低通滤波和高通滤波地震信号。选定频率使得该频率上的地震信号中,基本上只存在多次反射。产生具有单位值的掩码,这里高通滤波信号的幅度超出所选定的阈值。该掩码对于所有其他滤波信号具有单一值。将低通滤波信号与掩码值相乘。将相乘的信号加入高通滤波信号以产生多次反射的模型。从相干滤波的地震信号中减去多次反射模型。
本发明的其他的方面和优点根据以下说明和附上的权利要求将会变得明显。
附图说明
图1示出包括水层回响的海洋地震数据采集的一个例子。
图2是根据本发明方法的示例性实现方式的流程图。
图3示出具有可读存储介质的通用可编程计算机。
具体实施方式
图1示意性地示出本发明可用的海洋地震数据采集的例子。地震勘测船10沿着诸如湖或者海洋的水体18的表面18A移动。勘测船10包括以12显示的并统称为“记录系统”的设备,该设备包括装置(未单独显示),用于在选定时间启动地震能量源14,用于确定任何时候的船10、源14和多个地震传感器16中每一个的测地学位置(geodetic position),以及用于记录震传感器16产生的信号。记录系统12可能包括的各种装置在本领域为大家所熟知并且无需详述。
地震传感器16通常为水听器或者其他压力或压力梯度响应检测器。在其他的例子中可采用其他类型的地震传感器,例如质点运动响应检测器。地震传感器16沿着缆索(称为横幅)以定距离间隔的位置布置。为了例示的清楚,图1中的例子仅仅示出一个横幅,但是横幅的数目、横幅上检测器的数目以及有一个或多个横幅是否由勘探船10或者不同的船拖曳并不限于本发明的范围。
在水18的表面18A以下的选定深度拖曳地震能量源14和地震传感器16。在地震数据采集期间,源14(例如可以是气枪或气枪阵列)在选定时间启动并且地震传感器16响应由源14发射的地震能量而产生的信号由记录系统12中适当的装置记录。信号记录一般以地震能量源14的启动时间为索引。可以为每个地震传感器16做单独的信号记录或者选定的地震传感器16组可进行信号通信,使得由所有在该组中的检测器产生的信号组成对应的记录。通常,将有数千个单独的信号记录(被称为“信道”),对应于在以一定距离间隔的地震传感器组的每一个处所检测到的信号。
通常,在启动源14后,地震能量向外传播并且通过水18向下传播(如标记22所示),直到到达水底18B。地震能量继续向下通过地表直到到达一个或多个表面下声阻抗边界20。能量可在这些边界反射并且向上传播(如标记24所示),在那里其到达地震传感器16并被检测。该地震能量继续向上传播直到它到达水面18A。然后地震能量从水面18A反射并再一次向下传播(如标记26所示),在那里其可被地震传感器166检测并继续向下传播。该向下传播的地震能量,不管是从水面18A反射或者是从源14直接传播的,也可能从水底18B反射。水面18A和水底18B之间的反射可能延续若干次,在由检测器16检测的地震信号中产生假象,称作“水层回响”或“水层多次反射(water layermultiples)”。本发明的目的在于减少在地震信号中的“水层多次反射”的效应。
本发明的方法基于“初始反射”的带宽,即在表面下边界一次反射的地震能量(如图1中的24处所示),其受制于较低的频率(例如小于大约40Hz),而多次反射的带宽较宽。
参考图2,在31处,根据本发明方法的第一处理单元是将原始记录的地震信号分类成为公共迹属性集合(common raceattribute gathers)。与地震信号处理相关使用的“迹”表示一种记录,通常是关于地震信号幅度的记录时间,该地震信号由在接收器阵列中的特定地震接收器(或者这样的接收器组)检测。属性应当选择为使得可以将相干滤波应用到该集合迹上以增强所关心的地震信号和减少随机噪声的效应。一个公共属性集合的例子是公共的源到接收器偏移量集合。其他的例子包括公共的中点集合和公共的发射点集合。
如上所述,信号记录通常包括在源启动时多个单独的地震传感器的每一个或者这类检测器组的每一个在源启动时所做的一个信号记录。这种公共偏移量集合是地震信号记录的集合,其中当记录地震信号的时候,每个选定的单独信号记录表示沿着水面(图1中18A)源的位置和地震传感器之间的相同的距离。初级反射和多次反射在这种公共偏移量集合中作为相干事件出现,即信号中出现的高幅度事件将在一定程度上显现该高幅度事件的时间和该采集几何学(偏移量)之间的关系。
在32,一个被称作相干滤波的过程(一个例子是FX-反褶积)在公共偏移量集合上被执行以衰减随机噪声,即衰减在那些看起来与采集几何不具有任何关系的信号中的事件。在FX-反褶积中,输入信号的空间窗口(例如10个单独的传感器信号或者“迹”)和时间窗口(例如20毫秒)被傅里叶变换为频率位置(F-X)域。在该位置(X)维度中设计反褶积算子以预测该信号的相干部分。减去该相干部分将会留下不相干部分(随机噪声),然后该不相干的部分可以被反向变换并且从该原始信号中减去。然后对后继的窗进行相似的处理。较好地,在窗口之间有时间和/或空间的重叠。该相干滤波的结果将会在以下说明的过程中使用。
在33,将低通滤波器应用于地震信号记录。选定低通滤波器的截止频率为预计包含在原始反射中的最高频率。低通滤波器输出的信号包含基本上整个原始反射(例如图1中的24)和多次反射的低频部分(例如图1中的26)。该低通滤波信号可表示为Mlow。该低通滤波器可限定为具有0-30Hz的通带,从30Hz开始滑衰(roll off)并且40Hz的频率以及以上被强烈地衰减。根据该特定地震数据的频率组成,可为该低通滤波器的其他截止频率选定相似的滤波器特性曲线。
在34,将高通滤波器应用于信号记录(在32)。信号记录的高通滤波的输出基本上不包含原始反射,但是将包含多次反射的高频部分。高通滤波的记录可表示为Mhigh。该高通滤波器可限定为强烈地衰减0到30Hz之间的频率,从40Hz开始滑衰以及使40Hz以及更高的频率自由地通过。在33和34分别使用的低频滤波器和高通滤波器的过渡带(滑衰)应该选定为使得当它们的幅度频谱被求和时,对于所有频率,最终的幅度频谱等于单位值(unity)。这两种滤波器也应该为零相位滤波器,免得改变它们各自输入信号的相位。
在37,将在上文高通滤波之后的迹样本的大小(绝对值)与表示为阈值1的阈值相比。如果迹样本的大小大于阈值1,那么设置对应的迹输出样本等于单位值。如果该迹的大小小于或等于阈值1,那么设置迹的输出样本等于零。上述程序产生一种图(或者“掩码”),该图的样本值或者为1(表示多次子波的样本);或者为零(表示原始反射子波的样本或者表示代表无多次反射子波的样本)。因为包含在公共迹属性集合的迹中的样本值可以随每次地震勘测有很大的不同,所以可以计算阈值1的数值,例如计算为在每个迹集合中所有迹样本的平均绝对值。在实际的例子中,可以通过解释过程来选定阈值1的值。例如,可以对一个公共偏移量集合检查时间窗中的数据值的平均幅度[平均绝对值](作为初始值)和前述处理的结果,该数据值包含了任何多重反射刚刚到达之前(即,刚刚在两倍于海洋底部反射时间之前)的反射。上述步骤可以利用阈值1的逐步减小和逐步升高的值进行重复。视觉检查可能导致阈值1的数值的选择看起来对于整个数据集的处理来说是最优化的。如果阈值1过高,该幅度图只由零组成,并且所有低频率数据从该模型中被消除。如果阈值1设置太低,则幅度图全由1组成,并且所有低频数据被包含在该模型数据中,该模型数据会从原始数据自适应地减去(以下说明),因此衰减或者除去了该过程输出的原始反射。
在38,该图被逐个样本地与来自33的低通滤波信号的输出相乘。结果应该基本上除去了表示原始反射子波的样本并且仅仅留下多次反射的低频部份Mlow。注意,如果多次反射与原始反射一致,那么最终的地震子波将会等于两种子波的和。所以38的输出可以表示为Mlow+Poverlaps
在39,将来自38的相乘输出与来自32的高通滤波输出相加,产生全带宽多次反射Mlow+Mhigh(加上任何Poverlaps)。该Poverlaps子波的幅度通过相干滤波可以实质性地降低,如在40所示。可以采用的相干滤波的一个例子是FX-反褶积。假定Poverlaps幅度通过先前应用的掩码分段,那么该Poverlaps幅度应该在原始反射幅度存在的频率范围内,相对于相干滤波器(例如FX-反褶积)随机出现。任何原始反射能量将因此被衰减,从而将它保存在如下所述的自适应减去步骤中。
在41,该过程包括将幅度斜度(ramp)应用到每个迹以便当没有多次反射(即,当时间值小于两倍的水-海底反射时间时)不会执行适应性减去(以下说明)。斜度值可以从时间等于零(源启动时间)到t(=两次水底反射时间减去一个预选的滞后时间)期间取值为零。斜度值在两倍于水-底反射时间时可以线性增加到1.0。斜度值可以保持为1.0直到每个迹的结束。将斜度逐个样本地与输入的迹幅值相乘。在41执行的斜度输出是一个表示包含在原始的迹集合中的多次反射子波的模型的迹集合。
在42,从在32完成的FX-反褶积(相干滤波)自适应地减去来自41的上述多次反射的模型。该自适应减去的输出可以在43通过相干滤波器(例如FX-反褶积)处理。在45将这种滤波的输出与在43中对相干滤波的输入的选定片断(经过标度化)相加。由于在某些情况下可能不希望将FX-反褶积应用于短时间(浅的)数据,所以“标度”只限于应用F-X反褶积的部分。在45,该过程的输出为多次反射基本衰减的迹。
本发明的另一个方面涉及存储在计算机可读介质中的计算机程序。参考图3,根据图2说明的上述步骤可以以存储在计算机可读媒介(例如软盘64、只读光盘存储器60或者构成通用可编程计算机50的一部分的磁性硬盘驱动器66)上的计算机可读编码的方式实现。本领域已知的计算机50可以包括中央处理单元60、诸如键盘68之类的用户输入设备和诸如平板LCD显示器或者阴极射线管显示器之类的用户显示器52。计算机可以该包括适当的装置(分别如58和56所示)用于查询该软盘和只读光盘存储器。根据本发明的这方面,计算机可读介质包括在操作时可使得计算机50执行上述提出和根据图2说明的动作的逻辑。
本发明可以提供改进的从海洋地震数据中分离原始反射与多次反射的能力。
本发明已经根据有限的实施例进行说明,受益于本发明的本领域技术人员,将理解可以设计不离开在这里揭露的本发明范围的其他的实施例。因此,本发明的范围应该仅仅由附上的权利要求限定。

Claims (18)

1.一种用于从海洋地震信号中分离多次反射的方法,包括:
相干滤波该地震信号;
在选定的频率低通滤波该地震信号;
在该选定的频率高通滤波该地震信号,该频率的选择使得该频率之上的地震信号中基本上只存在多次反射;
产生具有单位值的掩码,这里该高通滤波信号的幅度超出所选定的阈值,该掩码对于所有其他的滤波信号具有单位值;
将该低通滤波信号与该掩码值相乘;
将该相乘的信号加入高通滤波信号以产生多次反射的模型;
从所述相干滤波的地震信号中减去该多次反射模型;以及
对减去的结果至少进行存储和显示中的一项。
2.如权利要求1所述的方法,进一步包括在相干滤波、低通滤波和高通滤波前将地震信号分类成公共迹属性集合。
3.如权利要求1所述的方法,其特征在于该相干滤波包括FX-反褶积。
4.如权利要求1所述的方法,进一步包括相干滤波所减去的结果。
5.如权利要求4所述的方法,其特征在于相干滤波所减去的结果包括FX-反褶积。
6.如权利要求1所述的方法,其特征在于所述减去包括适应性减去。
7.如权利要求1所述的方法,其特征在于所述频率大约为40Hz。
8.如权利要求1所述的方法,其特征在于所述低通滤波的滑衰部份和所述高通滤波的滑衰部分的联合幅度频谱基本上是单位的。
9.如权利要求1所述的方法,其特征在于所述低通滤波和所述高通滤波基本上是零相位的。
10.一种存储在计算机可读介质内的计算机程序,所述程序具有在操作时可以使得可编程计算机完成动作的逻辑,包括:
相干滤波该地震信号;
在选定的频率低通滤波该地震信号;
在该选定的频率高通滤波该地震信号,该频率的选择使得该频率之上的地震信号中基本上只存在多次反射;
产生具有单位值的掩码,这里该高通滤波信号的幅度超出所选定的阈值,该掩码对于所有其他的滤波信号具有单位值;
将该低通滤波信号与该掩码值相乘;
将该相乘的信号加入高通滤波信号以产生多次反射的模型;
从所述相干滤波的地震信号中减去该多次反射模型;以及
对减去的结果至少进行存储和显示中的一项。
11.如权利要求10所述的计算机程序,进一步包括在相干滤波、低通滤波和高通滤波前将地震信号分类成公共迹属性集合。
12.如权利要求10所述的计算机程序,其特征在于该相干滤波包括FX-反褶积。
13.如权利要求10所述的计算机程序,进一步包括相干滤波所减去的结果。
14.如权利要求13所述的计算机程序,其特征在于相干滤波所减去的结果包括FX-反褶积。
15.如权利要求10所述的计算机程序,其特征在于所述减去包括适应性减去。
16.如权利要求10所述的计算机程序,其特征在于所述频率大约为40Hz。
17.如权利要求10所述的计算机程序,其特征在于所述低通滤波的滑衰部份和所述高通滤波的滑衰部分的联合幅度频谱基本上是单位的。
18.如权利要求10所述的计算机程序,其特征在于所述低通滤波和所述高通滤波基本上是零相位的。
CN200910139598.1A 2008-06-30 2009-06-30 用于在地震数据中衰减多次反射的方法 Expired - Fee Related CN101620276B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US12/215862 2008-06-30
US12/215,862 2008-06-30
US12/215,862 US7675812B2 (en) 2008-06-30 2008-06-30 Method for attenuation of multiple reflections in seismic data

Publications (2)

Publication Number Publication Date
CN101620276A true CN101620276A (zh) 2010-01-06
CN101620276B CN101620276B (zh) 2015-03-25

Family

ID=41172445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910139598.1A Expired - Fee Related CN101620276B (zh) 2008-06-30 2009-06-30 用于在地震数据中衰减多次反射的方法

Country Status (12)

Country Link
US (1) US7675812B2 (zh)
EP (1) EP2141514B1 (zh)
CN (1) CN101620276B (zh)
AT (1) ATE555401T1 (zh)
AU (1) AU2009202480B2 (zh)
BR (1) BRPI0902020B1 (zh)
CA (1) CA2668376C (zh)
EA (1) EA017331B1 (zh)
EG (1) EG25387A (zh)
MX (1) MX2009007165A (zh)
MY (1) MY147741A (zh)
SG (1) SG158010A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103582825A (zh) * 2011-05-13 2014-02-12 沙特阿拉伯石油公司 同步源地震数据的变频滤波
CN104395779A (zh) * 2012-07-10 2015-03-04 雪佛龙美国公司 用于估计和减弱地震数据中的噪声的系统和方法
CN108919357A (zh) * 2018-05-16 2018-11-30 中国海洋石油集团有限公司 一种基于频谱重构的鬼波压制方法

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100274492A1 (en) * 2009-04-24 2010-10-28 Susanne Rentsch Determining attributes of seismic events
US8300498B2 (en) * 2009-06-30 2012-10-30 Pgs Geophysical As Method for dynamic aperture determination for three-dimensional surface-related multiple elimination
US8818730B2 (en) 2010-07-19 2014-08-26 Conocophillips Company Unique composite relatively adjusted pulse
WO2012074612A1 (en) 2010-12-01 2012-06-07 Exxonmobil Upstream Research Company Primary estimation on obc data and deep tow streamer data
SG193233A1 (en) * 2011-03-31 2013-10-30 Exxonmobil Upstream Res Co Method of wavelet estimation and multiple prediction in full wavefield inversion
AU2012290435B2 (en) * 2011-07-29 2014-11-20 Shell Internationale Research Maatschappij B.V. Method for increasing broadside sensitivity in seismic sensing system
US9075162B2 (en) 2011-11-10 2015-07-07 Pgs Geophysical As Method and system for separating seismic sources in marine simultaneous shooting acquisition
US20150032379A1 (en) * 2013-03-13 2015-01-29 Westerngeco L.L.C. Attenuation of multiple reflections
US9784869B2 (en) 2013-05-15 2017-10-10 Pgs Geophysical As Noise models by selection of transform coefficients
US9817143B2 (en) 2013-10-30 2017-11-14 Pgs Geophysical As Methods and systems for constraining multiples attenuation in seismic data
US11391857B2 (en) * 2013-12-30 2022-07-19 Pgs Geophysical As Methods and systems for attenuating residual acoustic energy in seismic data
US10670757B2 (en) * 2014-02-26 2020-06-02 Pgs Geophysical As Methods and systems for quantifying coherency and constraining coherency-based separation in simultaneous shooting acquisition
US9964656B2 (en) 2014-08-29 2018-05-08 Pgs Geophysical As Methods and systems to remove particle-motion-sensor noise from vertical-velocity data
US10444386B2 (en) 2014-08-29 2019-10-15 Pgs Geophysical As Methods and systems that determine a velocity wavefield from a measured pressure wavefield
WO2018071628A1 (en) * 2016-10-13 2018-04-19 Downunder Geosolutions (America) Llc Method for the attenuation of multiple refelections in shallow water settings
CN108732624A (zh) * 2018-05-29 2018-11-02 吉林大学 一种基于pca-emd的并行震源地震数据随机噪声压制方法
NO20210922A1 (en) 2019-02-21 2021-07-21 Halliburton Energy Services Inc Waveform processing utilizing an amplitude adaptive data mask
CN111929726B (zh) * 2019-09-04 2023-02-10 中国石油天然气股份有限公司 地震相干数据体处理方法及装置
GB2590177B (en) * 2019-11-21 2023-11-01 Cgg Services Sas Methods and devices performing adaptive subtraction of multiples

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1195204A (en) * 1966-06-15 1970-06-17 Seismograph Service England Methods for Use in Identifying Multiple Reflections.
US4486865A (en) 1980-09-02 1984-12-04 Mobil Oil Corporation Pressure and velocity detectors for seismic exploration
US4935903A (en) 1989-05-30 1990-06-19 Halliburton Geophysical Services, Inc. Reinforcement of surface seismic wavefields
EG19158A (en) 1989-08-25 1996-02-29 Halliburton Geophys Service System for attenuation of water-column reverberation
US5010526A (en) * 1990-07-30 1991-04-23 Mobil Oil Corporation Method of digital filtering
GB2273359B (en) * 1992-12-12 1997-01-15 Schlumberger Ltd Method for improving signal to noise ratio
US5572483A (en) * 1995-07-25 1996-11-05 Western Atlas International, Inc. Method of reducing noise in seismic signals by adaptive filtering of a noise reference
US5703833A (en) * 1995-11-13 1997-12-30 Mobil Oil Corporation One step inversion/separation scheme using a plurality of vibrator sources
US5621700A (en) 1996-05-20 1997-04-15 Schlumberger Technology Corporation, Geco-Prakla Div. Method for attenuation of reverberations using a pressure-velocity bottom cable
US5818795A (en) * 1996-10-30 1998-10-06 Pgs Tensor, Inc. Method of reduction of noise from seismic data traces
US6094620A (en) * 1999-01-08 2000-07-25 Exxonmobil Upstream Research Company Method for identifying and removing multiples from seismic reflection data
US7415401B2 (en) * 2000-08-31 2008-08-19 Exxonmobil Upstream Research Company Method for constructing 3-D geologic models by combining multiple frequency passbands
US6480440B2 (en) * 2001-03-07 2002-11-12 Westerngeco, L.L.C. Seismic receiver motion compensation
US7239578B2 (en) * 2005-03-03 2007-07-03 John M. Robinson Removal of noise from seismic data using radon transformations
BRPI0614103A2 (pt) * 2005-07-28 2012-11-20 Exxonmobil Upstream Res Co mÉtodo implementado por computador para reduzir o ruÍdo dependente de desvio de alta freqÜÊncia misturado com um sinal real em um sinal registrado por um receptor em uma prospecÇço eletromagnÉtica de fonte controlada em uma regiço subterrÂnea ao largo da costa e mÉtodo para produzir hidrocarbonetos a partir de uma regiço de sub-superfÍcie

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103582825A (zh) * 2011-05-13 2014-02-12 沙特阿拉伯石油公司 同步源地震数据的变频滤波
CN103582825B (zh) * 2011-05-13 2016-12-14 沙特阿拉伯石油公司 同步源地震数据的变频滤波
CN104395779A (zh) * 2012-07-10 2015-03-04 雪佛龙美国公司 用于估计和减弱地震数据中的噪声的系统和方法
US9612351B2 (en) 2012-07-10 2017-04-04 Chevron U.S.A. Inc. System and method for estimating and attenuating noise in seismic data
CN104395779B (zh) * 2012-07-10 2017-11-07 雪佛龙美国公司 用于估计和减弱地震数据中的噪声的系统和方法
CN108919357A (zh) * 2018-05-16 2018-11-30 中国海洋石油集团有限公司 一种基于频谱重构的鬼波压制方法
CN108919357B (zh) * 2018-05-16 2019-10-11 中国海洋石油集团有限公司 一种基于频谱重构的鬼波压制方法

Also Published As

Publication number Publication date
EG25387A (en) 2011-12-25
US7675812B2 (en) 2010-03-09
AU2009202480A1 (en) 2010-01-14
AU2009202480B2 (en) 2014-03-27
CA2668376C (en) 2014-09-09
CA2668376A1 (en) 2009-12-30
EA017331B1 (ru) 2012-11-30
MX2009007165A (es) 2010-02-17
BRPI0902020B1 (pt) 2019-10-15
BRPI0902020A2 (pt) 2010-04-13
EP2141514A2 (en) 2010-01-06
SG158010A1 (en) 2010-01-29
EP2141514A3 (en) 2011-02-23
EA200900650A1 (ru) 2010-04-30
CN101620276B (zh) 2015-03-25
MY147741A (en) 2013-01-15
ATE555401T1 (de) 2012-05-15
EP2141514B1 (en) 2012-04-25
US20090323470A1 (en) 2009-12-31

Similar Documents

Publication Publication Date Title
CN101620276B (zh) 用于在地震数据中衰减多次反射的方法
EP3339910B1 (en) Device and method for model-based deblending
US6738715B2 (en) Method for attenuating noise in seismic data
US6775618B1 (en) Method and system for reducing effects of sea surface ghost contamination in seismic data
CN103245969B (zh) 用于在震源虚反射去除之后确定源特征波形的方法和系统
CN101014881B (zh) 用于压力传感器和质点运动传感器记录的地震数据中水底多次波的衰减的系统
AU2018220021B2 (en) Systems and methods for removing acquisition related effects from seismic data
EP2787370B1 (en) Device and method for de-blending simultaneous shot data
US20050013194A1 (en) Method for seismic exploration utilizing motion sensor and pressure sensor data
US20060250890A1 (en) Method for deghosting and water layer multiple reflection attenuation in marine seismic data
NO20161660A1 (en) Marine surveys conducted with multiple source arrays
WO2014195508A2 (en) Systems and methods for de-noising seismic data
US9448315B2 (en) Device and method for denoising ocean bottom data
CN102103217A (zh) 用于衰减双传感器地震数据中的干扰噪声的方法
CN102124375A (zh) 使用双传感器记录系统衰减地震干扰噪声
Zhou et al. Multiple suppression by 2D filtering in the parabolic τ–p domain: a wave‐equation‐based method1
CN106199699B (zh) 利用传输矩阵法清除鬼波的方法
US8417458B2 (en) Removing ground-roll from geophysical data
US11391857B2 (en) Methods and systems for attenuating residual acoustic energy in seismic data
US20190113644A1 (en) Marine Seismic Use of a Harmonic Distorted Signal
Wang et al. The 3D converted wave seismic data processing in Lamadian area of Daqing oilfield

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150325

Termination date: 20210630