CN109283579A - 一种基于克希霍夫积分反偏移的层间多次波预测方法 - Google Patents
一种基于克希霍夫积分反偏移的层间多次波预测方法 Download PDFInfo
- Publication number
- CN109283579A CN109283579A CN201811323096.XA CN201811323096A CN109283579A CN 109283579 A CN109283579 A CN 109283579A CN 201811323096 A CN201811323096 A CN 201811323096A CN 109283579 A CN109283579 A CN 109283579A
- Authority
- CN
- China
- Prior art keywords
- formula
- interbed multiple
- wave
- migration
- indicates
- 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
Links
- 230000005012 migration Effects 0.000 title claims abstract description 47
- 238000013508 migration Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000003384 imaging method Methods 0.000 description 4
- 239000011229 interlayer Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
技术领域
本发明属于地震资料数据处理领域,是一种既可减少计算量又可保证精度的基于克希霍夫积分反偏移的层间多次波预测方法。
背景技术
由于地层当中存在强反射界面,地震波在界面之间经过多次反射形成层间多次波,若不加以去除将严重影响成像剖面的质量,甚至会误导后期地震资料的解释,因此预测与衰减层间多次波是地震资料数据处理的关键环节。
多次波压制主要是根据多次波与一次波性质上的差异,应用一定的数学处理方法分离或预测出多次波,并在此基础上衰减多次波。当前压制多次波的方法主要有两种:一种是滤波法,其主要是依据多次波与一次波视速度或周期性的差异,利用数学变换将多次波分离,达到压制多次波的效果,但是在多次波与一次波视速度相近的情况下,该方法无法达到理想的压制效果;另一种常用的多次波压制方法是预测减去法,该方法是依据多次波与一次波的关系预测出多次波,并在相应的炮集域或CMP域中减去相应的预测结果,以达到压制多次波的目的。
常规的多次波压制方法对于层间多次波无法起到良好的效果,其原因在于层间多次波视速度与一次波相近,并且周期性差异也不大,因此用常规的滤波方法难以识别与剔除层间多次波,而采用预测减去法计算量大,成本高,限制了其在实际资料处理中的应用。
发明内容
本发明要解决的技术问题在于提供一种基于克希霍夫积分反偏移的层间多次波预测方法,所述方法保证在较小的计算量下对多次波实现最有效的压制。总共分为两步实现对层间多次波的预测,第一步在成像剖面域预测一次波对应的多次波,第二步利用克希霍夫积分反偏移方法将其变换到炮集域。
本发明采取以下技术方案:
一种基于克希霍夫积分反偏移的层间多次波预测方法,其特征在于它具体包括以下步骤:
(1)给定一次波叠前时间剖面IMT,给定偏移速度v(x,t),利用式(1)对一次波叠前时间偏移记录进行“时——深”转换,得到叠前深度偏移剖面IMD;
式中,IMD表示叠前深度偏移剖面,IMT表示叠前时间偏移剖面,v(x,t)为地震波速,x表示横向距离,t表示时间,z表示深度;
(2)输入叠前深度偏移剖面IMD,利用式(2)得到波数域的层间多次波预测结果IMKm;
式中,x表示横向距离,kz为波数,z′1、z′2、z′3为深度,IMKm为波数域的层间多次波预测结果;
(3)利用式(3)对IMKm作反傅里叶变换,得到深度域层间多次波剖面IMDm;
式中x表示横向距离,z表示深度,kz为波数
(4)利用式(4)对深度域层间多次波剖面IMDm做“深——时”转换,最终得到层间多次波的叠前时间偏移剖面IMTm;
式中v(x,t)为地震波速,x表示横向距离,z表示深度,t表示时间;
(5)依据层间多次波的叠前时间偏移剖面IMTm,利用式(5)基于克希霍夫积分反偏移算法生成层间多次波炮集记录IDTm;
式中,IDTm为多次波炮集记录,x表示横向距离,t表示时间,WD(x,t)为积分权重因子,其可由公式(6)求得;
式中,v(x,t)为地震波速,r0、r分别为入射线、绕射线的传播距离,θ0与θ分别为入射线r0、绕射线r与面元法向矢量n的夹角。
本发明与现有技术相比的有益效果:
本发明提出了一种基于克希霍夫积分反偏移的层间多次波预测方法,该方法基于叠前时间偏移剖面的层间多次波预测有效减少了多次波预测的计算量,而通过克希霍夫反偏移创建层间多次波炮集记录其显著降低了实际问题的复杂程度,使层间多次波预测能够用于实际地震资料处理。
附图说明
图1基于一次波速度偏移成像的叠前时间偏移剖面;
图2仅包含层间多次波的叠前时间偏移剖面IMTm;
图3预测的层间多次波炮集记录与偏移速度谱示例(a)层间多次波炮集记录示例(炮号为110);(b)层间多次波记录的偏移速度谱示例(x=8250m);
图4原始记录的层间多次波特征分析(a)原始炮集记录示例(炮号为110);(b)偏移速度谱示例(x=8250m);
图5衰减多次波后的炮集记录与去除的层间多次波干扰示例(a)衰减多次波后的炮集记录(炮号为110);(b)去除的层间多次波干扰(炮号为110);
图6衰减层间多次波后炮集记录的叠前时间偏移剖面。
具体实施方式
下面通过实施例来对本发明的技术方案作进一步研究,但本发明的保护范围不受实施例任何形式上的限制。
本发明提出的一种基于克希霍夫积分反偏移的层间多次波预测方法,具体的实施过程主要分为以下两步:
(1)在成像剖面域预测一次波对应的多次波;(2)利用克希霍夫积分反偏移方法将其变换到炮集域。
采用海上二维拖缆观测方式模拟含有强层间多次波的炮集记录,其震源与检波器均位于海面(即深度为0)。模拟的数据集共含有200炮,炮间距为50m;每炮含有118道,其最小偏移距与道间距均为25m;各炮记录的时间采样率为4ms,记录长度达到5116ms。
下面将详细阐述本发明的具体实施过程:
(1)输入图1所示基于一次波速度偏移成像的原始叠前时间偏移剖面,利用式(1)对一次波叠前时间偏移记录进行“时——深”转换,得到叠前深度偏移剖面IMD;
(2)输入叠前深度偏移剖面IMD,利用公式(2)逆散射级数法进行层间多次波预测,得到波数域的层间多次波预测结果IMKm;
(3)利用式(3)对IMKm作反傅里叶变换,得到深度域层间多次波剖面IMDm;
(4)利用式(4)对深度域层间多次波剖面IMDm做“深——时”转换,最终得到层间多次波的叠前时间偏移剖面IMTm,所得结果如图2所示。通过对比可知,图2所示的多次波剖面中包含了丰富的层间多次波信息,其主要同相轴的旅行时与原始剖面中的基本一致,进一步验证了基于叠前时间偏移剖面进行层间多次波预测的有效性;
(5)输入层间多次波的叠前时间偏移剖面IMTm(见图2),利用时域克希霍夫积分反偏移方法生成层间多次波炮集记录。图3(a)与(b)分别给出了预测的层间多次波炮集记录及其偏移速度谱,为使深部反射清晰可见,对所示记录进行了球面扩散补偿处理。与图4所示的原始炮集记录相比,预测的记录中含有丰富的多次波信息,且强层间多次波同相轴的旅行时与原始记录中的基本一致;
图5显示了衰减多次波后的炮集记录以及减去的层间多次波干扰,经过对比可知,原始炮集记录(见图4)中的多次波同相轴受到了明显压制。对衰减多次波后的炮集记录进行克希霍夫叠前时间偏移处理,所得偏移剖面如图6所示,其中层间多次波同相轴已被完全去除,剖面的信噪比得到显著提升。
由此可见,基于叠前时间偏移剖面的层间多次波预测有效减少了多次波预测的计算量,而通过克希霍夫反偏移创建层间多次波炮集记录其显著降低了实际问题的复杂程度,使层间多次波预测能够用于实际地震资料处理,是一种较为理想的层间多次波预测方法。
Claims (1)
1.一种基于克希霍夫积分反偏移的层间多次波预测方法,其特征在于它具体包括以下步骤:
(1)给定一次波叠前时间剖面IMT,给定偏移速度v(x,t),利用式(1)对一次波叠前时间偏移记录进行“时——深”转换,得到叠前深度偏移剖面IMD;
式中,IMD表示叠前深度偏移剖面,IMT表示叠前时间偏移剖面,v(x,t)为地震波速,x表示横向距离,t表示时间,z表示深度;
(2)输入叠前深度偏移剖面IMD,利用式(2)得到波数域的层间多次波预测结果IMKm;
式中,x表示横向距离,kz为波数,z′1、z′2、z′3为深度,IMKm为波数域的层间多次波预测结果;
(3)利用式(3)对IMKm作反傅里叶变换,得到深度域层间多次波剖面IMDm;
式中x表示横向距离,z表示深度,kz为波数
(4)利用式(4)对深度域层间多次波剖面IMDm做“深——时”转换,最终得到层间多次波的叠前时间偏移剖面IMTm;
式中v(x,t)为地震波速,x表示横向距离,z表示深度,t表示时间;
(5)依据层间多次波的叠前时间偏移剖面IMTm,利用式(5)基于克希霍夫积分反偏移算法生成层间多次波炮集记录IDTm;
式中,IDTm为多次波炮集记录,x表示横向距离,t表示时间,WD(x,t)为积分权重因子,其可由公式(6)求得;
式中,v(x,t)为地震波速,r0、r分别为入射线、绕射线的传播距离,θ0与θ分别为入射线r0、绕射线r与面元法向矢量n的夹角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811323096.XA CN109283579B (zh) | 2018-11-08 | 2018-11-08 | 一种基于克希霍夫积分反偏移的层间多次波预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811323096.XA CN109283579B (zh) | 2018-11-08 | 2018-11-08 | 一种基于克希霍夫积分反偏移的层间多次波预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109283579A true CN109283579A (zh) | 2019-01-29 |
CN109283579B CN109283579B (zh) | 2019-07-19 |
Family
ID=65174770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811323096.XA Active CN109283579B (zh) | 2018-11-08 | 2018-11-08 | 一种基于克希霍夫积分反偏移的层间多次波预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109283579B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112578453A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 多次波压制方法及系统 |
CN114114419A (zh) * | 2020-08-31 | 2022-03-01 | 中国石油化工股份有限公司 | 层间多次波预测、压制方法和提高地震资料成像品质方法 |
CN116088054A (zh) * | 2023-03-07 | 2023-05-09 | 中国海洋大学 | 一种基于成像道集的剖面同相轴判别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1094339A2 (en) * | 1999-10-18 | 2001-04-25 | Petroleo Brasileiro S.A. Petrobras | A method for the supression of multiple reflections from marine seimsic data |
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN105510976A (zh) * | 2016-01-11 | 2016-04-20 | 中国海洋大学 | 一种多次波组合自适应衰减方法 |
CN105572735A (zh) * | 2014-10-17 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种提高叠前深度偏移成像精度的方法及装置 |
CN105676291A (zh) * | 2016-01-11 | 2016-06-15 | 中国海洋大学 | 一种基于同相轴优化追踪的多次波匹配衰减方法 |
-
2018
- 2018-11-08 CN CN201811323096.XA patent/CN109283579B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1094339A2 (en) * | 1999-10-18 | 2001-04-25 | Petroleo Brasileiro S.A. Petrobras | A method for the supression of multiple reflections from marine seimsic data |
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN105572735A (zh) * | 2014-10-17 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种提高叠前深度偏移成像精度的方法及装置 |
CN105510976A (zh) * | 2016-01-11 | 2016-04-20 | 中国海洋大学 | 一种多次波组合自适应衰减方法 |
CN105676291A (zh) * | 2016-01-11 | 2016-06-15 | 中国海洋大学 | 一种基于同相轴优化追踪的多次波匹配衰减方法 |
Non-Patent Citations (2)
Title |
---|
宋建国 等: "基于克希霍夫衍射理论的水层多次波预测方法", 《地球物理学进展》 * |
谭军 等: "基于同相轴优化追踪的多次波匹配衰减方法", 《中国石油大学学报(自然科学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112578453A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 多次波压制方法及系统 |
CN114114419A (zh) * | 2020-08-31 | 2022-03-01 | 中国石油化工股份有限公司 | 层间多次波预测、压制方法和提高地震资料成像品质方法 |
CN114114419B (zh) * | 2020-08-31 | 2024-04-09 | 中国石油化工股份有限公司 | 层间多次波预测、压制方法和提高地震资料成像品质方法 |
CN116088054A (zh) * | 2023-03-07 | 2023-05-09 | 中国海洋大学 | 一种基于成像道集的剖面同相轴判别方法 |
CN116088054B (zh) * | 2023-03-07 | 2023-06-13 | 中国海洋大学 | 一种基于成像道集的剖面同相轴判别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109283579B (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109283579B (zh) | 一种基于克希霍夫积分反偏移的层间多次波预测方法 | |
CN102590859B (zh) | 垂向各向异性介质准p波方程逆时偏移方法 | |
US20160320509A1 (en) | Suppressing near-surface scattered surface waves | |
CN103926622B (zh) | 一种基于l1范数多道匹配滤波压制多次波的方法 | |
CA2357604A1 (en) | Method for identifying and removing multiples from seismic reflection data | |
CN110058305A (zh) | 一种基于卷积神经网络的das地震数据降噪方法 | |
Xu et al. | Elastic model low-to intermediate-wavenumber inversion using reflection traveltime and waveform of multicomponent seismic data | |
CN106896409B (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
Yang et al. | An efficient and effective common reflection surface stacking approach using local similarity and plane-wave flattening | |
CN107894613A (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
CN107144880A (zh) | 一种地震波波场分离方法 | |
CN110658557B (zh) | 基于生成对抗网络的地震数据面波压制方法 | |
CN107144879A (zh) | 一种基于自适应滤波与小波变换结合的地震波降噪方法 | |
CN106932824A (zh) | 陆地地震勘探资料的降维自适应层间多次波压制方法 | |
CN106707342B (zh) | 共炮点道集多级面波压制方法和装置 | |
US6112155A (en) | Multiple suppression in geophysical data | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN104570114B (zh) | 一种基于波场分解的逆时偏移噪音压制方法 | |
CN111257938A (zh) | 基于小波互相关时移地震虚拟震源波场重构方法和系统 | |
CN113031072B (zh) | 虚同相轴层间的多次波压制方法、装置及设备 | |
CN115236733A (zh) | 基于深度学习的das-vsp数据背景噪声压制方法 | |
Zhang et al. | Seismic random noise attenuation using modified wavelet thresholding | |
Li et al. | Progressive multitask learning for high-resolution prediction of reservoir elastic parameters | |
CN112946733A (zh) | 海上立体观测系统联合压制多次波多缆处理方法及系统 | |
Guo et al. | Denoising and wavefield separation method for DAS VSP via deep learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |