一种基于聚类分析的扇区运行性能综合检测方法及系统的利记博彩app
【技术领域】
[0001] 本发明涉及航空领域，尤指一种空中交通管制扇区运行性能综合检测方法及系 统。
【背景技术】
[0002] 伴随着航空运输业的发展，为了保证各类飞行活动的安全和有序，空中交通管制 服务应运而生并不断得到发展完善，至20世纪80年代趋于成熟。现代空中交通管制服务 的主要内容是：空中交通管制员（简称为"管制员"，下同）依托现代通信、导航、监视技术， 对所辖航空器实施管理和控制，协调和指导其运动路径和模式，以防止空中航空器与航空 器相撞及在机场机动区内航空器与障碍物相撞，维护和加快空中交通的有序流动。空中交 通管制扇区（简称为"管制扇区"，下同）是空中交通管制（简称为"管制"，下同）的基本 空间单元。一般情况下，为航空器提供空中交通管制服务的空域被划设为若干管制扇区，每 个管制扇区对应一个管制员工作席位。管制扇区运行性能是管制扇区内航空器运行态势的 技术性指标提炼，既反映管制员对所辖管制扇区提供管制服务的质量和水准，又反映特定 管制空域使用效能。因此，对管制扇区运行性能的有效检测是调整管制运行策略、优化管制 空域结构的基础和前提。
[0003] 比如，公开号为CN104332073A的专利文献在2015-02-04公开了一种智能空中交 通管制系统，包括数据接收接口模块、数据处理中心、应急超控模块、电子飞行计划显示模 块、航班监控显示模块和航班控制指令发送模块。通过接收飞机准确的实时飞行信息，数据 处理中心即可展开飞行计划的精确计算，并优化整个飞行队列，如缩短飞机间距，合理调配 对应跑道的起降顺序，实时更改着陆角，曲线进近等，有效提高整个飞行队列的效率，加强 安全性。最终，整个飞行计划将显示在管制员的监控显示器上。航空管制员可以随时通过 显示器观察飞行队列起飞着陆的优先顺序，并通过更改飞机起飞/着陆航迹，航速等参数， 调整队列顺序。从而提尚机场通勤率，提尚管制员的工作效率。
[0004] 但即便如此，目前针对空中交通管制扇区运行性能的研究较少，大部分研究主要 体现在以下几个孤立方面：（1)空中交通流密度，分为战略和战术两层面，其中前者主要 体现为空域复杂性指标，后者主要体现为管制单元空中交通拥挤程度判定。目前，空中交 通流密度指标在应用上仍以管制单元的航空器架次统计作为主要呈现。（2)管制运行安 全性能，包括定量和定性两方面。定量方面，国际民航组织（ICAO)依据碰撞风险分析制定 的总的安全目标等级（TLS)是1.5X10-8次致命飞行事故/飞行小时，而我国民航空管系 统根据危险接近风险分析将事故征候万架次率作为关键安全指标。定性方面，ICAO推荐 采用威胁差错管理（Threat and Error Management，TEM)或日常运行安全监测（Normal Operations Safety Survey，N0SS)方法，实施定性的管制运行安全性能评价。国内学者围 绕人、机、环、管理等4类因素分别建立了安全风险评估指标体系，并开展了指标权重分析。 (3)管制运行效率性能，主要围绕航班延误指标方面。目前，国外航班延误统计指标涉及延 误架次率及延误时间。我国民航欠缺航班延误时间的细化统计，在航班延误原因界定、统计 指标设计、统计方法及流程等方面亟待改善。（4)管制员工作负荷，是管制扇区容量评估的 重要考量。国外学者从生理/行为特征、主观测评、工作细分的角度，分别提出了电击皮肤 的反应、心率、心电图、血压、体液等生理指标，设备操作次数、陆空通话时间记录等行为指 标；ATWIT技术、NASA - TLX量表、SWAT量表和MCH法等主观测评技术；DORATASK、MBB法、 RAMS法等衡量管制员工作时间的方法。国内学者发展了主观测评方法，提出了基于可拓学 的管制员工作负荷评价模型。
[0005] 例如，公开号为CN104636890A的专利文献在2015-05-20公开了一种空中交通管 制员工作负荷测量方法，包括：步骤A :确定管制负荷测量指标，该管制负荷测量指标包括 眼动指标和语音指标；步骤B :实时记录各眼动指标对应的眼动指标数据，以及各语音指标 对应的语音指标数据；步骤C :对记录的眼动指标数据进行因子分析，计算出眼动指标数据 的眼动综合因子；步骤D :以眼动综合因子和语音指标为输入因素，管制综合指标值为输出 因素，建立管制负荷回归模型。该方法能够实时、无干扰的测量管制员的综合指标，实用性 强。
[0006] 但该方法参考的数据比较局限，指标维度单一，不够全面、综合，该测量方法仅利 用局限的数据片面的对管制员工作负荷进行了测量，存在一定的局限性，预测可靠度不高。
【发明内容】
[0007] 本发明提供一种可以提高扇区性能的检测结果的可靠性的空中交通管制扇区运 行性能检测方法和系统。
[0008] 本发明的目的是通过以下技术方案来实现的：
[0009] 一种扇区运行性能综合检测方法，包括步骤：
[0010] 步骤1 :输入历史扇区运行性能指标数据，包括扇区通行性检测指标、扇区复杂性 检测指标、扇区安全性检测指标、扇区经济性检测指标和管制员工作负荷检测指标；
[0011] 步骤2 :基于聚类分析模型处理历史扇区运行性能指标数据，建立扇区运行性能 档案；
[0012] 步骤3 :输入当前扇区运行性能指标数据，对照扇区运行性能档案寻找最大相似 类，以此判断扇区运行性能等级，并计算得到最小相似度；
[0013] 步骤4 :根据最小相似度确定是否触发扇区异常运行响应告警。
[0014] 进一步的，所述步骤1中的所述扇区通行性检测指标为扇区流量、扇区航行里程、 扇区航行时间和扇区交通流密度；
[0015] 所述扇区复杂性检测指标为扇区航空器爬升次数、扇区航空器下降次数、扇区航 空器改速次数、扇区航空器改航次数；
[0016] 所述扇区安全性检测指标为扇区短期冲突告警频率和扇区最低安全高度告警频 率；
[0017] 所述扇区经济性检测指标为扇区饱和度、扇区排队长度、扇区航空器延误架次率、 扇区航空器延误时间、扇区航空器平均延误时间；
[0018] 所述管制员工作负荷检测指标为陆空通话信道占用率、陆空通话次数。
[0019] 进一步的，在所述步骤1与步骤2之间还包括对历史扇区运行性能指标数据进行 无量纲化处理，具体过程如下：
[0020] 令第i个样本的第j项指标的实际值为Xli ,，yu为无量纲化处理后的指标值，《 为第j项指标的均值，S]为第j项指标的标准差，则对于正向指标，
，对于逆 向指标，先取其倒数或取负使其正向化，再使用上述公式对其进行无量纲化处理。
[0021] 进一步的，所述步骤2具体包括如下过程：
[0022] 步骤2. 1 :根据预设的聚类数目k对历史扇区运行性能指标数据点进行分类，形成 k个初始类，并确定每一初始类中心结点；
[0023] 步骤2. 2 :计算历史扇区运行性能指标数据中的每个数据点到k个初始类中心结 点之间的欧几里德距离，并根据距离类中心结点最短的原则将所有的数据进行重新分类形 成k个类；
[0024] 步骤2. 3 :依次计算各类变量的均值，并以均值作为新的k个类的类中心结点；
[0025] 步骤2. 4 :计算新确定的类中心结点距上一个类中心结点之间的偏移量，当最大 偏移量大于设定值时，返回分类模块，否则聚类结束得到k个类，从而得到扇区正常运行性 能档案。
[0026] 进一步的，所述步骤3具体包括如下步骤：
[0027] 步骤3. 1 :实时数据输入与标准化：
[0028] 获取当前扇区运行性能的上述17项指标数据r = Irj, j = 1，2, . . .，17}，tj为无 量纲化处理后的第j项指标值，^为第j项指标的均值，S]为第j项指标的标准差，则对于 正向指标，
对于逆向指标，先取其倒数或取负使其正向化，再使用上述公式对 其进行无量纲化处理；
[0029] 其中均值和标准差从扇区运行性能档案中提取；
[0030] 步骤3. 2 :最大相似类寻找：
[0031] 计算检测数据无量纲化指标值t与扇区运行性能档案中各个类的类中心点(；的 相似度sim(t, Ck)，并求得与扇区运行性能档案中k个类中心点的最大相似度simm(t, Ck)及 其类Class (k)，其中sim(t, Ck)可记为：
[0033] 其中，t,为检测数据无量纲化后第j项指标值，c ki ]为第k个类中心点的第j项无 量纲化指标值；
[0034] 当前扇区运行性能等级即以扇区运行性能档案中的最大相似类确定。
[0035] 步骤3. 3 :最小相似度计算：
[0036] 计算t与所属Class (k)中各样本无量纲化数据的相似度，得到最大相似类中的最 小相似度min Sim。
[0037] 进一步的，步骤4 :若计算得到的min Sim与预设的阈值进行比较，min Sim超过预 先设定的相似度阈值，则判断扇区运行性能情况异常，并进行扇区异常运行响应告警。
[0038] 一种空中交通管制扇区运行性能检测系统，包括：
[0039] 扇区运行性能档案创建模块，根据聚类分析后的历史扇区运行性能指标数据创建 扇区运彳丁性能档案；
[0040] 检测模块：根据输入的当前扇区运行性能指标数据，对照扇区运行性能档案计算 得到最大相似类和最小相似度，并根据最小相似度判断扇区运行性能是否异常。
[0041] 进一步的，所述扇区运行性能档案创建模块包括聚类分析单元和创建单元，所述 聚类分析单元用于对输入的，经过无量纲化处理的历史扇区运行性能指标数据进行聚类分 析处理，所述创建模块用于根据聚类分析处理得到的结果创建扇区运行性能档案；
[0042] 所述检测模块包括相似度处理单