CN110865042A - 气体浓度检测方法、装置和*** - Google Patents
气体浓度检测方法、装置和*** Download PDFInfo
- Publication number
- CN110865042A CN110865042A CN201911172382.5A CN201911172382A CN110865042A CN 110865042 A CN110865042 A CN 110865042A CN 201911172382 A CN201911172382 A CN 201911172382A CN 110865042 A CN110865042 A CN 110865042A
- Authority
- CN
- China
- Prior art keywords
- frequency
- gas
- light intensity
- sub
- harmonic
- 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.)
- Withdrawn
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 22
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000010521 absorption reaction Methods 0.000 claims description 29
- 230000003287 optical effect Effects 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 10
- 239000007789 gas Substances 0.000 description 98
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 16
- 230000003321 amplification Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 238000003916 acid precipitation Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000004868 gas analysis Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000009965 odorless effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 description 1
- 238000005406 washing 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/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/59—Transmissivity
-
- 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
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/129—Using chemometrical methods
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
本公开涉及一种气体浓度检测方法、装置和***，其中，所述方法包括获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频；确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；基于所述二次谐波信号得到所述待检测浓度的气体的浓度。本公开实施例可以简单且方便的检测出气体的浓度。
Description
技术领域：
本公开涉及信号处理技术领域，尤其涉及一种气体浓度检测方法、装置、和***。
背景技术：
近些年，随着技术的进步和现代工业的快速发展，人们的生活水平不断提高，但同时也对环境产生了负面的影响，大气污染导致了酸雨、温室效应、雾霾等有害的自然现象，不仅会对社会的工业和农业造成直接或间接的经济损害，更是会对人们的健康产生致命的威胁。甲烷是沼气、天然气，以及煤矿和坑道气的主要成分，是一种无色无味的可燃性气体并且在空气中的含量达到一定程度会有发生***的危险。另外在大气中甲烷及二氧化碳等气体浓度所占比例逐渐的增大导致了地球表面气温的上升，从而造成了一系列自然灾害问题，如地球两端南极和北极冰川融化和海平面上升等。综上所述，气体浓度检测在众多领域被广泛的应用，具有十分重要的意义。
发明内容：
本公开提出了一种气体浓度检测方法、装置和***，能够简单方便的得到待检测气体的浓度。
根据本公开的一方面，提供了一种气体浓度检测方法，包括：
获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后得到的透射光的光强信号；
对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频；
确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
在一些可能的实施方式中，所述对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，包括：
利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
基于确定的每个子带的频率范围，确定所需频率对应的子带。
在一些可能的实施方式中，所述确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号，包括：
利用第一预设方式得到所述所需频率对应的子带的频域小波系数；
对所述频域小波系数执行傅里叶逆变换，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络；
利用所述包络得到透射光强信号中的二次谐波信号。
在一些可能的实施方式中，所述基于所述二次谐波信号得到所述待检测浓度的气体的浓度，包括：
获取所述入射光的光参数；
利用第二预设方式，基于所述二次谐波信号和所述光参数，得到所述气体浓度。
在一些可能的实施方式中，所述第二预设方式的表达式为：
其中，An(v0)为透射光强信号的n次谐波信号幅值，C为待检测气体的气体浓度，L为气体的吸收光程长，I0为入射光的光强，
为气体在频率v处吸收系数， V0为气体的中心频率。
根据本公开的第二方面，提供了一种气体浓度检测装置，包括：
获取模块，用于获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；
傅里叶变换模块，用于对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
谐波小波变换模块，用于对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频，以及确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
确定模块，用于基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
在一些可能的实施方式中，所述谐波小波变换模块还用于利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
基于确定的每个子带的频率范围，确定所需频率对应的子带。
在一些可能的实施方式中，所述谐波小波变换模块还用于，利用预设方式得到所述所需频率对应的子带的频域小波系数；
对所述频域小波系数执行傅里叶逆变换，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络；
利用所述包络得到透射光强信号中的二次谐波信号。
在一些可能的实施方式中，所述确定模块还用于获取所述入射光的光参数，利用预设方式，基于所述二次谐波信号和所述光参数，得到所述气体浓度。
根据本公开的第三方面，提供了一种气体浓度检测***，包括：
气体吸收池，用于存储待检测气体；
激光器，用于根据调制信号向所述气体吸收池发射入射光；
探测器，用于检测从气体吸收池透射出的透射光的透射光强信号；
处理器，用于执行如第一方面中任意一项所述的气体浓度检测方法，得到所述待检测气体的气体浓度。
在本公开实施例中，可以将透射光强信号进傅里叶变换，得到其频谱，然后选取频域小波系数，再将频域小波系数进行FFT逆变换，从而得到二次谐波分量，进而得到待检测气体的浓度。本公开实施例通过频域的滤波处理，可以滤掉高频电子噪声、低频环境变化，及激光器输出光强波动对测量结果的影响，提高测量***的抗干扰能力。与现有技术相比，该方法可简化算法，同时省略锁相放大技术中低通滤波参数设置的复杂性，避免了滤波器参数设置不匹配对二次谐波检测灵敏度的影响。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，而非限制本公开。
根据下面参考附图对示例性实施例的详细说明，本公开的其它特征及方面将变得清楚。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，这些附图示出了符合本公开的实施例，并与说明书一起用于说明本公开的技术方案。
图1示出根据本公开实施例的一种气体浓度检测方法的流程图；
图2示出根据本公开实施例的一种气体浓度检测方法中求取二次谐波信号的流程图；
图3示出根据本公开实施例的气体浓度检测***的示意图；
图4示出根据本公开实施例的一种气体浓度检测方法中步骤S30的流程图；
图5示出根据本公开实施例的气体浓度检测装置的框图。
具体实施方式
以下将参考附图详细说明本公开的各种示例性实施例、特征和方面。附图中相同的附图标记表示功能相同或相似的元件。尽管在附图中示出了实施例的各种方面，但是除非特别指出，不必按比例绘制附图。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在 B这三种情况。另外，本文中术语“至少一种”表示多种中的任意一种或多种中的至少两种的任意组合，例如，包括A、B、C中的至少一种，可以表示包括从A、 B和C构成的集合中选择的任意一个或多个元素。
另外，为了更好地说明本公开，在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解，没有某些具体细节，本公开同样可以实施。在一些实例中，对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述，以便于凸显本公开的主旨。
本公开实施例的气体浓度检测方法的执行主体可以是任意的数据处理装置，例如，气体浓度检测方法可以由终端设备或服务器或其它处理设备执行，其中，终端设备可以为用户设备(User Equipment，UE)、移动设备、用户终端、终端、蜂窝电话、无绳电话、个人数字处理(Personal Digital Assistant，PDA)、手持设备、计算设备、车载设备、可穿戴设备等。在一些可能的实现方式中，该气体浓度检测方法可以通过处理器调用存储器中存储的计算机可读指令的方式来实现。
图1示出根据本公开实施例的一种气体浓度检测方法的流程图，图2示出根据本公开实施例的一种气体浓度检测方法中求取二次谐波信号的流程图，如图1 和2所示，所述气体浓度检测方法包括：
S10：获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；
S20：对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
S30：对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频；
S40：确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
S50：基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
在一些可能的实施方式中，可以向充满待检测气体的容器(如图3所示的气体吸收池300)发射入射光，并通过探测器接收透射光，并获得透射光对应的光强信号。探测器可以为光电探测器，可以将光信号转换为电信号。图3示出根据本公开实施例的气体浓度检测***的示意图。其中，激光器200可以在激光驱动器100的驱动下，发射出入射光，该入射光的光强、波长等信息可以根据激光驱动器100的生成的驱动信号确定。本领域技术人员可以根据需求调整激光驱动器发射的驱动信号。本公开实施例中，待检测气体为甲烷气体，选用近红外半导体激光器，通过该激光器对近红外处1.65μm的甲烷吸收谱线进行测量示。为便于安装连接，激光光源使用中心波长为1653.7nm的单模蝶形封装DFB(nanoplus) 激光器。其工作温度范围是在-25-50℃，典型值28℃时的工作波长为1653.7nm，驱动电流最大不超过120mA，并且典型值为70mA，输出功率是5mW。
激光器在驱动信号的驱动下发出入射光，该入射光照射到气体吸收池300内的待检测浓度的气体，经过该气体洗手池300内的气体。发出的透射光被探测器 400接收，探测器400将该透射光的光信号转换为电信号形式的透射光强信号进行处理。其中，本公开实施例中的待检测气体可以为甲烷气体，在其他实施例中，也可以为其他气体，本公开对此不作具体限定。
其中，本公开实施例采用的吸收池为反射型吸收池，其参数如表1所示，激光光束(入射光)在吸收池内形成多次反射，在很小空间内使有效光程达到数十米，实现在小容积气体样品中对光线的高效吸收。该气体吸收池非常稳固，抗光束走离，安装简单，易于操作，非常适宜实验室痕量气体分析和工业过程气体检测。其主要突出特点如下所列：
(1)两种气室状态，可在开放状态或封闭状态间切换；
(2)气封闭状态下，气压可低至数十帕斯卡，有效降低谱线间重叠干扰。
表1 Herriot气体吸收池的参数
本公开实施例选用的探测器未400PDA20CS(-EC)InGaAs Switchable GainDetector型号的高精度放大探测器，具有0-70dB的放大设置，光谱响应度曲线范围宽(覆盖800至1700nm)、暗电流低、可靠性较高等特点。
激光调制参数对二次谐波信号影响复杂，调制参数包含扫描频率、调制系数和调制频率等三个参数。它们共同决定了激光光源的工作状态：扫描波长宽度和扫描时间。另外待测信号的特征也被调制参数所决定，不同的参数对应检测到的信号幅度和形状都可能会不同。因此调制参数确定并设定之后尽量不要在修改。本公开实施例使用的锯齿波扫描输出波长时需完全覆盖谐波信号，并且应尽可能减少引入不必要干扰信号的引入，所以扫描范围必须选取合适。本实验中激光器电流的中心值为70mA，因此应在其左右范围内设定扫描区间，经多次试验测试后，选定的65mA-85mA为最佳电流区间。检测信号的峰值和信噪比主要受扫描频率影响，信号的峰值和信噪比会随着扫描频率的升高而单调减小。因此选择较小的扫描频率可以在浓度测量中提高信号的信噪比，扫描频率选取的范围通常为几-十几Hz。本公开实施例经过多次测试后，选择2Hz的扫描频率。
波长调制可以抑制1/f噪声并提高检测***灵敏度。随着调制频率的增大检测器的1/f噪声跟着减小，但达到一定值时减小的不再明显，因此仪器的选择很重要。根据本实验装置的性能指标，通过多次的实验测试，本公开选取的调制频率是31.2kHz。
在一些可能的实施方式中，处理器500可以从探测器400接收透射光强信号，根据该透射光强信号得到气体吸收吃内的待检测气体的气体浓度。下面，对气体浓度的检测过程进行详细说明。
其中，首先可以对透射光强信号执行傅里叶变换处理(如快速傅里叶变换)，得到透射光强信号对应的频域信息，即得到透射光强信号的频谱。本公开实施例中，探测器输出的透射光强信号可以表示成：
I(v)＝I0(v)·exp[-α(v)CL]≈I0(v)·[-α(v)CL]；
其中，I(v)表示输出的透射光的光强信号(透射光强信号)，v为调制频率， I0(v)表示入射光的光强信号，α(v)表示吸收系数，C待检测气体浓度，L表示光程长度。
进而对透射光强信号执行傅里叶变换(FFT)可以得到频谱。其中，在得到频谱的情况下，可以进一步根据谐波小波包分解的方式，得到所需频率对应的子带。图4示出根据本公开实施例的一种气体浓度检测方法中步骤S30的流程图，其中，所述对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，包括：
S31：利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
S32：基于确定的每个子带的频率范围，确定所需频率对应的子带。
在一些可能的实施方式中，可以对得到的频谱进行谐波小波包分解得到多层的频率信号，其中可以根据采样频率和所需频带宽度(设定的子带宽度)来确定谐波小波包分解层数J，或者也可以通过预先配置好的层数信息确定分解层数 J以及设定的子带宽度B，其中可以通过S表示每个子带的子带号，其中第j层的子带个数为2j，对应的，第j层的每个子带的带宽为B＝fm/2j，其中信号的最高分析频率fm是采样频率fs的一半，因此第j层的每个子带的上限m和下限n 分别为：
其中，所需频率是指生成所述入射光的调制频率的二倍频，其中，生成入射光的调制频率是指激光驱动器对激光器进行驱动时采用的调制信号的调制频率。例如，调制频率可以约为31.2kHz，所需频率可以为62.4kHz，用正弦信号和锯齿波对激光器波长进行高频调制，调制后的光强信号对气体进行吸收被探测器接收，产生高频调制的电压信号(透射光强信号)。
在一些可能的实施方式中，在得到每个子带对应的频率范围的情况下，可以确定所需频率对应的子带，得到该对应的子带的子带号s。进一步的，可以得到该对应的子带的频率范围对应的频域小波系数。本公开实施例可以利用第一预设方式得到所述所需频率对应的子带的频域小波系数，其中，第一预设方式的表达式可以为：
Fm表示第s个子带的频域小波系数，aj，m表示第j层的频域小波展开系数。
表示第j层的所需频率对应的子带的频谱，t表示时间。基于上述，可以得到所需频率对应的子带的频域小波系数。
在得到上述频域小波系数的情况下，可以对该频域小波系数执行傅里叶逆变换，得到时域信号，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络，该包络即可以为二次谐波信号。
进一步地，可以利用第二预设方式基于二次谐波信号确定待检测气体的浓度，其中，第二预设方式的表达式为：
其中，An(v0)为透射光强信号的n次谐波信号幅值，C为待检测气体的气体浓度，L为气体的吸收光程长，I0为入射光的光强，
为气体在频率v处吸收系数， V0为气体的中心频率，其中，本公开实施例中得到的二次谐波信号中，n取值为 2。
下面举例说明本公开实施例的气体浓度检测的处理过程，气体吸收池对应的长光程吸收池(有效吸收光程14.5m)。***中激光光源使用中心频率为 1653.7nm的单模蝶形封装DFB(nanoplus)激光器。激光输出波长由温度和恒流源控制，激光波长随电流的变化率为0.011nm/mA。恒流源、生成三角波及正弦波的信号发生器芯片均封装在TDLAS控制器(PortInstruments，LLC，PCI-1D(Rev B)) 内，用以提供激光驱动器的调制信号。激光波长由频率为2Hz锯齿波在甲烷气体最强吸收峰1653.7nm附近扫描，同时，频率为31.2KHz的正弦波对DFB激光器输出波长进行调制。调制后的激光经过气体吸收池后到达光电探测器(PD，THORLABS PDA20CS)。由采集卡(NI Pxle-6366)将PD信号采集到处理器上。其中，将浓度为990ppm的甲烷标准气体充满气体池，压强20Kpa，温度为28℃，激光完成多次反射后进入到PD探测器中。PD信号即为调制后的透射光强信号，用谐波小波包方法检测二次谐波，并利用二次谐波信号确定检测到的气体浓度。
基于上述配置，本公开实施例的气体浓度检测方法可以将透射光强信号进行快速傅里叶变换，得到其频谱，然后选取频域小波系数，再将频域小波系数进行FFT逆变换，从而得到二次谐波分量，进而得到气体浓度。通过谐波小波包分析实现频域滤波处理，可以有效滤除高频电子噪声。另外，通过该方法可以提高了***自身的抗干扰能力。该方法的检测精度可达到与正交锁相放大算法同等的水平，并且克服了正交锁相放大器滤波器参数设置不准确导致检测失灵。相对传统锁相放大进行检测二次谐波技术，该方法更加的简化。
本领域技术人员可以理解，在具体实施方式的上述方法中，各步骤的撰写顺序并不意味着严格的执行顺序而对实施过程构成任何限定，各步骤的具体执行顺序应当以其功能和可能的内在逻辑确定。
可以理解，本公开提及的上述各个方法实施例，在不违背原理逻辑的情况下，均可以彼此相互结合形成结合后的实施例，限于篇幅，本公开不再赘述。
此外，本公开还提供了气体浓度检测装置、气体浓度检测***，上述均可用来实现本公开提供的任一种气体浓度检测方法，相应技术方案和描述和参见方法部分的相应记载，不再赘述。
图5示出根据本公开实施例的气体浓度检测装置的框图，如图5所示，所述气体浓度检测装置包括：
获取模块10，用于获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；
傅里叶变换模块20，用于对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
谐波小波变换模块30，用于对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频，以及确定所述所需频率对应的子带的频域小波系数，利用所述频域小波系数的傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
确定模块40，用于基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
所述谐波小波变换模块还用于利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
基于确定的每个子带的频率范围，确定所需频率对应的子带。
在一些可能的实施方式中，所述谐波小波变换模块还用于，利用预设方式得到所述所需频率对应的子带的频域小波系数；
对所述频域小波系数执行傅里叶逆变换，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络；
利用所述包络得到透射光强信号中的二次谐波信号。
在一些可能的实施方式中，所述确定模块还用于获取所述入射光的光参数，利用预设方式，基于所述二次谐波信号和所述光参数，得到所述气体浓度。
另外，本公开实施例还提供了一种气体浓度检测***，如图3所示为该***的框图。所述***可以包括：
气体吸收池，用于存储待检测气体；
激光器，用于根据调制信号向所述气体吸收池发射入射光；
探测器，用于检测从气体吸收池透射出的透射光的透射光强信号；
处理器，用于执行如第一方面中任意一项所述的气体浓度检测方法，得到所述待检测气体的气体浓度。
在一些实施例中，本公开实施例提供的装置具有的功能或包含的模块可以用于执行上文方法实施例描述的方法，其具体实现可以参照上文方法实施例的描述，为了简洁，这里不再赘述。
本公开实施例还提出一种计算机可读存储介质，其上存储有计算机程序指令，所述计算机程序指令被处理器执行时实现上述方法。计算机可读存储介质可以是非易失性计算机可读存储介质。
本公开实施例还提出一种电子设备，包括：处理器；用于存储处理器可执行指令的存储器；其中，所述处理器被配置为上述方法。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种气体浓度检测方法，其特征在于，包括：
获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；
对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频；
确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
2.根据权利要求1所述的方法，其特征在于，所述对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，包括：
利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
基于确定的每个子带的频率范围，确定所需频率对应的子带。
3.根据权利要求1所述的方法，其特征在于，所述确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号，包括：
利用第一预设方式得到所述所需频率对应的子带的频域小波系数；
对所述频域小波系数执行傅里叶逆变换，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络；
利用所述包络得到透射光强信号中的二次谐波信号。
4.根据权利要求1-3中任意一项所述的方法，其特征在于，所述基于所述二次谐波信号得到所述待检测浓度的气体的浓度，包括：
获取所述入射光的光参数；
利用第二预设方式，基于所述二次谐波信号和所述光参数，得到所述气体浓度。
5.根据权利要求4所述的方法，其特征在于，所述第二预设方式的表达式为：
其中，An(v0)为透射光强信号的n次谐波信号幅值，C为待检测气体的气体浓度，L为气体的吸收光程长，I0为入射光的光强，
为气体在频率v处吸收系数，V0为气体的中心频率。
6.一种气体浓度检测装置，其特征在于，包括：
获取模块，用于获取通过探测器检测的透射光强信号，所述透射光强信号为入射光经过待检测浓度的气体之后的光强信号；
傅里叶变换模块，用于对所述透射光强信号执行傅里叶变换，得到所述透射光强信号的频谱；
谐波小波变换模块，用于对所述频谱进行谐波小波包频带划分，确定所需频率对应的子带，所需频率为生成所述入射光的调制频率的二倍频，以及确定所述所需频率对应的子带的频域小波系数，对所述频域小波系数进行傅里叶逆变换，利用所述傅里叶逆变换的处理结果得到透射光强信号中的二次谐波信号；
确定模块，用于基于所述二次谐波信号得到所述待检测浓度的气体的浓度。
7.根据权利要求6所述的装置，其特征在于，所述谐波小波变换模块还用于利用谐波小波包的分解层数，以及设定的子带宽度，确定每个子带的频率范围；
基于确定的每个子带的频率范围，确定所需频率对应的子带。
8.根据权利要求6所述的装置，其特征在于，所述谐波小波变换模块还用于，利用预设方式得到所述所需频率对应的子带的频域小波系数；
对所述频域小波系数执行傅里叶逆变换，并利用希尔伯特变换得到所述傅里叶逆变换的处理结果的包络；
利用所述包络得到透射光强信号中的二次谐波信号。
9.根据权利要求6-8中任意一项所述的装置，其特征在于，所述确定模块还用于获取所述入射光的光参数，利用预设方式，基于所述二次谐波信号和所述光参数，得到所述气体浓度。
10.一种气体浓度检测***，其特征在于，包括：
气体吸收池，用于存储待检测气体；
激光器，用于根据调制信号向所述气体吸收池发射入射光；
探测器，用于检测从气体吸收池透射出的透射光的透射光强信号；
处理器，用于执行如权利要求1-5中任意一项所述的气体浓度检测方法，得到所述待检测气体的气体浓度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911172382.5A CN110865042A (zh) | 2019-11-26 | 2019-11-26 | 气体浓度检测方法、装置和*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911172382.5A CN110865042A (zh) | 2019-11-26 | 2019-11-26 | 气体浓度检测方法、装置和*** |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110865042A true CN110865042A (zh) | 2020-03-06 |
Family
ID=69656370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911172382.5A Withdrawn CN110865042A (zh) | 2019-11-26 | 2019-11-26 | 气体浓度检测方法、装置和*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110865042A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113390825A (zh) * | 2021-05-17 | 2021-09-14 | 西安理工大学 | 基于tdlas的时频域联合的气体浓度反演方法及装置 |
| CN113757570A (zh) * | 2021-07-28 | 2021-12-07 | 北京市燃气集团有限责任公司 | 一种天然气管道甲烷泄漏检测装置 |
| CN114705654A (zh) * | 2022-04-01 | 2022-07-05 | 华电智控（北京）技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
| WO2012115150A1 (ja) * | 2011-02-25 | 2012-08-30 | 三菱重工業株式会社 | 信号処理装置およびレーザ計測装置 |
| CN103033481A (zh) * | 2012-12-21 | 2013-04-10 | 重庆川仪自动化股份有限公司 | 基于fft的激光分析仪二次谐波滤波方法 |
| CN103886199A (zh) * | 2014-03-19 | 2014-06-25 | 中国人民解放军装备学院 | 调制光谱信号的谐波小波分析方法 |
| CN104237161A (zh) * | 2014-10-15 | 2014-12-24 | 中国科学院合肥物质科学研究院 | 一种煤炭自燃指标气体多组分实时在线远程监测装置及方法 |
| CN104458650A (zh) * | 2013-09-25 | 2015-03-25 | 西门子公司 | 用于测量试样气体的气体成分浓度的气体分析仪和方法 |
| CN107860726A (zh) * | 2017-10-19 | 2018-03-30 | 江汉大学 | 准分布式气体浓度检测二次谐波归一化方法 |
| CN110110607A (zh) * | 2019-04-12 | 2019-08-09 | 上海交通大学 | 基于谐波小波包提取狭窄带宽信号的方法 |
-
2019
- 2019-11-26 CN CN201911172382.5A patent/CN110865042A/zh not_active Withdrawn
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012115150A1 (ja) * | 2011-02-25 | 2012-08-30 | 三菱重工業株式会社 | 信号処理装置およびレーザ計測装置 |
| CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
| CN103033481A (zh) * | 2012-12-21 | 2013-04-10 | 重庆川仪自动化股份有限公司 | 基于fft的激光分析仪二次谐波滤波方法 |
| CN104458650A (zh) * | 2013-09-25 | 2015-03-25 | 西门子公司 | 用于测量试样气体的气体成分浓度的气体分析仪和方法 |
| CN103886199A (zh) * | 2014-03-19 | 2014-06-25 | 中国人民解放军装备学院 | 调制光谱信号的谐波小波分析方法 |
| CN104237161A (zh) * | 2014-10-15 | 2014-12-24 | 中国科学院合肥物质科学研究院 | 一种煤炭自燃指标气体多组分实时在线远程监测装置及方法 |
| CN107860726A (zh) * | 2017-10-19 | 2018-03-30 | 江汉大学 | 准分布式气体浓度检测二次谐波归一化方法 |
| CN110110607A (zh) * | 2019-04-12 | 2019-08-09 | 上海交通大学 | 基于谐波小波包提取狭窄带宽信号的方法 |
Non-Patent Citations (3)
| Title |
|---|
| 侯铁双 等: "基于谐波小波包变换的信号检测", 《西安邮电大学学报》 * |
| 张瑞峰 等: "可调谐激光遥测甲烷浓度的研究", 《电子测量技术》 * |
| 王吉元: "基于小波分析的虚拟甲烷浓度谐波检测仪", 《测控技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113390825A (zh) * | 2021-05-17 | 2021-09-14 | 西安理工大学 | 基于tdlas的时频域联合的气体浓度反演方法及装置 |
| CN113757570A (zh) * | 2021-07-28 | 2021-12-07 | 北京市燃气集团有限责任公司 | 一种天然气管道甲烷泄漏检测装置 |
| CN114705654A (zh) * | 2022-04-01 | 2022-07-05 | 华电智控（北京）技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
| CN114705654B (zh) * | 2022-04-01 | 2022-11-18 | 华电智控（北京）技术有限公司 | 一种兼容光学干涉效应的气体浓度测量方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110865042A (zh) | 气体浓度检测方法、装置和*** | |
| Zheng et al. | Performance improvement of a near-infrared CH4 detection device using wavelet-denoising-assisted wavelength modulation technique | |
| US8085403B2 (en) | Photoacoustic sensor | |
| CN102768197B (zh) | 用于以减小的压力相关性检测气体浓度的方法和设备 | |
| CN109102455B (zh) | 缺陷检测方法、检测图像生成方法、***及存储设备 | |
| CN101435773A (zh) | 基于准连续二极管激光器调制光谱气体监测方法和装置 | |
| Rabasović et al. | Low-cost, portable photoacoustic setup for solid samples | |
| CN110044824A (zh) | 一种基于石英音叉的双光谱气体检测装置及方法 | |
| CN101482501A (zh) | 温室气体的激光检测***及方法 | |
| CN109632687A (zh) | 四氟化碳气体浓度检测校正方法、装置、设备及存储介质 | |
| CN113390825B (zh) | 基于tdlas的时频域联合的气体浓度反演方法及装置 | |
| CN108267420A (zh) | 种子含水量检测***和方法 | |
| CN102809547A (zh) | 散射增强可调谐二极管激光痕量气体检测装置及检测方法 | |
| CN111521581B (zh) | 一氧化碳和甲烷组分判断及浓度检测方法、装置及应用 | |
| CN105158205A (zh) | 甲烷气团界面识别和可调量程激光遥测甲烷浓度的方法和装置 | |
| CN103115877A (zh) | 一种测量不同气体浓度的光学检测*** | |
| CN107560730A (zh) | 双腔式光声光谱仪 | |
| CN105044584B (zh) | 一种用于检测半导体器件的电荷及电场响应的*** | |
| CN111970051A (zh) | 一种电光相位调制器残余幅度调制系数的测试装置及方法 | |
| CN103411923A (zh) | 一种采用双路可调谐激光吸收光谱技术的归一化降噪方法 | |
| Rui-Feng et al. | Influence of laser intensity in second-harmonic detection with tunable diode laser multi-pass absorption spectroscopy | |
| CN110646384B (zh) | 一种半导体材料电阻率光学测量方法 | |
| Werner et al. | Frequency behaviour of light‐addressable potentiometric sensors | |
| CN114441477B (zh) | 基于宽吸收谱线气体的浓度检测方法及装置 | |
| CN115993346A (zh) | 基于tdlas及温度补偿的大气质量监测方法及*** |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200306 |
|
| WW01 | Invention patent application withdrawn after publication |