CN106340208A - 基于4d航迹运行的空中交通管制系统的工作方法 - Google Patents

Abstract

本发明涉及一种基于4D航迹运行的空中交通管制系统的工作方法，所述空中交通管制系统包括数据通信模块、监视数据融合模块、机载终端模块、管制终端模块，其中监视数据融合模块用于实现空管雷达监视数据与自动相关监视数据的融合，为管制终端模块提供实时航迹信息；管制终端模块包括飞行前无冲突4D航迹生成、飞行中短期4D航迹生成、实时飞行冲突监控与告警、飞行冲突解脱4D航迹优化这四个子模块；上述系统的空中交通工作方法，依靠管制终端模块，处理飞行计划数据并利用隐马尔科夫模型生成4D航迹，实现空域交通状况潜在的交通冲突的分析，以及采用模型预测控制理论方法提供最优解脱方案。本发明可有效防止飞行冲突，提高空中交通的安全性。

Description

基于4D航迹运行的空中交通管制系统的工作方法

本申请是申请号为：201510007981.7，发明创造名称为《空中交通管制系统的管制 方法》，申请日为：2015年1月7日的发明专利申请的分案申请。

技术领域

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

背景技术

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

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

基于航迹运行的空中交通运行方式需要在战略层面上对单航空器飞行航迹进行推算和优化，对多航空器构成的交通流实施协同和调整；在预战术层面上通过修正交通流中个别航空器的航迹以解决拥塞问题，并保证该交通流中所有航空器的运行效率；而在战术层面上预测冲突和优化解脱方案，将航空器间隔管理从固定的人工方式转变为考虑航空器性能、管制规则和环境等因素在内的可变的间隔控制方式，因此面向4D航迹的运行对空中交通管制提出了新的要求。

发明内容

本发明要解决的技术问题是在于克服现有技术的不足，提供一种基于4D航迹运行的基于4D航迹运行的空中交通管制系统的工作方法，可有效防止飞行冲突，提高空中交通的安全性。

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

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

飞行前无冲突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；

步骤C、实时飞行冲突监控与告警模块建立从航空器的连续动态到离散冲突逻辑的观测器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能违反空中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供及时的告警信息；

步骤D、飞行冲突解脱4D航迹优化模块在保证系统满足航空器性能和管制规则约束条件下，通过选择不同的解脱目标函数，采用模型预测控制理论方法，计算航空器冲突解脱4D航迹；并通过数据通信模块将航空器冲突解脱4D航迹发送给机载终端模块执行；

步骤E、机载终端模块接收并执行管制终端模块发布的4D航迹数据。

进一步的，步骤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秒。

进一步的，所述步骤D的具体实施过程如下：

步骤D1、对飞行冲突解脱过程建模：将冲突解脱航迹视为连续的三段光滑曲线，给定解脱航迹的起点和终点，依据航迹限制条件，建立包含加速度、爬升或下降率、转弯率的多变量最优冲突解脱模型；

步骤D2、对不同飞行条件下冲突解脱变量约束建模：其中t时刻需实施冲突解脱航空器k的变量约束可描述为：a_k(t)≤a_M、ω_k(t)≤ω_M、γ_k(t)≤γ_M，a_M、ω_M、γ_M分别为最大的加速度、转弯率和爬升或下降率；

步骤D3、设定航空器避撞规划的终止参考点位置P、避撞规划控制时域Θ、轨迹预测时域Υ；

步骤D4、在每一采样时刻，基于航空器当前的运行状态和历史位置观察序列，获取空域风场变量的数值；

步骤D5、设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个航空器赋予不同的权重以及融入实时风场变量滤波数值，得到各个航空器的避撞轨迹和避撞控制策略且各航空器在滚动规划间隔内仅实施其第一个优化控制策略；

步骤D6、在下一采样时刻，重复步骤D4至D5直至各航空器均到达其解脱终点。

更进一步的，步骤D3中：终止参考点位置P即为航空器的下一个航路点，避撞规划控制时域Θ为300秒，轨迹预测时域Υ为300秒；

步骤D4的具体过程如下：

D4.1)设定航空器的停靠位置为轨迹参考坐标原点；

D4.2)在航空器处于直线运行状态和匀速转弯运行状态时，构建空域风场线性滤波模型x(t+△t)＝F(t)x(t)+w(t)和z(t)＝H(t)x(t)+v(t)获取风场变量数值，其中△t表示采样间隔，x(t)表示t时刻的状态向量，z(t)表示t时刻的观测向量，F(t)和H(t)分别表示状态转移矩阵和输出测量矩阵，w(t)和v(t)分别表示系统噪声向量和测量噪声向量；在航空器处于变速转弯运行状态时，构建空域风场非线性滤波模型

x(t+△t)＝Ψ(t,x(t),u(t))+w(t)、z(t)＝Ω(t,x(t))+v(t)和u(t)＝[ω_a(t),γ_a(t)]^T，

其中Ψ(·)和Ω(·)分别表示状态转移矩阵和输出测量矩阵，ω_a(t)和γ_a(t)分别表示转弯率和加速率；

D4.3)根据所构建的滤波模型获取风场变量的数值；

步骤D5的具体过程如下：令

$d_{it}^2=\left|\right|P_i\left(t\right)-P_i^f|{|}_2^2={(x_{it}-x_i^f)}^2+{(y_{it}-y_i^f)}^2,$

其中表示t时刻航空器i当前所在位置P_i(t)和下一航路点P_i ^f间的距离的平方，P_i(t)＝(x_it,y_it)，那么t时刻航空器i的优先级指数可设定为：

$L_{it}=100\frac{d_{it}^{-2}}{\underset{i=1}{\overset{n_t}{\Σ}}{d}_{it}^{-2}},$

其中n_t表示t时刻空域内存在冲突的航空器数目，由优先级指数的含义可知，航空器距离其终点越近，其优先级越高；

设定优化指标

$\begin{array}{c}{J}^{*}(u_1^{\left(t\right)},u_1^{(t+\mathrm{\Δ}t)},...,u_1^{(t+p\mathrm{\Δ}t)},...,u_{n_t}^{\left(t\right)},u_{n_t}^{(t+p\mathrm{\Δ}t)},...,u_{n_t}^{(t+p\mathrm{\Δ}t)})=\underset{s=1}{\overset{p}{\mathrm{\Σ}}}\underset{i=1}{\overset{n_t}{\mathrm{\Σ}}}{L}_{it}\left|\right|P_i(t+s\mathrm{\Δ}t)-P_i^f|{|}_2^2\\ =\underset{s=1}{\overset{p}{\mathrm{\Σ}}}\underset{i=1}{\overset{n_t}{\mathrm{\Σ}}}{\left(P_i\right(t+s\mathrm{\Δ}t)-P_i^f)}^T{Q}_{it}\left(P_i\right(t+s\mathrm{\Δ}t)-P_i^f)\end{array}$

，其中i∈I(t)表示航空器代码且I(t)＝{1,2,…,n_t}，P_i(t+s△t)表示航空器在时刻(t+s△t)的位置向量，P_i ^f表示航空器i的下一航路点，u_i表示待优化的航空器i的最优控制序列，Q_it为正定对角矩阵，其对角元素为航空器i在t时刻的优先级指数L_it，并且

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

进一步的，所述步骤C的具体实施过程如下：

步骤C1、构造基于管制规则的冲突超曲面函数集：建立超曲面函数集用以反映系统的冲突状况，其中，冲突超曲面中与单一航空器相关的连续函数为第I型超曲面，与两架航空器相关的连续函数为第II型超曲面；

步骤C2、建立由航空器连续状态至离散冲突状态的观测器：需要根据管制规范建立观测器，观测系统系统穿越超曲面而产生的冲突事件，以便控制器做出相应的控制决策指令；观测器ξ用于观测系统中航空器位置的连续变化而产生冲突事件，称为第I型观测器，为第II型观测器；

步骤C3、设计从冲突到冲突解脱手段的离散监控器，该离散监控器可描述为函数其中S是观测器观测向量展成的空间，D是所有决策向量d展成的空间；当观测器的离散观测向量表明某一非期望的状态出现时，立刻发出相应的告警。

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

(2)本发明的基于4D航迹运行的空中交通管制系统的工作方法在航空器冲突解脱过程中，融入了高空风场的影响，所采用的滚动解脱轨迹规划方案能够根据高空内风场的变化及时调整解脱轨迹，提高了航空器冲突解脱的鲁棒性。

(3)本发明的基于4D航迹运行的空中交通管制系统的工作方法为航空器配置精确的空管间隔，严格控制航空器通过航路点的时间窗，降低了交通流无序性，提高了空域安全性。

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

(5)本发明的基于4D航迹运行的空中交通管制系统的工作方法不再局限于保持单个航空器之间的安全间隔，而是从宏观上对空域内的交通流实施有效控制，管制工作可以更多的转移到航空器起飞时刻、进场排序、恶劣天气改航等方面。

(6)本发明的基于4D航迹运行的空中交通管制系统的工作方法基于不同性能指标的航空器最优解脱航迹可以显著地提高航空器运行的经济性，以及空域的利用率。

附图说明

图1为本发明的空中交通管制系统的组成示意图；

图2为机载终端模块组成示意图；

图3为数据通信模块组成示意图；

图4为监视数据融合模块组成示意图；

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

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

图7为航空器航迹冲突监控与告警方法流程示意图；

图8为航空器解脱4D航迹优化方法流程示意图。

具体实施方式

(实施例1)

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

1.机载终端模块

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

如图2所示，其具体实施方案如下：

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

2.数据通信模块

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

如图3所示，其具体实施方案如下：

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

3.监视数据融合模块

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

如图4所示，其具体实施方案如下：

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

4.管制终端模块

管制终端模块104包括飞行前无冲突4D航迹生成、飞行中短期4D航迹生成、实时飞行冲突监控与告警、飞行冲突解脱4D航迹优化这四个子模块。

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

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

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

首先，进行航空器状态转移建模。航空器沿航迹飞行的过程表现为在航段之间动态切换过程，根据飞行计划中航空器的飞行高度剖面，建立单个航空器在不同航段转移的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轨迹。

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

首先，对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列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秒。

(3)实时飞行冲突监控与告警

当系统有可能出现违反安全状态集的状态时，通过控制器实施状态监控，对航空器实施有效的管制措施，避免飞行冲突的发生。

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

首先，构造基于管制规则的冲突超曲面函数集。空中交通管制约束的违反都可以视为被控对象(管制空域飞行的多架航空器)构成系统穿越超曲面而产生的事件，建立超曲面函数集用以反映系统的冲突状况。其中，冲突超曲面中与单一航空器相关的连续函数为第I型超曲面，而将与两架航空器相关的连续函数为第II型超曲面。

然后，建立由航空器连续状态至离散冲突状态的观测器。需要根据管制规范建立观测器，观测系统系统穿越超曲面而产生的冲突事件，以便控制器做出相应的控制决策指令。观测器ξ用于观测系统中航空器位置的连续变化而产生冲突事件，称为第I型观测器，为第II型观测器。

最后，设计从冲突到冲突解脱手段的离散监控器。当观测器的离散观测向量表明某一非期望的状态出现时，立刻发出相应的告警。该离散监控器可描述为函数其中S是观测器观测向量展成的空间，D是所有决策向量d展成的空间。

(4)飞行冲突解脱4D航迹优化

在保证使得系统满足控制规范的条件下，通过选择不同的解脱目标函数，采用最优控制理论方法，使得控制器给出的控制输入能达到最优。

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

步骤D1、对飞行冲突解脱过程建模：将冲突解脱航迹视为连续的三段光滑曲线，给定解脱航迹的起点和终点，依据航迹限制条件，建立包含加速度a_i(t)、爬升或下降率γ_i(t)、转弯率ω_i(t)的多变量最优冲突解脱模型。

步骤D2、对不同飞行条件下冲突解脱变量约束建模：其中t时刻需实施冲突解脱航空器k的变量约束可描述为：a_k(t)≤a_M、ω_k(t)≤ω_M、γ_k(t)≤γ_M，a_M、ω_M、γ_M分别为最大的加速度、转弯率和爬升或下降率。

步骤D3、设定航空器避撞规划的终止参考点位置P、避撞规划控制时域Θ、轨迹预测时域Υ。终止参考点位置P即为航空器的下一个航路点，避撞规划控制时域Θ为300秒，轨迹预测时域Υ为300秒。

步骤D4、在每一采样时刻t，基于航空器当前的运行状态和历史位置观察序列，获取空域风场变量的数值，其具体过程如下：

D4.1)设定航空器的停靠位置为轨迹参考坐标原点；

D4.2)在航空器处于直线运行状态和匀速转弯运行状态时，构建空域风场线性滤波模型x(t+△t)＝F(t)x(t)+w(t)和z(t)＝H(t)x(t)+v(t)获取风场变量数值，其中△t表示采样间隔，x(t)表示t时刻的状态向量，z(t)表示t时刻的观测向量，F(t)和H(t)分别表示状态转移矩阵和输出测量矩阵，w(t)和v(t)分别表示系统噪声向量和测量噪声向量；在航空器处于变速转弯运行状态时，构建空域风场非线性滤波模型

x(t+△t)＝Ψ(t,x(t),u(t))+w(t)、z(t)＝Ω(t,x(t))+v(t)和u(t)＝[ω_a(t),γ_a(t)]^T，其中Ψ(·)和Ω(·)分别表示状态转移矩阵和输出测量矩阵，ω_a(t)和γ_a(t)分别表示转弯率和加速率；

D4.3)根据所构建的滤波模型获取风场变量的数值。

步骤D5、设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个航空器赋予不同的权重以及融入实时风场变量滤波数值，得到各个航空器的避撞轨迹和避撞控制策略且各航空器在滚动规划间隔内仅实施其第一个优化控制策略，具体过程如下：令

其中表示t时刻航空器i当前所在位置P_i(t)和下一航路点P_i ^f间的距离的平方，P_i(t)＝(x_it,y_it)，那么t时刻航空器i的优先级指数可设定为：

$L_{it}=100\frac{d_{it}^{-2}}{\underset{i=1}{\overset{n_t}{\mathrm{\Σ}}}{d}_{it}^{-2}},$

其中n_t表示t时刻空域内存在冲突的航空器数目，由优先级指数的含义可知，航空器距离其下一航路点越近，其优先级越高。

设定优化指标

$\begin{array}{c}{J}^{*}(u_1^{\left(t\right)},u_1^{(t+\mathrm{\Δ}t)},...,u_1^{(t+p\mathrm{\Δ}t)},...,u_{n_t}^{\left(t\right)},u_{n_t}^{(t+p\mathrm{\Δ}t)},...,u_{n_t}^{(t+p\mathrm{\Δ}t)})=\underset{s=1}{\overset{p}{\mathrm{\Σ}}}\underset{i=1}{\overset{n_t}{\mathrm{\Σ}}}{L}_{it}\left|\right|P_i(t+s\mathrm{\Δ}t)-P_i^f|{|}_2^2\\ =\underset{s=1}{\overset{p}{\mathrm{\Σ}}}\underset{i=1}{\overset{n_t}{\mathrm{\Σ}}}{\left(P_i\right(t+s\mathrm{\Δ}t)-P_i^f)}^T{Q}_{it}\left(P_i\right(t+s\mathrm{\Δ}t)-P_i^f)\end{array}$

，其中i∈I(t)表示航空器代码且I(t)＝{1,2,...,n_t}，P_i(t+s△t)表示航空器在时刻(t+s△t)的位置向量，P_i ^f表示航空器i的下一航路点，u_i表示待优化的航空器i的最优控制序列，Q_it为正定对角矩阵，其对角元素为航空器i在t时刻的优先级指数L_it，并且

步骤D6、在下一采样时刻，重复步骤D4至D5直至各航空器均到达其解脱终点。

机载终端模块接收并执行管制终端模块发布的4D航迹数据。

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

Claims (2)

1.一种基于4D航迹运行的空中交通管制系统的工作方法，所述空中交通管制系统包括机载终端模块、数据通信模块、监视数据融合模块以及管制终端模块；监视数据融合模块用于实现空管雷达监视数据与自动相关监视数据的融合，为管制终端模块提供实时航迹信息；其特征在于：

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

飞行前无冲突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-yb，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；

步骤C、实时飞行冲突监控与告警模块建立从航空器的连续动态到离散冲突逻辑的观测器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能违反空中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供及时的告警信息；

步骤D、飞行冲突解脱4D航迹优化模块在保证系统满足航空器性能和管制规则约束条件下，通过选择不同的解脱目标函数，采用模型预测控制理论方法，计算航空器冲突解脱4D航迹；并通过数据通信模块将航空器冲突解脱4D航迹发送给机载终端模块执行；

步骤E、机载终端模块接收并执行管制终端模块发布的4D航迹数据；

所述步骤D的具体实施过程如下：

步骤D1、对飞行冲突解脱过程建模：将冲突解脱航迹视为连续的三段光滑曲线，给定解脱航迹的起点和终点，依据航迹限制条件，建立包含加速度、爬升或下降率、转弯率的多变量最优冲突解脱模型；

步骤D2、对不同飞行条件下冲突解脱变量约束建模：其中t时刻需实施冲突解脱航空器k的变量约束可描述为：a_k(t)≤a_M、ω_k(t)≤ω_M、γ_k(t)≤γ_M，a_M、ω_M、γ_M分别为最大的加速度、转弯率和爬升或下降率；

步骤D3、设定航空器避撞规划的终止参考点位置P、避撞规划控制时域Θ、轨迹预测时域Υ；

步骤D4、在每一采样时刻t，基于航空器当前的运行状态和历史位置观察序列，获取空域风场变量的数值；

步骤D5、设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个航空器赋予不同的权重以及融入实时风场变量滤波数值，得到各个航空器的避撞轨迹和避撞控制策略且各航空器在滚动规划间隔内仅实施其第一个优化控制策略；

步骤D6、在下一采样时刻，重复步骤D4至D5直至各航空器均到达其解脱终点。

2.根据权利要求1所述的基于4D航迹运行的空中交通管制系统的工作方法，其特征在于：步骤D3中：终止参考点位置P即为航空器的下一个航路点，避撞规划控制时域Θ为300秒，轨迹预测时域Υ为300秒。