冷原子荧光光谱中数字基线估计方法的研究
余鑫, 周伟*, 谢栋材, 肖峰, 李昕雨
成都理工大学核技术与自动化工程学院, 四川 成都 610059
*通讯作者 e-mail: zhouwei@cdut.edu.cn

作者简介: 余鑫, 1998年生,成都理工大学核技术与自动化工程学院硕士研究生 e-mail: 15608187550@163.com

摘要

水体中汞的存在形态有单质汞、 无机汞和有机汞三种。 其中, 甲基汞是主要的有机汞形态, 毒性远高于单质汞和无机汞。 测量水体中的甲基汞的方法有很多, 冷原子荧光光谱法是测量水体中甲基汞的推荐方法。 冷原子荧光光谱法是原子发射光谱法和原子吸收光谱法综合发展而来的元素分析方法。 经过多年的发展与完善, 是元素分析最常用的分析技术之一, 具有灵敏度高、 检出限低等特点, 被广泛应用于环境科学、 生命科学、 地质等领域。 由于检测仪器的背景噪声和色谱柱分离效果的影响, 冷原子荧光光谱出现基线漂移和信号拖尾等干扰因素, 严重影响冷原子荧光光谱数据的峰面积计算和痕量甲基汞的定量分析。 其中, 基线漂移是最主要干扰因素。 目前, 改进模拟器件参数和数字基线估计是解决基线漂移的的两种重要手段。 在改进模拟器件参数方面, 有激发光源使用空心阴极汞灯、 闭环控制的热阴极低压汞灯等, 但存在实验设备复杂、 成本高昂等缺陷; 在数字基线估计方面, 有最小二乘法、 差值拟合法等, 但存在基线估计不稳, 含量计算不准等缺点。 基于此, 提出了一种基于小波变换的数字基线估计方法。 首先, 分析甲基汞的冷原子荧光光谱微观信号和基线漂移现象, 建立冷原子荧光光谱信号和基线漂移数理模型; 其次, 根据冷原子荧光光谱信号模型的特点, 以小波变换为研究基础, 研究合适的母小波模型, 将母小波模型与基线漂移模型进行卷积, 卷积结果恒为零, 理论上证明了基线漂移现象经过小波变换后会被消除; 再次, 以100 pg标样甲基汞为例, 实验验证了小波变换能够有效地消除基线漂移的干扰和信号拖尾的影响; 最后, 在仪器相对标准差(RSD)为1.29%~3.40%的条件下, 对0, 10, 20, 50, 100, 500以及1 000 pg的标样甲基汞溶液进行5次重复实验, 分别建立小波变换前后峰面积平均值校准曲线, 校准曲线的相关系数(R2)由小波变换前的0.994提高到小波变换后的0.997。 实验结果表明, 该方法能够有效地消除测量仪器基线漂移和信号拖尾的影响, 提升了系统测量准确性。

关键词: 甲基汞; 冷原子荧光光谱; 数字基线估计; 小波变换
中图分类号:O433.4 文献标志码:A
The Study of Digital Baseline Estimation in CVAFS
YU Xin, ZHOU Wei*, XIE Dong-cai, XIAO Feng, LI Xin-yu
College of Nuclear Technology and Automation Engineering, Chengdu University of Technology, Chengdu 610059, China
*Corresponding author
Abstract

There are three forms of mercury in water: elemental mercury, inorganic mercury and organic mercury. Methylmercury is the main organic mercury form, and is much higher than that of elemental mercury and inorganic mercury. Cold vapor atomic fluorescence spectrometry(CVAFS) is the recommended method for measuring methylmercury in the water. CVAFS is an element analysis method developed from atomic emission and absorption spectrometry. After years of development and improvement, it has become one of the most commonly used technologies for element analysis. It is widely used in environmental protection, life science, geology and other fields with the characteristics of high sensitivity and low detection limit. However, affected by the background noise of the excitation light source, electronic components of the detection instrument and the separation effect of the chromatographic column, the signal of CVAFS will haveproblems such as baseline drift and signal tailing, which will seriously influence the peak area calculation of the CVAFS’s data and the quantitative analysis of trace methylmercury. Baseline drift is the most critical problem. At present, improving analogue device parameters and digital baseline estimation are two important ways to solve baseline drift. In terms of improving the parameters of the analogue device, there are hollow cathode mercury lamps and closed-loop controlled hot cathode low-pressure mercury lamps with disadvantages such as complex experimental equipment and high cost. The digital baseline estimation includesthe least square method, difference fitting method and so on, as all of them have weaknesses like unstable baseline estimation and inaccurate content calculation. Thus, a digital baseline estimation method based on wavelet transform was proposed. Firstly, by analyzing the microscopic signal of CVAFS and baseline drift of methylmercury, the mathematical model of the signal of CVAFS and baseline drift wasestablished. Secondly, according to the characteristics of the signal of the CVAFS model and wavelet transform, an appropriate mother wavelet model was established. The mother wavelet model was convoluted with the baseline drift model, and the convolution result was always zero. Theoretically, it indicated that the baseline drift wouldbe eliminated after wavelet transform. Thirdly, taking 100 pg standard sample methylmercury as an example, the experiments verified that wavelet transformation could eliminate baseline drift and solve the problem of signal tailing. Finally, under the condition that the relative standard deviation (RSD) of the instrument is 1.29%~3.40%, the experiments were repeated for 5 times for standard methylmercury solutions of 0, 10, 20, 50, 100, 500 and 1 000 pg, and the calibration curves of the average peak area before and after wavelet transform were established respectively. The calibration curve’s correlation coefficient (R2) is increased from 0.994 to 0.997 after the wavelet transform. The experimental results showed that this method could effectively eliminate the influence of baseline drift and signal tailing and improve the system’s measurement accuracy.

Keyword: Methylmercury; CVAFS; Digital baseline estimation; Wavelet transform
引言

甲基汞是水体中一种痕量剧毒的有机物, 当生物体饮用或食用被甲基汞污染的水、 肉类、 蔬菜或粮食等多种食品时, 就会发生中枢神经系统中毒, 造成严重的伤害[1]。 为此, 国家相关部门制定了水质甲基汞的标准[2]。 标准中推荐的方法就是冷原子荧光光谱法。 该方法具有灵敏度高、 检出限低等优点, 非常适合水质甲基汞的测量。 但是, 在测量仪器研制过程中, 电子元件、 电路设计以及处理方法等方面不足, 导致了基线偏移、 信号拖尾等技术瓶颈。 其中, 稳定基线是最关键的技术。

现有的仪器稳定基线一般采用改进模拟器件参数的方法, 如激发光源使用空心阴极汞灯[3], 以及闭环控制的热阴极低压汞灯[4]。 改进模拟器件参数的方法存在实验设备复杂, 成本高昂等缺陷, 一般由实验室自行搭建, 推广难度相对较大, 并且, 由于模拟器件的属性, 基线的稳定能力有限, 故有必要对数字基线进行处理。 目前没有发现国内外对冷原子荧光光谱进行数字基线估计的相关研究。 传统的数字基线估计的算法主要有最小二乘法[5]、 差值拟合法[6]、 平均阈值法[7]、 卡尔曼滤波法[8]等。 这些方法普遍应用在核信号数字分析中, 在一定程度上解决了基线估计的问题, 但存在估计基线不够准确等问题。

基于上述问题, 本文提出采用小波变换对甲基汞数字信号进行基线处理的方法。 建立母小波与甲基汞数字信号的小波变换数理模型, 验证了该方法的可行性, 通过不同浓度的标样甲基汞实验, 建立校准曲线, 实验结果表明该方法具有估计并自动扣除基线、 信号光滑去噪和准确性高等优势。

1 实验部分
1.1 检测系统

搭建甲基汞检测系统, 该系统如图1所示, 由吹扫捕集模块、 气相色谱模块、 原子化器模块、 冷原子荧光检测模块构成。 待衍生剂将样品中的甲基汞转化为易挥发的甲基汞衍生物后, 吹扫捕集模块将甲基汞衍生物从水样中吹扫分离, 经由捕集管进行富集; 气相色谱模块将不同形态的汞衍生物按时间分离; 原子化器模块是将甲基汞衍生物裂解为汞原子; 冷原子荧光检测模块将汞原子激发后, 采集汞的特征荧光信号并转换为电信号分析处理[9, 10, 11]

图1 甲基汞检测系统实物图Fig.1 Physical map of methylmercury detection system

重复实验发现, 由于激发光源的光强不是常量, 存在波动, 并且在电子学元器件温漂、 工频干扰以及白噪声等因素作用下, 基线会漂移。 基线漂移实测图如图2所示。

图2 系统基线漂移Fig.2 System baseline drift

1.2 荧光信号数理模型

在理想情况下, 多个汞原子可视为一个单位的荧光标准微球, 它的荧光强度与光电倍增管输出电压有一定的特性关系。 图3为基于冷原子荧光光谱的电压信号微观原理图, 横坐标为荧光标准微球进入激发区域的时间, 纵坐标为光电倍增管输出的电压信号。 将进入荧光池的荧光标准微球作为研究对象, 当荧光标准微球逐渐进入激发区域时, 外层一部分电子被激发, 在退激发时产生特征荧光被光电倍增管检测到并输出电压信号; 随着荧光标准微球进入激发区域的体积逐渐增大, 被激发的外层电子逐渐变多, 退激发产生的荧光信号逐渐增强, 光电倍增管的输出电压也逐渐变大; 当整个荧光标准微球进入到激发区域时, 被激发的外围电子数量达到最大, 产生的荧光信号也最强, 光电倍增管输出的电压达到峰值; 当荧光标准微球逐渐离开激发区域时, 被激发的外层电子逐渐减少, 荧光信号逐渐变弱, 电压信号逐渐降低。

图3 基于冷原子荧光光谱的电压信号微观原理图Fig.3 Microscopic schematic diagram of voltage signal based on CVAFS

因此, 理想情况, 汞单质在进入和远离激发区域时的电压信号关于时间的函数呈现高斯分布, 可以表示为

ν(t)=Ae-(t-tp)22δ2(1)

式(1)中A是电压信号的幅值, tp是达峰时间, δ 是电压信号的方差。

峰面积H可表示为

H=-+ν(t)dt=2πδA(2)

根据高斯信号的傅里叶变换特性, 傅里叶变换时移特性和频率标度可以推导出式(1)的傅里叶变换表达式

V(ω)=2πδAe-(δ2ω22+tp)(3)

1.3 小波变换推导

给定高斯函数[12]

h0(t)=-e-t22(4)

对式(4)二次求导, 再进行尺度变换和时移可以得到式(5)

h(t)=1μ1-(t-tc)22μ2e-(t-tc)22μ2(5)

式(5)中, μ是尺度因子, tc是中心时间。

式(5)满足

Ψ(0)=-+h(t)dt=0(6)

即式(5)满足母小波的容许条件[13], 也就是说, 式(5)可作小波变换的母小波。

式(5)的傅里叶变换表达式

H(ω)=2πμω2e-(μ2ω22+tc)(7)

式(3)与式(7)相乘

Y(ω)=V(ω)H(ω)=2πδAμω2e-(ω2(δ2+μ2)2+(tp+tc))(8)

对式(8)傅里叶反变换得到式(1)和式(4)的卷积形式

y(t)=2πδAμδ2+μ21-(t-tp-tc)2δ2+μ2e-(t-tp-tc)22(δ2+μ2)(9)

对于理想信号, 基线恒为零, 只需要提取信号幅值以及计算出方差, 就可以得到峰面积。 但是根据图2, 原始信号存在基线漂移。

不妨设基线漂移的表达式为

n(t)=kt+a(10)

式(5)与式(10)卷积

h(t)* n(t)=-+h(τ)n(t-τ)dτ=-+1μ1-(τ-tc)2μ2e-(τ-tc)22μ2(k(t-τ)+a)dτ=2π(kt+a+ktc)-2π|μ|2(kt+a+ktc)μ2=0(11)

由式(11)可知, 基线漂移现象经过小波变换后会被消除。

母小波式(5)满足容许条件, 由式(6)可得

-+h(t)dt(h(t)0)=--+h(t)dt(h(t)< 0)(12)

即母小波的时间轴上半部分面积等于下半部分面积, 同理, 式(9)也满足式(6)。

式(9)的离散形式可表示为

y(n)=2πδAμδ2+μ21-(n-tp-tc)2δ2+μ2·  e-(n-tp-tc)22(δ2+μ2)nN(N为自然数)0others(13)

因此, 在进行数字信号处理时, 基线已被扣除, 只需对式(13)正数部分累加

M=n=-+y(n) (y(n)> 0)(14)

可得到小波变换后的峰面积M

从上述推导可以看出, 该方法能有效地估计并自动扣除基线, 简化了计算峰面积的工作, 提高了峰面积的准确性。 理论上验证了该方法的可行性以及有效性。

2 结果与讨论

水体中含有不同形态的汞元素, 冷原子荧光光谱法测量甲基汞, 需要和色谱柱联用, 利用不同形态汞元素在色谱柱中的保留时间不同的特性, 可将甲基汞按时间差分离出来。 但是色谱柱分离效果并不能达到100%[14], 故实测的信号不仅存在基线漂移, 还会有拖尾的现象。

以100 pg标样甲基汞为例。 100 pg标样甲基汞单峰信号如图4所示。 可以看到, 约在1 200~1 600 s段存在拖尾的现象, 并且, 总体存在基线漂移现象, 信号不够平滑, 毛刺较多, 实测信号难以估计其基线。

图4 100 pg标样甲基汞单峰信号Fig.4 Single peak signal of 100 pg standard methylmercury

根据式(9)和式(11)的推论, 实测信号达峰时间tp在500~1 100 s之间。 多次实验发现, 式(5)母小波选取中心时间tc=1 000 s, 尺度因子μ=137.7时, 得到的小波变换波形最优。 母小波的时域图如图5所示。

图5 母小波时域图Fig.5 Mother wavelet time domain diagram

依据式(9), 100 pg标样甲基汞单峰信号与母小波进行卷积后得到图6。

图6 100 pg标样甲基汞单峰与母小波卷积结果Fig.6 100 pg standard methylmercury single peak and mother wavelet convolution result

由图6可看出, 原始信号经过小波变换后, 基线漂移被消除且自动扣除基线, 式(11)得证; 拖尾现象影响波形的对称性, 但不影响峰面积的计算, 同时, 该母小波对实测信号有光滑的作用, 消除了毛刺, 利于算出准确的峰面积。

在实验室条件下, 分别对0, 10, 20, 50, 100, 500以及1 000 pg的标样甲基汞溶液进行5次实验, 并计算峰面积。 结果如表1所示。 相对标准差(RSD)在1.29%~3.40%之间, 仪器精密度良好。

表1 不同含量标样甲基汞峰面积 Table 1 Peak area of methylmercury with different content standards

依次对5次重复实验的结果根据式(14)对变换的结果进行计算, 分别建立小波变换前后峰面积平均值校准曲线, 两条校准曲线的相关系数和表达式如图7所示。

图7 小波变换前后峰面积校准曲线Fig.7 Peak area calibration curve before and after wavelet transform

由图7可知, 标样甲基汞的峰面积校准曲线的相关系数(R2)由小波变换前的0.994提高到小波变换后的0.997, 说明小波变换能有效提高系统测量的准确性。 小波变换前后仪器相关系数情况如表2所示。

表2 小波变换前后相关系数对比 Table 2 Correlation coefficient comparison before and after wavelet transform
3 结论

在冷原子荧光光谱分析系统中引入小波分析方法进行基线的数字估计是一次成功的尝试。 系统中, 光谱信号的高斯模型建立, 母小波函数和尺度的选取是整个研究的关键技术要点。 得益于小波分析方法在多分辨率、 相对带宽恒定以及快速的奇异点检测等方面的优势, 系统实现荧光光谱基线的估计、 扣除以及浓度计算等数字化处理, 解决了系统基线漂移和信号拖尾等问题, 提升了系统分析的准确性和稳定性。 在本研究的基础上, 继续探讨提升前端激发光源稳定性的技术方法, 将更有利于进一步提升冷原子荧光光谱分析系统的综合技术指标, 这也是下一步需要继续研究的方向。

参考文献
[1] Wang Xinjin, Yan Mengling, Zhao Lina, et al. Oxidative Medicine and Cellular Longevity, 2016, 2016: 5137042. [本文引用:1]
[2] HJ 977-2018. Water Quality—Determination of Alkylmercury—Purge and Trap/Gas Chromatography Cold Vapor Atomic Fluorescence Spectrometry(水质 烷基汞的测定吹扫捕集/气相色谱-冷原子荧光光谱法). Ministry of Ecology and Environment of the People’s Republic of China(中华人民共和国生态环境部). Beijing: China Environmental Press(北京: 中国环境出版社), 2018. 11. [本文引用:1]
[3] WANG Xue-feng, LI Xiang, KANG Wen-gui(王雪枫, 李湘, 康文贵). Physical Testing and Chemical Analysis(Part B: Chemical Analysis)(理化检验-化学分册), 2020, 56(3): 344. [本文引用:1]
[4] QIN Yu-qian, DUAN Fa-jie, FU Xiao, et al(秦玉倩, 段发阶, 傅骁, ). Journal of Chongqing University(重庆大学学报), 2021, 44(5): 87. [本文引用:1]
[5] LI Yan-bo, LI Chao-hai, HE Zi-shu, et al(李彦波, 李朝海, 何子述, ). Nuclear Electronics & Detection Technology(核电子学与探测技术), 2014, 34(8): 991. [本文引用:1]
[6] WANG Xin, FAN Xian-guang, XU Ying-jie, et al(王昕, 范贤光, 许英杰, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2014, 34(8): 2117. [本文引用:1]
[7] CHEN Wei, WAN Wen-jie, YU Jie, et al(陈伟, 万文杰, 喻杰, ). Nuclear Techniques(核技术), 2019, 42(4): 58. [本文引用:1]
[8] Grover A, Lall B. IEEE Transactions on Instrumentation and Measurement, 2020, 69(9): 7306. [本文引用:1]
[9] He Qian, Yu Xueke, Li Yanbin, et al. Microchemical Journal, 2018, 141: 148. [本文引用:1]
[10] Kerstin Leopold, Michael Foulkes, Paul Worsfold. Analytica Chimica Acta, 2010, 663(2): 127. [本文引用:1]
[11] Jia Xiaoyu, Gong Dirong, Han Yi, et al. Talanta, 2012, 88: 724. [本文引用:1]
[12] Loska L, Janczyszyn J. Applied Radiation and Isotopes, 1997, 48(1): 127. [本文引用:1]
[13] Qin Zhangjian, Chen Chuan, Luo Junsong, et al. Nuclear Inst. and Methods in Physics Research, A, 2017, 858: 22. [本文引用:1]
[14] Arakawa S, Furuno M, Fukusaki E. Journal of Bioscience and Bioengineering, 2021, 132(1): 102. [本文引用:1]