[发明专利]一种多波长激光雷达气溶胶粒子谱分布反演方法及系统

说明书

本发明公开了一种多波长激光雷达气溶胶粒子谱分布反演方法及系统，该方法包括：获取光学数据和正则参数；光学数据包括气溶胶消光系数和气溶胶后向散射系数；确定核矩阵；根据核矩阵、光学数据以及正则参数，计算权系数；根据光学数据、核矩阵、权系数，计算评估误差；确定评估集合；判断评估误差是否在评估集合内，若否则调整正则参数，返回计算权系数步骤，若是则根据权系数，反演确定大气气溶胶粒子谱分布。与现有技术相比，本发明通过采用NCRE算法选取合适的正则参数，实现大气气溶胶粒子谱分布精确反演，且计算简单，适用范围广。

技术领域

本发明涉及激光雷达数据技术领域，特别涉及一种多波长激光雷达气溶胶粒子谱分布反演方法及系统。

背景技术

气溶胶是大气的成分之一，它对气候变化、环境污染、人类健康等有重要的影响而成为全球关注的焦点。特别近些年来，在全球各地越来越频繁发生的沙尘、雾霾现象，使得气溶胶的物理化学以及光学属性受到越来越多的关注。气溶胶消光系数、粒子谱分布等是气溶胶的基本光学和物理特性，它们对于研究气溶胶的辐射迫力和气溶胶对气候变化的影响有非常重要的作用。此外，气溶胶粒子谱分布也是重要的环境质量标准之一，气溶胶对人体健康的影响与气溶胶颗粒的粒径大小有关。因此，精确测量气溶胶粒子谱分布及其时空变化十分重要。

现阶段对气溶粒子谱观测分析的方法主要有直接采样分析法和遥感法。直接采样分析法主要采用粒子谱仪和光谱仪等直接获取气溶胶粒子谱相关数据并进行分析，但这种方法局限于点测量，无法实现大范围空间上的连续探测。遥感法采用太阳光度计、激光雷达等遥感探测工具对大气气溶胶进行探测，并使用专门的算法对遥感数据进行处理，从而获取气溶胶粒子谱分布的信息。其中，激光雷达作为主动遥感探测工具，可连续探测气溶胶高分辨率的时空分布，具有明显的优势。特别是多波长激光雷达技术，可以同时获取多个波长的激光雷达信号与气溶胶粒子相互作用的信息，对气溶胶粒子谱分布等进行探测，广受关注。

多波长激光雷达气溶胶粒径分布探测，是通过探测多个波长的激光与大气分子和气溶胶粒子相互作用产生的后向散射光信号，反演气溶胶消光系数和后向散射系数，进而利用获得的光学数据，采用相关的方法求解第一类Fredholm积分方程，从而得到气溶胶粒子谱分布。第一类Fredholm积分方程表示了气溶胶物理特性与光学特性之间的关系，但该方程不能直接进行求解，需应用一些特别的数学方法进行处理，否则会产生不适定反演问题。因此，采用何种方法求解第一类Fredholm积分方程，是实现气溶胶粒子谱分布高精度反演的先决条件。

目前，求解第一类Fredholm积分方程常用带约束条件的正则法进行求解，该方法不需要预先知道解的形式，可以采用与不适定问题相邻近的适定问题的解去逼近不适定问题的解，其反演的能力已经在相关研究中证明。但如何选取正则参数一直是亟待解决的技术问题，因为过大的正则参数会使反演结果过度平滑，丢失有用信息，而过小的正则参数会使反演结果产生振荡，去噪不充分。现有方法利用GCV(Generalized Cross-Validation)法确定正则参数，即选取GCV参数取最小值时的拉格朗日乘子作为正则参数。当拉格朗日乘子取值较小的时候，GCV参数随拉格朗日乘子的增大而减小，并出现强烈的抖动，说明GCV参数值不稳定；当拉格朗日乘子取某一范围值时，抖动消失，并且当拉格朗日乘子取全局最小值时，GCV参数值达到稳定；当拉格朗日乘子取值大于全局最小值时，GCV参数继续增加。该方法需计算GCV参数，并且需分析参数值随拉格朗日乘子的变化过程，因而存在计算复杂、只适合某些特定场合等不足。

发明内容

本发明的目的是提供一种多波长激光雷达气溶胶粒子谱分布反演方法及系统，通过采用NCRE算法选取合适的正则参数，实现大气气溶胶粒子谱分布精确反演，且计算简单，适用范围广。

为实现上述目的，本发明提供了如下方案：

一种多波长激光雷达气溶胶粒子谱分布反演方法，所述多波长激光雷达气溶胶粒子谱分布反演方法包括：