隧道及井巷工程超前探测模拟及其偏移技术研究.pdf

第 2 6卷增 1 2 0 0 7年 7月 岩石力学与工程学报 C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d E n g i n e e r i n g V o 1 ． 2 6 S u p p ． 1 J u ／ y，2 0 0 7 隧道及井巷工程超前探测模拟及其偏移技术研究 张平松 ～，刘盛东 ，吴健生 f 1 ．同济大学 海洋与地球科学学院，上海2 0 0 0 9 2 ；2 ． 安徽理工大学 资源与环境工程系，安徽 淮南2 3 2 0 0 1 摘要隧道及井巷工程中前方地质构造超前探测是一项技术难题，目前国内在基于反射波法探测方面开展了大量 的应用性研究工作，但较少进行探测方法布置及数据处理技术的研究，在实际工程探测中未取得较为理想的效果。 利用多炮检反射波观测系统，通过水槽超声波试验模拟隧道及井巷掘进过程中可能遇到的不同地质条件，进行反 射波超前数据采集，并实现了对数据绕射扫描偏移方法处理。探测模拟结果表明，多炮检反射波观测系统对不同 地质条件产生的波阻抗差异具有明显的反应，其偏移结果与实际界面位置相吻合，但多层介质探测时对后续界面 的分辨率降低 因此必须进一步完善三维多分量立体数据的采集与波动方程数据偏移方法技术的研究，增强对隧 道及井巷工程前方复杂地质构造的探测能力。 关键词。隧道工程；超前探测；超声波试验；绕射偏移；地质构造 中啊分类号l u 4 5 文献标识码A 文章缡号t 1 0 0 0 6 9 1 5 2 0 0 7 增 1 2 8 4 7 0 5 S TUDY ON DETECTI NG S I M ULATI ON AHEAD OF TUNNEL AND L ANEⅥ Y AND I TS MI GRAT I ON T E CHNI QUE S ZHANG Pi ng s o n g1 ～ ， L I U S h e n g d o n g ，W U J i a n s h e n g 1 f 1 ． S c h o o l o fO c e a na n d E a r t h S c i e n c e ，T o n g j i U n i v e r s i t y ， S h a n g h a i 2 0 0 0 9 2 ，C h i n a ；2 ． De p a r t me n t o f R e s o u r c e s a n dE n v i r o n m e n t a l E n g i n e e r i n g ，A n h u i U n i v e r s i t y o fS c i e n c e a n d T e c h n o l o g y ，Hu a i n a n ，A n h u i 2 3 2 0 0 1 ，C h i n a Abs t r a c t Th e d e t e c tio n a n d f o r e c a s tin g f o r g e o l o g i c a l s t r u c t ur e a h e a d of t h e tun ne l an d l a ne wa y a r e the o n e o f the mo s t d i ff i c u l t i s s u e s i n the wo r l d ． M u c h a p p l i c a t i o n wo r k i s p e rfo r me d wi th the s e i s mi c r e f l e c ti o n wa v e me tho d i n Ch i n a a t p r e s e n t ． Bu t the s t u d y o n d e t e c ti n g o b s e r v a t i o n s y s t e m i n 3 D s p a c e and d a t a p r o c e s s i n g t e c h n i q u e s i s n o t p e rf o rm e d c o n s i d e r a b l y ． T h e p e rf e c t e x p l o r a ti o n e ff e c t i s n o t g a i n e d i n the p r a c ti c al d e t ect i n g p r o j ect s ，and thi s s i tua ti o n i s n o t g o o d for the e x t e n d i n g a p p l i c a t i o n o f the a d v a n c e d f o r e c a s ti n g t e c h n i q u e s o f t u n n e l and l a n e wa y ． Th e diffe r e n t d e t e c ti n g c o n d i ti o n s tha t wi l l b e me t ah e a d o f the t u n n e l and l a n e wa y t h r o u g h u l t r a s o n i c e x p e r i me n t s i n the fl u me a r e p r e s e n t e d ．Th e f o r e c a s t i n g d a ta u s i n g mu l t i s h o t - r e c e i v e r s y s t e m o f the r e f l e c ti o n wa v e a r e c o l l e c t e d ，a n d the r e f r a c ti o n s c a n n i n g mi g r a ti o n for the c o l l e c ti n g da t a i s r e ali z ed ． Th e s i mu l a ti o n r e s u l t s s h o w tha t the mu l ti s h o t - r e c e i v e r s y s t e m h a s o b v i o u s r e a c ti o n o n the wa v e i mp e d an c e o f the i n t e rfa c e s i n d i ffe r e n t g e o l o g i c al c o n d i ti o n s ． T h e mi g r a ti o n r e s u l t s a r e i n a c c o r d an c e wi th tho s e o f p r a c ti c al i n t e rfa c e s ，b u t the d i s t i n g u i s hin g a b i l i t y for the b a c k i n t e rfa c e s wi l l b e u n s a ti s f a c t o r y wh e n mu l ti l a y e r me dia a r e d e t e c t ed ． S o the 3 D mu l t i c o mp o n e n t and s t e r e o s c o p i c d a t a c o l l e c ti n g me tho d s a l e p r o p o s e d，and wa v e e q u a ti o n mi g r a ti o n t e c h n i q u e s a r e d e v e l o p e d t O i mp r o v e e x p l o r a t i o n a b i l i t y for c o mp l e x g e o l o g i c a l s t r u c t u r e s ah e a d o f the t u n n e l an d l a n e wa y ． Ke y wo r d s tun n e l i n g e n g i n e e ri n g ；ah e a d d e t e c ti n g；u l tra s o n i c e x p e rime n t ；r e f r a c ti o n s c a n n i n g mi g r a ti o n g e o l o g i c al s t r u c t u r e 收藕 日期I 2 0 0 6 0 32 5 修目 日期I 2 0 0 6 0 51 5 作者筒介。 张平松 1 9 7 1 一 ，男，1 9 9 6 年毕业于安徽理工大学资源与环境工程系地质工程专业，现为博士研究生、副教授。主要从事地质工程、 工程物探方面 的教学与研 究工作 。E - ma i l p s z h a n g 7 1 1 6 3 ． t o m 维普资讯 岩石力学与工程学报 2 0 o 7年 1 引 言 在隧道及井巷工程施工过程 中，往往会因为地 质条件不详而遇到大量的影响安全生产的地质问 题，其中断层构造 、软弱层及其含水异常体是所要 解决的重点问题。对巷道前方地质异常体的位置及 特征进行探测与评价称为超前探测。 目前用于超前 探测的方法很多，主要有直接钻探法 、地震波法、 井巷 电阻率法和 电磁波法等【 l 4 ] 。但是在诸多的物 探方法中，利用反射地震波方法进行超前探测是主 要的，也是准确率最高的一种地球物理探测方法。 具体来说，依据反射波原理进行超前探测的方法主 要有隧道地震剖面系统t5 - 8 1 [1 T S P系统 、负视速度 法与隧道垂直地震剖面法 即 T V S P法 、水平声 波剖面 即 HS P法 、真反射成像法【 l 即 T I 法 、综合地震成像系 即 I S I S系统 、陆地声纳 法【 l 也叫高频地震反射法 以及隧道地震层析成像 法u 即 T S T法 等。其中多数方法是从国外直接引 进的，目前 由瑞士 Amb e r g公司完成的 T S P系统探 测准确率相对较高，主要应用在 隧道超前探测中， T R T和 I S I S方法国内还鲜见有相关应用文章 。总的 来说，目前国内开展地震波超前探测都是以技术应 用为主，除何振起等[ 1 6 - 1 8 ] 进行了负视速度法的理论 及处理算法研究外，对隧道及井巷中超前探测技术 原理研究较少，特别是对反射地震波法的空间数据 采集与处理技术没能有效地深入 。本文试图结合隧 道及井巷工程超前探测模型，进行超声波数据采集 与偏移技术研究，以求获得实际探测应用的基础内 容，为国内开展超前探测方法技术研究提供一些思 路。 2 超前探测模拟及数据采集 超前探测模拟试验是在水槽 中进行的，结合隧 道及井巷工程实际可能遇到的地质条件，模拟中分 别选取不 同材料代替前方 的高波速 、低波速 目标 体，进行超声波数据采集与处理 。 2 ． 1 模拟试验 图 1为多炮检反射波超前探测水槽试验布置 图。该观测系统是在 T S P方法基础上建立，采用多 激发多接收系统布置 。 试验利用 WS D一1 型超声波 仪进行数据采集，其发射探头作为激发源，三分量 图 1 多炮检反射波超前探测水槽试验布置图 F i g ． 1 S c h e mi c d i a g r a m o f e x p e r i me n t i n flu me u s i n g mu l ti-- s h o t --r e c e i v e r s y s t e m o f r e fl e c t i o n wa v e 加速度检波器进行数据接收。 发射探头按照设定的移动步距在中心线上滑动 并发射超声波 ，来模拟多炮激发。探测地质界面是 以方形玻璃、方形泡沫塑料 、球形泡沫塑料等材料 分别来进行模拟 。其中界面的角度 可按要求进行 调整， 目的是模拟隧道及井巷工程中前方 目标地质 体与巷道中心线夹角不同所产生的差异。图2所示 为多炮检反射波法隧道及井巷观测系统图。 图 2 多炮检反射波法隧道及井巷观测系统图 F i g ． 2 S c h e ma ti c dia gr am o f o b s e r v a t i o n s y s t e m u s i n g mu l ti - s h o t --r e c e i v e r s y s t e m o f r e fl e c ti o n wa v e a h e a d o f t u n n e l a n d r o a d wa y 2 ．2 数据采集 试验过程 中对坑 道掘进前方可能遇到 的多种 地质条件进行了模拟测试，这里仅选择一部分内容 进行说明。表 1为模拟试验 中探测布置有关参数及 说 明。 数据采集时，1 样 ， 接收点偏移距分别为 0 ．2 5 和 0 ． 1 5 II 1 ，发射探头移动步距为0 ．0 1 I I 1 ，从而获得 不同模拟条件下的多炮多接收点记录 。 3 数据处理及偏移 对于隧道及井巷 内反射地震波探测数据处理， 维普资讯 第2 6 卷 增1 张平 松 笠． 堕 丛 鳌 猩塑萱堡 垫丛基 垫 塑 墨 塑 表 1 不同探测条件及其布置参数 T a b l e 1 Di ff e r e n t d e t e c ti n g c o n d i ti o n s a n d i t s l a y o u t p a r a me t e r s 除采用负视速度法等传统方法外，可直接对数据进 行偏移处理[ 2 。偏移技术可分为两类 ，一类是射 线偏移法 ，以Hu y g e n s 原理为基础 ，另一类是波动 方程偏移法，以物理光学原理为基础。 目前从适用 角度来说，射线扫描偏移法较易实现，其偏移结果 能满足实际要求。本文应用 的是基于射线扫描的偏 移方法。 3 ． 1 绕射扫描偏移 绕射扫描偏移是建立在射线偏移的基础上使反 射波 自动归位到真实位置上的一种方法。根据 H u y g e n s原理 ，地下每一个反射点 P都可以看成是 一 个子波震源 ，进行绕射扫描偏移时，把每一个 网 格点看成是一个反射点，则其反射波或绕射波旅行 时为 式中v为地震波的速度 r r d ms ；h为点P的垂直深 度 m ；t 川 为扫描点P处第 f 炮第_『 个接收点的绕 射波旅行时 ms ，J1 ，2 ，3 ，⋯ ，m，且 m为参 与叠加的记录道。图 3为绕射扫描偏移网格图。这 样把记录道上t 时刻的振幅值a 与点P 的振幅值叠 加起来，作为点P的总振幅值 A，则有 A i ， 2 j l 当对 X h平面按 A x，A h划分的方格 网上每一 点p x ， 都进行计算，只要划分得足够细，总可 以 在所要求的精度上反映反射点的全部可能位置。这 样，使反射界面上的叠加扫描点P的总振幅A更加 增大，不在反射界面上的扫描点P的总振幅 A进一 步相对减小，既提高了信噪比，又把反射界面 自动 图3 绕射扫描偏移网格图 F i g ． 3 Me s h ma p o f r e f r a c ti o n s c a n n i n g mi g r a t i o n 偏移到其空间真实位置上去 。在进行偏移叠加时， 还要进行参数选择，如速度、扫描范围的选择等。 为了提高运算速度，对每一个炮检对进行振幅叠加 时要控 制其扫描范 围，对于横 向可选择炮检距 的 5 ～6倍范围，纵向选择设定的深度范 围进行扫描。 为了提高超前探测对前方地质构造的分辨率，扫描 时增加了每一炮检对 的扫描点与 X轴夹角控制 ，其 范 围为 4 5 。 ～6 0 。 。 3 ． 2 处理过程与实现 多炮检反射波探测数据处理具有一定的流程 ， 通常从原始记录开始，包括预处理、带通滤波、初 至波拾取、反射波提取、纵横波分离、速度分析、 扫描偏移 以及对各个反射界面 的提取过程。同时必 须对偏移范 围进行设定，反射波提取时上下行波分 离可采用二维滤波或 p域滤波实现，获得上行波。 目前对该数据的处理是通过自行编制的软件进行实 现的。 4 偏移结果及认识 数据采集时采用三分量传感器 ，可接 收来 自不 同方向的超声波。图4为多层介质条件下 1 接收点 X分量波形记录。按处理流程进行数据处理；图5为 提取 的上行波记录；图 6为该记录 的偏移结果。通 过对 5 种不同条件下的超前模拟数据进行偏移叠加 处理，提取测试中心线处的偏移结果进行对比 见 图7 ，其坐标原点为 1 接收点位置。 对比分析可知图 7中未放任何介质 的为水槽 测试区域背景数据偏移结果 见图7 a ，其后续的强 相位是 由水槽边界反射所致 ，这为后续 的对 比解释 提供了基础。 与背景值相比较， 绕射点测试结果 见 图 7 b 在 0 ． 4 5 5 m处具有明显的起跳 ，因此可 以对 ． ， f 维普资讯 2 8 5 0 岩石力学与工程学报 2 0 0 7焦 时间／ ms 图4 水槽两层介质超前探测 1 接收点X分量记录 F i g ． 4 Re c o r d o f x - d i r ecfi o n o f r e c e i v e r No ． 1 f o r t wo l a y e r 时间／ ms 图5 1 接收点 X分量上行波分离结果 F i g ． 5 Re s ult o f x - d i r e c t i o n o f u p wa r d wa v e r e c o r d o f rec e i v e r No ． 1 L ／ m 图 6 1 接收点X分量偏移结果图 F i g ． 6 M i g r a t i o n r e s u l t s o f x - d i r e c t i o n o f r e c e i v e r No ． 1 0 ． 0 0．2 0 ． 4 0 ． 6 0．8 1 ． 0 a ￡／ m 0 ． 0 0 ． 2 0 ．4 0 ． 6 0 ． 8 1 ． 0 0 ． 0 0 ． 2 b ￡／m c ￡／ m 0．0 0 ．2 0 ． 4 0 ． 6 0 ． 8 1 ． 0 d ￡／ m e 图7 水槽不同介质条件下偏移结果对比分析图 F i g ．7 Co mp a r i s o n o f mi g r a t i o n r e s u l ts wi t h d i ffe r e m me sa 绕 射点进 行分辨 。 比较 图 7 a 和 c 可 以看 出在 0 ． 5 3 0 m处存在强反射， 而 图 7 a 中相 同位置处没有 看到异常相位，因此可判断高速玻璃板的位置，但 在图 7 c 中存在后续的强相位，分析为受多次反射 波影响的结果。对于单层低速介质来说，试验模型 中泡沫板与绕射点是放在同一位置，虽然位置基本 相 同但其偏移结果 中的反射相位较绕射点的强 见 图 7 d 。图 7 e 为两层介质的偏移结果，与单层介 质模型偏移结果 图 7 c 和 d 相 比，0 ． 4 5 5 m处反射 特征清晰，而第 2层高速介质反射特征表现不够明 显 ，但与图 7 d 相 比，在 0 ． 5 3 0 m处有起跳反应 ， 同时后续水槽边界反射特征也削弱 。该试验结果表 明，利用多炮检观测系统进行超前反射波数据采集， 所模拟的前方不 同介质体在偏移结果图中特征表现 明显，可进行准确 的位置判断。 5 结论 通过水槽超声波模拟试验数据及其偏移结果表 明，在隧道及井巷 中可以利用多炮多接收反射波法 对前方地质构造及其异常进行探测与位置确定，且 能获得较好的效果，受篇幅限制本文对现场实际探 测应用未作具体说明。但与实际地质条件相 比，二 维数据采集方法 以及扫描法偏移处理还存在诸 多问 题，比如立体空间中单一测线数据采集对构造结构 面的归位效果较差 ，绕射扫描偏移对复杂结构探测 精度和成像时的波形特征都不够准确。因此进行三 维多分量立体数据采集并对数据进行波动方程叠前 偏移 ， 包括 K i r c h h o ff积分偏移、有限差分偏移和有 限元偏移等应用将会大大提高成像质量与效果【2 。 目前 T S P也仅仅实现 的是二维数据 的扫描偏移【 5 】 ， 因此其对 不 同地质体 的分辨 能力也存在一定 的问 题。同时加强对 多层结构面重叠 的特征分辨，对有 效探测距离的确定，以及与隧道和井巷中心线不同 夹角地质构造反应特征的研究程度，为实际探测应 用提供必要的基础资料，这将是反射波法超前探测 需要深入研究的方向。 参考文献 R e f e r e n c e s 【 1 】 赵永贵，刘浩，孙宇， 等．隧道地质超前预报研究进展叨．地 球物理学进展， 2 0 0 3 ， 1 8 3 4 6 04 6 4 ． Z H AO Y o n g g u i ， L I U H a o ， S U N Y u ,e t a l ，R e s e a r c h p r o g r e s s i n t u n n e l g e o l o g i c a l p r e d i c t io n 叨 ． P r o g r e s s i n G e o p h y s i c s ，2 0 0 3 ，1 8 3 4 6 0 4 6 4 ． i n C h i n e s e 维普资讯 第 2 6 卷 张平松，等．隧道及井巷工程超前探测模拟及其偏移技术研究 【 2 】 【 3 】 [ 4 】 【 5 】 [ 6 】 【 8 】 【 9 】 【 l 0 】 【 1 2 】 何继善， 柳建新．隧道超前探测方法技术与应用研究叨．工程地球 物理学报，2 0 0 4 ，l 4 2 9 3 2 9 8 ． H E J i s h a n ，L I U J i a n x i n ． A s t n d y o f t u n n e l a d v a n c e d d e t e c ti o n t e c h n o l o g y a n d i t s a p p l i c a ti o n [ J ] ， C h i n e s e J o u r n a l o f E n g i n e e r in gG e o p h y s i c s ， 2 0 0 4 ， l 4 2 9 32 9 8 ． i n C h i n e s e 钟宏伟， 赵凌．中国隧道工程超前预报技术现状分析[ J ] ．人民长 江 ， 2 0 0 4 ，3 5 9 1 51 9 ． Z H ON G Ho n g we i ， Z HA O L ing ． A n a l y s i s o f s ta tus q u o o f f o r wa r d p r e d i c t i o n t e c hn o l o g y o f Ch i n e s e t u n n e l w o r k s 们． Y ang t z e R i v e r ，2 0 0 4 ，3 5 9 1 51 9 ． i n C h in e s e 陈建峰．隧道施工地质超前预报技术比较叨．地下空间，2 0 0 3 ， 2 3 1 58 ， 1 6 ． C H E N J i a n f e n g ．C o mp a r i s o n o f t e c h n i q u e s o f g eol o g i c al p r e d i c ti o n o f tu n n e l s [ J ] ． Un d e r g r o u n d S p a c e ， 2 0 0 3 ， 2 3 1 5 8 ， 1 6 ． i n C h i n e s e A mb e r g Me a s u r i n g T e c hni q u e L t d ． ．O p e r a ti o n m a n u al 1 9 ． 6 1 [ R ] ． S wi s s Amb e r g M e a s u r i n g T e c h n iq u e Lt d ． ， 2 0 01 ． Amb e r g M e as u r i n g T ec h n i q u e Ltd．．Co mp a r i s o n o f v a r i o u s t u n n e l g eol o gic al for e c ast t ech n i q u e s [ R ] ．S wi s s A mb e r g Me asu r i n g Tech n i q u e Lt d ． ， 2 0 0 3 ． 刘志刚，刘秀峰． T s P 隧道地震勘探 在隧道隧洞超前预报中的应 用与发展[ J ] ．岩石力学与工程学报，2 0 0 3 ，2 2 8 1 3 9 91 4 0 2 ． L I U Z h i g ang ，L I U X i u f e n g ． T S P a p p l i c a ti o n and d e v e l o p m e n t i n tunn e l l e a d f o r e c ast [ J ] ．C h i n e s e J o u r n al o f R o c k Me c h a n i c s and E n g i n eer i n g ，2 0 0 3 ，2 2 8 1 3 9 91 4 0 2 ． i n C h i n e s e 戴前伟， 何刚， 冯德山．T S P一 2 0 3 在隧道超前预报 中的应用[ J ] ． 地球物理学进展， 2 0 0 5 ， 2 o 2 4 6 0 4 6 4 ． DA I Qi a n w e i ， HE G ang ， F E N G D e s h a n ． Ap p l i c a ti o n o f t h e T S P一2 0 3 s y s t e m i n g e o l o gic al a d v anc e d p r e d i c t i o n o f tunne l [ J ] ． P r o gre s s in G e o p h y s i c s ，2 0 0 5 ， 2 0 2 4 6 0 4 6 4 ． i n C h i n e s e 程增庆，吴奕峰，尚晶，等． 探测掌子面前方断层的 T V S P方 法们．石油地球物理勘探 ，1 9 9 4 ，2 9 增 2 1 5 8 1 64． C HE NG Ze n g q ing ， W U Yi f e n g， S HANG J i n g，e t a1． Th eTVS Pme tho dfo r d e t e c t i n g the f a u l t s a h e a d o f tu nn e l [ J ] ． G e o p h y s i c al E x p l o r a ti o n f o r P e t r o l e u m， 1 9 9 4 ，2 9 S u p p ． 2 1 1 5 81 64 ． in C h i n e s e BRUCKL E， CHWAT AL W ， DOLZ LMULLE R J ， e t a1． A s tud yo fthe a p p li c a ti o n o f VS Pt o e x p l o r a ti o nahe a d o f atunn e l [ J ] ． I me ma ti o n al J o u r n a l o f Ro c k Me c h ani c s and Mi n i n g S c i e n c e s ，2 0 0 1 ，3 8 6 8 3 38 41 ． I NAZAKI T， I S AHAIH， KAW AMURA S， e ta1． S t e p wi s ea p p l i c a ti o n o f h o r i z o n t a l s e i s mi c p r o f d i n g for tunn e l p r e d i c t i o n ah e a d o f the f a c e 们 ． Th e L e a d i n gE d g e ，1 9 9 9 ， 1 8 1 2 1 4 2 9 1 4 3 1 ． OTTO R， BUTI N E A， BR ET TE REBNER H， e t a1． Th e a p p l i c a ti o n f 1 3 】 【 l 4 】 【 l 5 】 【 1 6 1 【 l 7 】 【 1 8 】 【 1 9 】 f 2 l 】 【 2 2 】 o f T R T t r u e r e fl ecti o n t o mo g r a p h y a t the U n t e r wal d tu nn e l [ J]． G eop h y s i c s ，2 0 0 2 ，2 0 2 5 1 5 6 ． BORM G， GI ES E R， KLOS E C D， e t a1． I S I S int e gr a t ed s e i s mi c i ma g ing s y s t e m for the g e o l o gic al p r e d i c ti o n ahe a d o f u n d e r gro u n d c o n s t r u c ti o n [ C ] ／ ／ E A G E the 6 5 t I l C o n f e r e n c e and E x h i b i ti o n ． No rw a y 【 s ． n ． 】 ，2 0 0 3 ． 钟世航．陆地声纳法的原理及其在铁路地质勘测和隧道施工中的 应用[ J ] ．中国铁道科学，1 9 9 5 ，l 6 4 4 8 5 5 ． Z H ON G S h i h a n g ． Th e p r in c i p l e an d a p p li c a ti o n o f l an d s o n ar me thod i n r a j l wa y g eol o gic al s u r v e y and t u nne l c o n s t r u c ti o n [ J ] ．C h i n a R a i l wa y S c i e n c e ，1 9 9 5 ，1 6 4 4 8 5 5 ． i n C h in e s e ZHAO Y G， J ／ ANG H，ZHAO X P T u n n e l s e i s mi c to mo g r a p h y me thod for g e o l o g i c al p r e d i c ti o n and i ts a p p li c a ti o n [ J ] ．Ap p li ed Geop h y s i c s ，2 0 0 6 ，3 2 6 9 7 4 ． 何振起 ，李海，梁颜忠．利用地震反射法进行隧道施工地质超 前预报[ J ] ． 铁道工程学报， 2 0 0 0 ， 4 8 1 8 5 ． HE Z h e n q i ， L I H a i ， LI ANG Y a n z h o n g ．Geo l o g i c al s u pe r - l e a d i n g for e c as t d u r i n g tunn e l c o n s t r u c ti o n b y u t il i z ing s e i s mi c r e s p o n s e anal y s i s me thod [ J ] ． J o u r n al o f R a i l wa y E n g i n e e ri n g S o c i e t y ，2 0 0 0 ， 4 8 1 8 5 ． i n C h i n e s e 武思玉， 何振起．隧道超前预报叨．铁道知识，1 9 9 6 ， 5 1 4 ． WU S i y u ，HE Z li e n g q i ． F o r e c ast ahe a d o f t u n n e l [ J ] ． K n o w l edg e o f th e R a i l wa y ， 1 9 9 6 ， 5 1 4 , i n C h i n e s e 曾昭璜．隧道地震反射法超前预报 [ J ] ．地球物 理学报 ，1 9 9 4 ， 3 7 2 2 6 8 2 7 1 ． Z E N G Z h a o h u a n g ． P r e d i c ti o n ahe a d o f the t u n n e l f a c e b y th e s e i s mi c r e fl e c t i o n me tho d s [ J ] ． A c ta G e oph y s i c a S i n i c a ， 1 9 9 4 ，3 7 2 2 6 8 2 7 1 ． i n C h i n e s e 牟永光．地震勘探资料数字处理方法【 M】 ． 北 京 石油工业出版社， 1 9 8 1 ． MOU Y o n g g u ang ．Me tho d s o f d a ta p r o c e s s ing for s e i s mic e x p l o r a t i o