留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于近红外吸收光谱技术的高精度CO2检测系统的研制

李恒宽 朴亨 王鹏 姜炎坤 李峥 陈晨 曲娜 白晖峰 王彪 李美萱

李恒宽, 朴亨, 王鹏, 姜炎坤, 李峥, 陈晨, 曲娜, 白晖峰, 王彪, 李美萱. 基于近红外吸收光谱技术的高精度CO2检测系统的研制[J]. 红外与激光工程, 2023, 52(3): 20210828. doi: 10.3788/IRLA20210828
引用本文: 李恒宽, 朴亨, 王鹏, 姜炎坤, 李峥, 陈晨, 曲娜, 白晖峰, 王彪, 李美萱. 基于近红外吸收光谱技术的高精度CO2检测系统的研制[J]. 红外与激光工程, 2023, 52(3): 20210828. doi: 10.3788/IRLA20210828
Li Hengkuan, Piao Heng, Wang Peng, Jiang Yankun, Li Zheng, Chen Chen, Qu Na, Bai Huifeng, Wang Biao, Li Meixuan. Development of high precision CO2 detection system based on near infrared absorption spectroscopy[J]. Infrared and Laser Engineering, 2023, 52(3): 20210828. doi: 10.3788/IRLA20210828
Citation: Li Hengkuan, Piao Heng, Wang Peng, Jiang Yankun, Li Zheng, Chen Chen, Qu Na, Bai Huifeng, Wang Biao, Li Meixuan. Development of high precision CO2 detection system based on near infrared absorption spectroscopy[J]. Infrared and Laser Engineering, 2023, 52(3): 20210828. doi: 10.3788/IRLA20210828

基于近红外吸收光谱技术的高精度CO2检测系统的研制

doi: 10.3788/IRLA20210828
基金项目: 国家重点研发计划 (2018YFC1503802);吉林省科技发展计划项目(20200201050JC, 20220203016SF);吉林省发改委产业技术研究与开发项目 (2022C045-5);长春市科技发展计划项目 (21ZGG14)
详细信息
    作者简介:

    李恒宽,男,硕士生,主要从事红外光谱检测技术温控系统研究

    通讯作者: 李美萱,女,讲师,博士,主要从事多光谱关联成像技术、光学系统设计与制造方面的研究
  • 中图分类号: TN219

Development of high precision CO2 detection system based on near infrared absorption spectroscopy

Funds: National Key Research and Development Program of China (2018YFC1503802); Science and Technology Department of Jilin Province of China (20200201050JC, 20220203016SF); Industrial Technology Research and Development Project of Jilin Provincial Development and Reform Commission (2022C045-5); Science and Technology Department of Changchun City of China (21ZGG14)
  • 摘要: 为了准确测量地震断裂带溢出的痕量CO2气体浓度,文中采用可调谐半导体激光吸收光谱(TDLAS)技术,选取波数4978.202 cm−1作为CO2检测系统的吸收谱线,采用有效光程为40 m的多通池,以STM32作为主控和数据处理核心器件,研制了高精度CO2检测系统。针对系统中的探测器噪声与光学干涉条纹噪声,利用卡尔曼-小波分析算法滤波提升系统性能。实验表明,与滤波前相比,系统在50 ppmv CO2浓度下的二次谐波信噪比提升了2.06倍。在不同CO2浓度下(50、300、1 000、4 000、8 000 ppmv),系统误差为2.57%~2.66%。系统测量4 000 ppmv浓度下的CO2时检测精密度达到20.9 ppmv。利用Allan方差分析得出,积分时间在约61 s时对应的最低探测下限(MDL)为5.2 ppmv,实现了对CO2气体的高精度测量。结果表明,所设计的高精度CO2系统可以在气体检测领域为预测地震前兆提供良好前景。
  • 图  1  波数为4978.202 cm−1的CO2吸收谱线

    Figure  1.  CO2 absorption spectra with wave numbers of 4978.202 cm−1

    图  2  TDLAS痕量CO2检测系统原理图

    Figure  2.  Schematic diagram of TDLAS trace amount of CO2 detection system

    图  3  利用卡尔曼-小波分析算法滤波前后效果对比

    Figure  3.  Comparison of the second harmonic collected by the system before and after filtering by Kalman-wavelet analysis algorithm

    图  4  50~8000 ppmv CO2浓度标定图

    Figure  4.  Diagram of 50-8000 ppmv CO2 concentration calibration

    图  5  卡尔曼-小波分析算法滤波前后CO2误差

    Figure  5.  Error result of Kalman-wavelet analysis algorithm

    图  6  系统在4000 ppmv浓度下CO2 检测稳定性对比实验

    Figure  6.  Comparison experiment of CO2 detection stability of the system at 4000 ppmv concentration

    图  7  Allan方差探测下限图

    Figure  7.  Allan deviation method detection limit

  • [1] Chen Z, Li Y, Giovanni M. Spatial and temporal variations of CO2 emissions from the active fault zones in the capital area of China [J]. Applied Geochemistry, 2020, 112: 104489. doi:  10.1016/j.apgeochem.2019.104489
    [2] Zhou X C, Sun F X, Chen Z. Degassing of CO2, CH4, Rn and Hg in the rupture zones produced by Wenchuan Ms 8.0 earthquake [J]. Acta Petrologica Sinica, 2017, 33(1): 291-303. (in Chinese)
    [3] Xiong T, Gao Ming. Mainstream NDIR breathing CO2 monitoring system based on new light chamber structure [J]. Infrared and Laser Engineering, 2020, 49(6): 20190575. (in Chinese) doi:  10.3788/irla.30_2019-0575
    [4] Zhao A X, Tang X J, Liu J H. Spectral wavelength selection and dimension reduction using Elastic Net in spectroscopy analysis [J]. Infrared and Laser Engineering, 2014, 43(6): 1977-1981. (in Chinese) doi:  10.3969/j.issn.1007-2276.2014.06.050
    [5] Liang Y T, Tian F C. Research progress of coal mine gas detection technology in China [J]. Journal of China Coal Society, 2021, 46(6): 1701-1714.
    [6] Ren Q, Chen C, Wang Y. A prototype of pppbv-level mid-infrared CO2 sensor for potential application in deep-sea natural-gas-hydrate exploration [J]. IEEE Transactions on Instrumentation and Measurement, 2020, 69(9): 7200-7208.
    [7] Schaeffer S M, Bowling D R. Long-term field performance of a tunable diode laser absorption spectrometer for analysis of carbon isotopes of CO2 in forest air [J]. Atmospheric Chemistry & Physics, 2008, 8(17): 5263-5277.
    [8] Wang B, Fan X L, Huang S. Development of VCSEL based carbon dioxide detecting system using infrared spectroscopy [J]. Laser Journal, 2020, 41(8): 22-25. (in Chinese)
    [9] Wang C Y, Zhao S Q, Shi H W. Research on CO2 detection system based on TDLAS-WMS [J]. Transducer and Microsystem Technologies, 2021, 40(3): 52-57. (in Chinese)
    [10] Jia J W, Li W, Chai H, et al. Gas detection technology algorithm based on TDLAS [J]. Infrared and Laser Engineering, 2019, 48(5): 0517007. (in Chinese) doi:  10.3788/IRLA201948.0517007
    [11] Olver F W, Lozier D, Boisvert R F. NIST Handbook of Functions [M]. New York: International Statistical Review, 2011, 1: 131-132.
    [12] Kireev S V, Kondrashov A A. Kalman's method to improve accuracy of online 13C16O2 measurement in the exhaled human breath using tunable diode laser absorption spectroscopy [J]. Laser Physics Letters, 2011, 40(6): 992-996.
    [13] Zhou X, Jin X. Harmonic wavelet analysis of TDLAS signals [J]. Infrared and Laser Engineering, 2014, 43(6): 1722-1727. (in Chinese) doi:  10.3969/j.issn.1007-2276.2014.06.005
    [14] Chen C, Chen H D, Wang Y Z. Optical design and verification of multipass cell with two spherical mirrors using space equation method [J]. IEEE Transactions on Instrumentation and Measurement, 2021, 70: 1-8.
  • [1] 孙景旭, 费强, 谢虹波, 王硕, 谢新旺.  真空低温高灵敏度中红外辐射计 . 红外与激光工程, 2023, 52(11): 20230136-1-20230136-10. doi: 10.3788/IRLA20230136
    [2] 韦敏习, 尚万里, 侯立飞, 孙奥, 车兴森, 杨国洪.  高精度X射线柱面弯晶检测方法 . 红外与激光工程, 2022, 51(10): 20211121-1-20211121-5. doi: 10.3788/IRLA20211121
    [3] 王晨, 梁勖, 林颖, 方晓东.  MOPA结构准分子激光同步触发设计 . 红外与激光工程, 2021, 50(11): 20200516-1-20200516-6. doi: 10.3788/IRLA20200516
    [4] 侯月, 黄克谨, 于冠一, 张鹏泉.  基于红外TDLAS技术的高精度CO2同位素检测系统的研制 . 红外与激光工程, 2021, 50(4): 20200083-1-20200083-5. doi: 10.3788/IRLA20200083
    [5] 李周, 李铭扬, 余毅, 何锋赟, 蔡立华.  基于二次修正提高宽动态红外辐射测量精度 . 红外与激光工程, 2020, 49(10): 20200142-1-20200142-6. doi: 10.3788/IRLA20200142
    [6] 薛庆生, 卢继涛, 孙茜, 栾晓宁, 王福鹏.  海水吸光度特性高精度检测研究 . 红外与激光工程, 2020, 49(2): 0203010-0203010. doi: 10.3788/IRLA202049.0203010
    [7] 郝利丽, 王强, 唐红霞, 牟海维, 赵远.  线性和二次电光效应下小振幅相干耦合空间孤子对 . 红外与激光工程, 2019, 48(S1): 35-42. doi: 10.3788/IRLA201948.S106006
    [8] 刘鸿彬, 李铭, 舒嵘, 胡以华, 黄庚华.  少光子灵敏度精密激光测距方法及验证 . 红外与激光工程, 2019, 48(1): 106001-0106001(7). doi: 10.3788/IRLA201948.0106001
    [9] 熊浩西, 易仕和, 丁浩林, 徐席旺, 欧阳天赐.  三维无规则不均匀折射率场光线追迹新方法 . 红外与激光工程, 2019, 48(5): 503005-0503005(9). doi: 10.3788/IRLA201948.0503005
    [10] 卢伟业, 朱晓睿, 李越胜, 姚顺春, 卢志民, 曲艺, 饶雨舟, 李峥辉.  TDLAS直接吸收法和波长调制法在线测量CO2的比较 . 红外与激光工程, 2018, 47(7): 717002-0717002(6). doi: 10.3788/IRLA201847.0717002
    [11] 李志明, 王玺, 聂劲松, 胡瑜泽.  飞秒激光诱导硅表面高频周期结构 . 红外与激光工程, 2018, 47(1): 106003-0106003(6). doi: 10.3788/IRLA201847.0106003
    [12] 吴志豪, 施安存, 孙悦, 段靖远, 盖轶冰.  一种新型VLC高精度混合定位算法的实现 . 红外与激光工程, 2017, 46(11): 1126001-1126001(6). doi: 10.3788/IRLA201746.1126001
    [13] 秦侠格, 姬忠鹏, 徐映宇, 舒嵘.  双光路互参考高精度AOTF衍射效率测试方法及装置 . 红外与激光工程, 2017, 46(4): 417001-0417001(6). doi: 10.3788/IRLA201746.0417001
    [14] 范慧敏, 裘桢炜, 袁银麟, 康晴, 洪津.  通道式偏振遥感器偏振解析方向测量误差分析及验证 . 红外与激光工程, 2017, 46(12): 1217007-1217007(8). doi: 10.3788/IRLA201746.1217007
    [15] 伍越, 陈国柱, 沈咏, 刘曲, 邹宏新.  二次谐波产生的光学倍频腔自锁定系统 . 红外与激光工程, 2016, 45(6): 606002-0606002(7). doi: 10.3788/IRLA201645.0606002
    [16] 何韵, 许文渊.  一种高精度大动态范围的距离模拟脉冲发生方法研究 . 红外与激光工程, 2015, 44(11): 3408-3412.
    [17] 卢发铭, 张盛, 夏元钦, 陈德应.  双色飞秒强激光作用下的CO2分子高次谐波 . 红外与激光工程, 2014, 43(1): 77-80.
    [18] 周鑫, 金星.  谐波小波在TDLAS 信号分析中的应用 . 红外与激光工程, 2014, 43(6): 1722-1727.
    [19] 陈启梦员, 张国玉, 王哲, 王凌云, 高玉军.  甚高精度星敏感器测试用静态星模拟器设计 . 红外与激光工程, 2014, 43(6): 1830-1835.
    [20] 郭方, 练国富, 张宁.  复眼三维目标定位精度分析 . 红外与激光工程, 2014, 43(12): 4088-4093.
  • 加载中
图(7)
计量
  • 文章访问数:  190
  • HTML全文浏览量:  100
  • PDF下载量:  76
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-05
  • 修回日期:  2021-11-27
  • 录用日期:  2021-12-06
  • 刊出日期:  2023-03-25

基于近红外吸收光谱技术的高精度CO2检测系统的研制

doi: 10.3788/IRLA20210828
    作者简介:

    李恒宽,男,硕士生,主要从事红外光谱检测技术温控系统研究

    通讯作者: 李美萱,女,讲师,博士,主要从事多光谱关联成像技术、光学系统设计与制造方面的研究
基金项目:  国家重点研发计划 (2018YFC1503802);吉林省科技发展计划项目(20200201050JC, 20220203016SF);吉林省发改委产业技术研究与开发项目 (2022C045-5);长春市科技发展计划项目 (21ZGG14)
  • 中图分类号: TN219

摘要: 为了准确测量地震断裂带溢出的痕量CO2气体浓度,文中采用可调谐半导体激光吸收光谱(TDLAS)技术,选取波数4978.202 cm−1作为CO2检测系统的吸收谱线,采用有效光程为40 m的多通池,以STM32作为主控和数据处理核心器件,研制了高精度CO2检测系统。针对系统中的探测器噪声与光学干涉条纹噪声,利用卡尔曼-小波分析算法滤波提升系统性能。实验表明,与滤波前相比,系统在50 ppmv CO2浓度下的二次谐波信噪比提升了2.06倍。在不同CO2浓度下(50、300、1 000、4 000、8 000 ppmv),系统误差为2.57%~2.66%。系统测量4 000 ppmv浓度下的CO2时检测精密度达到20.9 ppmv。利用Allan方差分析得出,积分时间在约61 s时对应的最低探测下限(MDL)为5.2 ppmv,实现了对CO2气体的高精度测量。结果表明,所设计的高精度CO2系统可以在气体检测领域为预测地震前兆提供良好前景。

English Abstract

    • 每年中国有上千万人因自然灾害造成生命财产损失,地震所造成的损失占比约为50%以上。为了保证人民的生命财产安全,寻找准确分析地震前兆的方法仍然是急需解决的问题。因为地壳运动会将CO2等气体排出到地表,在地震前夕断裂带附近的CO2地表浓度会出现异常,对地震带附近CO2气体进行高精度测量,可以为分析地震前兆提供重要的帮助[1-2]

      目前,主要测量CO2浓度的方法有非分散性红外线技术(NDIR)分析技术、电化学技术、色谱分析技术等,但上述方法普遍存在易受其背景气体干扰,精度较低、无法实现实时监测等缺点[3-5]。可调谐半导体激光吸收光谱技术(TDLAS)技术由于其具有不受背景气体干扰、精度高、可实时监测等优点,近年来已经成为国内外的研究热点,在气体检测领域得到广泛应用[6]。2008年,Schaeffer等利用TDLAS技术对大气中CO2浓度进行测量,现场测量得到的最低探测下限(MDL)为4.5 ppmv。2020年,王彪采用新型VCSEL激光器,研制了一套CO2气体检测系统。系统最低探测下限达到9 ppmv,测量浓度为100 ppmv的CO2系统误差为7.3%。2021年,吕淑媛利用分布反馈式激光器(DFB) 作为光源搭建了基于TDLAS技术的CO2检测系统,其最低探测下限达到15 ppmv[7-9]。然而,以上研究均没有对于TDLAS系统中存在的噪声进行分析并利用算法对系统噪声进行抑制处理,仪器的精度存在较大的提升空间。

      文中首先介绍了TDLAS原理,然后介绍了CO2检测系统的各部分组成,最后通过实验验证了采用卡尔曼小波算法的TDLAS系统具有高准确度、高稳定性、低探测下限的优秀性能,为高准确度CO2系统预测地震前兆提供良好前景。

    • 根据Beer-Lambert定律,对于单一频率激光,其通过气体吸收后的光强可表示为:

      $$ {I_\lambda } = {I_0}\exp [ - \alpha (\upsilon ,P,T)cL] $$ (1)

      式中:$ {I_0} $为入射光强;$ {I_\lambda } $为出射光强;$ \alpha $$ \upsilon ,P,T $的函数;$ \upsilon $为单色光频率;$ P $,$ T $为单色光所通过介质的压强与温度;c为气室中分子数浓度;$ L $为光程。

      波长调制技术(WMS)在TDLAS技术基础上,引入高频正弦波作为调制信号,加入到激光器驱动信号中,产生的光信号通过多通池后经由锁相放大器对吸收谱线的二次谐波进行解调。当$ \alpha (\upsilon ,P,T)cL $远小于0.05时,公式(1)变为:

      $$ {I_\lambda } = {I_0}[ - \alpha (\upsilon ,P,T)cL] $$ (2)

      经过简省后,二次谐波系数与浓度关系式可表示为:

      $$ {I_{2f}} \propto {I_0}{\alpha _0}cL $$ (3)

      式中:$ {\alpha _0} $$ \upsilon ,P,T $为定值时$ \alpha (\upsilon ,P,T) $得到的常数。当$ L $为定值时,由公式(3)可得到$ {I_{2 f}} $信号与气体浓度成正相关关系[10]

    • CO2吸收谱线的选取如图1所示。

      图  1  波数为4978.202 cm−1的CO2吸收谱线

      Figure 1.  CO2 absorption spectra with wave numbers of 4978.202 cm−1

      为了实现痕量CO2的高精度测量,同时避免与水蒸气吸收谱线重叠与交叉干扰,系统选取了较强的CO2吸收谱线(4978.202 cm−1)。系统通过调节注入DFB激光器的电流使得其输出的中心波长在2 008 nm并作为CO2的检测光源。

      TDLAS痕量CO2检测系统原理图如图2所示。首先,将信号发生器产生的低频锯齿波信号和高频正弦波信号通过加法器叠加,提供给电流源用以扫描和调制DFB激光器波长。激光器输出的光束经过准直器后射入长光程多通池,出射光由光电探测器接收。

      为了提高CO2气体浓度的探测下限,该系统采用有效光程为40 m的自主研发柱面镜多通池。该多通池可在温度范围0~40 ℃及压强范围1.333~101.325 kPa的情况下稳定工作,以确保系统在现场测量过程中的可靠性。

      多通池温度控制方面,温度采集电路通过传感器实时采集多通池的工作温度,通过16位A/D模块AD7606将多通池温度模拟信号转变为数字信号发送给STM32 F103 ZET6芯片,芯片中采集的温度值与设定值进行比较,运用PID算法对温度进行控制,后通过单片机发送给D/A模块,并输出相应的PWM控制 TEC驱动芯片,从而控制TEC实现对被控对象的进行温度控制。温度控制系统实验室内控制精密度可达到0. 01 ℃,调节时间为20 min。

      多通池压力控制方面,待测气体经比例阀进入多通池,其次经过三通由压力传感器得到多通池中实时的压强值,再通过比例阀与气泵相接。压力传感器得到的实时压强值,以模拟量的形式传给压力控制电路。压力控制电路将实时压强值与目标压强值对比,改变驱动比例阀的PWM信号的占空比,改变比例阀的开口程度,从而实现对多通池内压强的控制。

      图  2  TDLAS痕量CO2检测系统原理图

      Figure 2.  Schematic diagram of TDLAS trace amount of CO2 detection system

      电源模块采用外部24 V直流电源输入,可输出两路12 V/100 mA、一路24 V/1250 mA以及两路12 V/2 000 mA对系统进行供电,电源纹波均小于50 mV,实现了系统电源供给一体化。

      采集和信号处理部分,系统测量采用STM32 F103 ZET6芯片进行控制与数据采集。在信号处理端,使用锁相放大器(LIA)解调光电探测器输出的信号,输出二次谐波信号由AD7606芯片转换为数字信号供STM32 F103 ZET6处理器采集与进行算法处理。

    • 探测器噪声和光学条纹干涉是影响TDLAS系统实际检测结果的重要因素。探测器噪声以热噪声和散粒噪声为主,均属于白噪声。系统的光学条纹干涉噪声受标准具长度、激光器温度、激光器驱动电流频率和调制幅度的影响,其在信号中的表达式推导如公式(4)。背景信号的二次函数可以展开为:

      $$ \begin{gathered} {S_{bg}}(\upsilon ,x) = \eta {I_1}\left(1 - \frac{F}{2}\right) + \eta \delta I\left(1 - \frac{F}{2}\right)\cos (\omega t - \phi ) + \\ \frac{{\eta {I_1}F}}{2}\cos (\varphi (\upsilon ,x)) + \frac{{\eta \delta IF}}{2}\cos (\omega t - \phi )\cos (\varphi (\upsilon ,x)) \\ \end{gathered} $$ (4)

      式中:$ {I_1} $为直流输出强度,$ \eta $为仪器响应参数;$ \delta I $为电流调制强度;x为标准具长度;$ \varphi (\upsilon ,x) $=$ 4\pi \upsilon {\text{x}} $$ F $为精细度系数;$ \phi $为频率调制和强度调制的相位差。二次谐波信号主要来源于公式(4)中的第三项:

      $$ \frac{{\eta {I_1}F}}{2}\cos (\varphi (\upsilon ,x)) = \frac{{\eta {I_1}F}}{2}\cos (4\pi \upsilon x + 4\pi {\delta _\upsilon }x\cos \omega t) $$ (5)

      式中:$ {\delta _\upsilon } $为频率调制幅度;定义关于调制电流归一化的频率调制幅度$ m = 2 x{\delta _\upsilon } $$ m $代入公式(5)并进行贝塞尔展开,最终背景信号二次谐波$ \;\mu (t) $可以表示为:

      $$ \mu (t) = \eta {I_1}F\cos (4\pi x{\upsilon _1}){J_2}(2\pi m) $$ (6)

      式中:$ {J_2} $为贝塞尔函数。通过上式可知,二次谐波背景信号中的光学条纹干涉会引起二次谐波波形以频率为$ {\upsilon _1} $的低频余弦形式波动[11]。这将在很大程度上影响系统精度。

    • 为了解决TDLAS系统中的探测器噪声与光学条纹干涉的影响,文中采用了卡尔曼-小波去噪的方法。卡尔曼滤波技术是利用公式(6)更新估计和公式(7)观测估计两个基本理论方程来对待采集数据完成估值推测的:

      $$ {x_{k + 1}} = {x_k} + {\varepsilon _k} $$ (7)
      $$ {z_k} = {x_k} + {\eta _k} $$ (8)

      式中:$ {x_{k + 1}} $为系统状态矩阵;$ {z_k} $为状态矩阵的实际观测量;$ {\varepsilon _k} $$ {\eta _k} $均属于系统探测器噪声引起的相互独立的随机变量。噪声观测估计方程将新测量的数据代入到预先的数据估计中从而得到更加接近于实际值的修正数据估计。利用卡尔曼滤波对于数据的预测性,可以有效的去除系统中的探测器噪声。二次谐波信号$ {\mu _{(t)}} $经过卡尔曼滤波后,对其做小波变换并表示为:

      $$ {W_{\psi u}} = (a,b) = \frac{1}{{\sqrt a }}\int_{ - \infty }^{ + \infty } {u(t)} \psi \left(\frac{{t - b}}{a}\right){\rm{d}}t $$ (9)

      式中:${\psi _{a,b}} = \dfrac{1}{{\sqrt a }}\psi \left(\dfrac{{t - b}}{a}\right)$为母小波时间平移b和尺度伸缩a的结果。由于小波变换可以分别在时域层面和频域层面进行对信号进行分解,可以根据光学条纹干涉在频域上频率为$ \upsilon $与时域上近似余弦波波形的特征,对其进行分离,以达到去除光学条纹干涉的目的。综上,利用卡尔曼-小波去噪方法可以有效去除系统中存在的探测器噪声与光学条纹干涉的影响,从而提升系统的稳定性与准确度[12-14]

    • 系统采用卡尔曼-小波分析方法消除探测器噪声以及光学条纹干涉的影响。图3为系统对浓度为50 ppmv CO2样本吸收线的二次谐波采用卡尔曼-小波分析去噪方法前后效果对比图。

      图  3  利用卡尔曼-小波分析算法滤波前后效果对比

      Figure 3.  Comparison of the second harmonic collected by the system before and after filtering by Kalman-wavelet analysis algorithm

      图3上图中蓝色波形是滤波前的气体浓度二次谐波波形,红色波形是利用滤波后的波形。SNR由7.5提升到了15.5,是原信噪比的2.06倍。拟合残差1$ \sigma $为0.38946 V。

    • 图4展示了气体浓度标定实验的结果。

      图  4  50~8000 ppmv CO2浓度标定图

      Figure 4.  Diagram of 50-8000 ppmv CO2 concentration calibration

      实验室配气系统配备了五种不同浓度的CO2气体(50、300、1000、4000、8000 ppmv),每个浓度值均采样500点以上并取平均值,在实验室25 ℃的环境下对该仪器的气体浓度进行了标定实验。如图4所示,CO2浓度标定线性度达到了0.997 7。

    • 实验针对五种不同浓度的CO2气体(50、300、1000、4000、8000 ppmv)进行测量,对该仪器应用卡尔曼-小波分析滤波前后的气体浓度检测准确度进行了测量与对比。

      图5为卡尔曼-小波分析算法滤波前CO2气体浓度误差对比图,黑色部分为滤波前的五种不同浓度CO2气体测量误差曲线,误差为4.47%~4.83%。红色部分为滤波后的误差曲线,误差为2.57%~2.66%。可以看出,使用卡尔曼-小波分析滤波算法可以明显减小系统误差,提高检测准确度。

      图  5  卡尔曼-小波分析算法滤波前后CO2误差

      Figure 5.  Error result of Kalman-wavelet analysis algorithm

    • 系统在4000 ppmv浓度条件下进行了稳定性对比实验,结果如图6所示。

      图  6  系统在4000 ppmv浓度下CO2 检测稳定性对比实验

      Figure 6.  Comparison experiment of CO2 detection stability of the system at 4000 ppmv concentration

      图6中红色部分为1000 s内采用卡尔曼-小波分析滤波前得到的浓度结果,波动范围为3 886.3~4 122.1 ppmv,标准差为35 ppmv。蓝色部分为系统滤波后探测到的浓度结果,波动范围为3 952.3~4 062.5 ppmv,标准差为20.9 ppmv。由图可知,系统采用卡尔曼-小波分析滤波后测得的浓度数据波动范围更小,得到数据的标准差由35 ppmv下降为20.9 ppmv。系统稳定性得到了较大提升。

    • 系统对浓度为300 ppmv的CO2气体样本进行测量,然后利用Allan方差来分析探测下限。

      图7所示,图7(a)与图7(b)中呈下降趋势的红色实线表示系统在白噪声占主导的区域呈现的预期的响应。图7(a)表示系统直接测量的浓度值得到的艾伦方差与积分时间$ \tau $的关系。由图可知,最优的积分时间约72 s时所对应的最低探测下限为7.2 ppmv。图7(b)表示系统通过卡尔曼-小波分析滤波后的浓度值得到的艾伦方差与积分时间t的关系。由图可知,最优的积分时间约61 s时所对应的最低探测下限为5.5 ppmv。实验利用卡尔曼-小波分析算法提升了系统的探测下限,使得最低探测下限由7.2 ppmv提升到5.5 ppmv。

      图  7  Allan方差探测下限图

      Figure 7.  Allan deviation method detection limit

    • 文中研制了面向地震监测的高精度CO2检测系统,并针对系统噪声利用卡尔曼-小波分析算法进行滤波,实验表明滤波后系统准确度、稳定性、最低探测下线指标均得到了提升。结合该方法的系统可以使地球化学气体测量具有更广泛的应用前景,为实现准确分析地震前兆提供重要的帮助。下一步将在目前的高精度CO2检测系统研制基础上加入主成分分析(PCA)算法,通过对二次谐波主特性进行提取的方式进一步提升系统准确度与稳定性。

参考文献 (14)

目录

    /

    返回文章
    返回