激光与光电子学进展, 2018, 55 (1): 013004, 网络出版: 2018-09-10   

基于PARAFAC和ART算法的油类污染物荧光检测 下载: 1065次

Fluorescence Detection of Oil Pollutants Based on PARAFAC and ART Algorithms
作者单位
1 华北理工大学电气工程学院, 河北 唐山 063210
2 燕山大学测试计量技术及仪器河北省重点实验室, 河北 秦皇岛 066004
摘要
采用平行因子分析(PARAFAC)和交替残差三线性(ART)算法,对石油类污染物进行测量与识别,重点对比分析了两种算法对油种鉴别的差异。在实验中,将以CCl4为溶剂的95号汽油、0号柴油与普通煤油溶液作为研究对象,以不同浓度的石油类物质混合液作为实验样本,利用F-7000荧光分光光度计对样本进行检测,以得到各样本的三维荧光数据。测量样本的组分数估计值预设为3时,采用PARAFAC算法得到的柴油、汽油和煤油样品的回收率分别为(95.60±3.60)%、(94.67±3.66)%和(95.49±4.49)%;ART算法无需预设组分数,其测量得到的柴油、汽油和煤油样本的回收率分别为(96.58±2.17)%、 (95.17±9.17)%和(95.90±8.90)%。结果表明:两种算法都可用于三组分石油类污染物的识别与测量,均能得到较高的回收率;ART算法因无需预先设定组分数而更具优势。
Abstract
Parallel factor analysis (PARAFAC) and alternating residual tri-linearization (ART) algorithms are used to measure and identify petroleum pollutants. The differences between the two algorithms in oil identification are emphatically compared and analyzed. The CCl4 solutions of No. 95 gasoline, No. 0 diesel and kerosene are used as the research objects. We take petroleum mixed solutions with different concentrations as samples to measure the three-dimensional fluorescence data of each sample by F-7000 fluorescence spectrometer. When PARAFAC algorithm is applied and the component number is set to 3, the recovery rates of diesel, gasoline and kerosene are (95.60±3.60)%, (94.67±3.66)% and (95.49±4.49)%, respectively. ART algorithm does not require a preset component number, and the recovery rates of diesel, gasoline and kerosene are (96.58±2.17)%, (95.17±9.17)% and (95.90±8.90)%, respectively. The results show that the two algorithms can be used for the measurement and identification of three kinds of petroleum pollutants, and high recovery rates can be obtained. ART algorithm does not require presetting component number, so it has more advantages.

1 引言

近年来,人类在生产活动中对海洋造成的污染危及动植物的生存,并间接损害人体健康。水污染的治理显得尤为迫切。由于海上石油泄漏污染事故频发,在水污染治理中,油类污染物的防治是重中之重[1-4]

原油及其炼制而成的石油产品中含有大量多环芳烃,具有较强的荧光特性,因此采用荧光分析法对石油类污染物进行测定与识别较为简便。在荧光分析法中,三维荧光光谱法因快速、灵敏的检测优点而在监测领域得到广泛应用[5-9]。平行因子分析(PARAFAC)算法专门用于解决三线性数据问题,较多地应用在化学分析方面[10-11]。近几年,有多位学者对此进行了大量实验研究[12-15]。PARAFAC算法需要预设组分数,基于此,杜文[16]提出了交替残差三线性(ART)算法,实现了无需预设组分数的三线性数据分析。本文采用上述两种算法,分别对三组分石油类污染物进行定性识别与定量测定,通过实验结果对比两种算法的优点[17-20]

2 基本原理

2.1 PARAFAC算法

PARAFAC算法广泛应用于三维及多维数据的分析中。三线性数据由某些三维荧光数据组成[21]。PARAFAC算法的最大优点是其分解结果唯一,但该算法对成分数敏感,需要预先设定组分数。一个三线性数据阵X的PARAFAC算法模型如图1所示,可以理解为由三个载荷矩阵ABC构成。

图 1. PARAFAC算法模型分解示意图

Fig. 1. Model decomposition schematic of PARAFAC algorithm

下载图片 查看所有图片

PARAFAC算法的数学表达式为

xijk=n=1Nainbjnckn+eijk,(1)

式中xijk为三维荧光阵内的元素,在石油类物质混合液三维荧光光谱分析中的实际意义为相对荧光强度(第k个样本第j个发射波长第i个激发波长处);ain为相对荧光强度(第n个组分第i个激发波长处);bjn为相对荧光强度(第n个组分第j个发射波长处);ckn为相对浓度(第k个样本第n个组分处);eijk为残差矩阵中的元素;N为混合液的组分数。(1)式对应到石油类物质混合液荧光数据阵X(I×J×K)中,样本数量为K,激发波长数量为I,发射波长数量为J

算法过程:确定模型因子数(即样品组分数)N,由X和初始化的BC估计载荷矩阵A;用得到的A的估计值和初始化的C,采用相同的方法估计B;由A的估计值和B的估计值估计C,具体过程见参考文献[ 10]。

2.2 ART算法

相比于传统的PARAFAC算法,新的二阶校正方法——ART算法——不需要预先估计组分数N,它在自身收敛过程中即可自行确定。ART算法是基于三线性模型发展而来的,三线性成分模型还可以用矩阵形式来描述:

X..k=k=1KcknanbTn+E..k,k=1,2,,K,(2)

式中X..kE..k分别为立方形矩阵XEK方向的第k层,an为与激发波长相关的荧光强度阵,bn为与发射波长相关的荧光强度阵。

将PARAFAC算法中的残差平方和sSSR按照组分方向写成:

sSSR=i=1Ij=1Jk=1Keijk2=k=1KX..k-n=1NcknanbTnF2(3)

D..k=X..k-n=1,nn0Nckn0anbTn,k=1,2,,K,(4)L=k=1KD..k-ckn0an0bTn02=k=1KtrDT..kD..k+aTn0an0bTn0bn0cTn0cn0-2k=1Kckn0aTn0D..kbn0,(5)

式中L为当组分数是n0时对应的损失函数的残差平方和,ckn0为第k个样本中组分数为n0时浓度阵里的元素,an0为样品组分数为n0的与激发波长相关的荧光强度阵,bn0为样品组分数为n0的与发射波长相关的荧光强度阵,cn0为组分数为n0的浓度阵。对an0求偏导,可得:

an0=k=1Kckn0D..kbn0an02cn02,(6)bn0=k=1Kckn0DT..kan0an02cn02(7)

λ=‖an02bn02cn04,则有

λan0=k=1Kckn0D..kk=1Kckn0DT..kan0(8)

这表明an0bn0分别是 k=1Kckn0D..k的左右奇异向量,即

an0,v,bn0=svdk=1Kckn0D..k,0,(9)

式中v为正交阵。

同理对ckn0求偏导,可得

ckn0=aTn0DT..kbn0an02bn02,(10)cn0=c1n0,c2n0,,ckn0,,cKn0T(11)

由此可以得到an0bn0cn0

ART算法的迭代过程如下:

1) 初始化ABC,由ABC和(4)式得到D..k

2) 初始化cn0,由(9)式和D..k得到an0bn0

① 应用(10)式和(11)式,由得到的an0bn0D..k计算cn0

② 应用(9)式,由得到的cn0D..k计算an0bn0

③ 重复步骤①和步骤②,直到cn0收敛为止。

3) 重复步骤1)和步骤2)得到an0+1bn0+1cn0+1

4) 重复步骤2)和步骤3),直到X收敛。

3 实验分析

3.1 样品配制

准备3个100 mL的容量瓶,分别加入0号柴油、95号汽油和普通煤油各1 g,经溶剂CCl4定容后,可得到10000 mg·L-1柴油、汽油和煤油的标准溶液。然后取适量标准溶液置于20个容量瓶中,之后加CCl4定容配制成不同浓度的样品,各样本的质量浓度如表1所示。其中,1~10号样品为校正样品,11~20号样品为预测样品。

3.2 实验设备

实验仪器采用可快速完成三维扫描的日立F-7000型荧光分光光度计。其参数设置如下:光电倍增管(PMT)电压为400 V,扫描速率为12000 nm·min-1,扫描步长为5 nm,入射狭缝为10 nm,出射狭缝为10 nm,激发波长、发射波长的扫描范围分别为250~430 nm和310~520 nm,设置发射起始波长始终超前20 nm于激发起始波长,以消除瑞利散射。响应时间为自动模式。仪器已经经过激发与发射校正处理。

表 1. 样品的质量浓度

Table 1. Mass concentration of samples

NumberMass concentration /(mg·L-1)NumberMass concentration /(mg·L-1)
No. 0 dieselNo. 95 gasolineKeroseneNo. 0 dieselNo. 95 gasolineKerosene
10.20.10.5110.30.30.3
20.10.50.2120.30.10.5
30.50.20.1130.50.30.1
40.50.50.5140.10.50.2
50.80.80.8150.40.30.2
61.01.01.0160.70.70.7
70.50.81.0170.60.90.8
80.81.00.5180.80.60.7
91.00.50.8190.90.80.6
100.81.00.8200.80.90.6

查看所有表

3.3 软 件

采用在MATLAB上编写的PARAFAC和ART算法程序进行处理,算法在MATLAB8.0及以上版本上实现。

4 结果与讨论

4.1 标准样品的三维与二维荧光光谱

在以上仪器设定条件下分别测量柴油、汽油、煤油标准溶液,得到其三维光谱如图2所示。从图2中不难发现,柴油、汽油、煤油标准溶液在其设定范围内都有强度不同的荧光产生,虽然荧光强度的变化范围较大,但每种油的荧光峰均有明显重叠。因此,当溶液中存在三种及以上石油类物质时,利用光谱方法难以定性或定量测定某种油的含量。实验所得三维荧光数据采用PARAFAC和ART算法分析处理。

图3是标准溶液的二维激发与发射光谱,用于石油类物质混合物的定性识别比对。

图 2. 标准溶液的三维荧光光谱。(a)柴油;(b)汽油;(c)煤油

Fig. 2. Three-dimensional fluorescence spectra of standard solutions. (a) Diesel; (b) gasoline; (c) kerosene

下载图片 查看所有图片

图 3. 柴油、汽油、煤油标准溶液的二维荧光光谱。(a)激发光谱;(b)发射光谱

Fig. 3. Two-dimensional fluorescence spectra of diesel, gasoline and kerosene standard solutions. (a) Excitation spectra; (b) emission spectra

下载图片 查看所有图片

4.2 混合样品的PARAFAC分析

利用F-7000型荧光分光光度计扫描样品,得到三维光谱。三维荧光光谱的高灵敏度使其自身难以抵御溶剂中杂质的干扰,故采用空白扣除法予以消除。

采用PARAFAC模型分析测量得到20×43×37的三维数据阵X。为确定模型的组分数,利用核一致诊断方法。组分数为2的核一致系数为77.9%;组分数为3的核一致系数为98.7%。组分数为2时,由于实验选取样品中汽油与煤油的荧光光谱形状相似,荧光峰位置接近,以致算法将汽油光谱与煤油光谱混淆,作为一个整体进行分解。因此,选定模型的实际组分数为3。

图4为因子数为3的三组分样品PARAFAC算法模型的分析结果。表2显示了三组分样品中柴油、汽油和煤油的预测浓度和回收率。

4.3 混合样品ART分析

采用ART算法分析测量得到20×43×37的三维数据阵X

ART算法处理分析数据时,其残差平方和sSSR随组分数N的变化曲线如图5所示。由图5可知,当N=3时,sSSR最小,说明ART方法不需要预先估计组分数N,它能在自身收敛过程中自行确定N值的大小。

图 4. PARAFAC算法模型估计的荧光光谱。(a)激发光谱;(b)发射光谱

Fig. 4. Fluorescence spectra estimated by PARAFAC algorithm model. (a) Excitation spectra; (b) emission spectra

下载图片 查看所有图片

表 2. 采用PARAFAC算法得到的混合溶液的预测质量浓度与回收率

Table 2. Predicted mass concentration and recovery rate of mixed solutions obtained by PARAFAC algorithm

NumberPredicted massconcentration /(mg·L-1)Recoveryrate /%Average recoveryrate /%
No. 0dieselNo. 95gasolineKeroseneNo. 0dieselNo. 95gasolineKeroseneNo. 0dieselNo. 95gasolineKerosene
110.2870.2790.29095.6793.0096.6795.60±3.6094.67±3.6695.49±4.49
120.2830.0920.47294.3392.0094.40
130.4790.2820.09195.8094.0091.00
140.0920.4850.19392.0097.0096.50
150.3760.2860.19294.0095.3396.00
160.6870.6650.67498.1495.0096.29
170.5750.8700.77895.8396.6797.25
180.7740.5640.67396.7594.0096.14
190.8570.7300.56695.2291.2594.33
200.7860.8850.57898.2598.3396.33

查看所有表

图6为ART算法对三维数据阵的分析结果。比对图6图3中三种石油类产品标准溶液的真实激发、发射光谱可以看到,谱图在荧光强度上发生了变化,但算法分辨的光谱图吻合得很好。三组分样品中柴油、汽油和煤油的预测质量浓度和回收率如表3所示。

图 5. (a) PARAFAC和(b) ART算法残差平方和随组分数的变化

Fig. 5. Residual square sum of (a) PARAFAC and (b) ART algorithms changing with component number

下载图片 查看所有图片

图 6. ART算法估计的荧光光谱。(a)激发光谱;(b)发射光谱

Fig. 6. Fluorescence spectra estimated by ART algorithm. (a) Excitation spectra; (b) emission spectra

下载图片 查看所有图片

表 3. ART算法得到的混合溶液的预测质量浓度与回收率

Table 3. Predicted mass concentration and recovery rate of mixed solutions obtained by ART algorithm

NumberPredicted massconcentration /(mg·L-1)Recoveryrate /%Averagerecovery rate /%
No. 0dieselNo. 95gasolineKeroseneNo. 0dieselNo. 95gasolineKeroseneNo. 0dieselNo. 95gasolineKerosene
110.2850.2930.26195.0097.6787.0096.58±2.1795.17±9.1795.90±8.90
120.2890.0860.48496.3386.0096.80
130.4870.2910.09597.4097.0095.00
140.0950.4700.20395.0094.00101.50
150.3850.2730.19096.2591.0095.00
160.6870.6730.65298.1496.1493.14
170.5870.8820.77297.8398.0096.50
180.7600.5790.68295.0096.5097.42
190.8650.7620.58996.1195.2598.17
200.7900.9010.59198.75100.1198.50

查看所有表

4.4 PARAFAC与ART算法分析结果的比较

与经典的PARAFAC算法相比,ART算法不需要预先估计组分数N,从图5中可以看出ART算法在处理分析数据时sSSRN的变化趋势,在N=3处出现了最小的sSSR,这是因为其计算过程中各组分是交替迭代收敛的,当分解组分数高于体系中的实际组分数时,分解残差急剧变大,从而识别出真实的组分数。因此,ART算法可推广到4组分以及更高组分的石油类物质混合物的同时定量检测中。PARAFAC算法得到的sSSRN的增大而减小,因而无法自行确定正确的N

图4图6图3的标准二维光谱进行比较,由实验结果可知,在对石油类物质混合物的定性分析上,PARAFAC和ART两种算法都能达到识别的效果,但ART算法分解得到的激发光谱比PARAFAC算法分解得到的激发光谱拟合得更好,可得到更多的光谱信息。因此,定性分析时选择ART算法较好。

PARAFAC和ART算法得到的预测样品的质量浓度和回收率分别见表2表3,分析表2表3可得到以下结论:1) PARAFAC算法得到柴油、汽油和煤油的平均回收率分别为(95.60±3.60)%、(94.67±3.66)%和(95.49±4.49)%,这表明PARAFAC算法能够测定三组分石油类物质混合液中各组分的浓度,但需预设组分数,且收敛速率缓慢;2) ART算法得到柴油、汽油和煤油的平均回收率分别为(96.58±2.17)%、(95.17±9.17)%和(95.90±8.90)%,这表明ART算法能够同时直接快速测量三组分石油类物质混合液中各组分的浓度,其预测浓度和真实浓度很接近,结果较为理想。

5 结论

分析石油类物质混合液时采用PARAFAC和ART算法均能达到预期目的,但ART算法的残差平方和sSSR随着组分数N的增大而增大,可自行确定正确的组分数而无需预设组分数N,且收敛速率快;两种算法在同时直接定量测定石油类污染物浓度时,在可接受误差范围内均能得到较好的回收率。

参考文献

[1] 高乐华. 我国海洋生态经济系统协调发展测度与优化机制研究[D]. 青岛: 中国海洋大学, 2012.

    Gao LH. Research on coordinative of development assessment and optimization of marine ecological-economic system in China[D]. Qingdao: Ocean University of China, 2012.

[2] 常江. 水中矿物油的FTIR-ATR光谱分析[D]. 秦皇岛: 燕山大学, 2014.

[3] 潘钊, 王玉田, 邵小青, 等. 三维荧光光谱及平行因子分析在石油类污染物检测分析中的应用[J]. 光谱学与光谱分析, 2012, 32(3): 714-718.

    Pan Z, Wang Y T, Shao X Q, et al. Application of PARAFAC method and 3-D fluorescence spectra in petroleum pollutant measurement and analysis[J]. Spectroscopy and Spectral Analysis, 2012, 32(3): 714-718.

[4] 周冰冰. 荧光光谱法水中矿物油的检测及研究[D]. 秦皇岛: 燕山大学, 2015.

    Zhou BB. Fluorescence detection and study of mineral oil in water[D]. Qinhuangdao: Yanshan University, 2015.

[5] 胡扬俊, 陈国庆, 朱纯, 等. 橙汁饮品的荧光光谱特性及模式识别[J]. 发光学报, 2013, 34(8): 1066-1072.

    Hu Y J, Chen G Q, Zhu C, et al. Fluorescence spectra characteristics and pattern recognition of orange juice beverages[J]. Chinese Journal of Luminescence, 2013, 34(8): 1066-1072.

[6] 金丹, 张玉钧, 李国刚, 等. 菲的三维荧光光谱特性研究[J]. 光谱学与光谱分析, 2009, 29(5): 1319-1322.

    Jin D, Zhang Y J, Li G G, et al. Study on three-dimensional fluorescence spectra of phenanthrene[J]. Spectroscopy and Spectral Analysis, 2009, 29(5): 1319-1322.

[7] 崔志成, 刘文清, 赵南京, 等. 水中油浓度快速测量方法研究[J]. 光谱学与光谱分析, 2008, 28(6): 1332-1335.

    Cui Z C, Liu W Q, Zhao N J, et al. Study on the method for rapid determination of oil concentration in water[J]. Spectroscopy and Spectral Analysis, 2008, 28(6): 1332-1335.

[8] 奚留华, 陈国庆, 朱焯炜, 等. 基于三维荧光光谱测定芝麻香型白酒的年份[J]. 激光与光电子学进展, 2016, 53(12): 123002.

    Xi L H, Chen G Q, Zhu Z W, et al. Determination of vintages of sesame flavor liquors by three-dimensional fluorescence spectroscopy[J]. Laser & Optoelectronics Progress, 2016, 53(12): 123002.

[9] 潘钊, 王玉田, 吴希军, 等. 胶束增敏荧光方法测定水中石油类成分及其含量[J]. 仪器仪表学报, 2013, 34(12): 2839-2845.

    Pan Z, Wang Y T, Wu X J, et al. Identification and concentration measurement of petroleum in water based on micelle sensitization fluorescence method[J]. Chinese Journal of Scientific Instrument, 2013, 34(12): 2839-2845.

[10] 杨丽丽, 王玉田, 鲁信琼. 三维荧光光谱结合二阶校正法用于石油类污染物的识别和检测[J]. 中国激光, 2013, 40(6): 0615002.

    Yang L L, Wang Y T, Lu X Q. Identification and measurement of petroleum pollutant by three-dimensional matrix fluorescence with second-order calibration methods[J]. Chinese Journal of Lasers, 2013, 40(6): 0615002.

[11] 刘佳, 吴海龙, 方冬梅, 等. 三维荧光二阶校正法用于血浆和尿液中柔红霉素的快速测定[J]. 分析化学, 2008, 36(3): 316-320.

    Liu J, Wu H L, Fang D M, et al. Rapid determination of daunomycin hydrochloride in plasma and urine using three-dimensional fluorescence coupled with second-order calibration[J]. Chinese Journal of Analytical Chemistry, 2008, 36(3): 316-320.

[12] 吕桂才, 赵卫红, 王江涛. 平行因子分析在赤潮藻滤液三维荧光光谱特征提取中的应用[J]. 分析化学, 2010, 38(8): 1144-1150.

    Lü G C, Zhao W H, Wang J T. Applications of parallel factor analysis in feature extraction of excitation-emission matrix spectrum of dissolved organic matter in red tide algae growth process[J]. Chinese Journal of Analytical Chemistry, 2010, 38(8): 1144-1150.

[13] 刘海龙, 吴希军, 田广军. 三维荧光光谱技术及平行因子分析法在绿茶分析及种类鉴别中的应用[J]. 中国激光, 2008, 35(5): 685-689.

    Liu H L, Wu X J, Tian G J. Three-dimensional fluorescence spectroscopy combined with parallel factor analysis as a complementary technique for green tea characterization[J]. Chinese Journal of Lasers, 2008, 35(5): 685-689.

[14] 郭卫东, 杨丽阳, 王福利, 等. 水库型河流溶解有机物三维荧光光谱的平行因子分析[J]. 光谱学与光谱分析, 2011, 31(2): 427-430.

    Guo W D, Yang L Y, Wang F L, et al. Parallel factor analysis for excitation emission matrix fluorescence spectroscopy of dissolved organic matter from a reservoir-type river[J]. Spectroscopy and Spectral Analysis, 2011, 31(2): 427-430.

[15] Christensen J H, Tomasi G. Practical aspects of chemometrics for oil spill fingerprinting[J]. Journal of Chromatography A, 2007, 1169(1/2): 1-22.

[16] 杜文. 多元校正及模型转移中的缺损数据重构和交替残差多线性方法研究[D]. 长沙: 湖南大学, 2016.

    DuW. Missing data recovery and alternating residual trilinearization applied to multivariate calibration and calibration transfer[D]. Changsha: Hunan University, 2016.

[17] 张燕. 化学计量学二阶校正法在荧光和色谱分析中的应用[D]. 长沙: 湖南大学, 2009.

    ZhangY. Applications of second-order calibration methods in fluorometric and chromatographic analysis[D]. Changsha: Hunan University, 2009.

[18] 刘佳. 化学计量学二阶校正法用于复杂体系抗癌药物的定量分析[D]. 长沙: 湖南大学, 2008.

    LiuJ. Second-order calibration methods applied to determination of antitumor drugs in complex systems[D]. Changsha: Hunan University, 2008.

[19] 邓琥, 尚丽平. 基于四组分同时测定的PARAFAC和SWATLD算法及比较[J]. 光谱学与光谱分析, 2009, 29(4): 1088-1092.

    Deng H, Shang L P. Comparison between PARAFAC and SWATLD based on four components simultaneous determination[J]. Spectroscopy and Spectral Analysis, 2009, 29(4): 1088-1092.

[20] 康超. 化学多维校正基础理论及其在复杂体系中的定量应用研究[D]. 长沙: 湖南大学, 2015.

    KangC. Chemical multi-way calibration fundamental theories and their applications to quantitative analysis in complex systems[D]. Changsha: Hunan University, 2015.

[21] 杜树新, 杜阳锋, 袁之报. 三维荧光光谱的特征区域选择方法[J]. 发光学报, 2012, 33(3): 341-345.

    Du S X, Du Y F, Yuan Z B. Characteristic region selection methods for three-dimensional fluorescence spectrometry[J]. Chinese Journal of Luminescence, 2012, 33(3): 341-345.

陈至坤, 弭阳, 沈小伟, 程朋飞. 基于PARAFAC和ART算法的油类污染物荧光检测[J]. 激光与光电子学进展, 2018, 55(1): 013004. Chen Zhikun, Mi Yang, Shen Xiaowei, Cheng Pengfei. Fluorescence Detection of Oil Pollutants Based on PARAFAC and ART Algorithms[J]. Laser & Optoelectronics Progress, 2018, 55(1): 013004.

本文已被 2 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!