水热金刚石压腔结合拉曼光谱技术进行Na2SO4-水体系溶解结晶动力学研究
王世霞1, 杨梦1, 吴嘉2, 郑海飞3
1. 上海理工大学理学院, 上海 200093
2. 中国石油大学(北京)油气资源与探测国家重点实验室, 北京 102249
3. 北京大学造山带与地壳演化教育部重点实验室, 北京 100871

作者简介: 王世霞, 女, 1983年生, 上海理工大学理学院讲师 e-mail: wangshixia@usst.edu.cn

摘要

成矿作用过程中, 温压条件改变导致矿物溶解重结晶, 溶液中溶质浓度发生变化。 从溶液中析出晶体的粒度同时存在着时间和空间的分布, 是复杂的动力学过程。 当前对矿物在流体中溶解重结晶动力学研究主要使用高压釜或活塞圆筒等封闭设备测定溶液溶质浓度的变化或固相矿物的形态变化, 降温淬火反应会影响样品的真实组成。 使用可进行原位观测的金刚石压腔结合拉曼光谱分析技术, 研究无水芒硝-饱和Na2SO4溶液随体系温度压力变化所出现的晶体溶解重结晶过程。 通过原位观测无水芒硝溶解、 结晶变化来探究硫酸钠晶体在不同温压条件下的溶解结晶动力学反应机制。 结果表明常温条件下无水芒硝拉曼位移分别位于449.9, 620.5, 632.9, 647.4, 993.3, 1 101.8, 1 132.2和1 153.1 cm-1。 随着体系温度的缓慢升高, 固相Na2SO4的晶形不断发生变化, 温度至193 ℃时无水芒硝(Na2SO4)完全溶解, 降温重结晶出现了新的1 196.5 cm-1拉曼特征峰, 重结晶晶体为芒硝(Na2SO4·10H2O); 金刚石原位观测结果显示迅速升温过程中无水芒硝发生部分溶解重结晶, 重结晶区域拉曼特征峰显示依然为无水芒硝。 拉曼光谱定量分析结果显示, 溶液中SO42-, H2O的拉曼谱峰面积比值参数更能反映SO42-浓度的变化, 体系达到溶解重结晶平衡状态时, SO42-/H2O峰面积比值AR为(0.016 6±0.000 4), 误差为2.4%。 应用Johnson-Mehl-Avrami-Kolmogorov(JMAK)模型结合溶液中SO42-/H2O峰面积比值变化对体系中固相无水Na2SO4的溶解重结晶过程进行动力学拟合, 计算得出无水硫酸钠在109 ℃条件下的溶解结晶反应的反应级数为1.266 7, 反应平衡常数为0.001 26。 综上所述, 水热金刚石压腔装置实验步骤少, 过程简便, 可避免由于淬火过程中退化交换作用等造成的误差, 应用水热金刚石压腔原位观测的装置优势结合拉曼光谱定量分析技术可实现高温高压条件下矿物在水溶液中溶解结晶动力学过程的原位观察和测定, 是一种高效的动力学研究手段。

关键词: HDAC; 高温高压; 拉曼光谱; Na2SO4-水体系; 溶解重结晶; 动力学
中图分类号:O657.3 文献标志码:A
Research on Dissolution Crystallization Kinetics of Na2SO4-H2O System Study Using Hydrothermal Diamond Anvil Cell and Raman Spectra
WANG Shi-xia1, YANG Meng1, WU Jia2, ZHENG Hai-fei3
1. School of Science, University of Shanghai for Science and Technology, Shanghai 200093, China
2. State Key Laboratory of Petroleum Resources and Prospecting-College of Geoscience, China University of Petroleum, Beijing102249, China
3. Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, Peking University, Beijing 100871, China
Abstract

During the mineralization process, the dissolution of primary mineral and the formation of secondary mineral could happen on the conditions of changing temperature and pressure. The dissolution and recrystallization of the minerals would cause the changing concentration of the solute, and the crystal which is from recrystallization depends on the reaction process. The process of dissolution and recrystallization is a complex dynamic process. At present, high-pressure autoclave and piston-cylinder are mainly used for the study of the kinetics on the dissolution and recrystallization of the minerals, whereas the cooling quenching reaction will affect the true composition of the sample. In this experiment, the process of the crystal recrystallization of thenardite saturated Na2SO4 solution with the change of temperature and pressure was traced by using the in-situ observation of diamond anvil cell with Raman spectroscopy. The dissolution and recrystallization kinetics of sodium sulfate crystals during different temperature and pressure conditions were investigated by in situ observation and spectrometry. The results showed that the Raman spectroscopy of thenardite at room temperature were 449.9, 620.5, 632.9, 647.4, 993.3, 1 101.8, 1 132.2 and 1 153.1 cm-1 respectively. The crystal shape of solid thenardite changed continuously with the slow increase of temperature, and thenardite was dissolved completely when system temperature reached to 193 ℃. The recrystallized crystal appeared with decreasing temperature rapidly, and the new 1 196.5 cm-1 Raman characteristic peak of recrystallized crystal showed the appeared crystal was mirabilite (Na2SO4·10H2O). In-situ observation of diamond showed that thenardite partially dissolveed and recrystallized during the rapid heating process, and the Raman characteristic peak of the recrystallized region was still thenardite. The process of dissolution and crystallization was controlled by diffusion. The Raman spectroscopy improved can be used for quantitative analysis. Compared with the parameter of peak intensity, area and SO42-/H2O intensity ratio in the system solution, the area ratio of SO42-/H2O in the solution reflected the change of SO42- concentration precisely in the solution during the reaction. The SO42-/H2O peak area ratio (AR) is (0.016 6±0.000 4), and the error is 2.4% when solution reaches the dissolved recrystallization equilibrium state. Based on Johnson-Mehl-Avrami-Kolmogorov (JMAK) model and the SO42-/H2O peak area ratio in the solution, the kinetics fitting of dissolution and recrystallization can be simulated. The results showed that the reaction order of dissolution and crystallization of anhydrous sodium sulfate is 1.266 7 and the equilibrium constant of the reaction is 0.001 26 at the temperature of 109 ℃. In summary, the device of hydrothermal diamond anvil cell is simple to operate, and can also avoid errors caused by degeneration and exchange in the quenching process. The advantage of in-situ observation of hydrothermal diamond anvil cell combined with Raman spectroscopy quantitative analysis can be applied to the kinetics of dissolution and crystallization of minerals in aqueous solution under high temperature and high pressure conditions. It is an efficient kinetic research method and is of great significance for studying rapid phase transition.

Keyword: HDAC; High pressure and high temperature; Raman spectroscopy; Dissolution and recrystallization; Kinetics researches
引 言

矿物在水溶液中溶解结晶不仅引起溶液相溶质浓度发生变化, 从溶液中析出晶体的粒度同时存在着时间和空间的分布。 对矿物在水溶液中溶解和结晶及其动力学研究的两种方法, 一是测定溶液中溶质浓度的变化[1, 2], 二是测定固相的变化[3]。 目前高温高压条件下的实验多采用高压釜或活塞圆筒等封闭设备来进行, 封闭设备不能进行原位观测, 降温淬火过程中伴随反应的发生会对样品的真实组分产生影响, 而固液两相在线实时变化信息的准确获取对于溶解结晶过程识别和参数计算具有较大的影响。 相对于传统实验方法的局限性, 水热金刚石压腔可以进行镜下原位观察, 原位拉曼光谱可定量分析溶液中离子的浓度[4, 5, 6]

芒硝为干涸盐湖的化学沉积物, 它从含钠和硫酸根的过饱和盐湖中晶出, 与石膏、 石盐、 无水芒硝、 泻利盐等共生。 应用水热金刚石压腔结合拉曼光谱技术研究无水芒硝在饱和水溶液中的溶解结晶过程, 通过原位观测无水芒硝溶解, 结晶变化, 对该反应中溶液中S O42-浓度变化进行测定, 来探究硫酸钠晶体在不同温度压力条件下的溶解结晶反应机制。

1 实验部分
1.1 装置

高温高压实验采用类似于Mao-Bell金刚石压腔[7, 8], 实验体系的温度由温控器测得并校正, 压力根据石英压标确定[9]。 拉曼光谱应用北京大学地球与空间科学学院Renishow100型激光拉曼光谱仪进行原位测定。

1.2 Na2SO4拉曼光谱

Na2SO4是无水芒硝的主要成分, 无水芒硝为斜方晶系, 晶体结构由[SO4]四面体和[NaO6]八面体组成, 集合体为粒状、 块状或粉末状, 易溶于水, 在潮湿空气中易水化, 逐渐变为粉末状的芒硝。 芒硝为十水硫酸钠 (Na2SO4· 10H2O), 单斜晶系, 晶体结构中[Na(H2O)6]八面体联结成锯齿状链, 链间以[SO4]和2个缓冲H2O分子以氢氧— 氢键相联结, 呈致密块状, 纤维状集合体。 图1显示为常温不同压力调节下无水硫酸钠和含水硫酸钠拉曼光谱图。 常温条件下无水硫酸钠拉曼位移分别为表征[S O42-]基团中S— O对称弯曲振动的449.9 cm-1, S— O变形弯曲振动的620.5, 632.9和647.4 cm-1, S— O对称伸缩振动的993.3 cm-1和S— O反对称伸缩振动1 101.8, 1 132.2和1 153.1 c m-1[10]。 其中S— O对称伸缩振动产生的993.3 cm-1拉曼谱峰峰型好, 强度大, 作为表征硫酸钠晶体的拉曼特征峰。 雷咪等认为1 150 cm-1为芒硝(Na2SO4· H2O)拉曼特征峰, 本实验测试结果显示1 150 cm-1为无水芒硝(Na2SO4)拉曼峰, 1 187 cm-1则为芒硝特征谱峰, 产生这一结果原因可能在于其测定样品为芒硝脱水后产物。 图1对比可得含水硫酸钠在3 396和3 519 cm-1出现O— H振动峰, 以及在1 187 cm-1处[S O42-]基团出现新的S— O反对称伸缩振动峰, 据矿物中基团所呈现的拉曼位移可对硫酸钠固相的成分进行区分。

图1 常温下无水硫酸钠和十水硫酸钠的拉曼光谱图Fig.1 Raman spectra of Na2SO4 and Na2SO4· 10H2O at room temperature

1.3 方法

选择无水Na2SO4晶体和饱和Na2SO4溶液共存态作为初始状态。 体系1研究缓慢升温Na2SO4晶体的溶解重结晶过程; 体系2研究快速升温Na2SO4晶体的溶解重结晶过程。 实验体系1缓慢升温使无水Na2SO4晶体完全溶解在饱和Na2SO4溶液中, 随后降温发生Na2SO4晶体重新结晶; 实验体系2快速升至109 ℃, 随温度升高 Na2SO4溶解度降低导致Na2SO4晶体发生重结晶。 实验过程中应用拉曼光谱仪测定体系溶液中S O42-和H2O拉曼谱峰强度和面积参数, 通过追踪体系溶液中S O42-/H2O拉曼谱峰参数变化趋势, 来探究硫酸钠晶体的溶解结晶反应机制。

2 结果与讨论
2.1 缓慢升温过程中Na2SO4晶体的溶解重结晶

首先在压腔中装入压标物质石英和样品无水硫酸钠晶体, 然后加入饱和Na2SO4溶液迅速密封。 常温条件下放置12 h后进行加温实验, 每次升温10~20 ℃, 升温之后等待10~20 min进行拉曼测试, 原位观测实验过程以及该过程中硫酸钠晶体的拉曼光谱图如图2和图3所示。

图2 不同温度压力下固相Na2SO4的拉曼光谱图Fig.2 Raman spectra of Na2SO4 at different temperature and pressure

图3 Na2SO4溶解结晶过程中固相的变化Fig.3 The change of solid phase in the process of dissolution and recrystallized of Na2SO4

硫酸钠溶液体系中S O42-和H2O的拉曼特征峰位置在980和3 430 cm-1附近, 随温度压力变化, 发生定向移动[6]表1所示为实验体系缓慢升温无水Na2SO4溶解降温重结晶过程溶液中S O42-和H2O的拉曼峰强度和面积参数。

将溶液中S O42-和H2O的拉曼峰强度、 面积及其比值与温度之间的关系作图, 如图4所示。 溶液中S O42-和H2O各自的拉曼峰强度及面积和温度压力之间没有很好的规律性, 但是两者的强度比值IR( ISO42-/ IH2O)在107 ℃之前随着温度的升高而逐渐降低, 在107~168 ℃之间随着温度的升高逐渐升高, 在168 ℃之后随温度升高IR值又逐渐降低, 241 ℃时的IR=0.594 8, 比初始IR比值0.648 4低。 溶液中S O42-和H2O的拉曼峰面积比值AR( ASO42-/ AH2O)在73 ℃之前随着温度的升高而逐渐降低, 之后随温度升高AR值逐渐升高, 241 ℃时的AR=2.65× 10-2, 大于初始AR值1.93× 10-2

表1 Na2SO4溶解重结晶过程中溶液S O42-和H2O的拉曼峰强度和面积参数 Table 1 Raman Intensity and area of S O42- and H2O during the dissolution of Na2SO4

图4 溶液中S$O^{2-}_{4}$和H2O的拉曼峰强度、 面积及其比值随温度变化图
(a): S$O^{2-}_{4}$拉曼峰强度随温度变化趋势图; (b): H2O拉曼峰强度随温度变化趋势图; (c): S$O^{2-}_{4}$拉曼峰面积随温度变化趋势图; (d): H2O拉曼峰面积随温度变化趋势图; (e): S$O^{2-}_{4}$/H2O拉曼峰强度比值随温度变化趋势图; (f): S$O^{2-}_{4}$/H2O拉曼峰面积比值随温度变化趋势图
Fig.4 The intensity and area of S$O^{2-}_{4}$ and H2O in the aqueous solution at different temperature and their ratios
(a): The intensity trend chart of S$O^{2-}_{4}$; (b): The intensity trend chart of H2O; (c): The area trend chart of S$O^{2-}_{4}$; (d): The area trend chart of H2O; (e): The intensity trend chart of S$O^{2-}_{4}$/H2O; (f): The area trend chart of S$O^{2-}_{4}$/H2O

通过和图3中的溶解现象对比, 无水Na2SO4溶解过程中溶液中S O42-和H2O的拉曼峰面积比值更能反映溶液中S O42-浓度的变化, 可作为追踪体系溶解结晶反应的参数进行定量研究。

2.2 快速升温至109 ℃过程中Na2SO4晶体的溶解重结晶

升温过程中, 无水硫酸钠的溶解度减小, 促使固相生长。 本实验常温条件下快速升温至109 ℃, 研究该升温过程中无水硫酸钠的生长速率。 实验过程如图5所示, 固相样品为无水硫酸钠。 图6所示整个实验过程中拉曼光谱未发生突变, 无水硫酸钠无相变。

图5 109 ℃条件下水溶液中固相Na2SO4随时间的变化Fig.5 The change of solid phase of Na2SO4 in the aqueous solution with time at 109 ℃

图6 109 ℃溶液中S O42-和H2O的拉曼峰参数随时间的变化Fig.6 Raman spectra of S O42- and H2O in the aqueous solution at 109 ℃ with time

由图6可以看出溶液中S O42-的拉曼峰位置随着时间逐渐降低, H2O的拉曼峰位置随着时间逐渐升高, S O42-/H2O峰面积比值随着时间呈现出逐渐降低的趋势, 这是由于压腔中的硫酸钠固相逐渐生长造成的。 反应时间1 800 s时溶液中S O42-和H2O的拉曼峰强度随着时间降低趋势变缓, 反应时间2 400 s时趋于稳定, 体系反应达到平衡, 故对体系中固相无水Na2SO4的溶解结晶过程的动力学计算仅采用2 400 s前的数据。 平衡时的峰面积比值AR取2 400 s之后平衡数据平均值(表2所示), 为(0.016 6± 0.000 4), 误差为2.4%。

研究对体系中固相无水Na2SO4的溶解结晶过程选取Johnson-Mehl-Avrami-Kolmogorov(JMAK)[11, 12, 13]模型进行拟合

ln[-ln(1-F)]=nlnk+nlnt

其中, F为反应进度; t为反应时间; n为反应级数; k为反应速率常数。 利用ln[-ln(1-F)]和lnt作图, 根据最小二乘法拟合得到的线性关系的斜率n和截距nlnk可得到反应的级数和反应速率。

根据FA=(AR-AR0)/(AReq-AR0)计算得到反应进度F, 由初始反应的0.126 6到反应平衡的1.113 9(表1所示)。

根据JMAK方程ln[-ln(1-F)]=nlnk+nlnt, 结合表2数据利用ln[-ln(1-F)]对lnt作图(图7所示), 由AR(S O42-/H2O)值得到的线性关系分别为ln[-ln(1-FA)]=1.269lnt-8.4743 R2=0.89可得到反应级数n为1.269, nlnk为-8.474 3, 反应速率常数k为0.0012 6。

图7 表2中数据利用ln[-ln(1-F)]-lnt作图Fig.7 Plot of ln[-ln(1-F)] vs lnt in Table 2

表2 109 ℃溶液中S O42-和H2O的拉曼峰面积和强度参数随时间的变化 Table 2 The intensity and area of Raman shifts of S O42- and H2O in the aqueous solution at 109 ℃ with time
3 结 论

(1)硫酸钠晶体在饱和溶液中重结晶受溶解过程控制。 逐步增温过程中无水硫酸钠晶体完全溶解, 降温结晶十水硫酸钠。 快速增温并保持恒温过程中, 无水硫酸钠晶体部分溶解至过饱和, 再结晶为无水硫酸钠。

(2)根据拉曼光谱定量分析结果, 溶液中S O42-和H2O的面积比值更能反映溶液中S O42-浓度的变化, 硫酸钠在109 ℃条件下溶解重结晶过程实验中反应级数为1.269, 反应平衡常数为0.001 26。

(3)利用金刚石压腔结合拉曼光谱技术可以实现高温高压条件下矿物在水溶液中溶解结晶动力学过程的原位观察和测量, 对于研究快速相变具有重要意义。

参考文献
[1] Bennett P C. Geochimica et Cosmochimica Acta, 1991, 55: 1781. [本文引用:1]
[2] Davis M C, Wesolowski D J, Rosenqvist J, et al. Geochimica et Cosmochimica Acta, 2011, 75: 401. [本文引用:1]
[3] Wang H M, Henderson G S, Brenan J M. Geochimica et Cosmochimica Acta, 2004, 68: 5197. [本文引用:1]
[4] LIU Chuan-jiang, ZHENG Hai-fei(刘川江郑海飞)Earth Science Frontiers(地学前缘), 2012, 19(4): 141. [本文引用:1]
[5] Liu C J, Zheng H F, Du J G, et al. Journal of Crystal Growth, 2015, 410: 39. [本文引用:1]
[6] TIAN Feng, ZHENG Hai-fei, SUN Qiang(田锋郑海飞孙樯)Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2015, 35(4): 924. [本文引用:2]
[7] Mao H K, Bell P M. Carnegie Institute Washington Yearbook, 1978, 77: 904. [本文引用:1]
[8] Bassett W A. European Journal of Mineralogy, 2003, 15: 773. [本文引用:1]
[9] Schmidt C, Ziemann M A. American Mineralogist, 2000, 85: 1725. [本文引用:1]
[10] LEI Mi, CHEN Long, HUANG Bi-sheng, et al(雷咪陈龙黄必胜)China Journal of Traditional Chinese Medicine and Pharmacy(中华中医药杂志), 2016, 31(7): 2811. [本文引用:1]
[11] Barmak K. Metallurgical and Materials Transactions A, 2010, 41A: 2711. [本文引用:1]
[12] Khanna Y P, Taylor T J. Polymer Engineering and Science, 1988, 28: 1042. [本文引用:1]
[13] Fournier J, Regnier E, Faure F, et al. Journal of Non-Crystalline Solids, 2018, 489: 77. [本文引用:1]