一种确定消光系数边界值的新算法 下载: 1143次
1 引言
激光雷达是一种将短激光脉冲以一个或多个波长发射到大气中,通过与大气相互作用后测量回波信号的主动式遥感仪器,可以探测气溶胶消光系数、后向散射系数等参数的时空分布。激光雷达已被证明是提供大气边界层(ABL)结构以及连续观测云和气溶胶属性的强大工具之一,也是大气环境的监测手段之一,可以进行连续监测,具有高时空分辨率[1-3]。沈吉等[4]利用拉曼散射激光雷达探测了平流层气溶胶的消光系数廓线,分析了平流层大气气溶胶消光系数的变化特征。项衍等[5]通过讨论气溶胶光学特性反演方法的流程,分析了影响反演结果准确度的各种不确定性因素。常用的大气消光系数反演算法有Collis斜率法[6]、Klett法[7]和Fernald法[8]。当假设大气为均匀分布时,可以使用Collis斜率法对消光系数进行反演。但实际上均匀分布的大气状况不会出现,所以大多采用Klett法或Fernald法求解消光系数。对于这两种算法,不同时间和不同天气条件下消光系数的边界值是一个不确定的值,故精确反演气溶胶消光系数的关键环节就是边界值的选取。1993年,Kovalev[9]提出了一种适用于双组分散射大气的激光雷达信号迭代反演方法,该方法对原始信号进行了变换,使得单组分大气的激光雷达方程求解成为可能。2010年,Marchant等[10]利用迭代最小二乘法从弹性激光雷达数据中得到了有用的信息,该方法能处理任意数量的通道,并在参考点以上和低信噪比的情况下生成稳定的解。2012年,熊兴隆等[11-12]分别通过不动点迭代法和改进的牛顿法求解消光系数的边界值。李红旭等[13]利用弦截法得到了消光系数的初始边界值,然后对边界值进行迭代修正。陈涛等[14]将激光雷达距离校正回波信号在一段范围内的积分值和此范围内的大气透过率作为判断依据,利用穷举法反复迭代确定合适的边界值,但该方法的取值步长固定,迭代速度较慢,运算量大。
本文提出了一种基于横向Steffensen型3阶方法反演消光系数边界值的新方法。首先选取参考点,结合激光雷达方程和Fernald反演算法构建关于消光系数边界值与消光系数均值的非线性方程,然后利用一种横向Steffensen型3阶方法进行求解,最终得到较准确的消光系数边界值。利用本文方法对仿真和实测回波信号进行消光系数的反演,结果表明,本文方法可以得到较准确的消光系数。本文方法不需求导,迭代次数较少,收敛速度可提高到3阶。
2 仪器介绍
本文使用的激光雷达为第二代CAMLTM CE 370-2后向散射激光雷达。对于一个低电平雷达信号的精确探测来说,可通过光子计数系统进行信号的采集,雷达的最大垂直分辨率为15 m。
表 1. CE 370-2系统的主要技术参数指标
Table 1. Main technical parameters of CE 370-2 system
|
3 大气消光系数反演算法
激光雷达发射的532 nm的激光处在可见光范围内,激光束通过大气中的分子和气溶胶时会受到散射和吸收,接收系统收集回波信号并进行处理,以获得所需信息。Fernald反演算法[8]将大气散射分为分子散射和气溶胶散射两部分,因此激光雷达接收到的距离
式中:
由激光的雷达方程可知,大气气溶胶消光系数和后向散射系数是不确定的,所以激光雷达比是一个必要的参数[15]。对于气溶胶,
根据Fernald反演算法,低于边界高度
式中:
4 消光系数边界值的确定方法
4.1 横向Steffensen型三阶方法求解非线性方程
求解形如
在
式中:
令
将(4)式中的
求解(5)式得到新的迭代公式为
其中,
令
将
4.2 消光系数边界值的确定算法
首先构建关于消光系数边界值的非线性方程[19],然后利用横向Steffensen型3阶方法得到迭代公式,求解此方程就可解得消光系数的边界值。算法流程图如
1) 构建非线性方程。首先在激光雷达最大有效探测距离范围内选取
式中:
2) 确定迭代公式。
3) 选取适当的迭代初值。
4) 用(12)式生成迭代序列{
5) 判断迭代是否停止,判断条件为|
迭代过程中当
5 实验验证及分析
采用CAMLTM CE 370-2后向散射激光雷达系统参数得到的仿真信号和实际测得的回波信号对本文方法进行验证。分别使用本文算法(算法一)、文献[
13]中提到的弦截法(算法二)和文献[
11]中的不动点迭代法(算法三)获得大气气溶胶消光系数的边界值,然后将其代入Fernald反演算法中求解大气消光系数,最后对这三种算法得到的结果进行比较。实验中边界高度
图 2. 仿真信号及反演结果。(a)激光雷达仿真信号;(b)算法一的反演结果;(c)算法二的反演结果;(d)算法三的反演结果
Fig. 2. Simulation signals and inversion results. (a) Lidar simulation signal; (b) inversion result of method 1; (c) inversion result of method 2; (d) inversion result of method 3
从
表 2. 三种反演算法的结果比较
Table 2. Inversion results of three inversion algorithms
|
2015年11月10日,当日空气质量为中度污染,空气质量指数为160,PM2.5的质量浓度为120 μg/m3。利用算法一对10日11:16—11:53实际测得的激光雷达回波信号进行反演,得到的消光系数边界值为0.00816 km-1,设置的迭代初值为0.2 km-1,迭代次数为3次;算法二得到的大气气溶胶消光系数边界值为0.00816 km-1,设置的迭代初值为0.2 km-1和0.3 km-1,迭代次数为5次;算法三得到的大气气溶胶消光系数边界值为0.00825 km-1,设置的迭代初值为0.2 km-1,迭代次数为82次。三种算法中设置
图 3. 2015年11月10日实测信号及反演结果。 (a)激光雷达回波信号;(b)算法一的反演结果;(c)算法二的反演结果;(d)算法三的反演结果
Fig. 3. Measured signals and inversion results on November 10, 2015. (a) Lidar echo signal; (b) inversion result of method 1; (c) inversion result of method 2; (d) inversion result of method 3
表 3. 三种反演算法的结果比较
Table 3. Inversion results of three inversion algorithms
|
表 4. 数值迭代结果
Table 4. Results of numerical iteration
|
由
6 结论
本文提出了一种基于横向Steffensen型3阶方法反演大气气溶胶消光系数边界值的新算法。利用横向Steffensen型3阶方法求解消光系数边界值的非线性方程,得到消光系数边界值。利用该方法计算得到的消光系数边界值与设置的迭代初始值无关,不需要求导运算,且能在一定条件下将收敛速度提高到3阶。通过对仿真信号和实际观测得到的回波信号进行反演,证实了本文方法的可靠性。将本文算法与利用弦截法和不动点迭代法反演得到的消光系数曲线进行比较,结果表明:本文算法的
[2] 梅亮. 沙氏大气激光雷达技术及其研究进展[J]. 激光与光电子学进展, 2018, 55(9): 090004.
[3] 虞历尧, 曹念文, 沈吉. Rayleigh-Raman-Mie激光雷达探测南京北郊对流层气溶胶光学特性[J]. 激光与光电子学进展, 2018, 55(11): 110102.
[4] 沈吉, 曹念文. 拉曼散射激光雷达反演平流层气溶胶消光系数廓线[J]. 中国激光, 2018, 45(6): 0609002.
[5] 项衍, 刘建国, 张天舒, 等. 激光雷达探测气溶胶光学特性的不确定性因素研究[J]. 激光与光电子学进展, 2018, 55(9): 092801.
[6] Collis R TH, Russell PB. Lidar measurement of particles and gases by elastic backscattering and differential absorption[M] ∥Hinkley E D. Laser monitoring of the atmosphere. Topics in applied physics. Berlin, Heidelberg: Springer, 1976, 14: 71- 151.
[8] Fernald F G. Analysis of atmospheric lidar observations: some comments[J]. Applied Optics, 1984, 23(5): 652-653.
[9] Kovalev V A. Lidar measurement of the vertical aerosol extinction profiles with range-dependent backscatter-to-extinction ratios[J]. Applied Optics, 1993, 32(30): 6053-6065.
[10] Marchant C C, Moon T K, Gunther J H. An iterative least square approach to elastic-lidar retrievals for well-characterized aerosols[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(5): 2430-2444.
[11] 熊兴隆, 蒋立辉, 冯帅, 等. 基于不动点原理的大气气溶胶消光系数边界值确定方法[J]. 光电子·激光, 2012, 23(2): 303-309.
Xiong X L, Jiang L H, Feng S, et al. Determination of the boundary value of atmospheric aerosol extinction coefficient based on fixed point principle[J]. Journal of Optoelectronics·Laser, 2012, 23(2): 303-309.
[12] 熊兴隆, 蒋立辉, 冯帅, 等. 改进的牛顿法确定大气消光系数边界值[J]. 红外与激光工程, 2012, 41(7): 1744-1749.
Xiong X L, Jiang L H, Feng S, et al. Using improved Newton method to determine the boundary value of atmospheric extinction coefficient[J]. Infrared and Laser Engineering, 2012, 41(7): 1744-1749.
[13] 李红旭, 常建华, 朱玲嬿, 等. 基于微脉冲激光雷达的能见度反演算法[J]. 电子测量与仪器学报, 2017, 31(10): 1555-1560.
Li H X, Chang J H, Zhu L Y, et al. Visibility inversion algorithm based on micro pulse lidar[J]. Journal of Electronic Measurement and Instrumentation, 2017, 31(10): 1555-1560.
[14] 陈涛, 吴德成, 刘博, 等. 低层大气中确定气溶胶后向散射系数边界值的新方法[J]. 光学学报, 2010, 30(6): 1531-1536.
[15] 宋秀瑜, 曹念文, 杨思鹏. 探究影响南京地区大气气溶胶光学特性反演的因素[J]. 激光与光电子学进展, 2017, 54(4): 040101.
[16] Sasano Y. Tropospheric aerosol extinction coefficient profiles derived from scanning lidar measurements over Tsukuba, Japan, from 1990 to 1993[J]. Applied Optics, 1996, 35(24): 4941-4952.
[17] 陈莎莎, 徐青山, 徐赤东, 等. 基于微脉冲激光雷达计算整层大气气溶胶光学厚度[J]. 光学学报, 2017, 37(7): 0701002.
[18] Candela V, Peris R. A class of third order iterative Kurchatov-Steffensen (derivative free) methods for solving nonlinear equations[J]. Applied Mathematics and Computation, 2019, 350: 93-104.
[19] 熊兴隆, 蒋立辉, 冯帅, 等. 非线性方程法确定低空探测机载激光雷达消光系数边界值[J]. 光电子·激光, 2012, 23(7): 1356-1362.
Xiong X L, Jiang L H, Feng S, et al. Constructing and solving the nonlinear equation of airborne lidar for determining the boundary value of the extinction coefficient for atmospheric aerosol in lower atmosphere[J]. Journal of Optoelectronics·Laser, 2012, 23(7): 1356-1362.
[20] 熊兴隆, 冯帅, 蒋立辉, 等. 一种新的大气消光系数边界值确定方法[J]. 光电子·激光, 2011, 22(11): 1699-1705.
Xiong X L, Feng S, Jiang L H, et al. A novel method for determining the boundary value of the atmospheric extinction coefficient[J]. Journal of Optoelectronics·Laser, 2011, 22(11): 1699-1705.
[21] 周军, 岳古明, 戚福第, 等. 大气气溶胶光学特性激光雷达探测[J]. 量子电子学报, 1998, 15(2): 140-148.
Zhou J, Yue G M, Qi F D, et al. Optical properties of aerosol derived from lidar measurements[J]. Chinese Journal of Quantum Electronics, 1998, 15(2): 140-148.
Article Outline
陈晓楠, 毕京平, 王凯欣, 韩冰, 柳云雷. 一种确定消光系数边界值的新算法[J]. 激光与光电子学进展, 2019, 56(24): 240102. Xiaonan Chen, Jingping Bi, Kaixin Wang, Bing Han, Yunlei Liu. Method for Determining Boundary Value of Extinction Coefficient[J]. Laser & Optoelectronics Progress, 2019, 56(24): 240102.