一种飞行器航迹数据驱动的终端区空域时空相关分析方法与流程

本发明属于飞行技术管控领域，具体涉及一种飞行器航迹数据驱动的终端空域运行模式分析方法。

背景技术：

在空中交通管理中，管制人员都是在对所管辖空域运行情况认知理解的基础上制定所有决策的。目前对空域运行状态的理解主要是基于静态的飞行计划，缺乏针对空管实际运行数据的深入分析，把握空域运行的变化规律和时空分布特征，导致在空域运行状态识别、空域流量协同水平判定、飞行调配等方面存在较大差距。

我们提出了一种基于实际航迹数据的终端区空域飞行流模式相关性分析方法，旨在从飞行流相关性的角度揭示空域运行的规律。

技术实现要素：

为了解决上述问题，本发明提供一种飞行器航迹数据驱动的终端空域运行模式分析方法，所述分析方法首先根据部分历史航空器航迹数据进行空间维度聚类，再根据聚类结果对所有航迹数据进行分类，得到终端区内每个机场的全部飞行流模式，最终从飞行流模式相关性的角度分析终端区空域的时空运行模式；

进一步地，所述方法包括以下步骤：

s1：输入d天内终端区内所有航空器的航迹数据，d为大于30的正整数；

s2：根据终端区空间范围剪裁航迹数据；

s3：选取m天的航迹数据进行空间维度聚类得到轨迹模式,m为d/4的整数部分；

s4：将剩余航迹数据按照聚类得到的轨迹模式进行分类，将d天内的所有航迹进行分类得到终端区的飞行流模式；

s5：计算所有飞行流模式彼此之间的时空相关性；

s6：根据时空相关性数值识别关键终端区空域运行模式；

进一步地，所述s2具体如下：根据终端区的水平边界，对原始的航迹数据进行剪裁，保留从终端区进出点至机场的航迹数据集合t，根据航迹数据序列的高度值，区分进场航迹t_in和离场航迹t_out；

进一步地，所述s3具体包括：

s31：分别从t_in和t_out集合中，选择m天航迹作为空间聚类的输入；

s32：基于密度的轨迹空间聚类算法得到终端区内每个机场的轨迹模式集合其中x表示机场集合，y表示某机场的进场轨迹模式集合，z表示某机场的离场轨迹模式集合；

进一步地，所述s2具体为：

s321：面向非均匀采样的航迹密度聚类算法，输入t_in和t_out，最小集合mint，边界距离eps；

s322：计算t_in内所有航迹之间的fréchet距离；

s323：对于尚未处理过的航迹p∈t_in，根据eps计算其fréchet距离意义下的eps临域n_f；

s324：若n_f包含的航迹数量小于mint，则将p标识为噪音；

s3251：否则，创建一个新的轨迹模式c^in_x_y，将p加入c^in_x_y；

s3252：对于每个尚未处理过的航迹p’∈n_f，根据eps计算其fréchet距离意义下的eps临域n’_f；

s3253：如果n’_f包含的航迹数量大于mint，n_f＝n_f+n’_f；

s3254：如果p’还未加入任何轨迹模式，将p’加入c^in_x_y；

进一步地，所述s4具体如下：

s41：将剩余的轨迹，通过聚类得到的轨迹模式，计算每个航迹与轨迹模式之间的eps距离；

s42：将d天的所有航迹均进行分类，分别得到每个机场进港、离港的飞行流模式矩阵wx,z，wx,y，n*p，

其中n表示终端区共有n种轨迹模式，p表示一天分为p个时段，其中wi*j为一个向量，包含j时段第i种轨迹模式的航迹编号；

进一步地，所述s5中计算飞行流模式的时空相关性方法如下：

通过相关函数ccf度量两个空间对象相关性，以ccf作为飞行流时空相关分析的度量指标，通过计算ccf，从每个飞行流模式的角度描述不同时间延迟下与其他飞行流模式之间关于流量特征的相关关系，反映metroplex空域内飞行流时空相关性的局部异质性；

进一步地，所述ccf的计算方法如下：

1)：在计算ccf之前，对于每个飞行流模式，统计每10min间隔内的航班数量作为流量；

2)：构建每个飞行流模式的流量特征时间序列，对于两个流量特征时间序列x和y，在时间延迟s下，x和y的相关函数ccf计算公式为：

式中：ρk为xt和yt+k这两个时间序列时间的相关程度，t为统计时间段；k为时间延迟，xt为x模式在时间段t内的航班数；yt+k为y模式在时间段t+k内的航班数，为x模式在实际用到的t-k个时间段内航班数的平均值；为y模式在实际用到的t-k个时间段内航班数的平均值；

本发明的有益效果如下：

本发明提出了一种基于实际航迹数据的终端区空域飞行流模式相关性分析方法，旨在从飞行流相关性的角度揭示空域运行的规律。拓展了对空域运行规律的认知角度。为空管调配决策提供更全面的信息参考。

附图说明

图1为本发明所述分析方法中整体框架图；

图2为本发明所述分析方法中经过裁剪后的某个终端区的所有航迹图；

图3为本发明所述分析方法中聚类后得到轨迹模式图；

图4为本发明所述分析方法中中所有航迹分类后取中间线结果图；

图5为本发明所述分析方法中飞行流模式经计算后正相关图例；

图6为本发明所述分析方法中飞行流模式经计算后负相关图例；

图7为本发明所述分析方法中飞行流模式经计算后不相关图例。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细描述。应当理解，此处所描述的具体实施例仅仅用于解释本发明，并不用于限定本发明。相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解，在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。

下面结合附图和具体实施例对本发明作进一步说明，但不作为对本发明的限定。下面为本发明的举出最佳实施例：

本发明提供一种飞行器航迹数据驱动的终端区空域时空相关分析方法，所述方法具体步骤如图1所示，包括以下：

s1：输入30天内终端区内所有航空器的航迹数据；

s2：根据终端区空间范围剪裁航迹数据；

s3：选取7天的航迹数据进行空间维度聚类得到轨迹模式；

s4：将剩余航迹数据按照聚类得到的轨迹模式进行分类，将一个月内的所有航迹进行分类得到终端区的飞行流模式，所有航迹分类后取中间线，结果如图4所示；

s5：计算所有飞行流模式彼此之间的时空相关性；

s6：根据时空相关性数值识别关键终端区空域运行模式。

所述s2具体如下：根据终端区的水平边界，对原始的航迹数据进行剪裁，保留从终端区进出点至机场的航迹数据集合t，根据航迹数据序列的高度值，区分进场航迹t_in和离场航迹t_out，经过裁剪后的某个终端区的所有航迹如图2中所示。

所述s3具体包括：

s31：分别从t_in和t_out集合中，选择7天航迹作为空间聚类的输入，聚类后得到轨迹模式结果如图3所示，；

s32：基于密度的轨迹空间聚类算法得到终端区内每个机场的轨迹模式集合其中x表示机场集合，y表示某机场的进场轨迹模式集合，z表示某机场的离场轨迹模式集合；

s321：面向非均匀采样的航迹密度聚类算法，输入t_in(t_out同理)，最小集合mint，边界距离eps；

s322：计算t_in内所有航迹之间的fréchet距离；

s323：对于尚未处理过的航迹p∈t_in，根据eps计算其fréchet距离意义下的eps临域n_f；

s324：若n_f包含的航迹数量小于mint，则将p标识为噪音；

s3251：否则，创建一个新的轨迹模式c^in_x_y，将p加入c^in_x_y；

s3252：对于每个尚未处理过的航迹p’∈n_f，根据eps计算其fréchet距离意义下的eps临域n’_f；

s3253：如果n’_f包含的航迹数量大于mint，则n_f＝n_f+n’_f；

s3254：如果p’还未加入任何轨迹模式，则将p’加入c^in_x_y，c^in_x_y表示的是一种新的轨迹模式，该轨迹模式为进场航迹in的集合，具体包含的数据有机场集合x，机场集合x种所有的进场轨迹模式集合y。

所述s4具体如下：剩余的轨迹，根据聚类得到的轨迹模式，计算每个航迹与轨迹模式之间的eps距离，将一个月内的所有航迹均进行分类；分别得到每个机场进港、离港的飞行流模式矩阵wx,z，wx,y，n*p，n表示终端区共有n种轨迹模式，p表示一天分为p个时段，其中wi*j为一个向量，包含了j时段第i种轨迹模式的航迹编号。

所述s5中计算飞行流模式的时空相关性具体如下：

相关函数(crosscorrelationfunction，ccf)用来度量两个空间对象相关性，在本发明中选择ccf作为飞行流时空相关分析的度量指标。作为一个局部度量指标，本发明中通过计算ccf，从每个飞行流模式的角度描述不同时间延迟下与其他飞行流模式之间关于流量特征的相关关系，从而反映metroplex空域内飞行流时空相关性的局部异质性。在计算ccf之前，对于每个飞行流模式，统计每10min间隔内的航班数量，即流量。构建每个飞行流模式的流量特征时间序列。对于两个流量特征时间序列x和y，在时间延迟s下，x和y的相关函数ccf计算公式为：

式中：ρk为xt和yt+k这两个时间序列时间的相关程度，即从终端区所有的的飞行流中选择的两个飞行流，t为统计时间段；k为时间延迟，xt为x模式在时间段t内的航班数；yt+k为y模式在时间段t+k内的航班数，为x模式在实际用到的t-k个时间段内航班数的平均值；为y模式在实际用到的t-k个时间段内航班数的平均值。各飞行流，图5、图6和图7分别为飞行流模式经计算后正相关图示、飞行流模式经计算后负相关图示和飞行流模式经计算后不相关图示。通过计算ccf，能够量化表示飞行流彼此之间的相关性数值，进而能够发现终端区内飞行流的整体相关关系。

以上所述的实施例，只是本发明较优选的具体实施方式的一种，本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。