0 引 言
随着现代工业及城市化建设的迅速发展, 化学气体带来的环境污染问题和安全问题日益严峻, 因此, 对人类生产生活中各种气体进行准确的识别和浓度监测具有重要意义
[1‐3]。金属氧化物半导体(Metal Oxide Semiconductor, MOS)气体传感器因具有灵敏度高、 响应速度快、 制作简单以及成本低廉等优点
[4], 被广泛应用于食品检测
[5]、 农业
[6]、 工业生产
[7]、 医疗诊断
[8]等领域。但MOS传感器的低选择性和交叉敏感问题限制了其发展及应用。提高选择性的常用方法是气体传感器阵列化, 通过不同敏感材料的多个传感器来构建电子鼻
[9‐11]。
然而, 传统电子鼻通常在恒温模式下工作, 仅靠单一维度的稳态响应电阻难以有效提高系统识别准确率。要获取更多气体分子特征信息需要增加传感器种类和数量, 这会导致制作成本、 体积和功耗大幅提高
[12]。尽管通过动态测试可以得到传感器响应恢复过程的更多特征信息, 但在实际现场检测中, 传感器长期暴露在复杂气体环境下, 动态响应过程难以获取。
因此, 研究人员采用温度调制技术
[13], 即利用气体传感器的温度依赖特性
[14], 通过施加周期性变化的加热电压, 获得传感器在不同温度下的动态响应, 从而构建虚拟传感阵列。与恒温测试的传感器静态响应信号相比, 温度调制模式下的信号蕴含着更丰富的气体响应信息, 结合模式识别技术可提高气体的识别率, 进而提高传感器的选择性。尤其随着微型电子机械系统(MEMS)气体传感器的成熟, 其功耗低、 热电响应速度快的特点为温度调制方法研究提供了有利的手段
[15]。目前的研究主要集中在气体传感器的温度调制模式优化
[16]和模式识别算法方面
[17]。
Krivetskiy等
[18]基于单个SnO
2传感器进行三角波调制, 结合人工神经网络对8种气体进行选择性检测, 分类识别错误率为13.2%; Meng等
[19]提出锯齿形矩形波温度调制方法, 有效提高了传感器的选择性。由于气体传感器的响应是热激发过程与表面反应过程共同作用的结果, 且与待测气体密切相关的表面反应速度远远低于热激发的变化速度
[20], 目前关于温度调制的研究报道获得的响应主要是一系列准稳态响应, 并未涉及快速温度调制模式下丰富的瞬态信息。
基于此, 本文采用快速脉冲温度调制技术, 在脉冲加热间隙的低温点进行信号检测, 提取表面反应过程相关信息, 可有效消除热激发对表面反应过程的影响, 为后续气体准确识别提供丰富的特异性信息。本文构建基于单个气体传感器的虚拟阵列并优化快速温度调制参数, 研究调制信号的处理及特征提取方法, 构建分类和浓度预测模型, 实现单个气体传感器对多种气体的准确识别和定量分析。该方法提高了单个传感器的选择性, 同时有效降低了功耗, 在高精度、 低成本、 低功耗现场检测方面具有应用潜力。
1 实 验
1.1 实验装置
气体传感器的测试装置如
图 1 所示, 主要由传感器、 密闭气室、 加热调理模块和信号采集模块组成。
传感器选用苏州慧闻纳米科技有限公司商用MEMS气体传感器系列的SMD1012氢气传感器, 并悬置于2 500 mL的玻璃气室中部。传感器加热端接入加热调理模块, 由单片机控制模块调制电压波形, 为传感器提供需要的工作温度, 电压调节范围为0~1.8 V; 测量端通过分压电阻形式接入信号采集模块, 由单片机进行数据采集读取, 并传输至上位机, 记录传感器电阻响应数据。
1.2 测试过程
选择H
2、 H
2S、 NH
3和C
2H
5OH作为目标测试气体。通过质量流量控制器(Mass Flow Controller, MFC)调节通入气室的标准气体和洁净空气的体积, 稀释获得设置浓度的H
2、 H
2S、 NH
3测试气体, C
2H
5OH气体则由蒸发法配制得到。在测试之前, 首先对传感器预热24 h以上, 以确保传感器响应的稳定。调制电压波形选择正弦调幅脉冲波形, 总周期时间为30 s, 调制电压波形及局部放大图如
图 2(a) 所示, 在该电压的驱动下, 传感器的工作温度经历高温-低温的循环变化过程。
首先进行正弦调幅脉冲波形参数优化, 研究不同频率和电压下的调制模式。根据传感器的常温工作特性, 选择1~1.8 V、 1.2~1.8 V、 1.4~1.8 V、 1.6~1.8 V 4种电压调制范围, 调制频率分别设计为10、 20、 30、 40和50 Hz, 共20种温度调制模式。分别在每种调制模式下进行4种气体测试, 样本气体浓度使用体积分数表示, H2测试浓度为20×10-6、 40×10-6、 60×10-6、 80×10-6和100×10-6; H2S测试浓度为0.1×10-6、 0.2×10-6、 0.3×10-6、 0.4×10-6和0.5×10-6, NH3测试浓度为20×10-6、 40×10-6、 60×10-6、 80×10-6和100×10-6; C2H5OH测试浓度为100×10-6、 200×10-6、 300×10-6、 400×10-6和500×10-6, 每一浓度重复性测试3次。
随后在优化的调制模式下进行气体定性定量分析, 对10种不同浓度的测试气体进行分类和浓度预测。具体目标测试气体浓度设置见
表 1。重复4组实验, 共获得160个响应信号样本, 将其中3组样本数据作为训练模型的训练集, 剩余1组样本数据作为模型预测和评估的测试集。
气体测试流程具体为: 首先将传感器置于充满干燥洁净空气的气室中清洗恢复20 min, 然后从低浓度开始依次通入设置浓度间隔相应体积的标准测试气体, 进行不同浓度的测量循环。每次通入测试气体后, 待样品在气室中扩散均匀、 气敏电阻响应稳定, 记录原始响应信号。
1.3 数据预处理与特征提取
在脉冲温度调制模式下, 气敏电阻的阻值会随工作温度周期性变化而变化, 同时气体与敏感材料的表面反应过程也会引起电阻的变化。因此, 呈现的响应图谱反映了工作温度与气体种类浓度的综合影响, 如
图 2(b) 所示。在此调制模式下, 电阻响应信号的上包络线是加热间隙(低温点)测得的膜电阻曲线, 主要反映了薄膜表面吸附状态的变化; 而下包络线是加热时(高温点)测得的, 该变化是由热激发过程与表面反应过程共同作用的结果
[20]。因此, 上下包络线蕴含了关于目标测试气体特征的重要信息, 可以从时域角度提取其上下包络线, 将周期振荡的原始电阻响应信号转化为两组低频响应信号。
本文首先基于极大极小值的寻峰算法, 自动识别提取响应图谱中的上下包络线, 并应用Savitzky-Golay滤波器进行降噪处理, 以进一步提高数据的可靠性和准确度。随后, 对上下包络线进行归一化处理, 减小基线偏移带来的误差影响。接着, 选用“db4”作为小波基, 对每个周期的气体响应数据进行3层离散小波变换(Discrete Wavelet Transform, DWT), 提取分解后的低频小波系数作为气体特征向量。
图 2(c) 显示了响应信号包络线经3层离散小波变换后的低频系数, 可以看出提取的小波系数数量仍较多。因此, 进一步选择主成分分析(Principal Component Analysis, PCA)以及线性判别分析(Linear Discriminant Analysis, LDA)对小波系数进行特征降维处理, 以提高后续算法的运行效率。
1.4 模式识别系统
基于上述多维响应数据, 本文采用模式识别算法对非线性的样本数据进行自动处理和判读, 实现对多种目标测试气体的定性及定量分析。设计的模式识别系统主要分为两个步骤: 利用分类模型对H
2、 H
2S、 NH
3和C
2H
5OH进行分类; 根据分类结果选择相应的浓度估计器来预测气体浓度, 该流程的详细实现如
图 3 所示。
首先, 选用支持向量机(Support Vector Machine, SVM)作为分类器, 构建四元分类模型。SVM通过最大化边缘距离实现结构风险最小化原则, 因而具备优秀的模型泛化能力和在小样本数据集上的分类性能, 是解决复杂分类问题的理想选择
[21]。在模型构建过程中, 首先对样本训练集进行预处理、 特征提取及特征降维, 然后利用网格搜索法并结合5折交叉验证, 优化分类模型的参数。在确定最优模型参数后, 使用完整的训练集样本数据重新训练分类模型, 并利用测试集来评估模型的性能。
随后, 引入了支持向量回归(Support Vector Regression, SVR)、 随机森林回归(Random Forest Regression, RF)和反向传播神经网络(Backpropagation Neural Network, BPNN) 3 种不同的回归策略, 针对不同种类气体分别构建浓度估计模型, 并采用均方根误差(Root Mean Square Error, RMSE)作为模型评价指标。该过程与分类模型的构建相似, 基于训练集数据, 采用网格搜索法结合5折交叉验证, 优化回归模型参数。通过分类模型的输出结果, 为测试集样本选择相应气体的浓度估计模型, 进而对模型的预测精度进行评估。
2 实验结果与分析
2.1 参数优化
利用正弦调幅脉冲波形对气体传感器进行温度调制, 在20种调制模式下, 分别对所列不同种类及浓度的气体进行测试, 研究不同的调制电压和调制频率参数对气体检测的影响。
图 4 为20种不同调制参数组合下, 4种气体样本响应数据提取的上包络线。
随着调制电压的增加, 即传感器低温检测点的温度升高, 不同气体的电阻响应包络线形状趋于相似, 响应特征之间的差异减小。这是因为通过扩展脉冲高温点至多个温度点并在固定的低温点进行检测的策略, 能够利用不同气体在不同温度下的吸附和脱附特性进行区分, 从而有效地消除了热激发对表面反应过程的干扰。然而, 随着低温点温度的升高, 检测温度范围变窄, 导致多种气体可能出现相似的表面反应特性, 进而降低了传感器的选择性。此外, 由于敏感薄膜与基底材料的膨胀系数不同, 传感器在高低温间周期性变化时, 过低的低温点电压会导致敏感薄膜产生剧烈的机械振动, 不仅影响薄膜与电极的接触性能, 还会降低传感器的稳定性和使用寿命。
随着调制频率的提高, 不同气体的电阻响应包络线之间的差异在减小。尽管微热板式传感器的热响应时间在毫秒级别, 但在较高频率下敏感电阻表面化学反应可能无法充分跟随温度变化的速率, 导致表面反应过程与温度变化之间产生延迟, 减小了不同气体间响应特征的区分度。
综合上述分析, 本文最终选择1.0 V、 10 Hz作为检测不同气体的最佳调制参数。
2.2 定性识别
首先, 初步评估了单个传感器对不同种类气体的响应特征差异。基于最佳调制参数下测得的样本集, 经预处理和DWT变换后, 选择非监督的PCA和有监督的LDA分别对特征向量进行降维处理。
图 5(a),
图 5(b) 展示了PCA前两个主成分构成的特征空间中样本的分布情况, 其中红色描边的点表示测试集的样本点, 3个椭圆区域分别代表
、
、
的置信区间。
结果显示, 基于上包络线特征的样本具有较好的分离度, 尽管H2和NH3样本在二维的投影空间中存在重叠区域, 但在三维或更高维度空间中, 仍能实现对气体样本的有效区分; 相比之下, 下包络线特征的样本点除了H2S外, 其它样本点高度重叠, 难以区分。这是因为上包络线特征与气体的表面反应过程密切相关, 其形状差异主要由测试气体的种类和浓度引起; 而下包络线特征则受热激发的影响较大, 与温度变化紧密相关, 不同气体特征之间的差异信息淹没在热调制信息中。
此外, 由于气体浓度差异的影响, 同一类别的样本点分布相对分散, 类内距离较大。得益于LDA最大化类间距离和最小化类内距离的特性
[22], 在
图 5(c),
图 5(d) 的LDA特征空间中, 样本点展现出更明显的可分性, 同时同类样本分布更加聚拢。因此, 本文选择上包络线特征并结合LDA构建分类模型, 在尽可能保留不同气体特征信息的同时, 减小气体浓度差异对检测分类性能的影响。
进一步利用测试集验证分类模型的性能。
图 6 为SVM对测试集预测结果的混淆矩阵, 对角线表示正确分类的样本, 结果表明SVM的分类准确率为100%, 能够完美识别4种目标测试气体。
2.3 定量分析
在准确识别出测试集样本的气体种类之后, 分别对每种气体的浓度进行预测。在
图 5(a) 和
图 5(b) 的PCA二维特征空间中, 同类气体样本点随浓度的变化沿特定矢量方向分布, 而相同浓度的同类气体样本点则垂直于该矢量方向分布, 尤其在下包络线特征空间中该现象更为明显。这表明对于同种气体, 响应包络线的差异主要是由浓度不同引起的幅度变化所致, 即使是相同浓度的样本, 也会因基线波动和测量误差而受到一定程度的影响。因此, 本文采用PCA对每种气体的下包络线特征进行单独处理, 此时的方差主要与浓度有关。各主成分与浓度呈单调相关, 将其作为浓度预测的输入向量, 并选择SVR、 RF和BPNN进行浓度预测。
以SVR为例, 4种目标气体的预测浓度与真实浓度关系如
图 7 所示, 其中对角线代表预测值与真实值一致。对于4种目标气体, 绝大多数样本点的预测浓度都紧邻对角线分布, 这部分误差主要是由不同组测量实验中传感器基线的偏移引起的幅值误差。此外, 预测浓度误差较大的点主要出现在浓度分布两端, 这是由于样本分布的两端数据相较中心区域更为稀缺。
表 2 展示了SVR、 RF和BPNN预测浓度的均方根误差RMSE比较, 可以看出3种算法的预测性能总体接近, 均受传感器基线偏移等因素影响, 但PCA结合SVR的回归模型稍优于其它两者, 对H
2、 H
2S、 NH
3和C
2H
5OH预测浓度的均方根误差分别为19.5×10
-6、 3.7×10
-6、 0.2×10
-6、 19×10
-6。
3 结 论
本文基于单个气体传感器, 采用温度调制技术构建了虚拟阵列传感系统。研究了动态调制信号的处理及特征提取方法, 优化了温度调制模式; 采用支持向量机实现了对H2、 H2S、 NH3和C2H5OH 4 种目标测试气体的准确识别, 对比了支持向量回归、 随机森林回归、 反向传播神经网络3种算法对气体浓度的预测性能, 最高预测误差分别为19.5×10-6、 3.7×10-6、 0.2×10-6、 19×10-6。本文方法有效提高了单个气体传感器的选择性, 同时降低了成本和功耗, 满足现场检测的快速响应及准确可靠的气体检测需求。
国家重点研发计划资助项目(2021YFB3201300)