隐伏煤矿区CSAMT法煤系地层勘查研究.pdf

第3 7 卷第1 期中国矿业大学学报 V 0 1 ．3 7N o ．1 2 0 0 8 年1 月J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g T e c h n o l o g y J a n ．2 0 0 8 隐伏煤矿区C S A M T 法煤系地层勘查研究 安志国1 ，底青云1 ，孑- a 5 韶华2 1 ．中国科学院地质与地球物理研究所，北京 1 0 0 0 2 9 2 ．福建省煤田地质局，福建沙县3 6 5 5 0 0 摘要对于复杂地形山区深部煤系地层的勘探，可控源音频大地电磁法 C S A M T 是比较有效的 一种方法．在福建山区应用C S A M T 法对隐伏矿区煤系地层进行了勘查研究．根据复杂山区数 据采集的特点，提出了针对该地区的地形改正和静态校正的数据处理方法，分析了原始曲线特征 并比对反演结果．结果表明，对有严重静态影响的原始数据采用有阀值的全剖面比较静态校正处 理后，反演模型的响应曲线和原始曲线的拟合误差均小于6 ％．煤系地层在C S A M T 反演电阻率 断面上表现为中低阻，和上覆介质有明显的电性差异．结合研究区地质资料，对C S A M T 反演结 果进行了地质解释，划分了地层并圈定了有希望的煤系范围． 关键词深部隐伏；煤系地层；C S A M T 法；数据处理解释 中图分类号P6 3 1 ．3 2 5文献标识码A文章编号1 0 0 0 - 1 9 6 4 2 0 0 8 0 1 0 1 1 8 0 7 S t u d yo fE x p l o r a t i o nf o rC o a lS t r a t u mw i t hC S A M T M e t h o d i nC o n c e a l e dM i n eA r e a A NZ h i g u 0 1 ，D IO i n g y u n l ，G U OS h a o h u a 2 1 ．I n s t i t u t eo fG e o l o g ya n dG e o p h y s i c s ，C h i n e s eA c a d e m yo fS c i e n c e s ，B e i j i n g1 0 0 0 2 9 ，C h i n a ； 2 ．C o a l f i e l dG e o l o g i cB r a n c ho fF u j i a nP r o v i n c e ，S h aC o u n t y ，F u j i a n3 6 5 5 0 0 ，C h i n a A b s t r a c t C o n t r o l l e ds o u r c ea u d i om a g n e t i ct e l l u r i c C S A M T i sa ne f f e c t i v ea p p r o a c ht od e t e c t d e e pc o a ls t r a t u mi nc o m p l i c a t e d t e r r a i nm o u n t a i n ．C S A M Tm e t h o dw a sa p p l i e df o rc o a ls t r a t u mi nS a n m i n ga r e a ，F u j i a np r o v i n c e ．B a s e do nt h ep r o p e r t yo fd a t aa c q u i s i t i o ni nc o m p l i c a t e d m o u n t a i n ，am e t h o do fd a t ap r o c e s s i n gi n f l u e n c e db yt e r r a i nc o r r e c t i o na n ds t a t i cc o r r e c t i o n w a sp r o p o s e d ，a n do r i g i n a lc u r v e sw e r ea n a l y z e da n dc o m p a r e dw i t hi n v e r t e dd a t a ．A l lp r o f i l e c o m p a r i s o nm e t h o dw a sa d o p t e dt oc a r r yo u ts t a t i cc o r r e c t i o n ，t h ef i t n e s s e r r o ro fr e s p o n s e c u r v e sb e t w e e no r i g i n a ld a t aa n di n v e r t e dd a t ai Sl e s st h a n6 ％．T h er e s u l t ss h o wt h a tt h ee l e c t r i cp r o p e r t yo fc o a ls t r a t u mi sm i d d l e l o wr e s i s t i v i t yi nc r o s ss e c t i o no fi n v e r t e dr e s i s t i v i t y ， a n dc o a ls t r a t u mh a sc l e a rd i f f e r e n c ec o m p a r e dw i t ho v e r b u r d e nm e d i u m ．B a s e do ng e o l o g i ci n f o r m a t i o no fc o a li nr e s e a r c ha r e a ，ag e o l o g i ci n t e r p r e t a t i o nh a sb e e nm a d eo nC S A M Tr e s u l t s ， t h es t r a t u m sa r ec a r v e du p ，a n dt h ep r o s p e c t i v ec o a li sp l o t t e do u t ． K e yw o r d s d e e pc o n c e a l e d ；c o a ls t r a t u m ；C S A M Tm e t h o d ；d a t ap r o c e s s i n ga n di n t e r p r e t a t i o n 近几十年来，福建露天和半覆盖煤矿分布区域 均已勘探、开发完毕，但是根据推覆构造等地质理 论研究，以及近几年的地质调查和相关的钻孔信息 推断，福建省内仍然存在深部的煤系地层．本文介 绍了隐伏煤矿区C S A M T 法煤系地层勘查的研究 成果，属现在倡导较多的“三下找煤”，即在红层下、 收稿晶期2 0 0 7 0 1 2 7 基金项目中国科学院知识创新重要方向项目 K Z C X 2 一Y w 一1 2 1 ；K Z C X 2 一Y w1 1 3 1 作者简介安志国 1 9 7 7 一 ，男，河南省安阳市人，博士研究生，从事电磁法数据处理和应用方面的研究 E r a i l 2 0 0 8 d r e a m f u t u r e 1 6 3 ．C O I T t T e l 0 1 0 8 2 9 9 8 6 2 5 万方数据 第1 期安志国等隐伏煤矿区C S A M T 法煤系地层勘查研究 老地层下和火山岩下找煤． 地震法口3 在地形起伏不大的情况下对煤层作 精细探测应用得较多，瞬变电磁法[ 23 在探测煤矿中 不良地质构造时能发挥很好的作用，重磁法[ 3 3 和7 能谱法[ 43 在煤矿普查、圈定靶区勘探中有着独特的 优势．此次福建煤系地层勘查的地区植被茂密，灌 木丛生，地形复杂，而且探O n , 0 目标体较深，最深达 15 0 0m ，上述地球物理方法难以发挥其优势．而在 此种条件下，C S A M T 法变频测深，探测深度大；横 向分辨率高；测量灵活；比值测量，地形影响较小等 优势凸显．因此，采用C S A M T 法在福建隐伏矿区 进行煤系地层勘查研究．本文首先简单介绍了 C S A M T 法的基本原理，重点陈述复杂山区数据采 集的特点和极端静态数据的处理解释方法，分析原 始资料比对反演结果，进行拟合分析，并结合本区 的地质资料，给出了以大田研究区为代表的 C S A M T 法煤系地层的电性结构特征并圈定了有 希望的煤系范围． 1方法原理 1 9 7 1 年G o l d s t e i n 口] ，1 9 7 5 年G o l d s t e i n 和 S t r a n g w a y [ 6 1 引入C S A M T 法，它是利用人工场源 发射变频电磁信号的一种测深技术．理论上讲，两 端接地的电偶极子向地下供交变电流，感应出反映 地下电性介质的电磁场，C S A M T 法通过测量感应 出来的电场和磁场，根据C a g n i a r d 电阻率公式计 算视电阻率，通过模型和视电阻率曲线的拟合反 演，得到各研究区块的电性深度资料，从而获得对 地电构造的认知．大勘探深度和高分辨率使得它在 寻找地下水睁1 0 ] 、勘查油气Ⅲ] 、矿产资源探测‘1 2 15 1 以及工程地质研究[ 1 6 。1 8 3 中应用颇多，此处不再赘 述． C S A M T 法测量交变电磁场的电场分量E 和 磁场分量H ，电场E 用不极化电极测量电位值获 得，而磁场H 用磁棒在相应地区的响应表示．根据 水平电偶极子理论和M a x w e l l 方程，可以计算E 和H 的3 个分量，用水平电场分量E 。和水平磁场 H 。的比值得出C a g n i a r d 电阻率口9 J ． E 一生掣． 3 C O S 2 0 2 ， 1 L z 一i 。 ＼u一／，＼1 ／ ￡7 【r 式中j 为供电电流强度；A B 为水平电偶极子长 度；0 为A B 和均匀半空问任意点M 的夹角． 一，．一 磐E i2012,S 1 n， 2 ，． ●i o o ，L Z ， 。4 7 r ，．。 式中1 0 。为水平层状模型中第1 层的电阻率． E 。刊川等 式中胁为磁导率；O 是A B 的中点；r 为0 到M 的 距离；0 9 为角频率． H ，一c ，㈨警√急 C O S0 s i n0 ． 4 H ，划㈨警√急 3 c o s 2 0 2 ． 5 H 一i 訾．S i ∽ 6 方程 ∞除以方程 5 ，简化得 舻刚1 倒E j 2 ．∽， 风一产。叫lH | ． 。” 因此，C a g n i a r d 电阻率为 ． 1} E 。l 2 B 一一5 f 舻。 阻抗相位P E 。‰。一H 。h 。． 趋肤深度D 一5 0 2 √手． C S A M T 法通常在电磁波传播的远场区采集 数据，常见的频率系列为2 1 3 4 2H z ，地表电阻率 一定的情况下，频率越低，探测深度越大．电极距的 大小决定C S A M T 法的横向分辨率．本次研究工 作为了满足探测深度的要求，数据采集采用自动发 射接收同步的4 F 3 0 频率系列，测量电极距为4 0 m ，时间同步由G P S 控制． 2 数据处理 野外采集到原始数据后，室内还要对这些数据 进行编辑处理，首先观察原始记录曲线形态，剔除 飞点、跳点，压制噪声，初步给出原始电阻率断面 图，检查评估静态、地形的影响程度，然后对经初步 圆滑处理后的数据进行静校、地改，进而进行反演． 对比反演结果和原始曲线进行拟合分析，若拟合差 大，则说明处理存在问题，需继续之前的静校、地改 和反演，直到拟合差足够满足反演结果的合理性要 求． 进行地质解释之前，还要分析各测深点的视电 阻率一频率曲线，结合研究的目标探测体来说明曲 线反映的地下介质的岩性变化．处理流程见图1 ． 下面就数据处理中涉及到的几个方面做一下介绍． 层 万方数据 1 2 0 中国矿业大学学报第3 7 卷 图1 数据处理流程 F i g ．1 D a t ap r o c e s s i n gf l o w 2 ．1 静态校正 本研究采用E M A P [ z o ] 法低通滤波和全剖面比 较校正处理[ 2 妇来实现所得数据的静态校正。 E M A P 滤波实现了极大限度地提取勘探区域的地 电信息，消除静态和地形影响、提高勘探精度、简化 数据处理方法、提高解释精度和分辨率． 汉宁窗滤波公式为 f 三f 1 c o s 垫1 I zI ≤詈 ， h z 百0 91 ， 、厶， 式中0 , 1 表示窗口宽度． 实际应用中，矗 z 离散成7 点滤波．针对本区 地质条件复杂、地表电性不均匀，还要采取全剖面 数据比较分析校正，结合岩性测试或者综合分析， 得出不受静态偏移的背景电性参数，再进行曲线平 移，选择合理的阀值进行静态校正，亦可取得较好 的校正效果．同时也会参照其他的静校正方法对比 研究以期消除、减小静态效应的影响，得到地下真 实的地电反映．图2 为静态校正前后对比图，其中 图2 a 为静校前，图2 b 为校正后消除了部分静态效 应，条带明显减少． 图2 静态校正前后对比图 F i g ．2C o m p a r i s o no fb e f o r ea n da f t e rs t a t i cc o r r e c t i o n 2 ．2 地形校正 首先对二维纯地形和二维地形下包含异常体 的各种变化条件下地形的影响进行正演模拟，模拟 结果可以直观形象展示出地形对C S A M T 法的视 电阻率影响特征和规律．然后用比值法 畸变张量 法 作C S A M T 地形校正，通过正演计算分别求出 纯地形舭和水平地形均匀半空问的视响应P o 来 求取校正因子，来对实测数据1 0 实进行校正． p 嵌一p 耍fp 她* p 。． 再者由于工区地形复杂，很难做到发射的电磁 波信号和各测点在同一高程平面上，这势必引起各 测点上数据的不一致性．为此需要详细测量各测点 的高程，，在二维地形影响正演计算结果的基础上根 据经验对山体或沟的形态可能的影响做了定性分 析，在反演结果解释时给予了考虑，具体如下 E ， E f C O S0 罢C O S0 一 V z 证V v c 。s 2 臼， 式中E 为水平电场；E 为斜坡方向的电场；V v 为 斜坡得到的电位；V ，为斜坡长度；0 为斜坡倾角；V ， 为斜坡在水平上的投影． 2 ．3曲线分析 在地质解释前，还要分析得到的视电阻率一频 率曲线，分析不含煤系地层和含煤系地层的差别， 结合相关的地质资料，就C S A M T 法反演结果给 出合理的地质解释． 从图3 a 可见，在视电阻率一频率曲线上，高频 段视电阻率值约为1 0 0f ／m ，随着频率降低，视电阻 率逐渐升高，大约在频率2 0 0H z 左右时出现过渡 区，电阻率变化不大，频率进一步变小时，电阻率曲 线又开始上升，这表明在接收频率范围内，测深曲 线没有反映存在低阻层，根据本区标本测线测定结 果，煤系地层呈中低阻显示，则可判断是正常地层 或完整岩层，没有煤系地层存在，该测深曲线是该 区典型的不含煤系地层的反映． 从图3 b 的曲线上可以看出，在视电阻率一频率 曲线上，高频段视电阻率值仍为1 0 0 ～2 0 0Q m ，但 随着频率降低，直到1 0 0 H z 左右视电阻率仍不是 很高，2 0 0 Q m 左右，以后曲线形态逐渐抬升，到 l O H z 左右到达峰值，随后视电阻率又开始变小，测 深曲线反映有低阻响应，说明存在低阻层，根据本 区岩性标本的测试，地球物理推测可能存在煤系地 层，而该视电阻率一测深曲线也是该区含煤系地层 的典型反映的显示． 万方数据 第1 期安志国等隐伏煤矿区C S A M q 、法煤系地层勘查研究 1 0 4 姜1 0 ， 斛 器| 0 2 察 1 0 ‘ 一． 9 0 ￡一相付 垂1 0 3 l 二囊了鬻警阻率蓥‘ ≥嗡、 9 0 基 } 频率／H z 频率／H z a 不含煤系地层 b 含煤系地层 图3 含煤系地层和不含煤系地层的电性响应曲线对比 F i g ．3C o m p a r i s o no fe l e c t r i cp r o p e r t yc u r v e sa b o u tn o n c o a ls t r a t u ma n dc o a ls t r a t u m 2 ．4 拟合分析 在保证采集数据质量良好的前提下，反演后的 数据必须和原始的视电阻率一频率数据进行拟合， 只有好的拟合结果可以对反演结果的可靠性做出 验证，也是给出合理的地质解释前提．得出反演结 果后，笔者对反演结果和原始数据进行了拟合对 比，表1 是含煤系地层的测深点的反演结果及其视 电阻率一频率曲线的拟合结果，从表1 可以看出，拟 合差都小于6 ％；也可从图4 拟合曲线图上看出， 原始视电阻率一频率曲线和反演后的电阻率一频率 曲线非常的接近，两者拟合得较好，较好地反映了 f d 靶 3 研究区地质和岩石电性特征 本次C S A M T 勘查研究工作集中在4 个试验 区块，满足三下找煤的要求，具体分布见图5 ． 各层的岩性变化，也说明了反演结果真实可靠． 2 5 表1原始数据和反演结果拟合表 T a b l e1F i t n e s so fo r i g i n a ld a t aa n di n v e r t e dd a t a 1 0 0 0 得 盛1 0 0 脚 察 1 0 5 一反损化视电阻率。原始视电阻翠 ＼“∥飞、吣⋯。 1 01 0 01 0 0 01 0 0 0 0 频率A I z 图4原始数据和反演数据拟合对比 F i g ．4 C u r v ec o m p a r i s o no fo r i g i n a ld a t aa n di n v e r t e dd a t a ，鹰潭、 i 。f N ，7浙扛省 江西省 ，≮| | 、j ＼ ‘ ，、、 』、 ，r ／福建省j 三明／ 、福州j 。 蒜舄沤／口 、I东 7 口 蓐 式验区 连城试驻口 ，1 0 海 瞄验区， 一■平潭 、‘’ b 泉州。卜’ 口 7 广东省j j 哈j ’ 湾 “金门 海 俞“ 潮州’ jL 1 f 謦{ 1 _ ，I ≯7 南 海 帮 北回归线 图5测区位置不意 F i g ．5 C S A M Ts u r v e ya r e a si nF u j i a np r o v i n c e 由于大田试验区测量剖面附近有钻孔验证，故 本文选择大田试验区作为典型代表，阐述解释 C S A M T 法勘查深部隐伏煤层的成果． 区内属中低山地貌，山势较高，地形较陡峭，平 均标高7 5 0m 以上．全年气候温暖、湿润、雨量充 沛，林木茂盛． 3 ．1 地质概况 试验区的出露地层中上元古界麻源群 P t m y 、加里东期花岗岩体 Y m 3 及上侏罗统长林 组 J 。c 、南园组 J 。n ，部分为上泥盆统桃子坑组 D 。t ；区内隐伏煤系地层为下二叠统童子岩组 P 2 t ，地层总厚度大于6 8 0m ． 童子岩组的盖层为印支期二层式推覆体．上层 为广平推覆体，其推覆体由麻源群 P t m y 和混合 花岗岩 T m 3 构成，下层为上泥盆桃子坑组一下石 炭林地组 D 。t z - C 。1 构成． 根据矿区周边的相关资料，总体为宽缓的复式 向斜．除了主体推覆断层以外，尚有N E 走向、S E 万方数据 中国矿业大学学报第3 7 卷 倾向，形成自西北往东南逐渐下落的台阶构造的正 断层，对下伏煤系地层的整体性有一定破坏作用， 煤系地层受其切割，往往保存不完整． 3 ．2 岩石电性特征 表2 为大田试验区的岩石标本测定表，采自勘 查剖面周围的相关钻孔． 表2岩石标本电性测定 T a b l e2 E l e c t r i cp r o p e r t yo fr o c ks a m p l e s 炭质泥岩、泥岩电阻率较低，约为1 0 0Q m 左 右；灰岩为15 0 0Q m 左右；砂砾岩、砂岩电阻率 23 0 8Q m ；而硅化岩、花岗岩、变粒岩这些火成岩 和老的变质岩电阻率一般较高，都在几千到几万 Q m ．上伏的老地层电阻率较高，而下伏的煤系地层 电阻率较低，具备开展地面电法工作的物性条件． 通过物性测定，结合多年来本区电法的工作成 果，综合本区及邻区资料，以及本次C S A M T 测得 的断面图结果，汇总给出本区煤系地层及几种常见 岩石电性统计结果 二叠系童子岩一文山组 P ，t ～P 。w 电阻率范 围5 0 ～8 0 0f l m ； 白垩系 K 电阻率范围1 0 0 ～10 0 0g l m ； 元古界麻源群 P t m y 电阻率范围10 0 0 ～ 2 00 0 00 m ； 下石炭一上泥盆系 C 。～D 。 电阻率范围 10 0 0 ～15 0 09 , m 花岗岩 7 m ～7 丌 电阻率范围9 0 0 ～3 80 0 0 Q m ； 二叠系栖霞组一石炭系 P ，q ～C 电阻率范围 6 0 0 ～20 0 0Q m ； 二叠系翠屏山组 P 。c p 电阻率范围4 0 0 ～ l2 0 0C Z m 上侏罗统南园组 J 。n 电阻率范围5 0 0 ～ 15 0 0D , m ；上侏罗统长林组 J 。c 电阻率范围2 0 0 ～8 0 0Q m ． 4C S A M T 反演结果和地质解释 前面提到，本文将大田工作区作为实例来解释 C S A M T 法在福建隐伏煤矿区勘查深部煤层的结 果．此剖面长度28 0 0m ，方位S E ，穿过煤层露头， 并且有钻孔验证，是本次资源勘查工作的重要剖 面，也是检验C S A M T 法探测隐伏煤系地层有效 性的一个典型实例．该剖面反演解释断面图见图 6 ． 02 8 05 6 0 8 4 0l1 2 01 4 0 01 6 8 0 1 9 6 02 2 4 02 5 2 0 2 8 0 0 水平距离佃 图6 大田测线电阻率反演解释断面 F i g ．6 I n v e r s i o na n di n t e r p r e t a t i o ns e c t i o no f s u r v e yl i n ei nD a t i a na r e a 从电阻率断面图上的电性结构特征看，横向上 可分为2 大段，第1 段为o ～11 2 0m 段，第2 段为 11 2 0 ～28 0 0m 段，下面分段论述． o ～11 2 0m 段，浅表层 标高4 0 0m 以浅 电 阻率较高，其电阻率可达5 0 0Q m 以上，尤其是o ～ 2 3 0m 段，其电阻率最高可达15 0 0Q m ，这个高阻 异常位于测线边部，地形又特别复杂，推测为干扰 异常．浅表层以下 标高4 0 0 ～一3 0 0m 电阻率较 低，为几十Q m ．结合地表露头和附近钻孑L 的地质 资料，推断浅表 标高4 0 0m 以浅 的中高阻层为 D 。一C ，l 和7 m 地层．其下部 标高4 0 0 ～一2 0 0m 的低阻层，有 6 0 0m 厚，推断为P ，t ～P ，w 煤系地 层，大体以4 0 。倾角向小号点倾斜．从老地层覆盖 于P 。新地层之上的接触关系，推断为残留推覆体． P l t ～P l w 地层以下 标高一2 0 0 ～- 5 0 0m ，电阻 率又有抬升，推断为P 。q 地层． 11 2 0m 号点的岩性边界推断为断层． 11 2 0 ～28 0 0m 段，电性特征分布在横向上不 均匀，成层性差，但仍可看出，3 0 0m 深度以浅 标 高4 0 0m 以上 为中高阻层 10 0 0n m 左右 ．3 0 0 ～7 0 0m 深 标高4 0 0 ～一3 0 0m 为中阻厚层，电 阻率约为8 0 0Q m ．7 0 0m 深 标高一3 0 0m 以下， 电阻率又变低，约为3 0 0 5 0 0Q m ．根据地质情况， 推测上部为J 。地层，中部为D 3 一C 。l 地层，下部的低 阻显示，有可能有局部残留煤系地层，深度较深． 万方数据 第1 期安志国等隐伏煤矿区C S A M T 法煤系地层勘查研究1 2 3 16 8 0m 附近有明显的低阻带，推断该低阻带 为一断层破碎带，向大号点倾斜． 图7 为大田区块测线9 0 0I T I 处的钻孔Z K 2 0 0 1 岩芯信息和该测深点的反演电阻率深度曲线对比 关系．在图6 中钻孔Z K 2 0 0 1 对应9 0 0r n 位置处， 根据电性特征解释地质情况，从浅到深依次为花岗 岩、D 。一C 。l ，P ，t ～P ，w 和P ，q 地层，反演解释同钻 孑L 岩芯对应较好，较好地验证了C S A M T 法的反 演结果的可靠性． 反演电阻率深度曲线图钻孔信息 电阻率／Q m 岩性实钻地层 层底厚 01 0 0 0 2 0 0 03 0 0 0 4 0 0 0 度／m l O O 花 岗Y m 31 1 4 岩 搋 D 3 t2 1 5 岩石 灰岩P f q2 4 0 2 0 0 I 煤 划3 0 0 隧 系 地 P t t I 4 8 7 4 0 0层 5 0 0 灰岩 P l q - P 】W 砂岩 泥岩 6 0 0 图7 测深点反演结果与钻孔信息对比 F i g ．7C o m p a r i s o nb e t w e e ni n v e r t e dr e s u l ta n d b o r e h o l ei n f o r m a t i o n 5结论 1 C S A M T 法对于复杂地形山区深部煤系地 层的勘探，是比较有效的一种方法．在复杂山区数 据采集时，采用地形改正和静态校正的数据处理方 法，提高了C S A M T 法煤系地层勘查准确度． 2 对有严重静态影响的原始数据采用有阀值 的全剖面比较静态校正处理后，反演模型的响应曲 线和原始曲线的拟合误差均小于6 ％． 针对复杂山区受地形、静态严重影响的数据处 理工作来进行研究，对于提高原始数据的可信度有 可取之处，也为C S A M T 方法研究积累了一些经 验． 致谢感谢福建省煤田地质局所给予的帮助和支 持，也感谢P r o f ．G a n q u a nX i e 和王妙月研究员给 予本稿中肯的建议和意见． 参考文献 [ 1 ] E 2 2 [ 3 ] [ 4 ] E 5 2 [ 6 ] [ 7 ] E 8 9 E 9 3 吴钦宝，陈同俊，陈风云．中国东部煤矿深部开采中 的地质勘探技术[ J ] ．地球物理学进展，2 0 0 5 ，2 0 2 3 7 0 3 7 3 ． W UQ i n b a o ，C H E NT o n g j u n ，C H E NF e n g y u n ． G e o l o g i c a le x p l o r a t i o nt e c h n i q u e s u s e di n d e e p c o a l m i n i n gi nC h i n a [ J ] ．P r o g r e s si nG e o p h y s i c s ， 2 0 0 5 ，2 0 2 3 7 0 3 7 3 ． 刘树才，刘志新，姜志海．瞬变电磁法在煤矿采区水 文勘探中的应用E J ] ．中国矿业大学学报，2 0 0 5 ，3 4 4 4 1 44 1 7 ． g I US h u c a i ，I 。I UZ h i x i n ，J A N GZ h i h a i ．A p p l i c a t i o n ’o fT E Mi nh y d r o g e o l o g i c a lp r o s p e c t i n go fm i n i n g d i s t r i c t [ J ] ．J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g ＆ T e c h n o l o g y ，2 0 0 5 ，3 4 4 4 1 44 17 ． 袁照令，沈博，李大明，等．高精度重磁勘探在推 覆体下找煤中的应用[ J ] ．物探与化探，1 9 9 9 ，2 3 4 2 8 3 2 8 9 ． Y U A NZ h a o l i n g ，S H E NB o ，L 1D a m i n g ，e ta 1 ． T h ea p p l i c a t i o no fh i g hr e s o l u t i o ng r a v i t ya n dm a g n e t i cs u r v e yi nc o a le x p l o r a t i o nu n d e rn a p p e [ J ] ．G e o p h y s i c a l G e o c h e m i c a le x p l o r a t i o n ，1 9 9 9 ，2 3 4 2 8 3 2 8 9 ． 申燕，邓爽岭，林建忠，等．贵州桐梓煤田地面7 能谱法找煤试验研究[ J ] ．新疆地质，2 0 0 6 ，2 4 2 2 1 5 2 1 7 ． S H E NY a n ，D E N GS h u a n g l i n g ，g i NJ i a n - z h o n g ，e t a 1 ．E x p l o r a t i o nf o rc o a lb yu s i n gu p p e rRr a ys p e c t r a a tt h eT o n g z ic o l l i e r yi nG u i z h o up r o v i n c e [ J ] ．X i n j i a n gG e o l o g y ，2 0 0 6 ，2 4 2 2 1 52 1 7 ． G 1 。D S T E I NMA ．M a g n e t o t e l l u r i ce x p e r i m e n t se m p l o y i n ga na r t i f i c i a ld i p o l es o u r c e [ D ] ．U n i v ．T o r o n t o 。1 9 7 1 ． G O I 。D S T E l NMA ，S T R A N G W A YDW ．A u d i o f r e q u e n c ym a g n e t o t e l l u r i cw i t h ag r o u n d e d d i p o l e s o u r c e [ J ] ．G e o p h y s i c s ，1 9 7 5 ，4 0 4 6 6 96 8 3 ． S A N D B E R GSK 。H H M A N NGW ．C o n t r o l l e d s o u r c ea u d i o f r e q u e n c ym a g n e t o t e I l u r i c si ng e o t h e r r e a le x p l o r a t i o n [ J ] ．G e o p h y s i c s ，1 9 8 2 ，4 7 1 1 0 0 1 16 ． B A T R E I 。I ．C ，JA C O B S NR ．D ．R e s u l t so facon t r o l l e d s o u r c ea u d i o f r e q u e n c ym a g n e t o t e l l u r i cs u r v e y a tt h eP u h i m a ut h e r m a la r e a [ J ] ．G e o p h y s i c s ，1 9 8 7 ， 5 2 4 6 6 5 - 6 7 7 ． 吴璐苹，石昆法，李荫槐，等．可控源音频大地电磁 法在地下水勘查中的应用研究[ J ] ．地球物理学报， 1 9 9 6 ，3 9 5 7 1 27 1 7 ． W U ，L u - p i n g ，S H IK u n - f a ，I ，IY i n h u a i ，e ta 1 ．A p 万方数据 1 2 4 中国矿业大学学报 第3 7 卷 p l i e ds t u d yo fC S A M Tm e t h o do ng r o u n dw a t e re x p l o r a t i o n [ J ] ．G e o p h y s i c aS i n i c a ，1 9 9 6 ，3 9 5 7 12 7 1 7 ． [ 1 0 ] W A N N A M A K E RPE ．T e n s o rC S A M Ts u r v e yO v e rt h es u l p h u rs p r i n g st h e r m a la r e a ，V a l l e sC a l d e r a ，N e wM e x i c o ，U ．S ．A ，P a r tI I m p l i c a t i o n sf o r s t r u c t u r eo ft h ew e s t e r nc a l d e r a [ J ] ．G e o p h y s i c s ， 1 9 9 7 ，6 2 2 4 5 卜4 6 5 ． [ 1 1 ] 王一新，钟慧智，姚治龙，等．可控源声频电磁测深 法在油气勘探中的应用[ J ] ．河南石油，】9 9 5 ，9 2 t 3 - 1 8 ． W A N GY i x i n ，Z H O N GH u i z h i ，Y A O ，Z h il o n g ， e ta 1 ．T h ea p p l i c a t i o no fC S A M Tm e t h o di ne x p l o r i n g { o ro i l g a s [ J ] ．H e n a nP e t r o l e u m ，1 9 9 5 ，9 2 1 3 1 8 ． [ 1 2 ] B A S O K U RAT ，R A S M U S S E NTM ，K A Y AC ， e ta 1 ．C o m p a r i s o no fi n d u c e dp o l a r i s a t i o na n dcon t r o l l e d s o u r c e a u d i o m a g n c t o t e l l u r