1.一种基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,包括步骤:S1:输入多波长激光雷达探测气溶胶后得到的光学特性参数,包括后向散射系数和消光系数,对所述光学特性参数进行归一化处理,得到与气溶胶负载无关的N个光学特性参数;
S2:建立气溶胶微物理特性参数与光学特性参数关联的查找表;
S3:计算所述光学特性参数与查找表中各气溶胶个体#j在光学特性域的距离d#j,并仅保留K个具有最小距离的气溶胶个体作为可行解的范围;
S4:根据Monte Carlo随机采样原理,按概率采样生成NMC个搜索执行顺序,在每个搜索执行顺序中,依次计算待反演个体与查找表中可行解的范围内气溶胶个体的距离,排除筛选后,取对于最后一个光学特性量有最小距离的气溶胶个体作为单次求解的备选解,共得到NMC个备选解;
S5:计算各备选解对应的物理特性参数,其中数浓度、表面积浓度和体积浓度的公式如下:式中, 和 分别为第#p个备选解在查找表中对应数浓度为1cm-3时的表面积浓度和体积浓度; 和 分别为第#p个备选解的数浓度、表面积浓度和体积浓度,λi对应不同激光雷达探测波长,n和m分别为后向散射系数和消光系数的探测波长数量,β*和α*分别为输入个例对应的后向散射系数和消光系数,β#p和α#p分别为第#p个备选解在查找表中对应的后向散射系数和消光系数;
平均NMC个备选解对应的物理特性参数,最终反演得到气溶胶微物理特性参数。
2.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,在步骤S1中,按下述方式对光学特性参数进行归一化处理得:Si,j=α(λj)/β(λi)
式中,i=1,...,n,j=1,...,m,β为后向散射系数和α为消光系数,n和m分别为后向散射系数和消光系数的探测波长数量。
3.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,所述步骤S2中,查找表中存储有对应数浓度为1cm-3时的气溶胶个体的多波长后向散射系数、消光系数、归一化后的光学特性参数及对应的气溶胶复折射率、有效粒径、表面积浓度、体积浓度;
后向散射系数:β(λ)=∫Kβ(r,λ,mr+imi)n(r)dr消光系数:α(λ)=∫Kα(r,λ,mr+imi)n(r)dr
有效粒径:reff=∫r3n(r)dr/∫r2n(r)dr
表面积浓度:st=4π∫r2n(r)dr
体积浓度:vt=4π/3∫r3n(r)dr
式中,λ为对应的激光探测波长,r为气溶胶的粒径,n(r)为气溶胶粒子的粒径分布,Kβ和Kα分别为后向散射系数和消光系数的核函数,mr+imi为气溶胶复折射率,i为虚数单位,mr和mi分别为复折射率的实部和虚部。
4.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,所述步骤S3中,距离d#j的计算公式如下:#j * #j *
d =∑G∈Ω|(G-G )/G|
式中,Ω为步骤S1中得到的N个光学特性参数集合,G*表示输入光学特性参数对应的归一化后的光学特性参数,G#j表示步骤S2中所构建查找表中第#j个个体的归一化后的光学特性参数。
5.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,所述步骤S4中,Monte Carlo随机采样的次数NMC小于全排列 且生成1~N的随机排列的采样方法可以有多种。
6.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,所述步骤S4中,所述的排除筛选为:按照距离 由小到大排列,排除距离较大的(1-ω)部分,ω∈(0,1)。
7.如权利要求1所述的基于多波长激光雷达反演气溶胶微物理特性的方法,其特征在于,所述步骤S5中,根据每个备选解#p,直接在所述的查找表获得对应的气溶胶复折射率有效粒径 求平均后得到最终的气溶胶复折射率和有效粒径。