一种空中交通管制系统的飞行冲突预警方法
【专利摘要】本发明涉及一种空中交通管制系统的飞行冲突预警方法，所述空中交通管制系统包括数据通信模块、监视数据融合模块、机载终端模块、管制终端模块，其中监视数据融合模块用于实现空管雷达监视数据与自动相关监视数据的融合，为管制终端模块提供实时航迹信息；管制终端模块包括飞行前无冲突4D航迹生成、飞行中短期4D航迹生成、实时飞行冲突监控与告警、这3个子模块；上述系统的飞行冲突预警方法，依靠管制终端模块，处理飞行计划数据并利用隐马尔科夫模型生成4D航迹，实现空域交通状况潜在的交通冲突的分析。本发明可有效预警飞行冲突，提高空中交通的安全性。
【专利说明】一种空中交通管制系统的飞行冲突预警方法
【技术领域】
[0001] 本发明涉及一种空中交通管制系统及方法，尤其涉及一种基于4D航迹运行的空 中交通管制系统对航空器轨迹进行预测并对飞行冲突预警的方法。
【背景技术】
[0002] 随着全球航空运输业快速发展与空域资源有限矛盾的日益突出，在空中交通流密 集的复杂空域，仍然采用飞行计划结合间隔调配的空中交通管理方式逐渐显示出其落后 性，具体表现在：（1)飞行计划并未为航空器配置精确的空管间隔，容易造成交通流战术管 理中的拥挤，降低空域安全性；(2)以飞行计划为中心的空管自动化系统对飞行剖面的推 算和航迹预测精度差，造成冲突化解能力差；(3)空中交通管制工作仍然侧重于保持单个 航空器之间的安全间隔，很难上升到对交通流进行战略性管理。对于航空器轨迹的预测并 从而并对飞行冲突预警显得尤为重要。
[0003] 4D航迹是以空间和时间形式，对某一航空器航迹中的各点空间位置（经度、炜度 和高度）和时间的精确描述，基于航迹的运行是指在4D航迹的航路点上使用"控制到达 时间"，即控制航空器通过特定航路点的"时间窗"。在高密度空域把基于4D航迹的运行 (TrajectorybasedOperation)作为基本运行机制之一，是未来对大流量、高密度、小间隔 条件下空域实施管理的一种有效手段，可以显著地减少航空器航迹的不确定性，提高空域 和机场资源的安全性与利用率。
[0004] 基于航迹运行的空中交通运行方式需要在战略层面上对单航空器飞行航迹进行 推算和优化，对多航空器构成的交通流实施协同和调整；在预战术层面上通过修正交通流 中个别航空器的航迹以解决拥塞问题，并保证该交通流中所有航空器的运行效率；而在战 术层面上预测冲突和优化解脱方案，则非常依赖于能否准确地对航空器的轨迹进行预测并 对飞行冲突进行预警，目前均不能准确实时地对航空器的轨迹进行预测，实时性上做的尤 为的差。
【发明内容】
[0005] 本发明要解决的技术问题是在于克服现有技术的不足，提供一种基于4D航迹运 行的空中交通管制系统的飞行冲突预警方法，可有效、准确、实时地预测航空器的轨迹并预 测飞行冲突。
[0006] 实现本发明目的的技术方案是提供一种空中交通管制系统的飞行冲突预警方法， 所述空中交通管制系统包括机载终端模块、数据通信模块、监视数据融合模块以及管制终 端模块；监视数据融合模块用于实现空管雷达监视数据与自动相关监视数据的融合，为管 制终端模块提供实时航迹信息；
[0007] 所述管制终端模块包括以下子模块：
[0008] 飞行前无冲突4D航迹生成模块，根据飞行计划和世界区域预报系统的预报数据， 建立航空器动力学模型，然后依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成航 空器无冲突4D航迹；
[0009] 飞行中短期4D航迹生成模块，依据监视数据融合模块提供的实时航迹信息，利用 隐马尔科夫模型，推测未来一定时间窗内的航空器4D轨迹；
[0010] 实时飞行冲突监控与告警模块，用于建立从航空器的连续动态到离散冲突逻辑的 观测器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能违 反空中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供及 时的告警信息；
[0011] 所述空中交通管制系统的飞行冲突预警方法包括如下几个步骤：
[0012] 步骤A、飞行前无冲突4D航迹生成模块根据飞行计划和世界区域预报系统的预报 数据，建立航空器动力学模型，并依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成 航空器无冲突4D航迹；
[0013] 步骤B、监视数据融合模块将空管雷达监视数据与自动相关监视数据进行融合， 生成航空器实时航迹信息并提供给管制终端模块；管制终端模块中的飞行中短期4D航迹 生成模块依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨 迹；所述依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨 迹的具体实施过程如下：
[0014] 步骤B6、对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列 X = [X1, X2,…，Xn]和y=Iiy1,y2，…，yn]，采用一阶差分方法对其进行处理获取新的航空器 离散位置序列ΛX= [ΛX1,Λχ2,…，ΛXn-J和Λy= [Λy"Λy2,…，ΛYn-J,其中 Δxb=χb+「xb，Δyb=yb+「yb(b= 1，2，…，n-1);
[0015] 步骤B7、对航空器轨迹数据聚类，对处理后新的航空器离散二维位置序列Λx和 Λy，通过设定聚类个数Μ'，采用遗传聚类算法分别对其进行聚类；
[0016] 步骤B8、对聚类后的航空器轨迹数据利用隐马尔科夫模型进行参数训练，通过将 处理后的航空器运行轨迹数据ΛX和Λy视为隐马尔科夫过程的显观测值，通过设定隐状 态数目Ν'和参数更新时段ζ'，依据最近的Τ'个位置观测值并采用B-W算法滚动获取最新 隐马尔科夫模型参数λ' ；
[0017] 步骤Β9、依据隐马尔科夫模型参数，采用Viterbi算法获取当前时刻观测值所对 应的隐状态q;
[0018] 步骤B10、通过设定预测时域h'，基于航空器当前时刻的隐状态q，获取未来时段 航空器的位置预测值O;
[0019] 步骤C、实时飞行冲突监控与告警模块建立从航空器的连续动态到离散冲突逻辑 的观测器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能 违反空中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供 及时的告警信息。
[0020] 进一步的，步骤B中，所述聚类个数M'的值为4,隐状态数目Ν'的值为3,参数更 新时段ζ'为30秒，Τ'为10,预测时域h'为300秒。
[0021] 进一步的，步骤B的B8具体是指：由于所获得的航迹序列数据长度是动态变化 的，为了实时跟踪航空器航迹的状态变化，有必要在初始航迹隐马尔科夫模型参数λ' = (π，Α，Β)的基础上对其重新调整，以便更精确地推测航空器在未来某时刻的位置；每隔时 段ζ'，依据最新获得的τ'个观测值（〇1，〇2，…，〇τ，）对航迹隐马尔科夫模型参数λ' =O，A,B)进行重新估计。
[0022]步骤B的BlO具体是指：每隔时段Φ，根据最新获得的隐马尔科夫模型参数λ' = (π，Α，Β)和最近H个历史观测值（〇1，〇2，…，〇Η)，基于航空器当前时刻的隐状态q，通过设定 预测时域h'，在时刻t获取航空器在未来时段h'的位置预测值0。
[0023] 更进一步的，时段Φ为4秒。
[0024] 进一步的，所述步骤A的航空器无冲突4D航迹按照以下方法生成：
[0025] 步骤Α1、进行航空器状态转移建模，根据飞行计划中航空器的飞行高度剖面，建立 单个航空器在不同航段转移的Petri网模型：E= (g，G，Pre，Post，m)为航空器阶段转移模 型，其中g表示飞行航段，G表示垂直剖面中飞行状态参数的转换点，Pre和Post分别表示 航段和航路点的前后向连接关系，^Z+表示航空器所处的飞行阶段；
[0026] 步骤A2、建立航空器全飞行剖面混杂系统模型如下，
[0027]vH=k(vcas,Mach,hp,tL0C),
[0028]vGS=μ(vCAS,Mach,hp,tL0C,vffS,a),
[0029] 其中veAS为校正空速，Mach为马赫数，hp为气压高度，a为风向预报与航路的夹 角，Vws为风速预报值，t为温度预报值，VA高度变化率，V<^为地速；
[0030] 步骤A3、采用混杂系统仿真的方式推测求解航迹：采用将时间细分的方法， 利用状态连续变化的特性递推求解任意时刻航空器在某一飞行阶段距参考点的航程 J(X) =J。+ 和高度= 其中Jtl为初始时刻航空器距参考点 的航程，Λτ为时间窗的数值，J(T)为τ时刻航空器距参考点的航程，Iitl为初始时刻航 空器距参考点的高度，1ι(τ)为τ时刻航空器距参考点的高度，由此可以推测得到单航空 器的4D航迹；
[0031] 步骤Α4、对多航空器耦合模型实施无冲突调配：根据两航空器预达交叉点的时 间，按照空中交通管制原则，对交叉点附近不满足间隔要求的航空器4D航迹进行二次规 划，得到无冲突4D航迹。
[0032] 进一步的，所述步骤B中监视数据融合模块将空管雷达监视数据与自动相关监视 数据进行融合，生成航空器实时航迹信息，具体按照以下方法：
[0033] 步骤Β1、将坐标单位和时间统一；
[0034] 步骤Β2、采用最邻近数据关联算法将属于同一个目标的点相关联，提取目标航迹； 步骤Β3、将分别从自动相关监视系统和空管雷达提取的航迹数据从不同的时空参
[0035] 考坐标系统变换、对准到管制终端统一的时空参考坐标系统；
[0036] 步骤M、计算两条航迹的相关系数，若相关系数小于某一预设阈值，则认为两条航 迹不相关；否则该两条航迹相关，可以进行融合；
[0037] 步骤Β5、对相关的航迹进行融合。
[0038] 更进一步的，所述步骤Β5中对相关的航迹进行融合，采用基于采样周期的加权平 均算法，其加权系数根据采样周期和信息精度确定，再利用加权平均算法将与之相关的自 动相关监视航迹和空管雷达航迹融合为系统航迹。
[0039] 进一步的，所述步骤C的具体实施过程如下：
[0040] 步骤CU构造基于管制规则的冲突超曲面函数集：建立超曲面函数集用以反映系 统的冲突状况，其中，冲突超曲面中与单一航空器相关的连续函数为第I型超曲 面，与两架航空器相关的连续函数\ :YXX^R为第II型超曲面；
[0041] 步骤C2、建立由航空器连续状态至离散冲突状态的观测器：需要根据管制规范建 立观测器，观测系统系统穿越超曲面而产生的冲突事件，以便控制器做出相应的控制决策 指令；观测器ξ用于观测系统中航空器位置的连续变化而产生冲突事件，称为 第I型观测器，4S为第II型观测器；
[0042] 步骤C3、设计从冲突到冲突解脱手段的离散监控器，该离散监控器可描述为函数 βD，其中S是观测器观测向量展成的空间，D是所有决策向量d展成的空间；当观测 器的离散观测向量表明某一非期望的状态出现时，立刻发出相应的告警。
[0043] 本发明具有积极的效果：（1)本发明的一种空中交通管制系统的飞行冲突预警方 法在航空器实时轨迹推测过程中，融入了随机因素的影响，所采用的滚动轨迹推测方案能 够及时提取外界随机因素的变化状况，提高了航空器轨迹推测的准确性。
[0044] (2)本发明的一种空中交通管制系统的飞行冲突预警方法对飞行冲突的预警效果 较好，可有效、准确、实时地预测航空器的轨迹并预测飞行冲突。
[0045] (3)本发明的一种空中交通管制系统的飞行冲突预警方法对飞行剖面的推算和航 迹预测精度高，进而使得冲突化解能力和自动化水平提高，降低了管制员的工作负荷。
【专利附图】
【附图说明】
[0046] 图1为飞行前无冲突4D航迹生成方法流程示意图；
[0047] 图2为飞行中短期4D航迹推测方法流程示意图；
[0048] 图3为航空器航迹冲突监控与告警方法流程示意图。
【具体实施方式】
[0049] (实施例1)
[0050] 本实施例的基于4D航迹运行的空中交通管制系统，包括机载终端模块101、数据 通信模块102、监视数据融合模块103以及管制终端模块104。以下对各部分的具体实施方 式分别进行详细描述。
[0051] 1.机载终端模块
[0052] 机载终端模块101是飞行员获取地面管制指令、参考4D航迹，以及输入飞行意图 的界面，同时还是采集当前航空器位置数据的接口。
[0053]其具体实施方案如下：
[0054] 机载终端模块101接收如下的信息输入：（I)ADS-B信息采集单元201通过机载 GPS采集的航空器位置向量、速度向量，以及本航空器的呼号，编码后通过信息及数据传递 给机载数据通信模块102 ; (2)航空器驾驶员需要将与地面管制指令不一致的飞行意图，通 过人机输入界面，以及约定的地面管制员可以识别的形式通过信息及数据传递给机载数据 通信模块102。另外机载终端模块101实现如下的信息输出：（1)通过终端显示屏幕，接收 和显示飞行员可以识别的飞行管制指令；（2)接收和显示地面管制终端飞行前生成的无冲 突4D航迹，以及当地面管制终端探测到冲突后计算的最优解脱4D航迹。
[0055] 2.数据通信模块
[0056] 数据通信模块102可实现空地双向数据通信，实现机载实时位置数据和飞行意图 数据单元202的下行传输和地面管制指令单元203,以及参考4D航迹单元204的上行传输。
[0057]其具体实施方案如下：
[0058] 下行数据通信：机载终端101通过机载二次雷达应答机将航空器识别标志和4D 位置信息，以及其他附加数据，如飞行意图、飞行速度、气象等信息传输给地面二次雷达 (SSR)，二次雷达接收后对数据报文进行解析，并传输给中央数据处理组件301解码，通过 指令航迹数据接口传输到管制终端104;上行数据通信：地面管制终端104通过指令航迹数 据接口，经中央数据处理组件301编码后，地面二次雷达的询问机将将地面管制指令或参 考4D航迹信息传递并显示在机载终端101。
[0059] 3.监视数据融合模块
[0060] 监视数据融合模块103实现空管雷达监视与自动相关监视ADS-B数据的融合，为 管制终端模块104中的飞行中短期4D航迹生成子模块和实时飞行冲突监控与告警子模块 提供实时航迹信息。
[0061] 其具体实施方案如下：
[0062] (1)在预处理阶段将坐标单位和时间统一，假设分别从ADS-B和空管雷达中提取 的数据是一系列离散点的坐标（如经度、炜度、海拔高度）、各点对应采集时间；（2)采用最 邻近数据关联算法将属于同一个目标的点相关联，提取目标航迹；(3)将分别从ADS-B和空 管雷达中提取的航迹数据从不同的时空参考坐标系统变换、对准到管制终端统一的时空参 考坐标系统；(4)计算两条航迹的相关系数，若相关系数小于某一预设阈值，则认为两条航 迹不相关，否则该两条航迹相关，可以进行融合；（5)对相关的航迹进行融合。由于ADS-B和 空管雷达的精度和采样周期不同，本系统采用基于采样周期的加权平均算法，其加权系数 根据采样周期和信息精度确定，再利用加权平均算法将与之相关的ADS-B航迹和空管雷达 航迹融合为系统航迹。
[0063] 4.管制终端模块
[0064] 管制终端模块104包括飞行前无冲突4D航迹生成、飞行中短期4D航迹生成、实时 飞行冲突监控与告警这3个子模块。
[0065] (1)飞行前无冲突4D航迹生成
[0066] 根据飞行数据处理系统（FDP)得到的飞行计划和世界区域预报系统（WAFS)发布 的风、温度的GRIB格点预报数据，对空中交通系统建立层次化的混杂系统模型，通过系统 在安全状态的演化，描述状态演化的时间轨迹，生成航空器航迹。
[0067] 如图1所示，其具体实施过程如下：
[0068] 首先，进行航空器状态转移建模。航空器沿航迹飞行的过程表现为在航段之间动 态切换过程，根据飞行计划中航空器的飞行高度剖面，建立单个航空器在不同航段转移的 Petri网模型：E= (g，G,Pre,Post,m)为航空器阶段转移模型，其中g表示飞行航段，G表 示垂直剖面中飞行状态参数（包括空速、高度、构型）的转换点，Pre和Post分别表示航段 和航路点的前后向连接关系，《I^Z1表示航空器所处的飞行阶段。
[0069]其次，建立航空器全飞行剖面混杂系统模型。航空器在单个航段内的飞行视为连 续过程，依据质点能量模型，推导航空器在不同的运行阶段同气象条件下的航空器动力学 方程，νΗ=κ(VCAS,Mach,hp,tLQC)，Vgs=μ(VCAS,Mach,hp,tLQC,vws,α)，其中VcAS为校正空速， Mach为马赫数，hp为气压高度，α为风向预报与航路的夹角，Vws为风速预报值，为温度 预报值，vH为高度变化率，V<^为地速。
[0070] 然后，采用混杂系统仿真的方式推测求解航迹。采用将时间细分的方法，利 用状态连续变化的特性递推求解任意时刻航空器在某一飞行阶段距参考点的航程 = +f)Λ和高度Α(τ) =Α。+fn?⑴也，其中Jtl为初始时刻航空器距参考点 JO JO 的航程，ΛT为时间窗的数值，J(T)为τ时刻航空器距参考点的航程，Iitl为初始时刻航 空器距参考点的高度，1ι(τ)为τ时刻航空器距参考点的高度，由此可以推测得到单航空 器的4D航迹。
[0071] 最后，对多航空器耦合模型实施无冲突调配。根据两航空器预达交叉点的时间，按 照空中交通管制原则，对交叉点附近不满足间隔要求的航空器4D航迹进行二次规划，得到 无冲突4D航迹。
[0072] (2)飞行中短期4D航迹生成
[0073] 依据管制雷达和自动相关监视系统ADS-B实施融合后获得航空器实时航迹数据， 利用隐马尔科夫模型，推测未来5分钟时间窗内的航空器4D轨迹。
[0074] 如图2所示，其具体实施过程如下：
[0075] 首先，对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列X =[XpX2,…，Xn]和y=Iiy^y2，…，yn]，采用一阶差分方法对其进行处理获取新的航空器 离散位置序列ΛX= [ΛX1,Λχ2,…，ΛXn-J和Λy= [Λy"Λy2,…，ΛYn-J,其中 AXb-Xb+l_Xb，AYb-yb+l_yb(b- 1，2，...，n_l)。
[0076] 其次，对航空器轨迹数据聚类。对处理后新的航空器离散二维位置序列ΛX和 Λy，通过设定聚类个数Μ'，采用遗传聚类算法分别对其进行聚类。
[0077] 然后，对聚类后的航空器轨迹数据利用隐马尔科夫模型进行参数训练。通过将处 理后的航空器运行轨迹数据ΛX和Λy视为隐马尔科夫过程的显观测值，通过设定隐状态 数目Ν'和参数更新时段ζ'，依据最近的Τ'个位置观测值并采用B-W算法滚动获取最新 隐马尔科夫模型参数λ'：由于所获得的航迹序列数据长度是动态变化的，为了实时跟踪 航空器航迹的状态变化，有必要在初始航迹隐马尔科夫模型参数λ' = (π，Α，Β)的基础上 对其重新调整，以便更精确地推测航空器在未来某时刻的位置。每隔时段ζ'，依据最新获 得的Τ'个观测值（〇1，〇2，…，〇Τ，）对航迹隐马尔科夫模型参数λ' = (π，Α，Β)进行重新估 计。
[0078] 再而，依据隐马尔科夫模型参数，采用Viterbi算法获取当前时刻观测值所对应 的隐状态q。
[0079] 最后，每隔时段tP，根据最新获得的隐马尔科夫模型参数λ' = (π，A，B)和最近 H个历史观测值（〇1，O2，…，〇H)，基于航空器当前时刻的隐状态q，通过设定预测时域h'，在 时刻t获取航空器在未来时段h'的位置预测值0。
[0080] 所述聚类个数M'的值为4,隐状态数目Ν'的值为3,参数更新时段ζ'为30秒， Τ'为10,预测时域h'为300秒，时段φ为4秒。
[0081] (3)实时飞行冲突监控与告警
[0082] 当系统有可能出现违反安全状态集的状态时，通过控制器实施状态监控，对航空 器实施有效的管制措施，避免飞行冲突的发生。
[0083] 如图3所示，其具体实施过程如下：
[0084] 首先，构造基于管制规则的冲突超曲面函数集。空中交通管制约束的违反都可以 视为被控对象（管制空域飞行的多架航空器）构成系统穿越超曲面而产生的事件，建立超 曲面函数集用以反映系统的冲突状况。其中，冲突超曲面中与单一航空器相关的连续函数 /I1:Z4i?为第〗型超曲面，而将与两架航空器相关的连续函数/?:ZxB为第π型超 曲面。
[0085] 然后，建立由航空器连续状态至离散冲突状态的观测器。需要根据管制规范建立 观测器，观测系统系统穿越超曲面而产生的冲突事件，以便控制器做出相应的控制决策指 令。观测器ξ用于观测系统中航空器位置的连续变化而产生冲突事件，称为第 I型观测器，ξ2 ?为第II型观测器。
[0086] 最后，设计从冲突到冲突解脱手段的离散监控器。当观测器的离散观测向量表明 某一非期望的状态出现时，立刻发出相应的告警。该离散监控器可描述为函数βJ 其中S是观测器观测向量展成的空间，D是所有决策向量d展成的空间。
[0087] 显然，上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的 实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其 它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发 明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。
【权利要求】
1. 一种空中交通管制系统的航空器轨迹预测方法，所述空中交通管制系统包括机载终 端模块、数据通信模块、监视数据融合模块以及管制终端模块；监视数据融合模块用于实现 空管雷达监视数据与自动相关监视数据的融合，为管制终端模块提供实时航迹信息；其特 征在于： 所述管制终端模块包括以下子模块： 飞行前无冲突4D航迹生成模块，根据飞行计划和世界区域预报系统的预报数据，建立 航空器动力学模型，然后依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成航空器 无冲突4D航迹； 飞行中短期4D航迹生成模块，依据监视数据融合模块提供的实时航迹信息，利用隐马 尔科夫模型，推测未来一定时间窗内的航空器4D轨迹； 实时飞行冲突监控与告警模块，用于建立从航空器的连续动态到离散冲突逻辑的观测 器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能违反空 中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供及时的 告警信息；所述空中交通管制系统的飞行冲突预警方法包括如下几个步骤： 步骤A、飞行前无冲突4D航迹生成模块根据飞行计划和世界区域预报系统的预报数 据，建立航空器动力学模型，并依据飞行冲突耦合点建立航迹冲突预调配理论模型，生成航 空器无冲突4D航迹； 步骤B、监视数据融合模块将空管雷达监视数据与自动相关监视数据进行融合，生成航 空器实时航迹信息并提供给管制终端模块；管制终端模块中的飞行中短期4D航迹生成模 块依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨迹；所 述依据航空器实时航迹信息和历史航迹信息推测未来一定时间窗内的航空器4D轨迹的具 体实施过程如下： 步骤B6、对航空器轨迹数据预处理，依据所获取的航空器原始离散二维位置序列X= [X1,x2, . . .，xn]和y=Iiy1,y2, . . .，yn]，采用一阶差分方法对其进行处理获取新的航空器离 散位置序列AX= [AX1,Ax2, ? --，Axn_J和Ay= [Ay"Ay2, ? --，Ayn_J,其中 Axb=xb+1-xb,Ayb=yb+1_yb(b=I, 2, ? ? ? ,n-l); 步骤B7、对航空器轨迹数据聚类，对处理后新的航空器离散二维位置序列Ax和Ay， 通过设定聚类个数M'，采用遗传聚类算法分别对其进行聚类； 步骤B8、对聚类后的航空器轨迹数据利用隐马尔科夫模型进行参数训练，通过将处理 后的航空器运行轨迹数据AX和Ay视为隐马尔科夫过程的显观测值，通过设定隐状态数 目N'和参数更新时段G'，依据最近的T'个位置观测值并采用B-W算法滚动获取最新隐马 尔科夫模型参数A' ； 步骤B9、依据隐马尔科夫模型参数，采用Viterbi算法获取当前时刻观测值所对应的 隐状态q; 步骤B10、通过设定预测时域h'，基于航空器当前时刻的隐状态q，获取未来时段航空 器的位置预测值〇 ; 步骤C、实时飞行冲突监控与告警模块建立从航空器的连续动态到离散冲突逻辑的观 测器，将空中交通系统的连续动态映射为离散观测值表达的冲突状态；当系统有可能违反 空中交通管制规则时，对空中交通混杂系统的混杂动态行为实施监控，为管制员提供及时 的告警信息。
2. 根据权利要求1所述的一种空中交通管制系统的飞行冲突预警方法，其特征在于： 步骤B中，所述聚类个数M'的值为4,隐状态数目N'的值为3,参数更新时段C为30秒， T'为10,预测时域h'为300秒。
3. 根据权利要求1或2所述的一种空中交通管制系统的飞行冲突预警方法，其特征在 于：步骤B的B8具体是指：由于所获得的航迹序列数据长度是动态变化的，为了实时跟踪 航空器航迹的状态变化，有必要在初始航迹隐马尔科夫模型参数X' =(IA，B)的基础上 对其重新调整，以便更精确地推测航空器在未来某时刻的位置；每隔时段G'，依据最新获 得的T'个观测值（〇1，〇2,. . .，〇T，）对航迹隐马尔科夫模型参数A' = (31，A，B)进行重新估 计； 步骤B的BlO具体是指：每隔时段？，根据最新获得的隐马尔科夫模型参数X' = (，A，B)和最近H个历史观测值（〇1，〇2, . . .，〇H)，基于航空器当前时刻的隐状态q，通过设 定预测时域h'，在时刻t获取航空器在未来时段h'的位置预测值0。
4. 根据权利要求1至3之一所述的一种空中交通管制系统的飞行冲突预警方法，其特 征在于：时段9为4秒。
5. 根据权利要求1至4之一所述的一种空中交通管制系统的飞行冲突预警方法，其特 征在于：所述步骤A的航空器无冲突4D航迹按照以下方法生成： 步骤A1、进行航空器状态转移建模，根据飞行计划中航空器的飞行高度剖面，建立单个 航空器在不同航段转移的Petri网模型：E= (g，G，Pre，Post，m)为航空器阶段转移模型， 其中g表示飞行航段，G表示垂直剖面中飞行状态参数的转换点，Pre和Post分别表示航段 和航路点的前后向连接关系，Z+表示航空器所处的飞行阶段； 步骤A2、建立航空器全飞行剖面混杂系统模型如下， vH-KCAS；Mach,hp,tL0C), vgs-t1CAS)Mach,hp,tL0C,vffS,a), 其中Vffls为校IH空速，Mach为马赫数，hp为气压高度，a为风向预报与航路的夹角，vws 为风速预报值，为温度预报值，V高度变化率，V<^为地速； 步骤A3、采用混杂系统仿真的方式推测求解航迹：采用将时间细分的方法，利 用状态连续变化的特性递推求解任意时刻航空器在某一飞行阶段距参考点的航程 卜人+ 和高度_) = &+JTb(T)?^其中Jci为初始时刻航空器距参考点 的航程，AT为时间窗的数值，J(T)为T时刻航空器距参考点的航程，Iltl为初始时刻航 空器距参考点的高度，h(T)为T时刻航空器距参考点的高度，由此可以推测得到单航空 器的4D航迹； 步骤A4、对多航空器耦合模型实施无冲突调配：根据两航空器预达交叉点的时间，按 照空中交通管制原则，对交叉点附近不满足间隔要求的航空器4D航迹进行二次规划，得到 无冲突4D航迹。
6. 根据权利要求1至5之一所述的一种空中交通管制系统的飞行冲突预警方法，其特 征在于：所述步骤B中监视数据融合模块将空管雷达监视数据与自动相关监视数据进行融 合，生成航空器实时航迹信息，具体按照以下方法： 步骤Bl、将坐标单位和时间统一； 步骤B2、采用最邻近数据关联算法将属于同一个目标的点相关联，提取目标航迹； 步骤B3、将分别从自动相关监视系统和空管雷达提取的航迹数据从不同的时空参考坐 标系统变换、对准到管制终端统一的时空参考坐标系统； 步骤B4、计算两条航迹的相关系数，若相关系数小于某一预设阈值，则认为两条航迹不 相关；否则该两条航迹相关，可以进行融合； 步骤B5、对相关的航迹进行融合。
7. 根据权利要求1至6之一所述的一种空中交通管制系统的飞行冲突预警方法，其特 征在于：所述步骤B5中对相关的航迹进行融合，采用基于采样周期的加权平均算法，其加 权系数根据采样周期和信息精度确定，再利用加权平均算法将与之相关的自动相关监视航 迹和空管雷达航迹融合为系统航迹。
8. 根据权利要求1至7之一所述的一种空中交通管制系统的飞行冲突预警方法，其特 征在于：所述步骤C的具体实施过程如下： 步骤C1、构造基于管制规则的冲突超曲面函数集：建立超曲面函数集用以反映系统的 冲突状况，其中，冲突超曲面中与单一航空器相关的连续函数hB&为第I型超曲面， 与两架航空器相关的连续函数心^〃为第II型超曲面； 步骤C2、建立由航空器连续状态至离散冲突状态的观测器：需要根据管制规范建立观 测器，观测系统系统穿越超曲面而产生的冲突事件，以便控制器做出相应的控制决策指令； 观测器I用于观测系统中航空器位置的连续变化而产生冲突事件，称S为第I型 观测器，^ : ZBS为第II型观测器； 步骤C3、设计从冲突到冲突解脱手段的离散监控器，该离散监控器可描述为函数pD，其中S是观测器观测向量展成的空间，D是所有决策向量d展成的空间；当观测 器的离散观测向量表明某一非期望的状态出现时，立刻发出相应的告警。
【文档编号】G08G5/04GK104504942SQ201510007919
【公开日】2015年4月8日 申请日期:2015年1月7日 优先权日:2015年1月7日
【发明者】韩云祥, 赵景波, 李广军 申请人:江苏理工学院