李建亭,刘雪敏,王学求,韩志轩,江瑶
(1.成都理工大学 地球科学学院,四川 成都 610059; 2.中国地质科学院 地球物理地球化学勘查研究所,河北 廊坊 065000)
通讯作者: 刘雪敏(1987-),女,博士,主要从事于勘查地球化学与矿床地球化学研究工作。Email: 451245437@qq.com
随着矿产勘查程度的提高,地表露头矿和浅部矿越来越少,寻找深部隐伏矿越来越受到矿产勘查界的重视,由此一系列找矿勘查方法应运而出,这主要包括电提取法、地气法、酶提取法、活动金属离子法、地球气纳微金属测量法、金属活动态测量法等多种非传统地球化学方法[1-6]。这些方法虽然原理不同,探测深度大(>100 m),探测信息微弱,但与矿化关系密切,因此,将它们统称为深穿透地球化学法[7]。上述方法主要通过在已知矿区开展试验,通过圈定的地表地球化学异常与深部隐伏矿体的空间对应关系来验证方法有效性,但这些深穿透地球化学方法并没有在未知的覆盖区开展广泛的找矿应用,因为目前国内外学术界在深穿透地球化学理论基础问题上还存在一定争议,特别是如何断定地表的地球化学异常就是来源于深部的隐伏矿体,有什么直接证据?要回答这个问题,实际上就是要示踪覆盖区地表圈定的地球化学异常的来源。
近年来,随着高分辨MC-ICP-MS技术的发展和同位素分析成本的降低,可借助同位素示踪技术来查证地表地球化学异常信息的来源。通常的做法是在覆盖区采集地表介质,如土壤、地下水、地表水、地气等,通过对比其同位素的含量和分布与矿石、围岩中同位素组成的关系,来示踪地表地球化学异常的来源。如汪明启等[8]、徐洋等[9]利用铅同位素分别在甘肃蛟龙掌铅锌矿区和山东邹平王家庄铜矿区成功示踪了异常区地气样品有深部矿石硫化物铅的加入,从而为地气用于覆盖区矿产勘查提供了可靠的理论基础。刘雪敏等[10]通过铅、硫同位素对新疆金窝子金矿床和内蒙古拜仁达坝—维拉斯托多金属矿床的异常区地气样品进行了示踪,提出了依托偏提取的思想,仅测量介质中活动态的铅、硫同位素的组成应更为合理。Saunders等[11]用Cu、S、Pb同位素在浅层低温Au-Ag矿区示踪了成矿流体的来源及其与构造的关系,为后期寻找深部的盲矿提供了依据。Mattew等[12]和Caritat等[13]分别利用矿区浅层地下水中的Pb、S、Sr、Pb同位素说明了地下水中的金属元素与深部矿体的关系。
本文在已知隐伏矿——罗卜岭斑岩型铜钼矿床上方采集了表层土壤以及典型钻孔中的矿石和围岩样品,分析6个微量元素(Cu、Mo、Ba、Pb、Zn、V)的含量变化以及硫、铅同位素的组成特征,来验证土壤金属活动态测量技术、微细粒级土壤全量测量技术在隐伏矿区的找矿效果,并根据铅、硫同位素的组成特征来识别地表地球化学异常的来源。
罗卜岭大型铜钼矿床位于福建省上杭县紫金山矿田的北东部(图1a),属于典型的斑岩型矿床。国内外对紫金山矿田的大地构造背景、地层、构造、岩浆活动已有大量文献发表[14-21],下文仅概述罗卜岭铜钼矿床的地质特征。
罗卜岭铜钼矿床位于宣和复背斜西南倾伏端东南翼,矿区出露的地层较为简单,仅在矿区南部出露新生界第四系地层和少量泥盆系天瓦岽组。该区大面积分布燕山期岩浆岩(图 1b),主要包括早白垩世四坊花岗闪长岩体和罗卜岭花岗闪长斑岩体, 也是该矿床的主要赋矿岩体,少量出露晚侏罗世五龙子中粒花岗岩体和才溪二长花岗岩体,在矿区深部还广泛发育早白垩世中寮似斑状花岗闪长岩体。矿区内断裂构造较为发育,以NE向和NW向断裂为主,近SN向次之,其中NE向断裂起控矿作用,NW向起控岩作用。另外,NE、NW向节理构造也十分发育,其中NE向节理为主要的容矿构造[23]。
1—早白垩世罗卜岭花岗闪长斑岩;2—早白垩世四坊花岗闪长岩; 3—早白垩世火山石帽山岩体;4—第四系冲积沉积物;5—晚侏罗世才溪岩体;6—中侏罗世紫金山复式岩体;7—白垩纪英安玢岩;8—白垩纪隐爆角砾岩;9—早石炭世岩体(林地组);10—晚泥盆世碎屑沉积岩体(天瓦岽组、桃子坑组);11—下震旦统楼子坝群;12—断层;13—264号测线; 14—经典矿床;15—矿化蚀变带(K-Phl:弱钾化-绢英岩化蚀变带,Chl-Phl:弱绿石化-绢英岩化蚀变带,Kl-Phy:高岭石化-黄铁绢英岩化,Ms-Di-Alu:白云岩化-地开石化-明矾石化);16—研究区;17—土壤采样点位;18-钻孔1—early Cretaceous Luoboling granodiorite porphyry;2—early Cretaceous Sifang granodiorite; 3—early Cretaceous Shimaoshan volcanic rocks;4—Quaternary alluvial sediments;5—late Jurassic Caixi monzogranite pluton;6—middle Jurassic Zijinshan granite batholith;7—Cretaceous dacitic porphyries;8—Cretaceous cryptoexplosive breccia pipes;9—early Cretaceous Lindi clasitic sediments;10—late Devonian clasitic sediments (Tianwadong & Taozikeng formations);11—Neoproterozoic Louziba metamorphosed clasitic sediments;12—fault;13—section view of exploration line 264;14—ore deposit;15—mineralization alteration zone (K-Phl: weakly potassium-sericitization alteration zone, Chl-Phl: weakly chlorogenic-sericitization alteration zone, Kl-Phy: kaolinitization-pyrite sericitization, Ms-Di-Alu: dolomitization-kazitritization-alum petrochemistry);16—the study area;17—soil sampling points;18—drill hole图1 福建紫金山矿田地质简图(a)及罗卜岭矿区地质简图(b)(修改自Zhong J等[22])Fig.1 Geological schematic diagram (a) of Zijinshan ore field in Fujian Province, geological schematic diagram (b) of Luoboling mining area (modified from Zhong J,et al.[22])
罗卜岭铜钼矿床主成矿元素为铜,伴生钼,在空间上呈“上铜下钼”分布。矿体平面上呈半圆弧展布,空间上呈马鞍状,以似层状、扁豆状和透镜状产出,主要赋存于中寮似斑状花岗闪长岩体外接触带的罗卜岭花岗闪长斑岩体和四坊花岗闪长岩体中,局部分布于五龙子中粒花岗岩体内,是典型的被岩浆岩覆盖的隐伏矿体。矿区热液蚀变强烈,按照蚀变矿物组合及其空间分布规律,罗卜岭铜钼矿床从内向外可划分出5个蚀变带:(弱)钾化—绢英岩化蚀变带(K-Phl)、(弱)绿泥石化—绢英岩化蚀变带(Chl-Phl)、高岭石化—黄铁绢英岩化蚀变带(Kl-Phy)、地开石—硅化蚀变带(Di-Q)和明矾石—地开石—硅化组合蚀变带(Di-Alu-Q),而铜钼矿体主要赋存于弱钾化—绢英岩化带内,部分位于绿泥石化—绢英岩化带内。主要的金属矿物有:黄铜矿、黄铁矿、辉钼矿,其次为铜蓝、斑铜矿、蓝辉铜矿、硫砷铜矿、方铅矿、闪锌矿等;非金属矿物以石英、白云母(绢云母)为主,其次为泥化蚀变矿物地开石、高岭石、蒙脱石、明矾石、叶蜡石,绿泥石—绢英岩化产物绿泥石、绢云母以及钾化蚀变矿物钾长石、黑云母组成。矿石结构主要有它形粒状结构、半自形—自形结构、填隙和充填结构、包含结构、交代残余结构和共结出结构。矿石构造较为简单,主要类型有浸染状、脉状、网脉状构造等[24-26]。
在罗卜岭斑岩型铜钼矿区地表采集了15件土壤样品、13件岩石(围岩和矿石)样品。具体采样方法如下:土壤样品沿264勘探线(位置见图1b)布置,该勘探线垂直于矿体走向,长约1.5 km,点距50~200 m不等,在矿区和围岩区共布置15个土壤采样点。每个采样点由3个子采样点组成,子采样点相距1.5~2 m,采集10~30 cm深的去掉腐殖层后的残坡积土并均匀混合为一件样品,每件样品采集1.5~2 kg。围岩、矿石采集于264勘探线ZKIX04、ZKVI03、ZKVII03钻孔内200 m以下的具有代表性的新鲜样品。
土壤样品风干、捏碎后,直接过200目(74 μm)孔径的筛子,仅对筛分后的微细粒级土壤进行研究。分析样品(全量、金属活动态)中测试元素(Cu、Mo、Ba、Pb、Zn、V)的含量变化以及全量样品的硫、铅同位素组成特征,其中,对于金属活动态样品仅分析了各元素的水提取态和黏土吸附态。采集的钻孔岩石样品在室内开展详细的矿物学研究后,挑选出矿石矿物较多的样品并粉碎至40~80目,在双倍镜下挑选出黄铁矿和黄铜矿。对围岩和挑选出的矿石(单矿物)样品分析了硫、铅同位素组成特征。测试的准确度和精确度利用实验室重复样和标准样品进行监测,分析结果可靠。具体分析实验室和测试方法见表1。
微细粒级土壤(<200目)中包含大量的粘土矿物、铁锰氧化物、盐类等物质,其吸附作用和可交换性能是金属元素的天然“捕获井”[27-28]。深部与成矿有关的物质在一种或多种营力(如离子扩散、电化学梯度、地下水溶解循环、地气流及地震泵等)的共同作用下,沿各种迁移通道卸载于地表,大部分被细粒级土壤捕获[29]。所以分析微细粒级土壤全量中的微量元素可以在一定程度上判别下伏是否存在隐伏矿体。如刘汉粮等[30]在甘肃柳元花牛山铅锌矿区107剖面线开展的微细粒级土壤全量测量实验,发现主成矿元素Pb、Zn的含量高值区与深部隐伏矿体相对应;韩志轩等[31]对江西通江岭铜矿区的19和20号剖面线进行了微细粒级土壤全量测量实验,发现Cu、Pb、Zn的含量高值区可以有效指示深部隐伏矿体。
表1 各介质的分析方法
在罗卜岭铜钼矿区264号测线上采集的15件微细粒级土壤全量样品的6个微量元素(Cu、Mo、Ba、Pb、Zn、V)统计参数见表2。由表2可知,6个元素的变异系数均低于1,如主成矿元素Cu、Mo的变异系数仅分别为0.63、0.85。这是因为一方面本矿床呈浸染状分布于斑岩体内,另一方面该矿床伴生较强的蚀变作用,与成矿有关的元素会进入蚀变带,造成成矿元素异常面积较大且分布均匀[33]。主成矿元素Cu、Mo的平均值明显高于福建省土壤背景值,富集系数均大于2,分别为4.44、2.79,显示出强富集特征;Ba的平均值与福建省土壤背景值相差不大,富集系数为1.13;V的平均值略低于福建省土壤背景值,富集系数为0.71。值得注意的是,伴生元素Pb、Zn富集特征具有明显的不同,Pb明显富集,Zn明显贫化,富集系数分别为3.17、0.35,这与一般热液型矿床中Pb和Zn密切共生的特征不符,可能Pb、Zn在地表土壤中已经发生了解耦,但是,这种解耦是由成矿过程中元素分异导致的,还是与表生成土过程中锌的强烈淋失相关,仍有待进一步研究。整体来看,矿区土壤微量元素总体有成矿元素变异系数低、富集系数高的特征。另外,Cu、Mo、Ba的极值比大于10,分别为10.1、24.28、10.24,反映出其在此剖面分布的均匀性较差,下方存在一定矿化,与已知事实相符。
图2为罗卜岭铜钼矿区264号勘探线微细粒级土壤6个微量元素的全量折线图。由图2可知,主成矿元素Cu、Mo从背景区到已知异常区含量明显升高,但最高值均出现在测线东南端背景区的14号点位,分别为120×10-6、46.14×10-6。另外,根据矿体形态,矿体右侧应该还有延伸。综合地表异常及矿体形态,推测264线南侧深部仍有成矿潜力,但成矿深度可能比已知矿体略浅。故本文将14~15号采样点位视为预测异常区。整体来看,主成矿元素Cu、Mo的含量高值区(含量大于120×10-6、10×10-6)与下伏隐伏矿体的分布较为吻合;Ba的含量变化趋势与Cu较为相似且大于300×10-6的范围与下伏隐伏矿体的分布高度吻合;V、Pb、Zn含量高值区的分布与下伏矿体的展布相关性较差,这可能与上文提及到的Pb、Zn在地表土壤的解耦有关,也可能与Pb、Zn活动性较强,易向矿体边缘富集有关[35]。值得注意的是,V、Pb、Zn的含量高值范围与赋矿岩体——罗卜岭花岗闪长斑岩接近地表的范围较为一致,这与赋矿岩体普遍发育矿化蚀变的地质事实相符。综上所述,地表微细粒级土壤全量中的Cu、Mo、Ba的含量高值区与深部隐伏矿体的展布相关性较强,对于预测深部矿体的规模有一定的指示意义,而V、Pb、Zn可以较为准确地圈定出接近地表矿化岩体的范围。
金属活动态测量的理论依据是深部与成矿有关的物质在一种或多种营力作用下,沿迁移通道到达地表并转化成各种活动态(超微细的亚微米至纳米级颗粒、胶体、离子和各种配合物),被可溶性盐类、黏土、氧化物和有机质所吸附或结合而存在于地表疏松介质中,在原介质元素含量的基础上形成活动态叠加含量,使用适当的提取剂将这些元素叠加含量提取出来,获取来自深部的矿化信息,从而达到寻找和评价隐伏矿的目的[7]。如宋雷鹰[36]在内蒙古哈如勒敖包矿区对成矿元素进行水提取态、黏土吸附态实验,发现贱金属元素如Cu、Pb、Zn、Ag等的水提取态对地质找矿具有较好的指示意义;杨刚刚等[37]在新疆东戈壁钼矿区进行的金属活动态测量剖面实验显示Mo、W的水提取态和Cu的铁锰氧化态在矿体上方有明显的连续高值异常区,并在主矿体上方达到最大,对深部隐伏矿体具有良好指示意义。
表2 罗卜岭铜钼矿区264号测线微细粒级土壤6个微量元素全量统计参数(n=15)
图2 福建罗卜岭铜钼矿区264号勘探线微细粒级土壤6个微量元素的全量折线图(修改自赵辰等[34])Fig.2 Line graphs of six trace elements’ concentration in the fine grained soils of the 264 exploration line in the Luoboling ore area (modified from Zhao C,et al[34])
图3 罗卜岭铜钼矿区264号勘探线微细粒级土壤6个微量元素的金属活动态折线图Fig.3 Line diagrams of active states for 6 trace elements in fine soils of 264 exploration line in Luoboling Cu-Mo mining area
罗卜岭铜钼矿区264号测线上采集的15件微细粒级土壤活动态样品的6个微量元素(Cu、Mo、Ba、Pb、Zn、V)统计参数见表3。由表3可知,主成矿元素Cu及伴生元素Pb、Zn的金属活动态(水提取态+黏土吸附态)的平均值分别占微细粒级土壤全量平均值的5.26%、5.48%、7.12%,说明水提取态和黏土吸附态是这些元素的重要存在形式。如活动态Cu的平均值为5.28×10-6,最大值可达17.37×10-6、活动态Pb平均值为6.06×10-6,最大值可达20.75×10-6,显然研究区局部受到了矿化影响。
图3为罗卜岭铜钼矿区264号勘探线微细粒级土壤6个微量元素的金属活动态折线图。由图3可知,6个微量元素的金属活动态在15号点位均存在较高异常,这与上文微细粒级土壤全量所预测的下伏矿体相对应。主成矿元素Cu在已知异常区上方仅12号点位存在异常,最高值(17.37×10-6)出现在15号点位;Ba与Mo在含量变化趋势上具有相似性,两者均在推测矿体上方14、15号点位显示出明显异常,最高值分别为2.28×10-6、278×10-9,但在已知异常区上方均未有异常显示;相较于土壤全量,活动态的V、Pb、Zn对于圈定接近地表的矿化岩体范围更为清晰。整体而言,本次各测试元素在隐伏矿体上方均有较好的指示和预测效果,V、Pb、Zn对于圈定接近地表的矿化岩体范围更为清晰,但相比于土壤全量,Cu、Mo、Ba的高值区与下伏隐伏矿体的展布相关性较差。这可能是因为本次研究采用的活动态提取方法用水提取态和黏土吸附来代表金属活动态,而本矿区处于中亚热带气候,雨量充沛,地表植被较多,为湿润气候森林覆盖区,测试元素可能在有机态和铁锰氧化物结合态中含量较多,所以是本次研究各测试元素的活动态与下伏隐伏矿体展布相关性较差的主要原因。建议未来针对不同的地貌区发展不同的活动态提取技术。
表3 罗卜岭铜钼矿区264号测线微细粒级土壤金属活动态含量统计参数(n=15)
属活动态=水提取态+黏土吸附态。
从理论上讲,可以利用地表土壤中的硫同位素组成来示踪和判别深部隐伏矿体。因为大多数金属矿床的含矿矿物为硫化物(如本矿床的黄铜矿、黄铁矿等),矿体周围的(S2)2-或S2-浓度高,可以在地气流的作用下迁移至地表土壤中,而自然界表生环境中存在着一类细菌(如Thiobacilli、Beggiatoa和Sulfolobus等)[38],能将低价态硫化物氧化成高价态的化合物(硫酸盐或亚硫酸盐)或者 (S2)2-/S2-直接被氧气氧化为亚硫酸盐或硫酸盐,形成的高价态化合物往往富集重的硫同位素,分馏约5%[39],故矿体上方土壤的硫同位素组成应与背景区存在一定的差异。
本次研究共获得15件微细粒级土壤全量样品、5件围岩(非赋矿围岩)样品、7件单矿物样品的硫同位素组成,此外还收集了前人在该矿床获得的硫同位素数据(表4),并绘制了硫同位素分布的频率直方图(图4)。由表和图可知:19件矿石样品的硫同位素分布比较集中,范围为-1.6‰~2.2‰,平均值为1.3‰,矿石硫同位素组成落在陨石硫同位素范围之内(0±3)‰,具有地幔来源特征[41]。其中1件辉钼矿硫同位素为2.2‰、1件闪锌矿硫同位素为-1.6‰、1件方铅矿硫同位素为-1.2‰、16件黄铁矿硫同位素组成范围为0.5‰~2.2‰,呈塔式分布,峰值在1‰~2‰之间;5件围岩硫同位素组成范围为6.2‰~16.6‰,分布范围较宽,主要集中在10‰~12‰之间; 5件背景区微细粒级土壤全量样品的硫同位素变化范围为-0.4‰~10.3‰,范围较宽,主要集中在4‰~7‰之间;10件异常区微细粒级土壤全量样品硫同位素变化范围最大,其范围为1.8‰~15.3‰,主要有两个集中范围:5‰~9‰、12‰~16‰。可以看出,从矿石→背景区土壤→围岩→异常区土壤,硫同位素的组成逐渐增大,且除矿石硫同位素组成较为集中外,其他介质中的硫同位素组成变化范围较宽,组成较为分散。
表4 罗卜岭铜钼矿区细粒级土壤、矿石、围岩样品的硫同位素数据(n=39)
图4 罗卜岭铜钼矿区矿石、围岩、土壤背景和土壤异常的硫同位素组成Fig.4 Sulfur isotopic composition diagrams of ores, rocks, soil background and soil anomalies in the Luoboling Cu-Mo mining area
本文利用硫同位素识别土壤中异常物质来源的基本思想是:对比异常区土壤样品和矿床周围不同介质(如矿石、围岩、土壤背景)的硫同位素组成特征,来探讨土壤中异常物质的可能来源。通过上文的讨论可以看出,研究区不同介质的硫同位素组成存在一定差异,这有利于使用硫同位素示踪法研究土壤中异常物质的来源。但首先要保证异常区与背景区的土壤有明显差异,这可以利用SPSS软件进行单因素方差分析来检验两者数据是否存在差异(当P值小于0.05时,存在明显差异)。方差分析结果如表5所示,土壤异常区与土壤背景区两组数据的P值为0.08,大于0.05,表明两组数据之间的δ34S 没有明显差别。从统计学的角度来看,这两组数据不能进一步分析。
由图4可知,土壤异常区硫同位素的组成集中区(5‰~9‰、12‰~16‰)与围岩(非赋矿围岩)硫同位素的组成集中区(10‰~12‰)特征相似,这与方差分析得出的结论基本一致,表明土壤异常区的硫同位素组成大多数信息继承自非赋矿围岩,而不能反映隐伏矿的硫同位素组成。另外,由表4、图5可知,5件背景区土壤样品δ34S的平均值为 5.20‰,10件异常区样品δ34S的平均值为9.77‰,说明异常区土壤确实明显受到深部矿体的(S2)2-/S2-运移至地表被氧化而富集重同位素的影响。本次研究的是土壤全量中硫同位素组成,受外界环境较大,若采用刘雪敏等[10]提出的依托偏提取的思想,仅测量土壤活动态中的硫同位素组成,可能更合理。
铅同位素之间的质量差相对较小,任何物理化学条件引起的铅同位素分馏均可以忽略不计,引起铅同位素组成变化的主要原因是放射性U和Th的衰变。方铅矿、闪锌矿、黄铁矿等硫化物中不含U、Th或含量极低,与矿物中铅的质量分数相比可忽略不计,且矿物形成后不再有放射性成因铅的明显加入。那么这些硫化物中的Pb同位素组成就可以反映原始热液中金属物质来源区的Pb体系及其初始Pb同位素组成特征,从而获得有关成矿物质来源的信息。铅同位素在成土作用过程中几乎不发生分馏[8],所以异常区土壤和背景区土壤208Pb/204Pb、207Pb/204Pb、206Pb/204Pb比值可能会出现明显差异,土壤的铅同位素组成能示踪异常物质的来源,可用于隐伏矿体的示踪。
图5 罗卜岭矿区264号勘探线土壤微细粒级全量中硫同位素组成分布Fig.5 Sulfur isotope distribution of total soil fine particles along exploration line 264 of Luoboling mining area
罗卜岭铜钼矿床各介质铅同位素数据(本次研究测试的数据和前人的数据)见表6。由表6可知,矿石208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的范围分别为38.67~39.215、15.595~15.674、18.417~18.764,变化范围较小,分布相对集中。在208Pb/204Pb-206Pb/204Pb(图6a)、207Pb/204Pb-206Pb/204Pb(图6b)图解中投点落于造山带附近,具有深部岩浆来源特征;围岩(非赋矿围岩)的3种铅同位素比值范围分别为38.721~39.854、15.643~15.709、18.525~19.86,相比矿石样品,208Pb/204Pb、206Pb/204Pb比值明显增高。在208Pb/204Pb-206Pb/204Pb(图6a)图解中投点落于造山带附近,具有深部岩浆来源特征,在207Pb/204Pb-206Pb/204Pb(图6b)图解中投点大多落于造山带附近,个别投点落在上地壳附近,亦显示出深部岩浆来源特征。土壤铅同位素整体介于矿石与围岩之间,相比矿石样品,具有更小变化范围,分布更为集中,其中异常区土壤208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的范围分别为38.694~38.804、15.626~15.655、18.48~18.607,背景区土壤208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的范围分别为38.787~38.846、15.639~15.661、18.555~18.645,异常区土壤3种铅同位素组成均略低于背景区土壤样品。在208Pb/204Pb-206Pb/204Pb(图6a)、207Pb/204Pb-206Pb/204Pb(图6b)图解中,土壤样品投点落于造山带附近,同样表现出具有深部岩浆来源特征。区别在于异常区土壤投点明显更接近矿石投点分布范围,而背景区土壤投点更靠近围岩投点分布范围。
表5 罗卜岭矿区不同介质硫同位素数据单因素方差分析P值
表6 罗卜岭铜钼矿区细粒级土壤、矿石、围岩样品的铅同位素数据(n=52)
首先利用SPSS软件对不同介质的3种铅同位素数据进行单因素方差分析,来检验数据间的差异性,结果见表7。异常区土壤和背景区土壤208Pb/204Pb、206Pb/204Pb的P值分别为0.01、0,均小于0.05,存在显著差异,表明两者2种铅同位素组成特征明显不同。从统计学的角度来看,这两组数据可以进一步分析。而207Pb/204Pb的P值为0.16,大于0.05,不存在显著差异,表明两者铅同位素组成特征相似。从统计学的角度来看,这两组数据不能进一步分析;矿石样品与围岩样品207Pb/204Pb的P值为0.2,大于0.05,不存在显著差异,表明两者铅同位素组成特征相似,围岩样品有矿石铅同位素的参与;异常区土壤样品与矿石样品208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的P值分别为0.17、0.07、0.66,均大于0.05,不存在显著差异,表明两者3种铅同位素组成特征相似,异常区土壤样品继承了矿石的铅同位素特征。异常区土壤样品与围岩样品208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的P值分别为0.02、0.01、0.01,均小于0.05,均存在显著差异,表明两者3种铅同位素组成特征明显不同;背景区土壤样品与矿石样品206Pb/204Pb的P值为0.01,小于0.05,存在显著差异,而两者208Pb/204Pb、207Pb/204Pb的P值分别为0.97、0.67,均大于0.05,不存在显著差异,背景区土壤样品可能部分继承了矿石铅同位素特征。背景区土壤样品与围岩样品208Pb/204Pb、207Pb/204Pb、206Pb/204Pb的P值分别为0.16、0.2、0.16,均大于0.05,不存在显著差异,表明两者3种铅同位素组成特征相似。
图6 罗卜岭铜钼矿区矿石、围岩、土壤背景和土壤异常208Pb/204Pb-206Pb/204Pb(a)、207Pb/204Pb-206Pb/204Pb(b)图解Fig.6 208Pb/204Pb vs.206Pb/204Pb(a)、207Pb/204Pb vs.206Pb/204Pb(b) composition diagrams of ores, wall rocks, soil background and soil anomalies in the Luoboling Cu-Mo mining area
表7 罗卜岭矿区不同介质铅同位素数据单因素方差分析P值
在207Pb/206Pb-206Pb/204Pb图解(图7)上,矿石、围岩、土壤样品投点近似呈一条直线分布且土壤样品投点介于矿石、围岩之间,异常区土壤样品投点更接近矿石分布范围而背景区土壤样品投点更接近围岩分布范围,表明该矿区土壤铅同位素的来源几乎没有外来干扰,而是继承了矿区矿石、围岩铅同位素特征。矿石、围岩部分样品投点几乎重合,这是因为本次采集样品除了中寮似斑状花岗闪长岩、才溪二长花岗岩等非赋矿围岩外,还有少量成矿母岩——罗卜岭花岗闪长斑岩。以上结果与方差分析得出的结论一致。
图7 罗卜岭矿区不同介质207Pb/204Pb-206Pb/204Pb图解Fig.7 Diagram of 207Pb/204Pb vs.206Pb/204Pb for various media in Luoboing mining area
从背景区土壤到异常区土壤,206Pb/204Pb、207Pb/204Pb、208Pb/204Pb值均有降低的趋势(图8),206Pb/204Pb值在异常区的最大值(18.540)比在背景区的最小值(18.555)还低,208Pb/204Pb值在异常区的最大值(38.788)与背景区的最小值(38.787)非常接近,但207Pb/204Pb值在异常区虽有降低的趋势,但这种趋势还不足以造成显著的差异。另外,由上文方差分析可知,背景区土壤样品与矿石样品208Pb/204Pb、207Pb/204Pb的P值大于0.05, 不存在显著差异,背景区土壤样品可能部分继承了矿石铅同位素特征,而206Pb/204Pb的P值小于0.05,二者不存在显著差异,背景区土壤样品没有矿石铅同位素特征。因此,综合来看,206Pb/204Pb在地表土壤中的变化与下伏隐伏矿体的分布吻合度最好,可以有效指示和示踪深部隐伏矿体。
图8 罗卜岭矿区264勘探线微细粒级土壤全量中铅同位素组成分布Fig.8 Lead isotope distribution of total soil fine particles along exploration line 264 of Luoboling mining area
本文在已知隐伏矿罗卜岭斑岩型铜钼矿床上方采集了表层土壤以及典型钻孔中的矿石和围岩样品,分析了6个微量元素(Cu、Mo、Ba、Pb、Zn、V)的含量变化以及硫、铅同位素的组成特征,来验证土壤金属活动态测量技术、微细粒级土壤全量测量技术在隐伏矿区的找矿效果,并根据铅、硫同位素的组成特征来识别地表矿致异常的来源。研究表明:
1) 地表微细粒级土壤全量中Cu、Ba、Mo的含量高值区与深部隐伏矿体的展布相关性较强,对于预测深部矿体的规模有一定的指示意义,而V、Pb、Zn可以较为准确地圈定出接近地表矿化岩体的范围。
2) 结合地表异常和矿体形态规模,土壤活动态、微细粒级土壤全量均显示出14、15号采样点下方极有可能存在着隐伏矿体。两种测量方法均可以根据V、Pb、Zn的含量变化较为准确地圈定出接近地表矿化岩体的范围。相比于微细粒级土壤全量,活动态的Cu、Ba、Mo的高值区与下伏隐伏矿体的展布相关性较差,这可能是因为本次研究采用的活动态提取方法仅包括金属的水溶态和黏土吸附态,但本矿区处于中亚热带气候,土壤富含大量的有机质和铁锰氧化物,活动态的金属可能主要以有机相态和铁锰氧化物相态存在,所以指示隐伏矿体分布的效果较差。建议针对不同的地貌特征发展不同的活动态提取技术。
3) 从矿石→背景区土壤→围岩→异常区土壤,硫同位素的组成逐渐增大,且除矿石硫同位素组成较为集中外,其他介质中的硫同位素组成变化范围较宽,组成较为分散。异常区土壤全量样品δ34S的平均值明显高于背景区,说明异常区土壤全量样品确实明显受到深部矿体中 (S2)2-/S2-运移至地表被氧化而富集重同位素的影响,但土壤异常区、围岩两者的硫同位素组成集中区特征相似,表明土壤异常区的硫同位素组成大多数信息继承自非赋矿围岩,而不能反映隐伏矿的硫同位素组成,建议测量活动态中的硫同位素组成更为合理。
4) 单因素方差分析显示异常区土壤全量样品206Pb/204Pb、207Pb/204Pb、207Pb/204Pb的比值与矿石不存在显著差异,异常区土壤样品继承了矿石的铅同位素特征,从而为微细粒级土壤全量测量用于覆盖区矿产勘查提供了理论基础。背景区土壤全量样品仅206Pb/204Pb的比值与矿石存在显著差异且206Pb/204Pb在地表土壤中的变化与下伏隐伏矿体的分布的吻合度最好,可以有效指示深部隐伏矿体。