CN107255627A - 一种基于级数展开的气体浓度测量方法及其检测装置 - Google Patents
一种基于级数展开的气体浓度测量方法及其检测装置 Download PDFInfo
- Publication number
- CN107255627A CN107255627A CN201710339154.7A CN201710339154A CN107255627A CN 107255627 A CN107255627 A CN 107255627A CN 201710339154 A CN201710339154 A CN 201710339154A CN 107255627 A CN107255627 A CN 107255627A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- phi
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
- G01N2021/396—Type of laser source
- G01N2021/399—Diode laser
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Optics & Photonics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于级数展开的气体浓度测量方法及其检测装置，根据高分辨率气体分子吸收光谱数据库、二极管激光器的光强调制系数和频率调制系数等参数快速拟合出某一浓度气体的归一化各次谐波信号，并与实测信号的扣除背景后归一化各次谐波信号进行最小方差拟合。当拟合方差小于设定阈值，则可认为该拟合浓度即为实测气体浓度，从而实现对气体浓度的测量。本发明可大大降低计算量；且利于嵌入式系统实现，促进仪器的小型化，数字化。
Description
技术领域
本发明属于气体传感和检测技术领域，具体涉及一种基于级数展开的气体浓度测量方法，属于激光吸收光谱领域，该方法可适用于一般的个人计算机和嵌入式系统中。
背景技术
可调谐半导体激光吸收光谱(Tunable Diode Laser Absorption Spectroscopy，TDLAS)技术具有高灵敏度、高选择性和快速响应等特点，可以实现对气体的高灵敏度、高分辨率及快速测量。由于它上述优点，TDLAS技术已经广泛应用于大气环境监测、燃烧诊断、等离子反应监测过程产物、烟道气体检测等。通常，TDLAS技术结合波长调制技术和长光程池吸收，以提高探测精度与灵敏度。波长调制技术是指分布反馈式二极管激光器(Distributed feedback laser，DFB)在低频扫描信号调制的基础上叠加一个高频正弦调制频率，通过在高频调制频率及高频调制频率的高次谐波(如二次谐波处)提取谐波信号的幅值信息，该幅值信息与待测气体的浓度、温度有关，由此实现对气体浓度、温度等参数的测量。目前主流的浓度反演算法都是通过锁相放大器提取谐波分量并用一次谐波对二次谐波进行归一化后，根据归一化后的谐波信号峰值来确定气体的浓度参数。该方法的确定在于实际测量过程中谐波信号峰值点会存在漂移、随机误差等影响，故测量精度存在较大误差。大量研究人员试图改进目前现有的缺点，因此发明一种改善从归一化后的谐波信号反演气体参数的方法是非常必要的。
发明内容
发明目的：为了解决现有技术存在的由归一化后的谐波信号峰值来确定气体的浓度参数时存在的漂移、随机误差等问题，本发明提供一种基于级数展开的气体浓度测量方法。
本发明的另一目的是提供一种基于级数展开的气体浓度检测装置。
技术方案：一种基于级数展开的气体浓度测量方法，包括以下步骤：
(1)通过DFB二极管激光器温度控制模块设定DFB二极管激光器的工作温度；
(2)通过DFB二极管激光器电流驱动模块设定DFB二极管激光器工作的中心电流以确保DFB二极管激光器静态工作点的原始入射光信号频率在吸收谱线中心附近；
(3)计算机控制DAC模块产生初始信号，并将初始信号传递至DFB二极管激光器电流驱动模块，经转换后形成注入电流注入二极管激光器，DFB二极管激光器受注入电流的调制发射激光；
(4)通过标准具，经光电探测器和ADC模块获得DFB二极管激光器发光频率与时间的函数关系，记为υ(t)；
(5)在待测气体未通入气体吸收池的情况下，光电探测器采集原始入射光强信号；
(6)将待测气体通入气体吸收池中，光电探测器采集透射光强信号，并将透射光强信号输入计算机，计算机对函数υ(t)、原始入射光强信号及透射光强信号进行处理得到待测气体浓度。
一种基于级数展开的气体浓度测量方法所用的检测装置，包括计算机、DAC信号发生器、DFB二极管激光器电流驱动模块、DFB激光二极管、DFB二极管激光器温度控制模块、气体吸收池、光电探测器、带通滤波放大器、ADC采样模块、标准具，所述计算机用户输出调制波形给DAC信号发生器；DAC信号发生器用于产生调制电压信号；DFB二极管激光器电流驱动模块用于将调制电压信号转换成注入电流注入DFB激光二极管；DFB二极管激光器温度控制模块用于设定DFB二极管激光器的工作温度；气体吸收池用于存放待测气体；光电探测器用于采集待测气体的光强信号；带通滤波放大器用于光强信号调理，ADC采样模块用于对光强信号采样并将采样信号输入计算机。
有益效果：本发明提供的一种基于级数展开的气体浓度测量方法及其检测装置，在计算仿真归一化二次谐波R2f/1f时，通过读取预存在系统中的归一化线型函数φ(υ(t))及其平方g(υ(t))＝φ(υ(t))·φ(υ(t))各次谐波分量系数φn和gn可以直接将浓度代入计算即可，无需将浓度值代入比尔朗伯公式得到仿真信号再通过数字锁相获取各次谐波，可大大降低计算量；且利于嵌入式系统实现，促进仪器的小型化，数字化；可以实现拟合法反演气体参数，实现快速、准确测量，具有方便、简单的优点。
附图说明
图1是一种基于级数展开的气体浓度测量方法的气体检测装置图；
图2(a)是吸收线型φ(υ(t))傅里叶级数展开直流分量系数φ0波形图；
图2(b)是吸收线型φ(υ(t))傅里叶级数展开一次谐波系数φ1波形图；
图2(c)是吸收线型φ(υ(t))傅里叶级数展开二次谐波系数φ2波形图；
图2(d)是吸收线型φ(υ(t))傅里叶级数展开三次谐波系数φ3波形图；
图3(a)是g(υ(t))＝φ(υ(t))·φ(υ(t))傅里叶级数展开直流分量系数g0波形图；
图3(b)是g(υ(t))＝φ(υ(t))·φ(υ(t))傅里叶级数展开直流分量系数g1波形图；
图3(c)是g(υ(t))＝φ(υ(t))·φ(υ(t))傅里叶级数展开直流分量系数g2波形图；
图3(d)是g(υ(t))＝φ(υ(t))·φ(υ(t))傅里叶级数展开直流分量系数g3波形图；
图4是根据傅里叶级数展开系数φn和gn逼近扣除背景后的归一化二次谐波信号的波形的拟合过程。
图1中有计算机1，DAC信号发生模块2，DFB二极管激光器电流驱动模块3，DFB激光二极管4，DFB二极管激光器温度控制模块5，气体吸收池6，光电探测器7，带通滤波放大电路8，ADC采样模块9，标准具10。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
基于级数展开的气体浓度测量方法及其检测装置，根据高分辨率气体分子吸收光谱(High Resolution Transmission，HITRAN)数据库、二极管激光器的光强调制系数和频率调制系数等参数快速拟合出某一浓度气体的归一化各次谐波信号，并与实测信号的扣除背景后归一化各次谐波信号进行最小方差拟合。当拟合方差小于设定阈值，则可认为该拟合浓度即为实测气体浓度，从而实现对气体浓度的测量。本发明可大大降低计算量；且利于嵌入式系统实现，促进仪器的小型化，数字化。
如图1所示，一种基于级数展开的气体浓度测量方法所用的检测装置，包括计算机、DAC信号发生器、DFB二极管激光器电流驱动模块、DFB激光二极管、DFB二极管激光器温度控制模块、气体吸收池、光电探测器、带通滤波放大器、ADC采样模块、标准具，所述计算机用户输出调制波形给DAC信号发生器；DAC信号发生器用于产生调制电压信号；DFB二极管激光器电流驱动模块用于将调制电压信号转换成注入电流注入DFB激光二极管；DFB二极管激光器温度控制模块用于设定DFB二极管激光器的工作温度；气体吸收池用于存放待测气体；光电探测器用于采集待测气体的光强信号；带通滤波放大器用于光强信号调理，ADC采样模块用于对光强信号采样并将采样信号输入计算机。
一种基于级数展开的气体浓度测量方法，包括以下步骤：
(1)通过DFB二极管激光器温度控制模块设定DFB二极管激光器的工作温度；
(2)通过DFB二极管激光器电流驱动模块设定DFB二极管激光器工作的中心电流以确保DFB二极管激光器静态工作点的原始入射光信号频率在吸收谱线中心附近；
(3)计算机控制DAC模块产生初始信号，并将初始信号传递至DFB二极管激光器电流驱动模块，经转换后形成注入电流注入二极管激光器，DFB二极管激光器受注入电流的调制发射激光；
所述步骤(3)中设所述初始信号为u＝us cos(ωst)+um cos(ωmt)，所述ωs、ωm、us、um分别为由计算机设定的扫描频率、调制频率、扫描电压系数、调制电压系数；转换后形成的注入电流为i(t)＝iavg+cos(2πfmt)，所述iavg为瞬时平均电流，可以为低频三角波、正弦波、锯齿波等，fm为调制频率。
(4)通过标准具，经光电探测器和ADC模块获得DFB二极管激光器发光频率与时间的函数关系，记为υ(t)；υ(t)＝υavg+kvcos(2πfmt)，所述υavg为瞬时平均电压，kv为光强调制系数。
(5)在待测气体未通入气体吸收池的情况下，光电探测器采集原始入射光强信号；
(6)将待测气体通入气体吸收池中，光电探测器采集透射光强信号，并将透射光强信号输入计算机，计算机对函数υ(t)、原始入射光强信号及透射光强信号进行处理得到待测气体浓度。
根据高分辨率气体分子吸收光谱(High Resolution Transmission，HITRAN)数据库、二极管激光器的光强调制系数和频率调制系数等参数快速拟合出某一浓度气体的归一化各次谐波信号，并与实测信号的扣除背景后归一化各次谐波信号进行最小方差拟合。当拟合方差小于设定阈值，则可认为该拟合浓度即为实测气体浓度，从而实现对气体浓度的测量，根据傅里叶级数展开系数φn和gn逼近扣除背景后的归一化二次谐波信号的波形的拟合过程见图4。具体处理过程包括以下步骤：
(61)通过对υ(t)进行级数展开与傅里叶分解推导从而得到归一化线型函数φ(υ(t))各次谐波分量系数φn以及φ(υ(t))的平方g(υ(t))＝φ(υ(t))·φ(υ(t))各次谐波分量系数gn；展开系数的波形图如图2和3所示。
(62)通过对原始入射光强拟合，确定原始入射光强信号以及频率随时间变化的傅里叶级数展开；
(63)采集透射光强信号，获取透射光强信号和扣除背景后的归一化二次谐波S2f/1f；对采集的透射光强信号进行平均滤波，经过数字锁相分别提取一次谐波分量和二次谐波分量，并进行扣除背景后归一化得到S2f/1f信号。
步骤(63)中将υ(t)代入到归一化线型函数(如Lorentz、Gauss或Voigt函数)中得到φ(υ(t))；由于归一化线型函数为偶函数，故将φ(υ(t))展开成傅里叶级数：
式中分别为n次谐波的分量系数，分别通过相敏解调求解：
φn＝(2-δn0)conv(φ(υ(t))*cos(2πnfmt),lps_coefficient)
式中conv表示卷积运算；lps_coefficient为低通滤波器；n＝0时，δn0为1；n>0时，δn0为1；
由于φ(υ(t))为fm周期函数，设g(υ(t))＝φ(υ(t))·φ(υ(t))也为fm周期函数，将g(υ(t))展开成傅里叶级数如下：
其中展开系数gn为：
式
(64)设原始入射光强信号为I0(t)＝Iavg+kIcos(2πfmt)，则由比尔-朗伯定律得到透射光强信号为：
It(t)＝I0exp(-PXLS(T)φ(υ))
式中P为气体总压，X为待测气体浓度，L为有效光程长，S(T)为线强，所述线强可通过查找HITRAN数据库获取；设A＝PLS(T)；可以将It展开成幂级数：
当PXLS(T)φ(υ)小于0.1时，
综上所述，
又因为I0(t)＝Iavg+kI cos(2πfmt)，故It＝I0exp(-PXLS(T)φ(υ))各次谐波系数分别为
式中除浓度X为未知数，其他参量均可在测量之前测得，归一化二次谐波信号R2f/1f为：
设定初始迭代浓度参数为X0，X0与φn和gn运算即可得到仿真归一化二次谐波R2f/1f；设定一个阈值，采用交替迭代法改变迭代值Xn，使得S2f/1f与R2f/1f的方差小于所述阈值，此时的Xn即为待测气体浓度。
Claims (8)
1.一种基于级数展开的气体浓度测量方法，其特征在于，包括以下步骤：
(1)通过DFB二极管激光器温度控制模块设定DFB二极管激光器的工作温度；
(2)通过DFB二极管激光器电流驱动模块设定DFB二极管激光器工作的中心电流以确保DFB二极管激光器静态工作点的原始入射光信号频率在吸收谱线中心附近；
(3)计算机控制DAC模块产生初始信号，并将初始信号传递至DFB二极管激光器电流驱动模块，经转换后形成注入电流注入二极管激光器，DFB二极管激光器受注入电流的调制发射激光；
(4)通过标准具，经光电探测器和ADC模块获得DFB二极管激光器发光频率与时间的函数关系，记为υ(t)；
(5)在待测气体未通入气体吸收池的情况下，光电探测器采集原始入射光强信号；
(6)将待测气体通入气体吸收池中，光电探测器采集透射光强信号，并将透射光强信号输入计算机，计算机对函数υ(t)、原始入射光强信号及透射光强信号进行处理得到待测气体浓度。
2.根据权利要求1所述的基于级数展开的气体浓度测量方法，其特征在于，所述步骤(3)中设所述初始信号为u＝uscos(ωst)+umcos(ωmt)，所述ωs、ωm、us、um分别为由计算机设定的扫描频率、调制频率、扫描电压系数、调制电压系数；转换后形成的注入电流为i(t)＝iavg+cos(2πfmt)，所述iavg为瞬时平均电流，fm为调制频率。
3.根据权利要求2所述的基于级数展开的气体浓度测量方法，其特征在于，所述步骤(4)中，υ(t)＝υavg+kvcos(2πfmt)，所述υavg为瞬时平均电压，kv为光强调制系数。
4.根据权利要求3所述的基于级数展开的气体浓度测量方法，其特征在于，所述步骤(6)计算机对信号进行处理的步骤包括：
(61)通过对υ(t)进行级数展开与傅里叶分解推导从而得到归一化线型函数φ(υ(t))各次谐波分量系数φn以及φ(υ(t))的平方g(υ(t))＝φ(υ(t))·φ(υ(t))各次谐波分量系数gn；
(62)通过对原始入射光强拟合，确定原始入射光强信号以及频率随时间变化的傅里叶级数展开；
(63)采集透射光强信号，获取透射光强信号和扣除背景后的归一化二次谐波S2f/1f；
(64)设定初始迭代浓度参数为X0，X0与φn和gn运算即可得到仿真归一化二次谐波R2f/1f；设定一个阈值，采用交替迭代法改变迭代值Xn，使得S2f/1f与R2f/1f的方差小于所述阈值，此时的Xn即为待测气体浓度。
5.根据权利要求4所述的基于级数展开的气体浓度测量方法，其特征在于，步骤(63)中将υ(t)代入到归一化线型函数中得到φ(υ(t))；由于归一化线型函数为偶函数，故将φ(υ(t))展开成傅里叶级数：
<mrow>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>&phi;</mi>
<mi>n</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;nf</mi>
<mi>m</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中分别为n次谐波的分量系数，分别通过相敏解调求解：
φn＝(2-δn0)conv(φ(υ(t))*cos(2πnfmt),lps_coefficient)
式中conv表示卷积运算；lps_coefficient为低通滤波器；n＝0时，δn0为1；n>0时，δn0为1；
由于φ(υ(t))为fm周期函数，设g(υ(t))＝φ(υ(t))·φ(υ(t))也为fm周期函数，将g(υ(t))展开成傅里叶级数如下：
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>g</mi>
<mi>n</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>&omega;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中展开系数gn为：
式中
6.根据权利要求4所述的基于级数展开的气体浓度测量方法，其特征在于，步骤(64)中归一化二次谐波信号R2f/1f的获取包括：
设原始入射光强信号为I0(t)＝Iavg+kIcos(2πfmt)，则由比尔-朗伯定律得到透射光强信号为：
It(t)＝I0exp(-PXLS(T)φ(υ))
式中P为气体总压，X为待测气体浓度，L为有效光程长，S(T)为线强，所述线强可通过查找HITRAN数据库获取；设A＝PLS(T)；可以将It展开成幂级数：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>t</mi>
<mo>=</mo>
<mi>I</mi>
<mn>0</mn>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>P</mi>
<mi>X</mi>
<mi>L</mi>
<mi>S</mi>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mi>&phi;</mi>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>I</mi>
<mn>0</mn>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>X</mi>
<mi>A</mi>
<mi>&phi;</mi>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>I</mi>
<mn>0</mn>
<mo>&CenterDot;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mrow>
<mi>n</mi>
<mo>!</mo>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mi>A</mi>
<mi>&phi;</mi>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
当PXLS(T)φ(υ)小于0.1时，
综上所述，
又因为
I0(t)＝Iavg+kI cos(2πfmt)，故It＝I0exp(-PXLS(T)φ(υ))各次谐波系数分别为
<mrow>
<msub>
<mi>G</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>></mo>
<mn>1</mn>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0.</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中除浓度X为未知数，其他参量均可在测量之前测得，归一化二次谐波信号R2f/1f为：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mn>2</mn>
<mi>f</mi>
<mo>/</mo>
<mn>1</mn>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>AX&phi;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<msub>
<mi>g</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
7.根据权利要求4所述的基于级数展开的气体浓度测量方法，其特征在于，步骤(63)中扣除背景后的归一化二次谐波S2f/1f的获取包括：对采集的透射光强信号进行平均滤波，经过数字锁相分别提取一次谐波分量和二次谐波分量，并进行扣除背景后归一化得到S2f/1f信号。
8.一种基于级数展开的气体浓度测量方法所用的检测装置，其特征在于：包括计算机、DAC信号发生器、DFB二极管激光器电流驱动模块、DFB激光二极管、DFB二极管激光器温度控制模块、气体吸收池、光电探测器、带通滤波放大器、ADC采样模块、标准具，所述计算机用户输出调制波形给DAC信号发生器；DAC信号发生器用于产生调制电压信号；DFB二极管激光器电流驱动模块用于将调制电压信号转换成注入电流注入DFB激光二极管；DFB二极管激光器温度控制模块用于设定DFB二极管激光器的工作温度；气体吸收池用于存放待测气体；光电探测器用于采集待测气体的光强信号；带通滤波放大器用于光强信号调理，ADC采样模块用于对光强信号采样并将采样信号输入计算机。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710339154.7A CN107255627B (zh) | 2017-05-15 | 2017-05-15 | 一种基于级数展开的气体浓度测量方法及其检测装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710339154.7A CN107255627B (zh) | 2017-05-15 | 2017-05-15 | 一种基于级数展开的气体浓度测量方法及其检测装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107255627A true CN107255627A (zh) | 2017-10-17 |
| CN107255627B CN107255627B (zh) | 2020-03-31 |
Family
ID=60027912
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710339154.7A Active CN107255627B (zh) | 2017-05-15 | 2017-05-15 | 一种基于级数展开的气体浓度测量方法及其检测装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107255627B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110779891A (zh) * | 2019-11-05 | 2020-02-11 | 北华航天工业学院 | 一种可调谐半导体激光监测二氧化碳浓度的方法 |
| CN115436320A (zh) * | 2022-09-05 | 2022-12-06 | 湖南五凌电力科技有限公司 | 一种混合气体背景下气体吸收线的修正方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101063655A (zh) * | 2006-04-27 | 2007-10-31 | Ir微系统股份有限公司 | 气体检测方法和气体检测装置 |
| CN101762557A (zh) * | 2009-10-28 | 2010-06-30 | 天津大学 | 基于有源内腔吸收光谱的气体种类识别方法及识别装置 |
| CN101936896A (zh) * | 2010-07-30 | 2011-01-05 | 合肥联合立体交通科技有限公司 | 一种乙醇气体浓度激光遥测装置 |
| CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
| CN103439291A (zh) * | 2013-08-09 | 2013-12-11 | 国家安全生产监督管理总局通信信息中心 | 一种基于tdlas的分布式光纤甲烷监测系统 |
| CN104749396A (zh) * | 2015-03-25 | 2015-07-01 | 中国科学院安徽光学精密机械研究所 | 一种基于波长调制-tdlas的高速气流流速检测系统 |
| CN105372205A (zh) * | 2015-11-25 | 2016-03-02 | 山西大学 | 基于S2f方法的免定标波长调制光谱气体检测方法 |
| CN105866068A (zh) * | 2016-03-28 | 2016-08-17 | 北京航空航天大学 | 一种基于递推最小二乘法的激光吸收光谱谐波信号峰值点提取方法 |
-
2017
- 2017-05-15 CN CN201710339154.7A patent/CN107255627B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101063655A (zh) * | 2006-04-27 | 2007-10-31 | Ir微系统股份有限公司 | 气体检测方法和气体检测装置 |
| CN101762557A (zh) * | 2009-10-28 | 2010-06-30 | 天津大学 | 基于有源内腔吸收光谱的气体种类识别方法及识别装置 |
| CN101936896A (zh) * | 2010-07-30 | 2011-01-05 | 合肥联合立体交通科技有限公司 | 一种乙醇气体浓度激光遥测装置 |
| CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
| CN103439291A (zh) * | 2013-08-09 | 2013-12-11 | 国家安全生产监督管理总局通信信息中心 | 一种基于tdlas的分布式光纤甲烷监测系统 |
| CN104749396A (zh) * | 2015-03-25 | 2015-07-01 | 中国科学院安徽光学精密机械研究所 | 一种基于波长调制-tdlas的高速气流流速检测系统 |
| CN105372205A (zh) * | 2015-11-25 | 2016-03-02 | 山西大学 | 基于S2f方法的免定标波长调制光谱气体检测方法 |
| CN105866068A (zh) * | 2016-03-28 | 2016-08-17 | 北京航空航天大学 | 一种基于递推最小二乘法的激光吸收光谱谐波信号峰值点提取方法 |
Non-Patent Citations (2)
| Title |
|---|
| 汪世美 等: "基于可调谐二极管激光吸收光谱遥测CH4浓度", 《光谱学与光谱分析》 * |
| 黄渐强 等: "用于气体检测的近红外半导体激光器温控系统", 《光子学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110779891A (zh) * | 2019-11-05 | 2020-02-11 | 北华航天工业学院 | 一种可调谐半导体激光监测二氧化碳浓度的方法 |
| CN110779891B (zh) * | 2019-11-05 | 2022-03-15 | 北华航天工业学院 | 一种可调谐半导体激光监测二氧化碳浓度的方法 |
| CN115436320A (zh) * | 2022-09-05 | 2022-12-06 | 湖南五凌电力科技有限公司 | 一种混合气体背景下气体吸收线的修正方法 |
| CN115436320B (zh) * | 2022-09-05 | 2024-03-15 | 湖南五凌电力科技有限公司 | 一种混合气体背景下气体吸收线的修正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107255627B (zh) | 2020-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102346137B (zh) | 气体浓度测量装置 | |
| CN106802288B (zh) | 基于可调谐激光和超连续谱激光的气体检测装置和方法 | |
| CN109991189B (zh) | 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法 | |
| CN105548075A (zh) | 一种玻璃药瓶内氧气含量的检测装置与方法 | |
| CN103616347B (zh) | 多气体谱线检测的最佳调制系数实现方法和装置 | |
| CN109696415B (zh) | 一种基于快速傅里叶变换的气体吸收率在线测量方法 | |
| CN103868885A (zh) | 基于复合多次谐波的气体浓度在线测量方法 | |
| CN103543124A (zh) | 一种基于软件锁相的可调激光吸收光谱气体检测方法 | |
| CN108801977B (zh) | 无标定痕量碳12和碳13二氧化碳气体探测装置及方法 | |
| CN105388120B (zh) | 基于wmrf模型的免定标波长调制光谱气体检测方法 | |
| CN106990071A (zh) | 一种全数字式气体检测仪 | |
| CN109270027A (zh) | 一种基于正弦调制时域拟合的气体吸收率在线测量方法 | |
| CN107037004A (zh) | 一种tdlas检测sf6电气设备中湿度的装置及方法 | |
| CN105158205B (zh) | 甲烷气团界面识别和可调量程激光遥测甲烷浓度的方法和装置 | |
| CN104568835A (zh) | 一种消除光纤传感系统中二次谐波剩余幅度调制的方法 | |
| CN109813639B (zh) | 一种基于红外光调制技术的颗粒物与气体浓度同步测量装置及其测量方法 | |
| CN110987870A (zh) | 基于波长调制光谱技术的实时监测气体浓度的系统和方法 | |
| CN109883991A (zh) | 碱金属原子蒸气同位素丰度的在线测量方法及装置 | |
| CN109724941A (zh) | 一种基于径向基神经网络的co2高温气体浓度检测方法 | |
| CN107255627A (zh) | 一种基于级数展开的气体浓度测量方法及其检测装置 | |
| CN204008440U (zh) | 用于差分气体检测的信号处理装置及其气体检测系统 | |
| CN110487748A (zh) | 一种基于直接激光吸收检测法的气体浓度的检测装置 | |
| CN105527247A (zh) | 一种基于正弦波调制的高灵敏度激光甲烷测量装置及方法 | |
| CN105866068A (zh) | 一种基于递推最小二乘法的激光吸收光谱谐波信号峰值点提取方法 | |
| CN107941467A (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 |