摘要:气体浓度检测一直是十分重要的工作,它与人类生活、环境变化和工业生产都息息相关,尤其是有毒有害气体的检测。随着光谱学的大力发展,可调谐半导体激光光谱技术(TDLAS)已被广泛应用于工业生产、环保监测、气象探空等领域,它具有灵敏度高、响应速度快、实时监测以及优秀的便携性等优点,成为了气体检测的重要技术之一。众所周知,传感器在出厂前需要标定,标定的好坏关系到传感器的测量精度,影响传感器的实际测量。对基于TDLAS技术的直接光谱标定算法进行了研究,推导了适用于大量程范围下的气体浓度与透射率对数的关系并进行了简化处理,提出了1/c拟合标定算法,它是用气体浓度倒数与透射率对数倒数的关系来进行拟合标定的算法,再通过拟合标定曲线反推气体浓度与透射率对数的关系。搭建了TDLAS水汽检测平台,以VaisalaHMT377在线湿度检测仪的水汽浓度值为标准值,对0.7%~50%VOL的水汽浓度进行了标定,比较了直接拟合算法和1/c拟合算法的结果和相对误差。实验结果表明:采用直接拟合算法对气体浓度与透射率对数曲线拟合,一次函数、二次函数和五次函数拟合相关系数为0.9468,0.9967和0.9999,均方根误差为0.0312,0.0078和0.0016,最大相对误差超过100%,存在严重误差;采用1/c拟合标定算法,气体浓度倒数与透射率对数倒数的拟合曲线呈一次函数线型,与理论相符,相关系数为0.9996,均方根误差为0.4901,对拟合曲线进行反推得到气体浓度与透射率对数拟合曲线,其相关系数为0.9999,均方根误差为0.0015,整个浓度范围内最大相对误差不超过4%,说明1/c拟合算法的有效性,提高了测量精度,扩大了标定量程范围。
气体浓度检测一直是十分重要的工作,它与人类生活、环境变化和工业生产都息息相关,尤其是有毒有害气体的检测[1,2]。近年来随着光谱技术的大力发展,可调谐半导体激光光谱技术(tunablediodelaserabsorptionspectroscopy,TDLAS)已逐渐发展成熟,它具有灵敏度高、响应速度快、实时监测以及优秀的便携性等优点,成为了气体检测的重要技术之一[3,4,5]。
标定是各类气体传感器出厂前必须要进行的一道工序,标定结果的好坏直接关系到传感器的最终测量精度[6]。基于TDLAS技术的直接吸收光谱标定算法主要是利用特定波长的激光穿过气体样品后,测量光强的透射率对数与标准浓度进行拟合,将实测的透射率对数代入拟合关系式来计算当前环境下的气体浓度[7,8]。根据Beer-Lambert定律可知,气体浓度与透射率对数是一次函数关系,标定时通常会采用一次函数进行拟合,但这种情况忽略了气体浓度对吸收系数的影响,尤其在气体浓度较高时,这会导致气体浓度与透射率对数的关系不是严格的一次函数关系,原有的标定算法会带来严重误差,影响测量精度。
本工作对TDLAS的直接吸收光谱标定算法进行研究,考虑了气体浓度对吸收系数的影响,推导了完整的气体浓度与透射率对数的关系式,提出了采用气体浓度倒数与透射率对数倒数的关系曲线进行拟合标定的算法;搭建了基于TDLAS水汽检测平台,对0.7%~50%浓度范围的水汽浓度进行了标定,对直接拟合标定算法和本算法作了比较分析。
1、理论分析
传统的TDLAS检测理论依据Beer-Lambert定律,
It=I0exp[−α(ν)cL] (1)Ιt=Ι0exp[-α(ν)cL] (1)
式(1)中,It为穿过待测气体后的透射光光强,I0为入射光强,α(ν)为吸收系数,c为待测气体的浓度,L为光吸收路径长度。
由式(1)可得,气体浓度与透射率对数Ratio的关系式
c=−1α(ν)Lln(ItI0) (2)Ratio=−ln(ItI0)=α(ν)cL (3)c=-1α(ν)Lln(ΙtΙ0) (2)Ratio=-ln(ΙtΙ0)=α(ν)cL (3)
根据式(2)可知,气体浓度与透射率对数是一次函数关系,但气体浓度较高时,吸收系数α(ν)会受到影响[9,10]。
α(ν)=NS⋅gL(ν)=NSγL/π(ν−ν0)2+γ2L=NSπγL[(ν−ν0)2/γ2L+1]=α0[(ν−ν0)2/γ2L+1] (4)α(ν)=ΝS⋅gL(ν)=ΝSγL/π(ν-ν0)2+γL2=ΝSπγL[(ν-ν0)2/γL2+1]=α0[(ν-ν0)2/γL2+1] (4)
其中,gL(ν)为洛伦兹线型函数,S为吸收谱线强度,N为单位体积内气体分子总数,α0=NS/πγY,为气体吸收峰中心位置(ν=ν0)的吸收系数,γL为气体吸收谱线半宽高,它与气体浓度大小有关,表达式如式(5)
γL(p,T)=pp0(T0T)n[(1−c)γair(p0,T0)+cγself(p0,T0)] (5)γL(p,Τ)=pp0(Τ0Τ)n[(1-c)γair(p0,Τ0)+cγself(p0,Τ0)] (5)
式(5)中,n为温度系数,p0和T0分别为标准气压和标准温度,γair是空气吸收谱线半高宽,γself是待测气体吸收谱线半高宽,由式(5)可知,γL是待测气体的线宽和背景空气的线宽共同决定。对于一般对痕量气体的含量进行检测,即待测气体的浓度c≪1,通常忽略掉后项,简化为γL(p,L)=γair(p0,T0)pp0(T0T)nairγL(p,L)=γair(p0,Τ0)pp0(Τ0Τ)nair。而对于待测气体浓度很高的环境,c≪1的条件不再成立,因此必须考虑c对线宽的影响。
根据式(3),当v=v0,即在吸收峰的中心位置,则
Ratio=α(ν0)cL=NScLπγL=Ks(T)S0N0pp0T0TcLπpp0(T0T)n[(1−c)γair(p0,T0)+cγself(p0,T0)] (6)Ratio=α(ν0)cL=ΝScLπγL=Κs(Τ)S0Ν0pp0Τ0ΤcLπpp0(Τ0Τ)n[(1-c)γair(p0,Τ0)+cγself(p0,Τ0)] (6)
式(6)中,Ks(T)是吸收谱线强度关于温度的修正系数,S0是标准气压温度下的吸收谱线强度,N0是标准气压温度下的单位体积分子数,对式(6)化简得到,气体浓度c的表达式。
c=[k0Lγair(p0,T0)Ratio(T0T)n−11Ks(T)+1−γself(p0,T0)γair(p0,T0)]−1 (7)c=[k0Lγair(p0,Τ0)Ratio(Τ0Τ)n-11Κs(Τ)+1-γself(p0,Τ0)γair(p0,Τ0)]-1 (7)
式(7)中,k0=S0N0/π,图1是气体浓度c与透射率对数Ratio的模拟关系曲线图。
图1气体浓度与透射率对数的模拟关系曲线图
Fig.1ThesimulationcurveofgasconcentrationvsRatio
根据图1可知,两者关系并非线性,导致标定过程变得复杂,因此对修正后的气体浓度的表达式进行简化。
令y=1c,x=[Ratio(T0T)n−11Ks(T)]−1y=1c,x=[Ratio(Τ0Τ)n-11Κs(Τ)]-1,则式(7)简化为
y=k0Lγair(p0,T0)x+1−γself(p0,T0)γair(p0,T0) (8)y=k0Lγair(p0,Τ0)x+1-γself(p0,Τ0)γair(p0,Τ0) (8)
式(8)即是考虑气体浓度对吸收影响后的气体浓度和透射率对数的关系式,其中气体浓度倒数与透射率对数倒数是一次线性函数关系。测量时,通过式(8)对气体浓度倒数与透射率对数倒数进行拟合,将透射率对数倒数代入式(8),求得气体浓度倒数,再进行一次求倒数运算即可得到实测浓度,本文将这种标点算法称为1/c拟合标定算法。这种标定算法大大简化了标定关系式的复杂程度,优化拟合的相关系数,扩大了标定量程,提高了测量精度。
2、实验部分
根据Beer-Lambert定律可知,气体吸收程度除了与待测气体浓度相关,吸收谱线强度和吸收路径长度对吸收峰的大小也有影响。为了尽可能地扩大标定量程范围,避免在低浓度下气体吸收峰信号太弱而无法检测以及高浓度下气体吸收峰过强而出现截顶的情况,本实验选取了中心波长为1368.597nm的水汽吸收峰,吸收强度为1.8×10-20cm-1·(molec·cm-2)-1,单次反射气室的光路长为55mm的开放式气室。图2是1365~1370nm的水汽吸收峰的光谱。
本实验搭建的TDLAS水汽检测系统包括激光器、系统控制单元和气室探头三个部分。激光器选用了波长1368nm的DFB半导体激光器,输出功率为3mW。系统控制单元通过单片机产生一个30Hz的低频锯齿波信号,经过电流驱动模块后输入到激光器中,通过对激光器的电流调制使出射光的波长覆盖选定的水汽吸收峰(1368.597nm),出射光经过准直镜后在气室进行反射吸收,光电探测器检测到反射吸收后的光信号并转化为电信号,最后对信号进行降噪放大处理,在PC端对信号进行处理分析。标定时,选用VaisalaHMT337在线湿度检测仪作为标准值,将TDLAS水汽检测系统的气室探头和VaisalaHMT337的湿度探头放入北京雅士林DHS-100恒温恒湿箱内,通过调节恒温恒湿箱来控制箱内的水汽浓度。整个实验装置测试图如图3所示。
图21365~1370nm范围内的水汽吸收峰的光谱
图3实验装置整体示意图
表1不同浓度下水汽吸收光谱透射率对数结果
3、结果与讨论
通过控制DHS-300恒温恒湿箱,设定18个不同水汽浓度测试点,浓度范围在0.7%~50%VOL,对水汽的直接吸收光谱进行了提取,对直接拟合和1/c拟合结果进行了比较。
3.1直接拟合
实验获取了18组不同水汽浓度下的透射率对数结果,如表1所示。
对表1中的数据进行拟合,选取了一次函数、二次函数以及五次函数,对三种拟合函数结果进行相关系数和均方根误差分析,其结果如表2所示,气体浓度和透射率对数曲线和对应的拟合曲线如图4所示。
表2三种拟合曲线的相关系数和均方根误差
图4气体浓度和透射率对数曲线、对应的拟合曲线和拟合曲线的相对误差
由图4所示,在大范围的浓度情况下气体浓度与透射率对数曲线不是一次函数关系,与理论分析相符合,其整体曲线趋势与图1中的模拟结果相近。表2的拟合结果也验证了一次函数拟合的结果不好,相关系数最小,均方根误差最大;二次函数和五次函数拟合的结果接近,相关系数都大于0.99,五次函数的相关系数达到了0.9999,均方根误差0.0016;从相对误差曲线上看,一次函数和二次函数拟合曲线在低浓度下的误差较大,超过了100%,测量完全失准。
3.21/c拟合
对表1中的气体浓度和透射率对数进行倒数处理,气体浓度倒数与透射率对数倒数曲线以及反推的气体浓度与透射率对数拟合曲线如图5所示。
图5中气体浓度倒数与透射率对数倒数曲线与理论推导的式(8)相符合,结果为一次函数线型,采用一次函数对曲线进行拟合,其拟合系数为0.9996,均方根误差为0.4901。将拟合后的曲线进行反推得到气体浓度与透射率对数曲线,其结果如图5所示,拟合系数为0.9999,均方根误差0.0015,在整个浓度范围内相对误差不超过4%。
图5气体浓度倒数与透射率对数倒数曲线以及反推后的拟合曲线
3.3两种算法的验证对比
为了对两种算法进行验证,实验选取了10%,20%,30%,40%,50%,60%和70%的标准浓度的水汽,对两种算法的结果进行了对比。其结果如表3所示。
表3多项式拟合和1/c算法的拟合结果
两种算法的相对误差结果曲线如图6所示,根据对比可知,二次拟合算法的相对误差超过了±25%,五次拟合的相对误差不超过±5%,采用1/c拟合标定算法的相对误差则不超±3%。从相对误差曲线趋势上看,多项式拟合算法的结果随着标准浓度范围的扩大,相对误差也在增加,而1/c拟合标定算法的相对误差结果则比较稳定,相对误差较小。多项式拟合只能针对标定浓度范围内的结果保证较好的相关性,不具有特定曲线模型,因此实测浓度超出标定浓度范围时,多项式拟合的误差就会急剧增大;1/c拟合方法则是从原理出发,因此对超出标定浓度范围的结果具有比较好的预测,误差结果较小。
图6多项式拟合和1/c算法的相对误差结果
4、结论
基于TDLAS原理分析了气体浓度与透射率对数的关系,并考虑了高浓度对线宽的影响,推导并简化了两者的关系,提出了利用气体浓度倒数与透射率对数倒数的关系来进行拟合标定的方法。实验结果表明:采用直接拟合算法对气体浓度与透射率对数曲线拟合,一次函数、二次函数和五次函数拟合相关系数为0.9468,0.9967和0.9999,均方根误差为0.0312,0.0078和0.0124,最大相对误差超过100%,存在严重误差;采用1/c拟合标定算法,气体浓度倒数与透射率对数倒数的拟合曲线呈一次函数线型,与理论相符,相关系数为0.9996,均方根误差为0.4901,对拟合曲线进行反推得到气体浓度与透射率对数拟合曲线,其相关系数为0.9999,均方根误差为0.0015,整个浓度范围内最大相对误差不超过4%,说明此方法有效,提高标定精度,扩大了标定量程范围。
参考文献:
[1]贾良权,刘文清,阚瑞峰,等.中国激光,2015,42(7):312.
[2]刘立富,张涵,温作乐,等.激光与光电子进展,2015,52(11):50.
[3]陈昊,鞠昱,韩立.光谱学与光谱分析,2018,38(6):16.
[5]JiangJ,ZhaoM,MaGM,etal.IEEESensorsJournal,2018,(99):1.
[6]贾军伟,李伟,柴昊,等.红外与激光工程,2019,48(5):202.
[7]张志荣,董凤忠,吴边,等,光电子·激光,2011,22(11):1691.
[8]胡雅君,赵学红,张锐,等.光学学报,2013,(11):296.
鞠昱,陈昊,韩立,常洋,张学健.基于TDLAS技术适用于大量程范围的气体浓度标定算法[J].光谱学与光谱分析,2020,40(12):3665-3669.
基金:国家自然科学基金项目(41705141);北京市科技计划课题(Z181100000118005)资助.
分享:
提出了一种利用数码相机进行观测的新型分光计方案及其调节方法,并研制了实验样机.利用数码相机替代了传统分光计中的望远镜部分,可观察并记录相关光学现象.调节分光计时,首先利用数码相机镜头对焦于无穷远作为标准进行平行光调节,然后利用反射狭缝像的位置调节载物台与旋转轴的垂直,最后以载物台为基准调节数码相机及平行光管与旋转轴垂直.对新型分光计的特性分析及实验结果表明:该新型分光计具有调节难度低、测量误差小、能同时观测多条光谱线的特点。
2020-12-08长波红外光谱(8~14μm)是介于中红外波段和太赫兹波之间的重要电磁辐射,对应着地球表面常温目标物体的辐射波段和地球“第三大气窗口”,相对于短波和中波红外辐射,长波红外辐射受大气散射影响较小。因此,长波红外辐射在红外夜视、资源探测、精确制导、安防报警等科研及国防领域具有极其重要的应用,为国家综合实力的重要体现。
2020-12-08气体浓度检测一直是十分重要的工作,它与人类生活、环境变化和工业生产都息息相关,尤其是有毒有害气体的检测。近年来随着光谱技术的大力发展,可调谐半导体激光光谱技术(tunablediodelaserabsorptionspectroscopy,TDLAS)已逐渐发展成熟,它具有灵敏度高、响应速度快、实时监测以及优秀的便携性等优点,成为了气体检测的重要技术之一。
2020-12-08随着科学技术的不断发展进步,军事和镜头相关的工商业领域都越来越多的使用红外光学系统。尤其在军事方面,光学系统发展迅速,导致军事方面对光学系统的性能有越来越高的要求。通常情况下,摄远物镜的系统长度小于其焦距[1],所以焦距相同时,把红外物镜设计成红外摄远物镜就可以很大程度上减少系统的制作成本。
2020-12-05光刻工艺是集成电路制造中最重要、最关键的工艺步骤之一。随着半导体技术的飞速发展,图形越来越密集,特征尺寸越来越小,对光刻工艺分辨率的要求越来越高。光刻工艺中一个重要的性能指标是每个图形的分辨率。在先进的半导体集成电路制造中,为获得高集成度器件分辨率很关键。
2020-11-20激光加工是一种非接触、无污染、无磨损的加工工艺,它包含激光切割、激光打孔、激光标刻等加工方式。激光打标技术已经在众多领域占有很大的比例,国家根据激光标刻行业制定了一系列的标准。工艺参数的设定在激光标刻过程中对其加工质量和效率起到了重要作用[1]。本文以激光标刻文字、条形码为研究对象,通过优化汉字笔画打标顺序和圆弧代替直线加工,能够激光标刻质量和效率。
2020-10-19在我们最近的工作中[11],将RF白噪声和波长调制被同时添加到OA-ICOS系统中,提出了一种痕量气体检测的新方法(RF-WM-OA-ICOS)。与未受RF噪声干扰的OA-ICOS相比,RF-WM-OA-ICOS的探测极限提高了约6倍。本文将进一步地详细研究RF噪声扰动对WM-OA-ICOS系统的二次谐波信号的影响,并基于该技术建立了一套高精密的甲烷气体测量装置。
2020-09-09自动随着公路建设的不断发展,对公路建设的质量提出了更高的要求,需要构建公路磨耗层质量自动监测模型,结合大数据信息处理方法,进行公路磨耗层质量自动监测的信息处理,提高公路磨耗层质量监测的自动化水平[1],从而提高公路建设的质量,提高对公路磨耗层质量信息健康管理水平[2],采用大数据信息管理和优化融合调度技术,进行公路磨耗层质量自动监测系统开发设计,在信息化管理平台下,进行公路磨耗层质量自动监测方法改进,对相关的公路磨耗层质量自动监测模型研究具有重要意义。
2020-08-10激光线性结构光传感器具有结构简单,体积小等特点,广泛用于焊接的焊缝追踪,轨迹规划等,称为焊接过程自动化的关键技术之一。而在结构光条定位中,为了提高定位精度,提出了亚像素的概念,是指在相机得到的像素和像素之间还有几十到十几个微米的距离。而通常我们计算出的坐标都是正整数,这表示的是我们是在对像素进行操作,而亚像素计算出来的坐标是实数,能够极大地提高算法的精度。
2020-08-10强度关联成像,又称为鬼成像或量子成像[1,2],是一种全新的成像体制,具有非局域性、抗干扰能力强和超分辨等优点,在遥感成像、显微成像、医学成像和超分辨成像等领域都具有巨大的应用潜力,目前已经成为研究热点[3,4].压缩感知理论首先在数学领域提出来[5],随后压缩成像作为压缩感知的一个研究方向迅速引起广大研究者的兴趣.由于多光谱成像数据是三维的,系统一次只能探测多个谱段内一个空间维度上的数据,一般需要扫描获得二维空间上的图像.
2020-08-10我要评论
期刊名称:光学技术
期刊人气:607
主管单位:中国兵器工业集团公司
主办单位:北京兵工学会,北京理工大学,中国北方光电工业总公司
出版地方:北京
专业分类:科学
国际刊号:1002-1582
国内刊号:11-1879/O4
邮发代号:2-830
创刊时间:1975年
发行周期:双月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:0.735
影响因子:0.645
影响因子:1.369
影响因子:0.874
影响因子:0.385
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!