Quantitative Determination of Sulfate in Hydrate and Water Systems by Raman Spectrometry
-
摘要: 硫酸根含量是天然气水合物孔隙水中重要的参数之一，硫酸根与水合物的赋存关系已成为天然气水合物地球化学异常研究的重要组成部分。为了更深入地了解硫酸根浓度异常产生的原因和机制，建立一种实时、快速的测定方法就显得尤为重要。本文利用自行设计的天然气水合物反应装置，建立了海水中硫酸根在高压低温环境下的拉曼光谱定量测定方法。硫酸根的浓度（ρ）在2.0~70.0 g/L范围内，与其拉曼光谱参数（R）有良好的线性关系，线性回归方程为R=0.0013ρ-0.000066，相关系数r = 0.9998。方法检出限为0.2 g/L，精密度（RSD，n=12）为2.5%~3.4%，回收率为102.1%~123.8%。测试结果表明，采用拉曼光谱法可在线检测天然气水合物-水体系中硫酸根的浓度变化；相对于离子色谱法，具有不破坏原样品和方便快捷等优点。通过对硫酸根浓度变化的初步分析，发现其浓度呈阶梯状升高，认为天然气水合物的形成是一个气体溶解-成核-生长反复进行的过程，直至整个反应体系达到平衡。
Abstract: Sulfate content is one of the important parameters in pore water of gas hydrate. The relationship between sulfate and gas hydrate has become an important method to study geochemical anomalies of gas hydrate. It is particularly important to establish a fast and real time method in order to determine and therefore understand the cause and mechanism of sulphate concentration anomalies. An explanation of a uniquely designed synthetic experiment for gas hydrate is described in this paper, which is a method for the quantitative determination of sulfate in seawater using Raman Spectrometry at high pressure and low temperature. The relationship between the sulfate content and its Raman spectrum parameters was in a good linear agreement within 2.0-70.0 g/L. The equation of linear regression is R = 0.0013ρ-0.000066. The related coefficient is r = 0.9998. The detection limit is 0.2 g/L. The recoveries are 102.1%-123.8% with precision of 2.5%-3.4% (RSD, n = 12). The test result shows that the changes of sulfate content in gas hydrate and water system on-line by Raman spectroscopy can be detected. Compared to ion chromatography, the original samples remain intact with this method, and it is convenient and rapid. By preliminary analysis of the changes of sulfate content, it was founded that sulfate content basically presented ladder form rise and determined that the formation of gas hydrate was a repeated process of gas dissolution, nucleation and growth until the entire reaction system reaches equilibrium. -
激光拉曼光谱分析作为一种非破坏性的分析方法，可以快速方便地对单个包裹体进行定性、半定量分析，现已成为流体包裹体研究的基本工具之一[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]。国内外学者对包裹体中常见的C-H-O-N-S体系的气相组分开展了比较广泛的定量研究[3, 4, 5, 6, 7]，取得了显著的成果。由于气相组分的拉曼定量分析与分子性质、温度、压力、仪器性能等诸多因素有关[3, 4, 29]，造成前人结果存在比较明显的差异，难以相互借用，如李维华等[5]与Wopenka等[30]测定的SO2的定量因子有近5倍的差别。因此在进行气相成分的定量分析之前，需要利用一系列混合气体标样对仪器进行标定。前人一般使用商用钢瓶装混合气进行仪器标定[3, 4, 5]，虽然上述标样易于购置、配比准确，却存在气体组成单一无法调节、费用高、需要经常更换钢瓶等缺点。如按10%的梯度对10%~90%的两种气体的混合物进行标定，需要购置9瓶钢瓶气轮换使用，并且钢瓶气一定的使用期限，超过期限需要重新购置。针对上述不足，本文提出了一种在线配置不同浓度和压力条件下混合气体标样的方法，以实现快速准确地对激光拉曼探针进行标定及测定气体拉曼定量因子的研究目的。
1. 在线标样制备装置和在线标样的制备
为了实现混合气体标样的制备，本次研究搭建了一套在线标样制备装置(图 1)。该装置可以同时接入三路钢瓶气体，每路钢瓶气分别连接一个减压阀用于控制气体的输出压力；利用带有刻度和活塞的体积转移器量取实验所需体积的气体并将量取的气体注入高压容器中进行混合；增压泵用于对高压容器中的混合气体进行增压；真空泵用于对装置进行抽真空；装置的输出端与石英毛细管相连接；管路中安装有真空表以及压力表用于监控系统的真空度以及线路中气体的压力值；线路中还设有两个排气孔用于排气及管路清洗。
实验所用的钢瓶气为高纯气体，浓度≥99.999%；毛细管规格为内径0.1 mm，外径0.3 mm，表面涂有一层聚酰亚胺保护膜，厚度约0.025 mm(美国Polymicro Technologies公司)。激光拉曼分析的仪器为Renishaw Invia型激光拉曼光谱仪(英国Renishaw公司)，使用Ar+激光器，波长为514 nm，光谱分辨率为2 cm-1。
在线混合气体标样制备的实验步骤如下。
(1) 打开阀门1~6、8、10，关闭阀门7、9、11，打开真空泵对管路、体积转移器及高压容器抽真空，待真空表读数≤10Pa时，关闭真空泵。
(2) 关闭阀门2~4、6、8、10，打开气瓶1的减压阀并调节至实验所需压力值，用体积转移器量取实验所需气体体积。
(3) 关闭阀门1、5、气瓶1的减压阀，打开阀门6、8，将体积转移器中的气体转移至高压容器中。
(4) 关闭阀门8，打开阀门1~6、8、10，对系统抽真空，待真空表读数≤10Pa时，关闭真空泵。
(5) 重复步骤(2)~(4)，量取实验所需体积及压力条件下的气体2并注入到高压容器中，使气体1和2充分混合。
(6) 关闭阀门6，打开阀门8、11，利用高压容器中的混合气体对管路进行清洗。
(7) 关闭阀门11，打开阀门9，打开电动增压泵，对高压容器中的气体进行增压，待达到实验所需的气体压力时，停止增压并进行激光拉曼分析，然后继续增压至下一个压力点并进行拉曼分析。
2. 结果与讨论
2.1 在线样品准确性验证
为了验证制样方法的准确性及重复性，将本研究制备的70% N2+30% CO2的在线标样与购置于大连大特气体公司生产的同等浓度的商用标样，在10 MPa条件下分别进行了激光拉曼分析。结果表明，本次研究制备的混合气体与商用钢瓶装标样具有相似的峰形(图 2)。利用英国Renishaw公司出品的Wire3.0软件对上述拉曼谱图进行了分析，结果表明本方法制备的混合气体与商用标样具有相似的CO2与N2的相对峰高以及相对峰面积值，其相对误差小于4%，并具有较好的重现性，能够满足实验要求。
![]()
图 2 商用标样与在线样品拉曼谱图Figure 2. The Raman spectra of commercial standard sample and on-line mixing sample2.2 CH4及CO2相对拉曼定量因子的测定
在测定单个包裹体气体组成方面，国内外多沿用“相对拉曼定量因子”的方法，即通常将N2的定量因子定为1.00，其他气体与N2进行比较，得到相对拉曼定量因子[3, 4]。本次研究分别对拉曼峰面积及峰高计算了相对拉曼定量因子，具体公式如下：
式中，Ag为气体g的拉曼峰面积；AN2为N2的拉曼峰面积；Cg为气体g的摩尔分数；CN2为N2的摩尔分数；Hg为气体g的拉曼峰高；HN2为N2的拉曼峰高；Fgr代表以峰面积为参考值时气体g相对于N2的拉曼定量因子；Ggr代表以峰高为参考值时气体g相对于N2的拉曼定量因子。
为了测定CO2以及CH4的相对拉曼定量因子，在室温、5 MPa和10 MPa压力条件下，分别制备了N2摩尔分数为30%、50%和70%的N2-CO2混合气体标样以及N2-CH4混合气体标样。
在上述标样的激光拉曼谱图(图 3)中能清晰地辨识出N2、CO2以及CH4的拉曼特征峰。气体的拉曼峰强度随浓度以及压力的增加而增加，信噪比随着压力由5 MPa增加到10 MPa增大约一倍。
![]()
图 3 N2-CO2以及N2-CH4在线混合气体拉曼谱图Figure 3. The Raman spectra of N2-CO2 and N2-CH4 on-line gas mixturesThis page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
虽然CO2在1286 cm-1附近以及1386 cm-1附近出现两个峰值，但是由于1286 cm-1附近的峰强度要小于1386 cm-1附近峰强度。因此本文仅针对CO2在1386 cm-1附近的峰计算了相对拉曼定量因子。
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
求得CH4和CO2相对拉曼定量因子之后，便可以对包裹体中CH4和CO2的相对含量进行计算，具体计算公式如下：
3. 地质样品应用
选取四川金沙岩孔剖面，震旦系的藻云岩样品进行应用研究。该样品溶洞发育，被后期亮晶白云石充填。溶洞充填的亮晶白云石中发育气液两相盐水包裹体。选取个体较大并且靠近样品表面的包裹体，对其气泡进行激光拉曼分析，结果表明包裹体的气泡主要由CH4和CO2组成(图 6)。
利用wire3.0对图 6中两个包裹体的拉曼相关参数进行求解，并分别利用公式(3) 和(4) 对包裹体a和b中的CH4和CO2摩尔浓度进行了计算，得到包裹体中CH4的摩尔分数为27.60%~31.63%，CO2的摩尔分数为68.37%~72.40%(表 1)。上述结果表明，利用本文所求得的拉曼定量因子F和G所得到计算的结果基本一致(两者的绝对偏差在2.5%以内)；包裹体a和b气相组成较接近，可能为同期捕获的产物。
表 1 包裹体样品分析结果Table 1. The analytical composition of gas in fluid inclusions包裹体 ACO2 HCO2 ACH4 HCH4 CCH4(%) CCO2(%) 据公式(3) 据公式(4) 据公式(3) 据公式(4) 包裹体a 3461.54 594.541 17891.2 4115.24 31.63 31.25 68.37 68.75 包裹体b 3137.87 732.481 14694.8 4251.27 29.54 27.60 70.46 72.40 4. 结语
本文利用自主搭建的在线标样制备装置，对N2-CH4以及N2-CO2进行在线混合增压，制备了N2摩尔浓度为30%、50%和70%，压力为5 MPa和10 MPa的N2-CH4以及N2-CO2混合气体在线标样。通过与商用混合钢瓶气体标样对比表明，该方法所使用的装置操作简单，制备的混合气体具有较高的准确性及重现性，能够方便、准确地对拉曼光谱仪进行标定，实现了不同压力和浓度条件下气体的相对拉曼定量因子的测定。通过对CH4及CO2的相对定量因子测定表明，气体压力在5~10 MPa的范围时，定量因子不受压力变化的影响，为固定值。地质样品应用表明，本方法可以方便、灵活、准确地按任意比例将两瓶及两瓶以上纯气体钢瓶样品进行混合及增压，弥补了商用钢瓶装混合气体标样费用高、气体组成单一固定等不足。
由于本次研究仅在5 MPa和10 MPa两个压力点进行了分析，因此对于相对定量因子在 < 5 MPa及 > 10 MPa压力条件下的变化规律还有待于进一步研究。另外由于缺乏已知气体组成的人工合成包裹体标样，对于本方法在包裹体应用中的误差范围还有待于进一步研究。
-
![]()
图 1 天然气水合物模拟实验离子参数实时检测装置
Figure 1. Apparatus for ion parameters real-time detection of gas hydrate simulation experiment
![]()
图 2 高压透明腔内液相样品的激光拉曼光谱测试
Figure 2. Laser Raman spectrum test in liquid samples in high transparent cavity
![]()
图 3 不同浓度的K2SO4溶液拉曼光谱图(1.0×105 Pa，25℃)
a—b中硫酸根特征峰的放大图；b—500~4000 cm-1范围内溶液光谱图。
Figure 3. Raman spectrums of different concentrations of K2SO4 solution (1.0×105 Pa, 25℃)
表 1 仪器工作条件
Table 1 Working parameters of the laser Raman spectrometry instrument
工作参数 设定条件 激光光源 532 nm 激光功率 80 mW 物镜 20倍 光栅刻线数 2400线/mm 拉曼位移范围 500~4000 cm-1 积分时间 1~10 s 
下载: 导出CSV
表 2 方法精密度实验
Table 2 Precision tests of the method
K2SO4标准样品的浓度(g/L) K2SO4标准溶液测定浓度 RSD(%) 分次测定值 平均值 5.0 g/L 5.6 5.3 5.7 5.2 5.8 5.5 5.5 3.4 5.8 5.5 5.5 5.4 5.7 5.4 35.0 g/L 37.4 36.0 36.7 37.9 36.4 36.5 36.4 2.5 36.5 36.1 36.0 36.6 36.5 34.1
下载: 导出CSV
表 3 方法回收率
Table 3 Recovery tests of the method
海水样品(50 ml) 加入量(g) 测定量(g) 回收量(g) 回收率(%) 平均回收率(%) RSD(%) 4.9 g/L 0.0500 0.3078 0.0613 122.6 123.8 23.6 0.2941 0.0476 95.1 0.3233 0.0768 153.6 0.1000 0.3521 0.1056 105.6 109.2 3.6 0.3552 0.1078 108.7 0.3599 0.1134 113.7 0.1500 0.3910 0.1445 96.3 102.1 5.7 0.3994 01529 101.9 0.4085 0.1620 108.0
下载: 导出CSV
表 4 不同温度下K2SO4溶液浓度
Table 4 Solution concentration of K2SO4 in different temperatures
样品浓度(g/L) 温度(℃) 测定浓度(g/L) 相对误差(%) 5.0 2 5.1 1.9 5.0 5 5.1 2.3 5.0 10 5.1 1.0 5.0 15 5.0 -0.2 5.0 20 4.9 -1.3 5.0 25 4.9 -1.6
下载: 导出CSV
表 5 不同压力下K2SO4溶液浓度
Table 5 Solution concentration of K2SO4 in different pressures
样品浓度(g/L) 压力(MPa) 测定浓度(g/L) 相对误差(%) 5.0 0 4.9 -1.2 5.0 2 5.0 0.9 5.0 4 5.1 1.0 5.0 6 5.0 -0.2 5.0 8 4.9 -1.3
下载: 导出CSV
-
计量
- 文章访问数: 1422
- HTML全文浏览量: 338
- PDF下载量: 9
京公网安备 11010202008159号