CN109118112B - 跨区域尾随间隔限制及离场时隙分配综合策略计算方法 - Google Patents

跨区域尾随间隔限制及离场时隙分配综合策略计算方法 Download PDF

Info

Publication number
CN109118112B
CN109118112B CN201811001806.7A CN201811001806A CN109118112B CN 109118112 B CN109118112 B CN 109118112B CN 201811001806 A CN201811001806 A CN 201811001806A CN 109118112 B CN109118112 B CN 109118112B
Authority
CN
China
Prior art keywords
time
flight
flow control
interval
ctot
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
CN201811001806.7A
Other languages
English (en)
Other versions
CN109118112A (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 University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201811001806.7A priority Critical patent/CN109118112B/zh
Publication of CN109118112A publication Critical patent/CN109118112A/zh
Application granted granted Critical
Publication of CN109118112B publication Critical patent/CN109118112B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G5/00Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/003Flight plan management

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及航空领域,具体而言涉及一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法,通过构建跨区域尾随间隔限制及离场时隙分配混合整数规划模型并进行求解,以计算出跨区域尾随间隔限制及离场时隙分配综合策略。针对航班在下游出区域流控点受到的约束生成MIT受控航班在上游进入区域流控点的尾随间隔限制策略以及区域内外起飞时隙受限航班的计算起飞时隙策略,基于实施管理中的运行限制和优化目标的分析,以及对计算起飞时隙和尾随间隔限制的实施特征分析,建立混合整数规划模型,进行求解,发现综合考虑尾随间隔限制及离场时隙分配比传统的先到先服务管理效率更高,减少总体航班延误成本。

Description

跨区域尾随间隔限制及离场时隙分配综合策略计算方法
技术领域
本发明涉及航空领域,具体而言,涉及一种区域尾随间隔限制及离场时隙分配综合策略计算方法及系统。
背景技术
下游出区域点受到外界流控时,为解决下游尾随间隔向上游放大传播,导致上游航班受到过大的流控限制的问题,建立尾随间隔限制与起飞离场时隙协同管理程序,对过出区域点受流控影响的航班分别采用尾随间隔和离场时隙分配策略实施管理。传统的先到先服务管理策略,按照计划到达下游点的先后顺序分配过点时隙,通常不会考虑到航班的延误成本,另外传统尾随间隔管理策略不考虑交通流的具体分布差异,因而计算上游尾随间隔过于粗放,导致航班产生大量不必要的延误。
发明内容
本发明的目的是提供一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法。
为了实现上述目的,本发明实施例采用的技术方案如下:
本发明实施例提供了一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法,包括构建用于跨区域尾随间隔限制及离场时隙分配的混合整数规划模型;以及对混合整数规划模型进行求解。
在本发明较佳的实施例中,构建所述混合整数规划模型的方法包括:
对混合整数规划模型涉及参数进行定义;
确定混合整数规划模型的决策变量;
定义混合整数规划模型的约束条件;以及
构建混合整数规划模型的目标函数。
在本发明较佳的实施例中,所述对混合整数规划模型涉及参数进行定义的方法包括:
Fctot:受CTOT影响的航班集合;
Fmit:受MIT影响的航班集合;
m:受CTOT和MIT影响的航班总数,m=|Fctot|+|Fmit|;
n:受CTOT影响的航班数量,n=|Fctot|;
Figure BDA0001781719640000021
受CTOT影响的航班i预计起飞时间i∈(1,...n),fi∈Fctot
Figure BDA0001781719640000022
受上游MIT影响的航班i预计过上游流控点B的时间i∈(n+1,...m),fi∈Fmit,使得
Figure BDA0001781719640000023
Figure BDA0001781719640000024
按升序进行排列;
Figure BDA0001781719640000025
航班i预计过下游流控点C的时间i∈(1,...m),fi∈Fctot∪Fmit
tbegin:上游流控的开始时间;
c1:地面延误成本;
c2:区域外空中延误成本;
c3:区域内空中延误成本;
a:空中可消耗延误占空中飞行时间百分比a∈[0,1];
b:下游流控限制时间间隔值;
d:单位流控时间间隔值,流控值通常为单位流控间隔值的倍数;
k:上游流控时段分为K个时间片,k∈(1,...K),表示第k个时间片;
h:上游流控从开始时间到结束时间的总时段长。
在本发明较佳的实施例中,所述确定决策变量的方法包括:
所述确定混合整数规划模型的决策变量的方法包括:
Figure BDA0001781719640000026
表示受CTOT影响的航班fi分配到的起飞时间,i∈(1,...n),fi∈Fctot
Figure BDA0001781719640000031
表示区域外航班fi受上游流控限制影响下的过点时间,i∈(n+1,...m),fi∈Fmit
Figure BDA0001781719640000032
表示航班f受下游流控限制影响下的过点时间,i∈(1,...m),fi∈FctotU Fmit
nk,k∈(1,...K):表示第k个时间片上游的流控间隔值为单位流控时间间隔值的倍数;
Figure BDA0001781719640000033
上式中,i∈(n+1,...m),k∈(1,...K)。
在本发明较佳的实施例中,所述定义混合整数规划模型的约束条件的方法包括:
步骤S1:受CTOT影响的航班分配的起飞时间大于等于预计起飞时间,即
Figure BDA0001781719640000034
步骤S2:航班受上游流控限制影响下的进入区域时间不早于预计进入时间,即
Figure BDA0001781719640000035
步骤S3:航班离开区域的时间不早于预计离开区域的时间,即
Figure BDA0001781719640000036
步骤S4:受CTOT影响的航班从起飞到出区域的飞行时间不小于预计飞行时间,同时不大于预计飞行时间加最大可消耗延误时间,即
Figure BDA0001781719640000037
步骤S5:区域外进入的航班在区域内的飞行时间不小于预计区域内飞行时间,同时不大于预计区域内飞行时间加最大可消耗延误时间,即
Figure BDA0001781719640000038
步骤S6:出区域航班之间的间隔时间不小于流控限制的间隔时间b,即
Figure BDA0001781719640000041
步骤S7:过上游流控点B的航班在唯一的一个时间片经过该上游流控点B,即
Figure BDA0001781719640000042
步骤S8:航班所属时间片与航班安排的过点时间之间满足的对应关系,即
Figure BDA0001781719640000043
步骤S9:航班过上游流控点B的时间为航班预计过点时间以及与前机满足流控间隔dnk的时间两者取大,即
Figure BDA0001781719640000044
fi,fj∈Fmit
在本发明较佳的实施例中,所述构建混合整数规划模型的目标函数的方法包括:
目标函数为最小化所有航班在地面及空中的总延误成本,目标函数中第一项表示受CTOT影响的航班在起飞机场的地面延误成本,第二项表示受MIT影响的航班分别在区域内和区域外的空中延误成本,即:
Figure BDA0001781719640000045
在本发明较佳的实施例中,所述对混合整数规划模型进行求解的方法,即
采用遗传算法与启发式算法相结合的混合算法对所述混合整数规划模型进行求解。
在本发明较佳的实施例中,所述采用遗传算法与启发式算法相结合的混合算法包括:
遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码;
混合编码后的航班处理顺序及流控间隔值作为启发式算法的输入,由启发式算法根据编码的航班顺序及上游流控间隔值计算航班起飞及进出点时隙,并返回适应度值给遗传算法;
遗传算法依据接收到的适应度值对混合整数规划模型进行求解。
在本发明较佳的实施例中,所述遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码的方法包括:
对所有出区域的航班处理顺序和上游流控间隔值分别编码;
假设有m个受CTOT和MIT影响的航班和k个时间片,则基于混合编码随机产生m+k个随机数,其中m作为航班处理顺序的随机数且采用十进制正整数排列编码方法,k作为各个时段流控间隔值的随机数且采用十进制随机正整数编码。
在本发明较佳的实施例中,遗传算法包括:
所述遗传算法中染色体对应的适应度函数为
Figure BDA0001781719640000051
Figure BDA0001781719640000052
其中y为中间变量,Z为混合整数规划模型中目标函数值,c4为对区内延误超出最大可消耗延误的惩罚系数,适应度函数值由启发式算法根据编码的航班顺序及上游流控间隔作为输入条件计算获得;
遗传算法选择算子:采用随机竞争对初始种群进行选择操作,即每次按轮盘赌选择机制选取一对个体,然后让这两个个体进行竞争,适应度高的个体被选中,重复计算,直到选满为止;
遗传算法交叉算子:对航班顺序采用部分映射杂交,即随机选择两个交叉点,子代基因中两个交叉点间的部分来自一个父代基因,其余部分来自另外一个父代基因;对上游流控间隔采用单点交叉,即随机选择一个交叉点,然后在该点互换两个父代染色体的部分基因;
遗传算法变异算子:航班序列部分采用交换基因值的变异方法,随机选择两个交换变异点,将其对应的基因进行互换;以及上游管制移交间隔适于采用基本位变异算子操作,随机选择变异点,取随机数来代替原有值。
相对于现有技术,本发明实施例具有以下有益效果:
本发明实施例提供了一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法,跨区域尾随间隔限制及离场时隙分配综合策略计算方法通过构建跨区域尾随间隔限制及离场时隙分配混合整数规划模型并进行求解,以计算出跨区域尾随间隔限制及离场时隙分配综合策略。针对航班在下游出区域流控点受到的约束生成MIT受控航班在上游进入区域流控点的尾随间隔限制策略以及区域内外起飞时隙受限航班的计算起飞时隙策略,基于实施管理中的运行限制和优化目标的分析,以及对计算起飞时隙和尾随间隔限制的实施特征分析,建立混合整数规划模型,进行求解,发现综合考虑尾随间隔限制及离场时隙分配比传统的先到先服务管理效率更高,可减少总体的航班延误成本。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1示出了本发明实施例提供的跨区域尾随间隔限制及离场时隙分配综合策略计算方法的流程图。
图2示出了本发明实施例提供的受限航班流及空域示意图。
图3示出了本发明实施例提供的遗传算法与启发式算法相结合的混合算法流程图。
图4示出了本发明实施例提供的遗传算法编码方式示例图。
图5示出了本发明实施例提供的遗传算法交叉算子示例图。
图6示出了本发明实施例提供的遗传算法变异算子示例图。
图7示出了本发明实施例提供的启发式算法流程图。
图8示出了本发明实施例提供的跨区域尾随间隔限制及离场时隙分配综合策略计算结果图。
图9示出了本发明实施例提供的跨区域尾随间隔限制及离场时隙分配综合策略实施效果比较图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。
实施例
请参阅图1,本发明实施例提供了一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法。跨区域尾随间隔限制及离场时隙分配综合策略计算方法通过构建跨区域尾随间隔限制及离场时隙分配混合整数规划模型并进行求解,以计算出跨区域尾随间隔限制及离场时隙分配综合策略。针对航班在下游出区域流控点受到的约束生成MIT受控航班在上游进入区域流控点的尾随间隔限制策略以及区域内外起飞时隙受限航班的计算起飞时隙策略,基于实施管理中的运行限制和优化目标的分析,以及对计算起飞时隙和尾随间隔限制的实施特征分析,建立混合整数规划模型,进行求解,发现综合考虑尾随间隔限制及离场时隙分配比传统的先到先服务管理效率更高,可减少总体的航班延误成本。
其中,构建所述混合整数规划模型的方法包括:
请参阅图2,采集出区域流控限制信息,包括下游流控开始结束时间、下游尾随间隔值,受限航班的预计起飞时间、预计进区域时间、预计出区域时间等信息;
对混合整数规划模型涉及参数进行定义;
确定混合整数规划模型的决策变量;
定义混合整数规划模型的约束条件;以及
构建混合整数规划模型的目标函数。
在本实施例中,所述对混合整数规划模型涉及参数进行定义的方法包括:根据采集出区域流控限制信息,包括下游流控开始结束时间、下游尾随间隔值,受限航班的预计起飞时间、预计进区域时间、预计出区域时间等信息对下列模型参数进行赋值:
Fctot:受CTOT(计算起飞时间,欧洲交通管制术语)影响的航班集合;
Fmit:受MIT(尾随间隔限制,交通管制术语)影响的航班集合;
m:受CTOT和MIT影响的航班总数,m=|Fctot|+|Fmit|;
n:受CTOT影响的航班数量,n=|Fctot|;
Figure BDA0001781719640000081
受CTOT影响的航班i预计起飞时间i∈(1,...n),fi∈Fctot
Figure BDA0001781719640000082
受上游MIT影响的航班i预计过上游流控点B的时间i∈(n+1,...m),fi∈Fmit,使得
Figure BDA0001781719640000083
Figure BDA0001781719640000084
按升序进行排列;
Figure BDA0001781719640000085
航班i预计过下游流控点C的时间i∈(1,...m),fi∈Fctot∪Fmit
tbegin:上游流控的开始时间;
c1:地面延误成本;
c2:区域外空中延误成本;
c3:区域内空中延误成本;
a:空中可消耗延误占空中飞行时间百分比a∈[0,1];
b:下游流控限制时间间隔值;
d:单位流控时间间隔值,流控值通常为单位流控间隔值的倍数;
k:上游流控时段分为K个时间片,k∈(1,...K),表示第k个时间片;
h:上游流控从开始时间到结束时间的总时段长。
在本实施例中,所述确定混合整数规划模型的决策变量的方法包括:
Figure BDA0001781719640000091
表示受CTOT影响的航班f分配到的起飞时间,i∈(1,...n),fi∈Fctot
Figure BDA0001781719640000092
表示区域外航班f受上游流控限制影响下的过点时间,i∈(n+1,...m),fi∈Fmit
Figure BDA0001781719640000093
表示航班fi受下游流控限制影响下的过点时间,i∈(1,...m),fi∈FctotUFmit
nk,k∈(1,...K):表示第k个时间片上游的流控间隔值为单位流控时间间隔值的倍数;
Figure BDA0001781719640000094
上式中,i∈(n+1,...m),k∈(1,...K)。
在本实施例中,所述定义混合整数规划模型的约束条件的方法包括:
步骤S1:受CTOT影响的航班分配的起飞时间大于等于预计起飞时间,即
Figure BDA0001781719640000095
步骤S2:航班受上游流控限制影响下的进入区域时间不早于预计进入时间,即
Figure BDA0001781719640000096
步骤S3:航班离开区域的时间不早于预计离开区域的时间,即
Figure BDA0001781719640000101
步骤S4:受CTOT影响的航班从起飞到出区域的飞行时间不小于预计飞行时间,同时不大于预计飞行时间加最大可消耗延误时间,即
Figure BDA0001781719640000102
步骤S5:区域外进入的航班在区域内的飞行时间不小于预计区域内飞行时间,同时不大于预计区域内飞行时间加最大可消耗延误时间,即
Figure BDA0001781719640000103
步骤S6:出区域航班之间的间隔时间不小于流控限制的间隔时间b,即
Figure BDA0001781719640000104
步骤S7:过上游流控点B的航班在唯一的一个时间片经过该上游流控点B,即
Figure BDA0001781719640000105
步骤S8:航班所属时间片与航班安排的过点时间之间满足的对应关系,即
Figure BDA0001781719640000106
步骤S9:航班过上游流控点B的时间为航班预计过点时间以及与前机满足流控间隔dnk的时间两者取大,即
Figure BDA0001781719640000107
fi,fj∈Fmit
在本实施例中,所述构建混合整数规划模型的目标函数的方法包括:
目标函数为最小化所有航班在地面及空中的总延误成本,目标函数中第一项表示受CTOT影响的航班在起飞机场的地面延误成本,第二项表示受MIT影响的航班分别在区域内和区域外的空中延误成本,即:
Figure BDA0001781719640000111
请继续参阅图1,所述进行求解的方法采用遗传算法与启发式算法相结合的混合算法对跨区域尾随间隔限制及离场时隙分配混合整数规划模型进行求解。
具体地,请参阅图3,所述遗传算法与启发式算法相结合的方法包括:
遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码;
混合编码后的航班处理顺序及流控间隔值作为启发式算法的输入,由启发式算法根据编码的航班顺序及上游流控间隔值计算航班起飞及进出点时隙,并返回适应度值给遗传算法(如图3中512,498,504等数值均为适应度值);
遗传算法依据接收到的适应度值对混合整数规划模型进行求解。
遗传算法编码方式示例如图4所示:所述遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码的方法包括:对所有出区域的航班处理顺序和上游流控间隔值分别编码。
假设有m个受CTOT和MIT影响的航班和k个时间片,则基于混合编码随机产生m+k个随机数,其中m可以作为航班处理顺序的随机数且采用十进制正整数排列编码方法,k可以作为各个时段流控间隔值的随机数且采用十进制随机正整数编码。
其中,启发式算法根据混合编码后的航班处理顺序及流控间隔值计算适应度值。
所述遗传算法中染色体对应的适应度函数为
Figure BDA0001781719640000112
Figure BDA0001781719640000121
其中y为中间变量,Z为混合整数规划模型中目标函数值,c4为对区内延误超出最大可消耗延误的惩罚系数,适应度函数值由启发式算法根据编码的航班顺序及上游流控间隔作为输入条件计算获得;
遗传算法选择算子:采用随机竞争对初始种群进行选择操作,即每次按轮盘赌选择机制选取一对个体,然后让这两个个体进行竞争,适应度高的个体被选中,重复计算,直到选满为止;
遗传算法交叉算子(如图5所示):对航班顺序采用部分映射杂交,即随机选择两个交叉点(如图5中第二到第四数字列),子代基因中两个交叉点间的部分来自一个父代基因,其余部分来自另外一个父代基因,对上游流控间隔采用单点交叉,即随机选择一个交叉点(如图5中第十一数字列),然后在该点互换两个父代染色体的部分基因。
遗传算法变异算子(如图6所示):航班序列部分采用交换基因值的变异方法,随机选择两个交换变异点(如图6中第二数字列和第四数字列),将其对应的基因进行互换,既加快了算法速度也可以保证变异后个体的可行性;上游管制移交间隔采用基本位变异算子操作,随机选择变异点(如图6中第十数字列),取随机数来代替原有值。
启发式算法流程如图7所示,算法输入、输出以及计算流程描述如下:
算法输入:航班排序序列,下游流控限制时间间隔值,上游分时段流控限制值;
算法输出:受CTOT影响的航班分配的起飞时间和出区域时间,MIT受控航班的进区域时间和出区域时间,上游尾随间隔值,适应度函数值;
计算流程:
对航班排序序列F中的顺序进行处理,使得MIT受控航班的顺序保持根据
Figure BDA0001781719640000131
的升序顺序,根据处理后的序列对航班进行依次处理;
对于受CTOT影响的航班,在出区域流控点C找到满足下游流控约束的时隙
Figure BDA0001781719640000132
减去从起飞机场到下游点的飞行时间
Figure BDA0001781719640000133
即得到航班分配的起飞时间
Figure BDA0001781719640000134
对于MIT受控航班,计算满足上游流控点B的流控约束下的上游过点时间
Figure BDA0001781719640000135
根据
Figure BDA0001781719640000136
计算到达出区域流控点C的最早到达时间
Figure BDA0001781719640000137
以这个时间为初始时间,找到航班满足在C点流控限制的出区域时间
Figure BDA0001781719640000138
对所有航班依次处理完后,根据
Figure BDA0001781719640000139
基于适应度函数计算公式计算适应度函数值(图3中“过点时间”指
Figure BDA00017817196400001310
三类过点时间)。
请参阅图8,针对不同交通流,传统的先到先服务策略与本方法所生成的跨区域尾随间隔限制及离场时隙分配综合策略结果对比见图9,横轴表示不满足下游流控限制的航班占航班总数的比例,纵轴为根据(先到先服务策略下目标函数值-本模型算法下目标函数值)/先到先服务策略下目标函数值计算出的比率值,结果对比显示随着不满足流控限制航班数量的增加,本模型算法计算出的策略结果比先到先服务策略优化的程度也增加。
综上所述,本发明实施例提供了一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法,跨区域尾随间隔限制及离场时隙分配综合策略计算方法通过构建跨区域尾随间隔限制及离场时隙分配混合整数规划模型并进行求解,以计算出跨区域尾随间隔限制及离场时隙分配综合策略。针对航班在下游出区域流控点受到的约束生成MIT受控航班在上游进入区域流控点的尾随间隔限制策略以及区域内外起飞时隙受限航班的计算起飞时隙策略,基于实施管理中的运行限制和优化目标的分析,以及对计算起飞时隙和尾随间隔限制的实施特征分析,建立混合整数规划模型,进行求解,发现综合考虑尾随间隔限制及离场时隙分配比传统的先到先服务管理效率更高,可减少总体的航班延误成本。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种跨区域尾随间隔限制及离场时隙分配综合策略计算方法,其特征在于,
构建用于跨区域尾随间隔限制及离场时隙分配的混合整数规划模型;以及
对混合整数规划模型进行求解;
构建所述混合整数规划模型的方法包括:
对混合整数规划模型涉及参数进行定义;
确定混合整数规划模型的决策变量;
定义混合整数规划模型的约束条件;以及
构建混合整数规划模型的目标函数;
所述对混合整数规划模型涉及参数进行定义的方法包括:
Fctot:受计算起飞时间CTOT影响的航班集合;
Fmit:受尾随间隔限制MIT影响的航班集合;
m:受CTOT和MIT影响的航班总数,m=|Fctot|+|Fmit|;
n:受CTOT影响的航班数量,n=|Fctot|;
Figure FDA0003097315220000011
受CTOT影响的航班fi预计起飞时间,i∈(1,...n),fi∈Fctot
Figure FDA0003097315220000012
受上游MIT影响的航班fi预计过上游流控点B的时间,i∈(n+1,...m),fi∈Fmit,使得
Figure FDA0003097315220000013
Figure FDA0003097315220000014
按升序进行排列;
Figure FDA0003097315220000015
航班fi预计过下游流控点C的时间,i∈(1,...m),fi∈FctotUFmit
tbegin:上游流控的开始时间;
c1:地面延误成本;
c2:区域外空中延误成本;
c3:区域内空中延误成本;
a:空中可消耗延误占空中飞行时间百分比,a∈[0,1];
b:下游流控限制时间间隔值;
d:单位流控时间间隔值;
k:上游流控时段分为K个时间片,k∈(1,...K),表示第k个时间片;
h:上游流控从开始时间到结束时间的总时段长;
所述确定混合整数规划模型的决策变量的方法包括:
Figure FDA0003097315220000021
表示受CTOT影响的航班fi分配到的起飞时间,i∈(1,...n),fi∈Fctot
Figure FDA0003097315220000022
表示区域外航班fi受上游流控限制影响下的过点时间,i∈(n+1,...m),fi∈Fmit
Figure FDA0003097315220000023
表示航班fi受下游流控限制影响下的过点时间,i∈(1,...m),fi∈FctotUFmit
nk,k∈(1,...K):表示第k个时间片上游的流控间隔值为单位流控时间间隔值的倍数;
Figure FDA0003097315220000024
上式中,i∈(n+1,...m),k∈(1,...K);
所述定义混合整数规划模型的约束条件的方法包括:
步骤S1:受CTOT影响的航班分配的起飞时间大于等于预计起飞时间,即
Figure FDA0003097315220000025
步骤S2:航班受上游流控限制影响下的进入区域时间不早于预计进入时间,即
Figure FDA0003097315220000026
步骤S3:航班离开区域的时间不早于预计离开区域的时间,即
Figure FDA0003097315220000027
步骤S4:受CTOT影响的航班从起飞到出区域的飞行时间不小于预计飞行时间,同时不大于预计飞行时间加最大可消耗延误时间,即
Figure FDA0003097315220000031
步骤S5:区域外进入的航班在区域内的飞行时间不小于预计区域内飞行时间,同时不大于预计区域内飞行时间加最大可消耗延误时间,即
Figure FDA0003097315220000032
步骤S6:出区域航班之间的间隔时间不小于流控限制的间隔时间b,即
Figure FDA0003097315220000033
步骤S7:过上游流控点B的航班在唯一的一个时间片经过该上游流控点B,即
Figure FDA0003097315220000034
步骤S8:航班所属时间片与航班安排的过点时间之间满足的对应关系,即
Figure FDA0003097315220000035
步骤S9:航班过上游流控点B的时间为航班预计过点时间以及与前机满足流控间隔dnk的时间两者取大,即
Figure FDA0003097315220000036
所述构建混合整数规划模型的目标函数的方法包括:
目标函数z为最小化所有航班在地面及空中的总延误成本,目标函数中第一项表示受CTOT影响的航班在起飞机场的地面延误成本,第二项表示受MIT影响的航班分别在区域内和区域外的空中延误成本,即
Figure FDA0003097315220000041
2.如权利要求1所述的跨区域尾随间隔限制及离场时隙分配综合策略计算方法,其特征在于,
所述对混合整数规划模型进行求解的方法,即
采用遗传算法与启发式算法相结合的混合算法对所述混合整数规划模型进行求解。
3.如权利要求2所述的跨区域尾随间隔限制及离场时隙分配综合策略计算方法,其特征在于,
所述采用遗传算法与启发式算法相结合的混合算法包括:
遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码;
混合编码后的航班处理顺序及流控间隔值作为启发式算法的输入,由启发式算法根据编码的航班顺序及上游流控间隔值计算航班起飞及进出点时隙,并返回适应度值给遗传算法;
遗传算法依据接收到的适应度值对混合整数规划模型进行求解。
4.如权利要求3所述的跨区域尾随间隔限制及离场时隙分配综合策略计算方法,其特征在于,
所述遗传算法对所有受外界流控的航班处理顺序及上游尾随间隔值进行混合式编码的方法包括:
对所有出区域的航班处理顺序和上游流控间隔值分别编码;
假设有m个受CTOT和MIT影响的航班和k个时间片,则基于混合编码随机产生m+k个随机数,其中m作为航班处理顺序的随机数且采用十进制正整数排列编码方法,k作为各个时段流控间隔值的随机数且采用十进制随机正整数编码。
5.如权利要求4所述的跨区域尾随间隔限制及离场时隙分配综合策略计算方法,其特征在于,
所述遗传算法还适于根据适应度值选择个体,适应度计算函数为
Figure FDA0003097315220000051
Figure FDA0003097315220000052
其中目标函数取值为z;目标函数取中间变量,以满足MIT受控航班在区内延误不超出最大可消耗延误;c4为对区内延误超出最大可消耗延误的惩罚系数;设计染色体对应的适应度函数为:
Figure FDA0003097315220000053
适应度函数值由启发式算法根据编码的航班顺序及上游流控间隔作为输入条件计算获得;
并且采用随机竞争对初始种群进行选择操作,即每次按轮盘赌选择机制选取一对个体,然后让这两个个体进行竞争,适应度高的个体被选中,重复计算,直到选满为止;
遗传算法交叉算子:对航班顺序采用部分映射杂交,即随机选择两个交叉点,子代基因中两个交叉点间的部分来自一个父代基因,其余部分来自另外一个父代基因;对上游流控间隔采用单点交叉,即随机选择一个交叉点,然后在该点互换两个父代染色体的部分基因;
遗传算法变异算子:航班序列部分采用交换基因值的变异方法,随机选择两个交换变异点,将其对应的基因进行互换;以及上游管制移交间隔适于采用基本位变异算子操作,随机选择变异点,取随机数来代替原有值。
CN201811001806.7A 2018-08-29 2018-08-29 跨区域尾随间隔限制及离场时隙分配综合策略计算方法 Active CN109118112B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811001806.7A CN109118112B (zh) 2018-08-29 2018-08-29 跨区域尾随间隔限制及离场时隙分配综合策略计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811001806.7A CN109118112B (zh) 2018-08-29 2018-08-29 跨区域尾随间隔限制及离场时隙分配综合策略计算方法

Publications (2)

Publication Number Publication Date
CN109118112A CN109118112A (zh) 2019-01-01
CN109118112B true CN109118112B (zh) 2021-08-24

Family

ID=64860557

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811001806.7A Active CN109118112B (zh) 2018-08-29 2018-08-29 跨区域尾随间隔限制及离场时隙分配综合策略计算方法

Country Status (1)

Country Link
CN (1) CN109118112B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110660274B (zh) * 2019-09-30 2020-09-25 中国电子科技集团公司第二十八研究所 一种机场容需平衡预测方法
CN110706520B (zh) * 2019-10-31 2021-08-31 中国电子科技集团公司第二十八研究所 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法
CN112233462B (zh) * 2020-09-30 2022-03-25 中国民用航空飞行学院 一种飞机着陆的动态时间间隔计算方法与系统
CN112396871B (zh) * 2020-10-21 2022-05-06 南京莱斯信息技术股份有限公司 基于航迹预测的进场延误分配吸纳方法
CN112634664B (zh) * 2020-12-15 2022-10-04 沈阳民航东北通信网络有限公司 基于空中交通流量管理的航班放行排序方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2891259A4 (en) * 2012-08-31 2016-10-05 Nokia Solutions & Networks Oy OPTIMIZATIONS FOR FREQUENT TRANSMISSION OF SMALL DATA
CN103778481B (zh) * 2014-01-23 2018-04-27 南京航空航天大学 一种针对多跑道航班进离场的动态调度方法
CN105094970B (zh) * 2015-07-13 2018-07-10 西安电子科技大学 一种求解分布式系统下可分任务多趟调度模型的方法

Also Published As

Publication number Publication date
CN109118112A (zh) 2019-01-01

Similar Documents

Publication Publication Date Title
CN109118112B (zh) 跨区域尾随间隔限制及离场时隙分配综合策略计算方法
CN109118111B (zh) 尾随间隔限制及起飞时隙分配综合策略管理决策支持系统
Feng et al. Heuristic hybrid game approach for fleet condition-based maintenance planning
CN109584638B (zh) 一种面向区域网络的先期航班时刻协同优化方法
CN107392355B (zh) 一种基于差分进化算法的多机场协同调度鲁棒优化方法
CN111754039B (zh) 纯电动公交线网综合集成优化设计的方法
CN112070355A (zh) 一种机场摆渡车的分配调度方法
CN113112874B (zh) 一种航路时隙与高度层协同优化分配方法
CN107679669A (zh) 一种基于超启发式方法的机场停机位调度方法及系统
CN108846596B (zh) 分时段终端区尾随间隔进场管理策略计算方法
Wei et al. A multiple objective optimization model for aircraft arrival and departure scheduling on multiple runways
CN110706520B (zh) 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法
CN114330971A (zh) 基于历史数据和中转旅客信息分配登机口的系统与方法
CN112862258B (zh) 一种考虑旅客偏好的有限理性的航班恢复方法
CN111260161A (zh) 一种众包任务下发的方法及设备
CN113257042A (zh) 一种基于时空分布特征的多机场终端区航班流生成方法
CN117540997A (zh) 基于航空公司航班消减偏好的航班预先飞行计划优化方法
CN112598328A (zh) 卫星厅模式下中转登机口多目标分配的优化方法及系统
CN110909946B (zh) 一种基于公路换乘的航班计划优化方法
CN108898891A (zh) 终端区尾随间隔进场管理决策支持系统
CN111785092A (zh) 一种面向航班延误的机场群航班编排优化方法
Hüsselmann et al. An improved solution methodology for the urban transit routing problem
CN114117928A (zh) 一种顾及旅客中转时间的机场登机口分配方法
Li et al. Hybrid optimization of arrival and departure flight schedules
Shidong et al. Multi-objectives optimization on flights landing sequence at busy airport

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