留言板

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

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

探测气溶胶-水汽的拉曼-米散射激光雷达系统

李路 邢昆明 赵明 邓迁 王邦新 庄鹏 施云

李路, 邢昆明, 赵明, 邓迁, 王邦新, 庄鹏, 施云. 探测气溶胶-水汽的拉曼-米散射激光雷达系统[J]. 红外与激光工程, 2023, 52(4): 20220484. doi: 10.3788/IRLA20220484
引用本文: 李路, 邢昆明, 赵明, 邓迁, 王邦新, 庄鹏, 施云. 探测气溶胶-水汽的拉曼-米散射激光雷达系统[J]. 红外与激光工程, 2023, 52(4): 20220484. doi: 10.3788/IRLA20220484
Li Lu, Xing Kunming, Zhao Ming, Deng Qian, Wang Bangxin, Zhuang Peng, Shi Yun. Raman-Mie scattering lidar system for detection of aerosol and water vapor in the atmosphere[J]. Infrared and Laser Engineering, 2023, 52(4): 20220484. doi: 10.3788/IRLA20220484
Citation: Li Lu, Xing Kunming, Zhao Ming, Deng Qian, Wang Bangxin, Zhuang Peng, Shi Yun. Raman-Mie scattering lidar system for detection of aerosol and water vapor in the atmosphere[J]. Infrared and Laser Engineering, 2023, 52(4): 20220484. doi: 10.3788/IRLA20220484

探测气溶胶-水汽的拉曼-米散射激光雷达系统

doi: 10.3788/IRLA20220484
基金项目: 安徽省教育厅自然科学研究重点项目(KJ2021 A0945);皖西学院高层次人才项目(WGKQ2022013,WGKQ2022015);中国科学院合肥物质科学研究院“十三五”规划重点支持项目(KP-2019-05);中国科学院战略性先导科技专项(XDA17040524);2019年安徽省科技重大专项(201903c08020013)
详细信息
    作者简介:

    李路,男,讲师,博士,主要从事大气激光雷达技术方面的研究

    通讯作者: 邢昆明,男,工程师,硕士,主要从事大气激光雷达光机结构优化设计方面的研究
  • 中图分类号: TN958.98

Raman-Mie scattering lidar system for detection of aerosol and water vapor in the atmosphere

Funds: Key Project of Natural Science Research Department of Education of Anhui Province (KJ2021A0945); High-Level Personnel Project of West Anhui University (WGKQ2022013, WGKQ2022015); Key Support Project of Hefei Institutes of Physical Science of the Chinese Academy of Sciences During the 13th Five-Year Plan (KP-2019-05); The Strategic Pilot Science and Technology Project of the Chinese Academy of Sciences (XDA17040524); Major Science Technology Project of Anhui Province in 2019 (201903C08020013).
  • 摘要: 设计和构建了发射波长为355 nm和532 nm的户外型全天时激光雷达系统,用于探测大气气溶胶和水汽。运用355 nm和532 nm的米散射、532 nm的偏振、氮气和水汽分子的拉曼激光雷达技术,用于对边界层结构、对流层气溶胶和云光学特性及其形态、水汽混合比进行连续探测研究。该系统结构紧凑,运输方便,具备远程操作、数据传输、一键式启动等功能。利用该系统对大气气溶胶和水汽进行探测,探测结果表明:在大气气溶胶的探测过程中,在重污染条件下混合层高度较干净天低,在0.5 km以下,而干净天在1 km左右;通过对消光系数、Angstrom指数和退偏振比分析可知,重污染条件下,底层大气气溶胶以球形粗粒子污染物为主,干净天底层大气气溶胶以球形细粒子污染物为主;在云层中,Ang-strom指数明显减小,且出现负值,说明云粒子半径较大。在水汽探测过程中,采用自标定方法获得系统的标定常数为121,与已标定的激光雷达系统对比,误差在±0.3 g/kg以内;连续探测结果表明可对夜晚5 km及白天混合层以内进行探测。该系统满足产品化的需求,可广泛运用于大气环境的监测领域中。
  • 图  1  大气气溶胶-水汽激光雷达产品结构示意图

    Figure  1.  Structural schematic of lidar system for atmospheric aerosol-water vapor

    图  2  分束片安装结构

    Figure  2.  Beam splitter mounting structure

    图  3  后继探测单元

    Figure  3.  Subsequent detection unit

    图  4  重度污染混合层高度

    Figure  4.  Mixed-layer height in the heavy pollution conditions

    图  5  重度污染天的消光系数(a)、Angstrom指数(b)及退偏振比(c)

    Figure  5.  Extinction coefficient (a), Angstrom index (b) and depolarization ratio (c) in the heavy pollution conditions

    图  6  干净天混合层高度

    Figure  6.  Mixed-layer height in the clean weather conditions

    图  7  干净天的消光系数(a)、Angstrom指数(b)及退偏振比(c)

    Figure  7.  Extinction coefficient (a), Angstrom index (b) and depolarization ratio (c) in the clean weather conditions

    图  8  水汽混合比垂直分布特征及绝对误差廓线图

    Figure  8.  Profile of vertical distribution feature of water vapor mixing ratio and absolute error

    图  9  水汽混合比时空分布

    Figure  9.  Spatiotemporal distribution of water vapor mixing ratio

    表  1  大气气溶胶-水汽激光雷达系统的主要参数、性能指标

    Table  1.   Main specifications, performance parameters of lidar system for atmospheric aerosol-water vapor

    Parameters
    Wavelength/nm532/355
    Frequency/Hz10
    Monopulse energy/mJ88/50
    Telescope diameter/mm400
    Field of view/mrad0.5-2(adjustable)
    Detection channel5 (355 nm, 532 nm S, 532 nm P, 386 nm and 407 nm)
    Filter bandwidth/nm1
    Capture cardA/D and photon counting
    Performance
    Measuring distance
    Spatial resolution
    Time resolution
    Operation temperature
    Storage temperature
    Way of working
    Power supply
    Volume
    0-12 km(atmosphere of clear sky)
    7.5 m
    15 min(adjustable depending on the measurement content)
    5- −40 ℃(−10- −40 ℃ under the condition of peripheral temperature control)
    5-−70 ℃(−30- −70 ℃ under the condition of peripheral temperature control)
    Manual, automatic and network control
    220 V/50 Hz/1500 W
    1.8 m×1.2 m×2 m(L×W×H
    下载: 导出CSV
  • [1] Han Yong, Wang Tijian, Rao Ruizhong, et al. Progress in the study of physic-optics characteristics of atmospheric aerosols [J]. Acta Physica Sinica, 2008, 57(11): 7396-7407. (in Chinese) doi:  10.7498/aps.57.7396
    [2] Wang Yufeng, Gao Fei, Zhu Chengxuan, et al. Raman Lidar for atmospheric temperature, humidity and aerosols up to tropo-sphere height [J]. Acta Optica Sinica, 2015, 35(3): 0328004. (in Chinese) doi:  10.3788/AOS201535.0328004
    [3] Whiteman D N. Examination of the traditional raman lidar technique. II. Evaluating the ratios for water vapor and aerosols [J]. Applied Optics, 2003, 42(15): 2593-608. doi:  10.1364/AO.42.002593
    [4] Shang Zhen, Xie Chenbo, Zhong Zhiqing, et al. Raman lidar for measurement of tropospheric water vapor [J]. Infrared and Laser Engineering, 2016, 45(12): 1211003. (in Chinese) doi:  10.3788/IRLA201645.1211003
    [5] Xie Chenbo, Zhou Jun, Yue Guming, et al. Mobile lidar system for measuring tropospheric aerosol and water vapor [J]. Infrared and Laser Engineering, 2007, 36(3): 365-368. (in Chinese) doi:  10.3969/j.issn.1007-2276.2007.03.021
    [6] Lv Weiyu, Yuan Ke′e, Wei Xu, et al. A mobile lidar system for aerosol and water vapor detection in troposphere with mobile lidar [J]. Infrared and Laser Engineering, 2016, 45(3): 0330001. (in Chinese) doi:  10.3788/irla201645.0330001
    [7] Tan Min, Wang Bangxin, Zhuang Peng, et al. Study on atmospheric temperature and water-vapor mixing ratio based on raman lidar [J]. Spectroscopy and Spectral Analysis, 2020, 40(5): 1397-1401. (in Chinese) doi:  10.3964/j.issn.1000-0593(2020)05-1397-05
    [8] 谭敏. 污染条件下气溶胶和水汽时空分布探测和数据分析研究[D]. 中国科学技术大学, 2018.

    Tan Min. Study on spatial and temporal distribution and data analysis of aerosol and water vapor under pollution conditions[D]. Hefei: University of Science and Technology of China, 2018. (in Chinese)
    [9] Di Huige, Hou Xiaolong, Zhao Hu, et al. Detections and analyses of aerosol optical properties under different weather conditions using multi-wavelength Mie lidar [J]. Acta Physica Sinica, 2014, 63(24): 244-251. (in Chinese) doi:  10.7498/aps.63.244206
    [10] Liu Dong, Tao Zongming, Wu Decheng, et al. Development of three-wavelength-raman- polarization lidar system and case study [J]. Acta Optica Sinica, 2013, 33(2): 0228001. (in Chinese) doi:  10.3788/AOS201333.0228001
    [11] Chi Ruli, Wu Decheng, Liu Bo, et al. Dual-wavelength Mie lidar observations of tropospheric aerosols [J]. Spectroscopy and Spectral Analysis, 2009, 29(6): 1468-1472. (in Chinese) doi:  10.3964/j.issn.1000-0593(2009)06-1468-05
    [12] Chen Tao, Wu Decheng, Liu Bo, et al. A new method for determining aerosol backscatter coefficient boundary value in the lower troposphere [J]. Acta Optica Sinica, 2010, 30(6): 1531-1536. (in Chinese) doi:  10.3788/AOS20103006.1531
    [13] Deng Qian, Wu Decheng, Kuang Zhiqiang, et al. 532 nm/660 nm dual wavelength lidar for self-calibration of water vapor mixing ratio [J]. Infrared and Laser Engineering, 2018, 47(12): 1230004. (in Chinese) doi:  10.3788/IRLA201847.1230004
    [14] Wu D, Wang Z, Liu D, et al. Independent calibration of water vapor Raman lidar by using additional elastic measurements at water vapor Raman wavelength [C]//EPJ Web of Conferences EDP Sciences, 2016, 119: 25007.
    [15] Deng Qian. Research on self-calibration method of Raman Lidar water vapor detection and development of all-solid-state system[D]. Hefei: University of Science and Technology of China(Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences), 2019. (in Chinese)
  • [1] 罗雄, 史悦, 范琪, 尹微, 彭涛, 赵培娥, 王柯, 周鼎富.  基于相干激光雷达气象多要素探测 . 红外与激光工程, 2023, 52(11): 20230138-1-20230138-10. doi: 10.3788/IRLA20230138
    [2] 王元祖, 孙东松, 韩於利, 郑俊, 赵一鸣.  大气模式对气溶胶光学参量反演及分类的影响研究 . 红外与激光工程, 2023, 52(1): 20220262-1-20220262-13. doi: 10.3788/IRLA20220262
    [3] 杨彬, 莫祖斯, 刘海姣, 卜令兵.  大气探测激光雷达突变信号处理方法研究(特邀) . 红外与激光工程, 2022, 51(1): 20211117-1-20211117-9. doi: 10.3788/IRLA20211117
    [4] 狄慧鸽, 华灯鑫.  底层大气探测激光雷达国内研究现状与进展(特邀) . 红外与激光工程, 2021, 50(3): 20210032-1-20210032-10. doi: 10.3788/IRLA20210032
    [5] 李珂, 刘秉义, 杨倩, 唐军武, 吴松华.  海洋激光雷达测量水体剖面偏振信号的仿真模拟 . 红外与激光工程, 2021, 50(6): 20211035-1-20211035-10. doi: 10.3788/IRLA20211035
    [6] 何秦, 郑硕, 秦凯, 胡明玉, 张亦舒.  基于车载激光雷达走航观测的石家庄及周边地区气溶胶空间分布特征 . 红外与激光工程, 2020, 49(S2): 20200353-20200353. doi: 10.3788/IRLA20200353
    [7] 赵曰峰, 高静, 潘杰, 王旭, 张玉容, 李辉, 王艳琪, 段孟君, 岳伟伟, 蔡阳健, 许化强, 王晶晶.  基于激光雷达的区域大气颗粒物探测 . 红外与激光工程, 2020, 49(S2): 20200324-20200324. doi: 10.3788/IRLA20200324
    [8] 陈鹏, 赵继广, 宋一铄, 王燊.  气溶胶环境下FMCW与脉冲激光探测性能对比 . 红外与激光工程, 2020, 49(6): 20190399-1-20190399-10. doi: 10.3788/IRLA20190399
    [9] 朱倩, 潘增新, 毛飞跃, 石瑞星, 臧琳, 卢昕.  中国区域CALIPSO和MERRA-2气溶胶三维参数对比验证 . 红外与激光工程, 2020, 49(S2): 20200350-20200350. doi: 10.3788/IRLA20200350
    [10] 滕曼, 庄鹏, 张站业, 李路, 姚雅伟.  大气气溶胶污染监测中应用的新型全天时户外型拉曼-米散射激光雷达系统 . 红外与激光工程, 2019, 48(7): 706001-0706001(7). doi: 10.3788/IRLA201948.0706001
    [11] 张青松, 侯再红, 谢晨波.  户外型探测臭氧和气溶胶激光雷达系统研制 . 红外与激光工程, 2019, 48(7): 706008-0706008(6). doi: 10.3788/IRLA201948.0706008
    [12] 陶宗明, 施奇兵, 谢晨波, 刘东, 张帅.  利用CCD和后向散射激光雷达精确探测近地面气溶胶消光系数廓线 . 红外与激光工程, 2019, 48(S1): 43-49. doi: 10.3788/IRLA201948.S106007
    [13] 陈超, 王章军, 宋小全, 张涛, 杜立彬, 孟祥谦, 刘兴涛, 李先欣, 李辉, 庄全风, 王秀芬.  扫描式气溶胶激光雷达研制与观测研究 . 红外与激光工程, 2018, 47(12): 1230009-1230009(7). doi: 10.3788/IRLA201847.1230009
    [14] 邓迁, 吴德成, 况志强, 刘东, 谢晨波, 王英俭.  用于水汽混合比自标定的532 nm/660 nm双波长激光雷达 . 红外与激光工程, 2018, 47(12): 1230004-1230004(5). doi: 10.3788/IRLA201847.1230004
    [15] 尹凯欣, 王海涛, 范承玉.  黑碳-氯化钠内混合气溶胶粒子有效吸收研究 . 红外与激光工程, 2017, 46(4): 411002-0411002(6). doi: 10.3788/IRLA201746.0411002
    [16] 陶宗明, 单会会, 张辉, 张连庆, 王申浩, 麻晓敏, 周浦城, 姚翎, 薛模根, 王邦新, 谢晨波, 刘东, 王英俭.  单波长发射五通道接收激光雷达系统研制 . 红外与激光工程, 2017, 46(10): 1030002-1030002(7). doi: 10.3788/IRLA201765.1030002
    [17] 刘秉义, 庄全风, 秦胜光, 吴松华, 刘金涛.  基于高光谱分辨率激光雷达的气溶胶分类方法研究 . 红外与激光工程, 2017, 46(4): 411001-0411001(13). doi: 10.3788/IRLA201746.0411001
    [18] 尚震, 谢晨波, 钟志庆, 王邦新, 王珍珠, 赵明, 谭敏, 刘东, 王英俭.  用于测量流层水汽的拉曼激光雷达 . 红外与激光工程, 2016, 45(12): 1211003-1211003(6). doi: 10.3788/IRLA201645.1211003
    [19] 孙贤明, 万隆, 王海华.  激光雷达探测水云退偏振比的敏感性研究 . 红外与激光工程, 2016, 45(9): 906001-0906001(5). doi: 10.3788/IRLA201645.0906001
    [20] 闫顺生.  提高拉曼探测气溶胶消光系数精度的玻-温模式 . 红外与激光工程, 2014, 43(9): 3015-3019.
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  304
  • HTML全文浏览量:  51
  • PDF下载量:  47
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-07-11
  • 修回日期:  2022-09-06
  • 网络出版日期:  2023-04-14
  • 刊出日期:  2023-04-25

探测气溶胶-水汽的拉曼-米散射激光雷达系统

doi: 10.3788/IRLA20220484
    作者简介:

    李路,男,讲师,博士,主要从事大气激光雷达技术方面的研究

    通讯作者: 邢昆明,男,工程师,硕士,主要从事大气激光雷达光机结构优化设计方面的研究
基金项目:  安徽省教育厅自然科学研究重点项目(KJ2021 A0945);皖西学院高层次人才项目(WGKQ2022013,WGKQ2022015);中国科学院合肥物质科学研究院“十三五”规划重点支持项目(KP-2019-05);中国科学院战略性先导科技专项(XDA17040524);2019年安徽省科技重大专项(201903c08020013)
  • 中图分类号: TN958.98

摘要: 设计和构建了发射波长为355 nm和532 nm的户外型全天时激光雷达系统,用于探测大气气溶胶和水汽。运用355 nm和532 nm的米散射、532 nm的偏振、氮气和水汽分子的拉曼激光雷达技术,用于对边界层结构、对流层气溶胶和云光学特性及其形态、水汽混合比进行连续探测研究。该系统结构紧凑,运输方便,具备远程操作、数据传输、一键式启动等功能。利用该系统对大气气溶胶和水汽进行探测,探测结果表明:在大气气溶胶的探测过程中,在重污染条件下混合层高度较干净天低,在0.5 km以下,而干净天在1 km左右;通过对消光系数、Angstrom指数和退偏振比分析可知,重污染条件下,底层大气气溶胶以球形粗粒子污染物为主,干净天底层大气气溶胶以球形细粒子污染物为主;在云层中,Ang-strom指数明显减小,且出现负值,说明云粒子半径较大。在水汽探测过程中,采用自标定方法获得系统的标定常数为121,与已标定的激光雷达系统对比,误差在±0.3 g/kg以内;连续探测结果表明可对夜晚5 km及白天混合层以内进行探测。该系统满足产品化的需求,可广泛运用于大气环境的监测领域中。

English Abstract

    • 随着社会经济的不断发展,环境问题日益严重,获得大气参数(如气溶胶、水汽、温度等)越来越受到环保和气象部门的重视。开展大气气溶胶和水汽的连续探测研究,可以提高天气预报的准确性,对研究大气污染气溶胶特性及扩散机理、云的形成、降雨等具有重要的科学意义。激光雷达作为一种主动遥感探测工具,具备传统大气探测方法无法比拟的时空分辨率及连续性,在大气参数测量中发挥着重要作用[1-4]

      2004年,谢晨波等[5]研制的国内首台车载式大气探测激光雷达系统,用于探测大气气溶胶和水汽,出射波长为355 nm (300 mJ)和532 nm (150 mJ)。2010年, 吕炜煜等[6]研制的车载激光雷达系统,用于探测对流层气溶胶和水汽,出射波长为355 nm(400 mJ)和1064 nm(400 mJ)。2014年,谭敏等[7-8]研制的温度、水汽和气溶胶拉曼激光雷达,用于探测大气温度、水汽和气溶胶,出射波长分别为355 nm (500 mJ)、532 nm (200 mJ)和1064 nm (300 mJ)。通过对以上激光雷达系统的分析可知,利用出射波长355 nm的氮气拉曼散射回波信号和水汽拉曼散射回波信号探测水汽混合比,能量在300 mJ以上;355、532、1064 nm的米散射用于探测大气气溶胶,能量在150 mJ以上。这种采用大能量的出射波长,可以有效提高白天的信噪比,实现水汽的白天探测,但所设计的激光雷达系统具有成本高、尺寸大、结构复杂等缺点,很难实现其市场化。

      随着激光雷达技术的不断成熟,激光雷达向着小型化、产品化及简单化方向发展,为了满足国家环保部门及时了解大气参数的测量要求,中国科学院安徽光学精密机械研究所中国科学院大气光学重点实验室研制了一套户外型全天时大气气溶胶-水汽激光雷达产品,用于气溶胶及水汽的长期探测。

    • 大气气溶胶-水汽激光雷达中三波长米散射激光雷达方程为:

      $$ \begin{split} \left[P(r,\;\lambda )-{P}_{b{g}'}(r,\;\lambda )\right]{r}^{2}=C\left[{\beta }_{{a}}(r,\;\lambda )+{\beta }_{m}(r,\;\lambda )\right]\cdot\\\mathrm{exp} \left \{-{{2}{\displaystyle {\int }_{{ 0}}^{{r}}\left[{\alpha }_{a}\text{(}{r}',\;\lambda \text{)}+{\alpha }_{m}{(}{r}',\;\lambda {)}\right]}\text{d}{r}'}\right \} \end{split}$$ (1)

      式中:$ P(r,\lambda ) $为雷达接收的高度为${{r}}$处的后向散射回波功率,km;${P_{{{bg'}}}}(r,\lambda )$为雷达接收的背景信号;$ \lambda $为激光波长,nm;C为激光雷达系统的常数;$\;{\beta }_{m}{(r,}\lambda )$${\alpha }_{m}{(r,}\lambda )$分别为空气分子的后向散射系数和消光系数,由大气模式计算获得;$\;{\beta }_{{a}}{(r,}\lambda )$${\alpha }_{a}{(r,}\lambda )$分别为大气气溶胶的后向散射系数,km−1·sr−1和消光系数,km−1

      反演大气气溶胶消光系数采用Fernald方法,该方法在反演过程中需要假设激光雷达比Sa(气溶胶消光后向散射比),受到气溶胶特性和波长的影响。文中假定激光雷达比Sa为常数,取532 nm和355 nm的Sa为50 sr。空气分子激光雷达比Sm(空气分子消光后向散射比)由大气模式获得,为$ {\text{8}}\pi {\text{/3}} $ sr[9-11]

      在利用Fernald方法反演的过程中,需要选取不含气溶胶粒子的洁净大气的高度作为标定高度rc,取${\beta _{{a}}}({r_c})/{\beta _m}({r_c})$为0.01。但激光雷达在实际探测过程中,由于受到污染或在有云的天气下,实际探测距离达不到干净的大气高度,这样选取的边界值会严重影响反演精度,特别对两波长的反演信号后续分析产生影响,边界高度的边界值的选择非常重要。根据文献[12],文中采用公式(2)作为Fernald方法反演边界值的判据,减小反演误差,利用迭代法可以找出准确的边界值。

      $$ \dfrac{\left|\dfrac{\left[P({r}_{c})-{P}_{b{g}^{\prime }}(r)\right]}{{\alpha }_{a}({r}_{{c}})+{\alpha }_{{m}}({r}_{c})}{r}_{c}{}^{2}\left\{2\tau \left[{\alpha }_{a}(r{}_{c})+{\alpha }_{m}({r}_{c})\right]+\dfrac{{\left\{2\tau \left[{\alpha }_{{a}}({r}_{c})+{\alpha }_{{m}}({r}_{c})\right]\right\}}^{2}}{2!}\right\}-2{\displaystyle {\int }_{{ r}_{0}}^{{r}_{c}}\left[P(r)-{P}_{b{g}^{\prime }}(r)\right]}{r}^{2}{\rm{d}}r\right|}{2{\displaystyle {\int }_{{ r}_{0}}^{{r}_{c}}\left[P(r)-{P}_{b{g}^{\prime }}(r)\right]}{r}^{2}{\rm{d}}r} \leqslant 5{\text{%}} $$ (2)

      其中,

      $$ \tau ({r}_{c},\lambda )={\displaystyle {\int }_{{ r}_{0}}^{{r}_{c}}\left[{\alpha }_{a}({r}',\lambda )+{\alpha }_{m}({r}',\lambda )\right]}{\rm{d}}{r}' $$ (3)

      式中:$ {r_c} $为边界点的高度;$ {r_0} $为激光雷达的回波信号全接收的最低高度。

      混合层高度是指大气气溶胶浓度梯度变化最大的高度,记为$ {h_{GM}} $,其表达式为:

      $$ {h_{GM}} = \min \frac{{{\rm{d}}[P(r){r^2}]}}{{{\rm{d}}r}} $$ (4)

      Angstrom指数A是度量气溶胶粒子尺度的一个重要参数,其表达式为:

      $$ {A_{{\lambda _1}/{\lambda _2}}} = {{ - }}\frac{{{\text{ln}}({\beta _{{{a}},{\lambda _1}}}(r)/{\beta _{{{a}},{\lambda _2}}}(r))}}{{\ln ({\lambda _1}/{\lambda _2})}} $$ (5)

      当细粒子模态(50~500 nm)和粗粒子模态(大于500 nm)共同影响光散射时,利用两个波长指数可以更深入地研究气溶胶粒子的半径[10]。文中研制的雷达系统有355 nm和532 nm两种波长,利用A532/355分析不同天气下的气溶胶粒子特性。

      偏振激光雷达的偏振特性可用于研究气溶胶粒子的形态特征,是探测冰晶云和沙尘气溶胶等非球形粒子的最有效的方法,文中系统设计的波长532 nm的退偏振比可以表示为:

      $$ \delta ({r},\lambda )=k(\lambda )\frac{{P}_{rs}(r)}{{P}_{rp}(r)} $$ (6)

      式中:$ {P_{rs}}(r) $$ {P_{rp}}(r) $分别为雷达接收的回波功率的平行分量和垂直分量;$ k(\lambda ) = {k_p}(\lambda )/{k_s}(\lambda ) $,为两通道的增益比,即标定因子,可通过实验获得。

      水汽拉曼激光雷达探测的水汽混合比(g/kg)可表达为:

      $$ w(r)={C}_{WV}\frac{T({\lambda }_{R,N},{r}_{0},r){P}_{R,WV}(r,{\lambda }_{R,WV})}{T({\lambda }_{R,WV},{r}_{0},r){P}_{R,N}(r,{\lambda }_{R,WV})} $$ (7)
      $$ \left\{\begin{array}{l}T({\lambda }_{R,N},{r}_{0},r)=\mathrm{exp}\left\{-{\displaystyle {\int }_{{ r}_{0}}^{r}\left[{\alpha }_{m}({r}^{\prime },{\lambda }_{R,N})+{\alpha }_{{a}}({r}^{\prime },{\lambda }_{R,N})\right]{\rm{d}}{r}^{\prime }}\right\}\\ T({\lambda }_{R,WV},{r}_{0},r)=\mathrm{exp}\left\{-{\displaystyle {\int }_{{ r}_{0}}^{r}\left[{\alpha }_{m}({r}^{\prime },{\lambda }_{R,WV})+{\alpha }_{{a}}({r}^{\prime },{\lambda }_{R,WV})\right]{\rm{d}}{r}^{\prime }}\right\}\end{array} \right.$$ (8)

      式中:$ {\lambda _{R,N}} $$ {\lambda _{R,WV}} $分别为N2分子和H2O分子的拉曼波长;${P}_{R,WV}(r,{\lambda }_{R,WV})$${P_{{{R}},N}}(r,{\lambda _{R,N}})$分别为高度为$ r $处的H2O分子和N2分子的拉曼回波信号;$ T({\lambda _{R,N}},{r_0},r) $/$ T({\lambda _{R,WV}},{r_0},r) $为大气透过率修正项,为雷达系统中高度$ {r_0} $$ r $处的散射体之间的大气在波长$ {\lambda _{R,N}} $$ {\lambda _{R,WV}} $的透过率的比值,可通过消光模式计算得到;$ {C_{WV}} $为水汽混合比的标定常数[13-15]

    • 为了满足在户外长期探测大气气溶胶-水汽的时空变化,设计及构建了一台户外型全天时激光雷达系统,该系统示意图及外观如图1所示,其外形尺寸为1.8 m×1.2 m×2 m (长×宽×高);表1为该系统主要性能参数,其主要由发射单元、接收单元、后继单元、数据采集单元及辅助控制单元组成。

      图  1  大气气溶胶-水汽激光雷达产品结构示意图

      Figure 1.  Structural schematic of lidar system for atmospheric aerosol-water vapor

      表 1  大气气溶胶-水汽激光雷达系统的主要参数、性能指标

      Table 1.  Main specifications, performance parameters of lidar system for atmospheric aerosol-water vapor

      Parameters
      Wavelength/nm532/355
      Frequency/Hz10
      Monopulse energy/mJ88/50
      Telescope diameter/mm400
      Field of view/mrad0.5-2(adjustable)
      Detection channel5 (355 nm, 532 nm S, 532 nm P, 386 nm and 407 nm)
      Filter bandwidth/nm1
      Capture cardA/D and photon counting
      Performance
      Measuring distance
      Spatial resolution
      Time resolution
      Operation temperature
      Storage temperature
      Way of working
      Power supply
      Volume
      0-12 km(atmosphere of clear sky)
      7.5 m
      15 min(adjustable depending on the measurement content)
      5- −40 ℃(−10- −40 ℃ under the condition of peripheral temperature control)
      5-−70 ℃(−30- −70 ℃ under the condition of peripheral temperature control)
      Manual, automatic and network control
      220 V/50 Hz/1500 W
      1.8 m×1.2 m×2 m(L×W×H

      发射单元中,激光器采用Nd:YAG激光器,输出波长为二倍频532 nm和三倍频355 nm,其能量分别为88 mJ和50 mJ,重复频率均为10 Hz。由于两波长从一个出光口射出,故采用高稳定的光学调整架及分束镜将两波长分离开,再利用反射镜和高稳定光学调整架将两波长反射到大气中,来调节发射光轴与接收光轴并保持高度的平行。由于出射光束小于13 mm,故所有的发射单元镜片采用1 in (1 in=2.54 cm)的圆形镜片。

      接收单元主要由望远镜、小孔光阑及目镜组成。望远镜采用卡塞格林望远镜,有效口径为400 mm,并在其系统焦点处安装具有可调节的小孔光阑,来实现在0.5~2 mrad范围内调节,适用于不同环境和弱信号探测的要求;通过小孔光阑的汇聚光经过焦距为100 mm的目镜将其准直为准平行光用于后继光学单元探测。

      后继光学单元光路设计优化是提高探测灵敏度的关键技术之一,其稳定性直接关系到探测信号是否完全被探测器所接收。首先对后继单元设计进行优化的原则:弱信号光程要短;分束片使用高性能的镀膜,增加反射和透射效率;通过分束片尽可能采用反射方式等,设计出的后继光路如图1所示;然后对后继单元进行结构设计,将所有分束片采用矩形,安装结构采用三角结构,基频远大于200 Hz,结构稳定性高,如图2(a)所示,并通过高精度的插槽方式固定于一体化的后继光学结构中,如图2(b)所示,其基频远大于200 Hz,稳定性高。

      图  2  分束片安装结构

      Figure 2.  Beam splitter mounting structure

      后继探测单元由5个探测通道组成,其中3个米散射探测通道分别用于探测355 nm出射波长与气溶胶粒子相互作用产生的米散射信号及532 nm的米散射型号垂直分量和水平分量,2个拉曼探测通道分别用于355 nm出射波长与氮气分子相互作用产生的386 nm拉曼散射信号和与水汽分子相互作用产生的407 nm拉曼信号。均采用探测器PMT,在每个探测器前安装1 nm的干涉滤光片来滤除背景光;以及安装会聚透镜,安装完成如图3所示。

      图  3  后继探测单元

      Figure 3.  Subsequent detection unit

      数据采集单元主要用于采集探测器转换的电信号,其采集系统采用模数转换A/D (采样精度为12 bits)和光子技术PC (最高计数率为250 MHz)相结合的数据采集器,空间分辨率为7.5 m。

      系统的辅助控制单元主要有以下几种功能:

      (1)利用工控机来确保发射信号与数据采集同步进行;

      (2)利用全自动反演软件,无需人为干预便可获得反演后的数据,如不同波长的距离订正回波信号、退偏振比以及水汽等数据获取,可快速直观地描述大气垂直分布的时空变化特征,同时,与自动测量程序相结合,可实现无人干预条件下的自动测量和数据反演;

      (3)利用设计控制软件可实现对雷达系统方舱进行温度、湿度的自动控制、一键启动、远程数据传输及应急保护等功能。

    • 2019年年初,合肥出现了严重的污染天,运用拉曼-米散射激光雷达对2019年1月份合肥地区上空大气进行了探测,一次探测累积时间为15 min。

    • 测量时间为2019年1月19日23:00,根据气象数据可知,此时AQI为237,PM2.5为187 μg/m3,为重度污染天。图4为实际探测的355 nm和532 nm (P和S)波长的距离平方信号(RSCS)(实线),及其导数(虚线),两波长回波信号随高度变化趋势相同。图5为根据图4的数据反演出的消光系数、Angstrom指数特性及532 nm的退偏振比。由于当时天气为重度污染,且出现低层薄云,影响了系统探测距离。由图4图5可以明显看出,分别在高度1.75 km、3 km以及3.6 km处出现薄云。

      图  4  重度污染混合层高度

      Figure 4.  Mixed-layer height in the heavy pollution conditions

      图4可知,重度污染天的混合层高度很低,对两波长回波信号进行一阶求导,计算出混合层高度分别为0.47 km (355 nm),0.36 km (532 nm),这是由于两波长所受气溶胶分布特征影响不同,从而导致两波长计算处理的混合层高度不同。

      图5可知,重度污染天消光系数较大,两波长的消光曲线具有良好的一致性。利用消光系数计算出的Angstrom指数在−2.5~3.8之间变化,退偏振比在0~0.11之间变化。1.25 km以下的污染层:Angstrom指数在0.8~1.8之间变化,说明在底层气溶胶存在较多的直径为1 μm的粗粒子,且在0.5~1.25 km基本保持不变,说明该区域气溶胶粒子均匀混合。退偏振比在0.03~0.11之间变化,随着高度不断增加而减小,但在0.7~1.25 km发生突变,峰值达到0.07,整个污染层主要以球形粒子污染物为主。在高度1.75 km、3 km以及3.6 km处云层中,两波长消光系数进入云层后数值大小发生反转,对应计算出的Angstrom指数出现较大的波动,并出现负值,说明在云层中粒子半径较大;对应的退偏振比都在0.05以下,属于水云。

      图  5  重度污染天的消光系数(a)、Angstrom指数(b)及退偏振比(c)

      Figure 5.  Extinction coefficient (a), Angstrom index (b) and depolarization ratio (c) in the heavy pollution conditions

    • 测量时间为2019年1月26日13:00,根据气象数据可知,此时AQI为24,PM2.5为14 μg/m3,为一级优天。利用激光雷达系统探测到的RSCS及其导数,如图6所示。根据探测数据反演得到消光系数、Angstrom指数特性及532 nm的退偏振比,如图7所示。由图6图7可知,在8 km处出现卷云。

      图6可知,355 nm和532 nm信号得到的混合层高度为1.1 km,这是由于在一级优天条件下,在混合层内气溶胶分布特征对两波长影响近似。

      图7可知,在一级优天的条件下,两波长在底层反演得到的消光系数远小于重度污染天的消光系数。根据Angstrom指数曲线可知,底层气溶胶的Angstrom指数比重度污染天要大,说明在一级优天条件下,底层大气中气溶胶粒子半径比重度污染天的小。根据退偏振比曲线可知,底层气溶胶在这两种天气下的退偏振比都在同一范围内(0~0.11),说明两种天气底层气溶胶以球形粒子为主。在1.7 km以下,Angstrom指数随着高度增加而呈现增加的趋势,说明气溶胶粒子半径随高度的增加而呈现减小的趋势。在1.7~3.9 km高度下,Angstrom指数随着高度增加而不断减小,最小值到达0.88,且退偏振比不断增加,最大值达到0.047,说明气溶胶粒子半径随高度的增加而增加,且在3.9 km处气溶胶层以直径为1 μm的粗粒子的球形粗粒子为主。在8 km处,卷云的两波长消光系数数值发生反转;对应的Angstrom指数波动较大,且出现负值,说明在云层中粒子半径较大;对应的退偏振比峰值到达0.2,介于水云和冰晶云之间,属于混合云。

      图  6  干净天混合层高度

      Figure 6.  Mixed-layer height in the clean weather conditions

      图  7  干净天的消光系数(a)、Angstrom指数(b)及退偏振比(c)

      Figure 7.  Extinction coefficient (a), Angstrom index (b) and depolarization ratio (c) in the clean weather conditions

    • 采用文献[15]研制的全固态紫外拉曼激光雷达系统探测的水汽混合比对文中研制的激光雷达进行标定,获得标定常数$ {C_{WV}} $,公式如下:

      $$ {C}_{WV}=\frac{{\displaystyle \sum _{n({r}_{0})}^{n({r}_{1})}\frac{{w}_{b}(r)}{\left[{P}_{R,WV}(r,{\lambda }_{R,WV})/{P}_{R,N}(r,{\lambda }_{R,N})\right]\cdot\left[T({\lambda }_{{R},N},{r}_{0},r)/T({\lambda }_{R,WV},{r}_{0},r)\right]}}}{\left[n(r)-n({r}_{0})+1\right]} $$ (9)

      式中:$ {w_b}(r) $为全固态紫外拉曼激光雷达系统获得的水汽混合比;$ n({r_0}) $$ n({r_1}) $分别为高度$ {r_0} $$ {r_1} $处对应的数据点数。

      利用全固态紫外拉曼激光雷达系统对文中所设计的拉曼-米激光雷达进行系统标定,标定时间为2020年11月20日01:45,根据气象数据可知,此时AQI为64,PM2.5为49 μg/m3,该天空气质量良好。由文献[15]提供的水汽混合比垂直分布如图8(a)(实线)所示,在5 km左右出现比较厚的云,使得5 km以上的数据迅速减小,甚至出现负值,这是由于信号没有穿透云层所致。综合考虑拉曼-米激光雷达的几何重叠因子的影响,文中选择0.5~3.5 km的水汽混合比对拉曼-米激光雷达系统进行标定,获得系统标定常数为121。

      图  8  水汽混合比垂直分布特征及绝对误差廓线图

      Figure 8.  Profile of vertical distribution feature of water vapor mixing ratio and absolute error

      拉曼-米激光雷达通过标定常数反演的水汽混合比,如图8(a)(虚线)所示,与全固态紫外拉曼激光雷达反演的水汽混合比,如图8(a)(实线)进行对比,两者探测结果大体一致,水汽结果能较好地吻合。对比两者的绝对误差,如图8(b)所示,绝大部分的绝对误差在±0.3 g/kg以内,但由于几何重叠因子的影响导致在0.5 km 以下误差较大。

    • 测量时间为2018年11月01日18:00到2018年11月02日16:00,根据气象数据可知,此时AQI平均为60,PM2.5平均为40.5 μg/m3,该天空气质量良好。图9给出了探测期间水汽混合比的连续时空分布,垂直分辨率为7.5 m,时间分辨率为15 min。

      图  9  水汽混合比时空分布

      Figure 9.  Spatiotemporal distribution of water vapor mixing ratio

      图9可知,夜晚测量探测高度可达5 km,白天由于背景太强,导致水汽混合比探测高度只有1.4~2 km。在夜晚探测期间,在2.5 km左右有明显的水汽层,随着时间的变化,峰值从1.5 g/kg增加到3 g/kg,由于信噪比低,该水汽层在白天探测不到。水汽主要分布在1.5 km以内,夜晚底层水汽混合比大部分在4~5 g/kg之间变化,白天水汽混合比大部分在3~4 g/kg之间变化,且夜晚底层水汽混合比要比白天高。由此表明,拉曼-米激光雷达系统可以实时、连续、精确探测大气中的水汽含量。

    • 研制了一款用于连续探测大气气溶胶和水汽的激光雷达系统。不同于传统的大气气溶胶和水汽激光雷达,该系统同时具备355 nm/532 nm波长米散射探测、532 nm退偏振探测、氮气和水汽分子拉曼探测的功能。两波长米散射探测功能可以精细探测大气边界层结构、大气气溶胶和云的消光特性及粗细粒子分布;532 nm退偏振探测功能可反映出大气气溶胶粒子和云粒子的形状特征,可识别出球形粒子(水云、污染型气溶胶和雾霾)和非球形粒子(沙尘与冰晶云等);氮气和水汽分子拉曼探测功能可获得水汽混合比的时空分布特征。该系统后继探测单元采用高稳定性的一体化结构,配置方舱具备防尘防水的恒温功能,可以直接在户外长期探测,有利于对局部气溶胶粒子和云粒子、水汽等物理参数间进行统计分析,已用于大气环境监测和科学等研究领域。

参考文献 (15)

目录

    /

    返回文章
    返回