CN101644782A - 基于极化滤波的多波多分量地震资料的去噪方法 - Google Patents

基于极化滤波的多波多分量地震资料的去噪方法 Download PDF

Info

Publication number
CN101644782A
CN101644782A CN200910167570A CN200910167570A CN101644782A CN 101644782 A CN101644782 A CN 101644782A CN 200910167570 A CN200910167570 A CN 200910167570A CN 200910167570 A CN200910167570 A CN 200910167570A CN 101644782 A CN101644782 A CN 101644782A
Authority
CN
China
Prior art keywords
parameter
ground roll
data
wave
horizontal
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
Application number
CN200910167570A
Other languages
English (en)
Other versions
CN101644782B (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.)
SHENGLI PETROLEUM ADMINISTRATION OF SINOPEC GROUP
Original Assignee
SHENGLI PETROLEUM ADMINISTRATION OF SINOPEC GROUP
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 SHENGLI PETROLEUM ADMINISTRATION OF SINOPEC GROUP filed Critical SHENGLI PETROLEUM ADMINISTRATION OF SINOPEC GROUP
Priority to CN2009101675709A priority Critical patent/CN101644782B/zh
Publication of CN101644782A publication Critical patent/CN101644782A/zh
Application granted granted Critical
Publication of CN101644782B publication Critical patent/CN101644782B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于极化滤波的多波多分量地震资料的去噪方法,该方法包括:A.根据所述多波多分量地震资料的纵波初始数据和转换波初始数据,获得面波初始水平参量和初始垂直参量;B.根据所述面波初始水平参量和初始垂直参量,利用极化滤波迭代方法计算面波目标水平参量和面波目标垂直参量;C.从所述纵波初始数据中减去所述面波目标垂直参量,及从所述转换波初始数据中减去所述面波目标水平参量,以得到去除噪音的数据。本发明所具有的有益效果:提高了转换波资料及纵波资料的信噪比。

Description

基于极化滤波的多波多分量地震资料的去噪方法
技术领域
本发明属于地震资料处理领域,尤其涉及一种基于极化滤波的多波多分量地震资料的去噪方法。
背景技术
地震勘探一直以来都是人们认识和了解地下介质的一种重要手段,多波多分量地震勘探是其中的一种方法。所谓多波多分量地震勘探,是区别与现有单一的纵波和横波勘探来说的。单一的纵波和横波勘探是用纵波和横波震源激发,用单分量(水平或者垂直)检波器来接收,根据接收的地震资料来探测地下的地质构造或岩性。而多波多分量地震勘探是分别用纵波和横波(沿测线方向偏振和垂直方向偏振)的震源激发,用三分量(一个垂直分量和两个水平分量)检波器来接收,从而得到多波多分量地震资料。多波多分量地震资料与上述单一波测量的地震资料相比具有信息量大及准确的优点。
在多波多分量地震资料中包含纵波资料和转换波资料,单一的纵波资料和转换波资料是我们想要的,但实际得到的纵波资料和转换波资料往往还含有对地震资料产生主要噪音影响的面波(如瑞雷面波),因此如何将面波去掉是人们比较关心的问题。常规去除面波的方法是使用带通滤波、F-K滤波、τ-p变换等方法。但带通滤波会将有效的低频成分滤掉,而F-K滤波及τ-p变换会产生一定的混波效应从而使有效波的保真度会受到影响。例如,在胜利油田垦71地区,由于面波的分布范围广泛,而且随偏移距的增加而增加,呈不规则分布,从而使地震资料的信噪比和分辨率较低。因此,常规的去噪方法无法有效地将地震资料中的面波去除。
发明内容
本发明实施例的目的在于,提供一种基于极化滤波的多波多分量地震资料的去噪方法,用以解决面波的分布范围广泛,而且随偏移距的增加而增加,呈不规则分布,从而使多波多分量地震资料的信噪比和分辨率较低的问题。
为解决上述技术问题,本发明实施例提供一种基于极化滤波的多波多分量地震资料的去噪方法,该方法包括以下步骤:
A、根据多波多分量地震资料的纵波初始数据和转换波初始数据,获得面波初始水平参量和初始垂直参量;
B、根据所述面波初始水平参量和初始垂直参量,利用极化滤波迭代方法计算面波目标水平参量和面波目标垂直参量;
C、从所述纵波初始数据中减去所述面波目标垂直参量,及从所述转换波初始数据中减去所述面波目标水平参量,以得到去除噪音的数据。
其中,所述极化滤波迭代方法包括以下步骤:
B1、利用极化椭圆参数求解方法求解得到多个椭圆;
B2、对所述求解的多个椭圆进行全局优选得到所需的椭圆(P),并将椭圆(P)的长轴(a)赋值给面波初始垂直参量及短轴的负值(-b)赋值给面波初始水平参量;
B3、根据所述面波初始水平参量和面波初始垂直参量分别计算转换波初始数据与面波初始水平参量差值(R0(N)-Rp′(N))的振幅包络(Erp(N))以及纵波初始数据与面波初始垂直参量差值(Z0(N)-Zp′(N))的振幅包络(Ezp(N));
B4、根据所述振幅包络(Erp(N))和(Ezp(N))利用互相关公式来计算此时面波水平参量的振幅包络(E′r)和面波垂直参量振幅包络(E′z);
B5、根据所述E′z和E′r来计算当所述(R0(N)-Rp′(N))与E′r的模量值(||(R0(N)-Rp(N)),E′r||)以及(Z0(N)-Zp′(N))与E′z的模量值(||(Z0(N)-Zp(N)),E′z||)为最小值时的面波水平参量和面波垂直参量,当所述两个模量值都达到最小值时计算得到的所述面波水平参量和面波垂直参量即为面波目标水平参量和目标垂直参量,否则将求得的所述面波水平参量和面波垂直参量分别赋值给面波初始水平参量和面波初始垂直参量后依次重复步骤B3、B4和B5。
所述极化椭圆参数求解方法是利用最小平方二乘法完成。
步骤A之前还可包括:获取多波多分量地震资料的纵波资料和转换波资料;在预设的频带范围内对获取的纵波资料和转换波资料进行低通滤波,以获得纵波初始数据和转换波初始数据。
步骤A中所述获得面波初始水平参量和初始垂直参量采用视速度差异、F-K滤波或者τ-p变换而实现。
本发明实施例具有以下有益效果:通过面波极化分析及椭圆拟合的方法,获得面波目标水平参量和目标垂直参量,然后从初始纵波数据和初始转换波数据中将面波目标水平参量和目标垂直参量减去,从而提高了转换波资料及纵波资料的信噪比。
附图说明
图1为本发明实施例提供的基于极化滤波的多波多分量地震资料去噪方法的一个较佳实施例的处理流程示意图;
图2为转换波资料单炮去噪前的效果示意图;
图3为对转换波资料进行滤波的效果对比示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本发明这种去噪方法,利用了面波在空间的极化轨迹是椭圆,而纵波及转换波在空间的极化轨迹为直线这一特性,通过三分量检波器得到的原始的纵波和转换波记录,对原始的纵波和转换波记录应用常规的压制面波方法(如F-K滤波及τ-p变换等)得到面波的初始水平参量及初始垂直参量,并把该初始水平参量及垂直参量作为初始值输入以加快极化滤波迭代法收敛速度,然后根据该极化滤波迭代法计算出面波目标水平参量和面波目标垂直参量,最后从所述纵波初始数据和转换波初始数据中减去所述面波目标水平参量和面波目标垂直参量,就可以有效的去除纵波及转换波中的面波,进而提高纵波及转换波资料的信噪比。
实施例一
图1为本发明实施例提供的基于极化滤波的多波多分量地震资料去噪方法的一个较佳实施例的处理流程示意图。该流程示意图包括以下步骤:
步骤101、获取多波多分量地震资料的纵波资料和转换波资料,并在预设的频带范围内对获取的纵波资料和转换波资料进行低通滤波,以获得纵波初始数据和转换波初始数据。
其中,该步骤的目的是将多波多分量地震资料的纵波资料和转换波资料中的高频成分滤除,并尽可能保留纵波资料和转换波资料中的面波信息。频带范围的设置由本领域技术人员根据实际情况来设置。
步骤102、根据所述纵波初始数据和转换波初始数据,获得面波初始水平参量和初始垂直参量。
本步骤中,获得面波初始水平参量和初始垂直参量的方式可以与现有技术相同,例如,根据所述纵波初始数据中包含的面波及纵波之间的视速度差异获得面波初始垂直参量,及所述转换波初始数据中包含的面波与转换波之间的视速度差异获得面波初始水平参量,或者采用F-K滤波及τ-p变换等。当然可以采用其他方式获得,只要能从纵波初始数据和转换波初始数据获得包含尽量多信息的面波初始水平参量和初始垂直参量即可。
步骤103、根据所述面波初始水平参量和初始垂直参量,利用极化滤波迭代方法计算出面波目标水平参量和面波目标垂直参量。
本步骤中,极化滤波迭代方法可以根据下面的步骤来实现:
B1、利用极化椭圆参数求解方法求解得到多个椭圆;
在本步骤中,可利用最小平方二乘法求解得到多个椭圆。当然也可以采用其它数学计算方法,只要能实现求解椭圆参数,并最终得到椭圆方程即可。
B2、对所述求解的多个椭圆进行全局优选得到所需的椭圆P,并将椭圆P的长轴a赋值给面波初始垂直参量及短轴的负值-b赋值给面波初始水平参量,即b=-Rp(N)、a=Zp(N);
在该步骤中,进行全局优选时可根据选择长短轴比值最大值时对应的椭圆。
B3、分别根据下面的公式(1)和(2)计算转换波初始数据与面波初始水平参量差值(R0(N)-Rp(N))的振幅包络Erp(N)以及纵波初始数据与面波初始垂直参量差值(Z0(N)-Zp(N))的振幅包络Ezp(N);
Erp(N)=E(R0(N)-Rp(N))        (1)
Ezp(N)=E(Z0(N)-Zp(N))        (2)
B4、根据步骤B3所求得的振幅包络Erp(N)和Ezp(N)利用下面的公式(3)和(4),当ERR和Epp达到最大时计算此时面波水平参量的振幅包络E′r和垂直参量振幅包络E′z
E RR ( i ) max = Σ k = - N / 2 N / 2 E rp ( k ) * E r ′ ( i + k ) , i = 1,2,3 , . . . , N - - - ( 3 )
E PP ( i ) max = Σ k = - N / 2 N / 2 E zp ( k ) * E p ′ ( i + k ) , i = 1,2,3 , . . . , N - - - ( 4 )
B5、根据步骤B4求得的E′r和E′z,利用下面的公式(5)和(6)分别来计算(R0(N)-Rp(N))与E′r的模量值||(R0(N)-Rp(N)),E′r||以及(Z0(N)-Zp(N))与E′z的模量值||(Z0(N)-Zp(N)),E′z||为最小值时的面波水平参量和面波垂直参量,当所述两个模量值都达到最小值时计算得到的所述面波水平参量和面波垂直参量即为面波目标水平参量和目标垂直参量分别为Rp′(N)、Zp′(N),否则将求得的所述面波水平参量和面波垂直参量分别赋值给面波初始水平参量和面波初始垂直参量后,依次重复步骤B3、B4和B5。
||(R0(N)-Rp(N)),E′r||min     (5)
||(Z0(N)-Zp(N)),E′z||min     (6)
步骤104、从所述纵波初始数据中减去所述面波目标垂直参量,及从所述转换波初始数据中减去所述面波目标水平参量,以得到去除噪音的数据。
下面将步骤103中使用极化滤波迭代法计算面波目标水平参量和垂直参量的过程进行如下详细描述:
一、极化椭圆参数的求解方法原理
椭圆曲线方程可为:
Ax2+Bxy+Cy2+Dx+Ey+F=0           (7)
其中,椭圆几何参数设置为:中心位置(X0,Y0),长短轴分别为(a,b),长轴a的转角为θ,A、B、C、D、E、F为椭圆曲线方程的系数,则椭圆几何参数与椭圆曲线方程的系数存在如下关系:
X 0 = BE - 2 CD 4 AC - B 2 Y 0 = BD - 2 AE 4 AC - B 2
Figure A20091016757000093
“±”取“+”时对应长轴a,取“-”时对应短轴b;
θ = 1 2 tan - 1 B A - C
从转换波初始数据R0(N)和纵波初始数据Z0(N)中任意选择5个点用于求取椭圆曲线方程的系数,即可以表示成如下的形式:
{ x 1 , x 2 , x 3 , x 4 , x 5 } ∈ C N 5 ( R N ) { y 1 , y 2 , y 3 , y 4 , y 5 } ∈ C N 5 ( Z N )
令MX=Y
其中,
Y = 1 1 1 1 1 X = A B C D E M = x 1 2 x 1 y 1 y 1 2 x 1 y 1 x 2 2 x 2 y 2 y 2 2 x 2 y 2 x 3 2 x 3 y 3 y 3 2 x 3 y 3 x 4 2 x 4 y 4 y 4 2 x 4 y 4 x 5 2 x 5 y 5 y 5 2 x 5 y 5
将上述两边左乘上M′,使得
M′MX=M′Y,其中M′为M的转置矩阵;
然后,可利用最小平方二乘法计算出椭圆曲线方程的系数A、B、C、D、E,当B2-4AC<0时,即可得到椭圆曲线方程的解。
二、极化滤波迭代方法原理
首先,在定周期时间范围N内对CN 5(R0(N),Z0(N))任意数据组合使用上述的极化椭圆参数的求解方法来得到椭圆参数,然后选择长短轴比值的最大值对应的椭圆P(a,b,θ,X0,Y0)。
然后,设转换波初始数据和纵波初始数据分别为R0(N)、Z0(N),面波初始水平参量和面波初始垂直参量分别为Rp(N)和Zp(N),并分别计算转换波初始数据与面波初始水平参量差值(R0(N)-Rp(N))的振幅包络Erp(N)以及纵波初始数据与面波初始垂直参量差值(Z0(N)-Zp(N))的振幅包络Ezp(N)。
然后,根据上述所求得的振幅包络Erp(N)和Ezp(N),利用下面的公式,当ERR和Epp达到最大时计算此时面波水平参量的振幅包络E′r和垂直参量振幅包络E′z
E RR ( i ) max = Σ k = - N / 2 N / 2 E rp ( k ) * E r ′ ( i + k )
E PP ( i ) max = Σ k = - N / 2 N / 2 E zp ( k ) * E p ′ ( i + k )
最后,分别计算当(R0(N)-Rp(N))与E′r的模量值||(R0(N)-Rp(N)),E′r||以及(Z0(N)-Zp(N))与E′z的模量值||(Z0(N)-Zp(N)),E′z||为最小值时的面波水平参量和面波垂直参量,当所述两个模量值都达到最小值时计算得到的所述面波水平参量和面波垂直参量即为面波目标水平参量和目标垂直参量分别为Rp′(N)、Zp′(N),否则将求得的所述面波水平参量和面波垂直参量分别赋值给面波初始水平参量和面波初始垂直参量后重新进行上面的计算过程。
由上述的实施例可见,本发明这种去除噪音的方法,利用了极化滤波迭代方法计算出了相对精确的面波目标水平参量和目标垂直参量,然后从转换波及纵波初始数据中减去面波目标水平参量和目标垂直参量即可有效的去除噪音,从而达到提高转换波资料和纵波资料的信噪比。
实施例二
根据实施例一提供的方法,本实施例给出了一个对在胜利油田垦71地区三维转换波实际资料去除面波的具体实例。如图2所示为转换波资料单炮去噪前的效果示意图,从图2中可以看到,由于采用单点数字检波器接收,转换波资料的噪音较纵波资料严重的多,而且低频、能量较强的面波发育范围广泛。图3为对转换波资料进行滤波的效果对比示意图,其中a图为使用常规滤波方法达到的滤波效果,b图为采用本发明实施例一提供的极化滤波迭代方法达到的滤波效果,从图3中可以看到采用本发明实施例一提供的方法,能有效去除转换波和纵波中的面波,从而提高了转换波资料及纵波资料的信噪比。
显然,上述实施例仅是本发明优选实施例而已,它并不限制本发明的保护范围,在本发明的保护范围内,所述领域技术人员还可以对本发明的方法做出各种改进及润饰,当然这些改进和润饰也视为本发明的保护范围。

Claims (5)

1、一种基于极化滤波的多波多分量地震资料的去噪方法,其特征在于,该方法包括以下步骤:
A、根据所述多波多分量地震资料的纵波初始数据和转换波初始数据,获得面波初始水平参量和初始垂直参量;
B、根据所述面波初始水平参量和初始垂直参量,利用极化滤波迭代方法计算面波目标水平参量和面波目标垂直参量;
C、从所述纵波初始数据中减去所述面波目标垂直参量,及从所述转换波初始数据中减去所述面波目标水平参量,以得到去除噪音的数据。
2、根据权利要求1所述的方法,其特征在于,所述极化滤波迭代方法包括:
B1、利用极化椭圆参数求解方法求解得到多个椭圆;
B2、对所述求解的多个椭圆进行全局优选得到所需的椭圆(P),并将椭圆(P)的长轴(a)赋值给面波初始垂直参量及短轴的负值(-b)赋值给面波初始水平参量;
B3、根据所述面波初始水平参量和面波初始垂直参量分别计算转换波初始数据与面波初始水平参量差值(R0(N)-R′p(N))的振幅包络(Erp(N))以及纵波初始数据与面波初始垂直参量差值(Z0(N)-Z′p(N))的振幅包络(Ezp(N));
B4、根据所述振幅包络Erp(N)和Ezp(N)利用互相关公式来分别计算此时面波水平参量的振幅包络(E′r)和面波垂直参量振幅包络(E′z);
B5、根据所述E′z和E′r来计算当所述(R0(N)-R′p(N))与E′r的模量值(||(R0(N)-Rp(N)),E′r||)以及(Z0(N)-Z′p(N))与E′z的模量值(||(Z0(N)-Zp(N)),E′z||)为最小值时的面波水平参量和面波垂直参量,当所述两个模量值都达到最小时计算得到的所述面波水平参量和面波垂直参量即为面波目标水平参量和目标垂直参量,否则将求得的所述面波水平参量和面波垂直参量分别赋值给面波初始水平参量和面波初始垂直参量后依次重复步骤B3、B4和B5。
3、根据权利要求2所述的方法,其特征在于,所述极化椭圆参数求解方法利用最小平方二乘法实现。
4、根据权利要求1所述的方法,其特征在于,步骤A之前还包括:
获取多波多分量地震资料的纵波资料和转换波资料;
在预设的频带范围内对获取的纵波资料和转换波资料进行低通滤波,以获得纵波初始数据和转换波初始数据。
5、根据权利要求1所述的方法,其特征在于,步骤A中所述获得面波初始水平参量和初始垂直参量采用视速度差异、F-K滤波或者τ-p变换而实现。
CN2009101675709A 2009-08-25 2009-08-25 基于极化滤波的多波多分量地震资料的去噪方法 Expired - Fee Related CN101644782B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101675709A CN101644782B (zh) 2009-08-25 2009-08-25 基于极化滤波的多波多分量地震资料的去噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101675709A CN101644782B (zh) 2009-08-25 2009-08-25 基于极化滤波的多波多分量地震资料的去噪方法

Publications (2)

Publication Number Publication Date
CN101644782A true CN101644782A (zh) 2010-02-10
CN101644782B CN101644782B (zh) 2011-12-07

Family

ID=41656745

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101675709A Expired - Fee Related CN101644782B (zh) 2009-08-25 2009-08-25 基于极化滤波的多波多分量地震资料的去噪方法

Country Status (1)

Country Link
CN (1) CN101644782B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915939A (zh) * 2010-08-17 2010-12-15 中国科学院地质与地球物理研究所 一种面波压制方法
CN102087367A (zh) * 2010-11-19 2011-06-08 中冶集团武汉勘察研究院有限公司 纵波-面波速度比值法划分基岩风化带的方法
CN102338886A (zh) * 2010-07-14 2012-02-01 中国石油天然气集团公司 一种有效衰减三分量地震记录中面波的极化滤波方法
CN103091714A (zh) * 2011-10-28 2013-05-08 中国石油化工股份有限公司 一种自适应面波衰减方法
CN103323877A (zh) * 2013-05-30 2013-09-25 吉林大学 一种基于地震勘探环境噪声指向性的去噪方法
CN104459769A (zh) * 2013-09-22 2015-03-25 中国石油化工股份有限公司 一种地震图像增强方法
CN104849752A (zh) * 2015-05-18 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 获取地震记录中的转换横波的高精度反褶积算子的方法
CN103217709B (zh) * 2012-01-18 2015-09-23 中国石油天然气集团公司 一种提高地震数据信噪比和分辨率的面波衰减方法
CN105652322A (zh) * 2016-01-07 2016-06-08 中国科学院地球化学研究所 多分量地震数据的t-f-k域极化滤波方法
CN106249288A (zh) * 2016-08-02 2016-12-21 中国地质大学(北京) 基于Shearlet域的极化滤波面波压制方法
CN106353791A (zh) * 2015-07-17 2017-01-25 中国石油化工股份有限公司 基于波形特征的多波多分量资料联合属性储层预测方法
CN107966734A (zh) * 2017-09-22 2018-04-27 中国地质大学(北京) 多分量地震数据的矢量去噪方法
CN112578428A (zh) * 2020-11-20 2021-03-30 中国矿业大学 一种基于巷道垂向虚拟测线的散射多波超前探测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100460893C (zh) * 2005-06-03 2009-02-11 中国石油大学(北京) 基于皮尔森体系独立分量分析地震去噪的方法
CN100349011C (zh) * 2005-06-03 2007-11-14 中国石油集团东方地球物理勘探有限责任公司 地震数据处理中压制与激发源无关的背景噪声的方法

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102338886A (zh) * 2010-07-14 2012-02-01 中国石油天然气集团公司 一种有效衰减三分量地震记录中面波的极化滤波方法
CN102338886B (zh) * 2010-07-14 2013-09-04 中国石油天然气集团公司 一种有效衰减三分量地震记录中面波的极化滤波方法
CN101915939B (zh) * 2010-08-17 2012-07-18 中国科学院地质与地球物理研究所 一种面波压制方法
CN101915939A (zh) * 2010-08-17 2010-12-15 中国科学院地质与地球物理研究所 一种面波压制方法
CN102087367B (zh) * 2010-11-19 2012-10-03 中冶集团武汉勘察研究院有限公司 纵波-面波速度比值法划分基岩风化带的方法
CN102087367A (zh) * 2010-11-19 2011-06-08 中冶集团武汉勘察研究院有限公司 纵波-面波速度比值法划分基岩风化带的方法
CN103091714A (zh) * 2011-10-28 2013-05-08 中国石油化工股份有限公司 一种自适应面波衰减方法
CN103091714B (zh) * 2011-10-28 2015-11-18 中国石油化工股份有限公司 一种自适应面波衰减方法
CN103217709B (zh) * 2012-01-18 2015-09-23 中国石油天然气集团公司 一种提高地震数据信噪比和分辨率的面波衰减方法
CN103323877A (zh) * 2013-05-30 2013-09-25 吉林大学 一种基于地震勘探环境噪声指向性的去噪方法
CN103323877B (zh) * 2013-05-30 2015-05-27 吉林大学 一种基于地震勘探环境噪声指向性的去噪方法
CN104459769B (zh) * 2013-09-22 2017-06-20 中国石油化工股份有限公司 一种地震图像增强方法
CN104459769A (zh) * 2013-09-22 2015-03-25 中国石油化工股份有限公司 一种地震图像增强方法
CN104849752A (zh) * 2015-05-18 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 获取地震记录中的转换横波的高精度反褶积算子的方法
CN106353791A (zh) * 2015-07-17 2017-01-25 中国石油化工股份有限公司 基于波形特征的多波多分量资料联合属性储层预测方法
CN105652322A (zh) * 2016-01-07 2016-06-08 中国科学院地球化学研究所 多分量地震数据的t-f-k域极化滤波方法
CN106249288A (zh) * 2016-08-02 2016-12-21 中国地质大学(北京) 基于Shearlet域的极化滤波面波压制方法
CN107966734A (zh) * 2017-09-22 2018-04-27 中国地质大学(北京) 多分量地震数据的矢量去噪方法
CN107966734B (zh) * 2017-09-22 2019-04-02 中国地质大学(北京) 多分量地震数据的矢量去噪方法
CN112578428A (zh) * 2020-11-20 2021-03-30 中国矿业大学 一种基于巷道垂向虚拟测线的散射多波超前探测方法

Also Published As

Publication number Publication date
CN101644782B (zh) 2011-12-07

Similar Documents

Publication Publication Date Title
CN101644782B (zh) 基于极化滤波的多波多分量地震资料的去噪方法
Fedi et al. Inherent power-law behavior of magnetic field power spectra from a Spector and Grant ensemble
CN102053261B (zh) 一种地震数据处理方法
CN102778693B (zh) 一种基于反射波层拉平提取并消除的绕射波分离处理方法
CN103376464A (zh) 一种地层品质因子反演方法
CN101598803B (zh) 一种直接得到转换波叠加剖面的方法
Wang et al. Up/down and P/S decompositions of elastic wavefields using complex seismic traces with applications to calculating Poynting vectors and angle-domain common-image gathers from reverse time migrations
CN106547020B (zh) 一种地震数据的保幅处理方法
US11209566B2 (en) Nonlinear signal comparison and high-resolution measurement of seismic or acoustic wave dispersion
CN107894613A (zh) 弹性波矢量成像方法、装置、存储介质及设备
CN103630932A (zh) 一种地震数据分形保幅方法
CN102262243B (zh) 一种滤波法可控震源地震数据谐波干扰压制方法
CN107884829A (zh) 一种联合压制浅海obc地震资料多次波的方法
CN103576197B (zh) 一种转换波角道集抽取方法
CN103675911A (zh) 一种基于纵波和转换横波联合反演截距和梯度的方法
CN105510975B (zh) 提高地震数据信噪比的方法及装置
CN106970417A (zh) 椭圆展开转换波速度分析方法与系统
CN102313903A (zh) 基于波动方程外推算子的vti介质中叠前时间偏移方法
CN102338888B (zh) 提高地震分辨率的可控震源资料相关方法
CN102053260B (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN113534243B (zh) 一种被动源Marchenko成像方法及系统
CN104635264A (zh) 叠前地震数据的处理方法及设备
CN102692643B (zh) 时变可控震源力信号反褶积方法
CN103777237A (zh) 一种基于空变加权镶边波数域滤波的地表高程平滑方法
CN104007466A (zh) 一种基于纵波振幅的无井约束叠前反演实现的储层与流体预测方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20111207

Termination date: 20190825

CF01 Termination of patent right due to non-payment of annual fee