CN103645502B - 一种曲波域中地震波衰减补偿方法 - Google Patents

一种曲波域中地震波衰减补偿方法 Download PDF

Info

Publication number
CN103645502B
CN103645502B CN201310674626.6A CN201310674626A CN103645502B CN 103645502 B CN103645502 B CN 103645502B CN 201310674626 A CN201310674626 A CN 201310674626A CN 103645502 B CN103645502 B CN 103645502B
Authority
CN
China
Prior art keywords
curvelet
signal
omega
compensation
seismic
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
Application number
CN201310674626.6A
Other languages
English (en)
Other versions
CN103645502A (zh
Inventor
赵伟
王德利
翁斌
郝振江
王小六
张云鹏
糜芳
朱振宇
丁继才
薛东川
刘永江
杨俊�
孙文博
彭利丽
孙佳林
孟大江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
Jilin University
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University, China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical Jilin University
Priority to CN201310674626.6A priority Critical patent/CN103645502B/zh
Publication of CN103645502A publication Critical patent/CN103645502A/zh
Application granted granted Critical
Publication of CN103645502B publication Critical patent/CN103645502B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种曲波域中地震波衰减补偿方法,其包括步骤:1)收集地震信号并进行预处理,得到地震的时域信号x(t);2)将时域信号x(t)做曲波变换,得到曲波域信号ψj,lj,Tk);3)将曲波域信号ψj,lj,Tk)进行平滑处理得到曲波系数Φj,lj,Tk);4)对曲波系数Φj,lj,Tk)进行阈值迭代处理,得到去除噪声后的信号Ψj,lj,Tk);5)逐点递推求每一频段补偿角度内去除噪声后的信号Ψj,lj,Tk)的深浅层能量比δj,llj,Tk);6)求取补偿角度内的所有频段补偿因子系数矩阵Dj,llj,Tk);7)将所有频段补偿角度内的补偿因子系数矩阵Dj,llj,Tk)的倒数加权,得到曲波域加权后的信号Cj,lj,Tk);8)对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿后的时域信号xx(t)。本发明可以广泛用于地震数据处理中。

Description

一种曲波域中地震波衰减补偿方法
技术领域
本发明涉及一种在石油勘探过程中使用的地震数据处理方法,特别是关于一种曲波域中地震波衰减补偿方法。
背景技术
在石油勘探过程中,地层对地震波的吸收造成了震源子波时变,使资料的分辨率降低,特别是深层信号,由于地层的低通滤波作用,高频信号衰减严重。因此,地表记录到的不同深度反射信号及其频带,已经不能反映地下的真实情况,影响了后续地震资料的处理和解释工作,也影响了石油勘探结果的准确性。为了提高地震资料的分辨率,就必须对地震波衰减进行吸收补偿处理。
目前,地震数据处理中采用的地震波衰减补偿方法主要有常规球面发散与吸收补偿方法(牟永光.地震资料数字处理方法.石油工业出版社.1981.3)、常规反Q滤波(Futterman W I.Dispersive body wave[J].Geophysics Res.1962.67:5279-5291)、常规谱白化和单道或多道反褶积(Robinson.E.A.Principles of digital Wienerfiltering.Geophysical prospecting[J].1967.15(3):311-333)等。常规球面发散与吸收补偿方法仅补偿地震波振幅随时间变化的衰减,而不能补偿随频率变化的衰减。常规反Q(Q为品质因子,表示地层吸收的量)滤波尽管可以补偿由于地震波衰减引起的振幅减弱与速度频散,但要求品质因子Q已知,而品质因子Q值不易求准,且反Q滤波降低了数据的信噪比。常规谱白化是在假设反射系数白噪的条件下,通过对分频数据进行自动增益来补偿振幅随时间和频率的吸收衰减,因此难以保持子波和振幅能量的信息,而且会相应增强炮间的干扰能量。单道或多道反褶积是基于时窗统计内平均条件下的时频补偿,难以较好的逐点补偿振幅衰减和频率吸收。
近年来,国内学者提出几种基于时频分析的地层吸收补偿方法,如基于短时傅立叶变换地层吸收补偿方法(白桦,李鲲鹏.基于时频分析的地层吸收补偿[J].石油地球物理勘探.1999,34(6):642-648)、基于小波包分解的地层吸收补偿方法(李鲲鹏,李衍达,张学工.基于小波包分解的地层吸收衰减[J].地球物理学报.2000.43(4):542-549)和基于广义S变换的地层吸收补偿方法(刘喜武,年静波,刘洪.基于广义S变换的吸收衰减补偿方法[J].石油物探.2006.45(1):9-14)等。这三种方法利用时频分析和地震波衰减特性求得地层吸收系数,可以一定程度上消除地震波衰减对地震记录的时变影响,改善地震资料的分辨率,但是上述几种方法对信号信噪比有一定要求,在补偿地震波衰减的同时也增强了噪声和随机干扰的能量。
综上所述,现有的技术对地震波衰减补偿的处理均有不足之处,如常规球面发散与吸收衰减补偿方法仅能补偿随时间变化的能量衰减,而不能补偿随频率的吸收衰减;常规反Q滤波方法Q值不易求准,且处理的同时降低了地震信号的信噪比;常规谱白化方法难以保持振幅;单道或多道反褶积方法不能点对点补偿吸收衰减变化;基于时频分析的吸收补偿方法对地震数据信噪比有一定要求,而且在补偿地震波衰减的同时相应提高了噪声和随机干扰的能量。
发明内容
针对上述问题,本发明的目的是提供一种受地震数据信噪比限制较小且只针对有效信号进行补偿,使得不同深度反射波的波形基本一致的曲波域中地震波衰减补偿方法。
为实现上述目的,本发明采取以下技术方案:一种曲波域中地震波衰减补偿方法,其包括以下步骤:
1)收集地震信号并进行预处理,得到地震的时域信号x(t),其中t为二维空间中时域信号x(t)的空间位置变量;
2)将时域信号x(t)做曲波变换,得到曲波域信号ψj,lj,Tk),并逐点补偿地震波衰减,其中,j是曲波变换的尺度因子且j=1,2,...,N;l是曲波变换的角度因子且l=1,2,...,N,ωj是频带,Tk是反射时间;
3)将曲波域信号ψj,lj,Tk)进行平滑处理得到曲波系数Φj,lj,Tk),减小反射波相互重叠的影响;
4)对曲波系数Φj,lj,Tk)进行阈值迭代处理,得到去除噪声后的信号Ψj,lj,Tk),减少随机噪声干扰;
5)逐点递推求每一频段补偿角度内去除噪声后的信号Ψj,lj,Tk)的深浅层能量比δj,llj,Tk),使高低频信号具有相同的深浅层能量比;
6)求取补偿角度内的所有频段补偿因子系数矩阵Dj,llj,Tk),使得补偿角度系数范围内的各频段深浅层能量的比值相同;
7)将所有频段补偿角度内的补偿因子系数矩阵Dj,llj,Tk)的倒数加权,得到曲波域加权后的信号Cj,lj,Tk),实现相应频段相应角度信号的分频定向补偿;
8)对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿后的时域信号xx(t)。
所述步骤1)中,时域信号x(t)通过以下方法获得:对叠前地震资料进行处理得到二维叠加剖面数据,将所得二维地震数据按照共中心点记录方式排列,得到地震的时域信号x(t)。
所述步骤2)中,曲波域信号ψj,lj,Tk)通过以下方法得到:
在二维空间中,设时域信号x(t)的空间位置变量是t,频率域变量是ω,对应的极坐标是(r,θ),半径窗是W(r),角时窗是V(t);半径窗W(r)在范围内支撑,角时窗V(t)在t∈(-1,1)范围内支撑,且满足容许性条件:
Σ - ∞ ∞ W 2 ( 2 j r ) = 1 r ∈ ( 3 4 , 3 2 )
Σ - ∞ ∞ V 2 ( t - l ) = 1 t ∈ ( - 1 2 , 1 2 )
对于每一个j≥1,利用傅立叶变换定义频率窗函数:
其中,是j/2的整数部分,极角θ范围是0~2π;频率窗函数Uj(r,θ)的支撑区间为一楔形区域,受半径窗W(r)和角时窗V(t)的支撑区间所限制;定义其中,频率窗函数Uj(r,θ)为Uj(ω)的极坐标形式,在j=1时,为时域信号x(t)在频率域对应的曲波变换“母”函数,经旋转平移得到其它尺度(j=2,3,...,N)和所有角度(l=1,2,...,N)上的曲波基函数其中平移参数k=(k1,k2)∈Z2表示位置,k1表示纵坐标的位置,k2表示横坐标的位置,Z2表示数据在二维空间上的分布区间的位置;
定义旋转角度序列范围在0<θl<2π;平移参数k=(k1,k2)∈Z2处的曲波变换为:
其中, t k j , l = R θ l - 1 ( k 1 · 2 - j , k 2 · 2 - j ) , R θ = cos θ sin θ - sin θ cos θ 表示θ弧度的旋转,是Rθ的逆,且时域信号x(t)做曲波变换后,转换为曲波域信号ψj,lj,Tk):
所述步骤3)中,曲波系数Φj,lj,Tk)通过以下方法得到:
采用Parzen窗作为一维平滑窗函数:
d ( τ ) = 1 - 6 ( | τ | / L ) 2 + 6 ( | τ | / L ) 3 | τ | ≤ L / 2 ; 2 ( 1 - | τ | / L ) 3 L / 2 ≤ | τ | ≤ L ; 0 | τ | > L .
曲波域信号ψj,lj,Tk)与平滑窗函数d(τ)相褶积,得到曲波系数Φj,lj,Tk):
Φ j , l ( ω j , T k ) = 1 T ∫ - T 2 + T k T 2 + T k | ψ j , l ( ω j , T k ) d ( τ - T k ) | dτ
其中,T是数据时间长度,且T=5s。
所述步骤4)中,去除噪声后的信号Ψj,lj,Tk)通过以下方法得到:
利用曲波域阈值迭代衰减随机噪声影响的基本原理描述为求解如下的优化反演问题:
P ϵ : x ~ = arg mi n x | | x | | 1 = Σ i = 1 N | x i | s . t . | | Ax - y | | 2 ≤ ϵ f ~ = A x ~
其中,Pε表示求解ε的优化反演问题,ε表示极小值,A=CT是反曲波变换,x是曲波系数Φj,lj,Tk),y是所述步骤1)中含噪声的地震的时域信号x(t),ε是一任意小量,是最终输出结果;
求解上述优化问题,反演一组满足经反曲波变换后与原始数据差的二范数小于某一噪声参量的l1最小的曲波系数Φj,lj,Tk),阈值迭代去噪声后的结果
Pε问题进一步优化为如下的优化反演问题:
P λ : x ~ λ = arg min x | | y - Ax | | 2 2 + λ | | x | | 1 f ~ λ = A x ~ λ
Pλ表示求解λ的优化反演问题;
求解Pλ问题中采用迭代阈值法,其迭代序列:
x=Tλ(x+AT(y-Ax))
Tλ(x):=sgn(x)·max(0,|x|-|λ|)
其中,Tλ是软阈值函数;
阈值迭代去除噪声后,输出的结果表示为去除噪声后的信号Ψj,lj,Tk)。
所述步骤5)中,深浅层能量比δj,llj,Tk)通过以下方法得到:
去除噪声后的信号Ψj,lj,Tk)在曲波域中的地震记录为Ψ(ω):
Ψ ( ω ) = Σ 0 K r k A k ( ω , T k ) e - jω T k ω ( ω )
其中,rk(k=0,1......K)是反射系数序列,K是地震道的道数,Tk是一个反射系数rk的反射时间,是Q滤波的振幅响应,其中,Q为品质因子,表示地层吸收的量,Qeq(Tk)为反射时间为Tk处的等效Q值,A0(ω,0)是初始时刻的等效振幅衰减,ω(ω)是曲波域中的地震子波;
经过所述步骤3)平滑处理后,地震记录上相邻反射面反射波互不重叠,则每一频段特定补偿角度内的深浅层能量比δj,llj,Tk):
δ j , ll ( ω j , T k ) = | Ψ j , ll ( ω j , T k ) | | Ψ j , ll ( ω j , 0 ) | = | r k , ll A 0 , ll ( ω j , 0 ) A k , ll ( ω j , T k ) ω ^ ( ω j ) | | A 0 , ll ( ω j , 0 ) ω ^ ( ω j ) | = | r k , ll A k , ll ( ω j , T k ) |
Ψj,llj,Tk)是每个频段补偿角度内的曲波系数矩阵,Ψj,llj,0)是每个频段补偿角度内的曲波系数矩阵的浅层系数值;是不同频段补偿角度系数的振幅响应,ll是曲波变换补偿角度系数,曲波变换补偿角度系数ll的选取原则为在时间域的角度范围只要不小于地震数据同相轴的角度范围即可,此角度系数一般选为±45°。
所述步骤6)中,补偿因子系数矩阵Dj,llj,Tk)通过以下方法得到:
补偿因子系数矩阵Dj,llj,Tk)是高频信号(j=3,4,5)特定补偿角度内的深浅层能量比δj,llj,Tk)同基准频带信号对应补偿角度内的深浅层能量比δjj,llljj,Tk)之比:
D j , . ll ( ω j , T k ) = δ j , ll ( ω j , T k ) δ jj , lll ( ω jj , T k ) = e - ω j T k Q eq - 1 ( T k ) e - ω jj T k Q eq - 1 ( T k )
其中,j=3,4,5,jj是基准频段的尺度值,lll是高频段曲波变换补偿角度系数ll对应的基准频段的对应角度系数值;不同尺度角度间的对应关系为曲波变换补偿角度系数ll的值是对方程右边取整数。
所述步骤7)中,曲波域加权后的信号Cj,lj,Tk)通过以下方法得到:
曲波域加权后的信号为Cj,llj,Tk),补偿角度系数之外的信号为Ψj,l-llj,Tk),曲波变换的角度因子l为曲波变换的尺度因子j上0~2π方向,经过分频定向补偿处理后的曲波域加权后的信号Cj,lj,Tk):
C j , ll ( ω j , T k ) = Ψ j , ll ( ω j , T k ) D j . ll ( ω j , T k )
Cj,lj,Tk)=Cj,llj,Tk)+Ψj,l-llj,Tk)。
所述步骤8)中,分频定向补偿后的时域信号xx(t)通过以下方法得到:
对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿处理后的时域信号xx(t);
由曲波变换原理得到反曲波变换的实现过程:
其中,M表示指标集,φj,l,k表示由指标(j,θ,k)确定的曲波族。
本发明由于采取以上技术方案,其具有以下优点:1、由于地震波振幅的衰减是逐点时变的,本发明通过曲波变换将时域信号变换得到相应的频域信号,因此能在曲波域逐点补偿地震波的能量衰减。2、本发明的曲波变换能多层次划分频带,将高频成分进一步细分,因此具有比常用的时频变换(如傅立叶变换和小波变换等)更高的时频分辨率,能够更好地补偿高频衰减,提高地震信号的分辨率和薄层解释的准确性。3、本发明采用阈值迭代法进行去噪,设定补偿角度提取有效信号进行处理,实现了对地震波衰减的定向补偿,即只定向补偿地震波有效信号的衰减,而未增强噪声和随机干扰的能量。4、本发明由于对曲波变换数据的信噪比没有太大要求,因此不需要在处理前进行去噪处理,不但可以减少有效信号预前去噪过程中的损失,而且还可以降低地震波衰减分频定向补偿处理的工作量。5、由于本申请人已经拥有有Matlab(Matrix&Laboratory,矩阵实验室)和Linux(一种网络操作系统)两种环境下的程序包,因此,地球物理数据的处理和解释人员都可以借助程序包进行本发明的地震波衰减分频定向补偿处理,具有较好的使用和推广价值。本发明可以广泛用于石油勘探过程中对各种地震波衰减进行吸收补偿处理的过程中。
附图说明
图1是本发明浅层数据处理前(a)后(b)对比示意图
图2是本发明中深层数据处理前(a)后(b)对比示意图
图3是本发明深层数据处理前(a)后(b)对比示意图
图4是本发明数据处理前(a)后(b)第50道地震信号的频谱分析示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明曲波域中地震波衰减补偿方法包括以下步骤:
1)收集地震信号并进行预处理,得到地震的时域信号x(t)
从磁带上读取地震信号,并对叠前地震资料进行常规处理,得到二维叠加剖面数据,将所得二维地震数据按照共中心点记录方式排列,此时所得信号为地震的时域信号x(t)。如图1(a)、图2(a)、图3(a)所示,分别为地震的浅、中、深层的时域信号数据,其中,各图的横轴CDP表示采用110道数据叠加而成,纵轴t表示时间,单位为秒,图1(a)截取的是地震波衰减分频定向补偿处理前后时间轴上0.7s到1.8s处的浅层数据;图2(a)截取的是地震波衰减分频定向补偿处理前后时间轴上2.0s到3.1s处的中深层数据;图3(a)截取的是地震波衰减分频定向补偿处理前后时间轴上3.4s到4.5s处的深层数据。
2)将时域信号x(t)做曲波变换(Curvelet变换),得到曲波域信号ψj,lj,Tk),并逐点补偿地震波衰减
地震波在地层中传播时,振幅总体趋势呈指数规律衰减,且在高频和低频中具有不同的衰减趋势,即地震波振幅的衰减是逐点时变的。因此,将时域信号x(t)变换成曲波域信号ψj,lj,Tk)才可以真正实现地震波衰减的逐点补偿,其中,j是曲波变换的尺度因子,即频率域变量为ω的频带数,范围是j=1,2,...,N;l是曲波变换的角度因子,表示曲波变换的尺度因子j内的曲波变换角度数,范围是l=1,2,...,N,ωj是不同频带,Tk是反射时间。
在二维空中,设时域信号x(t)的空间位置变量是t,频率域变量是ω,对应的极坐标是(r,θ),半径窗是W(r),角时窗是V(t)。半径窗W(r)在范围内支撑,角时窗V(t)在t∈(-1,1)范围内支撑,且满足容许性条件:
Σ - ∞ ∞ W 2 ( 2 j r ) = 1 r ∈ ( 3 4 , 3 2 ) - - - ( 1 )
Σ - ∞ ∞ V 2 ( t - l ) = 1 t ∈ ( - 1 2 , 1 2 ) - - - ( 2 )
对于每一个j≥1,利用傅立叶变换定义频率窗函数:
其中,是j/2的整数部分,极角θ范围是0~2π。频率窗函数Uj(r,θ)的支撑区间为一楔形区域,受半径窗W(r)和角时窗V(t)的支撑区间所限制。定义其中,频率窗函数Uj(r,θ)为Uj(ω)的极坐标形式,在j=1时,为时域信号x(t)在频率域对应的曲波变换“母”函数,经旋转平移得到其它尺度(j=2,3,...,N)、所有角度(l=1,2,...,N)上的曲波基函数其中平移参数k=(k1,k2)∈Z2表示位置,k1表示纵坐标的位置,k2表示横坐标的位置,Z2表示数据在二维空间上的分布区间的位置。
定义旋转角度序列范围在0<θl<2π。平移参数k=(k1,k2)∈Z2处的曲波变换为:
其中, t k j , l = R θ l - 1 ( k 1 · 2 - j , k 2 · 2 - j ) , R θ = cos θ sin θ - sin θ cos θ 表示θ弧度的旋转,是Rθ的逆,且时域信号x(t)做曲波变换后,转换为曲波域信号ψj,lj,Tk):
本实施例中,曲波变换的尺度因子j为5,曲波变换的角度因子l为16,即时域信号x(t)经曲波变换后,对应的曲波域信号ψj,lj,Tk)有5个频带,第一个频带有1个方向,方向间隔为2π;第二个频带共有16个方向,方向间隔为2π/16;第三个频带共有32个方向,方向间隔为2π/32;第四个频带共有32个方向,方向间隔为2π/32;第五个频带共有64个方向,方向间隔为2π/64。本实施例中考虑第二频带。
3)将曲波域信号ψj,lj,Tk)进行平滑处理得到曲波系数Φj,lj,Tk),减小反射波相互重叠的影响
实际地震记录中,相邻反射界面的反射波是有重叠的,重叠部分的地震信号会产生能量畸变和相位失真。实际数据处理中需要对信号进行平滑,减小反射波相互重叠的影响。
对信号的平滑需要通过加平滑窗的方式来解决,因此平滑窗的选取十分重要,本发明方法选择满足四阶累积量对称性质和归一化条件,且满足在所有频率内函数值不小于0的Parzen窗(ParzenWindow,又名核密度估计,Kernel Density Estimation)来作为一维平滑窗函数对曲波域信号ψj,lj,Tk)进行平滑,其中平滑窗函数Parzen窗的表达式为:
d ( τ ) = 1 - 6 ( | τ | / L ) 2 + 6 ( | τ | / L ) 3 | τ | ≤ L / 2 ; 2 ( 1 - | τ | / L ) 3 L / 2 ≤ | τ | ≤ L ; 0 | τ | > L . - - - ( 6 )
曲波域信号ψj,lj,Tk)与平滑窗函数d(τ)相褶积,可以得到曲波系数Φj,lj,Tk):
Φ j , l ( ω j , T k ) = 1 T ∫ - T 2 + T k T 2 + T k | ψ j , l ( ω j , T k ) d ( τ - T k ) | dτ - - - ( 7 )
其中,T是数据时间长度,且T=5s。
褶积后得到曲波系数Φj,lj,Tk)比单个曲波域信号ψj,lj,Tk)和平滑窗函数d(τ)都光滑。特别地,当平滑窗函数d(τ)为具有紧致集的光滑函数时,且曲波域信号ψj,lj,Tk)为局部可积时,它们的褶积曲波系数Φj,lj,Tk)也是光滑函数。利用这一性质,对原始曲波域信号ψj,lj,Tk)进行平滑可以有效修正地震信号在不同反射界面的反射波相互重叠处的能量畸变和相位失真,使反射波相互重叠的影响降至最低。
4)对曲波系数Φj,lj,Tk)进行阈值迭代处理,得到去除噪声后的信号Ψj,lj,Tk),减少随机噪声干扰
由步骤2)中的公式(3)可以导出曲波基函数曲波基函数的支撑区间在频率域满足长度≈宽度×宽度,这种“长条形”的支撑区间是方向性的一种体现。长条形的“各向异性”基使得曲波域信号ψj,lj,Tk)在时域信号x(t)与曲波基函数方向大致相同时取得较大的曲波系数Φj,lj,Tk),时域信号x(t)与曲波基函数方向垂直或近似垂直时取得较小的曲波系数Φj,lj,Tk)。
地震数据同相轴可以近似认为是曲线反射层,因此有效信号(即地震信号同相轴)曲波系数Φj,lj,Tk)较大,而没有确定方向的噪声和随机干扰则系数较小。为了衰减噪声的影响,保留有效信号,针对曲波系数Φj,lj,Tk)进行阈值迭代处理,即选取一个阈值λ,它决定了曲波系数Φj,lj,Tk)经过阀值迭代处理后保留的比例系数。略去较小系数,保留较大系数,即可以衰减噪声的影响,避免对噪声过多补偿而降低地震数据的信噪比。
利用曲波域阈值迭代衰减随机噪声影响的基本原理可以描述为求解如下的优化反演问题:
P ϵ : x ~ = arg mi n x | | x | | 1 = Σ i = 1 N | x i | s . t . | | Ax - y | | 2 ≤ ϵ f ~ = A x ~ - - - ( 8 )
其中,Pε表示求解ε的优化反演问题,ε表示极小值,A=CT是反曲波变换,x是曲波系数Φj,lj,Tk),y是步骤1)中含噪声的地震的时域信号x(t),ε是一任意小量,是最终输出结果。
求解上述优化问题,就是要反演一组满足经反曲波变换后与原始数据差的二范数小于某一噪声参量的l1最小的曲波系数Φj,lj,Tk),最终的输出结果即为阈值λ迭代去噪声后的结果
Pε问题可以进一步优化为如下的优化反演问题:
P λ : x ~ λ = arg min x | | y - Ax | | 2 2 + λ | | x | | 1 f ~ λ = A x ~ λ - - - ( 9 )
Pλ表示求解λ的优化反演问题。
求解Pλ问题中采用迭代阈值法,其迭代序列:
x=Tλ(x+AT(y-Ax)) (10)
Tλ(x):=sgn(x)·max(0,|x|-|λ|) (11)
其中,Tλ是软阈值函数。
计算过程中,阈值λ随着循环次数的增加而降低,不断迭代更新曲波系数Φj,lj,Tk),直到取得满意的结果为止。
阈值λ迭代去除噪声后,输出的结果表示为去除噪声后的信号Ψj,lj,Tk)。
5)逐点递推求每一频段补偿角度内去除噪声后的信号Ψj,lj,Tk)的深浅层能量比δj,llj,Tk),使高低频信号具有相同的深浅层能量比
如果地震波在地层中传播时,没有经过衰减,则不同深度的反射波具有相同的振幅谱。对地震记录做时频变换,不同频段的信号应具有相似的能量-时间分布关系,即所有频段内深层反射波与浅层反射波能量之比相同。
但是,由于地层的低通滤波作用,高频信号的衰减远远大于低频信号,使得不同频段内的地震信号具有不同的能量衰减趋势,如果对高频信号乘以一个补偿因子矩阵,使得高低频信号具有相同的深浅层能量比,即补偿了地震波的衰减。若要求取有效信号的补偿因子系数矩阵,则首先要求取曲波域所有频段特定补偿角度范围内的深浅层能量比。
去除噪声后的信号Ψj,lj,Tk)在曲波域中的地震记录为Ψ(ω):
Ψ ( ω ) = Σ 0 K r k A k ( ω , T k ) e - jω T k ω ( ω ) - - - ( 12 )
其中,rk(k=0,1......K)是反射系数序列,K是地震道的道数,Tk是一个反射系数rk的反射时间,是Q滤波的振幅响应,其中,Q为品质因子,表示地层吸收的量,Qeq(Tk)为反射时间为Tk处的等效Q值,A0(ω,0)是初始时刻的等效振幅衰减,ω(ω)是曲波域中的地震子波。
经过步骤3)平滑处理后,可以认为地震记录上相邻反射面反射波互不重叠,则每一频段特定补偿角度内的深浅层能量比δj,llj,Tk):
δ j , ll ( ω j , T k ) = | Ψ j , ll ( ω j , T k ) | | Ψ j , ll ( ω j , 0 ) | = | r k , ll A 0 , ll ( ω j , 0 ) A k , ll ( ω j , T k ) ω ^ ( ω j ) | | A 0 , ll ( ω j , 0 ) ω ^ ( ω j ) | = | r k , ll A k , ll ( ω j , T k ) | - - - ( 13 )
Ψj,llj,Tk)是每个频段补偿角度内的曲波系数矩阵,Ψj,llj,0)是每个频段补偿角度内的曲波系数矩阵的浅层系数值。是不同频段补偿角度系数的振幅响应,ll是曲波变换补偿角度系数,可以根据实际情况设定,应用曲波变换补偿角度系数ll可以实现对地震波衰减的定向补偿,即在此曲波变换补偿角度系数ll范围内的曲波域地震信号会进行处理,而此曲波变换补偿角度系数ll范围外保持不变。曲波变换补偿角度系数ll的选取原则为在时间域的角度范围只要不小于地震数据同相轴的角度范围即可。对于叠后同相轴水平或近水平的地震数据,此角度系数一般选为±45°。
本实施例中,第一频段为一个矩阵,故无角度,第二、三、四频段的补偿角度系数为±45°,第五频段的系数为±60°,这样可以在有效补偿中高频信号吸收衰减的同时,尽量多的补偿高频信号的衰减。
6)求取补偿角度内的所有频段补偿因子系数矩阵Dj,llj,Tk),使得补偿角度系数范围内的各频段深浅层能量的比值相同
由于大地具有低通滤波作用,低频成分衰减较慢,可以认为低频信号基本未衰减,故选取为基准频带。为使补偿角度系数范围内的各频段深浅层能量的比值相同,因此需要求出补偿因子系数矩阵Dj,llj,Tk)。
补偿因子系数矩阵Dj,llj,Tk)是高频信号(j=3,4,5)特定补偿角度内的深浅层能量比δj,llj,Tk)同基准频带信号对应补偿角度内的深浅层能量比δjj,llljj,Tk)之比:
D j , . ll ( ω j , T k ) = δ j , ll ( ω j , T k ) δ jj , lll ( ω jj , T k ) = e - ω j T k Q eq - 1 ( T k ) e - ω jj T k Q eq - 1 ( T k ) --- ( 14 )
其中,j=3,4,5,jj是基准频段的尺度值,lll是高频段曲波变换补偿角度系数ll对应的基准频段的对应角度系数值。不同尺度角度间的对应关系为曲波变换补偿角度系数ll的值是对方程右边取整数。
7)将所有频段补偿角度内的补偿因子系数矩阵Dj,llj,Tk)的倒数加权,得到曲波域加权后的信号Cj,lj,Tk),以实现相应频段相应角度信号的分频定向补偿
由步骤6)中的公式(14)可知,补偿因子系数矩阵Dj,llj,Tk)是高频信号(j=3,4,5)特定补偿角度内的深浅层能量比δj,llj,Tk)同基准频带信号对应补偿角度内的深浅层能量比δjj,llljj,Tk)之比,那么使各个频段相应补偿角度内每个频段的曲波系数矩阵Ψj,llj,Tk)倒数加权,使得各频段深浅层能量的比值相同,分频定向补偿功能即可实现。
加权后的信号为Cj,llj,Tk),补偿角度系数之外的信号为Ψj,l-llj,Tk),曲波变换的角度因子l为曲波变换的尺度因子j上0~2π方向,故经过分频定向补偿处理后的曲波域加权后的信号Cj,lj,Tk):
C j , ll ( ω j , T k ) = Ψ j , ll ( ω j , T k ) D j . ll ( ω j , T k ) - - - ( 15 )
Cj,lj,Tk)=Cj,llj,Tk)+Ψj,l-llj,Tk)(16)
8)对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿后的时域信号xx(t)
曲波域加权后的信号为Cj,lj,Tk),对其进行反曲波变换即可以得到分频定向补偿处理后的时域信号xx(t)。
由曲波变换原理可以得到反曲波变换的实现过程:
其中,M表示指标集,φj,l,k表示由指标(j,θ,k)确定的曲波族。
如图1所示,对比图1(a)和图1(b):CDP在1-110区域内,时间轴上0.9s-1.1s之间,可以看出同相轴(地震记录上各道振动相位相同的极值的连线,即图上能量较强的黑色部分)连续性明显增强,且比处理前更加清晰,即表示处理后的浅层数据分辨率有了明显提升。对比图1(a)和图1(b):CDP在1-110区域内,时间轴上1.3s-1.6s之间,可以看出处理后的数据除了同相轴连续性大大增强外,有效信号的能量也有少许增强(CDP在1-110区域内,时间轴上1.3s-1.4s之间,处理后同相轴的能量有一定加强)。处理前杂乱无章的信号经过处理后能量明显降低(如对比图1(a)和图1(b):CDP在1-110区域内,时间轴1.1s-1.2s之间),说明处理后随机噪声得到了明显衰减。
由图1(a)和图1(b)对比结果可知:曲波域地震波衰减分频定向补偿处理后浅层地震记录同相轴的分辨率有明显提升,连续性有很大增强,有效信号的能量有少许加强,随机噪声得到了明显衰减。
如图2所示,对比图2(a)和图2(b):CDP在1-110区域内,时间轴上2.4s-2.5s之间,可以看出处理后的数据除了同相轴连续性有较大增强外,一些原本模糊不清的同相轴也变得清晰起来,这说明处理后中深层数据的分辨率有了较大改善。处理后同相轴的能量相比处理前都得到了一定增强,即有效信号的能量衰减得到了较好的补偿(在CDP从1-110区域内,时间轴2.7s-2.8s之间,处理后同相轴的能量有较大加强),噪声和随机干扰也得到了较好的压制(在CDP从1-110区域内,时间轴2.7s-2.8s之间,杂乱无章的信号得到了较大的衰减)。
由图2(a)和图2(b)对比结果可知:处理后的中深层数据同相轴分辨率有较大改善,一些原本模糊的同相轴也变的清晰,噪声和随机干扰得到了较好的压制,有效信号的能量衰减得到了较好的补偿。
如图3所示,对比图3(a)和图3(b):CDP在1-110区域内,时间轴在3.6s-3.8s之间,可以看出处理后的数据同相轴的连续性和分辨率有较大增强外,原本模糊不清的同相轴也清晰可辨,噪声和随机干扰的能量也得到了较好的衰减,这说明处理后深层数据的分辨率和信噪比有了极大提高。由于分辨率的提高,许多地质薄层也有了清楚的表示(如两图CDP在1-110区域内,时间轴上4.0s-4.1s),可以看出处理后,原先混在一起的地质薄层得以分离开来,这对于后续的地质解释和储层圈定工作有较大的意义。
由图3(a)和图3(b)对比结果可知:处理后地震信号的分辨率和信噪比有极大的提高,模糊不清的同相轴变的清晰可辨,地质薄层也有了清楚的表示。
如图4所示,对比图4(a)和图4(b)可知:频率相同时,处理后的数据振幅普遍高于补偿前,尤其是中高频(40Hz-60Hz)处,处理后地震信号的振幅值有了明显的提升,这说明使用本发明可以衰减严重的中高频信号得到有效补偿,频带展宽。
如图1~4所示的结果可以看出:本发明可以使地震波衰减得到较好补偿,不同深度反射波的波形基本一致,中深层反射波的分辨率和信噪比有明显提高,高频加强,频带展宽。
综上所述,可以将本发明归纳为以下步骤:
1)收集地震信号并进行预处理,得到地震的时域信号x(t);
2)将地震的时域信号x(t)做曲波变换,得到曲波域信号ψj,lj,Tk),并逐点补偿地震波衰减;
3)将曲波域信号ψj,lj,Tk)进行平滑处理得到曲波系数Φj,lj,Tk),减小反射波相互重叠的影响;
4)对曲波系数Φj,lj,Tk)进行阈值迭代处理,得到去除噪声后的信号Ψj,lj,Tk),减少随机噪声干扰;
5)逐点递推求每一频段补偿角度内去除噪声后的信号Ψj,lj,Tk)的深浅层能量比δj,llj,Tk),使高低频信号具有相同的深浅层能量比;
6)求取补偿角度内的所有频段补偿因子系数矩阵Dj,llj,Tk),使得补偿角度系数范围内的各频段深浅层能量的比值相同;
7)将所有频段补偿角度内的补偿因子系数矩阵Dj,llj,Tk)的倒数加权,得到曲波域加权后的信号Cj,lj,Tk),实现相应频段相应角度信号的分频定向补偿;
8)对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿后的时域信号xx(t)。
由上述步骤可知,本发明只要确定曲波变换的尺度因子j、曲波变换的角度因子l,阈值λ和曲波变换补偿角度系数ll即可实现曲波域中地震波衰减补偿。
上述各实施例仅用于说明本发明,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (9)

1.一种曲波域中地震波衰减补偿方法,其包括以下步骤:
1)收集地震信号并进行预处理,得到地震的时域信号x(t),其中t为二维空间中时域信号x(t)的空间位置变量;
2)将时域信号x(t)做曲波变换,得到曲波域信号ψj,lj,Tk),并逐点补偿地震波衰减,其中,j是曲波变换的尺度因子且j=1,2,...,N;l是曲波变换的角度因子且l=1,2,...,N,ωj是频带,Tk是反射时间;
3)将曲波域信号ψj,lj,Tk)进行平滑处理得到曲波系数Φj,lj,Tk),减小反射波相互重叠的影响;曲波系数Φj,lj,Tk)通过以下方法得到:
采用Parzen窗作为一维平滑窗函数:
d ( τ ) = 1 - 6 ( | τ | / L ) 2 + 6 ( | τ | / L ) 3 | τ | ≤ L / 2 ; 2 ( 1 - | τ | / L ) 3 L / 2 ≤ | τ | ≤ L ; 0 | τ | > L .
曲波域信号ψj,lj,Tk)与平滑窗函数d(τ)相褶积,得到曲波系数Φj,lj,Tk):
Φ j , l ( ω j , T k ) = 1 T ∫ - T 2 + T k T 2 + T k | ψ j , l ( ω j , T k ) d ( τ - T k ) | d τ
其中,T是数据时间长度,且T=5s;
4)对曲波系数Φj,lj,Tk)进行阈值迭代处理,得到去除噪声后的信号Ψj,lj,Tk),减少随机噪声干扰;
5)逐点递推求每一频段补偿角度内去除噪声后的信号Ψj,lj,Tk)的深浅层能量比δj,llj,Tk),使高低频信号具有相同的深浅层能量比;
6)求取补偿角度内的所有频段补偿因子系数矩阵Dj,llj,Tk),使得补偿角度系数范围内的各频段深浅层能量的比值相同;
7)将所有频段补偿角度内的补偿因子系数矩阵Dj,llj,Tk)的倒数加权,得到曲波域加权后的信号Cj,lj,Tk),实现相应频段相应角度信号的分频定向补偿;
8)对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿后的时域信号xx(t)。
2.如权利要求1所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤1)中,时域信号x(t)通过以下方法获得:对叠前地震资料进行处理得到二维叠加剖面数据,将所得二维地震数据按照共中心点记录方式排列,得到地震的时域信号x(t)。
3.如权利要求1所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤2)中,曲波域信号ψj,lj,Tk)通过以下方法得到:
在二维空间中,设时域信号x(t)的空间位置变量是t,频率域变量是ω,对应的极坐标是(r,θ),半径窗是W(r),角时窗是V(t);半径窗W(r)在范围内支撑,角时窗V(t)在t∈(-1,1)范围内支撑,且满足容许性条件:
Σ - ∞ ∞ W 2 ( 2 j r ) = 1 r ∈ ( 3 4 , 3 2 )
Σ - ∞ ∞ V 2 ( t - l ) = 1 t ∈ ( - 1 2 , 1 2 )
对于每一个j≥1,利用傅立叶变换定义频率窗函数:
其中,是j/2的整数部分,极角θ范围是0~2π;频率窗函数Uj(r,θ)的支撑区间为一楔形区域,受半径窗W(r)和角时窗V(t)的支撑区间所限制;定义其中,频率窗函数Uj(r,θ)为Uj(ω)的极坐标形式,在j=1时,为时域信号x(t)在频率域对应的曲波变换“母”函数,经旋转平移得到其它尺度(j=2,3,...,N)和所有角度(l=1,2,...,N)上的曲波基函数其中平移参数k=(k1,k2)∈Z2表示位置,k1表示纵坐标的位置,k2表示横坐标的位置,Z2表示数据在二维空间上的分布区间的位置;
定义旋转角度序列范围在0<θl<2π;平移参数k=(k1,k2)∈Z2处的曲波变换为:
其中,表示θ弧度的旋转,是Rθ的逆,且时域信号x(t)做曲波变换后,转换为曲波域信号ψj,lj,Tk):
4.如权利要求1或2或3所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤4)中,去除噪声后的信号Ψj,lj,Tk)通过以下方法得到:
利用曲波域阈值迭代衰减随机噪声影响的基本原理描述为求解如下的优化反演问题:
P &epsiv; : x ~ = arg min x | | x | | 1 = &Sigma; i = 1 N | x i | s . t . | | A x - y | | 2 &le; &epsiv; f ~ = A x ~
其中,Ρε表示求解ε的优化反演问题,ε表示极小值,A=CT是反曲波变换,x是曲波系数Φj,lj,Tk),y是所述步骤1)中含噪声的地震的时域信号x(t),ε是一任意小量,是最终输出结果;
求解上述优化问题,反演一组满足经反曲波变换后与原始数据差的二范数小于某一噪声参量的最小的曲波系数Φj,lj,Tk),阈值迭代去噪声后的结果
Ρε问题进一步优化为如下的优化反演问题:
P &lambda; : x &lambda; = argmin x | | y - A x | | 2 2 + &lambda; | | x | | 1 f ~ &lambda; = A x ~ &lambda;
Ρλ表示求解λ的优化反演问题;
求解Ρλ问题中采用迭代阈值法,其迭代序列:
x=Tλ(x+AT(y-Ax))
Tλ(x):=sgn(x)·max(0,|x|-|λ|)
其中,Tλ是软阈值函数;
阈值迭代去除噪声后,输出的结果表示为去除噪声后的信号Ψj,lj,Tk)。
5.如权利要求1或2或3所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤5)中,深浅层能量比δj,llj,Tk)通过以下方法得到:
去除噪声后的信号Ψj,lj,Tk)在曲波域中的地震记录为Ψ(ω):
&Psi; ( &omega; ) = &Sigma; 0 K r k A k ( &omega; , T k ) e - j&omega;T k &omega; ( &omega; )
其中,rk(k=0,1......K)是反射系数序列,K是地震道的道数,Tk是一个反射系数rk的反射时间,是Q滤波的振幅响应,其中,Q为品质因子,表示地层吸收的量,Qeq(Tk)为反射时间为Tk处的等效Q值,A0(ω,0)是初始时刻的等效振幅衰减,ω(ω)是曲波域中的地震子波;
经过所述步骤3)平滑处理后,地震记录上相邻反射面反射波互不重叠,则每一频段特定补偿角度内的深浅层能量比δj,llj,Tk):
&delta; j , l l ( &omega; j , T k ) = | &Psi; j , l l ( &omega; j , T k ) | | &Psi; j , l l ( &omega; j , 0 ) | = | r k , l l A 0 , l l ( &omega; j , 0 ) A k , l l ( &omega; j , T k ) &omega; ^ ( &omega; j ) | | A 0 , l l ( &omega; j , 0 ) &omega; ^ ( &omega; j ) | = | r k , l l A k , l l ( &omega; j , T k ) |
Ψj,llj,Tk)是每个频段补偿角度内的曲波系数矩阵,Ψj,llj,0)是每个频段补偿角度内的曲波系数矩阵的浅层系数值;是不同频段补偿角度系数的振幅响应,ll是曲波变换补偿角度系数,曲波变换补偿角度系数ll的选取原则为在时间域的角度范围只要不小于地震数据同相轴的角度范围即可,此角度系数一般选为±45°。
6.如权利要求4所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤6)中,补偿因子系数矩阵Dj,llj,Tk)通过以下方法得到:
补偿因子系数矩阵Dj,llj,Tk)是高频信号(j=3,4,5)特定补偿角度内的深浅层能量比δj,llj,Tk)同基准频带信号对应补偿角度内的深浅层能量比δjj,llljj,Tk)之比:
D j , l l ( &omega; j , T k ) = &delta; j , l l ( &omega; j , T k ) &delta; j j , l l l ( &omega; j j , T k ) = e - &omega; j T k Q e q - 1 ( T k ) e - &omega; j j T k Q e q - 1 ( T k )
其中,j=3,4,5,jj是基准频段的尺度值,lll是高频段曲波变换补偿角度系数ll对应的基准频段的对应角度系数值;不同尺度角度间的对应关系为曲波变换补偿角度系数ll的值是对方程右边取整数。
7.如权利要求1或2或3或6所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤7)中,曲波域加权后的信号Cj,lj,Tk)通过以下方法得到:
曲波域加权后的信号为Cj,llj,Tk),补偿角度系数之外的信号为Ψj,l-llj,Tk),曲波变换的角度因子l为曲波变换的尺度因子j上0~2π方向,经过分频定向补偿处理后的曲波域加权后的信号Cj,lj,Tk):
C j , l l ( &omega; j , T k ) = &Psi; j , l l ( &omega; j , T k ) D j . l l ( &omega; j , T k )
Cj,lj,Tk)=Cj,llj,Tk)+Ψj,l-llj,Tk)。
8.如权利要求4所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤8)中,分频定向补偿后的时域信号xx(t)通过以下方法得到:
对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿处理后的时域信号xx(t);
由曲波变换原理得到反曲波变换的实现过程:
其中,M表示指标集,φj,l,k表示由指标(j,θ,k)确定的曲波族。
9.如权利要求5所述的一种曲波域中地震波衰减补偿方法,其特征在于:所述步骤8)中,分频定向补偿后的时域信号xx(t)通过以下方法得到:
对曲波域加权后的信号Cj,lj,Tk)进行反曲波变换,得到分频定向补偿处理后的时域信号xx(t);
由曲波变换原理得到反曲波变换的实现过程:
其中,M表示指标集,φj,l,k表示由指标(j,θ,k)确定的曲波族。
CN201310674626.6A 2013-12-11 2013-12-11 一种曲波域中地震波衰减补偿方法 Active CN103645502B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310674626.6A CN103645502B (zh) 2013-12-11 2013-12-11 一种曲波域中地震波衰减补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310674626.6A CN103645502B (zh) 2013-12-11 2013-12-11 一种曲波域中地震波衰减补偿方法

Publications (2)

Publication Number Publication Date
CN103645502A CN103645502A (zh) 2014-03-19
CN103645502B true CN103645502B (zh) 2016-08-17

Family

ID=50250750

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310674626.6A Active CN103645502B (zh) 2013-12-11 2013-12-11 一种曲波域中地震波衰减补偿方法

Country Status (1)

Country Link
CN (1) CN103645502B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105319593A (zh) * 2014-07-04 2016-02-10 中国石油化工股份有限公司 基于曲波变换和奇异值分解的联合去噪方法
CN104820242B (zh) * 2015-03-27 2017-06-06 中国石油化工股份有限公司 一种面向叠前反演的道集振幅分频补偿方法
CN104849757B (zh) * 2015-04-30 2018-02-27 长江大学 消除地震信号中随机噪声系统及方法
CN104932009B (zh) * 2015-05-29 2017-05-03 西北工业大学 补偿Morlet小波变换的复时‑频谱提高地震剖面分辨率的方法
CN105093329B (zh) * 2015-08-25 2017-07-07 中国石油天然气集团公司 地震数据的低频信息的自动补偿方法及装置
CN106910164A (zh) * 2015-12-23 2017-06-30 通用电气公司 一种对ct投影数据进行滤波的方法及装置
CN107238864A (zh) * 2016-03-28 2017-10-10 中国石油化工股份有限公司 基于优势频带约束的地震波能量时频补偿方法
CN107390267B (zh) * 2017-07-27 2019-03-01 西安交通大学 一种同步挤压变换域的地震资料衰减补偿方法
CN107589450B (zh) * 2017-09-01 2019-01-04 中国科学院地质与地球物理研究所 基于曲波变换和聚类的地震数据去噪方法和装置
CN109917459A (zh) * 2017-12-12 2019-06-21 中国石油天然气股份有限公司 一种压制地震噪声的方法、装置及系统
CN108519619B (zh) * 2018-04-11 2019-11-05 中国石油大学(华东) 基于特殊化曲波变换的时频分析技术
CN110967743B (zh) * 2018-09-28 2022-04-22 中国石油化工股份有限公司 一种分频迭代地震反演方法及系统
CN112698396B (zh) * 2019-10-23 2024-04-30 中国石油天然气集团有限公司 一种压制极强环境噪音的数据处理方法及系统
CN113219536B (zh) * 2021-06-25 2022-03-01 成都理工大学 依赖频率的纵横波衰减参数的叠前地震反演方法
CN114692676B (zh) * 2022-03-04 2023-06-20 华南理工大学 一种高频底座力天平复杂畸变信号的修正方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183787A (zh) * 2011-03-07 2011-09-14 中国海洋石油总公司 一种基于地震记录变子波模型提高地震资料分辨率的方法
CN103412329A (zh) * 2013-08-06 2013-11-27 中国海洋石油总公司 一种提高地震数据分辨率的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2168057B1 (en) * 2007-07-16 2019-09-18 ExxonMobil Upstream Research Company Geologic features from curvelet based seismic attributes
US8588028B2 (en) * 2009-09-22 2013-11-19 Westerngeco L.L.C. Slant-stack for time gather to angle gather conversion

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183787A (zh) * 2011-03-07 2011-09-14 中国海洋石油总公司 一种基于地震记录变子波模型提高地震资料分辨率的方法
CN103412329A (zh) * 2013-08-06 2013-11-27 中国海洋石油总公司 一种提高地震数据分辨率的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
An iterative curvelet thresholding algorithm for seismic random noise attenuation;Wang De-Li 等;《Applied Geophysics》;20101231;第7卷(第4期);第315-324页 *
Curvelet阈值迭代法在地震数据去噪和插值中的应用研究;仝中飞;《中国优秀硕士学位论文全文数据库·基础科学辑》;20090815(第8期);第A011-305页 *
基于Curvelet变换的地层吸收补偿方法研究;孙佳林;《中国优秀硕士学位论文全文数据库·基础科学辑》;20130815(第8期);第A011-154页 *
基于Curvelet变换的地震数据去噪方法研究;刘冰;《中国优秀硕士学位论文全文数据库·信息科学辑》;20081115(第11期);第I138-825页 *
稀疏变换反褶积和高精度子波提取方法研究;孟大江;《中国优秀硕士学位论文全文数据库·基础科学辑》;20121015(第10期);第A011-282页 *

Also Published As

Publication number Publication date
CN103645502A (zh) 2014-03-19

Similar Documents

Publication Publication Date Title
CN103645502B (zh) 一种曲波域中地震波衰减补偿方法
CA2707526C (en) Spectral shaping inversion and migration of seismic data
US9625593B2 (en) Seismic data processing
CN111208561A (zh) 基于时变子波与曲波变换约束的地震声波阻抗反演方法
CN110187388B (zh) 一种基于变分模态分解的稳定地震品质因子q估计方法
WO2007143355A2 (en) Diplet-based seismic processing
CN106680876B (zh) 一种地震数据联合去噪方法
CN102681014A (zh) 基于多项式拟合的规则线性干扰压制方法
Bakulin et al. Application of supergrouping to enhance 3D prestack seismic data from a desert environment
CN102305941A (zh) 由叠前时间偏移直接扫描确定地层叠加品质因子方法
CN106033125B (zh) 压制叠前大角度道集干涉的提频方法
CN111522062B (zh) 基于火山岩屏蔽量化分析的下伏地层振幅补偿法
CN112285778B (zh) 一种粘声TTI介质中纯qP波的逆时偏移成像方法
Banjade et al. Earthquake accelerogram denoising by wavelet-based variational mode decomposition
CN110780341B (zh) 一种各向异性地震成像方法
Nemeth et al. Separation of signal and coherent noise by migration filtering
Bakulin et al. Importance of phase guides from beamformed data for processing multi-channel data in highly scattering media
CN109143345B (zh) 基于模拟退火的品质因子q非线性反演方法及系统
CN110873900B (zh) 频率域叠前地震道q补偿方法及系统
US7889598B2 (en) Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity
Yan et al. Seismic spatially variant noise suppression method in the Tarim Basin based on FFDNet and transfer learning
Cheng et al. Q estimation based on the logarithmic spectral area double difference
CN112213775B (zh) 一种高覆盖次数叠前地震数据的保真提频方法
CN102998702B (zh) 保幅平面波叠前深度偏移方法
Bakulin et al. Seismic speckle as multiplicative noise explaining land reflections distorted by near-surface scattering

Legal Events

Date Code Title Description
PB01 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
CP01 Change in the name or title of a patent holder

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Co-patentee after: Jilin University

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC Research Institute

Patentee before: China National Offshore Oil Corporation

Co-patentee before: Jilin University

CP01 Change in the name or title of a patent holder
TR01 Transfer of patent right

Effective date of registration: 20200115

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC research institute limited liability company

Patentee before: China Offshore Oil Group Co., Ltd.

Co-patentee before: Jilin University

TR01 Transfer of patent right