CN113205707B - 一种航班流量时隙全局优化调度方法 - Google Patents

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

Info

Publication number
CN113205707B
CN113205707B CN202110494816.4A CN202110494816A CN113205707B CN 113205707 B CN113205707 B CN 113205707B CN 202110494816 A CN202110494816 A CN 202110494816A CN 113205707 B CN113205707 B CN 113205707B
Authority
CN
China
Prior art keywords
flight
time
flights
queue
cldt
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
Application number
CN202110494816.4A
Other languages
English (en)
Other versions
CN113205707A (zh
Inventor
刘成杰
靳学梅
苏祖辉
胥宝新
程先峰
席玉华
张海芹
张翰文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing LES Information Technology Co. Ltd
Original Assignee
Nanjing LES Information Technology Co. Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing LES Information Technology Co. Ltd filed Critical Nanjing LES Information Technology Co. Ltd
Priority to CN202110494816.4A priority Critical patent/CN113205707B/zh
Publication of CN113205707A publication Critical patent/CN113205707A/zh
Application granted granted Critical
Publication of CN113205707B publication Critical patent/CN113205707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G5/00Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/003Flight plan management
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G5/00Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/003Flight plan management
    • G08G5/0039Modification of a flight plan

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种航班流量时隙全局优化调度方法,步骤如下:确定受流控影响的航班集合范围;人工定义豁免航班;将非豁免航班分类处理;计算每个航班的最早计算落地时间的时间阀值;创建航班排序队列;根据用户输入的机场进场接收率分配虚拟的时隙资源;对每个航班分配计算落地时间。本发明方法采用了时隙资源属于航空公司规则,提高航空公司协同配合意愿,使其对航空公司具备对该时隙的拥有权,可用于时隙交换等。

Description

一种航班流量时隙全局优化调度方法
技术领域
本发明属于交通运输规划与管理领域,具体涉及民用航空空中交通管理(ATM)中的一种航班流量时隙全局优化调度软件方法。
背景技术
民航航班正常工作是世界性难题,也是中国民航当前面临的重要问题。随着国民经济快速发展和交通强国建设,民用航空运输领域日渐繁忙,国家空域内航班流量快速增长,空中交通管制员作为空中交通管理的核心主体,为应对流量增长,各部门都采取相应的技术手段增强各自服务保障能力,以提升服务质量和数量,辅助管制员的决策支持系统应运而生,流量管理系统/工具就其中一类辅助航班运行调度的工具。
目前,中国民航流量管理工具是7个地区级空中交通管理局(华北、东北、华东、中南、西南、西北、新疆)独立建设的工具,并未在全国层面统筹考虑设计航班排序调度算法,流量控制手段割裂低效,无形中存在可以利用时隙的浪费现象,全国航班运行的总体延误得不到缓解;存在的问题缺陷主要如下:
(1)流量管理措施层层转发,层层加码,时隙浪费。
(2)航空公司二次惩罚问题:航班由于机械故障或其他公司自身原因向民航空管即使提前通报延误信息,会自身延误基础上再延误。
(3)时隙资源不属航空公司:雷雨等不利条件下,对于航空公司积极响应民航空管提前取消航班,但这部分航班不再分配时隙,时隙资源本属于航空公司,公司可拿来交换给其他该公司的航班。
发明内容
针对于上述现有技术的不足,本发明的目的在于提供一种航班流量时隙全局优化调度方法,以解决现有技术中流控措施割裂层层转发加码导致时隙浪费和航空公司二次惩罚等问题。
为达到上述目的,本发明采用的技术方案如下:
本发明的一种航班流量时隙全局优化调度方法,步骤如下:
1)确定受流控影响的航班集合范围;
2)人工定义豁免航班;
3)将非豁免航班分类处理;
4)计算每个航班的最早计算落地时间的时间阀值;
5)创建航班排序队列Q;
6)根据用户输入的机场进场接收率分配虚拟的时隙资源;
7)对每个航班分配计算落地时间。
进一步地,所述步骤1)具体包括:确定包含在航班优化排序中的航班计划集合I,
Figure GDA0003579612610000011
I=A∪B
A={f∈F:t0≤ELDTf≤t1}
B={f∈F:ELDTf>t1 and(t0≤CLDTf≤t1 or t0≤SLDTf≤t1)}
式中,航班集合A表示航班的预计落地时间落入流控时间范围的航班,航班集合B表示预计落地时间晚于t1,且计划落地时间(SLDT)落入流控范围或已有计算落地时间(CLDT)在流控时间范围;F表示航班集合,f表示每架航班,ELDTf表示航班f预计落地时间,CLDTf表示航班f的计算落地时间,SLDTf表示航班f的时刻表落地时间,t0为流控开始时间,t1为流控结束时间。
进一步地,所述步骤2)具体包括:由流量人员定义航班集合E中的航班豁免,可直接指定某些航班豁免,或通过定义经过某些航路/航路点或某些机场起飞的航班集合豁免:
Figure GDA0003579612610000021
进一步地,所述步骤3)具体包括:已经被分配过计算落地时间(至少一次)的航班集合;非豁免但包含航班中除去已分配过时隙的航班集合;
F1={f∈I-E:f had a CLDT}
F2=(I-E)-F1
I=E∪F1∪F2
式中,F1表示受影响航班集合I中,除去步骤2中豁免航班和已拥有了计算落地时间时隙的航班,F2表示受影响航班集合I中,除去豁免航班和F1航班集合。
进一步地,所述步骤4)具体包括:
41)若航班是人工操作的豁免航班(如VIP或应急管理航班),则对于已分配时隙资源的航班,最早CLDT则取预计落地时间;对于未分配时隙资源的航班,则最早CLDT取预计落地时间或计划落地时间中的较小值;
对于已错过其预计撤轮档时间(EOBT)的航班,考虑最早CLDT的阀值,若当前时间>EOBT, 则最早CLDT=当前时间+VSP+空中飞行时间;
Figure GDA0003579612610000022
式中,CLDT为计算落地时间,slot表示航班所占用的时隙资源,EOBT表示预计离港时间, slotID表示时隙资源分配后的编号,SIBT表示计划到港时间,VIT表示航班在停机位和跑道头之间的可变滑行时间,VSP为系统可变参数;
42)若航班属于集合F1和F2,则:
当前时间current_time小于或等于SOBT,则最早CLDT=SIB T-VTT-vsp1;
当前时间current_time大于SOBT,则最早CLDT=当前时间+空中飞行时间+vsp2
Figure GDA0003579612610000031
式中,SOBT表示航班计划离港时间,vsp1、vsp2为用户定义适应性参数,VTT表示可变滑行时间,EET表示空中飞行所需时间。
进一步地,所述步骤5)具体包括:
51)根据航班的ELDT作为优先级,创建豁免航班的排序子队列Q(E);
52)创建已分配时隙的航班队列,根据已拥有的CLDT先后,创建航班排序子队列Q(F1);
53)创建剩下航班(即除去豁免航班集合和已有时隙航班集合,剩余受控航班)的子队列Q(F2),根据航班的SLDT初步排序;
54)将上述步骤51)、52)、53)中队列,重新排序合成一个带标识属性的排序总队列Q;
Figure GDA0003579612610000032
Figure GDA0003579612610000033
Figure GDA0003579612610000034
Figure GDA0003579612610000035
式中,
Figure GDA0003579612610000036
是算法创建的航班排序子队列,表示一架架航班,来自集合E,使用ELDT 作为优先级;
Figure GDA0003579612610000037
是算法创建的航班排序子队列,来自集合F1,使用CLDT作为优先级;
Figure GDA0003579612610000038
是算法创建的航班排序子队列,来自集合F2,使用SLDT作为优先级。
对航班队列Q进行重新编号排序变为q1,q2,...,q|I|
如某繁忙机场的进场接受率AAR设置为6架/小时,开始时间=1801,结束时间=2159,则时隙资源集合S={1810,1820,…,2140,2150}。
进一步地,所述步骤7)具体包括:根据步骤5)中的航班队列Q和步骤6)中的时隙资源slots对每个航班分配落地时隙;
计算方法采用双层循环计算,外层循环是流控措施的集合,内层循环是相关航班的待计算队列;
计算遍历队列中的每个航班f∈I,从剩余可用时隙集合S中,计算分配给该航班的时隙;
for f=q1,q2,...,q|I|
let t=min(t∈S:t≥earliest_CLDT and status(t)=open)
CLDTf=t
Figure GDA0003579612610000041
flight(t)=f
owner(t)=Airlinef
status(t)=filed
S=S-{t}
end for
式中,flight(t)表示时隙属于哪个航班,owner(t)表示时隙属于哪个航空公司,status(t)表示时隙是否被占用。
本发明的有益效果:
1、本发明方法减少了流量管理措施MIT从源头往下游层层转发,层层加码和时隙浪费等情况。
2、本发明方法采用时刻表时间(SOBT/SLDT)作为判断航班先到先服务的基准,避免对航空公司二次惩罚的现象。
3、本发明方法采用了时隙资源属于航空公司规则,提高航空公司协同配合意愿,使其对航空公司具备对该时隙的拥有权,可用于时隙交换等。
附图说明
图1为本发明方法的流程图。
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
术语解释:
Figure GDA0003579612610000042
Figure GDA0003579612610000051
参照图1所示,本发明的一种航班流量时隙全局优化调度方法,步骤如下:
1)确定受流控影响的航班集合范围;
确定包含在航班优化排序中的航班计划集合I,
Figure GDA0003579612610000052
I=A∪B
A={f∈F:t0≤ELDTf≤t1}
B={f∈F:ELDTf>t1 and(t0≤CLDTf≤t1 or t0≤SLDTf≤t1)}
式中,航班集合A表示航班的预计落地时间落入流控时间范围的航班,航班集合B表示预计落地时间晚于t1,且计划落地时间(SLDT)落入流控范围或已有计算落地时间(CLDT)在流控时间范围;F表示航班集合,f表示每架航班,ELDTf表示航班f预计落地时间,CLDTf表示航班f的计算落地时间,SLDTf表示航班f的时刻表落地时间,t0为流控开始时间,t1为流控结束时间。
2)人工定义豁免航班;
由流量人员定义航班集合E中的航班豁免,可直接指定某些航班豁免,或通过定义经过某些航路/航路点或某些机场起飞的航班集合豁免:
Figure GDA0003579612610000061
3)将非豁免(但包含在流控影响时间和空间范围内)航班分类处理;
已经被分配过计算落地时间(至少一次)的航班集合;非豁免但包含航班中除去已分配过时隙的航班集合;
F1={f∈I-E:f had a CLDT}
F2=(I-E)-F1
I=E∪F1∪F2
式中,F1表示受影响航班集合I中,除去步骤2中豁免航班和已拥有了计算落地时间时隙的航班,F2表示受影响航班集合I中,除去豁免航班和F1航班集合。
4)计算每个航班的最早计算落地时间的时间阀值;
41)若航班是人工操作的豁免航班(如VIP或应急管理航班),则对于已分配时隙资源的航班,最早CLDT则取预计落地时间;对于未分配时隙资源的航班,则最早CLDT取预计落地时间或计划落地时间中的较小值;
对于已错过其预计撤轮档时间(EOBT)的航班,考虑最早CLDT的阀值,如当前时间>EOBT, 则最早CLDT=当前时间+VSP+空中飞行时间;
Figure GDA0003579612610000062
式中,CLDT为计算落地时间,slot表示航班所占用的时隙资源,EOBT表示预计离港时间, slotID表示时隙资源分配后的编号,SIBT表示计划到港时间,VIT表示航班在停机位和跑道头之间的可变滑行时间,VSP为系统可变参数;
42)若航班属于集合F1和F2,则:
当前时间current_time小于或等于SOBT,则最早CLDT=SIBT -VTT-vsp1;
当前时间current_time大于SOBT,则最早CLDT=当前时间+空中飞行时间+vsp2
Figure GDA0003579612610000063
式中,SOBT表示航班计划离港时间,vsp1、vsp2为用户定义适应性参数,VTT表示可变滑行时间,EET表示空中飞行所需时间。
5)创建航班排序队列Q;
51)根据航班的ELDT作为优先级,创建豁免航班的排序子队列Q(E);
52)创建已分配时隙的航班队列,根据已拥有的CLDT先后,创建航班排序子队列Q(F1);
53)创建剩下航班(即除去豁免航班集合和已有时隙航班集合,剩余受控航班)的子队列Q(F2),根据航班的SLDT初步排序;
54)将上述步骤51)、52)、53)中队列,重新排序合成一个带标识属性的排序总队列Q;
Figure GDA0003579612610000071
Figure GDA0003579612610000072
Figure GDA0003579612610000073
Figure GDA0003579612610000074
式中,
Figure GDA0003579612610000075
是算法创建的航班排序子队列,表示一架架航班,来自集合E,使用ELDT 作为优先级;
Figure GDA0003579612610000076
是算法创建的航班排序子队列,来自集合F1,使用CLDT作为优先级;
Figure GDA0003579612610000077
是算法创建的航班排序子队列,来自集合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中,计算分配给该航班的时隙;
for f=q1,q2,...,q|I|
let t=min(t∈S:t≥earliest_CLDT and status(t)=open)
CLDTf=t
Figure GDA0003579612610000081
flight(t)=f
owner(t)=Airlinef
status(t)=filed
S=S-{t}
end for
式中,flight(t)表示时隙属于哪个航班,owner(t)表示时隙属于哪个航空公司,status(t)表示时隙是否被占用。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。

Claims (2)

1.一种航班流量时隙全局优化调度方法,其特征在于,步骤如下:
1)确定受流控影响的航班集合范围;
2)人工定义豁免航班;
3)将非豁免航班分类处理;
4)计算每个航班的最早计算落地时间的时间阈值;
5)创建航班排序队列Q;
6)根据用户输入的机场进场接收率分配虚拟的时隙资源;
7)对每个航班分配计算落地时间;
所述步骤1)具体包括:确定包含在航班优化排序中的航班计划集合I,
Figure FDA0003579612600000011
I=A∪B
A={f∈F:t0≤ELDTf≤t1}
B={f∈F:ELDTf>t1 and(t0≤CLDTf≤t1 or t0≤SLDTf≤t1)}
式中,航班集合A表示航班的预计落地时间落入流控时间范围的航班,航班集合B表示预计落地时间晚于t1,且计划落地时间落入流控范围或已有计算落地时间在流控时间范围;F表示航班集合,f表示每架航班,ELDTf表示航班f预计落地时间,CLDTf表示航班f的计算落地时间,SLDTf表示航班f的时刻表落地时间,t0为流控开始时间,t1为流控结束时间;
所述步骤3)具体包括:已经被分配过计算落地时间的航班集合;非豁免但包含航班中除去已分配过时隙的航班集合;
F1={f∈I-E:f had a CLDT}
F2=(I-E)-F1
I=E∪F1∪F2
式中,F1表示受影响航班集合I中,除去步骤2)中豁免航班和已拥有了计算落地时间时隙的航班,F2表示受影响航班集合I中,除去豁免航班和F1航班集合;
所述步骤4)具体包括:
41)若航班是人工操作的豁免航班,则对于已分配时隙资源的航班,最早CLDT则取预计落地时间;对于未分配时隙资源的航班,则最早CLDT取预计落地时间或计划落地时间中的较小值;
对于已错过其预计撤轮档时间的航班,考虑最早CLDT的阈值,若当前时间>EOBT,则最早CLDT=当前时间+VSP+空中飞行时间;
Figure FDA0003579612600000021
式中,CLDT为计算落地时间,slot表示航班所占用的时隙资源,EOBT表示预计离港时间,slotID表示时隙资源分配后的编号,SIBT表示计划到港时间,VIT表示航班在停机位和跑道头之间的可变滑行时间,VSP为系统可变参数;
42)若航班属于集合F1和F2,则:
当前时间current_time小于或等于SOBT,则最早CLDT=SIBT -VTT-vsp1;
当前时间current_time大于SOBT,则最早CLDT=当前时间+空中飞行时间+vsp2
Figure FDA0003579612600000022
式中,SOBT表示航班计划离港时间,vsp1、vsp2为用户定义适应性参数,VTT表示可变滑行时间,EET表示空中飞行所需时间;
所述步骤5)具体包括:
51)根据航班的ELDT作为优先级,创建豁免航班的排序子队列Q(E);
52)创建已分配时隙的航班队列,根据已拥有的CLDT先后,创建航班排序子队列Q(F1);
53)创建剩下航班的子队列Q(F2),根据航班的SLDT初步排序;
54)将上述步骤51)、52)、53)中队列,重新排序合成一个带标识属性的排序总队列Q;
Figure FDA0003579612600000023
Figure FDA0003579612600000024
Figure FDA0003579612600000025
Figure FDA0003579612600000026
式中,
Figure FDA0003579612600000027
是算法创建的航班排序子队列,表示一架架航班,来自集合E,使用ELDT作为优先级;
Figure FDA0003579612600000028
是算法创建的航班排序子队列,来自集合F1,使用CLDT作为优先级;
Figure FDA0003579612600000029
是算法创建的航班排序子队列,来自集合F2,使用SLDT作为优先级;
所述步骤7)具体包括:根据步骤5)中的航班队列Q和步骤6)中的时隙资源对每个航班分配落地时隙;
计算方法采用双层循环计算,外层循环是流控措施的集合,内层循环是相关航班的待计算队列;
计算遍历队列中的每个航班f∈I,从剩余可用时隙集合S中,计算分配给该航班的时隙。
2.根据权利要求1所述的航班流量时隙全局优化调度方法,其特征在于,所述步骤2)具体包括:由流量人员定义航班集合E中的航班豁免,可直接指定某些航班豁免,或通过定义经过某些航路/航路点或某些机场起飞的航班集合豁免:
Figure FDA0003579612600000031
CN202110494816.4A 2021-05-07 2021-05-07 一种航班流量时隙全局优化调度方法 Active CN113205707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110494816.4A CN113205707B (zh) 2021-05-07 2021-05-07 一种航班流量时隙全局优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110494816.4A CN113205707B (zh) 2021-05-07 2021-05-07 一种航班流量时隙全局优化调度方法

Publications (2)

Publication Number Publication Date
CN113205707A CN113205707A (zh) 2021-08-03
CN113205707B true CN113205707B (zh) 2022-06-07

Family

ID=77029282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110494816.4A Active CN113205707B (zh) 2021-05-07 2021-05-07 一种航班流量时隙全局优化调度方法

Country Status (1)

Country Link
CN (1) CN113205707B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114299762B (zh) * 2021-12-03 2022-11-08 南京莱斯信息技术股份有限公司 一种基于一致性监测的航班时隙自动优化方法及系统
CN115035748B (zh) * 2022-06-07 2024-04-02 成都民航空管科技发展有限公司 民航航班流量控制态势监控方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106709594A (zh) * 2016-11-18 2017-05-24 中国民航大学 一种基于排队论的航班滑出时间预测方法
CN107393348A (zh) * 2017-07-10 2017-11-24 南京航空航天大学 一种信息共享机制下的进离场航班协同排序方法
CN108460994B (zh) * 2018-01-23 2021-10-08 广州市中南民航空管通信网络科技有限公司 一种航班流控排序方法、电子设备及存储介质
CN109949617B (zh) * 2019-03-20 2020-11-27 中国电子科技集团公司第二十八研究所 一种基于动态规划的航班进场流量管理策略生成方法
CN112261109B (zh) * 2020-10-16 2021-08-24 中国民用航空华东地区空中交通管理局 一种基于区块链的多机场时隙交换系统和方法
CN112419791B (zh) * 2020-11-09 2022-05-10 南京莱斯信息技术股份有限公司 一种区域级繁忙终端航班排序调度方法

Also Published As

Publication number Publication date
CN113205707A (zh) 2021-08-03

Similar Documents

Publication Publication Date Title
CN103426331B (zh) 多机场协同放行系统航班排序决策方法
Filar et al. How airlines and airports recover from schedule perturbations: a survey
CN111160770B (zh) 一种进场航班动态协同排序方法
Pujet et al. Input-output modeling and control of the departure process of congested airports
CN113205707B (zh) 一种航班流量时隙全局优化调度方法
Bennell et al. Airport runway scheduling
AU2007324694B2 (en) Method and device for the control of air traffic management at an airport
US6584400B2 (en) Schedule activated management system for optimizing aircraft arrivals at congested airports
Anagnostakis et al. Runway operations planning and control sequencing and scheduling
CN107591034A (zh) 一种进港航班排序的实现方法
CN109711619B (zh) 考虑空地运行限制的多机场战略航班时刻协同优化方法
CN102737524A (zh) 用于飞行器轨迹管理的方法和系统
WO2018196700A1 (zh) 基于候补机制的航班分配系统和方法
CN109740871B (zh) 一种基于缓冲机制的航班时隙资源利用方法及对应系统
Mesgarpour et al. Models for aircraft landing optimization
CN113112874B (zh) 一种航路时隙与高度层协同优化分配方法
CN109544000A (zh) 面向航班正常性的航空公司排班计划优化方法和系统
CN116307542A (zh) 基于多计量点约束的进港排序方法及电子终端
Kwasiborska Development of an algorithm for determining the aircraft pushback sequence
CN114299762B (zh) 一种基于一致性监测的航班时隙自动优化方法及系统
CN110428665B (zh) 一种航路与机场时隙协同分配的随机双层规划方法
Çeçen et al. Optimising aircraft arrivals in terminal airspace by mixed integer linear programming model
CN117540997A (zh) 基于航空公司航班消减偏好的航班预先飞行计划优化方法
CN116129678B (zh) 基于航迹运行的进港航班时隙自主协商分配系统及方法
CN109598985B (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