基于航空公司间公平运行的航班计划分配方法及装置制造方法及图纸

本发明专利技术提供了基于航空公司间公平运行的航班计划分配方法及装置，该方法包括：首先，获取待规划的空域集合、航空公司集合、航班集合、以及待分配的航班计划集合，将获取的空域集合、航空公司集合、航班集合、以及航班计划集合输入至预先建立的航班计划管理模型，其中，该航班计划管理模型包括：飞机相遇概率模型和公平运行规则模型；然后，利用航班计划管理模型分别为航空公司集合中的各个航空公司确定相应的航班计划子集；最后，依次输出确定出的多个航班计划子集，并将各个航班计划子集分配至相应的航空公司。本发明专利技术实施例通过引入公平运行规则模型，使得为各个航空公司分配的航班计划更加合理公平。

【技术实现步骤摘要】
基于航空公司间公平运行的航班计划分配方法及装置
本专利技术涉及空中交通流量管理
，具体而言，涉及基于航空公司间公平运行的航班计划分配方法及装置。
技术介绍
航空公司是航空市场运营的主体，它的主营业务是航班。因此能否合理、有效地制定航班计划，对航空公司未来的发展和整体的经济效益有着决定性的影响，它是上承公司战略，下联公司收益的重要的中间环节。航班计划的编制与市场环境、旅客需求、航向、机型、机场条件、飞机性能等因素相关。一个好的航班飞行计划牵涉到中国民航和航空公司内部诸多章法制度，是一个十分复杂繁重的任务集合。国内许多学者对航班计划优化作了大量研究，取得了一定的成果，但真正被航空公司采纳并用到实际却很少见。主要原因是研究者通常采用简化的数学规划对航班计划进行优化，此类模型和算法一般只能对几种机型，小型航班网络的航班计划进行优化，对绝大多数航空公司来说是没有使用价值的。当前，相关技术中提供了一种航班计划的规划方法，该规划方法主要是：根据飞行安全、空中交通管制工作负荷来、飞行时间来对航班计划进行规划，即通过对成对飞机的航班计划进行冲突分析、对负责各个空域的管制人员的工作负荷进行分析、对各个航班的飞行时间进行分析来规划各个航空公司的航班计划。在实现本专利技术的过程中，专利技术人发现相关技术中至少存在以下问题：上述航班计划的规划方法中并未考虑各个航空公司间的公平性运行原则，因而，采用上述航班计划的规划方法得到的航班计划无法保证各个航空公司间公平运行。
技术实现思路
有鉴于此，本专利技术实施例的目的在于提供一种基于航空公司间公平运行的航班计划分配方法及装置，以使得为各个航空公司分配的航班计划更加合理公平。第一方面，本专利技术实施例提供了一种基于航空公司间公平运行的航班计划分配方法，该方法包括：获取待规划的空域集合、航空公司集合、航班集合、以及待分配的航班计划集合，其中，所述空域集合为s∈{1，...，S}，所述航空公司集合为所述航班集合为f∈{1，...，F}，所述航班计划集合为p∈{1，...，Pf}；将获取的所述空域集合、所述航空公司集合、所述航班集合、以及所述航班计划集合输入至预先建立的航班计划管理模型，其中，所述航班计划管理模型包括：飞机相遇概率模型和公平运行规则模型；利用所述航班计划管理模型分别为所述航空公司集合中的各个航空公司确定相应的航班计划子集；依次输出确定出的多个所述航班计划子集，并将各个所述航班计划子集分配至相应的航空公司。结合第一方面，本专利技术实施例提供了第一方面的第一种可能的实施方式，其中，利用所述航班计划管理模型分别为所述航空公司集合中的各个航空公司确定相应的航班计划子集包括：根据获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合、以及多个约束条件方程求解使目标函数最小化的解集；根据求得的所述使目标函数最小化的解集确定所述航空公司集合中的各个航空公司的航班计划子集；其中，所述目标函数公式为所述目标函数公式中的第一项表示为航班f选择航班计划p所需总消耗成本的总和，所述目标函数公式中的第二项表示对各个航空公司的加权平均协作低效率施加的惩罚，所述目标函数公式中的第三项表示对各个航空公司的加权平均协作不公平性施加的惩罚，所述目标函数公式中的第四项表示对各个航空公司的最大加权平均协作不公平性施加的惩罚，所述目标函数公式中的第五项表示对各个空域的工作负荷的总和，所述目标函数公式中的第六项表示对各个空域的工作负荷峰值大于预设工作负荷施加的惩罚，所述目标函数公式中的第七项表示交通管制员消除冲突航班计划(P，Q)冲突发生所需成本，cfp表示执行航班f的航班计划p∈Pf0所需总消耗成本，xfp表示航班f选择航班计划p二进制变量，如果为航班f选择航班计划p，则xfp等于1，否则xfp等于0，μD表示与获得的总加权协作低效率相关的惩罚成本系数，ωa表示每个航空公司a的加权系数，Ea(x)表示航空公司a的协作效率，μe表示与获得的总加权协作不公平性等级相关的惩罚成本系数，xe表示平均协作不公平性，表示与获得的最大加权协作不公平性相关的惩罚成本系数，表示最大加权不公平性，γs表示空域s的监控惩罚成本在每单位平均占有率的工作负荷，ωs表示空域s的平均航班数量，ψsn表示空域s中的监测工作负荷峰值大于预设工作负荷n时的惩罚成本，ysn表示附属于分段线性递增凸惩罚函数的凸组合权重，即空域s中的工作负荷峰值和平均值的差值所引起的惩罚，表示在最大可允许的监视工作负荷下同时占用空域s中的航班数量的最大值，表示消除冲突航班计划(P,Q)冲突发生所需的惩罚成本，zPQ表示选择冲突航班计划(P，Q)的二进制变量，如果选择冲突航班计划(P，Q)，则zPQ等于1。结合第一方面的第一种可能的实施方式，本专利技术实施例提供了第一方面的第二种可能的实施方式，其中，所述根据获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合、以及多个约束条件方程求解使目标函数最小化的解集包括：将获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合代入初始目标函数，得到最终目标函数；通过如下多个约束条件方程对所述最终目标函数中的参数进行约束，计算得到使所述最终目标函数最小化的解集；其中，多个所述约束条件方程包括：基础约束条件方程、第一约束条件方程、第二约束条件方程、第三约束条件方程、第四约束条件方程、第五约束条件方程，所述基础约束条件方程为所述第一约束条件方程为所述第一约束条件方程对应的约束条件为每个航班均对应唯一的航班计划；所述第二约束条件方程为所述第二约束条件方程对应的约束条件为空域s允许的航班数量小于所述空域s的最大航班数量；所述第三约束条件方程为所述第三约束条件方程对应的约束条件为对附属于分段线性递增凸惩罚函数的凸组合权重进行归一化；所述第四约束条件方程为xP+xQ≤1，所述第四约束条件方程对应的约束条件为排除所有存在致命性冲突的成对航班计划；所述第五约束条件方程为所述第五约束条件方程对应的约束条件为在任意时间点空域s内的冲突数量小于允许存在的同时冲突解决的最大冲突数量；其中，ve表示施加于平均协作不公平性xe上的最大值，即加权平均协作不公平性，ns表示空域s的最大航班数量，va表示航空公司a的协作公平性的变量，rs表示空域s中允许存在的同时冲突解决的最大冲突数量，Csi表示空域s占用航班计划的最大交迭集，Msk表示空域s中存在冲突的航班计划组的最大交迭集，Jsk(P)表示冲突曲线图Gsk中与航班计划P相邻的航班计划Q子集，Nsk表示包含入边集Msk的各航班计划，FC表示不可避免的相互冲突的航班计划(P，Q)的子集，A表示可解决的相互冲突的航班计划(P，Q)的子集，xP表示相互冲突的航班计划(P，Q)中航班计划P二进制变量，xQ表示相互冲突的航班计划(P，Q)中航班计划Q二进制变量。结合第一方面的第二种可能的实施方式，本专利技术实施例提供了第一方面的第三种可能的实施方式，其中，上述方法还包括：根据公式计算平均协作低效率，其中，xD表示平均协作低效率，ωa表示每个航空公司a的加权系数，Ea(x)表示航空公司a的协作效率；根据公式计算平均协作不公平性，其中，xe表示平均协作不公平性，ωa表示每个航空公司a的加权系数，va表示航空公司a的本文档来自技高网...

【技术保护点】
基于航空公司间公平运行的航班计划分配方法，其特征在于，所述方法包括：获取待规划的空域集合、航空公司集合、航班集合、以及待分配的航班计划集合，其中，所述空域集合为s∈{1，...，S}，所述航空公司集合为

【技术特征摘要】
1.基于航空公司间公平运行的航班计划分配方法，其特征在于，所述方法包括：获取待规划的空域集合、航空公司集合、航班集合、以及待分配的航班计划集合，其中，所述空域集合为s∈{1，...，S}，所述航空公司集合为所述航班集合为f∈{1，...，F}，所述航班计划集合为p∈{1，...，Pf}；将获取的所述空域集合、所述航空公司集合、所述航班集合、以及所述航班计划集合输入至预先建立的航班计划管理模型，其中，所述航班计划管理模型包括：飞机相遇概率模型和公平运行规则模型；利用所述航班计划管理模型分别为所述航空公司集合中的各个航空公司确定相应的航班计划子集；依次输出确定出的多个所述航班计划子集，并将各个所述航班计划子集分配至相应的航空公司。2.根据权利要求1所述的方法，其特征在于，利用所述航班计划管理模型分别为所述航空公司集合中的各个航空公司确定相应的航班计划子集包括：根据获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合、以及多个约束条件方程求解使目标函数最小化的解集；根据求得的所述使目标函数最小化的解集确定所述航空公司集合中的各个航空公司的航班计划子集；其中，所述目标函数公式为所述目标函数公式中的第一项表示为航班f选择航班计划p所需总消耗成本的总和，所述目标函数公式中的第二项表示对各个航空公司的加权平均协作低效率施加的惩罚，所述目标函数公式中的第三项表示对各个航空公司的加权平均协作不公平性施加的惩罚，所述目标函数公式中的第四项表示对各个航空公司的最大加权平均协作不公平性施加的惩罚，所述目标函数公式中的第五项表示对各个空域的工作负荷的总和，所述目标函数公式中的第六项表示对各个空域的工作负荷峰值大于预设工作负荷施加的惩罚，所述目标函数公式中的第七项表示交通管制员消除冲突航班计划(P，Q)冲突发生所需成本，cfp表示执行航班f的航班计划p∈Pf0所需总消耗成本，xfp表示航班f选择航班计划p二进制变量，如果为航班f选择航班计划p，则xfp等于1，否则xfp等于0，μD表示与获得的总加权协作低效率相关的惩罚成本系数，ωa表示每个航空公司a的加权系数，Ea(x)表示航空公司a的协作效率，μe表示与获得的总加权协作不公平性等级相关的惩罚成本系数，xe表示平均协作不公平性，表示与获得的最大加权协作不公平性相关的惩罚成本系数，表示最大加权不公平性，γs表示空域s的监控惩罚成本在每单位平均占有率的工作负荷，ωs表示空域s的平均航班数量，ψsn表示空域s中的监测工作负荷峰值大于预设工作负荷n时的惩罚成本，ysn表示附属于分段线性递增凸惩罚函数的凸组合权重，即空域s中的工作负荷峰值和平均值的差值所引起的惩罚，表示在最大可允许的监视工作负荷下同时占用空域s中的航班数量的最大值，表示消除冲突航班计划(P,Q)冲突发生所需的惩罚成本，zPQ表示选择冲突航班计划(P，Q)的二进制变量，如果选择冲突航班计划(P，Q)，则zPQ等于1。3.根据权利要求2所述的方法，其特征在于，所述根据获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合、以及多个约束条件方程求解使目标函数最小化的解集包括：将获取的所述空域集合、航空公司集合、航班集合、待分配的航班计划集合代入初始目标函数，得到最终目标函数；通过如下多个约束条件方程对所述最终目标函数中的参数进行约束，计算得到使所述最终目标函数最小化的解集；其中，多个所述约束条件方程包括：基础约束条件方程、第一约束条件方程、第二约束条件方程、第三约束条件方程、第四约束条件方程、第五约束条件方程，所述基础约束条件方程为所述第一约束条件方程为所述第一约束条件方程对应的约束条件为每个航班均对应唯一的航班计划；所述第二约束条件方程为所述第二约束条件方程对应的约束条件为空域s允许的航班数量小于所述空域s的最大航班数量；所述第三约束条件方程为所述第三约束条件方程对应的约束条件为对附属于分段线性递增凸惩罚函数的凸组合权重进行归一化；所述第四约束条件方程为xP+xQ≤1，所述第四约束条件方程对应的约束条件为排除所有存在致命性冲突的成对航班计划；所述第五约束条件方程为所述第五约束条件方程对应的约束条件为在任意时间点空域s内的冲突数量小于允许存在的同时冲突解决的最大冲突数量；其中，ve表示施加于平均协作不公平性xe上的最大值，即加权平均协作不公平性，ns表示空域s的最大航班数量，va表示航空公司a的协作公平性的变量，rs表示空域s中允许存在的同时冲突解决的最大冲突数量，Csi表示空域s占用航班计划的最大交迭集，Msk表示空域s中存在冲突的航班计划组的最大交迭集，Jsk(P)表示冲突曲线图Gsk中与航班计划P相邻的航班计划Q子集，Nsk表示包含入边集Msk的各航班计划，FC表示不可避免的相互冲突的航班计划(P，Q)的子集，A表示可解决的相互冲突的航班计划(P，Q)的子集，xP表示相互冲突的航班计划(P，Q)中航班计划P二进制变量，xQ表示相互冲突的航班计划(P，Q)中航班计划Q二进制变量。4.根据权利要求3所述的方法，其特征在于，所述方法还包括：根据公式计算平均协作低效率，其中，xD表示平均协作低效率，ωa表示每个航空公司a的加权系数，Ea(x)表示航空公司a的协作效率；根据公式计算平均协作不公平性，其中，xe表示平均协作不公平性，ωa表示每个航空公司a的加权系数，va表示航空公司a的协作公平性的变量；根据公式计算对各个航空公司的最大加权平均协作不公平性，其中，表示最大加权不公平性，ωa表示每个航空公司a的加权系数，表示航空公司a的协作公平性；其中，空域s的平均航班数量H表示考虑中的时间范围长度，单位为分钟，Ωs表示在某一时间范围子集内占用空域s的航班计划子集，表示航班f的航班计划p占用空域s的时间长度，单位为分钟，xfp表示航班f选择航班计划p二进制变量，如果为航班f选择航班计划p，则xfp等于1，否则xfp等于0；空域s的最大航班数量与平均航班数量的差值ns表示空域s的最大航班数量，ωs表示空域s的平均航班数量，表示在最大可允许的监视工作负荷下同时占用空域s中的航班数量的最大值，n表示在最大可允许的监视工作负荷下同时占用空域s中的航班数量，ysn表示附属于分段线性递增凸惩罚函数的凸组合权重，即空域s中的工作负荷峰值和平均值的差值所引起的惩罚；航空公司a的协作效率Dmax表示各个航空公司分配的航班计划的成本与各个航空公司单独优化的航班计划的成本之间的最大允许比例，Aa表示属于航空公司a的航班子集，Wf表示航空公司a附加于航班f∈Aa的相对优先权系数，表示每个航班f执行航班计划p所需总消耗成本的最小值，cfp表示执行航班f的航班计划p∈Pf0所需总消耗成本，xfp表示航班f选择航班计划p二进制变量，如果为航班f选择航班计划p，则xfp等于1，否则xfp等于0；航空公司a的协作公平性Ea(x)表示航空公司a的协作效率，ωa表示每个航空公司a的加权系数。5.根据权利要求2-4任一项所述的方法，其特征在于，所述方法还包括：根据公式cfp＝Ffp+Dfp+cf0，f∈{1，...，F}，p∈{1，...，Pf}计算执行航班f的航班计划p所需总消耗成本，其中，Ffp表示燃料消耗成本，Dfp表示航班抵达延误成本，cf0表示航班计划取消成本；通过如下公式计算执行航班f的航班计划p所需总消耗成本中的燃料消耗成本；其中，函数Ffuel-cost(·)是根据欧洲控制实验中心开发和维护的飞机性能数据库确定的；通过如下公式计算执行航班f的航班计划p所需总消耗成本中的航班抵达延误成本；其中，表示航班抵达延误时间，表示航班延误因子，lf表示航班f的乘客载荷估值，δ表示所有航空公司和航空公司各自的航班的每个乘客的平均延误成本；通过如下公式计算执行航班f的航班计划p所需总消耗成本中的航班计划取消成本；其中，Ffp表示执行航班f的航班计划p所需的燃料消耗成本，表示取消航班的预计延误价值，表示航班延误因子，lf表示航班f的乘客载荷估值，δ表示所有航空公司和航空公司各自的航班的每个乘客的平均延误成本。6.基于航空公司间公平运行的航班计划分配装置，其特征在...

【专利技术属性】
技术研发人员：陈志杰，朱永文，唐治理，王长春，付莹，刘晓丽，王家玲，
申请(专利权)人：中国人民解放军空军装备研究院雷达与电子对抗研究所，
类型：发明
国别省市：北京,11