CN110658159B - 一种波长调制气体浓度反演方法 - Google Patents
一种波长调制气体浓度反演方法 Download PDFInfo
- Publication number
- CN110658159B CN110658159B CN201810720907.3A CN201810720907A CN110658159B CN 110658159 B CN110658159 B CN 110658159B CN 201810720907 A CN201810720907 A CN 201810720907A CN 110658159 B CN110658159 B CN 110658159B
- Authority
- CN
- China
- Prior art keywords
- gas
- concentration
- signal
- theta
- gas concentration
- 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
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
本发明公开了一种波长调制气体浓度反演方法,包括下述步骤:(1)将相互正交的参考信号与探测信号相乘,并进行低通滤波后获得吸收光谱的相互正交的一组1f解调信号X1f、Y1f;(2)根据1f解调信号X1f、Y1f获得1f信号的幅角θ1f;(3)根据幅角θ1f获得被测气体浓度。本发明利用正交相位解调的1f信号X1f、Y1f,通过测量平面矢量(X1f,Y1f)幅角θ1f反演气体浓度,有效的屏蔽了光强抖动等噪声对气体浓度探测的影响,实现气体浓度的高灵敏度探测。
Description
技术领域
本发明属于可调谐激光气体吸收光谱探测领域,更具体地,涉及一种波长调制气体浓度反演方法。
背景技术
可调谐激光气体吸收光谱技术,利用气体分子对某个特定波长激光的吸收,通过波长调谐的方法获得该气体分子在该波长的吸收光谱测量其浓度及其信息。该技术广泛的应用于大气污染物监测、工业控制、安全生产等多个领域。
波长调制技术是用高频的调制信号控制激光的波长与强度,进而将气体吸收光谱信号的频率移至高频,再通过解调制探测,从而极大的降低低频噪声对吸收光谱信号的影响,提高系统信噪比。
通常使用1f归一化2f解调方法。使用频率为调制频率两倍,与调制信号同相位的参考信号,与探测到的光强信号相乘,再通过低通滤波的方法获得吸收光谱的X2r信号。为屏蔽解调相位误差带,通常使用频率为调制频率两倍,正交相位的两路参考信号,与探测到的光强信号相乘,再通过低通滤波的方法获得吸收光谱的X2f、Y2f信号。考虑到光强变化对结果的影响,通常采用1f信号对1f信号进行归一化处理。使用频率为调制频率,正交相位的两路参考信号,与探测到的光强信号相乘,在通过低通滤波的方法获得吸收光谱的X1f、Y1f信号,从而获得矢量(X1f,Y1f)的模R1f 最终获得单通道的1f归一化2f信号X2f/1f(X2f/1f=X2f/R1f)或正交解调的1f归一化的2f信号R2f/1f 使用获得的X2f/1f或R2f/1f有两种方法反演待测气体浓度。第一,采用标准浓度气体标定X2f/1f或R2f/1f目标气体吸收位置的峰值或峰峰值,建立气体浓度-峰值(或峰峰值)对应关系,再由测量光谱峰值(或峰峰值)反演气体浓度。第二,采用光谱数据库相应光谱吸收线的特征参数,根据比尔朗博气体吸收定律,结合测量环境条件(温度、压力等)与测量条件(激光光强、频率特征参数等)模拟特定浓度的单通道的1f归一化2f信号或正交解调的1f归一化的2f信号,与测量光谱进行拟合,残差最小的模拟单通道的1f归一化2f信号或正交解调的1f归一化的2f信号对应的目标气体浓度就是测量目标气体浓度。
为了获得高精度的测量结果,在传统解调方式的基础上,也发明了多种其他信号解调光谱反演方法。例如,德国西门子公司A.Hangauer发明了利用调制信号的多次谐波同时解调获得1f信号、2f信号、3f信号及4f信号,对他们同时进行拟合反演气体浓度。清华大学丁艳军发明了基于2次谐波与4次谐波的气体吸收谱线线宽和线型系数的测量等方法。
但在多种方法的应用过程中,激光光强抖动等噪声造成的气体浓度测量误差仍然是亟待解决的问题。
发明内容
针对现有技术的缺陷,本发明的目的在于提供一种波长调制气体浓度反演方法,旨在解决现有技术中激光光强抖动噪声导致气体浓度测量误差大的问题。
本发明提供了一种波长调制气体浓度反演方法,包括下述步骤:
(1)将相互正交的参考信号与探测信号相乘,并进行低通滤波后获得吸收光谱的相互正交的一组1f解调信号X1f、Y1f;其中,X1f、Y1f分别为坐标系中X轴和Y轴的1f解调信号;
(2)根据所述1f解调信号X1f、Y1f获得1f信号的幅角θ1f;
(3)根据所述幅角θ1f获得被测气体浓度。
其中,所述参考信号的频率、所述探测信号的频率均与调制信号的频率相同。
更进一步地,步骤(2)中幅角θ1f=arctan(Y1f/X1f),且幅角θ1f的峰峰值与气体吸光光谱积分吸光度成正比。
更进一步地,步骤(3)中通过标准浓度气体标定获得被测气体浓度。
更进一步地,标准浓度气体标定具体包括:
(1)在气体吸收池中依次通入包含待测气体浓度上、下限的一组已知浓度的不同浓度气体,测量幅角θ1f的峰峰值;
(2)建立气体浓度与幅角的对应关系模型;该对应关系模型可以是对应关系表,也可以是多项式拟合关系表达式;
(3)测量待测气体幅角θ1f的峰峰值,并根据待测气体幅角θ1f的峰峰值以及所述对应关系模型获得待测气体浓度。
更进一步地,通过在对应关系表中查表或内插求值或根据多项式拟合关系表达式计算获得待测气体浓度。
通过本发明所构思的以上技术方案,与现有技术相比,本发明中由于1f信号的幅值远高于2f信号,因此基于1f信号幅角的波长调制解调方法在原始信号噪声干扰较大的恶劣环境下仍能获得与光强无关的准确的解调结果。本发明利用正交相位解调的1f信号X1f、Y1f,通过测量平面矢量(X1f,Y1f)幅角θ1f反演气体浓度,有效的屏蔽了光强抖动等噪声对气体浓度探测的影响,实现气体浓度的高灵敏度探测。
附图说明
图1是本发明浓度反演流程示意图;
图2是本发明探测的原始光谱信号示意图;
图3是本发明1f解调信号X1f示意图;
图4是本发明1f解调信号Y1f示意图;
图5是本发明1f幅角信号θ1f示意图;
图6是不同积分吸光度下的θ1f数值计算结果;
图7是数值计算θ1r峰峰值与积分吸光度对应关系。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明具体涉及一种利用正交解调1f信号的幅角反演气体浓度的方法,实现高灵敏度气体浓度探测。本发明尤其适用于光强抖动等噪声明显的气体探测环境,其中1f信号是波长调制吸收光谱领域的通用术语,以调制频率同频率的参考信号解调获得的是1f信号(又称一次谐波),以调制频率2倍频的参考信号解调获得的是2f信号(又称二次谐波)。
本发明目的在于提供一种基于1f信号幅角特征的波长调制气体吸收光谱气体浓度反演方法,屏蔽光强抖动等噪声对探测信号幅值的影响,实现高灵敏度气体浓度探测。
根据本发明,采用与调制信号同频率的两路正交参考信号与探测信号相乘,再通过低通滤波的方法获得吸收光谱的X1f、Y1f信号,求Y1f除X1f结果的反正切,获得矢量(X1f,Y1f)的幅角θ1f(θ1f=arctan(Y1f/X1f))。θ1f是与吸收线型一阶导数相似的信号,θ1r信号的峰峰值与目标气体浓度成正比,从而可以通过标准浓度气体标定的方法由θ1f获得被测气体浓度。
下面结合附图与实施例对本发明作进一步说明。
图1所示的是本发明的浓度反演流程示意图;波长调制吸收光谱通过电流调制实现光强与频率同时调制,光强I0在直流光强i0的基础上调制响应的幅值为i1,频率v在基准频率v0的基础上调制响应的幅值为v1,调制频率相同为ω,存在相位差θ。光强I0,频率v可以写成时间t的函数;I0=i0+i1cos(ωt+θ),v=v0+v1 cos(ωt)。
气体吸光度τ(v)是气体浓度小压力P、光程L、线强S以及线型函数σ(v)的乘积,其中χPLS被称为积分吸光度,是与频率无关的函数,线型函数σ(v)在吸收光谱频率范围内的积分为1,气体吸光度τ(v)可以表达为傅里叶级数的形式,
τ(v0+v1 cos(ωt))=χPLSσ(v0+v1 cos(ωt))=∑Hk(v0,v1)cos(kωt)
其中,
I=I0 exp(-τ)
≈I0(1-τ)=i0+i1 cos(ωt+θ)-i0τ(v0+v1 cos(ωt))-i1 cos(ωt+θ)τ(v0+v1 cos(ωt))’
以与频率调制同频率,同相位的参考信号cosωt对原始探测信号进行解调获得X1f信号(如图3所示):
以与频率调制同频率,正交相位的参考信号sinωt对原始探测信号进行解调获得Y1f信号(如图4所示):
所以1f信号的幅角θ1f(如图5)可以表示为:
由于气体吸光度τ(ν)与积分吸光度、气体浓度都成正比,所以它的傅里叶级数H0、H1、H2......也都与积分吸光度、气体浓度成正比,因此θ1f是积分吸光度、气体浓度的函数。
假设i0:ii=10,频率与光强调制相位差θ=150°,对不同积分吸光度下的θ1f做数值计算结果如图6所示,不同积分吸光度下θ1f曲线在1、2位置的峰峰值有明显的区别。1、2位置的峰峰值与积分吸光度存在明显的单调递增关系(如图7所示)。而积分吸光度与气体浓度成正比。因此在气体压力、测量光程一定的情况下,可以通过θ1f曲线在1、2位置的峰峰值来测量气体浓度。
可以通过测量不同浓度目标待测气体光谱的θ1r信号,获得不同浓度下θ1f信号峰值1与谷值2的差值(峰峰值),建立浓度-峰峰值对应关系,对于未知浓度的目标待测气体通过测量θ1f信号的峰峰值,通过差值的方式获得气体浓度。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种波长调制气体浓度反演方法,其特征在于,包括下述步骤:
(1)将相互正交的参考信号与探测信号相乘,并进行低通滤波后获得吸收光谱的相互正交的一组1f解调信号X1f、Y1f;
(2)根据所述1f解调信号X1f、Y1f获得1f信号的幅角θ1f;幅角θ1f=arctan(Y1f/X1f),且幅角θ1f的峰峰值与气体吸光光谱积分吸光度成正比;
(3)根据所述幅角θ1f获得被测气体浓度。
2.如权利要求1所述的波长调制气体浓度反演方法,其特征在于,所述参考信号的频率、所述探测信号的频率均与调制信号的频率相同。
3.如权利要求1-2任一项所述的波长调制气体浓度反演方法,其特征在于,步骤(3)中通过标准浓度气体标定获得被测气体浓度。
4.如权利要求3所述的波长调制气体浓度反演方法,其特征在于,标准浓度气体标定具体包括:
(1)在气体吸收池中依次通入包含待测气体浓度上、下限的一组已知浓度的不同浓度气体,测量幅角θ1f的峰峰值;
(2)建立气体浓度与幅角的对应关系模型;
(3)测量待测气体幅角θ1f的峰峰值,并根据待测气体幅角θ1f的峰峰值以及所述对应关系模型获得待测气体浓度。
5.如权利要求4所述的波长调制气体浓度反演方法,其特征在于,通过在对应关系表中查表或内插求值或根据多项式拟合关系表达式计算获得待测气体浓度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810720907.3A CN110658159B (zh) | 2018-06-29 | 2018-06-29 | 一种波长调制气体浓度反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810720907.3A CN110658159B (zh) | 2018-06-29 | 2018-06-29 | 一种波长调制气体浓度反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110658159A CN110658159A (zh) | 2020-01-07 |
CN110658159B true CN110658159B (zh) | 2022-02-01 |
Family
ID=69028470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810720907.3A Active CN110658159B (zh) | 2018-06-29 | 2018-06-29 | 一种波长调制气体浓度反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110658159B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111351770A (zh) * | 2020-03-18 | 2020-06-30 | 中国科学院合肥物质科学研究院 | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 |
CN114062313B (zh) * | 2021-11-02 | 2023-07-14 | 中南大学 | 一种基于谐波特征三角形的气体浓度反演模型 |
CN114460038B (zh) * | 2021-12-31 | 2023-09-01 | 南京星空低碳科技中心(有限合伙) | 一种在线监测三氧化硫浓度的装置及方法 |
CN116952850B (zh) * | 2023-09-20 | 2024-01-05 | 安徽农业大学 | 一种基于非分散红外检测的动态减少抖动反演方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
CN102680020A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于波长调制光谱技术的气体参数在线测量方法 |
CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
CN103868885A (zh) * | 2014-03-27 | 2014-06-18 | 清华大学 | 基于复合多次谐波的气体浓度在线测量方法 |
US8896836B1 (en) * | 2012-04-23 | 2014-11-25 | Southwest Sciences Incorporated | Fluid properties measurements using wavelength modulation spectroscopy with first harmonic detection |
CN105372205A (zh) * | 2015-11-25 | 2016-03-02 | 山西大学 | 基于S2f方法的免定标波长调制光谱气体检测方法 |
CN107860726A (zh) * | 2017-10-19 | 2018-03-30 | 江汉大学 | 准分布式气体浓度检测二次谐波归一化方法 |
-
2018
- 2018-06-29 CN CN201810720907.3A patent/CN110658159B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590138A (zh) * | 2012-01-16 | 2012-07-18 | 清华大学 | 一种基于激光吸收光谱的气体浓度在线测量方法 |
US8896836B1 (en) * | 2012-04-23 | 2014-11-25 | Southwest Sciences Incorporated | Fluid properties measurements using wavelength modulation spectroscopy with first harmonic detection |
CN102680020A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于波长调制光谱技术的气体参数在线测量方法 |
CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
CN103868885A (zh) * | 2014-03-27 | 2014-06-18 | 清华大学 | 基于复合多次谐波的气体浓度在线测量方法 |
CN105372205A (zh) * | 2015-11-25 | 2016-03-02 | 山西大学 | 基于S2f方法的免定标波长调制光谱气体检测方法 |
CN107860726A (zh) * | 2017-10-19 | 2018-03-30 | 江汉大学 | 准分布式气体浓度检测二次谐波归一化方法 |
Non-Patent Citations (1)
Title |
---|
Design of Lock-In Amplifier in the Optical Fiber Methane Harmonic Detection;Wang Wen-qing等;《IEEE》;20141231;第367-370页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110658159A (zh) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110658159B (zh) | 一种波长调制气体浓度反演方法 | |
CN109100325B (zh) | 一种基于光谱吸收率二次谐波特征提取的气体浓度测量方法 | |
CN103323115B (zh) | 基于波长调制的气体吸收谱线线宽和线型系数的测量方法 | |
CN107063553B (zh) | 一种利用波长调制光谱测量气体压强和组分浓度的装置和方法 | |
CN107091818B (zh) | 一种多气室复杂组分气体分析系统及方法 | |
CN113447458B (zh) | 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法 | |
CN109470154B (zh) | 一种适用于光谱椭偏仪的薄膜厚度初值测量方法 | |
CN111537470B (zh) | 一种基于数字调制的tdlas气体浓度检测方法 | |
US20180149586A1 (en) | Gas analysis apparatus and gas analysis method | |
CN112255186B (zh) | 一种非均匀流场积分吸光度的计算方法和系统 | |
CN111351770A (zh) | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 | |
CN102519597B (zh) | 一种傅里叶变换光谱仪相位校正切趾方法 | |
CN107941467B (zh) | 直接获取分布反馈半导体激光源电流调制波长响应的方法 | |
CN109520941B (zh) | 在线光谱测量仪器的响应函数校正方法 | |
CN104729712B (zh) | 一种星载大气探测傅里叶变换光谱仪数据预处理方法 | |
CN103558182A (zh) | 一种用于激光气体在线分析仪确定气体浓度的方法 | |
US7751051B2 (en) | Method for cross interference correction for correlation spectroscopy | |
Wu et al. | TDLAS second harmonic demodulation based on Hilbert transform | |
CN115436320B (zh) | 一种混合气体背景下气体吸收线的修正方法 | |
CN111157115A (zh) | 一种水下布里渊散射光谱获取方法及装置 | |
CN113252601B (zh) | 基于波长调制光谱技术的多组分物质检测方法及装置 | |
CN113203468B (zh) | 一种光学干涉相位解调方法和声波测量系统 | |
CN116202975A (zh) | 一种水体参数预测方法、存储介质及终端设备 | |
CN204359684U (zh) | 一种高精度气体浓度检测装置 | |
CN103616566B (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 |