作者简介:药瑛(1990-),女,山西长治人,硕士研究生,主要从事实验地球化学研究.E-mail:yaoying@pku.edu.cn
CO2作为一种流体包裹体中常见的组分,对其中碳同位素的定量研究,能够提供很多地质历史时期地质流体源区的地球化学信息,并为地球碳循环的研究奠定基础。传统的质谱测试存在无法进行单个流体包裹体测试的缺点。拉曼光谱因其快速、无损、高精度的特点被广泛应用于地球科学各个研究领域,尤其在单一流体包裹体的研究中更为重要。当前,拉曼光谱在定性研究物质结构和成分方面已经十分成熟,相比之下,拉曼光谱定量研究仍显不足。结合前人的研究和实际的实验,根据对拉曼光谱基本原理的分析,提出激光拉曼光谱定量研究应当基于相对散射强度进行,且利用13CO2和12CO2的拉曼谱峰强度比(
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
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 (
地质流体在地质成岩及成矿过程中发挥着重要作用, 且以流体包裹体形式保存下来。对流体包裹体组成的研究, 是正确认识地质流体演化及流体与成矿作用关系的前提。拉曼光谱因具有与测试样品非接触、非破坏性、样品制备简单、操作方便、精度高等技术特点, 已经被普遍应用于流体包裹体的研究。流体包裹体以H2O-NaCl-CO2体系为代表[1, 2], 对这类流体包裹体的研究是研究更复杂流体包裹体的基础。
当前拉曼光谱主要应用于定性研究, 即根据振动频率进行分子类型的确定, 而定量研究明显不足。加强拉曼光谱定量化方面的研究是进一步发展激光拉曼技术所必需的。本文结合前人的研究方法和成果对应用拉曼光谱定量研究H2O-NaCl-CO2体系二氧化碳(CO2)的碳同位素组成进行了探讨。
碳与地球的形成和发展息息相关, 绝大部分(> 90%)的碳被储存在地球深部, 却又与地表碳有着千丝万缕的联系[3, 4]。研究地球内部的碳循环对地球演化史、矿床成因、全球气候变化[5]等问题的研究具有重要意义。在碳循环的研究中, 碳同位素组成是一个重要指标, 通常用δ 13CPDB来表征。地球碳循环的研究对象有单质碳和化合态碳, 其中金刚石作为典型的地球深部来源的单质碳被认为是研究深部碳循环的重要样品[6], 而化合态碳中CO2则是重要的研究对象之一, 其他还包括碳酸盐矿物、CH4和SiC等[7]。
CO2广泛存在于自然界中, 参与并影响了许多地质作用过程。地质作用过程中, CO2可以以气相、液相的形式被流体包裹体保存下来, 因此流体包裹体中CO2碳同位素组成的研究, 在研究地壳、地幔循环过程以及矿床、油气的形成等许多问题上都具有重要意义。
地幔流体的研究一直被认为是研究地球内部物质组成和演化的重要手段。地幔是地球碳的主要源区之一, 碳循环受地幔碳同位素组成和不均一分布的影响。结合来自地幔不同深度的流体包裹体样品中碳同位素的组成特征与其他的同位素信息能够很好地研究地幔流体的形成和演化。而且统计全球范围内地幔流体的同位素组成也有利于地幔物质的全球性对比[8]。
自1965年Roedder[9]在碱性玄武岩的超镁铁质捕掳体中发现大量CO2包裹体开始, 掀起了通过流体包裹体研究地幔流体的浪潮。研究证明, 富CO2的流体包裹体在地幔岩中普遍存在, 上地幔流体、下地壳麻粒岩化均具有富纯CO2流体的特征[10]。科学家们结合对不同地区的地幔岩石中的包裹体、金刚石、金伯利岩、碳酸盐、大洋玄武岩的碳同位素研究, 发现地幔流体中的δ 13C分布范围广, 但是大部分样品的δ 13C值集中分布在-5‰ 附近和-15‰ ~-25‰ [11]。这为研究地幔碳的多样性和原始地幔碳是否均一奠定了基础, 也为研究地球深部碳循环提供了手段。
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体系的研究很大程度上完善了我们对成矿作用的认识。
对流体包裹体中CO2碳同位素组成进行研究具有重要意义, 因此对其分析测试一直为研究者所关注, 并开展了大量的研究工作。流体包裹体中CO2碳同位素组成分析测试方法可大致分为破坏性分析测试方法和无损分析测试方法。
所谓破坏性分析测试方法是指通过破坏流体包裹体提取其中的待检测物质, 然后利用质谱仪进行包裹体内同位素含量以及比值的测定。针对群体流体包裹体进行分析测试的方法主要有机械破碎法、热爆法和分阶段加热法。近年来还发展了可进行单个流体包裹体的破坏性分析测试— — 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对单一包裹体进行的测试, 需要破坏样品, 因此不能进行重复测试和数据验证。
当前, 流体包裹体的无损分析测试主要应用红外光谱和激光拉曼光谱。理论上, 双原子分子简谐振动可表示为:
V=
式中: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的碳同位素组成的测定方面的研究和应用却非常缺乏。
1975年, Rosasco[30]利用拉曼光谱对天然流体包裹体进行了研究, 提出拉曼光谱在定性和定量研究方面应用的可能性。之后, 有越来越多的学者在该方面开展了工作, 显示了拉曼光谱在该领域的研究优势。近年来, 激光拉曼光谱已被广泛应用于流体包裹体的定性和定量研究中。对单个流体包裹体的定量研究集中在对其气液相组分、同位素化学组成分析及包裹体内的气体密度和压力、有机热成熟度等参数的测定[31]。
当前, 拉曼光谱定量研究过程中, 根据参考体系的类型可大致分为外标法和内标法。所谓内标法, 是指参考系与被测试对象在测试过程中, 拉曼光谱被单次测试同时记录下来, 从而保证了测试条件的一致。外标法则是在保证测试条件一致的条件下, 参考系与被测试对象先后测试完成。可以看出, 相比于外标法, 内标法更加准确快捷。
流体包裹体成分定量研究方面, 具有拉曼活性的分子, 如:CO2, CH4, N2等利用内标为参考系进行物质含量的定量研究。众多溶于水却在常温下无拉曼光谱活性的阴离子(Cl-, C
包裹体密度定量研究方面, 李淑玲[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中碳同位素组成方面的定量研究具有重要意义。
理论上, 分子振动拉曼活性的散射强度(I)可简化表示为:
式中:IL为入射激光强度, K为比例系数, σ 为拉曼散射截面, N为分子数。表面看来, 拉曼散射强度与拉曼分子振动活性的分子数目(浓度)具线性关系, 散射强度可直接应用于定量研究。然而, 入射光的强度与测试条件密切相关, 测试设备的光学配置、样品校准以及扫描时间和次数等都会影响入射激光强度进而影响散射强度。因此, 若应用散射强度进行定量研究必须消除测试条件对散射强度的影响。
理论上, 应用相对强度消除测试条件影响是最有效的方法:
式中:R为参考系, 其拉曼光谱在与样品完全相同的条件下测得; σ /σ R为散射截面比值, 可视为常数。因此, 相对拉曼强度I/IR就可以用于分子浓度的定量分析研究。
早在20世纪70年代就有科学家提出可以利用激光拉曼光谱对传统同位素的同位素比值进行测定, 却鲜有发展。如上文所述, 激光拉曼光谱在氢氧同位素比值测定方面已有不少研究成果, 然而对于CO2中碳同位素比值的测定研究却非常少, 对流体包裹体中CO2的定量研究则多集中在气体密度和压力上[13, 44]。
根据量子力学知识, 原子个数为N的线性分子具有3N-5个简正自由度。线性CO2分子有4个简正振动模式:ν 1 (
Irmer[45]利用拉曼光谱对天然12CO2气体和纯度为99%的13CO2气体进行测试, 得到了12C16O2, 13C16O2, 12C16O18O, 12C16O18O的拉曼散射谱峰位置。Afe等[46]肯定了用拉曼光谱测定碳同位素的方法。
1998年, 赫英等[47]对胜利油田火山岩中的单个包裹体进行碳同位素的分析, 探讨CO2气藏及来源问题, 并通过以CO2特征峰面积比表征碳同位素比值来确定碳同位素的比值, 但实验精确度并不高。Menneken等[48]利用拉曼光谱对包裹体中CO2碳同位素的组成进行测定, 实验证明在富13C的CO2中采用13CO2和12CO2的峰强比值来定量碳同位素含量, 精度可达± 4‰ , 研究中他们还认为CO2密度也是影响峰强比的一个重要因素。
近年来, 研究最为详细的是:Arakawa等[49]通过人工制备不同碳同位素组成的CO2流体测得13CO2和12CO2的拉曼特征峰位置, 实验发现
理论上来说, 利用拉曼光谱对CO2中碳同位素组成进行研究是可行的, 当前主要问题存在于测试数据的准确性和可靠性方面。根据拉曼光谱定量原理分析可以看出, 13CO2/12CO2比值可以用13CO2的1 370 cm-1与12CO2的1 388 cm-1的谱峰强度(
作者近期对含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)。
可以看出, 增加样品的曝光时间是一种有效提高拉曼光谱强度和光谱图信噪比的方式, 也是有效提高实验精度的方法。而如何能够避免拉曼谱峰过饱和的影响, 同时获得准确而强度较高的13CO2和12CO2(
如上文中所分析的, 只有能够在同一次测试中获得高强度的
初步实验研究发现, 在CO2的拉曼光谱图中, 5个典型的拉曼峰(
激光拉曼在单个流体包裹体中的CO2碳同位素比值的定量研究方面有可靠的理论依据, 前人在此基础上做了部分研究, 而且这种方法不仅在地球科学领域, 在医学、化学、生物学等其他领域的研究方面也有重要意义。激光拉曼的碳同位素组成研究在实验技术上已经相对成熟, 阻碍其广泛应用的是数据的准确性和可靠性, 解决这一问题的关键在于能够在单次测试中同时获得13CO2和12CO2的高拉曼散射强度(
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|
[42] |
|
[43] |
|
[44] |
|
[45] |
|
[46] |
|
[47] |
|
[48] |
|
[49] |
|