基于2.73 μm分布反馈式激光器同时在线测量呼出气体中的CO2和水汽
1 引言
呼出气体分析技术是一种无损、实时、高灵敏度的探测技术,能对某些疾病进行辅助诊断,并能检测人体的代谢状态,在医疗诊断领域有十分重要的作用。人体呼出气体中包含几千种气体分子,含量相对较多的有CO2、H2O、O2、N2等,绝大多数呼出气体的体积分数在10-9~10-12量级[1]。其中的一些呼出气体已被证实可以用作人体代谢监测、疾病诊断的呼吸生物标志物[2]。呼出气体中CO2的体积分数约为4%,可用来监测人体的代谢过程,常被作为人体氧化应激反应的生物标志物[3-4]。呼出气体分析除了可以追踪生物标志物外,还可以通过检测稳定同位素的方法来监测人体的健康状况[5],经常使用的示踪剂是
目前,对于呼出气体中CO2的检测最常用的是非分散红外(NDIR)技术[8-10],该技术通过光电探测器探测红外光的衰减量来反演得到气体的浓度,具有较高的精度和灵敏度,但容易受光源强度变化和背景漂移的影响,因此实验过程中需要经常对频率进行校准,并使用不同的滤光片。可调谐半导体激光吸收光谱技术(TDLAS)也常被应用于呼出气体中CO2的检测,已在工业和科学研究中被广泛应用[11-13]。此外,该技术具有速度快、灵敏度高等特点,可以更好地选择气体分子的特征吸收谱线。早在2002年,Roller等[14]就已利用TDLAS在5.2 μm附近对呼出气体中的CO2进行了测量。Hartmann等[15]在2 μm附近搭建了基于TDLAS的CO2检测装置,当气室程长为2.5 cm且采集频率为30 Hz时,该装置的探测灵敏度为3×10-4。然而如此小的程长会不可避免地导致出现光学干涉条纹,增加系统噪声,而且也会影响系统的稳定性,增加测量系统的复杂性[16];此外,很小的吸收程长并不利于痕量呼出气体的检测。Xiong等[17]利用TDLAS结合中空波导管(HWG)在2003 nm对呼出气体中的CO2进行了实时探测,当采集频率为25 Hz时,探测灵敏度为6×10-5。相比于传统的多通池,HWG气体样品的体积减小了几个数量级,并且有相对较大的吸收程长,但特别容易受到外界环境的干扰,不易调节[18]。Weidmann等[19]使用石英音叉增强型光声光谱技术同时对呼出气体中CO2和水汽进行了测量,得到噪声等效吸收灵敏度分别为10-8 cm-1·W-1·Hz-1/2和8×10-9 cm-1·W-1·Hz-1/2。音叉本身存在热噪声,这限制了探测灵敏度的提高。Bartlome等[20]利用中红外激光光谱仪对呼出水汽中的D2O进行了测量,在人体口服1.5 mL D2O后,发现呼出水汽中2H和1H同位素的比值增加。Stamyr等[21]通过红外腔衰荡光谱技术测量了人体呼出气体中HCN的背景水平,检测得到CO2和水汽体积分数的平均值分别为2.8%和2.0%。所查资料显示,对呼出气体中水汽的测量鲜有报道,原因是水汽在红外波段有较强的吸收带,容易对其他呼出气体的测量产生干扰,而且水汽的吸附性很强,难以消除。
本文采用波长调制光谱技术(WMS)搭建了可以同时检测呼出气体中CO2和水汽的装置,通过设置温度电流参数可以使2.73 μm分布反馈式(DFB)激光器同时覆盖CO2和水汽的吸收线,并在实验室条件下对呼吸循环过程中CO2和水汽含量的变化进行实时测量。此外,将WMS中的二次谐波(2
2 基本原理
在TDLAS中应用WMS可以有效提高探测灵敏度,这在吸收光谱测量中已被广泛采用[22-23]。在波长调制光谱测量中,将高频正弦信号叠加到DFB激光器的注入电流中,激光波长
式中:
式中:
通常取
式中:
式中:
3 实验与结果分析
3.1 实验装置
基于2.73 μm激光器的呼出气体检测装置如
3.2 谱线选择与参数优化
为了能在2.73 μm附近同时扫描到CO2和水汽的吸收线,根据HITRAN数据库[25],选取CO2和水汽的吸收线波数分别为3659.402 cm-1和3659.934 cm-1,将激光控制器的温度和电流分别设置为34 ℃和153.03 mA,通过电流扫描使激光波数的变化范围为3658.96~3660.97 cm-1,在该范围内,可产生干扰的其他呼出气体的含量是CO2和水汽的1/1010~1/104,可忽略其影响。在压力
由波长调制光谱理论可知,调制信号的各项参数均会对气体浓度的测量产生极大影响。为了得到最优的调制振幅、相位和频率,利用高纯CO2气体(体积分数≥99.999%)与高纯氮气(N2,体积分数≥99.999%)混合配制得到体积分数为4%的CO2,将其作为参考气体通入单通池中。不同压力下2
图 3. 不同压力下2f信号峰值随调制振幅的变化
Fig. 3. Variation of 2f signal peak height with modulation amplitude at different pressures
图 4. 2f信号峰值随(a)调制相位、(b)调制频率的变化
Fig. 4. Variations of 2f signal peak height with (a) modulation phase and (b) modulation frequency
3.3 浓度定标
为了得到2
图 5. 2f信号峰值与CO2体积分数之间的关系
Fig. 5. Relation between 2f signal peak height and volume fraction of CO2
图 6. (a) 2f信号峰值与CO2体积分数的拟合关系;(b)局部放大图
Fig. 6. (a) Fitted relation between 2f signal peak height and volume fraction of CO2; (b) partial enlarged detail
采用LI-610型便携式露点发生器配制体积分数为0.7%~7.3%的标准水汽。2
图 7. (a) 2f信号峰值与水汽体积分数之间的关系;(b)局部放大图
Fig. 7. (a) Relation between 2f signal peak height and volume fraction of water vapor; (b) partial enlarged detail
3.4 呼吸循环过程实时测量
对呼出气体中的CO2和水汽进行同时测量。在进行实验之前,利用真空泵将高纯氮气抽入单通池内,约1 min后将单通池抽至真空状态,控制单通池内的压力为21.33 kPa,按照
图 8. 实验测量得到的呼出气体中CO2和水汽的2f信号
Fig. 8. 2f signals of exhaled CO2 and water vapor obtained by experimental measurement
以每平均30个谱作为一个数据点,实际采集时间为184 s,时间间隔为0.92 s,对呼吸循环过程进行实时测量。
图 9. 2个志愿者呼吸过程中CO2体积分数的变化。(a)第1个志愿者;(b)第2个志愿者
Fig. 9. Variations of volume fraction of CO2 during respiration of two volunteers. (a) The first volunteer; (b) the second volunteer
图 10. 第1个志愿者呼吸过程中水汽体积分数的变化
Fig. 10. Variation of volume fraction of water vapor during respiration of the first volunteer
3.5 长时测量
为了得到该测量系统的测量精度,利用高纯氮气配制体积分数为4.35%的CO2,控制气体流速为20 mL/s,腔内压力为21.33 kPa。平均2000个光谱作为一个数据点,时间间隔为56.8 s,共采集117个连续光谱,测量时长约为1.8 h。
图 11. (a) CO2体积分数偏差随测量时间的分布;(b) CO2体积分数偏差的频数分布直方图以及高斯线型拟合
Fig. 11. (a) Distribution of CO2volume fraction deviation with measurement time; (b) frequency distribution histogram of CO2 volume fraction deviation and Gaussian lineshape fitting
Allan方差[27]可以获得测量系统的稳定性以及探测极限。实验中利用高纯氮气配制体积分数为4.35%的CO2,时间间隔为56.8 s,测量时长约为2 h,得到测量系统的Allan方差如
4 结论
本课题组开展了基于TDLAS的呼出气体中CO2和水汽的测量研究,使用中心波长为2.73 μm的DFB激光器作为探测光源,同时覆盖CO2和水汽的吸收线。利用WMS对数据进行处理,在有效光程为88 cm及积分时间为0.92 s时,CO2和水汽的探测灵敏度可分别达到4.33×10-3和1.37×10-4。该系统连续工作近2 h,稳定性良好;同时,CO2的探测精度在10-4量级,在最大积分时间为17 min时,CO2探测极限可以达到1.49×10-4。该测量系统具有结构简单、测量便捷的特点,可以用于临床呼出气体中CO2和水汽的检测。如果将单通池更换为有效程长更大的多通池,以及利用二次谐波比一次谐波(2
[2] Stacewicz T, Bielecki Z, Wojtas J, et al. Detection of disease markers in human breath with laser absorption spectroscopy[J]. Opto-Electronics Review, 2016, 24(2): 82-94.
[4] 姜琛昱, 孙美秀, 李迎新, 等. 激光光谱技术在呼出气体分析中的发展与未来[J]. 中国激光, 2018, 45(2): 0207015.
[8] Vogt JA, NahoussiN, FabinskiW, et al. Optimised NDIR technology for 13CO2 breath tests of i.e. drug/drug-interactions or gastric emptying for intensive care patients: new diagnostic opportunities [C]∥Dössel O, Schlegel W C. Proceedings of World Congress on Medical Physics and Biomedical Engineering, September 7-12, 2009, Munich, Germany. Heidelberg: Springer, 2009: 851- 854.
[11] 高彦伟, 张玉钧, 陈东, 等. 基于可调谐半导体激光吸收光谱的氧气浓度测量研究[J]. 光学学报, 2016, 36(3): 0330001.
[13] Patterson C S. McMillan L C, Longbottom C, et al. Portable optical spectroscopy for accurate analysis of ethane in exhaled breath[J]. Measurement Science and Technology, 2007, 18(5): 1459-1464.
[26] 周超, 刘宁武, 何天博, 等. 小波阈值去噪技术在呼出气体激光光谱诊断中的应用研究[J]. 中国激光, 2017, 44(11): 1111003.
吴涛, 张怀林, 孔维平, 何兴道, 陈忠平, 陈卫东. 基于2.73 μm分布反馈式激光器同时在线测量呼出气体中的CO2和水汽[J]. 光学学报, 2018, 38(12): 1230001. Tao Wu, Huailin Zhang, Weiping Kong, Xingdao He, Zhongping Chen, Weidong Chen. Simultaneous On-Line Measurement of Exhaled Carbon Dioxide and Water Vapor Based on 2.73 μm Distributed Feedback Laser[J]. Acta Optica Sinica, 2018, 38(12): 1230001.