光学学报, 2018, 38 (11): 1130005, 网络出版: 2019-05-09   

基于交替惩罚三线性分解的混合油液油种成分的检测 下载: 812次

Detection of Oil Species in Mixed Oil Based on Alternating Penalty Trilinear Decomposition
作者单位
1 燕山大学电气工程学院, 河北 秦皇岛 066004
2 燕山大学信息科学与工程学院, 河北 秦皇岛 066004
摘要
将三维荧光探测技术与交替惩罚三线性分解算法相结合,提出一种用于混合油液中油种成分检测的方法。首先利用FLS920型荧光光谱仪采集拟进行成分检测的混合油液(航空煤油和润滑油)在不同体积比配制条件下的20个样本的三维荧光光谱数据,利用Delaunay插值法对所获得的三维荧光光谱数据进行校正;然后利用核一致函数确定交替惩罚三线性分解算法解析三维荧光光谱数据时所需的组分数;最后利用均方根误差和相关系数矩阵,对交替惩罚三线性分解算法解析三维荧光光谱数据的效果进行评价。结果表明:在经交替惩罚三线性分解算法解析得到的均方根误差和相关系数矩阵中,非对角线上的元素值均满足所设阈值0.05和0.95的要求;在解决三维荧光光谱严重重叠的问题上,交替惩罚三线性分解算法优于平行因子算法,达到了对混合油液中油种成分快速检测的目的。
Abstract
Based on the combination of the three-dimensional fluorescence detection technique and the alternating penalty trilinear decomposition algorithm, a method is proposed to detect oil species in mixed oils. First, three-dimensional fluorescence spectra of 20 samples of the mixed oil (jet fuel and lube) with different volume ratios are obtained by the FLS920 fluorescence spectrometer, and calibrated by the Delaunay interpolation method. Then, the number of components required for the alternating penalty trilinear decomposition algorithm to analyze three-dimensional fluorescence spectra is determined by the core consistency function. Finally, the effectiveness of the alternating penalty trilinear decomposition algorithm in analyzing three-dimensional fluorescence spectra is evaluated by the root mean square error and the correlation coefficient matrix. The results show that both of the root mean square error and the non-diagonal elements in the correlation coefficient matrix meet the threshold requirements of 0.05 and 0.95, which are obtained by the alternating penalty trilinear decomposition algorithm. As for the solution of the serious overlapping problem of three-dimensional fluorescence spectra, the alternating penalty trilinear decomposition algorithm is superior to the parallel factor algorithm. The purpose of rapid detection of oil species in mixed oils can be achieved by the alternating penalty trilinear decomposition algorithm.

1 引言

近年来,随着全球经济的迅速发展,海洋上空的交通运输更加繁忙,旅客、货物的运输量与飞机的需求量呈现持续增长的态势。而海上恶劣的天气条件及飞机机械故障等因素可能会导致飞机的偶发性失事,造成严重的人员伤亡和财产损失,因此对海上失事飞机的搜救已成为世界各国不断加大投入研究的课题。众所周知,在飞机于海面失事的过程中,机体内的航空煤油与润滑油等油液会溢出机体,漂浮在海面上。如果可以对漂浮于海面上的疑似油液进行有效的检测和辨识,就可对该实际溢油源的性质加以判断,从而为下一步工作提供必要的技术支撑。因此,研究一种快速、高效的混合油液中油种成分的检测手段,对于海上搜救具有重要的现实意义。

目前,对混合油液中油种成分进行检测的方法主要有荧光光谱法[1]、红外光谱法[2]、气相色谱法[3]、气相色谱-质谱法[4]等。红外光谱法提供的信息较少,且检测精度易受环境和样本温度的影响[5];气相色谱法可提供精确的样本组分信息,但其所测组分容易受风化的影响,在实际使用中的检测误差较大[6];气相色谱-质谱法受风化程度的影响较小,但检测过程较繁琐,耗时较长,无法进行在线检测[7];荧光光谱法具有灵敏度高、选择性好、操作简便,以及可用于多组分混合物成分检测等优点[8]。鉴于此,本文采用荧光光谱法对航空煤油和润滑油的混合油液样本进行荧光特性分析。石油类物质的荧光光谱由多种多环芳烃的荧光光谱叠加而成,采用传统的化学方法检测其混合油液样本的油种成分比较困难[9],但是可利用二阶校正法具有的“数学分离”代替“物理或化学分离”的优势来加以解析[10]。二阶校正法与三维荧光光谱分析技术相结合可以实现多组分混合物的检测,该项技术具有广阔的应用前景。

本文利用FLS920型稳态荧光光谱仪测定20个航空煤油和润滑油混合油液样本的荧光光谱,并在保证不丢失荧光光谱信息的前提下进行激发和发射校正,以有效剔除实验仪器光谱特性带来的干扰;利用交替惩罚三线性分解(APTLD)算法和平行因子算法对航空煤油和润滑油混合油液样本检测的结果进行对比分析,以验证APTLD算法是否具有更优的检测能力。

2 混合油液中油种成分检测过程的实现

2.1 检测策略

快速、准确地检测混合油液中油种成分的关键在于获取被测混合油液样本的三维荧光光谱数据,并对其进行有效的分析和处理。处理过程的具体流程如图1所示,即分为5个步骤:样本荧光光谱数据的获取、原始荧光光谱数据的校正、组分数估计、三维荧光光谱数据算法处理、荧光光谱曲线的拟合和分析。

图 1. 混合油液中油种成分检测流程图

Fig. 1. Flow chart of oil species detection in mixed oil

下载图片 查看所有图片

2.2 检测基本原理

2.2.1 三线性分解模型

在进行荧光分析时,在选定的I个激发波长、J个发射波长范围内,对K个样本进行扫描,可以得到一个I×J×K型的三维响应数阵X。三维响应数阵X的分解实现过程可由三线性分解模型[11]表示,即

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

式中:i=1,2,…,I;j=1,2,…,J;k=1,2,…,K;N为体系中有荧光贡献的总组分数;xijk为三维响应数阵X中的元素(i, j, k),表示第k个样本在第i个激发波长、第j个发射波长时的荧光强度;ain为相对激发光谱阵AI×N中的元素(i, n);bjn为相对发射光谱阵BJ×N中的元素(j, n);ckn为相对浓度阵CK×N中的元素(k, n);eijk为三维残差阵EI×J×K中的元素(i, j,k)。对所得三维响应数阵X进行分解,可以得到具有清晰物理意义的唯一解ABC。三线性分解模型的图形表示如图2所示,其中G为三维核对角阵。

图 2. 三线性分解模型

Fig. 2. Trilinear decomposition model

下载图片 查看所有图片

2.2.2 核一致诊断法

核一致诊断法是基于三线性分解模型的化学秩估计方法。利用核一致函数计算经Tucker 3模型处理后的核心阵与理想情况下取正确组分数时的核心阵(对角元素为1,其他元素为0)之间的拟合程度,即核一致值,对体系的组分数进行估计。核一致函数[12]

1001-i=1Nj=1Nk=1N(gijk-tijk)2/i=1Nj=1Nk=1Ntijk2,(2)

式中:gijk为Tucker 3模型中核心阵的元素;tijk为理想情况下取正确组分数时核心阵的元素。在实际应用中,一般选取核一致值不小于最大值的60%作为体系的组分数估计结果。

2.2.3 APTLD算法

APTLD算法利用交替惩罚限制和交替最小二乘原理对3个不同的交替惩罚残差(AP残差)同时进行最小化,从而实现三线性分解[13]。APTLD算法的迭代过程如下。

1) 初始化载荷矩阵AB

2) 计算载荷矩阵

C=j=1JA+p(AT)+diag[1./diagm(BTB)]diag(bj)+pi=1IXTi(B+)Tdiag[1./diagm(ATA)]diag(ai)·j=1Jdiag(bj)ATA+pdiag[1./diagm(BTB)]diag(bj)+pj=1Jdiag[1./diagm(ATA)][diag(ai)]2+,(3)

式中:p为惩罚因子;bj为矩阵BJ×N的第j行元素矢量;Xi为三维响应数阵X的第i个水平面矩阵;ai为矩阵AI×N的第i行元素矢量;上标+为广义逆;diag为将矢量对角化;diagm为取矩阵对角元素化为一列向量。

3) 计算载荷矩阵A,并对其进行归一化。A的表达式为

A=k=1KXk{B+q(BT)+diag[1./diagm(CTC)]}diag(ck)+qk=1KXTj(C+)Tdiag(1./diagm(BTB))diag(bj)·k=1Kdiag(ck){BTB+qdiag[1./diagm(CTC)]}diag(ck)+qj=1Jdiag[1./diagm(BTB)][diag(bj)]2+,(4)

式中:Xk为三维响应数阵X的第k个水平面矩阵;q为惩罚因子;Xj为三维响应数阵X的第j个水平面矩阵;ck为矩阵CK×N的第k行元素矢量。

4) 计算载荷矩阵B,并对其进行归一化。B的表达式为

B=i=1IXi{C+r(CT)+diag[1./diagm(ATA)]}diag(ai)+rk=1KXTj(A+)Tdiag(1./diagm(CTC))diag(ck)·i=1Idiag(ai){CTC+rdiag[1./diagm(ATA)]}diag(ai)+rk=1Kdiag[1./diagm(CTC)][diag(ck)]2+,(5)

式中:r为惩罚因子。

5) 返回步骤3),依次迭代,直至满足收敛条件:

[SSSRm-SSSR(m-1)]/SSSR(m-1)10-6,(6)

式中:SSR为残差元素平方和;m为迭代次数。

6) 输出载荷矩阵ABC

2.2.4 评价标准

评价标准能反映所用算法分析结果的优劣,均方根误差(RMSE)和相关系数矩阵中非对角线元素值R的计算公式分别为

eRMSE=(ymixed-ypure)2/m',(7)R=C(1,2)/C(1,1)C(2,2),(8)

式中:ymixedypure分别为混合油液、纯油液经算法分解后的载荷矩阵BC中的光谱数据;m'为载荷矩阵BC中光谱数据的个数;R为相关系数矩阵中非对角线上的元素值;C(i, j)为矩阵(f,g)协方差矩阵中的元素,其中f为混合油液经算法分解后的载荷矩阵BC构成的矩阵,g为纯油液经算法分解后的载荷矩阵BC构成的矩阵。

3 实验

3.1 仪器设备

采用英国Edinburgh Instruments公司的FLS920型稳态荧光光谱仪测定实验样本的荧光光谱数据,光谱波长响应范围为200~900 nm,分辨率可达0.1 nm。该光谱仪配有制冷器,用来降低探测器的热噪声,液氮制冷范围为77~320 K,激发光源采用450 W氙灯,信噪比为6000∶1。所有实验均在室温下进行;激发与发射端狭缝的宽度为0.44 mm;设定激发波长范围为300~580 nm,发射波长范围为320~600 nm,步长均为10 nm;发射起始波长滞后激发起始波长20 nm,目的是避免瑞利散射的干扰。

3.2 样本与配制

实验样本为航空煤油、润滑油。在实验设计中,配制20个不同体积比的航空煤油和润滑油混合油液样本。按照表1所示的体积比,使用精密移液枪分别吸取相应体积的航空煤油和润滑油,并将它们置于烧杯中,经充分混合后,将混合油液移入比色皿中,待测定。

4 结果与讨论

4.1 光谱分析

在实验过程中,测量仪器的特性使得测量得到的荧光光谱皆为表观光谱,即未校正过的荧光光谱。因此,为了尽可能提高最终结果的精度,需要对光谱仪直接测量得到的光谱进行激发和发射校正。航空煤油和润滑油在校正前后的三维荧光光谱图和等高线图分别如图3图4所示,其中ExEm分别为激发波长和发射波长。由图3图4中航空煤油和润滑油的荧光光谱信息可知:航空煤油的最佳激发波长和发射波长分别为320 nm和340 nm,润滑油的最佳激发波长和发射波长分别为470 nm和530 nm;经激发和发射校正后,一些强度较弱的假峰得以去除,并使部分被掩盖的光谱信息凸显出来,明显消除了测量仪器本身对荧光光谱产生的干扰。

表 1. 混合油液中航空煤油和润滑油的配比

Table 1. Ratio of jet fuel to lube in mixed oil

SamplenumberVolume of jetfuel /mLVolume oflube /mLSamplenumberVolume ofjet fuel /mLVolume oflube /mL
13.000112.700.30
22.970.03122.400.60
32.940.06132.100.90
42.910.09141.801.20
52.880.12151.501.50
62.850.15161.201.80
72.820.18170.902.10
82.790.21180.602.40
92.760.24190.302.70
102.730.272003.00

查看所有表

图 3. 校正前后航空煤油的三维荧光光谱和等高线图。(a)校正前的三维荧光光谱;(b)校正前的等高线图;(c)校正后的三维荧光光谱;(d)校正后的等高线图

Fig. 3. Three-dimensional fluorescence spectra and contour maps of jet fuel before and after calibration. (a) Three-dimensional fluorescence spectrum before calibration; (b) contour map before calibration; (c) three-dimensional fluorescence spectrum after calibration; (d) contour map after calibration

下载图片 查看所有图片

4.2 组分数估计

在三线性分解模型建立过程中,组分数的准确估计是定性分析的关键。在本实验中,对所测得的前20个混合油液样本的荧光光谱数据进行激发和发射校正后,可构成29×29×20型三维响应数阵X,然后利用核一致函数对该三维响应数阵X进行组分数估计。航空煤油和润滑油混合油液样本的核一致值、残差平方和随组分数的变化如图5所示。当选取的组分数大于正确的组分数时,核一致值接近于0或为负数;当所选取的组分数小于或等于正确的组分数时,核一致值等于1或接近于1。由图5(a)可知:当组分数为1或2时,核一致值为1;当组分数大于2时,核一致值开始减小,偏离了三线性分解的模型。由图5(b)可知,组分数为2时的残差平方和显著小于组分数为1时的残差平方和,且与组分数为3时的残差平方和的差别较小。结合实验样本实际所含油种成分的情况,选择算法解析三维荧光光谱数据所需的组分数为2。

图 4. 校正前后润滑油的三维荧光光谱和等高线图。(a)校正前的三维荧光光谱;(b)校正前的等高线图;(c)校正后的三维荧光光谱;(d)校正后的等高线图

Fig. 4. Three-dimensional fluorescence spectra and contour maps of lube before and after calibration. (a) Three-dimensional fluorescence spectrum before calibration; (b) contour map before calibration; (c) three-dimensional fluorescence spectrum after calibration; (d) contour map after calibration

下载图片 查看所有图片

图 5. 三维响应数阵X。(a)核一致值;(b)残差平方和

Fig. 5. Three-dimensional response matrix X. (a) Core consistency value; (b) residual sum of squares

下载图片 查看所有图片

4.3 对比实验分析

4.3.1 基于APTLD算法的混合油液中油种成分检测

航空煤油和润滑油混合油液样本的荧光光谱出现了严重重叠,如图6所示的样本16的三维荧光光谱图,很难分辨出这两个油种的荧光特征峰。为了解决荧光光谱难分辨的问题,需要借助二阶校正法具有的“二阶优势”。利用4.2节中分析得到的组分数,对三维响应数阵X进行两因子分析,经APTLD算法解析得到的激发光谱和发射光谱如图7所示。其中,因子1为航空煤油组分,因子2为润滑油组分。由图7可知,混合油液样本解析后的荧光光谱与航空煤油、润滑油本身真实的荧光光谱的相似度都较高,说明APTLD算法对混合油液中各油种均具有良好的检测能力。

4.3.2 基于平行因子算法的混合油液中油种成分的检测

对三维响应数阵X进行两因子分析,经平行因子算法解析得到的激发和发射光谱谱图如图8所示,可见:在混合油液样本解析后的荧光光谱中,因子2与润滑油本身真实的荧光光谱的相似度较高,而因子1与航空煤油本身真实的荧光光谱的相似度较低。说明平行因子算法对混合油液中润滑油的检测能力与APTLD算法相差不大,但对航空煤油的检测能力低于APTLD算法。

图 6. 第16个样本的结果。(a)三维荧光光谱;(b)等高线图

Fig. 6. Result for the 16th sample. (a) Three-dimensional fluorescence spectrum; (b) contour map

下载图片 查看所有图片

图 7. APTLD算法解析得到的航空煤油和润滑油的结果。(a)激发光谱;(b)发射光谱

Fig. 7. Results for jet fuel and lube obtained by APTLD algorithm. (a) Excitation spectra; (b) emission spectra

下载图片 查看所有图片

图 8. 平行因子算法解析得到的航空煤油和润滑油的结果。(a)激发光谱;(b)发射光谱

Fig. 8. Results for jet fuel and lube obtained by parallel factor algorithm. (a) Excitation spectra; (b) emission spectra

下载图片 查看所有图片

4.3.3 两算法解析结果的对比分析

表2所示为经APTLD算法、平行因子算法解析后的航空煤油、润滑油纯油样本与其混合油液样本的激发、发射光谱曲线之间的RMSE及相关系数矩阵,其中hExhEm分别为航空煤油的激发、发射光谱的RMSE值,rExrEm分别为润滑油的激发、发射光谱的RMSE值。根据数据冗余度和精度的要求,设定RMSE和相关系数矩阵中非对角线元素的阈值分别为0.05和0.95。对于各光谱曲线之间的RMSE值,采用APTLD算法解析三维荧光光谱数据时,激发、发射光谱中航空煤油、润滑油对应光谱曲线的RMSE值均在阈值范围之内;采用平行因子算法解析三维荧光光谱数据时,激发、发射光谱中润滑油对应光谱曲线的RMSE值均在阈值范围之内,且比同种情况下经APTLD算法解析后的RMSE值大,而航空煤油对应光谱曲线的RMSE值均超出所设阈值。对于各光谱曲线之间的相关系数矩阵中的非对角线元素值,采用APTLD算法解析三维荧光光谱数据时,航空煤油、润滑油对应光谱曲线的数值均在阈值范围之内;而采用平行因子算法解析三维荧光光谱数据时,航空煤油、润滑油对应光谱曲线的数值均小于所设阈值。

由上述分析结果可知,在检测混合油液中各油种成分的效果方面,相对于平行因子算法,APTLD算法存在较大的优势,能够实现对混合油液中各油种成分快速检测的目的。

表 2. 纯油液与混合油液光谱曲线之间的RMSE与相关系数矩阵

Table 2. RMSE and correlation coefficient matrix between spectral curves of pure and mixed oil

SampleAPTLD algorithmParallel factor algorithm
RMSECorrelation coefficient matrixRMSECorrelation coefficient matrix
Jet fuelhEx=0.0223hEm=0.02861.0000 0.98260.9826 1.0000hEx=0.1232hEm=0.12491.0000 0.78620.7862 1.0000
LuberEx=0.0265rEm=0.03571.0000 0.97770.9777 1.0000rEx=0.0370rEm=0.04971.0000 0.81480.8148 1.0000

查看所有表

5 结论

本课题组利用三维荧光光谱技术和APTLD算法对混合油液中油种的成分进行了检测。研究结果表明,采用APTLD算法解析三维荧光光谱数据后,得到的RMSE和相关系数矩阵中非对角线上的元素值均满足所设阈值(0.05和0.95)的要求。在解决三维荧光光谱严重重叠的问题上,APTLD算法优于平行因子算法,能够实现对混合油液中油种成分的快速、准确检测。本研究提供了一种二阶校正法与荧光光谱分析技术相结合的方法,为判定实际溢油源的性质提供了必要的技术支持。

参考文献

[1] Yang R J, Dong G M, Sun X S, et al. Feasibility of the simultaneous determination of polycyclic aromatic hydrocarbons based on two-dimensional fluorescence correlation spectroscopy[J]. Spectrochimica Acta Part A, 2018, 190: 342-346.

[2] 孙通, 吴宜青, 李晓珍, 等. 基于近红外光谱和子窗口重排分析的山茶油掺假检测[J]. 光学学报, 2015, 35(6): 0630005.

    Sun T, Wu Y Q, Li X Z, et al. Discrimination of camellia oil adulteration by NIR spectra and subwindow permutation analysis[J]. Acta Optica Sinica, 2015, 35(6): 0630005.

[3] Bagur-González M G, Pérez-Castaño E, Sánchez-Viñas M, et al. . Using the liquid-chromatographic-fingerprint of sterols fraction to discriminate virgin olive from other edible oils[J]. Journal of Chromatography A, 2015, 1380: 64-70.

[4] Tseng W C, Chen P S, Huang S D. Optimization of two different dispersive liquid-liquid microextraction methods followed by gas chromatography-mass spectrometry determination for polycyclic aromatic hydrocarbons (PAHs) analysis in water[J]. Talanta, 2014, 120: 425-432.

[5] 韩云利. 海上溢油的油指纹鉴定研究[D]. 大连: 大连海事大学, 2008: 12- 15.

    Han YL. Study on the identification of oil spill fingerprints[D]. Dalian: Dalian Maritime University, 2008: 12- 15.

[6] 包木太, 文强, 崔文林, 等. 六种成品油的正构烷烃色谱指纹提取与鉴别[J]. 西安石油大学学报(自然科学版), 2007, 22(1): 87-90, 125.

    Bao M T, Wen Q, Cui W L, et al. Chromatogram fingerprint extraction of n-alkane of six product oils and identification of them[J]. Journal of Xi'an Shiyou University (Natural Science Edition), 2007, 22(1): 87-90, 125.

[7] 刘敏燕, 徐恒振. 水上溢油源快速鉴别技术[M]. 北京: 中国环境出版社, 2014: 3- 6.

    Liu MY, Xu HZ. Rapid identification technology for oil spill sources on water[M]. Beijing: China Environmental Press, 2014: 3- 6.

[8] 易莹, 周艳伟, Bridgeman J, 等. 三维荧光光谱法测定市政污水中有机污染物[J]. 环境监测管理与技术, 2017, 29(5): 37-39, 71.

    Yi Y, Zhou Y W, Bridgeman J, et al. Determination of organic pollutants in municipal sewage by three-dimensional fluorescence spectroscopy[J]. The Administration and Technique of Environmental Monitoring, 2017, 29(5): 37-39, 71.

[9] 程朋飞, 王玉田, 陈至坤, 等. 基于自加权交替三线性分解的荧光检测油类污染物[J]. 光谱学与光谱分析, 2016, 36(7): 2162-2168.

    Cheng P F, Wang Y T, Chen Z K, et al. The fluorescence detection of oil pollutants based on self-weighted alternating trilinear decomposition[J]. Spectroscopy and Spectral Analysis, 2016, 36(7): 2162-2168.

[10] 吴海龙, 李勇, 康超, 等. 三维荧光化学多维校正方法研究新进展[J]. 分析化学, 2015, 43(11): 1629-1637.

    Wu H L, Li Y, Kang C, et al. Research progress of three-dimensional fluorescence coupled with chemical multiway calibration[J]. Chinese Journal of Analytical Chemistry, 2015, 43(11): 1629-1637.

[11] Li Y N, Wu H L, Qing X D, et al. The maintenance of the second-order advantage: second-order calibration of excitation-emission matrix fluorescence for quantitative analysis of herbicide napropamide in various environmental samples[J]. Talanta, 2011, 85(1): 325-332.

[12] Bai X M, Liu T, Liu D L, et al. Simultaneous determination of α-asarone and β-asarone in Acorus tatarinowii using excitation-emission matrix fluorescence coupled with chemometrics methods[J]. Spectrochimica Acta Part A, 2018, 191: 195-202.

[13] 蔡晴, 吴海龙, 李元娜, 等. 三维荧光二阶校正法同时测定水果样中α-萘乙酸和吲哚-3-乙酸[J]. 化学通报, 2011, 74(1): 47-53.

    Cai Q, Wu H L, Li Y N, et al. Rapid determination of α-naphthylacetic acid and indole-3-acetic acid in fruits extract by using three-dimensional fluorescence coupled with second-order calibration method[J]. Chemistry, 2011, 74(1): 47-53.

孔德明, 张春祥, 崔耀耀, 李雨蒙, 王书涛. 基于交替惩罚三线性分解的混合油液油种成分的检测[J]. 光学学报, 2018, 38(11): 1130005. Deming Kong, Chunxiang Zhang, Yaoyao Cui, Yumeng Li, Shutao Wang. Detection of Oil Species in Mixed Oil Based on Alternating Penalty Trilinear Decomposition[J]. Acta Optica Sinica, 2018, 38(11): 1130005.

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

相关论文

加载中...

关于本站 Cookie 的使用提示

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