陷落柱地震波超前探测数值模拟与应用.pdf

第 2 4卷 3期 2 0 1 2年 3月 中国煤炭地质 C0AL GE0L0GY 0F CHD A Vo 1 ． 2 4 No ． 03 Ma r ． 2 01 2 d o i 1 0 ． 3 9 6 9 6 ．i s s n ． 1 6 7 4 - 1 8 0 3 ．2 0 1 2 ． 0 3 ． 1 5 文章编号 1 6 7 4 1 8 0 3 2 0 1 2 0 3 0 0 5 6 0 4 陷落柱地震波超前探测数值模拟与应用 王 勃 ， 刘盛东 ， 胡泽安 1 , 2 1 ． 中国矿业大学 资源与地球科学学院 ， 江苏 徐州2 2 1 1 1 6 2 ． 中国矿业大学 深部岩土力学与地下工程 国家重点实验 室． 江苏 徐州2 2 1 1 1 6 摘要 针对陷落柱与周 围煤体之间存在 明显波阻抗差异的特点， 建立巷道 内地震 波超前探测陷落柱数值模 型。采用 波场正演模拟方法 ， 分析巷道前方存在陷落柱 时的地震波场特征 ， 总结 出基于绕射偏移的陷落柱边界探测方法。在 新元矿南二正巷实 际超前探测中 ， 该方法准确 圈定 了巷道前方 陷落柱的范围。探测结果 和实际揭露情况对 比表 明 地震波超前探测技术可 以对巷道前方 1 2 0 m范围内的陷落柱构造进行有效预测。 关键词 地震波超前探 ； 陷落柱 ； 绕射偏移 ； 数值模拟 中 图分 类 号 P 3 6 1 ． 4 文 献 标 识 码 A S u b s i d e d Co l u mn S e i s mi c W a v e Ad v a n c e d P r o s p e c t i n g Nu me r i c a l S i mu l a t i o n a n d Ap p l i c a t i o n W a n g Bo 一．L i u S he n g d o n g a n d Hu Ze a n ， 1 ． S c h o o l o f Re s o u r c e a n d Ge o s e i e n c e s ， C UMT， Xu z h o u ， J i a n g s u 2 21 1 1 6 ； 2 ． S t a t e K e y L a b o r a t o r y f o r G e o me e h a n i e s ＆D e e p U n d e r g r o u n d E n g i n e e ri n g 。 X u z h o u ． J i a n g s u 2 2 1 1 1 6 Ab s t r a c t I n a l l u s i o n t o t h e o b v i o u s wa v e i mp e d a n c e d i ff e r e n c e e x i s t e d b e t we e n s u b s i d e d c o l u mn a n d p e rip h e r a l c o a l ma s s ．t h e n u me ri c a l mo d e l of s u b s i d e d c o l u mn r o a d wa y s e i s mi c w a v e a d v a n c e d p r o s p e c t i n g h a s b e e n mo d e l e d ．T h e a n a l y s i s u s i n g wa v e f i e l d f o r wa r d s i mu l a t i o n h a s f o u n d r o a d w a y f r o n t s u b s i d e n c e s e i s mi c w a v e f i e l d f e a t u r e s e x i s t e d ，a n d t h e n a s u b s i d e d c o l u mn b o u n d a r y d e t e c t i o n me t h o d b a s e d o n d i f f r a c t e d mi g r a t i o n s u mme d u p ． D u ri n g t h e a d v a n c e d p r o s p e c t i n g i n t h e S 2 ma i n r o a d wa y o f t h e Xi n y u a n c o a l mi n e ． u s i n g o f t h e me t h o d h a s c y c l e d o u t s u b s i d e d c o l u mn e x t e n t a h e a d o f r o a d w a y w i t h a c c u r a c y ． T h e c o n t r a s t o f p r o s p e c t i n g a n d p r a c t i c a l r e v e a l e d r e s u l t s d e mo n s t r a t e d t h a t t h e s e i s mi c w a v e a d v a n c e d p r o s p e c t i n g t e c h n o l o g y c a r l e f f e c t i v e l y p r e d i c t r o a d w a y f r o n t s u b s i d e d c o l u mn wi t h i n 1 2 0 m e x t e n t ． Ke y wo r d s s e i s mi c w a v e a d v a n c e d p r o s p e c t i n g ； s u b s i d e d c o l u mn ； d i ffr a c t e d mi g r a t i o n 0 引言 岩溶陷落柱是煤矿开采过程中常常遇到的一种 灾害性地质异常体 ，其普遍分布于我国华北石炭一 二叠纪煤 田中，是影响煤矿生产建设的特殊地质条 件之一。 岩溶陷落柱不仅破坏了煤层的连续性 ， 影响 煤矿 生产 ， 而且 还有 可能将 地下 水导 人工作 面 内 ， 给 煤矿安全生产造成极大威胁。 一 般情况下，探查岩溶陷落柱都是运用物探先 行 ， 物探 、 钻探 、 巷探相结合的技术手段进行综合探 查 ， 以确定其确切位置， 圈定范围 ， 查明形状及规模。 目前探测岩溶陷落柱的主要物探手段主要有地面地 震勘探 、 矿井无线电波透视技术 、 矿井地质雷达、 瑞利 波勘探技术、 地震波 C T和槽波地震勘探技术等 。 本文介绍地震波超前探测技术 ，结合 了巷道空 间特点 。采用 了巷道多次覆盖观测系统进行数据采 作者简介 王勃 1 9 8 4 一 ， 男 ， 安徽安庆人 ， 博士研究生 ， 主要从事矿 井 巷道超前探数值模拟与实际探测研究工作 。 收稿 日期 2 0 1 2 0 3 0 1 责任编辑 孙常长 集 ， 数据处理过程 中综合运用了波场分离 、 反射波提 取 、 叠前偏移成像等多种处理算法 ， 是一种多波多分 量全 空 间地震勘 探技 术 。 由于 井 中地震勘 探是 在巷 道 中激 发 和接收 地震 波 ， 既克服了复杂的表层干扰 ， 又离 目标体较近 ， 因 而 地震波 能量 和高 频成分 衰减 较少 ，可 获得高 精 度 和高分辨率的地震资料。为掌握井中地震的波场特 征 ， 本文采用波场正演模拟方法 ， 研究巷道内激发可 能产生的波场类型 包括有效波和干扰波 和传播特 点 ， 并在实践中验证 。实践证明， 地震波超前探测技 术 能有 效探 测掘 进 巷道 前 方 1 2 0 m 范 围 内断层 等地 质异常 ， 其探测成果可减少掘进事故 ， 保障矿井安全 高效生产 o _ 3 1 。 1 地震波超前探测技术及 系统原理 1 ． 1巷 道前方 高倾 角界 面负视 速 度时距 特征 研究反射波的时距曲线特征需首先建立观测系 统 。 由于受地下工程空间局限性的影响． 巷道地震探 测工作只能在巷道迎头附近有限区域展开。 因此 ． 结 3期 王 勃 ， 等 陷落柱地震波超前探测数值模拟与应用 5 7 合实际，坐标系的 轴原点位于巷道后方 ， 轴正 方向水平指向巷道前方 ；接收点分别设置在原点和 迎 头位置 ， 激 发点位 于原 点与 迎头 之间 。 设巷 道前 方 界面在 轴上出露点为 h， 0 ， 界面倾角为 0 ， 后置 接收点位于原点 O ， 0 ，任一激发点的坐标为 图 1 。 ．巷道 图 1 巷道 前 方 倾 斜 界 面 与 炮检 相 对 位 置 F i g u r e 1 Re l a t i v e po s i t i o n a m o n g r o a dwa y f r o n t s l a n t bo u n da r y p l a ne ，s ho t p o i n t a n d d e t e c t o r 根据反射定律和虚震源原理可得巷道前方界面 的时距 方程 如下 ．一R 一 v ／ x ％4 h h - x s i n 2 0 1 一 t ， 1 ’ 式 中 ￡ 『__后置接收点反射旅行时， S ； 任一激发点在 轴上的坐标 ， m； 前方界面和 轴 巷道轴线 的交点坐标 ， m ； 界面倾角，o ； 】 界 面上层 速度 ， m ／ s ； 由于巷道中仅能布置有限测线 ，因此仅能反映 出时距 曲线 的局 部 。图 2为巷 道前方 不 同倾 角 界面 在后 置接 收点 的时距 曲线 ，对 比不 同角度 的时距 曲 线可 以发现 ， 当界面为高倾角时 4 5 。 ～ 9 0 。 ， 由于最 小值点不在测线 内，测线只能反映出单调递减的反 射波 同相轴 ， 因此均表 现 出负视 速度 。 尤其 直立 界 面 0 9 0 。 时距 曲线 已简化为直线 ， 其斜率为一 1 ／ v 。 综上所述 ， 受巷道测线长度限制 ， 在前方界面的 反射波场中只能反映出局部反射波同相轴 该局部 同相轴和界面倾角大小相联系 ，当前方界面为高倾 角时 ， 时距曲线段表现出负视速度特征； 测线上前后 置接收点时距曲线特征相似 ．在高倾角界面时均具 有负视速度特征。 1 ． 2反射 波提 取 从以上分析可知，巷道条件下地震波场是异常 复杂 的。地震波超前探测的目标体主要集中在巷道 _～ z z ． s t 二 王 ≥ ～ ～1 5 a 一～～～ 一一 一 一 二 j ， ． ， 一，0 。 ＼√ ／／ g ＼ ＼ ， 、 、 』～ ． h ， 0 图 2 巷 遭 前 方 不 l司倾 角 界 面后 置 接 收 点 时 距 曲 线 图 Fi g u r e 2 Ti me g r a p h o f r o a d wa y f r o nt d i ffe r e nt d i p bo u n da r y p l a ne p o s t p o s i t i o n r e c e i v i ng p oi n t s 前方 ，而前方的地质信息完全体现在来 自巷道前方 反射波组中。 因此 ， 有效地提取巷道前方界面的反射 波成为地震波超前探测的关键技术之一。 从时距曲线规律看 出， 在线性观测系统条件下 ， 巷 道前 方界 面均 表现 出特殊 的负 视速 度特征 。本 文 基 于 此视 速 度差 异 规律 。 采 用 一 P滤波 方 法提 取 出 巷道前方界面反射波 ， 同时压制声波 、 面波等干扰波 以提高数据信噪 比。 2 超前探测数值模拟 地震波场数值模拟是研究复杂地区地震资料采 集 、 处理和解释的有效辅助手段。 为了有针对性地研 究巷道周边界面波场传播规律 ．通过合成地震记录 方式来研究地震超前探测技术对于前方陷落柱的探 测 效果【 。 2 ． 1 模 型 参数 模拟研究 中设计一个大小为 5 0 0 mx 5 0 0 m 的三 层数值模型。中间层为模拟煤层 ，上下为模拟顶底 板 ， 其 中模拟参数见表 1 。 网格划分为 0 ． 5 m x O ． 5 m。 巷 道前方存在一个陷落柱 ， 界面倾角为 9 0 。 。炮点和检 波点设置在巷道 同一侧帮煤层中部。 炮点在前方 ， 共 l 0炮， 后置 1 个检波器。 检波点坐标为 2 3 0 ， 2 5 0 ， 起 始炮 点坐标 为 2 6 0 ， 2 5 0 ， 炮距 5 m， 共 1 0炮 ， 末 端炮 点坐标为 3 0 5 ， 2 5 0 。 其空间位置关系见图 3 数值模 型 示意 图 。 58 中 国 煤 炭 地 质 第2 4 卷 表 1 巷道超前探测模 型参数 Ta b l e 1 Ro a dwa y a d v a n c e d pr o s p e c t i n g m o d e l pa r a m e t e r s 编号 名称 纵 ， 。位置 方向 ， m范围 z方 向 ／ m ／ m’ s ， Kg i n 图 3 数值 模 拟 模 型 示 意 图 Fi g u r e 3 A s c h e ma t i c d i a g r a m o f nu m e r i c a l s i m u l a tio n m o d e l 震 源类 型选用 点源 ，采用 波动方 程模 拟 的地震 波场包含了地震波传播的所有信息。为研究地震波 的传播机理和复杂地层 的解释提供了更多的佐证。 记录合成时选用雷克子波， 子波主频为 2 5 0 Hz 。 2 ． 2正 演模 拟超 前探 波场特征 图 4是正演模拟的波场快照。 图 4 a显示 的是 激发点 在煤层 激 发 。 地震 波 向右 传播时， 在陷落柱左界面产生反射波和透射波 。 由于 顶 底板 和煤 层存 在波 速差异 ，故 地震 波在煤层 中等 l盈鲤圃嘲 - ⋯‘ ’⋯ a 1 8 m s 盛魏 艄噶 ．． 时面呈剪 刀状 ， 并且 在顶 底板界 面上 产生转 换 波 。 图 4 b显示 地 震 波 向右传 播 至 右界 面 并 发 生反 射和透射 ， 反射波与转换波在陷落柱 内产生干涉； 地 震波向左传播至检波器位置 ， 由于受到巷道束缚 ． 信 号在其尾部干涉叠加形成槽波。 图4 c 显示陷落柱右侧界面产生的反射波 。 透过 左侧界面。 与左侧界面产生的反射波一起 ， 在巷道内 形成前后两组反射波组。多次波在陷落柱左侧界面 附近发 生干涉 ， 波 场变得 复杂 。 图 4 d显示 陷落柱左右界面产生的反射波在煤 层 内形成两组能量较强的反射波向左传播 。从波场 快照显示 的波的运动学和动力学特点．可以识别复 杂波场中几种主要类型的波 其中 1 为震源子波遇 到顶底板后产生的透射 P 波， 2为反射 P 波遇到顶底 板产生 的反射 P 波 。 3为波震波遇到顶底板后产生 的透 射 S波 ． 4为 反射 S波 遇 到顶 底 板 产 生 的二 次 转换 S波 ， 5为直达槽波 ， 6和 7为巷道陷落柱左界 面的两组反射 S波， 8为多次波和反射槽波。 巷 道 内 由于是点 源激 发 ， 只有 直达 P 波 、 直 达槽 波 、 反 射 P 波 、 反射 S波 以及 反射 槽波 和多 次 波 。在 分量上 ， 只有直达 P 波、 多次反射波 、 直达槽波 、 陷 落柱左侧界面反射 P 波 、陷落柱右侧界面反射 P波 以及反射槽波。 2 ． 3模拟地 震记 录和 处理 结果 观测 系 统参 数 为 检 波器 1个 ， 道 间距 l O re、 炮 ． 擘 ⋯t 黝 _ 蕊 ． 毒 第 b 2 4 m s ； 0 妇 爱 ⋯” 2 ’ ⋯一 ⋯ 目 l 一 -- i 骖 ，⋯ { ⋯ ⋯ ‘ h_ ⋯ ‘ l j c 3 6 ms d 5 2 m s 图 4 陷落 柱 震 波 超前 探 测 模 型 正 演 模 拟 波场 快 照 Fi g ur e 4 Fo r wa r d s i mu l a t i o n wa v e f i e l d s na p s h o t s o f s u b s i d e d c o l u m n s e i s mi c wa v e a dv a nc e d p r o s pe c tin g m o d e l 3 期 王 勃。 等 陷落柱地震 波超前探测数值模拟与应用 5 9 间距 5 m， 采样 时间 2 5 0 m s ， 采样间隔 0 ． 2 ms ， 合成地 震记 录 见 图 5 。 结合波场快照和波的时距 曲线特征 ，可得 图 5 a中①为直达 P 波 ， ②为顶底板多次反射波和直达 a合 成原始记录脉平分量 b预 处理后结果 图 5 数值模拟合成地震记录 Fi g ur e 5 Nu me r i c a l s i m u l a t i o n s y n t h e t i c s e i s m o g r a m s 槽波，③为巷道前方陷落柱左侧界面反射 P波， ④ 为巷道前方陷落柱右侧界面反射 P波 ， ⑤反射槽波 和多次反射、 折射波 。 图 5 b为采用 丁 一 P视速度滤波及切除技术等预 处理后的记录， 去除了直达波 ， 转换波 ， 直达槽波与 多次反射的干扰波。 以 2 m ／ m s的煤层 P波速度对 分量预处理后 信号进行叠前绕射深度偏移 ，并过滤相对弱能量异 常界面 ， 提取强能量反射界面 ， 结果见图 6所示 。在 方向上 1 1 0 m和 1 2 8 m处f 巷道最外侧检波器所在 位置为坐标原点 有两个强反射界面， 即为陷落柱前 后两个界面 ， 与模型中陷落柱位置吻合。 3 应用实例 新元井田位于沁水煤 田西北部 ，行政区划属山 西省晋中市寿阳县。 井 田内断层稀少 ， 没有岩浆岩侵 人的影响， 构造条件分类应属简单类略偏 中等。 现阶 l l 7厂 ／ | 图 6 超前探模型纵波深度偏移图 Fi g ur e 6 Co m p r e s s i v e wa v e d e p t h m i g r a t i o n g r a p h o f a d v a n c e d pr o s p e c t i ng m o d e l 段采区所采煤层为 3 煤 ， 3 煤为中灰 、特低硫 、 低 磷、 极易选的贫煤、 贫瘦煤 ， 属于山西组， 瓦斯含量较 高。其顶板为砂质泥岩 、 泥岩 ， 局部为中、 细粒砂岩 ， 底板为砂质泥岩， 局部为细、 粉砂岩。 2 0 1 0年 7月 1 3日对 新元矿 南二正巷 进前探 搴 暑 ● a． 旦 ● 一 - ． 一 日． “ 々 旦 一 一 ■ ● { 蕃 1 l I I } 地 震波超前探测 实测信号肜 量 预 处理结果 ①为陷落柱前界面反射纵波； ②为陷落柱后界面反射纵波预处理包 括道间平衡 、 道 内平衡 、 切除直达波 、 r P视速度滤波 图 7 实 测 信 号 与 处理 Fi g u r e 7 M e a s ur e d s i g n a l s a n d p r o c e s s i ng 测 ， 实际信号如图 7 。 地震探测数据 以 l - 8 m ／ ms r 纵波速度进行偏移 ． 结果显示当 日掘进头位置 2 0 3测点 向东 6 2 ． 6 m 前 1 0 0 m范 围内存在两处主要反射界面 图 8 ， 分别为 掘进头前方 3 0 ～ 3 4 m R 和 8 4 ～ 8 8 m R 。结合地质 资料 ， 预测 R 界面可能为陷落柱边界 ， 而 尺 界面可 0n ．． I6 0 ．． ． I钟 ． ．。 1 2 0 0 t 5 0 0 Dism l 、 0 ． _ - 一 、、 I 一 _ 邕 1 - I - I 。。 』 _ ， l ● 一 笤 0。 i 图 8 新元矿南二正巷超前探测纵波深度偏移 图 Fi g ur e 8 Co mpr e ssi v e wa v e de pt h m i g r a t i o n gr ap h o f a d v a nc e d p r os pe c t i n g i n S 2 ma i n r o a d wa y ． Xi n y ua n c o a i mi n e 能为陷落柱另一边界。 实际验证情况为 对应揭露一个陷落柱东边 界 ，岩 性显 著 变 化 ，位置 为 2 8 m； R 为 陷落 柱 西边 界 ， 位置为 8 4 m。 陷落柱呈椭圆形 ， 长轴 5 6 m， 陷落柱 内岩性杂乱 ， 棱角 明显 ， 以砂岩 、 砂质泥岩碎块充填 为主， 柱状边缘呈上大下小的倒八字形特征 ， 无水。 从对 比结果分析得出，地震波超前探测技术在 新元矿的应用较好 ， 探测陷落柱位置准确 ， 表明地震 波超前探测技术对巷道安全掘进起到很好 的指导作 用 ． 可 以有效 的指导 巷道 施工 。 下 转 第 6 5页 3期 王永明 利用测井 曲线进行煤 、 岩层追综对 比的方法及应用 6 5 出垂 向沉 积序 列 中不 同岩 性 的地层 ．而且 还 可 以有 效地 反 映 出上 述成 因标 志在垂 向上 的变化 。通 过全 区对比．不但可以直观的反映上述成因标志在空间 上 的变 化 特征 ． 而 且 在研 究地 质 构 造 、 水 工 环地 质 、 快速煤质分析 、岩石的弹性力学性质以及指导无芯 钻探等方面也十分重要。 参 考文 献 【 1 ]黄作华 ， 叶庆生． 煤 田测井综 合解 释【 M] ． 北京 煤炭工 业出版 社 ， 1 9 8 0 ． [ 2 ] 王惠廉． 综合地球物理测井[ M] ． 北京 高等教育 出版社 ， 1 9 8 4 ． 【 3 ] 尉 中良， 邹长春． 地球物理测井【 M】 ． 北京 高等教育 出版社 ， 2 0 0 4 ． 黄 智辉 ． 陈曜岑． 煤 田测 井资料在分 析沉积环境 及研究 聚煤带 中 的应用[ J ] ． 物探与化探 ， 1 9 8 5 ， l 9 ． [ 5 ]苗争鸣． 利用测井曲线物性特征确定鸳鸯湖矿 区石槽村井 田煤层 对 比的依据[ J ] ． 中国煤炭地质 ， 2 0 1 0 ， 2 2 1 ． 【 6 ] 罗忠文． 金沙龙凤煤矿测井 曲线特征及其在煤岩层对 比中的应 用 叫． 中国煤炭地质 ， 2 0 0 9 ， 2 1 5 ． [ 7 】胡菊艳 ， 郭鸿 ．测 井曲线在许疃井 田煤岩层对 比中的应用Ⅲ． 中国 煤炭地质 ， 2 0 0 8 ， 2 0 1 2 ． [ 8 ]余德志． 贵州 省晚三叠世含煤地层 特征及对 比I J 】 ． 中国煤炭地质 ， 2 0 0 9 ， 2 1 5 ． 【 9 ] 伍远学， 陆玉欢． 测井 曲线在含煤岩系的划分及对 比中的应 用一 以 贵州省盘县北部某井 田为例[ J ] ． 中国煤炭地质 ， 2 0 0 9 ， 2 1 5 ． [ 1 O ]韩绪 山， 谢波 ， 张心彬． 贵州省 金沙煤 田龙 潭组岩煤 层测井对 比 方法[ J ] ． 中国煤 田地质 ， 2 0 0 6 ， 1 8 3 ． 上 接第 5 5页 经矿业公 司依据地震 资料布设的 3个钻孔验证， 煤 层底板解释深度相对误差范 围为 1 ． 2 3 ％ 1 ． 4 6 ％． 符 合相关规范要求。对无煤 区边界所布设的一个验证 孔与解释成果也吻合良好 。 在地震条件极差的巨厚砾石覆盖 区，采用大吨 位可控震源可以取得高品质的地震资料 ，而中小型 可控震源由于多种因素限制难以达到勘探要求 。大 吨位可控震源能够很好地解决 山前戈壁等复杂工区 的激发问题，今后在此类困难地 区开展地震勘探时 可作借鉴 。 参考 文献 ⋯ 1 陶知非． 世纪之交论可控震源的发展与变化【 J 】 ． 物探装备 ， 2 0 0 0， 1 0 1 1 - 6 ． [ 2 】陈祖斌 ， 林君． 可控震源扫描控制技术研 究[ J ] ． 计算机工程与应用 ， 2 0 0 2 ， 3 8 7 2 4 1 2 4 3 ． [ 3 ] 张德 忠． 复杂地 表地 区地震 勘探实例【 M] ． 北 京 石油 工业出版社， 1 99 4 ． [ 4 】尹金章． 可控震 源技术在表 、 浅 层地震地质 条件极 复杂区的应用 [ J ] ． 中国煤 田地质 ， 1 9 9 9， 1 1 S 1 ． [ 5 ]张新生． 浅谈可控震源在煤 田地震勘探 中应用的几点认识[ J 】 ． 中国 煤 田地质 ， 2 0 0 6 ， 1 8 S 1 ． 【 6 】冯义德 ， 李建设 ， 可控 震源在古汉 山厚层砾石 区的地震勘探效 果 l J 】 ． 中国煤 田地质 ， 1 9 9 7， 9 3 ． 上 接第 5 9页 5 结论 地震波探测技术具有探测能力强、 效率高 、 成本 低等优点， 越来越多地被用于井下陷落柱的探测。 通 过数值模拟和现场探测应用对地震波超前探测陷落 柱 方法 可获得 如下 结论 。 ①巷道前方陷落柱的地震波响应在有限测线区 间上 ， 反射波同相轴为负视速度的特征。 ②地震波超前探测方法对于巷道陷落柱超前探 测具 有针 对性 和有效 性 。 ③ 现阶段 ， 地震波超前探测陷落柱 ， 主要是确定 陷落柱前后两个反射界面。对于陷落柱定性问题还 有待进一步探索 。 参考文献 【 l 】 刘 重举．陷落柱 的发育规 律及探 测技 术 『 J 1 ．煤 田地质 与勘探 ， 1 9 9 7 ，2 5 5 ． [ 2 】 胡运兵， 张宏敏， 吴燕清． 井 下多波多分 量地震反射 法探测 陷落柱 的应用研 究[ J l _ 矿业安全与环保， 2 0 0 8 ， 3 5 3 ． [ 3 ]Wa n g B o ，L i u S h e n g d o n g ，L i u J i n g ，e t c ．S t u d y o n a d v a n c e d p r e di c t i o n f o r mu l t i p l e d i s a s t e r s o u r c e s o f l a n e wa y u n d e r c o mp l i c a t e d g e o l o g i c a l c o n d i t i o n s [ J ] ． Mi n i n g S c i e n c e a n d T e c h n o l o g y ， 2 0 1 l ， 2 1 5 ． [ 4 】刘树才， 岳建 华， 刘志新． 煤 矿水 文物探技术与应用[ M] ． 江苏 徐州 中国矿业 大学 出版社． 2 0 0 5 ． 【 5 ]董守华， 王琦 ． 层析成像 在巷道无 线电波透视 法中的应用 【 J I． 中国 矿业大学学报 ， 2 0 0 3 。 3 2 5 5 7 9 ． 【 6 】 刘 志新． 无线 电波透视法在探测隐伏 导含水构造中的应用【 J 1 ． 工程 地球物理学报 ． 2 0 0 5 ， 2 f 3 1 1 8 5 ． [ 7 】 刘 盛东， 郭立全 ， 张平松． 巷道前方地 质构造 MS P法超前探 测技术 与应用研究【 J J ． 工程地球物理学报 ， 2 0 0 6 ， 6 3 4 2 8 4 4 2 ． [ 8 ] 张平松 ， 刘盛东 ， 曹煜． 坑道掘进立体电法超前预报技术研究 中 国煤炭地质 ， 2 0 0 9 ． 2 1 2 ． 【 9 ]张纪勋 ， 范志平． 瞬变电磁技术 在巷道超前探测 中的应用f J ] ． 中国 煤 田地质 ， 2 0 0 6 ， 1 8 3 ． [ 1 O ]刘盛 东， 郭立全． MS P技术及其在煤矿巷道小构造 探测 中的应 用 【 J J ． 安徽理工大学学报 自然科 学版 ， 2 0 0 7 ， 9 增刊 2 2 2 5 ． [ 1 1 】张平 松， 刘盛东， 吴健生． 隧道及井巷工程超前探测模拟及其偏移 技术研究f J ] ． 岩石力学与工程学报 ， 2 0 0 7 ， 2 6 增 1 2 0 4 8 2 0 5 1 ． [ 1 2 ]赵晶． MS P技术在任楼煤矿巷 道超前探测 中的应 用『J 1 ． 中国煤 炭 地质 ， 2 0 1 1 ， 2 3 1 0 ． [ 1 3 】刘忠远． 矿 井地震超前探测技 术在龙东矿 的应用 f J 1 ． 中国煤炭地 质 。 2 0 1 1 ． 2 3 1 ． [ 1 4 ]朱光 明， 李桂 花 ， 程建远． 煤矿 巷道内地震勘探 的数值模拟[ J ] _ 煤 炭学报． 2 0 0 8 ， 3 3 1 1 1 2 6 3 1 2 6 7 ．