蒙古国那仁陶勒盖金矿床原生晕特征及深部预测

吴涛涛，金 鑫，王庆双，姚 远，鲍庆中，周永恒，柴 璐

中国地质调查局沈阳地质调查中心，辽宁 沈阳 110034

0 引言

热液矿床原生晕具有明显的组分分带特征，利用原生晕分带研究能够对隐伏矿体进行预测，这种预测方法在危机矿山和小型矿山深部矿体预测方面已得到广泛应用，并取得了较好的找矿效果［1-7］.

那仁陶勒盖金矿床位于蒙古国中央省境内.该矿自20 世纪60 年代发现以来，累计探明黄金储量超过8 t，达到中型规模.该矿床的研究程度十分薄弱，苏联时期民主德国地质队在勘探过程中对该矿床的地质特征进行了较为系统的总结，近年来也有部分学者通过野外地质调查，对矿床地质特征及成因进行了探讨［8］.然而到目前为止，尚未有学者对那仁陶勒盖金矿床的元素地球化学特征及原生晕分带进行研究.随着近年来对该矿床的不断开发，矿区探明的资源量已基本开采完毕，查明矿床深部状况和寻找新的隐伏矿体成为矿山工作的重点.为此，笔者选择那仁陶勒盖金矿床最具代表性的①号主矿体作为研究对象，对其原生晕特征进行研究，计算原生晕轴向分带规律，建立矿体叠加晕理想模型，并在此基础上开展深部找矿预测，以期指导矿山勘查工作，同时弥补该矿床在原生晕分带研究方面的空白.

1 矿床地质概况

那仁陶勒盖金矿床位于蒙古国最重要的北肯特金成矿带［9-11］，区内分布有博洛、盖特苏尔特、苏尔特等大型金矿床和大量中小型金矿床和矿点.该区大地构造位置处在蒙古-鄂霍次克缝合带内，蒙古外贝加尔褶皱系之哈拉盆地地体中［8，12］（图1a）.

那仁陶勒盖金矿矿区地层出露简单（图1b），主要为下寒武统—上奥陶统哈拉群上段希尔古组（C3—O1h2），岩性主要为变质砂岩、类复理石变质杂砂岩，局部为粉砂质片岩、石英岩、泥质片岩等.矿区岩浆活动强烈，主要出露晚三叠世—早侏罗世哈扎尔花岗杂岩，岩性主要为花岗闪长岩，岩体侵入到希尔古组粉砂质片岩中，并在接触带附近发生强烈的接触交代变质作用.此外，矿区还零星出露闪长岩、石英斑岩等脉岩，脉岩走向大多以北东向和南东向为主、矿区构造较为简单，主要有北东向、北北东向、北东东向及北西向断裂构造，以北东向断裂构造及与其平行的断裂构造为主，该断裂构造控制矿区①号主矿体的走向.

图1 那仁陶勒盖金矿大地构造位置及矿区地质图（据文献［8］）Fig.1 Tectonic location and geological sketch map of Naran Tolgoi gold deposit（From Reference［8］）

矿区内共发现6 条含金矿体，走向均呈北东向（图1b）.其中①号矿体为矿区规模最大的一条矿体，为石英脉型，赋存于花岗闪长岩体和变质砂岩的接触带附近.矿体在地表断续出露，长度大于1 500 m，厚0.2～3.0 m，总体产状305°∠65°，倾角具有北陡南缓的特点；矿体沿走向和倾向均具舒缓波状，呈脉状、似脉状、透镜状、囊状产出，并具有向南西侧伏的趋势；矿体平均厚度0.95 m，平均品位8.92×10-6，局部品位和厚度变化较大.其余矿体均位于①号矿体两侧，均为石英脉型，走向与①号矿体近于平行，但成矿规模均较小.

通过野外地质调查，并结合前人研究成果［8］，将该矿床划分为3 个成矿阶段：石英-毒砂-黄铁矿阶段（Ⅰ），主要形成石英脉主体，在其内分布少量毒砂、黄铁矿等金属矿物；石英-多金属硫化物阶段（Ⅱ），主要形成含金的石英多金属矿物，包括硫砷铜矿、闪锌矿、砷黝铜矿、黄铜矿、方铅矿、自然金等脉石矿物和矿石矿物；碳酸盐阶段（Ⅲ），主要形成石英、绢云母、碳酸盐矿物.

2 样品采集与分析

选择那仁陶勒盖矿区11 和33 号勘探线控制的①号矿体1064、1010、964 和910 m 4 个中段及ZK11-1、ZK11-3 等钻孔的39 件矿石及近矿围岩进行系统采样，坑道和钻孔中采样以矿体为中心外延3～5 m 捡块取样（其中11 号勘探线采样位置见图2 所示）.样品经整理后送SGS-IMME Mongolia LLC 分析测试，分析的 主 要 元 素 为Cu、Zn、Pb、W、Mo、Ag、Sn、Bi、Sb、Hg、As、Au、Co、Ni、Nb、Ta、Li、Cs、Mn 共19 种.元素分析参照中国采矿行业标准DZ/T0130：2006 执行，其中As、Hg 元素分析先将样品溶于王水中，然后采用原子荧光光谱法分析；Cu、Zn、Pb、W、Mo、Ag、Sn、Bi、Sb、Co、Ni、Nb、Ta、Li、Cs、Mn 元素分析先采用4 酸消解法（HF ＋HNO3＋HClO4＋HCl）将样品进行溶解，然后使用ICP-MS 分析；Au 元素分析先将样品溶于王水中，然后使用ICP-MS 分析.

图2 那仁陶勒盖金矿11 号勘探线剖面图及原生晕采样位置示意图Fig.2 Profile of Naran Tolgoi gold deposit along No.11exploratory line with sampling location of primary halo

3 矿床原生地球化学异常特征

3.1 矿区围岩微量元素含量特征

那仁陶勒盖金矿床赋矿围岩主要为变质沉积岩（粉砂质片岩、板岩等）和花岗闪长岩，其微量元素含量特征如表1 所示.其特征如下：所有样品中浓集克拉克值大于1 的元素从大到小依次为Bi（47.29）→Au（9.67）→W（4.83）→Ag（4.18）→Cs（3.54）→Pb（2.12）→As（1.72）→Mo（1.38）→Li（1.14）→Sn（1.11）；两套围岩的Au、Ag、Bi 等元素浓集克拉克值均大于1，但在花岗闪长岩中Au、Ag 等主成矿元素浓集克拉克值明显高于变质沉积岩中的值.

表1 那仁陶勒盖金矿围岩中指示元素含量特征Table 1 Contents of indicator elements in the wall rock of Naran Tolgoi gold deposit

3.2 原生晕元素组合特征

为研究那仁陶勒盖金矿床中微量元素的相关性，笔者将39 件样品分为3 类：1）含矿石英脉样品；2）近矿（蚀变）花岗闪长岩样品；3）近矿（蚀变）砂岩、板岩等地层样品.在此基础上分别对这3 类样品进行聚类分析，聚类分析谱系图见图3.

图3 那仁陶勒盖金矿样品聚类分析谱系图Fig.3 Cluster analysis pedigree of Naran Tolgoi gold deposit samples

由图3 可以看出，在矿体中，与Au 关系密切的元素主要为Sb、Nb、Ta、As、W、Cs、Li 等；在花岗闪长岩中，与Au 关系较为密切的元素主要为As、Sb、W、Cs、Mo；在地层中，Au 显示出相对独立的地球化学行为，仅与Sb、Cu、Ni、Co 等元素有一定的相关性.将3 类样品元素的谱系图进一步对比可以看出，矿体和花岗闪长岩中的Au 与其他微量元素相关性有一定的相似性，而与地层中的差别较大，反映了该矿床的金成矿与岩体关系更密切的特点.

3.3 原生晕元素空间分布特征

那仁陶勒盖金矿床矿体组成较为简单，①号矿体为矿区规模最大的一条矿体，矿区11 号勘探线为矿山主要开采区域，钻孔控制矿体深度最深，且该矿体受构造控制明显，矿体与围岩界限清晰，为开展原生晕分带研究提供了良好的条件.为此，选择①号矿体11 号勘探线附近不同中段的5 件样品（表2）作为研究对象，对其垂向变化规律及轴向分带特征进行研究.

3.3.1 原生晕地球化学参数垂向变化规律

从①号矿体11 号勘探线中与Au 成矿密切相关的Au、As、Sb、Hg、Cu、Pb、Zn、Ag、W、Sn、Bi、Mo 等12种微量元素在不同中段的含量变化图（图4）中可以看出，从标高1064 m 到873 m，随着Au 含量的逐渐降低，Sb、Hg 在1010 m 中段出现峰值后呈明显减小的趋势，而As 则呈先增加后减小的趋势；Cu、Pb、Zn、Ag 在1010 m 中段出现峰值后呈迅速下降的趋势；W、Bi 在1010 m 中段出现峰值先下降后逐渐增加，Sn、Mo 与Cu、Pb、Zn、Ag 具有相似的变化规律.

图4 那仁陶勒盖金矿①号矿体11 号勘探线微量元素轴向变化曲线图Fig.4 Axial variation curves of trace elements in No.1 orebody of Naran Tolgoi gold deposit along No.11 exploratory line

3.3.2 原生晕轴向分带特征

元素轴向分带计算的方法众多［14-15］，不同学者依据不同原理对元素轴向分带规律进行了研究.本研究采用苏联学者C.B.格里戈良等人提出的分带指数法对研究区元素轴向分带进行计算［14］.其计算步骤为：首先整理各中段各元素含量及样品长度数据（表2），并计算线金属量，然后将线金属量标准化，使各元素的线金属量的最大值统一到同一个数量级，再计算分带指数，初步排出元素分带序列为（Au、Co、Ni、Nb）-（Cu、Zn、Pb、Mo、Ag、Bi、Sb、Hg）-（Sn、Ta、Mn）-As（W、Li、Cs），最后再计算同一个序列中各元素的变异性指数，最终得出分带序列（由浅至深）：Au-Nb-Co-Ni-Bi-Pb-Mo-Sb-Cu-Hg-Ag-Zn-Ta-Sn-As-W-Li-Cs.

表2 那仁陶勒盖金矿不同中段微量元素含量Table 2 Contents of trace elements in Naran Tolgoi gold deposit at different elevations

3.4 原生晕地球化学参数特征

根据原生晕轴向分带序列，选择前缘晕、近矿晕和尾晕特征元素以及由此建立的累加、累乘以及比值等，可作为参数反推金矿床轴（垂）向地球化学参数叠加结构的理想模型，从而评价矿体的深部含矿性，并进行矿体深部预测［1，16-18］.通常（Pb·Zn）/（Bi·Mo）和（Pb·Zn）/（Au·Ag）值可以作为判别矿体发育程度的指标，矿体附近比值通常会出现较大的波动和峰值；而（As＋Sb）/（Bi＋Mo）、100Sb/（Bi·Mo）从矿体头部到尾部一般呈现逐渐减小的规律，高比值区通常为矿体的上部［17］.为此，选择以上4 个指标作为分带性指数来指示矿体的剥蚀程度，进而对矿体深部远景进行地球化学评价（图5）.

由图5 可以看出，主矿体在地表至海拔873 m 附近，各曲线升降特点与李惠等总结的金矿床轴向地球化学参数叠加结构理想模型的同位叠加模型相似［1］.以海拔873 m 为转折点，近矿指标（Pb·Zn）/（Bi·Mo）和（Pb·Zn）/（Au·Ag） 随着海拔高度的降低而降低，而（As＋Sb）/（Bi＋Mo）指标由海拔873 m 的4.09 上升至710 m的15.23，100Sb/（Bi·Mo）由海拔873 m 的18.60 上升至710 m 的112.33，指标值均发生了反转，与李惠等总结的理想模型B 相似［1］，为深部盲矿前缘晕叠加所致.

4 深部矿体预测

前人通过对中国58 处金矿床原生晕轴向分带序列的概率统计，得出了中国金矿床原生晕综合轴向（垂直）分带序列从上到下为：B-I-As-Hg-F-Sb-Ba→Pb-Ag-Au-Zn-Cu→W-Bi-Mo-Mn-Ni-Cd-Co-V-Ti［1，16］.将那仁陶勒盖金矿①号矿体的元素分带特征与理论轴向分带序列对比发现，该序列显示出至少存在2 期叠加金成矿成晕的特点（图5）：一期为现在正在开采的金矿体，具体表现在Au 位于序列的头部，且Co、Bi、Ni等尾部晕位于序列的中部，说明正在开采的①号金矿体近地表位置即为矿体的中部，这与矿区地表即出露厚大矿体的事实相吻合，矿体至海拔873 m 矿体可能开始尖灭，这也可能是ZK11-2 钻孔在340 m 附近（海拔约810 m 处）的断层破碎带附近虽然广泛发育硅化、绢云母化、黄铁矿化等矿化蚀变现象，但金品位不高的原因.另一期金矿化处在①号金矿体的正下方，具体表现在As 显示反向分带的特征，出现在序列的尾部，且与W、Sn 等尾晕元素共同构成本次元素分带序列的尾部晕，显示出深部具有盲矿体或深部还存在第二个Au 富集地段的可能性，且该盲矿体的头部晕与①号矿体的尾部晕产生了叠加.

从分带指数法计算所得的轴向分带序列图（图5）可看出，主矿体在地表至海拔873 m 附近，各曲线升降特点与前人总结规律十分相似，较好的代表了矿体头晕、近矿晕和尾晕的特点.以海拔873 m 为转折点，（As＋Sb）/（Bi＋Mo）、100Sb/（Bi·Mo）这2 项指标发生反转，说明在710 m 以下可能构成1 个矿体前缘晕的特点，指示深部可能存在盲矿体.

图5 那仁陶勒盖金矿①号矿体叠加晕预测模型Fig.5 Prediction model for superposed halo of No.1 orebody in Naran Tolgoi gold deposit

在成矿地质条件方面，前人对同属北肯特金成矿带的苏尔特金矿床（图1a 中9 号矿床）测得蚀变矿物Ar-Ar 年龄为178 Ma，形成于早侏罗世，说明该地区金矿化时期主要为早侏罗世［19-20］，与那仁陶勒盖金矿床的赋矿岩体花岗闪长岩形成时代相近.本研究对矿区4 件远离矿体的花岗闪长岩样品测定Au 含量平均值为55.70×10-9，明显高于6 件地层样品中Au 含量平均值26.53×10-9（表1），远高于地壳平均丰度4.1×10-9［13］；花岗闪长岩中与金伴生的Ag 平均含量为0.56×10-6，也明显高于地层中Ag 含量平均值0.11×10-6（表1）和地壳平均丰度0.05×10-6［13］.在原生晕元素组合特征方面，矿体和花岗闪长岩中的Au 与其他微量元素相关性有一定的相似性，而与地层中的差别较大，反映了该矿床的金成矿与岩体关系更密切的特点.从以上分析可以看出，矿区晚三叠—早侏罗世花岗闪长岩不仅作为金矿体的赋矿围岩，它还可能为金成矿提供了物质来源.

此外，参考矿区施工的瞬变电磁剖面工作成果，在11 号勘探线西南方向900 m 处6 号勘探线和东北方向600 m 处23 号等勘探线PROTEM 反演模型断面推测矿体下部有断裂构造带向下延伸至海拔高度600 m以下，这与ZK11-1 和ZK11-3 钻孔在海拔约711 m处发现的断裂构造带相吻合.

基于上述研究，笔者认为在①号矿体下部还具有成矿潜力，推测盲矿体赋存位置在海拔710 m 以下的岩体与地层接触带附近.

5 结论

1）那仁陶勒盖金矿床中Au 与Sb、Nb、Ta、As、W、Cs、Li 等元素关系密切，与赋矿花岗闪长岩相似，结合区域上金矿床成矿时代多集中在早侏罗世以及岩体含Au 背景值较高的特点，认为该矿床金矿化与晚三叠—早侏罗世花岗闪长岩有关.

2）运用格里戈良分带指数法确定那仁陶勒盖金矿床元素轴向分带序列为：Au-Nb-Co-Ni-Bi-Pb-Mo-Sb-Cu-Hg-Ag-Zn-Ta-Sn-As-W-Li-Cs.与正常的原生晕序列对比，其前缘晕元素As 与尾晕元素W、Sn 在尾部叠加，显示出深部可能有盲矿体存在，且该盲矿体的头部晕与①号矿体的尾部晕产生了叠加.

3）综合元素轴向分带序列特征及原生晕地球化学参数特征，建立了那仁陶勒盖金矿①号矿体叠加晕理想模型，预测①号矿体深部存在盲矿体，为矿区深部矿体评价预测提供了依据.