CN103605157A - 衰减近地表散射波的方法 - Google Patents

衰减近地表散射波的方法 Download PDF

Info

Publication number
CN103605157A
CN103605157A CN201310479157.2A CN201310479157A CN103605157A CN 103605157 A CN103605157 A CN 103605157A CN 201310479157 A CN201310479157 A CN 201310479157A CN 103605157 A CN103605157 A CN 103605157A
Authority
CN
China
Prior art keywords
wave field
wave
big gun
near surface
gun collection
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
CN201310479157.2A
Other languages
English (en)
Other versions
CN103605157B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina 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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201310479157.2A priority Critical patent/CN103605157B/zh
Publication of CN103605157A publication Critical patent/CN103605157A/zh
Application granted granted Critical
Publication of CN103605157B publication Critical patent/CN103605157B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种衰减近地表散射波的方法,包括以下步骤:输入能量均衡的炮集波场和直达波场;根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理,获得时域近地表散射波场;用多道自适应匹配滤波技术将所述时域近地表散射波场从所述炮集波场中滤除。本发明可以预测近地表多种类型的散射波,包括侧面散射波,从而能有效预测并衰减复杂不均匀地表产生的近地表散射波,提高了复杂山地地震数据的体波信噪比,从而有利于获得高质量的地震勘探数据,因此,本发明适用于砾石区和山地等复杂近地表区地震资料处理。

Description

衰减近地表散射波的方法
技术领域
本发明涉及地球物理油气勘探中地震数据处理技术,尤其涉及一种衰减近地表散射波的方法。
背景技术
地震勘探是一种利用地下介质弹性和密度的差异,通过观测和分析大地对人工激发地震波的响应,推断地下构造形态和岩层性质的方法。地震勘探是钻探前勘测石油与天然气资源的重要手段,在煤田和工程地质勘察、区域地质研究和地壳研究等方面,地震勘探也得到广泛应用。地震勘探在分层的详细程度和勘查的精度上,都优于其他地球物理勘探方法。
地震勘探通常分为折射波地震勘探和反射波地震勘探,目前在石油和天然气资源勘查领域,主要采用反射波地震勘探,由地震数据采集、地震数据处理和地震资料解释三个阶段组成。在反射波地震勘探中,地下一次反射波是有效信号,地震数据处理的目的是提高一次反射波的信噪比和分辨率,而其他地震波都认为是噪声,所以噪声压制成为地震数据处理的首要环节。地震波从激发到接收,一般要两次经过近地表。当近地表由不均匀介质组成时,地震波要经历两次散射,这种散射就属于近地表散射。更严重的近地表散射是不均匀近地表对沿地表传播的地震波(如直达波、折射波和面波)的散射作用,由于这类波能量强,导致这种近地表散射影响比较明显。近地表散射波仅仅是由于近地表介质不均匀性造成的,与勘探目的层无关或关系不大,但是却严重影响了地震资料的品质,降低了地震资料的信噪比。
目前已经研发了多种噪声压制技术,主要分以下五类:(1)采集和处理相结合压制噪声,如采集上炮点组合和接收点组合;(2)利用频率分布特性压制噪声,如滤波法;(3)利用视速度差异压制噪声,如折射波和面波压制;(4)利用传播规律压制噪声,如压制多次波和切除直达波;(5)利用波形差异压制噪声,如调幅和调频法。这些方法常常把近地表散射波作为随机噪声进行处理,然而近地表散射波有其规律,而非随机噪声。为了压制近地表散射噪声,用逆散射方法衰减瑞雷散射波(瑞雷散射波实质是一种面波散射波),然而该方法反演计算量大,并且这种方法只能衰减瑞雷波散射波,而近地表散射波不仅包括瑞雷散射波,还包括直达散射波和折射散射波等。为了解决复杂近地表对地震波的影响,虚源法将检波器放置在复杂近地表下方的水平井中,但由于地形起伏剧烈的山峰、沟壑纵横的黄土塬和山前砾石区等因素严重影响了地震数据的质量,这种方法在山区实施起来非常困难。
发明内容
本发明的目的在于提供一种衰减近地表散射波的方法,以衰减复杂近地表散射波,提高复杂山地地震数据的体波信噪比。
为达到上述目的,本发明提供了一种衰减近地表散射波的方法,包括以下步骤:
输入能量均衡的炮集波场和直达波场;
根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理,获得时域近地表散射波场;
用多道自适应匹配滤波技术将所述时域近地表散射波场从所述炮集波场中滤除。
本发明首先根据直达波地震干涉测量技术对能量均衡的炮集波场和直达波场进行处理,获得时域近地表散射波场,其次,采用多道自适应匹配滤波技术将时域近地表散射波场从炮集波场中滤除,当用直达面波作为直达波场时,可以预测散射面波;当用初至波波场作为直达波场时,可以预测直达散射波和折射散射波,即本发明可以预测近地表多种类型的散射波,包括侧面散射波,从而能有效预测并衰减复杂不均匀地表产生的近地表散射波,提高了复杂山地地震数据的体波信噪比,从而有利于获得高质量的地震勘探数据,因此,本发明适用于砾石区和山地等复杂近地表区地震资料处理。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1为本发明实施例中衰减近地表散射波的方法的流程图;
图2为本发明实施例中x和y的对应关系图;
图3为本发明实施例中输入的能量均衡的炮集波场图;
图4为本发明实施例中输入的直达波场图;
图5为本发明实施例中采用直达波地震干涉测量技术处理后得到的时域近地表散射波场图;
图6为本发明实施例中采用多道自适应匹配滤波技术滤除掉时域近地表散射波后的炮集波场图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
下面结合附图,对本发明的具体实施方式作进一步的详细说明。
请参阅图1所示,本发明实施例的衰减近地表散射波的方法包括以下步骤:
步骤S1、输入能量均衡的炮集波场和直达波场。能量均衡的炮集波场和直达波场由以下步骤获得:
(1)、对输入的地震炮集数据进行预处理。具体的,首先输入地震炮集数据,其中每个炮集按道号顺序排列,其次,根据地震炮集数据中的测量信息设置地震道道头字,形成米桩号道字头;最后,通过逐炮检查共炮集记录剔除废炮废道,并用陷波技术消除单频工业干扰,再去除初至波之前的背景波场。其中:
在地震采集过程中,每放一炮,一般有多道接收,那么每炮就有多道记录,这些记录组合起来被称为共炮集,整个工区内所有共炮集组合起来就构成地震炮集数据。
米桩号就是一种以米为单位的桩号。由于野外施工时用大地坐标表示地震测线位置,而实际测线往往不是正南正北,若换算为米桩号后,即把测线旋转为正南或正东方向,便于处理计算。所以为了地震资料处理方便,常常把野外记录的炮点和接收点大地坐标转换为米桩号。
废炮废道是指报废的炮和报废的道。在施工中,有时会因为激发问题,导致整炮记录中不包含地下有效信息,这种记录就是废炮记录;有时因为接收问题,导致该道记录无效,这道记录就叫做废道记录。处理前必须人工去掉这些无效的记录。判断废炮废道的依据是:记录中是否包含地下有效信号,若没有,那就是废道,如果整炮记录所有道都是废道,那么该炮就是废炮。
(2)、将预处理后的地震炮集数据进行能量均衡处理,获取炮集波场。能量均衡处理既消除了接收因素差异影响,又消除了激发能量差异影响。
(3)、从炮集波场中分离直达波场。例如,从炮集波场中初至波旅行时开始,加一个时窗,从而分离出初至波波场,初至波波场由直达波和折射波组成。或者,结合直达面波分布区时窗对炮集波场进行倾角滤波处理,从而分离出直达面波波场。
图3和图4分别显示了该方法输入的能量均衡的炮集波场和直达波场,这里仅显示了地震炮集数据中一个炮集,本发明实施例中输入的数据体有2111炮,每炮1440道,有三条接收线,道长为8s,时间采样间隔为2ms。
步骤S2、根据直达波地震干涉测量技术对炮集波场和直达波场进行处理,获得时域近地表散射波场。具体包括以下步骤:
(a)、输入炮集波场和直达波场,分别进行傅里叶变换,从而得到频域炮集波场和频域直达波场。
(b)、根据公式 p s 1 ( r A , r B , ω ) = 1 NS Σ r S p 0 ( r A , r S , ω ) p * ( r B , r S , ω ) p ( r B , r S , ω ) p * ( r B , r S , ω ) 对频域炮集波场及频域直达波场的近源道进行反褶积型地震干涉测量计算,其中,ps1(rA,rB,ω)为近源道频域近地表散射波场,p0(rA,rS,ω)为频域直达波场,p(rB,rS,ω)为频域炮集波场,p*(rB,rS,ω)为p(rB,rS,ω)的复共轭,ω为圆频率,rS为震源位置,rA和rB为接收点位置,NS为叠加所用的炮数,rB为离rS最近的接收点位置,rS在rB和rA的单边位置上。
本发明实施例中,近源道是指炮检距小于600m的地震道,所用的炮数约为30炮,具体所用炮数可根据实际炮数而定,其中,炮检距是指炮点到检波点的距离。
(c)、根据公式 p s 2 ( r A , r B , ω ) = 1 NS Σ r S p 0 ( r A , r S , ω ) p ( r B , r S , ω ) p ( r B , r S , ω ) p * ( r B , r S , ω ) 对频域炮集波场及频域直达波场的远源道进行褶积型地震干涉测量计算,其中,ps2(rA,rB,ω)为远源道频域近地表散射波场,其与ps1(rA,rB,ω)共同构成一个完整的频域近地表散射波场,p0(rA,rS,ω)为频域直达波场,p(rB,rS,ω)为频域炮集波场,p*(rB,rS,ω)为p(rB,rS,ω)的复共轭,ω为圆频率,rS为震源位置,rA和rB为接收点位置,NS为叠加所用的炮数,rB为离rS最近的接收点位置,rS在rB和rA之间。
(d)、将频域近地表散射波场进行傅里叶反变换,获得对应的时域近地表散射波场。如图5所示,为本发明实施例利用直达波地震干涉测量技术得到的时域近地表散射波场图,对比图3和图5,可以清晰地看出,图5中炮集波场明显加强了这些线性散射波和双曲线侧面散射波。
步骤S3、用多道自适应匹配滤波技术将时域近地表散射波场从炮集波场中滤除。具体包括以下步骤:
(1)、由预定的第一道窗参数c1、第一时窗参数n1和滤波器长度参数nf,从炮集波场中提取数组y。
由预定的第二道窗参数c2、第二时窗参数n2和滤波器长度参数nf,从时域近地表散射波场中提取数组x;其中,c1=c2,nf为奇数,n2=n1+nf-1,且n1和n2居中对齐,即n2比n1在两侧各多出(nf-1)/2的长度,如图2所示。
(2)、采用最小二乘法求取一个匹配滤波器f,使得
Figure BDA0000395477630000051
其中,*表示褶积运算。
(3)、计算y-f*x。
(4)、由预定滑动参数滑动时窗并重复上述步骤(1)~(3)直至完成全部数据计算,然后对各次计算的结果取平均即可获得衰减时域近地表散射波后的炮集波场,其中,预定滑动参数为:道窗每次滑动一道,时窗每次滑动半个时窗。
图6为采用多道自适应匹配滤波技术滤除掉时域近地表散射波场后的炮集波场图。对比图5和图6可看出,图6明显的消除了包括线性散射波、混响面波和双曲线侧面散射波等在内的时域近地表散射波场。
在步骤S3后,将衰减近地表散射波后的炮集波场进行后续的常规地震资料处理(例如叠加和偏移等),就可获得最终地震成像剖面。
本发明实施例首先根据直达波地震干涉测量技术对能量均衡的炮集波场和直达波场进行处理,获得时域近地表散射波场,其次,采用多道自适应匹配滤波技术将时域近地表散射波场从炮集波场中滤除,当用直达面波作为直达波场时,可以预测散射面波;当用初至波波场作为直达波场时,可以预测直达散射波和折射散射波,即本发明可以预测近地表多种类型的散射波,包括侧面散射波,从而能有效预测并衰减复杂不均匀地表产生的近地表散射波,提高了复杂山地地震数据的体波信噪比,从而有利于获得高质量的地震勘探数据,因此,本发明实施例适用于砾石区和山地等复杂近地表区地震资料处理。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种衰减近地表散射波的方法,其特征在于,包括以下步骤:
输入能量均衡的炮集波场和直达波场;
根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理,获得时域近地表散射波场;
用多道自适应匹配滤波技术将所述时域近地表散射波场从所述炮集波场中滤除。
2.根据权利要求1所述的衰减近地表散射波的方法,其特征在于,所述根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理,获得时域近地表散射波场,具体包括以下步骤:
输入所述炮集波场和所述直达波场,分别进行傅里叶变换,从而得到频域炮集波场和频域直达波场;
根据公式 p s 1 ( r A , r B , ω ) = 1 NS Σ r S p 0 ( r A , r S , ω ) p * ( r B , r S , ω ) p ( r B , r S , ω ) p * ( r B , r S , ω ) 对所述频域炮集波场及所述频域直达波场的近源道进行反褶积型地震干涉测量计算,其中,ps1(rA,rB,ω)为近源道频域近地表散射波场,p0(rA,rS,ω)为频域直达波场,p(rB,rS,ω)为频域炮集波场,p*(rB,rS,ω)为p(rB,rS,ω)的复共轭,ω为圆频率,rS为震源位置,rA和rB为接收点位置,NS为叠加所用的炮数,rB为离rS最近的接收点位置,rS在rB和rA的单边位置上;
根据公式 p s 2 ( r A , r B , ω ) = 1 NS Σ r S p 0 ( r A , r S , ω ) p ( r B , r S , ω ) p ( r B , r S , ω ) p * ( r B , r S , ω ) 对所述频域炮集波场及所述频域直达波场的远源道进行褶积型地震干涉测量计算,其中,ps2(rA,rB,ω)为远源道频域近地表散射波场,其与所述ps1(rA,rB,ω)共同构成一个完整的频域近地表散射波场,p0(rA,rS,ω)为频域直达波场,p(rB,rS,ω)为频域炮集波场,p*(rB,rS,ω)为p(rB,rS,ω)的复共轭,ω为圆频率,rS为震源位置,rA和rB为接收点位置,NS为叠加所用的炮数,rB为离rS最近的接收点位置,rS在rB和rA之间;
将所述频域近地表散射波场进行傅里叶反变换,获得对应的时域近地表散射波场。
3.根据权利要求1所述的衰减近地表散射波的方法,其特征在于,所述用多道自适应匹配滤波技术将所述时域近地表散射波场从所述炮集波场中滤除,具体包括以下步骤:
(1)、由预定的第一道窗参数c1、第一时窗参数n1和滤波器长度参数nf,从所述炮集波场中提取数组y;由预定的第二道窗参数c2、第二时窗参数n2和所述滤波器长度参数nf,从所述时域近地表散射波场中提取数组x;其中,c1=c2,nf为奇数,n2=n1+nf-1,且n1和n2居中对齐;
(2)、采用最小二乘法求取一个匹配滤波器f,使得
Figure FDA0000395477620000021
其中,*表示褶积运算;
(3)、计算y-f*x;
(4)、由预定滑动参数滑动时窗并重复上述步骤(1)~(3)直至完成全部数据计算,然后对各次计算的结果取平均即可获得衰减时域近地表散射波后的炮集波场,其中,所述预定滑动参数为:道窗每次滑动一道,时窗每次滑动半个时窗。
4.根据权利要求1所述的衰减近地表散射波的方法,其特征在于,所述能量均衡的炮集波场和直达波场由以下步骤获得:
对输入的地震炮集数据进行预处理;
将预处理后的地震炮集数据进行能量均衡处理,获取所述炮集波场;
从所述炮集波场中分离直达波场。
5.根据权利要求4所述的衰减近地表散射波的方法,其特征在于,所述从炮集波场中分离直达波场,具体包括:
从所述炮集波场中初至波旅行时开始,加一个时窗,从而分离出初至波波场,所述初至波波场由直达波和折射波组成。
6.根据权利要求4所述的衰减近地表散射波的方法,其特征在于,所述从炮集波场中分离直达波场,具体包括:
结合直达面波分布区时窗对所述炮集波场进行倾角滤波处理,从而分离出直达面波波场。
7.根据权利要求4所述的衰减近地表散射波的方法,其特征在于,所述对输入的地震炮集数据进行预处理,具体包括:
输入地震炮集数据,其中每个炮集按道号顺序排列;
根据所述地震炮集数据中的测量信息设置地震道道头字,形成米桩号道字头;
通过逐炮检查共炮集记录剔除废炮废道,并用陷波技术消除单频工业干扰,再去除初至波之前的背景波场。
CN201310479157.2A 2013-10-14 2013-10-14 衰减近地表散射波的方法 Active CN103605157B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310479157.2A CN103605157B (zh) 2013-10-14 2013-10-14 衰减近地表散射波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310479157.2A CN103605157B (zh) 2013-10-14 2013-10-14 衰减近地表散射波的方法

Publications (2)

Publication Number Publication Date
CN103605157A true CN103605157A (zh) 2014-02-26
CN103605157B CN103605157B (zh) 2016-03-09

Family

ID=50123397

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310479157.2A Active CN103605157B (zh) 2013-10-14 2013-10-14 衰减近地表散射波的方法

Country Status (1)

Country Link
CN (1) CN103605157B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106814397A (zh) * 2016-12-21 2017-06-09 长江大学 一种多参数联合反演计算岩石散射衰减的方法
CN106950600A (zh) * 2017-02-16 2017-07-14 中国石油大学(华东) 一种近地表散射面波的去除方法
CN108983291A (zh) * 2018-04-28 2018-12-11 中国石油天然气股份有限公司 一种近地表散射波获取方法、装置及系统
CN109143362A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 基于共散射角道集的散射波分离方法
CN111458749A (zh) * 2020-04-13 2020-07-28 中国地质大学(北京) 应用于被动源地震勘探的面波与体波的分离方法及系统
CN111781642A (zh) * 2019-04-03 2020-10-16 中国石油天然气股份有限公司 地震数据层间多次波衰减方法及装置
CN112285790A (zh) * 2020-03-30 2021-01-29 中国科学院地质与地球物理研究所 一种电磁波场衰减系数的确定方法、装置以及介质
CN117214950A (zh) * 2023-09-25 2023-12-12 中海石油(中国)有限公司上海分公司 一种多次波压制方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4890264A (en) * 1988-03-21 1989-12-26 Atlantic Richfield Company Seismic exploration method and apparatus for cancelling non-uniformly distributed noise
CN101013161A (zh) * 2007-01-15 2007-08-08 中国石油大港油田勘探开发研究院 基于叠前波场模拟的地震勘探层位标定方法
GB2469211A (en) * 2007-12-14 2010-10-06 Shell Int Research Method of processing data obtained from seismic prospecting
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4890264A (en) * 1988-03-21 1989-12-26 Atlantic Richfield Company Seismic exploration method and apparatus for cancelling non-uniformly distributed noise
CN101013161A (zh) * 2007-01-15 2007-08-08 中国石油大港油田勘探开发研究院 基于叠前波场模拟的地震勘探层位标定方法
GB2469211A (en) * 2007-12-14 2010-10-06 Shell Int Research Method of processing data obtained from seismic prospecting
US20100268471A1 (en) * 2007-12-14 2010-10-21 Kiyashchenko Denis Method of processing data obtained from seismic prospecting
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
中国石化石油物探技术研究院: "《油气地球物理技术新进展--第80届SEG年会论文摘要》", 31 May 2012, 石油工业出版社 *
黄伟传 等: "反褶积干涉成像及其在随钻地震数据处理中的应用", 《地球物理学进展》 *
黄伟传 等: "地震干涉处理方法在随钻地震资料处理中的应用", 《石油地球物理勘探》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106814397A (zh) * 2016-12-21 2017-06-09 长江大学 一种多参数联合反演计算岩石散射衰减的方法
CN106814397B (zh) * 2016-12-21 2019-08-06 长江大学 一种多参数联合反演计算岩石散射衰减的方法
CN106950600A (zh) * 2017-02-16 2017-07-14 中国石油大学(华东) 一种近地表散射面波的去除方法
CN106950600B (zh) * 2017-02-16 2019-02-19 中国石油大学(华东) 一种近地表散射面波的去除方法
CN109143362A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 基于共散射角道集的散射波分离方法
CN108983291A (zh) * 2018-04-28 2018-12-11 中国石油天然气股份有限公司 一种近地表散射波获取方法、装置及系统
CN108983291B (zh) * 2018-04-28 2019-11-08 中国石油天然气股份有限公司 一种近地表散射波获取方法、装置及系统
CN111781642A (zh) * 2019-04-03 2020-10-16 中国石油天然气股份有限公司 地震数据层间多次波衰减方法及装置
CN112285790A (zh) * 2020-03-30 2021-01-29 中国科学院地质与地球物理研究所 一种电磁波场衰减系数的确定方法、装置以及介质
CN111458749A (zh) * 2020-04-13 2020-07-28 中国地质大学(北京) 应用于被动源地震勘探的面波与体波的分离方法及系统
CN117214950A (zh) * 2023-09-25 2023-12-12 中海石油(中国)有限公司上海分公司 一种多次波压制方法、装置、设备及存储介质
CN117214950B (zh) * 2023-09-25 2024-05-28 中海石油(中国)有限公司上海分公司 一种多次波压制方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN103605157B (zh) 2016-03-09

Similar Documents

Publication Publication Date Title
CN103605157B (zh) 衰减近地表散射波的方法
CN101852863B (zh) 一种利用高精度单道频谱分析技术处理地震数据的方法
CN102112894B (zh) 用地震表面波的波形评估土壤性质
CN101334483B (zh) 一种在地震数据处理中衰减瑞雷波散射噪声的方法
CN111158049B (zh) 一种基于散射积分法的地震逆时偏移成像方法
Reine et al. Robust prestack Q-determination using surface seismic data: Part 2—3D case study
CN103592682B (zh) 一种地震波场能量均衡方法及处理装置
CN102313900A (zh) 三维地震采集观测系统的激发位置确定方法
CN101852864B (zh) 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法
CN104360388A (zh) 一种三维地震观测系统评价方法
CN104330826A (zh) 一种去除复杂地表条件下多种噪音的方法
CN103713312B (zh) 一种虚源地震观测系统的设计方法
CN103454681A (zh) 评价三维地震观测系统成像效果的方法和设备
Gouveia et al. Deeper Vs profile constraining the dispersion curve with the ellipticity curve: A case study in Lower Tagus Valley, Portugal
CN104597485B (zh) 一种微小断层检测方法及断层检测装置
CN103513279B (zh) 一种基于地震波波动方程的照明分析计算方法及计算装置
CN103592684A (zh) 一种保持空间属性信息的海量地震数据压缩方法及装置
Kumari et al. Role of site effect for the evaluation of attenuation characteristics of P, S and coda waves in Kinnaur region, NW Himalaya
Banerjee et al. Determination of seismic wave attenuation for the Garhwal Himalayas, India
CN114415235B (zh) 一种基于正演模拟的非规则地震数据采集方法
Schwenk* et al. Correlation of the backscatter analysis of surface waves method (BASW) for anomaly detection
Askari et al. Estimation of phase and group velocities for multi-modal ground roll using the ‘phase shift’and ‘slant stack generalized S transform based’methods
CN104375180A (zh) 一种地震数据处理方法、装置及系统
Cooper et al. Improving seismic data quality by reprocessing and redesign of a 3D survey in an area of chaotically scattered source-generated noise
Almayahi et al. Fault-controlled contaminant plume migration: Inferences from SH-wave reflection and electrical resistivity experiments

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant