一种基于双波长激光雷达的边界层高度反演方法技术

本发明专利技术公开了一种基于双波长激光雷达的边界层高度反演方法。该方法通过利用双波长激光雷达355nm和532nm的后向散射信号构造大气粒子的特征序列。然后使用层次聚类分类器对特征序列进行分类，得到边界层之上的大气分子类和边界层之下的气溶胶类别序列。最后根据类别序列，通过查找类别序列的最大梯度处，检索出边界层高度的结果。该方法不依赖于气溶胶的垂直浓度廓线来查找边界层高度，而是根据大气粒子的二维分布，通过层次聚类的方法来确定边界层高度，与K均值、高斯混合模型等方法相比，精度提高了60％，避免了复杂气溶胶层次对查找边界层高度带来的影响。本发明专利技术方法具有执行简单，准确度高和适用性广泛的特点。准确度高和适用性广泛的特点。准确度高和适用性广泛的特点。

【技术实现步骤摘要】
一种基于双波长激光雷达的边界层高度反演方法


[0001]本专利技术属于大气探测
，具体涉及一种基于双波长激光雷达的边界层高度反演方法。

技术介绍

[0002]大气边界层是大气的最低层，大气边界层的高度与污染物的积累和扩散直接相关，也是大气化学模型和天气预报系统的重要输入参数。大气边界层高度的准确反演有利于研究污染物的时空分布，并进一步优化污染模拟。因此，准确反演大气边界层高度是非常重要的。
[0003]由于激光雷达具有很高的空间分辨率和连续稳定的工作能力，已成为探测边界层高度的主要手段。激光雷达系统可以使用气溶胶作为示踪剂来检测边界层高度，边界层高度可以从气溶胶浓度的垂直廓线推断。常规激光雷达边界层高度的算法包括梯度法，小波协方差变换法和理想剖面拟合法(Hayden,K.L.,et al.Steyn,D.G.,1997.The vertical chemical and meteorological structure of the boundary layer in the lower fraser valley during pacific 93.Atmos.Environ.31(14),2089
–
2105；Brooks,I.M.,2003.Finding boundary layer top:application of a wavelet covariance transform to lidar backscatter profiles.J.Atmos.Ocean.Technol.20,1092
–
1105；Steyn DG,et al.,The detection of mixed layer depth and entrainment zone thickness from Lidar backscatter profiles.Journal of Atmospheric and Oceanic Technology 1999；16:953
‑
9.)。这些算法都是基于气溶胶浓度的垂直梯度变化来识别边界层高度。这意味着，一旦气溶胶的垂直分布变得不均匀，气溶胶浓度的垂直梯度就会非常混乱，这些算法将难以准确地确定边界层高度(Wang,Z.,Cao,X.,Zhang,L.,Notholt,J.,Zhou,B.,Liu,R.,and Zhang,B.,2012.Lidar measurement of planetary boundary layer height and comparison with microwave profiling radiometer observation,Atmos.Meas.Tech.,5,1965
–
1972.)。之前提出基于粒子特性差异的边界层高度反演方法来克服这个限制，该方法需要基于激光雷达的色比和退偏比信号，通过计算大气粒子之间的粒子特性差异来获取边界层高度，可以有效的避免多层气溶胶的影响(马盈盈等，一种基于粒子特性差异的边界层高度反演方法:中国,108414997.B[P].2021
‑
07
‑
02.)。然而，有研究表明激光雷达观测的色比和退偏比信号存在一定的不确定性(Liu,B.,Ma,Y.,Gong,W.,Zhang,M.2017.Observations of aerosol color ratio and depolarization ratio over Wuhan.Atmospheric Pollution Research,8(6),1113
‑
1122.)，导致该算法不能有效的克服残留层的影响(Liu,B.,Ma,Y.,Gong,W.,Zhang,M.,&Yang,J.2019.Improved two
‑
wavelength Lidar algorithm for retrieving atmospheric boundary layer height.Journal of Quantitative Spectroscopy and Radiative Transfer,224,55
‑
61.)。这就意味着，当垂直气溶胶浓度较为稀薄时，会导致激光雷达反演的边界层结果失真。

技术实现思路

[0004]本专利技术的目的是为了克服
技术介绍
所述的不足，提出一种基于双波长激光雷达的边界层高度反演方法，利用355nm及532nm的后向散射信号以及层次聚类算法进行边界层高度的查找。用以克服色比和退偏比信号的不稳定性，改进激光雷达算法反演边界层高度不准确的问题，实现方式包括以下步骤：
[0005]步骤1，建立特征序列，包括通过双波长激光雷达数据可以得到355nm及532nm的后向散射信号，利用后向散射信号建立特征序列；
[0006]步骤2，特征序列分类，包括使用层次聚类算法对特征序列进行分类，得到边界层之上的大气分子类和边界层之下的气溶胶类别序列。
[0007]步骤3，查找边界层高度，包括根据类别序列，进行最大梯度检索，从而检索出边界层高度的结果。
[0008]而且，步骤1实现方式为，通过双波长激光雷达两个通道的回波信号计算出355nm及532nm的后向散射信号，将355nm的后向散射信号用A表示，532nm的后向散射信号用B表示，通过A和B信号一一对应构造二维特征序列。
[0009]而且，步骤2实现方式为，将特征序列输入到层次聚类模型中，首先把每个数据点当作一个单独的类，然后计算每个类两两之间的欧式距离，以第i个点和第j个点为例，
[0010][0011]d表示第j个点和第i个点之间的欧式距离，A(j)和A(i)分别表示第j个点和第i个点的355nm后向散射信号大小，B(j)和B(i)分别表示第j个点和第i个点的532nm后向散射信号大小。
[0012]然后，合并数据点中欧式距离最小的数据点为新的一类C。然后计算新类与当前各类之间的距离，使用均链接算法，含义为类C
m
和C
n
的距离等于两个簇类所有样本对的距离平均，数学表达式为，
[0013][0014]其中，D即为类间距离，N
m
和N
n
分别为类C
m
和C
n
中的样本个数，将类间距离最小的两个类再次合并为一个新的类，重复步骤2的计算过程，直到最后类别个数达到输出条件，即实现了对特征序列的分类。同时可以得到具有高度信息的类别序列，将类别序列中属于大气气溶胶的标注为1，属于大气分子的标注为2。
[0015]而且，步骤3实现方式为，对类别序列进行最大梯度查找，即可确定边界层高度所在的点，表达式为，
[0016]BLH＝max(diff(F))(3)
[0017]其中，BLH为反演出的大气边界层高度，F为由步骤2得到的类别序列，diff(F)表示对类别序列进行梯本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于双波长激光雷达的边界层高度反演方法，其特征在于：基于355nm和532nm的后向散射信号及层次聚类算法来获得边界层顶部，实现方式包括以下步骤：步骤1，利用双波长激光雷达数据建立特征序列X，用355nm和532nm后向散射信号构建特征序列X，双波长激光雷达355nm和532nm后向散射信号表示相邻粒子之间的消光能力的差异；步骤2，对特征序列X利用层次聚类算法进行分类，从而得到边界层之上的大气分子类和边界层之下的气溶胶类别序列F；步骤3，查找边界层高度，包括根据类别序列F，进行最大梯度检索，从而检索出边界层高度的结果。2.如权利要求1所述基于双波长激光雷达的边界层高度反演方法，其特征在于：步骤1实现方式为，使用双波长激光雷达两个通道的回波信号计算355nm和532nm后向散射信号，将355nm的后向散射信号用A表示，532nm的后向散射信号用B表示，得到表示相邻粒子间消光能力差异度大小的特征序列X。3.如权利要求1或2所述基于双波长激光雷达的边界层高度反演方法，其特征在于：步骤2实现方式为，利用层次聚类算法对特征序列X进行分类，得到边界层之上的大气分子类和边界层之下的气溶胶类别序列F；将特征序列X输入到层次聚类模型中，首先把每个数据点当作一个单独的类，然后计算每个类两两之间的欧式距离，d表...

【专利技术属性】
技术研发人员：朱忠敏，周向阳，杨威，徐文发，阳威，陈之射，
申请(专利权)人：武昌首义学院，
类型：发明
国别省市：