甘肃沃尔给楔卡金矿土壤地球化学特征及其评价

甘肃沃尔给楔卡金矿土壤地球化学特征及其评价
杜佰松;申俊峰;秦玉良;徐立为;聂潇;赵玉;牛刚;欧阳尔彪
【摘要】基于甘肃省沃尔给楔卡金矿1∶1万土壤地球化学测量数据,运用R型聚
类分析方法对测量数据进行分析,结果显示:8种元素可分为3类:第一类为Sb、As、Au、Pb、Ag中低温元素组合;第二类为Cu、Zn中高温元素组合;第三类为以Hg
为代表的挥发性元素.根据单元素异常分布特征及元素组合异常,结合地层、岩性、
构造等有利控矿因素综合圈定出3个找矿靶区,并运用系统核理论对其进行异常评价.Ⅰ、Ⅱ和Ⅲ号异常区系统的核分别为{Au,Pb}、{Au}及{Cu},系统稳定性强,其中
Ⅰ和Ⅱ号异常区经钻探和槽探等工程验证,矿化位置、规模与元素异常浓集中心位
置吻合较好,说明在这两个异常区内,土壤地球化学测量所显示的异常为矿致异常;Ⅲ号异常区系统的核为{Cu},并在该异常区内显示以Cu、Au为主要成矿元素的异常,这与野外踏勘中发现有孔雀石及铜蓝较为一致,建议展开对该异常区的相关工程查证.%Based on 1∶10 000 soil geochemical survey on the Woergei gold deposit in Gansu Province, the authors used the R-model factor analysis method to evaluate geochemical anomalies.Cluster analysis of the soil geochemical data obtained from the Woergei gold deposit shows that eight kinds of elements can be grouped into three categories, i.e.,(1)
Sb,As,Au,Pb,Ag;(2) Cu,Zn;(3) bined with favorable ore-controlling factors such as lithology and structure,the authors applied the system core theory to synthetic evaluation of main ore-forming elements and the abnormal element combination,and then selected three prospecting target areas.In No.Ⅰ,Ⅱ and Ⅲ anomaly target zones,the systems of nuclear are respectively {Au,Pb},{Au} and {Cu},the system stability is strong,and the
abnormal Ⅰ and Ⅱ areas were verified by drilling and pit engineering.The mineralization location,size and element anomaly concentrated heart position are in good agreement in the two abnormal zones,suggesting that soil geochemical survey is effective in the study area.However,No.Ⅲ anomaly target zone has the system of nuclear of {Cu},with the anomaly in the area displaying the main ore-forming elements Cu Au anomalies,the malachite and covellite were found in field investigation, which are relatively consistent with the anomaly target area,and hence drilling and pit engineering in the abnormal area should be conducted.
【期刊名称】《物探与化探》
【年(卷),期】2017(041)004
【总页数】7页(P641-647)
【关键词】沃尔给楔卡金矿;土壤地球化学测量;系统核理论
【作者】杜佰松;申俊峰;秦玉良;徐立为;聂潇;赵玉;牛刚;欧阳尔彪
【作者单位】中国地质大学(北京) 地球科学与资源学院,北京 100083;中国地质大学(北京) 地球科学与资源学院,北京 100083;甘肃鼎丰矿业有限责任公司,甘肃兰州730000;中国地质大学(北京) 地球科学与资源学院,北京 100083;中国地质大学(北京) 地球科学与资源学院,北京 100083;中国地质大学(北京) 地球科学与资源学院,北京 100083;甘肃鼎丰矿业有限责任公司,甘肃兰州 730000;甘肃鼎丰矿业有限责任公司,甘肃兰州 730000
【正文语种】中文
【中图分类】P632
位于中国腹地的西秦岭造山带是中国大陆中央造山带的重要组成部分之一，区内印支期花岗质岩体发育，金、铜等多金属矿床星罗棋布，多年来一直受矿业界和学术界所关注，是继胶东金矿之后，被认为可能是我国金矿找矿有希望取得重大突破的地区[1-5]。

研究区位于西秦岭造山带的西北缘，海拔 2 860～4 090 m，区内地
貌起伏相对平缓，属高寒半湿润气候，以化学风化为主，但元素的表生地球化学活性不强[6]，本区除南北局部陡崖地段基岩裸露外，其余区域浮土覆盖较广，植被
发育，较适宜开展土壤地球化学找矿工作。

前人对本区已经开展了地球化学方面的基础工作，并取得了一定的成果，如1988年甘肃省地矿局地球化学探矿队在该区开展了1∶20万化探扫面工作，发现了一系列以金为主的综合分散流异常；2011
年北京西域纵横能源科技有限公司开展了1∶5万水系沉积物地球化学测量工作，圈出2个综合异常区。

笔者基于前人的研究成果，为查证并进一步缩小找矿靶区，对研究区开展了1∶1万的土壤地球化学测量工作，对矿区元素的组合分布特征、单元素异常及其组合异常特征进行了研究，旨在为该区找矿工作提供科学依据。

研究区位于甘肃省合作市境内，大地构造位置处于西秦岭褶皱带北部断褶带与中部裂陷槽之间的过渡部位，属于夏河—礼县逆冲推覆构造带，新堡—力士山复背斜
北翼靠近背斜轴部位置。

区内主体构造格架由NWW向多层次逆冲推覆构造组成，以NW—SE向断裂构造为主[7-9]。

区域上出露地层主要为石炭系、二叠系、三叠系、新近系及第四系。

石炭系地层由下统巴都组(C1b)、中统下加岭组(C2x)、上统上加岭组(C3s)组成，主要岩性为灰岩、砂岩、砾岩及少量的火山岩。

石炭系地层
中丰富的碳酸盐岩为矿区内矽卡岩型矿床和矿化的形成提供了有利的地层条件。

二叠系地层主要出露下统大关山组(P1dg)，其上部岩性主要为深灰色含生物碎屑灰岩,下部为厚层状灰岩,底部为灰岩夹泥质灰岩，是区内主要的含矿层位；二叠系上
统地层为一套滨海—海陆交互相的碎屑岩建造。

三叠系地层主要出露下三叠统山
尕岭组(T1s)和中三叠统古浪堤组(T2g)，主要岩性为板岩、砂岩、砾屑灰岩及含有生物碎屑灰岩等。

新近系主要出露上新统临夏组一套干旱气候的内陆湖盆沉积。

第四系是一些坡积层、腐殖层。

本区构造较发育，断裂总体上以NWW向为主，主要断裂有：位于研究区以南的
力士山—卡加沙格逆断层、答浪沟逆断层、夏河—合作逆断层，经过研究区的西
土房—大槐沟逆断层和分布在研究区以北的铁沟逆断层。

这些断裂是本区的主要
控矿构造和容矿构造，控制着区域总体构造格架。

近年来,随着地质工作的不断深入,本区在Cu、Au找矿方面已取得重大突破，先后发现了德乌鲁、阿姨山、大侠等19处铜矿床(点),枣子沟、加甘滩、直合完干、早仁道、老豆、答浪沟、岗岔等38处金矿床(点)，构成了夏河—合作金铜多金属矿
集区[10-11]。

研究区内岩浆岩较为发育，在区域上与断裂构造共同形成显著的构造—岩浆岩带，除著名的夏河、阿姨山、达尔藏、德乌鲁、吉利、老豆，美武及冶力关等花岗岩类岩基之外，还出露大量小岩株和岩脉，岩性主要为花岗斑岩、花岗闪长岩、黑云母斑状二长花岗岩、花岗岩、石英闪长岩等中酸性侵入岩。

花岗质岩与本区的铜、金成矿时空关系密切，大部分金矿赋存于其岩体之中或附近。

研究区1∶1万土壤地球化学测量面积2 km2，按照100 m×40 m的网度采样，
共取样1 194件，样品采取率为92%以上，达到了1∶1万土壤地球化学测量的要求。

为提高样品的代表性，样品采集时，以采样点为中心，在两侧 10 m范围内按照一点多坑的原则进行采样。

采集深度为地表以下30～50 cm深处的B层或C层中的以粉砂质黏土、细沙质黏土为主细粒级物质。

样品采样过程中应有效地避开有机质、表层腐殖物及其他可能的污染源体，如遇有岩石露头、废石堆、农田、河床堆积等不能取样时弃点。

样品原始质量在500～1 000 g。

用60目不锈钢过筛后
混匀，放入纸袋并密封，以确保样品不受污染。

样品分析了Au、Ag、As、Sb、Cu、Pb、Zn、Hg，测试工作由中国冶金地质总
局一局测试中心完成。

其中Au采用泡沫塑料吸附—石墨炉原子吸收法测定，检出限为0.3×10-9，精密度(Dr,s，n=11)为4.25%；As采用氢化物发生—原子荧光
光谱法测定，检出限为1×10-6，Dr,s为1.83%；Hg采用冷蒸气—原子荧光光谱法测定，检出限为5×10-9，Dr,s为3.20%；Ag采用发射光谱法法测定，检出限为0.03×10-6，Dr,s为5.56%；Sb、Cu、Pb、Zn采用电感藕合等离子体质谱法测定，检出限分别为0.2×10-6、1.5×10-6、5×10-6、15×10-6，Dr,s分别为
1.56%、1.20%、
2.41%、
3.24%。

对研究区土壤各元素分布特征进行研究得知，各元素大致服从对数正态分布规律，因此可以选用传统的数理统计方法计算背景值和异常下限值。

为了避免特高含量对元素背景值及异常下限的影响，笔者采用平均值加2倍标准偏差对数据进行剔除，计算研究区内这8种元素的异常下限及指标参数(表1)。

从表中可知，该区Au、Zn含量与地壳元素丰度值相比基本持平，而Cu 则趋于贫化，低于地壳元素丰度值。

研究区土壤变异系数较大的元素为 Au、Hg，变异系数分别为47.69和38.23，表明这些元素发生了较强的分异作用，为本区有利于成矿的元素。

4.1 相关分析
运用SPSS统计分析软件对原始土壤地球化学数据进行相关分析，得到8种元素的相关性系数 (表2)。

结果显示，主成矿元素Au与伴生元素As、Sb、Pb、Ag呈
显著的正相关，其中Au与As的相关系数为0.76，相关性最好，与Sb、Pb的相关系数分别为0.71、0.57、0.46，说明Au与As、Sb、Pb、Ag可能来源于同一岩体或者同一成矿阶段元素组合；而Au与Zn、Hg几乎不相关，说明研究区Au、As、Sb、Pb、Ag与其他元素可能是分别来自同一岩体的不同部位或不同成矿阶
段的元素组合；Hg表现较为特殊，与研究区其他元素几乎不相关或呈微弱正相关，说明了矿区内可能存在多期次的热液流体活动。

4.2 R型聚类分析
为了直观地反映元素之间的组合规律，对研究区的原始土壤地球化学数据进行R
型聚类分析(图2)[13-15]，结果显示，在相关系数r=0.4的相似水平上，本区的元素可以分为3大类:第一类为 Sb、As、Au、Pb、Ag中低温热液元素组合；第二
类为Cu、Zn中高温热液元素组合，说明该矿区可能受到了另一期次的中高温热
液活动影响，叠加了Cu矿化；还有一类以具有挥发性、低温热液元素Hg为代表。

在第一类中，Au与As、Sb代表矿化元素组合，揭示了矿区金矿与前缘元素之间
的密切关系，剥蚀可能较浅，从而预示着矿体埋藏较深，深部具有良好找矿前景。

综合分析可知，与该矿床关系比较密切的元素为Au、As、Sb、Pb、Ag等，最为密切的是Au、As，因此将Au、As作为本次研究的最佳指示元素，以Au、As、Sb、Pb、Ag土壤地球化学异常组合来圈定找矿靶区，进行隐伏矿体定位预测。

5.1 元素组合特征
在土壤单元素异常分布图(图3)中，Au、As、Sb、Pb、Ag、Zn异常基本沿着黑
云母花岗岩与灰岩的接触带分布，其异常连续性好，分布范围较大，具有明显的异常中心，异常规模较大，吻合性较好，表明异常具有良好的找矿前景。

单元素及元素组合异常总体上呈WN—ES向串珠状展布，主要受矿区内WN—ES向的断裂
构造控制。

值得注意的是，Hg异常主要分布于远离花岗岩体的沉积岩系中，而在花岗岩体上方或其附近几乎不存在异常，这可能是Hg为易挥发元素，在热动力影响下易于逸散，因此Hg扩散较远，在外围显示异常，而在化探核心区几乎未显示异常。

因而可以根据Hg的异常分布特征判断热液流体的运移途径。

在矿区WS—EN方向上Hg异常呈零星分布，并且在西南方向上存在中高温热液元素Cu、Ag
异常区，而在矿区东北部存在大面积的As、Sb、Pb等中低温热液元素异常，推
断矿区从西南向东北成矿温度是逐渐降低的。

5.2 异常评价与验证
通过土壤地球化学测量工作，发现了3个土壤地球化学综合异常块段，由西向东
可分为Ⅰ、Ⅱ和Ⅲ号(图4)，其中Ⅰ和Ⅱ号异常强度相对较高，异常中心叠加明显，规模较大，是找矿最为有利的位置。

同时由综合异也可以看出，Ⅰ和Ⅱ号异常区主要出露于二叠系生物碎屑灰岩与燕山早期花岗质岩的接触带上或其附近，异常由Au、Ag、Cu、Pb、Zn、Sb、As组成。

众所周知，中酸性岩浆与碳酸盐类岩石
相遇时，容易产生有利的接触带构造，如果含矿热液沿着接触带或其附近发生交代作用，就有利于形成矽卡岩矿床[16]。

Ⅰ号异常区与北东向的断裂F3密切相关，
而Ⅱ号综合异常区与北西向的断裂F2密切相关，沿断裂两侧发育不同程度的蚀变，蚀变主要类型为褐铁矿化、绿泥石化、黄铁矿化及云英岩化等。

Ⅰ和Ⅱ号异常区的矿化位置与元素异常浓集中心位置吻合较好，说明在这两个异常区内，Ⅰ和Ⅱ号异常区土壤地球化学测量所显示的异常为有效的异常。

Ⅲ号异常区位于北东向的F3
断裂与北西向的F1断裂的交汇部位周围，相对于Ⅰ和Ⅱ号异常区，Ⅲ号异常区异常强度相对较低，异常主要由Au、Ag、Cu、Sb、As组成。

本次在研究矿区的化探异常时，应用系统核理论对其进行了元素组合异常的评价。

该理论自1994年由许进[17]提出后，周乐尧等[18]、杨利民等[19]把系统核理论
应用于化探异常评价中，并取得了一定效果，直到易桂花等[20]对该理论在化探异常评价中的稳定性进行分析，并提出了改进方法。

笔者在对研究区土壤地球化学原始数据相关性分析的基础上，将相关系数>0.4视为相关，并将相关的各个元素用
线连接起来，构建异常区的系统网络图(表3)。

依据杨利民[19]的理论计算方法，可以很容易地得到Ⅰ、Ⅱ和Ⅲ号异常区系统的核分别为{Au，Pb}、{Au}及{Cu}，系统核度均为1。

对于化探异常体系来说，一般
系统的核就是矿区的主成矿元素，并且系统核度越大，系统越稳定，系统的核越突出[18]。

依据易桂花等[20]针对系统核理论在化探异常评价中的稳定性分析而提出的连通系数、稳定系数，求出Ⅰ、Ⅱ和Ⅲ号异常区系统网络图的连通系数分别为
22/21、11/10、7/7，稳定系数分别为5/4、3/3、5/3。

通常连通系数越大，各
元素之间的连接边数就越多，系统网络图就越稳定。

一般情况下，连通系数>1，
便可认为系统可靠性大[20]。

在Ⅰ、Ⅱ号异常区内，野外踏勘过程中发现发育有大量黄铁矿、黄铜矿、孔雀石、斑铜矿及少量铜蓝等矿物，单样品化验结果显示Au 的含量最高可达30.8×10-6。

为了验证土壤测量方法的有效性和准确性，对研究区土壤地球化学异常区开展槽探和钻探验证。

在Ⅰ号异常区内，通过布设WTC042、WTC082、WTC122、
WTC123、WTC143、ZK100-1、ZK4-1控制工程揭露了Ⅰ、Ⅱ、Ⅲ号金矿体。

其中Ⅰ号矿体倾向171°，倾角55°，地表受WTC123控制，真厚度 2.48 m，金
含量2.62×10-6，深部受ZK100-1钻孔控制，在18.8～32.5 m深度发现3层厚1.1～1.5 m、金含量为2.19～10.47×10-6的矿体。

可以认为Ⅰ号异常区是由金
矿化所引起的矿致异常，是研究区重点找矿靶区。

在Ⅱ号异常区内，由探槽WTC182、EWTC2、EWTC3、EWTC4、EWTC6工程
揭露了Ⅳ号矿体。

矿体的品位较为稳定，一般金含量在(0.7～1.15)×10-6范围内，矿体的厚度不稳定，在0.4～2.02 m之间变化。

探槽样品的分析数据与化探异常
相一致，且异常高值点与金矿点较为吻合，可以认为Ⅱ号异常区也是由金矿化所引起的矿致异常区，这也与稳定性评价表(表3)的结果一致，由此可见利用系统核理
论开展化探异常评价可以取得很好的预期效果。

综上所述，可以认为Ⅱ号异常区为研究区区另一重点找矿靶区。

对于Ⅲ号异常区，根据系统稳定性评价(表3)，系统的核为{Cu}，系统稳定性较强，在该异常区显示
出铜异常，这与野外踏勘中发现有孔雀石及铜蓝较为一致，建议展开对该异常区的相关工程查证。

1) 研究区土壤地球化学数据R型聚类分析将元素分为3大类:第一类为 Sb、As、Au、Pb、Ag中低温热液元素组合；第二类为Cu、Zn中高温热液元素组合；还
有一类以具有挥发性、低温热液元素Hg为代表；说明了矿区内可能存在多期次的热液流体活动。

2) 在矿区的西南方向上存在中高温热液元素Cu、Ag异常区，而在矿区东北部存在大面积的As、Sb、Pb等中低温热液元素异常，推断矿区从西南向东北成矿温
度是逐渐降低的。

3) 矿区内Ⅰ、Ⅱ号异常区是以Au为主要成矿元素，As、Sb、Pb、Ag作为指示
元素的矿致异常，而Ⅲ号异常区系统的核为{Cu}，在该异常区内显示以Cu、Au
为主要成矿元素的异常，这与野外踏勘中发现有孔雀石及铜蓝较为一致，建议展开对该异常区的相关工程查证。

致谢：野外取样工作得到了沃尔给矿区有关工作人员的大力支持，在论文撰写过程中，得到了徐立为和刘圣强师兄的精心指导，同时也得到了匿名审稿老师对本文提出的宝贵建议和认真修改，在此一并表示衷心感谢！
【相关文献】
[1] 张国伟，张本仁，袁学诚，等.秦岭造山带与大陆动力学[M].北京:科学出版社，2001.
[2] 殷勇，殷先明.西秦岭北缘与埃达克岩和喜马拉雅型花岗岩有关的斑岩型铜—钼—金成矿作用[J].岩石学报，2009，25(5)：1239-1252.
[3] 鲍霖，申俊峰，曹卫东，等.甘肃岗岔金矿3号矿脉原生晕分带特征及其找矿意义[J].物探与化探，2014，38(4):660-666.
[4] 张琪，殷先明，殷勇，等.西秦岭与埃达克岩和喜马拉雅型花岗岩有关的金铜成矿及找矿问题[J].岩石学报，2009，25(12):3103-3122.
[5] 李金春，申俊峰，刘海明，等.西秦岭岗岔金矿赋矿围岩成岩时代及其地质意义[J].现代地质，2016，30(1):36-49.
[6] 徐家乐，冯治汉，牛洪斌，等.甘肃省区域化探方法技术和主要成果[J].甘肃地质学报，2002，11(2):67-73.
[7] 靳晓野，李建威，隋吉祥，等．西秦岭夏河—合作地区德乌鲁杂岩体的侵位时代、岩石成因及
构造意义[J].地球科学与环境学报，2013，35(3):20-38．
[8] 张国伟，张本仁，袁学诚，等.秦岭造山带与大陆动力学[M].北京:科学出版社，2001,1-855.
[9] 冯益民，曹宣铎，张二朋，等.西秦岭造山带结构造山过程及动力学[M].西安:西安地图出版社，2002.
[10]骆必继，张宏飞，肖尊奇.西秦岭印支早期美武岩体的岩石成因及其构造意义[J].地学前缘，2012，19 (3):199-213．
[11]金维浚，张旗，何登发，等．西秦岭埃达克岩的SHRIMP定年及其构造意义[J].岩石学报，2005，21(3):959-966．
[12]黎彤.化学元素的地球丰度[J].地球化学，1976，5(3):25-30．
[13]杨小峰，刘长垠，张泰然，等.地球化学找矿方法[M].北京:地质出版社，1997.
[14]于俊博，宋云涛，郭志娟，等.R型聚类分析在区域化探元素分组中的作用探讨[J].物探化探计算技术，2014，36(6):771-776.
[15]Hugh R R，杨学明，杨晓勇，等.岩石地球化学[M].合肥:中国科学技术大学出版社，2000.
[16]翟裕生，姚书振，蔡克勤.矿床学(第三版)[M].北京：地质出版社，2011.
[17]许进.系统核与核度理论及其应用[M].西安：西安电子科技大学出版社，1994.
[18]周乐尧，邱郁双.一种新的化探异常评价方法[J].地质与勘探，1998，34(6):40-43.
[19]杨利民，杨自安，罗铁良，等.青海两兰地区化探异常筛选评价的新方法研究[J].地质与勘探，2007，43 (5):77-81.
[20]易桂花，倪师军，张廷斌，等.系统核理论在化探异常评价中的稳定性分析及方法改进[J].成都理工大学学报：自然科学版，2015，42(1):115-121.。