VSP波场分离的径向中值滤波方法
VSP技术综述

VSP技术综述1前言垂直地震剖面技术(简称VSP技术)是一种垂直地震剖面是一种地面激发、井中接收的地震观测技术。
与地面地震相比, VSP技术中,地震波少经过一次地表低速带,其得到的地震资料的信噪比要好, 分辨率高, 波的运动学和动力学特征更加明显。
本文综合了一部分前人的研究成果,简要介绍了vsp技术的原理、采集和处理等方面的内容,并阐述了VSP的一些优缺点。
1.1研究目的及意义常规地震勘探是在地面激发地震波、地面布置检波器接收的一种勘探手段,这种勘探手段所得的剖面是常规地震剖面。
随着油田勘探开发难度的增大,常规勘探手段所得到的地震资料精度已经无法满足勘探的需要。
因此,出现了在地面激发、井中接收,利用直达波和反射波研究井旁构造和岩性的地震勘探方法。
这种方法就是垂直地震剖面法,简称为VSP(Vertical Seismic Profiling)方法。
VSP技术是一种检波器沿井孔放置,在地层内部接收地震波的方法。
与地面地震相比,VSP资料具有信噪比高、分辨率高、波的运动学和动力学特征明显等优点。
由于VSP观测系统中接收到的地震记录只穿过一次低降速带,地震波能量特别是高频成分相对于地面地震损失减少,具有更高的分辨率;VSP记录中既包含上行波,又包含下行波,波场信息丰富;VSP技术提供了地下地层结构同地面测量参数之间最直接的对应关系,可以为地面地震资料处理解释提供精确的时深转换及速度模型,可以可靠地识别地震反射层的地质层位,改善地面地震资料的解释效果,甚至可以利用VSP资料直接研究岩性和储层物性。
所以,VSP技术是一种很有前途的地震观测技术,研究VSP技术的理论及应用也有很重要的实际意义[3]。
1.2国内外研究及应用现状40 年代,一些前苏联科学家研制了体系完全的VSP野外采集系统及其相应的处理、解释理论,这使VSP 技术发展成为了一套完整的、独立的、新颖的观测体系。
在1973 年,加尔彼林院士出版了专著《垂直地震剖面》,这本书对前苏联十多年的研究工作做了很好的总结,为VSP 技术的发展奠定了坚实的基础[3]。
三维VSP技术应用和发展现状简述

实用文档能源地球物理方法----三维VSP技术应用和发展现状简述地测学院S0*******孟庆鑫前言 (1)一方法简介1.1 方法特点简介 (1)1.2 三维VSP方法 (1)1.2.1 3D VSP技术的应用现状 (1)1.2.2 3D VSP方法可解决的地质及地球物理问题 (2)1.3 几种相关方法简介 (2)1.3.1 逆VSP(R—VSP)技术 (2)1.3.2 多分量VSP技术 (2)1.3.3 套管井VSP技术 (2)1.3.4 高频VSP与地面地震资料综合应用 (3)二. 方法应用2.1 采集 (3)2.1.1 井下地震信号接收系统 (3)2.1.2 三维VSP观测系统 (3)2.1.2.1 概述 (3)2.1.2.2 具体实施 (4)2.2 处理 (4)2.2.1 三维VSP数据处理关键步骤 (4)2.2.2 主要处理流程和技术 (5)1.选择三维VSP资料处理实际工区范围 (5)2.炮点静校正采用直达波常规静校正方法 (5)3.三分量合成不同类型的地震波具有不同的偏振特性 (6)4.波场分离与反射波信号加强 (6)5.建立纵波与转换波地震数据体 (6)2.2.3 具体方法 (6)3D VSP波场分离处理 (7)3D VSP参数提取 (7)3D VSP成像 (8)三. 应用实例3.1 采集方法实例 (9)1. 工区概况 (9)2. 采集方法论证与设计 (10)3.2 处理方法实例 (11)3.3 对处理问题的认识和总结 (13)四、结束语 (14)三维VSP技术应用和发展现状简述前言垂直地震剖面(VSP)是一种井中地震观测技术。
苏联学者加里别林将VSP方法作为实际介质中地震波传播的实验研究方法提出来,该方法中包括早期用于测量地震速度和解决构造问题的地震测井和井中地质勘探。
加里别林对VSP方法最初的定义是“测量器件(震源或接收器)中的一种位于井中并在其中移动,而另一种位于地表或在另外一口井中进行测量的任意一种地震观测”。
地震勘探新方法作业题

地震勘探新方法作业题01综述1、写出5种与常规地面采集(地面激发、地面接收,主频20-40Hz)不同的地震勘探新方法新技术。
VSP:地面激发、井中接收(零偏、非零偏、Walkway、三维)井间地震:井中激发、井中接收时移地震/四维地震:多次采集随钻VSP:钻头激发多波多分量:纵波、横波激发(山地地震高分辨率采集高密度采集)2、写出地震勘探中5种解释新方法。
属性分析、地质统计学、反演:叠后反演、叠前反演(EI)、AVO、裂缝预测、信息融合技术、神经网络3、写出5种地震勘探基础理论新方法。
反演理论、小波变换、神经网络、模糊聚类、图形图像学、地震波模拟(数值模拟;物理模拟)、各向异性02 VSP1、什么是VSPVSP:垂直地震剖面,是一种井中地震观测技术。
也即在地面激发、井中放置检波器接收地震信号的一种地震观测技术。
2、VSP的采集方式(VSP的采集方式是指激发点、接收点的排列特点和空间分布特征)地面多次激发,井中三分量接收,激发-检波器提升-再激发-再提升。
3、VSP分为哪几种采集方式(三种)按激发点、接收点的分布特征可以将VSP的采集方式分为①常规VSP采集;②长排列资料采集;③三维VSP与三维地震联合采集4、零偏移距VSP有哪些应用求取各种速度、识别地面地震剖面上的多次波、标定地质层位、计算井旁的Q衰减因子等。
5、偏移距(非零偏)VSP有哪些应用查明井旁的地层构造细节、其作为二维观测可以作出一小段局部地震剖面,具有很高的垂向和横向分辨率描述井旁一定距离内的构造和岩性变化。
附:VSP应用:提取准确的速度及时深关系(零偏)标定地震地质层位(零偏)多次波的识别(零偏)提取反褶积因子预测井底下反射层的深度计算吸收衰减系数提取纵横波速度比及泊松比等参数6、在VSP中,什么是上行波和下行波。
直达波是上行波还是下行波,一次反射波是上行波还是下行波向下传播到达检波器的波/来自接收点上方向下传播的波称为下行波;向上传播到达检波器的波/来自接收点下方向上传播的波称为上行波。
非零偏VSP波场分离的量化质控

非零偏VSP波场分离的量化质控摘要:保幅的波场分离是vsp后续属性研究的关键步骤,也是vsp准确成像的基础处理流程。
文章引用参数反演的方法对非零偏vsp的波场进行分离,其中反演出的入射角和视慢度参数在波场传播过程中具有实际物理意义,利用方法本身的可逆性,将分离出的单纯波场对原始波场进行重构,再用重构波场与原始波场进行比较,从而达到对波场分离的保幅效果进行量化质控的目的。
应用该方法对实际资料进行了波场分离测试,并对原始资料进行了重构,验证了保幅的效果,实现了量化质控的振幅保真的单纯波场分离。
关键词: vsp 三分量波场分离量化质控一、引言在地震数据处理中,分离相互干扰的不同波场是一个很关键的问题,而vsp观测系统同时记录了地面激发产生的下行波和上行波,尤其在有偏移距的vsp观测系统下,能同时观测到下行p波、下行p-s波、上行反射p-p波和上行反射p-s波。
用于vsp波场分离方法很多,如f-k滤波、radon变换和中值滤波等,都属于标量波场分离方法。
非零偏vsp波场以一定角度传播到井中水平和垂直三分量检波器,这个角度随着偏移距离和传播时间的变化而变化,需要引入矢量波场分离方法对其进行分离。
esmersoy在1988年将反演的方法引入到两分量vsp的矢量波场分离中来,他通过最优化的方法估算出下行p波和sv波的视慢度和极化角参数,然后将参数代入方程求解波场。
随后这种技术被扩展到能同时分离上行p波和sv波并于下一年发表(leaney和esmersoy,1989)。
本文引用该反演方法,通过射线追踪正演约束反演参数,使反演方法更为稳定,参数更快收敛,提高程序运行效率,并着重利用参数反演方法本身的可逆性对分离出下行和上行p 波和sv波对原始波场进行重构,与原始输入波场进行相减,验证保幅效果。
二、原理参数反演波场分离方法假设各级井中三分量检波器接收到的波都是平面波,不同平面波场在时间域的时移,表现在频率域中就是线性的相位移动:然后再对这些频率成分进行逆fourier变换到时间域,就可以分离出假设的各个波场。
Vsp记录的波场分离方法综述

假定地震数据为 d ( x , t ) , 其拉冬变换为 m ( p , τ) , 则它们之间有如下的关系:
参考文献 [1] 朱光明. 垂直地展剖面法[M]. 石油工业出
版社, 1 9 8 8 . [ 2 ] 尧德中. 弹性波成像与纵横波分离[ D ] . 成
都地质学院博士论文, 1 9 9 1 . [ 3 ] 詹正彬, 等. 一种分离纵横波的方法[ J ] .
石油地球物理勘探, 1 9 9 0 , 2 5 ( 3 ) : 2 8 6 ~ 295. [ 4 ] 陆基孟. 地震勘探原理[ M]. 东营: 石油大 学出版社, 1 9 9 3 . [ 5 ] 涂齐崔. V S P 射线追踪模拟与处理方法 研究及应用[ D ] . 中国海洋大学研究生论 文, 2 0 0 6 .
4
科技资讯 SCIENCE & TECHNOLOGY INFORMATION
( 3 ) 在径向分量上, 直达 P 波和上行 P - S V 波近似在同一偏振平面内, 其视速度方 向相反。
( 4 ) 在切向分量上, P 波上行反射波和下 行转换波近似在同一偏振面内, 且具有相 反的视速度方向。
( 5 ) 根据各分量视速度方向上的差异, 利用 F - K 等方法将上述 4 种波场进行分 离, 得到上行 P 波和上行 P - S V 波。
(10)
(11) 为了方便, 将式( 1 1 ) 记为矩阵形式
(12)
称 L 为拉冬变换算子。为了克服离散 化和截断效应对拉冬变换的影响, 可用最 小平方反演算法求取双曲拉冬变换, 即建 立目标函数:
VSP地震勘探技术

9 or 13 Level (MLR)
Sample Rate: 1/4, 1/2, 1, 2 ms
Temperature Rating: 200 C
Pressure Rating: 20,000 psi
Max. Wireline Length 30,000 ft.
井下检波器应具条件:
(1)两端呈流线型尖端—避免管道波产生; (2)直径要小—避免井筒波,适应性更强; (3) 配备可伸张的推靠臂—便于检波器在井中移动;保证与井壁具有良
(1) 零井源距观测系统(零偏)
激发点与井口的水平距离d小
于150米的称之为零偏井源距 VSP。所谓的观测系统是指炮点
d
与接收点的相对位置关系.每激发
一次井下检波器由井底向上提升
一次 。
零井源距观测系统的作用:求 取地层速度、进行波场分析、制 作VSP地震道、预告未钻遇层位、 联结地面地震、测井曲线及地质 剖面 、为地面地震提供子波、处 理与解释的各种参数等。
井筒波的基本特征:强 度高,振幅不随深度衰减; 频谱宽,在高频范围内观 测时,沿流体柱方向有波 散;速度低(14001460),在记录上与横波 记录区重叠;可以有入射、 反射等多种类型。几种压 制方法:降低液面高度, 增加震源偏离观测井的距 离,压制高频低速波等。
2、电缆波 :电缆振动引起检波器振动 电缆波是一种因电缆振动引起的噪声。电缆波的速度与电缆结构有关。电缆波在记
(3)特殊复杂设计的VSP占用井场时间长,经 费开支大。与地面地震相比检波器组合级数少,叠 加次数低。
二、VSP资料采集
(一)、VSP野外采集装 备
VSP野外采集装备包括:井口震源、井下检波器、记录仪 器、电缆、参考检波器(近场检波器)
中值滤波法

中值滤波法
中值滤波法是一种基于排序统计理论的图像处理方法,其目的是去除图像中的噪声。
该方法的核心思想是将像素周围的像素值进行排序,然后取中间值作为该像素点的取值。
中值滤波法是一种非线性滤波方法,相对于线性滤波方法,它在去除噪声的同时可以保留更多的图像细节信息。
中值滤波法适用于各种类型的噪声,例如椒盐噪声、高斯噪声等。
中值滤波法的原理比较简单,但是计算复杂度较高,尤其是对于大尺寸的滤波器。
因此,在实际应用中,需要根据图像的噪声类型以及处理效果的要求来确定滤波器的尺寸。
总之,中值滤波法是一种简单有效的图像处理方法,在很多场合下都有着广泛的应用。
- 1 -。
VSP井地联合地震勘探技术应用——以华北束鹿泥灰岩勘探为例

VSP井地联合地震勘探技术应用——以华北束鹿泥灰岩勘探为例陈沅忠;张锐锋;唐传章;张彦斌;万学君;王刚【摘要】针对束鹿凹陷深部泥灰岩地震分辨率低、非均质性强的勘探难点,在该区开展了八方向Walkaway-VSP和全方位三维地震联合勘探技术应用研究.基于零井源距VSP和WalkawayVSP资料求取准确的地层速度、真振幅恢复TAR因子、品质因子Q及各向异性等地球物理参数,并应用到全方方位高密度三维地震数据处理,获得了更高精度的成像数据体.基于井驱数据处理结果提供的建议井位获得高产油气,表明VSP井地联合勘探技术在油气田勘探开发中具有良好的应用前景.【期刊名称】《石油地球物理勘探》【年(卷),期】2018(053)0z2【总页数】8页(P50-57)【关键词】VSP;速度;TAR因子;Q因子;各向异性;逆时偏移;井驱处理【作者】陈沅忠;张锐锋;唐传章;张彦斌;万学君;王刚【作者单位】东方地球物理公司新兴物探开发处,河北涿州072750;电子科技大学信息与通信工程学院,四川成都,611631;中国石油华北油田公司勘探事业部,河北任丘062550;中国石油华北油田公司勘探事业部,河北任丘062550;东方地球物理公司新兴物探开发处,河北涿州072750;东方地球物理公司研究院,河北涿州,072751;东方地球物理公司西南物探研究院,四川成都610041【正文语种】中文【中图分类】P6310 引言VSP采集是把检波器置于钻孔中,记录地表激发的地震信号。
相对于地面地震勘探[1-2],VSP的优点是把记录时间、记录波形与储层信息三者有机地结合起来,VSP观测到的反射波传播路径短、受近地表低速带和环境噪声影响均较小,因此VSP资料具有高分辨率、高信噪比等优势[3]。
随着井下仪器制造工艺的进步和多级井下检波器的应用及处理解释技术的发展,VSP技术应用越来越广泛,观测方式从零井源距VSP、非零井源距VSP向Walkaway-VSP(WVSP)、Walkaround-VSP、3D-VSP和井地联合勘探方向发展[4-6]。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
VSP波场分离的径向中值滤波方法聂爱兰【摘要】在VSP资料中,中值滤波法是一种简单有效的波场分离方法,主要适用于简单线性波场分离.提出并实现的径向中值滤波方法,除具有常规中值滤波方法的优势外,对具有发散状多方向线性同相轴轨迹的干扰波场,其波场分解可一次完成.通过井中地震(VSP)实例测试表明,该方法能够更有效地分离VSP上下行波场,以及P波和P-SV波分解,且在压制面波或提取面波处理方面具有潜在的应用前景.%Separating up going and down going wave field efficiently is the most fundamental and crucial part in VSP data processing. The median filtering is a simple but effective method for wave field separation, but conventional median filtering method is only suitable for simple linear wave field separation. In this paper, a radial median filtering method is proposed and realized and the new method amend the old one when radial linear interference wave field exists in data, in this situation, the new method can separate wave fields in just one time. VSP data processing test shows that this new method is more efficient in separating up going and down going wave field of VSP data, and also good in separating P and P-SV wave field of VSP data, and will have a good application in surface wave suppression and surface wave extraction.【期刊名称】《煤田地质与勘探》【年(卷),期】2011(039)005【总页数】4页(P69-71,75)【关键词】中值滤波;径向扫描;波场分离【作者】聂爱兰【作者单位】中煤科工集团西安研究院,陕西西安 710077【正文语种】中文【中图分类】P631井中VSP方法既能够高精度测定地层速度和获取精细的“时间-深度”标尺,亦能够利用VSP反子波算子及地层吸收等参量关系,高可信度地提高地面地震剖面的有效分辨率(包括消除多次波)。
VSP资料与地面地震资料相比,具有反射振幅强、信噪比高、初至波清晰可辨等优势。
因此,在复杂构造地震精细勘探中,VSP技术再度得到重视[1-6]。
在VSP资料中高质量地分离上行波和下行波是基础和关键内容。
中值滤波法是一种简单有效的波场分离方法,常规中值滤波法可分为5个步骤[1,7]:a.拾取各地震道下行直达波初至时间;b.依据各道初至时间做静态时移,将下行波排齐;c.沿深度轴按所选跨度(即道数)实施中值滤波运算,得到下行波;d.原始剖面与下行波剖面相减得到上行波剖面;e.对上行波和下行波剖面分别做反向时移,得到上、下行波剖面。
常规中值滤波方法主要适用于简单线性波场分离,本文提出的径向中值滤波方法,除具有常规中值滤波方法的性能外,对具有发散状多方向线性同相轴轨迹的干扰波场,其波场分解处理可一次性完成,并且可进一步对复杂的P波和P-SV波场进行精细分解,为多波参数分析提供可靠的波场剖面信息。
中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术。
它把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替[8]。
基本算法及步骤为:a.设有一组数(x1,x2, …, xn);b.对这N个数,按其数值大小排序;c.取重排序后的中间数值为输出:中值滤波具有算法简单,去噪明显等优点:a.中值滤波绝对阻止噪声峰值;b.中值滤波不改变阶跃函数在空间、时间上的位置;c.中值滤波对三角函数的平滑能力与滤波跨度有关;d.中值滤波对方波平滑能力与滤波跨度有关,当跨度足够长时,方波被完全平滑。
径向中值滤波法,是基于地震资料中各波场的视速度差异,沿强线性同相轴方向运用中值滤波法分离出线性干扰,进而利用“减去法”衰减强线性干扰,从而增强有效信号,得到期望的波场记录。
其实现步骤为(图1):a 对于地震道上的点(i, j)(i记录道号;j为采样点号),在记录上选取2m+1道(i为2m+1道的中心记录道号),形成该采样点的径向扫描道集;b.给定斜率扫描范围,在扫描道集里面逐点扫描,当某个扫描斜率的振幅值的和的绝对值最大时,则将该斜率定义为需分离的线性干扰的斜率;c.将该斜率记录下来的振幅数列进行中值滤波,其结果为该点的新记录;d.对所有输入道集记录逐点逐道重复b和c过程,分离出线性干扰记录;e.用输入道集记录减去被分离出的线性干扰记录,从而得到对线性干扰的衰减,完成径向中值滤波。
其中径向扫描斜率范围有正负,定义沿着顺序道时间增大时为正,反之为负。
对于VSP波场上、下行波分离,径向中值滤波最主要的处理参数一是时空域内径向扫描斜率取值范围;二是空间跨度(道数)。
时空域内径向扫描斜率取值范围,应视其下行波和上行波视速度值变化情况而定;径向空间跨度对处理效果有重要影响。
通过多组模型进行空间跨度测试分析,当跨度较小时,容易出现假频;当跨度较大时,对波场的信息都有损失。
一般道集径向空间跨度建议采用9~13道。
图2a为模拟的非零偏VSP垂直分量记录,观测井段为200~1 700 m,道距10 m,偏移距300 m。
对模型运用径向中值滤波法是以分离出下行波场得到上行波场为目标。
依据原始记录中下行波场的视速度和测试分析,径向扫描斜率取值范围为1~10 ms/道,空间跨度取值为11道效果最佳。
图2b为处理结果,下行波得到很好的分离,和图2a原始记录中相位和能量保持了一致;图2c的上行波场可为后续VSP处理提供可靠信息。
图3a为非零偏VSP原始Z分量记录, 观测井段1 500~4 560 m,道距20 m,偏移距2 000 m。
由该记录剖面可见,除常规目标信号P波(下行直达P波、上行反射P波)之外,与地层岩性密切相关的转换横波P-SV在此记录剖面中也具有可识别的品质特征(见图中标注)。
下面以该实测VSP记录为例,就径向中值滤波方法用于VSP记录上、下行波分离,以及P波和P-SV波分解处理步骤和测试结果作如下简述。
图3a中视下行波场为干扰,需被分离,径向扫描斜率范围取2~15 ms/道,空间跨度为11道。
图3b为对应的下行波分离结果。
图3c是图3a记录“减去”3b 记录所得的上行波记录剖面。
由图3处理结果可见,图3b中几乎没有上行波成份,下行P波和P-SV波能量在记录剖面上得到了很好的增强,并且与图3a相对比,两者的下行P波初至和下行P-SV波的相位特征具有高保真度及较好的一致性,这一特点相对其他二维滤波方法有一定优势(如F-K方法)。
图3b中下行P波和下行P-SV波,特别是图3c中的上行P波和上行P-SV波相互重叠,即便是不提取转换波信息(P-SV)而只考虑纵波(P)问题时,对下行波子波提取、上行波成像处理及解释都将造成一定的影响。
因此,将P波和P-SV波进一步分解是精细波场分离的目标之一。
对于上、下行波场中P波和P-SV波分离,可依据图3a所示Z分量记录中纵波和转换横波的上、下行波视速度(或时差)的取值区间(差异)进行分离。
图4是图3b对应的径向扫描斜率范围取8~15 ms/道、空间跨度为11道,下行P波和下行P-SV波的分离结果。
图5是图3c对应的径向扫描斜率范围取8-15 ms/道、空间跨度为11道,上行P波和上行P-SV波的分离结果。
由图4可见,下行P波初至和下行P-SV波能量得到了增强,并且两者相对比,下行P波初至和下行P-SV波的相位特征具有高度一致性。
图4b中下行P波初至被有效分离,下行P-SV波与下行P波初至的交点更清晰。
图5a中,浅层三组上行P波同相轴波组和图3b中的相位一致,并且能量得到了增强;图5b中上行P波被分离出去,上行P-SV波能量得到增强,极易识别。
分离出来的P波和P-SV波场记录都可以提取各自的波场信息,这为后续的VSP处理提供可靠资料。
理论和实际资料处理结果表明,径向中值滤波法除了能保持中值滤波特性之外,还具有如下特点:a.径向中值滤波不需拾取地震道下行直达波初至时间,也不需依据各道初至做静态时移将下行波排齐。
b.常规中值滤波法是沿深度轴实施中值滤波,是一维滤波运算;径向中值滤波是在道方向以及时间方向上进行中值滤波,是二维滤波运算。
c.常规中值滤波法最基础的工作是在拾取地震道下行直达波初至时间之后进行的,对资料的信噪比要求高;而径向中值滤波法只要地震资料线性干扰的视速度可识别就可以实施,相较而言,对地震资料的信噪比要求就低些,比常规中值滤波法有更广的适用性。
d.常规中值滤波法在需要消除某组下行波时,先将该组下行波排齐后再进行中值滤波,要消除多组下行波就需进行多次中值滤波运算;而径向中值滤波只需一次即可完成不同视速度的多组下行波分离。
e.反射地震资料中,面波一般都视为干扰,在处理时都对其进行压制;而径向中值滤波法可以将地震资料中的面波分离出来,分离出来的波场有利于研究地表层低降速带结构[9-10]。
径向中值滤波法是一种非线性滤波技术,不仅对VSP资料,对上下行波场分解处理可一次性完成,亦可进一步对复杂的P波和P–SV波场进行精细分解,为多波参数分析提供可靠的波场剖面信息;对叠前地面地震反射资料,可分离任意方向的线性干扰,信噪分离能力强、振幅保真度高,是反射资料高分辨处理的理想方法。
【相关文献】[1] 朱光明.垂直地震剖面方法[M].北京:石油工业出版社,1988.[2] 郭建.VSP技术应用现状及发展趋势[J].勘探地球物理进展,2004,27(1):1–8.[3] 涂齐崔.VSP射线追踪模拟与处理方法研究及应用[D].青岛:中国海洋大学,2006.[4] 于海超.非线性滤波技术在地震信号处理中的应用[D].大庆:大庆石油学院,2007.[5] 唐金良,曹辉.中值滤波在井间地震资料处理中的应用[J].石油物探,2005,44(1):47–50.[6] 夏洪瑞,董江伟.环形非对称中值滤波技术在可控震源资料去噪处理中的应用[J].石油物探,2009,48(5):510–514.[7] 李承楚.关于中值滤波的理论基础[J].石油地球物理勘探,1986,21(4):372–379.[8] 胡光瑞,郑志航,戚飞虎.二维数字信号处理Ⅱ[M].北京:科学出版社,1985.[9] 周熙蘘,王利国.利用瑞利面波调查表层结构的可行性探讨[J].石油地球物理勘探,2004,39(2):181–186.[10] 成锁.近地表调查中的面波方法与利用井资料提高分辨率的技术[D].长春:吉林大学,2008.。