扇形探测技术在超前探测中的应用研究.pdf

第 3 6卷 第 6期 2 0 0 7年 1 1月 中国矿业 大 学学报 J o u r n a l o f C h i n a Un i v e r s i t y o f M i n i n g Te c h n o l o g y Vo 1 ． 3 6 NO ． 6 NOV ． 2 00 7 文章 编 号 1 0 0 0 1 9 6 4 2 0 0 7 0 6 0 8 2 2 0 4 扇形探测技术在超前探测 中的应用研究 刘 志新 ，岳建 华，刘仰光 中国矿业 大学 资源 与地球科 学学 院 ， 江苏 徐州 2 2 1 1 1 6 摘要用扇形探测技术代替 了常规 的测量方法 以提 高矿 井瞬 变电磁法超前探 测资料 的解释精 度． 在数据 采 集过程 中， 调整 巷道 迎 头左侧 拐 角处 天 线的 法 线方 向 与左侧 巷 道 壁 的 夹 角 ， 以 6 0 。 ， 4 5 。 和 3 0 。 的 夹角进行 探测 ， 在 迎 头断面 处布置 2 ～ 3个 测 点 ， 调整 巷 道 迎 头右 侧 拐 角 处天 线 的 法 线 方 向与右侧巷 道壁 以 3 0 。 ， 4 5 。 和 6 0 。 的 夹 角进行 探 测． 通 过 调 整 天 线 的 法 线 与巷 道 底 板 的 夹 角， 探测巷道顶板、 顺层和巷道底板方向上介质的电性变化情况． 研 究结果表 明 在超前探测 中使 用扇形探测技术可提 高巷道超前探测资料的解释精度 ， 能更准确的反映出异常体的空间特征． 在 实际应 用 中探 测成果 图上 钻孔 中心 点位 置与 实 际揭 露钻 孔 中心 点位 置误差 不到 2 m． 关 键词 矿 井瞬 变 电磁 法 ；扇形 探测技 术 ； 超 前探 测 中图分类 号 P 6 3 l 文献 标识 码 A App l i c a t i o n o f S e c t o r De t e c t i o n Te c h no l o g y i n Ad v a n c e d De t e c t i o n LI U Z h i x i n，YUE J i a n h u a，LI U Ya n g g u a n g S c h o o l o f Re s o u r c e s a n d Ea r t h S c i e n c e ，Ch i n a Un i v e r s i t y o f Mi n i n g Te c h n o l o g y，Xu z h o u ，J i a n g s u 2 2 1 1 1 6 ，Ch i n a Ab s t r a c t S e c t o r d e t e c t i on t e c hn o l og y wa s u s e d t o e nh a nc e t h e e xp l a na t i o n a c c ur a c y of mi ne t r a ns i e nt e l e c t r o m a gn e t i c m e t ho d i n a d v a nc e d d e t e c t i on i n s t e a d o f c o n ve nt i on a l s u r ve y ． Du r i ng da t a a c q u i s i t i o n a n d d e t e c t i o n，t he a ng l e b e t we e n t h e n or ma l d i r e c t i on o f t he a nt e nn a a nd t he l e f t r o a d wa y i s a d j u s t e d t o 6 0 。 ， 4 5 。a n d 3 0 。 r e s p e c t i v e l y a t t h e h e a d o n l e f t c o r n e r o f r o a d wa y ， t wo or t h r e e s ur v e y po i n t s we r e l a i d ou t i n f r on t o f t he e xc a v a t e d s e c t i o n，a nd t he n t he a ng l e a t t h e h e a d o n r i g h t c o r n e r o f r o a d wa y i s a d j u s t e d t o t h e s a me d e g r e e t o t h e l e f t s i d e ．Th e v a r i a t i o n o f e l e c t r i c p r o pe r t i e s of m e di u m i n t he r oa dwa y r o of ，t h e be dd i ng a nd t he r o a d wa y f l o o r wa s d e t e c t e d b y a d j 『 u s t i n g t h e a n g l e b e t we e n t h e n o r ma l d i r e c t i o n o f t h e a n t e n n a a n d t h e r o a d wa y f l o o r ． The r e s ul t s s h ow t ha t t he da t a i nt e r p r e t a t i o n a c c ur a c y c a n be i mpr o v e d，a n d t he s pa t i a l f e a t ur e of a n oma l y bo d y c a n be r e f l e c t e d mor e e x ac t l y us i n g t h e s e c t or d e t e c t i on t e c hn o l _ 。 o g y i n a d v a n c e d d e t e c t i o n ．Th e p o s i t i o n e r r o r b e t we e n t h e b o r e c e n t e r o f d e t e c t i o n p l o t a n d t h a t o f d i s c l o s e d pl o t i s l e s s t ha n 2 m e t e r s ． Ke y wo r d s mi ne t r a n s i e n t e l e c t r o ma gn e t i c me t h o d；s e c t or de t e c t i on t e c hn o l og y；a dv a n c e d d e t e c t i on 目前用于巷道超前探测主要有 2种方法， 即直 接法 和 间接法． 直 接方 法 为 钻 探方 法 ， 通 过 对 钻 孔 岩芯变化情况的分析 ， 可确定前方是否存在地质异 常． 钻 探结 果 比较可靠 ， 但施 工周期 较长， 费用较 高， 对巷道的正常掘进生产影响较大， 另外 ， 钻孔探 测也只是一孔之见 ， 不能对前方采动影响范围内的 收稿 日期 2 0 0 7 0 41 0 基金项 目高等学 校博 士学科点专项科研基金项 目 2 0 0 5 0 2 9 0 5 0 1 ； 国家 自然科学基金项 目 4 0 6 7 4 0 7 4 作者简介 刘志新 1 9 7 6 一 ， 男 ， 安徽省萧县人 ， 讲师 ， 硕士 ， 从事地球探测与信息技术 方面的研究。 E - ma i l l i u z h x c u ret ． e d u ． c n T e l ； 1 3 8 5 2 1 5 0 8 3 8 维普资讯 第 6 期 刘志新等 扇形探测技术在超前探测 中的应 用研 究 8 2 3 情 况给 予评 价． 超 前探 测 的间接方 法 即采 用物 探方 法进行探测 ， 可用 于超前探测 的物探方法主要有 矿井 直流 电法 三点 三极 超前 探测 方法 、 矿 井 瞬变 电磁法、 瑞雷波法和矿井地质 雷达探测法． 其 中瑞 雷波 法 主要解决 地质 构造 界 面 的问题 ， 对构 造 的富 水性无 法进 行解 释 ； 矿井 地质 雷达 现在 主要 处 于研 究试 用 阶段 ， 主要 是 由 于其 探 测 深 度 较小 ， 不 能 满 足实 际生产 的需 要． 现 在 常用 于超 前探测 的物探 方 法为 矿井 直流 电法 三点 三极 超前 探测方 法 和矿 井 瞬变 电磁 法l_ 】 ] ． 1 三点 三极超前探 测方法技术 三 点 三极 超 前 探 测方 法 是 利 用矿 井 直 流 电法 工作 原 理 ， 施 工过程 中在 巷 道掘进 头 以一定 间隔布 置 A1 ， A2 ， A3这 3个供电电极 即 3点 ， 另一供电 电极 B在无穷远处， 测量 电极 MN 在巷 道内按箭 头所 示方 向以一定 间隔 移动 如 图 1 所 示 ． B 。 ． ．． 迎 ’ Ⅳ l I J 头 图 1 三点三极超前探测示意 图 Fi g．1 Ske t c h ma p o f a dv a nc e d de t e c t i o n us i n g t hr e e dot s t h r e e el e c t r o de me t ho d 在施工过程 中一次探测 3 组视 电阻率值 ， 目的 是 利用 同一 组 MN 测量 的 3组 视 电 阻率 值 进 行校 正 ， 消 除干扰 ， 提 高解 释 准 确 度 ． 测 量 电极 MN 的 间距 根 据地质 任 务和勘 探 的详 细程度 而定 ， 同时也 要 考虑信 噪 比 的大小． 在资料 处 理解 释时 主要 是利 用 曲线 比较法 确定 异常 点 ， 利 用几 何作 图法 确 定地 质异 常 的具体 位置 ． 三点三极超前探测方法技术和钻孑 L 探测相类 似 ， 只能反 映 出巷道 迎头 某一 方 向上 的地层 界 面异 常是否存在， 而不能获得巷道迎头前方在采动影响 范 围内 的构 造分 布情 况 ， 即不 能探 测 出巷 道 迎头 前 方全 方位 的 电性 变化 情况 ． 2 矿井瞬变 电磁法超前探测技术 2 ． 1 矿 井瞬 变 电磁 法原 理 在导电率为 、 导磁率为 ／1的巷道迎头表面放 置面积为 S的矩形发射天线 ， 并在天线 中供 以阶 跃脉冲电流 ， 一 {。I ㈩ 在电流断开之前 0时 ， 发射 电流在回线周 围 的大地和空间中建立起一个稳定 的磁场 ； t 一 0时 刻 ， 将 电流突然断开 ， 由该 电流产生的磁场也立即 消失． 一次磁场的这一剧烈变化通过介质传至发射 天线附近的巷道 围岩 当中， 并在围岩中激发出感应 电流以维持发射 电流断开之前存在 的磁场、 使空间 的磁场不会即刻消失． 由于介质的欧姆损耗 ， 这一 感 应 电流将 迅速 衰减 ， 由它产 生 的磁 场也 随之 迅速 衰 减 ， 这种 迅速衰 减 的磁 场又 在其周 围的介质 中感 应出新的强度更弱的涡流场 二次场 ． 这一过程继 续下去 ， 直至地层 的欧姆损耗将磁场能量消耗完毕 为止 如图 2所示 ， 伴随这一过程存在的电磁场便 是矿 井 瞬变 电磁 场． 图 2 T E M 探测方 向示 意图 Fi g． 2 Ske t c h ma p o f de t e c t i ng di r e c t i o n 由于矿井瞬变 电磁超前探测是在地下几百米 深 度 的井下 巷道 内进 行 ， 其 瞬变 电磁测量 装 置接 收 回线 中 的感 应 电位 为 巷 道周 围空 间有 效 探 测 范 围 内所有岩层导电性 的综合响应． 因此， 矿井瞬变电 磁视电阻率为全空间岩层导电性的综合反映 ， 视电 阻 率计 算公 式 。 C X ／A 0 2／3 一 C 6 ． 3 2 1 O - l Z S N ／ 。 V／ I -- 2 ／ 3 ～ ／ 3 ， 2 式 中 C为 全空 间 响 应 系数 ； S为 接 收 回线 线 圈 面 积 ； N 为 线 圈匝数 ； t 为二 次 场衰减 时 间 ； v／ i 为 接 收的归一化二次场电位场值． 在巷道 中迎头 附近 ， 通过 接收 线 圈测 量 二次 涡 流场的感应 电动势随时间的变化特征 ， 就可 以了解 巷 道迎 头 附近介 质 的 电磁 性质 变化 规律 ， 进 而确 定 电磁性质异常体 的规模、 产状和到巷道迎头的距离 等 ． 2 ． 2 超 前探 测施 工方 法技 术 由于受巷 道迎 头 空 间的限制 ， 矿 井 瞬变 电磁 法 的发 射 和接 收线 圈 的几 何 尺寸受 到 的一定 的制 约 ， 只能采用多匝小 回线的发射和接收装置形式 ， 即边 长为 2 ～3 m． 测点布置在巷道迎头空间位置 见如 图 3 ， 即从巷道迎头左侧开始 ， 首先使 发射、 接收 天线 的法 线垂 直巷 道左 侧 面进 行 测量 图 3中 的 1 号测点 ， 然后旋转天线 ， 使天线的法线方向与巷道 的左 侧分 别成 6 O 。 ， 4 5 。 和 3 O 。 的夹 角进 行 探测 图 3 中的 2 ， 3和 4号测点 ； 当天线 的法线方 向与巷道 维普资讯 8 2 4 中 国 矿 业 大 学 学 报 第 3 6卷 迎头界面垂直时 ， 根据其主迎头断面的宽度布置 2 ～ 3个测 点 图 3中的 5 ， 6和 7号测 点 ； 到巷 道 迎 头右侧时再旋转天线 ， 使法线方向与巷道右侧分别 成 3 O 。 ， 4 5 。 ， 6 O 。 和 9 O 。 的夹角 进行探 测 图 3中 的 8 ， 9 ， 1 0和 1 1 号测点 ． 即在多个角度采集数据， 从而 获得尽 可能 完整 的前 方空 间信 息 ， 故称 之 为扇形 测 深技术． 图 3 超前探测测点 布置示意 图 F i g ． 3 S k e t c h ma p o f a d v a n c e d d e t e c t i o n 在实际工作过程中对于每个发射点， 也可调整 天线 的法线 与巷道 底板 的夹 角大小 ， 以探测 巷道 顶 板 、 顺层 和底 板方 向 的围岩变 化情 况[ 8 ] ． 2 ． 3 资料处理 及成 果显 示技术 瞬变电磁法 的资料解释步骤是首先对采集到 的数据进行去噪处理， 根据晚期场或全期场公式 2 将 仪器 测量 的 电流归 一 化 值 转换 成不 同时 间 窗 口 对应的视电阻率值， 然后进行 时深转换处理， 得到 各测线视电阻率断面图． 在成果显示上主要有 2种 方法 第 1 种 为常用 的矩形成 图方法u “ ] ， 将 图 3 中的各测点视为在一条直测线上等间距分布 的点 ， 绘制的视电阻率等值线图中横坐标为测点标号 或 为相对 I 号 测点 的相 对 距 离 ， 纵 坐标 为 相 对 直 测 线 的探测深度； 第 2种为本文重点介绍的扇形成图 方法 ， 将 图 3中各 测点 测量 的不 同深度 上 的视 电阻 率值分布在实际平面位置， 将视电阻率等值线图绘 制成扇形方式， 其纵坐标为在迎头断面向正前方探 测 的实际距离、 横坐标对应 图 3中 I号和 I I号测 点 的探 测距 离加上 巷道 断面 的宽 度 ， 其 中巷道 的断 面宽度可根据正前方测点 的间隔进行等 比例的放 大． 3 实例 分析 煤矿资源普查与勘探初期 ， 在井田范围内施工 了许多勘探及水文观测钻孔 ， 其中有许多钻孔穿过 了煤层顶板或底板的含水层． 现在 ， 有些钻孔的准 确位置及封孔资料很难收集 ， 这给煤矿的安全生产 及地下水资源的合理开发利用带来了潜在的威胁． 由地 面已知地质资料可知， 某煤矿巷道掘进迎头附 近存在一个勘探钻孔， 该钻孔在煤层底板下的灰岩 中终止． 但是钻孔的孔斜及封孔资料不详 ， 如果该 钻孔封孔效果达不到要求， 并且在煤层附近的具体 位置不能准确确定 ， 当与煤层底板下的含水层导通 时 ， 在巷道掘进及工作面 回采的过程 中， 由于受采 动的影响或直接对钻孔 的揭露， 将会诱发涌水事故 的发生 ， 对煤 矿 的安全 生产 带来 较为严 重 的影 响． 在该 钻孔 附 近只 有一条 正在 掘进 的独头 巷道 ， 因此采用矿井瞬变 电磁法超前探测技术按图 3中 的测点布置方式进 行探测， 在探测过程调整 了发 射 、 接收天线与巷道底板 的夹角， 完成了向巷道顶 板 、 顺层 和 巷道底 板 的探测 任务 ． 查 明 了巷 道 迎头 正 前方 I 1 0 m 以 内、 前 方底 板 6 0 m 以上 、 顶板 6 O m 以下 一个 锥 形体 积 范 围 内地层 介质 的视 电阻 率 变 化情 况． 图 4为采用矩形成 图方法绘制 的视电阻率等 值线 图 图 中上面横 坐标 为在 巷道 内相 对于 I 号 测 点的相对距离 m ， 即先设计 各测 点 的间距 为 5 m； 下 面 横 坐标 为 巷 道 内布 置 的测 点 标 号 ， 纵 坐标 为 沿探测 方 向 的探 测距 离 m ． 实 际上 巷道 迎 头 断 面的宽度 3 ～6 m， 标 号为 2 ～1 O的距离最大为 6 m， 因此图 4中只提供 了相对的位置 ， 这给资料的 解释精度带来 了一定的误差． 在数据采集过程 中， 巷 道迎头 附 近标号 为 I ， 2 ， 1 0和 I I的 4个 测 点 位 置 已经利用工字钢或金属锚网支护好 ， 而金属体对 矿井瞬变 电磁法探 测 的数据会 产生较大干扰． 因 此 ， 这些测点 的资料对探测结果的解释又产生不利 的影响， 一般作为背景值进行处理． 相对距离／ m 姜 斟 赚 测点标号 图 4 超前探测视 电阻率断面 Fi g． 4 I s ol i ne s o f a ppa r e nt r e s i st i v i t y o f a d v a nc e d de t e c t i o n 从 图 4中可 以看 出在 巷道 迎头 前方 约 5 0 m 处 有一个视 电阻率值小 于 7 Q m 的狭长状封闭区 域 ， 可能是钻孔 的反映． 由于图 4中的横坐标是采 用相对的距离， 不能准确确定视电阻率值小于 7 Q 维普资讯 第 6期 刘 志新 等 扇形探测技术在超前探测 中的应用研究 8 2 5 m封闭区域在横 向上 的位置． 另外 ， 在两侧 由于 支护金属体的影响也 出现 了低阻区域． 图 5为 采 用 扇形 成 图方法 绘 制 的 视 电 阻率 等 值线图 该图中只利用图 3中的 2 ～1 O测点 即巷 道迎头 探测的数据， 将每个测点探测 的数据按实 际位置坐标绘在平面 图上 ， 横坐标为以巷道中点为 中心 点 向两侧 的实 际宽度 ， 纵 坐标 为相对 于巷 道 迎 头 的实际探测距离． 另外由于巷道迎头断面的宽度 一 般为 3 ～6 m， 为 了有 效利 用 5 ～7号 这 3个 测 点 的数据 ， 将巷道的宽度按一定 比例进行加宽． 这样 可准确计算出图中任 意一点到巷道迎头的具体距 离 ． 从 图 5中可 知 ， 在 巷 道 迎 头 正 前 方 5 5 m， 左 侧 3 0 m有一 个 圆形 的异 常 区域 ， 在视 电 阻率 值为 7 Q m的等值线 内又显示 出一个 更 低 的视 电 阻率 值即 6 ． 3 Q m对应等值线． 该异常区解释为钻孔 的 反映 ， 其 位 置与 地 面 钻孔 的位 置 基 本 一 致 ， 说 明 该钻 孔 的孔斜 度较 小． 根 据矿 井 瞬变 电磁法 在该 矿 区探测的经 验可知 ， 当视 电阻率 值小 于 4 Q m 时 ， 其富水性较强． 因此 ， 可判断该钻孔在煤层附近 富水 性较 弱 ， 没有 与 深 部 的 灰 岩 水 直 接 导 通 的 通 道 ， 钻 孔封 孔效 果较 好 ． 但 是 由 于受 钻 孔 对 地 层 破 坏的影响 ， 周围地层 的视 电阻率值呈低阻反 映， 随 着到钻 孔距 离 的 增 大 ， 地 层 的视 电 阻 率 值 逐 渐 稳 定． 该结果在煤层 回采过程 中被验证 ， 钻孔的实际 揭露位置位于当时探测迎头 正前方 5 3 ～5 4 m， 左 侧 3 1 ～3 2 m 范 围内 ， 钻 孔 中心 点 位 置与 探 测 成果 图 中的中心 点位 置误差 不 到 2 m_ E 避 宽度／ m 图 5 超前探测视 电阻率扇形等值线 Fi g． 5 Se c t o r i s o l i ne s of a p pa r e nt r e s i s t i v i t y of a d va n c e d d e t e c t i o n 另 外 ， 在 向巷 道迎 头前 上 方 顶 板方 向 和前 下 方 底板 方 向 探 测 的视 电阻率 等值 线 图 中， 在 与 图 5中同一位置存在一个类似椭圆形 的封闭区域 ， 视 电阻率值也为 7 Q 1T I ， 为钻孔 的反映． 对比图 4和 图 5可 以看 出 ， 图 5中钻孔 对应 的异 常 区域较 为 明 显 ， 与异常构造的实 际形状也较为一致 ； 图 4中的 异 常 区域 由于采 用 矩 形绘 图方 法 使 异 常 区域 产 生 变形 ， 给资料解释带来 了一定的难度 ， 另外图中采 用 相对 距 离 ， 给 异常 区域 位置 的确 定带来 较 大 的误 差 ． 4 结 论 三 点 三极 超前 探 测技 术 和矿 井 瞬变 电磁法 扇 形超前探测技术是 目前用 于巷道超前探测 的 2种 间接方法． 在扇形超前探测过程中， 首先调整巷道 迎 头左侧 拐 角处 天线 的法 线 方 向与 左 侧 巷道 壁 的 夹 角 ， 以 6 O 。 ， 4 5 。 和 3 O 。 的 夹 角进 行探 测 ； 然后 在 迎 头 断面处 布 置 2 ～ 3个 测 点 ； 最 后 调 整 巷道 迎 头 右 侧 拐角处 天线 的法 线 方 向与右侧 巷道 壁 的夹角 ， 以 3 O 。 ， 4 5 。 和 6 O 。 的夹角 进行 探测 ， 即完 成 一 次扇 形 超 前探 测． 另 外 ， 通过 调 整 天 线 的 法线 与巷 道 底 板 的 夹角 ， 可探测巷道顶板、 顺层和巷道底板方 向上介 质 的 电性 变化 情况 ， 这样 可 获得巷 道迎 头前 方 个 锥形数据体 ， 能够准确反映出探测范围内异常体的 空 间位 置 ． 应用 实 例证 明 在 超前 探 测 中使 用 扇 形 探测技术可提高巷道超前探测资料 的解释精度， 能 更准确的反映 出异常体 的空间特征． 实际应用表 明 ， 探 测成果 图中钻 孔 中心点 位置 与实 际揭露 钻孔 中心点 位置 误差 不 到 2 m． 参 考文 献 1- 1 ] 王齐 仁．地下地质 灾害 地球 物理 探测 研究 进展 [ J ] ． 地球 物理学进 展 ， 2 0 0 4 ， 1 9 3 4 9 7 5 0 3 ． W ANG Qi r e n ．P r o g r e s s i n t h e g e o p h y s i c a l e x p l o r a t i o n r e s e a r c h o n u n d e r g r o u n d g e o l o g i c a l h a z a r d s [ J ] ． Pr o g r e s s I n Ge o p h y s i c s ，2 0 0 4， 1 9 3 4 9 7 5 0 3 ． [ 2 ] 岳建华 ，刘树才 ， 刘 志新 ．巷道直流 电测深在探测陷 落柱 中的应用 [ J ] ． 中国矿业大 学学报 ，2 0 0 3 ， 3 2 5 4 79 481 ． YUE J J a n h u a ，L I U S h u - c a i ，LI U Z h i x i n ．Ap p l i c a t i on o f r oa d wa y DC e l e c t r i c a l s oun di n g i n de t e c t i n g c o l l a p s e c o l u mn s [ J ] ．J o u r n a l o f C h i n a Un i v e r s i t y o f Mi n i n g 8 L Te c h n o l o g y，2 0 0 3 ， 3 2 5 4 7 9 4 8 1 ． [ 3 ] 刘树才 ， 岳 建华 ，刘 江．西部 保水开采中的水文电 法勘探技术 [ J ] ．中国矿 业大 学 学报 ，2 0 0 4 ， 3 3 2 1 87 18 9． L I U S h u c a i ，YUE J i a n - h u a ，LI U J i a n g ．El e c t r i c a l m e t h od u s e d i n h ydr o l o gi c a l pr os pe c t i n g f o r wa t e r - k e e p i n g ～ mi n i n g i n w e s t e r n C h i n a E J ] ．J o u r n a l o f C h i n 3 Un i v e r s i t y o f Mi n i n g Te c h n o l o g y ， 2 0 0 4， 3 3 2 1 8 7 1 89 ． 下转 第 8 6 8页 维普资讯 8 6 8 中 国 矿 业 大 学 学 报 第 3 6卷 [ 4 ] [ 5 ] [ 6 ] 叶彦谦 ，王现．电子注 聚焦 理论 中所 出现 的非线 性方程 [ J ] ．应用数学学报 ，1 9 7 8 ， 1 1 1 3 4 1 ． YE Ya n q i a n，W ANG Xi a n．No nl i ne a r di f f e r e n t i al e q u a t i o n s i n e l e c t r o n b e a m f o c u s i n g t h e o r y[ J ] ．A c t a M a t h e ma t i c a e Ap pl i c a t a e Sc i ni c a，1 9 78， 1 1 1 3 41 ． 周钦德 ， 王怀 中．电子聚 焦理 论 中的一个 周期 解 问 题 [ J ] ．应用数学学报 ，1 9 8 8 ， 1 1 4 4 3 3 4 4 3 ． Z H0U Qi n g d e ，W ANG H u a i z h o n g ．A p e r i o d i c s o l u t i o n p r o b l e m i n e l e c t r o n b e a m f o c u s i n g t h e o r y[ J ] ． Ac t a M a t he ma t i c a e App l i c a t a e Sc i ni c a，1 9 88， 1 1 4 4 33 4 43 ． ZH ANG M e i r on g． Pe r i odi c s ol u t i on o f I i e na r d e q u a t i o n s wi t h s i n g u l a r f o r c e s o f r e p u l s i v e t y p e E J ] ． J o u r n a l o f Ma t h e ma t i c a l An a l y s i s a n d Ap p l i c a t i o n s ， 1 99 6 2 03 25 4 26 9． [ 7 2 [ 8 ] [ 9 ] [ 1 O ] T0RRES P J ．Ex i s t e n c e a n d u n i q u e n e s s o f e l l i p t i c s ol ut i o ns o f t he Br i l l o ui n e l e c t r on be a m f oc u s i ng sy s t e m [ J ] ． Ma t h e ma t i c a l Me t h o d s i n Th e Ap p l i e d S c i e nc e s ，2 00 0 23 1 1 39 - 1 1 4 3 ． ZH ANG M e i yu e，CH EN Ta i y o ng， LI U W e n - bi n． Exi s t e nc e of pos i t i v e pe r i od i c s ol u t i on f or t he e l e c t r o n b e a ms f o c u s i n g s y s t e m[ J ] ．Ma t h e m a t i c a l Me t h o d s i n Th e Appl i e d S c i e n c e s，2 00 5 28 7 7 9 78 8． 韩茂安 ， 顾 圣 士．非 线 性 系统理 论 和方 法 [ M] ．北 京 科 学出版社 ，2 0 0 4 ． LI U Chun - ge n． The s t a bi l i t y o f s u b h ar mo ni c s ol u t i o n s f o r h a mi l t o n i a n s y s t e ms [ J ] ．J o u r n a l o f Ma t h e ma t i c a l An a l y s i s a nd App l i c a t i on s，1 99 9 24 0 49 1 54 0． 责任编辑邓群 上接 第 8 2 5页 [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] 程久龙 ， 王 玉和．巷道 掘进 中 电阻 率 法超 前探 测原 理与应用[ J ] ．煤 田地质 与勘探 ，2 0 0 0 ， 2 8 1 6 0 6 2 ． CHENG J i u l o n g。W ANG Yu - h e ．Th e p r i n c i p l e a n d a pp l i c a t i o n of a dv a nc e s u r ve y i ng i n r o a dway e x c a v a t i o n b y r e s i s t i v i t y me t h o d [ J ] ．C o a l G e o l o g y a n d E x pl or a t i on，2 0 00， 2 81 60 62 ． 刘青雯．井下 电法超前探 测方 法及其应 用[ J ] ．煤 田 地质与勘探 ， 2 0 0 1 ， 2 9 5 6 0 6 2 ． LI U Qi n g we n ． Un d e r g r o u n d e l e c t r i c a l l e a d s u r v e y me t h o d a n d i t s a p p l i c a t i o n[ J ] ．C o a l Ge o l o g y a n d E x pl or a t i on，2 001， 29 5 6 0 62． 李貅．瞬变 电磁测 深的理论 与应用 [ M] ．西安 陕 西科学技术出版社 ， 2 0 0 2 ． 闫 述 ，陈明生 ，傅君眉．瞬变 电磁场的直接进域数 据分析[ J ] ．地球物理学报 ， 2 0 0 2 ， 4 5 2 2 7 5 ～ 2 8 4 ． YAN S h u，CHEN Mi n g s h e n g，F U J u n - me i ．Di r e c t t im e do ma i n nu me r i c a l a na l y s i s of t r a n s i e nt e l e c t r or ma g n e t ic f i e l d s [ J ] ．C h i n e s e J o u r n a l o f G e o p h y s i c s ， 20 02， 45 2 2 75 28 4． 刘志新 ，于景 邮．矿 井瞬变 电磁 法在 水文 钻孔 探测 中的应用[ ． J ] ．物探 与化探 ，2 0 0 6 ， 3 0 1 5 9 ～ 6 1 ． LI U Z h i x i n．YU J i n g c u n Th e a p p l i c a t i o n o f t h e [ 9 ] [ 1 O ] [ 1 1 ] mi n i n g t r a ns i e nt e l e c t r o ma gn e t i c me t ho d t o t he e xp l o r a t i o n o f h y d r o l o g i c a l b o r e h o l e [ J ] ．G e o p h y s i c a l a n d Ge o c h e m ic a l Ex pl or a t i o n，20 0 6， 3 0 1 59 - 61 ． 李貅 ， 郭文波 ．瞬变 电磁法 在煤 田矿 井 涌水通 道 勘察 中的应用 [ J ] ．西安 工程学 院学报 ，2 0 0 0 ， 2 2 3 35 38 ． I I Xi u，GUO W e n bo ． App l i c a t i on TEM t o e x pl or a t i o n o f w a t e r g u s h i n g c h a n n