胶北隆起中生代壳幔岩浆的混合反应是巨量金质来源的关键 *

田瑞聪 李大鹏 张文 田京祥 于晓卫 耿科 张岩

1. 中国地质大学(北京)地球科学与资源学院，北京 100083 2. 自然资源部金矿成矿过程与资源利用重点实验室，山东省金属矿产成矿地质过程与资源利用重点实验室，山东省地质科学研究院，济南 250013 3. 山东省地质调查院，济南 250014

胶东矿集区是世界著名的金矿资源基地，累计探明金矿资源量已达5480余吨，为继南非兰德金矿区和乌兹别克斯坦穆龙套金矿区后的世界第三大金矿区(陈玉民等, 2019)，深部找矿潜力巨大。与找矿和成矿作用有关的研究工作一直是热点，很多重大问题也得到了很好地解决，为深部找矿提供了很好的技术支持(邓军等, 2001; 庞绪成, 2005; 李建威等, 2010; 杨立强等, 2014; 宋明春等, 2015)，有些问题仍未解决或者存在争议，如金质来源和金元素丰度问题等。

许多学者开展了大量有关成矿物质来源与富集机制的专题研究(刘辅臣等, 1984; 刘连登等, 1984; 杨士望, 1986; 王鹤年, 1988; 孙丰月等, 1995; 杨忠芳, 1998; 孙景贵等, 1999; 邓军等, 1996, 1999; Yangetal., 2014; 陈玉民等, 2019; 田杰鹏, 2020; Tianetal., 2020; Wangetal., 2021)，取得了丰硕成果，主要有3种观点：①与早前寒武纪地质体有关，来源于地壳(杨士望, 1986; 吕古贤, 1998; 王义文等, 2002; 郭春影, 2009)；②来源于中生代花岗岩(Boyle, 1979; 王鹤年, 1988; Wangetal., 1998)。其中多数学者认为来源于玲珑花岗岩和郭家岭花岗岩(桑隆康等, 1992; Wangetal., 1998; 孙华山等, 2007; 于学峰等, 2012; 罗贤冬等, 2014; 田杰鹏等, 2016)，仅有少数学者认为来源与伟德山花岗岩也有密切关系(宋明春等, 2015; Songetal., 2014)；③来源于深部地幔(孙丰月等, 1995; 张连昌等, 2002; 毛景文等, 2002, 2005; 牛树银等, 2011; 丁正江, 2014; Goldfarb and Santosh, 2014; Li and Santosh, 2014; 朱日祥等, 2015; Wangetal., 2020; Dengetal., 2020a, b; Grovesetal., 2020)。有学者更强调了岩脉对金矿成矿的作用，如地幔来源的煌斑岩与金矿化有成因联系(刘辅臣等, 1984; 孙景贵, 1999; 罗镇宽等, 2001; 刘燊, 2004; 龚庆杰等, 2012; 祝德成等, 2018)或为金矿成矿提供了部分成矿物质；深源基性脉岩活动及与其相伴生的大规模地幔流体参与是胶东金矿形成的主要原因(王登红等, 1999; 毛景文等, 2005)。

现在越来越多研究者更倾向于来源于中生代花岗岩或地幔，其主要证据主要得益于同位素示踪技术的大量应用，元素地球化学方面很难拿出有说服力的证据，虽然在Au等成矿元素丰度方面做了大量的研究，但由于不同年代实验室、仪器、测试精度不同，以及单次数据量少和系统误差等原因，很难数据集成形成共性认识。有学者认为Au丰度较高，如黎彤(1992)认为玲珑花岗岩的地壳丰度为3.4×10-9，王鹤年等(1988)认为胶东群变质岩的原岩为玄武岩，岩石建造中Au丰度为地壳丰度的近5倍。然而，孙景贵(1999)认为，胶东变质岩和地层中金的丰度并非很高1.9×10-9～4.5×10-9，而与金矿有关的花岗岩金的量更低0.27×10-9～1.83×10-9。Yanetal.(1997)研究了中国东部上地壳化学组成，其中Au元素丰度为0.77×10-9，华北板块上地壳为0.74×10-9。21世纪随着分析仪器精度的提高，对Au元素丰度的研究也逐步深入，迟清华(2002)统计了中国东部地区酸性岩平均为0.54×10-9，片麻岩类为0.66×10-9，沉积岩类变化较大，同时给出了地壳、出露岩石地壳、火成岩、沉积岩和变质岩中的Au丰度，分别为1.0×10-9、0.80×10-9、0.90×10-9、0.60×10-9、1.0×10-9。史长义等(2005)给出了中国花岗岩Au元素丰度0.38×10-9。李瑞红等(2019)提出胶东地区和中国东部变质岩Au丰度相当，莱阳群和青山群沉积岩中Au丰度普遍较低，侵入岩Au丰度为1.37×10-9。陈玉民等(2019)认为胶东西北部玲珑花岗岩Au丰度为0.73×10-9，郭家岭花岗岩Au丰度为0.69×10-9，均低于华北地台大陆地壳丰度值(1.0×10-9)，胶东金矿集区存在金的地球化学亏损场。Wangetal.(2021)用高灵敏度Element XR同位素稀释法和/或聚氨酯泡沫预处理后的石墨炉原子吸收光谱法(GFAAS)，实现了低Au含量样品的可靠分析，得出胶东前寒武纪基底变质岩中金含量强烈亏损(Au<0.5×10-9)，很难成为胶东金矿床的主要金源。

上述Au丰度数据相差较大，对胶东地区Au丰度没有形成统一的认识，地质体不同的Au丰度说明了什么问题？与金质来源有无联系？基于这些问题，本文利用5529件Au等15种元素分析数据，对其进行统计分析，从不同地质体中的元素丰度变化、区域地球化学背景、元素变异系数、标准离差、正态分布等方面进行对照分析研究，试图寻找出胶东地区元素地球化学背景和可能的金矿成矿地球化学信息。根据不同地质单元特别是中生代不同期次侵入岩的元素丰度比较分析，探求与金矿成矿关系密切的地质体和可能的巨量金质来源或聚集迁移机制，为胶东地区巨量金质来源的研究提供一种新思路。

1 区域成矿地质背景

研究区位于中国东部，大地构造位置属华北板块之胶北隆起区和秦岭-大别-苏鲁造山带之胶南-威海隆起区(图1)。新太古代、古元古代基底变质岩系和中生代多期次多成因岩浆活动以及NE向断裂构造为主的构造格架，构成本区金矿的基本成矿地质背景(邓军等, 2000; 于学峰等, 2012; 宋明春等, 2013; 杨立强等, 2014; 田杰鹏等, 2016)。

图1 胶东地区地质矿产简图(据田杰鹏等, 2016修改)1-第四系；2-晚白垩世-古近纪王氏群砂岩；3-早白垩世青山群火山沉积岩；4-早白垩世莱阳群砂岩、粉砂岩；5-新元古代蓬莱群浅变质岩；6-古元古代芝罘群变质岩；7-古元古代粉子山群黑云变粒岩、大理岩等；8-古元古代荆山群黑云片岩、变粒岩等；9-早白垩世崂山花岗岩；10-早白垩世伟德山花岗岩；11-早白垩世郭家岭花岗岩；12-晚侏罗世玲珑花岗岩；13-中侏罗世文登花岗岩；14-晚三叠世花岗岩；15-南华纪花岗质片麻岩；16-中元古代变辉长岩；17-新太古代TTG；18-中生代脉岩；19-金矿床Fig.1 Geology and deposit sketch map of Jiaodong area (modified from Tian et al., 2016)1-Quaternary; 2-Late Cretaceous-Paleogene Wangshi Group sandstone; 3-Early Cretaceous Qingshan Group volcanic sedimentary rock; 4-Early Cretaceous Laiyang Group sandstone and siltstone; 5-Neoproterozoic Penglai Group epirock; 6-Paleoproterozoic Zhifu Group metamorphic rock; 7-Paleoproterozoic Fenzishan Group biotite granulite and marble; 8-Paleoproterozoic Jingshan Group biotite schist and granulite; 9-Early Cretaceous Laoshan granites; 10-Early Cretaceous Weideshan granites; 11-Early Cretaceous Guojialing granites; 12-Late Jurassic Linglong granites; 13-Middle Jurassic Wendeng granites; 14-Late Triassic granites; 15-Nanhuan granitic gneiss; 16-Mesoproterozoic metagabbro; 17-Neoarchean TTG; 18-Mesozoic dike rocks; 19-gold deposits

区内地层主要发育新太古代胶东岩群，古元古代荆山群、粉子山群，白垩纪莱阳群、青山群、王氏群。岩浆活动强烈，主要有新太古代马连庄序列变质基性-超基性岩、栖霞序列变质英云闪长岩、谭格庄序列变质奥长花岗岩、古元古代莱州序列变质基性-超基性岩、晚侏罗世玲珑花岗岩、早白垩世早期郭家岭花岗岩、早白垩世晚期伟德山花岗岩。侵入岩总体呈近东西或北东向展布的岩基、岩株、岩瘤状产出，具规模性的侵入体多聚集形成复式岩体，岩石类型齐全，中酸性岩规模大、分布广(张良等, 2014)。区内断裂构造发育，其中主要构造线是北北东向和北东向断裂构造，与金成矿休戚相关，控制了该区主要金矿床，北西和东西向断裂发育弱(宋明春等, 2015; 杨立强等, 2019)。以往地质勘查表明，胶东地区金矿床主要受三山岛断裂、焦家断裂、招远-平度断裂、台前-陡崖断裂、牟平-即墨断裂和金牛山断裂等6条NE向区域断裂构造带及其次级构造控制，共同构成了胶东金矿控矿断裂构造系统，控制着该区的主要金矿床(Dengetal., 2006; 张良等, 2013; 杨立强等, 2014; 图1)。

2 样品采集、分析方法与结果

2.1 样品采集

研究区共采集各类岩石样品5529件，其中未蚀变样品3888件，蚀变样品1641件，采样点分布于胶北隆起区，涉及胶东地区主要地质体。取样时，在采样点周围均匀敲取3～5块同一岩性岩石碎块组成一件样品。在地质体中采集无矿化且较新鲜岩石，在构造带、蚀变带、矿化带则尽量采集矿化较好的岩石。详细记录岩石新鲜程度和是否蚀变等，进行严格区分。

2.2 分析方法及精度

分析测试由山东省第四地质矿产勘查院实验测试中心承担。分析测试元素为Au、Ag、Hg、Cu、Zn、Sn、Sb、Pb、Bi、As、Nb、Ta、Mo、W、Co计15种元素。工作标准为《地质矿产实验室测试质量管理规范》(DZ/T 0130—2006)、《地球化学普查规范》(DZ/T 0011—2015)。样品在60℃以下烘干后再进行加工。样品加工过程中先用配备高铝瓷的颚式破碎机粗碎，再采用高铝瓷的无污染化探样品自动粉碎机加工，样品经粗碎、中碎、细碎至0.074mm，符合粒度要求的试样质量不少于加工前试样质量的90%，不缩分，不过筛，直接装入样品袋待检测(重量>300g)。样品制备后，测试样由业务室流转到分析室进行分析测试。各元素检出限符合国家、行业标准，一级标准物质/重复样合格率均为100%，各元素分析报出率均为100%。样品测试分析方法及检出限见表1。

表1 分析方法、检出限(×10-6，Au、Ag、Hg为×10-9)

分析质量监控包括5项内容，即国家Ⅰ级标准物质(GBW系列)监控、省级二级监控样监控、实验室内部检查、试样重复性密码分析、异常点的抽查检查。内检合格率均大于97%，质量符合要求。

2.3 分析结果及数据处理

5529 件样品分析数据采用SPSS软件进行了分类数据处理和统计分析，按照一次剔除大于平均值加3倍标准离差的数据后，有2624件未蚀变样品、1243件蚀变样品参与数据统计分析(表2)。元素丰度估计值的计算目前仍以算术平均值法为主。

图2 元素均值在各地质单元中的变化曲线图1-胶东岩群；2-荆山群；3-马连庄基性-超基性岩；4-栖霞英云闪长岩；5-谭格庄奥长花岗岩；6-莱州基性-超基性岩；7-玲珑花岗岩；8-郭家岭花岗岩；9-伟德山花岗岩；10-基性脉岩；11-中酸性脉岩Fig.2 A curve of the mean values of elements in each geological unit1-Jiaodong Group; 2-Jingshan Group; 3-Malianzhuang basic-ultrabasic rocks; 4-Qixia tonalite; 5-Tangezhuang trondhjemite; 6-Laizhou basic-ultrabasic rocks; 7-Linglong granite; 8-Guojialing granite; 9-Weideshan granite; 10-basic vein rock; 11-intermediate-acidic vein rock

3 讨论

3.1 胶北隆起不同地质体元素丰度特征

Au等15种元素丰度除Pb外，其他14元素均低于地壳丰度(黎彤, 1992)。中国花岗岩Au等元素丰度(史长义等, 2005)是大量实测数据统计结果，能很好的代表被测岩石的实际丰度(中值)，胶北中生代花岗岩中多数元素丰度与中国花岗岩大致相当或偏低，只有Au大于中国花岗岩丰度3倍，Bi元素明显偏低1倍。胶北地区中生代基性脉岩中Au、W丰度比地区平均丰度明显偏高1倍，其它元素大致相当或稍微偏高。构造岩中Au、Ag、Pb、W、As、Bi比岩石丰度值明显偏高，Co明显偏低，其它元素相当或者稍偏高，说明构造是成矿元素富集的有利场所。Cu、Zn、Co在偏酸性侵入岩中含量明显偏低，在基性岩和胶东岩群中明显偏高，说明其在地壳深部较浅部富集，Pb在中生代酸性侵入体中明显偏高。

图3 胶北隆起不同地质体Au丰度1-胶东岩群；2-荆山群；3-马连庄基性-超基性岩；4-栖霞英云闪长岩；5-谭格庄奥长花岗岩；6-莱州基性-超基性岩；7-玲珑花岗岩；8-郭家岭花岗岩；9-伟德山花岗岩；10-基性脉岩；11-中酸性脉岩；12-地壳丰度(黎彤, 1992)；13-中国花岗岩(史长义等, 2005)Fig.3 Au abundance of different geological bodies in Jiaobei uplift1-Jiaodong Group; 2-Jingshan Group; 3-Malianzhuang basic-ultrabasic rocks; 4-Qixia tonalite; 5-Tangezhuang trondhjemite; 6-Laizhou basic-ultrabasic rocks; 7-Linglong granite; 8-Guojialing granite; 9-Weideshan granite; 10-basic vein rock; 11-intermediate-acidic vein rock; 12-crustal abundance (Li, 1992); 13-Chinese granite (Shi et al., 2005)

不同元素在各地质体中均值变化情况可分为3类(图2)：早前寒武纪地质体、中生代花岗岩和中生代基性脉岩。Ag、Cu、Pb、Zn、Co在中生代花岗岩和早前寒武纪地质体中含量有明显不同，早前寒武纪地质体Ag、Cu、Zn、Co含量高，Pb含量低；中生代花岗岩Ag、Cu、Zn、Co含量低，Pb含量高，W、Sn、Mo、Nb、Ta、As、Sb、Bi、Hg元素则没有明显变化。中生代基性脉岩与中生代花岗岩Ag含量相当，Cu含量处于中生代花岗岩和早前寒武纪地质体之间，Pb含量处于两类地质体之间，其它元素含量差别不明显。因此，中生代地球化学环境继承了早前寒武纪的地球化学环境，但中生代又与早前寒武纪的地球化学环境有明显不同，中生代地壳是混合了大量地幔物质的产物。玲珑花岗岩Au、Ag、Cu含量较郭家岭花岗岩和伟德山花岗岩高，郭家岭花岗岩比玲珑花岗岩和伟德山花岗岩中Pb、Zn、Co含量高、Nb、Ta含量低；伟德山花岗岩Au、Ag、Cu、Pb、Zn、Co元素含量比玲珑花岗岩和郭家岭花岗岩最低，其他元素含量大致相当。

3.2 胶东地区Au元素地球化学背景

胶北隆起Au丰度的总平均值为1.31×10-9，与地壳估算值1.3×10-9(Rudnick and Gao, 2014)一致，侵入岩为1.37×10-9，中酸性侵入岩为1.22×10-9，中生代花岗岩为1.20×10-9(表2)。Au在新太古代谭格庄英云闪长岩、晚侏罗世玲珑花岗岩和中生代基性脉岩中含量明显偏高，分别为1.74×10-9、1.7×10-9、2.65×10-9；新太古代谭格庄奥长花岗岩和早白垩世伟德山花岗岩Au含量明显偏低，分别为0.74×10-9，0.75×10-9；而郭家岭花岗岩丰度处于中间为1.17×10-9(图3)。迟清华(2002)就中国东部火成岩、沉积岩和变质岩中Au含量原始数据经X±3S一次剔除后计算出的算术平均值分别为0.58×10-9、1.01×10-9和1.02×10-9，与本次数据处理方法相同，说明胶北地区变质岩和中国东部变质岩Au丰度相当，胶东地区莱阳群和青山群沉积岩中Au背景值普遍较低(李瑞红等, 2019)，胶北隆起侵入岩Au丰度(1.37×10-9)明显偏高，与胶莱盆地东北缘Au丰度(1.40×10-9)基本一致(李勇等, 2018)，是中国东部火成岩的2倍以上。另据史长义等(2005)的研究，中国东部花岗岩Au含量中值为0.38×10-9(为一次剔除X±2S)，本研究区中生代玲珑花岗岩、郭家岭花岗岩和伟德山花岗岩按照一次剔除X±2S后的中值平均为0.85×10-9，胶北隆起中生代花岗岩Au丰度是中国东部花岗岩的2倍以上。中生代基性脉岩是深部地幔基性岩浆信息的真实反映，其平均Au含量(2.65×10-9)远远高出其他地质体平均含量(1.31×10-9)，而同期中酸性脉岩的Au含量1.22×10-9，与胶北隆起中生代花岗岩平均含量接近，可以代表浅部地壳的Au含量。孙丰月等(1995)论证了金的上地幔源区的存在，由地幔C-H-O流体演化分异而来的煌斑岩中金平均丰度5.7×10-9，赫英等(2004)对胶东半岛地幔岩包体的研究认为胶东地区岩石圈地幔金平均丰度大致等于或略低于8.78×10-9。Wangetal.(2021)的新数据表明前寒武纪地壳基底中Au含量较低，<0.5×10-9，与其他地区的角闪岩相岩石相一致(Pattenetal., 2020; Pitcairnetal., 2015)。前寒武纪变质基底含金量低可以解释为从可能富金的地壳基底提取后的残余结果(Chenetal., 2019)。

2000年之前胶东地区Au丰度数据普遍偏高，特别是胶东岩群的Au丰度更高(Wangetal., 2021)，主要原因是仪器分析精度不高、样品数量偏少和统计时未剔除奇异值。由于不同实验室可能存在系统误差，因此，不同实验室的测试数据的分析对比，仅作参考，同时由于没有足够多的样品，Au丰度仅仅作为一个大致的参考值，但是同一实验室同一时期不同地质体的分析结果完全可以对比。综上，胶北隆起地壳和地幔均具有较高的Au地球化学背景，Au丰度在1.31×10-9左右是可信的，且地幔金丰度高于地壳。

3.3 壳幔岩浆混合反应析出巨量金质

郭家岭花岗岩和伟德山花岗岩均为壳幔混合成因(曲晓明等, 2000; 杨进辉和周新华, 2000; 杨进辉等, 2003; 郭敬辉等, 2005; 周建波等, 2003; Gossetal., 2010; 宋明春等, 2015)，从前述可知，地幔和地壳Au丰度均比郭家岭和伟德山花岗岩高，而混合后的郭家岭花岗岩和伟德山花岗岩Au含量均明显变低，最大的可能就是在地幔岩浆和地壳物质混合形成花岗岩浆演化过程中部分Au流失了。岩浆的结晶过程是一个结晶分异的过程，将岩浆中不能进入矿物晶格中的水分慢慢排出在矿物晶体之外，溶解在其中的不相容元素被带出，仅有很少形成流体包裹体被包裹在矿物中。据朱永峰和安芳(2010)的研究，花岗岩中93%以上的金赋存在硫化物中，在花岗质岩浆演化过程中，硅酸盐矿物的结晶或溶解不影响金的状态，热液就成了金活化迁移的主要载体，花岗岩浆的演化和结晶过程是一个明显的去金作用过程(孙景贵, 1999)，正常情况下，花岗岩中Au丰度要低于其他地质体。硫化物易溶于富水流体(Huang and Keppler, 2015)，铜与金等亲硫元素形成溶于水的络合物，形成富金流体(Sunetal., 2013)。从胶北隆起蚀变带中富含黄铁矿看，成矿热液含S、Fe特别高，说明壳幔岩浆反应演化过程中，从岩浆中析出大量的S、Fe，Au是亲S、Fe元素，Au较容易进入含S、Fe且含有大量Cl、CH4等高温强碱性热液中。很多研究者(朱日祥等, 2015; 杨立强等, 2014) 根据流体氢氧同位素落在岩浆水与大气降水之间，认为很可能与岩浆作用有关, 而氦-氩同位素落在地幔与地壳之间靠近地幔，正说明成矿流体以地幔来源为主，地壳为辅。根据前人碳氧同位素组成，大多数数据点投影在地幔多相体系和原始碳酸岩区外附近右上角花岗岩区内外，显示壳幔混源特征(杨立强等, 2014)，也证实成矿流体是壳幔岩浆混合反应来源。

表3 胶东地区金矿成矿同位素年龄

续表3Continued Table 3矿床矿物年龄(Ma)来源方法玲珑绢云母119.9±1.3Li et al., 2008Rb-Sr等时线玲珑绢云母121.4±2.4Li et al., 2008Rb-Sr等时线玲珑黄铁矿121.0±2.0Li et al., 2008Rb-Sr等时线玲珑钾长石117.5±1.0阳琼艳, 2013Ar-Ar坪大尹格庄绢云母130.71±0.70Yang et al., 2014Ar-Ar坪大尹格庄绢云母130.6±2.4Yang et al., 2014Ar-Ar等时线大尹格庄白云母128.38±0.69Yang et al., 2014Ar-Ar坪大尹格庄白云母130.6±2.4Yang et al., 2014Ar-Ar等时线大尹格庄绢云母132.74±0.83Yang et al., 2014Ar-Ar坪大尹格庄绢云母126.80±0.59Yang et al., 2014Ar-Ar坪大尹格庄绢云母127.7±2.4Yang et al., 2014Ar-Ar等时线夏甸独居石120.0±1.4Ma et al., 2013LA-MC-ICP-MS芝山石英119±5周燕等, 2007Rb-Sr等时线灵雀山石英116±23周燕等, 2007Rb-Sr等时线河西黄铁矿122.3±3.1侯明兰等, 2006Rb-Sr黑岚沟-大柳行黄铁矿117.8±6.5侯明兰等, 2006Rb-Sr庄子独居石119.0±3.1Li et al., 2018LA-ICPMS笏山独居石128.2±2.7Yang et al., 2018LA-MC-ICP-MS笏山独居石120.0±3.1Yang et al., 2018LA-MC-ICP-MS马家窑磷钇矿120±1.2冯凯等, 2017LA-MC-ICP-MS蓬家夼石英118.42±0.25张连昌等, 2002Ar-Ar坪蓬家夼石英117.03±0.13张连昌等, 2002Ar-Ar等时线蓬家夼石英120.53±0.49张连昌等, 2002Ar-Ar坪蓬家夼石英117.33±0.15张连昌等, 2002Ar-Ar等时线蓬家夼黑云母117.49±0.25张连昌等, 2002Ar-Ar坪蓬家夼黑云母116.83±0.36张连昌等, 2002Ar-Ar等时线蓬家夼绢云母119.5±0.3Li et al., 2006Ar-Ar坪蓬家夼绢云母119.7±0.2Li et al., 2006Ar-Ar坪蓬家夼绢云母119.1±0.2Li et al., 2006Ar-Ar坪蓬家夼绢云母119.6±0.4Li et al., 2006Ar-Ar坪蓬家夼绢云母120.8±0.4Li et al., 2006Ar-Ar坪蓬家夼绢云母120.3±0.5Li et al., 2006Ar-Ar坪蓬家夼绢云母120.9±0.4Li et al., 2006Ar-Ar坪西涝口绢云母116.5±0.61孙兴丽, 2014Ar-Ar坪西涝口绢云母116.62±0.75孙兴丽, 2014Ar-Ar等时线西涝口绢云母116.62±0.79孙兴丽, 2014Ar-Ar反等时线西涝口黑云母121.64±0.68孙兴丽, 2014Ar-Ar坪西涝口黑云母121.8±1孙兴丽, 2014Ar-Ar反等时线西涝口黑云母121.83±0.94孙兴丽, 2014Ar-Ar等时线辽上白云石104.8±5.1Li et al., 2019Sm-Nd等时线辽上黄铁矿105.5±9Li et al., 2019Rb-Sr等时线发云夼黄铁矿128.17±7.2张连昌等, 2002Rb-Sr等时线郭城石英116.2±2.4Tan et al., 2012Rb-Sr郭城黄铁矿119±2李杰等, 2020Rb-Sr等时线英格庄绢云母120.02±0.38陈炳翰, 2017Ar-Ar坪

图4 胶东地区金矿成矿年龄直方图Fig.4 Age histogram of gold deposits in Jiaodong area

图5 地质体元素标准离差1-胶东岩群；2-荆山群；3-马连庄基性-超基性岩；4-栖霞英云闪长岩；5-谭格庄奥长花岗岩；6-莱州基性-超基性岩；7-玲珑花岗岩；8-郭家岭花岗岩；9-伟德山花岗岩Fig.5 The standard deviation of the element in geological mass1-Jiaodong Group; 2-Jingshan Group; 3-Malianzhuang mafic and ultrabasic rocks; 4-Qixia tonalite; 5-Tangezhuang trondhjemite; 6-Laizhou mafic and ultrabasic rocks; 7-Linglong granite; 8-Guojialing granite; 9-Weideshan granite

图6 地质体元素变异系数1-胶东岩群；2-荆山群；3-马连庄基性-超基性岩；4-栖霞英云闪长岩；5-谭格庄奥长花岗岩；6-莱州基性-超基性岩；7-玲珑花岗岩；8-郭家岭花岗岩；9-伟德山花岗岩Fig.6 Variation coefficient of elements in geological bodies1-Jiaodong Group; 2-Jingshan Group; 3-Malianzhuang mafic and ultrabasic rocks; 4-Qixia tonalite; 5-Tangezhuang trondhjemite; 6-Laizhou mafic and ultrabasic rocks; 7-Linglong granite; 8-Guojialing granite; 9-Weideshan granite

虽然很多学者认为成矿物质来源于地幔，但是迄今未发现成矿热液直接来源于深部地幔的通道，壳幔岩浆的混合反应恰好解决了地幔物质怎么上升到地壳浅部的问题。含矿热液的上升依附于岩体的演化和上升，此现象在斑岩铜矿中更为明显。壳幔岩浆反应的机制、反应前后物质成分变化、反应的物理化学条件变化等还需深入研究。

胶北隆起发育大范围的中生代花岗岩，分布广泛，侵入持续时间长，从晚侏罗世到晚白垩世晚期，长达近100Myr，伟德山花岗岩的分布面积大致在2500km2以上，地幔、地壳Au含量分别为2.6×10-9和1.3×10-9，郭家岭和伟德山花岗岩Au含量分别为1.17×10-9和0.75×10-9。两种花岗岩均比地壳和地幔含量低，那么壳幔物质充分混合形成两类花岗岩浆过程中，大部分金可能被溶解有大量S、Cl等阴离子的H2O-CO2-CH4热液体系捕获溶解，形成含矿热液，含金热液逐渐汇聚向上运移，在岩体上方合适的构造部位成矿。如中生代两种花岗岩平均厚度按5km、平均值与壳幔平均值差近1×10-9估算，那么将有33600t金被析出。

3.4 中生代玲珑、郭家岭和伟德山三期花岗岩岩浆演化与成矿的耦合

胶东地区金矿总体成矿时间区间在132～105Ma之间(图4、表3)，郭家岭花岗岩的形成时间为130～125Ma(关康等, 1998; 杨进辉等, 2000; Wangetal., 1998; Yangetal., 2012; 罗贤冬等, 2014; Songetal., 2020)，伟德山花岗岩的形成时间为126～108Ma(周建波等, 2003; 郭敬辉等, 2005; 张田和张岳桥, 2008; Gossetal., 2010; 丁正江等, 2013; Songetal., 2020)，成矿时间恰好与两期花岗岩浆的冷凝结晶时间吻合。早期成矿时间与郭家岭花岗岩浆开始活动的时间相对应，即130Ma左右，如招远市大尹格庄金矿成矿年龄为130±4Ma(杨立强等, 2014)，显示该时期有明显的构造流体成矿活动，说明郭家岭花岗岩浆还未结晶时就对其上部地质构造和成矿有影响，郭家岭岩浆期的热液就已经开始上升成矿。而最晚的成矿时间为105Ma，比伟德山花岗岩最晚固结时间(108Ma)稍晚。这也间接说明郭家岭和伟德山花岗岩与成矿有一定的关系。

根据锆石年龄，玲珑花岗岩锆石结晶持续时间是24Myr，郭家岭花岗岩锆石结晶持续时间是5Myr，伟德山花岗岩锆石结晶持续时间是18Myr，从锆石结晶持续时间长短可以看出上述三类花岗岩岩浆冷凝结晶持续时间长短。这说明郭家岭花岗岩浆规模小，上升冷凝快，伟德山花岗岩浆规模大，上升冷却慢，能持续提供成矿热液。玲珑花岗岩中存在大量捕虏锆石，而郭家岭和伟德山花岗岩中存在很少量的捕虏锆石，说明壳源重熔的玲珑花岗岩熔融温度低，锆石等高温矿物没有被熔化，其持续时间长说明其规模大且上升缓慢。壳幔岩浆混合反应可能是岩浆结晶分异析出成矿物质的一种方式，也是地幔或下地壳成矿物质伴随岩浆上升到地壳浅部成矿的一种机制。大规模的地幔岩浆上升，才能携带大量深部成矿物质和热液一起上升。

胶东地区玲珑花岗岩分布范围广，赋存的金矿数量多。从成矿时间看，莱州市留村金矿成矿时间为151±3Ma(姜晓辉等, 2011)，而此时正是玲珑花岗岩侵位时间164～140Ma(郭敬辉等, 2005; Huetal., 1987; Maetal., 2013; Miaoetal., 1998; Wangetal., 1998; Yangetal., 2012; Zhangetal., 2010; Songetal., 2020)，矿区发育玲珑花岗岩体，周围未见其它中生代侵入岩，显然这次成矿与玲珑花岗岩关系密切。据徐述平等(2008)的研究，玲珑花岗岩成岩过程中，在岩体的边缘相Au、Pb、Bi出现了明显富集，Au元素的初步富集为后期大规模成矿作用做好了物质上的准备，具备一定的矿源层特征，同时有小规模的金矿成矿。留村金矿虽然是一个孤证，由于其远离中生代其它岩体，足以说明玲珑花岗岩与成矿关系密切，其它地方尚未发现较早的成矿，也可能是早期的成矿被后期大规模的成矿叠加改造而难测其早期成矿时间。

图7 部分元素的Q-Q正态分布图1～2-谭格庄奥长花岗岩；3～5-莱州基性-超基性岩；6-伟德山花岗岩. 纵坐标为期望的坐标，横坐标为观测值Fig.7 Q-Q normal distribution of the partial elements1～2-Tangezhuang trondhjemite-tonalite; 3～5-Laizhou mafic-ultramafic rocks; 6-Weideshan granite. The ordinate is standard of expectation, the abscissa is observed value

综上所述，中生代玲珑、郭家岭和伟德山三期规模较大的花岗岩的形成和演化是形成巨型金矿的基础条件，玲珑花岗岩呈大致层状体，其原始金丰度偏高，特别是边部，有小规模石英脉型金矿成矿，具有一定的初始矿源层特征，是较好的赋矿地质体。后期壳幔混合反应形成的郭家岭和伟德山花岗岩浆的演化，为大规模成矿提供了巨量的金质来源。似层状玲珑花岗岩中的金矿成矿与层控矿床类似，有层控矿床的特征，可称为“内生层控矿床”或“类层控矿床”(田杰鹏等, 2016)，是赋矿地质体为侵入岩的非标准层控矿床。

3.5 不同地质体对金矿成矿的有利度分析

从元素标准离差图上(图5)看出，15种元素离差在新太古代马连庄基性-超基性岩和早白垩世伟德山花岗岩中普遍偏低。从变异系数图上(图6)看出，15种元素中多数在新太古代马连庄基性-超基性岩、古元古代莱州基性-超基性岩、伟德山花岗岩中偏低，其中Au、Ag、Pb、Zn在上述地质体中均为较低变异系数。中生代花岗岩Au元素标准离差从高到低依次是玲珑花岗岩(2.38)、郭家岭花岗岩(0.75)、伟德山花岗岩(0.13，图5)；Au元素变异系数从高到低依次是玲珑花岗岩(1.40)、郭家岭花岗岩(0.65)、伟德山花岗岩(0.18)；Au元素含量极差从大到小依次为玲珑花岗岩(19.89)、郭家岭花岗岩(5.21)、伟德山花岗岩(0.78)。从标准离差和变异系数来看，对金矿成矿有利地质体依次为玲珑花岗岩、郭家岭花岗岩、伟德山花岗岩。

表4 服从正态分布的元素P值(置信度95%)

通过前述中生代花岗岩演化与金矿成矿耦合关系分析，花岗岩岩浆形成后有成矿元素进入热液并上升到上覆地质体中，就会增高上覆地质体中Au元素丰度，同时也使其标准离差和变异系数增大。郭家岭花岗岩的上覆围岩以玲珑花岗岩为主，其形成演化过程中使玲珑花岗岩Au元素丰度增高、标准离差和变异系数增大；同样，伟德山花岗岩的上覆围岩以玲珑花岗岩和郭家岭花岗岩为主，其形成演化过程中使玲珑花岗岩Au元素丰度进一步增高、标准离差和变异系数进一步增大，使郭家岭花岗岩Au元素丰度增高、标准离差和变异系数增大。即深部后期的花岗岩浆含矿热液(包括从地幔带上来的含矿热液)使上部早期的地质体又经历了一次地球化学过程，使得其中很多元素含量不均匀增加，造成了其标准离差和变异系数的增大，使上覆地质体更容易成矿。这从成矿热液来源的角度解释了胶北隆起Au元素标准离差和变异系数的大小与赋矿的正相关关系，因此，标准离差和变异系数是判断赋矿地质体的有效参数，是化探找矿的很好指标。

Vistelius (1960)认为，单一地球化学过程所形成的地质体，元素含量服从正态分布，由数个地球化学作用过程叠加所形成的复合地质体中元素含量偏离正态分布，并且多为正偏分布(其中有些服从对数正态分布)，他将这一结论称为“地球化学过程的基本定律”。一般来说，成矿作用总是出现在地质构造复杂，地球化学作用多次叠加的地区，因此，不服从正态分布的地质体，才具有找矿前景，地球化学作用叠加越多的地质体越是找矿有利的对象。

对元素进行正态性检验，是为了区别不同地质体对成矿的贡献大小。SPSS规定：当样本含量3≤n≤5000时，结果以Shapiro-Wilk(W检验)为准，当样本含量n>5000时，结果以Kolmogorov-Smirnov(D检验)为准，在SPSS软件统计结果中，P≥0.05时，认为是服从正态分布。同时再结合Q-Q图进一步确认，这是可以检验元素是否服从正态或对数正态分布的常用方法之一(成秋明, 2000)，如基本在直线附近，可以认为服从正态分布。各地质单元分析数据量<5000，采用SPSS软件处理时，置信度采用95%，采用W检验，大部分元素不服从正态分布，仅新太古代谭格庄变质奥长花岗岩的Pb、Co、古元古代莱州基性-超基性岩的Cu、Co、Mo、伟德山花岗岩的Zn的P值≥0.05(表4)，相应的在Q-Q图上均在直线附近(图7)。因此，上述元素在相应的地质体中可以认为是服从正态分布，而其他元素在各个地质单元中不服从正态分布，说明这3个地质体中的相应元素受后期地球化学影响较小。荆山群、新太古代栖霞英云闪长岩、新太古代马连庄基性-超基性岩、玲珑花岗岩、郭家岭花岗岩和构造岩无一元素服从正态分布，说明上述地质单元相对经历的地球化学影响较大，构造活动、岩浆活动都不同程度对其前期的地质体产生一定的地球化学影响。构造岩中也无一元素符合正态分布，其发生了热液蚀变，发生了大量元素的带入和带出，受地球化学影响是最大的。通过成矿元素和微量元素的标准离差、变异系数和正态性检验显示，中生代玲珑花岗岩和郭家岭花岗岩对金矿赋矿最为有利，伟德山花岗岩赋矿相对较差，但伟德山花岗岩提供了大量的金质，对成矿的贡献不容忽视。

4 结论

(1)胶北隆起壳幔岩浆混合反应是产生富含金质热液的主要形式，是巨量金质来源的主要机制，也是深部岩浆成矿物质上升的载体。首先壳幔混源的郭家岭花岗岩和伟德山花岗岩金丰度均低于地壳和地幔丰度的原因，可能是壳幔岩浆混合反应过程中，岩浆中的金质被高H2S-H2O-CO2-CH4的岩浆热液体系捕获带走所致。其次胶东金矿形成的主要成矿期为130～105Ma，与郭家岭(130～125Ma)和伟德山(126～108Ma)两期花岗岩岩浆演化侵位时间吻合，伴随该两期花岗岩岩浆演化和侵入全过程，其岩浆期含矿热液上升沉淀成矿。

(2)胶北隆起具有较高的Au背景场，Au丰度为1.31×10-9，中生代花岗岩丰度为1.20×10-9，晚侏罗世玲珑花岗岩为1.7×10-9，早白垩世郭家岭花岗岩为1.17×10-9，早白垩世伟德山花岗岩为0.75×10-9。胶北隆起中生代继承了早前寒武纪地球化学环境，但又与其有明显不同，中生代地壳是混合了大量地幔物质的产物。

(3)岩浆混合反应演化过程中形成的含矿热液上升到上覆早期地质体中，对部分元素丰度和地球化学参数会产生一定影响，造成元素丰度、标准离差和变异系数增大。胶北隆起玲珑花岗岩金丰度明显高于郭家岭和伟德山花岗岩，除其原始丰度高以外，与后期两种花岗岩的含矿热液上升“浸染”有很大关系。中生代玲珑花岗岩和郭家岭花岗岩对金矿赋矿有利，伟德山花岗岩赋矿相对较差，但对成矿提供了大量的金质和源源不断的热量。