CN115753678A - 气体分析方法 - Google Patents

气体分析方法 Download PDF

Info

Publication number
CN115753678A
CN115753678A CN202310031277.XA CN202310031277A CN115753678A CN 115753678 A CN115753678 A CN 115753678A CN 202310031277 A CN202310031277 A CN 202310031277A CN 115753678 A CN115753678 A CN 115753678A
Authority
CN
China
Prior art keywords
gas
spectrogram
waveform
measured
standard
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310031277.XA
Other languages
English (en)
Other versions
CN115753678B (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.)
Hangzhou Zetian Chunlai Technology Co ltd
Original Assignee
Hangzhou Zetian Chunlai Technology Co ltd
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 Hangzhou Zetian Chunlai Technology Co ltd filed Critical Hangzhou Zetian Chunlai Technology Co ltd
Priority to CN202310031277.XA priority Critical patent/CN115753678B/zh
Publication of CN115753678A publication Critical patent/CN115753678A/zh
Application granted granted Critical
Publication of CN115753678B publication Critical patent/CN115753678B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/20Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本申请公开了一种气体分析方法,方法包括:S1、构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系;S2、对第一标准气体进行检测,获得第一标准波形谱图;S3、对第二标准气体进行检测获得第二标准波形谱图,获得信号峰峰值;S4、将待测气体测量波形谱图进行对齐;S5、构建对齐后的待测气体测量波形谱图与第二标准波形谱图之间的函数关系,计算目标气体测量回归系数;S6、根据对齐后的待测气体测量波形谱图获得补偿信号幅值比例因子并推算干扰信号幅值比例因子,根据干扰信号幅值比例因子、信号峰峰值以及目标气体测量回归系数计算待测气体中目标气体浓度。利用本申请提供的方法可以排除背景气体的干扰,提高检测精度。

Description

气体分析方法
技术领域
本申请涉及气体检测技术领域,尤其涉及一种气体分析方法。
背景技术
基于激光吸收光谱技术可以实现NH3、CO2、CO、CH4、H2O、NO、H2S等气体组分测量,通常采用垂直腔面发射激光器或分布式反馈激光器作为可调谐激光光源,将频率发生电路产生的三角波信号和正弦波信号共同叠加在激光器驱动电流上,对半导体激光器输出波长进行扫描和调制。在接收端将调制的激光进行光电转换,并将输出电信号通过锁相放大器的同步检波得到二次谐波信号,通过二次谐波峰高等信息计算被测气体浓度值。
激光吸收光谱技术扫描范围非常小,远小于红外光谱法、紫外光谱法,因此,激光光谱通常未扫描到背景气体,测量不受背景气体干扰。但是在一些背景复杂、特殊工况下,被测气体吸收谱线与背景气体吸收谱线重合或部分重合,导致测量受到背景气体变化的干扰。
基于以上问题,有必要提出一种技术方案以排除背景气体对气体分析的干扰,提高气体测量的准确性。
发明内容
本申请的目的在于提供一种技术方案以排除背景气体对气体分析的干扰,提高气体测量的准确性。
基于以上目的,本申请提供一种气体分析方法,方法包括以下步骤:
S1、对含有不同浓度的背景气体的标准气体进行检测,构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系;
S2、通入第一标准气体,第一标准气体仅含定量浓度的目标气体,对第一标准气体进行检测,获得第一标准波形谱图;
S3、通入第二标准气体,第二标准气体包含定量浓度的目标气体和背景气体,对第二标准气体进行检测获得第二标准波形谱图,获得信号峰峰值;
S4、通入待测气体,待测气体包含目标气体,对待测气体进行检测,获得待测气体测量波形谱图,并将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐,形成对齐后的待测气体测量波形谱图;
S5、构建对齐后的待测气体测量波形谱图与第二标准波形谱图的吸收特征之间的函数关系,第二标准波形谱图的吸收特征的系数为目标气体测量回归系数,基于最小二乘法计算目标气体测量回归系数;
S6、根据对齐后的待测气体测量波形谱图获得补偿信号幅值比例因子,并根据补偿信号幅值比例因子获得干扰信号幅值比例因子,根据干扰信号幅值比例因子、信号峰峰值以及目标气体测量回归系数计算待测气体中目标气体浓度。
进一步的,步骤S1中包括:
S101、分别对L份仅含不同浓度背景气体的第三标准气体进行检测,获得L幅波形谱图,每幅波形谱图均包括表征背景气体补偿信号的第一部分以及表征背景气体干扰信号的第二部分;
S102、选择其中一个背景气体浓度所对应的波形谱图作为第三标准波形谱图,将剩余的L-1幅波形谱图作为背景气体测量谱,利用第三标准波形谱图第一部分表示任一背景气体测量谱的第一部分,并利用第三标准波形谱图的第二部分表示背景气体测量谱的第二部分;
S103、背景气体测量谱的第一部分作为背景气体补偿信号测量谱,背景气体测量谱的第二部分作为背景气体干扰信号测量谱,根据背景气体补偿信号测量谱和背景气体干扰信号测量谱,基于最小二乘法计算任一背景气体测量谱的补偿信号幅值比例因子和干扰信号幅值比例因子;
S104、对S103得到的多组补偿信号幅值比例因子和干扰信号幅值比例因子进行多项式拟合,获得补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系。
进一步的,步骤S102中,
背景气体补偿信号测量谱表示为:
Figure 100002_DEST_PATH_IMAGE001
背景气体干扰信号测量谱
Figure 100002_DEST_PATH_IMAGE002
式中,i表示采样点数,i =1, 2, …, n,j表示采样点数,j =1, 2, …, N,X i X j 分别表示第三标准波形谱图中第i和第j个采样点所对应的吸收特征,Y i '表示背景气体补偿信号测量谱的吸收特征,Z j 表示背景气体干扰信号测量谱的吸收特征;ce表示基线纵向偏移;回归系数d表示补偿信号幅值比例因子,回归系数f表示干扰信号幅值比例因子。
进一步的,补偿信号幅值比例因子表示为:
Figure 100002_DEST_PATH_IMAGE003
干扰信号幅值比例因子表示为:
Figure 100002_DEST_PATH_IMAGE004
式中,X i X j 分别表示第三标准波形谱图中第i和第j个采样点所对应的吸收特征,Y i '表示背景气体补偿信号测量谱的吸收特征,Z j 表示背景气体干扰信号测量谱的吸收特征,ij表示采样点数,i =1, 2, …, n,j =1, 2, …, N。
进一步的,采用四阶多项式拟合,获得背景气体补偿信号值与背景气体干扰信号值之间的函数关系,函数关系表示如下:
Figure 100002_DEST_PATH_IMAGE005
式中,a'b'c'd'e'代表各项拟合系数,d表示补偿信号幅值比例因子,f表示干扰信号幅值比例因子。
进一步的,步骤S4中,将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐包括如下步骤;
S401、分别确定第一标准波形谱图的采样范围起止点和待测气体波形谱图的采样范围起止点,待测气体波形谱图的采样范围起止点间距大于第一标准波形谱图的采样范围起止点间距;
S402、在待测气体测量波形谱图采样范围起止点之间,待测气体测量波形谱图中的采样点开始平移逐一与第一标准波形谱图做线性相关,然后比较求取最大线性相关系数;
S403、根据最大线性相关系数,计算待测气体测量波形谱图的波形最高点位置和偏移量,从而将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐,形成对齐后的待测气体测量波形谱图。
进一步的,最大线性相关系数表示为:
Figure 100002_DEST_PATH_IMAGE006
式中,r(XY)表示最大线性相关系数,X i 表示第一标准波形谱图中第i个采样点所对应的吸收特征值,Y i 表示待测气体测量波形谱图中第i个采样点所对应的吸收特征值,N表示第一标准波形谱图起止点范围内采样点总数。
进一步的,波形最高点位置表示为:
Figure 100002_DEST_PATH_IMAGE007
式中,
Figure 100002_DEST_PATH_IMAGE008
表示线性相关最大时待测气体测量波形谱图的终点,
Figure 100002_DEST_PATH_IMAGE009
表示线性相关最 大时待测气体测量波形谱图的起点。
进一步的,步骤S5中,对齐后的待测气体测量波形谱图与第二标准波形谱图的吸收特征之间的函数关系表示为:
Figure 100002_DEST_PATH_IMAGE010
式中,k表示采样点数,k =1, 2, …, m,a表示基线纵向偏移,b表示目标气体测量回归系数,X k 表示第二标准波形谱图中第k个采样点所对应的吸收特征;Y k 表示对齐后的待测气体测量波形谱图中第k个采样点所对应的吸收特征;
目标气体测量回归系数表示如下:
Figure 100002_DEST_PATH_IMAGE011
式中,k表示采样点数,m表示采样点总数,X k 表示第三标准波形谱图中第k个采样点所对应的吸收特征。
进一步的,待测气体中的目标气体浓度表示为:
Figure 100002_DEST_PATH_IMAGE012
式中,C为目标气体浓度;b为目标气体测量回归系数;V pp 为第二标准波形谱图的信号峰峰值;V pp '为第二标准波形谱图下归一化的信号峰峰值;b 0 为零点系数;f 0 为干扰信号幅值比例因子;K为量程校准系数;P为待测气体压力;L为光程;T表示待测气体温度;S(T)为被测气体吸收线强;g(T,P,C)为线型函数,根据待测气体压力及待测气体温度确定线性函数的系数。
根据以上说明,本申请提供一种气体分析方法,构建了补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系,并将待测气体波形谱图与第二标准波形谱图对比,基于最小二乘法计算目标气体测量回归系数。在对待测气体进行检测时,可以直接通过检测待测气体获得补偿信号幅值比例因子d,并通过检测获得的补偿信号幅值比例因子d推算出相应的干扰信号幅值比例因子f,利用推算得的干扰信号幅值比例因子f和目标气体测量回归系数对待测气体进行浓度计算,以排除背景气体的干扰,大大提高了检测精度。
附图说明
图1为本申请提供的气体分析方法的流程图;
图2为本申请提供的构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系的流程图;
图3为本申请提供的对仅含不同浓度背景气体的第三标准气体进行检测获得的波形谱图;
图4为本申请提供的待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐的流程图;
图5为本申请提供的待测气体检测结果示意图。
具体实施方式
以下将结合附图所示的具体实施方式对本申请进行详细描述,但这些实施方式并不限制本申请,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本申请的保护范围内。
如图1所示,本申请实施例提供一种气体分析方法,用于排除背景气体对气体分析的干扰,提高气体测量的准确性。本申请实施例提供的气体分析方法包括以下步骤:
S1、对含有不同浓度的背景气体的标准气体进行检测,构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系;
S2、通入第一标准气体,第一标准气体仅含定量浓度的目标气体,对第一标准气体进行检测,获得第一标准波形谱图;
S3、通入第二标准气体,第二标准气体包含定量浓度的目标气体和背景气体,对第二标准气体进行检测获得第二标准波形谱图,获得信号峰峰值;
S4、通入待测气体,待测气体包含目标气体,对待测气体进行检测,获得待测气体测量波形谱图,并将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐,形成对齐后的待测气体测量波形谱图;
S5、构建对齐后的待测气体测量波形谱图与第二标准波形谱图的吸收特征之间的函数关系,第二标准波形谱图的吸收特征的系数为目标气体测量回归系数,基于最小二乘法计算目标气体测量回归系数;
S6、根据对齐后的待测气体测量波形谱图获得补偿信号幅值比例因子,并根据补偿信号幅值比例因子获得干扰信号幅值比例因子,根据干扰信号幅值比例因子、信号峰峰值以及目标气体测量回归系数计算待测气体中目标气体浓度。
根据以上说明,本申请构建了不同背景气体浓度下补偿信号幅值比例因子与干扰信号幅值比例因子之间的函数关系,进而在对待测气体进行检测时,可以通过检测获得的补偿信号幅值比例因子推算出干扰信号幅值比例因子。并且,还将待测气体波形谱图与第二标准波形谱图比较,获得目标气体测量回归系数。本申请可以利用干扰信号幅值比例因子和目标气体测量回归系数排除背景气体对气体分析的干扰,从而提高了气体测量的准确性。
结合图2和图3,作为一种可选的实现方式,在步骤S1中,构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系具体可以包括以下步骤:
S101、分别对L份仅含不同浓度背景气体的第三标准气体进行检测,获得L幅波形谱图。每幅波形谱图均包括表征背景气体补偿信号的第一部分以及表征背景气体干扰信号的第二部分。
S102、选择其中一个背景气体浓度所对应的波形谱图作为第三标准波形谱图,将剩余的L-1幅波形谱图作为背景气体测量谱,
背景气体测量谱的第一部分作为背景气体补偿信号测量谱,利用第三标准波形谱图第一部分表示任一背景气体测量谱的第一部分,则可以有如下表示:
Figure DEST_PATH_IMAGE013
(1)
式中,i表示采样点数,i =1, 2, …, n,Y i '表示背景气体补偿信号测量谱的吸收特征,X i 表示第三标准波形谱图中第i个采样点所对应的吸收特征,c表示基线纵向偏移,回归系数d表示补偿信号幅值比例因子。
背景气体测量谱的第二部分作为背景气体干扰信号测量谱,利用第三标准波形谱图的第二部分表示背景气体测量谱的第二部分,则可以有如下表示:
Figure DEST_PATH_IMAGE014
(2)
式中,j表示采样点数,j =1, 2, …, NZ j 表示背景气体干扰信号测量谱的吸收特征,X j 表示第三标准波形谱图中第j个采样点所对应的吸收特征,e表示基线纵向偏移,回归系数f表示干扰信号幅值比例因子。
S103、作为一种可选的实现方式,可以基于最小二乘法计算补偿信号幅值比例因子和干扰信号幅值比例因子。
具体的,基于最小二乘法计算补偿信号幅值比例因子和干扰信号幅值比例因子包括以下步骤:
令绝对误差平方和最小,即:
Figure DEST_PATH_IMAGE015
(3)
根据微积分中的求极值方法,则公式(3)中c、d、e、f应当满足:
Figure DEST_PATH_IMAGE016
(4)
Figure DEST_PATH_IMAGE017
(5)
式中,n表示第三标准波形谱图第一部分中的总采样点数,N表示第三标准波形谱图第二部分中的总采样点数。
对公式(4)和(5)进行求解,则补偿信号幅值比例因子可以表示为:
Figure DEST_PATH_IMAGE018
(6)
式中,d表示补偿信号幅值比例因子,i表示采样点数,i =1, 2, …, nn表示第三标准波形谱图第一部分中的总采样点数,X i 表示第三标准波形谱图中第i个采样点所对应的吸收特征,Y i '表示背景气体补偿信号测量谱的吸收特征。
干扰信号幅值比例因子可以表示为:
Figure DEST_PATH_IMAGE019
(7)
式中,f表示干扰信号幅值比例因子,N表示第三标准波形谱图第二部分中的总采样点数,j表示采样点数,j =1, 2, …, NZ j 表示背景气体干扰信号测量谱的吸收特征,X j 表示第三标准波形谱图中第j个采样点所对应的吸收特征。
根据公式(6)和(7)可以得到每一浓度背景气体所对应的补偿信号幅值比例因子d以及干扰信号幅值比例因子f。通过对仅含不同浓度背景气体的第三标准气体进行检测,可以得到多组补偿信号幅值比例因子d以及干扰信号幅值比例因子f
S104、对获得的多组补偿信号幅值比例因子d和干扰信号幅值比例因子f进行多项式拟合,可以得到补偿信号幅值比例因子d以及干扰信号幅值比例因子f之间的函数关系。
作为一种可选的实现方式,本申请实施例采用四阶多项式拟合,则补偿信号幅值比例因子d与干扰信号幅值比例因子f之间的函数关系可以表示为:
Figure DEST_PATH_IMAGE020
(8)
式中,a'b'c'd'e'代表各项拟合系数,d表示补偿信号幅值比例因子,f表示干扰信号幅值比例因子。
根据以上说明,本申请对含不同浓度背景气体的第三标准气体进行检测,以获得每一背景气体浓度条件下的补偿信号幅值比例因子d和干扰信号幅值比例因子f,并对获得的多组补偿信号幅值比例因子和干扰信号幅值比例因子进行拟合,从而可以得到补偿信号幅值比例因子d与干扰信号幅值比例因子f之间的函数关系。
在对待测气体进行测量时,可以直接通过检测待测气体获得补偿信号幅值比例因子d,并通过检测获得的补偿信号幅值比例因子d推算出相应的干扰信号幅值比例因子f,利用推算得的干扰信号幅值比例因子f对待测气体进行浓度计算,可以排除背景气体对气体分析的干扰,从而提高了气体测量的准确性。
如图4所示,作为一种可选的实现方式,在步骤S4中,将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐,从而可以去除测量时吸收峰点位漂移的问题,提高测量精度。
具体的,在步骤S4中包括以下步骤:
S401、分别确定第一标准波形谱图的采样范围起止点和待测气体波形谱图的采样范围起止点,待测气体波形谱图的采样范围起止点间距大于第一标准波形谱图的采样范围起止点间距;
S402、在待测气体测量波形谱图采样范围起止点之间,待测气体测量波形谱图中的采样点开始平移逐一与第一标准波形谱图做线性相关,然后比较求取最大线性相关系数;
S403、根据最大线性相关系数,计算波形最高点位置和偏移量,从而将待测气体测量波形谱图的波形最高点位置与第一标准波形谱图对齐,形成对齐后的待测气体测量波形谱图。
其中,在步骤S401中,可以根据实际需求对采样范围起止点进行设置,从而在测量谱范围内找出与标准谱最匹配的一段,一般的,测量谱的起止点范围大于标准谱起止点范围。
作为一种可选的实现方式,在本申请实施例中,可以选择第一标准波形谱图中75至225范围内的采样点,并且,可以选择待测气体波形谱图中20至280范围内的采样点。
则,在待测气体测量波形谱图采样范围起止点之间,将待测气体测量波形谱图中的采样点开始平移逐一与第一标准波形谱图做线性相关,表示如下:
X75→Y20 X75→Y21…… X75→Y130
X76→Y21 X76→Y22…… X76→Y131
X77→Y22 X77→Y23…… X77→Y132
……………………
X225→Y170 X225→Y171…… X225→Y280
最大线性相关系数可以表示为:
Figure DEST_PATH_IMAGE021
(9)
式中,r(XY)表示最大线性相关系数,X i 表示所述第一标准波形谱图中第i个采样点所对应的吸收特征值,Y i 表示所述待测气体测量波形谱图中第i个采样点所对应的吸收特征值,N表示第一标准波形谱图起止点范围内采样点总数。例如,在本申请实施例中,第一标准波形谱图起止点定为75至225,则此时第一标准波形谱图起止点范围内采样点总数N=225-75+1=151。
根据最大线性相关系数,可以获得待测气体测量波形谱图的波形最高点位置X0,波形最高点位置可以表示为:
Figure DEST_PATH_IMAGE022
(10)
式中,
Figure 149183DEST_PATH_IMAGE008
表示线性相关最大时所述待测气体测量波形谱图的终点,
Figure 597613DEST_PATH_IMAGE009
表示线性相 关最大时所述待测气体测量波形谱图的起点。
则根据波形最高点位置可以求得待测气体测量波形谱图与第一标准波形谱图相比的偏移量,其偏移量可以表示为:
偏移量=X0-标准谱最高点位置
式中,X0表示待测气体测量波形谱图的波形最高点位置。
根据以上说明,在获得偏移量后可以将待测气体测量波形谱图与第一标准波形谱图对齐,形成对齐后的待测气体测量波形谱图,以去除测量时吸收峰点位漂移的问题,提高测量精度。
作为一种可选的实现方式,在步骤S5中,对齐后的待测气体测量波形谱图与第二标准波形谱图的吸收特征之间的函数关系表示为:
Figure DEST_PATH_IMAGE023
(11)
式中,k表示采样点数,k =1, 2, …, m,a表示基线纵向偏移,b表示目标气体测量回归系数,X k 表示第二标准波形谱图中第k个采样点所对应的吸收特征;Y k 表示对齐后的待测气体测量波形谱图中第k个采样点所对应的吸收特征;
根据最小二乘法,令绝对误差平方和最小,即:
Figure DEST_PATH_IMAGE024
(12)
式中,k表示采样点数,k =1, 2, …, m,a表示基线纵向偏移,b表示目标气体测量回归系数,X k 表示第二标准波形谱图中第k个采样点所对应的吸收特征;Y k 表示对齐后的待测气体测量波形谱图中第k个采样点所对应的吸收特征。
根据微积分中求极值的方法,a、b应当满足如下条件:
Figure DEST_PATH_IMAGE025
(13)
则公式(13)可以改写为:
Figure DEST_PATH_IMAGE026
(14)
对公式(14)进行求解,可得目标气体测量回归系数,目标气体测量回归系数表示如下:
Figure DEST_PATH_IMAGE027
(15)
式中,b表示目标气体测量回归系数,k表示采样点数,m表示采样点总数,X k 表示第三标准波形谱图中第k个采样点所对应的吸收特征,Y k 表示对齐后的待测气体测量波形谱图中第k个采样点所对应的吸收特征。
作为一种可选的实现方式,从第二标准波形谱图中读取的信号峰峰值包括目标气体的信号峰峰值和背景气体的信号峰峰值,则可以根据公式(16)计算待测气体中的目标气体浓度。待测气体中的目标气体浓度表示为:
Figure 195866DEST_PATH_IMAGE012
(16)
式中,C为目标气体浓度;b为目标气体测量回归系数;V pp 为第二标准波形谱图归一化下的目标气体的信号峰峰值;V pp '为第二标准波形谱图下归一化的背景气体的信号峰峰值;b 0 为零点系数;f 0 为干扰信号幅值比例因子; P为待测气体压力;L为光程;T表示待测气体温度;S(T)为被测气体吸收线强;g(T,P,C)为线型函数,根据待测气体压力及待测气体温度确定线性函数的系数。K为量程校准系数,可以是当通入第二标准气体时,根据已知的目标气体浓度,并结合公式(16)计算获得该量程校准系数K
如图5所示,本申请以含目标气体浓度为1000ppb的混合气作为检测对象进行测试,可以看到在使用本申请提供的算法排除背景气体干扰前测量值在1012ppb上下,且波动起伏较大,与1000ppb的测量真值相比偏差较大。而在使用本申请提供的气体分析方法后,测量值围绕1000ppb上下波动,波动起伏较小,可以认为使用本申请提供的气体分析方法测得的目标气体浓度与测量真值基本一致,大大提高了检测精度。
根据以上说明,本申请提供一种气体分析方法,构建了补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系,在对待测气体进行检测时,可以直接通过检测待测气体获得补偿信号幅值比例因子d,并通过检测获得的补偿信号幅值比例因子d推算出相应的干扰信号幅值比例因子f,并将待测气体波形谱图与第二标准波形谱图对比,基于最小二乘法计算目标气体测量回归系数,利用推算得的干扰信号幅值比例因子f和目标气体测量回归系数对待测气体进行浓度计算,以排除背景气体的干扰,大大提高了检测精度。
以上所揭露的仅为本申请的较佳实施例而已,然其并非用以限定本申请之权利范围,本领域普通技术人员可以理解:在不脱离本申请及所附的权利要求的精神和范围内,改变、修饰、替代、组合、简化,均应为等效的置换方式,仍属于申请所涵盖的范围。

Claims (10)

1.一种气体分析方法,其特征在于,所述方法包括以下步骤:
S1、对含有不同浓度的背景气体的标准气体进行检测,构建补偿信号幅值比例因子和干扰信号幅值比例因子之间的函数关系;
S2、通入第一标准气体,所述第一标准气体仅含定量浓度的目标气体,对所述第一标准气体进行检测,获得第一标准波形谱图;
S3、通入第二标准气体,所述第二标准气体包含定量浓度的目标气体和背景气体,对所述第二标准气体进行检测获得第二标准波形谱图,获得信号峰峰值;
S4、通入待测气体,所述待测气体包含所述目标气体,对所述待测气体进行检测,获得待测气体测量波形谱图,并将所述待测气体测量波形谱图的波形最高点位置与所述第一标准波形谱图对齐,形成对齐后的待测气体波形谱图;
S5、构建所述对齐后的待测气体测量波形谱图与所述第二标准波形谱图的吸收特征之间的函数关系,所述第二标准波形谱图的吸收特征的系数为目标气体测量回归系数,基于最小二乘法计算所述目标气体测量回归系数;
S6、根据所述对齐后的待测气体测量波形谱图获得补偿信号幅值比例因子,并根据所述补偿信号幅值比例因子获得所述干扰信号幅值比例因子,根据所述干扰信号幅值比例因子、所述信号峰峰值以及所述目标气体测量回归系数计算所述待测气体中目标气体浓度。
2.根据权利要求1所述的气体分析方法,其特征在于,
所述步骤S1中包括:
S101、分别对L份仅含不同浓度背景气体的第三标准气体进行检测,获得L幅波形谱图,每幅波形谱图均包括表征背景气体补偿信号的第一部分以及表征背景气体干扰信号的第二部分;
S102、选择其中一个背景气体浓度所对应的波形谱图作为第三标准波形谱图,将剩余的L-1幅波形谱图作为背景气体测量谱,利用第三标准波形谱图第一部分表示任一背景气体测量谱的第一部分,并利用第三标准波形谱图的第二部分表示背景气体测量谱的第二部分;
S103、所述背景气体测量谱的第一部分作为背景气体补偿信号测量谱,所述背景气体测量谱的第二部分作为背景气体干扰信号测量谱,根据所述背景气体补偿信号测量谱和所述背景气体干扰信号测量谱,基于最小二乘法计算任一背景气体测量谱的所述补偿信号幅值比例因子和所述干扰信号幅值比例因子;
S104、对S103得到的多组所述补偿信号幅值比例因子和所述干扰信号幅值比例因子进行多项式拟合,获得所述补偿信号幅值比例因子和所述干扰信号幅值比例因子之间的函数关系。
3.根据权利要求2所述的气体分析方法,其特征在于,所述步骤S102中,
所述背景气体补偿信号测量谱表示为:
Figure DEST_PATH_IMAGE001
所述背景气体干扰信号测量谱
Figure DEST_PATH_IMAGE002
式中,i表示采样点数,i =1, 2, …, n,j表示采样点数,j =1, 2, …, N,X i X j 分别表示第三标准波形谱图中第i和第j个采样点所对应的吸收特征,Y i '表示背景气体补偿信号测量谱的吸收特征,Z j 表示背景气体干扰信号测量谱的吸收特征;ce表示基线纵向偏移;回归系数d表示补偿信号幅值比例因子,回归系数f表示干扰信号幅值比例因子。
4.根据权利要求3所述的气体分析方法,其特征在于,
所述补偿信号幅值比例因子表示为:
Figure DEST_PATH_IMAGE003
所述干扰信号幅值比例因子表示为:
Figure DEST_PATH_IMAGE004
式中,X i X j 分别表示第三标准波形谱图中第i和第j个采样点所对应的吸收特征,Y i '表示背景气体补偿信号测量谱的吸收特征,Z j 表示背景气体干扰信号测量谱的吸收特征,ij表示采样点数,i =1, 2, …, n,j =1, 2, …, N。
5.根据权利要求4所述的气体分析方法,其特征在于,
采用四阶多项式拟合,获得所述背景气体补偿信号值与所述背景气体干扰信号值之间的函数关系,所述函数关系表示如下:
Figure DEST_PATH_IMAGE005
式中,a'b'c'd'e'代表各项拟合系数,d表示补偿信号幅值比例因子,f表示干扰信号幅值比例因子。
6.根据权利要求1所述的气体分析方法,其特征在于,
所述步骤S4中,将所述待测气体测量波形谱图的波形最高点位置与所述第一标准波形谱图对齐包括如下步骤;
S401、分别确定所述第一标准波形谱图的采样范围起止点和所述待测气体波形谱图的采样范围起止点,所述待测气体波形谱图的采样范围起止点间距大于所述第一标准波形谱图的采样范围起止点间距;
S402、在所述待测气体测量波形谱图采样范围起止点之间,待测气体测量波形谱图中的采样点开始平移逐一与第一标准波形谱图做线性相关,然后比较求取最大线性相关系数;
S403、根据所述最大线性相关系数,计算待测气体测量波形谱图的波形最高点位置和偏移量,从而将所述待测气体测量波形谱图的波形最高点位置与所述第一标准波形谱图对齐,形成对齐后的待测气体波形谱图。
7.根据权利要求6所述的气体分析方法,其特征在于,
所述最大线性相关系数表示为:
Figure DEST_PATH_IMAGE006
式中,r(XY)表示最大线性相关系数,X i 表示所述第一标准波形谱图中第i个采样点所对应的吸收特征值,Y i 表示所述待测气体测量波形谱图中第i个采样点所对应的吸收特征值,N表示第一标准波形谱图起止点范围内采样点总数。
8.根据权利要求7所述的气体分析方法,其特征在于,
所述波形最高点位置表示为:
Figure DEST_PATH_IMAGE007
式中,
Figure DEST_PATH_IMAGE008
表示线性相关最大时所述待测气体测量波形谱图的终点,
Figure DEST_PATH_IMAGE009
表示线性相关最 大时所述待测气体测量波形谱图的起点。
9.根据权利要求1所述的气体分析方法,其特征在于,
所述步骤S5中,所述对齐后的待测气体测量波形谱图与所述第二标准波形谱图的吸收特征之间的函数关系表示为:
Figure DEST_PATH_IMAGE010
式中,k表示采样点数,k =1, 2, …, m,a表示基线纵向偏移,b表示目标气体测量回归系数,X k 表示第二标准波形谱图中第k个采样点所对应的吸收特征;Y k 表示对齐后的待测气体测量波形谱图中第k个采样点所对应的吸收特征;
所述目标气体测量回归系数表示如下:
Figure DEST_PATH_IMAGE011
式中,k表示采样点数,m表示采样点总数,X k 表示第三标准波形谱图中第k个采样点所对应的吸收特征。
10.根据权利要求1所述的气体分析方法,其特征在于,
所述待测气体中的目标气体浓度表示为:
Figure DEST_PATH_IMAGE012
式中,C为目标气体浓度;b为目标气体测量回归系数;V pp 为第二标准波形谱图下归一化的目标气体的信号峰峰值;V pp '为第二标准波形谱图下归一化的背景气体的信号峰峰值;b 0 为零点系数;f 0 为干扰信号幅值比例因子;K为量程校准系数;P为待测气体压力;L为光程;T表示待测气体温度;S(T)为被测气体吸收线强;g(T,P,C)为线型函数,根据所述待测气体压力及待测气体温度确定所述线性函数的系数。
CN202310031277.XA 2023-01-10 2023-01-10 气体分析方法 Active CN115753678B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310031277.XA CN115753678B (zh) 2023-01-10 2023-01-10 气体分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310031277.XA CN115753678B (zh) 2023-01-10 2023-01-10 气体分析方法

Publications (2)

Publication Number Publication Date
CN115753678A true CN115753678A (zh) 2023-03-07
CN115753678B CN115753678B (zh) 2023-05-05

Family

ID=85348852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310031277.XA Active CN115753678B (zh) 2023-01-10 2023-01-10 气体分析方法

Country Status (1)

Country Link
CN (1) CN115753678B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115824995A (zh) * 2023-02-22 2023-03-21 天津市极光创新智能科技有限公司 一种用于气体分析的红外激光漫反射监测方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030022039A1 (en) * 2001-07-26 2003-01-30 Lloyd Greg A. Method of compensating a MOS gas sensor, method of manufacturing a MOS gas sensor, MOS gas sensor, and fuel cell system
US20080179530A1 (en) * 2007-01-30 2008-07-31 Xiang Liu Background Compensation By Multiple-Peak Measurements For Absorption Spectroscopy-Based Gas Sensing
US20140067282A1 (en) * 2012-09-03 2014-03-06 Sick Ag Method for the laser spectroscopy of gases
CN107110776A (zh) * 2014-11-11 2017-08-29 光谱传感器公司 具有复杂背景成分的样品气体中的目标分析物检测和量化
CN114813633A (zh) * 2022-04-12 2022-07-29 中国科学院合肥物质科学研究院 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法
CN114993987A (zh) * 2022-05-16 2022-09-02 北京航空航天大学 基于吸收光谱幅值调制的温度及气体浓度测量方法与系统
CN115372285A (zh) * 2022-10-26 2022-11-22 杭州泽天春来科技有限公司 尾气分析装置及方法
CN115436320A (zh) * 2022-09-05 2022-12-06 湖南五凌电力科技有限公司 一种混合气体背景下气体吸收线的修正方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030022039A1 (en) * 2001-07-26 2003-01-30 Lloyd Greg A. Method of compensating a MOS gas sensor, method of manufacturing a MOS gas sensor, MOS gas sensor, and fuel cell system
US20080179530A1 (en) * 2007-01-30 2008-07-31 Xiang Liu Background Compensation By Multiple-Peak Measurements For Absorption Spectroscopy-Based Gas Sensing
US20140067282A1 (en) * 2012-09-03 2014-03-06 Sick Ag Method for the laser spectroscopy of gases
CN107110776A (zh) * 2014-11-11 2017-08-29 光谱传感器公司 具有复杂背景成分的样品气体中的目标分析物检测和量化
CN114813633A (zh) * 2022-04-12 2022-07-29 中国科学院合肥物质科学研究院 一种激光多气体混叠光谱解调和非线性补偿的浓度检测方法
CN114993987A (zh) * 2022-05-16 2022-09-02 北京航空航天大学 基于吸收光谱幅值调制的温度及气体浓度测量方法与系统
CN115436320A (zh) * 2022-09-05 2022-12-06 湖南五凌电力科技有限公司 一种混合气体背景下气体吸收线的修正方法
CN115372285A (zh) * 2022-10-26 2022-11-22 杭州泽天春来科技有限公司 尾气分析装置及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
任双赞等: "一种基于双气室切换分时扫描的新型气体吸收光谱补偿方法", 《光谱学与光谱分析》 *
矫晓敏;尚丽平;: "应用TDLAS探测气体的误差分析与补偿研究" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115824995A (zh) * 2023-02-22 2023-03-21 天津市极光创新智能科技有限公司 一种用于气体分析的红外激光漫反射监测方法及系统

Also Published As

Publication number Publication date
CN115753678B (zh) 2023-05-05

Similar Documents

Publication Publication Date Title
JP5983779B2 (ja) ガス吸収分光装置及びガス吸収分光方法
CN101441165B (zh) 光谱仪的波长校准方法
CN103175822B (zh) 消除拉曼光谱仪台间差的方法
Kessler Jr et al. The deuteron binding energy and the neutron mass
US7511265B2 (en) Method and apparatus for accurate calibration of a reflectometer by using a relative reflectance measurement
US7282703B2 (en) Method and apparatus for accurate calibration of a reflectometer by using a relative reflectance measurement
CN105466560A (zh) 分光分析装置、以及分光分析装置的校正方法
JP6981817B2 (ja) 分光分析装置及び分光分析方法
CN115753678A (zh) 气体分析方法
US6615142B1 (en) Filtering to measure gas concentrations from spectral features
CN110927100B (zh) 用于测量气体通量的系统和测量气体通量的方法
US11193917B2 (en) Composition analysis apparatus and composition analysis method
RU2571185C2 (ru) Способ компенсации дрейфа амплитуды в спектрометре и спектрометр, осуществляющий указанный способ
CN113607686A (zh) 一种基于修正tdlas背景信号漂移的气体浓度检测方法、系统及计算机可读存储介质
CN112729544B (zh) 一种扫描吸收光谱的激光波长锁定系统与方法
CN115265785A (zh) 一种基于多项式拟合校正拉曼光谱漂移的方法
JP7440866B2 (ja) レーザ式ガス分析装置
CN114878515A (zh) 一种可调谐二极管激光吸收信号时间-频率转换方法
JP4146697B2 (ja) 温度計測方法および温度計測装置
CN113567602A (zh) 检测方法以及检测设备
CN111579526A (zh) 一种表征近红外仪器差异和校正的方法
KR102432166B1 (ko) 실리콘 단결정의 탄소 농도 측정 방법 및 장치
CN115684081B (zh) 激光气体分析系统
CN110865067B (zh) 使用对乙酰氨基酚快速校准拉曼光谱仪的方法
Swann et al. Accuracy limits for simple molecular absorption based wavelength references

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