赵理芳，李希元，李成立，赵俊康，包怡，陈中，刘和花，梁杰，王佩业
（1江苏省有色金属华东地质勘查局八一四队，江苏 镇江 212004；
2长江大学地球科学学院，湖北 武汉 430100；
3中国地质调查局沈阳地质调查中心，辽宁 沈阳 110000；
4紫金矿业集团股份有限公司矿产地质勘查院，福建 厦门 361000；
5江苏省地质局，江苏 南京 210018）
斑岩型矿床是指与中酸性浅成侵入岩（斑岩）有关的，以浸染状、细脉浸染状为主要矿化样式的岩浆热液矿床，矿种以Cu为主、伴生Au，具品位低、规模大的特征（叶天竺等，2014）。斑岩铜矿主要形成于聚合板块的活动大陆一侧，一般为典型的边缘构造岩浆活动带的陆缘弧和岛弧环境（张寿庭等，2011；
熊欣等，2014）。多宝山矿集区位于大兴安岭与小兴安岭交汇部位，中亚-蒙古斑岩铜矿带的东段，矿区位于新开岭隆起和古生代洋盆正负构造单元的边缘（杜琦等，1988），是中国重要的斑岩型铜钼矿成矿带，矿产较为丰富，以铜、金多金属矿产为主，包括多宝山铜钼矿床、争光岩金矿床、铜山铜矿床、三矿沟铁铜矿床等多金属矿床（石国明等，2018）。
多宝山矿集区的地质工作始于20世纪50年代，勘探研究于90年代进入高峰期，完成了区内1∶20万区调、重力、化探及1∶5万航测工作，并对重点成矿区（带）进行了1∶5万等大比例尺区域地质调查与矿产地质调查工作，这些工作比较系统地划分了区内的地层层序、岩浆活动的期次和隐伏岩体的范围，研究分析了多宝山矿集区的构造特征，对区内的矿产进行了成因及成矿规律研究，总结和划分了部分成矿远景区（杜琦等，1988；
2008；
石国明等，2018）。
矿集区植被覆盖，虽然岩石矿化时电阻率值会明显下降，形成局部低阻异常（范正国等，2004；
郑军，2018），但矿集区的采矿区、铜矿生产区、变电站、高压线、煤矿等干扰给电磁法的数据采集带来了较大的影响，此外，矿集区大量分布黄铁矿所造成的低阻异常，对于斑岩型矿床的深部勘查也带来了困扰。加上多宝山斑岩铜矿床具有多期成矿及叠加改造的特点（赵忠海等，2012；
刘宝山，2020；
刘宝山等，2020），造成了深部勘探整体研究程度低。
2017年—2020年中国地质科学院地质研究所研究牵头并联合国内等众多所高校、科研单位开展的《北方东部复合造山成矿系统深部结构与成矿过程》国家重点研发计划项目，在多宝山矿集区开展了大量的地球物理实物工作，获得了较高质量的野外数据。
由于多宝山斑岩型矿床的成矿地质体为花岗闪长斑岩，主要矿带位于斑岩体上盘（赵元艺等，1997），根据硫同位素组成的梯度变化及其他资料，可推断矿液进入矿带前首先流经了斑岩体（马德有，1984），因此，查明成矿斑岩体是斑岩铜矿找矿预测最关键的步骤（卿敏等，2019）。本文以三维广域电磁法和高精度重、磁数据为基础，结合地质、钻孔资料，对多宝山矿集区的隐伏斑岩体进行了研究。从已知到未知推测出斑岩体的位置，结合斑岩型矿床的矿化蚀变特征和岩石物性特征，根据电性特征推测了矿体位置，钻孔验证与解释相符。研究成果对于干扰较大的矿区勘探方法的选择和斑岩型矿床深部找矿具有很好的借鉴作用，同时对于多宝山、铜山、争光等矿区外围和深部找矿具有指导意义。
1.1 地质概况
多宝山矿集区的地层隶属兴安岭地层分区，区内出露地层由古至新为古生界、中生界及新生界。古生代地层主要发育在多宝山铜矿一带及裸河两岸，关鸟河两岸、东南部也有小面积出露，主要为奥陶系多宝山组、铜山组，其次为志留系及少量的泥盆系、二叠系；
中生界在测区南部广泛分布，主要位于三岔河村、清水河村以南，测区内中生界出露有三叠纪火山岩地层、白垩纪火山岩和沉积岩地层，以白垩系最为发育（石国明等，2018）；
新生界第四系主要分布于河谷两侧（图1）。奥陶系—泥盆系为海相地层，其中奥陶系火山岩发育，自晚二叠世开始为陆相沉积。中奥陶统属海相火山岩、火山碎屑岩和碳酸盐岩建造，构成了岛弧形火山-沉积岩系。多宝山组是由中-酸性角砾岩、熔岩及凝灰岩组成的海相火山岩，夹少量沉积岩层沉积岩，与成矿关系密切，是矿集区主要的赋矿层位（杜琦等，1988；
石国明等，2018）。
图1 多宝山矿集区地质图(据石国明等，2018）1—第四系；
2—白垩系火山沉积岩；
3—清水河组安山岩；
4—石炭系—二叠系陆相沉积；
5—志留系—泥盆系碎屑岩；
6—奥陶系上统沉积地层；
7—铜山组浅海相沉积；
8—多宝山组三段中酸性-酸性熔岩、火山碎屑岩；
9—多宝山组二段中性火山岩、火山碎屑岩；
10—多宝山组一段中性-中酸性火山碎屑岩；
11—燕山期二长花岗岩；
12—燕山期花岗闪长岩；
13—印支期二长花岗岩；
14—印支期花岗闪长岩；
15—华力西期花岗闪长岩；
16—加里东期花岗闪长岩；
17—加里东期花岗闪长斑岩；
18—闪长岩；
19—实测断层；
20—三维研究区；
21—广域电磁点；
22—分析剖面；
23—验证钻孔；
24—已知矿体Fig.1 Geological map of the Duobaoshan ore concentration area(after Shi et al.,2018)1—Quaternary;2—Cretaceous volcanic and sedimentary rocks;3—Andesite of Qingshuihe Formation;4—Carboniferous and Permian continental sediments;5—Silurian and Devonian clastic rocks;6—Upper Ordovician sedimentary rocks;7—Tongshan Formation shallow marine sediments;8—Intermediate acid-acid lava and pyroclastic rocks of the third member of Duobaoshan Formation;9—Intermediate volcanic rocks and pyroclastic rocks of the second member of Duobaoshan Formation;10—Intermediate to intermediate-acid pyroclastic rocks of the first member of Duobaoshan Formation;11—Yanshanian monzonitic granite;12—Yanshanian granodiorite;13—Indosinian monzonitic granite;14—Indosinian granodiorite;15—Variscan period granodiorite;16—Caledonian granodiorites;17—Caledonian granodiorite porphyry;18—Diorite;19—Measured fault;20—Three-dimensional study area;21—Wide field electromagnetic station;22—Study profile;23—Verification drill hole;24—Known ore bodies
多宝山矿集区岩浆活动有火山喷发和岩浆侵入2种形式，岩浆岩分布较广，岩石类型复杂，从超基性岩到酸性岩均有分布。矿集区曾发生加里东期、华力西期、印支期、燕山期4期岩浆侵入活动，以加里东期和华力西期为主。侵入岩多发育在NW向和NE向断裂交汇处，以酸性、中酸性花岗岩类为主，其次为中性闪长岩类及少量的基性辉长岩侵入体。与铜矿成因关系密切的侵入岩是加里东中期花岗闪长岩与花岗闪长斑岩，它们组成复式岩体。多宝山斑岩型铜矿床的形成，与加里东期浅成相花岗闪长斑岩密切相关，并因海西期的岩浆活动进一步富集成超大型铜矿床，严挌受NW向构造的控制，而燕山运动时遭受改造。燕山期变形可分为2期：早期为陆内造山带发育及坍塌，构造线方向北东东；
晚期为北东向断裂左行走滑。岩体中心部位比较厚，分枝多，形态复杂，北西延长方向上分枝少，形态规整，岩墙变薄，反映了岩浆是沿花岗闪长岩体的中心部位被动侵位的（赵元艺等，1997；
杜琦等，1988）。
1.2 矿化蚀变带特征
斑岩型铜矿一般均存在明显的蚀变带，其蚀变组合及其分带模式俗称“大白菜模式”，由内到外依次为：石英内核→钾化带（黑云母-钾长石带）→似千枚岩化带（绢云母-石英带)→泥化带→青磐岩化带，斑岩型铜（钼）矿化往往分布在钾化硅化带外的绢英岩带（绢云母-石英带）中（叶天竺等，2017）。
多宝山矿床具有典型的斑岩型铜矿床的蚀变分带特征，蚀变分带以花岗闪长斑岩体为核心，环绕着岩体呈北西向拉长的环带状分布，各带之间的界线不清，具有相互穿插、叠加现象，表现出多期次的蚀变特征，由内向外可划分为4个蚀变带：石英核、钾化带、绢云母化带、青磐岩化带（图2），多宝山铜（钼）矿的4个矿带均产于绢云母化带中（尹冰川等，1997；
杜琦等，2008）。
图2 多宝山铜矿床蚀变带分布示意图（据杜琦等，2008）Fig.2 The sketch map of alteration zonation of the Duobaoshan copper deposit(after Du et al.,2008)
从多宝山矿区地质图（图3）可见大砬子斑岩体（花岗闪长斑岩）呈北西向展布，与构造走向一致，矿体环斑岩体分布。斑岩体的地表出露面积较小，仅0.17 km2，延伸到100 m标高时，面积为0.18 km2，在深部，斑岩体则有明显膨大的趋势（赵元艺等，1997；
杜琦等，1988），隐伏斑岩体的范围可能远比地表更大，如D66勘探线剖面图显示在地表斑岩体未出露，但在深部仍有斑岩岩枝分布（图3），推测有较大的斑岩体隐藏在深部。
图3 多宝山矿区地质图（a）及D66勘探线剖面图（b，据杜琦等，2008）Fig.3 Geological map of the Duobaoshan mining area(a)and the cross-section of the No.66 profile line(b,after Du et al.,2008)
地质体的物性差异是引起重、磁、电异常的根本原因，掌握地质体的物性参数的特征及变化规律是做好重、磁、电资料地质解释工作的基础。
依据1∶20万区域地质图及1∶5万研究区地质图，设计了5条采集路线，在多宝山矿集区内100余处露头区采集岩石物性标本共计974块，采集物性样品涉及地层近20个、侵入岩体11个。此外还收集了研究区内具有代表性的6口钻孔的岩芯样品260块。对所有样品进行了甄别和加工，分别在现场与室内进行了密度、磁化率、电阻率物性参数测定，分组进行了正态统计分析，统计结果示于表1至表4。
从地层物性统计来看（表1），地层从新到古表现为密度加大的变化趋势，反映了压实作用对密度变化的影响。由于各时代地层沉积岩性的不同，局部地层的密度有所变化，同时各地层的岩性、密度特征也反映了地层沉积环境的变化。古生界多宝山组、铜山组含有大量的火山岩成分，均具有较大的密度。从地层磁性的特点可以看出，正常沉积岩为主的地层普遍为弱磁性或是无磁性，含有火山岩成分的地层具有磁性，磁性最强的地层为清水河组及甘河组，其次为龙江组，这3个强磁性层均与火山岩有关，多宝山组也因含有火山岩成分表现为中等的磁性特征。在电性方面，分布在多宝山矿集区的各地层的电阻率普遍较高，电阻率的平均值均在1000 Ω·m以上，多宝山组比相邻地层具有较高的电阻率。
表1 多宝山矿集区地层物性统计一览表Table 1 Physical property statistics of different stratum in the Duobaoshan ore concentration area
表2为岩浆岩的物性统计，可以看出从基性到中性到酸性，密度逐渐降低，磁性变化较大，整体呈降低的趋势，电阻率具有逐渐升高的特征，整体为高阻。玄武岩和安山岩为高密度、中强磁性、中高电阻率；
花岗闪长岩为中高密度、中磁性、高电阻率；
花岗岩为中低密度、弱磁性、高电阻率。
表2 岩浆岩物性统计一览表Table 2 Physical property statistics of magmatic rock samples
花岗闪长岩为本区的主要赋矿围岩，按照侵入的不同时期，对其物性进行了统计，见表3。从表3中可以看到加里东期和华力西期花岗闪长岩密度较为接近，密度都较高，达到2.75 kg/m3，印支期和燕山期花岗闪长岩密度较低，分别为2.65 kg/m3和2.69 kg/m3；
磁化率表现为时代越新，磁化率越大，可能与区内多期岩浆活动造成热液退磁有关；
电阻率则大都表现为高阻特征。
表3 不同时期花岗闪长岩标本物性统计Table 3 Physical property statistics of granodiorite samples of different ages
围岩经过矿化或蚀变后，岩石的物性均发生不同程度的变化（表4）。相比未蚀变的花岗闪长岩，蚀变、矿化的花岗闪长岩整体密度呈现增大的变化特征，推测为岩石在经过蚀变、矿化后，密度较大的金属矿物含量大量增加；
在磁性方面，矿化蚀变的岩石磁性弱于未蚀变的花岗闪长岩，这一特征说明岩石在矿化及蚀变的过程中，温度较高的流体热液对岩石、矿石原有磁性起到了退磁作用；
在电性方面，经过矿化、蚀变后的岩石、矿石，其导电性发生了显著的变化，蚀变后的岩石不同程度地表现为电阻率变低的特征。
表4 矿化蚀变物性变化统计表Table 4 Statistics table of physical property changes of mineralization or alteration
3.1 数据来源
1990年—1992年，黑龙江省物探大队在多宝山地区开展了1∶20万重力和1∶20万航磁测量工作；
1992年—1993年，地矿部遥感中心航测大队在多宝山地区完成了1∶5万航磁；
2013年中国地质科学院矿产资源研究所在多宝山矿集区进行了1∶1万大比例尺低空航磁勘查，其精度为±4.0 nT（杨生，2013）。
2017年—2020年，中国地质调查局沈阳地质调查中心为了探索综合地球物理在解决矿集区3000 m以浅深部透明化地质结构的作用及应用效果，在多宝山典型铜（金）矿集区选择了一面积为84 km2的长方形重点区域作为试验区开展了三维广域电磁法及重力勘探，勘探网度为150 m×150 m，外围还部署了点距为150 m的广域电磁法、重力剖面，共采集广域电磁法坐标点4614个、重力点5048个。广域电磁法Ⅰ级品率为95.6%，Ⅱ级品率4.4%（赵理芳等，2020），重力总精度为±0.078×10-5m/s2。
重力数据以2017年的高精度重力数据为主，外围采用1990年—1992年的重力数据作为扩边数据。磁力数据以2013年高精度航磁数据（比例尺1∶1万）为主，外围以较低精度1990年—1993年航磁数据（比例尺1∶5万、1∶20万）作为扩边数据。重力高精度数据范围与三维广域电磁数据范围一致，而航磁高精度数据范围大于研究范围，为尽可能地减少边缘影响，还是对不同精度的重、磁数据进行了扩边拼合处理。为保证数据精度和可靠性，拼合数据时，2017年的高精度重力数据保持不变，通过Oasis montaj软件与外围数据进行Grid knit拼合处理。航磁数据通过Oasis montaj软件进行了Microleveling校正，并同样以高精度数据为准，对不同的数据进行了Grid knit拼合处理。
广域电磁法（wide field electromagnetic method，简称WFEM）是中南大学何继善院士提出的（何继善，2010）。可在不限于“远区”的“广大区域”测量，颠覆了频率域电磁法只能在“远区”测量的思想，满足了“深地探测”战略3000 m以浅矿产资源勘探技术的要求（何继善，2019）。本次广域电磁法数据采集频率段为8192～0.015 625 Hz，为了减少场源的影响，全区只在测区的西北方向布设了唯一场源，场源AB平行于测线方向（SW-NE）布设，长度为1 km，收发距为30～11 km。
针对矿集区干扰源较多的问题，在矿集区进行了大地电磁测深（MT）和广域电磁（WFEM）采集试验，从采集的视电阻率曲线（图4a～d）可以看到,在厂房和采矿区附近干扰严重，而同点位的广域电磁（WFEM）点的数据质量明显优于大地电磁测深（MT）点。鉴于在多宝山矿集区广域电磁法数据能取得较高质量的数据，选择了在研究区施行三维广域电磁法。
图4 多宝山矿集区MT点与WFEM点视电阻率曲线对比图a.高压线干扰点；
b.厂房干扰点；
c.采矿干扰点；
d.一般干扰点Fig.4 Comparison of apparent resistivity curves of MT stations and WFEM stations in the Duobaoshan ore concentration areaa.Stations disturbed by high voltage wires;b.Stations disturbed by the plant;c.Stations disturbed by mining;d.General stations
3.2 重磁场特征分析
从布格重力异常平面图（图5a）上可以看到，矿集区重力高值异常呈一个“Y”型，“Y”的张口朝北东向，显示为出北西-南东的弧形特征，总体表现为西南高、中部高、北东低、东南低，重力高与重力低之间以密集梯级带相过渡。对照区域地质图结合物性资料分析，重力高异常主要由高密度的古生界或基底隆起引起；
东北部的重力低主要由低密度的中酸性岩所引起；
南部的重力低主要为中生界盆地造成。研究区的外缘南、北均存在明显的北东向重力异常梯度带，反映了存在大型断裂。尹冰川等（1997）认为，多宝山矿集区位于黑河-嫩江海西期拼接带及嫩江-呼玛走滑断裂所挟持的向南西突起的弧形构造带内，重力异常带展示了由这2条断裂所挟持的多宝山岛弧区。
图5 多宝山矿集区重力（a）与航磁（b）异常平面图Fig.5 Plan of gravity(a)and aeromagnetic(b)anomaly in the Duobaoshan ore concentration area
从已知矿床与布格重力异常叠合可以看到，多宝山铜矿、铜山铜矿和争光岩金矿均处于重力异常梯度带的交汇处。重力异常梯级带的交汇反映了断裂发育交汇，为构造活动地带，断裂为含矿热液流动提供通道，断裂交汇处压力较小，为成矿物质的沉淀提供了场所，为有利的成矿构造面。多宝山铜矿和铜山铜矿区的重力低反映了低密度的花岗闪长岩分布，重力高主要由多宝山组及下伏高密度基底地层引起，反映了花岗闪长岩和多宝山组与成矿作用关系密切；
争光岩金矿区，结合地表地质分析，与高密度闪长岩和多宝山组有关。
航磁化极异常（图5b）整体特征为：南侧以宽缓负磁异常为主，零星分布多个规模较小的正磁异常；
北侧分布多个规模较大的正磁异常，均与不同时期花岗闪长岩体的分布有关，根据物性测定资料和地表地质可知，年代越新的花岗闪长岩，其磁性越强。区内负磁异常主要反映了无磁性～弱磁性的地层及酸性侵入岩。从已知矿床与磁力异常叠合上看，多宝山铜矿、铜山铜矿位于低缓的负磁异常区，表明与斑岩铜矿有关的花岗闪长岩岩体具无磁或弱磁性；
争光岩金矿处于规模较小的条带状正磁异常与平缓负磁异常的接触部位，靠近负磁异常一侧，显示与金成矿关系密切的闪长岩体的外接触带的成矿条件更为有利。
3.3 三维反演
为获得反映深部地质结构的物性数据体，对重、磁、电进行了三维物性反演。利用加拿大哥伦比亚大学李耀国教授研发的MAG3D和GRAV3D软件完成了重、磁三维反演。重、磁采用上延3 km后的剩余场作为目标场，共119×144×81剖分单元，基本网格为200 m×200 m×100 m，外扩剖分单元遵从原则：平面上从里至外；
纵向由浅至深逐渐变稀疏。X、Y、Z各方向具体剖分如下：
X：800 400 200 113*200 200 400 800
Y：800 400 200 138*200 200 400 800
Z：52*100 26*200 3*4583.33
广域电磁数据采用中南大学戴世坤教授研发的“重磁电三维反演成像解释一体化系统”完成了电法三维反演，剖分网格为256×256×151，水平网格均匀剖分，垂向网格自地表至深部网格间距逐渐增大。
通过三维反演获取了反映地下三维物性结构特征的密度体、磁化率体、电阻率体（图6a～c），从图中可以看到，浅表多为高密度、弱磁、中阻，对应地表地质图，多为岩体与奥陶系，风化导致电阻率有所降低；
在中深部存在大片的高阻，其密度和磁性变化较大，依据标本物性特征，推测高密度及弱或无磁性体为奥陶系，低密度体为隐伏岩体，磁性越强其年代越新；
三维测区南部的低密度、中磁、中低阻体，依据地表地质和岩石标本物性特征，推测其为中生界火山岩沉积体。
图6 多宝山重、磁、电三维反演立体视图a.视密度；
b.视磁化率；
c.电阻率Fig.6 Stereoscopic view of 3D inversion density,magnetic susceptibility and resistivity of Duobaoshan gravity,magnetic and electromagnetic dataa.Apparent density;b.Apparent magnetic susceptibility;c.Resistivity
4.1 已知斑岩体异常特征
斑岩体在多宝山矿集区出露较小，其中面积最大的位于多宝山矿区，即大砬子花岗闪长斑岩岩体，其出露面积也仅0.16 km2。斑岩体与花岗闪长岩体的侵入接触界面复杂，形成犬牙交错的接触带，接触带处的花岗闪长岩明显碎裂，并有3～10 cm宽的热变质圈，使花岗闪长斑岩与花岗闪长岩体之间的界线变得不十分清晰。
为研究已知斑岩体的重、磁、电异常特征，选取了穿过大砬子花岗闪长斑岩岩体的AA"剖面线（位置见图1、图3），其外围为花岗闪长岩，花岗闪长斑岩位于花岗闪长岩的中心部位。在重、磁、电三维反演断面（图7）上，花岗闪长岩表现为低密度、弱磁、高阻的特征，花岗闪长斑岩位于花岗闪长岩的中间部位，为低密度、弱磁、中高阻的特征，花岗闪长斑岩周围具有密度、磁化率升高以及电阻率降低的特征，与岩石物性特征一致，推测为矿化、蚀变的影响，在矿化部位，电阻率显著降低。
图7 多宝山矿区AA"线重、磁、电反演异常特征图a.视密度断面；
b.视磁化率断面；
c.电阻率断面Fig.7 Characteristics of gravity,magnetic and electric inversion anomaly of AA"Line in Duobaoshan mining areaa.Apparent density section;b.Apparent susceptibility section;c.Resistivity section
4.2 隐伏斑岩体推测
多宝山矿区有相对较大面积的花岗闪长岩体和闪长斑岩体出露，而在铜山矿区地表未见花岗闪长岩体，作为斑岩型铜矿最重要直接的标志，斑岩体在铜山矿区一直未被地下普查和详查工程所揭露（王喜臣等，2007；
王乐等，2017）。仅在ZK1064-2、ZK1064-4钻孔发现4条花岗闪长斑岩岩枝，岩枝厚度小于1 m（石国明等，2018）。
根据矿床地质、地球化学及地球化学异常特征可得出，多宝山矿床、铜山矿床的成因类型均为斑岩型，区别在于铜山矿床的斑岩体隐藏于深部（赵元艺，1995）。由于铜山地区成矿地质体并未揭露，那么深部肯定还存在隐伏成矿地质体，存在着较大的找矿空间（叶天竺等，2014）。杜琦（2008）认为，铜山断层将铜山矿床矿体的顶部推覆上来。矿田已揭露的几处大型矿体产状陡倾，延深稳定，由顶部向下有明显的变宽现象，推测断层下盘的矿体资源量较大。但铜山断层也使成矿后的矿体发生了截切和错动，导致矿化中心不明（金山岩等，2014；
庞绪勇等，2017），给深部勘查带来了困难。王乐等（2017）根据铜山断层上、下盘蚀变带的类型、组合、强度对比分析，推测热液蚀变成矿中心与斑岩体顶部大致位于ZK1064-2钻孔正下方，垂直深度可能大于1400 m，斑岩体向西倾（图8）。
图8 铜山(T1064勘探线)推测斑岩体位置及热液蚀变矿化中心（据王乐等，2017）Fig.8 Inferred porphyry body location and hydrothermal alteration mineralization center in Tongshan(T1064 line)(after Wang et al.,2017)
从铜山铜矿BB"剖面线（T1064勘探线两侧延长）重、磁、电三维反演剖面的异常特征图（图9a～c）看到，在斑岩岩枝往深部延伸的方向，处于低密度、弱磁化率、中高阻区。类似于多宝山155线（图7）花岗闪长斑岩边缘的矿化、蚀变造成周缘密度、磁化率升高、电阻率降低的特征，初步推测存在斑岩体及矿体。由于电法相对于重、磁在纵向上分辨率更高，故推测的斑岩体位置以电法为主。从多宝山矿区AA"重、磁、电异常特征（图7a～c）可知，斑岩体为中高阻特征，由于斑岩型矿床矿化蚀变带呈环状，根据岩石标本物性统计可知矿化蚀变使岩石电阻率显著降低，故形成了环绕中高阻斑岩体的中低阻矿化蚀变带，在BB"线电法反演断面中可以看到，推测的中高阻斑岩体两侧为低阻特征（图9a～c）。推测出的BB"线花岗闪长斑岩的位置，位于海拔-900 m（深度1400 m）左右，向西南倾，与王乐等（2017）预测的斑岩体位置和倾向（图9a～c）较为吻合。从物探上证实了铜山矿床深部存在较大的斑岩体，也佐证了铜山铜矿其成因类型为斑岩型矿床。
图9 铜山（BB"线）剖面重、磁、电异常特征图a.视密度断面；
b.视磁化率断面；
c.电阻率断面Fig.9 Feature map of gravity,magnetic and electric anomalies in section of Tongshan(BB"line)a.Apparent density section;b.Apparent susceptibility section;c.Resistivity section
在多宝山、铜山、争光3个矿区的中深部都有与铜山斑岩体类似的电性异常特征（图10），即矿体对应着低阻、矿体之间为中高阻凸起，故可推测多宝山矿区、铜山矿区、争光矿区，均存在着隐伏的斑岩体，进一步证实了多宝山、铜山、争光属同一成矿系统——斑岩型铜矿系统。根据电性异常特征，中高阻（推测隐伏斑岩体）旁的中低阻（推测矿体）范围较广，故可判断多宝山铜矿外围可以扩大勘查，铜山铜矿矿体有进一步向下延伸的趋势，争光岩金矿在浅部为浅成低温热液型金矿床，多宝山组厚度较大，其下部有发现斑岩型铜矿的可能。由于早期的斑岩事件可能引起广泛的硫化作用，形成的溶液具有较高的溶解金和搬运大量Au(HS)2的能力，使其运移到浅部形成浅成低温热液型金矿床（陈革，1998）。铜山铜矿的上部围岩中发现了一定规模的金矿床，说明多宝山矿集区古生代构造岩浆活动带不但可以形成斑岩型铜矿床，还可以在铜矿上部形成热液型金矿床，形成“上金下铜（钼）”的成矿模式（谭成印等，2010）。矿集区北西端相对南东端剥蚀程度要高，成矿深度较浅的浅成低温热液型矿床很有可能已经被剥蚀（郝宇杰，2015）。
在矿区典型的电性断面（图10）中还可以发现，虽然矿体在中深部都为中低阻特征，但在浅部（海拔500 m以浅）矿体所处部位电阻率虽然仍低，但与浅部非矿部位比较，其电阻率相对较高，结合当地的地下水位高程为200～300 m的情况（李文睿，2016），推测地表由于含水及风化影响，电阻率普遍降低，而矿化区域风化、含水对其影响相对更小，导致矿化区电阻率反而比围岩高。
图10 典型矿区矿体、电阻率反演断面图及推测斑岩体Fig.10 Ore bodies,resistivity inversion section and inferred porphyry body of typical mining areas
从图10中DD"线可以看到铜山矿区已揭露的矿体位于推测斑岩体的右侧（北侧），根据斑岩型矿床的矿化蚀变对称分布的特征来看，推测斑岩体左侧（南侧）显示的低阻异常信息很有可能为隐伏矿体。
2020年，紫金矿业在铜山东南区域实施验证钻孔，该孔在孔深1150 m处揭露到一厚度为2 m左右的压扭性破碎带，其中可见厚20 cm的断层泥，和地表及坑道中所见到的铜山断层特征相似，穿过断层后即开始见到铜矿化体，和物探解译推测的位置也较为吻合。据矿方反馈，验证孔所揭露的矿化体铜品位多在0.40%～1.0%之间，局部见强硅钾化、黄铜矿化蚀变，平均品位1.62%，穿矿厚度30 m，单样品位最高4.33%。钻孔连续穿矿厚度600.16 m，平均品位0.684%，且因井故未穿透矿体。目前，下盘矿矿体走向控制长度超过700 m，最大垂深800 m，水平厚度30～300 m，控制以及推断的工业矿总资源量超过100万吨。
从电法三维反演断面（FF"线）及钻孔轨迹图（图11）中可以看到，矿体所处部位来看位于铜山断层的下盘（铜山断层浅部由地表地质、矿区钻井控制，深部根据电法反演断面推断），矿体处于中高阻背景下的中低阻，该测线与DD"测线临近（位置见图1），电性特征（图11）也与DD"测线（图10）相似，其中低阻体为图11的DD"线中推测斑岩体左侧（南侧）中低阻的延伸，应为同一矿体。
图11 FF"线三维反演电性断面及钻孔轨迹Fig.11 3D inversion electrical section and borehole trajectory of FF"line
由于矿化条件的差别，从多宝山矿区地质图（图3a）中可以看到，矿体不是围绕斑岩体绝对对称的，加上铜山断层的改造，在FF"线剖面上围绕中高阻的斑岩体对称的低阻异常并不明显，有待进一步研究分析。
（1）通过标本物性的统计分析，可以发现多宝山矿集区主要的成矿围岩-花岗闪长岩具有高阻、低密度、弱磁的特征，同时具有时代越新，磁性越强的特点；
经过矿化或蚀变后的岩体物性总体表现为密度增大、磁性减弱、电阻率显著变低的特征。这为应用重、磁、电资料解析奠定了重要的物性基础。
（2）由重磁场异常特征对照已知矿体的位置，可以发现多宝山矿集区花岗闪长岩和多宝山组与成矿作用密切相关，矿体主要位于重力异常和磁力异常的变化带，对应地质的构造交汇处及岩体的外接触带上。
（3）查明成矿斑岩体是斑岩铜矿找矿预测最关键步骤，通过从已知到未知，由浅及深，结合斑岩型矿床矿化、蚀变规律以及物性特征，对三维重、磁、电反演断面，特别是电法反演断面进行分析，可以有效地查明隐伏成矿斑岩体。
（4）多宝山斑岩型铜矿床、争光岩金矿床、铜山铜矿床等皆为岩浆期后铜金热液型矿床，多宝山、铜山、争光隐伏斑岩体的发现证实了多宝山矿集区的矿床类型为斑岩型铜矿，并为探明深部矿体指明了方向。多宝山铜矿外围可以扩大部署，铜山铜矿矿体有往下延伸的趋势，铜山矿区斑岩体南侧可能存在大型矿体，争光岩金矿下部也存在斑岩型铜矿的可能，多宝山矿集区具“上金下铜（钼）”的成矿特征。
（5）钻孔揭示了铜山矿区推测的斑岩体左侧（南侧）存在隐伏矿体，与综合解释较为吻合，预示着多宝山矿集区的铜资源储量前景巨大。同时说明通过重、磁、电等综合方法查明隐伏斑岩体的思路在斑岩型矿床的勘查较为有效。
致谢研究的过程中得到了中国地质调查局沈阳地质调查中心朱群研究员级高工、邵军研究员级高工、杨晓平高工、赵院冬高工和有色金属矿产地质调查中心杨生教授级高工、中南大学戴世坤教授的帮助和指导，撰写本文得到原中国冶金地质勘查总局三局真允庆教授级高工等专家的指导，中国科学院地质与地球物理研究所王乐老师提供了铜山勘探线热液蚀变成果矢量图（图8），在此一并致谢！
猜你喜欢铜山闪长岩宝山黑太阳金山(2022年1期)2022-02-13程宝山书法作品选鸭绿江(2021年17期)2021-10-13甘肃省永昌地区熊子山岩体锆石U-Pb年龄及其地质意义矿产与地质(2020年5期)2020-12-23关于RCA改性剂在闪长岩上面层的应用研究中国科技纵横(2020年12期)2020-11-28宝山再生能源利用中心概念展示馆现代装饰(2020年11期)2020-11-27吉林中部早三叠世高镁闪长岩地球化学特征、岩石成因及地质意义黑龙江科技大学学报(2020年5期)2020-11-20南太行綦村岩体中辉石的成因矿物学研究河北地质大学学报(2020年1期)2020-04-26———南京市江宁区铜山小学科技教育掠影">让乡村小学生科学素养落地生根———南京市江宁区铜山小学科技教育掠影江苏教育(2019年78期)2019-11-21非遗文化的保护与发展的实然和使然艺术评鉴(2018年17期)2018-12-29系统论视域下学校体育整体发展的个案研究体育教学(2015年12期)2015-12-26