多物探方法在胶东某金矿采空区勘查中的应用研究

第５５卷㊀第３期２０１９年５月地质与勘探ＧＥＯＬＯＧＹＡＮＤＥＸＰＬＯＲＡＴＩＯＮＶｏｌ.５５㊀Ｎｏ.３Ｍａｙꎬ２０１９ｄｏｉ:１０.１２１３４/ｊ.ｄｚｙｋｔ.２０１９.０３.０１２[收稿日期]２０１８－０６－１５ꎻ[改回日期]２０１９－０２－２２ꎻ[责任编辑]陈伟军ꎮ[基金项目]中国冶金地质总局山东局２０１６年度科技研发计划项目资助ꎮ[第一作者]张㊀扬(１９８４年－)ꎬ男ꎬ２０１４年毕业于中国地质大学(北京)ꎬ获硕士学位ꎬ工程师ꎬ注册安全工程师ꎬ国家(二级)安全评价师ꎬ主要从事地质勘查安全工程技术应用及安全管理工作ꎮＥ￣ｍａｉｌ:ｃｈｉｎａ＿ｚｙ１１１５＠１６３.ｃｏｍꎮ[通讯作者]王㊀璐(１９８６年－)ꎬ男ꎬ２０１０年毕业于吉林大学ꎬ工程师ꎬ主要从事地球物理勘查工作ꎮＥ￣ｍａｉｌ:６１６１０２６４３＠ｑｑ.ｃｏｍꎮ多物探方法在胶东某金矿采空区勘查中的应用研究张㊀扬１ꎬ王㊀璐１ꎬ刘智慧２ꎬ刘升台１ꎬ梅贞华１ꎬ杜利明１ꎬ林晓晖３ꎬ张建波４(１.中国冶金地质总局山东正元地质勘查院ꎬ山东济南㊀２５００００ꎻ２.中国冶金地质总局ꎬ北京㊀１０００２６ꎻ３.中国冶金地质总局二局ꎬ福建福州㊀３５００００ꎻ４.山东正元建设工程有限责任公司ꎬ山东济南㊀２５００００)[摘㊀要]近年来ꎬ采空区塌陷事故频发ꎬ采空区勘查重要性越来越高ꎬ为了更有效地开展采空区勘查工作ꎬ在胶东某金矿开展了采空区勘查试验ꎬ探索多物探方法在采空区勘查中的应用效果ꎮ在严重干扰条件下的矿区内开展了瞬变电磁法的多种装置类型㊁直流高密度电阻率法及浅层地震等多种物探方法选择试验ꎮ最终选择了高密度电法及等值反磁通瞬变电磁法进行了试验剖面ꎮ根据钻孔工程验证ꎬ证明了两种方法的有效性ꎬ同时对采空区充水条件的分析ꎬ使得推断结果更合理准确ꎬ对采空区勘查总结了一些经验ꎮ[关键词]㊀金矿㊀采空区㊀高密度电阻率法㊀等值反磁通瞬变电磁法[中图分类号]Ｐ６１８ ５１ꎻＰ６３１㊀㊀[文献标识码]Ａ㊀㊀[文章编号]０４９５－５３３１(２０１９)０３－０９ＺｈａｎｇＹａｎｇꎬＷａｎｇＬｕꎬＬｉｕＺｈｉｈｕｉꎬＬｉｕＳｈｅｎｇｔａｉꎬＭｅｉＺｈｅｎｈｕａꎬＤｕＬｉｍｉｎｇꎬＬｉｎＸｉａｏｈｕｉꎬＺｈａｎｇＪｉａｎｂｏ.Ａｐｐｌｉｃａｔｉｏｎｏｆｍｕｌｔｉ－ｇｅｏｐｈｙｓｉｃａｌｍｅｔｈｏｄｓｉｎｅｘｐｌｏｒａｔｉｏｎｏｆａｗｏｒｋｅｄ－ｏｕｔａｒｅａｏｆｇｏｌｄｍｉｎｉｎｇｉｎｔｈｅＪｉａｏｄｏｎｇｒｅｇｉｏｎ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ２０１９ꎬ５５(３):０８０９－０８１７.０㊀引言近年来ꎬ全国因采空区塌陷诱发的各种事故时有发生ꎬ造成了巨大的社会安全隐患和人民生命财产损失(张永双等ꎬ２０１７)ꎮ采空区作为影响目前矿山安全生产的主要危害源之一ꎬ存在巨大的安全隐患ꎬ研究如何有效解决采空区的勘查问题ꎬ对保护人民财产及人身安全具有重大意义ꎬ有着高度的社会价值(邓超文ꎬ２００７ꎻ范亮和荣毅ꎬ２０１１ꎻ朱亚军和王艳新ꎬ２０１２)ꎮ就目前的采空区探测技术而言ꎬ在欧美等相对发达国家主要以地球物理勘探为主要手段ꎻ国内目前还是主要依靠传统的钻探工程揭露ꎬ地球物理勘探技术只能起到一定的辅助作用ꎬ未能发挥出先进技术的优势(郭有刚等ꎬ２００４ꎻ吕惠进等ꎬ２００５ꎻ王华锋等ꎬ２０１３)ꎮ胶东地区作为重要的金矿资源集中地ꎬ矿山开采活动历史悠久ꎬ因采空区引发的地面㊁道路塌陷及房屋开裂事件频发(胡博文等ꎬ２０１５)ꎬ如何有效开展采空区勘查的问题较为紧迫ꎬ结合我单位在该地区的工作基础ꎬ本次在胶东某金矿开展采空区勘查综合物探方法的应用研究ꎬ通过在典型剖面上开展多种物探方法实验ꎬ对采空区勘查总结了一些经验ꎮ１㊀试验区概况试验区位于烟台蓬莱市某金矿矿山ꎬ矿山从１９９３年正式投产ꎬ已开采２５年ꎬ之前民采历史可追溯到明末清初ꎬ采空区问题较为严重ꎮ矿区内地形起伏较大ꎬ海拔９０~２８０ｍꎬ属沿海低山丘陵区(宋明春等ꎬ２００７)ꎮ区内采金业发达ꎬ中小型金矿山林立ꎬ电网密布ꎬ各种人文干扰复杂ꎬ开展地球物理工作的条件较为恶劣ꎮ工作区地层简单ꎬ仅见新生代第四纪ꎮ区内构造较为发育ꎬ区内主要的矿脉及矿体的分布基本由断裂构造控制ꎮ区内岩浆岩广泛分布ꎬ主要有新太古代栖霞超单元回龙夼单元㊁元古代玲珑超单元ꎬ中生代郭家岭超单元及中生代燕山晚期的各种脉岩ꎮ矿区出露的主体岩性有二长花岗岩㊁绢英岩化花岗质碎裂岩㊁断层泥㊁绢英岩和石９０８地质与勘探２０１９年英脉等及少量第四系松散堆积物等(宋国政等ꎬ２０１８)ꎮ矿区的主要隔水层为二长花岗岩ꎮ浅表部位有一层基岩裂隙水含水层ꎬ富水性一般较弱ꎬ伴随着深度的增加ꎬ风化程度越来越弱ꎬ逐渐的变化为相对的隔水层ꎮ位于开拓基准面以上ꎬ该含水层空间分布极不均匀ꎬ厚度变化较大ꎬ对深部开采不构成威胁ꎻ由绢英岩化花岗质碎裂岩㊁断层泥㊁绢英岩和含矿石英脉等组成的构造裂隙水含水层ꎬ是矿山开采地段主要的含水层ꎬ同时是直接充水的主要因素ꎮ因此ꎬ该矿床按照充水类型划分ꎬ属于构造裂隙水充水型的矿床(周彦章等ꎬ２００８ꎻ宋明春等ꎬ２０１４)ꎮ目前ꎬ断裂破碎带的构造裂隙水为影响矿坑充水的直接因素ꎬ大气降水渗透是影响矿坑充水的间接因素ꎮ近年来矿山新开采造成的采空区已全部实施了回填处理ꎬ目前危害较大的主要为责任灭失的老采空及旧的民采采坑(张强等ꎬ２００６)ꎮ花岗闪长岩㊁玲珑混合花岗岩及金矿石属高阻地质体ꎬ电阻率值在２０００Ω ｍ以上ꎬ其它各类岩矿石电阻率相近ꎬ在７３５~１８２０Ω ｍ之间ꎬ属中阻岩石ꎻ其中ꎬ碎裂岩为低阻岩石ꎬ明显低于围岩的电阻率ꎮ因此ꎬ胶东西部各类岩石之间存在明显的电阻率差异ꎬ其特征如下:郭家岭花岗闪长岩㊁玲珑岩体为高阻岩石ꎬ具有高阻电阻率异常特征ꎻ金矿石为高阻特征ꎬ由于矿体受断裂破碎带的控制ꎬ使岩石破碎呈现低阻特征ꎻ二者之间存在较为显著的电性特征差异ꎬ具体表现为:矿体所体现的电阻率值应该呈低阻中的相对高阻特征ꎻ断裂构造由于岩石破碎和断层泥充填ꎬ使得电阻率值明显低于围岩ꎬ因此在电性上呈低阻异常特征(曲少飞ꎬ２００８ꎻ董健等ꎬ２０１８)ꎮ从大柳行金矿收集的资料来看ꎬ矿区内分布的中生代石英脉具高阻特征ꎬ电阻率在９０００Ω ｍ以上ꎬ花岗岩类呈中阻特征ꎬ电阻率在３０００~４０００Ω ｍ之间ꎬ碎裂岩类呈低阻特征ꎬ电阻率小于１０００Ω ｍꎬ本区分布的各类岩性存在电性差异ꎬ具备开展高密度直流电阻率法㊁瞬变电磁法等地球物理勘查的物性条件(刘国兴ꎬ２００５ꎻ江玉乐等ꎬ２００７)ꎮ２㊀方法选择考虑到本次试验区属于重干扰区ꎬ电磁干扰等较为强烈ꎬ在试验正式开展前先进行了方法选择试验ꎬ寻找出能够适应场地条件的地球物理勘查方法或方法组合ꎮ２ １㊀瞬变电磁大定源装置试验本次瞬变电磁大定源装置测量采用美国ＺＯＮＧＥ公司生产的ＧＤＰ３２电法工作站ꎬ发射系统由３０千瓦发电机㊁发射机(ＧＧＴ－１０)和发射机控制器(ＸＭＴ－３２Ｓ)组成ꎬ接收机为ＧＤＰ３２ＩＩ接收机ꎬ采用ＴＥＭ－７０ｋ磁探头ꎮ发射框(Ｔｘ)布设为边长２００ｍ的正方形ꎮ接收ＴＥＭ探头(Ｒｘ)国产Ｔ－７Ｋ瞬变电磁探头ꎬ等效面积设为２５００ｍ２ꎬ平均发射电流１６Ａꎮ本区有十几家大小矿山处在生产状态ꎬ区内高压线林立ꎬ地下管线较多ꎬ属于电磁干扰严重地区ꎬ从采集到的原始数据电磁响应信号衰减异常ꎬ外部干扰信号能量较大ꎬ难以采集到足够的有效信息ꎬ数据可信度较差ꎬ不能满足正常观测的要求ꎮ２ ２㊀瞬变电磁中心回线装置本次瞬变电磁中心回线装置测量采用澳大利亚Ａｌｐｈａ公司生产的ＴｅｒｒａＴＥＭ瞬变电磁仪ꎬ具有多道㊁数字发射机和接收机一体化㊁自选时窗等功能ꎬ仪器的稳定性和抗干扰能力较好(金中国和邹林ꎬ２００２ꎻ周平和施俊法ꎬ２００７ꎻ孙怀凤等ꎬ２０１６)ꎮ发射框Ｔｘ边长使用了５０ˑ５０ｍ单匝ꎬ点距２０ｍꎮ接收线圈采用了０ ５ˑ０ ５ｍ(５００匝ꎬ等效面积１２５ｍ２)ꎮ发射电流均在９Ａ左右ꎬ叠加次数１６次ꎮ本次试验采集到的原始数据曲线图(图１)ꎮ为了进行对比ꎬ我们又在远离矿区ꎬ没有高压线干扰的地区进行了对比试验ꎬ采集到的原始数据曲线图(图２):从图１和图２２条电磁响应曲线看ꎬ在矿区内的测量结果与无干扰地区的测量结果相比ꎬ数据的稳定性差ꎬ电磁信号衰减趋势有较多跳点㊁畸变点ꎬ数据质量差ꎬ无法较好的解决电磁干扰的问题ꎬ无法满足工作要求ꎮ２ ３㊀瞬变电磁等值反磁通装置ＨＰＴＥＭ－０８等值反磁通瞬变电磁系统ꎬ它包括ＨＰＴＥＭ天线㊁仪器主机㊁１２Ｖ外接电池和操作ＰＣ组成ꎮ采集参数:等值反磁通瞬变电磁法选择半径为１ｍ左右的发射线圈(简称 微线圈 )ꎬ最大发射电流达３０Ａꎮ接收天线为内置ＨＰＴＥＭ天线ꎬ采用２４ｂｉｔꎬ６２５ｋｓｐｓ的采样率(采样间隔１ ６μｓ)ꎬ以多周期叠加采样取平均值作为实测记录信号ꎮ数据处理:采用曲线拟合的方法进行滤波ꎬ保留０１８第３期张㊀扬等:多物探方法在胶东某金矿采空区勘查中的应用研究图１㊀瞬变电磁响应曲线(Ｔｘ:５０ˑ５０ｍꎬ矿区内)Ｆｉｇ.１㊀Ｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｒｅｓｐｏｎｓｅｃｕｒｖｅ(Ｔｘ:５０ˑ５０ｍꎬｍｉｎｅ)１－原始观测数据ꎻ２ꎬ３ꎬ４－重复观测数据１－ｒａｗｄａｔａꎻ２ꎬ３ꎬ４－ｒｅｐｅａｔｅｄｍｅａｓｕｒｉｎｇｄａｔａ图２㊀瞬变电磁响应曲线(Ｔｘ:５０ˑ５０ｍꎬ无干扰地区)Ｆｉｇ.２㊀Ｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｒｅｓｐｏｎｓｅｃｕｒｖｅ(Ｔｘ:５０ˑ５０ｍꎬｎｏｎｄｉｓｔｕｒｂｉｎｇａｒｅａ)１－原始观测数据ꎻ２－重复观测数据１－ｒａｗｄａｔａꎻ２－ｒｅｐｅａｔｅｄｍｅａｓｕｒｉｎｇｄａｔａ数据的整体性ң使用５ＤＴＥＭ软件进行反演计算ꎬ对单点数据进行一维反演ꎬ同时对多点数据进行拟二维剖面联合反演(武军杰ꎬ２０１５)ꎮ该系统使用了等值反磁通的方法ꎬ可以起到消除收发线圈之间的耦合效应的作用ꎻ利用对偶中心耦合原理提高横向分辨率ꎮ它很好地解决了瞬变电磁浅层测量 盲区 以及抗电磁干扰弱的问题(李文尧ꎬ２０１０ꎻ韩进国等ꎬ２０１８)ꎮ零磁通面接收纯二次场ꎬ中心耦合优势ꎻ高密度小间隔采样ꎬ根据探测分辨率要求灵活抽道ꎮ抗干扰能力强ꎬ数据准确度更高ꎬ纵向分辨率更高ꎬ观察实测衰减曲线(图３)ꎬ数据质量满足工作需要ꎮ１１８地质与勘探２０１９年图３㊀等磁反磁通瞬变电磁法接收信号免受一次场干扰(矿㊀㊀㊀区内实测)Ｆｉｇ.３㊀ＯＣＴＥＭｒｅｃｅｉｖｉｎｇｓｉｇｎａｌｆｒｏｍｐｒｉｍａｒｙｆｉｅｌｄｉｎｔｅｒ￣㊀㊀㊀ｆｅｒｅｎｃｅ(ｍｅａｓｕｒｅｄｉｎｍｉｎｉｎｇａｒｅａ)２ ４㊀浅层地震(地震映像)本次浅层地震测量工作采用ＳＥ２４０４ＥＰ型数字工程地震探测仪ꎬ采用单道地震反射波法(又称为地震映像法)ꎮ测量前进行了波场调查ꎬ偏移距测试等工作ꎬ震源使用锤击震源(葛双成等ꎬ２００８)ꎮ受本区地下采矿活动影响ꎬ地下干扰震动剧烈且频繁ꎬ采集到的信号质量差ꎬ同相轴显示不明显ꎬ难以分辨有效信息与干扰信息ꎬ使用地震方法在本区效果不理想ꎮ２ ５㊀高密度电法试验本次高密度试验在已知民采空硐上布设高密度直流电阻率试验剖面一条ꎮ采集参数:本次试验剖面采用了温纳装置ꎬ电极排列方式为分布式ꎬ最大供电电压２５０Ｖꎮ电极距２０ｍꎬ最大隔离系数１６(祁民等ꎬ２００６ꎻ马德锡等ꎬ２００８)ꎮ数据处理:将传出的数据逐层检查ң对突变点进行三点圆滑ң加入地形点高程数据ң将数据整理成反演格式ң二维反演采用ＧＤＦ５ ０二维电法反演软件ꎬ把实测视电阻率数据和地形数据整理成反演数据文件ꎬ采用水平柱状体组合模型进行最小二乘反演ꎬ最终得到带地形的反演成果(董浩斌和王传雷ꎬ２００３ꎻ郭秀军等ꎬ２００４ꎻ杨振威ꎬ２０１２)ꎮ图４㊀试验剖面电阻率反演断面图Ｆｉｇ.４㊀Ｓｅｃｔｉｏｎａｌｍａｐｏｆｔｈｅｒｅｓｉｓｔｉｖｉｔｙｉｎｖｅｒｓｉｏｎｏｆｔｈｅｔｅｓｔｓｅｃｔｉｏｎ２１８第３期张㊀扬等:多物探方法在胶东某金矿采空区勘查中的应用研究在试验剖面中部点距１０ｍ和点距２０ｍ的拟断面等值线均呈现扭曲并收缩特点ꎬ反演断面图上(图４)同时呈现浅层电性特征反映为高阻ꎬ深部反映为次高阻区的变化ꎮ此处地表以下为埋深在４ｍ左右的民采平硐ꎬ硐口约３ˑ２ｍꎮ深部的次高阻区推测为基岩面ꎬ深部的低阻区推测为破碎带或含水空区ꎮ综合以上五种方法的试验效果ꎬ为了达到最好的探测效果ꎬ最终选择了使用高密度电法与等值反磁通瞬变电磁相结合的方法开展进一步的试验工作ꎮ既能够在浅部和深部同时得到可靠的高分辨率信息ꎬ又具有很好的抗电磁干扰能力(王启军等ꎬ２００９)ꎮ３㊀试验情况３ １㊀高密度电法试验剖面本次试验选择已知信息较多的地区ꎬ完成了高密度电法试验剖面１条ꎬ剖面沿北东方向布设ꎬ编号为１号测线ꎮ剖面成果解释如下(图５):图５㊀１线地质物探综合剖面图Ｆｉｇ.５㊀１ｌｉｎｅｇｅｏｌｏｇｉｃａｌｇｅｏｐｈｙｓｉｃａｌｃｏｍｐｒｅｈｅｎｓｉｖｅｐｒｏｆｉｌｅ１－矿脉ꎻ２－推断采空区ꎻ３－钻孔验证采空区ꎻ４－验证钻孔１－ｖｅｉｎｓꎻ２－ｉｎｆｅｒｒｅｄｇｏａｆꎻ３－ｄｒｉｌｌｉｎｇｖｅｒｉｆｉｅｄｇｏａｆꎻ４－ｖｅｒｉｆｙｄｒｉｌｌｉｎｇ在１线５５０~６３０ｍ㊁６４０~６９０ｍ和３３０~４１０ｍ处分布有ｃｋ１㊁ｃｋ２和ｃｋ３高阻和低阻异常ꎬ电阻率在３２００Ω ｍ~６４００Ω ｍ和４００Ω ｍ~８００Ω ｍ之间变化ꎬｃｋ１㊁ｃｋ２㊁ｃｋ３异常地表附近出露有多条矿化带并ｃｋ２南１００ｍ有已知民采井ＪＫ１３Ｃꎬ据此推测ｃｋ２异常为已知采空区引起ꎬｃｋ３为破碎带或含水空区引起ꎬｃｋ１异常为隐伏采空区引起ꎮｃｋ１㊁ｃｋ２异常顶部埋深一般２０~４０ｍ左右ꎬｃｋ３异常顶部埋深在５０ｍ左右ꎮ为了查验推断的电法异常ꎬ施工验证钻孔一个ꎬ孔深２６０ｍꎬ钻孔倾角７０ʎꎬ倾向北西ꎮ经钻探工程验证共揭露断层一处ꎬ采空区三处ꎬ在揭露的断层和采空区处实测拟断面等值线和反演断面等值线均具扭曲和收缩特点ꎬ所处深度与物探解译３１８地质与勘探２０１９年深度基本一致ꎮ其中断层位于３９ ５~４３ ６ｍꎬ岩性为弱褐铁矿化碎裂花岗岩ꎬ岩心采取率为８９％ꎬ该断层处在１线ｃｋ２已知采空区相对高阻异常上ꎬ电阻率值在３２００Ω ｍ~６４００Ω ｍꎬ钻探揭露断层深度与推断ｃｋ２异常顶部埋深在２０~４０ｍ之间相近ꎮ在１２０ ３~１２２ ３ｍ为采空区ꎬ岩性主要为绿泥石化和高岭土化ꎬ岩芯破碎ꎬ蚀变较强ꎬ岩心平均采取率为２５％ꎻ在１４０~１４１ ８ｍ为采空区ꎬ岩性为花岗岩ꎬ岩心破碎ꎬ呈棱角状ꎬ以上两处采空区处在１线ｃｋ１已知采空区相对高阻异常上ꎬ电阻率值在３２００Ω ｍ~６４００Ω ｍꎬ钻探揭露采空区与推断ｃｋ１异常相近ꎮ通过钻探验证工作得出以下结论:物探推断的采空区位置与钻探验证见采空区位置基本一致ꎬ证实了高密度直流电阻率法在本区采空区调查中方法有效ꎬ在金矿采空区调查中高密度直流电阻率法是一种可靠的技术ꎮ３ ２㊀瞬变电磁等值反磁通试验剖面完成瞬变电磁法剖面１条ꎬ剖面沿北东方向展布ꎬ编号为２号测线ꎮ剖面成果解释如下(图６):图６㊀２线地质物探综合剖面图Ｆｉｇ.６㊀ＧｅｏｌｏｇｉｃａｌｇｅｏｐｈｙｓｉｃａｌｃｏｍｐｒｅｈｅｎｓｉｖｅｐｒｏｆｉｌｅａｌｏｎｇｔｈｅｓｕｒｖｅｙｌｉｎｅＮｏ.２１－推断采空区ꎻ２－钻孔验证采空区ꎻ３－验证钻孔ꎻ４－强风化层界面１－ｉｎｆｅｒｒｅｄｇｏａｆꎻ２－ｄｒｉｌｌｉｎｇｖｅｒｉｆｉｅｄｇｏａｆꎻ３－ｖｅｒｉｆｙｄｒｉｌｌｉｎｇꎻ４－ｓｔｒｏｎｇｌｙｗｅａｔｈｅｒｉｎｇｌａｙｅｒｉｎｔｅｒｆａｃｅ共推断有４个采空影响范围或破碎岩体ꎬ其中心埋深在９０~１３８ｍ变化ꎬ电阻率在３００Ω ｍ~１１００Ω ｍ变化ꎬ４个异常影响范围大致分布在１２线的３７５~４３６ｍ㊁５８０~６４５ｍ㊁７６５~８２０ｍ㊁９４５~１０３５ｍꎬ表现为低阻异常ꎬ其中在地表４００~４４０ｍ为已知竖井大致范围ꎮ在剖面１０６０ｍ处施工了验证钻孔一个ꎬ深２４０ｍꎬ倾向北西ꎬ倾角７０ʎꎮ钻探工程共揭露采空区三处ꎬ在揭露的采空区处反演断面等值线具扭曲和收缩特点ꎬ所处深度与物探解译深度基本一致ꎮ其中在１６２ ６~１６４ ５ｍ为采空区ꎬ岩性主要为黑云母斜长花岗岩ꎻ１９７ ７~１９９ ４ｍ为采空区ꎬ岩性主要为伟晶岩ꎻ２０４ ２~２０６ １ｍ为采空区ꎬ岩性主要为含黄铁矿黑云母斜长花岗岩ꎮ岩心平均采取率为３０％ꎮ以上三处采空区均处在２线相对低阻异常上ꎬ电阻率值在１０００Ω ｍ左右ꎮ钻探工程揭露的采空区在埋深和位置上与推断４１８第３期张㊀扬等:多物探方法在胶东某金矿采空区勘查中的应用研究的采空区异常大致相同ꎬ通过钻探验证工作可知:物探推断采空区位置与钻探验证较为符合ꎬ证实等值反磁通瞬变电磁法在本区采空区调查中方法有效ꎬ并可以降低工频干扰ꎬ在本区采空区勘探中该法有较好的应用效果ꎮ３ ３㊀地下采空区充水分析由于试验区矿段－４９０ｍ标高以上矿体基本采空ꎬ矿山对采空区及时进行了尾砂胶结充填和高水固结充填ꎬ充填回水直接随采场回水井汇入排水系统ꎮ目前影响矿坑充水的直接因素是断裂破碎带ꎬ间接因素为大气降水渗透ꎬ现两矿段巷道已打通ꎬ－５００ｍ以上地下水已疏干ꎮ在本试验区范围内ꎬ深度２００ｍ以浅的地下采空区充填地下水的可能性不大ꎬ在此区域之外的民采空区有可能充填地下水ꎮ因此ꎬ对于试验区矿段范围以外的低阻区要关注它的含水性ꎮ４㊀结论依据不同地电条件分别选择不同的物探方法ꎬ钻探验证达到了预期的目的ꎬ高密度直流电阻率法和等值反磁通瞬变电磁法在本区采空区调查中方法有效ꎬ同时等值反磁通瞬变电磁法可以降低工频干扰ꎮ两种方法在本区采空区调查中均有较好的应用效果ꎮ通过对第四系覆盖和蓄水影响的综合分析ꎬ使得成果资料的解释更合理化ꎬ经验证后也证明了推断的准确性ꎬ这对采空区物探的资料解释是一个很好的启发ꎮ采空区一般多为场地较复杂的地区ꎬ异常的形成原因及影响因素较多ꎬ由于边界条件的限制不足ꎬ已知信息有限的情况下ꎬ地球物理勘查的成果本身就具有多解性ꎬ这对资料处理解释提出了更高的要求ꎮ[Ｒｅｆｅｒｅｎｃｅｓ]ＤｅｎｇＣｈａｏｗｅｎ.２００７.Ｔｈｅｐｒｉｎｃｉｐｌｅａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｅｌｅｃｔｒｉｃａｌｉｍａｇｉｎｇｓｕｒｖｅｙｓ[Ｊ].ＪｏｕｒｎａｌｏｆＳｈａｏｇｕａｎＣｏｌｌｅｇｅ－ＴｈｅＮａｔｕｒａｌＳｃｉｅｎｃｅｓꎬ２８(６):６５－６７(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＤｏｎｇＨａｏｂｉｎꎬＷａｎｇＣｈｕａｎｌｅｉ.２００３.Ｄｅｖｅｌｏｐｍｅｎｔａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｈｉｇｈ－ｒｅｓｉｓｔｉｖｉｔｙｅｌｅｃｔｒｉｃａｌｔｅｃｈｎｉｑｕｅ[Ｊ].ＥａｒｔｈＳｃｉｅｎｃｅＦｒｏｎｔｉｅｒｓꎬ１０(１):１７１－１７６(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＤｏｎｇＪｉａｎꎬＣｈｅｎＬｅｉꎬＺｈａｎｇＧｕｉｌｉ.２０１８.ＣｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｇｅｏｐｈｙｓｉｃａｌｆｉｅｌｄｓａｒｏｕｎｄｇｏｌｄｄｅｐｏｓｉｔｓｉｎｔｈｅＸｉａｄｉａｎｒｅｇｉｏｎꎬＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ５４(１):１３８－１４７(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＦａｎＬｉａｎｇꎬＱｉａｎＲｏｎｇｙｉ.２０１１.Ａｐｐｌｉｃａｔｉｏｎｏｆｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏ－ｍａｇｎｅｔ￣ｉｃｍｅｔｈｏｄｔｏｃｏａｌｍｉｎｅｇｏｂｓ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅ￣ｏｐｈｙｓｉｃｓꎬ８(１):２９－３３(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＧｅＳｈｕａｎｇｃｈｅｎｇꎬＬｉＸｉａｏｐｉｎｇꎬＳｈａｏＣｈａｎｇｙｕｎꎬＣｈｅｎＪｕｎꎬＣｈｅｎＸｉｎｇ.２００８.Ａｐｐｌｉｃａｔｉｏｎｏｆｓｅｉｓｍｉｃｒｅｆｒａｃｔｉｏｎａｎｄｒｅｓｉｓｔｉｖｉｔｙｆｏｒｅｘｐｌｏｒａｔｉｏｎｏｆｒｅｓｅｒｖｏｉｒｄａｍｓｉｔｅ[Ｊ].ＰｒｏｇｒｅｓｓｉｎＧｅｏｐｈｙｓｉｃｓꎬ２３(４):１２９９－１３０３(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＧｕｏＸｉｕｊｕｎꎬＪｉａＹｏｎｇｇａｎｇꎬＨａｎｇＸｉａｏｙｕꎬＮｉｕＪｉａｎｊｕｎꎬ２００４.Ａｐｐｌｉｃａｔｉｏｎｏｆｍｕｌｔｉ－ｅｌｅｃｔｏｄｅｓｅｌｅｃｔｒｉｃａｌｍｅｔｈｏｄｔｏｄｅｔｅｃｔｉｏｎｏｆｓｌｉｄｅ－ｆａｃｅｐｏ￣ｓｉｔｉｏｎ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＲｏｃｋＭｅｃｈａｎｉｅｓａｎｄＥｎｇｉｎｅｅｒꎬ２３(１０):１６６２－１６６９(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＧｕｏＹｏｕｇａｎｇꎬＬｉｕＦｅｎｇꎬＷａｎｇＢｉｎｗｕ.２０１１.ＴｈｅｅｘｐｌｏｒａｔｉｏｎｏｆＴＥＭｉｎｃｏａｌｍｉｎｅｄ－ｏｕｔａｒｅａ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｐｈｙｓ￣ｉｃｓꎬ７(２):１５１－１５４(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＨａｎＪｉｎｇｕｏꎬＣｈｅｎＱｉꎬＺｅｎｇＹｏｕｑｉａｎｇꎬＤｕＬｉｍｉｎｇꎬＷａｎｇＬｕꎬＺｈａｎｇＹａｎｇ.２０１８.Ａｐｐｌｉｃａｔｉｏｎｏｆｍｉｎｅｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄａｄｖａｎｃｅｄｐｒｅｄｉｃｔｉｏｎｏｆａｇｏｌｄｍｉｎｅ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ５４(２):３４４－３４７(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＨｕＢｏｗｅｎꎬＺｈａｎｇＦａｗａｎｇꎬＣｈｅｎＬｉꎬＺｈａｎｇＪｉｎ.２０１５.Ｓｔａｔｕｓａｎｄｐｒｏｓｐｅｃｔｅｖａｌｕａｔｉｏｎｏｆｍｉｎｅｇｅｏｌｏｇｉｃａｌｅｎｖｉｒｏｎｍｅｎｔｒｅｓｅａｒｃｈ[Ｊ].ＥａｒｔｈａｎｄＥｎｖｉｒｏｎｍｅｎｔꎬ４３(３):３７５－３７８(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＪｉａｎｇＹｕｌｅꎬＫａｎｇＷａｎｆｕꎬＺｈａｎｇＮａｎꎬＺｈａｎｇＺｈａｏｘｉａ.２００７.Ａｐｐｌｉｃａｔｉｏｎｏｆｈｉｇｈｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅｔｈｏｄｔｏｋａｒｓｔｅｘｐｌｏｒａｔｉｏｎ[Ｊ].ＪｏｕｒｎａｌｏｆＣｈｅｎｇｄｕＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ(Ｓｃｉｅｎｃｅ＆ＴｅｃｈｎｏｌｏｇｙＥｄｉｔｉｏｎ)ꎬ３４(４):４５２－４５５(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＪｉｎＺｈｏｎｇｇｕｏꎬＺｏｕＬｉｎ.２００２.Ａｐｐｌｉｃａｔｉｏｎｏｆｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏ－ｍａｇｎｅｔｉｃｍｅｔｈｏｄｓｔｏｐｒｏｓｐｅｃｔｉｎｇＭａｏａｏｃｈａｎｇＰｂ－ＺｎｏｒｅｄｉｓｔｒｉｃｔｉｎｎｏｒｔｈｗｅｓｔＧｕｉｚｈｏｕ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ(６):４８－５０(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＬｉＷｅｎｙａｏ.２０１０.Ａｃｃｕｒａｔｅｆｏｒｍｕｌａｏｆｓｅｌｆ－ｉｎｄｕｃｔａｎｃｅｃｏｅｆｆｉｃｉｅｎｔｏｆｒｅｃｔａｎｇｕｌａｒｃｏｉｌｉｎｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄｓ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ４６(１):１６０－１６４(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＬｉｕＧｕｏｘｉｎｇ.２００５.Ｐｒｉｎｃｉｐｌｅｓａｎｄｍｅｔｈｏｄｓｏｆｅｌｅｃｔｒｉｃｐｒｏｓｐｅｃｔｉｎｇ[Ｍ].Ｂｅｉｊｉｎｇ:ＧｅｏｌｏｇｉｃａｌＰｕｂｌｉｓｈｉｎｇＨｏｕｓｅ:６３－６９(ｉｎＣｈｉｎｅｓｅ).ＬüＨｕｉｊｉｎꎬＬｉｕＳｈａｏｈｕａꎬＬｉｕＢｏｇｅｎ.２００５.Ａｐｐｌｉｃａｔｉｏｎｏｆｒｅｓｉｓｔｉｖｉｔｙｔｏｍｏｇｒａｐｈｙｓｕｒｖｅｙｍｅｔｈｏｄｉｎｄｅｔｅｃｔｉｎｇｇｒｏｕｎｄｓｕｂｓｉｄｅｎｃｅ[Ｊ].Ｐｒｏ￣ｇｒｅｓｓｉｎＧｅｏｐｈｙｓｉｃｓꎬ２０(２):３８１－３８６(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＭａＤｅｘｉꎬＹｕＡｉｊｕｎꎬＧｅＬｉａｎｇｓｈｅｎｇꎬＦａｎＪｕｎｊｉｅꎬＺｈａｎｇＸｕｅｊｕｎꎬＷａｎｇＣｈｕｎｓｈｅｎｇ.２００８.Ａｐｐｌｉｃａｔｉｏｎｏｆｍｕｌｔｉ－ｅｌｅｃｔｒｏｄｅｅｌｅｃｔｒｉｃａｌｍｅｔｈｏｄｉｎｇｏｌｄｏｒｅｐｒｏｓｐｅｃｔｉｎｇ[Ｊ].ＧｅｏｌｏｇｙａｎｄＰｒｏｓｐｅｃｔｉｎｇꎬ４４(３):６５－６９(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＱｉＭｉｎꎬＺｈａｎｇＢａｏｌｉｎꎬＬｉａｎｇＧｕａｎｇｈｅꎬＣｈｅｎＹｏｕｍｉｎｇꎬＧｕａｎＧｕａｎｇ.２００６.Ｈｉｇｈ－ｒｅｓｏｌｕｔｉｏｎｐｒｅｄｉｃｔｉｏｎｏｆｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎｃｈａｒａｃｔｅｒｉｓ￣ｔｉｃｓｏｆｃｏｍｐｌｉｃａｔｅｄｕｎｄｅｒｇｒｏｕｎｄｃａｖｉｔｉｅｓ－Ｐｒｅｌｉｍｉｎａｒｙａｐｐｌｉｃａｔｉｏｎｏｆｈｉｇｈ－ｄｅｎｓｉｔｙｅｌｅｃｔｒｉｃａｌｔｅｃｈｎｉｑｕｅｉｎａｎａｒｅａｏｆＹａｎｇｑｕａｎꎬＳｈａｎｘｉ[Ｊ].ＰｒｏｇｒｅｓｓｉｎＧｅｏｐｈｙｓｉｃｓꎬ２１(１):２５６－２６２(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＱｕＳｈａｏｆｅｉ.２００８.Ｇｅｏｐｈｙｓｉｃａｌａｎｄｇｅｏｃｈｅｍｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄｇｅｏｌ￣ｏｇｙｏｆＤａｍｏｑｕｊｉａｇｏｌｄｄｅｐｏｓｉｔｉｎＬｏｎｇｋｏｕｏｆＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＬａｎｄａｎｄＲｅｓｏｕｒｃｅｓｉｎＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅꎬ２４(３):２６－３０(ｉｎＣｈｉ￣ｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＳｏｎｇＧｕｏｚｈｅｎｇꎬＬｉＳｈａｎꎬＹａｎＣｈｕｎｍｉｎｇꎬＢａｏＺｈｏｎｇｙｉꎬＦａｎＪｉａｍｅｎｇꎬＬｉｕＣａｉｊｉｅꎬＬｉｕＧｕｏｄｏｎｇꎬＧｕｏＺｈｉｆｅｎｇ.２０１８.ＧｅｏｌｏｇｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄｐｒｏｓｐｅｃｔｉｎｇｄｉｒｅｃｔｉｏｎｏｆＮｏ.ＩｍａｉｎｏｒｅｂｏｄｙｉｎｔｈｅＪｉａｏｊｉａｇｏｌｄ５１８地质与勘探２０１９年ｄｅｐｏｓｉｔ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ５４(２):２１９－２２９(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＳｏｎｇＭｉｎｇｃｈｕｎꎬＬｉＳａｎｚｈｏｎｇꎬＹｉＰｉｈｏｕꎬＣｕｉＳｕｎｘｕｅꎬＸｕＪｕｎｘｉａｎｇꎬＬüＧｕｘｉａｎꎬＳｏｎｇＹｉｎｇｘｉꎬＪｉａｎｇＨｏｎｇｌｉꎬＺｈｏｕＭｉｎｇｌｉｎｇꎬＺｈａｎｇＰｉｎｊｉａｎꎬＨｕａｎｇＴａｉｌｉｎｇꎬＬｉｕＣｈａｎｇｃｈｕｎꎬＬｉｕＤｉａｎｈａｏ.２０１４.ＣｌａｓｓｉｆｉｃａｔｉｏｎａｎｄｍｅｔａｌｌｏｇｅｎｉｃｔｈｅｏｒｙｏｆｔｈｅＪｉａｏｊｉａ－ＳｔｈｌｅｇｏｌｄｄｅｐｏｓｉｔｉｎＪｉａｏｄｏｎｇＰｅｎｉｎｓｕｌａꎬＣｈｉｎａ[Ｊ].ＪｏｕｒｎａｌｏｆＪｉｌｉｎＵｎｉｖｅｒｓｉｔｙ(ＥａｒｔｈＳｃｉｅｎｃｅＥｄｉ￣ｔｉｏｎ)ꎬ４４(１):８７－１０４(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＳｏｎｇＭｉｎｇｃｈｕｎꎬＹａｎｇＣｈｅｎｇｈａｉꎬＪｉａｏＸｉｕｍｅｉ.２００７.Ｐｒｏｂｅｉｎｔｏｔｈｅｄｉｖｉ￣ｓｏｎｏｆｍｅｔａｌｌｏｎｇｅｎｉｃｐｒｏｖｉｎｃｅｓａｎｄｔａｒｇｅｔａｒｅａｓｏｆｇｏｌｄｅｘｐｌｏｒａｔｉｏｎｉｎＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＣｏｎｔｒｉｂｕｔｉｏｎｓｔｏＧｅｏｌｏｇｙａｎｄＭｉｎｅｒａｌＲｅ￣ｓｏｕｒｃｅｓＲｅｓｅａｒｃｈꎬ２２(４):２４８－２５２ꎬ２５８(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＳｕｎＨｕａｉｆｅｎｇꎬＬｉＸｉｕꎬＬｕＸｕｓｈａｎꎬＬｉＳｈｕｃａｉꎬＲｅｎＢａｏｈｏｎｇ.２０１６.Ｔｒａｎ￣ｓｉｅｎｔｅｌｅｃ－ｔｒｏｍａｇｎｅｔｉｃｒｅｓｐｏｎｓｅｓｉｎｔｕｎｎｅｌｓｗｉｔｈｓｔｒｏｎｇｉｎｔｅｒｆｅｒｅｎｃｅｓａｎｄｔｈｅｃｏｒｒｅｃｔｉｎｇｍｅｔｈｏｄ:ＡＴＢＭｅｘａｍｐｌｅ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＧｅｏｐｈｙｓｉｃｓꎬ５９(１２):４７２０－４７３２(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＷａｎｇＱｉｊｕｎꎬＨｕＹａｌｉｎꎬＤｕＸｉｎｇｆｅｎｇꎬＬｉＳｈｉｊｕｎ.２００８.Ａｐｐｌｉｃａｔｉｏｎｏｆｔｈｅｈｉｇｈ－ｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅｔｈｏｄｉｎｐｒｏｊｅｃｔｅｘｐｌｏｒａｔｉｏｎ[Ｊ].ＰｒｏｇｒｅｓｓｉｎＧｅｏｐｈｙｓｉｃｓꎬ２３(４):１２９９－１３０３(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＷａｎｇＨｕａｆｅｎｇꎬＬｉｕＲｏｎｇｑｕａｎꎬＺｈｅｎｇＱｉａｎｇꎬＣｈｅｎＹｉｎｇｊｕｎꎬＧｕｏＤｏｎｇꎬＬｉｕＪｕｎꎬＹａｎＪｉａｙｏｎｇ.２０１３.Ａｐｐｌｉｃａｔｉｏｎｏｆｃｏｍｐｒｅｈｅｎｓｉｖｅｇｅｏｐｈｙｓ￣ｉｃａｌｐｒｏｓｐｅｃｔｉｎｇｍｅｔｈｏｄｓｔｏｇｏａｆｏｆｍｅｔａｌｏｒｅｍｉｎｅｓ:ＡｃａｓｅｓｔｕｄｙｏｆｔｈｅＷａｎｇｅｒｓｈａｎｍｉｎｅｇｏａｆꎬＪｉａｏｊｉａｇｏｌｄｍｉｎｅꎬＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ４９(３):４９６－５０４(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＷｕＪｕｎｊｉｅ.２０１５.ＡｐｐｌｉｃａｔｉｏｎｏｆＴＥＭｉｎｔｈｅＴｉａｎｙｕＣｕ－ＮｉｄｅｐｏｓｉｔｉｎＥａｓｔＴｉａｎＳｈａｎꎬＸｉｎｊｉａｎｇ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎꎬ５１(５):９７０－９７６(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＹａｎｇＺｈｅｎｗｅｉꎬＹａｎＪｉａｙｏｎｇꎬＬｉｕＹａｎꎬＷａｎｇＨｕａｆｅｎｇ.２０１２.Ｒｅｓｅａｒｃｈｐｒｏ￣ｇｒｅｓｓｅｓｏｆｔｈｅｈｉｇｈ－ｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅｔｈｏｄ[Ｊ].ＧｅｏｌｏｇｙａｎｄＥｘ￣ｐｌｏｒａｔｉｏｎꎬ４８(５):９６９－９７８(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＺｈａｎｇＱｉａｎｇꎬＸｕＲｕｆｅｎｇꎬＸｕＦａｎｇꎬＷａｎｇＪｉｕｈｕａꎬＸｉａｏＢｉｎｇｊｉａｎꎬＬｉｕＴｏｎｇꎬＹａｎｇＸｕｅｓｈｅｎｇ.２００６.Ｇｅｏｌｏｇｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｇｏｌｄｄｅ￣ｐｏｓｉｔｓａｎｄｏｒｅ－ｓｅａｒｃｈｉｎｇｄｉｒｅｃｔｉｏｎｉｎＤａｌｉｕｈａｎｇａｒｅａꎬＰｅｎｌａｉＣｉｔｙꎬＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＣｏｎｔｒｉｂｕｔｉｏｎｓｔｏＧｅｏｌｏｇｙａｎｄＭｉｎｅｒａｌＲｅ￣ｓｏｕｒｃｅｓＲｅｓｅａｒｃｈꎬ２１(４):２５８－２６１(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＺｈａｎｇＹｏｎｇｓｈｕａｎｇꎬＳｕｎＬｕꎬＹｉｎＸｉｕｌａｎꎬＭｅｎｇＨｕｉ.２０１７.ＰｒｏｇｒｅｓｓａｎｄｐｒｏｓｐｅｃｔｏｆｒｅｓｅａｒｃｈｏｎｅｎｖｉｒｏｎｍｅｎｔａｌｇｅｏｌｏｇｙｏｆＣｈｉｎａ:Ａｒｅｖｉｅｗ[Ｊ].ＧｅｏｌｏｇｙｉｎＣｈｉｎａꎬ４４(５):９０１－９１２(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).ＺｈｏｕＰｉｎｇꎬＳｈｉＪｕｎｆａ.２００７.Ｎｅｗｐｒｏｇｒｅｓｓｅｓａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄ(ＴＥＭ)ｉｎｄｅｅｐ－ｃｏｎｃｅａｌｅｄｏｒｅｆｉｎｄｉｎｇ[Ｊ].ＧｅｏｌｏｇｙａｎｄＰｒｏｓｐｅｃｔｉｎｇꎬ(６):６３－６９(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＺｈｏｕＹａｎｚｈａｎｇꎬＣｈｉＢａｏｍｉｎｇꎬＬｉｕＺｈｏｎｇｐｅｉ.２００８.Ｍｏｄｅｏｆｓｔｒｕｃｔｕｒａｌｃｏｎｔｒｏｌａｂｏｕｔｇｒｏｕｎｄｗａｔｅｒｐｏｕｒ－ｉｎｍｅｃｈａｎｉｓｍｉｎＸｉａｄｉａｎｇｏｌｄｄｅ￣ｐｏｓｉｔｓｉｎＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ[Ｊ].ＪｏｕｒｎａｌｏｆＪｉｌｉｎＵｎｖｉｅｒｓｉｔｙ(ＥａｒｔｈＳｃｉｅｎｃｅＥｄｉｔｉｏｎ)ꎬ３８(２):２５５－２６０(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂ￣ｓｔｒａｃｔ).ＺｈｕＹａｊｕｎꎬＷａｎｇＹａｎｘｉｎ.２０１２.Ｔｈｅａｐｐｌｉｃａｔｉｏｎｏｆｈｉｇｈｄｅｎｓｉｔｙｒｅｓｉｓ￣ｔｉｖｉｔｙａｎｄＴＥＭｍｅｔｈｏｄｔｏｕｎｄｅｒｇｒｏｕｎｄｋａｒｓｔｄｅｔｅｃｔｉｏｎ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｐｈｙｓｉｃｓꎬ９(６):７３８－７４２(ｉｎＣｈｉｎｅｓｅｗｉｔｈＥｎｇｌｉｓｈａｂｓｔｒａｃｔ).[附中文参考文献]邓超文.２００７.高密度电法的原理及工程应用[Ｊ].韶关学院学报－自然科学版ꎬ２８(６):６５－６７.董浩斌ꎬ王传雷.２００３.高密度电法的发展与应用[Ｊ].地学前缘ꎬ１０(１):１７１－１７６.董㊀健ꎬ陈㊀磊ꎬ张贵丽.２０１８.山东夏甸地区金矿床地球物理场特征[Ｊ].地质与勘探ꎬ５４(１):１３８－１４７.范㊀亮ꎬ荣㊀毅.２０１１.瞬变电磁法在煤矿采空区的应用研究[Ｊ].工程地球物理学报ꎬ８(１):２９－３３.葛双成ꎬ李小平ꎬ邵长云ꎬ陈㊀军ꎬ陈㊀星.２００８.地震折射和电阻率法在水库坝址勘察中的应用[Ｊ].地球物理学进展ꎬ２３(４):１２９９－１３０３.㊀郭秀军ꎬ贾永刚ꎬ黄潇雨ꎬ牛建军.２００４.利用高密度电阻率法确定滑坡面研究[Ｊ].岩石力学与工程学报ꎬ２３(１０):１６６２－１６６９.郭有刚ꎬ刘㊀峰ꎬ王斌武.２０１１.瞬变电磁法在煤矿采空区勘探中的应用[Ｊ].工程地球物理学报ꎬ７(２):１５１－１５４韩进国ꎬ陈㊀琦ꎬ增友强ꎬ杜利明ꎬ王㊀璐ꎬ张㊀扬.２０１８.井巷瞬变电磁法在某金矿超前预报中的应用[Ｊ].地质与勘探ꎬ５４(２):３４４－３４７.胡博文ꎬ张发旺ꎬ陈㊀立ꎬ张㊀瑾.２０１５.我国矿山地质环境评价方法研究现状及展望[Ｊ].地球与环境ꎬ４３(３):３７５－３７８.江玉乐ꎬ康万福ꎬ张㊀楠ꎬ张朝霞.２００７.高密度电法在岩溶勘察中的应用[Ｊ].成都理工大学学报(自科版)ꎬ３４(４):４５２－４５５.金中国ꎬ邹㊀林.２００２.瞬变电磁法在黔西北猫猫厂铅锌矿区找矿中的应用[Ｊ].地质与勘探ꎬ(６):４８－５０.李文尧.２０１０.瞬变电磁法矩形线圈自感的精确表达式[Ｊ].地质与勘探ꎬ４６(１):１６０－１６４.刘国兴.２００５.电法勘探原理与方法[Ｍ].北京:地质出版社:６３－６９.吕惠进ꎬ刘少华ꎬ刘伯根.２００５.高密度电阻率法在地面塌陷调查中的应用[Ｊ].地球物理学进展ꎬ２０(２):３８１－３８６.马德锡ꎬ于爱军ꎬ葛良胜ꎬ范俊杰ꎬ张学军ꎬ王春生.２００８.高密度电法在金矿勘查中的应用[Ｊ].地质与勘探ꎬ４４(３):６５－６９.祁㊀民ꎬ张宝林ꎬ梁光河ꎬ陈友明ꎬ管㊀刚.２００６.高分辨率预测地下复杂采空区的空间分布特征－高密度电法在山西阳泉某复杂采空区中的初步应用研究[Ｊ].地球物理学进展ꎬ２１(１):２５６－２６２.㊀曲少飞.２００８.山东龙口大磨曲家金矿地质及地球物理地球化学特征[Ｊ].山东国土资源ꎬ２４(３):２６－３０.宋国政ꎬ李㊀山ꎬ闫春明ꎬ鲍中义ꎬ范家盟ꎬ刘彩杰ꎬ刘国栋ꎬ郭志峰.２０１８.焦家金矿田Ⅰ号主矿体地质特征及找矿方向[Ｊ].地质与勘探ꎬ５４(２):２１９－２２９.宋明春ꎬ李三忠ꎬ伊丕厚ꎬ崔书学ꎬ徐军祥ꎬ吕古贤ꎬ宋英昕ꎬ姜洪利ꎬ周明岭ꎬ张丕建ꎬ黄太岭ꎬ刘长春ꎬ刘殿浩.２０１４.中国胶东焦家式金矿类型及成矿理论[Ｊ].吉林大学学报(地球科学版)ꎬ４４(１):８７－１０４.宋明春ꎬ杨承海ꎬ焦秀美.２００７.山东省金矿成矿区带划分及找矿方６１８第３期张㊀扬等:多物探方法在胶东某金矿采空区勘查中的应用研究向探讨[Ｊ].地质找矿论丛ꎬ２２(４):２４８－２５２ꎬ２５８.孙怀凤ꎬ李㊀貅ꎬ卢绪山ꎬ李术才ꎬ任宝宏.２０１６.隧道强干扰环境瞬变电磁响应规律与校正方法:以ＴＢＭ为例[Ｊ].地球物理学报ꎬ５９(１２):４７２０－４７３２.王启军ꎬ胡延林ꎬ都兴锋ꎬ李世军.２００９.高密度电阻率法在工程勘查中的应用[Ｊ].地球物理学进展ꎬ２４(１):３３５－３３９.王华锋ꎬ刘荣泉ꎬ郑㊀强ꎬ陈应军ꎬ郭㊀冬ꎬ刘㊀俊ꎬ严加永.２０１３.综合物探在金属矿采空区中的应用－以焦家金矿望儿山采空区为例[Ｊ].地质与勘探ꎬ４９(３):４９６－５０４.武军杰.２０１５.瞬变电磁法在新疆东天山天宇铜镍矿区的应用效果[Ｊ].地质与勘探ꎬ５１(５):９７０－９７６.杨振威ꎬ严加永ꎬ刘㊀彦ꎬ王华锋.２０１２.高密度电阻率法研究进展[Ｊ].地质与勘探ꎬ４８(５):９６９－９７８.张㊀强ꎬ徐汝峰ꎬ许㊀方ꎬ王久华ꎬ肖丙建ꎬ刘㊀同ꎬ杨学生.２００６.山东省蓬莱大柳行地区金矿地质特征及找矿方向[Ｊ].地质找矿论丛ꎬꎬ２１(４):２５８－２６１.张永双ꎬ孙㊀璐ꎬ殷秀兰ꎬ孟㊀晖.２０１７.中国环境地质研究主要进展与展望[Ｊ].中国地质ꎬ４４(５):９０１－９１２.周㊀平ꎬ施俊法.２００７.瞬变电磁法(ＴＥＭ)新进展及其在寻找深部隐伏矿中的应用[Ｊ].地质与勘探ꎬ(６):６３－６９.周彦章ꎬ迟宝明ꎬ刘中培.２００８.山东夏甸金矿床充水机理构造控制模式[Ｊ].吉林大学学报(地球科学版)ꎬ３８(２):２５５－２６０.朱亚军ꎬ王艳新.２０１２.高密度电法和瞬变电磁法在地下岩溶探测中的综合应用[Ｊ].工程地球物理学报ꎬ９(６):７３８－７４２.ＡｐｐｌｉｃａｔｉｏｎｏｆＭｕｌｔｉ－ＧｅｏｐｈｙｓｉｃａｌＭｅｔｈｏｄｓｉｎＥｘｐｌｏｒａｔｉｏｎｏｆａＷｏｒｋｅｄ－ＯｕｔＡｒｅａｏｆＧｏｌｄＭｉｎｉｎｇｉｎｔｈｅＪｉａｏｄｏｎｇＲｅｇｉｏｎＺＨＡＮＧＹａｎｇ１ꎬＷＡＮＧＬｕ１ꎬＬＩＵＺｈｉｈｕｉ２ꎬＬＩＵＳｈｅｎｇｔａｉ１ꎬＭＥＩＺｈｅｎｈｕａ１ꎬＤＵＬｉｍｉｎｇ１ꎬＬＩＮＸｉａｏｈｕｉ３ꎬＺＨＡＮＧＪｉａｎｂｏ４(１.ＳｈａｎｄｏｎｇＺｈｅｎｇｙｕａｎＧｅｏｌｏｇｉｃａｌＥｘｐｌｏｒａｔｉｏｎＩｎｓｔｉｔｕｔｅꎬＣｈｉｎａＭｅｔａｌｌｕｒｇｉｃａｌＧｅｏｌｏｇｙＢｕｒｅａｕꎬＪｉｎａｎꎬＳｈａｎｄｏｎｇ㊀２５００００ꎻ２.ＣｈｉｎａＭｅｔａｌｌｕｒｇｉｃａｌＧｅｏｌｏｇｙＢｕｒｅａｕꎬＢｅｉｊｉｎｇ㊀１０００２６ꎻ３.ＳｅｃｏｎｄＢｕｒｅａｕｏｆＣｈｉｎａＭｅｔａｌｌｕｒｇｉｃａｌＧｅｏｌｏｇｙꎬ㊀㊀㊀ＦｕｚｈｏｕꎬＦｕｊｉａｎ㊀３５００００ꎻ４.ＳｈａｎｄｏｎｇＺｈｅｎｇｙｕａｎＣｏｎｓｔｒｕｃｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇＣｏ.ＬｔｄꎬＪｉｎａｎꎬＳｈａｎｄｏｎｇ㊀２５００００)Ａｂｓｔｒａｃｔ:Ｉｎｒｅｃｅｎｔｙｅａｒｓꎬｃｏｌｌａｐｓｅａｃｃｉｄｅｎｔｓｏｆｔｅｎｏｃｃｕｒｉｎｗｏｒｋｅｄ－ｏｕｔａｒｅａｓｏｆｍｉｎｅｓꎬｉｎｄｉｃａｔｉｎｇｇｒｅａｔｉｍｐｏｒｔａｎｃｅｏｆｅｘｐｌｏｒａｔｉｏｎｔｏｔｈｅｓｅａｒｅａｓ.Ｔｈｉｓｐａｐｅｒｐｒｅｓｅｎｔｓａｎａｐｐｌｉｃａｔｉｏｎｏｆｍｕｌｔｉ－ｇｅｏｐｈｙｓｉｃａｌｍｅｔｈｏｄｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｔｈｅｗｏｒｋｅｄ－ｏｕｔａｒｅａｉｎａｇｏｌｄｍｉｎｅｏｆｔｈｅＪｉａｏｄｏｎｇｒｅｇｉｏｎ.Ｕｎｄｅｒｔｈｅｃｏｎｄｉｔｉｏｎｏｆｓｅｒｉｏｕｓｄｉｓｔｕｒｂａｎｃｅｗｉｔｈｉｎｔｈｅｍｉｎｉｎｇａｒｅａꎬｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｓｕｒｖｅｙｗｉｔｈｍａｎｙｄｅｖｉｃｅｔｙｐｅｓꎬｄｃｈｉｇｈｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅａｓｕｒｅｍｅｎｔａｎｄｓｈａｌｌｏｗｓｅｉｓｍｉｃｅｘｐｌｏｒａｔｉｏｎｗｅｒｅｃｏｎｄｕｃｔｅｄｓｕｃｃｅｓｓｉｖｅｌｙａｓｃｈｏｉｃｅｅｘｐｅｒｉｍｅｎｔｓ.Ｆｉｎａｌｌｙｔｈｅｈｉｇｈｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅｔｈｏｄａｎｄｔｈｅｅｑｕｉｖａｌｅｎｔｉｎｖｅｒｓｅｆｌｕｘｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄ(ＯＣＴＥＭ)ｗｅｒｅｃｈｏｓｅｎｉｎｔｈｅｔｅｓｔｓｅｃｔｉｏｎ.Ｄｒｉｌｌｉｎｇｖｅｒｉｆｉｃａｔｉｏｎｓｈｏｗｓｔｈｅｅｆｆｅｃｔｉｖｅｎｅｓｓｏｆｔｈｅｔｗｏｍｅｔｈｏｄｓꎬａｎｄｔｈｅａ￣ｎａｌｙｓｉｓｏｎｔｈｅｗａｔｅｒｆｉｌｌｉｎｇｃｏｎｄｉｔｉｏｎｍａｄｅｔｈｅｉｎｆｅｒｅｎｃｅｃｏｎｃｌｕｓｉｏｎｓｆｒｏｍｔｈｅｓｕｒｖｅｙｍｏｒｅｒｅａｓｏｎａｂｌｅａｎｄａｃｃｕｒａｔｅ.Ｋｅｙｗｏｒｄｓ:ｇｏｌｄｍｉｎｅꎬｗｏｒｋｅｄ－ｏｕｔａｒｅａꎬｈｉｇｈｄｅｎｓｉｔｙｒｅｓｉｓｔｉｖｉｔｙｍｅｔｈｏｄꎬＯＣＴＥＭ７１８。