地球物理学报  2015, Vol. 58 Issue (9): 3286-3305   PDF    
近地表反射和折射法的进展及应用
刘江平, 王莹莹, 刘震, 潘小康, 宗育泉    
中国地质大学(武汉)地球物理与空间信息学院, 武汉 430074
摘要:针对近地表物质非均质极强、各向异性明显及地形复杂等特点,系统阐述和讨论了近地表折射和反射法的国内外研究与应用进展,认为:① 综合利用纵、横波的优势,开展多波多分量联合勘探对提高浅层地震勘探的精度和分辨率具有重要作用;② 现在的浅层地震勘探主要是对地震剖面进行解释,容易忽略一些隐含的地质异常现象,属性提取技术是充分提取地震信息,进行全面综合解释的有效手段;③ 开展多层折射介质的观测系统和解释方法研究,尤其是折射层析成像研究,是提高多层折射介质成像精度的途径;④ 开展黏弹性、双相和各向异性介质的地震反射与折射波法研究是提高近地表地震勘探成像和物性参数提取精度的新思路.
关键词近地表     反射法     折射法     进展     应用    
Progress and application of near-surface reflection and refraction method
LIU Jiang-Ping, WANG Ying-Ying, LIU Zhen, PAN Xiao-Kang, ZONG Yu-Quan    
Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan 430074, China
Abstract: According to the strong heterogeneity, obvious anisotropy and complex terrain features of near-surface material, this paper systematically expounds and discusses the research and application progress of near-surface refraction and reflection method in domestic and abroad, finding that: taking advantage of P-wave and S-wave, and carrying out multi-wave and three-component exploration is important to improve the accuracy and resolution of shallow seismic exploration; the shallow seismic exploration interpretation now is mainly based on seismic profile, so it is easy to overlook some hidden geological anomalies, while attribute extraction technology is an effective means to fully get seismic information and conduct comprehensive interpretations; carrying out research on observing system and interpretation method of multilayer refraction media, especially refraction tomography, is an effective way to improve imaging precision of multilayer refraction media; studying seismic reflection and refraction method based on viscoelastic, two-phase and anisotropic media is an new idea to improve the accuracy of near-surface seismic imaging and physical parameters extraction.
Key words: Near-surface     Reflection method     Refraction method     Progress     Application    
1 引言

由于近地表物质具有非均质极强、各向异性明显、速度梯度大、泊松比高及地形复杂等特点,常规的一些适用于油气勘探等中深部地震勘探方法不再适用于浅层地震勘探.针对近地表地层特性,发展高分辨率、高信噪比的浅层地震勘探方法尤为重要.本文就近地表地震反射和折射法在数据采集、处理及解释方面的研究进展做简要阐述与讨论.

2 纵波反射

近地表地震反射波法主要基于地震波在弹性介质中的传播规律,由人工激发的地震波在向地下传播过程中,遇到弹性分界面时发生反射,通过检波器接收地震波反射信号,并由浅层工程地震仪记录,获得浅层地震反射资料.对采集的地震反射记录进行频谱分析、滤波、反褶积、速度分析、动校正叠加及时深转换等处理,以获得反映地下地质结构的地震剖面,达到近地表地质勘查的目的.

2.1 地震波的激发与观测

由于近地表地震地质条件复杂,浅层勘探的目标体多具有形态复杂、埋深浅、异常体小等特点,此时常规多次覆盖的结果多为一定范围内的平均效应,甚至会存在反射波相位畸变的问题,不利于反射波信号的识别;且常用的最佳时窗及最佳偏移距技术一般难以获得复杂目标体的有效反射波.为满足浅层地震勘探的需要,应采用合适的激发手段和合理的接收方式,即需要高频特性好、频带宽,能激发窄脉冲的震源,接收则可采用单道自激自收或点组合小排列的方式.

锤击震源具有灵活方便、频率高、成本低的优点,是近地表地震勘探中常用的地震波激发方式.由激震锤锤击置于锤击点的震板产生弹性波,要产生理想的地震波,锤击点应选在地层坚硬、地下水位浅的地方,且锤击时要保持激震锤与震板的瞬息接触.

炸药激发同样可以产生高频成分丰富的脉冲波,近地表勘探中,通常利用小药量激发产生高频脉冲来提高地震记录分辨率.针对巨厚黄土区及山前砾石堆积区,常规单炮地震波激发能量小、信噪比低的缺点,通常采用炮点组合激发的方式,而以往炮点组合参数的设计主要依赖于探区历史资料和野外实验,激发参数设计缺乏合理性.徐峰等(2011)基于相控理论,研究了组合震源地震波激发效率与排列长度、勘探深度的关系,通过对组合激发的时间或深度控制实现相位组合控制,使得发散传播的地震波集中于某个方向,激发定向地震波束,该方法在塔里木地区应用中对反射波信噪比改善具有明显的效果.等深等时激发的组合震源产生的波场向上传播时,由于传播距离短、能量强,且地表附近存在散射源,会在地震记录上形成强干扰,据此,刘福烈等(2014)提出了不等深延时组合激发方法,通过激发时间和深度控制各震源子波在空间上的叠加时差,在保证波场下传能量的同时有效衰减上传能量,压制虚反射,降低面波强度及散射噪声.

在浅层地震勘探中,常通过提高叠加次数来压制噪声,提高地震资料的信噪比,但由于多次覆盖具有低通滤波特性,覆盖次数过高会降低地震资料的分辨率.20世纪70年代发展起来的一种浅层高分辨率地震勘探方法——MINI-SOSIE法较好地解决了这个问题.将一个长时间的伪随机可控震源作为激发源,并通过参考道与记录道的相关运算得到地震记录,该方法产生的地震信号频带范围宽,且经过相关叠加后可得到较强的地震信号,由于随机噪声在随地震波传播的过程中相互干涉,在很大程度上压制了随机噪声,提高了地震记录的信噪比,该方法在一定程度上解决了常规地震勘探中为提高记录主频而降低记录信噪比之间的矛盾(杜兴忠和曹俊星,2008).城市活动断层探测与研究对合理进行城市规划建设及抗震设防具有重要作用,但由于城市探区多被第四纪松散沉积物覆盖,且地下电缆、管线等分布密集,给活动断层的探测带来巨大困难,针对城市的强干扰背景,吴子泉等(2005)通过研究及应用表明:通过可控震源,利用频率扫描方法和互相关处理技术,可有效压制城市的强背景干扰,对断层埋深进行精确定位.可控震源由于其信号频带宽且为无损激发,在地震勘探中得到广泛应用,但由于其信号能量较弱,导致可控震源记录信噪比较低.组合震源可增强信号能量,提高信噪比,但实际工作中组合台数受成本限制.姜弢等(20042008)提出了一种方向性震源——相控震源,令多个激震器按等间距排列,实际工作时,各个激震器依次激发,到达同一检波器的地震波信号会存在波程差,相控震源通过相位差来补偿波程差,使得经过反射的信号实现同相叠加,从而达到提高信噪比的目的.赵春蕾等(2013)对高精度组合式轻便小型可控震源进行了初步研究,通过野外试验表明该震源易于实现流动地震观测及震源组合,激发信号重复性好,通过相位控制可实现多台激发、信号同步叠加,适于特殊复杂的野外环境.近年来对可控震源在地震数据采集质量的研究主要集中于可控震源技术的发展,对激振器的研究相对较少,Wei等(2012)以液压激振器为例,讨论研究了激振器改进对提高采集信号的低频性能、高频性能及震源信号稳定性的重要性,为可控震源的发展提供了新方向.

依据工区近地表地震地质条件合理选择震源是获得高信噪比资料的关键.地下水位较浅,近地表土层较密实的地区,可采用锤击震源或小药量的炸药震源;近地表土层较薄,基岩面较浅时,可采用震源枪激发;而对土层疏松或有卵砾石层、地下水位较深的地区,宜采用可控震源(陈相府等,2005).近几年,针对浅层地震勘探的特点相继出现了几种新的激发方式.Rashed(2009)提出了用于浅层地震勘探的增压活塞震源,通过引燃封闭室内的气体产生爆炸推动活塞向下撞击汽缸底部,产生强大的冲击波,该震源轻便、成本低且激发能量大,通过试验表明:增压活塞震源可重复性高,易于控制,产生的地震波中含有大量的高频成分,是浅层地震勘探的理想震源.由Imhof等(2011)提出的气枪震源,比传统的锤击震源激发的能量穿透深度大,且避免了炸药震源对岩体的破坏,是一种新的地震波激发方式,该震源携带方便,利于高山地形勘探,在对拟建立的穿山公路隧道区开展地球物理探查,研究探区破碎程度和岩体硬度中取得了较好的效果.同年,Ge等(2011)提出了一种简单机械冲击系统(SMIS)激发震源,该地震波激发系统简单小巧且成本低、可重复性强、便于安装在浅孔的各个方向,并且可以产生清晰可重复的信号,在煤层试验中取得了很好的应用效果.

野外地震数据采集,经过多年的发展形成了较系统的观测方式.二维地震勘探中通常采用端点激发或中间激发的多次覆盖观测系统,复杂地区浅层三维地震常采用规则、非规则线束状观测系统及宽线观测系统,覆盖次数、偏移距大小及排列长度的选取与探测目标体的深度、大小及地表条件有关(沈阳等,2011).谢忠球和温佩琳(2000)结合最佳偏移距选择的问题及水平薄层模型三种等偏移距观测曲线的分析结论,提出了多偏移距等偏移观测方法,并通过实际应用表明,该观测方法可以提高地震记录分辨复杂目标的能力,对解决浅层及超浅层介质探测具有重要意义;王者江等(2005)提出了一种移动炮点观测系统,保持检波点不动而炮点依次向前小间距移动,达到所勘探的反射点的前移,在炮点移动一定次数后,所有检波器再整体前移,进行下一个排列的数据采集,该方法明显改善了野外施工效率;传统的基于地下水平层状介质假设的观测系统设计方法不适用于复杂构造系统,Liu等(2005)基于控制照明思想,提出了一种面向目标成像的观测系统成像方法,利用单程傅里叶有限差分波场传播算子将目标层的平面源延拓到地表,通过分析波场能量分布来确定炮点或检波点位置,朱金平等(2011)在此基础上发展了面向勘探目标成像的三维地震观测系统优化设计方法;周竹生等(2010)借鉴折射波勘探中常用的相遇观测系统,提出了一种固定单元排列接收的多次覆盖观测系统;针对沙丘地区地震数据采集存在的成本高、效率低及数据质量差等问题,Almalki和Alkhalifah(2012)提出了一种地面拖缆地震数据采集方法,实际应用表明,该方法工作效率高且采集到的信号频带宽,地震数据分辨率高,适用于沙漠戈壁覆盖区地震勘探.

2.2 反射处理解释及相关技术

地震波在传播过程中,受大地滤波作用的影响往往会引起高频成分的快速衰减及地震子波相位的畸变,且地表情况越复杂,这种高频衰减越快.为了提高地震资料的分辨率,必须进行近地表补偿.目前比较常用的近地表补偿方法主要有:谱白化、反Q滤波、多尺度联合分析及时频分析法.谱白化主要是通过展宽信号振幅谱来达到补偿频率衰减的目的,可以通过小波分解、子空间分解或Hlibert变换将信号分解,然后利用谱白化方法对分解的信号进行合理补偿,最后将处理后的信号合并得到近地表补偿后的信号,谱白化在一定程度上可以达到信号补偿的目的,但由于该方法是在分频后对各个频带数据进行增益,无法保证处理前后时间方向和空间各道上的相对振幅关系.反Q滤波是通过消除地震波在传播过程中的振幅衰减和相位畸变来实现信号补偿的,反Q滤波的准确性与Q值的准确求取密切相关.多尺度联合分析方法是指联合VSP资料、双井微地震测井资料或井间地震资料获取相应参数进行近地表补偿,由于VSP、测井资料等只有在特定的井位置处存在,限制了该方法的应用.当没有地层吸收衰减作用时,理论上对所有频率而言,同一频率的深浅层地震反射能量比应该相同,而实际上,在介质的吸收作用下,造成了各频率能量对时间的分布不同,时频分析方法就是通过引入时变因子来实现近地表补偿的.国内外对近地表补偿已开展过多项研究.Carpenter(1966)提出了一种通过滤波器相似性解计算平稳Q滤波器的方法; Hargreaves和Calvert(1991)最早开始了基于波场外推原理的近地表补偿研究; Shi等(2004)曾利用微测井资料进行低速层高频补偿;陆文凯等(2005)提出利用子空间分解技术进行谱白化处理获得高分辨率地震资料的方法;王建民等(2007)研究了在时间域利用Levinson算法及在频率域利用傅式变换求取近地表吸收补偿因子的方法;李子顺(2007)利用双井微地震测井资料计算滤波因子,并用其对地震资料进行确定性反褶积处理,有效恢复了近地表低速层、虚反射及近震源区导致的地震波衰减;张固澜(2010)丁进杰(2013)成功地将一种改进的广义S变化用于地层吸收衰减补偿;Wang(2011)研究了基于Futterman衰减模型,利用反演来实现衰减补偿的方法.

噪声会干扰地震记录中有效信号的处理与解释,常规地震勘探中,一般通过地震道叠加来压制随机噪声,但该方法耗时长及勘探成本较高.Jeng等(2009)讨论研究了自适应滤波在超浅层地震勘探随机噪声压制中的应用,联合利用最优二维中值滤波和二维自适应维纳滤波来压制反射记录中的随机噪声,该方法可有效实现叠前去噪及地滚波压制.大量研究表明,波束形成方法具有显著地压制噪声的能力(Hu and White,1998Panea and Drijkoningen,2008),在反射波法地震勘探中,针对目的层数据信噪比较低的情况,姜弢等(2012)基于波束形成原理,提出了一种基于接收阵列的时域地震波束形成方法,对常规采集到的原始地震记录通过分组及组内延时、叠加处理,得到定向地震记录,实现波束形成处理,并利用延时参数控制波束方向,然后对获得的定向地震波数据经过增益、动校正及水平叠加等处理后可获得较高信噪比的地震剖面,且可通过延时参数来选择性地提高不同深度目的层的信噪比.投影滤波是地震信号去噪中常用的方法,但其是基于地震记录水平一致性的假设前提的,Traonmilin和Gulunay(2013)对存在静校正位移的地震数据,提出了一种同时估计投影滤波算子和静校正位移的噪声压制方法,并通过合成记录和实际资料验证了该算法的去噪能力.传统的噪声压制技术有F-K滤波、线性拉冬变换、τ-p变换等,这些常用的噪声压制技术均是在单一尺度上进行的,Górszczyk等(2014)提出了一种基于Curvelet变换的多尺度噪声压制方法,利用地震信号与噪声在不同尺度、方向及空间位置上的差异,通过给定不同尺度上的阈值参数来实现噪声衰减,该方法可实现频谱重叠的有效信号和噪声的分离,在弯曲界面等复杂构造中仍具有较好的噪声压制效果.Pham等(2014)在近似模型的基础上,提出了一种基于稀疏表征的地震信号时变自适应滤波器设计框架,实现地震记录中多次反射噪声的去除.在山地不均匀地表和粗糙起伏地形的影响下,地震资料中往往存在较强的近地表散射波,降低了地震资料的信噪比,且其与体波在时间、空间和频率上重叠,干扰有效信号,甚至会产生复杂构造假象.吴希光等(2012a2012b2012c)对表层散射波进行了详细研究,总结了近地表散射波的基本性质、形态特征及一般规律,认为散射波是导致地震资料低信噪比的根本原因,并提出了通过增加检波器垂向组合来提高地震资料信噪比的方法.在Fleury等(2010)研究的基于地震干涉测量和扰动理论的散射波重建基础上,徐基祥(2014)将干涉测量技术用于分离近地表散射波,根据近地表散射波运动学规律,通过干涉测量法实现近地表散射波的相干加强,并利用多道匹配滤波技术从原炮集中消除近地表散射波,通过这种地震干涉测量技术可有效消除地震记录中近地表散射波的影响,提高资料品质,尤其适用于复杂地形及不均匀地表的山地勘探.

共中心点水平叠加是地震勘探的重要处理手段,但当地形起伏较大时,共中心点道集中的反射波同相轴已不再满足双曲线特征,而是一条高次曲线,并且会存在时间倒转现象(炮检距增大,地震波走时反而减少),此时采用常规基准面静校正和动校正已不能消除地形起伏对地震波走时的影响,难以达到同相叠加的目的,得到的水平叠加剖面上,会存在各界面反射同相轴连续性差,横向能量变化较大,及明显的同相轴畸变和错断现象,甚至会出现假同相轴特征.且地形起伏频率越大,同相轴畸变越严重,同相轴分叉、合并及错断现象越明显.据此,刘江平等(19981999)提出基于滑动基准面的静、动校正及速度分析方法,有效避免了常规水平基准面方法所存在的误差,该方法可以消除地形起伏影响,实现同相叠加,获得高质量的水平叠加剖面,提高复杂地形区勘探精度.之后,Liu等(2005)基于起伏地表下的时距曲线方程,讨论研究了非垂直射线路径下的动静校正量,并在此基础上进行速度分析,通过合成记录与实际应用表明,基于非垂直射线路径对地震数据进行校正后,可以准确提取叠加速度,该方法对复杂地形尤为适用.近年来,不少学者针对复杂地区的静校正问题进行了研究.何光明等(2006)从方法原理、方法特点及应用效果等方面对三种常用的静校正方法(折射静校正、无射线追踪层析静校正、非线性层析静校正)进行了对比分析,表明:折射静校正适用于近地表结构中存在明显折射层的情况,近地表结构复杂但不存在速度倒转和层尖灭时可选用无射线追踪层析静校正,非线性层析静校正可用于存在速度倒转、层尖灭等的复杂近地表条件,另外在不满足地表一致性要求时可采用波动方程基准面校正方法;谭昌勇等(2009)联合使用综合模型法野外一次静校正、综合寻优折射波静校正、综合寻优反射波静校正解决复杂地区的静校正问题,这种综合寻优静校正可避免初至拾取,是一种适用于复杂近地表的可靠快速的静校正方法;Barison等(2011)研究了波动方程基准面校正方法在复杂近地表条件下的应用及重要性;Haffinger和Verschuur(2012)提出了一种利用全波形方式直接从地震数据中估算基准面校正算子的方法,该方法不需要估计速度模型,且充分利用地震数据中的振幅、相位变化信息,提高了复杂近地表条件下基准面重建精度; Masjukov和Murzov(2013)假设经过短波长静校正后的地震波同相轴是局部水平的,提出了一种新的短波长静校正方法,该方法不需要进行速度分析,且比传统的基于共中心点叠加的剩余静校正方法更适于复杂近地表条件;Sun和Verschuur(2014)Sun等(2014)提出了一种自适应输入遗传算法用于解决复杂近地表问题,利用单程走时传播算子来描述地表和水平基准面之间的地震波传播,通过自适应输入遗传算法估计走时函数,并利用估计的走时对地震记录进行时间校正,在此基础上,通过基于走时算子的地表校正和速度模型更新交互实现来消除近地表影响,保留有效地质信息,该方法在随机选择和定期更新的子集里进行参数优化,计算速度快且精度高.

随着实际发展需要及勘探技术的提高,对工程勘探精度的要求越来越高,利用地震属性技术可以从地震资料中提取有效信息,揭示地震剖面中难以识别的地质异常现象,有效识别岩性和构造特征,提高覆盖层厚度、软弱夹层、隐伏构造及滑坡体等的勘探精度,为工程勘察提供更准确的地球物理信息.地震属性技术发展至今,已有多种属性参数可供选择.其中,振幅类属性主要从能量及振幅变化角度来研究地震波形特征和能量分布,可用于识别河道及三角洲砂体、划分岩层、确定不整合面、研究层序地层变化等,振幅属性的变化已被广泛用于地层研究中,可有效区分稳定地层和杂乱互层状地层;表征属性变化率的反射强度斜率属性和瞬时频率斜率属性可有效表征地层岩性纵向变化情况;能量半衰时属性则从能量角度研究能量的衰减情况,可用于判别裂缝及特殊岩性的变化等;曲率属性可用于描述微小断层,预测裂缝方向及分布;用于描述相邻道之间相似性的相关属性也可用于断层、裂缝的识别.由于属性自身的局限性和多解性,实际应用中,需要针对工程勘察的特点,结合已知地质背景及勘探目的进行属性优选与组合,剔除不相关属性,选择能够更合理、显著地反映浅层构造和岩性特征的地震属性参数进行综合解释,提高解释精度.付小明(2012)研究了一种通过统计方法研究地震数据体中的差异进行地震解释的方法,即利用灰度共生矩阵算法提取地震资料的角二阶矩(能量)、对比度、熵等纹理属性,并通过理论模型试算及工程实例应用表明,该方法可有效提高异常体识别精度,具有较高的可靠性和分辨率.Santana和Neto(2015)研究了地震属性技术在沉积物分类中的应用,从地震记录中提取地震属性,并对其进行统计分析,根据属性的一致性联合用于海底沉积物的地质分类评价,该方法可用于粗-中粒度的淤泥沉积物勘察.地震属性技术在浅层地震勘探解释中的应用具有较大发展前景.

地震反射法在工程建设、灾害防治及地质调查方面具有广泛应用.陈相府等(2005)将地震反射法用于采空区及老塌陷区勘测;石明等(2006)利用地震反射法与其他物探方法相结合进行堤防质量检测,通过实际应用表明,地震反射波法可以清楚地对堤身进行分层,确定各层厚度及获得重要的土力学分布参数;刘祖珉等(2006)将浅层地震反射法用于某重大工程备选场址的勘察与评价,以查明候选区的第四纪覆盖层及强风化层厚度、速度分布、有无隐伏断层及不良工程现象等,为工程设计提供了有效资料;高霞等(2007)曾利用振动冲击震源和工程地震仪开展海上反射波勘探,确定基岩面起伏形态及构造分布性质;谢冰等(2007)认为水电站进水口处,可将淤泥底面与进水口底板及淤泥顶面与上部水面视为两个波阻抗界面,采用合理的震源及观测系统,利用浅层地震反射波法对水电站进水口处的淤泥厚度进行了成功探测; Fradelizio等(2008)对地下水污染区开展浅层反射地震调查,研究古河道分布结构、边界及深度,为地下水污染防治提供依据;张凯和汪青松(2010)将反射波法用于采盐区塌陷坑及周边探查,确定塌陷原因并圈定安全范围,在此基础上对潜在的塌陷区域进行了成功预报;Pap and reou等(2011)研究了地震反射法在探测浅层掩埋体中的应用,并将其用于地下管道分布探测;由于采煤区附近往往存在采空区或采煤巷道,在该范围内施工时需先进行拟建场地勘探,曹江涛等(2014)采用浅层地震反射波法对该区进行工程勘察,确定地层结构,采煤巷道走向及深度、是否塌陷等,为拟建场地的施工设计提供了重要的地质信息.

2.3 特殊方法与技术

在地质条件复杂,地下水丰富,岩体松软,断裂带、造山带及溶洞发育的地方,隧道开挖过程中会存在坍塌、冒顶、涌水等潜在地质灾害的威胁,故隧道超前预报工作对工程的安全、质量等尤为重要.20世纪八、九十年代,世界各国陆续开始了对超前预报技术的研究,而中国则始于1990年初,其中反射波法是进行超前预报中最常用且准确率最高的一种方法(赵永贵,2007).由于隧道内观测空间有限,反射波孔径小,且隧道内不同方向均可形成反射波,为来自隧道掌子面前方的反射波资料判别与解释带来较大困难.最早发展起来的负视速度和水平声波剖面法利用走时反演方法推测隧道前方地质构造,仅适用于简单地质条件下的预报工作.隧道地震超前预报技术(Tunnel Seismic Prediction,TSP)利用偏移成像方法对采集的地震数据进行处理,通过对反射信号运动学及动力学特征的分析处理达到超前预报的目的,该方法在公路铁路系统上具有较多应用实例,但其观测方式及处理方法上的缺陷,造成了较多的漏报、误报事故.钟世航(2003)赵永贵和蒋辉(2009)对TSP技术存在的问题进行了详细的分析研究.隧道地震层析成像技术(Tunnel Seismic Tomographic,TST)则采用空间阵列观测方式,并利用对地下波场的方向滤波和回波提取、速度扫描、偏移成像等处理方法进行超前预报,该方法的速度分析及偏移定位精度明显提高,而且适用于复杂地质构造地区的预报工作,是现今发展较为完善的一种方法.肖启航和谢朝娟(2012)对TST技术进行了详细的分析研究并将其用于顶效隧道超前预报中,获取掌子面前方围岩波速分布及地质构造偏移图像,为隧道超前预报提供了可靠的资料.常旭等(2006)引用反射地震记录自激自收方式,提出了一种PMZO隧道观测方案,通过单道激发接收获取零偏移距地震数据,当地震道数达到一定数量(24道或48道)时,随着隧道向前开挖,每增加一个地震道,去掉一个最早的地震道,获得实时地震记录,并利用零偏移距逆时偏移算法进行成像,实现实时预报.隧道空间中的地震波场数值模拟是认识地震波传播规律,以及合理选择数据采集方式、资料处理解释方法的重要基础.鲁光银等(2011)基于一阶速度-应力弹性波方程,采用交错网格高阶有限差分进行隧道波场数值模拟,通过对垂直断层和软弱夹层模型的正演试算验证了该方法的有效性;刘江平等(2012)则将隧道作为空气介质建立地质模型,通过有限差分方法及吸收边界条件实现隧道空间的全波场二维模拟,并对具有不同界面倾角及异常体的地质模型进行正演模拟,分析波场特征,为观测系统的设计提供了更为可靠的理论依据;Tzavaras等(2012)对比分析了KPSDM(Kirchhoff Prestack Depth Migration)、FVM(Fresnel Volume Migration)及RIS(Reflection-Image-Spectroscopy)方法在隧道成像中的应用,研究表明三维RIS方法可以压制低频带散射噪声,提高高频带分辨率,而三维FVM方法可以降低反射体成像空间误差,联合利用FVM和RIS方法可有效实现隧道环境下地质构造的清晰成像;李枝文(2013)利用有限元分析软件建立隧道三维模型,对地震波在隧道中的传播及反射情况进行了分析研究; Cheng等(2014)在此基础上对不同的隧道地质模型进行了时间二阶空间四阶的有限差分正演模拟,并对获得的地震记录进行了基于P波和S波波动方程的二维逆时偏移成像,表明用于隧道勘探中的逆时偏移成像方法可实现对界面的准确成像,并能有效压制传统成像中容易出现的镜像假象,具有较高的精度;王朝令和刘争平(2012)王朝令等(2014)基于二维有限元数值方法,模拟研究复杂地表条件下隧道内的全波波场特征,并基于多个数值模型,先后研究了τ-p变换和F-K变换在波场分离中的有效性及影响因素,通过波场分离来提高隧道超前预报的精度.

90年代钟世航(1997)针对浅层勘探的大偏移距易形成宽角反射等问题,提出了单点测量,极小偏移距激发-接收的地震勘探方法——陆地声纳法,详细讨论研究了该方法适用的测线布置方式、地震资料处理流程及定性定量解释方法,并通过物理模拟和数值模拟研究常见探测对象的地震剖面特征,为实际地震勘探应用提供了依据.王荣等(2013)将陆地声纳法用于闹市区断层和岩溶探查,克服了各种振动噪声干扰,结合已有的钻孔资料成功确定岩层层面,并在断层、深部溶洞探测上取得了良好效果,为陆地声纳法用于闹市、公路上的地面探查及活断层探测提供了借鉴依据.

地震映像法是近十几年发展起来的用于探测浅部介质纵横向不均匀体的有效方法,它采集的主要是近震源处的弹性波场,以相同的小偏移距逐步移动测点来接收地震波信号,目前主要是利用波的动力学特征(振幅、相位、频率)解释勘探剖面下纵、横向不均匀体.其具有以下特点:(1)数据采集简单、工作效率高;(2)小偏移距小点距采集,可以得到能量较强的弹性波,且得到的数据信噪比高,抗干扰能力强;(3)数据处理过程简单,不需要进行校正处理,避免了动校正对浅层反射波的拉伸、畸变处理;(4)解释中除可以利用反射波信息外,还可以充分利用折射波、面波及绕射波等多种有效波联合解释.地震映像法在断层破碎带、岩溶勘察、滑坡体探测、堤坝隐患探测、地下管道及障碍物探测、垃圾坑勘察等中均已得到广泛应用,并逐渐用于海洋水域地质调查中.由于水中不传播剪切波,故水域映像法不会存在面波干扰,另外其在地震波激发接收系统上与陆地映像法有所不同,水域映像法通常选用能量大,频带宽的高频冲击震源,通过单点接收直接形成地震剖面,或采用多道走航式采集,经过数据处理获得地震映像剖面,但水域映像法易形成多次反射,影响地震资料品质,故多次波压制具有重要意义.李哲生等(1998)具体阐述了水域映像法的震源方式及数据采集方式要点,并将地震映像法用于深港西部通道深圳湾公路大桥海上勘探,成功探查水下浅层地层结构;王云安和刘章捷(2001)将地震反射波法和地震映射波法用于桥址区水域工程勘探,以查明该区的地质构造、岩性分布等地质条件,提供桥型设计方案选择的地质依据;针对水域映像中的多次波问题,张学强等(2006)将改进的预测反褶积方法用于多次波压制,明显提高了地震资料品质,取得了较好的效果;李小平(2009)将水域高密度地震映像法用于舟山大连岛工程勘察,并结合工程地质测绘,对研究区的海底地形、覆盖层厚度、地质构造及不良地质作用进行了有效勘察,为确定桥位方案提供了依据.但目前地震映像法只能进行定性的解释,对地层中不均匀体的形态、埋深等尚难以确定,充分利用观测到的波场特征进行定量解释及开展三分量测量是地震映像法今后研究和发展的方向.

井间地震技术是采用井中激发、井中接收,并对采集到的地震信号通过反演成像获得地质体结构及构造分布的一种方法,由于地震波传播距离短,更接近探测目标,且可以有效避开表面疏松地层及地面干扰源的影响,故具有高信噪比、高分辨率、解释成果直观的特点,广泛用于精细构造及地质灾害调查、隐伏地质体探测等工程地质勘察中.近年来,对井间地震的方法研究与应用取得了一系列进展.王建军等(2006)将井间地震技术用于解决桥墩基础断裂问题,准确地确定了桥墩基础断裂的位置和基本情况;井间地震正演模拟是研究井间地震波传播规律和波场特征的有效手段,李万万和裴正林(2008)基于波动方程交错网格有限差分法,进行了井间水平层状模型和井间低、高速夹层模型正演研究,表明井间地震波场复杂且广角反射波占优势,低速夹层内存在低频高能量导波,可将其作为探测地下裂隙破碎带等低速层的依据;Belina等(2009)发展了一种适用于井间地震的高分辨率全波形反演方法,通过低、高速异常体模型及随机模型实验表明,相比于传统的射线走时成像,该方法更适用于非均匀地质构造井间成像;左建军等(2011)进行了井间直达波走时和反射波走时层析成像研究,表明直达波走时层析反演对层速度比较敏感,而反射波走时层析反演能相对准确地确定反射界面形状和深度,利用直达波和反射波联合层析成像,可充分利用两者优点,减少多解性,提高反演精度和分辨率;谈顺佳等(2013)将井间地震用于岩溶勘察获得该区反演成果图,通过与钻孔资料对比,圈定出岩溶位置范围及发育分布形态,为工程设计提供了可靠的依据;郭淑君等(2014)将井间地震技术应用到唐山市环境工程地质勘察中,并利用弯曲射线追踪法和LSQR算法获得井间速度分布,对井控范围内的地层构造和地基注浆效果进行研究,准确探测了研究范围内的断裂破碎带和掩埋防空洞.根据VSP和井间地震数据采集方式的兼容性,通过合理设计观测系统,可以同时获取井间地震CT数据和施测孔的逆向VSP数据体,在不增加成本的情况下,提高勘探精度,刘永华等(2005)利用这种VSP与地震CT数据一体化观测系统进行三峡大坝工程物探检测,通过VSP剖面解释结果建立初始CT速度模型,并将VSP速度参数作为已知边界与约束条件,进行联合速度反演实现井间速度成像.VSP与井间地震CT联合勘探技术作为一种精细探测技术组合,在围岩精细分类、岩体完整性检测和复杂构造现象及病害地质体预测等方面具有广阔发展前景.

现场实时数据处理是提高浅层地震反射方法精度及推广应用的关键因素,Heilmann等(2013)提出了一种基于云计算的浅层地震数据实时成像分析技术,利用云基础设施强大的计算能力和不依赖于速度模型的并行网络共反射面元叠加成像方法,通过无线数据连接将现场采集到的地震数据实时上传到远程数据处理系统,并自动进行质量控制数据分析,对满足成像要求的地震数据进行预处理并实现快速高效的叠加、速度模型建立和叠前逆时偏移,这种云计算方法可以推动浅层反射在近地表地质调查中的广泛应用.

3 纵波折射

地震波向地下传播过程中,在遇到地层界面下层速度高于上层速度时,以入射角大于等于临界角传播到该界面的波会发生折射,地震波沿折射界面滑行,并在上覆介质中产生折射波.折射波法勘探就是利用检波器接收折射波,并经过相应处理解释获得近地表结构模型,完成工程勘察任务的.

3.1 折射波观测系统与正演模拟

折射波的特性决定了折射波法观测系统与常规观测系统不同,在折射波法勘探中,往往需要根据工区具体的地表条件和地质特征选择合理的观测系统进行数据采集.单边观测系统是一种单边一点放炮、全排列接收的数据采集方式,其只适用于折射界面较浅的观测条件;相遇观测系统则是一种测线两端放炮、全测线观测的方式,在两端点分别激发地震波获得正、反时距曲线,由于两折射波射线路径相同,其走时相等,两激发点构成互换点,相遇时距曲线的互换点特性是折射波法解释中的重要依据;追逐观测系统是在同一排列上接收来自同侧不同位置激发的折射波,其是由同侧不同点放炮的单边观测系统组成,该观测方式可以依据追逐曲线的平行性特点判断目的层有无穿透现象;追逐相遇观测系统兼容了相遇观测系统和追逐观测系统的优点,适用于需要同时观测深、浅折射界面或对同一折射界面重复观测的情况,在复杂地质条件下的观测效果较好.

折射初至正演对反演成像及沙漠、山地等具有复杂近地表结构地区的静校正研究具有重要的指导作用,折射波正演主要是通过射线追踪法来实现的,最先发展起来的射线追踪方法有试射法和弯曲法,但这两种方法难以处理速度变化较大的情况,且其采用的最小走时方法具有局限性,在此基础上进行改进,陆续提出了各类波前重建方法、慢度匹配法等,实现多值走时计算,提高计算效率(张钋等,2000).戴云等(1999)提出了一种基于Fermat原理的折射波正演方法,通过寻找固定炮点、检波点之间使得折射波传播时间最短的起始折射点及出射折射点位置,来确定相应炮点到检波点的折射波走时,这种折射波正演算法可以适用于任意地形和具有复杂折射界面的模型;徐涛等(2004)对传统的射线追踪方法进行了改进,提出了一种折射波逐段迭代射线追踪方法,利用推导出的透射中间点和折射点的一阶修正公式,通过逐段迭代最终求出满足费马原理的射线路径,该方法适用于复杂介质的折射波射线追踪;涂齐催和刘怀山(2006)利用线性走时插值射线追踪开展近地表模型正演计算,该方法追踪的射线路径不再是网格离散点的简单连线,而是在满足费马原理的精确点处与边界相交,并通过线性插值的方法得到交点走时,更接近于实际情况,适用于复杂结构及速度变化强烈的地区模型;针对线性插值中存在的累计误差问题,刘玲君等(2014)提出了一种边界线性走时插值射线追踪算法,通过在不同区域的边界上插入节点来确定射线的折射角度,以减小累积误差,提高精度.近年来,不少学者对各向异性介质中的射线追踪方法进行了研究.Zhou和Greenhalgh(2005)研究了最短路径射线追踪法在二维和三维各向异性介质中的应用;赵后越和张美根(2014)利用混合网格各向异性地震波最短路径射线追踪算法,实现各个节点初至走时的计算,该方法计算速度快,适于复杂各向异性介质折射波正演模拟研究.

3.2 折射解释及相关技术

地震波初至拾取技术在浅层折射波勘探中起着重要作用,初至拾取的精度直接影响近地表地下结构的探测结果,尤其是在地形复杂地区,研究一种精度高且稳定性强的初至拾取算法意义重大.Peraldi和Clement(1972)Gelchinsky和Shtivelman(1983)在假设各道脉冲形状一致的前提下分别提出了相关法及拐点校正法,但实际记录的波形在各种因素作用下肯定会发生变化;李辉峰等(2006)提出了一种基于图像边缘检测技术的初至自动拾取方法,对地震记录灰度化处理后,利用Kirsch微分算子进行边缘检测;张伟等(2009)在利用时窗能量比求取初至的基础上,提出了变换时窗能量比法,但由于其仅是通过单道的能量比的最大值获取初至时间,而有些情况初至波处的能量比仅为次极大值,因此存在固有的缺陷;Sabbione和Velis(2010)研究了分别用于海洋地震记录、炸药地震记录和可控震源记录的初至自动拾取方法,通过检测信号在能量比、分形维数等特定属性上的变化识别初至,并利用边缘保持平滑算子和初至拾取校正技术提高初至拾取精度及一致性;冯智慧等(2011)利用地震信号横向波形的相似性及互四阶累积量一维切片对高斯色噪声不敏感的特点,提出了一种基于高阶累积量一维切片的初至拾取方法,但初至拾取的准确性受参考道初至起跳时间的标定、搜索范围、运算窗口长度等的影响;徐钰等(2012)提出了多时窗能量比法,同时考虑能量比最大值和次极值,根据多时窗间的能量比特征可以在一定程度上避免后续波能量变化和初至波前强能量短持续噪声的影响,有效识别初至波,同时在算法中加入奇异点检测和相位域起跳点估算,来提高该算法的稳定性,该方法在准噶尔盆地中部某区的地震勘探应用取得了很好的效果;Kahrizi和Hashemi(2014)对传统的人工神经网络算法进行改进,并将其用于折射波初至拾取中,研究表明基于改进神经网络算法的初至拾取方法具有较好的抗噪性能;陈金焕等(2015)提出了一种在参考初至的基础上,利用二分法和改进的能量比法检测初至并去除不准确点的初至自动拾取算法,该方法初至拾取快速准确且精度较高.

在实际折射波法地震勘探中,由于远偏移距处的折射波能量较弱,而近偏移距处的折射波场在直达波等强能量同相轴的干扰下难以分辨,故发展一种有效提取弱信号,提高折射波信噪比的方法尤为重要.Dong等(2006)首次将地震干涉法用于折射波信号的提取中,并引入虚折射干涉法的概念,但该方法不能确定虚折射波场的激发位置和时间,且经过相关处理后的地震记录的偏移距减小;Mallinson等(2011)Mikesell和van Wijk(2012)针对上述问题发展了超虚拟折射干涉法,但其难以恢复近偏移距处的折射波场,邵宇蓝等(2014)在此基础上,通过正演模拟数据和实际资料对该方法的抗噪性进行了研究;乔宝平等(20142015)提出了逆虚折射干涉法,通过原始折射波场和虚折射波场的互相关求得,并充分利用超级虚折射法和逆虚折射法的互补性,通过两者叠加有效获取弱信号,保证折射波在不同偏移距下振幅恢复的一致性,明显提高了信噪比,同时利用频率域反褶积消除子波旁瓣的影响,提高折射波信号的分辨率,该方法可有效实现折射波的提取,提高初至拾取的精度.Nichols等(2010)通过研究表明将地震干涉法用于实际资料时产生的虚假能量中隐含有地下地质信息,利用地震干涉法形成虚炮点记录压制不相关噪声,产生虚折射信号,联合利用虚炮点记录、相关道集和实际炮点记录可确定地下不饱和砂层、饱和砂层的速度及地下水位的相对深度,该方法可用于复杂地质结构近地表特征的解释.

折射波勘探发展至今,已形成多种解释方法(Hawkins,1961Palmer,19812001王振东,1988de Franco,2005赵德亨等,2005).由哈莱斯于1958年提出的哈莱斯法,基于相遇时距曲线通过作图求出临界折射时间,进而求得基岩面深度进行折射层构造解释,该解释方法对地表相对平坦而基岩面起伏较大或基岩内存在速度变化的地区,具有较好的应用效果,但该方法作图复杂,应用受到限制;之后对其改进产生了共轭点法,从数学推导出发,应用计算机逻辑演绎自动成图,利用共轭点法勾画出的界面轮廓更接近模型构造形态,适于复杂基岩面形态的解释;时间项法是一种通过计算炮点及检波点的延迟时间并将延迟时间的坐标转换为深度的解释方法,该方法简单实用,但受倾角影响严重,因此比较适用于地层界面起伏不大或倾角影响可以忽略的地区;t0法又称加减法,利用相遇观测曲线,根据已定义的计算公式,通过加减法得到t0曲线和差数时距曲线,并用于求取界面速度,绘制折射界面,该方法具有速度快、成本低的优点,但其不适用于横向速度变化较大的地质结构;波前法是由Thornburgh根据惠更斯原理提出的由地表到时重建地下波前的方法,根据Hagedoorn成像原理,当某一点位于折射界面上时,这两个波前系统在此点的到时之和等于互换时间,据此可通过绘制波在介质中的时间场获得折射界面,实现折射界面成像,此方法的不足之处是需要先确定覆盖层的速度,而且当存在隐伏层或明显的变速层时,波前法难以达到理论精度;截距时间法基于斯奈尔定理,利用时距曲线上测量到的地震波视速度和截距时间确定多层倾斜模型,但该方法仅适用于倾角较小和地层平缓的地质结构解释;互换法通过将截距法的计算扩展到每一个检波器位置上,来获得简单的折射层结构和速度变换情况,但不适用于大构造解释;广义互换法首先确定最佳距离,在此基础上利用时深函数等求得折射层速度及深度,该方法可以处理解释不规则折射层,适用范围较广;时间褶积剖面法利用正向地震道和反向地震道褶积,消除长时差和强振幅变化,提高信噪比,通过褶积剖面反映折射层结构特征;多折射层成像方法是利用正反向地震道的褶积与互换道互相关得到褶积互相关记录剖面,对经过处理的地震道叠加得到多折射层成像剖面,该方法可有效提高地震信号信噪比,实现多反射层解释;除上述图解法和解析法外,包括层析成像在内的各种折射波反演方法逐渐得到广泛应用,折射波反演主要是利用折射波走时、波形等信息来确定近地表速度和厚度分布,研究地层结构,该方法精度较高,适用于相对复杂的近地表介质.

当近地表存在一些较大的高速埋伏体时,由于埋伏体对信号较强的衰减、散射作用,使得浅层折射波法难以准确获取反映地下基岩的地震信息,此种现象在我国某些沿海近海岸的潮间带区域尤为常见,给工程地震勘察带来巨大困难.薛荣俊等(2000)通过分析讨论地震折射波在低速层中经过三种不同形状(扁平椭球体、球体、长方体)的高速埋伏岩体的射线路径及时距曲线关系,给出了直接利用折射波时距曲线确定高速埋伏体的厚度及埋深的计算方法,为定量估算工区基岩上覆低速层内高速埋伏体的数量、位置及规模提供了可能;近地表异常体调查研究对工程安全建设具有重要作用,Engelsfeld等(20082011)基于程函方程数值模拟研究了地下空洞对地震记录初至到时的影响,并推导出利用折射波初至到时求取圆形空洞位置和大小的计算公式,在此基础上,利用数值模拟进一步研究了两层模型中存在圆形空洞、圆形异常体、凹边界层、凸边界层及垂直裂隙带时的折射波时距曲线特征,并总结这五种常见的近地表异常体的时距曲线变化规律,用于指导复杂地质模型分析;Harmankaya等(2013)提出了一种利用地震干涉得到的虚散射波走时来估计近地表散射体和衍射位置的方法,通过共炮点记录上的散射波互相关得到虚散射波,并利用虚散射波到时反演确定散射体和衍射位置,用于探测近地表结构.

在地震勘探中,有时会存在最终处理获得的反射波剖面混乱、同相轴难以辨别,反射资料无法解释,而单炮记录中的折射波信息清晰丰富的情况,对此,Taner等(1992)Coruh(1993)发展了共深度折射面元叠加的折射波地震勘探方法,并取得了初步的应用效果.冯太林等(2001)在折射波理论及共深度折射面元叠加成像理论的基础上,提出了一种适用于折射波发育区的共深度点叠加成像方法,并给出了折射波速度分析和动校正的具体步骤,该方法在模型及实例中均可以较好地实现对折射界面的成像,但折射波速度点拾取的精度会影响到折射波成像的正确性.

浅层折射波法在近地表勘探中具有广泛应用,主要用于解决工程、环境和水文地质问题,以及地震勘探静校正中的一些参数问题.其可以用于寻找隐伏断层、确定近地表厚度和速度结构、分析折射层地质年代、确定基岩面起伏和潜水面深度、判别裂隙发育程度及评价岩体质量等,也可以用于追踪浅层界面及没有稳定反射的界面,补充反射波法的不足.周鸿秋等(1991)曾将初至折射波法用于极浅层、小裂隙勘探中,并利用密点距、追逐相遇双重观测系统,及波场特征分析、自动化定量解释方法,提高微小裂隙的分辨率和定量解释的精度,在地质病害调查中取得了理想的效果;熊章强和方根显(2001)采用折射波法进行隧道工程地质调查,得到了有效的覆盖层及各风化层厚度以及基岩界面的埋藏深度,根据波速划分围岩类别,为隧道工程的设计和开挖提供了地质依据;王洪(2007)详细讨论研究了折射波法在堤坝病害探查中应注意的问题,并从折射波数据采集、资料处理解释等方面给出了一整套适用于堤坝病害探查的具体工作方法,在实际应用中,利用浅层折射波准确解释了界面深度和浸润线深度,并指出该坝体存在软弱层面且部分位置填筑土密实度差,为后续的坝体维护及坝体补强位置的选择提供了重要依据;针对部分煤田存在的煤层埋藏浅、上覆基岩薄、积沙厚且极不均匀、沙层含水量大等地质问题,张湖钢等(2009)将折射地震勘探方法应用于煤田勘查中,以查明采区内潜水面、基岩面的埋藏深度和起伏形态以及基岩面风化层厚度,掌握古河床、古冲沟的分布;Carvalho等(2009)利用折射波勘探数据获得研究区P波和S波速度,并结合SPT(St and ard Penetration Test,标准贯入试验)数据得到近地表岩土参数,开展研究区场地效应和微震研究; Asry等(2012)将折射法用于基岩面调查中,利用折射地震速度剖面研究近地表地层分布情况,确定残留土层厚度和基岩面深度.

4 层析成像与多波勘探

地震层析成像自20世纪80年代初被引入到应用地球物理领域后,在地震勘探中逐渐得到广泛应用.地震初至层析成像是一种利用地震初至波走时重建地下介质速度结构的地球物理勘探方法,通过初至波层析成像反演可以建立近地表速度模型.李福中等(2000)将初至层析反演用于新疆某工区三维地震资料,精确重建了风化层和低降速带速度结构模型,利用该速度结构对地震剖面进行静校正处理,取得了较好的效果;针对三维地震初至走时成像中存在的数据量大、计算速度慢等问题,Noble等(2010)将伴随状态法引入到走时不拟合函数梯度的计算中,该方法可有效实现并行计算,解决大数据量计算问题;对埋藏较浅的煤层,地表起伏引起的静校正量对煤田勘探影响较大,周国婷等(2012)将初至层析成像用于反演近地表速度结构,利用获得的速度模型求取炮点检波点的静校正量,通过该方法获得的地震剖面上基本消除了起伏地表和低降速带对煤层反射波的影响,提高了地震剖面精度和分辨率;山地地表条件复杂,速度横向变化较大,近地表建模困难,陈云峰和杨淑卿(2014)利用层析成像方法反演该区近地表速度模型,较准确地获得了近地表低降速带的速度变化情况,为静校正处理提供了可靠的速度资料;崔岩和王彦飞(2015)建立了初至波走时层析成像的Tikhonov正则化反演模型,在此基础上提出了一种层析成像的梯度优化算法,数值试验表明,该方法适用于速度剧烈变化的任意复杂构造模型,且求解稳定易于实现,具有一定应用前景.在全波形反演技术发展的基础上,Sheng等(2006)提出了一种初至波形反演成像方法,该方法相比走时层析成像,考虑了更一般的地震波传播效应且具有更高的成像精度;Yu等(2014)通过研究表明,利用加权初至波形反演可以提高浅层速度分布估计的精度,该方法与传统的初至波形反演相比,能减小地震记录中噪声及振幅变化的影响,提高近地表成像精度.

反射层析成像是利用反射波确定地下速度结构和界面深度的重要方法,随着浅层高分辨率地震勘探的发展,地震反射层析成像方法由单纯的反演速度结构转向速度结构和界面深度联合反演.Bishop等(1985)Zelt等(1992)分别推导了反射波走时对界面深度的偏导数,为速度深度联合反演奠定了理论基础;华标龙和刘福田(1995)推导了二维介质情况下,地震反射波走时反演界面深度的关系式,同时利用正交算子投影法解决速度和界面同时成像的问题,数值模拟试算取得了较好的效果;黄光远(2000)指出,波在空间的传播过程是一个随着时间和空间变换的频谱变化过程,即实际记录的波形信号与时间、空间及频率有关,同时介质的物理参数(如反射系数、波阻抗等)也与空间、频率有关,据此他提出了一种基于地震反射波谱的反射波成像方法,给出了反射系数与反射波谱间的关系,同时在反演公式中还考虑了覆盖层滤波作用及多次反射波作用的影响;周龙泉等(2006)在前人研究的基础上,基于波逆行原理推导了三维介质中地震波射线走时对界面偏导数的基本公式,并给出任意多个复杂界面情况下,反演所需走时对界面偏导数的系数矩阵,通过数值模拟计算验证了该方法的有效性,但只是基于反射波的反演成像,没有考虑折射波,因而虽然对界面起伏进行了相对准确的定位,但对高速体的上下边界定位精度不高;走时成像基于宽射线角度覆盖范围来确定速度异常体的空间位置,当射线角度覆盖面较窄时,传统的网格层析成像方法会产生模糊假象,成像精度降低,Zhou(2006)提出了一种多尺度可变形地层成像方法,通过直接反演速度界面几何结构来提高成像精度;地裂缝是近地表一种特殊的地质灾害现象,给人们的生产、生活及安全带来不便与威胁,探测研究地裂缝的形成机理、产状、分布等具有重要意义,燕利芳(2011)将反射波层析成像用于西安地裂缝勘探,初步获得了近地表处地裂缝的位置和走向,为地裂缝成功探测提供了一种新思路.受挤压作用的构造带地形复杂,地表条件变化剧烈,此时折射波能量在复杂的近地表结构下穿透深度有限,单一的折射波走时成像难以实现有效的速度估计,Boonyasiriwat等(2009)将时间域多尺度波形成像用于挤压构造带速度估计,联合利用反射波和折射波能量实现近地表及深部速度结构估计;魏剑平(2012)进行了初至波走时和反射波走时联合反演研究,并通过理论模型和实际资料试算,表明联合反演比利用单一波反演精度高,稳定性好.

综合利用地震纵波和横波信息的多波勘探,相比于单纯的纵波勘探具有更高的勘探精度,目前浅层多波勘探技术正受到越来越多学者的关注.在松散层勘探中,横波比纵波探测精度高且不受饱水度的影响,陈相府和安西峰(2007)将横波勘探用于第四纪松散层层序划分及厚度探测中并取得了较好的效果;朱德兵(2008)通过对土坝体介质参数分析,并结合各种弹性波的特点,表明横波具有能量强,反射波成分简单,对垂向薄层分辨能力强、易于分辨及抗干扰能力强的特点,在土坝体隐患探测中具有较大优势;Krawczyk等(2012)利用横波勘探研究沉积层中的小尺度构造,区分坑洞类型,实现地下侵蚀区高分辨率成像.但在实际应用中横波具有激发难度大、工作效率低的缺点,目前多采用纵波震源激发,三分量检波器接收的方法,通过对多分量记录进行转换波的速度分析得到横波信息,实现多波勘探.转换点的精确计算是利用转换波处理提取横波信息中重要的一步,水平反射界面可以通过激发接收点位置、偏移距等确定,而对于倾斜反射界面则较为复杂,刘杰和刘江平(2013)在证明了唯一P-SV波走时对应唯一转换点的基础上,基于Snell定律和射线几何关系建立方程并得出转换点坐标的表达式,利用搜索法实现倾斜反射界面上最佳转换点的准确定位.速度分析同样是多波勘探中的重要环节,李录明和罗省贤(1995)通过对P-SV时距曲线分析,提出了一种利用转换波记录直接提取横波信息的方法;邓志文(2011)提出了从转换波记录中提取纵、横波速度的速度谱迭代法和纵横波速度联合扫描法;徐浩(2013)详细给出了转换波横波速度提取的方法,并将转换波速度分析应用到多分量记录中,通过纵横波速度谱的层位对比解释,进一步说明了利用多分量地震记录进行横波速度提取的有效性.转换波静校正是三分量地震勘探资料处理中的难点之一,刘洋和魏修成(2008)研究了转换型折射波传播规律,并据此提出了PPS波的识别方法和基于PPS波的转换波静校正方法,在实际三分量地震资料的转换波静校正中,取得了较好的效果.Pugin等(2009)对浅层多分量地震勘探进行了初步研究,发现常用的垂直激发震源同样可以产生丰富的SV和SH波反射,且SV波频率较高,可以产生高分辨率的反射剖面,并尝试通过极化分析和坐标旋转来提高反射剖面质量,研究表明三分量地震数据中含有地层各向异性及地质结构不同方向变化信息,有助于更好地研究地下构造和地层特征;Woelz等(2009)利用二维布设的三分量勘探方法研究位于某考古遗址的淤积古港口的基底形态,基于波场极性及传播方位角特征进行波场分离,并利用分离的波场实现三维折射成像;Giustiniani等(2010)联合利用P波和S波研究浅层地质结构特征,通过共中心点反射地震处理技术和多折射层成像技术实现研究区的地震成像;刘保金等(2012)应用浅层三分量地震反射法探测隐伏活动断裂,综合利用获得的纵波、横波和转化波剖面特征进行地层分层和断层解释,与跨断层的钻孔联合地质剖面结果相吻合,准确性高.

浅层地震记录中通常含有反射波、折射波及回折波等多种波列信息,而在常规地震勘探中,通常只利用一种波列,造成信息浪费.由于浅层折射波勘探简单有效,通过反演、层析成像可以获得各层厚度及纵波速度,但当地层速度倒转、速度差异小及存在薄层(易形成隐伏层)时,折射波法勘探就不适用了.而反射波法勘探是以波阻抗差异为物性前提的,分层能力强,反射剖面形态可反映地下地层形态,有效探测隐伏构造,但其对浅层界面分辨能力差.针对工程勘察中遇到的日益复杂的地质问题,多波联合应用逐渐得到重视.孙渊(1993)曾在陕北黄土梁峁地区隧道勘探中用地震反射波法和折射波法查明黄土覆盖层与基岩风化面之间的地质构造和速度分布情况;张学强等(1999)基于波场特征及波速关系进行了地震记录中反射纵波、折射纵波和瑞雷面波的层位对比研究,为多波联合反演奠定了基础;刘江平等(2000)将地震反射波法和瞬态瑞雷面波法相结合,用于大堤防(隔)渗墙质量检测中,达到检测墙的连续性、确定顶底板埋深及厚度的目的;张建清等(2001)进行了浅层多波勘探技术研究;由于垃圾填埋区介质不均匀性强,速度低且变化剧烈,激发的地震波信噪比低,单一的反射技术难以达到勘察目的,De Iaco等(2003)在钻孔资料约束下,联合利用反射波和折射波数据研究调查区概况,得到垃圾填埋场基底和周围沉积物结构;李红星(2005)从隐伏层的识别和层位对比、近地表覆盖层速度和深度参数反演、反射波和折射波联合勘探及反射波和折射波联合成像等方面对多波联合应用进行了全面研究;Orl and o和Pelliccioni(2010)利用折射波、面波和反射P波及PS波地震数据研究探区冲击谷构造的地质背景和P波、S波速度剖面,实际应用表明,联合应用多种数据资料可有效降低多解性,提高解释精度;Mendes等(2014)联合利用地震波场中的折射波和反射波来研究近地表特征,通过初至成像得到速度模型,并将其用于近地表反射系数估计,该方法可有效提高近地表结构信息的连续性,减少勘探成本;针对我国西部黄土地区特殊地貌,李庆春等(2014)开展了路基下伏地质缺陷地震探测技术研究,数值试验表明:反射波成像和面波成像技术可以实现对凹陷、断层、地裂缝、地下低速体、软弱夹层及滑坡等典型地质缺陷的清晰成像,多种方法相结合可提高勘探精度,实现对异常体的立体探测及准确定位.

5 近地表折射和反射法的发展趋势

针对近地表非均质极强、各向异性明显、地形起伏大、勘探深度小和精度要求高等特点,为了有效服务于国民建设,高分辨率、高信噪比、高保真度和物性探测是近地表折射和反射法的发展方向.

(1)多波多分量地震勘探方法在油气勘探领域已进行了深入研究,但在工程勘察、活动断层探测等浅层地震勘探中研究应用较少,由于浅层横波勘探分辨率高,对地层结构成像效果好但不易激发,而纵波勘探深度大且易于激发,但其成像效果差,因此综合利用纵、横波的优势,开展多波多分量联合勘探对提高浅层地震勘探的精度和分辨率具有重要作用;

(2)现在的浅层地震勘探主要是对地震剖面进行解释,容易忽略一些隐含的地质异常现象,属性提取技术是充分提取地震信息,进行全面综合解释的有效手段,针对浅层工程勘探的特点合理选择属性参数或研究引入新属性参数是提高浅层地震勘探精度的重要手段;

(3)开展多层折射介质的观测系统和解释方法研究,尤其是折射层析成像研究,是提高多层折射介质成像精度的途径;

(4)由于近地表非均质极强、各向异性明显和裂隙孔隙发育,开展黏弹性、双相和各向异性介质的地震反射与折射波法研究尤为重要.

参考文献
[1] Almalki H, Alkhalifah T. 2012. Mapping the base of sand dunes using a new design of land-streamer for static correction applications. Journal of Petroleum Exploration and Production Technology, 2(2): 57-65.
[2] Asry Z, Samsudin A R, Yaacob W Z, et al. 2012. Geoelectrical resistivity imaging and refraction seismic investigations at Sg. Udang, Melaka. American Journal of Engineering and Applied Sciences,5(1): 93-97.
[3] Barison E, Brancatelli G, Nicolich R, et al. 2011. Wave equation datuming applied to marine OBS data and to land high resolution seismic profiling. Journal of Applied Geophysics, 73(3): 267-277.
[4] Belina F A, Ernst J R, Holliger K. 2009. Inversion of crosshole seismic data in heterogeneous environments: Comparison of waveform and ray-based approaches. Journal of Applied Geophysics, 68(1): 85-94.
[5] Bishop T N, Bube K P, Cutler R T, et al. 1985. Tomographic determination of velocity and depth in laterally varying media. Geophysics, 50(6): 903-923.
[6] Boonyasiriwat C, Valasek P, Zhu X H. 2009. Application of multiscale waveform tomography for high-resolution velocity estimation in complex geologic environments: Canadian Foothills synthetic data example. The Leading Edge, 28(4): 454-456.
[7] Cao J T, Qu W Q, Fang Y. 2014. Application of shallow seismic reflection wave method in engineering investigation of mountain area. Geotechnical Investigation & Surveying (in Chinese), (8): 84-88.
[8] Carpenter E W. 1996. Absorption of elastic wave sanoperatorfora constant, Qmechanism. AWRE Report, London.
[9] Carvalho J, Torres L, Castro R, et al. 2009. Seismic velocities and geotechnical data applied to the soil microzoning of western Algarve, Portugal. Journal of Applied Geophysics, 68(2): 249-258.
[10] Chang X, Liu Y K, Gui Z X. 2006. Zero-offset reverse time migration for prediction ahead of tunnel face. Chinese Journal of Geophysics (in Chinese), 49(5): 1482-1488.
[11] Chen J H, Cao Y S, Sun C L, et al. 2015. The algorithm for automatic first-breaks picking on seismic traces based on Dichotomy. Progress in Geophysics (in Chinese), 30(2): 688-694, doi: 10.6038/pg20150228.
[12] Chen X F, An X F, Wang G W. 2005. Application of high-resolution seismic detecting for undermined shallow strata. Progress in Geophysics (in Chinese), 20(2): 437-439.
[13] Chen X F, An X F, Wang G W. 2005. Application of high-resolution seismic detecting for undermined shallow strata. Progress in Geophysics (in Chinese), 20(2): 437-439.
[14] Chen X F, An X F. 2007. Application of seismic shear wave detection in dividing shallow strata of rock and soil. Progress in Geophysics (in Chinese), 22(5): 1655-1659.
[15] Chen Y F, Yang S Q. 2014. Building of near surface velocity model using tomography inversion method. Marine Geology Frontiers (in Chinese), 30(5): 47-51.
[16] Cheng F, Liu J P, Qu N N, et al. 2014. Two-dimensional pre-stack reverse time imaging based on tunnel space. Journal of Applied Geophysics, 104: 106-113.
[17] Coruh C. 1993. Composite refraction-reflection stack sections: tracing faults in the Atlantic Coastal plain sediments. // 63rd Ann. Internat. Mtg. Soc. Expl. Geophys., Expanded Abstracts, 1157-1160.
[18] Cui Y, Wang Y F. 2015. Tikhonov regularization and gradient descent algorithms for tomography using first-arrival seismic traveltimes. Chinese Journal of Geophysics (in Chinese), 58(4): 1367-1377, doi: 10.6038/cjg20150423.
[19] Dai Y, Zhang J Z, Zhong B S, et al. 1999. A refraction forward modelling algorithm by using Fermat theory. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 21(4): 323-326.
[20] De Iaco R, Green A G, Maurer H R, et al. 2003. A combined seismic reflection and refraction study of a landfill and its host sediments. Journal of Applied Geophysics, 52(4): 139-156.
[21] de Franco R. 2005. Multi-refractor imaging with stacked refraction convolution section. Geophysical Prospecting, 53(3): 335-348.
[22] Deng Z W. 2011. Approaches to mode-converted wave velocity analysis on 3C data. Oil Geophysical Prospecting (in Chinese), 46(4): 570-575.
[23] Ding J J, Dai Y S, Zhang Y N, et al. 2013. Improving seismic data resolution based on high frequency compensation. Progress in Geophysics (in Chinese), 28(6): 3214-3221, doi: 10.6038/pg20130648.
[24] Dong S Q, Sheng J M, Schuster J T. 2006. Theory and practice of refraction interferometry. // 76th SEG Annual Meeting, Expanded Abstracts, 3021-3025.
[25] Du X Z, Cao J X. 2008. Shallow subsurface seismic prospecting methods in high-noise condition. Progress in Geophysics (in Chinese), 23(5): 1616-1620.
[26] Engelsfeld T, umanovac F, Pavin N. 2008. Investigation of underground cavities in a two-layer model using the refraction seismic method. Near Surface Geophysics, 6(4): 221-231.
[27] Engelsfeld T, umanovac F, Krstic′ V. 2011. Classification of near-surface anomalies in the seismic refraction method according to the shape of the time-distance graph: A theoretical approach. Journal of Applied Geophysics, 74(1): 59-68.
[28] Feng T L, Zhang X G, Li Y D, et al. 2001. Research on methodology of stack imaging of refractive seismic recording. Chinese Journal of Geophysics (in Chinese), 44(1): 129-134.
[29] Feng Z H, Liu C, Feng X, et al. 2011. A method of automatic first arrival extraction based on one-dimensional slice of cross-fourth-order cumulants. Journal of Jilin University (Earth Science Edition) (in Chinese), 41(2): 559-564, 571.
[30] Fleury C, Snieder R, Larner K. 2010. General representation theorem for perturbed media and application to Green′s function retrieval for scattering problems. Geophysical Journal International, 183(3): 1648-1662.
[31] Fradelizio G L, Levander A, Zelt C A. 2008. Three-dimensional seismic-reflection imaging of a shallow buried paleochannel. Geophysics, 73(5): B85-B98.
[32] Fu X M. 2012. Seismic attributes and its application in engineering investigations [Master′s thesis] (in Chinese). Chengdu: Chengdu University of Technology.
[33] Górszczyk A, Adamczyk A, Malinowski M. 2014. Application of curvelet denoising to 2D and 3D seismic data-Practical considerations. Journal of Applied Geophysics, 105: 78-94.
[34] Gao X, Zhang Y W, Wei W B. 2007. Seismic geophysical exploration in Pingtan Strait Bridge, Fujian, East China Sea. Earth Science-Journal of China University of Geosciences (in Chinese), 32(4): 540-544.
[35] Ge M C, Hardy H R Jr, Wang H L, et al. 2011. Development of a simple mechanical impact system as a non-explosive seismic source. Geotechnical and Geological Engineering, 29(1): 137-142.
[36] Gelchinsky B, Shtivelman V. 1983. Automatic picking of first arrivals and parameterization of traveltime curves. Geophysical Prospecting, 31(6): 915-928.
[37] Giustiniani M, Tinivella U, Accaino F. 2010. P and S reflection and P refraction: An integration for characterising shallow subsurface. Journal of Applied Geophysics, 71(4): 149-156.
[38] Guo S J, Hou Y H, Li H L, et al. 2014. The application of seismic CT technology in urban geological survey. Computerized Tomography Theory and Applications (in Chinese), 23(1): 81-90.
[39] Haffinger P, Verschuur E. 2012. Estimation and application of near-surface full waveform redatuming operators. Geophysical Prospecting, 60(2): 270-280.
[40] Hargreaves N D, Calvert A J. 1991. Inverse Q-filtering by Fourier transform. Geophysics, 56(4): 519-527.
[41] Harmankaya U, Kaslilar A, Thorbecke J, et al. 2013. Locating near-surface scatterers using non-physical scattered waves resulting from seismic interferometry. Journal of Applied Geophysics, 91: 66-81.
[42] Hawkins L V. 1961. The reciprocal method of routine shallow seismic refraction investigations. Geophysics, 26(6): 808-819.
[43] He G M, He Z H, Huang D J, et al. 2006. The comparison of several static′s correction methods. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 28(4): 310-314.
[44] Heilmann Z, Deidda G P, Satta G, et al. 2013. Real-time imaging and data analysis for shallow seismic data using a cloud-computing portal. Near Surface Geophysics, 11(4): 407-421.
[45] Hu T, White R E. 1998. Robust multiple suppression using adaptive beamforming. Geophysical Prospecting, 46(3): 227-248.
[46] Hua B L, Liu F T. 1995. Joint reflection imaging of velocity and interface—theory and method. Chinese Journal of Geophysics (Acta Geophysica Sinica) (in Chinese), 38(6): 750-756.
[47] Huang G Y. 2001. New imaging method for seismic reflection wave and its theoretical basis. Science in China Series D: Earth Sciences (in Chinese), 44(2): 159-165.
[48] Imhof A L, Sánchez M, Calvo C, et al. 2011. Application of seismic refraction tomography for tunnel design in Santa Clara Mountain, San Juan, Argentina. Earth Sciences Research Journal, 15(2): 81-88.
[49] Jeng Y, Li Y W, Chen C S, et al. 2009. Adaptive filtering of random noise in near-surface seismic and ground-penetrating radar data. Journal of Applied Geophysics, 68(1): 36-46.
[50] Jiang T, Lin J, Chen Z B, et al. 2004. Seismic beam-forming in phased array of vibroseis. Journal of Jilin University (Information Science Edition) (in Chinese), 22(3): 181-184.
[51] Jiang T, Lin J, Yang D, et al. 2008. Analysis of directional seismic signal based on phased-array vibrator system. Chinese Journal of Geophysics (in Chinese), 51(5): 1551-1556.
[52] Jiang T, Lin J, Jia H Q, et al. 2012. Time-domain seismic beam-forming based on receiver arrays. Chinese Journal of Geophysics (in Chinese), 55(12): 4277-4287, doi: 10.6038/j.issn.0001-5733.2012.12.039 .
[53] Kahrizi A, Hashemi H. 2014. Neuron curve as a tool for performance evaluation of MLP and RBF architecture in first break picking of seismic data. Journal of Applied Geophysics, 108: 159-166.
[54] Krawczyk C M, Polom U, Trabs S. 2012. Sinkholes in the city of Hamburg-New urban shear-wave reflection seismic system enables high-resolution imaging of subrosion structures. Journal of Applied Geophysics, 78: 133-143.
[55] Li F Z, Xing G D, Bai X M, et al. 2000. Study of first break Tomographic inversion method for static corrections. Oil Geophysical Prospecting (in Chinese), 35(6): 710-718.
[56] Li H F, Zou Q, Jin W L. 2006. Method of automatic first breaks pick-up based on edge detection. Oil Geophysical Prospecting (in Chinese), 41(2): 150-155, 159.
[57] Li H X. 2005. Studies in the united application and simultaneous imaging of multi-wave [Master′s thesis] (in Chinese). Jilin: Jilin University.
[58] Li L M, Luo X X. 1995. Velocity analysis and interpretation of P-SV wave. Oil Geophysical Prospecting (in Chinese), 30(1): 66-74.
[59] Li Q C, Wu H, Zhou X M. 2014. Numerical modeling of elastic waves for geological defects under roadbed. Journal of Jilin University (Earth Science Edition) (in Chinese), 44(6): 2016-2025.
[60] Li W W, Pei Z L. 2008. Numerical simulation of elastic waves between wells. Geophysical and Geochemical Exploration (in Chinese), 32(2): 207-211.
[61] Li X P. 2009. Application of estuary high density seismic imaging to the island-mainland link project in Zhoushan. Journal of Geotechnical Investigation & Surveying (in Chinese), (11): 91-95.
[62] Li Z S, Yan Z H, Lin R. 1998. Seismic image method—one of new water engineering geophysical prospecting. Geotechnical Investigation and Surveying (in Chinese), (3): 70-73.
[63] Li Z S. 2007. Seismic wave attenuation and recovering method. Progress in Geophysics (in Chinese), 22(5): 1545-1551.
[64] Li Z W. 2013. The study of signal acquisition and imaging method on tunnel geological prediction [Master′s thesis] (in Chinese). Wuhan: Yangtze University.
[65] Liu B J, Zhao C B, Feng S Y, et al. 2012. Application of the three-component shallow seismic reflection method to probing buried active faults. Chinese Journal of Geophysics (in Chinese), 55(8): 2676-2686, doi: 10.6038/j.issn.0001-5733.2012.08.020.
[66] Liu F L, Xu F, Liu W G. 2014. Unequal-deep and time-delay source array. Oil Geophysical Prospecting (in Chinese), 49(2): 225-229.
[67] Liu J, Liu J P. 2013. A solution for calculating conversion of dip reflector based on the P-SV wave traveltime. Geophysical Prospecting for Petroleum (in Chinese), 52(6): 602-608, 616.
[68] Liu J P, Qian S H, Zhang S X. 1998. The stacking theoretic model calculation and analysis for the static and NMO correction of the seismic record in CMP gather at the slip reference plane. Geophysical & Geochemical Exploration (in Chinese), 22(1): 27-33.
[69] Liu J P, Yang Y Q, Li X Z. 1999. The method for extracting velocity under the complex topographic condition and its analysis. Geophysical & Geochemical Exploration (in Chinese), 23(4): 258-264.
[70] Liu J P, Zhang L Q, Zhang Y M, et al. 2000. The application of shallow seismic technology to the quality inspection of the leakage-proof embankment wall. Geophysical & Geochemical Exploration (in Chinese), 24(4): 302-306.
[71] Liu J P, Xia J H, Chen C, et al. 2005. Accurate elevation and normal moveout corrections of seismic reflection data on rugged topography. New Zealand Journal of Geology and Geophysics, 48(4): 707-716.
[72] Liu J P, Cheng F, Fan C Y, et al. 2012. Two-dimensional numerical simulation of tunnel-based seismic full-wave fields. Chinese Journal of Geotechnical Engineering (in Chinese), 34(9): 1705-1711.
[73] Liu L J, Xie Z H, Yang C. 2014. Ray tracing algorithm based on boundary linear traveltime interpolation. Journal of South China University of Technology (Natural Science Edition) (in Chinese), 42(5): 23-28, 35.
[74] Liu S W, Geng J H, Feng W. 2005. Controlled illumination and seismic acquisition geometry for target-oriented imaging. Applied Geophysics, 12(4): 230-234.
[75] Liu Y, Wei X C. 2008. Converted refraction wave propagation law and its application in converted wave static correction. Science in China (Series D: Earth Science) (in Chinese), 38(S1): 204-210.
[76] Liu Y H, Tian Z Y, Yu Z H. 2005. Combined prospecting method of engineering VSP and seismic CT and its application on geotechnology. Progress in Geophysics (in Chinese), 20(1): 267-272.
[77] Liu Z M, Lei W, Lu H Y, et al. 2006. The application of shallow seismic reflection in preselecting the location of a major engineering project. Journal of Chengdu University of Technology (Science & Technology Edition) (in Chinese), 33(4): 345-348.
[78] Lu G Y, Xiong Y, Zhu Z Q. 2011. Detection simulation ahead of tunnel face and reverse-time migration with reflection wave method. Journal of Central South University (Science and Technology) (in Chinese), 42(1): 136-141.
[79] Lu W K, Ding W L, Zhang S W, et al. 2005. A high resolution processing technique for 3_D seismic data based on signal subspace decomposition. Chinese J. Geophys. (in Chinese), 48(4): 896-901.
[80] Mallinson I, Bharadwaj P, Schuster G, et al. 2011. Enhanced refractor imaging by supervirtual interferometry. The Leading Edge, 30(5): 546-550.
[81] Masjukov A V, Murzov I S. 2013. Short-period static correction without normal move out. Geophysical Prospecting, 61(1): 92-103.
[82] Mendes M, Mari J L, Hayet M. 2014. Imaging geological structures up to the acquisition surface using a hybrid refraction-reflection seismic method. Oil & Gas Science and Technology-Rev. IFP Energies nouvelles, 69(2): 351-361.
[83] Mikesell D, van Wijk K. 2012. Seismic refraction interferometry with a semblance analysis on the crosscorrelation gather. Geophysics, 76(5): SA77-SA82.
[84] Nichols J, Mikesell D, van Wijk K. 2010. Application of the virtual refraction to near-surface characterization at the Boise Hydrogeophysical Research Site. Geophysical Prospecting, 58(6): 1011-1022.
[85] Noble M, Thierry P, Taillandier C, et al. 2010. High-performance 3D first-arrival traveltime tomography. The Leading Edge, 29(1): 86-93.
[86] Orlando L, Pelliccioni G. 2010. P and PS data to reduce the uncertainty in the reconstruction of near-surface alluvial deposits (Case study-Central Italy). Journal of Applied Geophysics, 72(1): 57-69.
[87] Palmer D. 1981. An introduction to the generalized reciprocal method of seismic refraction interpretation. Geophysics, 46(11): 1508-1518.
[88] Palmer D. 2001. Imaging refractors with the convolution section. Geophysics, 66(5): 1582-1589.
[89] Panea I, Drijkoningen G. 2008. The spatial data-adaptive minimum-variance distortionless-response beamformer on seismic single-sensor data. Geophysics, 73(5): Q29-Q42.
[90] Papandreou B, Brennan M J, Rustighi E. 2011. On the detection of objects buried at a shallow depth using seismic wave reflections. The Journal of the Acoustical Society of America, 129(3): 1366-1374.
[91] Peraldi R, Clement A. 1972. Digital processing of refraction data study of first arrivals. Geophysical Prospecting, 20(3): 529-548.
[92] Pham M Q, Duval L, Chaux C, et al. 2014. A primal-dual proximal algorithm for sparse template-based adaptive filtering: application to seismic multiple removal. IEEE Transactions on Signal Processing, 62(16): 4256-4269.
[93] Pugin A J M, Pullan S E, Hunter J A. 2009. Multicomponent high-resolution seismic reflection profiling. The Leading Edge, 28(10): 1248-1261.
[94] Qiao B P, Guo P, Wang P, et al. 2014. Effectively picking weak seismic signal near the surface based on reverse virtual refraction interferometry. Chinese Journal of Geophysics (in Chinese), 57(6): 1900-1909, doi: 10.6038/cjg20140621.
[95] Qiao B P, Guo P, Wang P, et al. 2015. Retrieval of super-virtual refraction by cross-correlation. Geophysical Prospecting, 63(3): 552-566.
[96] Rashed M A. 2009. GDP: A new source for shallow high-resolution seismic exploration. Journal of Applied Geophysics, 68(2): 243-248.
[97] Sabbione J I, Velis D. 2010. Automatic first-breaks picking: New strategies and algorithms. Geophysics, 75(4): V67-V76.
[98] Santana F R, Neto A A. 2015. Use of seismic attributes for sediment classification. The Journal of Engineering, 1-15.
[99] Shao Y L, Li T L, Liao Y, et al. 2014. Refraction wave SNR improvement based on SVI method. Marine Geology Frontiers (in Chinese), 30(5): 52-57, 69.
[100] Shen Y, Zhang T, Zhao M, et al. 2011. Discussion on complex areas of shallow seismic acquisition methods. Coal Technology (in Chinese), 30(6): 253.
[101] Sheng J M, Leeds A, Buddensiek M, et al. 2006. Early arrival waveform tomography on near-surface refraction data. Geophysics, 71(4): U47-U57.
[102] Shi M, Feng D S, Dai Q W. 2006. Application of integrated geophysical method in quality detection of dikes. Progress in Geophysics (in Chinese), 21(4): 1328-1331.
[103] Shi Z J, Tian G, Dong S X, et al. 2004. Attenuation compensation of low-velocity layers using micro-seismogram logs: case studies. Journal of Geophysics and Engineering, 1(3): 181-186.
[104] Sun Y. 1993. An application of the shallow synthetical seismic prospecting methods to the tunnel prospecting. Journal of Xi′an College of Geology (in Chinese), 15(2): 77-81.
[105] Sun Y M, Verschuur E, Vrolijk J W. 2014. Solving the complex near-surface problem using 3D data-driven near-surface layer replacement. Geophysical Prospecting, 62(3): 491-506.
[106] Sun Y M, Verschuur D J. 2014. A self-adjustable input genetic algorithm for the near-surface problem in geophysics. IEEE Transactions on Evolutionary Computation, 18(3): 309-325.
[107] Tan C Y, Yin C, Lin Y H, et al. 2009. Static corrections by integrative optimization in complex area. Oil Geophysical Prospecting (in Chinese), 44(3): 288-291.
[108] Tan S J, Zhang H, Wang R X. 2013. The cross-well seismic computerized tomography technology and it′s application in the cavern survey. Computerized Tomography Theory and Applications (in Chinese), 22(3): 439-446.
[109] Taner M T, Matsuoka T, Baysal E. 1992. Imaging with refractive seismic waves. // 62nd Ann. Internat, Mtg. Soc. Expl. Geophys., Expanded Abstracts, 1132-1135.
[110] Traonmilin Y, Gulunay N. 2013. Statics-preserving projection filtering. Geophysical Prospecting, 61(Suppl.1): 346-351.
[111] Tu Q C, Liu H S. 2006. The application of linear traveltime interpolation ray tracing to the calculation of first-break traveltime for the near-surface model. Geophysical and Geochemical Exploration (in Chinese), 30(2): 148-153.
[112] Tzavaras J, Buske S, Gross K. 2012. Three-dimensional seismic imaging of tunnels. International Journal of Rock Mechanics & Mining Sciences, 49: 12-20.
[113] Wang C L, Liu Z P. 2012. The modelling study of τ-p transform applied in the tunnel seismic prediction. Progress in Geophysics (in Chinese), 27(5): 2216-2225, doi: 10.6038/j.issn.1004-2903.2012.05.049.
[114] Wang C L, Liu Z P, Huang Y Y, et al. 2014. Modeling research of F-K transform applied in the tunnel earthquake prediction. Journal of Geodesy and Geodynamics (in Chinese), 34(4): 71-76, 82.
[115] Wang H. 2007. The application research of seismic refraction method in dike disease detection [Master′s thesis] (in Chinese). Changsha: Central South University.
[116] Wang J J, Liao Q T, Cao J W, et al. 2006. The application of the cross hole seismic computerized topography to the detection of fissures in bridge foundation. Geophysical and Geochemical Exploration (in Chinese), 30(2): 186-188.
[117] Wang J M, Chen S M, Su M X, et al. 2007. A study of the near surface high-frequency compensation technology in 3-D seismic exploration. Chinese Journal of Geophysics (in Chinese), 50(6): 1837-1843.
[118] Wang R, Zhong S H, Wang Z F. 2013. Application of Landsonar method on fault and karst exploration in busy town of Jinan. Progress in Geophysics (in Chinese), 28(6): 3313-3316, doi: 10.6038/pg20130661.
[119] Wang S D. 2011. Attenuation compensation method based on inversion. Applied Geophysics, 8(2): 150-157.
[120] Wang Y A, Liu Z J. 2001. Seismic geophysical exploration in Yangluo Yangtze bridge, Wuhan. Earth Science (in Chinese), 26(4): 388-390.
[121] Wang Z D. 1988. Shallow Seismic Exploration Technique (in Chinese). Beijing: Geological Publishing House.
[122] Wang Z J, He Q D, Tian G, et al. 2005. A study of engineering reflect seismic with unconventional geometry layout. Journal of Jilin University (Earth Science Edition) (in Chinese), 35(6): 796-800.
[123] Wei J P. 2012. High-precision joint velocity inversion with seismic first-break and reflection travel-time [Master′s thesis] (in Chinese). Xi′an: Chang′an University.
[124] Wei Z H, Hall M A, Phillips T F. 2012. Geophysical benefits from an improved seismic vibrator. Geophysical Prospecting, 60(3): 466-479.
[125] Woelz S, Rabbel W, Mueller C. 2009. Shear waves in near surface 3D media-SH-wavefield separation, refraction time migration and tomography. Journal of Applied Geophysics, 68(1): 104-116.
[126] Wu X G, Li Y L, Zhang M, et al. 2012a. Causes of and solutions to a low signal-to-noise ratio of seismic data in complex areas, Part Ⅰ: Scattering of surface wave resulting in low signal-to-noise ratio of seismic data. Natural Gas Industry (in Chinese), 32(1): 27-32.
[127] Wu X G, Li Z R, Li Z J, et al. 2012b. Causes and solutions of a low signal-to-noise ratio of seismic data in complex areas, Part Ⅱ: Basic property and shape features of surface scattered waves. Natural Gas Industry (in Chinese), 32(2): 38-42.
[128] Wu X G, Li Z R, Zhang X B, et al. 2012c. Causes of and solutions to a low signal-to-noise ratio of seismic data in complex areas, Part Ⅲ: Technical methods of surface scattered wave suppression. Natural Gas Industry (in Chinese), 32(3): 57-61.
[129] Wu Z Q, Liu Y S, Liu B J, et al. 2005. Application of geophysical method to the precise positioning of urban seismic mobile moving fault. Progress in Geophysics (in Chinese), 20(2): 528-533.
[130] Xiao Q H, Xie C J. 2012. Application of tunnel seismic tomography to tunnel prediction in karst area. Rock and Soil Mechanic (in Chinese), 33(5): 1416-1420.
[131] Xie B, Wang W, Wu C Y. 2007. Application of shallow earthquake reflection shooting in sludge thickness surveying in water intake of hydropower plant. Soil Engineering and Foundation (in Chinese), 21(2): 76-78.
[132] Xie Z Q, Wen P L. 2000. Shallow seismic reflection survey technique to complicated object. Journal of Central South University of Technology (in Chinese), 31(6): 477-480.
[133] Xiong Z Q, Fang G X. 2001. Application of refractive method on tunnel-engineering exploration. Hydrogeology and Engineering Geology (in Chinese), (1): 63-66.
[134] Xu F, Liu F L, Liang X H. 2011. A method for source array design based on the phased-array theory. Oil Geophysical Prospecting (in Chinese), 46(2): 170-175.
[135] Xu H. 2013. S-wave velocity extraction and case studies based on multi-component seismic records [Master′s thesis] (in Chinese). Wuhan: China University of Geosciences.
[136] Xu J X. 2014. Separating the near-surface seismic scattered wave using seismic interferometry method. Chinese Journal of Geophysics (in Chinese), 57(6): 1910-1923, doi: 10.6038/cjg20140622.
[137] Xu T, Xu G M, Gao E G, et al. 2004. Refraction ray-tracing in complex medium. Oil Geophysical Prospecting (in Chinese), 39(6): 690-693.
[138] Xu Y, Zeng W H, Song J G, et al. 2012. A new method for first-break picking in refraction exploration. Oil Geophysical Prospecting (in Chinese), 47(2): 218-224.
[139] Xue R J, Li J S, Tang Y. 2000. The calculation for the high speed mass of rock overlay by low speed medium in the tidal zone of the shallow sea. Journal of Ocean University of Qingdao (in Chinese), 30(4): 685-689.
[140] Yan L F. 2011. Reflected wave tomographic simulation and the application in Xi′an ground fissures [Master′s thesis] (in Chinese). Xi′an: Chang′an University.
[141] Yu H, Zhang D L, Wang X. 2014. Application of weighted early-arrival waveform inversion to shallow land data. Journal of Applied Geophysics, 102: 68-76.
[142] Zelt C A, Smith R B. 1992. Seismic traveltime inversion for 2-D crustal velocity structure. Geophys. J. Int., 108(1): 16-34.
[143] Zhang G L, Xiong X J, Rong J J, et al. 2010. Stratum absorption and attenuation compensation based on improved generalized S-transform. Oil Geophysical Prospecting (in Chinese), 45(4): 512-515.
[144] Zhang H G, He Z Y, Shan L, et al. 2009. Application and effect of shallow seismic refraction method in Jurassic coalfield loose beds, Northern Shaanxi. Coal Geology of China (in Chinese), 21(9): 55-57, 69.
[145] Zhang J Q, Zhang Z, Xia W L. 2001. Research on the new exploration technique of multiwave for shallow earthquake. Geotechnical Investigation & Surveying (in Chinese), (6): 65-68.
[146] Zhang K, Wang Q S. 2010. The shallow seismic reflection survey in the prediction of the collapse. Geophysical and Geochemical Exploration (in Chinese), 34(1): 85-88.
[147] Zhang P, Liu H, Li Y M. 2000. The situation and progress of ray tracing method research. Progress in Geophysics (in Chinese), 15(1): 36-45.
[148] Zhang W, Wang H, Li H C, et al. 2009. The application of statistic energy ratio of transform time windows to the pickup of seismic first arrival wave. Geophysical and Geochemical Exploration (in Chinese), 33(2): 178-180.
[149] Zhang X Q, Bao X Q, Liu J P, et al. 1999. Classification of seismic profiles and horizontal correlation of three shallow waves. Geophysical and Geochemical Exploration (in Chinese), 23(1): 7-13.
[150] Zhang X Q, Zhu X P, Yan Z, et al. 2006. Multiple wave suppression method for water seismic image source. // The 2nd International Conference on Environmental and Engineering Geophysics (in Chinese), 4-6.
[151] Zhao C L, Lu C, Hao T Y, et al. 2013. A study of the high-precision modular lightweight small vibrator. Chinese Journal of Geophysics (in Chinese), 56(11): 3690-3698, doi: 10.6038/cjg20131110.
[152] Zhao D H, Tian G, Wang B B. 2005. An overview of shallow seismic refraction methods. Global Geology (in Chinese), 24(2): 188-193.
[153] Zhao H Y, Zhang M G. 2014. Tracing seismic shortest path rays in anisotropic medium with rolling surface. Chinese Journal of Geophysics (in Chinese), 57(9): 2910-2917, doi: 10.6038/cjg20140916.
[154] Zhao Y G. 2007. Analysis and recommendation of tunnel prediction techniques at home and abroad. Progress in Geophysics (in Chinese), 22(4): 1344-1352.
[155] Zhao Y G, Jiang H. 2009. Theoretical problems and judgment of tunnel seismic prediction technology. //The 25th Annual Meeting of the Chinese Geophysical Society (in Chinese). Beijing: Science Press, 303-304.
[156] Zhong S H. 1997. Land sonar method and effectiveness. Geophysical & Geochemical Exploration (in Chinese), 21(3): 172-179.
[157] Zhong S H. 2003. Analysis on the failure reasons in survey forward from working face of tunnels by means of TSP. Chinese Journal of Rock Mechanics and Engineering (in Chinese), 22(1): 2443-2446.
[158] Zhou B, Greenhalgh S A. 2005. ‘Shortest path’ ray tracing for most general 2D/3D anisotropic media. Journal of Geophysics and Engineering, 2(1): 54-63.
[159] Zhou G T, Pan D M, Niu H, et al. 2012. Research and application of the tomographic static correction method. Geophysical and Geochemical Exploration (in Chinese), 36(5): 802-805, 841.
[160] Zhou H Q, Gu H M, Zhang X Q, et al. 1991. The application of optimization method of shallow layer refraction system to investigation of geological hazards. Earth Science—Journal of China University of Geosciences (in Chinese), 16(4): 445-452.
[161] Zhou H W. 2006. Multiscale deformable-layer tomography. Geophysics, 71(3): R11-R19.
[162] Zhou L Q, Liu F T, Chen X F. 2006. Simultaneous tomography of 3-D velocity structure and interface. Chinese Journal of Geophysics (in Chinese), 49(4): 1062-1067.
[163] Zhou Z S, Ni J X, Feng Y, et al. 2010. Research on the method of fixed unit array for shallow seismic reflection. Geotechnical Investigation & Surveying (in Chinese), (7): 84-88.
[164] Zhu D B. 2008. The advantage of S-wave in detecting hidden troubles in earth dam. Progress in Geophysics (in Chinese), 23(2): 612-616.
[165] Zhu J P, Dong L G, Cheng J B. 2011. Target-oriented 3D seismic optimal geometry design based on seismic illumination. Oil Geophysical Prospecting (in Chinese), 46(3): 339-348.
[166] Zuo J J, Lin S H, Kong Q F, et al. 2011. Imaging by a joint tomography of direct wave traveltime and reflection wave traveltime in crosswell seismic and its application. Oil Geophysical Prospecting (in Chinese), 46(2): 226-231.
[167] 曹江涛, 翟伟强, 方勇. 2014. 浅层地震反射波法在山区工程勘察中的应用. 工程勘察, (8): 84-88.
[168] 常旭, 刘伊克, 桂志先. 2006. 反射地震零偏移距逆时偏移方法用于隧道超前预报. 地球物理学报, 49(5): 1482-1488.
[169] 陈金焕, 曹永生, 孙成龙等. 2015. 基于二分法的地震波初至自动拾取算法. 地球物理学进展, 30(2): 688-694, doi: 10.6038/pg20150228.
[170] 陈相府, 安西峰, 王高伟. 2005. 浅层高分辨地震勘探在采空区勘测中的应用. 地球物理学进展, 20(2): 437-439.
[171] 陈相府, 安西峰. 2007. 地震横波勘探及其在浅层岩土分层中的应用. 地球物理学进展, 22(5): 1655-1659.
[172] 陈云峰, 杨淑卿. 2014. 利用层析成像方法反演近地表速度模型. 海洋地质前沿, 30(5): 47-51.
[173] 崔岩, 王彦飞. 2015. 基于初至波走时层析成像的Tikhonov正则化与梯度优化算法. 地球物理学报, 58(4): 1367-1377, doi: 10.6038/cjg20150423.
[174] 戴云, 张建中, 钟本善等. 1999. 用Fermat原理作折射波正演算法. 物探化探计算技术, 21(4): 323-326.
[175] 邓志文. 2011. 转换波速度分析方法. 石油地球物理勘探, 46(4): 570-575.
[176] 丁进杰, 戴永寿, 张亚南等. 2013. 基于高频补偿方法提高地震资料分辨率的初步研究. 地球物理学进展, 28(6): 3214-3221, doi: 10.6038/pg20130648.
[177] 杜兴忠, 曹俊星. 2008. 强噪声环境条件下的浅层地震勘探方法研究. 地球物理学进展, 23(5): 1616-1620.
[178] 冯太林, 张学工, 李衍达等. 2001. 折射波地震记录叠加成像方法研究. 地球物理学报, 44(1): 129-134.
[179] 冯智慧, 刘财, 冯晅等. 2011. 基于高阶累积量一维切片的地震信号初至自动拾取方法. 吉林大学学报(地球科学版), 41(2): 559-564, 571.
[180] 付小明. 2012. 地震属性的研究及在工程勘察中的应用[硕士论文]. 成都: 成都理工大学.
[181] 高霞, 张永旺, 魏文博. 2007. 平潭海峡大桥海上地震勘探. 地球科学(中国地质大学学报), 32(4): 540-544.
[182] 郭淑君, 侯延华, 李红立等. 2014. 地震CT技术在城市地质勘查中的应用研究. CT理论与应用研究, 23(1): 81-90.
[183] 何光明, 贺振华, 黄德济等. 2006. 几种静校正方法的比较研究. 物探化探计算技术, 28(4): 310-314.
[184] 华标龙, 刘福田. 1995. 界面和速度反射联合成像──理论与方法. 地球物理学报, 38(6): 750-756.
[185] 黄光远. 2000. 地震反射波成像新方法及其理论基础. 中国科学D辑: 地球科学, 30(6): 650-655.
[186] 姜弢, 林君, 陈祖斌等. 2004. 相控震源地震波定向技术. 吉林大学学报(信息科学版), 22(3): 181-184.
[187] 姜弢, 林君, 杨冬等. 2008. 相控震源定向地震波信号分析. 地球物理学报, 51(5): 1551-1556.
[188] 姜弢, 林君, 贾海青等. 2012. 基于接收阵列的时域地震波束形成方法. 地球物理学报, 55(12): 4277-4287, doi: 10.6038/j.issn.0001-5733.2012.12.039.
[189] 李福中, 邢国栋, 白旭明等. 2000. 初至波层析反演静校正方法研究. 石油地球物理勘探, 35(6): 710-718.
[190] 李辉峰. 邹强, 金文里. 2006. 基于边缘检测的初至波自动拾取方法. 石油地球物理勘探, 41(2): 150-155, 159.
[191] 李红星. 2005. 多波联合应用及同时成像研究[硕士论文]. 吉林: 吉林大学.
[192] 李录明, 罗省贤. 1995. P-SV转换波速度分析及解释方法. 石油地球物理勘探, 30(1): 66-74.
[193] 李庆春, 吴华, 周学明. 2014. 路基下伏地质缺陷弹性波数值模拟. 吉林大学学报 (地球科学版), 44(6): 2016-2025.
[194] 李万万, 裴正林. 2008. 井间地震弹性波传播特征数值模拟. 物探与化探, 32(2): 207-211.
[195] 李小平. 2009. 水域高密度地震映像法在舟山大陆连岛工程中的应用. 工程勘察, (11): 91-95.
[196] 李哲生, 颜志宏, 林榕. 1998. 水域工程物探勘察的新方法——地震映象法. 工程勘察, (3): 70-73.
[197] 李子顺. 2007. 地震波衰减规律及其恢复方法. 地球物理学进展, 22(5): 1545-1551.
[198] 李枝文. 2013. 隧道地质超前预报信号采集与成像方法研究[硕士论文]. 武汉: 长江大学.
[199] 刘保金, 赵成彬, 酆少英等. 2012. 应用三分量浅层地震反射方法探测隐伏活动断裂. 地球物理学报, 55(8): 2676-2686, doi: 10.6038/j.issn.0001-5733.2012.08.020.
[200] 刘福烈, 徐峰, 刘文革. 2014. 不等深延时激发震源组合. 石油地球物理勘探, 49(2): 225-229.
[201] 刘杰, 刘江平. 2013. 基于转换波走时计算倾斜界面转换点位置的搜索法. 石油物探, 52(6): 602-608, 616.
[202] 刘江平, 钱绍瑚, 张素新. 1998. 地震共中心点道集的滑动基准面静校正和动校正叠加理论模型计算和分析. 物探与化探, 22(1): 27-33.
[203] 刘江平, 杨永清, 李修忠. 1999. 复杂地形条件下地震速度提取方法及分析. 物探与化探, 23(4): 258-264.
[204] 刘江平, 张丽琴, 张友明等. 2000. 浅层地震技术在大堤防(隔)渗墙质量检测中的应用. 物探与化探, 24(4): 302-306.
[205] 刘江平, 程飞, 范承余等. 2012. 基于隧道空间全波场二维数值模拟与特征分析. 岩土工程学报, 34(9): 1705-1711.
[206] 刘玲君, 谢中华, 杨萃. 2014. 基于边界线性走时插值的射线追踪算法. 华南理工大学学报(自然科学版), 42(5): 23-28, 35.
[207] 刘洋, 魏修成. 2008. 转换型折射波传播规律及其在转换波静校正中的应用研究. 中国科学(D辑: 地球科学), 38(S1): 204-210.
[208] 刘永华, 田宗勇, 喻振华. 2005. 工程VSP与地震CT联合探测方法及其在岩土工程中的应用. 地球物理学进展, 20(1): 267-272.
[209] 刘祖珉, 雷宛, 鲁红英等. 2006. 浅层地震反射波法在某重大工程预选厂址中的应用. 成都理工大学学报, 33(4): 345-348.
[210] 鲁光银, 熊瑛, 朱自强. 2011. 隧道反射波超前探测有限差分正演模拟与偏移处理. 中南大学学报(自然科学版), 42(1): 136-141.
[211] 陆文凯, 丁文龙, 张善文等. 2005. 基于信号子空间分解的三维地震资料高分辨率处理方法. 地球物理学报, 48(4): 896-901.
[212] 乔宝平, 郭平, 王璞等. 2014. 基于逆虚折射干涉法有效提取近地表弱地震信号. 地球物理学报, 57(6): 1900-1909, doi: 10.6038/cjg20140621.
[213] 邵宇蓝, 李唐律, 廖仪等. 2014. 基于SVI方法的折射波信噪比改善效果研究. 海洋地质前沿, 30(5): 52-57, 69.
[214] 沈阳, 张涛, 赵民等. 2011. 复杂地区浅层地震勘探采集方法探讨. 煤炭技术, 30(6): 253.
[215] 石明, 冯德山, 戴前伟. 2006. 综合物探方法在堤防质量检测中的应用. 地球物理学进展, 21(4): 1328-1331.
[216] 孙渊. 1993. 浅层综合地震勘探方法在隧道工程勘探中的应用. 西安地质学院学报, 15(2): 77-81.
[217] 谭昌勇, 尹成, 林依华等. 2009. 复杂地区综合寻优静校正方法. 石油地球物理勘探, 44(3): 288-291.
[218] 谈顺佳, 张华, 王瑞雪. 2013. 井间地震CT技术在岩溶勘察中的应用. CT理论与应用研究, 22(3): 439-446.
[219] 涂齐催, 刘怀山. 2006. 利用线性走时插值射线追踪计算近地表模型初至波走时. 物探与化探, 30(2): 148-153.
[220] 王朝令, 刘争平. 2012. τ-p变换在隧道反射地震超前预报波场分离中应用的数值模拟研究. 地球物理学进展, 27(5): 2216-2225, doi: 10.6038/j.issn.1004-2903.2012.05.049.
[221] 王朝令, 刘争平, 黄云艳等. 2014. F-K变换在隧道反射地震预报波场分离中的数值模拟研究. 大地测量与地球动力学, 34(4): 71-76, 82.
[222] 王洪. 2007. 折射波法在堤坝病害探查中的应用研究[硕士论文]. 长沙: 中南大学.
[223] 王建军, 廖全涛, 曹建伟等. 2006. 应用井间地震CT探测某桥墩基础断裂. 物探与化探, 30(2): 186-188.
[224] 王建民, 陈树民, 苏茂鑫等. 2007. 近地表高频补偿技术在三维地震勘探中的应用研究. 地球物理学报, 50(6): 1837-1843.
[225] 王荣, 钟世航, 王泽峰. 2013. 陆地声纳法在济南闹市区作断层和岩溶探查. 地球物理学进展, 28(6): 3313-3316, doi: 10.6038/pg20130661.
[226] 王云安, 刘章捷. 2001. 武汉阳逻长江大桥水域工程地震勘探. 地球科学, 26(4): 388-390.
[227] 王云安, 刘章捷. 2001. 武汉阳逻长江大桥水域工程地震勘探. 地球科学, 26(4): 388-390.
[228] 王者江, 何樵登, 田钢等. 2005. 工程反射地震非常规方法野外施工的探讨. 吉林大学学报(地球科学版), 35(6): 796-800.
[229] 魏剑平. 2012. 高精度初至波和反射波联合速度反演[硕士论文]. 西安: 长安大学.
[230] 吴希光, 李亚林, 张孟等. 2012a. 复杂地区地震资料低信噪比的原因及对策(一): 表层散射波是导致地震资料低信噪比的根本原因. 天然气工业, 32(1): 27-32.
[231] 吴希光, 李志荣, 李正佳等. 2012b. 复杂地区地震资料低信噪比的原因及对策(二): 表层散射波的基本性质及形态特征. 天然气工业, 32(2): 38-42.
[232] 吴希光, 李志荣, 张晓斌等. 2012c. 复杂地区地震资料低信噪比的原因及对策(三): 表层散射波压制的技术方法. 天然气工业, 32(3): 57-61.
[233] 吴子泉, 刘元生, 刘保金等. 2005. 地球物理方法在城市地震活动断层精确定位中的应用. 地球物理学进展, 20(2): 528-533.
[234] 肖启航, 谢朝娟. 2012. TST技术在岩溶地区隧道超前预报中的应用. 岩土力学, 33(5): 1416-1420.
[235] 谢冰, 王伟, 吴昌勇. 2007. 浅层地震反射波法在探测水电站进水口淤泥厚度中的应用. 土木基础, 21(2): 76-78.
[236] 谢忠球, 温佩琳. 2000. 复杂目标浅层地震反射法观测技术. 中南工业大学学报, 31(6): 477-480.
[237] 熊章强, 方根显. 2001. 浅层折射波法在隧道工程地质调查中的应用. 水文地质工程地质, (1): 63-66.
[238] 徐峰, 刘福烈, 梁向豪. 2011. 基于相控理论的炮点组合设计技术. 石油地球物理勘探, 46(2): 170-175.
[239] 徐浩. 2013. 基于多分量地震记录的横波速度提取与实例研究[硕士论文]. 武汉: 中国地质大学.
[240] 徐基祥. 2014. 地震干涉测量法近地表散射波分离技术. 地球物理学报, 57(6): 1910-1923, doi: 10.6038/cjg20140622.
[241] 徐涛, 徐果明, 高尔根等. 2004. 复杂介质的折射波射线追踪. 石油地球物理勘探, 39(6): 690-693.
[242] 徐钰, 曾维辉, 宋建国等. 2012. 浅层折射波勘探中初至自动拾取新算法. 石油地球物理勘探, 47(2): 218-224.
[243] 薛荣俊, 李金山, 唐勇. 2000. 潮间带低速层中高速埋伏岩体计算模式的建立. 青岛海洋大学学报, 30(4): 685-689.
[244] 燕利芳. 2011. 反射波层析成像的数值模拟及在西安地裂缝中的应用研究[硕士论文]. 西安: 长安大学.
[245] 张固澜, 熊晓军, 容娇君等. 2010. 基于改进的广义S变换的地层吸收衰减补偿. 石油地球物理勘探, 45(4): 512-515.
[246] 张湖钢, 贺中银, 单力等. 2009. 浅层折射法在陕北侏罗纪煤田松散层调查中的应用及效果. 中国煤炭地质, 21(9): 55-57, 69.
[247] 张建清, 张智, 夏望麟. 2001. 浅层地震多波勘探新技术研究. 工程勘察, (6): 65-68.
[248] 张凯, 汪青松. 2010. 浅层地震反射波法在塌陷勘察中的超前预报. 物探与化探, 34(1): 85-88.
[249] 张钋, 刘洪, 李幼铭. 2000. 射线追踪方法的发展现状. 地球物理学进展, 15(1): 36-45.
[250] 张伟, 王海, 李洪臣等. 2009. 用变换时窗统计能量比法拾取地震初至波. 物探与化探, 33(2): 178-180.
[251] 张学强, 包军强, 刘江平等. 1999. 地震剖面分类与浅层三种波的层位对比. 物探与化探, 23(1): 7-13.
[252] 张学强, 朱翔鹏, 严哲等. 2006. 水上地震映像震源多次波的压制.// 第二届环境与工程地球物理国际会议论文集, 4-6.
[253] 赵春蕾, 卢川, 郝天珧等. 2013. 高精度组合式轻便小型可控震源的研究. 地球物理学报, 56(11): 3690-3698, doi: 10.6038/cjg20131110.
[254] 赵德亨, 田钢, 王帮兵. 2005. 浅层地震折射波法综述. 世界地质, 24(2): 188-193.
[255] 赵后越, 张美根. 2014. 起伏地表条件下各向异性地震波最短路径射线追踪. 地球物理学报, 57(9): 2910-2917, doi: 10.6038/cjg20140916.
[256] 赵永贵. 2007. 国内外隧道超前预报技术评析与推介. 地球物理学进展, 22(4): 1344-1352.
[257] 赵永贵, 蒋辉. 2009. 隧道超前预报技术中的理论问题与评价方法. //中国地球物理学会第25届年会. 北京: 科学出版社, 303-304.
[258] 钟世航. 1997. 陆地声纳法及其应用效果. 物探与化探, 21(3): 172-179.
[259] 钟世航. 2003. TSP作隧道掌子面前方地质预报几例失误原因分析. 岩石力学与工程学报, 22(1): 2443-2446.
[260] 周国婷, 潘冬明, 牛欢等. 2012. 层析静校正方法研究与应用. 物探与化探, 36(5): 802-805, 841.
[261] 周鸿秋, 顾汉明, 张学强等. 1991. 浅层折射波法优化系统在将军崖地质病害调查中的应用. 地球科学——中国地质大学学报, 16(4): 445-452.
[262] 周龙泉, 刘福田, 陈晓非. 2006. 三维介质中速度结构和界面的联合成像. 地球物理学报, 49(4): 1062-1067.
[263] 周竹生, 倪进鑫, 丰赟等. 2010. 浅层地震反射中的固定单元排列法研究. 工程勘察, (7): 84-88.
[264] 朱德兵. 2008. 土坝体隐患弹性波探测中的横波优势. 地球物理学进展, 23(2): 612-616.
[265] 朱金平, 董良国, 程玖兵. 2011. 基于地震照明、面向勘探目标的三维观测系统优化设计. 石油地球物理勘探, 46(3): 339-348.
[266] 左建军, 林松辉, 孔庆丰等. 2011. 井间地震直达波和反射波联合层析成像及应用. 石油地球物理勘探, 46(2): 226-231.