应用于流体包裹体CO2碳同位素组成的拉曼光谱定量研究探讨
药瑛, 孙樯*
造山带与地壳演化教育部重点实验室,北京大学地球与空间科学学院,北京 100871
*通信作者:孙樯(1970-),男,山东泰安人,副教授,主要从事地球化学研究.E-mail:qiangsun@pku.edu.cn

作者简介:药瑛(1990-),女,山西长治人,硕士研究生,主要从事实验地球化学研究.E-mail:yaoying@pku.edu.cn

摘要

CO2作为一种流体包裹体中常见的组分,对其中碳同位素的定量研究,能够提供很多地质历史时期地质流体源区的地球化学信息,并为地球碳循环的研究奠定基础。传统的质谱测试存在无法进行单个流体包裹体测试的缺点。拉曼光谱因其快速、无损、高精度的特点被广泛应用于地球科学各个研究领域,尤其在单一流体包裹体的研究中更为重要。当前,拉曼光谱在定性研究物质结构和成分方面已经十分成熟,相比之下,拉曼光谱定量研究仍显不足。结合前人的研究和实际的实验,根据对拉曼光谱基本原理的分析,提出激光拉曼光谱定量研究应当基于相对散射强度进行,且利用13CO212CO2的拉曼谱峰强度比(I+13/I+12)表征流体包裹体中CO2的碳同位素组成是可行的。而影响该方法应用的是其较低的准确度,解决这一问题的关键在于如何在单次实验中同时获得准确而强度较高的I+13I+12的值。

关键词: 拉曼光谱; 定量; CO2 碳同位素
中图分类号:P597 文献标志码:A 文章编号:1001-8166(2016)10-1032-09
Raman Quantitative Measurements for Carbon Isotopic Composition in CO2-Rich Fluid Inclusion: A Preliminary Study
Yao Ying, Sun Qiang*
Key Laboratory of Orogenic Belt and Crustal Evolution (Ministry of Education), School of Earth and Space Sciences, Peking University, Beijing 100871, China

First author:Yao Ying(1990-),female,Changzhi City,Shanxi Province,Master student. Research areas include experimental geochemistry.E-mail:yaoying@pku.edu.cn

*Corresponding author:Sun Qiang(1970-),male,Taian City,Shandong Province,Associate Professor. Research areas include geochemistry.E-mail:qiangsun@pku.edu.cn

Abstract

Carbon dioxide is a common and important component in fluid inclusions. Because carbon isotopic ratio of CO2in fluid inclusion can provide geochemical characteristics of source rocks, many works have been conducted to measure the carbon isotopic ratio (13C/12C) of CO2. In general, carbon isotopic composition (13C/12C) of CO2 fluid inclusions are measured by mass spectroscopy. However, mass spectroscopy is a destructive analytical technique, and can not be applied to measure single inclusion. Raman spectroscopy is an efficient non-contact and non-destructive method, and has been widely employed in many research fields. In fact, due to the difference of carbon atom mass, the Femi resonance of13CO2 is lower than that of12CO2, so they can be identified in Raman spectrum. In principle, the carbon isotopic composition (13CO2/12CO2) in CO2 is closely related to the Raman intensity ratio between13CO2 and12CO2 (I+13/I+12). Therefore, Raman spectroscopy can theoretically be utilized to measure the carbon isotopic composition (13CO2/12CO2) of CO2. However, this method is seriously hampered by the poor measurement accuracy and precision. In this study, based on the theoretical analysis, and in combination with previous Raman studies on carbon isotopic composition in CO2, a preliminary research on Raman quantitative measurement for carbon isotopic composition in CO2 was carried out.

Keyword: Raman spectroscopy; Quantitative analysis; Carbon isotope.
1 引言

地质流体在地质成岩及成矿过程中发挥着重要作用, 且以流体包裹体形式保存下来。对流体包裹体组成的研究, 是正确认识地质流体演化及流体与成矿作用关系的前提。拉曼光谱因具有与测试样品非接触、非破坏性、样品制备简单、操作方便、精度高等技术特点, 已经被普遍应用于流体包裹体的研究。流体包裹体以H2O-NaCl-CO2体系为代表[1, 2], 对这类流体包裹体的研究是研究更复杂流体包裹体的基础。

当前拉曼光谱主要应用于定性研究, 即根据振动频率进行分子类型的确定, 而定量研究明显不足。加强拉曼光谱定量化方面的研究是进一步发展激光拉曼技术所必需的。本文结合前人的研究方法和成果对应用拉曼光谱定量研究H2O-NaCl-CO2体系二氧化碳(CO2)的碳同位素组成进行了探讨。

2 碳同位素研究的地质意义

碳与地球的形成和发展息息相关, 绝大部分(> 90%)的碳被储存在地球深部, 却又与地表碳有着千丝万缕的联系[3, 4]。研究地球内部的碳循环对地球演化史、矿床成因、全球气候变化[5]等问题的研究具有重要意义。在碳循环的研究中, 碳同位素组成是一个重要指标, 通常用δ 13CPDB来表征。地球碳循环的研究对象有单质碳和化合态碳, 其中金刚石作为典型的地球深部来源的单质碳被认为是研究深部碳循环的重要样品[6], 而化合态碳中CO2则是重要的研究对象之一, 其他还包括碳酸盐矿物、CH4和SiC等[7]

CO2广泛存在于自然界中, 参与并影响了许多地质作用过程。地质作用过程中, CO2可以以气相、液相的形式被流体包裹体保存下来, 因此流体包裹体中CO2碳同位素组成的研究, 在研究地壳、地幔循环过程以及矿床、油气的形成等许多问题上都具有重要意义。

2.1 地幔流体时空分布及演化研究

地幔流体的研究一直被认为是研究地球内部物质组成和演化的重要手段。地幔是地球碳的主要源区之一, 碳循环受地幔碳同位素组成和不均一分布的影响。结合来自地幔不同深度的流体包裹体样品中碳同位素的组成特征与其他的同位素信息能够很好地研究地幔流体的形成和演化。而且统计全球范围内地幔流体的同位素组成也有利于地幔物质的全球性对比[8]

自1965年Roedder[9]在碱性玄武岩的超镁铁质捕掳体中发现大量CO2包裹体开始, 掀起了通过流体包裹体研究地幔流体的浪潮。研究证明, 富CO2的流体包裹体在地幔岩中普遍存在, 上地幔流体、下地壳麻粒岩化均具有富纯CO2流体的特征[10]。科学家们结合对不同地区的地幔岩石中的包裹体、金刚石、金伯利岩、碳酸盐、大洋玄武岩的碳同位素研究, 发现地幔流体中的δ 13C分布范围广, 但是大部分样品的δ 13C值集中分布在-5‰ 附近和-15‰ ~-25‰ [11]。这为研究地幔碳的多样性和原始地幔碳是否均一奠定了基础, 也为研究地球深部碳循环提供了手段。

2.2 成矿流体来源研究

CO2流体包裹体与造山型金矿成矿作用有密切联系[12, 13]。研究表明脉状金矿中常发育较多的富CO2包裹体(CO2-H2O体系); 造山带金矿流体具有低盐度、富CO2的特点, 成矿作用与CO2-H2O流体具有密切联系[14] , CO2在金矿形成过程中能够缓冲流体pH值, 进而促进Au富集成矿[15]。Klein等[16]对Carará 金矿含金石英脉中的富CO2包裹体进行研究, 结合δ 13CPDB和流体包裹体的密度数据, 推断出成矿流体中的CO2来源于紫苏花岗岩岩浆和同时代麻粒岩相变质作用。卢焕章等[14]、李吉君等[17]和王勇等[18]通过对矿床中CO2流体包裹体中δ 13C值的测定, 结合前人对地幔CO2流体的研究, 在推断成矿流体来源, 反演地质历史时期元素的运移, 探讨成矿原理方面获得重要成果。

Mumm等[19]根据加纳Ashanti金矿中极富CO2流体的成矿作用, 认为纯CO2或CO2含量远大于H2O的流体是一种新的成矿流体类型。虽然Mumm等的观点引起了争论, 但也不能否认对成矿流体中CO2-H2O体系的研究很大程度上完善了我们对成矿作用的认识。

2.3 油气藏成因研究

油气藏中常含富CO2的流体包裹体, CO2存在与否能指示石油形成所处的阶段; 碳同位素的不同组成可以指示油气藏成因。石油形成会经历不同的阶段, 即未成熟、低成熟、高成熟、过成熟阶段, 在不同阶段包裹体的成分类型不尽相同。未成熟阶段包裹体成分以H2O和CO2为主, 到高成熟阶段CO2消失[20]。戴金星等[21]总结在气藏中CO2含量≥ 60%时, CO2为无机成因; 含量< 15%时则为有机成因。δ 13 CCO2> -8‰ 时, CO2为无机成因, 当δ 13 CCO2< -10‰ 为有机成因。因此, δ 13 CCO2的值可以指示气藏的成因以及在成藏过程中有机因素和无机因素的影响过程。

3 流体包裹体中CO2碳同位素组成研究方法

对流体包裹体中CO2碳同位素组成进行研究具有重要意义, 因此对其分析测试一直为研究者所关注, 并开展了大量的研究工作。流体包裹体中CO2碳同位素组成分析测试方法可大致分为破坏性分析测试方法和无损分析测试方法。

3.1 破坏性分析测试方法

所谓破坏性分析测试方法是指通过破坏流体包裹体提取其中的待检测物质, 然后利用质谱仪进行包裹体内同位素含量以及比值的测定。针对群体流体包裹体进行分析测试的方法主要有机械破碎法、热爆法和分阶段加热法。近年来还发展了可进行单个流体包裹体的破坏性分析测试— — LA-ICP-MS方法。

机械破碎法是将样品去气后利用如真空球磨等机械方法粉碎样品, 以释放包裹体中的气体。热爆法和分阶段加热法都是利用石英管式加热炉加热流体包裹体使之爆裂, 在真空条件下收集释放的CO2。李洪伟等[22]对传统的热爆法进行了改进, 利用单个密封石英试管装样, 马弗炉批量加热以使包裹体爆裂得到气体, 既避免了传统方法样品的相互污染同时也提高了实验效率。分阶段加热法较传统热爆法更加优越, 因为对不同温度阶段收集的气体进行质谱测试, 能区分不同来源的气体, 也避免保存在矿物粒间的气体对包裹体成分测试的干扰[23~25]。方解石中的含CO2包裹体由于寄主矿物加热可分解产生CO2而不适合加热, 可利用CaCO3-CO2平衡分馏方程来获得包裹体中CO2碳同位素组成[26]

在进行流体包裹体CO2碳同位素破坏性分析测试研究中, 对群体包裹体进行的测试混入多期次的包裹体, 所测得的结果只能代表地质流体中同位素组成的平均值, 因而掩盖了不同期次的地质活动信息。此外, 破坏性分析测试不可避免地会受到测试时大气环境的影响, 这都可能导致实验误差。因此, 对单一包裹体进行非破坏性分析测试, 日益受到重视。

近年来, LA-ICP-MS被应用于对单个流体包裹体进行破坏性分析测试。该方法利用激光剥蚀来提取单个包裹体中的待测物质, 然后进行质谱分析, 获得的是离子含量比值。尽管LA-ICP-MS方法可以进行单个流体包裹体的原位分析, 但也具有明显的缺点:①激光剥蚀过程可能造成包裹体中流体成分的改变而产生误差, 其中包括激光的高能量溶蚀宿主矿物, 从而污染包裹体中的流体, 且激光导致宿主矿物破裂造成部分流体泄露[27]; ②该方法对于CO2的检测限为10-9mol, 而一般典型的CO2包裹体中CO2含量只有10-11mol, 因此其较高的检测限无法用于测试低含量同位素的组成; ③应用LA-ICP-MS对单一包裹体进行的测试, 需要破坏样品, 因此不能进行重复测试和数据验证。

3.2 无损分析测试方法

当前, 流体包裹体的无损分析测试主要应用红外光谱和激光拉曼光谱。理论上, 双原子分子简谐振动可表示为:

V= 12πkμ

式中:V为分子振动频率, k为键力常数, μ 为折合质量。可以看出, 对于分子振动具拉曼或红外活性的分子, 其不同同位素组成由于质量数的差异而对应不同的谱峰位置。因此, 分子光谱可以定性判断同位素的组成, 这也为定量分析测试提供了基础。

红外光谱测试可以用于检测单个流体包裹体的成分。李相贤等[28]利用傅里叶红外光谱法测试标准CO2(南京特种气配厂制), 较同位素质谱法(Isotope Ratio Mass Spectrometry, IRMS)所得数据, 平均差异仅为0.25%。这说明红外光谱法可以应用于单个流体包裹体中CO2的碳同位素组成的测定。但是, 通常情况下可用于红外光谱测试的流体包裹体, 其大小不低于25 μ m, 而许多天然流体包裹体的尺寸要小于25 μ m, 难以进行红外光谱测试。虽然同步辐射傅里叶转换红外光谱(Synchrotron Fourier Transform Infrared, SFTIR)的测试光斑可以缩小到3 μ m[29], 但是, 实际应用中受同步辐射的实验设备开放时间限制。

相比之下, 激光拉曼光谱的工作波长在可见光范围内, 其光斑可小至1 μ m, 且样品前期处理简单快捷, 因而成为进行单一流体包裹体无损分析测试的理想方法。然而其在单个流体包裹体中CO2的碳同位素组成的测定方面的研究和应用却非常缺乏。

4 激光拉曼光谱的定量研究

1975年, Rosasco[30]利用拉曼光谱对天然流体包裹体进行了研究, 提出拉曼光谱在定性和定量研究方面应用的可能性。之后, 有越来越多的学者在该方面开展了工作, 显示了拉曼光谱在该领域的研究优势。近年来, 激光拉曼光谱已被广泛应用于流体包裹体的定性和定量研究中。对单个流体包裹体的定量研究集中在对其气液相组分、同位素化学组成分析及包裹体内的气体密度和压力、有机热成熟度等参数的测定[31]

当前, 拉曼光谱定量研究过程中, 根据参考体系的类型可大致分为外标法和内标法。所谓内标法, 是指参考系与被测试对象在测试过程中, 拉曼光谱被单次测试同时记录下来, 从而保证了测试条件的一致。外标法则是在保证测试条件一致的条件下, 参考系与被测试对象先后测试完成。可以看出, 相比于外标法, 内标法更加准确快捷。

流体包裹体成分定量研究方面, 具有拉曼活性的分子, 如:CO2, CH4, N2等利用内标为参考系进行物质含量的定量研究。众多溶于水却在常温下无拉曼光谱活性的阴离子(Cl-, C O32-, S O42-等)以及阳离子(Na+, K+, Ca2+, Mg2+等), 通过以水的拉曼光谱图中O-H振动为参考系被定量研究。研究表明, 在含单一阴离子的溶液中, 离子浓度和水的O-H振动峰强度具良好的线性关系[32]; 常见的氯盐溶液中阳离子浓度和水O-H振动去卷积后的2个峰强比(低频/高频)之间也具有良好的线性关系[33]; 据此可以测定流体包裹体液相中各离子的浓度。Wang 等[34]进一步限定了孙樯等[35]提出的利用Cl-浓度与Δ K(定义的水拉曼峰比值相关参数)之间的线性关系定量Cl-的方法的适用范围。

包裹体密度定量研究方面, 李淑玲[36]指出在酸性条件下制得的CO2气体, 其拉曼光谱强度与浓度在一定压力范围内具线性关系。另外, Rosso等[37]研究表明在均一气相或液相CO2中, CO2费米共振峰的间距X与CO2的密度间的线性关系在一定密度范围内(≤ 1.2 g/cm3)可以用来定量CO2密度。

包裹体盐度的定量分析中, Mernagh等[38]提出的基于拉曼光谱的“ 频移参数法” 被广泛应用, 该方法通过水拉曼峰的频移参数与包裹体盐度之间的线性关系来计算包裹体的盐度[39]

包裹体压力定量研究中, 李晶等[40]在金刚石压腔(DAC)的实验条件下, 定量研究得到了NaCl-CaCl2-H2O系统中水Δ V(O-H)(实验条件下O-H振动峰较常温常压下振动峰的偏移量)与温度、压力、盐度之间的定量关系式。流体包裹体中同位素的定量分析中, C, H, N, O同位素都能够用拉曼光谱来分析[41~43], 其中对CO2中碳同位素组成的定量研究较为缺乏。因此对拉曼光谱在CO2中碳同位素组成方面的定量研究具有重要意义。

5 拉曼光谱CO2碳同位素组成的研究
5.1 拉曼散射定量研究原理

理论上, 分子振动拉曼活性的散射强度(I)可简化表示为:

I=KNσIL

式中:IL为入射激光强度, K为比例系数, σ 为拉曼散射截面, N为分子数。表面看来, 拉曼散射强度与拉曼分子振动活性的分子数目(浓度)具线性关系, 散射强度可直接应用于定量研究。然而, 入射光的强度与测试条件密切相关, 测试设备的光学配置、样品校准以及扫描时间和次数等都会影响入射激光强度进而影响散射强度。因此, 若应用散射强度进行定量研究必须消除测试条件对散射强度的影响。

理论上, 应用相对强度消除测试条件影响是最有效的方法:

I* =IIR=NNR×σσR

式中:R为参考系, 其拉曼光谱在与样品完全相同的条件下测得; σ /σ R为散射截面比值, 可视为常数。因此, 相对拉曼强度I/IR就可以用于分子浓度的定量分析研究。

5.2 拉曼光谱CO2碳同位素组成定量研究进展及存在问题

早在20世纪70年代就有科学家提出可以利用激光拉曼光谱对传统同位素的同位素比值进行测定, 却鲜有发展。如上文所述, 激光拉曼光谱在氢氧同位素比值测定方面已有不少研究成果, 然而对于CO2中碳同位素比值的测定研究却非常少, 对流体包裹体中CO2的定量研究则多集中在气体密度和压力上[13, 44]

根据量子力学知识, 原子个数为N的线性分子具有3N-5个简正自由度。线性CO2分子有4个简正振动模式:ν 1 ( g+), 位于1 339.5 cm-1的对称伸缩振动; ν 2 (∑ u), 位于667.3 cm-1的二度兼并弯曲振动; ν 3 ( u+), 位于2 349.1 cm-1反对称伸缩振动。其中, ν 2 和ν 3 具红外活性, 而ν 1 具拉曼活性。2ν 2 ( g+) 和ν 1 ( g+)间的费米共振导致拉曼光谱中分别位于1 285 cm-1 ( ν-12)和1 388 cm-1 ( ν+12)的费米共振谱峰对。此外, 在1 265 cm-1和1 410 cm-1处存在对称的热峰。由于同位素效应, 13CO2的费米振动峰出现在1 370 cm-1( ν+13)和1 265 cm-1( ν-13)处, 但由于12CO2的热峰影响掩盖了其在1 265 cm-1处的振动峰, 所以在谱图中只能看到13CO2在1 370 cm-1处的峰, 且十分微弱, 在大气压下不易被检测出来。因此, 13CO2在1 370 cm-1处的拉曼峰强( I+13)被用于碳同位素组成的定量研究。室温下CO2气相拉曼光谱如图1所示。

图1 室温(25 ℃)下CO2(g)拉曼光谱Fig.1 Raman spectra of CO2(g) at ambient temperature(25 ℃)

Irmer[45]利用拉曼光谱对天然12CO2气体和纯度为99%的13CO2气体进行测试, 得到了12C16O2, 13C16O2, 12C16O18O, 12C16O18O的拉曼散射谱峰位置。Afe等[46]肯定了用拉曼光谱测定碳同位素的方法。

1998年, 赫英等[47]对胜利油田火山岩中的单个包裹体进行碳同位素的分析, 探讨CO2气藏及来源问题, 并通过以CO2特征峰面积比表征碳同位素比值来确定碳同位素的比值, 但实验精确度并不高。Menneken等[48]利用拉曼光谱对包裹体中CO2碳同位素的组成进行测定, 实验证明在富13C的CO2中采用13CO212CO2的峰强比值来定量碳同位素含量, 精度可达± 4‰ , 研究中他们还认为CO2密度也是影响峰强比的一个重要因素。

近年来, 研究最为详细的是:Arakawa等[49]通过人工制备不同碳同位素组成的CO2流体测得13CO212CO2的拉曼特征峰位置, 实验发现 I+13/ I+12值与碳同位素含量比值间存在线性关系, 以此规律来确定碳同位素的比值。研究表明, 拉曼光谱对流体包裹体中CO2中碳同位素比值测定的精确度为5‰ , 考虑到流体密度的差异以及曲线拟合过程造成的误差, 最后所得的实验误差为20‰ 。Arakawa等还利用拉曼光谱对近日本海沟火山岩中的地幔捕虏体中流体包裹体CO2碳同位素组成(区分碳来源)进行了研究, 验证了拉曼光谱对CO2碳同位素比值的测定的准确性。应当指出, Arakawa等实验中合成的流体包裹体中13CO212CO2的含量近乎相等(多数合成包裹体中13CO2含量大于12CO2), 这有利于提高测试数据的准确性, 但这样的配比与天然包裹体中CO2碳同位素组成相差较大, 因为自然界中的13C含量要远远小于12C的含量。也就是说上述实验条件与实际相差甚远。

理论上来说, 利用拉曼光谱对CO2中碳同位素组成进行研究是可行的, 当前主要问题存在于测试数据的准确性和可靠性方面。根据拉曼光谱定量原理分析可以看出, 13CO2/12CO2比值可以用13CO2的1 370 cm-112CO2的1 388 cm-1的谱峰强度( I+13I+12)比值来表征。在天然流体包裹体中, 13CO2含量远远小于12CO2的含量, 因此13CO2谱峰强度较弱, 弱峰会造成较大的数据处理误差, 也就是说13CO2的1 370 cm-1谱峰的强度越强, 在数据处理的过程中误差就越小, 越有利于改进数据的准确性。在拉曼光谱测试过程中, 增加样品的曝光时间是一个增强拉曼散射强度的方法, 然而增强13CO2的1 370 cm-1谱峰强度的同时往往造成12CO2在1 285 cm-1(或1 388 cm-1)处谱峰强度达到过饱和, 从而无法准确获得12CO2在1 285 cm-1(或1 388 cm-1)的谱峰强度, 以至于无法进行准确的谱峰强度的比值计算(这也许是Arakawa等合成的流体包裹体中13CO212CO2含量近乎相等的原因)。

作者近期对含CO2的人造流体包裹体(H2O-NaCl-CO2体系)进行了拉曼光谱测试, 体系密封于石英毛细管中, 包含气相CO2和5% NaCl水溶液两相。实验设备使用雷尼绍公司激光拉曼仪, 激光514.5 nm, 光栅密度1 800/mm, 显微镜使用50倍物镜, 对CO2气相同一个测试点分别进行单次10~350 s时长的扫描, 以获得信噪比高的光谱图(图2)。

实验发现, 扫描时间小于70 s时13CO2在1 370 cm-1处的峰强很弱(图2a), 数据精确度低; 当扫描时间增长到80 s之后, 13CO2峰略微能与背景噪声区分开(图2b); 但是当扫描时间超过100 s之后, 12CO2在1 388 cm-1处的振动峰达到过饱和, 无法进行峰强的测定(图2c, d)。

可以看出, 增加样品的曝光时间是一种有效提高拉曼光谱强度和光谱图信噪比的方式, 也是有效提高实验精度的方法。而如何能够避免拉曼谱峰过饱和的影响, 同时获得准确而强度较高的13CO212CO2( I+13I+12)的拉曼谱峰强度成为解决问题的核心。解决了这一问题, 激光拉曼光谱法应用于富CO2流体包裹体中碳同位素的质谱研究将能得到广泛应用。

图2 不同扫描时间下流体包裹体中CO2气相拉曼光谱Fig.2 Raman spectra of CO2(g) in fluid inclusion with different scan time

5.3 拉曼光谱CO2碳同位素组成测定方法建议

如上文中所分析的, 只有能够在同一次测试中获得高强度的 I+13I+12的值, 才能保证利用拉曼光谱进行碳同位素组成测定的准确度。也就是说在增长样品曝光时间的同时能够获得强度已经过饱和的 I+12是解决问题的关键。在CO2拉曼谱图中寻找一个合适的参数作为中间值来估算过饱和的 ν+12振动的拉曼峰强( I+12)将是一个有效的办法。

初步实验研究发现, 在CO2的拉曼光谱图中, 5个典型的拉曼峰( ν-12, ν+12, ν+13, 2个热峰), 只有12CO2在1 388 cm-1处的振动峰在增加扫描时长的时候最容易达到过饱和, 而其余4个峰仍然保持正常的峰形。在对人造流体包裹体的拉曼实验中, 利用PeakFit 4.0对所CO2获拉曼光谱图进行平滑, 去基线和高斯— 洛伦兹曲线拟合的处理后获得的谱峰数据中发现, 12CO2的2个费米共振峰的拉曼强度成良好的线性关系: I-12=2.35191× I+12-397.3413(R2=0.99993)。因此依据这样的关系, 当12CO2在1 388 cm-1处的振动峰拉曼强度( I+12)达到过饱和时, 12CO2在1 285 cm-1处的费米共振峰拉曼强度( I-12)可以做中间值估算 I+12的值, 最终能够计算得到准确的 I+13/ I+12的值用于计算碳同位素的组成δ 13C。

6 结 论

激光拉曼在单个流体包裹体中的CO2碳同位素比值的定量研究方面有可靠的理论依据, 前人在此基础上做了部分研究, 而且这种方法不仅在地球科学领域, 在医学、化学、生物学等其他领域的研究方面也有重要意义。激光拉曼的碳同位素组成研究在实验技术上已经相对成熟, 阻碍其广泛应用的是数据的准确性和可靠性, 解决这一问题的关键在于能够在单次测试中同时获得13CO212CO2的高拉曼散射强度( I+13I+12), 寻找一个合适的参考值作为中间值估算过饱和的 I+12将是一个有效的方法, 只要解决了这一问题, 激光拉曼在CO2碳同位素的质谱研究上将得到广泛应用。

The authors have declared that no competing interests exist.

参考文献
[1] Azbej T, Severs M J, Rusk B G, et al. In situ quantitative analysis of individual H2O-CO2 fluid inclusions by laser Raman spectroscopy[J]. Chemical Geology, 2007, 237(3): 255-263. [本文引用:1]
[2] Chen J, Zheng H, Xiao W, et al. Raman spectroscopic study of CO2 -NaCl-H2O mixtures in synthetic fluid inclusions at high temperatures 1[J]. Geochimica et Cosmochimica Acta, 2004, 68(6): 1 355-1 360. [本文引用:1]
[3] Zhang Zhou, Zhang Hongfu. Diamond and deep carbon cycle[J]. Earth Science Frontier, 2011, 18(3): 268-283.
[张舟, 张宏福. 金刚石与深部碳循环[J]. 地学前缘, 2011, 18(3): 268-283. ] [本文引用:1]
[4] Javoy M. The major volatile elements of the Earth: Their origin, behavior, and fate[J]. Geophysical Research Letters, 1997, 24(2): 177-180. [本文引用:1]
[5] Sleep N H, Zahnle K. Carbon dioxide cycling and implications for climate on ancient Earth[J]. Journal of Geophysical Research Atmospheres, 2001, 106(E1): 1 373-1 399. [本文引用:1]
[6] Zhao Xin, Shi Guanghai, Zhang Ji. Review of lithospheric diamonds and their mineral inclusions[J]. Advances in Earth Science, 2015, 30(3): 310-322.
[赵欣, 施光海, 张骥. 岩石圈地幔中的金刚石及其矿物包裹体的研究进展[J]. 地球科学进展, 2015, 30(3): 310-322. ] [本文引用:1]
[7] Deines P. The carbon isotope geochemistry of mantle xenoliths[J]. Earth-Science Reviews, 2002, 58(3/4): 247-278. [本文引用:1]
[8] Wang Xianbin, Wu Maobing, Zhang Mingjie. Stable isotope geochemistry of mantle fluids[J]. Geology-Geochemistry, 2000, 28(3): 69-73.
[王先彬, 吴茂炳, 张铭杰. 地幔流体的稳定同位素地球化学综述[J]. 地质地球化学, 2000, 28(3): 69-73. ] [本文引用:1]
[9] Roedder E. Liquid CO2 inclusions in olivine-bearing nodules and phenocrysts from basalts[J]. American Mineralogist, 1965, 50: 1 746-1 782. [本文引用:1]
[10] Xu Jiuhua, Ding Rufu, Xie Yuling, et al. Pure CO2 fluids of the Sarekoubu gold deposit in southern margin of Altai Mountains[J]. Chinese Science Bulletin, 2005, 50(4): 380-386.
[徐九华, 丁汝福, 谢玉玲, . 阿尔泰山南缘萨热阔布金矿床的纯CO2流体[J]. 科学通报, 2005, 50(4): 380-386. ] [本文引用:1]
[11] Chu Xuelei. Carbon isotopes in mantle[J]. Advances in Earth Science, 1996, 11(5): 446-452.
[储雪蕾. 地幔的碳同位素[J]. 地球科学进展, 1996, 11(5): 446-452. ] [本文引用:1]
[12] Chen Yanjing, Ni Pei, Fan Hongrui, et al. Diagnostic fluid inclusions of different types hydrothermal gold deposits[J]. Acta Petrologica Sinica, 2007, 23(9): 2 085-2 108.
[陈衍景, 倪培, 范洪瑞, . 不同类型热液金矿系统的流体包裹体特征[J]. 矿物学报, 2007, 23(9): 2 085-2 108. ] [本文引用:1]
[13] Mernagh T P, Bastrakov E N, Zaw K, et al. Comparison of fluid inclusion data and mineralization processes for Australian orogenic gold and intrusionrelated gold systems[J]. Acta Petrologica Sinica, 2007, 23(1): 21-32. [本文引用:2]
[14] Lu Huanzhang. Role of CO2 fluid in the formation of gold deposits: Fluid inclusion evidences[J]. Geochemical, 2008, 37(4): 321-328.
[卢焕章. CO2流体与金矿化: 流体包裹体的证据[J]. 地球化学, 2008, 37(4): 321-328. ] [本文引用:2]
[15] Phillips G N, Evans K A. Role of CO2 in the formation of gold deposits[J]. Nature, 2004, 429(6 994): 860-863. [本文引用:1]
[16] Klein E L, Fuzikawa K. Origin of the CO2-only fluid inclusions in the Palaeoproterozoic Cararvein-quartz gold deposit, Ipitinga Auriferous district, SE-Guiana Shield, Brazil: Implications for orogenic gold mineralisation[J]. Ore Geology Reviews, 2010, 37(1): 31-40. [本文引用:1]
[17] Li Jijun, Cao Xicheng, Wang Weiming, et al. Analysis and application of components and isotopic compositions of fluid inclusions in volcanic reservoirs of the Xujiaweizi fault depression[J]. Oil & Gas Geology, 2014, 35(1): 124-130.
[李吉君, 曹喜铖, 王伟明, . 徐家围子断陷火山岩包裹体组分和同位素组成分析及应用[J]. 石油与天然气地质, 2014, 35(1): 124-130. ] [本文引用:1]
[18] Wang Yong, Hou Zengqian, Mo Xuanxue, et al. Stable isotope characteristics and origin of ore-forming fluids in copper-gold polymetallic deposits within strike-slip pull-apart basin of Weishan-Yongping continental collision orogenic belt, Yunnan Province, China[J]. Mineral Deposits, 2006, 25(1): 60-70.
[王勇, 侯增谦, 莫宣学, . 云南巍山—永平碰撞造山带走滑拉分盆地铜金多金属矿成矿流体系统: 稳定同位素特征及热液来源[J]. 矿床地质, 2006, 25(1): 60-70. ] [本文引用:1]
[19] Mumm A S, Oberthür T, Vetter U, et al. High CO2 content of fluid inclusions in gold mineralisations in the Ashanti Belt, Ghana: A new category of ore forming fluids?[J]. Mineralium Deposita, 1997, 32(2): 107-118. [本文引用:1]
[20] Xie Yihan, Fan Hongrui, Wang Yinglan. Hydrocarbon inclusions indicate hydrocarbon maturity and evolution[J]. Marine Origin Petroleum Geology, 1999, 4(3): 32-35.
[谢奕汉, 范宏瑞, 王英兰. 石油包裹体与石油成熟度及油气演化[J]. 海相油气地质, 1999, 4(3): 32-35. ] [本文引用:1]
[21] Dai Jinxing, Dai Chunsen, Song Yan, et al. Inorganic genetic carbon dioxide gas accumulation and their characteristics in east part of China[J]. China offshore Oil and Gas(Geology), 1994, 8(4): 215-222.
[戴金星, 戴春森, 宋岩, . 中国东部无机成因的二氧化碳气藏及其特征[J]. 中国海上油气(地质), 1994, 8(4): 215-222. ] [本文引用:1]
[22] Li Hongwei, Feng Lianjun, Chen Jian. A rapid method for determination of carbon isotope of inclusions by sealed quartz tube decrepitation[J]. Chinese Journal of Analytical Chemistry, 2014, 42(1): 127-130.
[李洪伟, 冯连君, 陈健. 密封石英管爆裂法快速分析包裹体中CO2碳同位素[J]. 分析化学, 2014, 42(1): 127-130. ] [本文引用:1]
[23] Fan Qicheng, Liu Ruoxin, Lin Zhuoran, et al. Preliminary study of carbon isotope in mantle CO2 fluid inclusion from eastern China[J]. Geochimica, 1996, 25(3): 264-269.
[樊祺诚, 刘若新, 林卓然, . 中国东部地幔CO2流体包裹体的碳同位素初步研究[J]. 地球化学, 1996, 25(3): 264-269. ] [本文引用:1]
[24] Nadeau S, Pineau F, Javoy M, et al. Carbon concentrations and isotopic ratios in fluid-inclusion-bearing upper-mantle xenoliths along the northwestern margin of North America[J]. Chemical Geology, 1990, 81(4): 271-297. [本文引用:1]
[25] Pineau F, Mathez E A. Carbon isotopes in xenoliths from the Hualalai Volcano, Hawaii, and the generation of isotopic variability[J]. Geochimica et Cosmochimica Acta, 1990, 54(1): 217-227. [本文引用:1]
[26] Liu Deliang, Sun Xianru, Li Zhensheng, et al. Carbon and oxygen isotope analysis of fluid inclusions in ordovician carbonate veins of Ordos Basin[J]. Acta Petrolei Sinica, 2007, 28(3): 68-74.
[刘德良, 孙先如, 李振生, . 鄂尔多斯盆地奥陶系碳酸盐岩脉流体包裹体碳氧同位素分析[J]. 石油学报, 2007, 28(3): 68-74. ] [本文引用:1]
[27] Ma Lichun, Tang Qingfeng, Zhang Qi, et al. Advances in chemical composition determination of individual inclusions trapped within evaporite minerals[J]. Advances in Earth Science, 2014, 29(4): 475-481.
[马黎春, 汤庆峰, 张琪, . 蒸发岩矿物单个流体包裹体成分测定方法研究进展[J]. 地球科学进展, 2014, 29(4): 475-481. ] [本文引用:1]
[28] Li Xiangxian, Gao Minguang, Xu Liang, et al. Carbon isotope ratio analysis in CO2 based on fourier transform infrared spectroscopy[J]. Acta Physica Sinica, 2013, 62(3): 7-15.
[李相贤, 高闽光, 徐亮, . 基于傅里叶变换红外光谱法CO2气体碳同位素比检测研究[J]. 物理学报, 2013, 62(3): 7-15. ] [本文引用:1]
[29] Nieuwoudt M K, Simpson M P, Tobin M, et al. Synchrotron FTIR microscopy of synthetic and natural CO2-H2O fluid inclusions[J]. Vibrational Spectroscopy, 2014, 75: 136-148, doi: DOI:10.1016/j.vibspec.2014.08.003. [本文引用:1]
[30] Rosasco G J. Laser-excited raman spectroscopy for nondestructive partial analysis of individual phases in fluid inclusions in minerals[J]. Science, 1975, 190(4 214): 557-560. [本文引用:1]
[31] Xi Binbin, Shi Weijun, Jiang Hong, et al. Synthesis of N2 -CH4 and N2 -CO2 gas mixtures as an online stand ard and determination of their Raman quantification factors of CH4 and CO2[J]. Rock and Mineral Analysis, 2014, 33(5): 655-660.
[席斌斌, 施伟军, 蒋宏, . N2-CH4(CO2)混合气体在线标样制备及其拉曼定量因子测定[J]. 岩矿测试, 2014, 33(5): 655-660. ] [本文引用:1]
[32] Chen Yong, Zhou Yaoqi, Zhang Dagang. Trying to determine the salinity of hydrous liquid by laser Raman Spectroscopy[J]. Chinese Journal of Light Scattering, 2003, 14(4): 216-223.
[陈勇, 周瑶琪, 章大港. 几种盐水溶液拉曼工作曲线的绘制[J]. 光散射学报, 2003, 14(4): 216-223. ] [本文引用:1]
[33] Yang Dan, Xu Wenyi. Gauss-lorentz quantitative research on O-H stretching raman spectra of water in common chlorine salt solution[J]. Spectroscopy and Spectral Analysis, 2013, 33(2): 376-382.
[杨丹, 徐文艺. 常见氯盐溶液水的拉曼包络线高斯—洛仑兹去卷积分峰定量研究[J]. 光谱学与光谱分析, 2013, 33(2): 376-382. ] [本文引用:1]
[34] Wang Xiaolin, Hu Wenxuan, Chou I-Ming. Raman spectroscopic characterization on the OH stretching band s in NaCl-Na2CO3-Na2SO4-CO2-H2O systems: Implications for the measurement of chloride concentrations in fluid inclusions[J]. Journal of Geochemical Exploration, 2013, 132(3): 111-119. [本文引用:1]
[35] Sun Qiang, Zhao Lei, Li Nuo, et al. Raman spectroscopic study for the determination of Cl- concentration (molarity scale) in aqueous solutions: Application to fluid inclusions[J]. Chemical Geology, 2010, 272(1/4): 55-61. [本文引用:1]
[36] Li Shuling. Study on Raman spectra of CO2-H2O system[J]. Rock and Mineral Analysis, 1998, 17(3): 177-180.
[李淑玲. 二氧化碳—水体系的拉曼光谱研究[J]. 岩矿测试, 1998, 17(3): 177-180. ] [本文引用:1]
[37] Rosso K M, Bodnar R J. Microthermometric and Raman spectroscopic detection limits of CO2 in fluid inclusions and the Raman spectroscopic characterization of CO2[J]. Geochimica et Cosmochimica Acta, 1995, 59(19): 3 961-3 975. [本文引用:1]
[38] Mernagh T P, Wilde A R. The use of the laser Raman microprobe for the determination of salinity in fluid inclusions[J]. Geochimica et Cosmochimica Acta, 1989, 53(4): 765-771. [本文引用:1]
[39] Ding Junying, Ni Pei, Rao Bing, et al. Evaluation of the laser Raman microprobe method for the determination of salinity in a single fluid inclusion by using synthetic fluid inclusions[J]. Geological Review, 2004, 50(2): 203-209.
[丁俊英, 倪培, 饶冰, . 显微激光拉曼光谱测定单个包裹体盐度的实验研究[J]. 地质论评, 2004, 50(2): 203-209. ] [本文引用:1]
[40] Li Jing, Mao Shide, Zheng Haifei. Raman spectroscopic quantitative study of NaCl-CaCl2-H2O system at high temperatures and pressures[J]. Spectroscopy and Spectral Analysis, 2014, 34(7): 1 747-1 753. [本文引用:1]
[41] Bloom S D, Harney R C, Milanovich F P. Stable isotope ratio measurements in nitrogen and oxygen using Raman scattering[J]. Applied Spectroscopy, 1976, 30(1): 64-67. [本文引用:1]
[42] Huang Gang, Cao Xiaohua, Long Xinggui. Application of laser raman spectroscopy to hydrogen isotope analsis[J]. Spectroscopy and Spectral Analysis, 2011, 31(1): 30-33.
[黄刚, 曹小华, 龙兴贵. 激光拉曼光谱分析氢同位素的应用[J]. 光谱学与光谱分析, 2011, 31(1): 30-33. ] [本文引用:1]
[43] Wang Shixia, Zheng Haifei. Research on the experiment of hydrogen isotope fractionation using diamond anvil cell and Raman spectra[J]. Spectroscopy and Spectral Analysis, 2011, 31(3): 230-236.
[王世霞, 郑海飞. 金刚石压腔结合拉曼光谱技术进行氢同位素分馏的实验研究[J]. 光谱学与光谱分析, 2011, 31(3): 230-236. ] [本文引用:1]
[44] Wang X, Chou I M, Hu W, et al. Raman spectroscopic measurements of CO2, density: Experimental calibration with High-Pressure Optical Cell (HPOC) and Fused Silica Capillary Capsule (FSCC) with application to fluid inclusion observations[J]. Geochimica et Cosmochimica Acta, 2011, 75(14): 4 080-4 093. [本文引用:1]
[45] Irmer G. Isotopes of C and O in CO2: A Raman study using gas stand ards and natural fluid inclusions[J]. Acta Universitatis Carolinae-Geologica, 2002, 46(1): 35-36. [本文引用:1]
[46] Afe O, Jaeger F, Wagner G. Detection of 13CO2 and 12CO2 using diode laser driven Raman scattering[J]. Geophysical Research Abstracts, 2007, 9: 3 071-3 072. [本文引用:1]
[47] He Ying, Zhu Xingguo, Xu Peicang. The isotopic compositions of CO2 for single fluid inclusions in volcanic rocks from Shengli Oil Field[J]. Chinese Journal of Geology, 1998, 33(3): 380-383.
[赫英, 朱兴国, 徐培苍. 胜利油田火山岩单个包裹体碳同位素[J]. 地质科学, 1998, 33(3): 380-383. ] [本文引用:1]
[48] Menneken M, Geisler T. An evaluation of the potential of using Raman spectroscopy to determine the carbon isotope composition of CO2 inclusions[J]. Journal of Geochemical Exploration, 2009, 101(1): 70. [本文引用:1]
[49] Arakawa M, Yamamoto J, Kagi H. Developing micro-raman mass spectrometry for measuring carbon isotopic composition of carbon dioxide[J]. Applied Spectroscopy, 2007, 61(7): 701-705. [本文引用:1]