CN101630019B - 一种校正预测的海底多次波大时差的方法 - Google Patents
一种校正预测的海底多次波大时差的方法 Download PDFInfo
- Publication number
- CN101630019B CN101630019B CN2008101167054A CN200810116705A CN101630019B CN 101630019 B CN101630019 B CN 101630019B CN 2008101167054 A CN2008101167054 A CN 2008101167054A CN 200810116705 A CN200810116705 A CN 200810116705A CN 101630019 B CN101630019 B CN 101630019B
- Authority
- CN
- China
- Prior art keywords
- prediction
- collection
- time difference
- bottom echo
- big gun
- 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.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是地震资料处理中校正预测的海底多次波大时差的方法。用多次波预测算子生成预测的海底反射波和多次波模型道集,自动拾取海底反射波的旅行时进行中值滤波,计算原始炮集和预测的海底反射波的旅行时差,对多次波模型道集做时差校正,应用时差校正后的多次波模型道集与原始数据进行数据匹配,从原始数据中减去匹配后的多次波道集,输出压制多次波后的地震剖面。本发明是在常规数据匹配之前先对大的时差进行校正的方法,因而能够适应存在拖缆漂移的采集数据,提高自适应减去法的效率,改善压制多次波的效果。
Description
技术领域
本发明涉及地球物理勘探技术,属于地震资料处理中数据匹配的技术范畴,是一种校正预测的海底多次波大时差的方法。
背景技术
在海上地震勘探中,由于作为自由界面的海水面和硬的海底面都是很强的反射界面,其中海水面的反射系数可达0.9,因此在地震记录中不可避免的会记录到海底和海水面产生的各种各样的多次波。多次波的存在会干扰人们对有效波的识别,影响后续的速度分析、偏移成像、解释的真实性和可靠性,因此必须加以压制。
压制多次波的方法一般分为两大类:一类是基于有效波和多次波之间差异的滤波方法,简称为滤波方法;另一类是基于波动方程的预测减去法,即首先利用原始地震数据,根据波动原理,预测出多次波,然后对预测的多次波与原始数据中的多次波进行匹配,从原始数据中减去匹配后的多次波就可以得到一次有效波。
在基于波动方程的预测减去方法中,根据波动原理利用原始地震数据预测的多次波,与原始数据中的多次波之间会存在振幅、相位和旅行时差异,因此需要对预测的多次波以原始数据中的多次波为目标进行匹配,然后用原始炮集减去匹配后的预测的多次波就达到消除或压制多次波的目的。
数据匹配一般要求匹配滤波算子的长度较短,因为长的匹配滤波算子要求目标数据,包括预测的多次波和原始地震数据的长度也要长,而当目标数据较长时,就有可能包含不期望参与数据匹配的信号和干扰波,进而影响匹配效果。然而在实际海上拖缆地震数据采集中,受海流的影响,拖缆漂移是不可避免的,因此一方面由拖缆漂移的地震数据得到的海底高程会存在误差,另一方面当把拖缆漂移的地震道投影到测线时,也不能真实地反映地震道的位置,这样就有可能导致预测的海底多次波与实际资料中的海底多次波存在较大的时差,这种时差客观上需要很长的匹配滤波算子来进行处理。因此,为了减少匹配算子的长度,提高匹配的效率与精度,在数据匹配之前,必须消除预测的海底多次波大的时差。
发明内容
本发明目的在于提供一种能够适应存在拖缆漂移的采集数据,提高自适应减去法的效率,改善压制多次波的效果的校正预测的海底多次波大时差的方法。
本发明的具体实施步骤如下:
1)采集炮集地震记录;
2)用多次波预测算子对输入原始炮集地震记录进行运算,生成预测的海底反射波和多次波模型道集;
步骤2)中所述的输入的原始炮集记录中要保留直达波;
3)利用原始炮集和预测的海底反射波和多次波模型道集自动拾取海底反射波的旅行时;
步骤3)中拾取海底反射波的旅行时是拾取海底反射波同相轴的初至对应的时间;
4)对拾取的海底反射波旅行时进行中值滤波:
步骤4)中所述的中值滤波是指,将某点为中点的窗口中的所有点的值由小到大进行排序,中间值作为该点的输出值。
步骤4)中所述的中值滤波窗口的长度取7-11个采样点。
5)计算原始炮集和预测的海底反射波的旅行时差,对多次波模型道集做时差校正;
步骤5)所述的计算原始炮集和预测的海底反射波的旅行时差是利用原始炮集中确定的海底反射波旅行时减去预测的海底反射波的旅行时。
步骤5)所述的时差校正逐道进行,按照多次波模型道与原始数据的海底反射时差整道上移或下移,向上移动时地震道尾部所缺部分充零,向下移动时地震道头部所缺部分充零。
6)应用时差校正后的多次波模型道集与原始数据进行数据匹配;
7)从原始数据中减去匹配后的多次波道集;
8)输出压制多次波后的地震剖面。
本发明不是采用增加匹配滤波算子长度的方法来校正匹配数据间的大的时差,而是在常规数据匹配之前先对大的时差进行校正的方法,因而能够适应存在拖缆漂移的采集数据,提高自适应减去法的效率,改善压制多次波的效果。
附图说明
图1原始炮集(左)及预测的海底反射波和海底多次波道集(右)。
图2原始炮集(左)及预测的海底反射波和海底多次波道集(右)中自动拾取的海底反射波的旅行时曲线。
图3原始炮集中的海底反射波及预测的海底反射波的旅行时之差。
图4对预测的海底反射波和海底多次波道集进行大的时差校正后的结果。
图5对预测的海底多次波进行大的时差校正和数据匹配后的结果。
图6原始炮集。
图7是原始炮集及大的时差校正后多次波自适应减去后的结果。
具体实施方案
在预测海底或自由表面多次波时,人们往往要把直达波切除,因为多次波预测算子会把直达波转换成海底反射波,这样在多次波自适应减去时,海底反射波也会被消除。
本发明的不切除原始炮集中的直达波,这样多次波预测算子就会把直达波转换成海底反射波,即预测的海底反射波;自动拾取原始炮集中和预测的海底反射波的旅行时,计算两者之间的时差;设原始炮集中和预测的海底多次波之间的时差与原始炮集中和预测的海底反射波之间的时差是相同的,这样就可以利用原始炮集中和预测的海底反射波旅行时之差,对预测的海底多次波进行大的时差校正。
本发明的具体实施方式为:
以下是本发明具体实例:
1)取得原始炮集地震记录。
2)按照步骤(2)对原始地震数据行处理得到预测的海底反射波和多次波模型道集
图1原始炮集及预测的海底反射波和海底多次波道集。可以看出原始炮集中的海底反射波和预测的海底反射波有较大的旅行时误差。
3)按照步骤3)-步骤4)分别在原始炮集、预测的海底反射波和海底多次波道集上拾取海底反射波的旅行时。
图2是分别在原始炮集及预测的海底反射波和海底多次波道集中自动拾取的海底反射波的旅行时曲线。
4)按照步骤5)计算原始炮集中的海底反射波和预测的海底反射波的旅行时误差,然后对预测的海底反射波和海底多次波道集进行大的时差校正。
图3是原始炮集的海底反射波及预测的海底反射波的旅行时之差,可以看出误差比较大,一般在-100.0ms-100.0ms之间。
图4是对预测的海底反射波和海底多次波道集进行大的时差校正后的结果。可以看出,经过校正,预测的多次波与原始数据中的多次波已经不存在大的旅行时误差。
5)按照步骤6)-步骤8)对时差校正后的预测多次波进行匹配,然后从原始数据中减去匹配后的多次波就可以得到有效波。
图5是对预测的海底多次波进行大的时差校正和数据匹配后的结果。
图6是原始炮集。
图7是原始炮集及大的时差校正后多次波自适应减去后的结果。可以看出,多次波的能量得到了很好的压制。
Claims (3)
1.一种校正预测的海底多次波大时差的方法,其特征在于实施步骤如下:
1)采集炮集地震记录;
2)用多次波预测算子对输入原始炮集地震记录进行通常的运算,生成预测的海底反射波和多次波模型道集,所述的输入原始炮集地震记录中要保留直达波;
3)利用原始炮集和预测的海底反射波和多次波模型道集自动拾取海底反射波的旅行时;其中所述拾取海底反射波的旅行时是拾取海底反射波同相轴的初至对应的时间;
4)对拾取的海底反射波旅行时进行中值滤波:
5)计算原始炮集和预测的海底反射波的旅行时差,对多次波模型道集做时差校正;
所述的计算原始炮集和预测的海底反射波的旅行时差是利用原始炮集中确定的海底反射波旅行时减去预测的海底反射波的旅行时;
所述的时差校正逐道进行,按照多次波模型道集与原始数据的海底反射时差整道上移或下移,向上移动时地震道尾部所缺部分充零,向下移动时地震道头部所缺部分充零;
6)应用时差校正后的多次波模型道集与原始数据进行数据匹配;
7)从原始数据中减去匹配后的多次波模型道集;
8)输出压制多次波后的地震剖面。
2.根据权利要求1所述的校正预测的海底多次波大时差的方法,其特征在于步骤4)中所述的中值滤波是指,将某点为中点的窗口中的所有点的值由小到大进行排序,中间值作为该点的输出值。
3.根据权利要求1所述的校正预测的海底多次波大时差的方法,其特征在于步骤4)中所述的中值滤波窗口的长度取7-11个采样点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101167054A CN101630019B (zh) | 2008-07-16 | 2008-07-16 | 一种校正预测的海底多次波大时差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101167054A CN101630019B (zh) | 2008-07-16 | 2008-07-16 | 一种校正预测的海底多次波大时差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101630019A CN101630019A (zh) | 2010-01-20 |
CN101630019B true CN101630019B (zh) | 2011-12-07 |
Family
ID=41575191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101167054A Active CN101630019B (zh) | 2008-07-16 | 2008-07-16 | 一种校正预测的海底多次波大时差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101630019B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9086505B2 (en) * | 2011-11-21 | 2015-07-21 | Cggveritas Services Sa | Device and method for computing depth velocity variations |
CN103576198B (zh) * | 2012-08-02 | 2016-04-06 | 中国石油天然气集团公司 | 一种二维海上地震资料自由表面多次波预测方法 |
CN103105623B (zh) * | 2012-12-13 | 2013-08-21 | 东北石油大学 | 一种地震勘探中的数据波形处理方法 |
US9588244B2 (en) * | 2013-07-25 | 2017-03-07 | Chevron U.S.A. Inc. | Predicting interbed multiples in seismic data using beam decomposition |
CN103984023B (zh) * | 2014-05-15 | 2017-01-25 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN104502974A (zh) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | 一种压制多次波的组合方法及装置 |
CN104834010B (zh) * | 2015-05-13 | 2017-09-12 | 中国石油集团东方地球物理勘探有限责任公司 | 两界面之间的多次波的叠后预测方法和装置 |
CN105445793B (zh) * | 2015-11-30 | 2018-04-06 | 中国石油天然气集团公司 | 一种确定坏道数据的方法及装置 |
CN106855640B (zh) * | 2017-02-03 | 2019-05-07 | 中国石油天然气集团公司 | 一种海底电缆地震数据初至计算方法及装置 |
CN110967755B (zh) * | 2018-09-30 | 2021-09-17 | 中国石油化工股份有限公司 | 基于单道数据的能量一致性匹配方法及*** |
CN111751881A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油天然气集团有限公司 | 一种海上采集地震数据旅行时的校正方法、装置及*** |
CN111781646B (zh) * | 2020-06-29 | 2022-04-15 | 广州海洋地质调查局 | 一种保幅保真线性相干噪音压制方法及处理终端 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1073236C (zh) * | 1994-09-02 | 2001-10-17 | 埃克森美孚上游研究公司 | 处理具有多次反射噪声的地震数据的方法 |
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
-
2008
- 2008-07-16 CN CN2008101167054A patent/CN101630019B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1073236C (zh) * | 1994-09-02 | 2001-10-17 | 埃克森美孚上游研究公司 | 处理具有多次反射噪声的地震数据的方法 |
US6418380B1 (en) * | 1997-10-10 | 2002-07-09 | Compagnie Generale De Geophysique | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration |
Non-Patent Citations (2)
Title |
---|
J.J.Robert等.压制多次波的策略.《石油物探译丛》.1992,(第5期),28-34. * |
左黄金.多次波去除方法研究.《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑(月刊)》.2006,(第03期),4-40. * |
Also Published As
Publication number | Publication date |
---|---|
CN101630019A (zh) | 2010-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101630019B (zh) | 一种校正预测的海底多次波大时差的方法 | |
US11105942B2 (en) | Generative adversarial network seismic data processor | |
US10317554B2 (en) | Noise attenuation via thresholding in a transform domain | |
CN103645497B (zh) | 一种基于出射角的多分量波场分离方法 | |
US20080137480A1 (en) | Method of Building a Subsurface Velocity Model | |
MX2014003060A (es) | Sistemas y metodos de filtracion de dominio de frecuencia y discriminacion de dominio de espacio-tiempo de datos sismicos. | |
WO2014209728A1 (en) | Seismic wavefield deghosting and noise attenuation | |
CN113625337B (zh) | 一种极浅水高精度地震资料快速成像方法 | |
NO20141519A1 (no) | Fremgangsmåter og systemer for optimalisering av generering av seismiske avbildninger | |
CN102928873A (zh) | 基于四维能量聚焦的地面微地震定位方法 | |
CN101609160A (zh) | 一种抗背景噪音强的地震数据初至波自动拾取方法 | |
US10877172B2 (en) | Prediction and subtraction of multiple diffractions | |
US20160161619A1 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
CN104536045A (zh) | 一种基于子波处理的鬼波压制方法 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
MX2014001334A (es) | Sistemas y metodos para detectar ruido del oleaje en un recolector sismico. | |
Wu et al. | Side-scan sonar and sub-bottom profiler surveying | |
CN103675900A (zh) | 一种确定转换波叠前时间偏移最佳速度剖面的方法 | |
US8982664B2 (en) | Optimal source signature determination | |
CN116660996A (zh) | 一种基于深度学习的漂流式浅海局部地声参数预测方法 | |
CN105372644B (zh) | 一种基于动态重修正的自适应波束形成方法及*** | |
Feng et al. | Data-domain wave equation reflection traveltime tomography | |
WO2017024112A1 (en) | Seismic wavefield deghosting | |
Wu et al. | Offset VSP P wave and converted S wave salt proximity migration study |
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 |