Characteristics of soil geochemical anomalies and metallogenic prediction in Liukuaidi South of Mulan County, Heilongjiang Province
-
摘要:
黑龙江木兰县六块地南大地构造位置属于小兴安岭-张广才岭岩浆弧，成矿区带为伊春-延寿多金属成矿亚带的重要组成部分，成矿地质条件优越。为查明该区各成矿元素和矿（化）体的分布特征，实现找矿突破，开展了1：2万土壤地球化学测量及成矿预测工作。利用12种成矿元素进行相关性分析、聚类分析、因子分析、单元素异常分析和组合元素异常分析，圈定了元素异常区，优选了找矿靶区。结果表明，区内Au、As、Sb、Mo元素变异系数值高，成矿潜力较强；As、Sb元素的相关性较好，相关系数为0.73，R型聚类分析将成矿元素分为4类，因子分析将分析元素分为4组；结合区内成矿地质条件和土壤地球化学异常特征，共圈出单元素异常113处，组合异常3处；测区主攻矿种为Au和Zn，Au元素主要找矿靶区为Au-5和Au-7异常区，Zn元素主要找矿靶区为Zn-6异常区，As、Sb、Mo找矿前景非常好，可作为第一找矿目标元素。
Abstract:The study area Liukuaidi South is located in the central and eastern part of Heilongjiang Province.The tectonic position belongs to the Xiaoxinganling-Zhangguangcailing magmatic arc.The metallogenic zone is an important part of the Yichun-Yanshou polymetallic metallogenic subbelt, and the metallogenic geological conditions are superior.To find out the distribution characteristics of ore-forming elements and ore-forming bodies in the study area and achieve the ore-prospecting breakthrough, a 1:20000 soil geochemical survey and ore-forming prediction were carried out.The correlation analysis, cluster analysis, factor analysis, single element anomaly analysis, and combined element anomaly analysis of 12 metallogenic elements are used to delineate the element anomaly area and optimize the prospecting target area.The results show that the variation coefficients of Au, As, Sb, and Mo are high and the metallogenic potential is strong.The correlation coefficient of As and Sb elements was 0.73.The metallogenic elements were divided into four groups by R-type cluster analysis and four groups by factor analysis.Combined with the ore-forming geological conditions and soil geochemical anomaly characteristics, 113 single-element anomalies and 3 combined anomalies were identified.The main ore types in the test area are Au and Zn, the main Au prospecting target areas are Au-5 and Au-7 abnormal areas, and the main Zn prospecting target areas are Zn-6 abnormal areas.Sb and Mo have very good prospecting prospects and can be used as the first prospecting target elements.
-
加甘滩金矿床位于甘肃省甘南州夏河-合作-临潭成矿带上，20世纪90年代以来，在该地区相继发现了以早子沟特大型金矿、加甘滩特大型金矿、地南大型金矿等为代表的大批金矿床(点)，探明总资源储量逾450 t，夏河—临潭地区已成为西秦岭地区最重要的金矿富集区。加甘滩金矿自1997年发现以来至2016年，基于大规模钻探、槽探等勘查工作，累计查明金资源量153.7 t，矿床平均品位2.76 g/t，成为中国西北地区备案资源量最大的金矿床。但加甘滩金矿找矿突破时间较短，研究程度极薄弱，只有部分勘查单位根据勘查实践从矿床地质特征、围岩蚀变、控矿因素、找矿标志、金的赋存状态等方面进行了基础性研究，并对矿床成因进行了初步探讨(张江苏等，2007；田向盛，2014；2016；代文军等，2016；白云等，2017)，而在成矿物质来源、热液流体运移沉淀、成矿时代等方面的研究几乎为空白。目前一般将西秦岭地区的主要金矿床划分为造山型和类卡林型2类(Mao et al., 2002；陈衍景等, 2004；Deng et al., 2016；Goldfarb et al., 2019；Qiu et al., 2020；Yu et al., 2020a, b)，或者造山型、卡林型和类卡林型3类(Liu et al., 2015)。近年，一些学者认为，夏河—合作地区一些金矿床的成矿作用与岩浆热液有关，因而将区内的早子沟、录斗艘、以地南等金矿床划分为与岩浆岩有关的金矿床(Jin et al., 2017；Sui et al., 2017；2018)。流体包裹体是研究成矿作用的“探针”(陈衍景等，2007)，封存了成矿流体的活动记录。稳定同位素通常是地质流体和某些热液矿物的主要组分，可以作为确定物质来源的示踪剂，在矿床研究中具有重要作用。本文在详细矿床地质调查研究的基础上，通过对加甘滩金矿床主成矿阶段的石英中流体包裹体、氢-氧同位素测试分析，结合矿石中黄铁矿、毒砂硫、铅同位素特征，揭示成矿流体来源，探讨矿床成因和成矿机制。
1. 区域地质背景
加甘滩金矿床所在的夏河-合作成矿带位于西秦岭造山带的西北部(图 1-a)，秦岭-祁连-昆仑中央造山带西延，处于中秦岭早古生代陆缘褶皱带和南秦岭陆表海盆的结合部位。其基本构造格架定型于三叠纪以来的印支造山运动，地质构造复杂，岩浆活动频发，内生成矿作用强烈。区域出露地层从老到新依次有泥盆系(D)、石炭系(C)、二叠系(P)、三叠系(T)，其中泥盆系—二叠系分布在夏河-临潭区域断裂以北，呈背斜构造，主要岩性组合为灰绿色石英砂岩、长石石英砂岩、炭质板岩、钙质板岩、粉砂岩、生物碎屑灰岩、灰质砾岩等，主体为一套滨浅海相陆源碎屑沉积建造。三叠系主要分布在夏河-临潭区域断裂以南，岩性为灰绿色—深灰色长石杂砂岩、泥质板岩、粉砂质板岩，局部夹泥晶灰岩，为一套斜坡相浊积岩建造。研究区北部侵入岩发育，夏河二长花岗岩-老豆花岗闪长岩-美武花岗闪长岩等主要侵入体呈NW向弧形展布。喷出岩主要分布在甘南州东北部录斗艘金矿一带，岩性组合为凝灰岩、火山角砾岩、流纹岩、英安岩、安山岩等。夏河-临潭NW向区域断裂带(组)控制着区域矿床(点)的分布，NE向断裂及其与NW向断裂的交会部位往往是赋矿的有利空间，著名的矿床有加甘滩金矿(特大型)、早子沟金矿(特大型)、以地南铜金矿(大型)、格娄昂金矿(大型)、早仁道金矿、德乌鲁金铜矿、下拉地铅锌矿等。
![]()
1—第四系；2—上三叠统华日组；3—中上三叠统；4—中下三叠统；5—二叠系；6—石炭系；7—泥盆系；8—晚三叠世花岗闪长岩；9—砂岩；10—板岩；11—砂岩夹粉砂岩；12—砂岩夹板岩；13—闪长岩脉；14—花岗闪长岩脉；15—断裂；16—正断层；17—逆断层；18—俯冲方向；19—断层运动方向；20—总应力方向；21—剖面位置；22—金矿体；23—金矿床(点)；24—铜矿床(点)；25—锑矿床(点)；26—砷矿床(点)；27—铅锌矿床(点)；AMS—阿尼玛卿古缝合带；CBS—柴达木北缘古缝合带；NQL—北祁连缝合带；MLS—勉略古缝合带；SDS—商丹古缝合带；GL—甘孜-理塘缝合带Figure 1. Sketch geological map of the geotectonic location (a), Xiahe-Hezuo regional geology (b) and the Jiagantan gold deposit(c)(Li et al., 2019; 2020; Liu et al., 2018)2. 矿床地质特征
2.1 矿区地质
矿区出露地层为中下三叠统和第四系，中下三叠统总体呈NNW—NW向展布，从上到下可分为4个岩性段，依次为：①灰色—浅灰绿色长石石英砂岩夹浅灰色粉砂质板岩、灰色粉砂岩段；②浅灰色粉砂质板岩段；③灰色—浅灰绿色长石石英砂岩夹浅灰色粉砂质板岩段；④灰绿色岩屑砂岩段。地层变形较强，揉皱、褶皱、节理、软沉积变形等均较发育。NW向断裂为区域夏河-临潭断裂带的组成部分，EW向断裂规模较小，形成较晚，切割NW向断裂。其中F1断裂是矿区主控断裂，产状210°~250°∠40°~60°，断层产状上缓下陡。断层破碎带最宽处约350 m，最窄处约180 m。F1断层控制着金矿化的分布、规模及产状变化(白云等，2017)，断层产状具有波状变化，控制的矿体形态随之膨大缩小、尖灭再现、分支复合。褶皱构造组合较复杂，矿区北西至中部地层走向为NW—SE向，倾角陡倾、较缓均有，倾向亦见反向变化，产状不协调处均为断层接触，表明矿区内存在规模较大的次一级NW向褶皱，但被NW向断裂严重破坏。矿区南东部位则发育相对连续的褶皱组合，表现为宽向斜，紧背斜，轴迹为近NS向，向N倾伏(倾伏角二十余度)，该组褶皱与矿区中部至北西部的构造形迹在形态和方位上明显不一致，二者呈断层接触。围岩蚀变主要有毒砂化、黄铁矿化、辉锑矿化、硅化、绢云母化、赤铁矿化、褐铁矿化、碳酸盐化等。
2.2 矿床地质
矿区共圈出金矿体316条，均受F1及其次级断裂控制，形成庞大的NW向矿体群。主矿体共16条，产出标高3310~2400 m，延长520~1560 m，厚度2.76~12.82 m，Au平均品位2.13×10-6~3.65×10-6，矿体呈脉状、板状、分叉状等，倾向220°~245°，倾角35°~45°(图 2)。
![]()
图 2 加甘滩金矿0线勘查线剖面图1—金矿体；2—锑矿体；3—钻孔Figure 2. No. 0 prospecting profile of the Jiagantan gold deposit加甘滩金矿以金为主，伴生少量锑、砷。矿石类型以石英脉型、构造蚀变岩型为主，见少量辉锑矿化板岩型、毒砂-黄铁矿-白云石化长石砂岩型等。矿石中金属矿物主要为毒砂、黄铁矿、辉锑矿、自然金，氧化矿石中见赤铁矿、褐铁矿等，非金属矿物有长石、石英、白云母、方解石等。矿石构造有角砾状、浸染状、细脉浸染状、纹层状、块状等，矿石结构主要有碎裂结构、假象交代结构、交代结构、胶状结构、纤维状结构等。
矿石矿物主要有黄铁矿、毒砂、自然金、赤铁矿、褐铁矿、辉锑矿、磁铁矿等。脉石矿物主要有石英、长石、绢云母、少量方解石等(图 3)。根据矿物嵌布共生组合关系，将加甘滩金矿的成矿过程划分为2期4阶段：热液期，包括黄铁矿-毒砂阶段(早期硅化，普遍发育浸染状黄铁矿-毒砂)、石英-多金属硫化物阶段(石英-黄铁矿-毒砂-自然金-磁铁矿-辉锑矿等)和石英-碳酸盐阶段(形成无矿纯白色石英-方解石脉)；表生期，黄铁矿、毒砂、辉锑矿等硫化物在(近)地表发生氧化，形成赤铁矿、褐铁矿、锑华等。
![]()
图 3 加甘滩金矿气液两相包裹体镜下照片L—液态；V—气态Figure 3. Microscopic picture of gas-liquid two-phase inclusions in the Jiagantan gold deposit3. 流体包裹体特征
本次研究的包裹体样品采自热液成矿期主成矿阶段的含矿石英脉，流体包裹体样品测试分析在核工业地质分析测试研究中心完成，冷热台型号为Linkam THMS600型，检测精度为±0.1℃，测温范围-196~600℃。
加甘滩金矿区石英中流体包裹体数量不多，孤立分布，偶见成群分布。形态可见椭球状、长条状、三角形及不规则状，大小4~20 μm。包裹体类型主要为气液两相，可进一步分为富气相和富液相，其中富液相两相包裹体占比较大，广泛发育，气液比一般介于10%~20%之间，偶见单相液态包裹体(图 3)。
![]()
图版Ⅰa~h.加甘滩金矿床典型矿石照片。a. 采坑；b.错断的石英大脉；c. 矿石中发育石英细脉；d. 毒砂交代黄铁矿；e、f. 他形黄铁矿；g. 细粒黄铁矿和毒砂；h. 被交代溶蚀的黄铁矿。Py—黄铁矿；Apy—毒砂图版Ⅰ.在岩相学研究的基础上，对加甘滩金矿床热液期第二阶段石英中的流体包裹体进行显微测温，获得42个均一温度和冰点温度(表 1)，均一温度为132.7~365.9℃(图 4-a)，平均值为248.7℃，整体显示为中温特征。冰点温度变化范围为-0.3~-5.1℃，根据H2O-NaCl体系盐度-冰点公式(Potter et al., 1978)，计算其盐度区间为0.53%~7.85% NaCl，主要集中在0~6.0% NaCl(图 4-b)，平均值为3.78%NaCl，显示成矿流体具有低盐度的性质。
表 1 加甘滩金矿矿石中石英流体包裹体均一温度、冰点与盐度Table 1. Homogenization temperature, freezing point and salinity of Jiagantan gold deposit序号 矿物 包体类型 大小/ μm 气液比 均一温度/℃ 冰点温度/℃ 盐度/%NaCl 1 石英 V+L 6×6 20 329.3 -1.7 2.90 2 石英 V+L 9×6 35 342.5 -2.1 3.54 3 石英 V+L 5×4 10 262.7 -2.9 4.78 4 石英 V+L 8×3 25 361.5 -1 1.74 5 石英 V+L 8×2 12 279.8 -2.1 3.54 6 石英 V+L 6×3 20 340.9 -2.9 4.78 7 石英 V+L 8×5 8 214.7 -0.9 1.57 8 石英 V+L 6×4 28 355.2 -1.8 3.06 9 石英 V+L 8×4 8 224.9 -3.2 5.23 10 石英 V+L 9× 8 8 163.7 -1.3 2.24 11 石英 V+L 16×12 5 173.6 -2.4 4.01 12 石英 V+L 10×6 10 181.7 -2.8 4.63 13 石英 V+L 9×5 8 167.4 -0.3 0.53 14 石英 V+L 6×4 5 176.2 -1.4 2.40 15 石英 V+L 6×5 3 180.3 -2.8 4.63 16 石英 V+L 10×3 12 216.6 -3.5 5.66 17 石英 V+L 7×2 l5 326.3 -1.3 2.24 18 石英 V+L 8×4 10 349.6 -3.7 5.95 19 石英 V+L 5×4 15 278.2 -5.1 7.85 20 石英 V+L 7×2 22 292.3 -2.7 4.47 21 石英 V+L 6×2 28 264.7 -0.6 1.05 22 石英 V+L 8×4 12 266.4 -2.3 3.85 23 石英 V+L 8×6 25 328.3 -0.7 1.22 24 石英 V+L 6×2 32 365.9 -1.4 2.40 25 石英 V+L 9×2 l8 286.7 -0.9 1.57 26 石英 V+L 9×5 17 279.9 -2.6 4.32 27 石英 V+L 16×12 15 259.2 -2 3.38 28 石英 V+L 8×5 9 250.7 -2.9 4.78 29 石英 V+L 8×7 15 344.6 -1.3 2.24 30 石英 V+L 14×7 8 201.3 -2.4 4.01 31 石英 V+L 8×6 8 219.3 -3.6 5.81 32 石英 V+L 10×6 10 201.7 -2 3.38 33 石英 V+L 5×3 1 261.3 -1.8 3.06 34 石英 V+L 9×5 12 273.4 -1.7 2.90 35 石英 V+L 8×3 10 247.4 -2.8 4.63 36 石英 V+L 9×6 15 263.9 -4.1 6.52 37 石英 V+L 8×5 5 132.7 -2.1 3.54 38 石英 V+L 11×8 12 179.4 -2.9 4.78 39 石英 V+L 12×4 7 135.2 -1.9 3.22 40 石英 V+L 10×5 10 149.1 -3.2 5.23 41 石英 V+L 6×4 8 146.2 -3.3 5.37 42 石英 V+L 8×7 8 169.7 -3.7 5.95 注：V—气态；L—液态 ![]()
图 4 加甘滩金矿流体包裹体均一温度直方图(a)及盐度直方图(b)Figure 4. Homogeneous temperature histogram (a) and salinity histogram (b) of fluid inclusions in Jiagantan gold deposit由盐度-均一温度-密度图可知，流体包裹体的密度介于0.6~0.95 g/cm3之间，集中分布在0.80 g/cm3附近(图 5)，表明成矿流体具有中低密度的性质。根据邵洁莲(1988)提出的流体压力经验公式：P1＝P0T1/T0(T0=374+920ω，P0=219+2620ω。式中：T0代表初始温度，T1代表均一温度，P0代表初始压力，ω代表成矿溶液盐度，温度单位为℃，压力单位为105 Pa，计算得出加甘滩矿区压力变化范围为101.6×105~313.4×105 Pa，平均值为189.5×105 Pa，估算深度变化范围为0.4~1.2 km，平均值为0.7 km。
![]()
图 5 加甘滩金矿流体包裹体均一温度-盐度- 密度(Ahmad，1980)Figure 5. Homogeneous temperature-salinity-density diagram of fluid inclusions in Jiagantan gold deposit4. 稳定同位素特征
4.1 氢-氧同位素
氢、氧同位素与包裹体样品为同一样品，分析测试在核工业地质分析测试研究中心完成，所用气体同位素质谱仪型号为MAT-253型，测试结果以V-SMOW为标准，精度±1‰。本次共测试加甘滩金矿氢-氧同位素样品9件。分析结果见表 2。
表 2 加甘滩金矿氢-氧同位素分析结果表Table 2. H-O isotopic analysis results of Jiagantan gold deposit样品原号 岩石/ 矿物 δDV-SMOW/‰ δ18OV-SMOW/‰ δ18OH2O/‰ JZK1408-6 石英 -99.3 20.1 11.32 JZK1408-15 石英 -99.9 21.3 12.52 JZK809-5 石英 -98 19.2 10.42 JZK1805-3 石英 -97.5 20.4 11.62 JZK1609-6 石英 -94.1 21.6 12.82 JZK1609-7 石英 -93.7 21.2 12.42 JZK1609-11 石英 -94.5 21.3 12.52 JGT2-2 石英 -101.2 20.3 11.52 JGT2-3 石英 -93.2 22.6 13.82 氧同位素按照分馏公式(Clayton，1972)：1000lnα石英-H2O=3.34×106/T2-3.31，1000lnα石英－H2O=δ18O石英-δ18OH2O进行换算，温度T取石英-多金属硫化物阶段的平均温度252.5℃。加甘滩金矿矿石中石英的δ18OH2O值在10.42‰~13.82‰之间，平均值为12.11‰；δD值为-101.2‰~-93.2‰，平均值为-96.82‰。在δ18OH2O-δD图解(图 6)中，所有样品点均落在变质水下方区域，靠近原生岩浆水的范围，表明成矿流体组成较复杂，可能既有岩浆水，也有变质水的参与。
![]()
图 6 不同类型水的δ18O-δD图解(据杨学明等，2000)Figure 6. Diagram of δ18O-δD for different types of water4.2 硫同位素
硫同位素测试的样品采自主成矿阶段的金矿石，分别测试黄铁矿、毒砂的硫同位素组成。测试由核工业北京地质研究院完成，仪器为Delta v plus气体同位素质谱计，分析精度优于±2‰，结果以V-CDT为标准。本次研究选择了加甘滩金矿热液期第2成矿阶段2件黄铁矿、1件毒砂样品进行了δ34S值的测定。测试结果见表 3。
表 3 加甘滩金矿硫、铅同位素分析结果Table 3. S and Pb isotopic analysis results of Jiagantan gold deposit样品号 样品名称 S Pb δ34SV-CDT/‰ 208Pb/204Pb 207Pb/204Pb 206Pb/204Pb Δβ Δγ JZK1408-20 黄铁矿 -12.5 38.35 15.59 18.17 17.03 23.38 JZK809-7 黄铁矿 -7.5 38.65 15.68 18.55 22.77 31.44 JZK1609-11 毒砂 -13.4 38.50 15.59 18.41 17.16 27.41 加甘滩金矿床含硫矿物中未见含SO42-等离子硫酸盐，矿区均为硫化物，成矿流体的δ34S值可用硫化物δ34S的平均值代表(Ohmoto et al., 1979)。根据硫同位素分析测试结果，加甘滩金矿黄铁矿δ34S值的变化范围为-12.5‰~-7.5‰，极差为5‰，平均值为-10.1‰；毒砂样品δ34S值为-13.4‰，明显低于黄铁矿。综合来看，δ34S值的范围为-13.4‰~-7.5‰，极差为5.9‰，平均值为-11.13‰，说明加甘滩金矿硫化物δ34S值变化范围较窄，硫同位素组成较稳定。
4.3 铅同位素
铅同位素测试的样品采自加甘滩金矿热液期第二成矿阶段的金矿石，分别测试黄铁矿、毒砂的铅同位素组成。测试仪器为Phoenix热表面电离质谱仪，实验条件为温度20℃，相对湿度20%~30%，精度优于0.005%。测试结果见表 3。
由测试数据可以看出，加甘滩金矿黄铁矿和毒砂铅同位素组成相似，均一程度较高。206Pb/204Pb值为18.14~18.55，平均值为18.38；207Pb/204Pb值为15.59~15.68，平均值为15.62；208Pb/204Pb值为38.35~38.65，平均值为38.50。
为进一步确定加甘滩金矿床矿石的铅来源，将硫化物的铅同位素值进行了铅同位素构造模式和构造环境投图(图 7)。在铅同位素构造模式图207Pb/204Pb-206Pb/204Pb(Zartman et al., 1981)(图 7-a)上，1个毒砂、1个黄铁矿样品点落在上地壳和造山带之间，1个黄铁矿样品点落在上地壳演化线附近。在铅同位素208Pb/204Pb-206Pb/204Pb模式图(Zartman et al., 1981)(图 7-b)上，所有样品的铅同位素值投点落在下地壳演化线和地幔演化线之间，并靠近地幔一侧。在207Pb/204Pb-206Pb/204Pb、208Pb/204Pb-206Pb/204Pb构造环境判别图(Zartman et al., 1981)(图 7-c、d)上，样品点均落在造山带区域，表明加甘滩金矿硫化物中的铅主要来源于造山带铅，具有混合铅的特点，为地幔、地壳相互作用铅同位素均一化的结果(Doe et al., 1979)。
![]()
图 7 铅同位素构造模式图和构造环境判别图(Zartman et al., 1981；A、B、C、D为各区域中样品相对集中区)Figure 7. Pb isotopic structural model map and structural environment discrimination map朱炳泉(1998)认为，Th、Pb的变化及Th、Pb与U-Pb同位素组成的相互关系对于地质过程与物质来源能提供更丰富的信息，为突出这种变化关系，将铅同位素表示成与同时代地幔的相对偏差，直观地表示为Δβ-Δγ成因分类图解。本次研究依据加甘滩金矿样品铅的相对偏差Δβ、Δγ(表 3)，将其放入Δβ-Δγ成因分类图解(图 8)。从图 8可以看出，铅分布较集中，均落入上地壳与地幔混合俯冲带铅范围，与岩浆作用紧密相关。
![]()
图 8 铅同位素△β-△γ成因分类图解(朱炳泉，1998)1—地幔源铅；2—上地壳铅；3a—岩浆作用；3b—沉积作用；4—化学作用铅；5—热水作用铅；6—中深变质铅；7—下地壳铅；8—造山带铅；9—古老页岩上地壳铅；10—退变质铅Figure 8. △β-△γ diagram of genetic classification lead isotopes5. 讨论
5.1 矿床成因
硫同位素组成的差异往往代表来源的差异，硫同位素来源组成分为4类(Ohmoto，1972；1986)：①δ34S值在0附近(-5‰~+5‰)，为地幔、地壳深部物质演化的结果；②δ34S值较大(大于+20‰)，一般认为来源于海相沉积地层或海水；③δ34S值在+5‰~+15‰之间，一般认为为混合源；④δ34S值为较大负值，为开放条件下的有机还原硫。本次研究的加甘滩金矿床中毒砂、黄铁矿δ34S值呈现明显的负值，变化范围较窄，介于上述条件①和④之间。
加甘滩金矿与西秦岭地区主要金矿及全球主要岩石类型硫同位素对比(图 9)揭示，加甘滩金矿及其周围的录斗艘金矿、早子沟金矿的δ34S值变化范围与自然界花岗岩相似，显示成矿流体可能与岩浆活动有关。研究区侵入岩中As、Sb等元素的含量(As=15×10-6~20×10-6，Sb=1.7×10-6~2.1×10-6)与该区二叠系、三叠系沉积岩的含量(As=20×10-6~26×10-6，Sb=1.3×10-6~2.2×10-6)非常接近①，暗示侵入岩与地层岩石发生了强烈的同化混染，相应地发生了元素的交换，一些典型岩浆成因矿床因此出现地层来源硫同位素组成特征(Muntean，2011)。故认为，加甘滩金矿床硫主要来自岩浆或岩浆热液，可能在后期成矿过程中因环境变化引起硫同位素分馏，或有一部分地层硫的加入。
![]()
图 9 加甘滩金矿与西秦岭地区典型金矿床、主要岩石类型硫同位素分布对比图(底图据靳晓野，2013)Figure 9. Comparison of sulfur isotope distribution between Jiagantan gold deposit and typical gold deposits and main rock types of West Qinling area从流体包裹体测试分析结果可以看出，加甘滩金矿成矿流体具有中低温、低盐度的特征，介于陈衍景等(2007)提出的改造热液和变质热液之间，其矿床类型对应卡林型和造山型金矿。而硫同位素组成显示为卡林型，并有岩浆来源特点。加甘滩金矿铅同位素主要为造山带铅，铅的来源与岩浆作用紧密相关，看起来加甘滩似乎又是岩浆热液型金矿。但是，加甘滩金矿流体包裹体主要为气液两相包裹体，并未发现大量富CO2的流体包裹体，不同于典型的造山型金矿和岩浆热液型金矿。铅同位素具有明显的岩浆来源特点，金的赋存形式既有粗粒金，又有显微—次显微金，距加甘滩不远的完安囊金矿石英脉中的明金甚至有大于1 cm者，区别于典型卡林型金矿(微细粒浸染型)。
叶天竺等(2017)研究认为，中低温岩浆热液型金矿床多产于造山运动的中晚期及陆块活化期的中酸性、酸性和偏碱性岩浆活动地区，是由岩浆结晶分异产生的气水溶液形成的金矿床，成矿主阶段温度大致为300±50℃。根据矿床与成矿地质体的空间距离及成矿主阶段的温度，将中低温岩浆热液型金矿床划分为中温热液型金矿床和远成低温热液型金矿床2个亚类。通过与中温热液型金矿床、远成低温热液型金矿床典型特征对比发现，加甘滩金矿主要受张性断裂控制，成矿主阶段包裹体峰值温度在250℃左右，硅化发育，辉锑矿等低温矿物常见于成矿作用主阶段；元素组合为Au-As-Sb-Hg-Ag，As、Sb含量较高，Au与As呈正相关关系，氢-氧、硫、铅等同位素显示岩浆来源特征，并且在加木龙—加甘滩一带遥感解译有环形构造，暗示该区存在热隆构造。近年来，加甘滩金矿东部的拉古河金矿钻探出深部岩体也证实了这一猜想。基于以上分析，加甘滩金矿床成矿特征与远成低温热液型金矿床更吻合，为与深部岩浆结晶分异作用有关的远离成矿地质体的远成低温岩浆热液型金矿床。
5.2 矿床形成机制
综合前文流体包裹体、硫铅稳定同位素特征等研究，加甘滩金矿成矿物质来源与岩浆作用关系密切，成矿物质可能来自深部隐伏岩体。野外观察表明，加甘滩金矿矿体主要受控于断裂构造，矿体与围岩的空间关系主要取决于断裂构造的产状。区域上，早子沟金矿、以地南金矿等产于岩浆活动较强区域的金矿体也均为断裂控制，且控矿断裂穿过岩体(脉)，说明成矿热液作用滞后于岩脉侵位作用。
不少学者(Meng et al., 2005；Zhang et al., 2007；Li et al., 2013)研究认为，华北板块与扬子板块在秦岭地区可能于中三叠世末期发生初始碰撞，中三叠世西秦岭地区的洋壳仍处于俯冲阶段，导致西秦岭印支早期(250~230 Ma)岛弧岩浆岩的侵位，并形成夏河-临潭弧岩浆活动带(李康宁等，2019；2020a)。岩石固结过程中由于收缩应力形成原生节理，并在后期构造运动影响下，形成各类褶曲、断层等变形构造。深部成矿物质伴随着构造-岩浆活动，沿构造有利部位先侵入到一定位置。热泉、下渗水及各种形式的封存水在热能的驱动下对流循环，萃取岩浆中的成矿物质形成含金的地下热液。含矿热液沿断裂破碎带上升，运移到距侵入体较远的部位，在开放的断裂和裂隙系统中充填、交代，同时围岩遭受广泛的酸淋滤，析出大量Fe2+加入溶液中，成矿溶液中的金络合物被还原成独立金，并富集成矿，在地层中形成脉状的远成岩浆热液型金矿床(李康宁等，2020b)。
6. 结论
(1) 加甘滩金矿流体包裹体主要为气液两相包裹体，其中富液相两相包裹体最发育，气液比一般介于10%~20%之间。成矿流体均一温度为132.7~365.9℃，平均值为248.67℃，盐度为0.53%~7.85% NaCl，平均值为3.78% NaCl。成矿流体具有中低温、低盐度特征。
(2) 石英中δ18O H2O值为10.42‰~13.82‰，δD值为-101.2‰~-93.2‰。成矿主阶段黄铁矿和毒砂δ34S值的范围为-13.4‰~-7.5‰，极差为5.9‰，平均值为-11.13‰。黄铁矿和毒砂铅同位素组成相似， 206Pb/204Pb值为18.14~18.55，平均值为18.38；207Pb/204Pb值为15.59~15.68，平均值为15.62；208Pb/204Pb值为38.35~38.65，平均值为38.50。成矿流体组成较复杂，可能既有岩浆水，也有变质水的参与。成矿物质来源可能与深部隐伏岩体和浅变质地层有关。
(3) 加甘滩金矿的形成与印支期洋壳俯冲产生的岩浆作用密切相关，但成矿作用稍滞后于岩浆作用(250~230 Ma)，矿床成因类型应该为与岩浆热液有关的金矿床。
-
![]()
图 1 六块地南区域地质图
1—现代河流冲积层；2—唐家屯组上段；3—唐家屯组下段；4—黑龙宫组；5—花岗斑岩；6—石英正长岩；7—正长花岗岩；8—石英闪长岩；9—角闪石岩；10—石英脉；11—闪长玢岩脉；12—安山岩脉；13—角岩化；14—不整合界线；15—推断性质不明断层；16—地质界线；17—岩相变化界线；18—层理产状；19—岩金矿矿(化)点
Figure 1. Geological map of the Liukuaidi South
![]()
图 3 六块地南土壤地球化学元素聚类分析谱系图
Figure 3. Cluster analysis of soil geochemical elements of the Liukuaidi South
![]()
图 4 测区组合异常图
C2t2—唐家屯组上段；C2t1—唐家屯组下段；D1hl—黑龙宫组；γπJ1—花岗斑岩；ξγJ1—正长花岗岩；δοJ1—石英闪长岩；δμ—闪长玢岩脉
Figure 4. The combination anomaly map of the survey area
![]()
图 8 六块地南查证区TC1内金矿体素描图(部分)
1—腐殖土层；2—残坡积层；3—破碎带；4—流纹质熔结凝灰岩；5—化学样取样位置及编号；6—原岩光谱取样位置及编号；7—标本取样位置及编号；8—薄片取样位置及编号；9—光片取样位置及编号；10—矿体；11—矿化体；12—云英岩化
Figure 8. Sketch of Au orebody in TC1 of the Liukuaidi South inspection areas(part)
![]()
图 9 六块地南查证区TC2内金矿体素描图(部分)
1—腐殖土层；2—残坡积层；3—流纹质熔结凝灰岩；4—化学样取样位置及编号；5—原岩光谱取样位置及编号；6—标本取样位置及编号；7—薄片取样位置及编号；8—矿体；9—矿化体；10—硅化
Figure 9. Sketch of Au orebody in TC2 of the Liukuaidi South inspection areas(part)
表 1 六块地南土壤地球化学元素特征参数
Table 1 Characteristic parameters of soil geochemical elements in the Liukuaidi South
元素 均值 中位数 极大值 极小值 标准差 省均值 变异系数 浓集系数 异常下限 异常个数 Au 1.38 1.07 47.20 0.54 1.76 1.00 1.27 1.38 2.12 11 Ag 0.07 0.07 0.34 0.03 0.03 0.09 0.43 0.79 0.11 18 Cu 14.07 13.37 50.15 6.07 3.86 18.50 0.27 0.76 19.58 5 Pb 27.55 26.67 72.73 14.26 5.78 22.30 0.21 1.24 36.64 10 Zn 63.26 60.54 417.58 32.90 17.33 71.30 0.27 0.89 82.28 11 As 16.67 12.24 377.63 5.15 19.22 9.50 1.15 1.75 27.20 11 Sb 1.48 0.68 159.03 0.30 5.99 0.61 4.05 2.42 2.96 3 Bi 0.25 0.25 0.55 0.10 0.05 0.37 0.20 0.69 0.34 10 W 2.04 1.97 12.89 0.60 0.57 1.85 0.28 1.10 2.60 7 Mo 1.29 1.75 14.95 0.50 0.82 1.48 0.64 0.87 1.85 12 Sn 3.00 2.96 5.40 1.55 0.44 2.40 0.15 1.25 3.80 5 Hg 0.03 0.03 0.14 0.01 0.01 0.05 0.31 0.62 0.04 10 注：Au单位为10-9；其他元素单位为10-6 
下载: 导出CSV
表 2 六块地南土壤地球化学测试数据成矿元素相关系数矩阵
Table 2 Mineralization element correlation coefficient matrix for soil geochemical test data of the Liukuaidi South
元素 Au Ag Cu Pb Zn As Sb Bi W Mo Sn Hg Au 1.00 Ag 0.37 1.00 Cu 0.67 0.51 1.00 Pb 0.31 0.46 0.42 1.00 Zn 0.41 0.47 0.60 0.42 1.00 As 0.08 0.29 0.16 0.22 0.10 1.00 Sb 0.01 0.11 0.07 0.07 -0.03 0.73 1.00 Bi 0.44 0.33 0.73 0.51 0.49 0.04 -0.003 1.00 W 0.78 0.40 0.83 0.36 0.49 0.03 0.07 0.60 1.00 Mo 0.11 0.26 0.04 0.24 0.14 0.51 0.21 0.06 -0.04 1.00 Sn 0.10 0.13 0.36 0.43 0.28 -0.01 -0.02 0.72 0.23 0.03 1.00 Hg 0.59 0.44 0.79 0.44 0.49 0.17 0.21 0.78 0.80 0.06 0.47 1.00
下载: 导出CSV
表 3 因子分析总方差解释
Table 3 Factor analysis of the total variance interpretation table
成分 初始特征值 提取平方和载入 旋转平方和载入 合计 方差的% 累积% 合计 方差的% 累积% 合计 方差的% 累积% 1 5.124 42.698 42.698 5.124 42.698 42.698 3.755 31.291 31.291 2 2.051 17.091 59.789 2.051 17.091 59.789 2.166 18.046 49.337 3 1.254 10.448 70.237 1.254 10.448 70.237 1.814 15.118 64.455 4 1.043 8.689 78.926 1.043 8.689 78.926 1.737 14.471 78.926 5 0.667 5.558 84.484 6 0.546 4.553 89.037 7 0.451 3.758 92.795 8 0.263 2.191 94.986 9 0.236 1.968 96.954 10 0.144 1.198 98.152 11 0.134 1.114 99.266 12 0.088 0.734 100.00
下载: 导出CSV
表 4 旋转前后的因子载荷矩阵
Table 4 Factor load matrix before and after rotation
元素 旋转后的成分矩阵因子 元素 旋转前的成分矩阵因子 成分 成分 1 2 3 4 1 2 3 4 Au 0.854 -0.042 -0.007 0.154 Au 0.714 -0.102 -0.482 -0.043 Ag 0.447 0.046 0.054 0.665 Ag 0.613 0.252 -0.105 -0.445 Cu 0.852 0.321 0.055 0.169 Cu 0.899 -0.107 -0.186 0.073 Pb 0.236 0.511 0.020 0.568 Pb 0.631 0.137 0.376 -0.285 Zn 0.528 0.244 -0.139 0.474 Zn 0.686 -0.052 -0.011 -0.330 As 0.027 -0.013 0.863 0.378 As 0.240 0.900 -0.003 0.144 Sb 0.064 0.010 0.943 -0.035 Sb 0.148 0.766 -0.076 0.528 Bi 0.532 0.769 -0.012 0.102 Bi 0.829 -0.212 0.349 0.175 W 0.934 0.170 0.018 0.018 W 0.836 -0.197 -0.378 0.144 Mo -0.117 -0.017 0.344 0.730 Mo 0.186 0.652 0.154 -0.426 Sn 0.071 0.940 -0.021 0.040 Sn 0.521 -0.186 0.734 0.214 Hg 0.766 0.484 0.190 0.062 Hg 0.888 -0.055 -0.029 0.263
下载: 导出CSV
表 5 六块地土壤异常下限、异常强度分级及单元素异常个数
Table 5 Lower limit of soil anomaly, classification of anomaly intensity and number of single element anomaly in Liukuaidi
元素 异常下限(T) 异常强度分级 单元素异常个数 一级 二级 三级 Au 2.12 2.12~4.24 4.24~8.48 >8.48 11 Ag 0.11 0.11~0.22 0.22~0.44 >0.44 18 Cu 19.58 19.58~39.16 39.16~78.32 >78.32 5 Pb 36.64 36.64~73.28 73.28~146.56 >146.56 10 Zn 82.28 82.28~164.56 164.56~329.12 >329.12 11 As 27.20 27.20~54.40 54.40~108.80 >108.80 11 Sb 2.96 2.96~5.92 5.92~11.84 >11.84 3 Bi 0.34 0.34~0.68 0.68~1.36 >1.36 10 W 2.60 2.60~5.20 5.20~10.40 >10.40 7 Mo 1.85 1.85~3.70 3.70~7.40 >7.40 12 Sn 3.80 3.80~7.60 7.60~15.20 >15.20 5 Hg 0.04 0.04~0.08 0.08~0.16 >0.16 10 注：Au含量单位为10-9, 其余为10-6
下载: 导出CSV
-
Ahrens L H. The lognormal distribution of the elements: a fundamental law of geochemistry and its subsidiary[J]. Geochemica et Cosmochimica Acta, 1954, 5: 49-73. doi: 10.1016/0016-7037(54)90040-X
Zhang G B, Yang Y C, Ye J W, et al. Geology, geochemistry, and genesis of the hot-spring-type Sipingshan gold deposit, eastern Heilongjiang Province, Northeast China[J]. International Geology Review, 2013, 55(4): 482-495. doi: 10.1080/00206814.2012.727572
Zuo R G, Wang J. Fractal/multifractal modeling of geochemical data: A review[J]. Journal of Geochemical Exploration, 2016, 164: 33-41. doi: 10.1016/j.gexplo.2015.04.010
蔡朝阳, 韩秀梅, 吴国学, 等. 统计分析软件SPSS在化探数据处理中的应用[J]. 矿床地质, 2010, 29(S1): 635-636. 曹园园, 李新虎. 地球化学综合异常的圈定及找矿效果[J]. 物探与化探, 2017, 41(1): 58-64. 陈会军, 付俊彧, 钱程, 等. 东北地区前中生代花岗岩类年龄与时空分布[J]. 地质通报, 2021, 40(6): 827-844. http://dzhtb.cgs.cn/cn/article/id/20210602 迟清华, 鄢明才. 应用地球化学元素丰度数据手册[M]. 北京: 地质出版社, 2007: 1-148. 董一博, 焦建刚, 刘凯, 等. 土壤地球化学测量在南秦岭夏家店金矿刘家峡测区的应用[J]. 地质与勘探, 2019, 55(3): 1202-1213. 高艳芳, 李俊英, 陈军威, 等. 地球化学勘查数据迭代处理的可视化及结果分析[J]. 物探与化探, 2016, 40(5): 1021-1025. 高艳芳, 柳青青, 王玮, 等. 正态分布和对数变换在化探数据处理中存在问题的讨论[J]. 物探化探计算技术, 2017, 39(3): 404-410. 蒋敬业, 程建萍, 祈世华, 等. 应用地球化学[M]. 武汉: 中国地质大学出版社, 2006: 1-340. 蒋永健, 魏俊浩, 周京仁, 等. 勘查地球化学新方法在矿产勘查中的应用及其地质效果[J]. 物探与化探, 2010, 34(2): 134-138. 李飞. 博茨瓦纳弗朗西斯敦地区土壤地球化学特征及靶区优选[D]. 成都理工大学硕士学位论文, 2015. 李申, 赵润东, 宋岳庭, 等. 相关性分析、聚类分析、因子分析的可靠性对比——以嘎拉勒和躬琼左波两条剖面为例[J]. 地质与勘探, 2018, 54(3): 574-583. doi: 10.3969/j.issn.0495-5331.2018.03.013 林强, 葛文春, 孙德有, 等. 中国东北地区中生代火山岩的大地构造意义[J]. 地质科学, 1998, (2): 3-13. 刘亚洲, 申维, 李姗. 小兴安岭-张广才岭地区成矿金属元素分布特征及找矿靶区预测[J]. 矿物岩石地球化学通报, 2018, 37(5): 922-931. 马晓阳. 我国东北森林沼泽区化探异常查证方法技术研究[D]. 中国地质大学(北京) 博士学位论文, 2006. 彭中山, 刘君. 西藏某铅多金属矿土壤地球化学特征及找矿[J]. 四川地质学报, 2014, 34(3): 399-403. 邵军, 王世称, 张炯飞, 等. 大兴安岭原始森林覆盖区化探异常查证方法研究与实践[J]. 地质与勘探, 2004, (2): 66-70. 史长义. 研究微量元素区域分散与富集规律的新方法[J]. 物探与化探, 1994, (3): 219-227. 孙景贵, 张勇, 邢树文, 等. 兴蒙造山带东缘内生钼矿床的成因类型、成矿年代及成矿动力学背景[J]. 岩石学报, 2012, 28(4): 1317-1332. 谭红艳. 黑龙江小兴安岭-张广才岭成矿带成矿系列及找矿远景评价[D]. 中国地质大学(北京) 博士学位论文, 2013. 王建, 朱立新, 马生明, 等. 冀北地区龙头山铅多金属矿床的发现及地物化综合找矿模型的建立[J]. 地质力学学报, 2019, 25(1): 9-18. 王磊, 韩润生, 王加昇. 地球化学勘查的新技术及发展趋势[J]. 物探与化探, 2015, 39(4): 686-690. 王秋印, 程华生, 朱梅花, 等. 土壤地球化学测量在内蒙古东乌旗查干楚鲁银铅锌多金属矿区的找矿效果[J]. 物探与化探, 2009, 33(6): 652-656. 王子正, 唐菊兴, 郎新海. 土壤地球化学测量在隐伏矿体勘探中的应用[J]. 西部探矿工程, 2007, (3): 92-95. 吴彦彬. 福建紫金山罗卜岭矿区土壤地球化学特征及找矿预测[D]. 桂林理工大学硕士学位论文, 2021. 肖霞, 倪师军, 白涛, 等. 西藏妥坝地区水系沉积物地球化学特征及找矿靶区[J]. 矿物岩石地球化学通报, 2016, 35(3): 570-575. 杨龙坤, 罗先熔, 文美兰, 等. 地电提取测量法在黑龙江金厂外围区寻找隐伏金矿的应用[J]. 桂林理工大学学报, 2015, 35(4): 809-816. 杨继权, 王秀琴, 刘殿生, 等. 黑龙江省大地构造单元划分及特征[J]. 世界地质, 2007, (4): 426-434, 452. doi: 10.3969/j.issn.1004-5589.2007.04.007 尹志刚, 张志浩, 周小刚, 等. 黑龙江省木兰县双合屯地区晚古生代火山岩时代厘定及其对构造环境的约束[J]. 地质通报, 2022, 41(12): 2146-2157. doi: 10.12097/j.issn.1671-2552.2022.12.007 于俊博, 宋云涛, 郭志娟, 等. R型聚类分析在区域化探元素分组中的作用探讨[J]. 物探化探计算技术, 2014, 36(6): 771-776. 袁和, 罗先熔, 李武毅, 等. 西藏邦卓玛地区土壤地球化学特征及找矿预测[J]. 地质与勘探, 2017, 53(3): 472-481. 臧金生, 王东晓, 赵瑞强. 化探异常定量评价[J]. 物探与化探, 2014, 38(6): 1114-1118. 赵欣怡, 罗先熔, 杨笑笑, 等. 河南洛宁石龙山金多金属矿预查区土壤地球化学特征及找矿远景分析[J]. 矿物岩石地球化学通报, 2020, 39(4): 1-11. 张国宾, 杨言辰, 梁冰, 等. 黑龙江东部完达山地区水系沉积物地球化学特征及成矿预测[J]. 中国地质, 2017, 44(3): 588-603. 中华人民共和国国土资源部. 地球化学普查规范(1: 50000, DZ/T 0011-2015) [S]. 北京: 中国标准出版社, 2015. -
期刊类型引用(7)
1. 焦阳，冯俊环. 西秦岭地区猪婆沟金矿成矿物质来源及矿床成因分析. 西北地质. 2024(01): 219-229 . 
百度学术
2. 袁臻. 多组分气体地球化学测量技术在加甘滩金矿应用研究. 世界有色金属. 2024(09): 145-147 . 百度学术
3. 张陇军. 甘肃西和县冉湾金矿成矿地质条件与找矿潜力研究. 世界有色金属. 2024(10): 106-108 . 百度学术
4. 张陇军. 甘肃两当县太阳寺地区金多金属矿成矿地质条件与资源潜力分析. 世界有色金属. 2024(12): 99-101 . 百度学术
5. LIU Jiawang，WANG Shuzhi，WANG Tao，CAO Chong，QIU Junting，XIE Fenquan. A Comparative Study of Ore-bearing Dikes and a Barren Dike in the Zaozigou Gold Deposit, Western Qinling: Zircon Morphology, U-Pb Chronology and Geochemical Constraints. Acta Geologica Sinica(English Edition). 2024(06): 1489-1500 . 必应学术
6. 樊婷，任子杰，张乾坤，何宇豪，宋昱晗，刘志，向璐. 双极端对比检测法定量分析高纯石英砂气液包裹体含量. 岩矿测试. 2024(06): 892-900 . 百度学术
7. 第鹏飞，汤庆艳，刘东晓，王玉玺，梁自兴，宋宏，刘聪，黎卓明，张燕. 西秦岭甘南地区金矿床黄铁矿微量元素地球化学特征及意义——以加甘滩和早子沟金矿为例. 稀土. 2023(04): 140-154 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 801
- HTML全文浏览量: 407
- PDF下载量: 177
- 被引次数: 7
