-
飞机上下翼面存在压力差,其经过的区域的空气受到扰动,最终呈现为一对涡旋的形式,其强度用涡环量来衡量。涡环量与飞机的重量成正比和飞机的速度成反比,初始涡环量Γ0用下式表示[1]:
$ {\Gamma _0} = \frac{{Mg}}{{\rho SBv}} = \frac{{Mg}}{{\rho {b_0}v}} $
(1) 式中,M, B分别表示飞机的质量和翼展,g为重力加速度,ρ为周围空气的密度,S为载荷系数,v为航空器速度,b0为初始涡核间距。
在模拟飞机尾涡剖面径向速度的模型中,Hallock-Burnham(HB)模型因使用实际数据进行修正且表述较为简单进而获得了广泛的应用,HB模型表示见下[1]:
$ {v_{\rm t}}(r) = \frac{{{\Gamma _0}}}{{2{\rm \pi} r}}\frac{{{r^2}}}{{{r^2} + {r_{\rm c}}^2}} = \frac{{{\Gamma _0}r}}{{2{\rm \pi} ({r^2} + {r_{\rm c}}^2)}} $
(2) 式中, r为距涡核的距离,rc为涡核半径(通常设为0.052b0),vt(r)表示切向速度。
激光雷达通过对大气气溶胶粒子的探测,在晴空条件下能对航空器尾涡进行较为精确的刻画。由于飞机尾涡的影响面是对飞机飞过垂直平面的空气扰动,因此, 常规的探测方式是对飞机的飞行轨迹的垂直剖面进行距离高度指示器(range height indicator, RHI)模式的探测,具体如图 1所示。
图 1中,O1,O2,O分别为左、右涡心和激光雷达探测点; α1,α2,θ分别表示雷达RHI扫描时左、右涡和探测空间点俯仰角; d1,d2,d分别为激光雷达点距离左、右涡心和探测空间点的距离; r1,r2分别表示探测空间点距离左、右涡心的距离。
激光雷达处于RHI模式探测时,方位角固定,改变俯仰角对激光雷达接收回波做快速傅里叶变换,得到固定距离门的径向风速速度,负值为接近激光雷达方向的速度,正值为远离激光雷达方向的速度,数值的大小即为速度的值。当雷达从设定的扫描俯仰角扫描一周时,即得到一个切面的径向速度风场:
$ {V_{\rm r}} = \left\{ {{v_{\rm r}}({\rho _i},{\theta _j})} \right\},(i \in [1,n],j \in [1,m]) $
(3) 式中,ρi表示距离激光雷达的径向距离, θj表示激光雷达的俯仰角角度, vr表示该极坐标对应的径向速度, n和m分别是激光雷达探测扫描所得到的俯仰角离散值的个数和径向距离门的个数,实地探测采用的激光雷达,距离门个数n=56。
以全球飞机类型数量最多的空客A320为例,其机型参量如表 1所示。
Table 1. Parameters of Airbus A320 and environment
parameter value M 78000kg B 37.57m v 69.96m/s g 9.81m/s2 ρ 1.16kg/m3 S π/4 将表 1中的机型参量和环境参量代入(1)式和(2)式,设定飞机左右涡的涡心距离地面的高度为50m,涡点连线的中点距离激光雷达的水平距离为200m,背景风场为静风。图 2a为依据这些参量所计算出的切向速率图。静风条件下HB模型得出的尾涡为左右对称的涡旋结构,最大切向速度值为15.30m/s。表 2中列出了实地探测中使用激光雷达的特征参量。
Table 2. Main parameters of the LiDAR
parameters value detection range 45m~6000m range resolution 15m angle resolution 0.1° laser source 1.55μm detection accuracy 0.1m/s 受限于激光雷达探测原理,脉冲激光雷达扫描得到的速度为径向速度,相对切向速度,径向速度损失了垂直于雷达扫描径向上的速度分量, 且有限的径向分辨率和角度分辨率,所得到的径向速度为(3)式所示的离散数据。根据现场探测试验的实际参量,假设脉冲多普勒激光雷达探测精度良好,径向分辨率为15m,角度分辨率为0.2°,根据图 1中描述的位置关系和设置的参量,计算出图 2a的切向速度场激光雷达探测到的径向速度场, 如图 2b所示。此时最大径向速度值为7.87m/s,相对15.30m/s的最大切向速度值损失了接近一半,因此,使用传统尾涡模型在匹配激光雷达所探测的径向速度场时存在较大的误差。此外,根据HB模型计算结果可得,空客A320的尾流在静风条件下的径向风速场影响区域为长100m、高70m的矩形区域。
-
多普勒测风激光雷达每次扫描获取的为离散数据,即固定方向固定距离的径向速度,如(3)式所示。如将径向速度场直接当作特征用于分类,同样的一组径向风速,不同的排列(角度、距离)可能会有截然不同的分类结果,因此在对激光雷达径向风速数据识别分类时,需先进行特征提取,提取中其中能表征飞机尾涡的特征。
本团队于2018年8月至10月在国内某机场使用激光雷达进行飞机尾流的探测实验,图 3展示为现场探测激光雷达和飞机的位置。作者从探测数据中选取了包含了具有一定背景风速较为典型空客A320的尾涡数据和非尾涡数据合计1273组,其中643组数据存在尾涡,另外630组数据不存在尾涡。图 4为实地连续探测的一组空客A320的径向速度图。尾流影响范围为长100m、高70m的矩形区域,这与HB模型计算的结果基本一致。整个背景风场的风向大体趋势是远离雷达的方向,其中红色和绿色区域相接的地方是探测出的尾涡区域,可以观察到尾涡受背景风场等因素的影响,朝右下移动。对比理论模型计算出的模拟探测图(如图 2b所示),在实际探测的数据中背景风场具有一定的风速,风场存在非均匀现象,此时使用波形匹配的方法计算出的波形和标准涡相差较大,对于图 4b,使用图像识别的方法特征已经不是很明显,变化色标的范围,效果会更差。
考虑到背景风场为影响飞机尾涡的重要因素之一,在非均匀背景风场下,初始涡和标准涡有较大差异,因此本文中选择平均背景风速作为重要的参考特征,用于表征背景风场对于飞机尾涡的影响,如下式所示:
$ \overline {v_{\rm r}} = \frac{1}{{nm}}\sum\limits_{i = 1}^n {\sum\limits_{j = 1}^m {{v_{\rm r}}({\rho _i},{\theta _j})} } $
(4) 图 5分别为多普勒激光雷达RHI扫描模式下,同一径向风速场的径向风速图、距离速度极差图和仰角-速度极差图。如图所示,真实的背景风场存在非均匀性,传统的方法对径向风速场按照俯仰角或距离提取多普勒速度极差[12],从而判断极差最大的位置为涡核中心,存在误判的可能。以距离为基准的最大正负速度包络曲线处代表的区域并不是尾涡区域,以扫描俯仰角为基准的最大正负速度包络曲线也不是尾涡区域。
鉴于此,本文中将完整径向速度场区域内的径向速度极差作为特征参量,即:
$ \begin{array}{c} \Delta v({\rho _1},{\rho _2},{\theta _1},{\theta _2}) = {\rm \max} \left( {\left\{ {{v_r}({\rho _i},{\theta _j})} \right\}} \right) - \\ {\rm \min} \left\{ {{v_r}({\rho _i},{\theta _j})} \right\} \end{array} $
(5) 式中,∣ρi cos θj-x0∣≤w/2, ∣ρi sin θj-y0∣≤h/2, w为设定矩形区域的宽度,h表示滑动窗口的长度,根据第1节中的计算和第2节中的实地探测的结果,设置w=100m,h=70m,x0为滑动窗口的中心横坐标,y0表示滑动窗口的中心纵坐标。
-
KNN分类算法是经典的统计学分类方法之一[13],其原理是通过从训练数据集中选取最接近待分类数据的k个实例,根据其参量和类别预测样本的分类,如图 6所示。图 6a中正方形和正六边形分别代表已分类的数据集,五角星代表待分类的数据[14-15]。
训练数据集T为KNN输入:
$ \mathit{\boldsymbol{T}} = [({\mathit{\boldsymbol{x}}_1},{\mathit{\boldsymbol{y}}_1}),({\mathit{\boldsymbol{x}}_2},{\mathit{\boldsymbol{y}}_2}), \cdots ,({\mathit{\boldsymbol{x}}_N},{\mathit{\boldsymbol{y}}_N})] $
(6) 式中,xi∈X⊆Rn为实例的特征向量,yi∈Y⊆{c1, c2, …, cK}表示实例的类别,对于径向风速场的类别,分为两类:存在尾涡和不存在尾涡,N为数据集的大小。
对于单组数据的特征xi包含了的特征表示如下:
$ {\mathit{\boldsymbol{x}}_i} = ({\mathit{\boldsymbol{x}}_{i,(1)}},{\mathit{\boldsymbol{x}}_{i,(2)}}, \cdots {\mathit{\boldsymbol{x}}_{i,(l)}}) $
(7) 式中, 下标i表示第i组数据, 下标(l)表示数据的第l维特征。此处l=2,对应第2节中所提出的平均背景风速和矩形区域速度极差,作为KNN分类的特征参量。影响KNN分类结果的因素主要包括已分类数据集、k值的选择、距离度量方式和分类决策规则。
如图 6b所示, k值分别为1, 3和5时,对于五角星的待分类点的预测结果分别为正六边形、正六边形和正方形,k值的选择对于预测分类的结果影响极大。选择较大的k值会减少特征较为相似的数据对预测分类结果的影响,较小的k值则会增加异常数据对于预测分类结果的影响,因此k值应依据数据的特点进行选取。在具体应用中,k值一般取一个比较小的数值,采用交叉验证法来选取最优k值。使用交叉验证的方法来确定k值的大小,通过将数据70%作为KNN的训练的训练数据集,故训练集的大小N=891,剩下部分作为测试集。
KNN采用闵可夫斯基(Minkowski)距离Lp来度量数据间的相似度:
$ d({\mathit{\boldsymbol{x}}_i},{\mathit{\boldsymbol{x}}_j}) = {L_p} = p\sqrt {\sum\limits_{l = 1}^n {{{\left| {{\mathit{\boldsymbol{x}}_{i,(l)}} - {\mathit{\boldsymbol{x}}_{j,(l)}}} \right|}^p}} } $
(8) 式中,xi表示测试点第i维特征值,xi, (l)表示样本点第l维特征值,p为常数值。
在实际应用过程中,KNN一般选用p=2时的闵可夫斯基距离,即欧氏距离,也是本文中所采用的距离度量方式,即:
$ d({\mathit{\boldsymbol{x}}_i},{\mathit{\boldsymbol{x}}_j}) = {L_2} = \sqrt {\sum\limits_{l = 1}^n {{{\left( {{\mathit{\boldsymbol{x}}_{i,(l)}} - {\mathit{\boldsymbol{x}}_{j,(l)}}} \right)}^2}} } $
(9) 分类决策采用多数表决法,即由输入实例的k个邻近的训练实例中多数类决定输入实例的类。
-
为评估不同k值的kNN分类方法的性能,采用准确率(accuracy, ACC)A、正预测值(positive predicted value,PPV)VPPV和真正率(true positive rate, TPR)RTPR作为评估指标,根据径向风场是否存在尾涡,已经KNN分类预测是否存在尾涡可得出该分类的混淆矩阵,如表 3所示[16]。
Table 3. Confusion matrix
condition positive condition negative predicted condition positive true positive(TP) false positive(FP) predicted condition negative false negative(FN) true negative(TN) 准确率A表示正确识别存在尾涡和不存在尾涡的径向风速场数据占所有数据的比率:
$ A = \frac{{\rm TP + TN}}{{\rm TP + TN + FP + FN}} $
(10) 正预测值VPPV实际存在的尾涡的径向风速场数据占识别存在的尾涡的径向风速场数据的比率:
$ {V_{\rm PPV}} = \frac{{\rm TP}}{{\rm TP + FP}} $
(11) 真正率(true positive rate, TPR)RTPR为所有存在尾涡的径向风场数据中,成功识别的比率:
$ {R_{\rm TPR}} = \frac{{\rm TP}}{{\rm TP + FN}} $
(12) 真负率(true negative rate, TNR)RTNR为所有不存在尾涡的径向风场数据中,成功识别的比率:
$ {R_{\rm TNR}} = \frac{{\rm TN}}{{\rm TN + FP}} $
(13) 考虑实际应用中VPPV和RTPR有时差距过大,这里再额外纳入数值F1来综合评估KNN模型的分类性能:
$ {F_1} = \frac{2}{{1/{V_{\rm PPV}} + 1/{R_{\rm TPR}}}} $
(14) 根据第3节中提到的交叉验证法,在实地探测的1273组径向速度场数据中,随机选取30%作为测试集数据,其余70%作为训练集数据, 测试得出KNN分类方法不同k值下的性能如图 7所示。
训练集性能指标随着k值的增大,总体上呈现的为下降的趋势,而测试集性能指标随着k值的增大,性能指标总体均为先增大后减小,不同k值会显著影响该方法的性能。当k=11时,分类的结果如表 3所示。依据分类结果计算分类器的性能,如表 4所示。训练集的平均性能指标达到了0.816,其中A=0.796,测试集的平均性能指标达到了0.729,其中A=0.717,按照样本数量对训练集和测试集的加权A达到0.772,效果较好,因此本文中选择的是11NN分类方法。
Table 4. Results of 11NN classifier
data set TP FN FN TN train set 340 48 134 369 test set 84 23 85 190 Table 5. Performance of 11NN classifier
data set ACC PPV TPR TNR F1 average train set 0.796 0.876 0.717 0.855 0.834 0.816 test set 0.717 0.785 0.497 0.892 0.750 0.729 作为评估预测与分类性能的指标,受试工作者曲线(receiver operating characteristic, ROC)已被广泛地用于医学以及机器学习等领域。当正负数据发生剧烈变化时候,其它指标也随时会发生较大的变化,而ROC曲线却能接近不变, 因此ROC曲线可全面客观地评价模型的性能, AUC(area under curve)是ROC曲线下面积,用于定量的评价分类型性能,能够在一定程度上衡量模型的鲁棒性。本文中提出的KNN的方法很大程度上受限于测试的数据,同样的模型并不一定具有通用性,因此使用ROC和AUC来评价11NN分类方法的通用性能。
图 8a是11NN分类器中不同阈值下的实际存在尾涡和实际不存在尾涡的数据概率分布。据此能得到对应的11NN分类器所得到的ROC曲线,如图 8b所示,其远离TPR=FPR的虚线,相应的AUC的值为0.855,说明11NN分类器在正负数据发生变化的同时,该分类器的性能不会发生较大的变化,说明11NN分类器具备了较好的鲁棒性,根据上面得出11NN的性能参量,说明本文中提出的方法能够很好地实现对尾涡的预测与分类。
基于k最近邻的激光雷达飞机尾涡识别
Identification of aircraft wake vortex based on k-nearest neighbor
-
摘要: 为了对脉冲多普勒激光雷达探测径向速度场进行识别,建立了基于k最近邻(KNN)方法的分类模型。结合飞机尾涡Hallock-Burnham模型和脉冲多普勒激光雷达特性,提取脉冲多普勒激光雷达探测径向速度场的特征参量。基于现有测试数据,在非均匀背景风场下利用KNN方法对飞机尾涡进行识别。结果表明,以准确率(ACC)和ROC曲线下的面积(AUC)作为评估标准,所提出的方法对尾流识别所获得的ACC和AUC分别为0.772和0.855。该方法可有效地识别飞机尾涡并具备较好的鲁棒性。Abstract: To identify aircraft wake vortex by pulsed doppler LiDAR's characteristics, a classification model based on k-nearest neighbor (KNN) was established in this paper. This approach by combining Hallock-Burnham model with pulsed doppler LiDAR's characteristics to extract the feature parameters of radial velocity of wind field was pursued. Based on the test dataset, the KNN was employed to identify aircraft wake vortex in the context of nonuniform wind field. The performance of the proposed method was evaluated in terms of the accuracy (ACC) and the area under ROC curve (AUC). The ACC and the AUC of our technique on test dataset are 0.772 and 0.855, respectively. Experimental results are presented to illustrate the validity and robustness of the proposed approach to aircraft wake vortex.
-
Key words:
- laser technique /
- Doppler LiDAR /
- k-nearest neighbor /
- feature extraction /
- wake vortex identification
-
Table 1. Parameters of Airbus A320 and environment
parameter value M 78000kg B 37.57m v 69.96m/s g 9.81m/s2 ρ 1.16kg/m3 S π/4 Table 2. Main parameters of the LiDAR
parameters value detection range 45m~6000m range resolution 15m angle resolution 0.1° laser source 1.55μm detection accuracy 0.1m/s Table 3. Confusion matrix
condition positive condition negative predicted condition positive true positive(TP) false positive(FP) predicted condition negative false negative(FN) true negative(TN) Table 4. Results of 11NN classifier
data set TP FN FN TN train set 340 48 134 369 test set 84 23 85 190 Table 5. Performance of 11NN classifier
data set ACC PPV TPR TNR F1 average train set 0.796 0.876 0.717 0.855 0.834 0.816 test set 0.717 0.785 0.497 0.892 0.750 0.729 -
[1] MURPHY B, O'CALLAGHAN J, FOX M, et al. Overview of the structures investigation for the american airlines flight 587 investigation[C]//46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. New York, USA: IEEE, 2013: 3369-3373. [2] GERZ T, HOLZÄPFEL F, DARRACQ D. Commercial aircraft wake vortices[J]. Progress in Aerospace Sciences, 2002, 38(3): 181-208. doi: 10.1016/S0376-0421(02)00004-0 [3] CORNELIUS J O, RUTISHAUSER D K. Enhanced airport capacity through safe, dynamics reductions in aircraft separation//Aircraft VOrtex Spacing System (AVOSS).New York, USA: IEEE, 2001: 77-103. [4] HOLZÄPFEL F, GERZ T, FRECH M, et al. The wake vortex prediction and monitoring system WSVBS Part Ⅰ: Design[J]. Air Traffic Control Quarterly, 2008, 17(9):301-322. [5] GERZ T, HOLZÄPFEL F, GERLING W, et al. The wake vortex prediction and monitoring system WSVBS(Part Ⅱ: Performance and ATC integration at Frankfurt airport)[J]. Air Traffic Control Quarterly, 2009, 17(4): 323-346. doi: 10.2514/atcq.17.4.323 [6] LIANG H J, DENG W X, LIANG Y A, et al. Review of aerocraft wake vortex separation dynamic reduction technology[J]. Journal of Ordnance Equipment Engineering, 2018, 39(12):15-19(in Chinese). [7] HANNON S M, THOMSON J A. Aircraft wake vortex detection and measurement with pulsed solid-state coherent laser radar[J]. Journal of Modern Optics, 1994, 41(11):2175-2196. doi: 10.1080/09500349414552031 [8] HARRIS M, YOUNG R I, KÖPP F, et al. Wake vortex detection and monitoring[J]. Aerospace Science and Technology, 2002, 6(5):325-331. doi: 10.1016/S1270-9638(02)01171-9 [9] XU Sh L, HU Y H, WU Y H. Identification of aircraft wake vortex based on Doppler spectrum features[J]. Journal of Optoelectronics·Laser, 2011, 22(12): 1826-1830(in Chinese). [10] PAN W J, ZHANG Q Y, ZHANG Q, et al. Identification method of aircraft wake vortex based on Doppler lidar[J]. Laser Technology, 2019, 43(2):233-237(in Chinese). [11] PAN W J, DUAN Y J, ZHANG Q, et al. Research on aircraft wake vortex recognition using AlexNet[J]. Opto-Electronic Engineering, 2019, 46(7): 190082(in Chinese). [12] ZHAO L Y, GU R P, WEI Z Q. Calculation of characteristics parameters of dynarnic wake vortex based on lidar echo[J]. Journal of Wuhan University of Science and Technology, 2018, 41(5):388-394(in Chinese). [13] LI H. Statistical learning method [M]. 2nd ed. Beijng: Tsinghua University Press, 2012:49-58(in Chinese). [14] COVER T, HART P. Nearest neighbor pattern classification[J]. IEEE Transactions on Information Theory, 2003, 13(1):21-27. [15] FAWCETT T. An introduction to ROC analysis[J]. Pattern Recognition Letters, 2006, 27(8):861-874. doi: 10.1016/j.patrec.2005.10.010 [16] RUPPERT D. The elements of statistical learning: Data mining, inference, and prediction[J]. Journal of the American Statistical Association, 2004, 99(466):567-567. doi: 10.1198/jasa.2004.s339