作者简介: 车少敏, 女, 1996年生, 西安交通大学能源与动力工程学院硕士研究生 e-mail: cheshaomin@foxmail.com
空气被夹带进入发动机润滑油系统形成油气两相流后, 会严重影响润滑系统的正常运行。因此实现发动机润滑油系统中油气两相流含气率的准确和快速测量具有十分重要的意义。基于光谱匹配方法, 对发动机润滑油系统在弯管处的油气两相流含气率进行了测量研究。首先使用近红外、 可见光和紫外光谱仪获取油气两相流在2个流速、 5个温度工况下0.10%~1.00%含气率范围(间隔0.06%)的吸光度数据并进行分析, 确认油气两相流在三个波段下均有吸收, 且吸收强度与含气率有关; 其次提出将数据预处理结合光谱相似性测度方法并应用于弯管含气率光谱分析中, 显著降低了基于原始光谱的含气率预测的最大相对误差。在近红外波段通过使用数据增强方法如中心化(center)和标准化(autoscaling)并结合光谱角余弦方法, 得出新润滑油的含气率最大相对误差相对于原始光谱得到的数值从48%降低到36%。利用该方法分别在三个波段进行含气率预测, 试验工况包括2个流速、 5个温度, 分析了两相流的温度和流速对含气率预测的影响。温度30.0 ℃, 流速5.1 m·min-1下, 含气率信息在紫外波段(193.5~413.8 nm)与光谱向量的方向或形状差异的联系更为紧密, 含气率预测的最大相对误差仅为6%。在近红外和可见光波段, 流速或温度一定时, 最大相对误差随着温度或流速的增加均有减小趋势。紫外波段下温度对含气率预测效果没有特定的影响规律, 随着两相流的流速增加, 含气率预测的最大相对误差有增大趋势。结果表明对于透光性较好的新润滑油, 使用紫外光谱仪进行含气率数据的采集, 含气率的预测选用标准化预处理结合光谱角余弦方法, 得到的含气率预测最大相对误差最小。
When air is entrained into the engine lubricating oil system to form an oil and gas two-phase flow, it will seriously affect the normal operation of the lubrication system. Therefore, it is very important to realize the accurate and rapid measurement of the void fraction of the oil-gas two-phase flow in the engine lubricating oil system. This paper determines the void fraction in the engine lubrication system at the elbow based on the spectral matching method. Firstly, the absorption data of oil-gas two-phase flow were obtained at two flow rates and five temperature conditions, covering a gas content range of 0.10% to 1.00% (with an interval of 0.06%). This was accomplished by utilizing near-infrared, visible, and ultraviolet spectrophotometers, followed by rigorous analysis. It was ascertained that the oil-gas two-phase flow demonstrates absorption across all three wavelengths, with the intensity of absorption being correlated to the gas content. Secondly, the data preprocessing method combined with spectral similarity measure is proposed and applied to the gas content spectral analysis of bent pipe, significantly reducing the maximum relative error of gas content prediction based on the original spectrum. Using data enhancement methods such as center and autoscaling combined with Spectral Angle Cosine methods in the near-infrared spectrum, the maximum relative error of the gas content of the new lubricating oil was reduced from 48% to 36% relative to the original spectrum. The method predicts the void fraction of gas-oil two-phase flow in three bands respectively. The experimental conditions include two flow rates and five temperatures, and the influence of the temperature and flow rate of the two-phase flow on the void fraction prediction is analyzed. At the temperature of 30.0 ℃ and the flow rate of 5.1 m·min-1, the gas content information in the ultraviolet band (193.5~413.8 nm) is more closely related to the spectral feature of the direction or shape difference of the spectral vector and the maximum relative error of the gas content prediction is only 6%. In the near-infrared and visible bands, the maximum relative error decreases with the increase of temperature or velocity when the flow rate or temperature is constant. There is no specific effect of temperature on the prediction of gas content in the ultraviolet band. With the increase of the two-phase flow rate, the maximum relative error of gas content prediction tends to increase. The results show that for new lubricating oil with good light transmittance, the gas content data is collected by the ultraviolet spectrometer, and the maximum relative error of gas content prediction is minimum by using standardized pretreatment combined with spectral Angle cosine method.
弯管为发动机润滑油系统中常见的一种润滑油输送管道结构。研究航空发动机润滑油系统在弯管处的油气两相流含气率检测问题, 可以为航空发动机测试系统提供理论依据。含气率可以反映出两相流的流动结构和流动规律, 其测定对于工业过程的安全监控以及正常运行是十分必要的。因此对于含气率的测量方法, 要求具有准确、 快速和操作简便的特点。
光谱法检测含气率的基本原理是Beer-Lambert定律, 作为一种非侵入式方法, 有不易受流型干扰、 快速灵敏的优点。光谱相似性测度分析法也叫光谱匹配法, 是一种通过对比不同光谱曲线之间的光谱特征, 如形状、 方向的相似性和距离大小等的分析方法, 包括欧式距离、 马氏距离、 光谱角余弦、 相关系数、 信息熵等特征参数[1]。张长兴等[2]将光谱匹配法与盲元检测相结合, 提高热红外高光谱成像仪的盲元检测精度。张从征等[3]利用光谱匹配法, 实现多个不同波段LED光源对目标光谱的拟合。李宏达等[4]利用该方法, 实现对土壤有机碳的含量估算。胡鹏飞等[5]通过光谱匹配法, 利用人造光源实现太阳光的模拟。王宝明[6]研究了多种化合物在红外光下的光谱特征, 得到了各个化合物的光谱带宽度和位置, 并对红外辐射材料的光谱特征匹配问题进行了探究。
但是单一的光谱匹配方法只考虑了光谱某一方面的特征, 在样品光谱的单个特征不明显时, 存在光谱识别精度低的问题, 因此不适用于复杂的光谱识别。针对这种问题, 近年来对光谱匹配的综合测度方法有以下研究。安斌等[7]将光谱角匹配法(spectral angle mapping, SAM)和欧氏距离法相结合, 提出了一种兼顾光谱亮度与光谱方向的距离测度方法, 结果表明该方法的效果比传统分类方法好。闻兵工等[1]在地物识别中综合考虑光谱的形状以及曲线信息熵特征, 提出了一种新的光谱匹配方法, 结果证明该方法是可行的。孔祥兵等[8]将光谱向量的大小, 形状和信息量这三种光谱特征进行融合, 提出了一种新的匹配方法(spectral pan-similarity measure, SPM)并成功应用于矿物识别。方圣辉等和张浚哲等[9, 10]提出了一种新方法, 可以适用于不同的光谱影像以及不同的匹配目的, 并且分类识别精度更高。魏祥坡等[11]采用光谱角余弦和光谱相关系数两者的加权平均, 对地物进行分类。王展等[12]提出一种采用多元散射校正结合增广拉格朗日的方法减弱光谱数据中的散射噪声影响, 再使用光谱角余弦测度的新方法, 应用于矿物和壁画颜料的识别中。刘益岑等[13]将光谱匹配法与最小距离匹配相结合, 应用于绝缘子污秽类型检测。杨文宗等[14]利用光对视觉的感知特性, 提出融合光谱分析与带通能量积分的方法, 用于壁画虚拟修复。郭业才等[15]将光谱匹配与特征融合相结合, 实现对高光谱图像降维, 从而提取目标的特征的算法, 用于高光谱视频目标跟踪。齐海超等[16]利用多光谱图像技术实现对马赛克的去除。Riahi等[17]结合光谱角匹配方法与主成分分析和匹配滤波等方法, 利用遥感卫星图像, 实现对矿物的识别与检测。Paul等[18]提出一种检测光谱特征细微变化的算法, 实现对食品新鲜程度的检测。Aulich等[19]针对三结太阳能电池的所有三个子电池系统, 改进光谱表征方法, 并进一步验证该方法在匹配条件之外的相关性。Harati等[20]将光谱匹配技术用于多阻尼比下的频谱匹配, 并在地震记录等方面取得良好结果。Challagulla等[21]利用光谱匹配的原理, 对地震的动参数进行研究, 以减小地震造成的经济损失。目前还没有见到将光谱相似性测度技术应用于油气两相流含气率的光谱检测的报道。
光谱相似性测度技术是基于对比不同光谱曲线之间的形状、 方向以及距离等特征的方法, 因此不需要考虑光谱强度和含气率之间的线性关系, 而是充分利用丰富的光谱特征信息进行光谱的相似度分析。如果能够将该方法应用于油气两相流含气率的检测, 预计会得到更好的光谱识别结果以及含气率预测结果。基于上述现状, 本工作使用近红外、 可见光和紫外光纤光谱仪获取油气两相流含气率吸光度数据, 并提出将数据挖掘技术结合光谱相似性测度方法, 用于发动机润滑油系统弯管处的油气两相流含气率的测量中, 找到一个合适的光谱相似性测度方法对含气率的大小进行准确、 快速的预测。重点是不同波段下体现含气率大小信息最为明显的光谱特征, 为航空发动机测试系统提供理论依据。
光谱吸光度法检测弯管中油气两相流的含气率大小是基于Beer-Lambert定律。该定律描述了从光源发出的一束特定波段范围的光, 在完全覆盖并经过具有一定厚度的均匀待测物后光的强度会发生变化, 这是由于一部分光被待测物吸收, 其中待测物的浓度c、 待测物的厚度l和待测物的吸收强度A之间的关系可以用式(1)描述
式(1)中: A为待测物的吸收强度; I为出射光强; I0为入射光强; K为待测物质的吸收系数(与待测物有关); l为待测物的厚度, cm; c为待测物的浓度, g· L-1。
在本实验中, 待测物的厚度l为弯管处的垂直厚度, 且固定不变; K与待测物有关, 当待研究的润滑油确定后, 油气两相流的吸收系数K也为定值。因此在相同的条件下, 油气两相流的吸收强度A仅与油气两相流含气率的大小有关。以新润滑油为实验对象, 将连续改变的不同波段范围的光, 经过不同含气率大小的油气两相流, 油气两相流会吸收某些波段范围的光, 使其强度发生变化, 则从光谱仪获取的不同吸光度数据能反映出油气两相流中不同含气率大小的信息。
搭建的弯管处的油气两相流含气率测量实验系统如图1所示。
含气率的大小可以通过式(2)计算
式(2)中: ε 为体积含气率, 取流通池与流通回路中含气率的大小相同; Vt为整个流通管道的体积大小, (mL); V为流通管道中所含润滑油的体积, (mL)。
选取新润滑油作为实验对象, 分别在近红外、 可见光和紫外波段下, 采集油气两相流含气率在0.10%~1.00%范围内, 间隔0.06%一共16组吸光度数据。实验对象为新润滑油时, 采集工况为2个流速: 5.1和7.7 m· min-1, 5个温度: 室温(21.0 ℃)、 30.0 ℃、 35.0 ℃、 40.0 ℃和43.6 ℃共10个工况; 每个波段范围下各取16组数据, 组成数据集进行数据分析。图2、 图3、 图4分别为以新润滑油为实验对象时, 油气两相流在室温(21.0 ℃), 流速5.1 m· min-1工况下, 含气率为0.10%、 0.28%、 0.46%、 0.58%和0.82%时的近红外波段、 可见光波段和紫外光波段的原始光谱吸光度图。其中图(a)和图(b)分别为油气两相流的吸收谱和含气率与各个波长吸光度相关性的R2图。
![]() | 图2 新润滑油在近红外波段的原始吸收光谱和含气率与波长的相关图 (a): 吸收光谱; (b): R2Fig.2 Origin absorbance spectra of new oil in near-infrared band and R2 (a): Absorption spectra; (b): R2 |
![]() | 图3 新润滑油在可见光波段的原始吸收光谱和含气率与波长的相关图 (a): 吸收光谱; (b): R2Fig.3 Origin absorbance spectra of new oil in visible band and R2 (a): Absorption spectra; (b): R2 |
![]() | 图4 新润滑油在紫外光波段的原始吸收光谱和含气率与波长的相关图 (a): 吸收光谱; (b): R2Fig.4 Origin absorbance spectra of new oil in ultraviolet band and R2 (a): Absorption spectra; (b): R2 |
分析可知, 以新润滑油为实验对象时, 不同含气率下的油气两相流在近红外、 可见光和紫外波段均有特定波段范围的吸收, 并且吸收强度与含气率的大小有关, 但是吸收强度和含气率大小之间均没有直观的线性相关性, 这是由于在实验过程中, 原始光谱信号会引入较多的噪声信息和光散射信息, 而散射会扭曲油气两相流的吸收强度和含气率之间的线性相关性。因此, 采用基于光谱相似性测度的光谱分析技术, 充分挖掘并利用丰富的光谱特征信息, 进行油气两相流含气率大小的预测。
由于采集到的原始光谱信息量巨大, 并且存在较多的噪声信息和光散射等信息, 使得吸收强度和含气率大小之间没有直观的线性相关性, 不能仅仅使用传统的基于Beer-Lambert定律的方法分析油气两相流的光谱数据, 需要采用基于光谱相似性测度的光谱分析技术, 充分挖掘并利用丰富的光谱特征信息, 根据某一特定的光谱特征以及对应的光谱相似性度量函数判断未知含气率大小的光谱与光谱库中各个已知含气率大小的光谱之间的相似性, 选取出与未知含气率大小的光谱相似性最大的参考光谱, 将该参考光谱的含气率值作为待求含气率的大小, 通过设置函数阈值去除预测坏值, 最终实现油气两相流含气率大小的预测。同时, 为了充分挖掘并利用丰富的光谱特征信息, 需要根据实际的光谱数据选择合适的数据降维方法以及数据预处理方法, 以便得到更好的含气率预测效果。
由实验得到的吸收光谱可以看出, 0.10%~1.00%的边缘范围含气率大小可能由于实验操作不当而导致采集到的吸光度数据有误, 因此为了得到更好的含气率预测效果, 从0.10%~1.00%范围的16组含气率吸光度数据中删除0.10%、 0.16%、 0.94%和1.00%含气率的吸光度数据, 用剩下的12组含气率吸光度数据进行光谱相似性测度分析。选取的预处理方法包括: 归一化、 中心化(center), 标准化(autoscaling), S-G卷积平滑, 多元散射校正(MSC), 标准正态变量变换(SNV), 一阶导数和二阶导数。选取主成分分析法即PCA(principal component analysis), 提取有效的光谱特征信息, 去掉冗余信息。
光谱相似性测度即根据某一特定的光谱特征及对应的度量函数判定测试光谱与参考光谱之间的相似性并根据相似性的大小对测试光谱的类别属性进行划分[1]。本工作选取的光谱相似性测度方法包括以下四种:
2.2.1 欧氏距离
欧式距离是最常见且较易于理解的一种距离计算方法。在欧几里得空间中连接两个被测样本之间的直线距离即欧氏距离, 其计算公式为式(3)
式(3)中: xi为向量x的第i个元素; yi为向量y的第i个元素; n为向量的维数。
需要注意的是, 欧氏距离的计算结果与样本的尺度有关。当样本各个维度方向的尺度不一致时, 欧氏距离的计算效果较差, 一般先使用归一化方法处理原始实验数据; 同时, 欧氏距离的计算结果与样本维度也有关, 当维数过高时, 欧氏距离的作用会越小。因此对原始数据进行预处理后, 一般使用PCA方法减小样本维度大小, 以获得较好的欧氏距离的计算效果。测试光谱与参考光谱之间的欧式距离越小, 其相似性越高。
2.2.2 马氏距离
马氏距离(Mahalanobis distance)是基于总体样本基础上的, 计算该点与这个总体样本分布之间的距离。并且该方法与尺度无关, 保证原始数据的各个维度尺度保持一致。马氏距离的计算公式为式(4)
式(4)中: X为样本向量; Y为样本向量; S为协方差矩阵。
在计算马氏距离之前, 需要保证数据的样本个数大于数据维数。一般使用PCA方法对原始数据进行降维处理, 保证总体数据协方差矩阵的逆矩阵存在, 进而进行马氏距离的计算。测试光谱与参考光谱之间的马氏距离越小, 其相似性越高。
2.2.3 光谱角余弦
光谱角余弦通常用来描述光谱向量在方向或形状上的差异。两个光谱向量之间的夹角余弦值计算公式为
式(5)中: θ 为光谱向量之间的夹角即光谱角, θ ∈ [0, π /2]; Xk为测试光谱中第k个波长的吸光度值; Yk为参考光谱中第k个波长的吸光度值; m为波长数。
光谱向量之间的夹角余弦值为[0, 1], 且夹角余弦值越大, 光谱之间的相似性越高。
2.2.4 光谱相关系数
光谱相关系数用于描述光谱向量之间的相关程度大小。光谱向量的吸光度均值用式(6)和式(7)计算
则光谱向量之间的相关系数rij的计算公式为
式(8)中:
光谱向量之间相关系数绝对值的大小取[0, 1], 并且该绝对值越大光谱之间越相似。
2.2.5 含气率大小的预测
使用光谱相似性度量函数, 分别求取测试样本集D(k× m)中每一个样本Di(i=1, 2, …, k)与其对应的参考光谱集Ri中各个样本的相似性测度值, 并设置合理的相似性测度值的阈值, 最后将与测试样本Di相似性最高的参考光谱所对应的含气率值, 作为样本Di的含气率大小。
图5所示为应用光谱相似性测度分析进行含气率大小预测的流程图。
采用光谱相似性测度方法得到油气两相流的含气率大小的预测值后, 需要利用一些指标对预测的效果进行评价, 选取的评价指标为相对误差, 计算公式为
式(9)中: δ 为相对误差; ε p为预测含气率值; ε 为计算含气率值即含气率的真实大小。
相对误差δ 越小, 表明含气率的预测效果越好。本工作中由于含气率的采样间隔为0.06%。故最小的δ 为6%。
含气率大小测定的不确定度来源主要包括以下三个方面: 重复测量引起的A类不确定度UArel、 测量整个管道体积时使用量筒量取润滑油引起的不确定度UB1rel和制备含气率样品时使用移液枪移除部分润滑油引起的不确定度UB2rel。选择室温21.0 ℃下0.10%~1.00%含气率范围中的某个含气率大小重复测量10次, 使用的500 mL量筒的最大允许误差为± 2.5 mL, 使用的1 mL移液枪的最大允许误差为± 5× 10-3 mL, 经计算, 含气率的测定结果ε r可表示为
在近红外光下, 分别使用欧氏距离、 马氏距离、 光谱角余弦、 光谱相关系数等四种光谱相似性测度方法进行含气率预测的结果如表所示。从表1中可以看出, 在近红外波段, 油气两相流的含气率大小信息主要体现在光谱向量的方向或形状差异这一光谱特征。
![]() | 表1 近红外波段下基于不同方法的含气率预测结果 Table 1 The prediction of void fraction based on different methods in near-infrared band |
可见光波段下基于四种光谱相似性测度方法的含气率预测对比结果如表2所示。可以看出, 在该波段以新润滑油为实验对象时, 油气两相流的含气率大小信息所表现的光谱特征为光谱之间的直线距离大小。
![]() | 表2 可见光波段下基于不同方法的含气率预测结果 Table 2 The prediction of void fraction based on different methods in visible band |
紫外波段下基于四种光谱相似性测度方法的含气率预测对比结果如表3所示。从表中可以看出, 在紫外波段, 油气两相流的含气率大小信息主要体现在光谱向量的方向或形状差异这一光谱特征上。并且使用光谱角余弦方法的含气率预测最大相对误差比近红外和可见光波段的都小, 这表明油气两相流的含气率信息在紫外光波段范围(193.5~413.8 nm)与光谱向量的方向或形状差异这一光谱特征的联系更为紧密。
![]() | 表3 紫外波段下基于不同方法的含气率预测结果 Table 3 The prediction of void fraction based on different methods in ultraviolet band |
为了研究光谱相似性测度方法对于不同温度以及流速下的油气两相流含气率预测的适应性, 选取5个温度: 室温(21.0 ℃)、 30.0 ℃、 35.0 ℃、 40.0 ℃和43.6 ℃, 2个流速: 5.1和7.7 m· min-1共10个工况下的油气两相流的吸光度数据进行光谱分析。
依据3.1节分析结果, 对近红外波段下均采用标准化预处理结合光谱角余弦方法, 可见光波段下均采用归一化结合欧氏距离方法, 紫外波段下均采用标准化预处理结合光谱角余弦方法, 分别对10个工况的含气率大小进行预测并对比, 其结果如表4, 表5和表6所示。
![]() | 表4 近红外波段不同工况下的含气率预测效果对比 Table 4 The comparison of void fraction prediction based on different conditions in near-infrared band |
![]() | 表5 可见光波段不同工况下的含气率预测效果对比 Table 5 The comparison of void fraction prediction based on different conditions in visible band |
![]() | 表6 紫外波段不同工况下的含气率预测效果对比 Table 6 The comparison of void fraction prediction based on different conditions in ultraviolet band |
可以看出, 在近红外波段, 两相流的流速固定不变时, 随着温度的升高, 含气率预测的最大相对误差在减小; 其中温度为最高温度43.6 ℃时, 流速7.7 m· min-1下的含气率预测最大相对误差为24%, 小于流速5.1 m· min-1下的30%。同样的, 除了室温21.0 ℃, 在其他4个温度, 随着两相流流速增加, 最大相对误差均在减小; 在可见光波段, 流速为5.1 m· min-1时, 温度30.0和40.0 ℃下的含气率预测最大相对误差为18%, 小于其他三个温度的最大相对误差; 流速为7.7 m· min-1时, 温度43.6 ℃下的含气率预测最大相对误差为12%, 小于其他四个温度的最大相对误差。除了温度30.0和40.0 ℃, 在其他3个温度, 随着两相流流速增加, 含气率预测的最大相对误差均在减小; 紫外波段下流速为5.1 m· min-1不变时, 温度30.0 ℃的含气率预测的最大相对误差为6%, 小于其他四个温度的最大相对误差; 流速为7.7 m· min-1不变时, 温度30.0和43.6 ℃的含气率预测的最大相对误差为18%, 小于其他三个温度的最大相对误差。除了43.6 ℃, 在其他4个温度, 随着两相流流速增加, 含气率预测的最大相对误差均在增大。由此可见, 温度对三个波段下的油气两相流含气率预测效果均有影响。尤其是在室温(21.0 ℃)时, 预测效果相对较差。
(1)弯管处的油气两相流在0.10%~1.00%含气率范围下, 两相流的流型主要为泡状流, 气泡大小均匀且较小。随着含气率增加, 气泡数量增多, 形成少数较大尺寸的气泡。
(2)弯管处的油气两相流在近红外、 可见光和紫外波段均有吸收, 吸收强度与含气率的大小有关。然而原始光谱数据由于引入较多无用信息和光散射信息, 使得波长变量和吸光度之间的相关性降低。
(3)油气两相流的吸光度数据中存在的主要问题为数据中心不明显和散射噪声信息较多。在近红外波段, 对于新润滑油, 通过使用数据增强方法如中心化和标准化以及多元散射校正(MSC)方法, 使得含气率预测的最大相对误差相对于原始光谱从48%降低到36%。
(4)在紫外波段(193.5~413.8 nm), 使用光谱角余弦方法的含气率预测效果相较于欧氏距离、 马氏距离、 光谱相关系数方法的效果较好, 最大相对误差仅为6%, 在近红外波段(1 000.6~2 473.2 nm)预测的最大相对误差为36%; 这表明油气两相流的含气率信息在紫外波段范围与光谱特征的联系更为紧密, 并且油气两相流的含气率大小信息主要体现在光谱向量的方向或形状差异这一光谱特征。在近红外和可见光波段, 流速一定时, 含气率预测的最大相对误差随着温度的升高而有减小的趋势; 并且温度一定时, 随着流速增加, 含气率预测的最大相对误差有减小的趋势; 紫外波段下温度对含气率没有特定的影响规律, 随着两相流流速的增加, 含气率预测的最大相对误差有增大的趋势。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|