CN114778475B - 一种可变光程ndir宽量程气体检测方法和系统 - Google Patents
一种可变光程ndir宽量程气体检测方法和系统 Download PDFInfo
- Publication number
- CN114778475B CN114778475B CN202210229949.3A CN202210229949A CN114778475B CN 114778475 B CN114778475 B CN 114778475B CN 202210229949 A CN202210229949 A CN 202210229949A CN 114778475 B CN114778475 B CN 114778475B
- Authority
- CN
- China
- Prior art keywords
- optical path
- gas
- light source
- concentration
- measured
- 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
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 151
- 238000001514 detection method Methods 0.000 title claims abstract description 41
- 238000001745 non-dispersive infrared spectroscopy Methods 0.000 title claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000010521 absorption reaction Methods 0.000 claims abstract description 18
- 230000007613 environmental effect Effects 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 16
- 230000014509 gene expression Effects 0.000 claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims description 24
- 230000004044 response Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000004445 quantitative analysis Methods 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 238000009423 ventilation Methods 0.000 claims description 5
- 230000003750 conditioning effect Effects 0.000 claims description 4
- 230000008901 benefit Effects 0.000 abstract description 5
- 239000007789 gas Substances 0.000 description 134
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 34
- 229910002092 carbon dioxide Inorganic materials 0.000 description 17
- 239000001569 carbon dioxide Substances 0.000 description 16
- 238000000691 measurement method Methods 0.000 description 10
- 239000000523 sample Substances 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000013480 data collection Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 3
- 238000000525 cavity enhanced absorption spectroscopy Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000004867 photoacoustic spectroscopy Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 239000012535 impurity Substances 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 description 2
- 206010010071 Coma Diseases 0.000 description 1
- 206010021143 Hypoxia Diseases 0.000 description 1
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000004847 absorption spectroscopy Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000228 dispersive infrared spectroscopy Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000840 electrochemical analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004817 gas chromatography Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000007954 hypoxia Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000005245 sintering Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/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/01—Arrangements or apparatus for facilitating the optical investigation
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种可变光程NDIR宽量程气体检测方法和系统，仅用一个通道，通过改变光程进行多次测量获得包含多个影响因素参量的多组数据，从而解析出气体浓度，得到气体浓度C和光程差ΔL的关系式，可以避免环境温湿度等因素的干扰。该检测系统包括：光源，用于发射红外光，光谱范围包含被测气体的特征吸收波长；光源驱动模块，用于给光源提供稳定的供电电压以及调制光源；气室，为单通道可变光程气室；探测器模块，用于获取待测气体在不同的气体浓度和光程下的原始电压数据；数据采集处理模块，用于对原始电压进行数据采集并进行数字信号处理。本发明测量范围更大、系统体积更小、成本更低、稳定性更高，为气体检测提供了更加简单准确的方法。
Description
技术领域
本发明涉及气体检测技术领域，具体地，涉及一种可变光程NDIR宽量程气体检测方法和系统。
背景技术
人类的生存离不开气体环境，且各种气体在整个地球生态环境中都有着广泛的分布。以二氧化碳为例，环境中二氧化碳浓度分布范围很宽，根据WMO温室气体公报，全球大气环境中二氧化碳平均摩尔分数为410.5ppm±0.2ppm，人体呼出气体中二氧化碳含量在4％左右(40000pm)，工业二氧化碳排放浓度更高，如大型钢铁厂烧结机废气可达4-7％，当环境中二氧化碳浓度高于5％时，会导致人严重缺氧，昏迷甚至死亡。现有的气体传感技术主要包括气敏法、质谱法、气相色谱法、电化学分析方法和各种光学方法，如色散红外光谱(NDIR)、傅里叶变换红外光谱(FTIR)、光声光谱(PAS)、腔增强吸收光谱(CEAS)以及可调谐二极管吸收光谱(TDLAS)技术等。非光学方法主要通过待测气体的物理化学特性，例如因化学反应产生的电阻、电容或电位的变化等，优势在于工艺简单，成本低，然而他们的长期稳定性以及与杂质气体的交叉敏感性较差。相反，光学方法具有特异性好的优势，但其中FTIR扫描速度慢，易受水汽及杂质气体干扰；PAS容易受其它气体影响，且系统较为复杂，成本较高；CEAS对光路要求高，调试过程复杂，稳定性差；TDLAS需对光源进行精确调制，系统复杂，成本较高。而NDIR具有选择性好，灵敏度高，响应速度快，寿命长且易维修等优势，因而在气体检测领域被广泛应用。
现有技术中的基于NDIR的气体技术大多采用具有参考通道和测量通道的双通道设计，通过测量通道和参考通道的比值以降低部分噪声和基线漂移，提高测量精度，但是参考通道与被测通道采用固定式的光程，其检测浓度范围较小，气室结构比较复杂，调试过程较繁杂。特别是在有些测量中的多通道设计为非共光路系统，抗干扰能力较弱，且测量通道与参考通道中还包含有不同相的噪声，这部分噪声无法通过比值消除，而且两个通道采用的探测器存在响应不一致的问题。此外，因为待测气体浓度测量会受环境因素(比如温度、湿度、气压等)干扰，测量通道和参考通道所测量的波长不一致，无法通过参考通道测量消除环境因素的干扰，所以双通道测量方法受温湿度等环境因素影响较大，需要定期校准零点或增设补偿算法。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种可变光程NDIR宽量程气体检测方法和系统，通过单通道变光程测量消除环境温湿度干扰，并且扩大气体测量范围。
本发明是通过以下技术方案实现的：
根据本发明的第一方面，提供一种可变光程NDIR宽量程气体检测方法，包括：
基于朗伯比尔定律，修正环境参数，建立传感器输出电压与气体浓度及光程差的函数关系；
根据所述函数关系，建立分析模型，并确定所述分析模型的常量；
利用所述分析模型对待测气体的浓度进行定量分析。
进一步地，所述建立传感器输出电压与气体浓度及光程差的函数关系，包括：
选择光谱范围包含被测气体的特征吸收波长的光源，在光源照射下，气体浓度C和光程差L以及传感器输出电压U之间的关系表达式S1为：
其中，K0为光电转换的系数，与传感器型号以及后续调理电路的放大倍数有关；ε是被测气体的摩尔吸收系数；S(L1)为光程长度为L1的损耗函数，S(L2)为光程长度为L2的损耗函数；f(T,R)为环境参数温度T和湿度R对测量结果的干扰函数；L1和L2分别是光程变化前后的光程长度；U1和U2为光程长度为L1和L2时探测器的输出电压。
进一步地，所述根据所述函数关系，建立分析模型，并确定所述分析模型的常量，包括：
在不含被测气体的条件下，测量不同光程下的传感器输出电压，并采用最小二乘法进行优化拟合，从而确定函数S(L)；
建立被测气体在所需测量浓度范围的一系列梯度气体浓度的样本，测量不同光程下，不同气体浓度对应的传感器电压，得到作为建模依据的原始电压数据；
对所述原始电压数据进行处理后，进行多项式优化拟合，从而确定表达式S1中的参数K0和ε，从而建立接收传感器输出电压-气体浓度和光程差的三维曲面关系模型，得到分析模型。
进一步地，所述利用所述分析模型对待测气体的浓度进行定量分析，包括：
确定合适的光程长度L1和L2；
根据所选取的光程长度L1和L2，测量对应光程下的传感器输出电压；
将光程差ΔL＝L2-L1和所述传感器输出电压代入已预先通过定标建立的所述分析模型，获得待测气体的浓度C。
进一步地，所述确定合适的光程长度L1和L2，包括：
以光程最小处Lmin能够测量到的最大气体浓度Cmax的一半为浓度阈值Cth，所述浓度阈值Cth对应的传感器输出电压为电压阈值Uth，首先测量被测气体在光程取最小L1＝Lmin时的传感器响应U1，若测量结果小于电压阈值Uth，则预判断其浓度范围为偏大，反之偏小；
根据预判断结果，若浓度范围偏大，L2选取大于L1小于光程中位数值，反之选取大于光程中位数值，其中光程中位数值为(Lmax-Lmin)/2。
根据本发明的第二方面，提供一种可变光程NDIR宽量程气体检测系统，用于实现上述的方法，包括：
光源，用于发射红外光,所述红外光的光谱范围包含被测气体的特征吸收波长；
光源驱动模块，用于给光源提供稳定的供电电压以及调制光源；
气室，所述气室为单通道可变光程气室；
探测器模块，用于获取待测气体在不同的气体浓度和光程下的原始电压数据；
数据采集处理模块，用于对所述原始电压进行数据采集并进行数字信号处理。
进一步地，所述光源为带有抛物反射面的宽带光源。
进一步地，所述气室为直射型气室，所述气室呈空心圆柱形，所述气室的内径被构造为能够放置轴向对准的光源和探测器，所述气室在空心圆柱的主体部分开设有通气槽。
进一步地，所述气室的表面开设有多个与外界连通的限位孔，相邻两个所述限位孔之间的间隔距离相同。
进一步地，所述探测器为热释电探测器、热电堆探测器和光电二极管中的任一种，所述探测器的前方设有滤光片，所述滤光片的光谱特性与待测气体的红外吸收峰波段匹配。
与现有技术相比，本发明具有如下至少之一的有益效果：
(1)本发明仅用一个通道，通过改变光程进行多次测量来获得包含了多个影响因素参量(如湿度、温度、探测器响应度等)的多组数据，从而解析出气体的浓度，由于温度和湿度等环境因素的影响在解析运算过程中抵消了，所以使用本发明变换光程的测量方法可以避免环境温湿度等因素的干扰。
(2)因为气室采用的是单通道，变光程进行的多次测量属于共光路测量，抗干扰能力强，噪声相同，能够得到比较彻底的消除，并且只需用一个探测器，克服了多个探测器响应不均一的问题，还具有结构简单、量程宽的特点。
(3)本发明的可变光程NDIR宽量程气体检测系统和方法，二氧化碳气体定量分析模型决定系数达0.9929，预测的均方根误差为0.0136，证明本发明的检测系统和方法能显著提高NDIR气体检测方法的量程范围；而且，除了适用于二氧化碳气体，本发明的气体检测系统和方法还适用于其他的气体，具有广泛的应用范围。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明实施例中单通道可变光程测量方法的示意图；
图2为本发明实施例中可变光程NDIR宽量程气体检测系统的结构框图；
图3为本发明实施例中可变光程NDIR宽量程CO2气体检测系统的示意图；
图4为本发明实施例中原始测量需要经过数据处理流程图；
图5为本发明实施例中固定光程下的传感器响应曲线；
图6为本发明实施例中7种气体浓度(0-6000ppm)以及对应的9种光程长度(1-9cm)的实际数据的C-ΔL三维曲面模型；
图7为本发明实施例中温度对测量系统的干扰影响。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明实施例中的可变光程NDIR宽量程气体检测方法，其设计原理如下：
在气体测量中，不同的气体分子具有各自不同的特征红外吸收峰，在光程一定的情况下，吸光度与浓度满足朗伯比尔定律如表达式(1)：
式中A是待测气体的吸光度，I0是入射光强，It是透射光强，ε是被测气体的摩尔吸收系数，C是气体浓度，L是光程长度。基于朗伯比尔定律，可选取两种不同的光程长度L1、L2作为测量自变量：
式中，I0为光源原始光强，I1、I2是经过不同的光程长度L1、L2吸收后透射的光强，经过探测器的光电转换，探测器输出电压U满足：
U＝K×I (4)
综合以上各式，可得
表达式(4)中K为光电转换的系数，与探测器型号以及后续调理电路的放大倍数有关。进一步简化可得出探测器输出电压与气体浓度及光程差的函数关系，满足表达式(5)，其中K1，K2为比例系数，与气体的摩尔吸收系数ε以及光电转换系数K有关，可看作常数。
本发明实施例提供一种可变光程NDIR宽量程气体检测方法，该方法基于上述设计原理实现，参阅图1，该方法包括：
步骤一：基于朗伯比尔定律，修正环境参数，建立传感器输出电压与气体浓度及光程差的函数关系，该步骤目的为修正环境参数，构建实际测量模型。
表达式(5)为不考虑环境温湿度等因素影响以及光在传播路径上损耗的理想化模型，在实际测量中需要考虑这些问题，现对模型进行环境因素干扰的修正。假设温度T和湿度R对待测气体浓度的干扰用函数f(T,R)表示，如表达式(6)；且在不同光程处，光源发出的光到达探测器过程存在不同程度的损耗，假设损耗函数为S(L)，可得到达探测器的光强I'和修正后的探测器输出的电压U'，如表达式(7)和(8)。
C'＝f(T,R)+C (6)
I'＝S(L)I (7)
U'＝K×I'＝K×S(L)×I (8)
那么修正后的模型为式(9)，
该表达式即为在光谱范围包含被测气体的特征吸收波长的光源照射下，气体浓度C和光程差L1-L2以及传感器输出电压U之间的关系表达式。式中，K0为光电转换的系数，与传感器型号以及后续调理电路的放大倍数有关；ε是被测气体的摩尔吸收系数；S(L1)为光程长度为L1的损耗函数，S(L2)为光程长度为L2的损耗函数；f(T,R)为环境参数温度T和湿度R对测量结果的干扰函数；L1和L2分别是光程变化前后的光程长度；U1和U2为光程长度为L1和L2时探测器的输出电压。由表达式(9)可知，温度和湿度的影响在运算过程中抵消了，所以使用变换光程的测量方法可以避免环境温湿度等因素的干扰，无需考虑环境的温湿度干扰，且只需要测量在两种不同的光程下的探测器电压就可以获得待测气体的浓度信息，而传统的双通道测量方法需要定期对温湿度进行漂移标定和校准。所以，变换光程的测量方法只需要更少的通道，而且不容易受环境的影响。
步骤二：建立分析数据模型，确定常量，具体地，根据步骤一中的函数关系，建立分析模型，并确定分析模型的常量。
基于公式(9)，建立分析模型，并确定分析模型的常量，包括：
首先，在不含被测气体的条件(即被测气体浓度为0ppm)下，测量不同光程下的传感器输出电压，并采用最小二乘法进行优化拟合，从而确定公式(9)中的函数S(L)。然后，建立被测气体在所需测量浓度范围的系列梯度气体浓度的样本，测量不同光程下，不同气体浓度对应的传感器电压，得到作为建模依的原始电压数据据；最后，对原始电压数据采用平滑、数字滤波算法等进行处理后，进行多项式优化拟合，从而确定公式(9)中的参数K0,ε。此外，以光程最小处Lmin能够测量到的最大气体浓度Cmax的一半定义为浓度阈值Cth，浓度阈值Cth对应的传感器输出电压定义为电压阈值Uth。
由此建立气体浓度和光程差为自变量，传感器输出电压为因变量的三维曲面关系模型，得到分析模型。
步骤三：利用分析模型对待测气体的浓度进行定量分析。具体地，采用具体的测试装置进行测量，将上述的分析模型应用于气体浓度定量，并基于以上三维曲面模型分析待测气体物质的浓度。
首先确定合适的L1和L2。以光程最小处Lmin能够测量到的最大气体浓度Cmax的一半为浓度阈值Cth，该浓度阈值Cth对应的传感器输出电压为电压阈值Uth。具体的，首先测量被测气体在光程取最小L1＝Lmin时的传感器响应U1，若测量结果小于电压阈值Uth，则预判断其浓度范围为偏大(>Cth)，反之偏小(<Cth)。根据预判断结果，若浓度范围偏大，L2选取大于L1小于光程中位数值，反之选取大于光程中位数值，其中光程中位数值为(Lmax-Lmin)/2。
然后根据所选取L2，测量此光程下传感器输出电压U2。
再将光程差ΔL＝L2-L1和传感器输出电压代入已预先通过定标建立的分析模型，获得待测气体的浓度C。
本发明实施例还提供一种可变光程NDIR宽量程气体检测系统，参阅图2，该系统包括：光源，用于发射红外光,红外光的光谱范围包含被测气体的特征吸收波长；光源驱动模块，用于给光源提供稳定的供电电压以及调制光源；气室，该气室为单通道可变光程气室；探测器模块，用于获取待测气体在不同的气体浓度和光程下的原始电压数据；数据采集处理模块，用于对原始电压进行数据采集并进行数字信号处理。
光源选用宽带光源，其光谱范围包含被测气体的特征吸收波长，为了提高能量利用率，在一些优选的实施例中，光源为带有抛物反射面的宽带光源,采用抛物反射面将反方向传播的光反射到样品通道中。光源驱动模块包括PWM发生器和光源驱动电路，PWM发生器用于产生固定占空比的方波信号传输至电源驱动电路，光源驱动电路用于驱动光源光强按类正弦曲线变化。
在一些优选的实施例中，单通道可变光程气室为直射型气室，结构简单，气室呈空心圆柱形，气室的内径被构造为能够放置轴向对准的光源和探测器，气室在空心圆柱的主体部分开设有通气槽，以保证气室内外气体能够自由交换。气室的圆筒式外壳设计减少光的发散，同时能够固定光源与探头位置。
为实现在与外界交换气体同时便于自由调整光程，气室的表面开设有多个与外界连通的限位孔，相邻两个限位孔之间的间隔距离相同。限位孔除了同通气槽一样具有通气功能外，还具有定位的功能，通过人工操作移动光源到限定位置来实现变化不同的光程。
在一些优选的实施例中，探测器模块包括探头及其驱动电路，探测器为热释电探测器、热电堆探测器和光电二极管中的任一种，探测器的前方设有滤光片，滤光片的光谱特性与待测气体的红外吸收峰波段匹配。
为实现对原始电压的数据采集以及数字信号处理，在一些优选的实施例中，数据采集处理模块包括滤波放大电路、AD转换和数据采集电路。
本实施例中系统还包括电源模块，电源模块通过电压转换电路共分为两部分，一部分用于驱动光源驱动模块，另一部分通过用于驱动探测器模块。
本发明一优选实施例中，CO2气体检测系统的结构图参照图3所示，采用的测量系统包括带有抛物反射面的宽带光源、直射型具有限位孔(与外界连通)的气室、热释电探测器、相应的调制电路模块、不同浓度的标准样品气体、气体密封罐以及计算机。气室材质为ABS材料，外形为空心圆柱，长度选取120mm，内径9.2mm，每个限位孔之间的距离为10mm，气室轴线与光源、探测器轴线完全重合，且光源可沿气室内壁自由移动。12V电源模块通过电压转换电路共分为两部分，一部分降压为0-3.3V，负责驱动光源驱动模块，包括PWM发生器以及光源驱动电路；另一部分通过电压转换电路转化为-5～+5V，驱动探测器模块，包括探头和放大滤波电路；其中，PWM发生器用于产生5Hz的方波，同时方波作为驱动信号供给光源驱动电路驱动光源光强产生周期性类正弦变化。探头满足电压供电要求-5～+5V即可探测到在不同的气体浓度和光程长度下的峰值电压信号强度，且在输出端产生峰值不同的但周期性与光源相同的类正弦电压变化曲线，此原始电压信号经过放大滤波电路后可传至PC端的数字信号处理模块进行数字信号处理。
根据上述原理，本发明实施例中的可变光程NDIR宽量程气体检测系统，仅用一个通道，通过改变光程进行多次测量，获得包含多个影响因素参量的多组数据，从而解析出气体的浓度；同时气室采用的是单通道，变光程进行的多次测量属于共光路测量，抗干扰能力强，噪声相同，能够得到彻底的消除，并且只用一个探测器，克服了探测器响应不均一的问题，还具有结构简单，特别是量程宽的特点。相对于现有技术中的双通道探测系统，本发明测量范围更大、系统体积更小、成本更低、稳定性更高，从而为二氧化碳等气体的检测提供了一种更简单准确的方法。
为了更好说明本发明上述的超宽范围的可变光程NDIR宽量程气体检测系统和方法，以下结合具体的应用实施例来进行描述，以便了解实现的细节，但是应当理解的是，以下实施例并不用于限定本发明。具体的，以二氧化碳的浓度测量为例。
本应用实施例中，二氧化碳的浓度定量分析方法按照以下步骤进行：
第一步，建立气体测量系统。
本步骤中采用的可变光程NDIR宽量程气体检测系统，参阅图3。为了测量二氧化碳的浓度，测量系统的具体结构参数选择为：
光源选择涵盖CO2的红外特征吸收峰4.26μm，自带抛物面式发射面型的宽带光源，减小光的散射；光源驱动模块包括PWM发生器以及光源驱动电路，PWM发生器产生固定占空比为50％的方波信号传输至光源驱动电路，驱动光源光强按类正弦曲线变化。探头选取热释电探测器，搭载宽测量范围二氧化碳滤光片(NBP4.2μm)，该探头具有低温漂以及抑制机械震荡干扰结构。气室，外形为空心圆柱，长度选取120mm，内径9.2mm，气室轴线与光源、探测器轴线完全重合，且光源可沿气室内壁自由移动气体测量系统搭建完成后，进行气体检测。
第二步，为了确定二氧化碳的C-ΔL模型中的参数，进行以下操作：
本发明实施例选取M种不同的浓度的标准气体样品和N种不同的光程长度作为定标自变量，使用单通道可变光程的NDIR宽量程高精度气体检测平台，获得探测器模块输出的原始电压信号。不同浓度测量前需进行洗气操作，即用氮气(纯度达99.99％)冲扫整个气体密封罐(包含气室在内)，然后用真空泵将容器抽至真空，重复以上操作3-4次，接着开始气体测量。对传感器输出的M*N组原始输出电压数据进行平滑、数字滤波算法，提取原始数据中电压峰值强度作为传感器输出电压。以气体浓度和光程长度作为自变量，以传传感器输出电压作为因变量，构建三维曲面定标模型，进行定量分析。
第三步，根据建立的定标模型获取待测气体的浓度，如图4所示的流程图进行浓度测量。
首先测量光程取最小L1＝Lmin时的传感器响应U1，测量结果小于电压阈值Uth预判断其浓度范围为偏大(>Cth)，反之偏小(<Cth)。根据预判结果，若浓度范围偏大，L2选取大于L1小于光程中位数值(Lmax-Lmin)/2即可，反之选取大于光程中位数值。然后根据所选取L2测量此光程下传感器输出电压U2。最后将光程差ΔL＝L2-L1和传感器输出电压代入已预先通过定标建立的气体测量模型，即可获得待测气体的浓度C。
实施例结果分析：
传统固定光程的测量方法中，光程会限制测量气体的浓度检测范围，如图5显示。图5为固定光程L分别为2、4、8cm时的传感器响应变化曲线，可知在长光程下测量范围有限，但是拟合曲线斜率搭说明测量精度高，反之短光程测量范围大，但测量精度小。
二氧化碳气体浓度分别选取(0，500，1000，2000，3000，4000，6000ppm)共七种，光程选取1，2，3，4，5，6，7，8，9cm共9种光程，进行原始电压数据采集，然后对原始电压数据进行处理。对于原始电压数据进行处理包括Goertzel数字滤波、归一化，获得传感器输出电压。利用所有样品的最终传感器输出电压，和对应的样品的浓度与光程，建立定量分析模型即三维曲面定标模型。以气体浓度和光程长度作为自变量，传感器响应作为因变量构建三维曲面定标模型。本实施例中，以7种气体浓度(0，496，1001，2002，3000，4004和6000ppm)作为x轴，光程差(ΔL＝1，2，3，4，5，6，7，8cm)作为y轴，相对应的传感器输出电压作为z轴构建三维曲面模型，并进行三维曲面拟合，拟合二元五次曲面方程的形式为：
z＝p00+p10x+p01y+p20x2+p11xy+p02y2+p21x2y+p12xy2+p03y3+p22x2y2+p13xy3+p04y4+p23x2y3+p14xy4+p05y5
p00＝0.6097；p10＝-0.06183；p01＝-0.201；p20＝0.006908；
p11＝-0.01503；p02＝-0.2439；p21＝0.004374；p12＝0.0387；
p03＝0.09449；p22＝-0.002187；p13＝-0.001685；
p04＝0.133；p23＝-0.001175；p14＝-0.01037；p05＝-0.05822；
其中，x为气体浓度，y为光程差，拟合曲面模型的决定系数R2与均方根误差RMSE的计算公式如式(10)和(11)所示，式中Cnum为定标浓度个数，Lnum为定标测量中光程数目，计算得R2可达0.9929，RMSE为0.0136，如图6所示，证明本发明实施例的检测系统和方法能显著提高NDIR气体检测方法的量程范围。
此外，为了验证模型式(9)的建立过程中消除了环境温湿度的干扰，选取不同温度的实验环境进行对照测量。实验过程中控制环境温度分别为25℃、35℃、45℃、55℃，待测气体选取浓度为2000ppm的标样气体，光程选取L1＝1cm，L2＝6cm，结果如附图7所示。图(7)中横坐标代表测量温度，a和b曲线纵坐标代表固定光程下的传感器原始输出值，说明随着温度增大，在单光程下传感器响应输出呈现递减趋势。采用经过本发明实施例提出的可变光程法进行对比测量方法后，结果如图7曲线(c)所示，纵坐标代表经过本发明实施例测量方法后的传感器响应，几乎不随温度的变化发出显著波动，即消除环境温度的干扰。
实验结果充分证明，本发明实施例提出的方法，相比于单光程测量方法，能获取更宽量程的待测气体信息，且消除环境温度和湿度的干扰(如图7所示)，能显著扩大非色散红外技术的测量范围和测量精度。除了适用于二氧化碳气体，本发明的气体检测系统和方法还适用于其他的气体，具有广泛的应用范围。
本发明上述实施例中的单通道可变光程的NDIR宽量程高精度检测方法和系统，仅用一个通道，通过改变光程进行多次测量来获得包含了多个影响因素参量(如湿度、温度、探测器响应度等)的多组数据，从而解析出气体的浓度，得到气体浓度C和光程差ΔL的关系式。由于温度和湿度等环境因素的影响在解析运算过程中被抵消了，所以使用变换光程的测量方法可以避免环境温湿度等因素的干扰，因为气室采用的是单通道，变光程进行的多次测量属于共光路测量，抗干扰能力强，噪声相同，能够得到彻底的消除，并且只用一个探测器，克服了探测器响应不均一的问题，还具有结构简单，特别是量程宽的特点。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。上述各优选特征在互不冲突的情况下，可以任意组合使用。
Claims (8)
1.一种可变光程NDIR宽量程气体检测方法，其特征在于，包括：
基于朗伯比尔定律，修正环境参数，建立传感器输出电压与气体浓度及光程差的函数关系，包括：
选择光谱范围包含被测气体的特征吸收波长的光源，在光源照射下，气体浓度C和光程差L以及传感器输出电压U之间的关系表达式S1为：
其中，K0为光电转换的系数，与传感器型号以及后续调理电路的放大倍数有关；ε是被测气体的摩尔吸收系数；S(L1)为光程长度为L1的损耗函数，S(L2)为光程长度为L2的损耗函数；f(T,R)为环境参数温度T和湿度R对测量结果的干扰函数；L1和L2分别是光程变化前后的光程长度；U1和U2为光程长度为L1和L2时探测器的输出电压；
根据所述函数关系，建立分析模型，并确定所述分析模型的常量，包括：
在不含被测气体的条件下，测量不同光程下的传感器输出电压，并采用最小二乘法进行优化拟合，从而确定函数S(L)；
建立被测气体在所需测量浓度范围的一系列梯度气体浓度的样本，测量不同光程下，不同气体浓度对应的传感器电压，得到作为建模依据的原始电压数据；
对所述原始电压数据进行处理后，进行多项式优化拟合，从而确定表达式S1中的参数K0和ε，从而建立接收传感器输出电压-气体浓度和光程差的三维曲面关系模型，得到分析模型；
利用所述分析模型对待测气体的浓度进行定量分析。
2.根据权利要求1所述的可变光程NDIR宽量程气体检测方法，其特征在于，所述利用所述分析模型对待测气体的浓度进行定量分析，包括：
确定合适的光程长度L1和L2；
根据所选取的光程长度L1和L2，测量对应光程下的传感器输出电压；
将光程差ΔL＝L2-L1和所述传感器输出电压代入已预先通过定标建立的所述分析模型，获得待测气体的浓度C。
3.根据权利要求2所述的可变光程NDIR宽量程气体检测方法，其特征在于，所述确定合适的光程长度L1和L2，包括：
以光程最小处Lmin能够测量到的最大气体浓度Cmax的一半为浓度阈值Cth，所述浓度阈值Cth对应的传感器输出电压为电压阈值Uth，首先测量被测气体在光程取最小L1＝Lmin时的传感器响应U1，若测量结果小于电压阈值Uth，则预判断其浓度范围为偏大，反之偏小；
根据预判断结果，若浓度范围偏大，L2选取大于L1小于光程中位数值，反之选取大于光程中位数值，其中光程中位数值为(Lmax-Lmin)/2。
4.一种可变光程NDIR宽量程气体检测系统，用于实现权利要求1-3任一项所述的方法，其特征在于，包括：
光源，用于发射红外光，所述红外光的光谱范围包含被测气体的特征吸收波长；
光源驱动模块，用于给光源提供稳定的供电电压以及调制光源；
气室，所述气室为单通道可变光程气室；
探测器模块，用于获取待测气体在不同的气体浓度和光程下的原始电压数据；
数据采集处理模块，用于对所述原始电压进行数据采集并进行数字信号处理。
5.根据权利要求4所述的可变光程NDIR宽量程气体检测系统，其特征在于，所述光源为带有抛物反射面的宽带光源。
6.根据权利要求4所述的可变光程NDIR宽量程气体检测系统，其特征在于，所述气室为直射型气室，所述气室呈空心圆柱形，所述气室的内径被构造为能够放置轴向对准的光源和探测器，所述气室在空心圆柱的主体部分开设有通气槽。
7.根据权利要求6所述的可变光程NDIR宽量程气体检测系统，其特征在于，所述气室的表面开设有多个与外界连通的限位孔，相邻两个所述限位孔之间的间隔距离相同。
8.根据权利要求4所述的可变光程NDIR宽量程气体检测系统，其特征在于，所述探测器为热释电探测器、热电堆探测器和光电二极管中的任一种，所述探测器的前方设有滤光片，所述滤光片的光谱特性与待测气体的红外吸收峰波段匹配。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210229949.3A CN114778475B (zh) | 2022-03-10 | 2022-03-10 | 一种可变光程ndir宽量程气体检测方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210229949.3A CN114778475B (zh) | 2022-03-10 | 2022-03-10 | 一种可变光程ndir宽量程气体检测方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114778475A CN114778475A (zh) | 2022-07-22 |
| CN114778475B true CN114778475B (zh) | 2024-08-20 |
Family
ID=82424079
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210229949.3A Active CN114778475B (zh) | 2022-03-10 | 2022-03-10 | 一种可变光程ndir宽量程气体检测方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114778475B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116183537A (zh) * | 2022-12-30 | 2023-05-30 | 上海交通大学 | 基于差分消元的抗干扰ndir混合气体检测方法和系统 |
| CN116202966A (zh) * | 2023-04-07 | 2023-06-02 | 安徽大学 | 一种光谱吸收式气体检测装置及方法 |
| CN118858216B (zh) * | 2024-01-12 | 2024-12-20 | 北京宝隆泓瑞科技有限公司 | 一种气体浓度检测方法、装置、存储介质及电子设备 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110056784A (zh) * | 2019-05-24 | 2019-07-26 | 昆仑杰信（北京）科技有限责任公司 | 一种带压气体全天候泄露监测与即时报警光电探测系统 |
| CN113406036A (zh) * | 2021-06-17 | 2021-09-17 | 桂林电子科技大学 | 基于光谱重构技术的便携式温室气体检测系统 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030215791A1 (en) * | 2002-05-20 | 2003-11-20 | Applied Spectral Imaging Ltd. | Method of and system for multiplexed analysis by spectral imaging |
| US20120078532A1 (en) * | 2010-09-24 | 2012-03-29 | David Edward Forsyth | Non-dispersive infrared sensor measurement system and method |
-
2022
- 2022-03-10 CN CN202210229949.3A patent/CN114778475B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110056784A (zh) * | 2019-05-24 | 2019-07-26 | 昆仑杰信（北京）科技有限责任公司 | 一种带压气体全天候泄露监测与即时报警光电探测系统 |
| CN113406036A (zh) * | 2021-06-17 | 2021-09-17 | 桂林电子科技大学 | 基于光谱重构技术的便携式温室气体检测系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114778475A (zh) | 2022-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114778475B (zh) | 一种可变光程ndir宽量程气体检测方法和系统 | |
| CN110044824B (zh) | 一种基于石英音叉的双光谱气体检测装置及方法 | |
| CN109991189B (zh) | 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法 | |
| US20160132617A1 (en) | Target Analyte Detection and Quantification in Sample Gases With Complex Background Compositions | |
| Cui et al. | Photonic sensing of the atmosphere by absorption spectroscopy | |
| CN106483088A (zh) | 一种基于紫外光调制的气体浓度测量装置及方法 | |
| CN109085133B (zh) | 一种基于实时反射率修正的离轴积分腔大气ch4浓度的测量装置及其测量方法 | |
| CN109813639B (zh) | 一种基于红外光调制技术的颗粒物与气体浓度同步测量装置及其测量方法 | |
| Shen et al. | Methane near-infrared laser remote detection under non-cooperative target condition based on harmonic waveform recognition | |
| CN103852446B (zh) | 一种基于光腔衰荡光谱技术的血液成分识别与分析仪器 | |
| Sun et al. | Stack emission monitoring using non-dispersive infrared spectroscopy with an optimized nonlinear absorption cross interference correction algorithm | |
| Matveev et al. | Midinfrared light-emitting diodes based on А3В5 heterostructures in gas-analyzer-equipment engineering: potential and applications in 2014–2018 | |
| Hu et al. | Triple-range adjustable plane-concave multipass gas cell for high-precision trace gas detection based on TDLAS | |
| Sun et al. | High accuracy wide range CO2 detection method based on difference optical path NDIR | |
| CN106596444A (zh) | 一种基于紫外宽带吸收光谱的氧气浓度测量系统及测量方法 | |
| Dang et al. | A compact self-calibration atmospheric CH4 sensor based on NIR laser and new type multi-pass gas cell | |
| CN114526831A (zh) | 一种露点霜点温度传感器及其测量方法 | |
| Xu et al. | Micro quartz crystal tuning fork based photothermal spectroscopy for trace gas detection | |
| CN206740638U (zh) | 一种并联式气体光谱分析双气室 | |
| CN101281124B (zh) | 宽带腔增强吸收光谱大气环境光电监测系统 | |
| CN206772801U (zh) | 一种串联式气体光谱分析双气室 | |
| Xu et al. | Design of nitrogen oxide detection system based on non-dispersive infrared technology | |
| Yu et al. | A temperature-compensated CO2 detection system based on non-dispersive infrared spectral technology | |
| CN202092947U (zh) | 烟气气体含量激光在线检测系统的光轴调节机构 | |
| He et al. | Research on high-precision oxygen partial pressure sensor based on TDLAS technology |
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 |