CN111398216A - 气体浓度检测的多频调制方法、气体浓度检测方法及系统 - Google Patents
气体浓度检测的多频调制方法、气体浓度检测方法及系统 Download PDFInfo
- Publication number
- CN111398216A CN111398216A CN202010430751.2A CN202010430751A CN111398216A CN 111398216 A CN111398216 A CN 111398216A CN 202010430751 A CN202010430751 A CN 202010430751A CN 111398216 A CN111398216 A CN 111398216A
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- modulation signal
- signal
- gas concentration
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000010521 absorption reaction Methods 0.000 claims abstract description 40
- 239000011521 glass Substances 0.000 claims description 10
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 43
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 7
- 239000001301 oxygen Substances 0.000 description 7
- 229910052760 oxygen Inorganic materials 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 5
- 230000033228 biological regulation Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 241001520299 Phascolarctos cinereus Species 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/01—Arrangements or apparatus for facilitating the optical investigation
-
- 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/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- 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)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种气体浓度检测的多频调制方法、气体浓度检测方法及系统,该多频调制方法包括:步骤S1:调节基础调制信号的幅值使得多频调制的气体浓度检测系统输出的二次谐波吸收曲线中峰值点横坐标与背景干涉条纹包络的极小值点横坐标相同;步骤S2:基于抖动调制信号的幅值调节区间、频率调节区间分别调节所述抖动调制信号的幅值和频率。其中,本发明在原单频调制系统中加入抖动调制信号构成多频调制系统,且通过调节基础调制信号的幅值以及抖动调制信号的幅值和频率可以有效地抑制背景干涉条纹噪声,有利于提高信噪比和系统气体浓度检测精度。
Description
技术领域
本发明属于气体检测技术领域,具体涉及一种应用于气体浓度检测的多频调制方法、气体浓度检测方法及系统。
背景技术
基于可调谐激光的TDLAS技术是近年来发展起来的一种新型光谱检测技术。其采用非接触式的测量方式,可以对待测样本进行精确选择分析,具有检测灵敏度高、分辨率高、实时性好等优点。因此,TDLAS技术已成为气体在线检测的主要手段。
激光在光学元件表面发生多次反射会引起干涉条纹,表现为一种周期性的低频噪声混叠在谐波信号中,对气体浓度特征值检测造成极大影响。干涉条纹受到驱动电流、激光器温度、标准具长度等多种因素的影响。为了能够有效抑制干涉条纹,科研工作者们提出了一些可行的方法。郭心骞等人提出了一种基于经验模态分解的多次去噪方法来降低多光程吸收光谱中的干涉噪声。D.S.Bomse等人提出了一种自适应奇异值分解算法。Alan Fried等人介绍了一种采用压强调制技术的半导体激光器,使多程吸收池中产生的干涉条纹最小化。此外,其他学者还提出了诸如平衡式气体检测法,傅里叶分析法,数字信号处理等方法。这些方法虽然对抑制干涉条纹有一定效果,但是会增加系统软件算法和硬件结构的复杂程度,影响系统实时性,且无法适用于不同检测对象,适应性不强。
发明内容
本发明的目的是提供一种气体浓度检测的多频调制方法、气体浓度检测方法及系统,所述方法在不增加系统复杂度的情况下,从TDLAS-WMS系统调制信号的角度考虑,即在原有的调制信号的基础上添加抖动调制信号,无需新增复杂的噪声抑制算法,利用所述方法可以有效抑制背景干涉条纹噪声,有利于提高信噪比和系统气体浓度检测精度。
一方面,本发明提供一种气体浓度检测的多频调制方法,包括如下步骤:
步骤S1:调节基础调制信号的幅值使得多频调制的气体浓度检测系统输出的二次谐波吸收曲线中的峰值点横坐标与背景干涉条纹包络的极小值点横坐标相同;
其中,在基于TDLAS-WMS单频调制的气体浓度检测系统中加入抖动调制信号得到多频调制的气体浓度检测系统,所述多频调制的气体浓度检测系统输出的二次谐波信号由二次谐波吸收信号叠加背景干涉条纹构成;加入的抖动调制信号的频率初始值和幅值初始值为预设值;
步骤S2:基于抖动调制信号的幅值调节区间、频率调节区间分别调节所述抖动调制信号的幅值和频率。
本发明为了降低背景干涉条纹的影响,在原调制系统中加入了抖动调制信号构成多频调制系统。且研究发现背景干涉条纹是受基础调制信号的幅值m1以及抖动调制信号的幅值m2影响,为了抑制背景干涉条纹噪声,可以通过调节基础调制信号的幅值m1以及抖动调制信号的幅值m2抑制背景干涉条纹。由于基础调制信号的幅值m1不仅会影响背景干涉条纹噪声水平大小,还会影响二次谐波吸收峰值信号,所以先调节m1使二次谐波的峰值点横坐标与背景干涉条纹包络的极小值点的横坐标相同,如此调节m1既可以使吸收信号的值较大,又可以一定程度上降低背景干涉条纹噪声;然后再调节抖动调制信号的幅值m2和频率ω2。其中,调节抖动调制信号的幅值m2和频率ω2时可以是在幅值调节区间、频率调节区间任意取值,也可以按照本申请下文阐述的方法进行取值,优选按照下文方法进行取值,可以更加有效地抑制背景干涉条纹噪声。
进一步优选,所述抖动调制信号的幅值理论最优值至少位于所述幅值调节区间,所述抖动调制信号的幅值理论最优值由零阶贝塞尔函数项达到第一零点处时得到。
进一步优选,所述抖动调制信号的幅值理论最优值如下:
Md=0.383/2lξ
式中,Md为抖动调制信号的幅值理论最优值,l是标准具长度,ξ是激光器频率-电流调谐系数。零阶贝塞尔函数项达到第一零点处时约满足M2=0.383,而M2=2lm2ξ,因此,得到抖动调制信号的幅值m2的幅值理论最优值Md的上述计算公式。
进一步优选,所述幅值调节区间为:[5,Md+30]μA,Md为抖动调制信号的幅值理论最优值。
进一步优选,基于抖动调制信号的幅值调节区间、频率调节区间分别调节所述抖动调制信号的幅值和频率的过程如下:
在所述幅值调节区间的范围以ΔMc等间距调节所述抖动调制信号的幅值,并每次调节幅值后采集N个二次谐波吸收曲线的峰值,再比较不同幅值下二次谐波吸收曲线的峰值的噪声水平,最后选取噪声水平最低时对应的幅值作为抖动调制信号的幅值,N为正整数;
在所述频率调节区间的范围以Δωc等间距调节所述抖动调制信号的频率,并每次调节频率后采集N个二次谐波吸收曲线的峰值,再比较不同幅值下二次谐波吸收曲线的峰值的噪声水平,最后选取噪声水平最低时对应的频率作为抖动调制信号的频率。
其中,用于衡量噪声水平的参数可以是二次谐波峰值的标准差也可以是Allan方差或者是其他参数,本发明对此不进行具体的限定。且为了排除环境因素等干扰,还可以将上述过程分别重复K1、K2次。
进一步优选,所述频率调节区间为:[0.2ω1,0.4ω1],ω1是基础调制信号的频率。
抖动调制信号频率ω2并无理论最优值,当ω2取ω1、2ω1、3ω1等ω1的整数倍时,往往会使二次谐波信号发生畸变;当ω2<0.2ω1时,背景干涉条纹噪声会变大;又考虑到系统硬件负担,ω2的值不宜过大,因此,ω2一般取0.2ω1-0.4ω1。
进一步优选,步骤S2中先调节所述抖动调制信号的幅值,再调节抖动调制信号的频率。
另一方面,本发明体用的一种气体浓度检测的多频调制方法,利用上述多频调制方法确定的基础调制信号的幅值、抖动调制信号的幅值及频率的气体浓度检测系统检测其他气体对象的浓度时,所述多频调制方法包括如下步骤:
首先,搭建利用检测前一类气体时按照所述多频调制方法调制的气体浓度检测系统,并用于检测其他气体对象的浓度;
其次,调节基础调制信号的幅值使得所述气体浓度检测系统输出的二次谐波吸收曲线中的峰值点横坐标与当前背景干涉条纹包络的极小值点横坐标相同。
此外,本发明还提供一种气体浓度检测方法,包括如下步骤:
首先,利用上述多频调制方法调节基础调制信号的幅值、抖动调制信号的幅值及频率;
其次,基于调制后的气体浓度检测系统检测气体浓度。
最后,本发明还提供一种气体浓度检测系统,包括:锁相放大器、激光驱动器、温控模块、激光器、玻璃瓶以及光电探测器,其中,激光驱动器和光电探测器均与所述锁相放大器连接,所述温控模块内置于所述激光驱动器,所述激光驱动器以及温控模块与所述激光器连接,所述玻璃瓶设置于所述激光器的激光发射端一侧并置于激光器和光电探测器之间,玻璃瓶内装有待检测浓度的气体;
所述锁相放大器输出驱动电流信号至激光驱动器,所述驱动电流信号由低频锯齿波、基础调制信号以及抖动调制信号组成,所述基础调制信号的幅值、抖动调制信号的幅值及频率根据上述多频调制方法确定;
激光驱动器从锁相放大器接收驱动电流信号后生成控制电流并与温控模块共同控制激光器发射调制光;
激光器发射的调制光经过玻璃瓶后被光电探测器接收,所述光电探测器将光信号转换为吸收电信号,并输送至锁相放大器进行解调,锁相放大器解调后输出二次谐波信号。
进一步优选,该气体浓度检测系统还包括与锁相放大器连接的数字示波器或数字处理器。其中,数字示波器用于观察二次谐波信号,数字处理器可以用于进行数据处理。
有益效果
一方面,本发明提供的上述多频调制方法是在原有的单频调制系统中加入一个抖动调制信号构成多频调制系统,从而并没有增加多余的硬件元件也没有增加多余的噪声抑制算法,避免了使用复杂的数据处理方式和增设系统元件而导致的系统复杂度的增加。在引入抖动调制信号后,通过调节基础调制信号的幅值、抖动调制信号的幅值及频率可以有效地降低背景条纹干涉噪声,是从更加“本质”的角度切入,抑制了背景干涉条纹噪声。
另一方面,影响多频调制法中m2,ω2的因素是标准具长度,激光器温度,控制电流平均值等系统硬件参数,而与检测对象无关。当硬件系统固定后,这些参数是不变的。所以将m2,ω2调节至最优后,即使面对不同的检测对象,也不用再做调整,此时的系统的背景干涉条纹噪声已经达到最小水平。而传统的噪声抑制算法需要改变的参数多,实现步骤较为复杂,当改变检测对象时,还需要对相应参数做出调整。
附图说明
图1为TDLAS-WMS多频调制法的玻璃瓶内氧气含量检测系统示意图;
图2为基础调制信号幅值m1的调节示意图;
图3为二次谐波峰值数据标准差随抖动调制信号幅值m2的变化示意图;
图4为二次谐波峰值数据标准差随抖动调制信号频率ω2的变化示意图;
图5为不同调制系统下的背景干涉条纹噪声;
图6为不同调制系统下的二次谐波吸收信号。
具体实施方式
下面将结合实施例对本发明做进一步的说明。
如图1所示为一种气体浓度检测系统的架构示意图,本实施例中玻璃瓶中装有氧气,本实施例中以氧气检测为例,其他可行的实施例中,气体检测对象可以是其他气体。
本实施例中气体浓度检测系统包括:锁相放大器、激光驱动器、温控模块、激光器、玻璃瓶、光电探测器、数字示波器以及数字处理器。其中,激光驱动器、光电探测器、数字示波器以及数字处理器均与锁相放大器连接,温控模块内置于激光驱动器,激光驱动器以及温控模块与所述激光器连接,玻璃瓶设置于激光器的激光发射端一侧并置于激光器和光电探测器之间。
锁相放大器输出驱动电流信号至激光驱动器,其驱动电流信号由低频锯齿波、高频正弦波(基础调制信号)和抖动调制信号组成。激光驱动器接受驱动电流信号后将其转换成控制电流并与温控模块共同作用,使激光器工作在电流调制模式下,并使激光中心波长维持在760.8nm处。激光器发出的调制光经过开放光路,穿过含有氧气的玻璃瓶,被光电探测器接收并转换成吸收电信号,送至锁相放大器中进行解调。锁相放大器输出的二次谐波信号最终被送入数字处理器进行数据处理;同时也可以通过数字示波器观察二次谐波信号,进行参数调节。其中,基础调制信号的幅值、抖动调制信号的幅值及频率按照下述的多频调制方法来确定。
二次谐波信号是由二次谐波吸收信号叠加了背景干涉条纹得到,降低背景干涉条纹带来的噪声可以抑制二次谐波基线漂移,提高信噪比和系统检测精度。本发明在原单频调制系统中加入抖动调制信号构成了多频调制系统,加入抖动调制信号后,激光驱动器的控制电流表示为:
i(t)=i0+m1 cosω1t+m2 cosω2t
式中,i(t)表示当前t时刻激光驱动器的控制电流,i0是激光器驱动器输出的平均电流,m1是基础调制信号的幅值,ω1是基础调制信号的频率,它们的初始值设置与原系统单频调制系统的相同;ω2是抖动调制信号的频率,其初始值设置为0.5ω1;m2是抖动调制信号的幅值,其初始值设置为0.5m1。其他可行的实施例中,ω2与m2的初始值可以是其他设置,本发明对此不进行具体的限定。
同时,加入抖动调制信号后,背景干涉条纹可以近似表示为如下:
I0Fcos(4πlv0)J2(2πM1)J0(2πM2)
其中,M1=2lm1ξ,M2=2lm2ξ,I0是激光器的平均输出光强,F是标准具精细度系数,l是标准具长度,ξ是激光器频率-电流调谐系数,v0是激光器平均输出频率,J2(2πM1)是二阶贝塞尔函数项,J0(2πM2)是零阶贝塞尔函数项,M1和M2是为了方便表示定义的两个变量。
从上述背景干涉条纹的表达式可知,调节m2使J0为零可以抑制背景干涉条纹。由于m1还会影响二次谐波吸收峰值信号,所以先调节m1使二次谐波的峰值点横坐标与背景干涉条纹包络的极小值点横坐标相等,如图2所示,选取此时m1的值。对于m2,其理论上的最优值应使J0达到第一零点处,此时M2=0.383,m2的幅值理论最优值Md可通过0.383/2lξ计算获得。然而上述背景干涉条纹公式是近似表达式,在实际应用中,最佳选择是调节抖动调制信号的幅值m2使其在幅值理论最优值的一定范围内变换,同时分析不同的m2值下的二次谐波峰值的噪声水平,进而找到m2的实验最优值。抖动调制信号的频率ω2并无理论最优值,当ω2取ω1、2ω1、3ω1等ω1的整数倍时,往往会使二次谐波信号发生畸变;当ω2<0.2ω1时,背景干涉条纹噪声会变大;又考虑到系统硬件负担,ω2的值不宜过大,因此,ω2一般的取值范围:0.2ω1-0.4ω1。本发明实施例中优选在频率调节范围内调节ω2使其在该范围内变化,同时分析不同ω2值下的二次谐波峰值的误差水平找到ω2的实验最优值。本发明实施例中m2和ω2的调节过程如下:
在幅值调节区间的范围[5,Md+30]μA以ΔMc=1μA等间距调高抖动调制信号的幅值m2,并每次调节幅值m2后采集500个二次谐波吸收曲线的峰值,再比较不同幅值下二次谐波的峰值的噪声水平,最后选取噪声水平最低时对应的幅值作为抖动调制信号的幅值m2,为了排除环境因素等干扰,优选该过程重复K1次,K1=5,其他可行的实施例中,K1可以设置为其他值。本实施例中噪声水平以二次谐波峰值的标准差,其可行的实施例中选择其他参数来表示噪声水平,譬如Allan方差。
在频率调节区间[0.2ω1,0.4ω1]的范围以Δωc=0.01ω1等间距调高抖动调制信号的频率ω2,并每次调节频率ω2后采集500个二次谐波吸收曲线的峰值,再比较不同频率下二次谐波的峰值的噪声水平,最后选取噪声水平最低时对应的频率作为抖动调制信号的频率。同理,为了排除环境因素等干扰,优选该过程重复K2次,K2=2,其他可行的实施例中,K2可以设置为其他值。本实施例中噪声水平以二次谐波峰值的标准差,其可行的实施例中选择其他参数来表示噪声水平,譬如Allan方差。
基于上述理论性描述,本发明实施例提供的一种气体浓度检测的多频调制方法包括如下步骤:
首选,调节基础调制信号的幅值使得多频调制的气体浓度检测系统输出的二次谐波吸收曲线中的峰值点横坐标与背景干涉条纹包络的极小值点横坐标相同;
然后,基于抖动调制信号的幅值调节区间调节所述抖动调制信号的幅值,本实施例中按照上述描述的过程执行,其他可行的实施例中,也可以在幅值调节区间选择幅值。
最后,基于抖动调制信号的频率调节区间调节抖动调制信号的频率,本实施例中按照上述描述的过程执行,其他可行的实施例中,也可以在频率调节区间选择频率。
需要说明的是,本实施例中先调节基础调制信号的幅值,再调节抖动调制信号的幅值,最后调节抖动调制信号的频率。其他可行的实施例中,可以先调节基础调制信号的幅值,至于抖动调制信号的幅值以及抖动调制信号的频率并无顺序要求。
由于影响多频调制法中m2,ω2的因素是标准具长度,激光器温度,控制电流平均值等系统硬件参数,而与检测对象无关。当硬件系统固定后,这些参数是不变的。所以将m2,ω2调节至最优后,即使面对不同的检测对象,也不用再做调整,此时的系统的背景干涉条纹噪声已经达到最小水平。因此,切换不同检测对象时,多频调制方法包括:
首先,搭建前一类气体对象确定的气体浓度检测系统,用于检测新气体对象的浓度;
然后,调节基础调制信号的幅值使得气体浓度检测系统输出的二次谐波吸收曲线中的峰值点横坐标与当前背景干涉条纹包络的极小值点横坐标相同。
从上述可知,在系统硬件平台不变的情况下,切换不同检测对象时,m2,ω2不需要再调整,只需要针对不同检测对象调节基础调制信号的幅值m1。
此外,本发明提供的多频调制方法是应用于氧气检测,本发明提供一种基于上述多频调制方法的气体浓度检测方法,其包括:
首先,利用上述多频调制方法调节系统中的基础调制信号的幅值、抖动调制信号的幅值及频率;
然后,基于调制后的气体浓度检测系统检测气体浓度。需要说明的是,该步骤是利用已知的浓度检测方法来实现。
为了实验本发明所述方法,实验内容如下:
实验中涉及的主要参数有:激光器中心频率v1为13143.2cm-1,氧气谱线中心频率v0为13142.6cm-1,大气压强P为1atm,激光器温度T为296K,开放光程长度L为30cm,玻璃瓶内氧气浓度X为21%,激光器频率-电流调谐系数ξ为0.6cm-1/mA,标准具长度l为20cm,精细度系数F为0.08。扫描锯齿波频率f为25Hz,扫描锯齿波幅值m为10mA,扫描锯齿波直流偏置k为2mA,基础调制信号的频率ω1为2500Hz,基础调制信号的幅值m1初始值设为0.14mA,抖动调制信号的频率ω2初始值设为0.5ω1=1250Hz,抖动调制信号的幅值m2初始值设为0.5m1=70μA。现在开始对m1、m2、ω2进行参数优化调节。
仔细微调m1的值,同时观察示波器中二次谐波吸收峰值点相对于背景干涉条纹包络极小值点的位置变化,若远离,则反向调节m1的值,直到二次谐波吸收峰与背景干涉条纹包络极小值点位置重合,如图2所示,此时m1=0.12mA。固定该值,然后调节m2。
通过公式m2=0.383/2lξ和实验参数,计算得到m2理论最优值为16μA。调节m2,使其在[5μA,46μA]之间以1μA的间隔增加,同时在每个m2值下测量500个二次谐波吸收峰值点,计算并记录这500个数据点的标准差。为排除环境波动带来的影响,将这一过程重复5次,分别对应图3中的实验1-实验5。图3显示了记录的二次谐波峰值数据标准差随m2的变化。在图中,实验3的曲线在m2=18μA时取得最小值,其余4个实验的曲线均在m2=17μA时取得最小值。这两个值都与理论最优值十分接近,根据多数原则,选择m2=17μA,固定该值。
最后调节ω2,使其在[500Hz,1000Hz]之间以25Hz的间隔增加,同时在每个ω2值下测量500个二次谐波吸收峰值点,计算并记录这500个数据点的标准差。为排除环境波动带来的影响,将这一过程重复2次,分别对应图4中的实验1-实验2。图4显示了记录的二次谐波峰值数据标准差随ω2的变化。图中两个实验的曲线重合度很高,在500Hz、550Hz、575Hz、600Hz、750Hz、800Hz、900Hz、950Hz、1000Hz等处,二次谐波信号都有不错的稳定性表现。考虑硬件负担,一般选择低频率,所以ω2=500Hz。
经过上述步骤,已完成多频调制系统的参数优化。现对比单频调制系统和多频调制系统的背景干涉条纹噪声水平。图5显示了无吸收气体条件下的背景干涉条纹2f信号。可以看出单频调制下的背景干涉条纹噪声要显著大于多频调制下的背景干涉条纹噪声,且会发生明显的基线漂移。通过计算得到:单频调制下的背景干涉条纹噪声均值为5.07,方差为2.85,方差均值比为56.2%;多频调制下的背景干涉条纹噪声均值为0.084,方差为0.0048,方差均值比为5.7%。图6显示,多频调制后的2f吸收信号与理想情况下的2f吸收信号拟合度很高。单频系统下的2f吸收信号已发生畸变和基线漂移,对系统测量精度和准确度造成很大影响。通过计算得到:单频系统下的2f吸收信号信噪比(SNR)为3.122dB,多频调制后的2f吸收信号信噪比为36.979dB,信噪比提升了一个数量级。综上所述,验证了本方法可以有效抑制系统干涉条纹噪声和基线漂移,提高系统信噪比和稳定性。
需要强调的是,本发明所述的实例是说明性的,而不是限定性的,因此本发明不限于具体实施方式中所述的实例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,不脱离本发明宗旨和范围的,不论是修改还是替换,同样属于本发明的保护范围。
Claims (10)
1.一种气体浓度检测的多频调制方法,其特征在于:包括如下步骤:
步骤S1:调节基础调制信号的幅值使得多频调制的气体浓度检测系统输出的二次谐波吸收曲线峰值点横坐标与背景干涉条纹包络的极小值点横坐标相同;
其中,在基于TDLAS-WMS单频调制的气体浓度检测系统中加入抖动调制信号得到多频调制的气体浓度检测系统,所述多频调制的气体浓度检测系统输出的二次谐波信号由二次谐波吸收信号叠加背景干涉条纹构成;加入的抖动调制信号的频率初始值和幅值初始值为预设值;
步骤S2:基于抖动调制信号的幅值调节区间、频率调节区间分别调节所述抖动调制信号的幅值和频率。
2.根据权利要求1所述的方法,其特征在于:所述抖动调制信号的幅值理论最优值至少位于所述幅值调节区间内,所述抖动调制信号的幅值理论最优值由零阶贝塞尔函数项达到第一零点处时得到。
3.根据权利要求2所述的方法,其特征在于:所述抖动调制信号的幅值理论最优值如下:
Md=0.383/2lξ
式中,Md为抖动调制信号的幅值理论最优值,l是标准具长度,ξ是激光器频率-电流调谐系数。
4.根据权利要求2所述的方法,其特征在于:所述幅值调节区间为:[5,Md+30]μA,Md为抖动调制信号的幅值理论最优值。
5.根据权利要求1所述的方法,其特征在于:基于抖动调制信号的幅值调节区间、频率调节区间分别调节所述抖动调制信号的幅值和频率的过程如下:
在所述幅值调节区间的范围以ΔMc等间距调节所述抖动调制信号的幅值,并每次调节幅值后采集N个二次谐波吸收曲线的峰值,再比较不同幅值下二次谐波吸收曲线的峰值的噪声水平,最后选取噪声水平最低时对应的幅值作为抖动调制信号的幅值,N为正整数;
在所述频率调节区间的范围以Δωc等间距调节所述抖动调制信号的频率,并每次调节频率后采集N个二次谐波吸收曲线的峰值,再比较不同频率下二次谐波吸收曲线的峰值的噪声水平,最后选取噪声水平最低时对应的频率作为抖动调制信号的频率。
6.根据权利要求1所述的方法,其特征在于:所述频率调节区间为:[0.2ω1,0.4ω1],ω1是基础调制信号的频率。
7.根据权利要求1所述的方法,其特征在于:步骤S2中先调节所述抖动调制信号的幅值,再调节抖动调制信号的频率。
8.一种气体浓度检测的多频调制方法,其特征在于:利用权利要求1所述的方法确定的基础调制信号的幅值、抖动调制信号的幅值及频率的气体浓度检测系统检测其他气体对象的浓度时,所述多频调制方法包括如下步骤:
首先,搭建权利要求1确定的所述气体浓度检测系统,用于检测其他气体对象的浓度;
然后,调节基础调制信号的幅值使得所述气体浓度检测系统输出的二次谐波吸收曲线中的峰值点横坐标与当前背景干涉条纹包络的极小值点横坐标相同。
9.一种气体浓度检测方法,其特征在于:包括如下步骤:
首先,利用权利要求1或者权利要求8所述方法调节基础调制信号的幅值、抖动调制信号的幅值及频率;
其次,基于调制后的气体浓度检测系统检测气体浓度。
10.一种气体浓度检测系统,其特征在于:包括:锁相放大器、激光驱动器、温控模块、激光器、玻璃瓶以及光电探测器,其中,激光驱动器和光电探测器均与所述锁相放大器连接,所述温控模块内置于所述激光驱动器,所述激光驱动器以及温控模块均与所述激光器连接,所述玻璃瓶设置于所述激光器的激光发射端一侧并置于激光器和光电探测器之间,所述玻璃瓶内装有待检测浓度的气体;
所述锁相放大器输出驱动电流信号至激光驱动器,所述驱动电流信号由低频锯齿波、基础调制信号以及抖动调制信号组成,所述基础调制信号的幅值、抖动调制信号的幅值及频率根据权利要求1或者权利要求8所述方法确定;
激光驱动器从锁相放大器接收驱动电流信号后生成控制电流并与温控模块共同控制激光器发射调制光;
激光器发射的调制光经过玻璃瓶后被光电探测器接收,所述光电探测器将光信号转换为吸收电信号,并输送至锁相放大器进行解调,锁相放大器解调后输出二次谐波信号。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010430751.2A CN111398216A (zh) | 2020-05-20 | 2020-05-20 | 气体浓度检测的多频调制方法、气体浓度检测方法及系统 |
PCT/CN2020/101417 WO2021232562A1 (zh) | 2020-05-20 | 2020-07-10 | 气体浓度检测的多频调制方法、气体浓度检测方法及系统和装置 |
GB2214653.4A GB2608554A (en) | 2020-05-20 | 2020-07-10 | Multi-frequency modulation method for gas concentration measurement, gas concentration measurement method and system, and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010430751.2A CN111398216A (zh) | 2020-05-20 | 2020-05-20 | 气体浓度检测的多频调制方法、气体浓度检测方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111398216A true CN111398216A (zh) | 2020-07-10 |
Family
ID=71435679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010430751.2A Pending CN111398216A (zh) | 2020-05-20 | 2020-05-20 | 气体浓度检测的多频调制方法、气体浓度检测方法及系统 |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN111398216A (zh) |
GB (1) | GB2608554A (zh) |
WO (1) | WO2021232562A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112986180A (zh) * | 2021-02-06 | 2021-06-18 | 复旦大学 | 一种光谱式气体传感数据处理方法及系统 |
CN114199820A (zh) * | 2021-12-07 | 2022-03-18 | 北京华亘安邦科技有限公司 | 一种气体浓度的检测方法和装置 |
CN114705654A (zh) * | 2022-04-01 | 2022-07-05 | 华电智控(北京)技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114942236B (zh) * | 2022-05-07 | 2024-07-05 | 山西大学 | 一种消除光强抖动对激光气体遥测影响的方法 |
CN115773996A (zh) * | 2022-11-29 | 2023-03-10 | 中国电子科技集团公司第四十八研究所 | 一种高可靠性激光氧分压检测装置及其使用方法和修正方法 |
CN115791698B (zh) * | 2023-02-02 | 2023-05-23 | 合肥清芯传感科技有限公司 | 气体浓度检测设备及其控制方法、控制装置和存储介质 |
CN115963081B (zh) * | 2023-03-16 | 2023-08-01 | 杭州泽天春来科技有限公司 | 一种气体分析仪 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026830A (ja) * | 2010-07-22 | 2012-02-09 | Shimadzu Corp | ガス濃度測定装置 |
CN104237161A (zh) * | 2014-10-15 | 2014-12-24 | 中国科学院合肥物质科学研究院 | 一种煤炭自燃指标气体多组分实时在线远程监测装置及方法 |
CN105548075A (zh) * | 2016-01-08 | 2016-05-04 | 楚天科技股份有限公司 | 一种玻璃药瓶内氧气含量的检测装置与方法 |
CN107991268A (zh) * | 2018-01-11 | 2018-05-04 | 中国石油大学(华东) | 基于多频调制和多频解调的多气体组份检测方法及装置 |
CN110044844A (zh) * | 2019-05-29 | 2019-07-23 | 中南大学 | 一种基于多次谐波联合分析的气体浓度反演方法及系统 |
CN110160976A (zh) * | 2019-06-06 | 2019-08-23 | 中南大学 | 基于二次谐波干涉包络补偿的气体浓度检测方法和系统 |
CN210005431U (zh) * | 2019-05-22 | 2020-01-31 | 中南大学 | 一种玻璃容器内气体浓度测量装置及玻璃容器质检设备 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011013126A (ja) * | 2009-07-03 | 2011-01-20 | Shimadzu Corp | ガス濃度測定装置 |
JP2016070686A (ja) * | 2014-09-26 | 2016-05-09 | 神栄テクノロジー株式会社 | Tdlas法による濃度測定装置。 |
CN107247034B (zh) * | 2017-05-24 | 2019-11-12 | 东南大学 | 一种基于吸收光谱技术的双频率波长调制方法 |
CN108872147B (zh) * | 2018-06-22 | 2020-07-24 | 中国科学院长春光学精密机械与物理研究所 | 一种基于tdlas技术的激光发生装置及气体浓度检测装置 |
-
2020
- 2020-05-20 CN CN202010430751.2A patent/CN111398216A/zh active Pending
- 2020-07-10 WO PCT/CN2020/101417 patent/WO2021232562A1/zh active Application Filing
- 2020-07-10 GB GB2214653.4A patent/GB2608554A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026830A (ja) * | 2010-07-22 | 2012-02-09 | Shimadzu Corp | ガス濃度測定装置 |
CN104237161A (zh) * | 2014-10-15 | 2014-12-24 | 中国科学院合肥物质科学研究院 | 一种煤炭自燃指标气体多组分实时在线远程监测装置及方法 |
CN105548075A (zh) * | 2016-01-08 | 2016-05-04 | 楚天科技股份有限公司 | 一种玻璃药瓶内氧气含量的检测装置与方法 |
CN107991268A (zh) * | 2018-01-11 | 2018-05-04 | 中国石油大学(华东) | 基于多频调制和多频解调的多气体组份检测方法及装置 |
CN210005431U (zh) * | 2019-05-22 | 2020-01-31 | 中南大学 | 一种玻璃容器内气体浓度测量装置及玻璃容器质检设备 |
CN110044844A (zh) * | 2019-05-29 | 2019-07-23 | 中南大学 | 一种基于多次谐波联合分析的气体浓度反演方法及系统 |
CN110160976A (zh) * | 2019-06-06 | 2019-08-23 | 中南大学 | 基于二次谐波干涉包络补偿的气体浓度检测方法和系统 |
Non-Patent Citations (5)
Title |
---|
LIANG,WK 等: ""Developing CH4 detection limit at lambda=1.654 mu m by suppressing optical interference fringes inwavelength modulation spectroscopy"", 《SENSORS AND ACTUATORS B-CHEMICAL》 * |
LIU,YS 等: ""A New Method for Second Harmonic Baseline Correction and Noise Elimination on Residual Oxygen Detection in Packaged Xilin Bottle"", 《SPECTROSCOPY AND SPECTRAL ANALYSIS》 * |
朱高峰 等: ""WMS的玻璃药瓶内氧气含量检测及其标定方法改进"", 《光谱学与光谱分析》 * |
朱高峰 等: ""基于波长调制光谱技术的氧气浓度检测及其标定方法"", 《分析试验室》 * |
梁文科: ""TDLAS系统中的噪声抑制研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112986180A (zh) * | 2021-02-06 | 2021-06-18 | 复旦大学 | 一种光谱式气体传感数据处理方法及系统 |
CN114199820A (zh) * | 2021-12-07 | 2022-03-18 | 北京华亘安邦科技有限公司 | 一种气体浓度的检测方法和装置 |
CN114199820B (zh) * | 2021-12-07 | 2024-06-07 | 北京华亘安邦科技有限公司 | 一种气体浓度的检测方法和装置 |
CN114705654A (zh) * | 2022-04-01 | 2022-07-05 | 华电智控(北京)技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
CN114705654B (zh) * | 2022-04-01 | 2022-11-18 | 华电智控(北京)技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2021232562A1 (zh) | 2021-11-25 |
GB202214653D0 (en) | 2022-11-16 |
GB2608554A (en) | 2023-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111398216A (zh) | 气体浓度检测的多频调制方法、气体浓度检测方法及系统 | |
US5267019A (en) | Method and apparatus for reducing fringe interference in laser spectroscopy | |
JP5528497B2 (ja) | 圧力依存性を低下させてガス濃度を検出するための方法および装置 | |
KR100186272B1 (ko) | 가스의 적외선 분광 분석방법 및 이것에 사용되는 장치 | |
CN111537470B (zh) | 一种基于数字调制的tdlas气体浓度检测方法 | |
KR100873550B1 (ko) | 가스 검출 방법 및 가스 검출 장치 | |
CN113252204B (zh) | 一种基于幅度调制的激光吸收光谱温度测量方法 | |
US5636035A (en) | Method and apparatus for dual modulation laser spectroscopy | |
CN1908624A (zh) | 气体检测方法和气体检测装置 | |
US6552792B1 (en) | Wavelength modulated photoacoustic spectrometer | |
CN108872146B (zh) | 一种基于tdlas的激光驱动装置和基于tdlas的激光发生装置 | |
US5930000A (en) | Line-locked diode laser for gas spectroscopy | |
CN109696415A (zh) | 一种基于快速傅里叶变换的气体吸收率在线测量方法 | |
CN110160976B (zh) | 基于二次谐波干涉包络补偿的气体浓度检测方法和系统 | |
CN110987870A (zh) | 基于波长调制光谱技术的实时监测气体浓度的系统和方法 | |
US8891085B2 (en) | Gas analyzer | |
US7414727B2 (en) | Gas detection method and gas detection device | |
JP2001210625A (ja) | エッチング深さの検出方法 | |
US9459209B2 (en) | Gas analysis device | |
US5742399A (en) | Method for stabilizing the wavelength in a laser spectrometer system | |
CN115275772B (zh) | 一种基于tdlas技术的特定时刻激光波长控制方法及装置 | |
EP0768524A2 (en) | Method for stabilizing the wavelength in a laser spectrometer system | |
US11650155B2 (en) | Gas analysis system and gas analysis method | |
Zhu et al. | Improved Savitzky-Golay filtering algorithm for measuring a pharmaceutical vial’s oxygen content based on wavelength modulation spectroscopy | |
CN112504981A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200710 |