CN1073236C - 处理具有多次反射噪声的地震数据的方法 - Google Patents
处理具有多次反射噪声的地震数据的方法 Download PDFInfo
- Publication number
- CN1073236C CN1073236C CN95116803A CN95116803A CN1073236C CN 1073236 C CN1073236 C CN 1073236C CN 95116803 A CN95116803 A CN 95116803A CN 95116803 A CN95116803 A CN 95116803A CN 1073236 C CN1073236 C CN 1073236C
- Authority
- CN
- China
- Prior art keywords
- repeatedly
- reflection wave
- ripple
- extremely
- wave
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000012545 processing Methods 0.000 title abstract description 14
- 238000012937 correction Methods 0.000 claims description 81
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 238000013480 data collection Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 101100400452 Caenorhabditis elegans map-2 gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/368—Inverse filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
处理包含多次反射波地震数据的方法,对推测的多次反射波生成波至确定双向旅行时,对记录的地震记录线进行正常时差校正操作。将正常时差校正应用于该数据。一旦正常时差校正操作和双向旅行时移动已被用于记录线,反向正常时差校正操作利用扩充速度函数而被用于记录线。结果记录线现在只是包括多次反射波的记录线。最后,多次反射波噪声可以利用定标因数和剩余时间移动从原始数据中扣除。从而克服了现有技术的限制,并实现较小的计算消耗。
Description
本发明涉及处理地震数据的一种方法。另一方面,本发明涉及处理具有多次反射噪声的地震数据的方法,还有一方面,本发明涉及在地震数据中隔离多次反向噪声的方法。
地震测定通常包含领先一个或多个地震源在地表面所形成的地震脉冲。由于不同地下结构层的界面上具有不同的弹性,这些地震脉冲被反向和/或折射一小部分传送到地下。
检测器如地震仪,小型地震仪或水中地震检波器根据检测的地震波的反射和/或折射产生地震的模拟电信号或地震记录线信号。来自检测器的地震模拟电信号或地震记录线信号被记录。或者,来自检测器的地震模拟电信号或地震记录线信号在记录之前被采样和数字化,对两种方式记录的地震数据依次被处理和分析以确定地下结构层的性质和结构。
根据记录数据形成地震剖面图,地震剖面图是一个地震图,它沿地震轮廓线图描绘地球剖面的地下分层,地震剖面图是地质学家研究确定地球地面下结构层特性的重要工具。而且,在地震采样阵列转换成地震剖面图之前,这些地震采样阵列亦可由地质学家解释,对地震数据必须处理,以减少噪声造成的退化。
地震分析通常包括对到达时间、到达振幅、到达速率、频率的特性的研究及从目标层位反射特性的研究。任何变化或异常的特性都特别引起人们兴趣。
多次反射能量,通常称为“多次反射波”(“multiples”)被称为地球物理学现象,这些多次反射波通常定义为已经反射多于一次的地震能量。多次反射波通常在一层内形成,这一层由很多不同密度的层定界。例如,由水底和水面对水定界。因此,水底趋于作为多次反射波的共源,虽然多次反射波亦能出现在陆地数据中。
当通过水旅行的地震信号的一部分被反射离开水底时,水底多次反射波的发生就出现。然后,这个反射信号在水底中向水面旅行,在水面经反射离开水面返回水底。在水底和水面之间的反射可以发生一次或多次。因为地震接收器仅测量信号的到达时间,而不在意旅行路径。而地震信号已经受迂回旅行路径,所以在地震记录中显得比它们在地面下所处实际位置更深。
传统的现有技术方法从地震数据中去除多次反射波,例如,预测的解卷积通常已被利用,以适当处理多次反射波的周期性。即在第一次旅行后,在记录线中的多次反射波将在每n毫秒重复一次。例如,在预测的解卷积中,来自地震记录线的较早部分的信息用来预测和解卷积记录线的后面部分的信息。
还有一种传统的技术方法,例如正常时差校正(normalmoveout),这些方法要求起始和多次反射波之间有一明显的速度差异。例如,对相同到达时间长路径的多次反射波旅行平均速度比起始反射波速度低,它们显示出更大的正常时差校正,且随共深点(common-depth-point)迭层而被减弱。
遗憾地,这些周期性的或明显的速度差别的前提限制了应用和/或限制了预测解卷积或正常时差校正技术的效率。
用来消除多次反射波作用的另一种技术方法。包括波方程方法,这种方法产生更严格的预测多次旅行时的处理方法,波方程方法能处理复杂的几何结构且无需地面下速度的知识。遗憾的是涉及的波方程方法的计算代价严重地限制了它的使用。
因此,在现有的技术中,需要处理包括多次反射波的地震数据的替换方法。
另外,在本技术领域中需要一种改进处理包括多次反射波的地震数据的方法,以克服技术方法的限制。
在现有技术中还需要一种改进处理包括多次反射波噪声的地震数据的方法,这种方法同数据的周期性无关。
在现有的技术中甚至还需要一种改进处理包括多次反射波噪声的地震数据的方法,这方法同起始波至和它的多次反射波之间的速度差无关。
在现有的技术中甚至还有另外要求,以改进处理包括多次反射波噪声的地震数据的方法,这些方法并无过分的计算耗费。
在本领域中的这些要求及其它要求,对熟知此技术的人来说,通过对本专利说明的回顾将会很清楚。
本发明的一个目标是为处理具有多次反射波的地震数据提供一种可选择的方法。
本发明的另一目标是为处理具有多次反射波的地震数据提供一种方法,以克服现有技术的限制。
本发明的再一个目标是为处理具有多次反射波的地震数据提供一种方法,这种方法并不需要在初始波至和它的多次反射波之间有周期性。
本发明甚至还有一个目标是为处理具有多次反射波的地震数据提供一种方法,这种方法并无过分的计算耗费。
本发明还有一个目标是为处理具有多次反射波的地震数据提供一种方法,这方法并不依靠初始波至和它的多次反射波之间速度差。
本发明的这些目标及其它目标对熟知地震数据处理技术的人,通过阅读本专利说明将会很清楚的。
为实现本发明的上述目的,本发明提供了一种处理地面下的地震数据的方法,其中具有多次反射波生成区,此区域至少在一个面同多次反射波生成波至邻接,且具有初始目标波至,其中地震数据包括由多次反射波生成波至反射出并在多次反射波生成区旅行的地震波生成的多次反射能量,而其中多次反射波生成波至的零偏移双向地震波旅行时、多次反射波生成区中地震波的速度、以及初始波至的均方根地震波速度是公知的,该方法包括:
(a)利用初始目标波至的均方根地震波速度,通过应用地震数据的正常时差校正操作产生偏移校正数据集;
(b)通过按时间下移偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的零偏移双向旅行时;
(c)利用对步骤(b)的时间移位数据集的反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向正常时差校正操作利用一个速度函数,此速度函数是初始目标波至的均方根地震波速度及多次生成区中地震波速度的函数;且
(d)用从地震数据中减去多次反射波反射能数据集产生校正数据集。
本发明还提供了一种处理地面下的地震数据的方法,所说地面下具有的多次反射波生成区至少一面同多次反射波生成波至邻接,且具有初始目标波至,其中地震数据包括由多次反射波生成波至反射出并在多次反射波生成区旅行的地震波生成的多次反射波反射能量;该方法包括:
(a)确定针对多次反射波生成波至的零偏置双向地震波旅行时,确定多次反射波生成区内地震波的速度及针对初始波至确定均方根地震波速度;
(b)利用对初始波至的均方根地震波速度通过应用正常时差校正操作产生一个偏移校正数据集;
(c)通过按时间下移的偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的零偏双向旅行时;
(d)以步骤(c)的时间移动数据集应用反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向动样正操作利用一个速度函数,该速度函数是初始目标波至均方根地震波速度及多次反射波生成区中的地震波速度的函数;且
(e)用从地震数据中减去多次反射波反射能数据集产生校正数据集。
本发明还提供了对地震数据产生多次反射波反射能量数据集的方法,该地震数据表示地下具有多次反射波生成波至,初始目标波至,其中的地震数据包括由多次反射波生成波至生成的多次反射波反射能量,其中针对多次反射波生成波至的双向零偏移地震能量旅行时,针对多次反射波生成波至的地震能量速度是已知的,这方法包括:
(a)利用初始目标波至的正常时差校正速度,通过应用地震数据的正常时差校正操作产生偏移校正数据集;
(b)通过按时间下移的偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的双向零偏移旅行时;且
(c)对步骤(b)的时间移动数据集的反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向正常时差校正操作的速度是初始目标波至的正常时差校正速度和多次反射波生成波至地震能量速度的函数。
根据本发明的优选实施例,提供了一种从包含多次反射波的地震数据中分离多次反射数据的计算简单而特殊的方法,在利用计算机时该方法是优先采用的。这种方法通常包括这样一步骤,测定零偏移双向(two-way)旅行时以检测多次反射产生的波至。这种方法通常要包括这样一步骤,通过一个正常时差校正操作的应用产生一个偏移校正数据集,并利用初始波至速度记录地震记录线。一旦正常时差校正已用于数据,下一步骤就包括向下移动偏移校正数据集中的所有波至,产生一时间移动数据集。正像前面所定的那样,时间移动量等于推测的多次反射波生成波至的双向旅行时。一旦正常时差校正操作和双向旅行时移动已被用于数据,该方法包含通过把反向正常时差校正操作用于时移数据集,利用扩充速度函数产生一多次反射波反射能量数据集。此函数是原来的初始速度函数和推测的多次反射波生成波至的双向旅行时的间隔速度的函数。结果数据集现在仅包括多次反射波的数据集。
根据本发明的另一实施例,其中包括提供处理多次反射波地震数据的计算简单而特殊的一种方法。这在利用计算机时是优选的方法。这个方法包括从地震数据中扣除多次反射波反射能量数据集。在优选实施例中,多次反射波反射能量数据是这样进行换算,使扣除的多次反射波反射能量数据为最小。
图1是表示本发明的方法的简图,图中显示了起始波至PE,水底WB,水面WS,源S和接收器R。
图2是一附图,它显示了关于正常时差校正操作的地表面上在源S、接收器G(中点为M)之间单层几何关系。
图3是一附图,它显示了关于正常时差校正操作的地表面上在源S、接收器G(中点M)之间多层几何关系。
图4是一附图,它显示了与图2或3的几何图形有关的综合共中心点(GMP)的记录线集。在图2或3中,对一平坦反射层的旅行时曲线是一双曲线,该曲线的顶点在零偏移点。
图5(a)和5(b)是一附图,它描述了非零偏移旅行时转成零偏移旅行时的图。图5(a)为正常时差校正之前的记录线集,图5(b)为正常时差校正之后的记录线集。
图6(a)和6(b)是一附图,它显示了利用计算机对正常时差校正传输曲线的操作,对这些情况下,在t(x)值上的振幅A并不需要落在整数采样位置上。
图7(a)是例1中综合以上数据的记录线集。
图7(b)是显示了在正常时差校正操作之后图7(a)的数据。
图7(c)是一记录线集,它显示了通过向水底推移0.2秒双向旅行时以后,图7(b)的正常时差校正数据。
图7(d)是一记录线集,它显示了用反向NMO方程,用图8的扩充速率函数,在反向正常时差校正操作应用之后图7(c)的数据。
图8是对例1的数据扩充速率函数的图。
图9(a)到9(c)是记录线集,它们分别显示了图7(a)的输入数据,图7(d)的噪声数据及利用本发明的方法之后的数据。
图9(d)是应用NMO及重叠以后,图9(a)的记录线数据及类似数据的重叠剖面。而图9(e)是依靠本发明的方法除去多次反射波以后图9(d)的记录线数据及类似数据的重叠剖面图。
在本发明的方法的实施中,可用任何适当方法获得地震数据,例如许多熟知的地震测定技术的任意一种。通常,声波源如在陆地上的甘油炸药或振动器或海上空气枪所提供的声波向下传输到地里再从地面下反射界面向上反射。这些反射波被诸如地震仪,小型地震仪或水中地震检波器这些检测器接收,并产生地震信号或地震记录线的模拟电信号。所收集到的数据然后用计算机经过各种地震处理技术处理。
本发明的方法可以参照图1作简要解释。图1是本发明的方法的图示表示。显示了初始波至PE,水底WB,水面WS,源S和接收器R。
在图1中,推测的多次反射波形成波至为水底WB,接收器R,源S有距离为S的偏移。旅行的路径P是从源S旅行到初始波至PE再到接收器R。符号a表示一正常时差校正操作,其中旅行路径P用偏移C校正。初始的旅行时现在同在零偏移时采集的数据相同,符号d是双向旅行时,它指从水面WS到推测的多次反射波生成波至,即水底WB。符号b表示通过双向旅行时d对路径P移动零偏移旅行时,反向正常时差校正然后应用到移动的零偏移旅行时b,以产生具有偏移c的路径p’,该路径表示一预测的水底多次反射波,这预测水底多次反射波被从地震数据中扣除。
处理包含多次反射波地震数据的本发明的实施中,通常需要某些有关初始波至的信息,尤其是那些形成任何多次反射的推测。一旦一个多次反射波生成的初始波至被确定,本发明的方法包括这一步骤,以确定双向零偏移旅行时的推测多次反射波生成波至。例如,在海上地震数据处理中,至少推测的多次反射波生成初始波至中一个将在水底形成,因此,需要确定到水底的双向旅行时。此双向旅行时到一推测多次反射波生成初始波至利用地震技术中所熟知的技巧很容易确定。
本发明的实施亦包括这一步骤,把正常时差校正操作(通常称“NMO”)应用到利用初始速度的地震记录线。沿接收器的线记录的地震数据的正常时差校正通常认为由于射点到接收器距离偏移变化,改变了反射到达时间而引起。正常时差校正是速度的函数,对较小范围的偏移,下沉随反射时间降低。
本发明的实施中,动态校正还可用于本技术领域任何已知技艺的数据。例如,正常时差校正可以用于由Ozdogan Yilmaz在《地震数据处理》中第157至162页之间所述的地震数据,该书于1987年由勘探地理学家协会(Society of Exploration Geophysicists)出版。在此仅作为参考。
关于正常时差校正的概念可进一步参考图2的说明。图2中显示了单一水平层的简单情况下关于正常时差校正操作在地面上对源S、接收器G中点M之间的几何关系。
在所给中点位置M上,从发射位置S沿射线路径到一深度点D,然后再返回到接收器位置G的旅行时是t(x)。用勾股定理作为偏移函数的旅行时方程是
这里x是源和接收器位置之间距离(偏移)
v是反射界面上面介质的地震速度
t(o)是沿垂直路径MD2次旅行时。
如图2所示,深度点D到表面的垂直投影,(沿法线到的反射层)同中点M重合。这仅仅适合当反射层是在点P水平向时的情况。方程1描述了双向时间面同零偏移点为顶点双曲线关系。
多次反射波在多层地下的结构层中形成,典型的例子是由水底多次反射波被畸变的初始波至。对多层正常时差校正如下说明。
图3所示水平向等速度层组成的地下结构,地表面上关于正常时差校正操作在源S,接收器G,中点M之间的几何关系。在图3中,这些层的厚度依靠双向零偏移时间定义。这N层具有间隔速度V1,V2…Vn。
对射线路径从源S到深度点D,返回到接收器R的旅行时方程为
t2(x)=c0+c1x2+c2x4+c3x6+… (方程2)
这里C0=t2(0)
c1=1/Vrms 2
c2,c3…是取决于层厚和间隔速度的复杂函数。rms(均方根)速度Vrms下到位于深度点D所在反射层。定义为
这里Δti是通过第i层的垂直双向时间
当形成的偏移量近似值相对比是微小深度时,方程3各项可以如下项,
t2(x)=t2(o)+x2/Vrms 2 (方程4)
比较方程1和方程4表明为正常时差校正所需的速度,对于偏移同深度比较小时,水平分层地面下的速度等于均方根速度。
现参考图4,其中显示图2的几何结构,按其中心点(CMP)采集的记录线集的例子。图4还表示一共深点(CDP)记录线集,因为所有的同每个源接收器对有关的射线路径从相同的地面下深度点D反射。所有在此共中心点记录线包括来自同深度点的反射,在给定偏移x的双向时间t(x)和双向零偏移时间t(o)之间差称为正常时差校正(“NMO”),从方程4知,当偏移x和双向时间t(x)和t(o)已知时,均方根速度可计算出来。
一旦正常时差校正速度被估算出,旅行时能被校正以除去偏移的影响。现在参考图5(a)和图5(b),正常时差校正包括非零偏移旅行时t(x)到零偏移旅行时t(o)的变换,对图5(a)显示了正常时差校正之前的记录线集。图5(b)显示了正常时差校正之后的记录线集。然后,对NMO校正的记录线集中记录线累加以得到一特殊CMP位置上一迭合记录线集,其说明如下:
这种正常时差校正操作的应用最适宜利用计算机。然而,当用计算机时,对旅行时曲线t(x)的时间值A没有落在相应于x的时间值的整数位置时会产生某些问题。
现在参考图6(a)和6(b),其中显示利用计算机对正常时差校正所用技术。按所给的整数值t(o),速度和偏移,用方程4计算t(x),在那些例子中,作为旅行时曲线t(x)的时间值并不发生于相应于x时间值的输入的整数取样位置上,t(x)的每边的整数值(用实点表示整数值)被用来内插到整数值之间以计算t(x)的时间值。然后,这计算的时间值映射成相应偏移上整数的样值t(o)(用A’表示)。作为一般规则,正常时差校正随着偏移而增加,随着深度而减少。
一旦正常时差校正用于地震数据中,对于推测多次反射生成波至的双向旅行时就被确定,在本发明的实施例中的下一步是通过等于双向零偏移旅行时的时间量,按时间向下移动记录线中所有波至,对预测的多次反射波生成波至。
一旦通过双向零偏移旅行时对预测的多次反射波生成波至已按时间向下移动数据。一个反正常时差校正操作被用于利用扩充的速度函数记录线中。这个函数是起动速度函和对预测的多次反射波生成波至间隔速度的函数。
反向正常时差校正可以从方程4中确定的t2(x)得到,方程4中所用扩充速度函数为均方根速度,而时间值已由双向旅行时移到多次反射波生成波至。
作为时间t的函数扩充速度函数如下所示
Vextended(t)={(V2t+Vm 2tm)/(t+t(o))}1/2………(方程5)
这里Vm是多次反射生成层内的速度
tm是多次反射生成层内的旅行时
V是对关心的起始波至的均方根速度
t(o)是对多次反射生在波至的零偏移双向旅行时。
利用方程5的扩充速度函数和时间移动地震数据,可得到如下反向正常时差校正方程
t2(x)=(t(o)+tm(x))2+(x2/vextended 2(t))(方程6)
在反向正常时差校正操作应用后,可得到仅包含多次反射波的多次反射波记录线集。下一步,在消除多次反射噪声时是使原始记录线和如下所有时间深度的多次记录线之间能量差最小。
在原始记录线和所有时间深度的多次记录线之间能量差E亦称为“被扣除的记录线能”(“subtracted ttace energy”)并示于下面方程7中。其中Di为原始数据记录线集
Mi为多次反射波记录线集
所有下标i为时间指数
本发明的方法假设多次反射波,用小于1的反射系数,从界面反射。因此,利用振幅定标因数α。而且,剩余时间移动τ。是在最佳化中考虑用确定双向旅行时去推测多次反射波生成波至的误差。
对所给τ,α可以对方向7对α求导,设置方程7的导数为0得到
对α求解方向8*,得到方程9
α=(ΣDiMi)/(∑Mi 2) (方程9)。
一旦α和τ相结合,最小被扣除的记录线能E被确定,多次反射波记录线集被从原始数据记录线集中减去得到一无多次反射波的记录线集。
例1
本例子利用50米记录线间隔的综合海上数据,在本例中,双向零偏移水底时间为0.2秒,第一初始波至在1.8秒时,第一个形成的浅层,多次反射波是在2.0秒时。
图7(a)是本例的输入数据的记录线集,图7(b)显示在应用了正常时差校正后图7(a)的数据集。图7(c)显示在用0.2秒双向旅行时移动数据到水底以后图7(b)的数据。
用上面方程2和水底时间的间隔速度确定扩充的速度函数。对本例所确定的扩充速度函数,显示在图8中。
图7(d)显示用上面方程6及在图8的扩充速度函数,在应用反向正常时差校正操作之后的图7(c)的数据显示的仅是生成的多次反射波。
图9(a)到9(c)分别显示图7(a)的输入数据,图7(d)的噪声数据,及在应用本发明的方法后的数据。噪声数据用常数0.8调节,它用最小被扣除的记录线能来确定。
图9(d)是在应用NMO和迭合以后图9(a)的记录线集数据和相似数据的迭合剖面。而图9(e)是用本发明的方法,移动多次反射波以后图9(c)的记录线集数据和类似数据的迭合剖面。在图9(d)和图9(e)之间最值得注意的差别是在图9(d)中在2.0秒处存在水底多次反射,而在图9(e)处没有这种情况。
在此所给的说明想描述本发明的优选实施例。对于在本技术领域具有一般技艺的人可以在相当于本发明权利要求范围内对本发明的细节做各种变化,但并不脱离本发明的精神。因此,这意味着所有这些变化包含在本发明的权利要求范围之中。
Claims (14)
1.一种处理地面下的地震数据的方法,其中具有多次反射波生成区,此区域至少在一个面同多次反射波生成波至邻接,且具有初始目标波至,其中地震数据包括由多次反射波生成波至反射出并在多次反射波生成区旅行的地震波生成的多次反射能量,而其中多次反射波生成波至的零偏移双向地震波旅行时、多次反射波生成区中地震波的速度、以及初始波至的均方根地震波速度是公知的,该方法包括:
(a)利用初始目标波至的均方根地震波速度,通过应用地震数据的正常时差校正操作产生偏移校正数据集;
(b)通过按时间下移偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的零偏移双向旅行时;
(c)利用对步骤(b)的时间移位数据集的反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向正常时差校正操作利用一个速度函数,这速度函数是初始目标波至的均方根地震波速度及多次生成区中地震波速度的函数;且
(d)用从地震数据中减去多次反射波反射能数据集产生校正数据集。
2.权利要求1的方法,用步骤(c)中的速度函数v(t)表示如下:
v(t)=((v2t+vm 2tm)/(t+t(o))}1/2这里vm是多次反射波生成区中地震波的速度
tm是多次反射波生成区中地震波的旅行时
v是感兴趣的初始波至的均方根速度
t(o)是多次反射波生成波至的零偏移旅行时。
3.权利要求1的方法中,多次反射波生成区域是水,多次反射波生成波至是水底。
4.权利要求1的方法中,步骤(d)的校正数据集是通过从地震数据中减去多次反射波反射能数据集生成,其中多次反射波反射能数据集是通过振幅定标因数和剩余时间移动来调节。
5.权利要求4的方法中振幅定标因数和剩余移动被选择为多次反射波反射能数据集和地震数据之间相减能量的最小平方,方程如下: 这里E是被扣除的能量
α是振幅定标因数
τ是剩余时间移动
D是地震数据
M是多次反射能数据集
i是时间指数。
6.权利要求1的方法中,多次反射波生成区是水,多次反射波生成波至是水底,用在步骤(c)中的速度函数v(t)如下所示:
v(t)={(v2t+vm 2tm)/(t+t(o))}1/2这里vm是多次反射波生成区中地震波速度
tm是多次反射波生成区中地震波速度
v是感兴趣的初始波至的均方根速度
t(o)是多次反射波生成波至的零偏移旅行时;而步骤(d)的校正集是从地震数据中减去反射能量数据集产生的,其中多次反射波反射能量数据集通过振幅定标因数和剩余时间移动来调节,振幅定标因数和剩余时间移动选择成多次反射波反射能量数据集和地震数据之间相减能量的最小平方,方程如下: 这里E是被扣除的能量
α是振幅定标因数
τ是剩余时间移动
D是地震数据
M是多次反射波反射能数据集
i是时间指数。
7.处理地面下的地震数据的方法,所说地面下具有的多次反射波生成区至少一面同多次反射波生成波至邻接,且具有初始目标波至,其中地震数据包括由多次反射波生成波至反射出并在多次反射波生成区旅行的地震波生成的多次反射波反射能量;这方法包括:
(a)确定针对多次反射波生成波至的零偏置双向地震波旅行时,确定多次反射波生成区内地震波的速度及针对初始波至确定均方根地震波速度;
(b)利用对初始波至的均方根地震波速度通过应用正常时差校正操作产生一个偏移校正数据集;
(c)通过按时间下移的偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的零偏双向旅行时;
(d)以步骤(c)的时间移动数据集应用反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向正常时差校正操作利用一个速度函数,这速度函数是初始目标波至均方根地震波速度及多次反射波生成区中的地震波速度的函数;且
(e)用从地震数据中减去多次反射波反射能数据集产生校正数据集。
8.权利要求7的方法中,用在步骤(d)中的速度函数v(t)表示如下
v(t)={(v2t+vm 2tm)/(t+t(o))}1/2这里vm是多次反射波生成区中地震波的速度
tm是多次反射波生成区中地震波的旅行时
v是感兴趣的初始波至的均方根速度
t(o)是多次反射波生成波至的零偏移旅行时。
9.权利要求7的方法中多次反射波生成区域是水,而多次反射波生成波至是水底。
10.权利要求7的方法中步骤(e)的校正数据集通过从地震数据中减去多次反射波反射能数据集生成,其中多次反射波反射能量数据集用振幅定标因数和剩余时间移动来调节。
11.权利要求10的方法中振幅定标因数和剩余时间移动被选择为多次反射波反射能数据集和地震数据之间相减能量的最小平方,方程如下: 这里E是被扣除的能量
α是振幅定标因数
τ是剩余时间移动
D是地震数据
M是多次反射波反射能数据集
i是时间指数
12.对地震数据产生多次反射波反射能量数据集的方法,该地震数据表示地下具有多次反射波生成波至,初始目标波至,其中的地震数据包括由多次反射波生成波至生成的多次反射波反射能量,其中针对多次反射波生成波至的双向零偏移地震能量旅行时,针对多次反射波生成波至的地震能量速度是已知的,这方法包括:
(a)利用初始目标波至的正常时差校正速度,通过应用地震数据的正常时差校正操作产生偏移校正数据集;
(b)通过按时间下移的偏移校正数据集中的所有波至,产生时移数据集,所移时间量等于对多次反射波生成波至的双向零偏移旅行时;且
(c)对步骤(b)的时移数据集的反向正常时差校正操作产生一多次反射波反射能量数据集,其中反向正常时差校正操作的速度是初始目标波至的正常时差校正速度和多次反射波生成波至地震能量速度的函数。
13.权利要求12的方法,用在步骤(d)中的速度函数v(t)表示如下:
v(t)={(v2t+vm 2tm)/(t+t(o))}1/2这里vm是多次反射波生成区中地震波的速度
tm是多次反射波生成区中地震波的旅行时
v是感兴趣的初始波至的均方根速度
t(o)是多次反射波生成波至的零偏移旅行时。
14.权利要求12的方法中多次反射波生成区域是水,而多次反射波生成波至是水底。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US300474 | 1994-09-02 | ||
US08/300,474 US5521881A (en) | 1994-09-02 | 1994-09-02 | Method of processing seismic data having multiple reflection noise |
US300,474 | 1994-09-02 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1118441A CN1118441A (zh) | 1996-03-13 |
CN1073236C true CN1073236C (zh) | 2001-10-17 |
Family
ID=23159244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN95116803A Expired - Fee Related CN1073236C (zh) | 1994-09-02 | 1995-08-31 | 处理具有多次反射噪声的地震数据的方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US5521881A (zh) |
CN (1) | CN1073236C (zh) |
AU (1) | AU687590B2 (zh) |
GB (1) | GB2292801B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630019B (zh) * | 2008-07-16 | 2011-12-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种校正预测的海底多次波大时差的方法 |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6427124B1 (en) | 1997-01-24 | 2002-07-30 | Baker Hughes Incorporated | Semblance processing for an acoustic measurement-while-drilling system for imaging of formation boundaries |
US6023443A (en) * | 1997-01-24 | 2000-02-08 | Baker Hughes Incorporated | Semblance processing for an acoustic measurement-while-drilling system for imaging of formation boundaries |
US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
GB9906456D0 (en) * | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
US6442490B1 (en) * | 2000-06-30 | 2002-08-27 | Pgs Americas, Inc. | Vertical cable time processing |
US6906982B2 (en) * | 2002-04-16 | 2005-06-14 | Shell Oil Company | Marine time-lapse seismic surveying |
US6636810B1 (en) * | 2002-05-24 | 2003-10-21 | Westerngeco, L.L.C. | High-resolution Radon transform for processing seismic data |
US6873913B2 (en) * | 2002-09-20 | 2005-03-29 | Paradigm Geophysical Corporation | Method of reconstructing seismic records to obtain high resolution signals |
EP1620752A2 (en) * | 2003-05-02 | 2006-02-01 | WesternGeco, L.L.C. | Method for computing and using timing errors that occur in multiples predicted by multiple prediction algorithm |
US7181347B2 (en) * | 2003-09-23 | 2007-02-20 | Westerngeco, L.L.C. | Method for the 3-D prediction of free-surface multiples |
WO2005052639A1 (en) * | 2003-10-28 | 2005-06-09 | Western Geco, Llc | A method for estimating porosity and saturation in a subsurface reservoir |
US7791980B2 (en) * | 2004-05-21 | 2010-09-07 | Westerngeco L.L.C. | Interpolation and extrapolation method for seismic recordings |
CN100349005C (zh) * | 2004-12-29 | 2007-11-14 | 中国石油天然气集团公司 | 地震数据多炮组合叠前深度偏移方法 |
CN100349009C (zh) * | 2004-12-29 | 2007-11-14 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN101285894B (zh) * | 2008-05-30 | 2011-02-09 | 中国科学院地质与地球物理研究所 | 起伏地表下采集的地震资料的直接叠前时间偏移方法 |
CN101334483B (zh) * | 2008-06-13 | 2011-01-26 | 徐基祥 | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 |
CN102081168B (zh) * | 2009-12-01 | 2012-07-18 | 中国石油天然气集团公司 | 一种提高地震数据处理中面波衰减性能的方法 |
US9014985B2 (en) * | 2010-04-30 | 2015-04-21 | Chevron U.S.A. Inc. | System and method for compensating time and offset varying near-surface effects in seismic data background |
CN102305945B (zh) * | 2011-06-20 | 2013-09-25 | 电子科技大学 | 一种线性噪声消除方法 |
WO2014139585A1 (en) * | 2013-03-15 | 2014-09-18 | Fmc Kongsberg Subsea As | Method for determining a position of a water/cement boundary between pipes in a hydrocarbon well |
WO2014140363A2 (en) * | 2013-03-15 | 2014-09-18 | Fmc Kongsberg Subsea As | Method for determining a position of a water/cement boundary between pipes in a hydrocarbon well |
US20160334526A1 (en) * | 2014-01-13 | 2016-11-17 | Cgg Services Sa | Method and apparatus for modeling and separation of primaries and internal multiples using solutions of the two-way wave equation |
CN110618457B (zh) * | 2018-06-20 | 2021-02-26 | 中国石油化工股份有限公司 | 一种可控震源地震数据的联合去噪方法及其装置 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4189704A (en) * | 1975-09-29 | 1980-02-19 | Geophysical Systems Corp. | Method and apparatus for determination of optimum velocity in processing seismic signals from low energy sources |
US4319347A (en) * | 1976-03-08 | 1982-03-09 | Western Geophysical Co. Of America | Seismic method and system of improved resolution and discrimination |
DE2658954C2 (de) * | 1976-12-24 | 1979-04-12 | Deutsche Texaco Ag, 2000 Hamburg | Verfahren zur Unterdrückung von Störwellenzügen in seismischen Daten |
US4209854A (en) * | 1978-05-04 | 1980-06-24 | Mobil Oil Corporation | Method and system for enhancing primary seismic reflection signals in common depth point sets |
US4887243A (en) * | 1982-12-27 | 1989-12-12 | Mobil Oil Corporation | Removal of surface multiples |
US4937794A (en) * | 1985-05-06 | 1990-06-26 | Western Atlas International, Inc. | Seismic noise suppression method |
US4802146A (en) * | 1985-05-23 | 1989-01-31 | Mobil Oil Corporation | Method for moveout correction and stacking velocity estimation of offset VSP data |
US4922362A (en) * | 1988-03-04 | 1990-05-01 | Schlumberger Technology Corporation | Methods for deconvolution of unknown source signatures from unknown waveform data |
US4907205A (en) * | 1988-04-29 | 1990-03-06 | Amoco Corporation | Method for attenuating multiple reflection events in seismic data |
FR2653900B1 (fr) * | 1989-10-30 | 1992-02-28 | Total Petroles | Procede perfectionne de prospection sismique marine. |
US4992993A (en) * | 1990-06-18 | 1991-02-12 | Western Atlas International, Inc. | Correction for variable water-column velocity in seismic data |
US5150331A (en) * | 1991-03-25 | 1992-09-22 | Amoco Corporation | Method for enhancing seismic data |
US5404296A (en) * | 1991-09-25 | 1995-04-04 | Tinking Machines Corporation | Massively parallel computer arrangement for analyzing seismic data pursuant to pre-stack depth migration methodology |
US5365492A (en) * | 1993-08-04 | 1994-11-15 | Western Atlas International, Inc. | Method for reverberation suppression |
-
1994
- 1994-09-02 US US08/300,474 patent/US5521881A/en not_active Expired - Lifetime
-
1995
- 1995-08-24 GB GB9517386A patent/GB2292801B/en not_active Expired - Fee Related
- 1995-08-31 CN CN95116803A patent/CN1073236C/zh not_active Expired - Fee Related
- 1995-08-31 AU AU30392/95A patent/AU687590B2/en not_active Ceased
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630019B (zh) * | 2008-07-16 | 2011-12-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种校正预测的海底多次波大时差的方法 |
Also Published As
Publication number | Publication date |
---|---|
GB9517386D0 (en) | 1995-10-25 |
AU687590B2 (en) | 1998-02-26 |
US5521881A (en) | 1996-05-28 |
CN1118441A (zh) | 1996-03-13 |
GB2292801B (en) | 1998-02-18 |
GB2292801A (en) | 1996-03-06 |
AU3039295A (en) | 1996-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1073236C (zh) | 处理具有多次反射噪声的地震数据的方法 | |
CN1271420C (zh) | 海洋地震数据处理方法和地震测量方法 | |
CN1181358C (zh) | 用于局部浪高确定的方法和在海洋地震信号中的声传感器 | |
CN103477246B (zh) | 具有减少的拖影和/或增大的失真容限的海洋振动器扫描 | |
CN1065630C (zh) | 地震激振器的串接扫描法 | |
CN1275049C (zh) | 海洋地震勘测 | |
CN1275050C (zh) | 一种地震数据处理方法和装置 | |
CN1015672B (zh) | 消除地震数据中的多普勒相位色散的方法 | |
CN1013226B (zh) | 多分量地震数据深度成象法 | |
CN1011350B (zh) | 垂直地震剖面的上、下行地震事件的区分法 | |
CN1246623A (zh) | 海洋地震学测量中船舶移动影响的校正方法及装置 | |
CN1205084A (zh) | 源分离的高保真度震动源地震探矿法 | |
CN1797037A (zh) | 一种地震波波阻抗反演的方法 | |
CN1404581A (zh) | 估算地震介质特性的系统和方法 | |
EA026344B1 (ru) | Система и способ получения и обработки сейсмических данных о полях упругих волн | |
CN104533396A (zh) | 一种远探测声波的处理方法 | |
Pratt | Characterizing and imaging sedimentary strata using depth‐converted spectral ratios: An example from the Atlantic Coastal Plain of the Eastern United States | |
MX2014010606A (es) | Extraccion de datos de cizallamiento sv a partir de datos marinos de la onda p. | |
CN101251602B (zh) | 海洋组合可控震源短时分段扫描方法 | |
Pan et al. | Reflection signals and wellbore scattering waves in acoustic logging while drilling | |
CN87103519A (zh) | 具有极高分辨率的在水平钻孔中的地震勘探方法 | |
CN1656388A (zh) | 处理地震数据的方法和装置 | |
CN1013000B (zh) | 从已获得的压缩波反射数据估算切变波反射数据的方法 | |
Gamal et al. | Validity of the refraction microtremors (ReMi) method for determining shear wave velocities for different soil types in Egypt | |
Gómez de la Peña | The origin and tectono-sedimentary structure of the Alboran Basin |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent for invention or patent application | ||
CB02 | Change of applicant information |
Applicant after: Exxon Production Research Co. Applicant before: Exxon Production Research Co. |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: EXXON PRODUCTION RESEARCH CO. TO: EXXONMOBIL UPSTREAM RESEARCH COMPANY |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |