孔德明,董 瑞,崔耀耀,王书涛
1.燕山大学电气工程学院,河北 秦皇岛 066004 2.燕山大学信息科学与工程学院,河北 秦皇岛 066004 3.Department of Telecommunications and Information Processing,Ghent University,B-9000 Ghent,Belgium
随着生产技术的提高,人类对石油的需求量也越来越大,然而石油从开采到使用过程中,不可避免地会有石油发生泄漏[1]。其中,海洋石油的开采、含油废水的排放以及海上交通事故是造成海面溢油的主要原因。海面溢油会对海洋环境造成污染,并且溢油中的多环芳烃类化合物是一种有毒物质,会危及海洋生物的生命安全[2]。因此,对海面溢油进行有效检测并准确分析其成分具有重要的实际意义。但由于海面的状况比较复杂,不同海域的自然状况不同。特别地,不同海域中海水的盐度会有明显的区别[3],所以在不同的海水盐度条件下,准确、快速地识别海面溢油,对保护海洋生态环境具有重要意义。
石油是一种具有荧光效应的有机混合物,其所含的成分比较复杂,并且不同类型的石油中所含的成分及其含量不同。在众多检测石油的方法中,三维荧光光谱法是一种有效检测石油的方法[4]。杨丽丽等利用平行因子分析和自加权交替三线性分解分别实现了对混合油的分类[5]。陈至坤等利用平行因子分析对0#柴油、95#汽油和煤油的三维荧光光谱矩阵进行解析,实现了3种石油浓度的预测[6]。但是平行因子分析对组分敏感,收敛速度较慢,容易产生局部最优解,从而导致识别效果不理想[7]。三阶校正在保证二阶校正“二阶优势”的基础上引入新一维度的数据,可以提高方法的分辨能力,能够得到更好的分析结果[8-9]。
本文利用三维荧光光谱技术结合交替加权残差约束四线性分解(alternating weighted residue constraint quadrilinear decomposition,AWRCQLD)算法,并通过引入海水盐度作为新一维度数据,实现了不同盐度条件下对航空煤油和润滑油的混合油定性及定量分析。
三阶校正是在二阶校正的基础上发展起来的一种处理四维数据的模型。其中,四维数据是在三维数据基础上增加新的一维信息,从而提高校正模型分辨数据的能力,使得分析结果更加准确。所以三阶校正可以在保留“二阶优势”的基础上更加精确地对数据进行解析,使得三阶校正模型具有“三阶优势”[10-13]。四线性分解的原理如图1所示。
图1 四线性分解模型图Fig.1 Quadrilinear decomposition model
在L个不同实验条件下采集K个实验样本的I个激发波长和J个发射波长下荧光强度,会得到一个大小为I×J×K×L的四维矩阵X,将X中的每一个元素记作Xijkl。四线性分解模型如式(1)
(1)
其中,i=1,2,…,I,j=1,2,…,J,k=1,2,…,K,l=1,2,…,L;N为实验样本中具有荧光效应的组分数;ain,bjn,ckn,dln,eijkl分别是荧光强度矩阵、相对激发光谱矩阵、相对发射光谱矩阵、第四维度矩阵及四维残差矩阵中的元素。
AWRCQLD算法采用四线性模型的伪完全拉伸矩阵形式来设计算法的目标函数,通过使用交替加权的残差函数作为四线性模型误差的约束项,最小化四个新的基于最小二乘的目标函数[13]。通过把加权残差函数作为约束项构造的四个新的目标函数如式(2)—式(5)
(2)
(3)
(4)
(5)
其中,λA,λB,λC和λD为约束项的系数。通过选择约束项的系数来平衡每个目标函数中两个部分的权重,可以拟合四线性模型的损失函数。AWRCQLD算法以交替的方式最小化四个目标函数,具体的方法如式(6)—式(9)
(6)
(7)
(8)
(9)
AWRCQLD算法的步骤如下:
(1)随机初始化矩阵A,B,C;
(2)使用式(9)计算矩阵D;
(3)使用式(6)计算矩阵A并将其归一化;
(4)使用式(7)计算矩阵B并将其归一化;
(5)使用式(8)计算矩阵C并将其归一化;
(6)使用式(9)计算D;
采用纯净水和人工海盐配制了3种盐度的模拟海水,并以此为溶剂来配制油类实验样本。具体配制步骤如下:(1)利用精密电子秤称取3份适量的人工海盐并将其分别移入3个烧杯中,分别加入适量的纯净水,用玻璃棒搅拌使其充分溶解,最后分别移入3个容量瓶中定容,配成3种盐度分别为0,10和20的模拟海水;(2)为了使石油能够充分溶于模拟海水溶剂,向上述配制的模拟海水溶剂中加入适量十二烷基硫酸钠(SDS),得到3种盐度条件下浓度为0.1 mol·L-1的SDS溶剂;(3)利用精密电子秤分别称取0.1 g的航空煤油和润滑油并移入两个烧杯中,然后加入浓度为0的SDS溶剂溶解。将溶解后的航空煤油和润滑油分别移入两个10 mL容量瓶中并用相应的SDS溶剂对其进行定容,得到浓度为10 mg·mL-1的航空煤油和润滑油溶液,此为一级储备液;(4)使用移液枪分别移取1 000 μL的航空煤油和润滑油的一级储备液于两个10 mL的容量瓶中,用盐度为0的SDS溶剂对其进行稀释并定容,得到浓度为1 mg·mL-1的航空煤油和润滑油溶液,此为二级储备液;(5)用移液枪分别移取不同体积的航空煤油和润滑油的二级储备液于10 mL的容量瓶中混合,并加入盐度为0的SDS溶剂对其进行定容,得到11组混合油溶液的实验样本,每组实验样本中航空煤油和润滑油的具体浓度如表1所示,其中,C1—C8为校正样本,T1—T3为预测样本;(6)分别利用盐度为10和20的SDS溶剂按照上述步骤(3)—(5)再次配制混合油实验样本,共计得到33组混合油实验样本。
利用FLS920荧光光谱仪采集样本的荧光光谱数据。设定其激发与发射狭缝宽度为0.10 mm,发射波长范围为280∶5∶500 nm,激发波长范围为260∶10∶480 nm,为了避免瑞利散射的干扰,设置发射波长滞后激发波长20 nm。
表1 实验样本配制浓度(mg·mL-1)Table 1 Experimental sample preparation concentration (mg·mL-1)
对采集的光谱数据进行预处理,利用扣除空白法去除溶液中的拉曼散射。以T3为例,图2(a)—(f)是T3分别在0,10和20盐度下的三维荧光光谱图和等高线图。其中,Em为发射波长,Ex为激发波长。由图2可知,3种盐度下主要荧光峰形相似,当海水的盐度变化时,航空煤油和润滑油的混合油的荧光强度也发生变化。因此,可利用不同盐度条件下的三维光谱数据矩阵构建四维光谱数据矩阵。
图2 样本T3在不同海水盐度条件下的荧光光谱图Fig.2 Fluorescence spectra of sample T3 at different seawater salinity
将盐度作为新的一维数据引入,用来构建四维光谱数据矩阵,所得到四维光谱数据矩阵X的大小为23×45×11×3,其中23表示激发波长数,45表示发射波长数,11表示实验样本的总数(8个校正样本,3个预测样本),3表示3种不同的盐度条件。
首先,使用核一致诊断法对溶液的组分数进行估计,如图3所示。当组分数小于等于2时,核一致值大于75%;当组分数大于2时,核一致值迅速降低至0。因此,选择2作为四线性分解的最佳组分数。
图3 核一致曲线Fig.3 Curves of core consistency value
然后,利用AWRCQLD算法对数据矩阵X进行组分数为2的四线性分解,其解析结果如图4所示。图4(a)为航空煤油、润滑油及解析后的混合油的相对发射光谱,图4(b)为航空煤油、润滑油及解析后的混合油的相对激发光谱,图4(c)为第四维度的解析结果。由图4(a)和(b)可知,混合油中的factor1与润滑油的光谱曲线重合度良好,因此可以判定混合油中的factor1为润滑油;混合油中的factor2与航空煤油的光谱曲线重合度良好,因此可以判定混合油中的factor2为航空煤油。润滑油的荧光主峰位于λem/λex=350/290~300 nm处,航空煤油的荧光主峰位于λem/λex=340~350/280~300 nm处。由图4(c)可知,当盐度增加时,润滑油的荧光强度先增大后减小,在盐度为10时,润滑油的荧光强度最大;当盐度增加时,航空煤油的荧光强度先减小后增大,在盐度为10时,航空煤油的荧光强度最小。原因可能是海水中的离子会与溶液中的有机物形成络合物,同时海水中的胶体粒子会吸附溶液中的离子;随着盐度的增加,胶体的稳定性可能会受到影响,会导致胶体聚集、再分散与溶解,而航空煤油和润滑油的成分有区别,可能会因此受到不同影响。综上可知,AWRCQLD算法能够对航空煤油和润滑油的混合油进行良好的定性分析。
图4 定性分析结果(a):相对发射光谱图;(b):相对激发光谱图;(c):第四维度分解图Fig.4 Qualitative analysis results(a):Relative emission spectra;(b):Relative excitation spectra; (c):Fourth dimensional decomposition graph
经过AWRCQLD算法对四维光谱数据矩阵X四线性分解后,得到相应预测样本的相对荧光光谱矩阵和浓度矩阵,根据所得到的相对荧光光谱矩阵和浓度矩阵对实验样本进行回归分析,如图5所示。图5(a)为航空煤油的回归曲线,图5(b)为润滑油的回归曲线。表2所示为三个预测样本中润滑油和航空煤油的预测浓度、回收率及均方根误差(root mean square error,RMSE)。由表2可知,航空煤油和润滑油的回收率结果较为理想,航空煤油的回收率范围为100.2%~109%,均方根误差为0.002 1 mg·mL-1;润滑油的回收率范围为91.8%~109.3%,均方根误差为0.004 8 mg·mL-1。因此,AWRCQLD定量分析结果良好。
表2 对航空煤油和润滑油定量分析的结果(mg·mL-1)Table 2 Quantitative analysis results of jet fuel and lube (mg·mL-1)
表3所示为实验样本的回归曲线方程及航空煤油和润滑油的浓度与荧光强度的线性相关系数,其线性相关系数分别为0.893和0.838 5,一般认为,当两个变量相关系数大于等于0.8时,表明两者高度相关,所以航空煤油和润滑油的浓度与其荧光强度在0.01~0.1 mg·mL-1的浓度范围内线性相关。
通过引入盐度维的方式构建了四维光谱数据阵,并利用三阶校正方法中的AWRCQLD对航空煤油和润滑油的混合油进行分析,得到合理的光谱分辩和浓度预测。比较0,10和20盐度下航空煤油和润滑油的荧光强度,发现海水盐度对航空煤油和润滑油的荧光强度有影响。在0~20盐度范围内,航空煤油和润滑油的荧光强度随着海水盐度的增加出现不同的变化趋势。航空煤油的荧光强度先减小后增大,润滑油的荧光强度先增大后减小。经过AWRCQLD算法解析后航空煤油的回收率为100.2%~109%,均方根误差为0.002 1 mg·mL-1;润滑油的回收率为91.8%~109.3%,均方根误差为0.004 8 mg·mL-1。本文将海水盐度作为新的一维信息引入,实现了对不同海水盐度条件下航空煤油和润滑油混合油的定性和定量分析,为不同海水盐度条件下海面溢油的检测提供了参考。
图5 航空煤油和润滑油的浓度回归曲线(a):航空煤油浓度回归曲线;(b):润滑油浓度回归曲线Fig.5 Concentration regression curve of jet fuel and lube(a):Concentration regression curve of jet fuel; (b):Concentration regression curve of lube
表3 回归曲线方程及相关系数Table 3 Regression curve equation and correlation coefficient
注:x表示实验样本的浓度,y为实验样本的荧光强度