一种航班流量时隙全局优化调度方法与流程

本发明属于交通运输规划与管理领域，具体涉及民用航空空中交通管理(atm)中的一种航班流量时隙全局优化调度软件方法。

背景技术：

民航航班正常工作是世界性难题，也是中国民航当前面临的重要问题。随着国民经济快速发展和交通强国建设，民用航空运输领域日渐繁忙，国家空域内航班流量快速增长，空中交通管制员作为空中交通管理的核心主体，为应对流量增长，各部门都采取相应的技术手段增强各自服务保障能力，以提升服务质量和数量，辅助管制员的决策支持系统应运而生，流量管理系统/工具就其中一类辅助航班运行调度的工具。

目前，中国民航流量管理工具是7个地区级空中交通管理局(华北局、东北局、华东局、中南局、西南局、西北局、新疆局)独立建设的工具，并未在全国层面统筹考虑设计航班排序调度算法，流量控制手段割裂低效，无形中存在可以利用时隙的浪费现象，全国航班运行的总体延误得不到缓解；存在的问题缺陷主要如下：

(1)流量管理措施层层转发，层层加码，时隙浪费。

(2)航空公司二次惩罚问题：航班由于机械故障或其他公司自身原因向民航空管即使提前通报延误信息，会自身延误基础上再延误。

(3)时隙资源不属航空公司：雷雨等不利条件下，对于航空公司积极响应民航空管提前取消航班，但这部分航班不再分配时隙，时隙资源本属于航空公司，公司可拿来交换给其他该公司的航班。

技术实现要素：

针对于上述现有技术的不足，本发明的目的在于提供一种航班流量时隙全局优化调度方法，以解决现有技术中流控措施割裂层层转发加码导致时隙浪费和航空公司二次惩罚等问题。

为达到上述目的，本发明采用的技术方案如下：

本发明的一种航班流量时隙全局优化调度方法，步骤如下：

1)确定受流控影响的航班集合范围；

2)人工定义豁免航班；

3)将非豁免航班分类处理；

4)计算每个航班的最早计算落地时间的时间阀值；

5)创建航班排序队列q；

6)根据用户输入的机场进场接收率分配虚拟的时隙资源；

7)对每个航班分配计算落地时间。

进一步地，所述步骤1)具体包括：确定包含在航班优化排序中的航班计划集合i，i＝a∪b

a＝{f∈f:t0≤eldtf≤t1}

b＝{f∈f:eldtf＞t1and(t0≤cldtf≤t1ort0≤sldtf≤t1}

式中，航班集合a表示航班的预计落地时间落入流控时间范围的航班，航班集合b表示预计落地时间晚于t1，且计划落地时间(sldt)落入流控范围或已有计算落地时间(cldt)在流控时间范围；f表示航班集合，f表示每架航班，eldtf表示航班f预计落地时间，cldtf表示航班f的计算落地时间，sldtf表示航班f的时刻表落地时间，t0为流控开始时间，t1为流控结束时间。

进一步地，所述步骤2)具体包括：由流量人员定义航班集合e中的航班豁免，可直接指定某些航班豁免，或通过定义经过某些航路/航路点或某些机场起飞的航班集合豁免：

进一步地，所述步骤3)具体包括：已经被分配过计算落地时间(至少一次)的航班集合；非豁免但包含航班中除去已分配过时隙的航班集合；

f1＝{f∈i-e:fhadacldt}

f2＝(i-e)-f1

i＝e∪f1∪f2

式中，f1表示受影响航班集合i中，除去步骤2中豁免航班和已拥有了计算落地时间时隙的航班，f2表示受影响航班集合i中，除去豁免航班和f1航班集合。

进一步地，所述步骤4)具体包括：

41)若航班是人工操作的豁免航班(如vip或应急管理航班)，则对于已分配时隙资源的航班，最早cldt则取预计落地时间；对于未分配时隙资源的航班，则最早cldt取预计落地时间或计划落地时间中的较小值；

对于已错过其预计撤轮档时间(eobt)的航班，考虑最早cldt的阀值，若当前时间>eobt,则最早cldt＝当前时间+vsp+空中飞行时间；

式中，cldt为计算落地时间，slot表示航班所占用的时隙资源，eobt表示预计离港时间，slotid表示时隙资源分配后的编号，sibt表示计划到港时间，vit表示航班在停机位和跑道头之间的可变滑行时间，vsp为系统可变参数；

42)若航班属于集合f1和f2，则：

当前时间current_time小于或等于sobt，则最早cldt＝sidt-vtt-vsp1；

当前时间current_time大于sobt，则最早cldt＝当前时间+空中飞行时间+vsp2

式中，sobt表示航班计划离港时间，vsp1、vsp2为用户定义适应性参数，vtt表示可变滑行时间，eet表示空中飞行所需时间。

进一步地，所述步骤5)具体包括：

51)根据航班的eldt作为优先级，创建豁免航班的排序子队列q(e)；

52)创建已分配时隙的航班队列，根据已拥有的cldt先后，创建航班排序子队列q(f1)；

53)创建剩下航班(即除去豁免航班集合和已有时隙航班集合，剩余受控航班)的子队列q(f2)，根据航班的sldt初步排序；

54)将上述步骤51)、52)、53)中队列，重新排序合成一个带标识属性的排序总队列q；

式中，是算法创建的航班排序子队列，表示一架架航班，来自集合e，使用eldt作为优先级；是算法创建的航班排序子队列，来自集合f1，使用cldt作为优先级；是算法创建的航班排序子队列，来自集合f2，使用sldt作为优先级。

对航班队列q进行重新编号排序变为q1,q2,...,q|i|。

如某繁忙机场的进场接受率aar设置为6架/小时，开始时间＝1801，结束时间＝2159，则时隙资源集合s＝{1810，1820，…,2140,2150}。

进一步地，所述步骤7)具体包括：根据步骤5)中的航班队列q和步骤6)中的时隙资源slots对每个航班分配落地时隙；

计算方法采用双层循环计算，外层循环是流控措施的集合，内层循环是相关航班的待计算队列；

计算遍历队列中的每个航班f∈i，从剩余可用时隙集合s中，计算分配给该航班的时隙；

式中，flight(t)表示时隙属于哪个航班，owner(t)表示时隙属于哪个航空公司，status(t)表示时隙是否被占用。

本发明的有益效果：

1、本发明方法减少了流量管理措施mit从源头往下游层层转发，层层加码和时隙浪费等情况。

2、本发明方法采用时刻表时间(sobt/sldt)作为判断航班先到先服务的基准，避免对航空公司二次惩罚的现象。

3、本发明方法采用了时隙资源属于航空公司规则，提高航空公司协同配合意愿，使其对航空公司具备对该时隙的拥有权，可用于时隙交换等。

附图说明

图1为本发明方法的流程图。

具体实施方式

为了便于本领域技术人员的理解，下面结合实施例与附图对本发明作进一步的说明，实施方式提及的内容并非对本发明的限定。

术语解释：

参照图1所示，本发明的一种航班流量时隙全局优化调度方法，步骤如下：

1)确定受流控影响的航班集合范围；

确定包含在航班优化排序中的航班计划集合i，i＝a∪b

a＝{f∈f:t0≤eldtf≤t1}

b＝{f∈f:eldtf＞t1and(t0≤cldtf≤t1ort0≤sldtf≤t1}

式中，航班集合a表示航班的预计落地时间落入流控时间范围的航班，航班集合b表示预计落地时间晚于t1，且计划落地时间(sldt)落入流控范围或已有计算落地时间(cldt)在流控时间范围；f表示航班集合，f表示每架航班，eldtf表示航班f预计落地时间，cldtf表示航班f的计算落地时间，sldtf表示航班f的时刻表落地时间，t0为流控开始时间，t1为流控结束时间。

2)人工定义豁免航班；

由流量人员定义航班集合e中的航班豁免，可直接指定某些航班豁免，或通过定义经过某些航路/航路点或某些机场起飞的航班集合豁免：

3)将非豁免(但包含在流控影响时间和空间范围内)航班分类处理；

已经被分配过计算落地时间(至少一次)的航班集合；非豁免但包含航班中除去已分配过时隙的航班集合；

f1＝{f∈i-e:fhadacldt}

f2＝(i-e)-f1

i＝e∪f1∪f2

式中，f1表示受影响航班集合i中，除去步骤2中豁免航班和已拥有了计算落地时间时隙的航班，f2表示受影响航班集合i中，除去豁免航班和f1航班集合。

4)计算每个航班的最早计算落地时间的时间阀值；

41)若航班是人工操作的豁免航班(如vip或应急管理航班)，则对于已分配时隙资源的航班，最早cldt则取预计落地时间；对于未分配时隙资源的航班，则最早cldt取预计落地时间或计划落地时间中的较小值；

对于已错过其预计撤轮档时间(eobt)的航班，考虑最早cldt的阀值，如当前时间>eobt,则最早cldt＝当前时间+vsp+空中飞行时间；

式中，cldt为计算落地时间，slot表示航班所占用的时隙资源，eobt表示预计离港时间，slotid表示时隙资源分配后的编号，sibt表示计划到港时间，vit表示航班在停机位和跑道头之间的可变滑行时间，vsp为系统可变参数；

42)若航班属于集合f1和f2，则：

当前时间current_time小于或等于sobt，则最早cldt＝sidt-vtt-vsp1；

当前时间current_time大于sobt，则最早cldt＝当前时间+空中飞行时间+vsp2

式中，sobt表示航班计划离港时间，vsp1、vsp2为用户定义适应性参数，vtt表示可变滑行时间，eet表示空中飞行所需时间。

5)创建航班排序队列q；

51)根据航班的eldt作为优先级，创建豁免航班的排序子队列q(e)；

52)创建已分配时隙的航班队列，根据已拥有的cldt先后，创建航班排序子队列q(f1)；

53)创建剩下航班(即除去豁免航班集合和已有时隙航班集合，剩余受控航班)的子队列q(f2)，根据航班的sldt初步排序；

54)将上述步骤51)、52)、53)中队列，重新排序合成一个带标识属性的排序总队列q；

式中，是算法创建的航班排序子队列，表示一架架航班，来自集合e，使用eldt作为优先级；是算法创建的航班排序子队列，来自集合f1，使用cldt作为优先级；是算法创建的航班排序子队列，来自集合f2，使用sldt作为优先级。

对航班队列q进行重新编号排序变为q1,q2,...,q|i|。

如某繁忙机场的进场接受率aar设置为6架/小时，开始时间＝1801，结束时间＝2159，则时隙资源集合s＝{1810，1820，…,2140,2150}。

6)根据用户输入的机场进场接收率分配虚拟的时隙资源；

7)对每个航班分配计算落地时间；

根据步骤5)中的航班队列q和步骤6)中的时隙资源slots对每个航班分配落地时隙；

计算方法采用双层循环计算，外层循环是流控措施的集合，内层循环是相关航班的待计算队列；

计算遍历队列中的每个航班f∈i，从剩余可用时隙集合s中，计算分配给该航班的时隙；

式中，flight(t)表示时隙属于哪个航班，owner(t)表示时隙属于哪个航空公司，status(t)表示时隙是否被占用。

本发明具体应用途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进，这些改进也应视为本发明的保护范围。