LIBS技术的沉积物中重金属元素Zn的分析研究
张超, 朱林, 郭金家*, 李楠, 田野, 郑荣儿
中国海洋大学信息科学与工程学院, 山东 青岛 266100
*通讯联系人 e-mail: opticsc@ouc.edu.cn

作者简介: 张 超, 1995年生, 中国海洋大学信息科学与工程学院硕士研究生 e-mail: zhangc_ouc@163.com

摘要

随着经济的发展, 海洋、 河流污染日益严重, 沉积物是判断海洋、 河流重金属污染程度的重要窗口, 实际调查中也迫切需要能对沉积物污染程度快速检测的技术。 激光诱导击穿光谱技术(LIBS)具有实时、 样品预处理简单、 多种元素同时检测等优点, 是一种可用于在线分析沉积物中重金属元素的有效手段。 应用实验室研发的便携式LIBS装置开展了沉积物中Zn元素的检测, 该装置采用Nd:YAG脉冲激光器作为激发光源, 用单通道微型光纤光谱仪进行分光检测, 装置整体尺寸为371 mm×192 mm×294 mm, 可方便的进行现场测量。 为实现LIBS技术对沉积物的定量分析, 首先在实验室对沉积物的LIBS数据进行预处理, 对八种沉积物样品进行干燥、 过筛、 混匀和压片处理后进行测量, 通过对比采用四种不同标准化方法处理后光谱的相对标准偏差(RSD), 确定最佳光谱标准化方法为结合等离子体特征的标准化, 采用该方法后的RSD降至未标准化处理光谱的五分之一; 之后分别采用单变量分析方法和带动量因子的自适应学习率三层BPNN算法进行定量分析, BPNN算法的训练相关系数R2达到0.99, 预测误差为3.2%, 相比单变量定量分析方法预测结果改善了一倍。 同时, 针对现场快速分析, 对比了不同干燥时间下沉积物光谱强度变化, 12 min后信号强度随干燥时间的变化不大, 因此确定现场探测对沉积物的干燥时间为12 min。 利用该装置在胶州湾海域六个不同站位进行了现场沉积物采集和测量试验, 在甲板上对现场采集的沉积物样品快速烘干后进行LIBS探测, 并采用上述实验室建立的定量方法对Zn元素进行浓度预测, 与带回实验室采用原子吸收分光光度法(AAS)的检测结果相比, LIBS探测结果的平均误差为8.32%。 验证了该装置应用于现场沉积物快速测量的可行性, 同时证明了结合等离子体特征的标准化方法以及BPNN用于沉积物中Zn元素浓度预测的适用性。

关键词: 激光诱导击穿光谱(LIBS); 便携式装置; 沉积物; 重金属元素Zn; 定量分析
中图分类号:O433.4 文献标志码:A
Laser-Induced Breakdown Spectroscopy for Heavy Metal Analysis of Zn of Ocean Sediments
ZHANG Chao, ZHU Lin, GUO Jin-jia*, LI Nan, TIAN Ye, ZHENG Rong-er
College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China
*Corresponding author
Abstract

With the rapid development of economy, the intensely deepening pollutions in the ocean and rivers have become a major issue. As a key indicator of environmental assessment, the sediments can offer important information on heavy metal pollution. Infield investigation, rapid and real-time detection method for sediments is lacking for the time being. With the advantages of fast and real-time analysis, simple sample pretreatment, multiple element detection, Laser-induced Breakdown Spectroscopy (LIBS) can meet the requirements of field sediment probe. In this paper, a portable LIBS system was developed and used for the element of Zn detection from sediment samples.The system was equipped with a Nd:YAG pulses laser as the source and a miniature optical fiber spectrometer. With a size of 371 mm×192 mm×294 mm, the system was convenient for field measurement. To quantify the Zn in sediment samples, the LIBS spectra of sediments obtained after the samples pretreatment shouldbe standardized in the first step. The best method was standardized with characteristic plasma parameters through a comparison with some typical methods to reduce the relative standard deviation (RSD) to its fifth of original data. A backpropagation artificial neural network (BPNN) with the driving factorwas adopted as the calibration model for quantitation. An R2 value over 0.99 and a twice better predictive result were achieved with the single variable analysis. Also, we found that the spectral intensity of sediments was not significantly affected by drying longer than 12 minutes. Therefore, LIBS detection can be carried out after a 12-minute drying to the samples in the field investigation. The field experiment was made at 6 different sampling sites in Jiaozhou Bay using the portable LIBS system. A LIBS detection can be carried out after a 12-minute drying to the samples in the field investigation. The concentration prediction of Zn was got by using of the above methods. The results showed their coincidencein comparison with the results using Atomic Absorption Spectroscopy (AAS) method with the average relative deviation of 8.32%. The obtained results proved the performance of the portable LIBS system in field sediment sample analysis and showed the applicability for predicting Zn concentration using the standardization with characteristic plasma parameters and BPNN combined approach.

Keyword: Laser-induced breakdown spectroscopy (LIBS); Portable system; Sediments; Heavy metal Zn; Quantitative analysis
引言

随着经济的发展, 运输、 采矿等工业排污与人类生活排污等活动造成了水体污染, 水中的金属污染不断加重, 经过一定后污染物最终进入沉积物中[1, 2]。 常用于沉积物中重金属测量的方法很多, 包括传统的实验室方法, 如: 原子吸收光谱法(atomic absorption spectroscopy, AAS)[3]、 电感耦合等离子体质谱法(inductively coupled plasma mass spectrometry, ICP-MS)[4], 以及可用于现场探测的激光诱导击穿光谱(laser induced breakdown spectroscopy, LIBS)[5]等。

LIBS作为一种实时快速的光谱探测技术, 具有多组分同时探测、 无需或只需简单样品预处理等突出优势, 已被广泛应用到工业[6]、 环境[7]、 地质[8]、 海洋[9]等领域。 在对沉积物这种成分复杂、 含水量高的样品探测分析中, 样品的预处理不可或缺, 因为较大的形态差异和较高的含水量会严重影响LIBS探测的稳定性和灵敏度[10]。 Dukki Han等对北冰洋的沉积物样品进行了冷冻干燥、 过筛和研磨压片处理, 应用LIBS技术获得的样品中Al和Cu等金属元素含量与ICP-MS检测结果相近[5]。 Lazic等采用单色仪和增强型CCD探测器搭建了一套船基LIBS系统, 应用于南极科考调查中的海洋沉积物的探测分析[11], 他们对样品进行烘干48 h后探测。 Elmer等使用J200(ASI)系统对河流沉积物进行LIBS探测, 该系统采用了高性能ICCD探测器配置, 样品经过48 h干燥后再加入粘合剂压片处理, LIBS系统对其中的Cr, Cu和Pb含量预测结果与传统的ICP-OES方法测量结果有很好的一致性[12]。 然而, 高性能配置的系统体积大、 样品脱水烘干时间长等特点使LIBS技术用于沉积物的现场快速检测分析受到了一定限制。

针对沉积物的现场快速检测需求, 搭建了一套基于微型光纤光谱仪的便携式LIBS装置, 采用标准化方法以改善LIBS数据质量, 通过不同定量分析方法的比对以寻求预测沉积物中Zn元素浓度的最佳方案, 并通过现场试验进行评估。

1 实验部分
1.1 基于微型光纤光谱仪的LIBS装置

便携式LIBS装置如图1所示, 整体尺寸为371 mm× 192 mm× 294 mm, 重量为10 kg。 激发光源选用1 064 nm半导体泵浦固体激光器(MONTFORT, M‐NANO), 单脉冲能量最高为40 mJ, 频率为10 Hz。 激光经过透镜L1和L2(f1=-38.1 mm, f2=75 mm)扩束, 由反射镜M1反射后透过二向色镜M2由组合透镜L3(f3=25.4 mm)聚焦至样品表面。 LIBS信号原路返回, 经二向色镜M2反射, 通过透镜L4(f4=50.8 mm)耦合进入微型光纤光谱仪(Avantes, AvaSpec-ULS2048)。 光谱仪的光谱探测范围是300~550 nm, 探测延时最低为1.28 μ s, 光谱分辨率为0.29 nm。 实验样品置于旋转平台, 保证相邻击穿点之间没有重叠区域。 我们选择Ca Ⅰ 422 nm和Zn Ⅰ 468 nm的光谱线对探测延时和激光能量进行评估, 实验结果表明最佳实验条件为激光能量40 mJ、 探测延时为1.28 μ s。 对每个样本进行100次光谱测量, 每次光谱数据为5个激光脉冲激发光谱的平均, 采集门宽为1.05 ms。

图1 便携式LIBS装置图
(a): 系统光路图; (b): 系统内部结构图
Fig.1 Schematic diagram of the portable LIBS instrument
(a): Optical layout; (b): Internal structure charts of the instrument

1.2 样品

便携式LIBS装置进行标定的实验样品来源于小清河流域八个不同地点的沉积物, 其特点为颗粒大小不一、 杂质众多, 含水率较高, 其中主要金属元素包括Ca, Fe和Mn等。 将样品一分为二, 分别采用AAS和LIBS检测Zn的浓度。 AAS检测结果如表1所示。 LIBS实验中统一对样品进行干燥、 过筛、 混匀和压片处理, 处理后的样品为直径13 mm、 重300 mg的圆形薄片。 沉积物的典型LIBS光谱如图2所示。

图2 沉积物典型LIBS光谱图Fig.2 A typical sediment LIBS spectrum

表1 实验样品中Zn浓度的AAS检测结果 Table 1 AAS detection results of the concentration of Zn in experimental samples
2 结果与讨论

图3为本文设计的基于便携式LIBS装置的沉积物中Zn元素浓度预测流程图, 通过对比建立最佳定量分析方法, 实现对沉积物中Zn元素浓度的准确预测。

图3 基于便携式LIBS装置的沉积物中Zn元素浓度预测流程图Fig.3 Flow chart of concentration prediction of Zn in sediment using portable LIBS system

2.1 沉积物光谱标准化方案选择

在LIBS探测实验中, 激光通量不稳定、 基体效应等因素都会导致光谱信号的波动, 从而影响激光烧蚀量与光谱强度的相关性[13], 进而影响定量分析的准确性, 因此需要对LIBS光谱进行标准化处理。 针对沉积物中的重金属元素Zn, 分别进行了min-max标准化、 Z-score标准化、 全谱标准化以及结合等离子体特征的标准化。 结合等离子体特征的标准化主要利用标准等离子体温度、 电子数密度和目标元素的电子数密度, 对目标元素的谱线强度进行修正[13], 公式为

IZn=a1I468nm+a2(Isum)+a51-a3(Isum-Imean)-a4(FWHMCa-FWHMmean)

其中a1, a2, a3, a4, a5为拟合系数, Isum为Zn元素多条谱线和, Imean为谱线强度平均值, FWHMmean为谱线半高全宽平均值。

图4给出了不同标准化方法下Zn Ⅰ 468 nm特征峰的相对标准偏差(relative standard deviation, RSD)结果, 其中结合等离子体特征的标准化方法的RSD优于其他方法, 平均RSD由38.7%降至7.2%, 提高了沉积物LIBS光谱的可重复性。

图4 不同标准化方法Zn元素的RSD对比图Fig.4 The RSD results of Zn with different standardized methods

2.2 沉积物光谱定标方案选择

图2表明Zn的三条特征谱线468.01, 472.22和481.05 nm激发上能级相同, 因此使用三条特征峰的谱峰面积之和建立定标曲线, 并对X10#样品的Zn浓度进行预测。 图5(a)对比了对光谱进行结合等离子体特征的标准化处理前后的定标曲线, 可以看出, 标准化对定标结果有一定的改善, 但单变量定量分析拟合曲线线性相关度不高, 说明单纯的线性关系不能很好地减弱基体效应的影响。 故采用带动量因子的自适应学习率的三层BPNN(back propagation neural network)算法, 将标准化后的光谱作为神经网络输入, 选择的隐藏神经元个数为150~170, 迭代次数为5 000次, 学习速率为0.01, 动力参数为0.85~0.95, 输出层神经元个数为1, 即为所预测元素浓度。 训练结果如图5(b)所示。

图5 沉积物中Zn定标曲线
(a): 单变量方法对沉积物中Zn的定标曲线; (b): BPNN算法对沉积物中Zn的定标曲线
Fig.5 Calibration curve for Zn in sediment
(a): Calibration curve for Zn with univariate quantitative method; (b): Calibration curve for Zn with BP neural network

三种定标方法的定量分析结果对比见表2。 采用等离子体特征的标准化后光谱的重复性有提高, 定量预测的误差由10.5%降低至6.7%; BPNN算法的定标曲线的相关系数R2提高至0.998, 预测相对误差为3.2%, 相较于单变量定量分析方法改善了一倍。 由此可见, 采用结合等离子体特征的标准化和BPNN算法的定量效果最佳, 从光谱中提取的有效信息更多, 可以更大程度地校正基体效应, 以提高定量分析的精度。

表2 三种定标方法结果 Table 2 Quantitative result of three analysis methods
2.3 样品烘干时间选择

沉积物中水分的存在会大大降低LIBS烧蚀效率进而影响信号强度[14], 而长时间脱水处理又限制了现场检测的快速性。 为了确定合适现场探测的样品烘干时间, 探究了烘干时间对LIBS信号的影响, 实验结果如图6所示。 在烘干时间低于12 min时, LIBS信号强度增加迅速, 12 min后LIBS信号随烘干时间的变化不明显。 因此在现场试验中选择沉积物样品的烘干时间为12 min。

图6 Ca和Zn的谱线强度随样品烘干时间的变化曲线Fig.6 The emission intensities of Ca Ⅰ 422 nm and Zn Ⅰ 468 nm as a function of sample’ s drying time

2.4 现场试验

2019年3月我们应用便携式LIBS装置在胶州湾海域进行了现场探测试验, 船基试验现场照片和试验站点位置分别如图7(a)和图7(b)所示。 所选择的站点包括3个河口站点(海泊河口S1、 李村河口S2、 大沽河口S6)和1个近岸站点(S3)和2个离岸较远站点(S4和S5)。 图7(c)展示的是利用沉积物取样器获得的试验站位沉积物原样和经过制备后的分析样品图片。 在每个取样站位获得的沉积物样品均一分为二, 分别用于现场LIBS探测和实验室AAS比对检测。 现场试验中将采集的沉积物置于直径为3 cm的圆盘内, 去除明显的颗粒物后将表面抹平, 然后烘干12 min达到表面干燥状态, 便可进行LIBS现场探测。 图7(d)是沉积物现场LIBS探测获得的典型光谱图。 利用2.2建立的结合等离子体的标准化与BPNN算法对各站位沉积物样品(J01#— J06#)中重金属Zn元素含量进行预测, 预测结果以及实验室AAS的比对检测结果如表3所示。

图7 便携式LIBS装置胶州湾现场试验照片与站点位置图
(a): 便携式LIBS装置现场试验照片; (b): 现场沉积物采样站点图; (c): 沉积物样品与现场制备照片; (d): 现场典型光谱图
Fig.7 Field measurements pictures and site positions in Jiaozhou Bay
(a): Field measurement picture of the portable LIBS system; (b): Site map for sediments sampling; (c): Picture of a typical sediment sample; (d): A typical sediment LIBS spectrum on site

表3 现场LIBS探测结果与实验室AAS测量结果对比 Table 3 Comparisons of LIBS measurements on site and AAS results in lab

表3中可以看出, LIBS预测结果与AAS检测结果有较好的一致性。 两者均表明两个河口站位样品(J01#, J02#)中Zn元素含量明显高于其他区域; 大沽河口站位样品(J06#)的Zn含量的明显偏低。 前者是由于河流上游工厂居多, 从而加重了河口区的重金属污染所致。 后者是因为大沽河口站位的沉积物属于粉砂质黏土, 砂砾对重金属的吸附作用较弱[15], 从而造成沉积物中金属元素含量较低。 从离岸较远的两个站位样品(J04#, J05#)中测得的Zn含量预测值均低于离岸较近的站位样品(J03#)中的预测结果, 这也与常理认知相吻合。 以AAS方法为标准, LIBS测量结果平均相对误差8.32%, 最大相对误差12.77%, LIBS测量结果的平均均方根误差(root mean square error, RMSE)为6.92, 最大为9.36。 对比结果证明了所搭建的LIBS装置及建立的定量分析方法有望成为用于沉积物成分现场快速分析的有效手段。

3 结论

针对沉积物现场快速检测的需求, 使用便携式LIBS装置在实验室中利用样品X01#— X07#建立定标曲线, 对比三种不同的定量方法得出浓度预测最佳方案: 通过结合等离子体特征的标准化方法减小光谱波动, 与原始光谱相比, RSD由38%下降至7%; 利用带动量因子的自适应学习率的三层BPNN算法对沉积物中的Zn元素进行浓度预测, 神经网络训练的相关系数达到0.99以上, 预测误差为3.2%。 为了实现在现场的快速检测样品, 烘干时间选择12 min。 应用该系统在胶州湾开展现场试验, 结果表明, 河口处的沉积物Zn元素浓度整体高于远离河口区域, 符合实际情况; LIBS与AAS两者所测得的浓度具有较好的一致性, 平均误差为8.3%。 证明了本定量方法的有效性, 也验证了LIBS技术用于沉积物现场快速测量的可行性。

参考文献
[1] Baturin G N. Doklady Earth Sciences, 2012, 444(2): 729. [本文引用:1]
[2] Meng M, Shi J B, Yun Z J, et al. Marine Pollution Bulletin, 2014, 88(1-2): 347. [本文引用:1]
[3] DONG Wei-feng, BAO Mu-tai, CAI Ying, et al(董炜峰, 包木太, 蔡莹, ). Marine Environmental Science(海洋环境科学), 2013, 32(1): 137. [本文引用:1]
[4] Bastami K D, Bagheri H, Haghparast S, et al. Marine Pollution Bulletin, 2012, 64(12): 2877. [本文引用:1]
[5] Han D, Joe Y J, Ryu J S, et al. Spectrochimica Acta Part B: Atomic Spectroscopy, 2018, 146: 84. [本文引用:2]
[6] Boué-Bigne F. Spectrochimica Acta Part B: Atomic Spectroscopy, 2014, 96: 21. [本文引用:1]
[7] Pand hija S, Rai N K, Rai A K, et al. Applied Physics B: Lasers and Optics, 2010, 98(1): 231. [本文引用:1]
[8] Khajehzadeh N, Kauppinen T K. IFAC-PapersOnLine, 2015, 48(17): 119. [本文引用:1]
[9] Guo J J, Lu Y, Cheng K, et al. Applied Optics, 2017, 56(29): 8196. [本文引用:1]
[10] LIU Jin-qing, ZHANG Yong, YIN Ping, et al(刘金庆, 张勇, 印萍, ). Marine Geology & Quternary Geology(海洋地质与第四纪地质), 2016, 36(1): 137. [本文引用:1]
[11] Barbini R, Colao F, Lazic V, et al. Spectrochimica Acta Part B: Atomic Spectroscopy, 2002, 57(7): 1203. [本文引用:1]
[12] Austria E S, Nuesca G M, Lamorena R B. Environmental Science and Pollution Research, 2018, 25(2): 1. [本文引用:1]
[13] Li L Z, Wang Z, Yuan T B, et al. Journal of Analytical Atomic Spectrometry, 2011, 26(11): 2274. [本文引用:2]
[14] Michel A P M, Sonnichsen F. Spectrochimica Acta Part B: Atomic Spectroscopy, 2016, 125: 177. [本文引用:1]
[15] CHU Ya, LIU Song-yu, CAI Guo-jun, et al(储亚, 刘松玉, 蔡国军, ). Chinese Journal of Underground Space and Engineering(地下空间与工程学报), 2014, 10(6): 1312. [本文引用:1]