井地瞬变电磁联合探测在2000m深部找矿中的应用.pdf

井地瞬变电磁联合探测在 2 000 m 深部找矿中的应用 武军杰1,2，刘彬1,2，智庆全1,2，邓晓红1,2，王兴春1,2，杨毅1,2，刘东明1,2，邱金柱3 1. 中国地质科学院地球物理地球化学勘查研究所，河北 廊坊 065000；2. 自然资源部地球物理电磁法探测 技术重点实验室，河北 廊坊 065000；3. 辽宁省有色地质一三队有限责任公司，辽宁 丹东 118000 摘要 深部矿产探测是一个难题，井中瞬变电磁方法Transient Electromagnetic , TEM 通过传 感器沿钻孔深入地下，更接近于深部目标体，从而提高见矿率和找矿效果。在国外井中瞬变电磁被 广泛地应用于深部找矿勘查中，并取得了许多重要的找矿成果，而在国内受限于仪器、理论研究以 及解释技术发展较慢等诸多原因，在实际深部找矿应用方面较少。以辽宁青城子白云金矿区 2 000 m 深部找矿作为典型案例，开展井中 TEM 和地面 TEM 的探测，尝试将地面和井中 TEM 资料进行综 合解释，以期更合理地解释钻孔周围的有利成矿部位。实测结果表明井中 TEM 原始剖面曲线在 1 4001 650 m 深度段出现典型负异常特征，显示钻孔外局部低阻异常体的存在；此外，井中 TEM 原始曲线晚期道中出现极化现象，是由于断裂和破碎带中富集的与金密切相关的黄铁矿化引起；地 面 TEM 反演结果清楚显示了测线下方典型分布情况，发现了隐伏断裂深部延伸较大，成矿条件良好； 综合分析解释，TEM L52 线所推测的隐伏断裂深部仍有良好的金富集条件；经钻孔资料证实，地面 和井中 TEM 综合解释结果正确可靠。研究成果表明，井地瞬变电磁法联合探测在深部找矿中具有 较好应用前景。 关键词瞬变电磁法；井中 TEM；地面 TEM；深部矿产；白云金矿；青城子矿集区；极化现象 中图分类号P631 文献标志码A 文章编号1001-1986202207-0070-09 Application of surface and borehole TEM joint exploration in 2 000 m deep mineral exploration WU Junjie1,2, LIU Bin1,2, ZHI Qingquan1,2, DENG Xiaohong1,2, WANG Xingchun1,2, YANG Yi1,2, LIU Dongming1,2, QIU Jinzhu3 1. Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000, China; 2. Key Laboratory of Geophysical Electromagnetic Probing Technologies, Ministry of Natural Resources, Langfang 065000, China; 3. No.103 Branch Non-ferrous Geological of Liaoning Province Co., Ltd, Dandong 118000, China Abstract Deep mineral exploration is a challenge. The borehole Transient Electromagnetic TEM goes deep underground through the sensor along the borehole, and gets closer to the deep target, so as to improve the ore discovery rate and prospecting effect. TEM is widely used in deep prospecting and exploration in foreign countries, and has achieved many important results. However, in China, due to the slow development of instruments, theoretical research and interpretation technology, it is less used in practical deep prospecting. This paper introduces a typical case about the 2 000 m deep prospecting in Baiyun Gold Deposit in Qingchengzi, Liaoning Province. The borehole and surface TEM exploration is carried out, and the TEM data are tried to be comprehensively interpreted in order to explain the favorable metallogenic sites around the borehole more reasonably. The measured results show that the original TEM profile curve in the borehole has typical negative anomaly characteristics at the depth range of 1 400-1 650 m, indicating the exist- ence of local conductive anomalies outside the borehole. In addition, polarization in the middle and late stages is caused by the mineralization of gold-related pyrite concentrated in fracture zones. The inversion results of surface TEM clearly show the typical distribution under the survey line, and the large extension of deep concealed faults, showing good 收稿日期 2021-11-30；修回日期 2022-04-15 基金项目 国家重点研发计划项目2018YFC0603803；中央级公益性科研院所基本科研业务费专项项目AS2020Y01 第一作者 武军杰，1979 年生，男，河北栾城人，博士，正高级工程师，从事瞬变电磁方法技术研究. E-mail 第 50 卷 第 7 期煤田地质与勘探Vol. 50 No.7 2022 年 7 月COAL GEOLOGY borehole TEM; surface TEM; deep mineral resources; Baiyun Gold Deposit; Qingchengzi ore concentration area; polarization phenomenon 辽东青城子矿区位于华北克拉通东北部，是我国 重要的金、银、铅、锌等多金属矿产集中地，区内先后 发现多个大中型铅锌矿、白云、小佟家堡子、林家三 道沟等多个金矿床以及若干小型矿床和矿点1-2。其 中白云金矿位于青城子矿田北部，是最重要的大型金 矿床之一，深部找矿潜力巨大3-6。前人基础性地质研 究表明，白云金矿区经历过多期构造运动，地质结构复 杂，深部地质结构研究还存在薄弱之处6。近几年矿 区开展的航磁、航空瞬变电磁、15 万重力、音频大地 电磁等面积性工作对区域性地层、构造以及岩体分布 的进一步认识起到关键作用7-8，并预测了深部成矿有 利区9。矿区及外围重点地段开展的地面电磁法对于 探测含矿地层和控矿构造的空间形态和深部延伸情况 也获得了良好效果10-13，白云矿区深部找矿亟待突破。 深部矿产探测一直都是一个难题，而瞬变电磁 法Transient Electromagnetic , TEM 是一种矿 区深部找矿重要的地球物理找矿方法。井中 TEM 法 一般是指地井瞬变电磁法，是在地面发射，钻孔中接 收的一种瞬变电磁测量装置14-15。深部目标体受地面 发射框激发产生的异常场被井中传感器所测量，而由 于传感器在钻孔中测量，受近地表的干扰小，能够更接 近于深部目标体，获得比地面 TEM 更可靠的深部信 息16。对于规模小，埋深大的高电导目标矿体，或者存 在浅部地质干扰的情况下，其优势就更加突出17。因 此，利用井中 TEM 方法能够获得受地面发射框激发 的钻孔周围数百米范围内的有用地质信息，从而提高 见矿率和找矿效果18-19。 20 世纪，在加拿大、澳大利亚等国，井中瞬变电磁 被广泛地应用于深部找矿勘查中，并取得了许多重要 的找矿成果。在加拿大 Sudbury 铜镍矿，采用地井瞬 变电磁法先后发现了 1 280 m 深处的 Linsley 矿体和 2 400 m 以深的 Victoria 富铜镍矿床20。另外，采用 地井瞬变电磁法还在已开采的矿床之下 1 2001 500 m 深处发现了一处高品位的底板矿床，为矿山资源开发 作出重大贡献21-22。国外在资料解释方面较为成熟， 发展了一维反演、等效电流环以及三维反演等解释方 法23-27。 虽然国内受制于仪器设备、数据采集以及解释技 术等各项技术，总体发展较晚，研究程度较低，但是近 年来，国内相关研究学者对于井中瞬变电磁法的研究 保持了较大的热度，主要包括三维正演方法、不同介 质条件下不同形态异常体的三分量响应特征的分 析28-33，以及数据处理和应用方面的研究34-41。而在 实际应用案例方面，由于国内金属矿产钻孔普遍在 1 000 m 以浅，多数井中 TEM 成果反映深度也多在几 百米35,37，深井中的 TEM 成果报道更少。 以上研究多是针对团块状目标矿体如块状铜镍 硫化物矿体直接找矿的探测技术及解释方法，而现有 方法技术在白云金矿区并不完全适用。白云矿区蚀变 岩型和石英脉型金矿体本身的规模都不足以引起电磁 场的有效变化，因此，TEM 深部找矿的目标并非金矿 体本身，而是金的深部有利富集空间，包括地层中片岩 和大理岩的分界面、破碎带以及断裂的由陡变缓部 位3，这加大了井中 TEM 解释的难度。而充分利用地 面和井中 TEM 信息开展综合解释无疑会获得更合理 结果42。因此，为在白云金矿区 5002 000 m 第二找 矿空间中寻找深部金矿，在重点区开展了地面瞬变电 磁法的深部探测，并充分利用在瞬变电磁异常区实施 的 2 000 m 深孔开展了井中 TEM 测量，进而尝试将地 面和井中 TEM 资料进行综合解释，以期使这 2 种装 置优势互补，能够同时在纵向和横向上获取可靠数据， 从而更合理地解释钻孔周围的目标地质体。 1 矿区地质与地球物理概况 1.1 地质概况 白云金矿区内出露地层主要为元古界辽河群上部， 由古元古代大石桥组三段、盖县组和第四系组成。大 石桥组三段下部由中粒大理岩和白云岩大理岩组成， 中部由薄层含石墨中粒大理岩和硅质带状大理岩组成， 上部由条纹透闪石片岩和硅质条带大理岩组成。盖县 组下部为黑云母浅粒岩、变粒岩、透闪石变粒岩夹薄 层细粒大理岩，上部主要为矽线石云母片岩，黑云母片 岩、夹黑云母变粒岩、浅粒岩、透闪石变粒岩。第四 系地层由砾石、砂、黏土等现代冲洪积层组成。 研究区构造经历过多期地质运动，形成了复杂的 褶皱和系列断层。白云金矿区主体以东西向构造为主， 发育有多个褶皱和韧脆性推覆体。东西向构造从北 向南依次为白云山背斜、阳沟向斜、阳沟石湖沟背斜、 第 7 期武军杰等井地瞬变电磁联合探测在 2 000 m 深部找矿中的应用 71 姚家沟天桥沟李家堡子倒转向斜、苏家堡子翻转背 斜和顾家堡子倒转向斜5-6。此外，区内还有 NW、NE、 NS 向断裂分布。这些断层大多被石英斑岩、花岗斑 岩、闪长玢岩、二长斑岩和煌斑岩等脉岩充填。这些 岩脉大多已矿化和蚀变3。 1.2 地球物理特征 以往资料显示，片岩具有中高电阻率，大理岩呈高 阻，而透辉透闪片岩和含石墨大理岩为低阻显示。其 余岩性单元中石英斑岩、硅钾蚀变岩平均电阻率为高 阻，含金蚀变岩和硅化带一般表现为中高阻，区内岩石 的电阻率受石墨、破碎带较大影响。结合物性搜集结 果和钻孔综合测井情况，总结了工作区地层电性情况， 其中盖县组分为两层以片岩为主，总体呈高阻、低极 化特征。大石桥组三段上层呈低阻、低极化特征，其 中低电阻主要由透辉透闪片岩引起；中部由薄层含石 墨大理岩和硅质条带大理岩组成，呈低阻特征，局部破 碎带呈高极化特征；下层主要为白云石大理岩，呈高阻 特征。 2 野外数据采集 2.1 采集参数 为探测白云矿区深部矿体，笔者团队在白云矿区 姚家堡子一带开展了地面定源回线瞬变电磁法测量， 完成多条 TEM 测线。其中，L52 线是较为典型的一条 剖面，反演结果显示深部仍为低阻，存在大规模断裂破 碎带，有利于深部成矿。当地矿山企业在 TEM 异常 区实施了 2 000 m 钻探，在深部见矿。之后在钻孔中 开展了井中瞬变电磁法，进一步探测深部地质结构的 空间分布。TEM 地面测线和钻孔位置如图 1 所示。 地面及井中野外数据采集参数见表 1。其中，地 面点距 50 m，钻孔中测量基本点距为 5 m，一般会在异 常段加密至 12 m。 2.2 仪器设备 野外数据采集中使用加拿大 Crone PEM 系统，该 姚家堡子 顾家堡子 倒转向斜 白云山村 L52 600 0 ZK1 N 荒甸子 姚家沟 顾家沟 地质界线 盖县组浅粒岩 盖县组变粒岩 盖县组片岩 第四系 大石桥组大理岩 变粒岩互组 大石桥组透辉透 闪片岩 大石桥组 角闪片岩 大石桥组硅质条 带大理岩 大石桥组大理岩 TEM 测点TEM 发射框 倒转向斜轴 0 100 200 m 向斜轴 倒转背斜轴 背斜轴 石英斑岩 闪长斑岩 花岗斑岩 图 1 白云矿区地面及井中 TEM 野外布置 Fig.1 Field layout of surface and borehole TEM in Baiyun mining area 72 煤田地质与勘探第 50 卷 系统包括 4.8 kW 的发射系统和接收系统。发射系统 由 5 kW 发电机、整流器、发射机和回线框组成。发 电机产生的 220 V 交流电经整流器转换为直流电压， 发射机根据采集参数在回路中发射双极梯形波。发射 系统和接收系统的同步方式是石英钟。接收系统包 括 PEM 接收机、地面探头和井下探头。井中瞬变电 磁系统照片如图 2 所示。 发射机 整流器 发电机 电缆、 绞车 接收器 探 头 图 2 Crone 井中 TEM 系统 Fig.2 Borehole TEM system of Crone 3 结果及讨论 3.1 地面 TEM 反演结果 地面 TEM 测线 L52 位置如图 1 所示，图 3 为地 面定源回线瞬变电磁法 L52 线反演和解释结果。由 图 3 可知，测线下方的电性分布总体分为五层，300 m 以浅呈高阻，为盖县组片岩的反映；300500 m 呈低阻， 为透辉透闪片岩及含石墨大理岩的反映；5001 400 m 总体呈相对高阻，而高阻特征为南部薄，北部厚； 1 4001 800 m 为低阻，推断为破碎带；1 800 m 以下高 阻为白云大理岩反映，属大石桥组三段下部。 点号 深度/m 0 200 400 600 800 1 000 1 200 1 400 1 600 1 800 2 000 0100 m 推断断裂 100200300400500 N 600 s/m 230 220 140 90 82 72 62 54 40 图 3 地面瞬变电磁 L52 线反演断面图 Fig.3 Inversion section of surface TEM L52 3.2 井中 TEM 曲线特征 为了解不同时间道瞬变响应沿钻孔的变化特征， 将 31 个时间道分为 4 组分别绘制。图 4 为钻孔 ZK1 中实施的井中 TEM 实测剖面曲线图，测量深度 950 1 700 m。图中可以看出，4 组剖面曲线在 9501 300 m 段仅在个别段如 9701 020、1 100、1 140、1 220 m 等 出现微弱异常，为岩性互层中含炭质薄低阻层的反映。 而第 4 组2531 道晚期道剖面曲线在 1 4001 650 m 段出现明显负异常图 5，为典型的孔旁异常特征，显 示该段孔旁一定距离存在低阻体。该段钻孔编录显示 为煌斑岩、破碎带、蚀变带、含石墨大理岩等，为深部 断裂特征。 为进一步分析该段井中 TEM 的瞬变响应特征， 表 1 野外瞬变电磁数据采集参数 Table 1 Parameters of field data collection for TEM 项目地面TEM井中TEM 点距/m 505 发射机功率/kW 4.84.8 发射框/mm 400800300300 发射电流/A 2020 发射时间/ms 5020 关断时间/s 500500 探头类型地面井中 探头等效面积/m2 3 8506 000 接收信号感应电动势感应电动势 叠加次数 256512256512 同步方式石英钟电缆 第 7 期武军杰等井地瞬变电磁联合探测在 2 000 m 深部找矿中的应用 73 选择 1 0001 600 m 段的原始衰减曲线进行对比图 6。 由图 6 可知，井中 TEM 衰减曲线衰减规律与地面曲 线衰减规律不同。在早、中期道中，瞬变响应幅值逐 渐增大之后再衰减，这主要由测点深度、发射框与钻 孔的相对位置关系以及发射框与下方介质耦合关系综 合导致。1 000、1 100、1 200 和 1 300 m 晚期道曲线 衰减正常，而 1 400 m 开始在晚期道出现反号现象，而 且随着深度增加，反号的幅值逐渐增大。该现象为典 型的极化现象，为深部破碎带蚀变带中局部富集的石 墨化、黄铁矿化导致。而前期研究表明，该矿区金成 矿与黄铁矿化密切相关，因此，钻孔中出现的极化现象 在一定程度上指示了一定范围内金的富集。而钻孔编 录也证实了这一结论。 3.3 TEM 井地联合解释结果 为更详细了解钻孔周围地质结构以及研究矿区深 部金的有利富集部位，尝试将地面 TEM 结果与井中 TEM 测量结果进行综合对比分析。测线 L52 和钻孔 ZK1 的相对位置关系如图 1 所示，钻孔与地面测线直 线距离约为 400 m。考虑到两者同处顾家堡子倒转向 斜北翼的同一构造区，地质背景相似，且地面和井中 TEM 的地面发射框亦有叠合之处，对于深部 2 000 m 目标体的激发也有相似的激发条件，因此，尝试将地 面 L52 线 TEM 反演结果和钻孔 ZK1 井中 TEM 实测 结果进行联合解释，以分析 2 000 m 深部地质结构及 成矿有利部位。 感应电动势/nV 深度/m 9005 1 000 1 100 1 2002531 道 1 300 1 400 1 500 1 600 1 700 05101520 感应电动势/nV 深度/m 9000 1 000 1 100 1 2001724 道 1 300 1 400 1 500 1 600 1 700 4080 120 160 200 感应电动势/nV 深度/m 9000 1 000 1 100 1 200916 道 1 300 1 400 1 500 1 600 1 700 100200300 感应电动势/nV 深度/m 9000 1 000 1 100 1 20018 道 1 300 1 400 1 500 1 600 1 700 100200300 图 4 钻孔 ZK1 井中 TEM 剖面曲线 Fig.4 ZK1 borehole TEM plot curves 感应电动势/nV 时间道 24 25 26 27 28 29 30 31 深度/m 1 3004 1 400 1 500 1 600 1 700 04810226 图 5 钻孔 ZK1 井中 TEM 剖面曲线1 3001 700 m Fig.5 TEM plot curves of ZK1 borehole at depths of 1 300 m to 1 700 m 74 煤田地质与勘探第 50 卷 图 7 为地面和井中 TEM 综合解释图，综合了 L52 线瞬变电磁反演结果、井中 TEM 曲线、视电阻率曲线 以及钻孔编录资料。 首先，地面 TEM 反演结果中纵向电性分布与视 电阻率测井曲线基本一致，地表 300 m 以浅为高阻， 300500 m 为低阻，5001 400 m 整体呈高阻含局部低 阻薄层，而 1 4001 750 m 呈低阻。其中 300 m 处的 高低阻界面为片岩和大理岩分界面。是浅部金的有利 成矿层位。 其次，在井中瞬变电磁 1 4001 650 m 段出现的负 异常，是由 0200 号点下方 1 4001 800 m 段低阻团 块引起的，其异常中心距离钻孔约 100 m。在点号 200500 号下方不同深度多处出现明显的电阻率不连 续之处，推断为隐伏断裂，如图 7 所示。该隐伏断裂形 态呈上陡下缓的特征，而这种形态为深部金及其伴生 的黄铁矿化的局部富集提供了条件，而正是富集的黄 铁矿化引起了井中 TEM 响应晚期出现极化效应。目 前在 1 4001 500 m 处已经钻探见矿，还可推断深部 1 6002 000 m 仍是深部成矿有利区。 图 7 中文字部分是根据钻孔编录资料总结的大致 时间/ms 1 000 m 102 100 101 102 103 101100101102 感应电动势/nV 时间/ms 1 300 m 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 500 m 虚线为负值虚线为负值虚线为负值 102 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 400 m 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 550 m 102 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 450 m 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 600 m 102 102 101 100 101 102 101100101102 感应电动势/nV 时间/ms 1 100 m 102 100 101 102 103 101100101102 感应电动势/nV 时间/ms 1 200 m 102 100 101 102 103 101100101102 感应电动势/nV 图 6 不同深度测点井中 TEM 衰减曲线对比 Fig.6 Comparison of borehole TEM decay curves at different depths 第 7 期武军杰等井地瞬变电磁联合探测在 2 000 m 深部找矿中的应用 75 岩性情况，地面和井中 TEM 结果与钻孔编录资料基 本吻合，也证实了地面 TEM 的探测深度和反映的电 性结构的可靠性。 点号 深度/m 0 200 400 600 800 1 000 1 200 1 400 1 600 1 800 2 000 0100 m 推断断裂 白云石大理岩 含石墨大理岩、煌斑岩、 蚀变带、黄铁矿化 煌斑岩、蚀变带 片岩、变粒岩、 大理岩、破碎带 片岩 ZK01 ZK01 透辉透闪片岩 含石墨大理若 钻孔 金矿化 地井 TEM 曲线 视电阻率 测井曲线 100200300400500 N 600 s/m 230 220 140 90 82 72 62 54 40 图 7 地面和井中瞬变电磁综合解释结果 Fig.7 Comprehensive interpretation result of surface and borehole TEM 4 结 论 a. 井中 TEM 勘探，由于接收探头沿钻孔深入地 下，获得更加可靠的深部地质体及目标矿体的信息，包 括局部低阻异常体引起的剖面曲线异常特征，以及断 裂中富集的与金密切相关的黄铁矿化引起的激电效应， 因此，井中 TEM 能够在深部金矿探测中起到重要 作用。 b. 辽东白云矿区 2 000 m 深孔地面和井中 TEM 综合解释结果表明，地面 TEM 能够大致反映总体的 地下电性结构，井中 TEM 深入地下可以获取更为可 靠的纵向电性分布情况，两者优势互补可以提供更加 合理的解释成果。 c. 井中 TEM 方法由于充分利用钻孔深入地下，因 此，可以获得高纵向分辨率的可靠深部信息。本次仅 是将地面和井中联合进行了综合解释，建议下一步深 入研究地面和井中数据的联合反演方法，获得更加合 理可靠的地下电性分布情况，从而更加广泛应用于金 属矿、煤炭、地下水等领域深部勘探中。 致谢致谢 辽宁省有色地质一三队有限责任公司 刘福兴总工、王伟正高级工程师以及李生辉高级工程 师提供了白云矿区地质资料，并且为地面、地-井 TEM 野外施工提供诸多帮助，在此表示感谢。 参考文献References LIU Jun，LIU Fuxing，LI Shenghui，et al. ation of the Baiy- un gold deposit，Liaodong gold province，NE ChinaConstraints from zircon UPb age，fluid inclusion，and CHOPbHe iso- topesJ. Ore Geology Reviews，2018，104686706. 1 SUN Guotao，ZENG Qingdong，LI Taiyang，et al. Ore genesis of the Baiyun gold deposit in Liaoning Province，NE China Con- straints from fluid inclusions and zircon U-Pb agesJ. Arabian Journal of Geosciences，2019，129117. 2 王伟，刘福兴，郭强，等. 白云金矿床60号脉带原生晕分带特征 及深部含矿性评价J. 矿床地质，2021，403523538. WANG Wei，LIU Fuxing，GUO Qiang，et al. Characteristics of primary halo and uation of deep mineralization along No.60 vein zone of Baiyun gold depositJ. Mineral Deposits， 2021， 403523538. 3 曾庆栋，陈仁义，杨进辉，等. 辽东地区金矿床类型、成矿特征 及找矿潜力J. 岩石学报，2019，35719391963. ZENG Qingdong，CHEN Renyi，YANG Jinhui，et al. The metal- logenic characteristics and exploring ore potential of the gold de- posits in eastern Liaoning ProvinceJ. Acta Petrologica Sinica， 2019，35719391963. 4 刘国平，艾永富. 辽宁白云金矿床某些基本问题探讨J. 矿床 地质，1999，183219225. LIU Guoping，AI Yongfu. A discussion on some major problems of the Baiyun gold deposit，eastern LiaoningJ. Mineral Deposits， 1999，183219225. 5 张拴宏，胡国辉，李建锋，等. 辽东白云小佟家堡子矿集区控 矿构造及成矿有利区预测J. 地球科学，2020，45113885 3899. ZHANG Shuanhong，HU Guohui，LI Jianfeng，et al. Orecon- trolling structures and metallogenic favorable area prediction in BaiyunXiaotongjiabuzi ore concentration area，eastern Liaon- ing ProvinceJ. Earth Science，2020，451138853899. 6 程莎莎，彭莉红，孙栋华，等. 青城子矿集区深部地质构造探测 及找矿意义J. 物探与化探，2021，454859868. CHENG Shasha，PENG Lihong，SUN Donghua，et al. The ex- ploration of deep geological structure in the Qingchengzi ore concentration area and its prospecting significanceJ. Geophys- ical Geochemical Exploration，2021，454859868. 7 ZHAO Weijun，SHA Deming，HUAN Hengfei，et al. Prelimin- ary structure framework and concealed intrusive rocks in the Qingchengzi ore field in northeast China disclosed by largescale twodimensional audiomagnetotelluric soundingC//The Soci- ety of Exploration Geophysicist and the Chinese Geophysical So- ciety GEM 2019. Xian，May 19-22，2019. 8 底青云，薛国强，雷达，等. 华北克拉通金矿综合地球物理探测9 76 煤田地质与勘探第 50 卷 研究进展以辽东地区为例J. 中国科学地球科学，2021， 51915241535. DI Qingyun，XUE Guoqiang，LEI Da，et al. Summary of techno- logy for a comprehensive geophysical exploration of gold mine in North China CratonJ. Science China Earth Sciences，2021， 51915241535. 武军杰，智庆全，邓晓红，等. 辽东白云金矿区深部地质结构的 瞬变电磁法探测J. 地球科学，2020，451140274037. WU Junjie，ZHI Qingquan，DENG Xiaohong，et al. Exploration of deep geological structure of Baiyun gold deposit in eastern Liaoning Province with TEMJ. Earth Science， 2020， 4511 40274037. 10 王兴春，邓晓红，陈晓东，等. 基于高温超导的瞬变电磁法在青 城子矿集区的应用J. 地球科学，2021，46518711880. WANG Xingchun， DENG Xiaohong， CHEN Xiaodong， et al. Application effect of TEM based on high temperature supercon- ducting sensor in Qingchengzi oreconcentrated areaJ. Earth Science，2021，46518711880. 11 WU Junjie，CHEN Xiaodong，YANG Yi，et al. Application of TEM based on HTS SQUID magnetometer in deep geological structure exploration in the Baiyun gold deposit，NE ChinaJ. Journal of Earth Science，2021，32117. 12 WU Junjie，ZHI Qingquan，DENG Xiaohong，et al. Deep gold exploration with SQUID TEM in the Qingchengzi orefield，east- ern Liaoning，northeast ChinaJ. Minerals，2022，12102. 13 BARNETT C. Simple inversion of timedomain electromagnetic dataJ. Geophysics，1984，497925933. 14 武军杰，李貅，智庆全，等. 电性源地井瞬变电磁全域视电阻 率定义J. 地球物理学报，2017，60415951605. WU Junjie，LI Xiu，ZHI Qingquan，et al. Full field apparent res- istivity definition of borehole TEM with electric sourceJ. Chinese Journal of Geophysics，2017，60415951605. 15 EADIE T，STALTARI G. Introduction to downhole electromag- netic sJ. Exploration Geophysics，1987，183247351. 16 THOMAS L. Short noteA simple interpretation aid for electro- magnetic anomaliesJ. Exploration Geophysics， 1987， 183 349351. 17 HUGHES N A，RAVENHURST W. Three component DHEM surveying at BalcoomaJ. Exploration Geophysics，1996，273 7789. 18 ADAMSON M，RAY B，HUIZI A. Using low frequency ground and downhole TDEM to explore for massive sulfide mineralisa- tion in the Carajs mineral ProvinceJ. ASEG Ex