CN116088054A - 一种基于成像道集的剖面同相轴判别方法 - Google Patents
一种基于成像道集的剖面同相轴判别方法 Download PDFInfo
- Publication number
- CN116088054A CN116088054A CN202310206536.8A CN202310206536A CN116088054A CN 116088054 A CN116088054 A CN 116088054A CN 202310206536 A CN202310206536 A CN 202310206536A CN 116088054 A CN116088054 A CN 116088054A
- Authority
- CN
- China
- Prior art keywords
- phase axis
- profile
- slope
- section
- imaging
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Abstract
本发明涉及一种基于成像道集的剖面同相轴判别方法,属于地震数据处理领域,本发明方法第一步基于给定的偏移速度在叠前时间偏移剖面上手动拾取出同相轴,第二步利用克希霍夫积分反偏移方法将所拾同相轴从叠前时间偏移剖面转换到成像道集中,第三步在成像道集中计算出所拾同相轴的斜率,第四步通过比较给定斜率阈值和成像道集中所拾同相轴的斜率判断其为有效信息还是干扰波信息。本发明方法快速、精确地实现了叠前时间偏移剖面同相轴的判断,显著降低了实际问题的复杂程度,使叠前时间偏移剖面同相轴的判别能够用于实际地震资料处理。
Description
技术领域
本发明属于地震数据处理领域,具体的涉及一种基于成像道集的剖面同相轴判别方法。
背景技术
由于地层当中存在强反射界面,地震波在界面之间经过多次反射会形成多次波,尽管使用多次波压制方法可以将其去除,但仍有部分多次波会与一次波共同存在于叠前时间偏移剖面中,形成“真假难辨”的剖面同相轴,这会严重影响剖面质量,甚至会误导后期地震资料的解释,因此叠前时间偏移剖面同相轴判别是地震资料数据处理的重要环节。
发明内容
本发明要解决的技术问题在于提供一种基于成像道集的剖面同相轴判别方法,本发明方法保证在较小的计算量下对成像道集域的叠前时间偏移剖面同相轴实现快速精确的判断。总共分为四步实现对剖面同相轴的判别,第一步基于给定的偏移速度在叠前时间偏移剖面上手动拾取出同相轴,第二步利用克希霍夫积分反偏移方法将所拾同相轴从叠前时间偏移剖面转换到成像道集中,第三步在成像道集中计算出所拾同相轴的斜率,第四步通过比较给定斜率阈值和成像道集中所拾同相轴的斜率判断其为有效信息还是干扰波信息。
本发明采取以下技术方案:
一种基于成像道集的剖面同相轴判别方法,采用克希霍夫积分反偏移方法实现叠前时间偏移剖面同相轴判别过程,具体步骤如下:
第一步、地震数据经过处理后得到的叠前时间偏移剖面,基于给定的偏移速度在叠前时间偏移剖面上,手动拾取出剖面同相轴;
第二步、将叠前时间偏移剖面上拾取的同相轴进行克希霍夫积分反偏移处理:
(1)
式中,
x为水平坐标,
t为时间,
I ma
(x,t)表示经过反偏移得到的成像道集中的同相轴,
sec(
x,t)为叠前时间偏移剖面上拾取的同相轴,
ds为积分的面元,
W D (
x,
t)为克希霍夫积分的权重因子,其计算式为:
(2)
式中,
r 0、
r分别为入射射线、绕射射线的传播距离,
θ 0与
θ分别为
r 0、
r与面元法向矢量
n的夹角,
v(
x,t)为偏移速度;
第三步、通过式(1)、(2)将叠前时间偏移剖面上拾取的同相轴转换到成像道集中,通过式(3)计算得到其在成像道集中的斜率:
S Ima =(
t max -
t min )/ (
x max -
x min ) (3)
式中,
S Ima 为式(2)得到的
I ma
(x,t)的斜率,
x min 和
x max 分别为式(2)得到的
I ma
(x,t)的最小和最大
x坐标,
t min 和
t max 分别为式(2)得到的
I ma
(x,t)的最小和最大
t坐标;
第四步、由第三步得到成像道集中同相轴的斜率
S Ima ,基于给定的斜率阈值
S threshold 来判断其在叠前时间偏移剖面中是否为有效信息,若
S Ima <=
S threshold ,所拾的同相轴
sec(
x,t)为有效信息,否则为干扰波信息。
进一步,所述第四步中的
S threshold 取值为0.55。
本发明的原理:由于底下的地质构造形态多变,因此地震剖面上的同相轴难以有一个科学合理的判别准则。考虑到地震成像剖面多是由克希霍夫积分叠前时间偏移成像技术得到,在成像的过程中应用的是有效波的速度场,因此在成像道集中,有效信息的同相轴应该近似水平。而由于多次波的速度和有效信息的速度差异较大,其在成像道集中的同相轴应该具有一定斜率。因此将剖面中的同相轴映射至成像道集域即可通过同相轴的斜率判断其是否为有效信息。
本发明与现有技术相比的有益效果:
本发明提出了一种剖面同相轴判别方法,基于成像道集的剖面同相轴判别快速、精确地实现了叠前时间偏移剖面同相轴的判断,而通过克希霍夫积分反偏移创建的叠前时间偏移剖面所拾同相轴的成像道集显著降低了实际问题的复杂程度,使叠前时间偏移剖面同相轴的判别能够用于实际地震资料处理。
附图说明
图1基于叠前时间偏移剖面拾取的一次波同相轴;
图2使用克希霍夫积分反偏移得到的图1所拾同相轴的成像道集;
图3基于叠前时间偏移剖面拾取的多次波同相轴;
图4使用克希霍夫积分反偏移得到的图3所拾同相轴的成像道集。
具体实施方式
面通过实施例结合附图来对本发明的技术方案做进一步解释,但本发明的保护范围不受实施例任何形式上的限制。
本发明提出的剖面同相轴判别方法,具体的实施过程主要分为以下四步:
(1)基于给定的偏移速度和输入的叠前时间偏移剖面,手动拾取剖面同相轴;(2)使用克希霍夫积分反偏移算法生成所拾同相轴的对应成像道集;(3)在成像道集中计算所拾同相轴的斜率;(4)比较成像道集中所拾同相轴的斜率和给定斜率的阈值,判断其为有效信息还是干扰波信息。
下面将详细阐述本发明的具体实施过程:
图1、图3均为我国黄海某工区的地震数据经过处理后得到的叠前时间偏移剖面,图1为基于叠前时间偏移剖面拾取的一次波同相轴(有效信息),图3为基于叠前时间偏移剖面拾取的多次波同相轴(干扰波信息)。
(1)基于叠前时间偏移剖面(如图1、图3所示),手动拾取出剖面上的同相轴,其结果分别如图1箭头指向的横线、图3箭头指向的横线;
(2)将叠前时间偏移剖面上拾取的同相轴进行克希霍夫积分反偏移处理:
(1)
式中,
x为水平坐标,
t为时间,
I ma
(x,t)表示经过反偏移得到的成像道集中的同相轴,
sec(
x,t)为叠前时间偏移剖面上拾取的同相轴,
ds为积分的面元,
W D (
x,
t)为克希霍夫积分的权重因子,其计算式为:
(2)
式中,
r 0、
r分别为入射射线、绕射射线的传播距离,
θ 0与
θ分别为
r 0、
r与面元法向矢量
n的夹角,
v(
x,t)为偏移速度。
将图1、图3所示的叠前时间偏移剖面上拾取的同相轴进行克希霍夫积分反偏移处理,即得到其对应的成像道集(如图2、图4所示);
(3)通过式(1)、(2)将叠前时间偏移剖面上拾取的同相轴转换到成像道集中,通过式(3)计算得到其在成像道集中的斜率:
S Ima =(
t max -
t min )/ (
x max -
x min ) (3)
式中,
S Ima 为式(2)得到的
I ma
(x,t)的斜率,
x min 和
x max 分别为式(2)得到的
I ma
(x,t)的最小和最大
x坐标,
t min 和
t max 分别为式(2)得到的
I ma
(x,t)的最小和最大
t坐标;
计算图2、图4所示的成像道集中的同相轴的斜率。得到图2所示的成像道集中的同相轴的
t max 为1310ms,
t min 为1300ms,
x max 为1230m,
x min 为1210m,则同相轴斜率为0.5;图4所示的成像道集中的同相轴的的
t max 为1400ms,
t min 为1300ms,
x max 为540m,
x min 为516m,则同相轴斜率为4.16。
(4)由第三步得到成像道集中同相轴的斜率
S Ima ,基于给定的斜率阈值
S threshold 来判断其在叠前时间偏移剖面中是否为有效信息,若
S Ima <=
S threshold ,所拾的同相轴
sec(
x,t)为有效信息,否则为干扰波信息。
将图2、图4所示成像道集中的同相轴斜率的计算结果与给定斜率的阈值(取0.55)相比较,显然图2的同相轴斜率值小于给定斜率的阈值,图4的同相轴斜率值大于给定斜率的阈值。即可判断:图1拾取的剖面同相轴为有效信息,图3拾取的剖面同相轴为干扰波信息。该判别结果与实际资料一致,即验证了本发明的有效性。
由此可见,基于成像道集的剖面同相轴判别快速、精确地实现了叠前时间偏移剖面同相轴的判断,而通过克希霍夫积分反偏移创建的叠前时间偏移剖面所拾同相轴的成像道集显著降低了实际问题的复杂程度,使叠前时间偏移剖面同相轴的判别能够用于实际地震资料处理,是一种较为理想的剖面同相轴判别方法。
Claims (2)
1.一种基于成像道集的剖面同相轴判别方法,其特征在于,所述方法采用克希霍夫积分反偏移方法实现叠前时间偏移剖面同相轴判别过程,具体步骤如下:
第一步、地震数据经过处理后得到的叠前时间偏移剖面,基于给定的偏移速度在叠前时间偏移剖面上,手动拾取出剖面同相轴;
第二步、将叠前时间偏移剖面上拾取的同相轴进行克希霍夫积分反偏移处理:
(1)
式中,x为水平坐标,t为时间,I ma (x,t)表示经过反偏移得到的成像道集中的同相轴,sec(x,t)为叠前时间偏移剖面上拾取的同相轴,ds为积分的面元,W D (x,t)为克希霍夫积分的权重因子,其计算式为:
(2)
式中,r 0、r分别为入射射线、绕射射线的传播距离,θ 0与θ分别为r 0、r与面元法向矢量n的夹角,v(x,t)为偏移速度;
第三步、通过式(1)、(2)将叠前时间偏移剖面上拾取的同相轴转换到成像道集中,通过式(3)计算得到其在成像道集中的斜率:
S Ima =( t max -t min )/ ( x max -x min ) (3)
式中,S Ima 为式(2)得到的I ma (x,t)的斜率,x min 和x max 分别为式(2)得到的I ma (x,t)的最小和最大x坐标,t min 和t max 分别为式(2)得到的I ma (x,t)的最小和最大t坐标;
第四步、由第三步得到成像道集中同相轴的斜率S Ima ,基于给定的斜率阈值S threshold 来判断其在叠前时间偏移剖面中是否为有效信息,若S Ima ≤S threshold ,所拾的同相轴sec(x,t)为有效信息,否则为干扰波信息。
2.根据权利要求1所述的一种基于成像道集的剖面同相轴判别方法,其特征在于,所述第四步中的S threshold 取值为0.55。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310206536.8A CN116088054B (zh) | 2023-03-07 | 2023-03-07 | 一种基于成像道集的剖面同相轴判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310206536.8A CN116088054B (zh) | 2023-03-07 | 2023-03-07 | 一种基于成像道集的剖面同相轴判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116088054A true CN116088054A (zh) | 2023-05-09 |
CN116088054B CN116088054B (zh) | 2023-06-13 |
Family
ID=86204638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310206536.8A Active CN116088054B (zh) | 2023-03-07 | 2023-03-07 | 一种基于成像道集的剖面同相轴判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116088054B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120290210A1 (en) * | 2011-05-06 | 2012-11-15 | Philip Cheung | Systems and methods for investigating a formation surrounding a borehole |
CN103777242A (zh) * | 2012-10-24 | 2014-05-07 | 中国石油化工股份有限公司 | 一种深度聚焦和道集同相轴拉平联合的速度判别方法 |
CN104280766A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | 一种利用局部数据同相轴斜率的直接偏移方法 |
US20150078124A1 (en) * | 2013-09-19 | 2015-03-19 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
CN105223607A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 基于反稳相滤波的绕射波成像方法 |
WO2016041185A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种高效叠前时间偏移速度分析方法 |
CN109283579A (zh) * | 2018-11-08 | 2019-01-29 | 中海油田服务股份有限公司 | 一种基于克希霍夫积分反偏移的层间多次波预测方法 |
CN111856577A (zh) * | 2020-07-28 | 2020-10-30 | 中国石油大学(北京) | 一种降低逆时偏移地表炮检距道集计算量的方法 |
CN113031063A (zh) * | 2021-04-09 | 2021-06-25 | 中国海洋大学 | 一种基于成像道集相关加权的逆时偏移成像方法 |
CN114859410A (zh) * | 2022-04-22 | 2022-08-05 | 西安石油大学 | 基于差分进化算法的b-cosfire滤波器的地震同相轴提取方法 |
-
2023
- 2023-03-07 CN CN202310206536.8A patent/CN116088054B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120290210A1 (en) * | 2011-05-06 | 2012-11-15 | Philip Cheung | Systems and methods for investigating a formation surrounding a borehole |
CN103777242A (zh) * | 2012-10-24 | 2014-05-07 | 中国石油化工股份有限公司 | 一种深度聚焦和道集同相轴拉平联合的速度判别方法 |
CN104280766A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | 一种利用局部数据同相轴斜率的直接偏移方法 |
US20150078124A1 (en) * | 2013-09-19 | 2015-03-19 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
CN105223607A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 基于反稳相滤波的绕射波成像方法 |
WO2016041185A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种高效叠前时间偏移速度分析方法 |
CN109283579A (zh) * | 2018-11-08 | 2019-01-29 | 中海油田服务股份有限公司 | 一种基于克希霍夫积分反偏移的层间多次波预测方法 |
CN111856577A (zh) * | 2020-07-28 | 2020-10-30 | 中国石油大学(北京) | 一种降低逆时偏移地表炮检距道集计算量的方法 |
CN113031063A (zh) * | 2021-04-09 | 2021-06-25 | 中国海洋大学 | 一种基于成像道集相关加权的逆时偏移成像方法 |
CN114859410A (zh) * | 2022-04-22 | 2022-08-05 | 西安石油大学 | 基于差分进化算法的b-cosfire滤波器的地震同相轴提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116088054B (zh) | 2023-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110058305A (zh) | 一种基于卷积神经网络的das地震数据降噪方法 | |
CN101630015B (zh) | 一种提高初至波拾取精度和效率的方法 | |
CN111666944B (zh) | 一种红外弱小目标检测方法及装置 | |
CN106646598A (zh) | 一种fast‑aic法微地震信号拾取方法 | |
CN108983284A (zh) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 | |
CN115097442B (zh) | 基于毫米波雷达的水面环境地图构建方法 | |
CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
CN107219554A (zh) | 陆地地震资料的剩余静校正量的自动获取方法 | |
CN105629300B (zh) | 提高复杂构造偏移数据信噪比的方法 | |
CN115905805A (zh) | 一种基于全局信息判别gan的das数据多尺度噪声消减方法 | |
CN107991707A (zh) | 一种基于shear let域内峰度特性的井中微地震初至波拾取方法 | |
CN108983287A (zh) | 一种基于凸集投影算法的曲波变换反假频地震数据重建方法 | |
CN116088054B (zh) | 一种基于成像道集的剖面同相轴判别方法 | |
CN110389377A (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
CN109581481B (zh) | 一种便携式高频可控震源地震信号谐波干扰消除方法 | |
CN102338885B (zh) | 三分量vsp资料初至时间自动拾取方法 | |
CN110850474A (zh) | 一种海洋地震资料震源鬼波压制方法和系统 | |
CN105676169B (zh) | 一种脉冲声目标的定向方法及装置 | |
CN112558159A (zh) | 一种声波测井初至拾取方法 | |
JP3263752B2 (ja) | 埋設物探知に用いる受信信号の処理方法 | |
CN111898478B (zh) | 一种基于多尺度形态学的初至自动拾取方法 | |
CN104280766B (zh) | 一种利用局部数据同相轴斜率的直接偏移方法 | |
JPH0572332A (ja) | 不可視物体探査方法 | |
CN111624668A (zh) | 一种用于频分电法的谐波校正方法 | |
CN111650645A (zh) | 一种变偏移距vsp弯线校正处理方法和装置 |
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 |