类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用.pdf

第 59 卷 第 4 期2023 年 7 月地质与勘探GEOLOGY AND EXPLORATIONVol.59 No.4July，2023类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用顾观文1，2，王顺吉1，武 晔1，2，许志河1，2 （1.防灾科技学院地球科学学院，河北廊坊 065201；2.河北省地震动力学重点实验室，河北廊坊 065201）摘 要 时间域激发极化法是应用于金矿勘查的一种重要地球物理勘探手段，以往主要采用中梯扫面和激电测深组合的方法获取激电异常。中梯扫面受限于单一固定供电极距，仅能反映某一深度的激电异常分布；常规对称四极激电测深由于需要对固定测点布设不同极距的供电电极，工作效率较低，并且不适合三维探测。近些年随着大功率、阵列观测技术发展起来的类中梯激电装置可在观测区布设多组不同极距的供电电极，兼顾了激电测深与激电剖面二者的特点，为常年干旱的戈壁滩地区金属矿深部勘查提供了一种可能的探测手段。本文基于类中梯装置大功率三维激电测量技术在甘肃花牛山金矿区开展时间域激发极化法实验工作。实验结果表明：类中梯装置大功率三维激电观测能快捷高效获得研究区不同深度的激电异常；采用与类中梯装置相适应的三维反演方法对激电数据进行三维反演，可获取研究区地下三维精细电性结构；基于矿区的三维电性结构并结合相关地质资料构建的三维地质-地球物理模型，与矿区已知金矿点位置吻合，并发现了三处潜在成矿点。类中梯装置大功率三维激电测量技术在供电困难的戈壁滩地区金属矿勘查中应用效果明显。关键词三维地质-地球物理模型 类中梯装置 大功率 时间域激发极化法 花牛山金矿 甘肃中图分类号P631.3 文献标识码A 文章编号0495-5331（2023）04-0872-11Gu Guanwen,Wang Shunji,Wu Ye,Xu Zhihe.Application of the similar central gradient array high-power 3D induced polarization method in the Huaniushan gold exploration of Gansu ProvinceJ.Geology and Exploration,2023,59(4):0872-0882.0 引言时 间 域 激 发 极 化 法（Timedomain induced polarization method，简称 TDIP）已广泛应用于金属矿和非金属矿勘查、地下水寻找、油气开发、地热和城市地质调查等领域，并取得了明显的应用效果（余官娣，1987；Yatini and Laesanpura，2013；Yatini et al.，2016；孙仁斌等，2017；孙二峰和郭峰伟，2019；薛胜利和凌丹丹，2019；张刚和吕庆田，2019；戴国水，2020；庞军，2021；张文雨和胡浩，2021；罗敏玄等，2022；张参辉等，2022）。随着目前新一轮找矿工作的开展，对勘探目标的深度和精度提出了更高的要求（黄俊革等，2004；Ren and Tang，2010，2014；Pan and Tang，2014）。传统的中间梯度装置激电法虽然具有工作效率高的优点，但受限于采用单一固定长度的供电极距，主要用于激电快速扫面工作，仅反映地下某一深度的电性分布特征，无法刻画不同深度电性结构的变化，对地下介质的电阻率和极化率的分布显示不够全面和精细；传统的对称四极激电测深法由于需要对固定测点布设不同的供电电极距，工作效率较低，并且不适合三维探测。而近些年发展起来的类中梯激电装置可在观测区布设多组不同极距的供电电极，兼顾了激电测深与激电剖面二者的特点，与大功率供电、三维阵列观测顾观文收稿日期 2022-11-01；改回日期 2023-05-04；责任编辑 陈伟军。基金项目 中央高校创新团队项目（编号：ZY20215108）、国家重大科学仪器设备开发专项（编号：2011YQ050060）、河北省高等学校科学技术研究项目（编号：ZC2022056，ZC2021213）和河北省硕士在读研究生创新能力培养项目（编号：CXZZSS2023184）联合资助。第一作者 顾观文（1975年-），男，2020年毕业于吉林大学，获博士学位，教授，现主要从事电法勘探方法、理论及应用研究工作。Email：sun_。doi:10.12134/j.dzykt.2023.04.013872顾观文等：类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用第 4 期技术相结合，能快捷高效地获得研究区不同深度的激电异常，其空间分辨率高，探测深度大（行英弟和朱永余，1995；张慧等，2008；吴小平等，2015；陈伟军等，2016；Wang et al，2017；曹创华等，2018）。花牛山矿区是甘肃北山地区一个重要的金银铅锌多金属重晶石成矿带，近年来属于热点研究区域，许多学者利用不同的方法对其研究（聂凤军等，2002；杨建国等，2010a，2010b；云龙等，2021；孙雨等，2022；王二腾等，2022），大多以地质和化探手段为主，兼以少量的重力、磁法物探测量工作。该金矿区位于甘肃、新疆交界部位的戈壁滩区，属大陆性半沙漠气候，常年干旱，几乎没有地表水，近地表电阻率很大，导致在该地区开展电法测量工作较为困难。常规的激电法供电功率有限，输出电流难以穿透高阻层，无法获取矿区深部电阻率和极化率信息。为了解决常规激电法在戈壁滩供电困难、探测深度小和空间分辨率低等问题，本文采用大功率多组极距的类中梯装置在花牛山金矿区开展三维激电测量工作，通过在观测区布设不同极距的供电电极，并结合阵列观测技术，获得不同深度的电性特征信息。在此基础上对采集的时间域激电数据进行三维反演，获取研究区地下三维精细电性结构，结合矿区的地质和钻探资料，构建了三维地质-地球物理模型，解译了金矿的容矿有利层位和空间分布。1 类中梯装置三维激电测量方法中间梯度装置激电法（激发极化法）因其工作效率高、可以旁测等特点，在有色金属矿床勘探领域发挥了重要的作用（umi，1961；李金铭，2005；马一行等，2019；卢卯等，2020）。传统的中间梯度装置激电法采用单一固定长度的供电极距，主要用于激电快速扫面工作，只能反映地下某一深度的电性分布特征，无法刻画电性参数随深度的变化。近些年随着大功率阵列激电勘探设备和三维激电采集技术的发展，为了获得地下精细的三维激电信息，在传统中间梯度装置激电法的基础上发展形成了类中梯装置的三维激电测量法。类中梯装置激电法的基本原理和常用的时间域激电法几乎相同，仍然是以地下岩石、矿物之间的导电性差异为物理基础，研究人工供电形成地下稳定电场的时空分布及其变化规律特征以解决金属矿的勘查问题。类中梯装置激电法与传统中梯激电法相比，其优势在于可在观测区域两边布设不同极距的供电电极。根据目标体走向及工区工作条件，可以选择垂直于构造走向（横向）、平行于构造走向（纵向）或平行与垂直组合（横向、纵向组合）的方式布设供电电极。通过改变供电电极的极距，能够获取到地下不同深度的电性信息，兼顾了激电剖面测量和激电测深测量的特点，特别是适合于大功率供电、阵列式三维激电测量。图1为类中梯装置三维观测示意图，图中Ah1、Ah2、Ah3、Ahn、Bh1、Bh2、Bh3、Bhn为横向供电电极，Av1、Av2、Av3、Avn、Bv1、Bv2、Bv3、Bvn为纵向供电电极，M1、N1、M2、N2、Mn-1、Nn-1、Mn、Nn为测量电极，中间黑色圆点为观测区域的测量点，即测量电极MN的中点，棕绿色部分是异常体。类中梯装置三维激电测量获取的两个主要参数为视电阻率（s）和视极化率（s）。视电阻率s的计算式为：s=K(A，B，M，N)UMNI（1）式（1）中s为视电阻率（m）；A、B 为供电电极，M、N为测量电极；K(A，B，M，N)为装置系数，与供电电极 A、B 和测量电极 M、N 的位置有关，无量纲；UMN为测量电极MN间的电位差（V），I为发射电流（A）。时间域激发极化法测量取得的原始数据为全波形数据，在一个周期内具有两个完整的充放电过程，选取合适的采样时刻点来计算充电达饱和的总场电位差VMN和断电后不同时刻的二次场电位差V2 MN。在断去供电电流后观测和研究不同时刻的激电二次场的时间特性，得到多个时刻的s。计算视极化率s的公式如下：s=V2MNVMN 100%（2）式（2）中s是一个无量纲的值，VMN为总场电位差（V），V2 MN为二次场电位差（V）。2 工作区地质与地球物理特征2.1 地质特征地质特征2.1.1 地层地层花牛山矿区的矿体大多集中于震旦系洗肠井群（ZX）和中奥陶统花牛山群（O2hm）。（1）ZX主要是在花牛山矿田的北部，近似东西向展布，其主要分布在一矿区、二矿区和花牛山矿田的东部，区域内有二岩组（ZXb）、三岩组（ZXc）、四岩组（ZXd）露出。873地质与勘探2023 年ZXb：有酸性斑岩脉和石英脉的浅变质、浅海相的泥质岩，岩性比较单一。ZXc：主要由含有细粒状石英脉的大理岩夹绢云千枚岩和粉砂质板岩等岩性的浅海相碳酸盐岩和泥质岩构成。大理岩和千枚岩主要分布在一、二矿区及东部地区，是矿区铅锌矿的主要赋矿层位。四岩组（ZXd）：分布于一、二矿区北部，主要为板岩和角岩等。（2）O2hm花牛山群被印支期和华力西中期的中酸性岩体侵蚀了顶部和底部，大致以东西向分布在三矿区，岩性主要为黑云母角岩、斑点板岩、变粉砂岩、变玄武岩、变安山岩，局部夹大理岩。2.1.2 矿床地质特征矿床地质特征花牛山金矿点位于花牛山金银铅锌矿田东北部印支期花岗岩北接触带，产于一个次级同斜倾伏次级背斜构造内，背斜规模较小，轴向近东西，向南陡倾，向东倾伏。核部为第二岩性段（Jxpb3）粉砂泥质板岩，已透辉石角岩化，两翼为第三岩性段（Jxpc3）大理岩，含金磁黄铁矿体产于背斜两翼两种岩性接触带部位的钙硅质角岩中，大理岩对成矿有圈闭作用，如图2所示。2.1.3 岩浆岩岩浆岩区内发育有火山岩和侵入岩。火山岩主要赋存在洗肠井组中，主要岩性为玄武岩。可能的火山口位置在花牛山东（图2）。侵入岩体的岩性有钾长-二长花岗岩、花岗岩、斑状花岗岩、闪长岩、石英闪长岩、花风闪长岩、蛇纹石化超基性岩、煌斑岩脉等。其中主要岩体为花牛山岩体，其总出露范围大约为78 km2，东西长约40 km，南北宽约为几千米至8 km不等，由早、晚两期的侵入体所构成。早期侵入岩为花牛山岩片的主岩体，总出露范围约为54 km2，与其成矿作用关系也较为紧密。其主要岩性有钾长花岗岩与二长花岗岩，均具有中细粒花岗结构、块状构造。钾长花岗岩与二长花岗岩之间为渐变接触，并没有明确的岩相划分带，因此二者除了在钾长岩和斜长岩的成分方面有差异之外，其它岩相学、矿物学和地球化学性质也大致相同。晚期侵入岩呈椭圆状岩株侵于早期钾长花岗岩中，出露范围大约为24 km2，在部分区域还可见有0.31 m的冷凝边和早期二长花岗岩残留体。岩性主要是钾长花岗图1 类中梯装置三维激电法观测示意图Fig.1 Schematic diagram of the three-dimensional（3D）induced polarization method of observation in similar central gradient array 874顾观文等：类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用第 4 期石，肉红色，有粗粒的似斑状结构、块状构造。对于区内侵入岩体的时代，杨建国等（2010a）认为是华力西期至印支期构造-岩浆的结果，不存在燕山期岩浆活动。聂凤军等（2002）测定钾长-二长花岗岩40Ar/39Ar 同位素年龄为为 194.251.96 Ma，认为花牛山岩体属燕山期构造-岩浆活动的产物。2.1.4 构造构造区内构造发育情况具有如下特征：（1）褶皱花牛山矿区主体为一复式向斜构造，由4个次级向斜和3个次级背斜组成。褶皱轴向大体呈近东西、北东东走向。向斜核部由洗肠井组第三岩性段组成，花牛山主峰一带则由玄武岩夹炭质千枚岩、硅质岩及大理岩透镜体构成，两翼主要由第二岩性段等组成。（2）断裂矿区断裂发育，在花西滩-花西山-三矿区的断裂走向约为北西西向，在花黑滩-一矿带的断裂走向为北北西方向，在花牛山金矿断裂走向以北西向为主。2.2 地球物理特征地球物理特征花牛山金矿区电性参数测量结果如表1，矿区内岩、矿石电性参数存在明显的差异。图2 甘肃花牛山金矿地质简图（据杨建国等，2010a）Fig.2 Geological sketch of the Huaniushan gold mine in Gansu Province（modified from Yang et al.，2010a）1-中奥陶统花牛山群；2-蓟县系平头山组第三岩性段；3-蓟县系平头山组第二岩性段；4-蓟县系平头山组第一岩性段；5-印支期二长花岗岩；6-印支期钾长花岗岩；7-华力西期斑状花岗岩；8-华力西期闪长岩；9-华力西期石英闪长岩；10-华力西期花岗闪长岩；11-华力西期超基性岩；12-变砂岩或砂板岩；13-大理岩化灰岩；14-玄武岩；15-千枚岩或千枚状板岩；16-英安岩；17-断裂；18-古火山机构；19-氧化铅锌矿体或铁锰帽；20-金矿床；21-钼矿体；22-花牛山金银铅锌矿；23-花牛山金矿1-Middle Ordovician Huaniushan Group；2-third lithologic section of the Jixian System Pingtoushan Formation；3-second lithologic section of the Jixian System Pingtoushan Formation；4-first lithologic section of the Jixian System Pingtoushan Formation；5-Indosinian monzonitic granite；6-Indosinian moyite；7-Variscan porphyritic granite；8-Variscan diorite；9-Variscan quartz diorite；10-Variscan granodiorite；11-Variscan ultrabasic rock；12-meta sandstone or silt slate；13-marbleized limestone；14-basalt；15-phyllite or phyllite slate；16-dacite；17-fault；18-ancient volcanic edifice；19-lead-zinc oxide ore body or iron-manganese cap；20-gold deposits；21-molybdenum ore body；22-Huaniushan gold-silver-lead-zinc mine；23-Huaniushan gold mine875地质与勘探2023 年表1 花牛山金矿区岩（矿）石电性参数统计表Table 1 Statistics of electrical parameters for rocks（ores）in the Huaniushan gold mine序号123456岩石名称结晶灰岩钾长花岗岩含金磁黄铁矿石矿化角岩角岩化砂板岩条带状结晶灰岩样品数202020194014视极化率(%)最小值1.211.351.2410.410.460.74最大值5.212.8775.4750.9341.986.13算术平均值3.221.7350.8237.932.823.38视电阻率(m)最小值609.962934.2044.9496.26545.622054.25最大值25309.0716322.6713698.0564042.63179730.0029221.54算术平均值7368.165914.972297.956732.9028782.6014251.68注：数据由西安地质调查中心2015年提供。3 数据采集及成果初步分析本次数据采集工作使用的仪器设备为中国地质科学院地球物理地球化学勘查研究所自主研制、重庆地质仪器厂生产的大深度三维电磁探测系统。该系统可进行TDIP/FDIP、SIP、CSAMT、AMT等方法技术的测量，具有大功率供电、抗干扰接收、多参量提取等特点。大功率多功能电磁发射系统由EM-T160发射机、整流电源和200 kW柴油发电机组构成。系统最高发射电压可达1600 V，最高功率160 kW，可实现超大功率供电，保证在大深度探测时的信号强度以及足够的信噪比。接收系统采用多台套接收设备完成三维电磁测量，每个接收机有三个电通道，可以测量多个方向、多个测点的电信号。采用大功率多组极距的类中梯装置在该矿区开展三维激电测量，以获得金矿区深部精细电性结构。TDIP测量按照1：10，000物探规则测网要求布设，共布设13条测线，每条测线长度为960 m，线距100 m，点距20 m，测点数共637个。大功率多组极距的类中梯装置激电测量如图3所示，在工作区南北中心线（L720线）沿东西两端布置了 4组供电电极，供电电极与测线互相平行，供电电极距分别为2000 m、3000 m、4000 m和5000 m，接收电极 MN距离为40 m。供电电压近1000 V，供电电流4.215 A，测区观测点的电位差范围为35 mV（旁测线L660和L780）500 mV（主测线L720），满足测区信号强度的要求。不同供电极距视极化率平面等值线显示（图4），在测区南部大片花岗岩出露区，时域激电观测的视极化率值普遍较低，数值上一般3%，测网北部即金矿区视极化率数值不均匀，蓟县系千枚岩粉砂岩和结晶灰岩分布区的视极化率与南部花岗岩的视极化率值较接近，高值异常出现在千枚岩粉砂岩与灰岩接触带附近（测线L740L770），视极化率值大于4%，中心区域达到6%以上。高视极化异常形态上呈矩形，沿测线方向宽度约300 m，垂直测线方向即南北方向长度约 600 m。供电极距 AB=2000 m、3000 m和4000 m的视极化率高值异常形态及延伸范围（图4a、b和c）较吻合，当AB极距增加到5000 m时，视极化率异常中心区域有向南部（测线L720L740）迁移的趋势（图4d）。不同供电极距视电阻率平面等值线显示（图5），测区基本上为中、高阻电性特征，测区大部分区域视电阻率值普遍大于2000 m。局部低阻异常出现在测区中北部附近，即测线L720L760线的300500 m范围。低阻异常视电阻率值小于2000 m，形态上呈长方形，走向近南北，规模150 m600 m。不同供电极距反映的低阻区形态、位置基本一致。4 金矿区激电数据三维反演解释采用不完全高斯-牛顿法对获取的激电数据进行三维反演（吴小平等，2015），根据正则化理论（Tikhonov，1950），定义三维反演问题的目标函数（Tikhonov and Arsenin，1977）为：=12d+2m（3）式（3）中的d和m分别为数据拟合项和模型约束项，为正则化因子。调节的大小，可改变数据拟合项和模型约束项在目标函数中所占的权重，876顾观文等：类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用第 4 期值的选取范围为110，值越大，模型约束项在目标函数中所占权重越大，反之，数据拟合项所占权重大，一般可通过实验的方法来确定最佳的值（吴小平和汪彤彤，2001）。反演输入数据为观测点的电位差和视极化率，正则化因子=5，迭代反演5次，获得金矿区地下三维电阻率和极化率电性结构。基于获取的三维电性结构，并结合该地区的地质资料和钻探资料，构建金矿区的三维地质-地球物理模型（图6）。从图6可以看出，矿区存在4处孤立“低阻、高极化”异常。其中，在矿区西部边界存在2处“低阻、高极化”异常（异常范围未封闭，编号、），中北部“低阻、高极化”异常（编号）位于花牛山金矿采区、紧邻出露的花岗岩体的东侧，呈纺锤状，走向近南北，中心埋深约160 m；南部号异常为囊状孤立异常，号异常为扁豆状孤立异常，地表出露岩性为印支期花岗岩。结合矿区地质资料，号异常为已知似层状含金磁黄铁矿体的反映，据钻孔资料（钻孔位于号异常体上方，编号ZK1-10）和其他已知资料，该金矿体赋存于大理岩化灰岩与长英质角岩之间，走向近北西，倾向南东，在矿体与南北向花岗岩枝交汇部位，因后期叠加改造形成厚大矿体。金矿体中硫化物发育，是引起“低阻、高极化”异常的主要因素；号“低阻、高极化”异常位于南北向花岗岩枝西侧，与号图3 花牛山金矿区地质图与类中梯装置激电测量布置图Fig.3 Geological map of the Huaniushan gold mining area and the layout of the induced polarization survey of the similar central gradient array 1-第四系冲洪积砂砾层；2-蓟县系平头山组第三岩性段；3-蓟县系平头山组第二岩性段；4-印支期二长花岗岩；5-印支期钾长花岗岩；6-华力西期花岗闪长岩；7-石英脉/花岗岩脉；8-花岗斑岩脉/斜长花岗斑岩脉；9-供电电极；10-TDIP 测线1-Quaternary alluvial gravel layer；2-third lithologic section of the Jixian System Pingtoushan Formation；3-second lithologic section of the Jixian System Pingtoushan Formation；4-Indosinian monzonitic granite；5-Indosinian moyite；6-Variscan granodiorite；7-quartz vein/granite vein；8-granitic porphyry vein/plagiogranite porphyry vein；9-power supply electrode；10-TDIP survey lines877地质与勘探2023 年图4 花牛山金矿区TDIP不同供电极距视极化率异常平面等值线图Fig.4 Contours showing anomalies of apparent polarization at different power supply pole distances of TDIP in the Huaniushan gold mining areaa-AB=2000 m；b-AB=3000 m；c-AB=4000 m；d-AB=5000 m异常成矿条件相似。南部孤立“低阻、高极化”异常（、号异常）位于大片花岗岩体出露区、印支期不同侵入期次花岗岩分带线附近。在东西走向延伸方向，有蓟县系大理岩和长英质质角岩出露，在东南角“低阻、高极化”异常区地表花岗岩与大理岩接触带，见有少量含黄钾铁钒化含金阳起石矽卡岩氧化铁帽，故分析认为、号“低阻、高极化”异常区地表出露的花岗岩体中可能有蓟县系地层“捕虏体”，该“捕虏体”具有低电阻率电性特征，推测由隐伏多金属矿化物矿化引起。878顾观文等：类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用第 4 期图5 花牛山金矿区TDIP不同供电极距视电阻率异常平面等值线图Fig.5 Contours of apparent resistivity anomalies at different power supply pole distances of TDIP in the Huaniushan gold mining areaa-AB=2000 m；b-AB=3000 m；c-AB=4000 m；d-AB=5000 m5 结论本文通过在花牛山金矿区开展基于类中梯装置的大功率三维激电法实验工作，获得结论如下认识：（1）在获取花牛山金矿区地下精细三维电性结构的基础上，结合相关地质和钻探资料，构建的三维地质-地球物理模型清晰地展示了矿致异常赋存的有利部位，其中一处局部异常与矿区已知金矿点位置吻合，此外，发现了三处潜在成矿点，为矿区潜在矿产资源的评价提供了依据。（2）相比于传统的中梯装置只能获得某一深度的电性变化，在测区布设不同极距的供电电极，并结合三维阵列测量技术的类中梯观测装置，能够反映地下不同深度的电性特征。（3）基于大功率多组极距的类中梯装置三维激电测量技术，能够在长期干旱、戈壁滩等恶劣环境下快捷高效获得测区三维电阻率及极化率信息，且具有空间分辨率高、探测深度大等优点。879地质与勘探2023 年图6 花牛山金矿区三维地质-地球物理模型Fig.6 3D geological-geophysical models of the Huaniushan gold mining areaa-基于金矿区极化率信息构建的三维地质-地球物理模型；b-基于金矿区电阻率信息构建的三维地质-地球物理模型；1-大理岩化灰岩；2-印支期钾长花岗岩；3-印支期二长花岗岩；4-华力西期斑状花岗岩；5-角岩；6-花岗斑岩脉/斜长花岗斑岩脉；7-石英脉/花岗岩脉；8-异常体编号；9-金矿脉；10-钻孔a-3D geological-geophysical model based on polarization of gold mining area；b-3D geological-geophysical model based on resistivity of gold mining area；1-marbleized limestone；2-Indosinian moyite；3-Indosinian monzonitic granite；4-Variscan porphyritic granite；5-hornstone；6-granitic porphyry vein/plagiogranite porphyry vein；7-quartz vein/granite vein；8-abnormal body number；9-gold ore vein；10-drill hole and number致谢：感谢西安地质调查中心王凯教授级高工、杨建国研究员给予本文工作的指导和帮助。ReferencesChen Weijun，Hong Wanhua，Hao Qingqing，Chu Shaoxiong，Xing Baoshan.2016.Application of high-power induced polarization and continuous conductivity sounding for prospecting in the ChifengChaoyang gold concentration regionJ Geology and Exploration，52（1）：152-158（in Chinese with English abstract）.Cao Chuanghua，Deng Zhuan，Kang Fangping，Liu Jianxin，Tong Xiaozhong.2018.Three-dimensional forward and inspiration of abnormal cuboids under intermediate gradient way of IP method J.Computing Techniques for Geophysical and Geochemical Exploration，40（4）：459-467（in Chinese with English abstract）.Dai Guoshui.2020.Application of time-domain induced polarization method in groundwater resources exploration J.China Resources Comprehensive Utilization，38（8）：24-26（in Chinese with English abstract）.Huang Junge，Ruan Baiyao，Bao Guangshu.2004.Resistivity inversion on 3-D section based on FEMJ.Journal of Central South University of Technology（Natural Science），（2）：295-299（in Chinese with English abstract）.Li Jinming.2005.Geoelectric field and electrical explorationM.Beijing：Geological Publishing House：217-228（in Chinese）.Lu Mao，Huang Kai，Jiang Haimin，Zeng Guodao，Pu Qinglong，Tang Yongyong.2020.Application of the induced polarization method in mineral prospecting of the Xixiangou lead-zinc deposit in the northwestern Guizhou，China J.Acta Mineralogica Sinica，40（4）：502-509（in Chinese with English abstract）.Luo Minxuan，Li Xin，Wu Feng，Nie Xiaoli，Huang Jinyan，Yang Guolin.2022.An application effect of double frequency IP method in search and time domain IP method in Tibet J.China Manganese Industry，40（2）：45-49（in Chinese with English abstract）.Ma Yihang，Yan Tingjie，He Peng，Liu Zujian，Su Hang，Guo Shuo.2019.Exploration in coverage area：A case study of the Changtuxili Mn-Ag-Pb-Zn polymetallic ore depos it，Inner Mongolia J.Geophysical and Geochemical Exploration，43（4）：709-717（in Chinese with English abstract）.Nie Fengjun，Jiang Sihong，Liu Yan，Chen Wen，Liu Xinyu.2022.（40）Ar/（39）Ar isotopic age dating on K-feldspar separates from eastern Huaniushan granite，Gansu Province，and its geological significanceJ.Chinese Journal of Geology，（4）：415-422（in Chinese with English abstract）.Pan Kejia，Tang Jingtian.2014.2.5-D and 3-D DC resistivity modelling using an extrapolation cascadic multigrid methodJ.Geophysical Journal International，197（3）：1459-1470.Pang Jun.2021.Study on the application of induced polarization method in cryptoexplosive breccia tube type lead-zinc depositJ.World Nonferrous Metals，（3）：211-212（in Chinese with English abstract）.Ren Zhengyong，Tang Jingtian.2014.A goal-oriented adaptive finiteelement approach for multi-electrode resistivity systemJ.Geophysical Journal International，199（1）：136-145.Ren Zhengyong，Tang Jingtian.2010.3D direct current resistivity modeling with unstructured mesh by adaptive finite-element methodJ.Geophysics，75（1）：H7-H17.umi F.1961.The induced polarization method in ore investigation J.Geophysical Prospecting，9（3）：459-477.880顾观文等：类中梯装置的大功率三维激电法在甘肃花牛山金矿勘查中的应用第 4 期Sun Erfeng，Guo Fengwei.2019.Application of induced polarization method in urban geological surveyJ.China Energy and Environmental Protection，41（7）：124-128（in Chinese with English abstract）.Sun Renbin，Chu Lixia，Zhao Yijun，Wang Ning，Liu Qiang2017.Application of the time-domain induced polarization method to a Zn-W polymetallic deposit in Xianghuang Banner，Inner Mongolia J Geology and Exploration，53（3）：519-527（in Chinese with English abstract）Sun Yu，Liu Jiajun，Zhao Yingjun，Zhai Degao，Liu Zhenjiang，Zhang Yafeng，Zhang Fangfang，Tian Feng，Qin Kai.2022.Alteration mineral mapping based on the GF-5 hyperspectral data and its geological applicationAn example of the Huaniushan area in Guazhou County of Gansu Province J.Geology in China，49（2）：558-574（in Chinese with English abstract）.Tikhonov A N.1950.On determining electric characteristics of the deep layers of the Earths crust J.Doklady Akademii nauk，73（2）：295-297.Tikhonov A N，Arsenin V Y.1977.Solution of Ill-Posed Problems M .New York：John Wiley and S