胶-辽-吉造山带辽河群石墨矿碳同位素特征及成因分析*

朱建江 刘福来 刘福兴 石闯 王舫 徐文涛

1.自然资源部深地动力学重点实验室，中国地质科学院地质研究所，北京 1000372.辽宁省有色地质一〇三队有限责任公司，丹东 118008

石墨由于其特殊的结构特征，在军事及工业材料上有着广泛的应用，随着新能源汽车及石墨烯高科技材料的出现，石墨资源的地位越来越受到社会的关注(Crossley, 2000; Novoselovetal., 2004; Singhetal., 2011; Luqueetal., 2014; 肖克炎等, 2016)。中国是世界最大的石墨生产国，也是国际市场的最大出口国，因此石墨矿床的形成、分布规律和勘查远景是我国矿产资源研究的重要问题之一(陈毓川等, 2010; 王登红等, 2013; 李超等, 2015)。全球石墨矿床分布相对集中，主要分布在中国、印度、巴西、捷克、加拿大和墨西哥等少数几个国家。据美国地质调查局数据，截至2019年，全球已探明的天然石墨储量约为3亿吨，其中中国储量约为7300万吨，占全球的1/4(http://minerals.usgs.gov)。中国石墨储量虽然位居世界第一，但石墨矿床分布广泛，品位参差不齐，研究程度较低。针对中国石墨矿产资源的勘探开发研究现状，深入探讨石墨矿床的基本特征、成因机制以及成矿分布规律，对于我国石墨矿床的勘探以及石墨产业的发展具有重要意义(陈毓川等, 2010; 王登红等, 2013)。

中国石墨矿床具有一定的时空分布特征，在时间上表现为石墨主要形成于元古代和太古代，少量形成于古生代和中、新生代；在空间上石墨矿主要分布在古老陆块边缘，受老地块制约，形成了分布相对集中“东多西少”的空间格局(李超等, 2015; 肖克炎等, 2016)。华北克拉通周缘的晶质石墨矿床占全国储量的74%，主要位于西北缘丰镇活动带、南缘晋豫活动带和东缘胶-辽-吉活动带，形成华北克拉通周缘三条石墨成矿带(李超等, 2015; 肖克炎等, 2016; Zhongetal., 2019)。尽管这些石墨矿从20世纪50年代就开始开采，但对于这些石墨矿的成因机制、物质来源、形成时代和分布规律仍不是很清楚，对这些问题的解答可以为石墨矿床资源形成规律及勘察提供依据。

陈衍景等(2000)对其中的代表性地层(西北缘的丰镇群、南缘的水滴沟群和东缘的荆山群和辽河群)进行了总结，表明这些含石墨矿的地层沉积时代为古元古代(2300～2050Ma)，经历了后期19～18亿年的变质事件。岩相学特征表明以上地区的石墨矿床的矿石类型基本一致，主要为石墨片麻岩、石墨透辉岩、石墨大理岩和混合岩化石墨片麻岩等4种类型，碳同位素表明石墨矿具有有机和无机碳两种来源，以有机为主。近年来，有研究者对华北克拉通西北缘的丰镇石墨成矿带进行了详细的研究(Yangetal., 2014; Zhongetal., 2019)，提出该地区的石墨矿沉积于活动大陆边缘环境，在后期碰撞造山过程中有机质变质并逐渐富集成石墨矿。相对于丰镇石墨成矿带，胶-辽-吉石墨成矿带的研究比较薄弱。兰心俨(1981)对该带中的南墅石墨矿床进行了碳同位素方面的研究，认为石墨的碳主要来自于有机质沉积，很少无机碳的加入。李凯月等(2018)对胶北荆山群张舍石墨矿研究表明其石墨碳质来源以有机质为主，并与无机碳混合，产生了均一化，导致碳同位素变重；此外，石墨拉曼光谱峰特征指示张舍石墨矿经历了高角闪岩相-麻粒岩相的变质作用。Wangetal.(2020b)对胶东两处石墨矿(刘各庄石墨矿和大梁子口石墨矿)碳同位素研究表明刘各庄石墨矿碳源主要为沉积的有机物，而大梁子口石墨矿形成过程还有无机碳的参与。然而到目前为止，对于胶-辽-吉带北部的辽河群石墨矿床的碳来源、形成时代、成矿机制和分布规律仍不是很清楚。

本文以胶-辽-吉造山带中的辽河群石墨矿(甜水乡马沟石墨矿)为研究对象，对其进行了详细的岩相学、微量元素地球化学、拉曼光谱学和碳同位素等方面的研究，限定了石墨矿的碳来源和成矿机制，并对其形成时代以及对地球古环境的启示进行了探讨。

图1 华北克拉通东部陆块胶-辽-吉造山带地质简图(据Zhao et al., 2005修改)Fig.1 Simplified geological map of the Jiao-Liao-Ji belt in the Eastern Block, North China Craton (after Zhao et al., 2005)

1 地质背景及样品岩石学特征

1.1 区域地质概况

胶-辽-吉带位于华北克拉通的东部，呈北东-南西向展布，长约1000km、宽约50～300km，向东延伸到朝鲜半岛境内，向西可能穿越郯庐断裂延伸至徐州-蚌埠一带(图1，Zhaoetal., 2005, 2012; Caietal., 2020)。该带经历了十分复杂的构造演化过程，具有非常复杂的物质组成，记录了多期岩浆-变质事件。北部的龙岗地块主要由大量分布的新太古代英云闪长岩-奥长花岗岩-花岗闪长岩(TTG)(白瑾, 1993; 万渝生等, 2001; Wanetal., 2005)及少量变质表壳岩(鞍山群, 辽宁省地质矿产局, 1989; 翟明国等, 1990)组成。南部狼林地块主要由一系列新太古代闪长岩-英云闪长岩-花岗闪长岩侵入体组成(辽宁省地质矿产局, 1989; 陆孝平等, 2004; Zhaoetal., 2006)，最近的同位素年代学资料表明，狼林地块主要由古元古代(1.8～1.9Ga) 岩石组成，与辽-吉岩系类似，因此提出狼林地块可能属于胶-辽-吉造山带的一部分(吴福元等, 2016)。胶-辽-吉带主要由古元古代巨量变质火山-沉积岩系、多期具有不同成因的花岗质岩石(A 型花岗质片麻岩、碱性花岗岩、钙碱性花岗岩及环斑花岗岩等)以及大量变基性侵入体(辉长岩和辉绿岩等)组成(Lietal., 2004; Luoetal., 2004, 2008;陆孝平等, 2004; Luetal., 2006; Li and Zhao, 2007; 刘福来等, 2015; Wangetal., 2017a, b; Xu and Liu, 2019)。该带古元古代变质火山-沉积岩系主要包括吉南地区的吉安群和老岭群、朝鲜半岛的摩天岭群、辽东地区的南辽河群和北辽河群、胶北地区的荆山群和粉子山群(辽宁省地质矿产局, 1989; Zhaoetal., 2005; 赵国春, 2009)。近年来的研究表明胶-辽-吉带可以向南延伸至蚌埠地区的五河群(Guo and Li, 2009; Caietal., 2020)。有研究者根据岩石组合、岩浆作用、构造变质作用的差异性，将胶-辽-吉带分为由南辽河群、吉安群和荆山群组成的南部带(贺高品和叶慧文, 1998a, b; Lietal., 2005; Zhaoetal., 2005, 2012; Zhouetal., 2008)以及由北辽河群、老岭群和粉子山群组成的北部带(李三忠等, 2001; Zhaoetal., 2005; Luetal., 2006; Luoetal., 2008)。

图2 胶-辽-吉带甜水乡地区(辽河群石墨矿)地质简图(据辽宁省第一区域地质测量队，1975(1)辽宁省第一区域地质测量队. 1975. 1:20万辽阳幅地质图修改)

1.2 辽河群变质火山-沉积岩系

辽河群是由斋藤林次于1938年所建立的“辽河系”演变而来(转引自辽宁省地质矿产局1989)，主要出露于辽东南地区，不整合于由太古宙鞍山群和连山关花岗岩组成的变质基底之上，自下而上可以划分为浪子山组、里尔峪组、高家峪组、大石桥组和盖县组，其上部被中元古代榆树砬子组不整合覆盖(辽宁省地质矿产局, 1989)。辽河群在空间上被盖县-析木城-塔子岭-叆阳界线分为北辽河群和南辽河群(王惠初等, 2011, 2015)，它们分别可以和胶北地区的粉子山群和荆山群对比(Zhaoetal., 2005, 2012)。北辽河群主要由石英岩、板岩、片岩、千枚岩和大理岩等岩石组成，这些变沉积岩石大多数记录了顺时针的变质P-T轨迹(贺高品和叶慧文, 1998a, b; Lietal., 2005; Zhaoetal., 2005, 2012; Zhouetal., 2008)。南辽河群的主要组成岩石为含石墨黑云母片麻岩、夕线-石榴黑云母片麻岩、石榴石-十字石片岩和大理岩，其大多记录了逆时针的变质P-T轨迹(贺高品和叶慧文, 1998a, b; 李三忠等, 2001; Zhaoetal., 2005; Luetal., 2006; Luoetal., 2008)。相对于北辽河群，南辽河群缺失浪子山组，其中的里尔峪组和高家峪组经历了更高的变质作用以及更明显的混合岩化作用。最近的研究发现在南辽河群和吉安群中的岩石也存在有顺时针的变质P-T轨迹(Caietal., 2017; Liuetal., 2019)；此外，同位素年代学研究表明南辽河群和北辽河群形成于相同的构造环境，并且具有相同的物质来源(Wangetal., 2017a, 2020a)。

辽河群浪子山组的主体岩性为砾岩、石英岩、含石榴石云母片岩、含石墨白云石英片岩、石墨白云(二云)长英质粒状岩石和千枚岩等。里尔峪组和高家峪组主要由变质火山岩、片岩、含石榴石磁铁矿云母片岩、斜长角闪岩以及少量的大理岩组成，其中里尔峪组中发育大量硼矿和黄铁矿，高家峪组中发育有磷矿和石墨矿(刘福来等, 2015; Tianetal., 2017)。大石桥组是一套以碳酸盐岩为主的地层，自下而上共分为三个岩性段，一段岩性主要为方解石大理岩和白云石大理岩，夹透闪大理岩和透闪岩；二段主要由石榴十字石云母石英片岩、钙质黑云变粒岩和条带状大理岩夹白云质大理岩等组成；三段主要岩性为厚层菱镁矿和白云质大理岩，夹薄层千枚岩和板岩，大石桥组中发育有大量的滑石、菱镁矿和岫岩玉等矿床(陈从喜, 2000; Chen and Cai, 2000; 蒋少涌等, 2004)。位于辽河群上部的盖县组的岩石组合主要为变质砂岩、千枚岩和变质粉砂岩(刘福来等, 2015; Tianetal., 2017)。

1.3 辽河群石墨矿产出层位及岩石学特征

辽河群石墨矿主要产于高家峪组中。北辽河群高家峪组底部的主要岩石组成为二云母片岩、含石榴石二云母石英片岩、黑云母片岩和千枚岩等；中部以白云石大理岩、透辉-透闪石大理岩以及碳质方解石大理岩为主；上部的岩性主要为黑色碳质泥沙质板岩，夹有含碳质石英方解大理岩及含碳质凝灰岩等。南辽河群高家峪组以含石墨为其主要特征，主要由含石墨黑云片岩-片麻岩、含石墨透闪长英质粒状岩石、含石墨大理岩、石墨透闪岩、夕线黑云斜长片岩-片麻岩、斜长角闪岩、二云母片岩等岩石组成。本文研究的石墨矿样品采自于北辽河群的甜水乡马沟地区(高家峪组)，取自两口钻井岩芯(ZK3-7和ZK7-4；由辽宁省有色地质局一〇三队提供)，岩芯采样区地质简图及位置见图2，代表性的石墨钻孔岩芯手标本照片见图3a, b。

含石墨矿的主要岩石类型有：含石墨变质杂砂岩、含石墨黑云母长英质片麻岩和含石墨透闪大理岩。其岩相学特征如下：含石墨变质杂砂岩，主要由石英、长石、透闪石、石墨，以及少量黑云母和黄铁矿等组成，其中石墨主要呈片状、条带状，少量为粒状(图3c, d)。含石墨黑云母长英质片麻岩主要由石英、长石、黑云母和石墨组成，还含有少量的黄铁矿，其中长石部分绢云母化，石墨以粒状为主，少量片状(图3e, f)。含石墨透闪大理岩中的主要矿物为方解石、透闪石、石墨，以及少量的石英和黄铁矿，石墨主要呈粒状(图3g, h)。

2 分析方法

2.1 全岩主量、微量元素分析

本文挑选了代表性的含石墨岩石样品进行全岩地球化学测试分析。首先手标本剥掉表面的风化层，用去离子水洗净，打碎成细小颗粒后，选择没有裂隙和脉体的部分，再次用去离子水洗净，然后在玛瑙研钵中研磨成 200 目以下的粉末；接着在澳实分析检测(广州)有限公司进行全岩主量和微量元素的测试分析。主量元素采用XRF(X射线荧光光谱仪PANalytical Axios)方法测试，烧失量(LOI)通过样品在1000℃下灼烧后的质量差得到，主量元素的分析误差(1σ)小于1%。对于微量元素分析步骤如下：首先将大约0.2g样品和0.9g偏硼酸锂混合均匀，1000℃下熔化，冷却之后溶解在100mL 4% HNO3/2% HCl溶液中；然后再使用电感耦合等离子体质谱仪(Perkin Elmer Elan 9000)对溶液进行测试。标样选用GSR-1和GSR-3，测试结果表明大多数微量元素的实验相对误差小于10%。

[5]芦千文、张益：《对现代农业产业化联合体发展的调查与思考——安徽省宿州市为例》，《农业经济与管理》2017年第2期。

图4 辽河群石墨矿log(Fe2O3/K2O)-log(SiO2/Al2O3)分类图解(a, 据Herron, 1988)和Th/Sc-Zr/Sc判别图解(b, 据McLennan et al., 1993)Fig.4 Log(Fe2O3/K2O) vs. log(SiO2/Al2O3) (a, after Herron, 1988) and Th/Sc vs. Zr/Sc (b, after McLennan et al., 1993) classification diagrams of the graphite deposits in the Liaohe Group

2.2 拉曼光谱学分析

本文中的拉曼光谱学分析测试工作是在中国地质科学院地质研究所显微激光拉曼实验室完成的。拉曼光谱仪型号为Horiba LabRAM HR Evolution，搭载Olympus BX41显微镜，使用100倍物镜，光源为double Nd: YAG laser，波长532nm，强度80mW，光谱分辨率为±0.59cm-1，激光束斑为1～2μm，信号采集时间为5～10s，每次测试前用单晶硅进行校正。测试石墨的拉曼光谱时，为了防止激光发热对石墨拉曼谱峰的影响，我们将激光强度设为1mW，测试过程参照Beyssacetal. (2003)的步骤。

2.3 石墨碳稳定同位素分析

本文中的石墨碳同位素分析是在中国科学院地质与地球物理研究所稳定同位素地球化学实验室测试完成的，使用的仪器为ThermoFisher 253质谱仪，搭配 GasBench II系统。测试岩石样品中石墨的碳同位素过程如下：首先称量一定量的样品粉末(200目)和6N HCl反应24h，确保样品中的碳酸岩矿物被完全反应掉。再经过中和、干燥之后将样品放入锡胶囊中，将装有样品的锡胶囊封住放入EA自动送样机中灼烧，将产生的CO2导入质谱仪中测试碳同位素值，标样为已知碳同位素的尿素(IVA)以及石墨(GBW04407)，测量误差为 0.15‰。其中碳同位素用南卡罗莱纳州白垩纪皮迪建造中的箭石进行标准化(VPDB，δ13Ccarb)。

3 测试结果

3.1 含石墨样品全岩地球化学特征

图5 辽河群石墨矿原始地幔标准化微量元素蛛网图(a)和球粒陨石标准化稀土元素配分图(b)(标准化值据Sun and McDonough, 1989)上地壳组分引自Rudnick and Gao (2003)；PAAS-太古宙后澳大利亚平均页岩(Taylor and McLennan, 1985)；E-MORB-富集洋中脊玄武岩；N-MORB-正常洋中脊玄武岩Fig.5 Primitive mantle-normalized trace element diagram (a) and chondrite-normalized REE patterns (b) for the graphite deposits in the Liaohe Group (normalization values after Sun and McDonough, 1989)The data of upper continental crust are from Rudnick and Gao (2003); PAAS-post-Archean Australian average shales (Taylor and McLennan, 1985); E-MORB-enriched mid-ocean ridge basalt; N-MORB-Normal mid-ocean ridge basalt

辽河群含石墨岩石样品具有相似的微量元素特征(图5)，但不同微量元素的含量具有一定的差异性。其中大离子亲石元素(LILEs)含量Rb (22.7×10-6～170×10-6)、Ba(121×10-6～734 ×10-6)、Sr(51.7×10-6～190.5×10-6)具有较宽的变化范围，Rb (87.1×10-6)、Ba(251.9 ×10-6)、Sr(94.1×10-6)的平均值均小于太古宙后澳大利亚平均页岩(PAAS；Taylor and McLennan, 1985)中Rb、Ba和Sr的含量(表1)，表明辽河群含石墨岩石在后期变质作用过程中，大离子亲石元素发生了元素迁移作用。这些样品的高场强元素含量(HFSEs：Nb=5.9×10-6～11.1×10-6, Ta=0.55×10-6～0.95×10-6, Zr=109×10-6～176×10-6, Hf=2.7×10-6～4.6×10-6)变化范围较窄，和PAAS以及上地壳成分特征相似，表明它们在后期变质作用过程中，高场强元素未发生明显的迁移。在原始地幔标准化微量元素图解中(图5a)，这些样品具有Th、U元素正异常和Nb、Ta、Sr、Ti元素负异常的特征，和上地壳组分的微量元素特征相似(Rudnick and Gao, 2003)。在球粒陨石标准化稀土元素配分图中(图5b)，这些样品具有轻稀土(LREE)富集[(La/Yb)N= 4.9～10.3]和Eu负异常(Eu/Eu*=0.51～0.78)的特征，和PAAS和上地壳的稀土元素特征相似(Taylor and McLennan, 1985; Rudnick and Gao, 2003)。

3.2 石墨拉曼光谱学特征

碳质物质的拉曼光谱峰可以分为两个区域，一级谱峰区(1100～1800cm-1)和二级谱峰区(2500～3100cm-1)(Tuinstra and Koenig, 1970; Nemanich and Solin, 1979)。石墨拉曼光谱峰在一级谱峰区的主峰为G峰 (≈1580cm-1)，

表1 辽河群含石墨岩石样品主量(wt%)、微量(×10-6)元素数据

表2 辽河群含石墨岩石样品石墨碳同位素测试结果

石墨结晶度越高，G峰越尖锐(Beyssacetal., 2002)。此外，石墨拉曼光谱峰在一级谱峰区还存在D1(≈1350cm-1)和D2(≈1620cm-1)两个缺陷峰，这两个峰的面积随着石墨结晶度的升高而逐渐降低(Nemanich and Solin, 1979; Beyssacetal., 2002)。在二级谱峰区，石墨拉曼光谱峰的位置主要位于2700cm-1附近，为S1峰，随着石墨结晶度的升高，S1谱峰对称度降低，其谱峰可以进一步分解为两个小谱峰(Wopenka and Pasteris, 1993; Beyssacetal., 2003; Reich and Thomsen, 2004)。

我们对辽河群含石墨岩石样品中的石墨进行了详细的显微拉曼光谱学分析，代表性结果如图6所示。辽河群石墨矿中的石墨整体都具有非常尖锐的G峰和微弱的D1、D2缺陷峰，二级谱峰区的S1峰都显示不对称的特征，并且可以分解为两个小的谱峰，表明辽河群石墨矿中的石墨都具有比较高的结晶度(Beyssacetal., 2002)。此外，不同含石墨岩石样品中的石墨拉曼光谱峰具有一定的差异性(图6)，而且同一样品不同石墨颗粒的拉曼光谱峰也存在一定的差异(图6a, b)。峰面积比值R2[R2=D1/(G+D1+D2)]位于0.03～0.20区域内，峰强度比值R1 (R1= D1/G)为0.03～0.28。

3.3 石墨碳同位素特征

辽河群含石墨岩石样品(钻井岩芯ZK3-7、ZK7-4)石墨碳同位素测试结果列在表2以及投在图7中。这些样品的石墨碳同位素具有比较宽的变化范围(δ13CPDB=-16.49‰～-25.93‰)，和华北克拉通周缘其它地区(佳木斯、华北克拉通东南缘、华北克拉通孔兹岩带和内蒙古中部)以及印度南部Kerala孔兹岩带的石墨矿具有相似的碳同位素特征(图8)。其中钻井岩芯ZK3-7的样品石墨碳同位素值δ13CPDB(‰)位于-16.49‰～-25.72‰区间，石墨碳同位素的值随着样品的埋藏深度先降低后升高。取自钻井岩芯ZK7-4的样品石墨碳同位素值和ZK3-7相似(δ13CPDB=-17.39‰～-25.93‰)，其石墨碳同位素的值也具有随着深度先降低后升高的特征。不同含石墨岩石样品具有相似的石墨碳同位素特征(表2和图7)。

4 讨论

4.1 石墨结晶度和形成温度

有机物转变为石墨的过程主要受温度控制，随着变质温度的升高，有机碳的结晶度逐渐升高(Landis, 1971; Wadaetal., 1994; Nishimuraetal., 2000; Beyssacetal., 2002)。石墨的拉曼光谱峰对石墨结构的变化非常敏感(Wopenka and Pasteris, 1993)，其中G峰(1580cm-1)主要和石墨晶体中的E2G2震动模式(在石墨芳香烃结构平面内的震动)有关(Beyssacetal., 2002)。此外，由于在石墨晶体的层间存在杂环原子(例如O、H、N)或结构缺陷，在石墨拉曼光谱的一级谱峰区常出现另外两个缺陷峰D1和D2峰，这两个缺陷峰强度越高，石墨的结晶度越低(Wopenka and Pasteris, 1993; Beyssacetal., 2003)。石墨拉曼光谱在二级谱峰区的谱峰主要受谐波和组合衍射影响，石墨晶体由二维到三维转变过程中会使得二级谱峰区S1峰分解为两个小峰，反应其具有高结晶度特征(Wopenka and Pasteris, 1993; Beyssacetal., 2003; Reich and Thomsen, 2004)。由于石墨的形成是不可逆的，其结晶度不受退变质过程的影响，因而许多学者用石墨拉曼谱峰特征来计算其经历的峰期变质温度(Beyssacetal., 2002, 2019; Busemannetal., 2007; Aoyaetal., 2010; Hilchie and Jamieson, 2014; 李凯月等, 2018)。

图6 辽河群不同含石墨岩石样品中石墨拉曼光谱学特征(a、b) 含石墨透闪大理岩(ZK3-7H4)；(c)含石墨黑云母长英质片麻岩(ZK7-4H31)；(d)含石墨变质杂砂岩(ZK7-4H55)Fig.6 Raman spectra of graphite in different graphite-bearing rock samples from the Liaohe Group(a, b) graphite-bearing tremolite marble (ZK3-7H4); (c) graphite-bearing biotite feldspar gneisses (ZK7-4H31); (d) graphite-bearing meta-greywacke (ZK7-4H55)

Beyssacetal. (2002)通过对比分析大量的有机物变质成因石墨的拉曼光谱学特征，发现可以通过峰面积比R2 [R2=D1/(G+D1+D2)]和峰强度比R1(R1=D1/G)两个参数对石墨的结晶度进行定量化分析，而且可以根据R2比值计算石墨结晶温度，其计算方程为T(℃) =-445×R2+641(误差范围±50℃，温度范围330～650℃)。我们对辽河群石墨矿中不同含石墨岩石样品中的石墨进行了详细的拉曼光谱学分析，通过谱峰分解计算得出这些样品的R2比值位于0.03～0.2 区间内，对应的石墨形成温度为551～627℃。贺高品和叶慧文(1998a, b)采用传统矿物对温压计对南、北辽河群的P-T条件进行了估算，得出其峰期变质温度压力条件分别为610～670℃/0.60～0.68GPa和600～640℃/0.64～0.73GPa，和我们计算的石墨拉曼光谱峰记录的峰期温度相似。然而，最近Liuetal. (2019)在南辽河群的三家子地区发现了含堇青石的泥质麻粒岩，相平衡模拟表明其峰期变质温度可以达到790～840℃。尽管到目前为止在北辽河群还没有发现麻粒岩相的变质记录，本文研究区的石墨矿是否和华北克拉通孔兹岩带中的石墨矿一样经历了麻粒岩相的变质？我们通过石墨拉曼光谱学研究认为这种可能性不大。原因如下，石墨拉曼光谱峰主要受温度控制，其拉曼光谱参数对温度很敏感，如果含石墨的岩石样品达到了麻粒岩相的变质作用，石墨的拉曼光谱峰将记录这一变质温度条件，石墨拉曼光谱峰一级谱峰区D1、D2缺陷峰将消失，对应的R2比值为零(Chopinetal., 1991; Beyssacetal., 2002)。而辽河群石墨矿中没有发现D1和D2缺陷峰消失的石墨拉曼峰，其最高结晶度对应的温度大约为627℃，记录了高角闪岩相的变质作用。

4.2 石墨成因分析：有机成因VS.无机成因

石墨是地球表面稳定存在的单质碳，是地表碳储库的重要组成部分(Mackenzieetal., 2004)。岩石中的石墨主要有两种形成机制：(1)有机物经过变质作用转变成石墨；(2)在一定的温度压力条件下从碳过饱和的C-H-O流体中沉淀结晶形成石墨(Luqueetal., 1998, 2009；Zhuetal., 2020)。

稳定碳同位素地球化学分析可以示踪碳的来源，常用来研究金刚石的成因(Thomassotetal., 2007; Walteretal., 2011)以及早期(太古宙)生命演化(Schidlowski, 2001; Uenoetal., 2002; Van Zuilenetal., 2003; Papineauetal., 2010; Leplandetal., 2011)，也是研究石墨成因的有效手段(Luqueetal., 2012; Yangetal., 2014; Zhongetal., 2019)。由于碳同位素的分馏作用，原始地球的碳储库发生碳同位素分馏，轻碳同位素(12C)倾向于富集在有机物中，碳酸盐则富集重碳同位素(13C)(Javoyetal., 1986; Chackoetal., 2001; Luqueetal., 2012)，因此会形成具有不同碳同位素值的碳储库。形成石墨的碳主要有三个来源：有机物、碳酸盐和地幔碳。其中有机物的碳同位素(δ13CPDB‰)多位于-17‰～-40‰范围内，平均为-26‰～-28‰(Schidlowski, 1987, 2001; Hoefs, 2009)。典型海洋碳酸盐往往富集重碳同位素，其δ13CPDB值位于-2‰～+4‰区间内(Sharp, 2007)。地幔碳同位素的值则位于-7‰附近(Hahn-Weinheimer and Hirner, 1981; Weisetal., 1981)。

图7 辽河群含石墨岩石样品石墨碳同位素(a)岩芯钻孔ZK3-7；(b)岩芯钻孔ZK7-4；纵坐标左侧为岩芯编号，右侧为取样深度Fig.7 Carbon isotope of graphite in different graphite-bearing rock samples from the Liaohe Group(a) Borehole ZK3-7; (b) Borehole ZK7-4; The left side of the ordinate is the sample number, and the right side is the sampling depth

由于地球有机物相对无机碳(碳酸盐和地幔碳)富集轻同位素，许多学者认为变沉积岩中富集12C碳同位素的石墨是有机物在高级变质作用过程中变质形成的(Landis, 1971; Grew, 1974; Santosh and Wada, 1992; Radhikaetal., 1995; Dissanayakeetal., 2000)。前人对印度南部Kerala孔兹岩带中的石墨进行了详细的碳同位素研究，发现呈浸染状分布在层状变泥质岩中的石墨具有轻的碳同位素组成(δ13CPDB=-34.3‰～-17.5‰; 图8)，认为这些石墨是有机物变质形成的，具有有机碳的同位素特征(Santosh and Wada, 1992, 1993a, b; Radhikaetal., 1995; Radhika and Santosh, 1996)。然而Kerala孔兹岩带剪切带以及伟晶花岗岩中的石墨都记录了更重的碳同位素特征，分别为-8.2‰～-12.4‰和 -10‰～-15.1‰(图8；Radhika and Santosh, 1996)，被认为是从含CO2流体中沉淀结晶出来的。华北克拉通西北部的孔兹岩带中的石墨矿具有和印度南部Kerala孔兹岩带石墨相似的碳同位素特征，在变质沉积岩中石墨具有轻的碳同位素值(δ13CPDB=-25.3‰～-25.7‰)，而在石英脉体(δ13CPDB=-19.1‰～-20.9‰)和长英质浅色体(δ13CPDB=-15.8‰～-16.8‰)中的石墨则分别具有更重的碳同位素值(图8和图9d；Yangetal., 2014)，Yangetal. (2014)认为变质沉积岩中的石墨是由有机物变质形成的，而对于石墨脉体以及长英质浅色体中的石墨，则是多种碳源汇聚的结果，是从有机物脱挥发分形成的含碳流体和外来的含CO2流体混合形成的流体中沉淀出来的。陈衍景等(2000)总结了华北克拉通周缘不同地区的石墨矿碳同位素特征(图8和图9a-c)，包括佳木斯(δ13CPDB=-17.0‰～-26.5‰)、内蒙古中部(δ13CPDB=-11.4‰～-28.9‰)和华北克拉通东南缘(δ13CPDB=-14.7‰～-26.8‰)，它们的碳同位素值都具有很宽的变化范围，陈衍景等(2000)提出混合岩化过程中无机碳的加入是导致变沉积岩中石墨碳同位素变重的主要原因。从以上讨论我们可以看出不管是印度Kerala孔兹岩带中的石墨，还是华北克拉通周缘的石墨，其碳同位素值都具有很宽的变化范围，具有复杂的成因过程，容易造成石墨成因的多解性。图10总结了不同碳源形成石墨的碳同位素特征，如果石墨的碳来源于富12C的有机物，则形成的石墨也具有轻的碳同位素特征；如果石墨的碳来源于富13C的碳(碳酸盐或地幔碳)，则形成的石墨具有重的碳同位素特征；此外，如果富12C以及富13C的两种流体混合，则会形成具有中间碳同位素值的石墨(Crespoetal., 2004; Luqueetal., 2012)。

辽河群石墨矿的碳同位素值具有很宽的变化范围，从 -16.49‰到-25.93‰，和华北克拉通周缘其它地区的石墨矿相似。在图8中，这些石墨矿岩石的石墨碳同位素值主要位于生物成因区域内，然而单独依靠石墨碳同位素并不能很好的限定石墨是通过有机物变质形成还是从含碳流体中沉淀形成的。因为流体沉淀形成石墨的碳同位素值也可以位于生物成因区域内(富集轻碳同位素)，例如英国Borrowdala(Barrenecheaetal., 2009; Luqueetal., 2009; Ortegaetal., 2010)、西班牙Huelma(Barrenecheaetal., 1997)、美国New Hampshire(Rumble III and Hoering, 1986; Rumble IIIetal., 1986)和Black Hills(Nabeleketal., 2003; Huff and Nabelek, 2007)地区的石墨矿都是通过流体沉淀形成的，然而它们的碳同位素值大多位于生物成因区域内(图8)。因此，对于辽河群石墨矿的成因机制及碳源问题，需要进一步结合岩相学特征和沉积构造环境来限定。

图8 辽河群石墨矿和其它地区的石墨矿石墨碳同位素特征对比数据来源：华北克拉通周缘其它地区(陈衍景等，2000；Yang et al., 2014); 印度南部Kerala孔兹岩带(Santosh and Wada, 1992, 1993a, b; Radhika et al., 1995; Radhika and Santosh, 1996); 英国Borrowdala(Barrenechea et al., 2009; Luque et al., 2009; Ortega et al., 2010); 西班牙Huelma(Barrenechea et al., 1997); 美国New Hampshire(Rumble III and Hoering, 1986; Rumble III et al., 1986); Black Hills地区(Nabelek et al., 2003; Huff and Nabelek, 2007)Fig.8 Carbon isotope composition of graphites from the graphite deposit of the Liaohe Group and other regionsData sources: the North China Craton (Chen et al., 2000; Yang et al., 2014); Kerala Khondalite Belt, Southern India (Santosh and Wada, 1992, 1993a, b; Radhika et al., 1995; Radhika and Santosh, 1996); Borrowdala (UK) (Barrenechea et al., 2009; Luque et al., 2009; Ortega et al., 2010); Huelma (Spain) (Barrenechea et al., 1997); New Hampshire (USA) (Rumble III and Hoering, 1986; Rumble III et al., 1986); and Black Hills (USA) (Nabelek et al., 2003; Huff and Nabelek, 2007)

辽河群含石墨岩石的微量元素数据表明，这些样品具有Th、U元素正异常和Nb、Ta、Sr、Ti元素负异常的特征，此外，球粒陨石标准化稀土元素图解表明它们具有轻稀土富集和Eu负异常的特征，这些地球化学特征和PAAS以及上地壳的微量元素特征相似。Th/Sc-Zr/Sc图解可以用来很好的指示沉积物的再循环过程(McLennanetal., 1993; Asieduetal., 2000)，沉积物再循环过程中Zr/Sc的比值相对Th/Sc会发生显著的增加。在Th/Sc-Zr/Sc图解中(图4b)，辽河群含石墨样品的Th/Sc和Zr/Sc比值显示一定的正相关性，落在源区组分变化线上，没有出现Zr的显著增加现象，表明这些样品没有经历后期沉积物的再循环过程。Lietal. (2015)对辽河群变沉积岩进行了详细的地球化学研究，根据化学蚀变指数(CIA; Nesbitt and Young, 1982)以及成分变化指数(ICV; Coxetal., 1995)，表明辽河群变沉积岩具有低的成熟度，未经历长途搬运，是快速堆积的产物。前人通过碎屑锆石和地球化学研究，认为辽河群变沉积岩(辽河群石墨矿的原岩)主要形成于弧后盆地活动大陆边缘环境(Lietal., 2015; Wangetal., 2017a, 2020a)。岩相学研究表明辽河群石墨矿的石墨主要以浸染状分布在辽河群变沉积岩中，因此我们认为在辽河群沉积物快速堆积过程中，还伴随着大量有机物的加入。沉积物快速堆积且未经历长途搬运，为有机物的埋藏保存提供了良好的环境。在后期的变质作用过程中，有机物发生变质分解并逐渐形成具有高结晶度的石墨。岩相学研究表明含石墨岩石样品中普遍含有一定量的黄铁矿，表明其沉积环境具有一定的还原性。有机物在还原环境下分解的产物主要以CH4为主，而CH4富集轻的碳同位素(12C)，因此会使残余的有机碳同位素逐渐变重(Rumble III and Hoering, 1986; Wadaetal., 1994)。这可以用来解释辽河群石墨矿具有比较宽泛的碳同位素特征，并且具有逐渐变重的趋势。然而，有研究者用多组分流体的瑞利分馏来解释具有中间碳同位素特征的石墨成因(图11; Ray, 2009; Luqueetal., 2012; Yangetal., 2014)，这种模型可以很好的解释华北克拉通孔兹岩带中石英脉以及浅色体中具有较重碳同位素特征的石墨成因(Yangetal., 2014)。但是，岩相学研究表明辽河群石墨矿大多以浸染状分布在层状地层中，没有流体沉淀结晶石墨的特征，因此不太可能是从多组分流体中通过瑞利分馏形成的。因此，我们认为辽河群石墨矿是由有机物经过变质作用形成的，其变质程度可以达到高角闪岩相，低于华北克拉通孔兹岩带中石墨的变质程度(麻粒岩相；Yangetal., 2014)；有机物在变质作用过程中发生不同程度的CH4分解作用，导致石墨碳同位素具有向富13C同位素演化的趋势。

图9 华北克拉通周缘石墨矿石墨碳同位素直方图(a-c)佳木斯、内蒙古、华北克拉通东南缘(陈衍景等，2000);(d)华北克拉通西北缘(Yang et al., 2014)Fig.9 Histograms assembling the carbon isotope composition of graphites from Jiamusi Block, Inner Mongolia, southeastern margin of North China Craton (a-c, after Chen et al., 2000) and northwestern margin of North China Craton (d, after Yang et al., 2014)

图10 不同碳源形成不同碳同位素值石墨示意图(据Crespo et al., 2004修改)Fig.10 Schematic diagram depicting the carbon isotope ratio in graphite derived from different carbon sources (modified after Crespo et al., 2004)

图11 在600℃温度下，单组分(CO2，黄线)和多组分(CO2+CH4，蓝线)含碳流体瑞利分馏形成石墨碳同位素演化曲线(据Luque et al., 2012修改)Fig.11 Evolution curves of δ13C values for graphite precipitated from single component (CO2, yellow lines) and multicomponent (CO2+CH4, blue lines) during Rayleigh fractionation at 600℃ (modified after Luque et al., 2012)

4.3 辽河群石墨矿形成时代

前人对胶-辽-吉带内辽河群原岩时代开展了大量的同位素年代学研究(辽宁省地质矿产局, 1989; Luoetal., 2004, 2008; Luetal., 2006; 李壮等, 2015; 刘福来等, 2015; 王惠初等, 2015; Wangetal., 2017a, 2020a; 刘平华等, 2017)。早期基于K-Ar、Rb-Sr和Sm-Nd全岩等时线或单颗粒锆石蒸发定年研究，前人认为辽河群形成于2.3～1.9Ga(张秋生等, 1988; 辽宁省地质矿产局, 1989; 白瑾, 1993)，然而由于胶-辽-吉带辽河群在后期经历了多期变质-深熔事件叠加，会引起K-Ar、Rb-Sr和Sm-Nd同位素体系的开放，导致上述研究方法不能很好的制约辽河群的沉积时代。随后，许多学者使用LA-ICP-MS锆石定年方法对辽河群碎屑锆石进行了详细的研究，以其中岩浆锆石的最小年龄代表沉积的最老年龄，而最大变质年龄则代表沉积的最新年龄，通过这一方法得出辽河群的沉积时代为2.05～1.95Ga (Luoetal., 2004, 2008; 刘福来等, 2015)。近年来，Wangetal. (2020a)对北辽河群进行了大量碎屑锆石年代学的研究工作，年龄数据表明北辽河群里尔峪组的年龄谱峰为2.17Ga，高家峪组和大石桥组都具有双年龄谱峰，分别为2.51Ga/2.17Ga和2.51Ga/2.18Ga。此外，在北辽河群发育大量切穿沉积地层的变质基性岩脉(Mengetal., 2014; Xuetal., 2018)，年代学研究表明其侵位时间为2.13Ga。因此，辽河群的沉积时代应早于2.13Ga。综上我们认为位于北辽河群高家峪组中的石墨矿沉积时代为2.13～2.17Ga。根据前述讨论(章节4.2部分)，我们认为在2.13～2.17Ga，有大量的有机物沉积在弧后盆地活动大陆边缘环境，在后期的弧-陆碰撞造山以及后碰撞作用过程中(2000～1895Ma; 1875～1850Ma; Xu and Liu, 2019)，有机物逐渐发生变质形成石墨并聚集形成石墨矿床。

4.4 对地球古环境的启示

年代学研究表明，辽河群石墨矿沉积时代大约在2.13～2.17Ga，和华北克拉通西北缘孔兹岩带石墨矿的沉积时代(2.07～2.15Ga; Lietal., 2011; Zhangetal., 2014)相似。陈衍景等(2000)总结了华北克拉通周缘石墨矿特征，认为这些含石墨岩石的沉积时代可能为2.05～2.3Ga。这些年代学数据表明华北克拉通周缘巨量石墨矿床大部分都形成于古元古代层侵纪(2.3～2.05Ga)。与此同时，在世界各克拉通层侵纪碳酸盐地层序列中普遍发现了碳同位素(δ13Ccarb)正异常现象(Baker and Fallick, 1989a, b; Karhu, 1993; Karhu and Holland, 1996; Melezhik and Fallick, 1996; Aharon, 2005; Bekkeretal., 2006; Frauensteinetal., 2009; Maheshwarietal., 2010；Tangetal., 2011; Pufahl and Hiatt, 2012; Martinetal., 2013)。这一全球性δ13Ccarb正漂移事件又称为拉玛岗地-瓦图里事件(Lomagundi-Jatuli Event, LJE; Baker and Fallick, 1989a,b)。有许多研究者认为LJE事件和有机物的埋藏有关(Karhu and Holland, 1996; Bekkeretal., 2008)。根据公式δ13Cin=δ13Ccarb(1-forg)+δ13Corgforg(δ13Cin、δ13Ccarb、δ13Corg分别代表进入大气中的碳同位素值、沉积碳酸盐碳同位素值和有机物碳同位素值，forg则代表有机物埋藏的比例)，在δ13Cin和δ13Corg值不变的情况下，如果有机物埋藏量增多，相应的碳酸盐矿物δ13Ccarb会升高(Eguchietal., 2020)。此外，根据光合作用CO2+H2O=CH2O+O2，有机物的埋藏量控制着地表氧含量的变化，有机物埋藏量的增加会导致大气氧含量逐渐升高(Karhu and Holland, 1996)。大量研究表明地球表生环境首次大规模充氧事件(大氧化事件，GOE)发生时间大约在2.3～2.0Ga(Bekkeretal., 2006; Lyonsetal., 2014; Luoetal., 2016)，和LJE事件的时间相一致。在这一时期，地球各圈层性质发生了全球性突变，例如发育全球性硅铁建造，红层、沉积磷矿和蒸发岩普遍发育，叠层石及碳酸盐岩大量出现，发育大量石墨矿床等(陈衍景, 1990； Tangetal., 2016)。华北克拉通周缘古元古代(2.3～2.05Ga)大量石墨矿床的沉积很有可能是GOE和LJE事件的沉积响应。

关于地球古元古代早期大氧化事件的成因，前人已经做了大量的研究(Kastingetal., 1993; Ohmoto, 1996; Holland, 1999, 2002, 2006, 2009; Ohmotoetal., 2004, 2006; Canfield, 2005; Kump and Barley, 2007; Kump, 2008; Lyonsetal., 2014; Leeetal., 2016; Luoetal., 2016; Duncan and Dasgupta, 2017; Eguchietal., 2020)。传统简单模型认为这时期蓝藻细菌的出现可以导致大气氧含量的快速增加，然而Brocksetal. (1999)在2.7～2.8Ga沉积地层中也发现了蓝藻细菌，表明GOE事件不能简单的用蓝藻细菌的出现来解释。Kastingetal. (1993)认为大气H2的逃逸可以改变地球的氧化还原状态，使地幔逐渐氧化，进而使火山喷出氧化性气体含量增加。Leeetal. (2016)认为2.7～2.5Ga大陆壳形成后，地壳岩石由铁镁质向硅酸质的转变使地壳吸附氧的能力降低，可以促使大气氧含量的增加；大陆壳形成后风化作用可以为海洋生物提供大量的营养物质，促使光合作用产氧量的增加。此外，2.5Ga左右板块俯冲的初始启动可以使大量有机物被俯冲带入地幔深部并固存下来，进而促使大气氧含量的增加(Duncan and Dasgupta, 2017)。最近，Eguchietal. (2020)提出，2.3～2.0Ga期间地球构造转变导致火山喷出CO2含量的大量增加，可以使地表沉积的碳酸盐和有机物含量都显著增加，进而导致大气氧含量的升高。综合来看，不管是哪种模式，大气氧含量的升高都伴随着大量有机物的沉积。华北克拉通周缘古元古代早期(2.3～2.05Ga)沉积的大量有机物可能记录了这一时期大气充氧过程，这些有机物在后期经历了多期变质事件的改造而转变成了具有高结晶度的石墨，进一步对其进行沉积构造环境方面的研究对于揭示大氧化事件的成因仍具有重要的意义。

5 结论

根据对胶-辽-吉造山带辽河群石墨矿(甜水乡马沟石墨矿)详细的岩相学、地球化学、拉曼光谱学以及碳同位素等方面的研究，本文得到如下几点认识：

(1)辽河群石墨矿主要岩石类型为含石墨变质杂砂岩、含石墨黑云母长英质片麻岩和含石墨透闪大理岩，地球化学特征表明它们具有低的成熟度，未经历长途搬运，是快速堆积的产物，主要沉积于活动大陆边缘环境。

(2)辽河群石墨矿中的石墨具有非常高的结晶度，其拉曼峰面积比值R2=0.03～0.2，对应的石墨形成温度为551～627℃，记录了高角闪岩相的变质作用。

(3)辽河群石墨矿碳同位素值具有很宽的变化范围(δ13CPDB=-16.49‰～-25.93‰)，是有机物在变质过程中脱CH4造成的。其沉积时代为2.13～2.17Ga，在后期的弧-陆碰撞造山以及后碰撞作用过程中(2000～1895Ma; 1875～1850Ma)，有机物逐渐发生变质形成石墨并聚集形成石墨矿床。

(4)华北克拉通周缘大量石墨矿的沉积时代(2.3～2.05Ga)和大氧化事件发生的时间一致，可能是地球早期大氧化事件的沉积响应。

致谢感谢辽宁省有色地质一〇三队在野外工作中提供的大力支持和帮助；感谢中国地质科学院地质研究所显微激光拉曼实验室张聪副研究员在石墨拉曼光谱测试过程中提供的帮助；感谢中国科学院地质与地球物理研究所稳定同位素分析实验室李洪伟老师在石墨碳同位素测试过程中提供的帮助；感谢审稿专家汤好书研究员和李旭平教授对本文提出了宝贵的修改意见。