留言板

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

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

新型光室结构的主流式NDIR呼吸CO2监测系统

熊涛 高明 DesGibson DavidHutson

熊涛, 高明, DesGibson, DavidHutson. 新型光室结构的主流式NDIR呼吸CO2监测系统[J]. 红外与激光工程, 2020, 49(6): 20190575. doi: 10.3788/IRLA20190575
引用本文: 熊涛, 高明, DesGibson, DavidHutson. 新型光室结构的主流式NDIR呼吸CO2监测系统[J]. 红外与激光工程, 2020, 49(6): 20190575. doi: 10.3788/IRLA20190575
Xiong Tao, Gao Ming, Des Gibson, David Hutson. Mainstream NDIR breathing CO2 monitoring system based on new light chamber structure[J]. Infrared and Laser Engineering, 2020, 49(6): 20190575. doi: 10.3788/IRLA20190575
Citation: Xiong Tao, Gao Ming, Des Gibson, David Hutson. Mainstream NDIR breathing CO2 monitoring system based on new light chamber structure[J]. Infrared and Laser Engineering, 2020, 49(6): 20190575. doi: 10.3788/IRLA20190575

新型光室结构的主流式NDIR呼吸CO2监测系统

doi: 10.3788/IRLA20190575
基金项目: 陆军装备预研项目;陕西省重点实验室项目(17JS052)
详细信息
    作者简介:

    熊涛(1995-),男,硕士生,主要从事光学设计理论及技术方面的研究。Email:421505973@qq.com

    高明(1964-),男,教授,主要从事光学设计、光电测试技术方面的研究。Email:minggao1964@163.com

  • 中图分类号: O657.3

Mainstream NDIR breathing CO2 monitoring system based on new light chamber structure

  • 摘要: 针对主流式非分散红外(NDIR)呼吸CO2监测系统灵敏度及信噪比低下的问题,提出利用复合抛物面聚光器(CPC)提升监测系统性能的方法。利用ZEMAX模拟光室并优化CPC和常用的圆管聚光器以及锥形聚光器,采用单光源双光路的主流式一体化方法进行监测系统的气室设计,以STM32F100为主控芯片实现硬件系统控制,对具有不同光室结构的监测系统进行CO2浓度标定和实时人体呼吸检测实验,获取其CO2浓度与系统输出信号的对应关系和CO2波形图。结果表明,CPC光室其模拟光学效率为4.3%,比常用的聚光器最高可提升89.77倍;装有CPC的监测系统其灵敏度与信噪比为8.940 7、24.65,比装有常用聚光器的系统的灵敏度最高可提升3.811倍、信噪比提升1.926倍;并且安装CPC的系统运行稳定,响应快,能够实时的反应出被测者的呼吸CO2波形图。
  • 图  1  CPC纵剖面[12]

    Figure  1.  CPC longitudinal cross section[12]

    图  2  (a)CPC聚光系统的优化结果;(b)铝制CPC

    Figure  2.  (a) Optimization results of the CPC focusing system; (b) Aluminum CPC

    图  3  (a)安装了CPC的信号采集系统(气体腔室);(b)呼吸CO2监测系统;(c)监测系统主要部分示意图

    Figure  3.  (a) Signal acquisition system (gas chamber) with CPC; (b) Respiratory CO2 monitoring system; (c) Schematic diagram of the main parts of the monitoring system

    图  4  (a)圆管聚光器;(b)锥形聚光器;(c)复合抛物面聚光器在ZEMAX中的光线追迹图及其加工实物图;(d)腔室光电导红外探测器

    Figure  4.  (a) Straight cylinder concentrator; (b) Cone concentrator; (c) compound parabolic concentrator ray trace using ZEMAX, and samples of them; (d) Photoconductive infrared detector of the chamber

    图  5  (a) CMOS探测器(Point Grey Research, Inc)(左)以及红外光源EMIRS200(Axetris)(右);(b)CPC验证实验环境

    Figure  5.  (a) CMOS detector(Point Grey Research, Inc) (left) and IR source EMIRS200(Axetris)(right); (b) Environment of CPC verification experiment

    图  6  标定实验环境示意图

    Figure  6.  Calibration experiment environment diagram

    图  7  在CPC与探测器不同距离时,(a) ZEMAX模拟获得的探测器上光斑图;(b)验证实验中由CMOS探测器捕获的光斑图;(c) CPC与探测器不同距离时模拟结果(银色)与实验结果(橘色)的归一化的光斑截面光亮度分布比较

    Figure  7.  When the CPC and the detector are at different distances, (a) Spot pattern on the detector obtained by ZEMAX simulation; (b)Spot pattern captured by the CMOS detector in the verification experiment;(c) Comparison of the normalized light spot cross-section brightness distribution between the simulation results (silver) and experimental results (orange) at different distances between the CPC and the detector

    图  8  (a)圆管聚光器,(b)锥形聚光器,(c)复合抛物面聚光器在信号通道探测器上的光线分布图

    Figure  8.  (a)Straight cylinder concentrator, (b)cone concentrator , (c)light distribution of compound parabolic concentrator on the signal channel detector

    图  9  10 000~73 000 ppm CO2浓度范围内,三种结构对应的系统的输出信号与CO2浓度的关系

    Figure  9.  Relationship between the output signals of three systems corresponding to the CO2 concentration in the range of 10 000 ppm to 73 000 ppm CO2 concentration

    图  10  CO2波形图

    Figure  10.  CO2 waveform

    表  1  ZEMAX输入参数

    Table  1.   Input parameters in ZEMAX

    ParametersValue
    Power of the source/W1
    Type of IR sourceThermal infrared emitters (Lambertian source)
    Detector size/mm27.2×5.6
    Number of analyzing rays10 000 000
    IR source size/mm22.1×1.8
    Distance from light source to CPC/mm0.5
    Distance from CPC to detector/mm20
    下载: 导出CSV

    表  2  CPC优化后的参数

    Table  2.   Optimized parameters of CPC

    ParametersValue
    Radial aperture/mm1.79
    Angle/(°)22.86
    Length/mm9.99
    下载: 导出CSV

    表  3  三种结构的光学效率模拟结果

    Table  3.   Optical efficiency simulation results of three structures

    TypeLight efficiency
    Straight cylinder concentrator0.047%
    Cone concentrator1.2%
    Compound parabolic concentrator4.3%
    下载: 导出CSV

    表  4  系统输出值与CO2浓度的三阶多项式拟合系数及校准系数

    Table  4.   Third-order polynomial fitting coefficient and calibration coefficient of system output and CO2 concentration

    TypeabcdR2
    Straight cylinder0.000 70–0.148 710.49424.9430.969 9
    Cone0.001 5–0.314 822.86652.9630.978 1
    CPC0.001 1–0.262 324.03368.7790.985 3
    下载: 导出CSV

    表  5  监测系统分别安装三种结构后的信噪比

    Table  5.   Signal-to-noise ratio of the monitoring system after installing three structures

    TypeSNR
    Straight cylinder12.8
    Cone17.2
    CPC24.65
    下载: 导出CSV
  • [1] Lu Weiye, Zhu Xiaorui, Li Yuesheng, et al. Comparison of direct absorption and wavelength modulation methods for online measurement of CO2 by TDLAS [J]. Infrared and Laser Engineering, 2018, 47(7): 0717002. (in Chinese) doi:  10.3788/IRLA201847.0717002
    [2] Lv Shuyuan, Du Shaoyong. Detection system of CO2 using hollow-core crystal fiber based on LabVIEW [J]. Infrared and Laser Engineering, 2018, 47(11): 1117002. (in Chinese) doi:  10.3788/IRLA201847.1117002
    [3] Vaessen P H B, Knape J T A. An evaluation of the Integrated Pulmonary Index (IPI) for the detection of respiratory events in propofol sedated patients undergoing upper gastrointestinal endoscopy [J]. Journal of Hepatology and Gastrointestinal Disorders, 2016, 2(2): 117−128.
    [4] Cinar O, Acar Y A, Arziman İ, et al. Can mainstream end-tidal carbon dioxide measurement accurately predict the arterial carbon dioxide level of patients with acute dyspnea in ED [J]. The American Journal of Emergency Medicine, 2012, 30(2): 358−361. doi:  10.1016/j.ajem.2010.12.014
    [5] Singh B S, Gilbert U, Singh S, et al. Sidestream microstream end tidal carbon dioxide measurements and blood gas correlations in neonatal intensive care unit [J]. Pediatric Pulmonology, 2013, 48(3): 250−256. doi:  10.1002/ppul.22593
    [6] Mayrwöger J, Hauer P, Reichl W, et al. Modeling of infrared gas sensors using a ray tracing approach [J]. IEEE Sensors Journal, 2010, 10(11): 1691−1698. doi:  10.1109/JSEN.2010.2046033
    [7] Han J H, Han S W, Kim S M, et al. High detection performance of NDIR CO2 sensor using stair-tapered reflector [J]. IEEE Sensors Journal, 2013, 13(8): 3090−3097. doi:  10.1109/JSEN.2013.2262268
    [8] Zhang G, Wu X. A novel CO2 gas analyzer based on IR absorption [J]. Optics and Lasers in Engineering, 2004, 42(2): 219−231. doi:  10.1016/j.optlaseng.2003.08.001
    [9] Ashraf S, Mattsson C G, Thungström G, et al. Evaluation of a CO2 sensitive thermopile with an integrated multilayered infrared absorber by using a long path length NDIR platform[C]//2018 IEEE International Instrumentation and Measurement Technology Conference (I2MTC). IEEE, 2018: 1−6.
    [10] Tian M, Su Y, Zheng H, et al. A review on the recent research progress in the compound parabolic concentrator (CPC) for solar energy applications [J]. Renewable and Sustainable Energy Reviews, 2018, 82: 1272−1296. doi:  10.1016/j.rser.2017.09.050
    [11] Welford W T. High Collection Nonimaging Optics[M]. US:Elsevier, 2012.
    [12] Wang Le, Zhang Shusheng, Zhai Jing. Modeling of LED reflector cup based on compound parabolic concentrator [J]. Laser & Optoelectronics Progress, 2012(10): 24.
    [13] Vincent T A, Gardner J W. A low cost MEMS based NDIR system for the monitoring of carbon dioxide in breath analysis at ppm levels [J]. Sensors and Actuators B: Chemical, 2016, 236: 954−964. doi:  10.1016/j.snb.2016.04.016
    [14] Lin Zhe. Research and development of non-dispersive infrared CO2 gas concentration analyzer[D]. Shenyang: Shenyang Ligong University, 2008.(in Chinese)
  • [1] 朱子任, 白进周, 付景静, 苏新军, 叶静涵, 刘羽, 杨银辉, 黄文武, 李辉, 郑义军, 谭荣清.  紧凑型双路TEA CO2激光器谐振腔热稳定性及补偿 . 红外与激光工程, 2023, 52(10): 20230020-1-20230020-10. doi: 10.3788/IRLA20230020
    [2] 李恒宽, 朴亨, 王鹏, 姜炎坤, 李峥, 陈晨, 曲娜, 白晖峰, 王彪, 李美萱.  基于近红外吸收光谱技术的高精度CO2检测系统的研制 . 红外与激光工程, 2023, 52(3): 20210828-1-20210828-7. doi: 10.3788/IRLA20210828
    [3] 刘灿, 靳长明, 谭秋林, 张磊, 张雅楠.  具有温度补偿的非分光红外CO2传感器设计 . 红外与激光工程, 2022, 51(3): 20210746-1-20210746-7. doi: 10.3788/IRLA20210746
    [4] 侯月, 黄克谨, 于冠一, 张鹏泉.  基于红外TDLAS技术的高精度CO2同位素检测系统的研制 . 红外与激光工程, 2021, 50(4): 20200083-1-20200083-5. doi: 10.3788/IRLA20200083
    [5] 郑义军, 朱子任, 谭荣清, 石佳俊, 黄文武, 孙科, 李志永, 钟艳红, 李辉, 刘峻曦, 刁伟伦.  双光路快速调谐脉冲CO2激光器 . 红外与激光工程, 2020, 49(1): 0105001-0105001(6). doi: 10.3788/IRLA202049.0105001
    [6] 杨杰, 黄见, 苑克娥, 孙培育, 刘琳琳, 胡顺星.  小型化大气二氧化碳垂直廓线测量系统标定 . 红外与激光工程, 2019, 48(5): 517004-0517004(7). doi: 10.3788/IRLA201948.0517004
    [7] 左维康, 朱亚丹, 邱敏, 刘继桥, 陈卫标.  APD探测器的性能对激光雷达反演CO2浓度误差影响研究 . 红外与激光工程, 2018, 47(4): 406002-0406002(7). doi: 10.3788/IRLA201847.0406002
    [8] 邵明振, 薛向尧, 张文豹, 王光, 马亚坤, 张晨.  高功率TEA CO2激光器9R波段输出谱线 . 红外与激光工程, 2018, 47(7): 706001-0706001(5). doi: 10.3788/IRLA201847.0706001
    [9] Yasha Saxena, Chulhong Kim, Yao Junjie.  从实验室到市场:光声成像的产业化进程 . 红外与激光工程, 2017, 46(11): 1103001-1103001(14). doi: 10.3788/IRLA201746.1103001
    [10] 潘其坤, 陈飞, 石宁宁, 李殿军, 谢冀江, 何洋, 郭劲.  声光调Q CO2激光器波长调谐理论分析与实验研究 . 红外与激光工程, 2017, 46(7): 705002-0705002(5). doi: 10.3788/IRLA201746.0705002
    [11] 张阔, 陈飞, 李若斓, 杨贵龙.  大功率CO2激光器输出窗口热性能分析 . 红外与激光工程, 2017, 46(2): 205005-0205005(6). doi: 10.3788/IRLA201645.0205005
    [12] 张钊, 孙景旭, 刘则洵, 李宪圣, 任建岳.  离轴抛物面镜像面扫描太赫兹光学设计 . 红外与激光工程, 2016, 45(7): 703003-0703003(7). doi: 10.3788/IRLA201645.0703003
    [13] 于德洋, 郭立红, 陈飞, 孟范江, 杨贵龙, 邵明振.  高功率TEA CO2激光器控制系统设计 . 红外与激光工程, 2016, 45(7): 705002-0705002(6). doi: 10.3788/IRLA201645.0705002
    [14] 黄见, 胡顺星, 曹开法, 苑克娥, 时东锋.  可用于探测大气CO2垂直廓线的无线传感系统 . 红外与激光工程, 2016, 45(4): 417004-0417004(6). doi: 10.3788/IRLA201645.0417004
    [15] 张明军, 高文英, 牛泉云, 袁兴起.  聚光光伏系统菲涅耳聚光器性能分析与仿真 . 红外与激光工程, 2015, 44(8): 2411-2416.
    [16] 杨英, 冯庆.  锐钛矿型TiO2(101)面吸附CO2分子的光学气敏传感机理 . 红外与激光工程, 2015, 44(7): 2116-2121.
    [17] 何洋, 陆君, 陈飞, 李殿军, 郑长彬, 谢冀江, 郭劲.  声光调Q CO2激光器的实验研究 . 红外与激光工程, 2015, 44(8): 2280-2285.
    [18] 杨皓聿, 田爱玲, 刘丙才.  大口径抛物面镜子孔径拼接测量 . 红外与激光工程, 2014, 43(4): 1296-1300.
    [19] 郑义军, 谭荣清, 王东雷, 张阔海, 黄文武, 刘世明, 李能文, 孙科, 卢远添, 刁伟伦.  新型高重复频率脉冲CO2 激光器 . 红外与激光工程, 2013, 42(10): 2697-2701.
    [20] 吴杰, 廖光, 毛宏霞, 朱希娟, 董雁冰.  高超声速流场中CO2的红外辐射研究 . 红外与激光工程, 2013, 42(5): 1113-1116.
  • 加载中
图(10) / 表(5)
计量
  • 文章访问数:  748
  • HTML全文浏览量:  226
  • PDF下载量:  56
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-05
  • 修回日期:  2020-04-20
  • 刊出日期:  2020-07-01

新型光室结构的主流式NDIR呼吸CO2监测系统

doi: 10.3788/IRLA20190575
    作者简介:

    熊涛(1995-),男,硕士生,主要从事光学设计理论及技术方面的研究。Email:421505973@qq.com

    高明(1964-),男,教授,主要从事光学设计、光电测试技术方面的研究。Email:minggao1964@163.com

基金项目:  陆军装备预研项目;陕西省重点实验室项目(17JS052)
  • 中图分类号: O657.3

摘要: 针对主流式非分散红外(NDIR)呼吸CO2监测系统灵敏度及信噪比低下的问题,提出利用复合抛物面聚光器(CPC)提升监测系统性能的方法。利用ZEMAX模拟光室并优化CPC和常用的圆管聚光器以及锥形聚光器,采用单光源双光路的主流式一体化方法进行监测系统的气室设计,以STM32F100为主控芯片实现硬件系统控制,对具有不同光室结构的监测系统进行CO2浓度标定和实时人体呼吸检测实验,获取其CO2浓度与系统输出信号的对应关系和CO2波形图。结果表明,CPC光室其模拟光学效率为4.3%,比常用的聚光器最高可提升89.77倍;装有CPC的监测系统其灵敏度与信噪比为8.940 7、24.65,比装有常用聚光器的系统的灵敏度最高可提升3.811倍、信噪比提升1.926倍;并且安装CPC的系统运行稳定,响应快,能够实时的反应出被测者的呼吸CO2波形图。

English Abstract

    • 目前,国内外针对CO2气体有多种检测方法,其在遥感探测、农业、医学、工业等领域有着广泛的应用[1-2],常用的呼吸CO2浓度监测系统主要是根据NDIR原理设计的[3],如Cinar等人[4]在2012年利用EMMA主流式急诊CO2计对162位病人进行临床实验,通过测定病人潮气末CO2值(ETCO2)来预测病人的肺泡CO2浓度,从而判定病人是否呼吸困难。Singh等人[5]在2007年至2011年间利用Nellcor微型旁流式CO2监测仪通过监测34名ICU低出生体重(VLBW)新生儿和48名非ICU新生儿的ETCO2来预测血液CO2浓度(PCO2)。根据监测模块的测量位置,监测系统大致有两种实现方法:旁流式、主流式呼吸监测法。近年来,通过改变光学结构来提高性能的方法被广泛地用于旁流式CO2监测系统,Johann等人[6]在2011年为了提升旁流式NDIR CO2监测系统的性能,利用TracePro设计模拟了螺旋管聚光器(spiral)与椭圆管聚光器(ellipsoid),其模拟结果表明,装有spiral的系统其灵敏度比之ellipsoid系统提升了10%;Han等人[7]在2013年的报道中,对旁流式CO2监测系统中台阶锥形聚光器(stair-tapered)和锥形聚光器(cone)进行了模拟和实验验证,该团队发现stair-tapered系统的输出电压灵敏度的模拟及实验结果分别是cone系统的2倍和1.27倍;张等人[8]在2004年提出的旁流式CO2分析仪中,采用圆管聚光器进行光室设计,其样机在0~3%测量范围内的精度可达±0.028%。Ashraf等人[9]在2018年的研究中,通过在旁流式CO2监测仪光室内安装特制反射镜,来增加光线在气体中的光程从而提升系统探测下限值,实验表明,最终该监测仪可实现几个ppm至1 300 ppm的CO2浓度测量。然而,由于旁流式与主流式监测系统光室结构的不同,因此其光室设计思路存在一定差异,但是,鲜有关于通过优化光室结构的方法提升主流式CO2监测系统性能的报道。

      为了实现高信噪比、高灵敏度并且实时监测CO2浓度的目的,文中设计了一款结合CPC的NDIR主流式CO2监测系统。CPC的参数借助ZEMAX软件进行了优化,基于优化结果加工出模型之后,进行了CPC聚光效果验证实验。搭建了主流式CO2监测系统,对该系统先后安装三种不同结构的光室以对比它们对监测系统灵敏度以及信噪比的提升程度,并且对性能最佳的系统进行了实时人体呼吸CO2监测实验。

    • 呼吸CO2监测主要基于Beer-Lambert定律[7]

      $$I = {I_0}\exp ( - kcl)$$ (1)

      式中:$I_0$是入射光强度;$I$是接收到的光强度;$c$是呼吸管中的CO2浓度;l是光线在CO2气体中走过的光程;k是CO2对波长为4.26 μm的光线的吸收系数,其值为10−18 cm−1/(molecule × cm−2),它比CO2在近红外波段的吸收系数大四五个数量级。Beer-Lambert定律描述了一种通过样品对特征吸收带的强度的吸收来测量组分浓度的方法。

      在目标频带中,特征吸收频率对应于检测物质。红外光通过样品室被气体吸收一定比例的光能量,并且此比例与样品气体浓度有关,符合Beer-Lambert定律。

      从公式(1)中可以推知,相同气体浓度下,${I_0}$越大,探测器输出电压的幅值越高。因此,从理论上而言,通过提升光学模块的光能传输效率来提升监测系统的性能是可行的。

      CPC是一类具有抛物线型纵剖面的反射型聚光器(见图1),其在太阳能收集领域有着极为广泛的应用[10],其特征在于呈阶梯状的传输效率曲线,可以有效地收集从0°到最大角度(称为接收角)的光[11]

      图  1  CPC纵剖面[12]

      Figure 1.  CPC longitudinal cross section[12]

      CPC是由如下参数定义的:出光孔径半径$a'$;长度$L$;接收角${\theta _{acc}}$;抛物线焦点$f$。这五个参数可完全定义一个三维CPC,它们之间有如下基本关系:

      $$f = a'\left( {1 + \sin {\theta _{acc}}} \right)$$ (2)
      $$a' = a \cdot \sin {\theta _{acc}}$$ (3)
      $$L = (a + a') \cdot {\rm{ctg}}{\theta _{acc}}$$ (4)

      图1中右半部分抛物线坐标解析式$z = z\left( x \right)$为:

      $$z = d \cdot w\left( x \right) - {{e }}\cdot \sqrt {w\left( x \right)} + g$$ (5)

      其中$d$$e$$g$为常数:

      $$d = \frac{{\cos {\theta _{acc}}}}{{16f \cdot {{\sin }^2}{\theta _{acc}}}}$$ (6)
      $${{e}} = \frac{1}{2} + \frac{{2b \cdot \cos {\theta _{acc}}}}{{16f \cdot {{\sin }^2}{\theta _{acc}}}} + \frac{{a' \cdot \cos {\theta _{acc}}}}{{4f \cdot \sin {\theta _{acc}}}}$$ (7)
      $$\begin{split} g =\;& \frac{b}{2} + \frac{{{b^2} \cdot \cos {\theta _{acc}}}}{{16f \cdot {{\sin }^2}{\theta _{acc}}}} + \frac{{a{'^2} \cdot \cos {\theta _{acc}}}}{{4f}} + \cdots \\ &\cdots + \frac{{a'b \cdot \cos {\theta _{acc}}}}{{4f \cdot \sin {\theta _{acc}}}} + a' \cdot \sin {\theta _{acc}} - f \cdot \cos {\theta _{acc}} \end{split} $$ (8)

      其中$b$$c$是常数:

      $$b = 4f \cdot \cos {\theta _{acc}} - 2a' \cdot \sin {\theta _{acc}}$$ (9)
      $$\begin{split} c = \; & {\left( {a'} \right)^2} \cdot \sin {\theta _{acc}} - 4{f^2} \cdot \sin {\theta _{acc}} + \cdots \\ & \cdots + 4fa' \cdot \left( {1 - \cos {\theta _{acc}}} \right) \end{split} $$ (10)

      公式$w\left( x \right)$定义如下:

      $$w\left( x \right) = {b^2} - 4 \cdot \left( {4f \cdot x + c} \right) \cdot \sin {\theta _{acc}}$$ (11)
    • 搭建主流式呼吸CO2监测系统,利用ZEMAX模拟优化CPC、常用的圆管聚光器(Straight cylinder)与圆锥聚光器(Cone)并加工出实物。分别结合这三款聚光器构造监测系统的光学部分,以STM32F100为主控芯片实现硬件系统控制,利用该芯片进行信号的调制解调以及数据映射等信号处理。通过一系列CO2浓度标定以及实时人体呼吸监测实验获得安装不同聚光器时的监测系统其灵敏度、信噪比以及实时CO2波形图。

    • 由于人体呼吸气体中CO2的浓度基本处于400~50 000 ppm之间,根据Beer-Lambert定律选择呼吸管道直径为20 mm[13],并通过ZEMAX软件建立光源及探测器模型,将CPC的最佳聚焦距离优化至20 mm,该光路系统主要输入参数见表1,优化之后的CPC参数见表2

      表 1  ZEMAX输入参数

      Table 1.  Input parameters in ZEMAX

      ParametersValue
      Power of the source/W1
      Type of IR sourceThermal infrared emitters (Lambertian source)
      Detector size/mm27.2×5.6
      Number of analyzing rays10 000 000
      IR source size/mm22.1×1.8
      Distance from light source to CPC/mm0.5
      Distance from CPC to detector/mm20

      表 2  CPC优化后的参数

      Table 2.  Optimized parameters of CPC

      ParametersValue
      Radial aperture/mm1.79
      Angle/(°)22.86
      Length/mm9.99

      基于表2中CPC优化结果以及1.2节的CPC理论模型,计算出CPC的纵剖面坐标解析式:

      $$\begin{split} z\left( x \right) = \;& - 2.375x - 3.313\sqrt { - 15.521x + 71.846} + \cdots \\ & \cdots + 26.238 \end{split} $$ (12)

      将公式(12)导入数控机床加工出如图2(b)所示的CPC。

      图  2  (a)CPC聚光系统的优化结果;(b)铝制CPC

      Figure 2.  (a) Optimization results of the CPC focusing system; (b) Aluminum CPC

    • 基于上文的模拟结果,建立了如图3(a)所示的3D模型。从图3(a)中可见,安装在光源之后的CPC,可将光源发出的光线聚焦到位于呼吸管正对面的探测平面上,同时,它们都被铝制支撑结构固定在呼吸管的进、出光口。

      图  3  (a)安装了CPC的信号采集系统(气体腔室);(b)呼吸CO2监测系统;(c)监测系统主要部分示意图

      Figure 3.  (a) Signal acquisition system (gas chamber) with CPC; (b) Respiratory CO2 monitoring system; (c) Schematic diagram of the main parts of the monitoring system

      设计的呼吸CO2监测系统样机主要包括四部分:(1)由呼吸管、CPC、LED红外光源和PbSe探测器组成的信号采集系统;(2)信号调制及放大电路;(3) STM32F100 (STMicroelectronics Inc.) 微控制器,集成了高性能的ARM®Cortex®-M3数据处理模块以及用于存储测量数据的高速嵌入式存储器;(4)PC和样机之间的数据通信端口(USART接口)。样机中使用的红外光源(Gas sensing solutions Inc.)可以辐射出3.3~5.3 μm波段的红外光。此外,探测器选用的是四通道光电导红外探测器(MBXP-15-4 from OPTO DIODE,USA),每个通道有效探测面的尺寸皆为1 mm×1 mm。文中仅使用该探测器中心波长为4.26 μm滤波器(半功率宽为70 nm)的CO2通道,和用以补偿温度和湿度等环境变化的安装有中心波长为3.85 μm滤波器(半功率宽为100 nm)的参考通道。

      图3(c)所示,气体信号通道与参考通道的信号通过调理、放大电路之后分别进入ARM微控制器的两个模数转换通道并进行数字信号的解调、降噪等处理环节。ARM微控制器以及光源驱动电路使该红外LED产生约1 kHz的高频红外光辐射。此外,微控制器能同时对信号进行采样并且实时计算CO2浓度值,该值通过USART端口传输到PC端的用户界面(GUI)以便进行数据可视化以及相关诊断算法的执行。

      为了方便对比,利用ZEMAX设计了针对20 mm探测距离而优化的两款常用的聚光器:圆管聚光器,圆锥聚光器。优化设计时,将它们的外形尺寸限制在同一水平。

      图4所示为三类聚光器的模拟图、加工样品实物(其材料皆为铝材),以及光电导红外探测器。图中聚光器底部的中心为红外LED光源,这三个聚光器之间相互独立,并且易拆卸,便于监测系统切换不同的腔室光学结构。

      图  4  (a)圆管聚光器;(b)锥形聚光器;(c)复合抛物面聚光器在ZEMAX中的光线追迹图及其加工实物图;(d)腔室光电导红外探测器

      Figure 4.  (a) Straight cylinder concentrator; (b) Cone concentrator; (c) compound parabolic concentrator ray trace using ZEMAX, and samples of them; (d) Photoconductive infrared detector of the chamber

    • 在验证实验中(实验是在23 ℃室温,标准大气压环境下的暗室进行的),采用的是图5(a)所示的黑体辐射光源,它的辐射光波段为1 ~14 μm,该CMOS探测器(像元尺寸为4.5 μm,感光面像素数为1 600×1 200)可以探测0.3 ~1.1 μm波段的光信号。如图5(b)所示,为该验证实验的实验环境,波形发生器的方波信号进入电源驱动的光源调制电路,该电路输出的光源调制信号驱动光源发出光线,光线经由CPC汇聚被探测器接收,探测器接收到的信号可在PC端查看储存。

      图  5  (a) CMOS探测器(Point Grey Research, Inc)(左)以及红外光源EMIRS200(Axetris)(右);(b)CPC验证实验环境

      Figure 5.  (a) CMOS detector(Point Grey Research, Inc) (left) and IR source EMIRS200(Axetris)(right); (b) Environment of CPC verification experiment

    • 实验中使用的气体为BOC公司生产的CO2与N2图6所示为标定实验的实验环境(标定实验是在室温21 ℃,标准大气压下进行的)。通过多气体控制台内的N2流量控制仪(MFC)和CO2流量控制仪调节气体流量,两种气体在后段气体管道内充分混合后进入冷却系统(内有除湿机构),管道末端的排风扇将冷却后的混合气体引入监测系统的光学腔室内,光源出射的光线穿过混合气体(被吸收部分光能量)后被探测器接收,探测器的输出信号通过信号放大电路后进入ARM微控制器进行模数转换和信号解调,解调后的信号被输进示波器中,得到CO2波形的实时变化图形。

      图  6  标定实验环境示意图

      Figure 6.  Calibration experiment environment diagram

      分别将三种聚光器安装进监测系统的光室内,通过多气体控制台产生不同浓度的CO2混合气体,对监测系统进行浓度标定实验。

    • 图7(a)为CPC光路系统中探测器模拟结果,图7(b)为验证实验中探测器接收的信号,图7(c)为模拟结果与实验验证结果之间的比较。结果表明,当CPC与探测器距离为8 mm时,光斑呈光圈状,这是由于光源发出的光大部分都直接经由CPC的内壁反射,偏离最佳聚焦距离的探测器接收到这些被反射的光后显示为光圈,而探测器接收的从光源直接出射的光线数量相较经过反射的光线数量比较少,因此,探测器上显示为一个亮圈。当探测器逐渐向最佳聚焦距离靠近时,亮圈逐渐向中心聚拢,当探测器处于最佳聚焦距离(20 mm)时,亮圈聚拢为一个亮斑,当距离为大于20 mm时,亮斑逐渐发散。因此,CPC聚光效果满足设计需求。

      图  7  在CPC与探测器不同距离时,(a) ZEMAX模拟获得的探测器上光斑图;(b)验证实验中由CMOS探测器捕获的光斑图;(c) CPC与探测器不同距离时模拟结果(银色)与实验结果(橘色)的归一化的光斑截面光亮度分布比较

      Figure 7.  When the CPC and the detector are at different distances, (a) Spot pattern on the detector obtained by ZEMAX simulation; (b)Spot pattern captured by the CMOS detector in the verification experiment;(c) Comparison of the normalized light spot cross-section brightness distribution between the simulation results (silver) and experimental results (orange) at different distances between the CPC and the detector

      由于CPC、支撑结构以及呼吸管存在一定的加工误差;光源、CPC与探测器之间存在一定的不同轴现象;并且模拟光源与实际光源存在一定的偏差,因此验证实验的结果与模拟得到的结果很难完全拟合。但值得一提的是,在CPC与探测器的距离为20 mm时它们都表明CPC具有超越其他距离下的良好聚光性能。这对文中研究的主流式CO2监测仪的信号采集系统来说是最为重要的一点。

    • 表3为三种不同结构的聚光器光学效率模拟结果,该结果显示CPC的光学效率是锥形聚光器的3.58倍,是圆管聚光器的89.77倍。图8是探测器接收到的光线分布图,从图中亦可以看出CPC后的探测器接收到的光线数量明显大于其他两种光学结构后的探测器所接收到的。

      表 3  三种结构的光学效率模拟结果

      Table 3.  Optical efficiency simulation results of three structures

      TypeLight efficiency
      Straight cylinder concentrator0.047%
      Cone concentrator1.2%
      Compound parabolic concentrator4.3%

      图  8  (a)圆管聚光器,(b)锥形聚光器,(c)复合抛物面聚光器在信号通道探测器上的光线分布图

      Figure 8.  (a)Straight cylinder concentrator, (b)cone concentrator , (c)light distribution of compound parabolic concentrator on the signal channel detector

      图8表3的模拟结果都表明,相较于其他两种常用聚光器,CPC具有最佳的聚光性能,其抛物面型的内壁可以将相当一部分光线反射至目标探测面,由于目标探测面尺寸微小且离聚光器有一定的距离,因此锥形聚光器其锥形内壁只有小部分环带可以将光线反射至探测面。而圆管聚光器由于其内壁平行于主光轴,因此经由内壁反射至探测面光线的数量受限于聚光器与探测面的距离大小,由于文中研究的主流式CO2监测系统的该距离大小为20 mm,因此,当光学腔室结构为圆管聚光器时,经由内壁反射的光线损失了绝大部分,探测面接收到的光线多为光源直接出射的光。

    • 由于监测系统对光源辐射的非均匀响应特性、红外光源辐射能量不稳定以及红外传感器放大电路非线性等因素的影响,使得利用Beer-Lambert定律对CO2浓度标定结果存在一定的偏差,因此在文中,整个监测系统被视为一个整体,通过对标定实验获取监测系统数据,并将它和的CO2与N2的混合气体浓度进行最小二乘拟合,从而确定监测系统处理最终值与CO2浓度的对应关系,计算出当前时刻人体呼吸CO2浓度值。

      根据林喆[14]针对此问题建立的检测模型以及实验结果显示,在校准温度下,三阶多项式拟合为最小二乘建模时的最佳模型,其数学表达式如下:

      $$V = a \cdot {C^3} + b \cdot {C^2} + c \cdot C + d$$ (13)

      式中:V为从监测系统获得的实验结果;C为CO2浓度;R2为拟合校准系数。由实验得到的监测系统数据处理最终值与CO2浓度的拟合关系系数和拟合校准系数见表4

      表 4  系统输出值与CO2浓度的三阶多项式拟合系数及校准系数

      Table 4.  Third-order polynomial fitting coefficient and calibration coefficient of system output and CO2 concentration

      TypeabcdR2
      Straight cylinder0.000 70–0.148 710.49424.9430.969 9
      Cone0.001 5–0.314 822.86652.9630.978 1
      CPC0.001 1–0.262 324.03368.7790.985 3
    • 正常人体潮气末CO2浓度一般为50 000 ppm,基于10 000 ~73 000 ppm浓度范围的CO2浓度标定结果可得到图9所示线性拟合曲线。

      图  9  10 000~73 000 ppm CO2浓度范围内,三种结构对应的系统的输出信号与CO2浓度的关系

      Figure 9.  Relationship between the output signals of three systems corresponding to the CO2 concentration in the range of 10 000 ppm to 73 000 ppm CO2 concentration

      图9显示了在10 000~73 000 ppm的所有浓度标定实验数据,可以看出三种结构对相同CO2浓度的响应程度并不一致,由于它们的线性拟合曲线线性度较高,因此,其斜率直接反应出了当前结构对CO2浓度响应的灵敏度。在浓度标定实验中,CtCt-1分别表示不同的浓度的CO2混合气体;VtVt-1分别表示对应的样机输出信号值。灵敏度定义为:

      $$\kappa = \frac{{\left\| {{C_t} - {C_{t - 1}}} \right\|}}{{\left\| {{V_t} - {V_{t - 1}}} \right\|}}$$ (14)

      由公式(14)可得不同光学腔结构的监测系统灵敏度分别为:圆管聚光器:${\kappa _{{\rm{cyl}}}} = 2.346\;2$,锥形聚光器:${\kappa _{{\rm{cone}}}} = $5.454 8,CPC:${\kappa _{CPC}} = 8.940\;7$,其量纲为1 000 ppm/mV。安装CPC光学腔的样机其灵敏度是安装圆管聚光器光学腔监测系统的3.811倍,是安装锥形聚光器光学腔监测系统的1.639倍。这表明,对于相同浓度的CO2,安装CPC的系统相对其他两种结构而言会返回更高的输出电压,因此,它可以对更微小的CO2浓度变化做出响应。

    • 进行浓度标定实验时,由于样机系统中存在1/f噪声,白噪声,以及高斯噪声,它们的存在会直接影响到样机可测量的最低CO2浓度值。因此,确定样机的信噪比参数是十分必要的。文中利用公式(15)和(16)计算混合噪声的RMS值。

      $$Signa{l_{RMS}} = \sqrt {\frac{1}{N}\left( {\sum\limits_{n = i}^N {Signal_i^2} } \right)} $$ (15)
      $$Nois{e_{RMS}} = \sqrt {\frac{1}{N}{{\left( {\sum\limits_{n = i}^N {Signa{l_i} - Signa{l_{avg}}} } \right)}^2}} $$ (16)

      利用公式(17)计算分别安装三种光学腔结构的样机其信噪比。

      $$SNR = \frac{{Signa{l_{RMS}}}}{{Nois{e_{RMS}}}}$$ (17)

      表5为监测系统安装三种光学腔结构后的信噪比,可以得出,安装CPC后系统的信噪比是安装圆管聚光器系统的1.926倍,是安装锥形聚光器系统的1.433倍。由于安装CPC的系统具有更高的灵敏度,且其信噪比也是最高的,因此,相对其他两种结构而言,它可以测量出更低的CO2浓度值。

      表 5  监测系统分别安装三种结构后的信噪比

      Table 5.  Signal-to-noise ratio of the monitoring system after installing three structures

      TypeSNR
      Straight cylinder12.8
      Cone17.2
      CPC24.65
    • 对具有最佳性能的安装CPC光室的样机进行呼吸CO2浓度测试,当对人体进行呼吸监测时,管道适配器中潮湿的水蒸气(文中在冷却系统中有除湿机构,因此短时间内的测量基本不受水蒸气影响)和快速变化的呼吸气流会对CO2浓度测量产生极大的干扰,导致最后CO2浓度出现偏差。

      为了验证该样机系统的实用效果,对两位健康测试者进行呼吸CO2浓度的实时监测实验。实验前使测试者放松状态,避免有意识的控制呼吸节奏。对两位被测者进行的呼吸CO2监测,实验表明,该样机具有良好的响应速度,信噪比和灵敏度。图10为其中一位被测者的呼吸二氧化碳波形图,该图清晰的显示了呼吸CO2浓度变化规律以及EtCO2等特征信息。

      图  10  CO2波形图

      Figure 10.  CO2 waveform

    • 文中提出利用CPC提升主流式呼吸CO2监测系统性能,对CPC进行了优化、模拟和验证实验,结果表明模拟结果与实验结果具有良好的一致性;搭建了监测系统,并对圆管聚光器、锥形聚光器、CPC分别进行了模拟和对比试验,模拟结果与实验结果表明,比之圆管聚光器,CPC可使光室的光学效率提升89.77倍,使系统的灵敏度提升3.811倍,信噪比提升1.926倍;比之锥形聚光器,可使光室光学效率提升3.58倍,使系统灵敏度提升1.639倍,信噪比提升1.433倍。光室内安装CPC可使样机分辨更微小的CO2浓度变化,测量更低的CO2浓度值。对监测系统进行的实时人体呼吸CO2浓度测试结果表明,该样机具备良好的响应速度、信噪比以及灵敏度,这对进一步研发微型化高性能主流式呼吸CO2监测仪具有很好的参考价值。

参考文献 (14)

目录

    /

    返回文章
    返回