CN111351770A - 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 - Google Patents
一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 Download PDFInfo
- Publication number
- CN111351770A CN111351770A CN202010192673.7A CN202010192673A CN111351770A CN 111351770 A CN111351770 A CN 111351770A CN 202010192673 A CN202010192673 A CN 202010192673A CN 111351770 A CN111351770 A CN 111351770A
- Authority
- CN
- China
- Prior art keywords
- signal
- absorption
- background
- phase angle
- 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.)
- Pending
Links
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 80
- 230000003287 optical effect Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000000862 absorption spectrum Methods 0.000 title claims abstract description 16
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 13
- 238000001914 filtration Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012886 linear function Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 24
- 238000001514 detection method Methods 0.000 abstract description 9
- 230000007547 defect Effects 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 42
- 238000005516 engineering process Methods 0.000 description 5
- 238000004847 absorption spectroscopy Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 238000002168 optical frequency-domain reflectometry Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000001285 laser absorption spectroscopy Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000013598 vector Substances 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/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
Abstract
Description
技术领域
本发明涉及可调谐激光气体吸收光谱探测技术领域,具体涉及一种波长调制吸收光谱吸收光程与气体浓度同步反演方法。
背景技术
随着红外二极管激光技术的飞速发展,可调谐二极管激光吸收光谱技术(TDLAS)在温室气体和大气污染监测、燃烧诊断、工业过程控制等领域得到了广泛的应用。在激光吸收光谱中,吸光度通常与气体浓度和吸收路径长度的乘积成正比。因此,吸收路径长度的测量误差直接决定了反演气体浓度的精度。在单通或多通气池的测量配置中,吸收路径是明确定义的,或者可以通过已知气体浓度的参考测量进行校准。然而,多通气池的吸收路径长度也会因机械不稳定性而变化,导致长期测量的不确定性增加。此外,在许多应用中,如气体在散射介质吸收光谱(GASMAS)中的应用、气体浓度的开放路径遥感等,吸收路径长度也可能是未知的或很难确定的。因此,在测量气体吸收信号时测量吸收路径长度一直是人们关注的问题,相关研究人员在这一领域投入了大量的工作。
意大利的Gianfrani等人使用高精度的迈克尔逊干涉仪测量了吸收路径长度。Du等人使用光频域反射计(OFDR)测量多通气池的吸收路径长度。最近,Lou等人演示了一种同时测量气体吸收光谱和吸收路径长度的方法,方法是在一个多气电池中使用光调频连续波(FMCW)干涉术。这种方法虽然可以达到一定的精度,但是测量范围受到二极管激光器的线宽的限制,例如,对于一个典型的10mhz的线宽,测量范围为10米。在GASMAS技术中,激光束在多孔的散射介质中发生了明显的散射,使得吸收路径长度难以测量。Mei等人使用频域光子迁移(FDPM)技术和FMCW技术测量了散射介质中的气体吸收和路径长度,但这需要额外的光电子器件来进行路径长度测量。在开放路径遥感应用中,吸收路径长度一般采用单独的GPS或测距模块进行测量,但这在很多应用中是不方便的。
发明内容
本发明的目的在于提供一种波长调制吸收光谱吸收光程与气体浓度同步反演方法,该方法能够解决现有技术中存在的不足,实现吸收光程和气体浓度的快速、高灵敏度同步反演。
为实现上述目的,本发明采用了以下技术方案:
一种波长调制吸收光谱吸收光程与气体浓度同步反演方法,该方法包括以下步骤:
(1)获取吸收信号It和背景信号It,0。
(4)根据一次谐波的X通道解调信号和Y通道解调信号,求得吸收信号的一次谐波相位角θ1f。
(6)根据吸收信号的一次谐波相位角和背景信号的一次谐波相位角,求得扣除背景的一次谐波相位角θ1f-0。
(8)对扣除背景的一次谐波相位角θ1f-0进行差分计算,得到扣除背景的一次谐波差分相位角δθ1f-0,并采用以下公式对扣除背景的差分相位角δθ1f-0进行最小二乘法线型拟合,求得待测气体的浓度;
其中,压强表示P,L表示吸收光程,T表示温度,S(T)表示吸收线强,表示激光光强零阶傅里叶系数,表示激光光强一阶偶数傅里叶系数,表示激光光强一阶奇数傅里叶系数,a表示调制深度,v表示频率,线型函数φ为一种已知的函数(由高斯线型与洛伦兹线型卷积得到),对实验测得线型按上式进行最小二乘拟合,可得到上式中的唯一变量气体浓度X值。
进一步的,步骤(1)中所述的“吸收信号It和背景信号It,0”,其采用信号采集模块获取,所述信号采集模块包括可调谐窄带激光器、激光器控制器、分束镜、反射镜、光电探测器一和光电探测器二。
在激光器控制器的控制下,激光器发出波长调制扫描信号,波长调制扫描信号经分束镜分成两束,一束直接被光电探测器一吸收,得到背景信号It,0,另一束穿过待测流场,再经反射镜反射后被光电探测器二吸收,得到吸收信号It。
由以上技术方案可知,本发明利用正交相位解调的一次谐波信号X1f、Y1f,通过测量平面矢量(X1f、Y1f)一次谐波相位角的基线值及其差分信号的线型,同步反演吸收光程与气体浓度,实现吸收光程与气体浓度的高灵敏度探测。因为扣除背景的一次谐波相位角及其差分信号几乎是无背景的,且与解调相位无关,因此,本发明能够快速同步反演吸收光程与气体浓度,实现高灵敏度吸收光程与气体浓度探测,尤其适用于遥测等噪声明显的气体探测环境。
附图说明
图1是本发明的工作原理示意图;
图2是本发明中信号采集模块的工作原理示意图;
图3是本发明中的探测的原始光谱信号示意图;
图4是本发明中一次谐波相位角的示意图;
图5是本发明中扣除背景的一次谐波相位角差分信号测量及拟合示意图。
其中:
1、激光器控制器,2、可调谐窄带激光器,3、分束镜,4、光电探测器一,5、光电探测器二,6、反射镜。
具体实施方式
下面结合附图对本发明做进一步说明:
如图1所示的一种波长调制吸收光谱吸收光程与气体浓度同步反演方法,该方法包括以下步骤:
(1)获取吸收信号It和背景信号It,0。
(2)采用正弦信号sin(2πfmt)分别与吸收信号It、背景信号It,0相乘,并对相乘后的结果分别进行低通滤波(低通滤波截止频率为扫描频率的5倍),得到一次谐波的X通道解调信号X1f及X通道背景解调信号
(8)对扣除背景的一次谐波相位角θ1f-0进行差分计算(即求其导数),得到扣除背景的一次谐波差分相位角δθ1f-0,δθ1f-0如图5中的“测量光谱”所示。
采用以下公式对扣除背景的差分相位角δθ1f-0进行最小二乘法线型拟合,得到如图5中所示“拟合光谱”,由拟合结果求得待测气体的浓度;
其中,压强表示P,L表示吸收光程,T表示温度,S(T)表示吸收线强,表示激光光强零阶傅里叶系数,表示激光光强一阶偶数傅里叶系数,表示激光光强一阶奇数傅里叶系数,a表示调制深度,v表示频率,线型函数φ为一种已知的函数(由高斯线型与洛伦兹线型卷积得到),对实验测得线型按上式进行最小二乘拟合,可得到上式中的唯一变量气体浓度X值。
步骤(7)和步骤(8)中涉及的吸收光程计算公式与气体浓度计算公式中的各变量在具体实施过程中容易获取,且这两个公式减弱了环境影响,提高了测量精度。
进一步的,步骤(1)中所述的“吸收信号It和背景信号It,0”,其采用信号采集模块获取,所述信号采集模块包括可调谐窄带激光器2、激光器控制器1、分束镜3、反射镜6、光电探测器一4和光电探测器二5。不同于现有技术,本发明设置了两个光电探测器(光电探测器一4和光电探测器二5),还设置了一个分束镜3。分束镜3用于对激光分束,得到两束相干光,其中一束用于获取吸收信号,另一束作为参考信号,这两束信号中就包含了气体浓度与吸收光程的信息,本发明利用这些信息进行气体浓度与吸收光程同步反演。
在激光器控制器1的控制下,可调谐窄带激光器2发出波长调制扫描信号,波长调制扫描信号经分束镜分成两束,一束直接被光电探测器一4吸收,得到背景信号It,0,另一束穿过待测流场,再经反射镜6反射后被光电探测器二5吸收,得到吸收信号It。
波长调制扫描信号的扫描频率fs为100Hz,调制频率fm为50kHz。光电探测器一获取的信号为模拟信号,该模拟信号采用数据采集卡采集并保存,模拟信号转换为数字信号后,得到It,0,如图3中的“2背景信号”所示。穿过待测流场的波长调制信号,分子受光的作用,发生受激吸收,产生吸收信号,该过程中的激光光强的变化满足Beer-Lambert定律,激光穿过待测流场后,再经远处反射镜反射,由光电探测器二获得吸收信号,并用数据采集卡采集并保存,得到It,如图3中的“1吸收信号”所示。
本发明用相互正交的参考信号与探测信号相乘,低通滤波后获得吸收光谱的相互正交的1f解调信号X1f、Y1f;根据1f解调信号X1f、Y1f获得一次谐波信号的幅角θ1f;根据扣除背景的一次谐波相位角基线值获得吸收光程;根据扣除背景的一次谐波相位角θ1f-0的差分信号的线型,获得被测气体的浓度。
现有的基于吸收光谱理论的气体传感器,大多需要借助于其他技术进行吸收光程的获取,比如激光测距仪或GPS技术,这将大大增加仪器的体积与成本。目前,也有一些技术可以同时测量气体吸收光谱和吸收光程,比如光调频连续波干涉术,但这种方法的测距范围受二极管激光器的线宽限制,例如对于一个典型的10MHz线宽的激光器,测量范围仅为10米。而本发明提出的波长调制吸收光谱吸收光程与气体浓度同步反演方法,可以同步实现气体浓度测量与吸收光程测量,在实际应用中,不仅可以大大减少仪器的数量和体积,降低测量成本,还可以保证气体传感所必须的吸收光程的测量精度,同时测量距离远大于10米,实现了远距离高精度测量。
目前大多数吸收光谱法气体传感器未考虑相位角信息,本发明使用相位角信息进行气体浓度与吸收光程同步反演是已有专利与文献未被报道的,特别是使用相位角信息进行距离的测量。也可以说,本发明是在已有的吸收光谱法气体传感器技术的基础上,对所获得的信号进行的进一步挖掘,得到了更多有用的信息——吸收光程。并且因为扣除背景的一次谐波相位角及其差分信号几乎是无背景的,且与解调相位无关,因此,本发明能够快速同步反演吸收光程与气体浓度,实现高灵敏度吸收光程与气体浓度探测,尤其适用于遥测等噪声明显的气体探测环境。
实验发现,受大气温度、压力变化引起的空气扰动,以及激光器漂移和数据采集系统的误差影响,在光程700米条件下,采用本发明所述的方法对大气CO2和H2O进行吸收光程与浓度同步反演,吸收光程测量精度为0.01m,气体浓度测量精度为0.1ppm。本发明实现的距离测量精度达到了常用的激光测距仪测量精度,气体浓度测量精度与同类吸收光谱法气体传感器测量精度一致。
以上所述的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案作出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (2)
1.一种波长调制吸收光谱吸收光程与气体浓度同步反演方法,其特征在于:该方法包括以下步骤:
(1)获取吸收信号It和背景信号It,0;
(4)根据一次谐波的X通道解调信号和Y通道解调信号,求得吸收信号的一次谐波相位角θ1f;
(6)根据吸收信号的一次谐波相位角和背景信号的一次谐波相位角,求得扣除背景的一次谐波相位角θ1f-0;
(8)对扣除背景的一次谐波相位角θ1f-0进行差分计算,得到扣除背景的一次谐波差分相位角δθ1f-0,并采用以下公式对扣除背景的差分相位角δθ1f-0进行最小二乘法线型拟合,求得待测气体的浓度;
2.根据权利要求1所述的一种波长调制吸收光谱吸收光程与气体浓度同步反演方法,其特征在于:步骤(1)中所述的“吸收信号It和背景信号It,0”,其采用信号采集模块获取,所述信号采集模块包括可调谐窄带激光器、激光器控制器、分束镜、反射镜、光电探测器一和光电探测器二;
在激光器控制器的控制下,激光器发出波长调制扫描信号,波长调制扫描信号经分束镜分成两束,一束直接被光电探测器一吸收,得到背景信号It,0,另一束穿过待测流场,再经反射镜反射后被光电探测器二吸收,得到吸收信号It。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010192673.7A CN111351770A (zh) | 2020-03-18 | 2020-03-18 | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010192673.7A CN111351770A (zh) | 2020-03-18 | 2020-03-18 | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111351770A true CN111351770A (zh) | 2020-06-30 |
Family
ID=71192970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010192673.7A Pending CN111351770A (zh) | 2020-03-18 | 2020-03-18 | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111351770A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112798550A (zh) * | 2021-04-14 | 2021-05-14 | 四川大学 | 一种宽测量范围的激光吸收光谱燃烧诊断方法 |
CN113281299A (zh) * | 2021-03-26 | 2021-08-20 | 安徽波汇智能科技有限公司 | 一种利用高次谐波的气体检测方法 |
CN113295642A (zh) * | 2021-05-17 | 2021-08-24 | 中国科学院合肥物质科学研究院 | 一种氨气分子吸收线参数的中红外光谱测量系统及方法 |
CN113358160A (zh) * | 2021-07-08 | 2021-09-07 | 天津大学 | 一种大气数据测量方法及系统 |
CN113607686A (zh) * | 2021-08-06 | 2021-11-05 | 广东电网有限责任公司广州供电局 | 一种基于修正tdlas背景信号漂移的气体浓度检测方法、系统及计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
CN107247034A (zh) * | 2017-05-24 | 2017-10-13 | 东南大学 | 一种基于吸收光谱技术的双频率波长调制方法 |
CN108931498A (zh) * | 2018-05-23 | 2018-12-04 | 哈尔滨工业大学 | 多通池内气体吸收光谱和吸收光程同步测量的装置和方法 |
CN110658159A (zh) * | 2018-06-29 | 2020-01-07 | 中国科学院合肥物质科学研究院 | 一种波长调制气体浓度反演方法 |
-
2020
- 2020-03-18 CN CN202010192673.7A patent/CN111351770A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680428A (zh) * | 2012-05-16 | 2012-09-19 | 清华大学 | 一种基于一次谐波信号的气体温度和浓度在线测量方法 |
CN107247034A (zh) * | 2017-05-24 | 2017-10-13 | 东南大学 | 一种基于吸收光谱技术的双频率波长调制方法 |
CN108931498A (zh) * | 2018-05-23 | 2018-12-04 | 哈尔滨工业大学 | 多通池内气体吸收光谱和吸收光程同步测量的装置和方法 |
CN110658159A (zh) * | 2018-06-29 | 2020-01-07 | 中国科学院合肥物质科学研究院 | 一种波长调制气体浓度反演方法 |
Non-Patent Citations (2)
Title |
---|
CHENGUANG YANG 等: "Simultaneous measurement of gas absorption and path length by employing the first harmonic phase angle method in wavelength modulation spectroscopy", 《OPTICS EXPRESS》 * |
CHENGUANG YANG 等: "Wavelength modulation spectroscopy by employing the first harmonic phase angle method", 《OPTICS EXPRESS》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281299A (zh) * | 2021-03-26 | 2021-08-20 | 安徽波汇智能科技有限公司 | 一种利用高次谐波的气体检测方法 |
CN112798550A (zh) * | 2021-04-14 | 2021-05-14 | 四川大学 | 一种宽测量范围的激光吸收光谱燃烧诊断方法 |
CN112798550B (zh) * | 2021-04-14 | 2021-07-13 | 四川大学 | 一种宽测量范围的激光吸收光谱燃烧诊断方法 |
CN113295642A (zh) * | 2021-05-17 | 2021-08-24 | 中国科学院合肥物质科学研究院 | 一种氨气分子吸收线参数的中红外光谱测量系统及方法 |
CN113358160A (zh) * | 2021-07-08 | 2021-09-07 | 天津大学 | 一种大气数据测量方法及系统 |
CN113358160B (zh) * | 2021-07-08 | 2022-05-17 | 天津大学 | 一种大气数据测量方法及系统 |
CN113607686A (zh) * | 2021-08-06 | 2021-11-05 | 广东电网有限责任公司广州供电局 | 一种基于修正tdlas背景信号漂移的气体浓度检测方法、系统及计算机可读存储介质 |
CN113607686B (zh) * | 2021-08-06 | 2023-11-21 | 广东电网有限责任公司广州供电局 | 一种基于修正tdlas背景信号漂移的气体浓度检测方法、系统及计算机可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111351770A (zh) | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 | |
US7407817B2 (en) | Surface plasmon resonance sensors and method for detecting samples using the same | |
CN108981953B (zh) | 一种基于干涉调制原理的激光吸收光谱测温方法和系统 | |
US6130439A (en) | Instrument for measuring the refractive index of a fluid | |
CN109991189B (zh) | 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法 | |
CN108760681A (zh) | 一种基于波形分解的路径平均温度测量系统与方法 | |
CN101532885B (zh) | 温度测量装置和温度测量方法 | |
CN103954589B (zh) | 一种光学材料折射率的精密测量装置及方法 | |
CN111829982A (zh) | 多组分气体浓度检测及温度误差修正方法 | |
Guo et al. | A portable laser-based sensor for detecting H2S in domestic natural gas | |
CN114384045B (zh) | 一种痕量气体浓度和路径长度实时检测系统与方法 | |
CN110207733B (zh) | 基于扫频激光的光纤干涉仪臂长差测量装置及方法 | |
CN103364371A (zh) | 同轴式光热干涉的大气气溶胶吸收系数差分测量新方法 | |
CN112378616A (zh) | 基于波长调制吸收光谱的高速流场多参数测量系统及方法 | |
CN114993987B (zh) | 基于吸收光谱幅值调制的温度及气体浓度测量方法与系统 | |
CN110044847A (zh) | 一种不受光源漂移影响的全内反射式折射率传感方法 | |
CN113447458A (zh) | 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法 | |
CN105865500B (zh) | 一种光纤激光干涉型传感器的移相解调法 | |
Zhang et al. | Water vapor detection system based on scanning spectra | |
CN103278475B (zh) | 透明介质折射率的测量装置及测量方法 | |
CA2870177A1 (en) | High resolution calibration-free infrared spectrometer | |
JP2012173218A (ja) | 干渉計及び測定方法 | |
Yeh et al. | Real-time measurement of glucose concentration using position sensing detector | |
CN106247953B (zh) | 一种同时测量相位和古斯-汉欣位移的方法及装置 | |
CN115343233B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200630 |