CN115436320A - 一种混合气体背景下气体吸收线的修正方法 - Google Patents

一种混合气体背景下气体吸收线的修正方法 Download PDF

Info

Publication number
CN115436320A
CN115436320A CN202211078916.XA CN202211078916A CN115436320A CN 115436320 A CN115436320 A CN 115436320A CN 202211078916 A CN202211078916 A CN 202211078916A CN 115436320 A CN115436320 A CN 115436320A
Authority
CN
China
Prior art keywords
gas
absorption
coefficient
line
absorption line
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
CN202211078916.XA
Other languages
English (en)
Other versions
CN115436320B (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.)
Hunan Wuling Power Technology Co Ltd
Wuling Power Corp Ltd
Original Assignee
Hunan Wuling Power Technology Co Ltd
Wuling Power Corp 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 Hunan Wuling Power Technology Co Ltd, Wuling Power Corp Ltd filed Critical Hunan Wuling Power Technology Co Ltd
Priority to CN202211078916.XA priority Critical patent/CN115436320B/zh
Publication of CN115436320A publication Critical patent/CN115436320A/zh
Application granted granted Critical
Publication of CN115436320B publication Critical patent/CN115436320B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/39Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/39Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
    • G01N2021/396Type of laser source
    • G01N2021/399Diode laser
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/12Circuits of general importance; Signal processing
    • G01N2201/121Correction signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/12Circuits of general importance; Signal processing
    • G01N2201/127Calibration; base line adjustment; drift compensation

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

本申请提供一种混合气体背景下气体吸收线的修正方法,包括:S1:将输出中心频率υn的激光器通过浓度为c1的待测气体,获得待测气体对称的吸收线的二次谐波,S2:将输出中心频率υn的激光器通过包含浓度为c2的待测气体n的混合气体中,获得浓度为c2的待测气体在混合气体背景下的实测二次谐波信号Mn(υ);S3:建立待测气体在混合气体背景下拟合吸收线的二次谐波信号模型
Figure DDA0003832912680000011
S4:根据二次谐波信号模型
Figure DDA0003832912680000014
和实测二次谐波信号Mn(υ)获得混合气体背景分离后的待测气体吸收线
Figure DDA0003832912680000012
S5:根据待测气体吸收线
Figure DDA0003832912680000013
和标准吸收二次谐波An(υ)获得谱线形变参数;S6:根据横向伸缩系数和平移系数完成实时吸收线修正。可有效解决混合气体背景对待测气体吸收谱线干扰问题。

Description

一种混合气体背景下气体吸收线的修正方法
技术领域
本发明属于变压器油中溶解气体检测技术领域,具体涉及一种混合气体背景下气体吸收线的修正方法。
背景技术
相较于现有的基于物理/化学的变压器油中溶解气体检测方法具有稳定性差、需定期更换和校准会增加操作和维护的工作量等不足之处,可调谐二极管激光吸收光谱技术(Tunable Diode Laser Absorption Spectroscopy,TDLAS)具有线宽窄、波长扫描快、可在室温下工作等优点在变压器油中溶解气体检测方面具有广阔应用前景。然而在实际测量系统中,激光器驱动电流与温度的漂移等会造成激光器波长偏移,而且由于变压器油中溶解气体种类多,混合气体中的其他分子之间的碰撞和发射也会影响待测气体吸收线的形状,如不进行吸收线修正,会严重影响气体测量准确性,特别是对于低浓度气体测量准确性,从而影响变压器运行状态分析和早期故障诊断。
现有的解决方案中没有考虑到气体吸收线形变来源于温度、光强幅度调制等多因素的影响耦合,仅通过改变激光器驱动参数或者补偿某一影响量难以达到吸收线整体复原的效果,且针对动态变化的工况条件适应性差。特别地,考虑到变压器油中溶解气体检测气体种类多(包括但不限于一氧化碳、二氧化碳、甲烷、乙烷、乙烯、乙炔),混合气体中的其他分子之间的碰撞和发射也会影响待测气体吸收线的形状。
发明内容
本发明要解决的技术问题是提供一种混合气体背景下气体吸收线的修正方法,以解决现有技术存在的上述至少一个问题。
基于上述目的,本申请中一个或多个实施例提供了一种混合气体背景下气体吸收线的修正方法,内容包括以下步骤:
S1:将输出中心频率υn的激光器通过浓度为c1的待测气体n,获得待测气体关于中心频率υn对称的吸收线的二次谐波,即标准吸收二次谐波An(υ),υ为该激光器的瞬时频率;
S2:将输出中心频率υn的激光器通过包含浓度为c2的待测气体n的混合气体中,获得浓度为c2的待测气体在混合气体背景下的实测二次谐波信号Mn(υ),其中c2<c1
S3:建立待测气体在混合气体背景下拟合吸收线的二次谐波信号模型
Figure BDA0003832912660000011
S4:根据二次谐波信号模型
Figure BDA0003832912660000012
和实测二次谐波信号Mn(υ)获得混合气体背景分离后的待测气体吸收线
Figure BDA0003832912660000013
S5:根据待测气体吸收线
Figure BDA0003832912660000014
和标准吸收二次谐波An(υ)获得谱线形变参数;
S6:根据横向伸缩系数和平移系数完成实时吸收线修正,得到修正后的气体吸收线A′n(v)。
基于本发明的上述技术方案,还可以作出如下改进:
可选的,所述步骤S3中的二次谐波信号模型
Figure BDA0003832912660000021
表示为:
Figure BDA0003832912660000022
其中,N为混合气体中气体种类数,
Figure BDA0003832912660000023
分别为混合气体中第i种气体的吸收系数αi(υ)的一阶、二阶和三阶导数,xi,yi,zi分别为拟合系数。
可选的,第i种气体的吸收系数αi(υ)的表示为αi(υ)=SiΦ(ν)Ng,其中Si为气体的吸收线强,Ng为气体的分子数密度,
Figure BDA0003832912660000024
为高斯线型函数,高斯线型函数表示为
Figure BDA0003832912660000025
其中
Figure BDA0003832912660000026
为高斯线型函数的半高半宽,T为气体绝对温度,m为分子质量,υi为第i种气体在该激光器输出频率范围内的吸收线频率。
可选的,步骤S4中获得混合气体背景分离后的待测气体吸收线
Figure BDA0003832912660000027
包括:建立二次谐波信号模型
Figure BDA0003832912660000028
和实测二次谐波信号Mn(υ)的误差公式
Figure BDA0003832912660000029
Figure BDA00038329126600000210
其中υ1、υ2分别表示输出中心频率υn的激光器瞬时频率的范围,且满足:υ1<υ2;应用最小二乘法求得误差E最小时
Figure BDA00038329126600000211
的参数;根据求得的参数,得到待测气体n的气体吸收线
Figure BDA00038329126600000212
可表示为:
Figure BDA00038329126600000213
可选的,步骤S5中的谱线形变参数包括横向伸缩系数和平移系数,其中,横向伸缩系数通过计算
Figure BDA00038329126600000214
A n(υ)吸收峰峰值到谷值的距离得到;平移系数通过吸收线
Figure BDA00038329126600000215
A n(υ)吸收峰峰值位置的差值得到;
其中,定义气体吸收线
Figure BDA00038329126600000216
的吸收峰峰值位置为p*,平移系数为s1,An(υ)吸收峰峰值位置的位置p,则平移系数s1=p*-p;
定义气体吸收线
Figure BDA00038329126600000217
的吸收峰左波谷位置为
Figure BDA00038329126600000218
右波谷位置为
Figure BDA00038329126600000219
横向伸缩系数为sR,伸缩系数为sL,A n(υ)吸收峰左波谷位置为t1,右波谷位置为t2,则伸缩系数
Figure BDA00038329126600000220
横向伸缩系数
Figure BDA00038329126600000221
可选的,所述的横向伸缩系数和平移系数完成实时吸收线修正包括:根据横向伸缩系数,采用线性插值完成气体吸收线的修正;根据平移系数,横向平移吸收线完成气体吸收线的修正;
其中,根据平移系数s1,若s1>0,则将待测气体吸收线
Figure BDA00038329126600000222
向左平移|s1|个点,若s1<0,则将待测气体吸收线
Figure BDA00038329126600000223
向右平移|s1|个点,其中|s1|表示s1的绝对值,若s1=0,则不操作;
根据伸缩系数sL,若sL>1,则对吸收线
Figure BDA0003832912660000031
进行压缩插值,插值系数为1/|sL|;若sL<1,则进行拉伸插值,插值系数为1/|sL|;若sL=1,则不操作;
根据伸缩系数sR,若sR>1,则对吸收线
Figure BDA0003832912660000032
进行压缩插值,插值系数为1/|sR|;若sR<1,则进行拉伸插值,插值系数为1/|sR|;若sR=1,则不操作。
本发明的有益效果是,本发明提供了一种混合气体背景下气体吸收线的修正方法,解决了应用TDLAS技术实现变压器油中溶解气体检测时谱线形变影响检测精度的问题,从而有效提高对变压器油中溶解气体检测准确性、可靠性。可有效解决混合气体背景对待测气体吸收谱线干扰问题,对实现单个激光器下多种气体检测具有重要意义,从而有效简化基于TDLAS的变压器油中溶解气体检测系统结构,降低系统成本;
还可以修正温度等环境干扰对吸收谱线的影响,有效提高基于TDLAS的变压器油中溶解气体检测系统的测量精度和抗干扰能力。
附图说明
图1为本发明实施例的一种混合气体背景下气体吸收线的修正方法的步骤流程示意图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本公开进一步详细说明。
需要说明的是,除非另外定义,本申请中一个或多个实施例使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本申请中一个或多个实施例中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的组成部分。“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接,而是可以包括电性的连接,不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
本申请中一个或多个实施例的一一种混合气体背景下气体吸收线的修正方法,包括以下步骤:
S1:将输出中心频率υn的激光器通过浓度为c1的待测气体n,获得待测气体关于中心频率υn对称的吸收线的二次谐波,即标准吸收二次谐波An(υ),υ为该激光器的瞬时频率;
S2:将输出中心频率υn的激光器通过包含浓度为c2的待测气体n的混合气体中,获得浓度为c2的待测气体在混合气体背景下的实测二次谐波信号Mn(υ),其中c2<c1
S3:建立待测气体在混合气体背景下拟合吸收线的二次谐波信号模型
Figure BDA0003832912660000033
S4:根据二次谐波信号模型
Figure BDA0003832912660000041
和实测二次谐波信号Mn(υ)获得混合气体背景分离后的待测气体吸收线
Figure BDA0003832912660000042
S5:根据待测气体吸收线
Figure BDA0003832912660000043
和标准吸收二次谐波An(υ)获得谱线形变参数;
S6:根据横向伸缩系数和平移系数完成实时吸收线修正,得到修正后的气体吸收线A′n(v)。
可以理解的是,在本实施例中,解决了应用TDLAS技术实现变压器油中溶解气体检测时谱线形变影响检测精度的问题,从而有效提高对变压器油中溶解气体检测准确性、可靠性。可有效解决混合气体背景对待测气体吸收谱线干扰问题,对实现单个激光器下多种气体检测具有重要意义,从而有效简化基于TDLAS的变压器油中溶解气体检测系统结构,降低系统成本;并能够修正温度等环境干扰对吸收谱线的影响,有效提高基于TDLAS的变压器油中溶解气体检测系统的测量精度和抗干扰能力。
根据本实施例的技术方案,基于得到的修正后的气体吸收线A′n(v)实现待测气体浓度特征值的计算。
在一种可选的实施例方式中,所述步骤S3中的二次谐波信号模型
Figure BDA0003832912660000044
表示为:
Figure BDA0003832912660000045
其中,N为混合气体中气体种类数,
Figure BDA0003832912660000046
分别为混合气体中第i种气体的吸收系数αi(υ)的一阶、二阶和三阶导数,xi,yi,zi分别为拟合系数。
可以理解的是,在本实施例中,建立的二次谐波信号模型
Figure BDA0003832912660000047
中涵盖了多种混合气体的存在的情况,并基于混合气体的吸收系数的一阶、二阶和三阶导数,以及拟合系数进行综合计算,使得结果更加精准可靠。
在一种可选的实施例方式中,第i种气体的吸收系数αi(υ)的表示为αi(υ)=SiΦ(ν)Ng,其中Si为气体的吸收线强,Ng为气体的分子数密度,
Figure BDA0003832912660000048
为高斯线型函数,高斯线型函数表示为
Figure BDA0003832912660000049
其中
Figure BDA00038329126600000410
为高斯线型函数的半高半宽,T为气体绝对温度,m为分子质量,υi为第i种气体在该激光器输出频率范围内的吸收线频率。
可以理解的是,在本实施例中,于现有技术不同的是,本实施例中所采用的技术方案中考虑的对气体吸收线的影响因素更多,基于不同类型的影响因素,采用对应的方式进行综合考虑,并基于此对气体吸收线进行进一步修正。
在一种可选的实施例方式中,步骤S4中获得混合气体背景分离后的待测气体吸收线
Figure BDA00038329126600000411
(υ)包括:建立二次谐波信号模型
Figure BDA00038329126600000412
和实测二次谐波信号Mn(υ)的误差公式
Figure BDA00038329126600000413
Figure BDA0003832912660000051
其中υ1、υ2分别表示输出中心频率υn的激光器瞬时频率的范围,且满足:υ1<υ2;应用最小二乘法求得误差E最小时
Figure BDA0003832912660000052
的参数;根据求得的参数,得到待测气体n的气体吸收线
Figure BDA0003832912660000053
可表示为:
Figure BDA0003832912660000054
Figure BDA0003832912660000055
可以理解的是,在本实施例中,步骤S4中根据二次谐波信号模型
Figure BDA0003832912660000056
和实测二次谐波信号Mn(υ)获得混合气体背景分离后的待测气体吸收线
Figure BDA0003832912660000057
具体的,建立
Figure BDA0003832912660000058
Mn(υ)的误差公式,在误差最小的情况下确定
Figure BDA0003832912660000059
的参数,并据此确定气体吸收线的表达式。
在一种可选的实施例方式中,步骤S5中的谱线形变参数包括横向伸缩系数和平移系数,其中,横向伸缩系数通过计算
Figure BDA00038329126600000510
A n(υ)吸收峰峰值到谷值的距离得到;平移系数通过吸收线
Figure BDA00038329126600000511
A n(υ)吸收峰峰值位置的差值得到;
其中,定义气体吸收线
Figure BDA00038329126600000512
的吸收峰峰值位置为p*,平移系数为s1,An(υ)吸收峰峰值位置的位置p,则平移系数s1=p*-p;
定义气体吸收线
Figure BDA00038329126600000513
的吸收峰左波谷位置为
Figure BDA00038329126600000514
右波谷位置为
Figure BDA00038329126600000515
横向伸缩系数为sR,伸缩系数为sL,A n(υ)吸收峰左波谷位置为t1,右波谷位置为t2,则伸缩系数
Figure BDA00038329126600000516
横向伸缩系数
Figure BDA00038329126600000517
在一种可选的实施例方式中,所述的横向伸缩系数和平移系数完成实时吸收线修正包括:根据横向伸缩系数,采用线性插值完成气体吸收线的修正;根据平移系数,横向平移吸收线完成气体吸收线的修正;
其中,根据平移系数s1,若s1>0,则将待测气体吸收线
Figure BDA00038329126600000518
向左平移|s1|个点,若s1<0,则将待测气体吸收线
Figure BDA00038329126600000519
向右平移|s1|个点,其中|s1|表示s1的绝对值,若s1=0,则不操作;
根据伸缩系数sL,若sL>1,则对吸收线
Figure BDA00038329126600000520
进行压缩插值,插值系数为1/|sL|;若sL<1,则进行拉伸插值,插值系数为1/|sL|;若sL=1,则不操作;
根据伸缩系数sR,若sR>1,则对吸收线
Figure BDA00038329126600000521
进行压缩插值,插值系数为1/|sR|;若sR<1,则进行拉伸插值,插值系数为1/|sR|;若sR=1,则不操作。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的系统。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令系统的制造品,该指令系统实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。

Claims (6)

1.一种混合气体背景下气体吸收线的修正方法,其特征是,其包括以下步骤:
S1:将输出中心频率υn的激光器通过浓度为c1的待测气体n,获得待测气体关于中心频率υn对称的吸收线的二次谐波,即标准吸收二次谐波An(υ),υ为该激光器的瞬时频率;
S2:将输出中心频率υn的激光器通过包含浓度为c2的待测气体n的混合气体中,获得浓度为c2的待测气体在混合气体背景下的实测二次谐波信号Mn(υ),其中c2<c1
S3:建立待测气体在混合气体背景下拟合吸收线的二次谐波信号模型
Figure FDA0003832912650000011
S4:根据二次谐波信号模型
Figure FDA0003832912650000012
和实测二次谐波信号Mn(υ)获得混合气体背景分离后的待测气体吸收线
Figure FDA0003832912650000013
S5:根据待测气体吸收线
Figure FDA0003832912650000014
和标准吸收二次谐波An(υ)获得谱线形变参数;
S6:根据横向伸缩系数和平移系数完成实时吸收线修正,得到修正后的气体吸收线A′n(ν)。
2.如权利要求1所述的一种混合气体背景下气体吸收线的修正方法,其特征是,所述步骤S3中的二次谐波信号模型
Figure FDA0003832912650000015
表示为:
Figure FDA0003832912650000016
其中,N为混合气体中气体种类数,
Figure FDA0003832912650000017
分别为混合气体中第i种气体的吸收系数αi(υ)的一阶、二阶和三阶导数,xi,yi,zi分别为拟合系数。
3.如权利要求2所述的一种混合气体背景下气体吸收线的修正方法,其特征是,第i种气体的吸收系数αi(υ)的表示为αi(υ)=SiΦ(ν)Ng,其中Si为气体的吸收线强,Ng为气体的分子数密度,
Figure FDA0003832912650000018
为高斯线型函数,高斯线型函数表示为
Figure FDA0003832912650000019
Figure FDA00038329126500000110
其中
Figure FDA00038329126500000111
为高斯线型函数的半高半宽,T为气体绝对温度,m为分子质量,υi为第i种气体在该激光器输出频率范围内的吸收线频率。
4.如权利要求3所述的一种混合气体背景下气体吸收线的修正方法,其特征是,步骤S4中获得混合气体背景分离后的待测气体吸收线
Figure FDA00038329126500000112
包括:建立二次谐波信号模型
Figure FDA00038329126500000113
和实测二次谐波信号Mn(υ)的误差公式
Figure FDA00038329126500000114
其中υ1、υ2分别表示输出中心频率υn的激光器瞬时频率的范围,且满足:υ1<υ2;应用最小二乘法求得误差E最小时
Figure FDA00038329126500000115
的参数;根据求得的参数,得到待测气体n的气体吸收线
Figure FDA00038329126500000116
可表示为:
Figure FDA00038329126500000117
5.如权利要求4所述的一种混合气体背景下气体吸收线的修正方法,其特征是,步骤S5中的谱线形变参数包括横向伸缩系数和平移系数,其中,横向伸缩系数通过计算
Figure FDA00038329126500000118
An(υ)吸收峰峰值到谷值的距离得到;平移系数通过吸收线
Figure FDA00038329126500000119
An(υ)吸收峰峰值位置的差值得到;
其中,定义气体吸收线
Figure FDA0003832912650000021
的吸收峰峰值位置为p*,平移系数为s1,An(υ)吸收峰峰值位置的位置p,则平移系数s1=p*-p;
定义气体吸收线
Figure FDA0003832912650000022
的吸收峰左波谷位置为
Figure FDA0003832912650000023
右波谷位置为
Figure FDA0003832912650000024
横向伸缩系数为sR,伸缩系数为sL,An(υ)吸收峰左波谷位置为t1,右波谷位置为t2,则伸缩系数
Figure FDA0003832912650000025
横向伸缩系数
Figure FDA0003832912650000026
6.如权利要求5所述的一种混合气体背景下气体吸收线的修正方法,其特征是,所述的横向伸缩系数和平移系数完成实时吸收线修正包括:根据横向伸缩系数,采用线性插值完成气体吸收线的修正;根据平移系数,横向平移吸收线完成气体吸收线的修正;
其中,根据平移系数s1,若s1>0,则将待测气体吸收线
Figure FDA0003832912650000027
向左平移|s1|个点,若s1<0,则将待测气体吸收线
Figure FDA0003832912650000028
向右平移|s1|个点,其中|s1|表示s1的绝对值,若s1=0,则不操作;
根据伸缩系数sL,若sL>1,则对吸收线
Figure FDA0003832912650000029
进行压缩插值,插值系数为1/|sL|;若sL<1,则进行拉伸插值,插值系数为1/|sL|;若sL=1,则不操作;
根据伸缩系数sR,若sR>1,则对吸收线
Figure FDA00038329126500000210
进行压缩插值,插值系数为1/|sR|;若sR<1,则进行拉伸插值,插值系数为1/|sR|;若sR=1,则不操作。
CN202211078916.XA 2022-09-05 2022-09-05 一种混合气体背景下气体吸收线的修正方法 Active CN115436320B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211078916.XA CN115436320B (zh) 2022-09-05 2022-09-05 一种混合气体背景下气体吸收线的修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211078916.XA CN115436320B (zh) 2022-09-05 2022-09-05 一种混合气体背景下气体吸收线的修正方法

Publications (2)

Publication Number Publication Date
CN115436320A true CN115436320A (zh) 2022-12-06
CN115436320B CN115436320B (zh) 2024-03-15

Family

ID=84246280

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211078916.XA Active CN115436320B (zh) 2022-09-05 2022-09-05 一种混合气体背景下气体吸收线的修正方法

Country Status (1)

Country Link
CN (1) CN115436320B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115753678A (zh) * 2023-01-10 2023-03-07 杭州泽天春来科技有限公司 气体分析方法
CN117347301A (zh) * 2023-09-28 2024-01-05 湖南五凌电力科技有限公司 混合气体检测装置、方法、电子设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103323115A (zh) * 2013-06-24 2013-09-25 清华大学 基于波长调制的气体吸收谱线线宽和线型系数的测量方法
WO2014112502A1 (ja) * 2013-01-16 2014-07-24 横河電機株式会社 レーザガス分析装置
CN107255627A (zh) * 2017-05-15 2017-10-17 东南大学 一种基于级数展开的气体浓度测量方法及其检测装置
CN109696415A (zh) * 2019-01-15 2019-04-30 清华大学 一种基于快速傅里叶变换的气体吸收率在线测量方法
US20200088576A1 (en) * 2016-12-23 2020-03-19 Siemens Aktiengesellschaft Method for Correcting a Wavelength and Tuning Range of a Laser Spectrometer
CN111521581A (zh) * 2020-04-24 2020-08-11 北京航天控制仪器研究所 一氧化碳和甲烷组分判断及浓度检测方法、装置及应用
CN114942235A (zh) * 2022-03-30 2022-08-26 北京遥测技术研究所 一种用于复杂背景环境提取有效吸收信息的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014112502A1 (ja) * 2013-01-16 2014-07-24 横河電機株式会社 レーザガス分析装置
CN103323115A (zh) * 2013-06-24 2013-09-25 清华大学 基于波长调制的气体吸收谱线线宽和线型系数的测量方法
US20200088576A1 (en) * 2016-12-23 2020-03-19 Siemens Aktiengesellschaft Method for Correcting a Wavelength and Tuning Range of a Laser Spectrometer
CN107255627A (zh) * 2017-05-15 2017-10-17 东南大学 一种基于级数展开的气体浓度测量方法及其检测装置
CN109696415A (zh) * 2019-01-15 2019-04-30 清华大学 一种基于快速傅里叶变换的气体吸收率在线测量方法
CN111521581A (zh) * 2020-04-24 2020-08-11 北京航天控制仪器研究所 一氧化碳和甲烷组分判断及浓度检测方法、装置及应用
CN114942235A (zh) * 2022-03-30 2022-08-26 北京遥测技术研究所 一种用于复杂背景环境提取有效吸收信息的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115753678A (zh) * 2023-01-10 2023-03-07 杭州泽天春来科技有限公司 气体分析方法
CN117347301A (zh) * 2023-09-28 2024-01-05 湖南五凌电力科技有限公司 混合气体检测装置、方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN115436320B (zh) 2024-03-15

Similar Documents

Publication Publication Date Title
CN115436320A (zh) 一种混合气体背景下气体吸收线的修正方法
CN105300923B (zh) 一种近红外光谱分析仪在线应用时无测点温度补偿模型修正方法
CN104764869A (zh) 一种基于多维特征量的变压器气体故障诊断和报警方法
JP6981817B2 (ja) 分光分析装置及び分光分析方法
CN107703097B (zh) 利用近红外光谱仪构建快速预测原油性质的模型的方法
JP2013113664A (ja) レーザガス分析装置
CN101226147A (zh) 一种气体浓度定量分析仪
CN103868884B (zh) 一种基于调制系数的气体吸收率在线测量方法
CN111521581B (zh) 一氧化碳和甲烷组分判断及浓度检测方法、装置及应用
CN102252982A (zh) 一种用于激光气体分析仪的波长漂移补偿方法
CN110658159A (zh) 一种波长调制气体浓度反演方法
CN105675522A (zh) 一种原油含水率检测装置
CN107941467B (zh) 直接获取分布反馈半导体激光源电流调制波长响应的方法
CN114993944A (zh) 一种甲烷和一氧化碳的共检测方法、装置及设备
CN117330882B (zh) 一种用于滤波器的自动化测试方法及系统
CN103558182A (zh) 一种用于激光气体在线分析仪确定气体浓度的方法
CN111220571B (zh) 一种基于幅值离散度的二次谐波信号拟合方法及系统
CN104458630B (zh) 一种紫外差分气体分析仪的数据处理方法及系统
CN112014344A (zh) 一种污水在线监测方法
CN109521002B (zh) 一种固体燃料颗粒流的燃料特性测量方法
CN114047161B (zh) 激光器发射波长的自诊断方法
US9970867B2 (en) Method of determining the concentration of a gas component and spectrometer therefor
CN108896511A (zh) 一种智能化修复光谱分析仪谱图形变的方法
CN110567932B (zh) 气体拉曼光谱的处理方法及处理装置
CN106124420A (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
GR01 Patent grant
GR01 Patent grant