激光雷达测量大气气溶胶的数据处理研究

赵云珊 陈鹏举 杨 昭 李 强

（连云港师范高等专科学校物理与电子工程系 江苏 连云港 222006）

0 引言

气溶胶粒子是指悬浮在大气中的直径为0.001-10μm的固体或液体粒子，其质量仅占整个大气质量的十亿分之一，但却对大气中的气候变化、云量及云的寿命、光学特性的改变等有着十分重要的影响，通过气溶胶表面的化学反应能影响大气的组成成份，进而影响全球的气候环境。由于目前缺乏全球范围内空间和时间分布资料，大气气溶胶和云的具体作用机制还存在很多不确定性，这大大限制了人类对气候系统及全球气候变化趋势的认识。

激光雷达是一种主动式光学遥感设备，它以激光为光源，通过探测激光与大气相互作用的辐射信号来遥感大气。由于使用探测束的波长较短且定向性强，激光雷达具有很高的空间、时间分辨能力及探测灵敏度，目前已经成为探测大气气溶胶物理-光学性质的一个重要手段。通过激光光波与大气的相互作用，会产生包含气体原子、分子、大气气溶胶粒子和云等有关信息的辐射信号，利用相应的反演方法则可以获得气体原子、分子、大气气溶胶粒子和云等大气成份的信息。

1 激光雷达系统

表1 激光雷达系统光学结构参数Table 1 Parameters of mie-scattering lidar system

本系统激光发射与接收单元不同轴，导致在一定范围内发射激光束只能逐渐进入接收望远镜的视场，因此接收信号强度必须受系统重叠因子 Y（r）的修正[3]。

根据激光雷达方程，接收的信号光子数为[4]

其中，η是探测器量子效率；λ是激光波长（nm)；E0是激光发射脉冲能量(J)；h是普朗克常数；c是光速（m/s）；Ar是望远镜有效接收面积（m2）；Y（r）是重叠因子；Tsys是系统光学总透过率；Δr是激光雷达距离分辨率(m)；βb（r）是距离 处大气后向散射系数(km-1Sr-1)；σr是距离 r处的大气消光系数(km-1)。对于确定的激光雷达系统及其参数，在水平激光发射和均匀大气情况下，（1）式可以简化为

一般远场激光束波面为高斯分布。根据表1中的参数可得出接收相对光子数的分布曲线。但是，经过40倍率扩束镜的出射光束发散角不一定压缩了40倍，需要调整扩束镜筒长才能达到最佳效果。过大的发散角会对Y（r）及相对接收光子数的分布产生影响，如图2和3所示。图中可以看出，随着发散角的增大，重叠因子及接收相对光子数分布形状变化也较大。当扩束倍率为10时，其相对误差达到60%以上。

图2 不同扩束倍率下的重叠因子

2 数据处理方法

在过去的几十年，已经出现了多种利用激光雷达回波信号反演粒子消光系数和后向散射系数的方法。常用的方法主要包括斜率法，Klett方法[5]，Fernald 方法[6]，以及线性迭代方法[7]等。

2.1 斜率法

由激光雷达方程可得：

假设某个区域内（rtop＞r＞rbottom）大气均匀分布，即后向散射系数和消光均为常数，对(3)式两边取对数，则激光雷达方程变为线性回归方程，假定从上到下对信号进行积分：

甲状腺乳头状癌为发病率较高的内分泌系统恶性肿瘤，分化好，恶性程度较低。其发病与性别、年龄、碘的摄入等多因素有显著的相关性[13-15]。随着对其研究的不断推进，我国甲状腺癌的诊治逐渐规范，甲状腺癌诊断方法由单纯使用超声发展为影像定位-细胞病理-分子靶标多层次早期诊断体系，诊断过程较为系统，具有科学性[16]。TPC恶性程度较低，主要治疗方案以外科手术为主。孟晓敏等对甲状腺乳头状癌原发灶外科治疗进行meta分析得出：TPC原发灶外科治疗术后生存率较为可观，且甲状腺部分切除术损伤喉返神经或甲状旁腺发生率较低[17]。

2.2 Klett方法

斜率法仅仅适合于均匀大气的探测。对于非均匀大气，σ(r),β(r)不再为定值，斜率法将不再适用。

Klett假定 σ(r),β(r)存在下述关系：

对(4)式进行距离微分，结合(5)式得：

上式为一阶伯努利方程。 假设 X(r)＝ln（β′（r）），则由(6)式得：

由于(7)式分母中的两项之差可以很小甚至为零，故其解很不稳定，经常会产生严重发散的结果。假定j为1，Klett提出了更加稳定的解的形式：

式(7)为前向积分式，σ(r0)为其边界条件，式(8)为后向积分式，σ（rM）为其边界条件。利用(7)、(8)可以求得不同高度处的大气消光系数。边界条件σ(r0)和σ（rM）可以假设某一小段范围内大气均匀利用斜率法计算得到。

2.3 Fernald 方法

当激光在大气中传播时，要受到大气分子和气溶胶粒子（云粒子）两部分同时作用。因此激光雷达方程中σ(r),β(r)实际上包括两部分，即分子贡献部分和气溶胶（云）贡献部分，故有：

对于校准的激光雷达方程，得到的前向反演解为：

后向反演解为：

其中，Sm为大气分子的后向散射比，而Sp为粒子的后向散射比，

其中，rc为参考高度，在此高度上气溶胶散射足够小可以忽略。可通过选取近乎不含气溶胶的清洁大气层所在的高度来确定。

斜率法和Klett方法获得的是包含大气分子和大气气溶胶粒子或云粒子的大气总消光系数，而Fernald方法将大气分子与大气气溶胶或云的消光系数分开考虑来求解激光雷达方程，是激光雷达方程各种反演方法具有代表性也是最常用的一种方法。

2.4 线性迭代法

线性迭代方法考虑了大气分子和大气气溶胶粒子或云粒子两种大气成份，同时可以对后向散射比进行调整。该方法将大气（云/气溶胶）分为等厚的N份，垂直厚度为Δr＝(rt-rb)/N，距离l从1→N。具体迭代公式如下：

迭代过程持续进行，直到β值收敛于一点，迭代过程结束。

3 对流层云和气溶胶后向散射系数反演

本实验中主要利用Fernald方法进行数据处理和后向散射系数的反演工作。下图给出了激光雷达系统在2011年11月12日测量得到的衰减的后向散射系数进行20km水平平均获得的廓线，其中2-4km范围内存在一层气溶胶层，其后向散射系数相对于其他高度明显偏大。

图3 对流层气溶胶衰减的后向散射系数

上面仅给出了单一特征存在时其后向散射系数的反演结果，当多种特征同时存在时，反演算法变得非常复杂，需要首先反演上层特征的后向散射系数、消光系数及其光学厚度，并对激光雷达廓线进行修正，这样才能够准确获得下层特征的光学特性，这也是今后研究的重点工作。

［1］Yang Zhao,Li Qiang.Study about atmosphere extinction coefficient based on 1064nm Mie-scattering lidar[J].Chinese J.Laser Technology,2006，30（2）：170-173.（in Chinese）.

［2］American National Standards Institute,“American national standard for the safe use of lasers”[J].ANSI Z136.1-1986,pp.1-96,1986.

［3］Zhong Zhiqing,Zhou Jun.Calculation of detect signal on micropulse lidar[J].Chinese J.Quantum Electron,2003，20（5）：618-622.（in Chinese）.

［4］Yang Chunhu,Sun Dongsong,Li Hongjing.Photon counting applied to imaging lidar[J].Chinese J.Infrared and Laser Engineering.2005，34（5）：517-520（in Chinese）.

［5］J.D.Klett.Stable Analytical Inversion Solution for Processing Lidar Returns,Appl[J].Opt.,1981,Vol.20:211-220.

［6］F.G.Fernald.Analysis of atmospheric lidar observations:some comments,Appl[J].Opt.,1984,Vol.23:652-653.

［7］C.M.R.Platt,S.A.Young,P.J.Manson et al.The Optical Properties of Equatorial Cirrus from Observations in the ARM Pilot Radiation Observation Experiment,J.[J].Atmos.Sci.,1998,Vol.55:1977-1996.