1.一种基于激光雷达测量的气溶胶分类方法,其特征在于,包括:S1.对大气层高度进行分层形成多个高度层列,0‑4km高度每200m一层,4‑8km高度每
500m一层,每层高度假定气溶胶均匀分布;
S2.根据激光雷达系统采集装置的设定采样率确定激光雷达的分辨率,对分割好的每个高度层列的大气雷达信号进行单独存储;
S3.对每个高度层列的激光雷达方程,利用多波长雷达的拉曼通道和米散射通道,联立每个空间分辨率下的激光雷达方程,联立求解得355nm和532nm波长对应位置的后向散射系数β和消光系数α;
S4.利用正则化方法对每个分辨率高度上的粒子谱分布进行反演,利用得到的粒子谱反推激光雷达比,并通过迭代算法将LR误差降低至合理区间之内,得到对应的1064nm波长的后向散射系数;
S5.计算参数LR355、LR532、kβ(355,532)、kβ(355,1064)、kβ(532,1064)、kα(532,1064)、CR(532/355)和δ532,并将这些参数设置为一个8维向量,用该8维向量表示该层大气的气溶胶光学参数;LR355 是355nm波长的激光雷达比,LR532 是532nm波长的激光雷达比,kβ(355,
532)是355nm和532nm波长下的埃指数后向散射比,kβ(355,1064)是355nm和1064nm波长下的埃指数后向散射比,kβ(532,1064)是532nm和1064nm波长下的埃指数后向散射比,kα(532,1064)是532nm和1064nm波长下的埃指数消光比,CR(532/355)是532nm和355nm波长下的色比,δ532 是532nm波长下的退偏比;
S6.选取数据库的标准参量,将六种粒子的标准向量x列出;
S7.运用余弦相似度算法计算观测向量与六种气溶胶标准参量的相似度,并得到最大值;
S8.将该最大值对应的高度层气溶胶种类归类为相似度值最大的气溶胶。
2.根据权利要求1所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,所述拉曼通道对应的氮气拉曼散射激光雷达方程如下所示:;
米散射通道对应的米散射激光雷达方程如下所示:
;
式中, 表示激光雷达系统发射的 波长的激光在 高度处的氮气拉曼散射后向散射的信号的能量大小, 表示激光雷达系统发射的 波长的激光在 高度处的米散射后向散射的信号的能量大小,表示高度,表示米散射通道的系统常数, 表示拉曼通道的系统常数, (r)表示在 高度处的激光氮气分子的拉曼效应的后向散射系数,表示积分当中位于 高度区间内的中间变量,亦指代高度, 表示发射波长为的激光对应 高度的气溶胶消光系数, 表示发射波长为 的激光对应 高度的大气分子消光系数, 表示发射波长 对应的氮气拉曼散射 的激光对应 高度的气溶胶消光系数, 表示发射波长 对应的氮气拉曼散射 的激光对应 高度的大气分子消光系数, 表示发射波长为 的激光对应 高度的气溶胶后向散射系数, 表示发射波长为 的激光对应 高度的大气分子后向散射系数。
3.根据权利要求2所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,对拉曼散射激光雷达方程进行取对数和求导,直接获得气溶胶的消光系数 :;
式中, 表示在 高度处的氮气分子数密度;
根据经验值选择对流层上部高度 处的气溶胶粒子浓度为0,可推导出气溶胶的后向散射系数 :;
联立得到气溶胶的消光系数 和后向散射系数 ,并得到激光雷达比LR, , 为发射波长 高度 对应的LR;
通过上述过程,得到355nm波长的不同高度处的消光系数 、后向散射系数、激光雷达比 和532nm波长的不同高度处的消光系数 、后向散射系数 、激光雷达比 。
4.根据权利要求3所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,
1064nm波长的后向散射系数β使用米散射计算公式进行计算,使用Tikhonov正则化反演粒子谱分布,根据算得的粒子谱分布通过米散射理论计算出新的LR,并将计算得到的新LR同先前假设的LR进行迭代计算,直至误差小于3%即可认为所得LR为理想LR,通过该理想LR计算得到1064nm下的后向散射系数。
5.根据权利要求4所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,设定对流层上层某高度 处无气溶胶,8km以上的空间认为不存在气溶胶粒子,每一层的高度都认为只有一种气溶胶粒子种类,得到28个剖面,每个剖面按照一种气溶胶参数进行计算;
将光学参数进行扩展,包括埃指数消光比、埃指数后向散射比、色比、退偏比:埃指数后向散射比: ;
埃指数消光比: ;
色比: ;
退偏比: ;
上式中, 是 波长的垂直接收通道的后向散射系数, 是 波长的平行接收通道信号的后向散射系数。
6.根据权利要求5所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,将所述8维向量表述为:=( LR355,LR532,kβ(355,532),kβ(355,1064),kβ(532,1064),kα(355,532),CR(532/
355),δ532)。
7.根据权利要求6所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,把每种参考气溶胶的向量 的值表述出来,利用已有的数据对其进行赋值,赋值的气溶胶粒子的分类标准为清洁大陆型、污染大陆型、海洋型、烟雾型、沙尘型、火山灰型六种纯气溶胶粒子。
8.根据权利要求7所述的一种基于激光雷达测量的气溶胶分类方法,其特征在于,通过每个剖面计算得到向量 ,并将 与六种标准气溶胶粒子的 进行距离计算,通过余弦相似度算法计算二者的相似性: ;
其中A、B分别代表的是每个气溶胶的参考向量和每层气溶胶的计算得到的真实向量,通过真实值与六个参考值进行计算,得到相似度最大的组合为目标类气溶胶,认为该层剖面的气溶胶以目标类气溶胶为主。