搜索

文章查询

x

留言板

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

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

基于连续量子级联激光器的1103.4 cm–1处NH3混叠吸收光谱特性研究

李梦琪 张玉钧 何莹 尤坤 范博强 余冬琪 谢皓 雷博恩 李潇毅 刘建国 刘文清

基于连续量子级联激光器的1103.4 cm–1处NH3混叠吸收光谱特性研究

李梦琪, 张玉钧, 何莹, 尤坤, 范博强, 余冬琪, 谢皓, 雷博恩, 李潇毅, 刘建国, 刘文清
PDF
HTML
导出引用
导出核心图
  • 由于NH3在大气气溶胶化学中具有重要作用, 所以快速和精确反演NH3浓度对环境问题非常重要. 本文以9.05 μm的室温连续量子级联激光器(quantum cascade laser, QCL)作为光源, 采用波长扫描直接吸收可调谐二极管激光吸收光谱(tunable diode laser absorption spectroscopy, TDLAS)技术, 研究了QCL在1103.4 cm–1的光谱特性, 获得了激光器控制的温度电流与波长的关系. 设计了QCL二级温控的低压实验平台, 测量氨气在1103.4 cm–1处的6条混叠吸收线, 在降低压强的情况下谱线展宽变小, 使混叠光谱分离, 由此计算各条吸收线的线强, 进一步对测量不确定度进行分析. 针对混叠严重的光谱提出了低压分离单光谱精确反演气体浓度的方法, 并进行了实验验证. 通过与HITRAN数据库进行结果对比, 得出氨气在1103.4 cm–1的实验测量线强值与数据库偏差为2.71%—4.71%, 实验测量线强值的不确定度在2.42%—8.92%, 极低压条件下反演浓度与实际值的偏差在1%—3%.
      通信作者: 张玉钧, yjzhang@aiofm.ac.cn
    • 基金项目: 国家级-国家自然科学基金青年项目(41805124)
    [1]

    何莹 2017 博士学位论文(合肥: 中国科学技术大学)

    He Y 2017 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [2]

    Gu B J, Ju X T, Chang J, Ge Y, Vitousek P M 2015 Proc. Natl. Acad. Sci. U.S.A. 112 8792

    [3]

    Fowler D, Coyle M, Skiba U, Sutton M A, Voss M 2013 Philos. Trans. R. Soc. B-Biol. Sci. 368 0164

    [4]

    Vitousek P M, Menge D N L, Reed S C, Cleveland C C 2013 Philos. Trans. R.Soc. B-Biol. Sci. 368 0119

    [5]

    Bhattacharyya P, Roy K S, Neogi S, Adhya T K, Rao K S, Manna M C 2012 Soil Tillage Res. 124 119

    [6]

    王飞, 黄群星, 李宁, 严建华, 池涌, 岑可法 2007 物理学报 56 3867

    Wang F, Huang Q X, Li N, Yan J H, Chi Y, Cen K F 2007 Acta Phys. Sin. 56 3867

    [7]

    Tao L, Sun K, Miller D J, Khan M A, Zondlo M A 2012 Opt. Lett. 37 1358

    [8]

    陈玖英, 刘建国, 何亚柏, 王辽, 江强, 许振宇, 姚路, 袁松, 阮俊, 何俊锋, 戴云海, 阚瑞峰 2013 物理学报 62 224206

    Chen J Y, Liu J G, He Y B, Wang L, Jiang Q, Xu Z Y, Yao L, Yuan S, Ruan J, He J F, Dai Y H, Kan R F 2013 Acta Phys. Sin. 62 224206

    [9]

    王立明, 张玉钧, 李宏斌, 周毅, 尤坤, 何莹, 刘文清 2012 中国光学快报 10 74

    Wang L M, Zhang Y J, Li H B, Zhou Y, You K, He Y, Liu W Q 2012 Chin. Opt. Lett. 10 74

    [10]

    Webber M E, Baer D S, Hanson R K 2001 Appl. Opt. 40 2031

    [11]

    Xu L H, Liu Z, Yakovlev 2004 Infrared Phys. Technol. 45 31

    [12]

    Jia H, Zhao W, Cai T 2009 J. Quant.Spectrosc.Radiat. Transfer 110 347

    [13]

    Sur R, Spearrin R M, Peng W Y, Strand C L, Jeffffries J B, Enns G M, Hanson R K 2016 J. Quant.Spectrosc.Radiat. Transfer 175 90

    [14]

    Romh J E, Cacciani P, Taher F 2016 J. Mol.Spectrosc. 326 122

    [15]

    Yang S, Li J, Wang R 2017 Applied Optics and Photonics China Beijing, China, October 24, 2017 p61

    [16]

    You K, Zhang Y J, Wang L M, Li H B, He Y 2013 Adv. Mater. Res. 760 84

    [17]

    鲁一冰, 刘文清, 张玉钧, 张恺, 何莹, 尤坤, 李潇毅, 刘国华, 唐七星, 范博强, 余冬琪, 李梦琪 2019 光谱学与光谱分析 39 2657

    Lu Y B, Liu W Q, Zhang Y J, Zhang K, He Y, You K, Li X Y, Liu G H, Tang Q X, Fan B Q, Yu D Q, Li M Q 2019 Spectrosc. Spect. Anal. 39 2657

    [18]

    费业泰 2015 误差理论与数据处理 (北京: 机械工业出版社) 第83页

    Fei Y T 2015 Error Theory and Data Processing (Beijing: Machinery Industry Press) p83 (in Chinese)

    [19]

    Xin Z 2005 Pap. Reg. Sci. 8418 3

    [20]

    PogányA, Klein A, Ebert V 2015 J. Quant. Spectrosc. Radiat. Transfer 165 108

    [21]

    GoldensteinC S, Hanson R K 2015 J. Quant. Spectrosc. Radiat. Transfer 152 127

    [22]

    GoldensteinC S, Jeffries J B, Hanson R K 2013 J. Quant. Spectrosc. Radiat. Transfer 130 100

    [23]

    聂伟, 阚瑞峰, 许振宇, 杨晨光, 陈兵, 夏晖晖, 魏敏, 陈祥, 姚路, 李杭, 范雪丽, 胡佳屹 2017 物理学报 66 054207

    Nie W, Kan R F, Xu Z Y, Yang C G, Chen B, Xia H H, Wei M, Chen X, Yao L, Li H, Fan X L, Hu J Y 2017 Acta Phys.Sin. 66 054207

    [24]

    林洁丽 2001 博士学位论文 (武汉: 中国科学院研究生院 (武汉物理与数学研究所))

    Lin J L 2001 Ph. D. Dissertation (Wuhan: Chinese Academy of Sciences (Wuhan Institute of Physics and Mathematics)) (in Chinese)

  • 图 1  QCL特性测量实验装置简图

    Fig. 1.  QCL characteristic measurement experimental device diagram.

    图 2  实验气路图

    Fig. 2.  Experimental gas path diagram.

    图 3  (a)未加二级温控时信号漂移(b)加二级温控后的信号稳定输出

    Fig. 3.  (a) Signal drift without secondary temperature control; (b) signal stable output after adding temperature control

    图 4  QCL的温度、电流与波长的关系

    Fig. 4.  The temperature, current and wavelength of QCL

    图 5  (a)常压下氨气特征吸收光谱信号(b)不同压力下氨气特征吸收光谱信号

    Fig. 5.  (a) Characteristic absorption spectrum signal of ammonia gas under normal pressure; (b) characteristic absorption spectrum signal of ammonia gas under different pressure

    图 6  不同压力下各条吸收线拟合结果

    Fig. 6.  Fitting results of various absorption lines under different pressures.

    图 7  实验测量的各吸收线的光谱线强值

    Fig. 7.  Spectral line strength values of each absorption line measured by experiment.

    表 1  NH3吸收线及主要参数

    Table 1.  NH3 absorption line and main parameters.

    同位素
    Isotopologue
    波数
    ν/cm–1
    线强
    S/(cm–1/(molecule cm–2))
    空气展宽
    γair/cm–1·atm–1
    自展宽
    γself/cm–1·atm–1
    爱因斯坦系数
    A/s–1
    14NH31103.430477.119 × 10200.08750.4527.183
    14NH31103.434326.131 × 10200.09330.5215.263
    14NH31103.441221.514 × 10–190.08180.3888.654
    14NH31103.469787.756 × 10200.07630.3289.694
    14NH31103.479518.172 × 10200.09940.5952.877
    14NH31103.485757.824 × 10200.0710.27410.32
    下载: 导出CSV

    表 2  实验测量线强与HITRAN数据库对比分析

    Table 2.  Comparison of experimental measurement line strength and HITRAN database.

    v0/cm–1SH/cm–2atm–1SM/cm–2atm–1E/%
    1103.430471.783391.704394.43
    1103.434321.535891.463624.71
    1103.441223.792753.69012.71
    1103.469781.942971.878253.33
    1103.479512.047181.970893.73
    1103.485751.960011.883833.89
    下载: 导出CSV

    表 3  混叠光谱线强测量不确定度

    Table 3.  Uncertainty in the measurement of the intensity of the overlapping spectral line.

    积分吸光度的
    不确定度 ΔA/%
    压强的不
    确定度 ΔP/%
    浓度的不
    确定度 Δm/%
    光程的不
    确定度 ΔL/%
    温度的不
    确定度 ΔT/%
    线强的不
    确定度 ΔS/%
    min2.420.150.20.20.12.44
    max8.928.93
    下载: 导出CSV

    表 4  实验测量线强反演浓度与HITRAN数据库线强反演浓度作对比分析

    Table 4.  Contrast analysis of experimentally measured line strong inversion concentration and HITRAN database line strong inversion concentration.

    实验气压/
    torr
    HITRAN数据库线强
    反演浓度 m1/ppm
    实验测量线强反演
    浓度 m2/ppm
    m1与实际浓度
    偏差/%
    m2与实际浓度
    偏差/%
    3.2104.98102.554.982.55
    8.5103.89101.943.891.94
    12.6102.90101.402.901.40
    19.7103.08101.073.081.07
    25.4103.22101.003.221.00
    31.8106.23102.886.232.88
    90.7108.53104.918.534.91
    760108.86104.938.864.93
    下载: 导出CSV

    表 5  不同压力、标气浓度的比较结果

    Table 5.  Comparison results of different pressure and standard gas concentrations.

    压力P/torr
    反演浓度
    m1, m2/ppm
    标气浓度m/ppm
    3.28.512.619.725.4760
    151.68m1156.68155.32154.73154.79154.81159.62
    m2154.12153.54153.05152.77152.70155.48
    202.00m1207.13206.67205.14205.45205.86210.97
    m2205.25204.73204.29203.88203.82207.06
    250.34m1256.08255.03254.26254.45254.71260.25
    m2254.37253.79253.25252.91252.86256.61
    413.00m1420.47419.73418.18418.55418.94423.77
    m2419.65418.33417.61416.87416.76420.51
    下载: 导出CSV
  • [1]

    何莹 2017 博士学位论文(合肥: 中国科学技术大学)

    He Y 2017 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [2]

    Gu B J, Ju X T, Chang J, Ge Y, Vitousek P M 2015 Proc. Natl. Acad. Sci. U.S.A. 112 8792

    [3]

    Fowler D, Coyle M, Skiba U, Sutton M A, Voss M 2013 Philos. Trans. R. Soc. B-Biol. Sci. 368 0164

    [4]

    Vitousek P M, Menge D N L, Reed S C, Cleveland C C 2013 Philos. Trans. R.Soc. B-Biol. Sci. 368 0119

    [5]

    Bhattacharyya P, Roy K S, Neogi S, Adhya T K, Rao K S, Manna M C 2012 Soil Tillage Res. 124 119

    [6]

    王飞, 黄群星, 李宁, 严建华, 池涌, 岑可法 2007 物理学报 56 3867

    Wang F, Huang Q X, Li N, Yan J H, Chi Y, Cen K F 2007 Acta Phys. Sin. 56 3867

    [7]

    Tao L, Sun K, Miller D J, Khan M A, Zondlo M A 2012 Opt. Lett. 37 1358

    [8]

    陈玖英, 刘建国, 何亚柏, 王辽, 江强, 许振宇, 姚路, 袁松, 阮俊, 何俊锋, 戴云海, 阚瑞峰 2013 物理学报 62 224206

    Chen J Y, Liu J G, He Y B, Wang L, Jiang Q, Xu Z Y, Yao L, Yuan S, Ruan J, He J F, Dai Y H, Kan R F 2013 Acta Phys. Sin. 62 224206

    [9]

    王立明, 张玉钧, 李宏斌, 周毅, 尤坤, 何莹, 刘文清 2012 中国光学快报 10 74

    Wang L M, Zhang Y J, Li H B, Zhou Y, You K, He Y, Liu W Q 2012 Chin. Opt. Lett. 10 74

    [10]

    Webber M E, Baer D S, Hanson R K 2001 Appl. Opt. 40 2031

    [11]

    Xu L H, Liu Z, Yakovlev 2004 Infrared Phys. Technol. 45 31

    [12]

    Jia H, Zhao W, Cai T 2009 J. Quant.Spectrosc.Radiat. Transfer 110 347

    [13]

    Sur R, Spearrin R M, Peng W Y, Strand C L, Jeffffries J B, Enns G M, Hanson R K 2016 J. Quant.Spectrosc.Radiat. Transfer 175 90

    [14]

    Romh J E, Cacciani P, Taher F 2016 J. Mol.Spectrosc. 326 122

    [15]

    Yang S, Li J, Wang R 2017 Applied Optics and Photonics China Beijing, China, October 24, 2017 p61

    [16]

    You K, Zhang Y J, Wang L M, Li H B, He Y 2013 Adv. Mater. Res. 760 84

    [17]

    鲁一冰, 刘文清, 张玉钧, 张恺, 何莹, 尤坤, 李潇毅, 刘国华, 唐七星, 范博强, 余冬琪, 李梦琪 2019 光谱学与光谱分析 39 2657

    Lu Y B, Liu W Q, Zhang Y J, Zhang K, He Y, You K, Li X Y, Liu G H, Tang Q X, Fan B Q, Yu D Q, Li M Q 2019 Spectrosc. Spect. Anal. 39 2657

    [18]

    费业泰 2015 误差理论与数据处理 (北京: 机械工业出版社) 第83页

    Fei Y T 2015 Error Theory and Data Processing (Beijing: Machinery Industry Press) p83 (in Chinese)

    [19]

    Xin Z 2005 Pap. Reg. Sci. 8418 3

    [20]

    PogányA, Klein A, Ebert V 2015 J. Quant. Spectrosc. Radiat. Transfer 165 108

    [21]

    GoldensteinC S, Hanson R K 2015 J. Quant. Spectrosc. Radiat. Transfer 152 127

    [22]

    GoldensteinC S, Jeffries J B, Hanson R K 2013 J. Quant. Spectrosc. Radiat. Transfer 130 100

    [23]

    聂伟, 阚瑞峰, 许振宇, 杨晨光, 陈兵, 夏晖晖, 魏敏, 陈祥, 姚路, 李杭, 范雪丽, 胡佳屹 2017 物理学报 66 054207

    Nie W, Kan R F, Xu Z Y, Yang C G, Chen B, Xia H H, Wei M, Chen X, Yao L, Li H, Fan X L, Hu J Y 2017 Acta Phys.Sin. 66 054207

    [24]

    林洁丽 2001 博士学位论文 (武汉: 中国科学院研究生院 (武汉物理与数学研究所))

    Lin J L 2001 Ph. D. Dissertation (Wuhan: Chinese Academy of Sciences (Wuhan Institute of Physics and Mathematics)) (in Chinese)

  • [1] 李宁, 翁春生. 基于多波长激光吸收光谱技术的气体浓度与温度二维分布遗传模拟退火重建研究. 物理学报, 2010, 59(10): 6914-6920. doi: 10.7498/aps.59.6914
    [2] 凌六一, 谢品华, 林攀攀, 黄友锐, 秦敏, 段俊, 胡仁志, 吴丰成. 基于O2-O2吸收的非相干宽带腔增强吸收光谱浓度反演方法研究. 物理学报, 2015, 64(13): 130705. doi: 10.7498/aps.64.130705
    [3] 王传位, 李宁, 黄孝龙, 翁春生. 基于多角度投影激光吸收光谱技术的两段式速度分布流场测试方法. 物理学报, 2019, 68(24): 247801. doi: 10.7498/aps.68.20191223
    [4] 聂伟, 阚瑞峰, 许振宇, 杨晨光, 陈兵, 夏晖晖, 魏敏, 陈祥, 姚路, 李杭, 范雪丽, 胡佳屹. 66116618 cm-1之间氨气光谱线强的测量. 物理学报, 2017, 66(5): 054207. doi: 10.7498/aps.66.054207
    [5] 林桂江, 周志文, 赖虹凯, 李 成, 陈松岩, 余金中. Si/SiGe量子级联激光器的能带设计. 物理学报, 2007, 56(7): 4137-4142. doi: 10.7498/aps.56.4137
    [6] 周康, 黎华, 万文坚, 李子平, 曹俊诚. 太赫兹量子级联激光器频率梳的色散. 物理学报, 2019, 68(10): 109501. doi: 10.7498/aps.68.20190217
    [7] 徐刚毅, 李爱珍. 量子级联激光器有源核中界面声子的特性研究. 物理学报, 2007, 56(1): 500-506. doi: 10.7498/aps.56.500
    [8] 黎华, 韩英军, 谭智勇, 张戎, 曹俊诚. 半绝缘等离子体波导太赫兹量子级联激光器工艺研究. 物理学报, 2010, 59(3): 2169-2172. doi: 10.7498/aps.59.2169
    [9] 谭智勇, 陈镇, 韩英军, 张戎, 黎华, 郭旭光, 曹俊诚. 基于太赫兹量子级联激光器的无线信号传输的实现. 物理学报, 2012, 61(9): 098701. doi: 10.7498/aps.61.098701
    [10] 万文坚, 尹嵘, 谭智勇, 王丰, 韩英军, 曹俊诚. 2.9THz束缚态向连续态跃迁量子级联激光器研制. 物理学报, 2013, 62(21): 210701. doi: 10.7498/aps.62.210701
    [11] 周超, 张磊, 李劲松. 基于单个量子级联激光器的大气多组分测量方法. 物理学报, 2017, 66(9): 094203. doi: 10.7498/aps.66.094203
    [12] 汤媛媛, 刘文清, 阚瑞峰, 张玉钧, 刘建国, 许振宇, 束小文, 张帅, 何莹, 耿辉, 崔益本. 基于室温脉冲量子级联激光器的NO气体检测中的光谱处理方法研究. 物理学报, 2010, 59(4): 2364-2368. doi: 10.7498/aps.59.2364
    [13] 朱永浩, 黎华, 万文坚, 周涛, 曹俊诚. 三阶分布反馈太赫兹量子级联激光器的远场分布特性. 物理学报, 2017, 66(9): 099501. doi: 10.7498/aps.66.099501
    [14] 李金锋, 万婷, 王腾飞, 周文辉, 莘杰, 陈长水. 太赫兹量子级联激光器中有源区上激发态电子向高能级泄漏的研究. 物理学报, 2019, 68(2): 021101. doi: 10.7498/aps.68.20181882
    [15] 常俊, 黎华, 韩英军, 谭智勇, 曹俊诚. 太赫兹量子级联激光器材料生长及表征. 物理学报, 2009, 58(10): 7083-7087. doi: 10.7498/aps.58.7083
    [16] 马欲飞, 何应, 于欣, 于光, 张静波, 孙锐. 基于中红外量子级联激光器和石英增强光声光谱的CO超高灵敏度检测研究. 物理学报, 2016, 65(6): 060701. doi: 10.7498/aps.65.060701
    [17] 周彧, 曹渊, 朱公栋, 刘锟, 谈图, 王利军, 高晓明. 基于7.6 m量子级联激光的光声光谱探测N2O气体. 物理学报, 2018, 67(8): 084201. doi: 10.7498/aps.67.20172696
    [18] 管林强, 邓昊, 姚路, 聂伟, 许振宇, 李想, 臧益鹏, 胡迈, 范雪丽, 杨晨光, 阚瑞峰. 基于可调谐激光吸收光谱技术的二硫化碳中红外光谱参数测量. 物理学报, 2019, 68(8): 084204. doi: 10.7498/aps.68.20182140
    [19] 聂伟, 阚瑞峰, 许振宇, 姚路, 夏晖晖, 彭于权, 张步强, 何亚柏. 基于TDLAS技术的水汽低温吸收光谱参数测量. 物理学报, 2017, 66(20): 204204. doi: 10.7498/aps.66.204204
    [20] 李宁, 吕晓静, 翁春生. 基于光强与吸收率非线性同步拟合的吸收光谱测量方法. 物理学报, 2018, 67(5): 057801. doi: 10.7498/aps.67.20171905
  • 引用本文:
    Citation:
计量
  • 文章访问数:  341
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-03
  • 修回日期:  2020-01-14
  • 刊出日期:  2020-04-01

基于连续量子级联激光器的1103.4 cm–1处NH3混叠吸收光谱特性研究

  • 1. 中国科学院环境光学与技术重点实验室, 安徽光学精密机械研究所, 合肥 230031
  • 2. 中国科学技术大学, 合肥 230026
  • 通信作者: 张玉钧, yjzhang@aiofm.ac.cn
    基金项目: 国家级-国家自然科学基金青年项目(41805124)

摘要: 由于NH3在大气气溶胶化学中具有重要作用, 所以快速和精确反演NH3浓度对环境问题非常重要. 本文以9.05 μm的室温连续量子级联激光器(quantum cascade laser, QCL)作为光源, 采用波长扫描直接吸收可调谐二极管激光吸收光谱(tunable diode laser absorption spectroscopy, TDLAS)技术, 研究了QCL在1103.4 cm–1的光谱特性, 获得了激光器控制的温度电流与波长的关系. 设计了QCL二级温控的低压实验平台, 测量氨气在1103.4 cm–1处的6条混叠吸收线, 在降低压强的情况下谱线展宽变小, 使混叠光谱分离, 由此计算各条吸收线的线强, 进一步对测量不确定度进行分析. 针对混叠严重的光谱提出了低压分离单光谱精确反演气体浓度的方法, 并进行了实验验证. 通过与HITRAN数据库进行结果对比, 得出氨气在1103.4 cm–1的实验测量线强值与数据库偏差为2.71%—4.71%, 实验测量线强值的不确定度在2.42%—8.92%, 极低压条件下反演浓度与实际值的偏差在1%—3%.

English Abstract

    • 近20年来, 我国的氨气排放量居全球之首[1]. 氨气是大气中仅次于N2O的碱性微量气体, 会影响区域空气的质量和大气能见度, 并导致一系列的环境问题[2-5]. 同时, 氨气极易溶于水, 并具有腐蚀性, 吸入过量氨会对人体的器官功能产生影响, 引起中毒及器官损伤, 危害人类身体健康[6]. 因此, 氨气排放特性对研究全球氮循环、气候变化和大气污染情况、保护人类健康具有重要意义.

      新兴的热电(thermoelectrics, TE)冷却分布反馈(distributed-feedback,DFB)连续波(continuous wave, CW)量子级联激光器(quantum cascade laser, QCL)能在室温下长期工作, 可用于在中红外光谱区域构建紧凑、灵敏的痕量气体传感器[7]. 可调谐二极管激光吸收光谱(tunable diode laser absorption spectroscopy, TDLAS)技术利用可调谐半导体激光器的波长可调谐特性扫描气体分子的单根或多根特征吸收谱线, 获得待测气体的特征吸收光谱, 从而对待测气体进行定量分析, 具有光谱分辨率较高、装置简单、便于操作、数据处理较简单等优点, 比较适用于各种波段较强吸收线的线参数测量[8,9]. 将这种技术与有效的中红外光源相结合, 利用分子固有振动能级在中红外波段的基频吸收, 可以检测大气中10–9甚至10–20量级的痕量气体.

      近年来, 国内外发表了多篇关于利用TDLAS技术进行NH3光谱参数测量的报道. Webber等[10]利用TDLAS技术测量了适用于环境监测和燃烧产物监测的6667 cm–1附近的6条NH3较强吸收线的线强并与其他文献中的测量值进行了对比; Xu等[11]利用外腔可调谐二极管激光器测量了6460—6522 cm–1的NH3光谱参数并与傅里叶变换红外光谱技术对比, 得出前者灵敏度更高的结论; Jia等[12]利用TDLAS技术测量了6667 cm–1附近20条NH3强吸收线并得到了它们的的线强、自展宽和空气展宽系数; Sur等[13]利TDLAS技术, 采用2个量子级联激光器, 测量了962 cm–1附近的NH3的v2谱带的Q支跃迁, 得到了吸收线的线强及温度依赖线宽系数; Romh等[14]在室温下测量得到了6369—6578 cm–1波长范围的氨气同位素(15NH3) 的吸收光谱参数; Yang等[15]基于TLAS的氨传感器, 使用5 m HWG作为气室, 利用直接吸收光谱技术, 对位于1046.4 cm–1处的无干扰NH3吸收线进行检测, 并通过HITRAN数据库的模拟模型拟合测得的谱线获得了NH3浓度.

      在之前的报道中, 鲜见关于NH3在1103.4 cm–1处吸收线测量的报道. 而1103.4 cm–1范围的部分吸收线强度较强, 其他气体吸收线对其干扰影响较小, 适用于NH3高灵敏监测. 因此, 本文设计了二级温控的低压光谱检测平台, 测量了NH3在1103.4 cm–1波段范围的吸收光谱, 通过混叠光谱分离获得了多条吸收线的线强参数, 提出利用低压分离单光谱离进行浓度精确反演的新方法, 并进行了方法验证.

    • 根据Lambert-Beer定律, 当红外光通过一段气体介质时, 由于气体吸收, 光强会发生衰减,

      ${I_t}\left( v \right) = {I_0}\left( v \right)\exp \left( { - {k_v}L} \right),$

      式中, ${I_0}\left( v \right)$为入射光强(单位: mW), ${I_t}\left( v \right)$为透射光强(单位: mW), $v$为入射光的频率(单位: cm–1), $L$为有效吸收光程(单位: cm), ${k_v}$为光谱吸收系数[16](单位: cm–1). 对于单一气体的单一跃迁而言,

      ${k_v} = S\left( T \right)\phi \left( v \right)Pm,$

      其中, $S\left( T \right)$为温度$T$时的吸收线强(单位: cm–2·atm–1, 1 atm=101325 Pa), $P$为气体压强(单位: atm), $m$为气体浓度, $\phi \left( v \right)$为吸收线的归一化线型函数(单位: cm).

      吸光度可用于描述待测气体的吸收特性, 定义吸光度${\alpha _v}$

      ${\alpha _v} = \ln \left[ {\frac{{{I_0}\left( v \right)}}{{{I_t}\left( v \right)}}} \right] = {k_v}L = S\left( T \right)\phi \left( v \right)PmL.$

      对吸光度${\alpha _v}$进行积分得到积分吸光度$A$

      $A = \int {{a_v}} {\rm{d}}v = S\left( T \right)PmL.$

      由此可见, 单一跃迁的积分吸光度面积将随压力呈线性变化[17]. 因此, 将吸光度信号进行归一化处理, 对不同压力下的测量结果(积分吸光度与光程和浓度乘积的商)进行线性拟合, 通过计算斜率即可得到线强数据, 如公式(5)所示.

      ${S_T} = \frac{A}{{PmL}}.$

      针对(5)式得到的实验线强数据, 采用不确定度反映其可信赖程度. 不确定度是指由于测量误差的存在, 对被测量值的不能肯定的程度. 它是测量结果质量的指标. 不确定度越小, 所述结果与被测量的真值越接近, 从而质量越高, 水平越高, 其使用价值越高. 假使(5)式中各变量相互独立, 则实验线强的不确定度计算如下[18]:

      $\begin{split}\Delta S =\;& \left[{{\left( {\frac{{\partial S}}{{\partial A}}} \right)}^2}\Delta _A^2 + {{\left( {\frac{{\partial S}}{{\partial P}}} \right)}^2}\Delta _P^2+ {{\left( {\frac{{\partial S}}{{\partial m}}} \right)}^2}\Delta _m^2 \right. \\ & \left.+ {{\left( {\frac{{\partial S}}{{\partial L}}} \right)}^2}\Delta _L^2 + {{\left( {\frac{{\partial S}}{{\partial T}}} \right)}^2}\Delta _T^2 \right]^{{1/2}}.\\[-18pt]\end{split}$

      通过锯齿扫描信号改变激光器的注入电流来对激光器的波长进行调谐, 激光通过测量的气体介质被探测器检测, 得到信号谱线. 通过对谱线上没有气体吸收的区域进行低阶的多项式拟合, 可以得到近似的初始激光强度, 由这2个强度能够得到随时间变化的吸光度, 再由公式(7)可反演气体浓度. 吸收谱线的线型主要分为洛伦兹线型、高斯线型和Voigt线型3种, 在常温低压状态下(10—50 torr), 吸收线为Voigt线型[19-22],

      $m = \frac{A}{{{S_T}PL}}.$

      由HITRAN数据库知, 在大气压力下, 1103.4 cm–1附近相邻且吸收强度高的6条NH3吸收线及主要参数如表1所示:

      同位素
      Isotopologue
      波数
      ν/cm–1
      线强
      S/(cm–1/(molecule cm–2))
      空气展宽
      γair/cm–1·atm–1
      自展宽
      γself/cm–1·atm–1
      爱因斯坦系数
      A/s–1
      14NH31103.430477.119 × 10200.08750.4527.183
      14NH31103.434326.131 × 10200.09330.5215.263
      14NH31103.441221.514 × 10–190.08180.3888.654
      14NH31103.469787.756 × 10200.07630.3289.694
      14NH31103.479518.172 × 10200.09940.5952.877
      14NH31103.485757.824 × 10200.0710.27410.32

      表 1  NH3吸收线及主要参数

      Table 1.  NH3 absorption line and main parameters.

      表1可知, NH3在1103.4 cm–1附近相邻2个吸收线的中心距离小于它们的展宽系数, 6根吸收线重叠形成混叠光谱. 为实现谱线分离, 便于浓度反演和光谱线强参数测量, 本文使用降压方法, 测量氨气在低压条件下的吸收情况.

    • 实验使用AdTech optics的连续DFB QCL激光器作为光源进行测量, 其具有标准的高热负荷(high heat load, HHL)封装, 边模抑制比为30 dB, 并外加底座LSM-HL0111. 温度控制器LDT-5500B和电流源LDX-3232分别用于控制QCL的温度和电流输出, 函数发生器(Tektronics, AFG3022B)向电流驱动器提供锯齿波扫描信号调谐激光稳定输出. 激光束发射穿过气池中心, 由中红外检测器(MCT-13-1.0)采集光信号并放大处理后输出. 来自检测器的信号由示波器(SIGLENT, SDS1102CNL)记录. 实验配置如图1所示.

      图  1  QCL特性测量实验装置简图

      Figure 1.  QCL characteristic measurement experimental device diagram.

      实验气路如图2所示, 由电容式薄膜真空计CDG-800(采用陶瓷膜片增加抗腐蚀能力), 有效光程长度为15 cm的单程不锈钢气体池(两端窗片材料是硒化锌, 带3°斜角以防止干涉现象的产生, 并内喷特氟龙涂料以减少氨气的吸附性), 气体流量计MF5706, 真空泵JH-VYB-110和Swagelok低压耐腐蚀气阀和光滑的不锈钢管道连接构成. 气路其中一端通过流量计连接标准气瓶, 控制匀速缓慢通气; 另一端经过单向阀连接真空泵, 抽气后静置20 min再读数. 实验前采用SWAGELOK检漏液(MS-SNOOP-8OZ)对气路进行检测, 发现无明显泄漏后开始实验.

      图  2  实验气路图

      Figure 2.  Experimental gas path diagram.

      当激光器长期工作时, 由于QCL激光器的实际温度调谐范围小(大约在25—35 ℃), 当设定温度值逼近温度调谐范围边界时, 激光器的温度无法长时间稳定在设定值(图3(a)). 因此需要优化局部散热及控温功能. 为了扩大调谐范围, 提高系统测量的长期稳定性, 利用LTC1923脉宽调制器为主芯片设计了二次温度控制电路模块. 加二级温控功能后, 激光器温度可调谐范围增大到15—42 ℃, 且能长期稳定在0.01 ℃范围内不变, 激光器输出波长无明显漂移(14 d后), 如图3(b)所示.

      图  3  (a)未加二级温控时信号漂移(b)加二级温控后的信号稳定输出

      Figure 3.  (a) Signal drift without secondary temperature control; (b) signal stable output after adding temperature control

      图4显示了氨气浓度为100 ppm, 锯齿波扫描信号频率为50 Hz, 峰峰值为0.1 V, 占空比为9: 1的常压情况下, 不同电流和工作温度下QCL的中心波长位置.

      图  4  QCL的温度、电流与波长的关系

      Figure 4.  The temperature, current and wavelength of QCL

      图4可知, QCL的中心波长随温度和电流的增大而线性增大, 温度调谐系数为0.83 nm/K, 电流调谐系数为84 nm/A. 由于实验需要测定1103.4 cm–1附近的吸收线, 故应设定温度为33.8 ℃, 电流为400 mA左右.

      常压下(760 torr)氨气在1103.4 cm–1附近的吸收情况如图5(a)所示. 此时吸收谱线已连接成片, 无法分辨. 为了获得各条吸收谱线, 调节锯齿波扫描信号为1 kHz, 峰峰值为0.060 V, 占空比为9∶1, 利用真空泵不断降低压强至25.4 Torr时, 相邻吸收谱线逐渐分离至清晰可见, 图5(b)显示了压强分别为25.4, 19.7, 12.6, 8.5 和3.2 Torr下的氨气特征吸收光谱信号.

      图  5  (a)常压下氨气特征吸收光谱信号(b)不同压力下氨气特征吸收光谱信号

      Figure 5.  (a) Characteristic absorption spectrum signal of ammonia gas under normal pressure; (b) characteristic absorption spectrum signal of ammonia gas under different pressure

      图5可见, 相临6条吸收谱线会随着压力减小逐渐分离开来, 且吸收谱线的展宽随压强降低而减小. 通过提取吸光度信号, 进行时域-频域转换, 再进行Voigt线型的多线拟合得到了各条吸收线的拟合结果(图6).

      图  6  不同压力下各条吸收线拟合结果

      Figure 6.  Fitting results of various absorption lines under different pressures.

      图6可以看出, 线型拟合的一致性较好, 随着压力的降低, 混叠光谱分离得愈明显, 各条独立光谱的获取准确性越高, 以测量得到的低压吸光度信号计算线型拟合残差, 在3.2 torr的压力之下计算的拟合残差约为2.42%. 同理对实验中一组压力条件下的拟合残差进行计算, 整体在2.42%— 8.92%. 压力越低, 光谱越明晰, 多线拟合的整体残差越小.

    • 根据公式(2)计算A/mLP的线性拟合结果, 斜率代表不同位置处的吸收线强, 以1103.43047和1103.48575 cm–1为例(图7).

      图  7  实验测量的各吸收线的光谱线强值

      Figure 7.  Spectral line strength values of each absorption line measured by experiment.

      计算上述实验得到的6条氨气吸收线强, 并将其与HITRAN数据库进行对比(表2). 低压实验所测得的线强数据与HITRAN数据库中线强数据偏差在2.71%—4.71%; 1103.43047和1103.43432 cm–1这2条吸收线在极低压条件下仍有部分交叠, 造成了多线拟合时偏差相对较大.

      v0/cm–1SH/cm–2atm–1SM/cm–2atm–1E/%
      1103.430471.783391.704394.43
      1103.434321.535891.463624.71
      1103.441223.792753.69012.71
      1103.469781.942971.878253.33
      1103.479512.047181.970893.73
      1103.485751.960011.883833.89

      表 2  实验测量线强与HITRAN数据库对比分析

      Table 2.  Comparison of experimental measurement line strength and HITRAN database.

      表中, v0为无压力频移时的中心频率(源自Hitran); SH为Hitran数据库线强值, SM为测量得到的线强值; E为测量值和Hitran值的相对偏差.

      根据不确定度的公式(6)可知, 测量中线强的不确定度和吸光度、压力、浓度、光程、温度有关, 因此分析各个参数, 得到实验中的不确定度与误差分析如下: 抽气过程中由于氨气的吸附性, 真空泵无法将气体池抽至绝对真空, 这样引入的浓度值误差为10/5000 = 0.2%[23]; 根据各实验器材的说明参数可得电容式薄膜真空计测量值的不确定度为读数的0.15% (忽略气路泄露的影响), 温度测量的不确定度为0.1%, 实验中气池有效测量光程的不确定度为0.2% (这里忽略由于充入氨气后导致气体池内折射率的变化、窗片的折射以及低压下微小形变导致的光程的微小变化); 混叠吸收线提取积分吸光度不确定度的大致范围在2.42%—8.92% (忽略时域频域的转换误差).

      根据以上线强计算和不确定度与误差分析, 得到了1103.4 cm–1附近的混叠的6条氨气吸收线的线强测量不确定度(表3). 由公式(7)可计算得到实验线强范围为$2.44\% \leqslant \Delta S \leqslant 8.93\% $, 大于一般单根谱线线强(< 5%), 这是光谱混叠的结果; 线强测量的不确定度主要和混叠吸光度的分离和提取有关.

      积分吸光度的
      不确定度 ΔA/%
      压强的不
      确定度 ΔP/%
      浓度的不
      确定度 Δm/%
      光程的不
      确定度 ΔL/%
      温度的不
      确定度 ΔT/%
      线强的不
      确定度 ΔS/%
      min2.420.150.20.20.12.44
      max8.928.93

      表 3  混叠光谱线强测量不确定度

      Table 3.  Uncertainty in the measurement of the intensity of the overlapping spectral line.

      其中, $\Delta A$是计算得到的积分吸光度的不确定度, $\Delta P$是实验压强的不确定度, $\Delta m$是实验浓度的不确定度, $\Delta L$是实验有效光程的不确定度, $\Delta T$是实验温度的不确定度.

    • 当吸收段有光谱重叠时, 通常将该吸收段所有吸收线线强之和作为${S_T}$, 重叠的吸光度信号直接线型拟合作为总吸光度, 然而实验数据和现有研究表明, 这种近似方法包含了原本无“吸收”区域, 具有不可忽视的相对误差. 研究表明[24], 将2条部分重叠的谱线视为单条谱线时所得到的压致展宽线宽与其实际值的相对误差$\varepsilon $, 与实际线宽和中心位置距离有关, 在中心位置距离小于20 MHz时, 谱线重叠部分越多, 这种误差就越大, 因此本文提出了采用降低压强分离单光谱的算法, 依次计算单一吸收线的积分吸光度反演气体浓度, 再求平均值((8)式). 根据实验比较, 发现此种方法具有可行性.

      低压条件下, $n$条吸收线全部分离后, 气体浓度可以表示为

      $m = \frac{{\rm{1}}}{{nPL}}\left( {\frac{{{A_1}}}{{{S_1}}} + \frac{{{A_2}}}{{{S_2}}} + \cdots \frac{{{A_n}}}{{{S_n}}}} \right).$

      其中, ${A_{\rm{1}}}, {A_2}, \cdots, {A_n}$是每条吸收线对应的积分吸光度; ${S_{\rm{1}}}, {S_2}, \cdots, {S_n}$是每条吸收线线强. 以实验测量的100 ppm氨气在不同气压下的特征光谱结合HITRAN数据库的线强和实验测量线强, 利用不分离混叠谱线的方法和低压分离单光谱反演浓度的方法分别由(7)和(8)式反演常压和其他较低压力下的气体浓度, 并比较相对偏差, 结果如表4所示.

      实验气压/
      torr
      HITRAN数据库线强
      反演浓度 m1/ppm
      实验测量线强反演
      浓度 m2/ppm
      m1与实际浓度
      偏差/%
      m2与实际浓度
      偏差/%
      3.2104.98102.554.982.55
      8.5103.89101.943.891.94
      12.6102.90101.402.901.40
      19.7103.08101.073.081.07
      25.4103.22101.003.221.00
      31.8106.23102.886.232.88
      90.7108.53104.918.534.91
      760108.86104.938.864.93

      表 4  实验测量线强反演浓度与HITRAN数据库线强反演浓度作对比分析

      Table 4.  Contrast analysis of experimentally measured line strong inversion concentration and HITRAN database line strong inversion concentration.

      表4可知, 极低压条件下(3.2—25.4 Torr)HITRAN数据库线强反演浓度与实际值偏差约在3%—5%, 实验测量线强反演浓度与实际值偏差约在1%—3%, 略小于HITRAN数据库的计算值. 此时, 计算值与实际值的偏差随压力的升高而降低; 当压力从31.8 Torr开始逐渐向常压升高时, 计算值与实际值的偏差随压力的升高而增大, 计算值最接近实际值时的压力值在19.7—25.4 Torr之间; 当压力值小于90.7 Torr时, 使用低压分离单光谱反演浓度方法计算的偏差值比常压的不分离混叠谱线方法的偏差值小.

      选择不同浓度的标气重复上述实验步骤, 实验结果如表5所示:

      压力P/torr
      反演浓度
      m1, m2/ppm
      标气浓度m/ppm
      3.28.512.619.725.4760
      151.68m1156.68155.32154.73154.79154.81159.62
      m2154.12153.54153.05152.77152.70155.48
      202.00m1207.13206.67205.14205.45205.86210.97
      m2205.25204.73204.29203.88203.82207.06
      250.34m1256.08255.03254.26254.45254.71260.25
      m2254.37253.79253.25252.91252.86256.61
      413.00m1420.47419.73418.18418.55418.94423.77
      m2419.65418.33417.61416.87416.76420.51

      表 5  不同压力、标气浓度的比较结果

      Table 5.  Comparison results of different pressure and standard gas concentrations.

      显然, 对不同浓度气体, 低压条件下的计算偏差都明显小于常压下的偏差. 因此, 相比于常用的不分离混叠谱线的方法, 低压分离单光谱谱线反演浓度方法准确性更高.

    • 本文基于TDLAS技术, 对中心波长为9.05 µm的室温连续QCL激光器的调谐特性进行研究, 搭建了低压实验平台, 通过降低压强实现了氨气1103.4 cm–1处6条混叠光谱的高灵敏检测, 计算了各条吸收线的线强并分析了实验线强的不确定度, 进行了浓度反演. 分析结果表明, 实验线强值与HITRAN数据库中的值相差2.71%—4.71%, 线强测量的不确定度主要和混叠吸光度的分离和提取有关, 约在2.42%—8.92%, 极低压条件下实验线强反演浓度计算偏差约在1%—3%, 而HITRAN数据库中的线强值反演浓度计算偏差约在3%—5%, 以上结果均表明本次实验数据可靠. 使用不分离混叠谱线的方法与低压分离单光谱谱线方法分别反演浓度并进行对比, 发现低压分离单光谱谱线反演浓度值与原浓度偏差更小, 准确性更高, 这对于提高环境监测等领域中气体浓度反演精度具有重要意义.

参考文献 (24)

目录

    /

    返回文章
    返回