CN113362604B - 一种基于扇区动态交通特征的管制员工作负荷评估方法 - Google Patents
一种基于扇区动态交通特征的管制员工作负荷评估方法 Download PDFInfo
- Publication number
- CN113362604B CN113362604B CN202110812169.7A CN202110812169A CN113362604B CN 113362604 B CN113362604 B CN 113362604B CN 202110812169 A CN202110812169 A CN 202110812169A CN 113362604 B CN113362604 B CN 113362604B
- Authority
- CN
- China
- Prior art keywords
- sector
- sample
- index
- traffic
- classification
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供了一种基于扇区动态交通特征的管制员工作负荷评估方法,包括构建扇区动态交通特征度量指标集合,并以等长时间片为单位构建分时交通动态交通样本数据集;采用高斯混合聚类算法进行交通样本聚类,实现对样本特征的识别与分类;针对每个分类,采用熵权法对样本包含的各个交通指标进行客观赋权,实现交通特征度量向管制员工作负荷度量的演进;基于随机森林算法建立扇区交通运行场景分类模型,以聚类结果为分类标签,对待评估样本进行分类决策,并依托分类结果进行管制员工作负荷度量。本方法消除了现有管制员工作负荷评估方法中人为因素对于评估结果干扰性大的问题,提升了评估结果的客观性与精确。
Description
技术领域
本发明涉及一种基于扇区动态交通特征的管制员工作负荷评估方法。
背景技术
随着空中交通流量不断增长,扇区交通复杂性不断提高,如何客观度量管制员工作负荷是当下的热点问题。准确合理的评估管制员的工作负荷对于识别运行潜在安全风险,提升管制员工作效率具有重要意义。
目前,管制员工作负荷评估方法大致分有三类:一是利用管制员生理反应,通过测量管制员通话、心率等行为或心理特征,评估管制员工作负荷;二是通过调查问卷等主观询问的方式由管制员自评工作负荷;三是通过管制模拟机的训练数据统计管制员工作负荷。其中测量管制员生理特征的方法会干扰管制员的管制行为,且管制员个体间差异性较大,操作难度大;调查问卷等主观询问的方法无法量化度量管制员的工作负荷,各负荷分级之间界限模糊。通过管制模拟机评估的方式评估周期长,成本较高。此外,上述方法缺乏针对特定场景的战术级管制员工作负荷评估能力,因此导致评估方法的应用领域较为狭隘。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种基于扇区动态交通特征的管制员工作负荷评估方法,旨在依托扇区交通运行特征,实现对不同运行场景下管制员工作负荷的客观度量。所述方法包括构建扇区动态交通特征度量指标集合,并以等长时间片为单位构建分时交通动态交通样本数据集;采用高斯混合聚类算法进行交通样本聚类,实现对样本特征的识别与分类;针对每个分类,采用熵权法对样本包含的各个交通指标进行客观赋权,实现交通特征度量向管制员工作负荷度量的演进;基于随机森林算法建立扇区交通运行场景分类模型,以聚类结果为分类标签,对待评估样本进行分类决策,并依托分类结果进行管制员工作负荷度量。本方法消除了现有管制员工作负荷评估方法中人为因素对于评估结果干扰性大的问题,提升了评估结果的客观性与精确。
本发明方法包括如下步骤:
步骤1,构建扇区动态交通特征度量指标集合;
步骤2,构建动态交通样本数据集;
步骤3,基于自适应高斯混合聚类算法进行交通特征分类;
步骤4,基于熵权法进行扇区动态交通特征指标赋权;
步骤5,基于随机森林进行管制员工作负荷评估。
步骤1包括:构建扇区动态交通特征度量指标集合W={L1,L2,L3},其中L1为管制员思考负荷相关联的交通特征指标集合, L1={Lflightlevel,LAreaspecial,LNumspecial,Lmix,Lroute,Lborder},L2为管制员通话负荷相关联的交通特征指标集合,L2={LACtime,LACdis,Lc,Ld,Lh},L3为管制员操作负荷相关联的交通特征指标集合,L3={LACnum,LACover};
其中各指标说明如下:
Lflightlevel:扇区所用高度层个数,指统计时段内,航空器占用飞行高度层的个数;
LAreaspecial:特殊空域面积,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区域面积之和;
LNumspecial:特殊空域个数,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区的个数;
Lmix:机队混合指数,指统计时段内,扇区内运行航空器机型种类数量;
Lroute:扇区路径数量,指统计时段内,扇区进离扇区飞行路径数量之和;
Lborder:扇区边界接近数量,指统计时段内,扇区内飞行的航空器距离扇区边界距离小于Z1海里的事件次数;
LACtime:航空器飞行时间,指统计时段内,航空器在扇区内飞行时间的平均值,起飞航空器的飞行时间从起飞开始至离开扇区,着陆航空器的飞行时间从进入扇区开始至降落;
LACdis :航空器扇区内平均飞行距离,指统计时段内,航空器在扇区飞行的水平距离的平均值;
Lc:爬升航空器数量,指统计时段内,在扇区飞行过程中航空器出现爬升动作的次数之和;
Ld:下降航空器数量,指统计时段内,在扇区飞行过程中航空器出现下降动作的次数之和;
Lh:巡航航空器数量,指统计时段内,在扇区飞行过程中一直保持平飞动作航空器数量之和;
LACnum:航空器数量,指统计时段内,进离扇区的航空器数量,所述进离扇区的航空器包括起降航空器和飞越航空器;
LACover:飞越航空器数量,指统计时段内,飞越扇区的航空器数量。
步骤2包括:根据采集的交通数据,以等长时间片为单位进行统计计算,形成动态交通样本数据集T:
其中N为样本总量,M为交通特征度量指标数量,xNM为第N个样本的第M个指标的样本值,XN为第N个样本的指标横向量,YM为第M个指标的样本竖向量。
步骤3包括:步骤3.1,对动态交通样本数据集T进行MAX-MIN(最大值-最小值) 标准化处理,对于全体M个交通度量指标,针对第m个指标,第n个样本,标准化处理后的指标为:
其中,max(Ym)是动态交通样本数据集T中第m个指标的最大值,min(Ym)是动态交通样本数据集T第m个指标的最小值,并设置聚类类别个数为Kinitial=2,最大聚类类别个数K;
步骤3.2,如果Kinitial≤K,则表明尚未完成各个分类数的聚类,执行步骤3.3;如果Kinitial>K,则表明已完成各个分类数的聚类,执行步骤3.7;
步骤3.3,以Kinitial为分类类别数,并以Kinitial为初始参数构建高斯混合模型:
式中:P(x)为满足高斯混合模型的概率函数,x为动态交通样本数据集T中的随机向量,λi为第i个高斯分布的混合系数,并且满足:λi>0,并且p(x|ui,Σi)为服从均值ui和协方差Σi的第i个高斯分布的概率密度函数:
步骤3.4,对于动态交通样本数据集T中第n个样本xn,判断所属分类,n≤N;
步骤3.5,计算聚类类别数为Kinitial时聚类结果的轮廓系数其中s(n)为第n个样本的轮廓系数,且a(n)为样本xn到同簇其他样本xm的平均距离,b(n)为平均不相似度,指样本xn到其他簇中样本xm'的平均距离,其中,|Cq| 为样本xn的同簇样本数量,d(n,m)为样本xn到样本xm之间的距离,|Cp|为第p簇的样本数量,第p簇与样本xn非同簇,d(n,m')为样本xn到样本xm'之间的距离;将轮廓系数S(Kinitial)添加至集合R中;
步骤3.6,迭代更新分类数量Kinitial=Kinitial+1,执行步骤3.3;
步骤3.7,选择集合R中的最大值,其对应的Kinitial为最佳聚类类别数Koptimal;以最佳聚类数Koptimal对应的聚类结果作为最终聚类结果,完成交通特征分类。
步骤3.4包括:
步骤3.4.1,根据贝叶斯定理,求出初始参数下样本xn由第j个高斯成分组成的后验概率γnj,其中:
步骤3.4.2,根据后验概率γnj和初始设定的分类类别,利用最大似然估计更新{λi,ui,Σi}后,再次计算后验概率γij;
步骤3.4.3,判断结果是否收敛,如果收敛则执行步骤3.4.4,否则执行步骤3.4.2;
步骤3.4.4,选择后验概率γij最大的类别作为样本xn所属类别。
步骤4包括:
步骤4.1,设置聚类类别为k,初始化k=1,对于第k个类别,存在M个评价指标,有Nk个样本数据,则经过归一化后得到非负矩阵:
步骤4.2,根据非负矩阵计算概率矩阵,对于非负矩阵的第g行第m列的元素xgm,其概率pgm为:
步骤4.3,计算第k类别下,第m个指标的信息熵ekm和信息效用值dkm,最终计算权重ωkm,其中:
dkm=1-ekm,
步骤4.4,计算第k类别各指标权重后,继续计算下一类别指标权重,令k=k+1,判断k<=Kinitial是否成立,如果不成立则表明所有类别下的指标权重计算完毕,执行步骤4.5,否则继续计算第k类别下各指标权重;
步骤4.5,得到指标权重矩阵V,并根据指标权重矩阵计算管制员工作负荷Lnk:
步骤4.6,得到各个样本的管制员工作负荷指标后,选择最大值和最小值,利用步骤3.1所描述的最大值-最小值标准化处理方法,将管制员工作负荷指标归一化,归一化后的管制员工作负荷指标最大值为1,最小值为0。
步骤5包括:
步骤5.1,根据聚类结果对样本数据进行标记,得到标记后的动态交通样本数据集T,在动态交通样本数据集T中,随机抽取一部分作为测试集Ttest,剩余的作为训练集 Ttrain,其中训练集样本个数为Ntrain;
步骤5.2,选择决策树最大个数n_estimators和最大特征数m_features作为可调节参数,根据样本个数和指标个数设置参数调优范围为:n_estimators={10,11,...,20},m_features={1,2,...,M},其中n_estimators为基决策树个数,m_features为最大特征数;
步骤5.3,利用网格搜索方法,对n_estimators和m_features进行遍历;
步骤5.4,当m_features=d,n_estimators=r时,其中,d为决策树最大个数参数值,r为最大特征数参数值,有放回抽取Ntrain次样本,获得与Ntrain大小相同的数据集随机选择nt个特征值,重复r次,建立r棵决策树,其中nt<M;
步骤5.5,计算每棵决策树的结果,通过多数投票法,对所有叶子节点的输出结果进行累计处理,取预测概率总和最大的那一类作为分类结果;
步骤5.6,利用十折交叉验证计算分类准确性Pdr,并存储Pdr;
步骤5.7,判断遍历是否结束,如果没有则按照网格将n_estimators和m_features移至下一网格处,重新调整参数执行步骤5.4,如果结束则执行步骤5.8;
步骤5.8,选择Pdr最大的参数作为最优参数,以此建立基于随机森林的扇区交通运行场景分类模型。
步骤5.9,根据扇区交通运行场景分类模型,确定已采集的指定统计时间段内的动态交通样本数据的类别,并提取该类别对应的指标权重,利用指标值和指标权重加权求和的方法计算统计时间段内的管制员工作负荷指数。
利用利用步骤3.1所描述的最大值-最小值标准化处理方法,将得到的管制员工作负荷指数归一化,归一化后的数值越靠近1,则表明该时间段内的管制员工作负荷强度越大,反之,归一化之后的数值越接近0.则表明该时间段内的管制员工作负荷强度越小,从而实现管制员工作负荷的评估。
有益效果:本发明提出了一种基于扇区动态交通特征的管制员工作负荷评估方法,通过提取交通运行特征、交通域向负荷域的特征演进、决策模型构建实现了客观的管制员工作负荷评估,为动态、实时、精细的管制员工作负荷评估提供了理论支撑。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/ 或其他方面的优点将会变得更加清楚。
图1是本发明方法流程图。
具体实施方式
本发明提出一种基于扇区动态交通特征的管制员工作负荷评估方法,能够依据扇区交通运行特征,客观度量扇区管制员工作负荷。
本发明针对扇区结构和航班流分布,构建扇区交通运行特征度量指标集合以及样本集,利用高斯混合聚类算法对样本进行聚类,实现交通运行特征提取;利用熵权法对实现交通特征度量向管制员工作负荷度量的演进;利用决策森林模型对待评估样本进行分类,从而计算对应的管制员工作负荷。
技术方案:本发明是这样实现的,一种基于扇区动态交通特征的管制员工作负荷评估方法,所述方法包括以下步骤,流程图如图1所示。
步骤1:构建扇区动态交通特征度量指标集合
管制员工作负荷是在给扇区内航空器提供各类管制服务过程中管制员所承担的心理和生理压力,管制员工作的主要内容包括指挥航空器运行、维持航空器安全间隔、为飞行员提供情报、与其他管制员协调等,因此,国内外相关研究机构通常将管制员工作负荷细分为思考负荷、通讯负荷和操作负荷三大类。
其中,管制员的思考负荷来源于对扇区结构和扇区交通流态势的感知,以及在感知的基础上进行冲突解脱所作出的决策。航空器在扇区运行过程中,管制员会根据扇区内航空器占用的高度层个数、特殊空域面积、特殊空域个数、机队混合指数、路径数量、扇区边界接近数量进行大致判断。
管制员的通话负荷来源于保障航空器正常运行的基础通话和进行冲突解脱的附加通话。正常运行情况下,航空器在扇区内的飞行时间和飞行距离决定了管制员的基础通话数量,冲突解脱过程中,飞行员根据指令采取爬升、下降、巡航等操作,因此爬升航空器数量、下降航空器数量、巡航航空器数量等指标与管制员附加通话相关;
管制员的操作负荷来源于管制过程中管制员记录信息和协调移交。扇区内飞行的航空器数量、飞越扇区的航空器数量等因素都会影响管制员的操作负荷。
综合以上分析,本发明基于管制员产生负荷的工作流程,梳理扇区动态交通特征中与管制员思考、通讯和操作负荷相关联的因素,以此构建扇区动态交通特征度量指标集合。包括:扇区所用高度层个数、特殊空域面积、特殊空域个数、扇区路径数量、航空器数量、航空器飞行时间、机队混合指数、航空器扇区内平均飞行距离、爬升航空器数量、下降航空器数量、巡航航空器数量、飞越航空器数量、扇区边界接近数量。
扇区动态交通特征指标集合表示为W={L1,L2,L3},其中L1为管制员思考负荷相关联的交通特征指标集合,L1={Lflightlevel,LAreaspecial,LNumspecial,Lmix,Lroute,Lborder},L2为管制员通话负荷相关联的交通特征指标集合,L2={LACtime,LACdis,Lc,Ld,Lh},L3为管制员操作负荷相关联的交通特征指标集合,L3={LACnum,LACover}
其中各指标说明如下:
Lflightlevel:扇区所用高度层个数,指统计时段内,航空器占用飞行高度层的个数;
LAreaspecial:特殊空域面积,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区域面积之和;
LNumspecial:特殊空域个数,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区的个数;
Lmix:机队混合指数,指统计时段内,扇区内运行航空器机型种类数量;
Lroute:扇区路径数量,指统计时段内,扇区进离扇区飞行路径数量之和;
Lborder:扇区边界接近数量,指统计时段内,扇区内飞行的航空器距离扇区边界距离小于Z1(一般取值为10)海里的事件次数;
LACtime:航空器飞行时间,指统计时段内,航空器在扇区内飞行时间的平均值,起飞航空器的飞行时间从起飞开始至离开扇区,着陆航空器的飞行时间从进入扇区开始至降落;
LACdis :航空器扇区内平均飞行距离,指统计时段内,航空器在扇区飞行的水平距离的平均值;
Lc:爬升航空器数量,指统计时段内,在扇区飞行过程中航空器出现爬升动作的次数之和;
Ld:下降航空器数量,指统计时段内,在扇区飞行过程中航空器出现下降动作的次数之和;
Lh:巡航航空器数量,指统计时段内,在扇区飞行过程中一直保持平飞动作航空器数量之和;
LACnum:航空器数量,指统计时段内,进离扇区的航空器数量,所述进离扇区的航空器包括起降航空器和飞越航空器;
LACover:飞越航空器数量,指统计时段内,飞越扇区的航空器数量;
步骤2:构建交通样本集
根据国际民航组织9971文件中描述,扇区管制员工作负荷通常以15分钟或60分钟为单位进行度量。因此,根据采集的交通数据,以等长时间片为单位进行统计计算,形成动态交通动态交通样本数据集T:
其中N为样本总量,M为交通特征度量指标数量,xNM为第N个样本的第M个指标的样本值,XN为第N个样本的指标横向量,YM为第M个指标的样本竖向量。
步骤3:基于自适应高斯混合聚类算法进行交通特征分类
影响管制员工作负荷的因素包含步骤1提出的各个交通度量指标,但不同的交通流时空分布会体现出不同的管制员工作负荷特征,因此需首先将复杂、多样的交通流态势提取为典型态势场景,实现交通特征的分类。本发明采用自适应的高斯混合聚类算法进行交通特征分类。
自适应的高斯混合聚类算法步骤如下所示:
步骤3.1,对动态交通样本数据集T进行MAX-MIN(最大值-最小值)标准化处理,对于全体M个交通度量指标,针对第m个指标,第n个样本,标准化处理后的指标为:其中,max(Ym)是动态交通样本数据集T中第m个指标的最大值,min(Ym)是动态交通样本数据集T第m个指标的最小值,并设置聚类类别个数为 Kinitial=2,最大聚类类别个数K;
步骤3.2,如果Kinitial≤K,则表明尚未完成各个分类数的聚类,执行步骤3.3;如果Kinitial>K,则表明已完成各个分类数的聚类,执行步骤3.7;
步骤3.3,以Kinitial为分类类别数,并以Kinitial为初始参数构建高斯混合模型:式中:x为动态交通样本数据集T中的随机向量,λi为混合系数,并且满足:λi>0,并且p(x|ui,Σi)为服从均值ui和协方差Σi的第i个高斯分布的概率密度函数
步骤3.4,对于动态交通样本数据集T中第n个样本xn,判断所属分类,n≤N,具体方法为:
步骤3.4.1,根据贝叶斯定理,求出初始参数下样本xn由第j个高斯成分组成的后验概率γnj,其中:
步骤3.4.2,根据后验概率γnj和初始设定的分类类别利用最大似然估计更新{λi,ui,Σi}后,再次计算后验概率γij;
步骤3.4.3,判断结果是否收敛,如果收敛则执行步骤3.4.4,否则执行步骤3.4.2;
步骤3.4.4,选择后验概率γij最大的类别作为样本xn所属类别;
步骤3.5,计算聚类类别数为Kinitial时聚类结果的轮廓系数其中a(n)为样本xn到同簇其他样本xm的平均距离,b(n)为平均不相似度,指样本xn到其他簇中样本xm'的平均距离,其中,|Cq|为样本xn的同簇样本数量,d(n,m)为样本xn到样本xm之间的距离,|Cp|为样本xn的同簇样本数量,d(n,m')为样本xn到样本xm'之间的距离;将轮廓系数S(Kinitial)添加至集合R中;
步骤3.6,迭代更新分类数量Kinitial=Kinitial+1,执行步骤3.3;
步骤3.7,选择集合R中的最大值,其对应的Kinitial为最佳聚类类别数Koptimal;以最佳聚类数Koptimal对应的聚类结果作为最终聚类结果。
步骤4:基于熵权法进行扇区动态交通特征指标赋权
在完成交通特征的分类后,需要对各类交通特征所能表达的管制员工作负荷进行提取,本发明采用熵权法对扇区动态交通特征进行负荷权重赋权。
熵权法能够根据信息熵的原理,利用指标的熵值来判断指标的离散程度,熵值小,指标的离散程度大,对管制员工作负荷的影响也会越大,其具体步骤如下:
步骤4.1,设置聚类类别为k,初始化k=1,对于第k个类别,存在M个评价指标,有Nk个样本数据,则经过归一化后得到非负矩阵:
步骤4.2,根据非负矩阵计算概率矩阵,对于非负矩阵的第g行第m列的元素xgm,其概率pgm为:
步骤4.3,计算第k类别下,第m个指标的信息熵ekm和信息效用值dkm,最终计算权重ωkm,其中:
dkm=1-ekm
根据熵权法计算各个交通流态势下各个管制员工作负荷测量指标的权重后,便可计算出每个类别下的管制员工作负荷值;
步骤4.4,计算第k类别各指标权重后,继续计算下一类别指标权重,令k=k+1,判断k<=Kinitial是否成立,如果不成立则表明所有类别下的指标权重计算完毕,执行步骤4.5,否则继续计算第k类别下各指标权重;
步骤4.5,得到指标权重矩阵V,并根据指标权重矩阵计算管制员工作负荷指标Lnk:
步骤4.6,得到各个样本的管制员工作负荷指标后,选择最大值和最小值,利用步骤3.1所描述的最大值-最小值标准化处理方法,将管制员工作负荷指标归一化,归一化后的管制员工作负荷指标最大值为1,最小值为0。
步骤5:基于随机森林分类模型进行管制员工作负荷评估
随机森林是一种能够利用多棵决策树对样本进行训练并进行预测的分类模型,通过投票等方式集成多棵决策树的结果最终得到更为准确分类。随机森林具有准确率高、适用于大数据集、能够处理高维特征等优点。
在高斯混合聚类算法的结果基础上,利用随机森林分类模型,对指定扇区的交通特征数据进行训练,得出给予随机森林的交通特征分类决策模型,当面对特定的待决策对象时,可以根据随机森林模型将动态交通样本进行分类,进而通过对应的管制负荷权重计算管制员工作负荷。算法具体步骤为:
步骤5.1,根据聚类结果对样本数据进行标记,得到标记后的动态交通样本数据集T,在动态交通样本数据集T中,随机抽取20%作为测试集Ttest,80%作为训练集Ttrain,其中训练集样本个数为Ntrain;
步骤5.2,选择决策树最大个数n_estimators和最大特征数m_features作为可调节参数,根据样本个数和指标个数设置参数调优范围为:n_estimators={10,11,...,20},m_features={1,2,...,M},其中n_estimators为基决策树个数,m_features为最大特征数;
步骤5.3,利用网格搜索方法,对n_estimators和m_features进行遍历;
步骤5.4,当m_features=d,n_estimators=r时,其中,d为决策树最大个数参数值,r为最大特征数参数值,有放回抽取Ntrain次样本,获得与Ntrain大小相同的数据集随机选择nt个特征值,重复r次,建立r棵决策树。
步骤5.5,计算每棵决策树的结果,通过多数投票法,对所有叶子节点的输出结果进行累计处理,取预测概率总和最大的那一类作为分类结果;
步骤5.6,利用十折交叉验证计算分类准确性Pdr,并存储Pdr;
步骤5.7,判断遍历是否结束,如果没有则按照网格将参数移至下一网格处,重新调整参数执行步骤5.4,如果结束则执行步骤5.8;
步骤5.8,选择Pdr最大的参数作为最优参数,以此建立基于随机森林的扇区交通运行场景分类模型。
步骤5.9,根据扇区交通运行场景分类模型,确定已采集的指定统计时间段内的动态交通样本数据的类别,并提取该类别对应的指标权重,利用指标值和指标权重加权求和的方法计算统计时间段内的管制员工作负荷指数。
利用利用步骤3.1所描述的最大值-最小值标准化处理方法,将得到的管制员工作负荷指数归一化,归一化后的数值越靠近1,则表明该时间段内的管制员工作负荷强度越大,反之,归一化之后的数值越接近0.则表明该时间段内的管制员工作负荷强度越小,从而实现管制员工作负荷的评估。
本发明提供了一种基于扇区动态交通特征的管制员工作负荷评估方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (1)
1.一种基于扇区动态交通特征的管制员工作负荷评估方法,其特征在于,包括如下步骤:
步骤1,构建扇区动态交通特征度量指标集合;
步骤2,构建动态交通样本数据集;
步骤3,基于自适应高斯混合聚类算法进行交通特征分类;
步骤4,基于熵权法进行扇区动态交通特征指标赋权;
步骤5,基于随机森林进行管制员工作负荷评估;
步骤1包括:构建扇区动态交通特征度量指标集合W={L1,L2,L3},其中L1为管制员思考负荷相关联的交通特征指标集合,L1={Lflightlevel,LAreaspecial,LNumspecial,Lmix,Lroute,Lborder},L2为管制员通话负荷相关联的交通特征指标集合,L2={LACtime,LACdis,Lc,Ld,Lh},L3为管制员操作负荷相关联的交通特征指标集合,L3={LACnum,LACover};
其中各指标说明如下:
Lflightlevel:扇区所用高度层个数,指统计时段内,航空器占用飞行高度层的个数;
LAreaspecial:特殊空域面积,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区域面积之和;
LNumspecial:特殊空域个数,指统计时段内,不能正常使用的危险区、禁止区、限制区和警告区的个数;
Lmix:机队混合指数,指统计时段内,扇区内运行航空器机型种类数量;
Lroute:扇区路径数量,指统计时段内,扇区进离扇区飞行路径数量之和;
Lborder:扇区边界接近数量,指统计时段内,扇区内飞行的航空器距离扇区边界距离小于Z1海里的事件次数;
LACtime:航空器飞行时间,指统计时段内,航空器在扇区内飞行时间的平均值,起飞航空器的飞行时间从起飞开始至离开扇区,着陆航空器的飞行时间从进入扇区开始至降落;
LACdis :航空器扇区内平均飞行距离,指统计时段内,航空器在扇区飞行的水平距离的平均值;
Lc:爬升航空器数量,指统计时段内,在扇区飞行过程中航空器出现爬升动作的次数之和;
Ld:下降航空器数量,指统计时段内,在扇区飞行过程中航空器出现下降动作的次数之和;
Lh:巡航航空器数量,指统计时段内,在扇区飞行过程中一直保持平飞动作航空器数量之和;
LACnum:航空器数量,指统计时段内,进离扇区的航空器数量,所述进离扇区的航空器包括起降航空器和飞越航空器;
LACover:飞越航空器数量,指统计时段内,飞越扇区的航空器数量;
步骤2包括:根据采集的交通数据,以等长时间片为单位进行统计计算,形成动态交通样本数据集T:
其中N为样本总量,M为交通特征度量指标数量,xNM为第N个样本的第M个指标的样本值,XN为第N个样本的指标横向量,YM为第M个指标的样本竖向量;
步骤3包括:步骤3.1,对动态交通样本数据集T进行MAX-MIN(最大值-最小值)标准化处理,对于全体M个交通度量指标,针对第m个指标,第n个样本,标准化处理后的指标为:
其中,max(Ym)是动态交通样本数据集T中第m个指标的最大值,min(Ym)是动态交通样本数据集T第m个指标的最小值,并设置聚类类别个数为Kinitial=2,最大聚类类别个数K;
步骤3.2,如果Kinitial≤K,则表明尚未完成各个分类数的聚类,执行步骤3.3;如果Kinitial>K,则表明已完成各个分类数的聚类,执行步骤3.7;
步骤3.3,以Kinitial为分类类别数,并以Kinitial为初始参数构建高斯混合模型:
式中:P(x)为满足高斯混合模型的概率函数,x为动态交通样本数据集T中的随机向量,λi为第i个高斯分布的混合系数,并且满足:并且p(x|ui,Σi)为服从均值ui和协方差Σi的第i个高斯分布的概率密度函数:
步骤3.4,对于动态交通样本数据集T中第n个样本xn,判断所属分类,n≤N;
步骤3.5,计算聚类类别数为Kinitial时聚类结果的轮廓系数其中s(n)为第n个样本的轮廓系数,且a(n)为样本xn到同簇其他样本xm的平均距离,b(n)为平均不相似度,指样本xn到其他簇中样本xm'的平均距离,其中,|Cq|为样本xn的同簇样本数量,d(n,m)为样本xn到样本xm之间的距离,|Cp|为第p簇的样本数量,第p簇与样本xn非同簇,d(n,m')为样本xn到样本xm'之间的距离;将轮廓系数S(Kinitial)添加至集合R中;
步骤3.6,迭代更新分类数量Kinitial=Kinitial+1,执行步骤3.3;
步骤3.7,选择集合R中的最大值,其对应的Kinitial为最佳聚类类别数Koptimal;以最佳聚类数Koptimal对应的聚类结果作为最终聚类结果,完成交通特征分类;
步骤3.4包括:
步骤3.4.1,根据贝叶斯定理,求出初始参数下样本xn由第j个高斯成分组成的后验概率γnj,其中:
步骤3.4.2,根据后验概率γnj和初始设定的分类类别,利用最大似然估计更新{λi,ui,Σi}后,再次计算后验概率γij;
步骤3.4.3,判断结果是否收敛,如果收敛则执行步骤3.4.4,否则执行步骤3.4.2;
步骤3.4.4,选择后验概率γij最大的类别作为样本xn所属类别;
步骤4包括:
步骤4.1,设置聚类类别为k,初始化k=1,对于第k个类别,存在M个评价指标,有Nk个样本数据,则经过归一化后得到非负矩阵:
步骤4.2,根据非负矩阵计算概率矩阵,对于非负矩阵的第g行第m列的元素xgm,其概率pgm为:
步骤4.3,计算第k类别下,第m个指标的信息熵ekm和信息效用值dkm,最终计算权重ωkm,其中:
dkm=1-ekm,
步骤4.4,计算第k类别各指标权重后,继续计算下一类别指标权重,令k=k+1,判断k<=Kinitial是否成立,如果不成立则表明所有类别下的指标权重计算完毕,执行步骤4.5,否则继续计算第k类别下各指标权重;
步骤4.5,得到指标权重矩阵V,并根据指标权重矩阵计算管制员工作负荷Lnk:
步骤5包括:
步骤5.1,根据聚类结果对样本数据进行标记,得到标记后的动态交通样本数据集T,在动态交通样本数据集T中,随机抽取一部分作为测试集Ttest,剩余的作为训练集Ttrain,其中训练集样本个数为Ntrain;
步骤5.2,选择决策树最大个数n_estimators和最大特征数m_features作为可调节参数,根据样本个数和指标个数设置参数调优范围为:n_estimators={10,11,...,20},m_features={1,2,...,M},其中n_estimators为基决策树个数,m_features为最大特征数;
步骤5.3,利用网格搜索方法,对n_estimators和m_features进行遍历;
步骤5.4,当m_features=d,n_estimators=r时,其中,d为决策树最大个数参数值,r为最大特征数参数值,有放回抽取Ntrain次样本,获得与Ntrain大小相同的数据集随机选择nt个特征值,重复r次,建立r棵决策树,其中nt<M;
步骤5.5,计算每棵决策树的结果,通过多数投票法,对所有叶子节点的输出结果进行累计处理,取预测概率总和最大的那一类作为分类结果;
步骤5.6,利用十折交叉验证计算分类准确性Pdr,并存储Pdr;
步骤5.7,判断遍历是否结束,如果没有则按照网格将n_estimators和m_features移至下一网格处,重新调整参数执行步骤5.4,如果结束则执行步骤5.8;
步骤5.8,选择Pdr最大的参数作为最优参数,以此建立基于随机森林的扇区交通运行场景分类模型;
步骤5.9,根据扇区交通运行场景分类模型,确定已采集的指定统计时间段内的动态交通样本数据的类别,并提取该类别对应的指标权重,利用指标值和指标权重加权求和的方法计算统计时间段内的管制员工作负荷指数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110812169.7A CN113362604B (zh) | 2021-07-19 | 2021-07-19 | 一种基于扇区动态交通特征的管制员工作负荷评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110812169.7A CN113362604B (zh) | 2021-07-19 | 2021-07-19 | 一种基于扇区动态交通特征的管制员工作负荷评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113362604A CN113362604A (zh) | 2021-09-07 |
CN113362604B true CN113362604B (zh) | 2022-05-17 |
Family
ID=77539768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110812169.7A Active CN113362604B (zh) | 2021-07-19 | 2021-07-19 | 一种基于扇区动态交通特征的管制员工作负荷评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113362604B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113838311B (zh) * | 2021-09-16 | 2023-09-05 | 民航数据通信有限责任公司 | 一种管制扇区的最优性能航班量分析方法及装置 |
CN113822564A (zh) * | 2021-09-16 | 2021-12-21 | 民航数据通信有限责任公司 | 一种空域仿真分析的飞行计划最小样本量确认方法及装置 |
CN114202921B (zh) * | 2021-12-09 | 2023-08-04 | 洛阳师范学院 | 交通流量的多尺度符号动态熵分析方法 |
CN116523397B (zh) * | 2023-04-25 | 2024-03-08 | 长安大学 | 基于熵权法和gmm聚类算法的城市交通网络弹性评估方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103226900B (zh) * | 2013-03-21 | 2015-10-28 | 北京工业大学 | 一种基于加权图模型的空域扇区划分方法 |
CN104835354B (zh) * | 2015-05-20 | 2017-03-01 | 青岛民航空管实业发展有限公司 | 一种管制负荷管理系统及其管制员工作负荷的评估方法 |
CN105160201A (zh) * | 2015-09-30 | 2015-12-16 | 成都民航空管科技发展有限公司 | 一种基于gabp神经网络的管制员工作负荷的预测方法及系统 |
CN108090613A (zh) * | 2017-12-18 | 2018-05-29 | 南京航空航天大学 | 一种进近管制扇区管制员工作负荷预测方法 |
CN110060513A (zh) * | 2019-01-24 | 2019-07-26 | 中国民用航空飞行学院 | 基于历史轨迹数据的空中交通管制员工作负荷评估方法 |
-
2021
- 2021-07-19 CN CN202110812169.7A patent/CN113362604B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113362604A (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113362604B (zh) | 一种基于扇区动态交通特征的管制员工作负荷评估方法 | |
CN106023655B (zh) | 一种扇区空中交通拥堵状态监测方法 | |
CN103489005B (zh) | 一种基于多分类器融合的高分辨率遥感影像分类方法 | |
CN105956790B (zh) | 低空飞行态势安全性评估指标及其评估方法 | |
CN110008442A (zh) | 一种基于飞行员综合能力评价的飞行员训练控制方法 | |
Kuravsky et al. | Mathematical foundations of flight crew diagnostics based on videooculography data | |
CN104239712B (zh) | 一种雷达抗干扰效能实时评估方法 | |
CN111626366B (zh) | 基于运行特征的区域扇区场景相似识别方法 | |
CN106251065A (zh) | 一种基于眼动行为指标体系的管制效率评估方法 | |
CN109062245B (zh) | 一种无人机地面站系统软件的可靠性智能分配方法 | |
CN113344423B (zh) | 一种基于机器学习的飞行员场景适用性诊断方法及系统 | |
CN115828140A (zh) | 邻域互信息与随机森林相融合故障检测方法、系统及应用 | |
CN110082738B (zh) | 基于高斯混合和张量循环神经网络的雷达目标识别方法 | |
CN111340388A (zh) | 一种基于航班qar数据的飞行员操作品质评价方法 | |
CN115064009A (zh) | 一种终端区无人机与有人机冲突风险等级划分方法 | |
CN115544714A (zh) | 一种基于飞行器编队的时序动态对抗威胁评估方法 | |
CN114398769B (zh) | 一种无人直升机飞行控制系统评分自动获取方法 | |
CN108986554B (zh) | 一种基于模糊综合评判的空域扇区拥挤程度动态识别方法 | |
CN113222229B (zh) | 一种基于机器学习的非合作无人机轨迹预测方法 | |
CN115345408A (zh) | 一种基于熵权法的终端区进场效率评估方法 | |
CN115188223A (zh) | 一种基于机场天气实况的飞机起降安全指数确定方法 | |
Arts et al. | Trajectory Based Flight Phase Identification with Machine Learning for Digital Twins | |
Dong et al. | Evaluation for Trainee Pilot Workload Management Competency During Approach Phase Based on Flight Training Data | |
CN112488146A (zh) | 基于深度学习的飞机危险识别与趋势预测方法 | |
CN118313733B (zh) | 一种飞行动作质量评估方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |