CN108801977B - 无标定痕量碳12和碳13二氧化碳气体探测装置及方法 - Google Patents

无标定痕量碳12和碳13二氧化碳气体探测装置及方法 Download PDF

Info

Publication number
CN108801977B
CN108801977B CN201810359053.0A CN201810359053A CN108801977B CN 108801977 B CN108801977 B CN 108801977B CN 201810359053 A CN201810359053 A CN 201810359053A CN 108801977 B CN108801977 B CN 108801977B
Authority
CN
China
Prior art keywords
gas
pressure
harmonic
measured
temperature
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
CN201810359053.0A
Other languages
English (en)
Other versions
CN108801977A (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201810359053.0A priority Critical patent/CN108801977B/zh
Publication of CN108801977A publication Critical patent/CN108801977A/zh
Application granted granted Critical
Publication of CN108801977B publication Critical patent/CN108801977B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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

本发明无标定痕量碳12和碳13二氧化碳气体探测装置及方法，ICL带间级联激光器发射出的红外光经过反射镜折射到多通池中；压力传感器和温度传感器分别将采集到的压强信息和温度信息传输给DSP数据采集模块，红外光从多通池出射之后由光电探测器进行接收，光电探测器将红外光的光功率信号转换为电信号，同时将电信号发送给锁相放大器；锁相放大器接收光电探测器向其发送的电信号并提取电信号中的一次谐波和二次谐波，并将提取的一次谐波和二次谐波发送给DSP数据采集模块；本发明方法是用DSP软补偿的方法取代了复杂的硬件装置，使得系统的使用寿命更长，使其适用于更广泛的应用环境。

Description

无标定痕量碳12和碳13二氧化碳气体探测装置及方法
技术领域
本发明是一种无标定的对痕量气体进行化学探测的装置及方法，尤其是对12CO2/13CO2痕量气体探测的装置及方法。
背景技术
红外检测吸收光谱技术(LAS)是一个新兴的学科，在大气的环境检测，气体泄漏，工业生产控制等领域都有着诸多的应用，它对气体的探测有很好的分辨，并且更加的灵敏，相比于传统的技术，红外检测系统的体积更加的小巧，稳定性强，响应速度快，所以具有更好的检测实时性。近年来半导体激光技术迅速发展，尤其是带间级联激光器(ICL)的发展配合中红外吸收光谱技术和多反射气体吸收气室，对波长调制技术的高灵敏度的气体检测应用提供了一种新的可能。其中波长调制技术(TDLAS-WMS)是随着半导体激光器发展起来的一种红外光谱技术，随着半导体激光器的不断完善和发展，波长调制技术在对矿下有毒气体的检测，大气的污染等都有着非常广泛的应用前景，并且波长调制技术和谐波检测技术的相结合，可以实现对气体分子的高灵敏度检测。波长调制技术就是在直接吸收法的基础之上提取吸收线的一次谐波和二次谐波，进行进一步的分析和浓度信息的提取。原理是利用半导体激光器的电流调谐特性，使得半导体激光器的中心波长在待测气体的吸收波段内，通过改变注入半导体激光器的电流，不仅实现激光器的扫描和调制，而且将引起激光器光功率的变化，也就是说波长调制还伴随着强度调制，可以有效的降低检测中的1/f噪声。
(S.Neethu,R.Verma,S.S.Kamble,et al.Validation of wavelengthmodulation spectroscopy techniques for oxygen concentrationmeasurement.Sensors and Actuators B:Chemical,vol.192,2014)S.Neethu团队利用中心波长为760nm(25℃，70mA)的可调谐分布反馈激光器，通过光程为56cm的多通池(Analytical Corporation,USA)，之后光束打到光功率计上(Newport,USA)，最后由锁相放大器提取出吸收信号的一次谐波和二次谐波，通过二次谐波与一次谐波的比值提取出氧气的气体浓度信息。
(Chunguang Li,Lei Dong,Chuantao Zheng,et al.Compact TDLAS basedoptical sensor for ppb-level ethane detection by use of a 3.34m room-temperature CW interband cascade laser.Sensors and Actuators B,vol.232,2016)，研究小组对中心波长为3.4μm的CW(连续波长调制)ICL(带间级联)激光器进行电流的调协和调制，并通过光程为57.6m的多通池，吸收之后打到光电探测器上，最后由锁相放大器提取一次谐波和二次谐波得到浓度信息。
(Kyle Owen·Aamir Farooq.A calibration-free ammonia breath sensorusing a quantum cascade laser with WMS 2f/1f,Appl.Phys.B 2014,116:371–383)此小组用的是波长为9063nm的量子级联(QCL)激光器作为光源，通过光程为76.45m的多通池，通过光电探测器进行光电转换，最后由锁相放大器提取一次谐波与二次谐波的浓度信息。
上述方法都是利用TDLAS-WMS-2f/1f的检测方法对浓度信息进行提取，此方法测量精度高，并消除了光功率误差带来的干扰，但都有着一定的缺陷，根据比尔朗伯定律知道气体的温度和压强对气体的吸收线有着直接的影响，所以现有技术中都用温度控制器和压力控制器对温度和压强进行了控制和标定，但由于长时间的使用温度控制器和压力控制器都会不同程度的产生飘移，则会对系统产生不同程度的误差。
发明内容
本发明的目的在于针对传统的TDLAS-WMS气体探测方法由于温控和压控的长期使用，会使其产生误差和飘移的不足，提供无标定痕量碳12和碳13二氧化碳气体探测装置及方法。
本发明的目的通过以下技术方案实现的：无标定痕量碳12和碳13二氧化碳气体探测装置，其特征在于：包括ICL带间级联激光器、反射镜、多通池、压力传感器、光电探测器、温度传感器、锁相放大器、DSP数据采集模块及ICL激光器电流驱动装置，
所述ICL带间级联激光器与ICL激光器电流驱动装置连接，ICL激光器电流驱动装置用于驱动ICL带间级联激光器发射出红外光；
所述反射镜位于ICL带间级联激光器出射光路上，ICL带间级联激光器发射的红外光经过反射镜折射进入到多通池；
所述多通池内充满待测气体，多通池的进气口与配气装置相连接,并在多通池的进气口与配气装置连接处设置有压力传感器，在多通池侧壁上设置有温度传感器；
所述压力传感器用于检测多通池内部气压，压力传感器与DSP数据采集模块通信连接，其将采集到的压强信息传输给DSP数据采集模块；
所述温度传感器用于检测多通池内部待测气体温度，温度传感器与DSP数据采集模块通信连接，其将采集到的温度信息传输给DSP数据采集模块；
所述光电探测器用于接收从多通池出射的红外光，并将红外光的光功率信号转换为电信号，同时将电信号发送给锁相放大器；
所述锁相放大器用于产生驱动电流，驱动电流是由调制电流和调谐电流组成，锁相放大器将其产生的调制电流和调谐电流发送给ICL激光器电流驱动装置，同时锁相放大器接收光电探测器向其发送的电信号并提取电信号中的一次谐波和二次谐波，并将提取的一次谐波和二次谐波发送给DSP数据采集模块；
所述ICL激光器电流驱动装置接收锁相放大器向其发送的调制电流和调谐电流并进行叠加。
无标定痕量碳12和碳13二氧化碳气体探测方法，其特征在于：该方法采用上述的探测装置进行无标定痕量12CO2/13CO2气体探测，包括以下步骤：
A、ICL带间级联激光器在ICL激光器电流驱动装置的驱动下发射出红外光，红外光经过反射镜,最后折射到多通池中；
B、红外光经过折射后进入到多通池，多通池里面充满了CO2和N2的混合气体，压力传感器和温度传感器分别采集多通池内部压强信息和温度信息并传输给DSP数据采集模块，红外光从多通池出射之后由光电探测器进行接收，光电探测器将红外光的光功率信号转换为电信号，同时将电信号发送给锁相放大器；
C、锁相放大器接收光电探测器发送的电信号并进行分析和处理，预先设置锁相放大器锁定信号的频率，锁相放大器提取电信号中的一次谐波和二次谐波，并将提取的一次谐波和二次谐波发送给DSP数据采集模块；
D、DSP数据采集模块同时接收锁相放大器向其发送的一次谐波和二次谐波，压力传感器向其发送的压强信息及温度传感器向其发送的温度信息，经DSP数据采集模块处理后得到待测气体浓度，实现无标定痕量12CO2/13CO2气体的探测。
步骤D所述待测气体浓度的获得过程如下：
根据比尔朗伯定律，有
Figure BDA0001635476270000041
其中I(λ)是气体吸收之后的光强，I0(λ)是初始光强，
Figure BDA0001635476270000042
为吸收系数，C是待测气体浓度，L是吸收路径的长度，即总光程，λ是吸收波长，将比尔朗伯定律展开则有
Figure BDA0001635476270000043
分别提取出一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，
I1f＝nI0 (2)
Figure BDA0001635476270000044
利用
Figure BDA0001635476270000045
算法对待测气体浓度C进行求解，有
Figure BDA0001635476270000046
其中
Figure BDA0001635476270000047
为单一气体在吸收线中心的吸收截面，n为光强调制系数，J为调制深度；
N＝C/P (5)
其中P是待测气体压强，N是待测气体体积分数比；
Figure BDA0001635476270000048
的大小与温度、压强有关，
A、当压强大于1atm时，
Figure BDA0001635476270000049
满足Lorenz线形，则
Figure BDA00016354762700000410
的表达式有
Figure BDA00016354762700000411
谱线强度S的表达式为：
Figure BDA0001635476270000051
其中S0是温度为T0时的吸收谱线强度，T0取273K，k是玻尔兹曼常数，E″是气体的低能态能量,T是待测气体温度；
ΔνL是Lorenz线形下的谱线半宽，表达式为：
Figure BDA0001635476270000052
其中Δν0是标准状态下的谱线半宽，标准状态是压强为P0及温度为T0，P0为1个标准大气压，T0为273K，P是待测气体压强，T是待测气体温度；
B、当压强小于0.01atm时，
Figure BDA0001635476270000057
满足Gaussion线形，则
Figure BDA0001635476270000058
的表达式有
Figure BDA0001635476270000053
ΔνD是Gaussion线形下的谱线半宽，表达式为：
Figure BDA0001635476270000054
其中M为待测气体分子的相对分子质量，ν0为吸收波数；
C、当压强在0.01atm～1atm时，
Figure BDA0001635476270000059
满足Voigt线形，则
Figure BDA00016354762700000510
的表达式有
Figure BDA0001635476270000055
其中Voigt线形下的谱线半宽Δν和t的表达是分别为：
Figure BDA0001635476270000056
Figure BDA0001635476270000061
其中u是传输路径的光程，若多通池3内的压强恒定，则u大小等于L，t为与吸收分子数目成正比的一个参数；
根据压力传感器4测到的待测气体压强P数值，判断所测压强大小所在的范围，将待测气体压强P和待测气体温度T代入到其对应的表达式中，便可求得谱线半宽和谱线强度S,将谱线半宽和谱线强度S代入到对应的表达式中，便可得到
Figure BDA0001635476270000063
最后将
Figure BDA0001635476270000062
代入到式(4)中，便可求出待测气体浓度C，通过公式(5)既可求出待测气体的体积分数比N，上述方法即为无标定痕量12CO2/13CO2气体探测的方法。
通过上述设计方案，与现有技术相比本发明可以带来如下有益效果：本发明提出的无标定气体探测方法相比于传统的TDLAS-WMS气体探测方法，少了对气体的温度和压强的控制，传统的TDLAS-WMS气体探测方法由于温控和压控的长期使用，会使其产生误差和飘移，只能手动对其调整和重新标定，这样就减小了系统的使用寿命，而本发明中的温度传感器和压力传感器的体积更小，设置更简单，并用DSP软补偿的方法取代了复杂的硬件装置，使得系统的使用寿命更长，使其适用于更广泛的应用环境；并且传统的TDLAS-WMS气体探测方法只能探测单一的气体浓度信息，而本发明可以同时探测到12CO2/13CO2两种气体的浓度信息。
附图说明
图1本发明无标定痕量碳12和碳13二氧化碳气体探测装置的结构框图。
图2本发明实施例中12CO2和13CO2的选定吸收线示意图，在气体温度和压强分别为0℃和20torr的条件下测定。
图中各标记如下：1-ICL带间级联激光器、2-反射镜、3-多通池、4-压力传感器、5-光电探测器、6-温度传感器、7-锁相放大器、8-DSP数据采集模块、9-ICL激光器电流驱动装置。
具体实施方式
为了避免混淆本发明的实质，公知的方法、过程、流程及元件并没有详细的叙述。
下面结合附图和实施例对本发明做进一步的详细说明：
随着中红外激光器的发展，中红外激光气体探测技术也日益趋向成熟，其中TDLAS-WMS技术更是得到了广泛的发展与应用，本发明提供的装置就是在此技术基础之上提供了一种无标定的气体探测方法，把气体的温度和压强从硬补偿变为了软补偿，并对12CO2/13CO2进行测量。
对于气体探测第一步的工作就是进行吸收线的选取，由图2所示，12CO2的吸收线的吸收波数为2315.4/cm，13CO2吸收线的吸收波数为2315.2/cm,之所以要选取这两条吸收谱线，是因为两条吸收谱线的距离挨的非常近，在ICL带间级联激光器1的谱线调谐范围之内，而且没有产生谱线的交叠干扰，并且水线的干扰也可以忽略不计，并根据实际工程经验，采用TDLAS技术可探测到的最小气体吸收线强度为0.001，又由图可知在标定气体温度和压强分别为0℃、20torr的条件下90ppbv体积分数比的12CO2和90ppbv体积分数比13CO2的吸收强度分别为0.008和0.002，所以根据换算，12CO2和13CO2的探测下限分别为11ppbv和45ppbv左右，本发明提供的方案不对温度和压强进行标定，而是通过测得的气体温度和压强并根据与吸收线强度的关系对计算公式内的参数进行补偿，从而得到待测气体浓度，具体的补偿方法由下文可知。
在吸收线选线之后，就大致确定了中心波长的范围，最终选择了中心波长为4319nm(2315.4/cm)的ICL带间级联激光器1，并作为整个系统的光源，ICL带间级联激光器1的驱动电流是由调谐电流和调制电流组成的，调谐电流是低频的锯齿波，调制电流是高频正弦波，由于驱动电流的大小与ICL带间级联激光器1的输出频率成线性关系，所以调谐电流起到了一个扫频的作用，而调制电流则是为了降低1/f噪声的干扰，即闪烁噪声的干扰。
如图1所示，ICL带间级联激光器1发出红外光后经过反光镜2折射后进入到多通池3的通光孔，多通池3内充满着待测气体，红外光在多通池3内经过多次反射之后从多通池3的出射口射出，并且多通池3内的气体温度和压强会由温度传感器6和压力控制器4传给DSP数据采集模块8进行数据处理；红外光进入到光电探测器5后，光电探测器5会把红外光的光功率信号转换为电信号，并把转换完成的电信号传给锁相放大器7进行分析和处理，锁相放大器7会根据之前的设置提取出一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)(谐波的基频就是驱动电流中高频调制信号的频率)，并且此装置的这个驱动电流(锯齿波调谐电流与正弦波调制电流)也是由锁相放大器7提供的，并由ICL激光器电流驱动装置9进行叠加。
锁相放大器7提取出的一次谐波的光强I1f(λ)、二次谐波的光强I2f(λ)随压力传感器4检测到的压强信息及温度传感器6检测到的温度信息一并传给DSP数据采集模块8进行数据的分析和处理，根据比尔朗伯定律，有
Figure BDA0001635476270000081
其中I(λ)是气体吸收之后的光强，I0(λ)是初始光强，
Figure BDA0001635476270000082
为吸收系数，C是待测气体浓度，L是吸收路径的长度，即总光程，λ是吸收波长，将比尔朗伯定律展开则有
Figure BDA0001635476270000083
分别提取出一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，
I1f＝nI0 (2)
Figure BDA0001635476270000084
利用
Figure BDA0001635476270000085
算法对待测气体浓度C进行求解，有
Figure BDA0001635476270000086
其中
Figure BDA0001635476270000087
为单一气体在吸收线中心的吸收截面，n为光强调制系数，J为调制深度；
N＝C/P (5)
其中P是待测气体压强，N是待测气体体积分数比。
因为
Figure BDA0001635476270000091
的大小与温度、压强有关，并且在不同的压强下，
Figure BDA0001635476270000092
的表达形式不一样所以需要分情况进行讨论：
A、当压强大于1atm时，
Figure BDA0001635476270000093
满足Lorenz线形，则
Figure BDA0001635476270000094
的表达式有
Figure BDA0001635476270000095
谱线强度S是只与温度有关的参数，它的关系表达式为：
Figure BDA0001635476270000096
其中S0是温度为T0时的吸收谱线强度，T0取273K，k是玻尔兹曼常数，E″是气体的低能态能量,T是待测气体温度；
ΔνL是Lorenz线形下的谱线半宽，表达式为：
Figure BDA0001635476270000097
其中Δν0是标准状态下的谱线半宽，标准状态是压强为P0及温度为T0，P0为1个大气压，T0为273K，P是待测气体压强，T是待测气体温度；
B、当压强小于0.01atm时，
Figure BDA0001635476270000098
满足Gaussion线形，则
Figure BDA0001635476270000099
的表达式有
Figure BDA00016354762700000910
ΔνD是Gaussion线形下的谱线半宽，表达式为：
Figure BDA00016354762700000911
其中M为待测气体分子的相对分子质量，ν0为吸收波数；
C、当压强在0.01atm～1atm时，
Figure BDA00016354762700000912
满足Voigt线形，则
Figure BDA00016354762700000913
的表达式有
Figure BDA0001635476270000101
其中Voigt线形下的谱线半宽Δν和t的表达是分别为：
Figure BDA0001635476270000102
Figure BDA0001635476270000103
其中u是传输路径的光程，若多通池3内的压强恒定，则u大小等于L，t为与吸收分子数目成正比的一个参数；
将待测气体压强P和待测气体温度T代入到其对应的表达式中，便可求得三种情况下的谱线半宽(Δν,ΔνL或ΔνD)和谱线强度S,将谱线半宽和谱线强度S代入到对应的表达式中，便可得到
Figure BDA0001635476270000104
最后将
Figure BDA0001635476270000105
代入到式(4)中，便可求出待测气体浓度C，通过公式(5)既可求出待测气体的体积分数比N，上述方法即为无标定痕量12CO2/13CO2气体探测的方法。
实施例1
ICL带间级联激光器1是nanoplus公司生产的4319nm的ICL带间级联激光器1，驱动电路的调谐电流是45mA～50mA并且频率为10Hz的锯齿波，调制电流是频率为5000Hz正弦波，并且通过ICL激光器电流驱动装置9把调谐电流和调制电流进行叠加，ICL带间级联激光器1受到电流驱动后发出的红外光经过反光镜2折射到多通池3中，此多通池3的单路光程定为15cm，总光程L为24m，所以红外光需要在多通池3内经过130次～140次反射后射出，反射镜2是镀金的凸面镜，折射率能达到99.99％，在多通池3的玻璃壁上装着一个型号为pt1000的温度传感器6，来感知多通池3内气体的温度，并将气体的温度信息传输给DSP数据采集模块8，在多通池3的进气口前端连接着一个Bronkhorst公司生产的IQ+FLOW系列的压力传感器4,并把压力信息传输给DSP数据采集模块8进行分析和处理；红外光从多通池3射出之后，射入到光电探测器5上，光电探测器5是vigo公司生产的中红外碲镉汞探测器，具有很高的光电转化率，从光电探测器5发出的电信号传输给锁相放大器7，锁相放大器7是上海昕虹公司生产的数字锁相放大器，通过上位机对锁相放大器7进行参数的设置，提取出输入电信号的一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并把一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)传输给DSP数据采集模块8，DSP数据采集模块8在接收温度信息、压强信息、一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并根据以上公式便可求得待测气体的体积分数比，其中配气装置配得的气体(12CO2,13CO2和N2的混合气体)的体积分数比分别是8ppmv的12CO2和96ppbv的13CO2；由公式(4)有：
Figure BDA0001635476270000111
其中测得的二次谐波的光强I2f(λ)峰值大小和一次谐波的光强I1f(λ)的中心值大小分别为920和8，J是一个设定参数，设定值的大小为3.43，总光程L为24m,光强调制系数n的值为0.63，测得的温度T为300K,将温度参数代入到公式(7)中
Figure BDA0001635476270000112
其中T0是273K,玻尔兹曼常数k的大小为1.38×10-23，基于Hitran数据库进行仿真，气体的低能态能量E″的大小为1.387×10-23，并在273K的温度下S0的大小为1.87×10-20，将以上参数代入到公式(6)中得到谱线强度S的大小为1.78×10-20。又测得气体的压强为800torr,经判断其满足条件具体实施方式A,由公式(8)可知，谱线半宽ΔνL：
Figure BDA0001635476270000121
其中T0和T分别为273K和300K,P和P0分别为800torr和20torr,基于Hitran数据库进行仿真，在273K,20torr的条件下的Δν0为0.03，将其带入到公式(8)中有ΔνL的值为1.23，将求得的S和ΔνL分别带入到公式(6)中
Figure BDA0001635476270000122
可得到
Figure BDA0001635476270000123
的大小为4.6×10-21，并将
Figure BDA0001635476270000124
代入到公式(4)中
Figure BDA0001635476270000125
求得的待测气体浓度C的大小为1.91×1020，单位是(分子数/m3)，又通过公式(5)算出体积比分数为8.06ppmv。
以上方法是对12CO2进行的探测，用同样的方法检测此待测气体13CO2的浓度，体积分数比为96.12ppbv。
实施例2
ICL带间级联激光器1是nanoplus公司生产的4319nm的ICL带间级联激光器1，驱动电路的调谐电流是45mA～50mA并且频率为10Hz的锯齿波，调制电流是频率为5000Hz正弦波，并且通过ICL激光器电流驱动装置9把调谐电流和调制电流进行叠加，ICL带间级联激光器1受到电流驱动后发出的红外光经过反光镜2折射到多通池3中，此多通池3的单路光程定为15cm，总光程L为24m,所以红外光需要在多通池3内经过130次～140次反射才能射出，反射镜2是镀金的凸面镜，折射率能达到99.99％，在多通池3的玻璃壁上装着一个型号为pt1000的温度传感器6，来感知多通池3内气体的温度，并将气体的温度信息传输给DSP数据采集模块8，在多通池3的进气口前端连接着一个Bronkhorst公司生产的IQ+FLOW系列的压力传感器4,并把压力信息传输给DSP数据采集模块8进行分析和处理；红外光从多通池3射出之后，射入到光电探测器5上，光电探测器5是vigo公司生产的中红外碲镉汞探测器，具有很高的光电转化率，从光电探测器5发出的电信号传输给锁相放大器7，锁相放大器7是上海昕虹公司生产的数字锁相放大器，通过上位机对锁相放大器7进行参数的设置，提取出输入电信号的一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并把一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)传输给DSP数据采集模块8，DSP数据采集模块8在接收温度信息、压强信息、一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并根据以上公式便可求得待测气体的体积分数比，其中配气装置配得的气体(12CO2,13CO2和N2的混合气体)的体积分数比分别是0.7ppmv的12CO2和8ppbv的13CO2；由公式(4)有：
Figure BDA0001635476270000131
其中测得的二次谐波的光强I2f(λ)峰值大小和一次谐波的光强I1f(λ)的中心值大小分别为920和8，J是一个设定参数，设定值的大小为3.43，总光程L为24m,光强调制系数n的值为0.63，测得的温度T为300K,将温度参数代入到公式(7)中
Figure BDA0001635476270000132
其中T0是273K,玻尔兹曼常数k的大小为1.38×10-23，基于Hitran数据库进行仿真，得到E″的大小为1.387×10-23，并在273K的温度下S0的大小为1.87×10-20，将以上参数代入到公式(7)中得到S的大小为1.78×10-20。又测得气体的压强为5torr,经判断其满足条件具体实施方式B,由公式(10)可知，谱线半宽ΔνD有：
Figure BDA0001635476270000141
其中T是300K，M是44，ν0为2315.4/cm，将其带入到公式(10)中有ΔνD的值为0.0022，将求得的S和ΔνD分别带入到公式(9)中
Figure BDA0001635476270000142
可得到
Figure BDA0001635476270000143
的大小为1.75×10-18，并将
Figure BDA0001635476270000144
代入到公式(4)中
Figure BDA0001635476270000145
求得的待测气体浓度C的大小为1.01×1016，单位是(分子数/m3),通过公式(5)算出体积分数比为0.71ppmv。
以上方法是对12CO2进行的探测，用同样的方法检测此待测气体13CO2的浓度，体积分数比为8.13ppbv。
实施例3
ICL带间级联激光器1是nanoplus公司生产的4319nm的ICL带间级联激光器1，驱动电路的调谐电流是45mA～50mA并且频率为10Hz的锯齿波，调制电流是频率为5000Hz正弦波，并且通过ICL激光器电流驱动装置9把调谐电流和调制电流进行叠加，ICL带间级联激光器1受到电流驱动后发出的红外光经过反光镜2折射到多通池3中，此多通池3的单路光程定为15cm，总光程L为24m,所以红外光需要在多通池3内经过130次～140次反射才能射出，反射镜2是镀金的凸面镜，折射率能达到99.99％，在多通池3的玻璃壁上装着一个型号为pt1000的温度传感器6，来感知池内气体的温度，并将气体的温度信息传输给DSP数据采集模块8，在多通池3的进气口前端连接着一个Bronkhorst公司生产的IQ+FLOW系列的压力传感器4,并把压力信息传输给DSP数据采集模块8进行分析和处理；光从多通池3射出之后，射入到光电探测器5上，光电探测器5是vigo公司生产的中红外碲镉汞探测器，具有很高的光电转化率，从光电探测器5发出的电信号传输给锁相放大器7，锁相放大器7是上海昕虹公司生产的数字锁相放大器，通过上位机对锁相放大器7进行参数的设置，提取出输入电信号的一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并把一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)传输给DSP数据采集模块8，DSP数据采集模块8在接收温度信息、压强信息、一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，并根据以上公式便可求得待测气体的体积分数比，其中配气装置配得的气体(12CO2,13CO2和N2的混合气体)的体积分数比分别是0.95ppmv的12CO2和10ppbv的13CO2；由公式(4)有：
Figure BDA0001635476270000151
其中测得的二次谐波的光强I2f(λ)峰值大小和一次谐波的光强I1f(λ)的中心值大小分别为920和8，J是一个设定参数，设定值的大小为3.43，总光程L为24m，光强调制系数n的值为0.63，测得的温度T为300K,将温度参数代入到公式(7)中，
Figure BDA0001635476270000152
其中T0是273K,玻尔兹曼常数k的大小为1.38×10-23，基于Hitran数据库进行仿真，E″的大小为1.387×10-23，并在273K的温度下S0的大小为1.87×10-20，将以上参数代入到公式(7)中得到S的大小为1.78×10-20。又测得气体的压强为100torr,经判断其满足条件具体实施方式C,由公式(12)可知，谱线半宽Δν有：
Figure BDA0001635476270000161
此公式中的ΔνL和ΔνD，是将测得的温度和压强分别带入到公式(8)和(10)中得到的，整理后得到的Δν为582.79，将得到的Δν代入到公式(11)中
Figure BDA0001635476270000162
可得到
Figure BDA0001635476270000163
的大小为1.54×10-19，并将
Figure BDA0001635476270000164
代入到公式(4)中
Figure BDA0001635476270000165
求得的待测气体浓度C的大小为3.00×1017，单位是(分子数/m3),通过公式(5)算出体积分数比为0.96ppmv。
以上方法是对12CO2进行的探测，用同样的方法检测此待测气体13CO2的浓度，体积比分数为10.18ppbv。

Claims (1)

1.无标定痕量碳12和碳13二氧化碳气体探测方法，
采用无标定痕量碳12和碳13二氧化碳气体探测装置，包括ICL带间级联激光器(1)、反射镜(2)、多通池(3)、压力传感器(4)、光电探测器(5)、温度传感器(6)、锁相放大器(7)、DSP数据采集模块(8)及ICL激光器电流驱动装置(9)，
所述ICL带间级联激光器(1)与ICL激光器电流驱动装置(9)连接，ICL激光器电流驱动装置(9)用于驱动ICL带间级联激光器(1)发射出红外光；
所述反射镜(2)位于ICL带间级联激光器(1)出射光路上，ICL带间级联激光器(1)发射的红外光经过反射镜(2)折射进入到多通池(3)；
所述多通池(3)内充满待测气体，多通池(3)的进气口与配气装置相连接,并在多通池(3)的进气口与配气装置连接处设置有压力传感器(4)，在多通池(3)侧壁上设置有温度传感器(6)；
所述压力传感器(4)用于检测多通池(3)内部气压，压力传感器(4)与DSP数据采集模块(8)通信连接，其将采集到的压强信息传输给DSP数据采集模块(8)；
所述温度传感器(6)用于检测多通池(3)内部待测气体温度，温度传感器(6)与DSP数据采集模块(8)通信连接，其将采集到的温度信息传输给DSP数据采集模块(8)；
所述光电探测器(5)用于接收从多通池(3)出射的红外光，并将红外光的光功率信号转换为电信号，同时将电信号发送给锁相放大器(7)；
所述锁相放大器(7)用于产生驱动电流，驱动电流是由调制电流和调谐电流组成，锁相放大器(7)将其产生的调制电流和调谐电流发送给ICL激光器电流驱动装置(9)，同时锁相放大器(7)接收光电探测器(5)向其发送的电信号并提取电信号中的一次谐波和二次谐波，并将提取的一次谐波和二次谐波发送给DSP数据采集模块(8)；
所述ICL激光器电流驱动装置(9)接收锁相放大器(7)向其发送的调制电流和调谐电流并进行叠加；
其特征在于：包括以下步骤：
A、ICL带间级联激光器(1)在ICL激光器电流驱动装置(9)的驱动下发射出红外光，红外光经过反射镜(2),最后折射到多通池(3)中；
B、红外光经过折射后进入到多通池(3)，多通池(3)里面充满了CO2和N2的混合气体，压力传感器(4)和温度传感器(6)分别采集多通池(3)内部压强信息和温度信息并传输给DSP数据采集模块(8)，红外光从多通池(3)出射之后由光电探测器(5)进行接收，光电探测器(5)将红外光的光功率信号转换为电信号，同时将电信号发送给锁相放大器(7)；
C、锁相放大器(7)接收光电探测器(5)发送的电信号并进行分析和处理，预先设置锁相放大器(7)锁定信号的频率，锁相放大器(7)提取电信号中的一次谐波和二次谐波，并将提取的一次谐波和二次谐波发送给DSP数据采集模块(8)；
D、DSP数据采集模块(8)同时接收锁相放大器(7)向其发送的一次谐波和二次谐波，压力传感器(4)向其发送的压强信息及温度传感器(6)向其发送的温度信息，经DSP数据采集模块(8)处理后得到待测气体浓度，实现无标定痕量12CO2/13CO2气体的探测；
步骤D所述待测气体浓度的获得过程如下：
根据比尔朗伯定律，有
Figure FDA0002595395180000021
其中I(λ)是气体吸收之后的光强，I0(λ)是初始光强，
Figure FDA0002595395180000022
为吸收系数，C是待测气体浓度，L是吸收路径的长度，即总光程，λ是吸收波长，将比尔朗伯定律展开则有
Figure FDA0002595395180000023
分别提取出一次谐波的光强I1f(λ)和二次谐波的光强I2f(λ)，
I1f(λ)＝nIo (2)
Figure FDA0002595395180000024
利用
Figure FDA0002595395180000025
算法对待测气体浓度C进行求解，有
Figure FDA0002595395180000031
其中
Figure FDA0002595395180000032
为单一气体在吸收线中心的吸收截面，n为光强调制系数，J为调制深度；
N＝C/P (5)
其中P是待测气体压强，N是待测气体体积分数比；
Figure FDA0002595395180000033
的大小与温度、压强有关，
当压强大于1atm时，
Figure FDA0002595395180000034
满足Lorenz线形，则
Figure FDA0002595395180000035
的表达式有
Figure FDA0002595395180000036
谱线强度S的表达式为：
Figure FDA0002595395180000037
其中S0是温度为T0时的吸收谱线强度，T0取273K，k是玻尔兹曼常数，E″是气体的低能态能量,T是待测气体温度；
ΔνL是Lorenz线形下的谱线半宽，表达式为：
Figure FDA0002595395180000038
其中Δν0是标准状态下的谱线半宽，标准状态是压强为P0及温度为T0，P0为1个标准大气压，T0为273K，P是待测气体压强，T是待测气体温度；
当压强小于0.01atm时，
Figure FDA0002595395180000039
满足Gaussion线形，则
Figure FDA00025953951800000310
的表达式有
Figure FDA00025953951800000311
ΔνD是Gaussion线形下的谱线半宽，表达式为：
Figure FDA0002595395180000041
其中M为待测气体分子的相对分子质量，ν0为吸收波数；
当压强在0.01atm～1atm时，
Figure FDA0002595395180000042
满足Voigt线形，则
Figure FDA0002595395180000043
的表达式有
Figure FDA0002595395180000044
其中Voigt线形下的谱线半宽Δν和t的表达是分别为：
Figure FDA0002595395180000045
Figure FDA0002595395180000046
其中u是传输路径的光程，若多通池(3)内的压强恒定，则u大小等于L，t为与吸收分子数目成正比的一个参数；
根据压力传感器(4)测到的待测气体压强P数值，判断所测压强大小所在的范围，将待测气体压强P和待测气体温度T代入到其对应的表达式中，便可求得谱线半宽和谱线强度S,将谱线半宽和谱线强度S代入到对应的表达式中，便可得到
Figure FDA0002595395180000047
最后将
Figure FDA0002595395180000048
代入到式(4)中，便可求出待测气体浓度C，通过公式(5)即可求出待测气体的体积分数比N，上述方法即为无标定痕量12CO2/13CO2气体探测的方法。
CN201810359053.0A 2018-04-20 2018-04-20 无标定痕量碳12和碳13二氧化碳气体探测装置及方法 Active CN108801977B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810359053.0A CN108801977B (zh) 2018-04-20 2018-04-20 无标定痕量碳12和碳13二氧化碳气体探测装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810359053.0A CN108801977B (zh) 2018-04-20 2018-04-20 无标定痕量碳12和碳13二氧化碳气体探测装置及方法

Publications (2)

Publication Number Publication Date
CN108801977A CN108801977A (zh) 2018-11-13
CN108801977B true CN108801977B (zh) 2020-09-18

Family

ID=64093380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810359053.0A Active CN108801977B (zh) 2018-04-20 2018-04-20 无标定痕量碳12和碳13二氧化碳气体探测装置及方法

Country Status (1)

Country Link
CN (1) CN108801977B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109444068B (zh) * 2018-12-29 2020-12-15 郎溪杰博电器科技有限公司 一种红外碳硫分析仪的模糊预测控制分析系统
CN111504945B (zh) * 2020-06-08 2023-06-13 朗思传感科技(深圳)有限公司 一种光纤光热气体传感装置及方法
CN112834562B (zh) * 2021-01-04 2022-04-12 吉林大学 一种热导式混合气体中氦气浓度检测的装置及其方法
CN113533207B (zh) * 2021-07-27 2022-08-09 吉林大学 一种基于tdlas技术的高准确度检测装置及修正方法
CN115494020B (zh) * 2022-09-30 2024-02-23 湖南五凌电力科技有限公司 痕量气体检测方法、装置及系统
CN117664901A (zh) * 2023-12-14 2024-03-08 深圳市诺安智能股份有限公司 一种基于可调谐滤光片的多气体传感器及气体检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104126106A (zh) * 2011-12-20 2014-10-29 小利兰·斯坦福大学托管委员会 一种用于气体感测的免校准扫描波长调制光谱的方法
CN106654844A (zh) * 2016-12-30 2017-05-10 中国科学院合肥物质科学研究院 一种基于室温qcl激光器的同位素探测在线锁频装置和方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104126106A (zh) * 2011-12-20 2014-10-29 小利兰·斯坦福大学托管委员会 一种用于气体感测的免校准扫描波长调制光谱的方法
CN106654844A (zh) * 2016-12-30 2017-05-10 中国科学院合肥物质科学研究院 一种基于室温qcl激光器的同位素探测在线锁频装置和方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ICL-based TDLAS sensor for real-time breath gas analysis of carbon monoxide isotopes;RAMIN GHORBANI 等;《OPTICS EXPRESS》;20170523;第25卷(第11期);全文 *
免标定的TDLAS氧气在线测量系统;康信文;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20180415;第7-10、17-18、37页 *
激光气体检测及温压补偿研究;黄云彪;《仪器仪表用户》;20150831;第22卷(第4期);"1 TDLAS 浓度测量原理"、"3 浓度信号的温压修正"部分 *

Also Published As

Publication number Publication date
CN108801977A (zh) 2018-11-13

Similar Documents

Publication Publication Date Title
CN108801977B (zh) 无标定痕量碳12和碳13二氧化碳气体探测装置及方法
CN104903703B (zh) 气体吸收分光装置以及气体吸收分光方法
CN109490216B (zh) 一种免校准的激光光声光谱微量气体检测仪器及方法
US11346774B2 (en) Infra red spectrometer
CN109991189B (zh) 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法
CN104730021A (zh) 一种非分光红外气体传感器的标定方法
Wang et al. Enhancing off-axis integrated cavity output spectroscopy (OA-ICOS) with radio frequency white noise for gas sensing
CN103616347B (zh) 多气体谱线检测的最佳调制系数实现方法和装置
Wang et al. Laser frequency locking and intensity normalization in wavelength modulation spectroscopy for sensitive gas sensing
CN109813639B (zh) 一种基于红外光调制技术的颗粒物与气体浓度同步测量装置及其测量方法
CN102353633A (zh) 烟气气体含量激光在线检测方法及系统
CN114945816A (zh) 用于对空芯光纤中的气体进行实时、高分辨率分子分析的受激拉曼光谱仪
CN110927100A (zh) 用于测量气体通量的系统和测量气体通量的方法
CN112748087A (zh) 基于温度调谐吸收光谱技术的多组分检测装置及方法
CN110320178A (zh) 一种基于波长调制光谱技术的气体检测系统及其检测方法
Song et al. Interband cascade laser-based ppbv-level mid-infrared methane detection using two digital lock-in amplifier schemes
CN112067582A (zh) 基于腔增强吸收光谱技术探测水汽稳定同位素装置及方法
US7064329B2 (en) Amplifier-enhanced optical analysis system and method
Li et al. Mid-infrared telemetry sensor based calibration gas cell for CO detection using a laser wavelength locking technique
Lu et al. Bias-immune absorption spectrum extraction by downsampling the wavelength modulation laser signal for precise measurement of dynamic gas temperature
CN107389609A (zh) 一种基于多模激光光谱技术的一氧化碳气体浓度检测方法
CN114397273B (zh) 基于二次-四次谐波联用的气体浓度测量装置及测量方法
CN207636485U (zh) 波长捷变的可调谐半导体激光吸收光谱气体检测装置
CN108181261A (zh) 基于太赫兹时域光谱检测混合气体各组分含量的装置
CN113533249A (zh) 一种基于开放光路和波长调制的n2o检测系统及方法

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Yanzhang

Inventor after: Ren Qiang

Inventor after: Chen Chen

Inventor after: Zheng Chuantao

Inventor after: Liu Zhiwei

Inventor after: Xie Hongtao

Inventor before: Chen Chen

Inventor before: Ren Qiang

Inventor before: Zheng Chuantao

Inventor before: Wang Yanzhang

Inventor before: Liu Zhiwei

Inventor before: Xie Hongtao

GR01 Patent grant
GR01 Patent grant