CN108981953B - 一种基于干涉调制原理的激光吸收光谱测温方法和*** - Google Patents

一种基于干涉调制原理的激光吸收光谱测温方法和*** Download PDF

Info

Publication number
CN108981953B
CN108981953B CN201811235909.XA CN201811235909A CN108981953B CN 108981953 B CN108981953 B CN 108981953B CN 201811235909 A CN201811235909 A CN 201811235909A CN 108981953 B CN108981953 B CN 108981953B
Authority
CN
China
Prior art keywords
data
absorption
absorption spectrum
laser
fitting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811235909.XA
Other languages
English (en)
Other versions
CN108981953A (zh
Inventor
徐立军
曹章
李远青
邱爽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201811235909.XA priority Critical patent/CN108981953B/zh
Publication of CN108981953A publication Critical patent/CN108981953A/zh
Application granted granted Critical
Publication of CN108981953B publication Critical patent/CN108981953B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明提供一种基于干涉调制原理的激光吸收光谱测温方法和***,所用元件包括激光发生单元,马赫曾德光纤干涉仪,光电探测单元和数据处理单元等。其中,激光信号产生后通过马赫曾德干涉仪,然后穿过待测气体,由光电探测器采集数据并上传至数据处理单元处理。从测量数据中提取吸收光强数据和相对波数与吸收光强数据的映射,拟合得到基线和吸收光谱参数的初始值,并通过迭代对其进行修正,根据修正的吸收光谱参数解算积分吸收率,由多条谱线的积分吸收率解算出被测气体温度。本发明光学***结构简单,无可动部件且可抑制背景辐射噪声,整体效果突出、适用范围广、可靠性高,具有广泛的应用前景。

Description

一种基于干涉调制原理的激光吸收光谱测温方法和***
(一)技术领域
本发明提出一种基于干涉调制原理的激光吸收光谱测温方法和***,特别涉及一种在测量***中复用为标准具和光强调制器的马赫曾德光纤干涉仪,将激光调制技术应用于温度测量领域。
(二)背景技术
基于激光吸收光谱的气体参数测量方法作为一种快速、灵敏度高的非接触测量方法得到了飞速的发展,在温度和物质浓度测量领域是一种有效的方法,被广泛应用在燃烧诊断、燃烧控制等方面,应用潜力巨大,前景可观。
典型的单光路吸收光谱的测量方法分为直接吸收法(DAS,Direct AbsorptionSpectroscopy)和波长调制法(WMS,Wavelength Modulation Spectroscopy)。两条吸收谱线的单光路激光吸收光谱技术用来获取流场沿激光路径的平均温度和组分浓度。如X.Liu和J.B.Jeffries在2006年于《应用物理B》(Applied Physics B)第82卷第3期469-478页发表的题为《用于测量燃气轮机排气温度的可调二极管激光传感器的研制》(Development ofa tunable diode laser sensor for measurements of gas turbine exhausttemperature)一文中,在用于发电机的大型工业燃气轮机中验证,沿激光路径的平均温度由近红外区两根水吸收谱线测量吸光度的比值转换得到,并正在开发多线策略来解决不均匀温度分布的问题。DAS方法广泛应用于低温低压测量。然而,如M.A.Bolshov和Yu.A.Kuritsyn在2015年于《光谱化学学报B原子光谱学》(Spectrochimica Acta Part BAtomic Spectroscopy)第106期45-66页发表的题为《作为燃烧诊断技术使用的可调谐二极管激光光谱技术》(Tunable diode laser spectroscopy as a technique forcombustion diagnostics)一文中,DAS方法存在诸如光散射引起的基线偏移、测试室的热结构的发射、光学标准具效应等问题。对于高温环境的测量过程,尤其是高温条件的动态过程测量中,红外波段的辐射能量(密度)随温度升高而增加,在温度达到1000K时,由于热辐射引起的背景噪声无法作为微小噪声被忽略。
在WMS方法中,使用高频正弦信号叠加在慢速斜坡扫描波形上,调制激光的波长和强度。如Gregory B.Rieker和Jay B.Jeffries在2009年于《应用光学》(Applied Optics)第48卷第29期5546页发表的题为《在恶劣环境中测量气体温度和浓度的免校准波长调制法吸收光谱技术》(Calibration-free wavelength-modulation spectroscopy formeasurements of gas temperature and concentration in harsh environments)一文中,对比于直接吸收方法,在超音速气流中二氧化碳测量中,WMS-2f方法在信噪比方面提升了4倍。如Alexander Klein和Oliver Witzel在2014年于《传感器》(Sensors)第14卷第11期21497-21513页发表的题为《快速、分时复用的直接吸收法和波长调制法激光吸收光谱技术》(Rapid,time-division multiplexed,direct absorption-and wavelengthmodulation-spectroscopy)一文中,表明在获取实验测量后,WMS的最终结果需要通过相当复杂的计算建模和拟合算法获得。如Guojie Tu和Fengzhong Dong在2015年于《IEEE传感器期刊》(IEEE Sensors Journal)第15卷第6期3535-3542页发表的题为《常压条件下可调谐二极管激光吸收光谱***的随机噪声和长期漂移分析》(Analysis of Random Noise andLong-Term Drift for Tunable Diode Laser Absorption Spectroscopy System atAtmospheric Pressure)一文中,指出在严酷的工业环境中,压力和温度的不同给绝对吸收光谱半宽带来了不确定因素,从而降低了长期传感器部署的适用性。在早期的TDLAS实现中,通常采用外部光学斩波器或铌酸锂振幅调制器。如James R.P.Bain和Michael Lengden在2015年于《IEEE传感器期刊》(IEEE sensors journal)第16卷第3期1-1页发表的题为《基于可调谐二极管激光技术采用平衡检测外调幅法恢复绝对吸收线形状的方法》(Recoveryof Absolute Absorption Line Shapes in Tunable Diode Laser Spectroscopy usingExternal Amplitude Modulation with Balanced Detection)一文中,指出典型的光学斩波器的调制频率一般可达到kHz,铌酸锂振幅调制器的调制频率可达到MHz,但需使激光达到偏振状态而损失一部分光强,且需要的光路结构与传统的直接吸收法相比更复杂。
如Y Deguchi和M Noda在2002年于《细菌学期刊》(Journal of Bacteriology)第190卷第7期2637-2641页发表的题为《工业应用的温度和物种浓度监测使用的激光诊断技术》(Industrial applications of temperature and species concentrationmonitoring using laser diagnostics)一文中,干涉仪无论在DAS方法还是WMS方法中都作为标准具来确保光学频率的同步。
基于以上背景,本文发明了一种基于干涉调制的激光吸收光谱测温方法与***,通过单一光路双谱线对均匀温度进行测量,通过一条光路同时获取标准具数据和吸收谱信息。在相对简化***结构的基础上,实现了光强的调制,一定程度上对背景噪声有抑制作用,并实现了温度的测量。
(三)发明内容
基于激光吸收光谱技术的气体测温仪的核心技术是如何从激光信号中准确地提取出吸收光谱信息。一般的基于激光吸收光谱技术的气体测温***,均需要单独一条光路作为标准具使用,以提取相对波数随时间的变化关系,光学***较为复杂。
本发明提出一种基于干涉调制的吸收光谱测温***和基于迭代的光谱重建方法,马赫曾德光纤干涉仪作为光强调制器件的同时复用为标准具,不再需要单独一条光路作为标准具使用,简化了光学***并可抑制背景噪声。
所用元件包括:激光发生单元、马赫曾德光纤干涉仪、光电探测单元和数据处理单元等。
本发明采用的技术方案是:激光经过马赫曾德光纤干涉仪之后,通过准直镜变为平行的空间光,穿过被测气体后由光电探测器采集。通过对采集测量信号进行解调并拟合得到光谱数据和基线数据的参数初值,而后通过迭代过程对其进行修正,得到光谱拟合参数的最终值。通过光谱数据的拟合参数求得其积分吸收率,并解算得到被测气体的温度值。
通过采集数据分别解算输出激光相对波数数据和吸收光强数据过程中,采用的解算方法步骤如下:
步骤一,激光发生模块产生的激光信号首先经过马赫曾德干涉仪得到光强调制波长扫描的测量信号,而后穿过被测对象,最终被光电探测器采集,采集信号在原理上形式如式(1)所示,
Figure GDA0002406341760000021
其中,IMZI(t)是光电探测器采集到随时间变化的测量数据,I0(t)是由激光器产生的激光信号未经过马赫曾德干涉仪和被测对象前的光强数据,vFSR是干涉仪的自由光谱范围(FSR,Free spectral range),K1、K2、φ0是由马赫曾德光纤干涉仪自身结构决定的参数,当使用的干涉仪固定时这些参数均为常量,v(t)是输出激光相对波数数据,α(v(t))为被测对象的吸收谱;
步骤二,根据采集数据解算输出激光相对波数数据v(t),经过马赫曾德干涉仪和被测对象的激光信号如(1)式所示,激光相对波数数据可以通过光强数据解算得到,在激光信号强度两次极值点之间,波数变化为一个自由光谱范围;为了在有限的采样点数下更加准确提取极值点的位置,提取相对波数数据时采用如下方法,对IMZI(t)信号进行高通滤波,其截止频率为激光信号扫描频率;记录滤波后信号的过零点,每两个相邻的过零点的中点位置即为极值点位置;两个极值点所在数据点所代表的相对波数相差一个自由光谱范围,由此可得到数据激光的相对波数数据v(t);
步骤三,对采集数据进行解调,得到吸收光强数据:将激光的相对波数数据v(t)代入(1),获得相对波数v与测量数据的映射关系,如(2)式
Figure GDA0002406341760000031
根据马赫曾德干涉仪的自由光谱范围vFSR,可以给出如(3)式所述信号
Figure GDA0002406341760000032
将(2)(3)联立,进行正交解调后可以得到吸收光强数据I(t),如(4)式所示
Figure GDA0002406341760000033
其中I0(t)是由激光器产生的激光信号在进入马赫曾德干涉仪前的光强信号,Ibase(t)是与I0(t)成正比的信号,在数据处理过程中为理想情况下的基线数据,α(v(t))为被测对象的吸收谱。
在解算过程通过拟合的方式得到基线参数和吸收光谱参数,并通过迭代修正拟合结果,进而解算出被测气体的温度值,采用的解算方法包括如下步骤:
步骤一,确定在不同过程中拟合所使用的数学模型,根据吸收光强数据,通过最小二乘拟合的方式估算基线的拟合参数,采用三次多项式对基线进行拟合,得到基线的拟合参数组ParaBase0,使用的模型如(5)所示
Figure GDA0002406341760000041
其中,拟合参数组ParaBase0为拟合的多项式参数。而后,基于比尔朗伯定律得到吸收率的数据α(t),如(6)所示
Figure GDA0002406341760000042
代入激光的相对波数数据v(t),得到吸收谱数据α(v),并在此基础上进行拟合,采用的模型如(7)所示,
Figure GDA0002406341760000043
其中,A为吸收谱的积分吸收率,φv(v)为Voigt线型展宽函数。由于Voigt线型函数计算包含卷积分的计算,因此无法直接得到其解析解,只能通过数值计算来进行求解,为了快速得到谱线中心处的吸收系数和Voigt线型半高宽,采用了如(8)所示简化模型,
Figure GDA0002406341760000044
其中,v0为吸收谱中心波数,ΔvC为碰撞展宽谱线半宽,ΔvD为多普勒展宽谱线半宽,ΔvV为Voigt线型函数谱线半宽,可由经验公式(9)求得,φv(v0)为Voigt线型函数在谱线中心频率处取值,可由(10)、(11)和(12)求得。
Figure GDA0002406341760000045
Figure GDA0002406341760000046
Figure GDA0002406341760000047
Figure GDA0002406341760000051
因此,以简化的吸收谱模型(7)对吸收谱数据进行拟合时,吸收谱参数ParaSpec0需要包括A、ΔvC、ΔvD和v0四个参数。
将α(v)作为原始数据,采用吸收谱拟合模型(7)进行拟合,得到吸收光谱的拟合参数ParaSpec0
步骤二,将获得的基线和吸收光谱的拟合参数作为初始值,根据吸收光强数据I(t)与使用基线的拟合参数ParaBasei-1计算得到的基线数据Ibasefitting(i-1)(t)相除并根据(6)求得吸收光谱数据αi(t);
步骤三,将吸收光谱数据αi(t)和激光的相对波数数据v(t),得到吸收光谱与激光的相对波数之间的映射αi(v),吸收谱模型(7)对吸收谱数据进行拟合,得到吸收光谱参数ParaSpeci
步骤四,以吸收光强数据为拟合目标,根据吸收光谱参数ParaSpeci和相对波数对吸收光强数据的映射v(t),对基线进行拟合得到参数ParaBasei,和拟合的吸收光强数据Ii(t);
步骤五,判断吸收光强数据I(t)与其拟合数据Ii(t)的残差平方Si和与上一次残差平方和的相对变化,当变化量大于0.1%时,返回步骤二,否则,输出修正后吸收光谱的拟合参数ParaSpeci和基线的拟合参数ParaBasei
步骤六,根据吸收光谱参数得到积分吸收率,通过多条吸收谱的积分吸收率进行解算得到被测气体温度。
本发明的优点在于:本发明中通过复用马赫曾德光纤干涉仪作为标准具和光强调制器,测量过程只需要一条光路,化简了基于可调谐吸收光谱技术采用直接吸收法测温的***结构;本发明中通过复用马赫曾德光纤干涉仪作为标准具和光强调制器,在直接吸收法测温的基础上进行了光强调制,在一定程度上对背景噪声有抑制作用;在解算算法中,以采用基线拟合方法得到的光谱参数和基线参数作为初值,通过迭代的方式对拟合结果进行修正,得到更加准确的光谱参数,减小了***误差。
(四)附图说明
图1是单一路径激光吸收光谱测温***结构图;
图2是是解调后吸收信号随时间的变化;
图3是相对波数随时间的关系;
图4仿真数据和对比迭代前后结果的关系。
(五)具体实施方式
在本实施例中,给定一个均匀温度的被测对象气体,通过数值***明了该方法的有效性。
下面结合附图对本发明作进一步的说明:
在图1中,可调谐半导体激光器控制模块(1)、可调谐半导体激光器(2)共同组成激光发生单元,通过光纤与马赫曾德光纤干涉仪(3)相连,通过准直镜(4)将光纤内的传输光转变成准直光(平行光),穿过被测气体(5)后,通过光电探测器(6),由数据采集模块(7)对测量数据进行采集和上传,通过数据处理单元(8)进行后续数据处理得到最终的解算结果。
仿真过程中,仿真的环境设置如下,吸收目标为水分子,摩尔浓度0.5%,温度298K,压强1atm,吸收距离80cm。使用的马赫曾德光纤干涉仪的参数K1=K2=0.5,vFSR=0.01cm-1,选定的吸收谱的中心波数分别为7183.89cm-1、7185.59cm-1
步骤一,由激光发生模块发出激光,通过马赫曾德光纤干涉仪后通过准直镜变为空间光,穿过待测区域,后经探测器接收,得到一路吸收信号,信号形式如(1)所示,***结构如图1所示;
Figure GDA0002406341760000061
步骤二,根据采集数据重建相对波数与吸收光强数据间的映射关系,首先对IMZI(t)信号高通滤波,其截止频率为激光信号扫描频率;而后,记录其过零点,每两个相邻的过零点的中点位置即为极值点位置;两个极值点所在数据点所代表的相对波数相差一个自由光谱范围(FSR),据此,我们可以得到相对波数与吸收光强数据间的映射关系v(t),结果如图2所示。
步骤三,根据探测信号解算出光强吸收数据,首先,根据相对波数与吸收光强数据间的映射关系v(t),进行坐标变换,得到如(2)所示
Figure GDA0002406341760000062
根据马赫曾德干涉仪的自由光谱范围vFSR,可以给出参考信号的形式如(3)式所述
Figure GDA0002406341760000063
进行解调后可以得到与传统的直接吸收法测量信号相同的吸收光强数据,如(4)式所示
I(t)=I0(t)·α(v(t)) (4)
其中I0(t)是由激光器产生的激光信号在进入马赫曾德干涉仪前的光强信号,α(v(t))为被测对象的吸收谱,结果如图3所示。
步骤四,通过迭代过程对基线参数和吸收光谱参数进行解算和修正。首先通过基线拟合的方式得到基线和吸收谱的拟合参数的估算值,而后以估算的基线和吸收光谱的拟合参数作为初始值。其中,中心波数为7183.89cm-1的光谱拟合参数初始值ParaSpec7183-0的ΔvC、ΔvD、A和v0分别是0.1860、0.1860、0.0113和0.0838,基线拟合参数初始值ParaBase7183-0=[0,0,0.0003,3.9044];中心波数为7185.59cm-1的光谱拟合参数初始值ParaSpec7185-0的ΔvC、ΔvD、A和v0分别是0.5162、0.5162、0.0973和0.1952,基线拟合参数初始值ParaBase7185-0=[0,0,0.0003,3.8880];
步骤五,开始迭代过程,根据吸收光强数据与使用基线拟合参数ParaBasei计算得到的基线数据相除并取负对数得到吸收光谱数据;之后,对吸收光谱数据进行Voigt函数拟合,得到吸收光谱参数ParaSpeci;其次,以吸收光强数据为拟合目标,根据吸收光谱参数ParaSpeci和相对波数对吸收光强数据的映射v(t),对基线进行拟合得到参数ParaBasei,和拟合的吸收光强数据Ii(t);最后,判断吸收光强数据I(t)与其拟合数据Ii(t)的残差平方Si和与上一次残差平方和的相对变化,当变化量大于0.1%时,返回继续迭代过程的第一步继续进行拟合,否则,认为此过程已收敛,输出修正后吸收光谱的拟合参数ParaSpeci和基线的拟合参数ParaBasei。其中,中心波数为7183.89cm-1的吸收谱拟合参数迭代经历了186次,光谱拟合参数ParaSpec7183-186的ΔvC、ΔvD、A和v0分别是0.0117、0.1055、0.0047和0.0841,基线拟合参数值ParaBase7183-186=[-1.0848e-14,2.1269e-10,0.0003,3.8971],中心波数为7185.59cm-1的吸收谱拟合参数迭代经历了290次,光谱拟合参数值ParaSpec7185-290的ΔvC、ΔvD、A和v0分别是0.0716、0.1928、0.0240和0.1980,基线拟合参数值ParaBase7185-290=[-6.10e-15,9.21e-11,2.95e-4,3.8837];
根据吸收光谱参数得到积分吸收率,通过多条吸收谱的积分吸收率进行解算得到被测气体温度,解算得到的吸收谱初始值与通过迭代修正后的吸收谱的对比如图4所示。设定温度为298K,解算温度为299.9K。
以上对本发明及其实施方式的描述,并不局限于此,附图中所示仅是本发明的实施方式之一。在不脱离本发明创造宗旨的情况下,不经创造地设计出与该技术方案类似的结构或实施例,均属本发明保护范围。

Claims (4)

1.基于干涉调制原理的激光吸收光谱测温方法,包括激光发生单元、马赫曾德光纤干涉仪、准直镜、光电探测器、数据采集模块和数据处理单元,其特征在于,激光信号由激光发生单元产生后经过马赫曾德光纤干涉仪完成强度调制,穿过被测对象后由光电探测器探测,由数据采集模块采集数据,并将其上传至数据处理单元,数据处理过程中,通过采集数据分别解算出吸收光强数据和相对波数与吸收光强数据之间的映射关系,而后通过拟合的方式根据吸收光强数据中近似无吸收的部分解算出基线参数的估计值,根据相对波数与吸收光强数据间的映射关系解算得吸收光谱参数的估计值,并通过迭代的方式对基线和吸收光谱参数的估计值进行修正,根据修正后的吸收光谱参数得到对应谱线在被测目标环境下的积分吸收率,结合多条谱线积分吸收率解算出被测气体的温度值。
2.根据权利要求1所述的基于干涉调制原理的激光吸收光谱测温方法,其特征在于:通过采集数据分别解算出吸收光强数据和相对波数与吸收光强数据之间的映射关系,包括以下步骤,
步骤一,激光发生单元产生的激光信号首先经过马赫曾德干涉仪得到光强调制波长扫描的测量信号,而后穿过被测对象,最终被光电探测器采集,采集信号在原理上形式如式(1)所示,
Figure FDA0002406341750000011
其中,IMZI(t)是光电探测器采集到随时间变化的测量数据,I0(t)是由激光发生单元产生的激光信号未经过马赫曾德干涉仪和被测对象前的光强数据,vFSR是干涉仪的自由光谱范围(FSR,Free spectral range),K1、K2、φ0是由马赫曾德光纤干涉仪自身结构决定的参数,当使用的干涉仪固定时这些参数均为常量,v(t)是输出激光相对波数数据,α(v(t))为被测对象的吸收谱;
步骤二,根据采集数据解算输出激光相对波数数据v(t),经过马赫曾德干涉仪和被测对象的激光信号如(1)式所示,激光相对波数数据可以通过光强数据解算得到,在激光信号强度两次极值点之间,波数变化为一个自由光谱范围;为了在有限的采样点数下更加准确提取极值点的位置,提取相对波数数据时采用如下方法,对IMZI(t)信号进行高通滤波,其截止频率为激光信号扫描频率;记录滤波后信号的过零点,每两个相邻的过零点的中点位置即为极值点位置;两个极值点所在数据点所代表的相对波数相差一个自由光谱范围,由此可得到数据激光的相对波数数据v(t);
步骤三,对采集数据进行解调,得到吸收光强数据:将激光的相对波数数据v(t)代入(1),获得相对波数v与测量数据的映射关系,如(2)式
Figure FDA0002406341750000012
式中I0(v)是随波数变化的未经过马赫曾德干涉仪和被测对象前的激光光强,是将相对波数数据v(t)代入激光初始光强数据I0(t)的结果;
根据马赫曾德干涉仪的自由光谱范围vFSR,可以给出如(3)式所述信号
Figure FDA0002406341750000021
将(2)(3)联立,进行正交解调后可以得到吸收光强数据I(t),如(4)式所示
Figure FDA0002406341750000022
其中I0(t)是由激光发生单元产生的激光信号在进入马赫曾德干涉仪前的光强信号,Ibase(t)是与I0(t)成正比的信号,在数据处理过程中为理想情况下的基线数据,α(v(t))为被测对象的吸收谱。
3.根据权利要求1所述的基于干涉调制原理的激光吸收光谱测温方法,其特征在于:解算过程通过拟合的方式得到基线参数和吸收光谱参数,并通过迭代修正拟合结果,进而解算出被测气体的温度值,解算过程包括如下步骤,
步骤一,根据吸收光强数据,通过最小二乘拟合的方式估算基线的拟合参数,采用三次多项式对基线进行拟合,得到基线的拟合参数组ParaBase0,使用的模型如(5)所示
Figure FDA0002406341750000023
其中,t为时间,拟合参数组ParaBase0为拟合的多项式参数,ParaBase0(j)为基线拟合参数组中的第j个元素,而后,基于比尔朗伯定律得到吸收率的数据α(t),如(6)所示
Figure FDA0002406341750000024
代入激光的相对波数数据v(t),得到吸收谱数据α(v),并在此基础上进行拟合,采用的模型如(7)所示,
Figure FDA0002406341750000025
其中,A为吸收谱的积分吸收率,φv(v)为Voigt线型函数,由于Voigt线型函数计算包含卷积分的计算,因此无法直接得到其解析解,只能通过数值计算来进行求解,为了快速得到谱线中心处的吸收系数和Voigt线型半高宽,采用了如(8)所示简化模型,
Figure FDA0002406341750000031
其中,v为吸收谱波数点,v0为吸收谱中心波数,ΔvC为碰撞展宽谱线半宽,ΔvD为多普勒展宽谱线半宽,ΔvV为Voigt线型函数谱线半宽,可由经验公式(9)求得,φv(v0)为Voigt线型函数在谱线中心频率处取值,可由(10)、(11)和(12)求得:
Figure FDA0002406341750000032
Figure FDA0002406341750000033
Figure FDA0002406341750000034
Figure FDA0002406341750000035
因此,以简化的吸收谱模型(7)对吸收谱数据进行拟合时,吸收谱参数ParaSpec0需要包括A、ΔvC、ΔvD和v0四个参数;
将α(v)作为原始数据,采用吸收谱拟合模型(7)进行拟合,得到吸收光谱的拟合参数ParaSpec0
步骤二,将获得的基线和吸收光谱的拟合参数作为初始值,根据吸收光强数据I(t)与第(i-1)次迭代中计算得到的基线数据Ibasefitting(i-1)(t)相除,并根据(6)求得第i次迭代中的吸收光谱数据αi(t),其中第(i-1)次迭代中的基线数据Ibasefitting(i-1)(t)是使用第(i-1)次迭代得到的基线的拟合参数ParaBasei-1根据式(5)计算得到的;
步骤三,将吸收光谱数据αi(t)和激光的相对波数数据v(t),得到吸收光谱与激光的相对波数之间的映射αi(v),吸收谱模型(7)对吸收谱数据进行拟合,得到第i次迭代中拟合得到的吸收光谱参数ParaSpeci
步骤四,以吸收光强数据为拟合目标,根据第i次迭代中拟合得到的吸收光谱参数ParaSpeci和相对波数对吸收光强数据的映射v(t),对基线进行拟合得到第i次迭代拟合得到的基线拟合参数ParaBasei,和第i次迭代拟合的吸收光强数据Ii(t);
步骤五,判断吸收光强数据I(t)与其第i次迭代拟合数据Ii(t)的残差平方Si和与上一次残差平方和的相对变化,当变化量大于0.1%时,返回步骤二,否则,输出修正后吸收光谱的拟合参数ParaSpeci和基线的拟合参数ParaBasei
步骤六,根据吸收光谱参数得到积分吸收率,通过多条吸收谱的积分吸收率进行解算得到被测气体温度。
4.基于干涉调制原理的激光吸收光谱测温***,包括激光发生单元、马赫曾德光纤干涉仪、准直镜、光电探测器、数据采集模块和数据处理单元,其特征在于,激光信号由激光发生单元产生后经过马赫曾德光纤干涉仪完成强度调制,穿过被测对象后由光电探测器探测,由数据采集模块采集数据,并将其上传至数据处理单元,数据处理过程中,通过采集数据分别解算出吸收光强数据和相对波数与吸收光强数据之间的映射关系,而后通过拟合的方式根据吸收光强数据中近似无吸收的部分解算出基线参数的估计值,根据相对波数与吸收光强数据间的映射关系解算得吸收光谱参数的估计值,并通过迭代的方式对基线和吸收光谱参数的估计值进行修正,根据修正后的吸收光谱参数得到对应谱线在被测目标环境下的积分吸收率,结合多条谱线积分吸收率解算出被测气体的温度值。
CN201811235909.XA 2018-10-24 2018-10-24 一种基于干涉调制原理的激光吸收光谱测温方法和*** Active CN108981953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811235909.XA CN108981953B (zh) 2018-10-24 2018-10-24 一种基于干涉调制原理的激光吸收光谱测温方法和***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811235909.XA CN108981953B (zh) 2018-10-24 2018-10-24 一种基于干涉调制原理的激光吸收光谱测温方法和***

Publications (2)

Publication Number Publication Date
CN108981953A CN108981953A (zh) 2018-12-11
CN108981953B true CN108981953B (zh) 2020-07-31

Family

ID=64543704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811235909.XA Active CN108981953B (zh) 2018-10-24 2018-10-24 一种基于干涉调制原理的激光吸收光谱测温方法和***

Country Status (1)

Country Link
CN (1) CN108981953B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110031198B (zh) * 2019-04-25 2020-11-10 上海工程技术大学 一种考虑能量的应变载荷谱编制方法
CN111239072B (zh) * 2019-11-25 2023-05-12 中国航空工业集团公司北京长城计量测试技术研究所 一种精确测量燃烧气体温度的方法
CN113252204B (zh) * 2021-04-01 2022-04-01 北京航空航天大学 一种基于幅度调制的激光吸收光谱温度测量方法
CN112798550B (zh) * 2021-04-14 2021-07-13 四川大学 一种宽测量范围的激光吸收光谱燃烧诊断方法
CN113390826B (zh) * 2021-05-28 2022-11-01 中国科学院力学研究所 一种用于视线非均匀吸收光谱测量的积分吸光度测量误差估计方法
CN113176222B (zh) * 2021-06-09 2022-07-29 复旦大学 一种基于直接吸收谱的气体浓度反演方法
CN113358160B (zh) * 2021-07-08 2022-05-17 天津大学 一种大气数据测量方法及***
CN114993987B (zh) * 2022-05-16 2024-04-30 北京航空航天大学 基于吸收光谱幅值调制的温度及气体浓度测量方法与***
CN117747008B (zh) * 2024-02-07 2024-05-03 山东省科学院激光研究所 一种气体激光吸收光谱的基线拟合及降噪方法、***

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102680020A (zh) * 2012-05-16 2012-09-19 清华大学 一种基于波长调制光谱技术的气体参数在线测量方法
CN103557886A (zh) * 2013-10-29 2014-02-05 东南大学 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时测量的方法
JP2014085301A (ja) * 2012-10-26 2014-05-12 Nippon Telegr & Teleph Corp <Ntt> 波長掃引光源用光強度計測装置
CN107219190A (zh) * 2017-05-05 2017-09-29 东南大学 一种基于三角波调制的免标定气体参数测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9766483B2 (en) * 2016-01-21 2017-09-19 Sumitomo Electric Industries, Ltd. Optical transceiver implementing erbium doped fiber amplifier

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102680020A (zh) * 2012-05-16 2012-09-19 清华大学 一种基于波长调制光谱技术的气体参数在线测量方法
JP2014085301A (ja) * 2012-10-26 2014-05-12 Nippon Telegr & Teleph Corp <Ntt> 波長掃引光源用光強度計測装置
CN103557886A (zh) * 2013-10-29 2014-02-05 东南大学 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时测量的方法
CN107219190A (zh) * 2017-05-05 2017-09-29 东南大学 一种基于三角波调制的免标定气体参数测量方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Rapid, Time-Division Multiplexed, Direct Absorption- and Wavelength Modulation-Spectroscopy;Alexander Klein 等;《Sensors》;20141114;第14卷(第11期);21497-21513 *
二极管激光吸收传感器测量超声速流场的温度和速度;王广宇 等;《光学学报》;20130930;第33卷(第9期);0912009-1至0912009-7 *
基于波长调制谐波信号主峰拟合的气体浓度测量方法;杨阳 等;《北京航空航天大学学报》;20171130;第43卷(第11期);2187-2192 *
多谱线 TDLAS 测温中谱线选择的实验研究;陈亚婧 等;《电子测量与仪器学报》;20171231;第31卷(第12期);1948-1952 *

Also Published As

Publication number Publication date
CN108981953A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
CN108981953B (zh) 一种基于干涉调制原理的激光吸收光谱测温方法和***
CN111077109B (zh) 一种基于双光梳光谱技术的温度和浓度测量***和方法
Huang et al. Frequency-division multiplexing and main peak scanning WMS method for TDLAS tomography in flame monitoring
CN111077110B (zh) 一种基于双光梳光谱的温度场和浓度场测量***和方法
CN109100325B (zh) 一种基于光谱吸收率二次谐波特征提取的气体浓度测量方法
JP5983779B2 (ja) ガス吸収分光装置及びガス吸収分光方法
CN108760681A (zh) 一种基于波形分解的路径平均温度测量***与方法
CN109991189B (zh) 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法
CN113252204B (zh) 一种基于幅度调制的激光吸收光谱温度测量方法
CN111829982A (zh) 多组分气体浓度检测及温度误差修正方法
CN112378616B (zh) 基于波长调制吸收光谱的高速流场多参数测量***及方法
CN109696415B (zh) 一种基于快速傅里叶变换的气体吸收率在线测量方法
CN109100044B (zh) 基于单光路多光谱的气体温度概率密度分布拟合重建方法
CN108562237B (zh) 一种采用hcn气室在光频域反射传感***中进行光谱校准的装置和方法
CN110514622B (zh) 一种基于直方图信息的二维温度和浓度重建***与方法
CN114034653B (zh) 基于深度学习的波长调制吸收光谱层析重建***
CN114993987B (zh) 基于吸收光谱幅值调制的温度及气体浓度测量方法与***
CN113447458A (zh) 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法
CN111351770A (zh) 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法
CN109270027A (zh) 一种基于正弦调制时域拟合的气体吸收率在线测量方法
CN111089661A (zh) 一种基于激光吸收光谱的温度快速提取方法
Xu et al. FPGA-based real-time implementation of temperature measurement via tunable diode laser absorption spectroscopy
CN114544547A (zh) 一种用于波长调制光谱信号在线递推解调的方法与***
CN103411923B (zh) 一种采用双路可调谐激光吸收光谱技术的归一化降噪方法
CN108333143B (zh) 一种基于可调谐激光吸收光谱的水汽浓度测量修正方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant