CN112444881B - 一种鬼波压制方法 - Google Patents
一种鬼波压制方法 Download PDFInfo
- Publication number
- CN112444881B CN112444881B CN201910818531.4A CN201910818531A CN112444881B CN 112444881 B CN112444881 B CN 112444881B CN 201910818531 A CN201910818531 A CN 201910818531A CN 112444881 B CN112444881 B CN 112444881B
- Authority
- CN
- China
- Prior art keywords
- ghost
- wave field
- wave
- free
- ghost 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.)
- Active
Links
- 238000003825 pressing Methods 0.000 title claims description 23
- 230000007613 environmental effect Effects 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000014509 gene expression Effects 0.000 claims description 52
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000005056 compaction Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 description 10
- 230000001629 suppression Effects 0.000 description 9
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 5
- 125000004122 cyclic group Chemical group 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000004613 tight binding model Methods 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种鬼波压制方法,通过原始数据得到环境因素影响下的初始无鬼波波场,由于初始无鬼波波场受到环境因素的影响和制约,初始无鬼波波场中鬼波却除得并不准确,本发明从初始无鬼波波场开始反复叠代,通过构建约束对象并判断约束对象是否小于或等于预设值,判断对象大于预设值则生成下一级无鬼波波场,直到构建的判断对象小于预设值为止,此时,将待处理无鬼波波场作为目标无鬼波波场,在最终得到的目标无鬼波波场中,鬼波已基本被消除。该方法从初始无鬼波波场开始反复叠代,消除了环境因素的影响和制约,能够快速准确地消除原始地震记录中的鬼波,从而提高了低频能量,很好地保证地数据的拓频效果,提高了地震资料的分辨率和精度。
Description
技术领域
本发明涉及海洋地震宽频采集数据处理领域,尤其涉及一种鬼波压制方法。
背景技术
海洋地震勘探中,检波器除接收到海底一次反射波外,还接收到与自由表面相关的鬼波。鬼波作为一种特殊的噪声,影响海底一次反射的波形和振幅,甚至产生虚假同相轴,降低地震勘探精度,鬼波造成的陷波现象和低频端能量损失导致地震信号有效频带宽度不足。
对于常规拖缆鬼波压制,多种鬼波压制技术被提出,如滤波方法,根据实际采集信号是一次波与鬼波算子褶积的假设,通过反滤波算法实现鬼波压制,但由于鬼波延迟时间难以精确计算,导致滤波方法计算精度较低;逆散射级数法(ISS)和基于格林函数理论的鬼波压制算法从地震信号波动理论出发,基于高阶近似和波场散射理论压制鬼波,但是该技术计算量大。但是,常规采集资料原始频宽宽度不足,经过鬼波压制处理后的地震资料频带宽度往往也难以达到宽频信号的要求,且常规采集资料处理存在计算效率、稳定性和精度难以兼顾的问题。
因此,上下缆、双检、犁式电缆、斜缆等宽频采集技术逐渐得到推广和应用,其中斜缆采集技术作为一种拥有连续变化深度检波器的采集技术,可以获得连续变化的频带宽度,相比常规采集,斜缆采集技术可以获得宽频带的海洋地震原始资料,数据说服力强,尤其受到地球物理学者的重视。
对于斜缆采集资料鬼波压制,地球物理学者们也提出了诸多处理技术,但传统的的斜缆鬼波压制技术对地下假设信息依赖性强,稳定性差,受采集环境等制约因素明显,由于受环境制约,斜缆上的检波器可能与预设位置存在偏差,并不刚好处于预设位置处,数据拓频效果难以保证,不利于技术普及推广,制约海洋海洋宽频地震勘探的发展。
发明内容
为解决现有技术中受环境因素制约,鬼波压制的数据拓频效果难以保证的技术问题,本发明提供一种鬼波压制方法,具体方案如下:
一种鬼波压制方法,其特征在于,包括以下步骤:
S1:获取原始数据,根据所述原始数据获取环境因素影响下的初始无鬼波波场,以所述初始无鬼波波场作为待处理无鬼波波场;
S2:根据所述待处理无鬼波波场构建判断对象,并判断所述判断对象是否小于或等于预设值;
若是,则将所述待处理无鬼波波场作为目标无鬼波波场并结束;
若否,则根据所述待处理无鬼波波场获取所述下一级无鬼波波场;
S3:更新所述待处理无鬼波波场,将所述下一级无鬼波波场作为所述待处理无鬼波波场并返回步骤S2。
进一步的,在步骤S2中,根据所述待处理无鬼波波场构建判断对象的操作步骤包括:
所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取与所述待处理无鬼波波场对应的鬼波波场;
根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间;
根据所述待处理无鬼波波场、与所述待处理无鬼波波场对应的鬼波波场以及与所述待处理无鬼波波场对应的鬼波波场延时时间构建所述判断对象。
进一步的,所述判断对象的表达式为:
其中:
Pupi(ω,x,y)为待处理无鬼波波场在频率域上的表达式;
Pdi(ω,x,y)为与所述待处理无鬼波波场对应的鬼波波场在频率域上的表达式;
ω为频率;
x,y为用于获取所述原始地震记录的检波器的坐标;
j为虚数单位;
e为自然常数;
Ti为与所述待处理无鬼波波场对应的鬼波波场延时时间;
O为判断对象。
进一步的,所述预设值的取值范围为大于0且小于或等于0.1。
进一步的,步骤S2中,根据所述待处理无鬼波波场获取所述下一级无鬼波波场的操作步骤包括:
所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间;
根据所述鬼波延时时间获取下一级无鬼波波场对应的第一鬼波算子;
根据下一级无鬼波波场对应的第一鬼波算子和原始地震记录获取下一级无鬼波波场。
进一步的,根据所述鬼波延时时间获取下一级无鬼波波场对应的第一鬼波算子的操作步骤中,下一级无鬼波波场对应的第一鬼波算子满足以下表达式:
其中:
Ti为与所述待处理无鬼波波场对应的鬼波波场延时时间;
j为虚数单位;
ω为频率;
e为自然常数;
Gp(i+1)为下一级无鬼波波场对应的第一鬼波算子。
进一步的,根据下一级无鬼波波场对应的第一鬼波算子和原始地震记录获取下一级无鬼波波场的操作步骤中,下一级无鬼波波场满足以下表达式:
其中:
Pup(i+1)(ω,x,y)为下一级无鬼波波场在频率域上的表达式;
P(ω,x,y)为原始地震记录在频率域上的表达式;
ω为频率;
x,y为用于获取所述原始地震记录的检波器的坐标;
为Gp(i+1)的倒数。
进一步的,步骤S1中,根据所述原始数据获取环境因素影响下的初始无鬼波波场的操作步骤包括:
所述原始数据包括原始地震记录和用于计算初始无鬼波波场对应的第二鬼波算子的计算数据;根据所述原始地震记录获取所述原始地震记录的镜像数据;根据所述计算数据获取所述初始无鬼波波场对应的第二鬼波算子,根据所述第二鬼波算子获取第一鬼波算子;
根据所述原始地震记录、所述原始地震记录的镜像数据、所述初始无鬼波波场对应的第一鬼波算子以及所述初始无鬼波波场对应的第二鬼波算子获所述取初始无鬼波波场。
进一步的,所述初始无鬼波波场满足以下表达式:
其中:
Pup1(ω,x,y)为所述初始无鬼波波场在频率域上的表达式;
P2(ω,x,y)为P(ω,x,y)的平方,P(ω,x,y)为所述原始地震记录在频率域上的表达式;
M2(ω,x,y)为M(ω,x,y)的平方,M(ω,x,y)为所述原始地震记录的镜像数据在频率域上的表达式;
为GP的平方,Gp为初始无鬼波波场对应的第一鬼波算子;
为Gm的平方,Gm为初始无鬼波波场对应的第二鬼波算子;
进一步的,所述初始无鬼波波场对应的第二鬼波算子Gm满足以下表达式:
其中:
e为自然常数;
j为虚数单位;
ω为频率;
为环境影响因子;
z为用于获取原始地震记录的检波器的深度;
c为水流速度。
进一步的,环境影响因子的取值范围为大于0且小于或等于0.1。
与现有技术相比,本发明提供一种鬼波压制方法,该鬼波压制方法通过原始数据得到环境因素影响下的初始无鬼波波场,由于初始无鬼波波场受到环境因素的影响和制约,初始无鬼波波场中鬼波却除得并不准确,本发明从初始无鬼波波场开始反复叠代,通过构建约束对象并判断约束对象是否小于或等于预设值,判断对象大于预设值则生成下一级无鬼波波场,直到构建的判断对象小于预设值为止,此时,将待处理无鬼波波场作为目标无鬼波波场,在最终得到的目标无鬼波波场中,鬼波已基本被消除。该方法从初始无鬼波波场开始反复叠代,消除了环境因素的影响和制约,能够快速准确地消除原始地震记录中的鬼波,从而提高了低频能量,很好地保证地数据的拓频效果,提高了地震资料的分辨率和精度。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:
图1为采用本发明的鬼波压制方法对模拟资料去鬼波前后单炮记录对比,包括图1a、图1b和图1c,其中图1a为鬼波压制前的数据,图1b为鬼波压制后的数据,图1c为鬼波数据;
图2为采用本发明的鬼波压制方法对模拟资料去鬼波前后频谱对比,包括图2a和图2b,其中图2a为鬼波压制前的数据,图2b为鬼波压制后的数据。
图3为采用本发明的鬼波压制方法对斜缆实际采集资料去鬼波前后单炮记录对比,包括图3a和图3b,其中图3a为鬼波压制前的数据,图3b为鬼波压制后的数据;
图4为采用本发明的鬼波压制方法对斜缆实际采集资料去鬼波前后叠加剖面对比,包括图4a和图4b,其中图4a为鬼波压制前的数据,图4b为鬼波压制后的数据;
图5为采用本发明的鬼波压制方法对斜缆实际采集资料去去鬼波前后频谱对比。
在附图中,相同的标识对象采用相同的附图标记,附图并未按实际比例绘制。
具体实施方式
下面将结合附图对本发明作进一步的说明。
本实施例提供一种鬼波压制方法,该方法包括以下步骤:
S1:获取原始数据,根据所述原始数据获取环境因素影响下的初始无鬼波波场,以所述初始无鬼波波场作为待处理无鬼波波场;
S2:根据所述待处理无鬼波波场构建判断对象,并判断所述判断对象是否小于或等于预设值;
若是,则将所述待处理无鬼波波场作为目标无鬼波波场并结束;
若否,则根据所述待处理无鬼波波场获取所述下一级无鬼波波场;
S3:更新所述待处理无鬼波波场,将所述下一级无鬼波波场作为所述待处理无鬼波波场并返回步骤S2。
与现有技术相比,该鬼波压制方法通过原始数据得到环境因素影响下的初始无鬼波波场,由于初始无鬼波波场受到环境因素的影响和制约,初始无鬼波波场中鬼波却除得并不准确,本发明从初始无鬼波波场开始反复叠代,通过构建约束对象并判断约束对象是否小于或等于预设值,判断对象大于预设值则生成下一级无鬼波波场,直到构建的判断对象小于预设值为止,此时,将待处理无鬼波波场作为目标无鬼波波场,在最终得到的目标无鬼波波场中,鬼波已基本被消除。该方法从初始无鬼波波场开始反复叠代,消除了环境因素的影响和制约,能够快速准确地消除原始地震记录中的鬼波,从而提高了低频能量,很好地保证地数据的拓频效果,提高了地震资料的分辨率和精度。
本实施例中的除后文提到的理想状态下的完全无鬼波波场外,其余无鬼波波场如初始无鬼波波场以及叠代出的各级无鬼波波场等并非完全准确地去除了鬼波,只是叠代次数越多,相应得到的无鬼波波场就越接近于理想状态下的完全无鬼波波场,也就意味着鬼波去除得越准确,最终得到的目标无鬼波波场已基本上准确地将鬼波去除。理想状态下的完全无鬼波波场则是鬼波被全部准确地去除后的波场。
本实施例的鬼波压制方法具体步骤如下。
获取原始数据,在海洋斜缆地震采集中,假设斜缆上某一检波器获取到的时间域原始地震记录为p(t,x,y),根据原始地震记录p(t,x,y)可得到该原始地震记录p(t,x,y)的镜像数据m(t,x,y),将时间域的原始地震记录p(t,x,y)和镜像数据m(t,x,y)转换到频率域,得到原始地震记录在频率域上的表达式P(ω,x,y)以及镜像数据在频率域上的表达式M(ω,x,y)。
理想状态下,假设斜缆不受环境因素的影响,斜缆的检波器处于预设位置处,检波器获取到的原始地震记录与没有环境因素影响的理想状态下的完全无鬼波波场之间的关系在频率域上的表达式为:
P(ω,x,y)=GpkPup(ω,x,y) (1)
式(1)中:
P(ω,x,y)为原始地震记录在频率域上的表达式;
Gpk为没有环境因素影响的理想状态下的完全无鬼波波场对应的第一鬼波算子;第一鬼波算子用于表征原始地震记录与相应的无鬼波波场之间的关系;
Pup(ω,x,y)为没有环境因素影响的理想状态下的完全无鬼波波场在频率域上的表达式;
ω为频率;
x,y为用于获取所述原始地震记录的检波器的坐标。
检波器获取到的原始地震记录的镜像数据与没有环境因素影响的理想状态下的完全无鬼波波场之间的关系在频率域上的表达式为:
M(ω,x,y)=GmkPup(ω,x,y) (2)
式(2)中:
M(ω,x,y)为原始地震记录的镜像数据在频率域上的表达式;
Gmk为没有环境因素影响的理想状态下的完全无鬼波波场对应的第二鬼波算子;第二鬼波算子用于表征原始地震记录的镜像数据与相应的无鬼波波场之间的关系;
其余参数含义与式(1)相同。
根据镜像原理:
式(3)中,*表示复数的共轭形式,即Gmk和Gpk为共轭复数,各参数与前式中的参数含义相同。
对于每一级无鬼波波场,相应的第一鬼波算子和第二鬼波算子为共轭复数。
完全无鬼波波场对应的第二鬼波算子Gmk的获取方式为,获取的原始数据还包括用于计算初始无鬼波波场对应的第二鬼波算子的计算数据,所述计算数据包括水流速度以及用于获取原始地震记录的检波器的深度,通过以下表达式计算完全无鬼波波场对应的第二鬼波算子Gmk:
Gmk=(1-e-2jz2πω/c) (4)
e为自然常数;
j为虚数单位;
z为用于获取原始地震记录的检波器的深度;
c为水流速度;
其余参数与前述各表达式中的对应的参数的含义相同。
在实际勘探采集过程中,由于采集过程中检波器深度受采环境影响不一定完全沉放于预设位置,即检波器不一定在预设深度,则考虑环境因素影响后,在环境因素影响下的初始无鬼波波场对应的第二鬼波算子的表达式为:
式(5)中:
Gm为初始无鬼波波场对应的第二鬼波算子;
为环境影响因子;环境影响因子通常由工作人员根据经验取值,优选的,其取值范围为,/>
其余参数与前述表达式中相应参数的含义相同。
对于初始无鬼波波场,初始无鬼波波场对应的第一鬼波算子和第二鬼波算子为也同样为共轭复数,即:
式(6)中,Gp为初始无鬼波波场对应的第一鬼波算子,利用初始无鬼波波场对应的第二鬼波算子Gm结合式(6)即可得到初始无鬼波波场对应的第一鬼波算子Gp。
综合所述原始地震记录、所述原始地震记录的镜像数据、所述初始无鬼波波场对应的第一鬼波算子以及所述初始无鬼波波场对应的第二鬼波算子获所述取初始无鬼波波场,初始无鬼波波场满足以下表达式:
式(7)中:
Pup1(ω,x,y)为所述初始无鬼波波场在频率域上的表达式;
P2(ω,x,y)为P(ω,x,y)的平方,P(ω,x,y)为所述原始地震记录在频率域上的表达式;
M2(ω,x,y)为M(ω,x,y)的平方,M(ω,x,y)为所述原始地震记录的镜像数据在频率域上的表达式;
为GP的平方,Gp为初始无鬼波波场对应的第一鬼波算子;
为Gm的平方,Gm为初始无鬼波波场对应的第二鬼波算子。
将式(5)代入式(7)后的得到的初始无鬼波波场的表达式为:
式(8)中各参数与前述表达式中相应参数的含义相同。
Pup1(ω,x,y)为循环叠代的初始无鬼波波场,以初始无鬼波波场Pup1(ω,x,y)作为待处理无鬼波波场,开始下述循环叠代的步骤。
根据所述待处理无鬼波波场构建判断对象。根据所述待处理无鬼波波场构建判断对象的操作步骤包括:所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取与所述待处理无鬼波波场对应的鬼波波场;
与所述待处理无鬼波波场对应的鬼波波场的表达式为:
Pdi(ω,x,y)=P(ω,x,y)-Pupi(ω,x,y) (9)
式(9)中:
Pdi(ω,x,y)为与所述待处理无鬼波波场对应的鬼波波场在频率域上的表达式;
Pupi(ω,x,y)为待处理无鬼波波场在频率域上的表达式;
下标中的i用于标示叠代次数;
其余参数含义与前述表达式对应参数的含义相同。
根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间Ti;根据所述待处理无鬼波波场、与所述待处理无鬼波波场对应的鬼波波场以及与所述待处理无鬼波波场对应的鬼波波场延时时间构建所述判断对象。所述判断对象的表达式为:
式(10)中:
O为判断对象,表征了Pup1(ω,x,y)和的差异程度;
Ti为与所述待处理无鬼波波场对应的鬼波波场延时时间;
其余各参数的含义与前述表达式相应参数的含义相同。
针对根据待处理无鬼波波场构建的判断对象进行判断,判断该判断对象是否小于或等于预设值,即判断式(10)中的判断对象O是否小于或等于预设值,预设值通常由工作人员根据经验设定,预设值的预付范围为大于0且小于或等于0.1。
若式(10)中的判断对象O小于或等于预设值,则将待处理无鬼波波场Pupi(ω,x,y)作为目标无鬼波波场并结束。
若式(10)中的判断对象O大于预设值,则根据所述待处理无鬼波波场获取所述下一级无鬼波波场。根据所述待处理无鬼波波场获取所述下一级无鬼波波场的操作步骤包括:所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间,与所述待处理无鬼波波场对应的鬼波波场延时时间在前述步骤中已获取。根据所述鬼波延时时间获取下一级无鬼波波场对应的第一鬼波算子。下一级无鬼波波场对应的第一鬼波算子满足以下表达式:
式(11)中:
Gp(i+1)为下一级无鬼波波场对应的第一鬼波算子;
其余参数的含义与前述表达式中相应参数的含义相同。
根据下一级无鬼波波场对应的第一鬼波算子和原始地震记录获取下一级无鬼波波场。
下一级无鬼波波场满足以下表达式:
式(12)中:
Pup(i+1)(ω,x,y)为下一级无鬼波波场在频率域上的表达式;
为Gp(i+1)的倒数;
其余各参数所含义与前述各表达式中相应参数所含义相同。
更新待处理无鬼波波场,将下一级无鬼波波场Pup(i+1)(ω,x,y)作为待处理无鬼波波场重复上述循环叠代步骤,直到式(10)的判断对象小于或等于预设值为止。
例如待处理无鬼波波场为初始无鬼波波场Pup1(ω,x,y),即下标i=1,表示用于第一次叠代的参数,与初始无鬼波波场Pup1(ω,x,y)对应的鬼波波场Pd1(ω,x,y)则为:
Pd1(ω,x,y)=P(ω,x,y)-Pup1(ω,x,y) (13)
式(13)中各参数含义与前述相应参数含义相同。
根据初始无鬼波波场Pup1(ω,x,y)和P(ω,x,y)可得到初始无鬼波波场Pup1(ω,x,y)对应的鬼波波场延时时间T1,此时,判断对象可表示为:
若式(11)中的判断对象O小于或等于预设值,则将Pup1(ω,x,y)作为目标无鬼波波场并结束,或若式(11)中的判断对象O大于预设值,则根据初始无鬼波波场Pup1(ω,x,y)对应的鬼波波场延时时间T1计算第二级无鬼波波场对应的第一鬼波算子GP(2),即:
式(15)中各参数含义与前述相应参数含义相同。
第二级无鬼波波场对应的第一鬼波算子GP(2)和原始地震记录P(ω,x,y)可得到第二级无鬼波场Pup(2)(ω,x,y):
式(16)中各参数含义与前述相应参数含义相同。
更新待处理无鬼波波场,将第二级无鬼波场Pup(2)(ω,x,y)作为待处理无鬼波波场,重复循环叠代步骤,得到第三级无鬼波场Pup(3)(ω,x,y)和第四级无鬼波场Pup(4)(ω,x,y)等,直到待处理无鬼波波场构建的判断对象小于或等于预设值为止。最终得到目标无鬼波波场。
本实施例通过原始地震记录、原始地震记录的镜像数据等求解初始无鬼波波场,通过循环叠代进行差异匹配和约束反演,根据差异大小,即根据构建的判断对象是否小于或等于预设值,来进行第一鬼波算子的更新,得到最优的第一鬼波算子及相应的目标无鬼波波场,从而实现斜缆数据的去鬼波处理。本实施例的方法无需镜像偏移,大大减少了计算量,通过迭代约束实现精确鬼波算子求解,从而快速准确地将鬼波去除。在最终得到的目标无鬼波波场中,鬼波已基本被消除,通过该方法能够快速准确地消除原始地震记录中的鬼波,从而提高了低频能量,很好地保证地数据的拓频效果,提高了地震资料的分辨率和精度。
为了验证基于约束反演迭代的斜缆数据鬼波压制技术的实用效果,首先采用本实施例的鬼波压制方法对海洋斜缆模拟资料进行处理。模拟的斜缆单炮记录如图1a所示,可见海底一次反射后紧随一个与之相位相反的同相轴,即鬼波。通过对模拟的单炮记录进行去鬼波处理,图1b为去鬼波后的斜缆单炮记录,图1c为去除的鬼波,与原始单炮记录对比,可以看到,鬼波同相轴得到很好的去除。
从图2频谱分析可见,图2a为去鬼波压制前的频谱,鬼波压制前的地震记录存在明显的陷频现象,随着炮检距的增大,第一个陷波点频率在逐渐减小。图2b为采用本实施例的鬼波压制方法进行鬼波压制后的频谱,鬼波引起的陷波现象得到了很好的压制,频谱的能量得到了补偿,频带得到了有效拓宽。
对某地区实际资料进行处理,图3为鬼波压制前后单炮记录对比,从图3a和图3b的对比可见,采用本实施例的鬼波压制方法进行鬼波压制后,使得紧随一次波之后的鬼波能量得到了很好压制,地震记录的分辨率得到了有效提升。
图4为鬼波压制前后的叠加剖面对比,从图4a和图4b的对比可见,采用本实施例的鬼波压制方法进行鬼波压制后,使得地震剖面的分辨率得到了有效提升,深层能量明显增强,成像更清晰,目的层成像更为清晰。
通过频谱分析,如图5所示,可以看到,采用本实施例的鬼波压制方法进行鬼波压制后,低频端能量得到了明显提升,数据的频带宽度得到了有效拓宽。
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以对其中部分或者全部技术特征进行等同替换。尤其是,只要不存在逻辑或结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。
Claims (10)
1.一种鬼波压制方法,其特征在于,包括以下步骤:
S1:获取原始数据,根据所述原始数据获取环境因素影响下的初始无鬼波波场,以所述初始无鬼波波场作为待处理无鬼波波场;
S2:根据所述待处理无鬼波波场构建判断对象,并判断所述判断对象是否小于或等于预设值;
若是,则将所述待处理无鬼波波场作为目标无鬼波波场并结束;
若否,则根据所述待处理无鬼波波场获取下一级无鬼波波场;
S3:更新所述待处理无鬼波波场,将所述下一级无鬼波波场作为所述待处理无鬼波波场并返回步骤S2;
在步骤S2中,根据所述待处理无鬼波波场构建判断对象的操作步骤包括:
所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取与所述待处理无鬼波波场对应的鬼波波场;
根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间;
根据所述待处理无鬼波波场、与所述待处理无鬼波波场对应的鬼波波场以及与所述待处理无鬼波波场对应的鬼波波场延时时间构建所述判断对象。
2.根据权利要求1所述的鬼波压制方法,其特征在于,所述判断对象的表达式为:
其中:
Pupi(ω,x,y)为待处理无鬼波波场在频率域上的表达式;
Pdi(ω,x,y)为与所述待处理无鬼波波场对应的鬼波波场在频率域上的表达式;
ω为频率;
x,y为用于获取所述原始地震记录的检波器的坐标;
j为虚数单位;
e为自然常数;
Ti为与所述待处理无鬼波波场对应的鬼波波场延时时间;
O为判断对象。
3.根据权利要求2所述的鬼波压制方法,其特征在于,所述预设值的取值范围为大于0且小于或等于0.1。
4.根据权利要求1-3中任一项所述的鬼波压制方法,其特征在于,步骤S2中,根据所述待处理无鬼波波场获取所述下一级无鬼波波场的操作步骤包括:
所述原始数据包括原始地震记录,根据所述待处理无鬼波波场和所述原始地震记录获取获取与所述待处理无鬼波波场对应的鬼波波场延时时间;
根据所述鬼波延时时间获取下一级无鬼波波场对应的第一鬼波算子;
根据下一级无鬼波波场对应的第一鬼波算子和原始地震记录获取下一级无鬼波波场。
5.根据权利要求4所述的鬼波压制方法,其特征在于,根据所述鬼波延时时间获取下一级无鬼波波场对应的第一鬼波算子的操作步骤中,下一级无鬼波波场对应的第一鬼波算子满足以下表达式:
其中:
Ti为与所述待处理无鬼波波场对应的鬼波波场延时时间;
j为虚数单位;
ω为频率;
e为自然常数;
Gp(i+1)为下一级无鬼波波场对应的第一鬼波算子。
6.根据权利要求5所述的鬼波压制方法,其特征在于,根据下一级无鬼波波场对应的第一鬼波算子和原始地震记录获取下一级无鬼波波场的操作步骤中,下一级无鬼波波场满足以下表达式:
其中:
Pup(i+1)(ω,x,y)为下一级无鬼波波场在频率域上的表达式;
P(ω,x,y)为原始地震记录在频率域上的表达式;
ω为频率;
x,y为用于获取所述原始地震记录的检波器的坐标;
为Gp(i+1)的倒数。
7.根据权利要求1-3中任一项所述的鬼波压制方法,其特征在于,步骤S1中,根据所述原始数据获取环境因素影响下的初始无鬼波波场的操作步骤包括:
所述原始数据包括原始地震记录和用于计算初始无鬼波波场对应的第二鬼波算子的计算数据;根据所述原始地震记录获取所述原始地震记录的镜像数据;根据所述计算数据获取所述初始无鬼波波场对应的第二鬼波算子,根据所述第二鬼波算子获取第一鬼波算子;
根据所述原始地震记录、所述原始地震记录的镜像数据、所述初始无鬼波波场对应的第一鬼波算子以及所述初始无鬼波波场对应的第二鬼波算子获所述初始无鬼波波场。
8.根据权利要求7所述的鬼波压制方法,其特征在于,所述初始无鬼波波场满足以下表达式:
其中:
Pup1(ω,x,y)为所述初始无鬼波波场在频率域上的表达式;
P2(ω,x,y)为P(ω,x,y)的平方,P(ω,x,y)为所述原始地震记录在频率域上的表达式;
M2(ω,x,y)为M(ω,x,y)的平方,M(ω,x,y)为所述原始地震记录的镜像数据在频率域上的表达式;
为GP的平方,Gp为初始无鬼波波场对应的第一鬼波算子;
为Gm的平方,Gm为初始无鬼波波场对应的第二鬼波算子;
9.根据权利要求8所述的鬼波压制方法,其特征在于,所述初始无鬼波波场对应的第二鬼波算子Gm满足以下表达式:
其中:
e为自然常数;
j为虚数单位;
ω为频率;
为环境影响因子;
z为用于获取原始地震记录的检波器的深度;
c为水流速度。
10.根据权利要求9所述的鬼波压制方法,其特征在于,环境影响因子的取值范围为大于0且小于或等于0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818531.4A CN112444881B (zh) | 2019-08-30 | 2019-08-30 | 一种鬼波压制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818531.4A CN112444881B (zh) | 2019-08-30 | 2019-08-30 | 一种鬼波压制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112444881A CN112444881A (zh) | 2021-03-05 |
CN112444881B true CN112444881B (zh) | 2024-04-09 |
Family
ID=74734190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910818531.4A Active CN112444881B (zh) | 2019-08-30 | 2019-08-30 | 一种鬼波压制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112444881B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514889B (zh) * | 2021-07-13 | 2022-06-21 | 中山大学 | 一种提升海洋深反射地震数据中低频信号能量的处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536045A (zh) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | 一种基于子波处理的鬼波压制方法 |
WO2016063123A2 (en) * | 2014-10-23 | 2016-04-28 | Cgg Services Sa | Systems and methods for deghosting seismic data using migration of sparse arrays |
CN109143363A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋拖缆双检采集鬼波压制方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10788597B2 (en) * | 2017-12-11 | 2020-09-29 | Saudi Arabian Oil Company | Generating a reflectivity model of subsurface structures |
-
2019
- 2019-08-30 CN CN201910818531.4A patent/CN112444881B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016063123A2 (en) * | 2014-10-23 | 2016-04-28 | Cgg Services Sa | Systems and methods for deghosting seismic data using migration of sparse arrays |
CN104536045A (zh) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | 一种基于子波处理的鬼波压制方法 |
CN109143363A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋拖缆双检采集鬼波压制方法及系统 |
Non-Patent Citations (2)
Title |
---|
基于波场延拓和反演的变深度缆地震数据鬼波压制方法;王建花;王艳冬;刘国昌;;石油地球物理勘探(05);第885-893页 * |
拉东域基于最小二乘反演的斜缆数据鬼波压制方法;张威;韩立国;李洪建;叶林;张齐;;世界地质(02);第595-601页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112444881A (zh) | 2021-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101948509B1 (ko) | 지구 물리학적 데이터의 반복 반전의 아티팩트 감소 | |
CN108845317B (zh) | 一种基于分层介质格林函数的频域逆时偏移方法 | |
CN101551465B (zh) | 一种自适应识别和消除地震勘探单频干扰的方法 | |
CN103424777A (zh) | 一种提高地震成像分辨率的方法 | |
CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
CN113805234B (zh) | 在被动源地震数据中增强面波的处理方法 | |
CN112444881B (zh) | 一种鬼波压制方法 | |
CN112593922A (zh) | 一种阵列声波测井评价固井二界面胶结质量的方法及装置 | |
CN111045076A (zh) | 多模式瑞雷波频散曲线并行联合反演方法 | |
CN113640881A (zh) | 多偏移距二维横向高分辨率瞬态面波探测方法 | |
CN110703332B (zh) | 一种鬼波压制方法 | |
CN110398778B (zh) | 一种基于等效采样的航空电磁浅层数据相关辨识方法 | |
CN110737022B (zh) | 一种可控震源激发地震资料黑三角区噪音的压制方法 | |
CN112630830A (zh) | 一种基于高斯加权的反射波全波形反演方法及系统 | |
CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
CN115712146A (zh) | 基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法 | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
CN112379427B (zh) | 面波噪声压制方法及装置 | |
CN113064205B (zh) | 菲涅尔带约束的浅水多次波衰减方法 | |
CN113138419B (zh) | 提取下行子波和衰减参数的方法、装置 | |
CN112327361B (zh) | 基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法 | |
CN115993650B (zh) | 一种基于棱柱波的地震干涉成像方法 | |
CN111722287B (zh) | 一种基于渐进数据同化方法的震相特征识别波形反演方法 | |
CN117434589A (zh) | 基于补偿生成可控震源扫描信号的方法和装置 | |
CN117890974A (zh) | 基于非线性聚束滤波的二维散射波成像方法、系统及设备 |
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 |