基于4D航迹运行的航空器轨迹预测方法与流程

技术领域

本发明涉及一种空中交通管制系统及方法，尤其涉及一种基于4D航迹运行的空中交通管制系统对航空器轨迹进行预测的方法。

背景技术：

随着全球航空运输业快速发展与空域资源有限矛盾的日益突出，在空中交通流密集的复杂空域，仍然采用飞行计划结合间隔调配的空中交通管理方式逐渐显示出其落后性，具体表现在：(1)飞行计划并未为航空器配置精确的空管间隔，容易造成交通流战术管理中的拥挤，降低空域安全性；(2)以飞行计划为中心的空管自动化系统对飞行剖面的推算和航迹预测精度差，造成冲突化解能力差；(3)空中交通管制工作仍然侧重于保持单个航空器之间的安全间隔，很难上升到对交通流进行战略性管理。对于航空器轨迹的预测显得尤为重要。

4D航迹是以空间和时间形式，对某一航空器航迹中的各点空间位置(经度、纬度和高度)和时间的精确描述，基于航迹的运行是指在4D航迹的航路点上使用“控制到达时间”，即控制航空器通过特定航路点的“时间窗”。在高密度空域把基于4D航迹的运行(Trajectory based Operation)作为基本运行机制之一，是未来对大流量、高密度、小间隔条件下空域实施管理的一种有效手段，可以显著地减少航空器航迹的不确定性，提高空域和机场资源的安全性与利用率。

基于航迹运行的空中交通运行方式需要在战略层面上对单航空器飞行航迹进行推算和优化，对多航空器构成的交通流实施协同和调整；在预战术层面上通过修正交通流中个别航空器的航迹以解决拥塞问题，并保证该交通流中所有航空器的运行效率；而在战术层面上预测冲突和优化解脱方案，则非常依赖于能否准确地对航空器的轨迹进行预测，目前均不能准确实时地对航空器的轨迹进行预测，实时性上做的尤为的差。

技术实现要素：

本发明要解决的技术问题是在于克服现有技术的不足，提供一种基于4D航迹运行的空中交通管制系统的基于4D航迹运行的航空器轨迹预测方法，可有效、准确、实时地预测航空器的轨迹。

实现本发明目的的技术方案是提供一种基于4D航迹运行的航空器轨迹预测方法由空中交通管制系统实施，所述空中交通管制系统包括机载终端模块、数据通信模块、监视数据融合模块以及管制终端模块；监视数据融合模块用于实现空管雷达监视数据与自动相关监视数据的融合，为管制终端模块提供实时航迹信息；

所述管制终端模块包括以下子模块：

飞行前无冲突4D航迹生成模块，根据飞行计划和世界区域预报系统的预报数据，建立航空器动力学模型，然后依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成航空器无冲突4D航迹；

飞行中短期4D航迹生成模块，依据监视数据融合模块提供的实时航迹信息，利用隐马尔科夫模型，推测未来一定时间窗内的航空器4D轨迹；

所述空中交通管制系统的基于4D航迹运行的航空器轨迹预测方法包括如下几个步骤：

步骤A、飞行前无冲突4D航迹生成模块根据飞行计划和世界区域预报系统的预报数据，建立航空器动力学模型，并依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成航空器无冲突4D航迹；

步骤B、监视数据融合模块将空管雷达监视数据与自动相关监视数据进行融合，生成航空器实时航迹信息并提供给管制终端模块；管制终端模块中的飞行中短期4D航迹生成模块依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨迹；所述依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨迹的具体实施过程如下：

步骤B6、对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列x＝[x₁,x₂,…,x_n]和y＝[y₁,y₂,…,y_n]，采用一阶差分方法对其进行处理获取新的航空器离散位置序列△x＝[△x₁,△x₂,…,△x_n-1]和△y＝[△y₁,△y₂,…,△y_n-1],其中△x_b＝x_b+1-x_b,△y_b＝y_b+1-y_b(b＝1,2,…,n-1)；

步骤B7、对航空器轨迹数据聚类，对处理后新的航空器离散二维位置序列△x和△y，通过设定聚类个数M'，采用K-means聚类算法分别对其进行聚类；

步骤B8、对聚类后的航空器轨迹数据利用隐马尔科夫模型进行参数训练，通过将处理后的航空器运行轨迹数据△x和△y视为隐马尔科夫过程的显观测值，通过设定隐状态数目N'和参数更新时段ζ'，依据最近的T'个位置观测值并采用B-W算法滚动获取最新隐马尔科夫模型参数λ'；

步骤B9、依据隐马尔科夫模型参数，采用Viterbi算法获取当前时刻观测值所对应的隐状态q；

步骤B10、通过设定预测时域h'，基于航空器当前时刻的隐状态q，获取未来时段航空器的位置预测值O。

进一步的，步骤B中，所述聚类个数M'的值为4，隐状态数目N'的值为3，参数更新时段ζ'为30秒，T'为10，预测时域h'为300秒。

进一步的，步骤B的B8具体是指：由于所获得的航迹序列数据长度是动态变化的，为了实时跟踪航空器航迹的状态变化，有必要在初始航迹隐马尔科夫模型参数λ'＝(π,A,B)的基础上对其重新调整，以便更精确地推测航空器在未来某时刻的位置；每隔时段ζ'，依据最新获得的T'个观测值(o₁,o₂,…,o_T')对航迹隐马尔科夫模型参数λ'＝(π,A,B)进行重新估计；

步骤B的B10具体是指：每隔时段根据最新获得的隐马尔科夫模型参数λ'＝(π,A,B)和最近H个历史观测值(o₁,o₂,…,o_H)，基于航空器当前时刻的隐状态q，通过设定预测时域h'，在时刻t获取航空器在未来时段h'的位置预测值O。

更进一步的，时段为4秒。

进一步的，所述步骤A的航空器无冲突4D航迹按照以下方法生成：

步骤A1、进行航空器状态转移建模，根据飞行计划中航空器的飞行高度剖面，建立单个航空器在不同航段转移的Petri网模型：E＝(g,G,Pre,Post,m)为航空器阶段转移模型，其中g表示飞行航段，G表示垂直剖面中飞行状态参数的转换点，Pre和Post分别表示航段和航路点的前后向连接关系，表示航空器所处的飞行阶段；

步骤A2、建立航空器全飞行剖面混杂系统模型如下，

v_H＝κ(v_CAS,Mach,h_p,t_LOC)，

v_GS＝μ(v_CAS,Mach,h_p,t_LOC,v_WS,α)，

其中v_CAS为校正空速，Mach为马赫数，h_p为气压高度，α为风向预报与航路的夹角，v_WS为风速预报值，t_LOC为温度预报值，v_H为高度变化率，v_GS为地速；

步骤A3、采用混杂系统仿真的方式推测求解航迹：采用将时间细分的方法，利用状态连续变化的特性递推求解任意时刻航空器在某一飞行阶段距参考点的航程和高度其中J₀为初始时刻航空器距参考点的航程，△τ为时间窗的数值，J(τ)为τ时刻航空器距参考点的航程，h₀为初始时刻航空器距参考点的高度，h(τ)为τ时刻航空器距参考点的高度，由此可以推测得到单航空器的4D航迹；

步骤A4、对多航空器耦合模型实施无冲突调配：根据两航空器预达交叉点的时间，按照空中交通管制原则，对交叉点附近不满足间隔要求的航空器4D航迹进行二次规划，得到无冲突4D航迹。

进一步的，所述步骤B中监视数据融合模块将空管雷达监视数据与自动相关监视数据进行融合，生成航空器实时航迹信息，具体按照以下方法：

步骤B1、将坐标单位和时间统一；

步骤B2、采用最邻近数据关联算法将属于同一个目标的点相关联，提取目标航迹；步骤B3、将分别从自动相关监视系统和空管雷达提取的航迹数据从不同的时空参

考坐标系统变换、对准到管制终端统一的时空参考坐标系统；

步骤B4、计算两条航迹的相关系数，若相关系数小于某一预设阈值，则认为两条航迹不相关；否则该两条航迹相关，可以进行融合；

步骤B5、对相关的航迹进行融合。

更进一步的，所述步骤B5中对相关的航迹进行融合，采用基于采样周期的加权平均算法，其加权系数根据采样周期和信息精度确定，再利用加权平均算法将与之相关的自动相关监视航迹和空管雷达航迹融合为系统航迹。

本发明具有积极的效果：(1)本发明的空中交通管制系统的基于4D航迹运行的航空器轨迹预测方法在航空器实时轨迹推测过程中，融入了随机因素的影响，所采用的滚动轨迹推测方案能够及时提取外界随机因素的变化状况，提高了航空器轨迹推测的准确性。

(2)本发明的空中交通管制系统的基于4D航迹运行的航空器轨迹预测方法对飞行剖面的推算和航迹预测精度高，进而使得冲突化解能力和自动化水平提高，降低了管制员的工作负荷。

附图说明

图1为飞行前无冲突4D航迹生成方法流程示意图；

图2为飞行中短期4D航迹推测方法流程示意图。

具体实施方式

(实施例1)

本实施例的基于4D航迹运行的空中交通管制系统，包括机载终端模块101、数据通信模块102、监视数据融合模块103以及管制终端模块104。以下对各部分的具体实施方式分别进行详细描述。

1.机载终端模块

机载终端模块101是飞行员获取地面管制指令、参考4D航迹，以及输入飞行意图的界面，同时还是采集当前航空器位置数据的接口。

其具体实施方案如下：

机载终端模块101接收如下的信息输入：(1)ADS-B信息采集单元201通过机载GPS采集的航空器位置向量、速度向量，以及本航空器的呼号，编码后通过信息及数据传递给机载数据通信模块102；(2)航空器驾驶员需要将与地面管制指令不一致的飞行意图，通过人机输入界面，以及约定的地面管制员可以识别的形式通过信息及数据传递给机载数据通信模块102。另外机载终端模块101实现如下的信息输出：(1)通过终端显示屏幕，接收和显示飞行员可以识别的飞行管制指令；(2)接收和显示地面管制终端飞行前生成的无冲突4D航迹，以及当地面管制终端探测到冲突后计算的最优解脱4D航迹。

2.数据通信模块

数据通信模块102可实现空地双向数据通信，实现机载实时位置数据和飞行意图数据单元202的下行传输和地面管制指令单元203，以及参考4D航迹单元204的上行传输。

其具体实施方案如下：

下行数据通信：机载终端101通过机载二次雷达应答机将航空器识别标志和4D位置信息，以及其他附加数据，如飞行意图、飞行速度、气象等信息传输给地面二次雷达(SSR)，二次雷达接收后对数据报文进行解析，并传输给中央数据处理组件301解码，通过指令航迹数据接口传输到管制终端104；上行数据通信：地面管制终端104通过指令航迹数据接口，经中央数据处理组件301编码后，地面二次雷达的询问机将将地面管制指令或参考4D航迹信息传递并显示在机载终端101。

3.监视数据融合模块

监视数据融合模块103实现空管雷达监视与自动相关监视ADS-B数据的融合，为管制终端模块104中的飞行中短期4D航迹生成子模块和实时飞行冲突监控与告警子模块提供实时航迹信息。

其具体实施方案如下：

(1)在预处理阶段将坐标单位和时间统一，假设分别从ADS-B和空管雷达中提取的数据是一系列离散点的坐标(如经度、纬度、海拔高度)、各点对应采集时间；(2)采用最邻近数据关联算法将属于同一个目标的点相关联，提取目标航迹；(3)将分别从ADS-B和空管雷达中提取的航迹数据从不同的时空参考坐标系统变换、对准到管制终端统一的时空参考坐标系统；(4)计算两条航迹的相关系数，若相关系数小于某一预设阈值，则认为两条航迹不相关，否则该两条航迹相关，可以进行融合；(5)对相关的航迹进行融合。由于ADS-B和空管雷达的精度和采样周期不同，本系统采用基于采样周期的加权平均算法，其加权系数根据采样周期和信息精度确定，再利用加权平均算法将与之相关的ADS-B航迹和空管雷达航迹融合为系统航迹。

4.管制终端模块

管制终端模块104包括飞行前无冲突4D航迹生成、飞行中短期4D航迹生成这2个子模块。

(1)飞行前无冲突4D航迹生成

根据飞行数据处理系统(FDP)得到的飞行计划和世界区域预报系统(WAFS)发布的风、温度的GRIB格点预报数据，对空中交通系统建立层次化的混杂系统模型，通过系统在安全状态的演化，描述状态演化的时间轨迹，生成航空器航迹。

如图1所示，其具体实施过程如下：

首先，进行航空器状态转移建模。航空器沿航迹飞行的过程表现为在航段之间动态切换过程，根据飞行计划中航空器的飞行高度剖面，建立单个航空器在不同航段转移的Petri网模型：E＝(g,G,Pre,Post,m)为航空器阶段转移模型，其中g表示飞行航段，G表示垂直剖面中飞行状态参数(包括空速、高度、构型)的转换点，Pre和Post分别表示航段和航路点的前后向连接关系，表示航空器所处的飞行阶段。

其次，建立航空器全飞行剖面混杂系统模型。航空器在单个航段内的飞行视为连续过程，依据质点能量模型，推导航空器在不同的运行阶段同气象条件下的航空器动力学方程，v_H＝κ(v_CAS,Mach,h_p,t_LOC)，v_GS＝μ(v_CAS,Mach,h_p,t_LOC,v_WS,α)，其中v_CAS为校正空速，Mach为马赫数，h_p为气压高度，α为风向预报与航路的夹角，v_WS为风速预报值，t_LOC为温度预报值，v_H为高度变化率，v_GS为地速。

然后，采用混杂系统仿真的方式推测求解航迹。采用将时间细分的方法，利用状态连续变化的特性递推求解任意时刻航空器在某一飞行阶段距参考点的航程和高度其中J₀为初始时刻航空器距参考点的航程，△τ为时间窗的数值，J(τ)为τ时刻航空器距参考点的航程，h₀为初始时刻航空器距参考点的高度，h(τ)为τ时刻航空器距参考点的高度，由此可以推测得到单航空器的4D航迹。

最后，对多航空器耦合模型实施无冲突调配。根据两航空器预达交叉点的时间，按照空中交通管制原则，对交叉点附近不满足间隔要求的航空器4D航迹进行二次规划，得到无冲突4D航迹。

(2)飞行中短期4D航迹生成

依据管制雷达和自动相关监视系统ADS-B实施融合后获得航空器实时航迹数据，利用隐马尔科夫模型，推测未来5分钟时间窗内的航空器4D轨迹。

如图2所示，其具体实施过程如下：

首先，对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列x＝[x₁,x₂,…,x_n]和y＝[y₁,y₂,…,y_n]，采用一阶差分方法对其进行处理获取新的航空器离散位置序列△x＝[△x₁,△x₂,…,△x_n-1]和△y＝[△y₁,△y₂,…,△y_n-1],其中△x_b＝x_b+1-x_b,△y_b＝y_b+1-y_b(b＝1,2,…,n-1)。

其次，对航空器轨迹数据聚类。对处理后新的航空器离散二维位置序列△x和△y，通过设定聚类个数M'，采用K-means聚类算法分别对其进行聚类。

然后，对聚类后的航空器轨迹数据利用隐马尔科夫模型进行参数训练。通过将处理后的航空器运行轨迹数据△x和△y视为隐马尔科夫过程的显观测值，通过设定隐状态数目N'和参数更新时段ζ'，依据最近的T'个位置观测值并采用B-W算法滚动获取最新隐马尔科夫模型参数λ'：由于所获得的航迹序列数据长度是动态变化的，为了实时跟踪航空器航迹的状态变化，有必要在初始航迹隐马尔科夫模型参数λ'＝(π,A,B)的基础上对其重新调整，以便更精确地推测航空器在未来某时刻的位置。每隔时段ζ'，依据最新获得的T'个观测值(o₁,o₂,…,o_T')对航迹隐马尔科夫模型参数λ'＝(π,A,B)进行重新估计。

再而，依据隐马尔科夫模型参数，采用Viterbi算法获取当前时刻观测值所对应的隐状态q。

最后，每隔时段根据最新获得的隐马尔科夫模型参数λ'＝(π,A,B)和最近H个历史观测值(o₁,o₂,…,o_H)，基于航空器当前时刻的隐状态q，通过设定预测时域h'，在时刻t获取航空器在未来时段h'的位置预测值O。

所述聚类个数M'的值为4，隐状态数目N'的值为3，参数更新时段ζ'为30秒，T'为10，预测时域h'为300秒，时段为4秒。

显然，上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。