CN103413462A - 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法 - Google Patents

一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法 Download PDF

Info

Publication number
CN103413462A
CN103413462A CN2013103013259A CN201310301325A CN103413462A CN 103413462 A CN103413462 A CN 103413462A CN 2013103013259 A CN2013103013259 A CN 2013103013259A CN 201310301325 A CN201310301325 A CN 201310301325A CN 103413462 A CN103413462 A CN 103413462A
Authority
CN
China
Prior art keywords
flight
time
sector
path
air traffic
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.)
Granted
Application number
CN2013103013259A
Other languages
English (en)
Other versions
CN103413462B (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201310301325.9A priority Critical patent/CN103413462B/zh
Publication of CN103413462A publication Critical patent/CN103413462A/zh
Application granted granted Critical
Publication of CN103413462B publication Critical patent/CN103413462B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,包括:考虑地面等待、改航、空中调速和空中等待等空中流量调控措施,建立一个多目标空中交通网络流优化模型;基于时空分离思想,提出路径时间分配系统算法,通过分别对航班进行路径搜索和时间优化操作,对多目标模型进行求解;同时设计预选算子嵌入路径时间分配系统算法中,提取有效路径解,减小搜索空间,提高运行效率。本发明可以同时最小化空域拥挤和航班延误,且运行效率高,符合空中交通网络流优化的实时性要求。

Description

一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法
技术领域
本发明涉及一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,属于空中交通流量管理领域。
背景技术
随着民航事业的飞速发展,空中交通流量急剧增加,空域拥挤与大规模航班延误现象也日趋频繁。因此,空中交通流量管理(Air Traffic Flow Management;以下简称:ATFM),即统一组织、全局优化空中交通流量的运行,最大限度地提高航路网络的利用率和飞行流量的运行效率,缓解空域拥挤和航班延误,已成为空中交通管理领域的一个热点问题。
空中交通网络流优化(Air Traffic Network Flow Optimization;以下简称:ATNFO)是ATFM领域中的典型的优化问题,具体通过合理调配航路网络中每一航班的起飞降落时间和飞行路径,以最大限度缓解空中交通拥堵程度和航班的总体飞行延误。现有工作主要是将空中交通流量优化问题建模为单目标数学规划模型,包括:0-1整数规划模型,以机场容量和扇区容量为约束,通过优化航班的起降时间最小化航班延误;混合0-1整数规划模型,将航班延误代价分为地面等待代价和空中延误代价,在求解过程中均衡航班各类延误代价。然而,缓解空中交通拥堵与减少航班延误在实际中往往是相互冲突的,因此ATNFO问题实质上是一个多目标空中交通网络流优化(Multi-objective Air Traffic Network FlowOptimization;以下简称:MATNFO)问题。并且,在实际中空中交通拥挤很难用线性函数来表示,故整数规划模型无法用来准确刻画MATNFO问题。
现有MATNFO问题研究以最小化航班延误和最小化交通拥挤度为目标建模,从随机优化的角度使用多目标遗传算法(Multi-Objective Genetic Algorithm;以下简称:MOGA)求解得出最优的起飞时间和航班路径。该模型仅支持地面等待和改航两种空中流量调控措施,以缓解交通拥挤和减少航班延误,并没有考虑空中调速和空中等待对缓解拥挤和延误的作用。
发明内容
本发明的目的在于解决上述问题,建立一个更加实用的MATNFO模型,支持地面等待、改航、空中调速和空中等待等流量调控措施,并提供一种适用于求解MATNFO模型的路径时间分配(Route and Time-slot Assignment;以下简称:RTA)系统算法,获取MATNFO问题最优解,同时满足实时性要求。
本发明的一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,包括以下几个步骤:
步骤1、建立一个通用的双目标空中交通网络流优化数学模型,同时优化空域拥挤和航班延误;
步骤2、基于时空间分离的问题求解思想,提出一个路径时间分配系统算法,通过分别对航班进行路径搜索和时间优化操作,对双目标空中交通流量优化问题求最优解;同时,设计预选算子提取有效路径解,将预选算子嵌入路径时间分配系统算法中,以减小搜索空间,提高问题求解效率。
本发明提供的空中交通网络流优化方法,与现有技术相比,取得了如下技术效果:
(1)建立了更加符合实际的模型,支持地面等待、改航、空中调速和空中等待等空中流量调控措施,可以确保飞行的安全和高效;
(2)RTA是一种包括路径搜索模块、时间优化模块和预选算子的系统算法,可以同时最小化空域拥挤和航班延误,其中预选算子可以减小解的搜索空间,相比现有的MOGA、非支配排序遗传算法—Ⅱ(Non-dominated Sorting Genetic Algorithm II;以下简称:NSGA-II)和基于分解的多目标遗传算法(Multi-objective Evolutionary Algorithm basedon Decomposition;以下简称:MOEA/D)提高了解的质量的同时也提高了运行效率,符合MATNFO实时性要求。
附图说明
图1为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例中空域航路网络和扇区划分示意图;
图2为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例的流程图;
图3为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例中预选算子提取有效路径解流程图;
图4为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例中时间优化模块流程图;
图5为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例中时间优化模块中时间优化示意图;
图6为本发明综合考虑空域拥堵和航班延误的空中交通网络流量优化方法实施例中预选算子流程图;
图7为本发明的方法流程图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
图1为本发明的航路网络和扇区划分示意图。航路网络中节点由机场和航路点构成,边由航路段构成,所有节点依次按阿拉伯数字编号,如图1中数字1到17所示。空域共划分为4个扇区,依次表示为S1,S2,S3,S4。当航班在航路网络中飞行时,其路径可以由所经过航路点或所经过扇区表示,故图1中航班f的路径可以表示成航路网络节点路径
rf=(2,16,14,13,5),或扇区路径sf=(S2,S1,S4)。
本发明的一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,流程如图7所示,包括以下几个步骤:
步骤1、建立一个通用的双目标空中交通网络流优化数学模型,同时优化空域拥挤和航班延误;
本发明中问题描述为:给定航路网络、相应的扇区结构以及航班初始飞行计划,优化航路网络中飞行的每一航班的起飞降落时间和飞行路径,最大限度减少航路网络的交通拥挤程度和航班的总体飞行延误。下面将对本发明中MATNFO数学模型的进行描述,该模型的建立基于如下假设:
(1)所有航班的起飞降落时间和飞行路径在一定的范围内可选;
(2)同一起飞机场—降落机场(Original–Destination;以下简称:OD)对之间的航班,可选路径集合相同。
本发明中MATNFO数学模型的参数如表一所示:
表1 参数含义
Figure BDA00003528954000041
本发明同时最小化空域拥挤和航班延误,即目标函数为:
(1)最小化空域拥挤(Airspace Congestion;以下简称:AC)
Figure BDA00003528954000042
其中,
Figure BDA00003528954000043
表示扇区Sk在时段T内总工作负荷,表示扇区Sk在时段T内最拥挤时刻的工作负荷,即最大工作负荷;
Figure BDA00003528954000045
给扇区总工作负荷和最大工作负荷分配不同的权重。
W S k ( t ) = W mo S k ( t ) + W co S k ( t ) + W cf S k ( t ) , 为监视负荷、协调负荷及冲突负荷之和。
(2)最小化航班延误(Total Flight Delay;以下简称:TFD)
min TFD = Σ f ∈ F | c f g · | τ f d - τ f 0 d | + c f a · ( τ f a - τ f d - t f min ) | - - - ( 2 )
其中,航班飞行延误成本由两部分构成:航班由于延迟起飞所产生的地面等待成本
Figure BDA00003528954000051
Figure BDA00003528954000052
为单位时间地面等待成本;航班由于空中等待、减速调速或改航所产生的空中延误成本
Figure BDA00003528954000053
Figure BDA00003528954000054
为单位时间空中飞行成本,为航班f在网络中沿着所选飞行路径到达目的机场的最短飞行时间。
相对于已有双目标数学建模工作中仅考虑航班地面等待和改航带来的延误代价,本发明引入了航班降落时间f∈F的优化,可更全面考虑航班由于地面等待、改航、空中调速和空中等待等空中流量调控措施带来的延误代价,从而使得我们的数学模型更符合于实际空中交通流量管理问题。
模型中约束条件:
(1)rf∈Rf,其中
Figure BDA00003528954000057
表示航班f从起飞机场到降落机场的可选路径集合。|Rf|是可选路径集合的数目;
(2) τ f d ∈ [ T ‾ f d , T ‾ f d ] τ f a ∈ [ T ‾ f a , T ‾ f a ] , 其中
Figure BDA000035289540000510
Figure BDA000035289540000511
为航班f的起飞、降落时间, [ T ‾ f d , T ‾ f d ]
Figure BDA000035289540000513
分别为航班有效的起飞、降落时间范围;
(3)
Figure BDA000035289540000514
其中L(rf)是航路rf的长度,[vmin,vmax]是飞机飞行速度的变化范围。
综上,本发明中通用的双目标空中交通网络流优化数学模型表示为:
目标函数:
Figure BDA000035289540000515
min TFD = Σ f ∈ F | c f g · | τ f d - τ f 0 d | + c f a · ( τ f a - τ f d - t f min ) |
约束条件:
rf∈Rf,f∈F
τ f d ∈ [ T ‾ f d , T ‾ f d ] , f ∈ F
τ f a ∈ [ T ‾ f a , T ‾ f a ] , f ∈ F
v f = L ( r f ) τ f a - τ f d ∈ [ v min , v max ] , f ∈ F
步骤2、基于时空间分离的问题求解思想,提出一个路径时间分配系统算法,通过分别对航班进行路径搜索和时间优化操作,对双目标空中交通流量优化问题求最优解,同时,设计预选算子提取有效路径解,将预选算子嵌入路径时间分配系统算法中,以减小搜索空间,提高问题求解效率;
流程如图2所示,具体包括如下步骤:
步骤201,初始化航班路径选择概率模型P={P1,P2,...,Pf,...,P|F|}为均匀分布,根据飞行计划初始化种群Pop(g),初始化非支配解集EPop为空,设定RTA演化代数(Generation;以下简称:g)g=0;其中|F|为总航班数,
Figure BDA00003528954000061
f=1,2,...,|F|表示航班f可选路径集rf中各路径(即
Figure BDA00003528954000062
)被航班f选择飞行的概率向量,
|Rf|为航班f可选路径规模,且
Figure BDA00003528954000063
步骤202,根据航班路径选择概率模型P生成航班当前最优飞行路径,更新Pop(g);
步骤203,令种群SPop(g)=Φ,设计预选算子从种群Pop(g)中提取有效路径解放入种群SPop(g)中;
步骤204,将种群SPop(g)中每一个体进行时间优化,即利用时间平移和时间缩放两个操作对每一个体中时间参数进行优化;其中时间平移操作用来优化航班起飞时间以均衡各扇区负荷,减少空域拥挤;时间缩放操作用来优化航班降落时间以减少航班延误;
步骤205,更新种群SPop(g)为SPop(g)=SPop(g-1)∪EPop,计算种群SPop(g)中的每一个体目标函数值,并根据目标函数值对SPop(g)中个体进行非支配排序;清空非支配排序解集EPop,将种群SPop(g)中前50%的个体复制到EPop中,统计EPop中个体路径解情况,根据公式(3)对航班路径选择概率模型P进行更新;
P f ( r f = r f k ) = 1 NP R Σ i = 1 NP R find ( r f = r f k ) - - - ( 3 )
其中,NPR为更新航班路径选择概率模型P的个体数目,即EPop中个体数目;find为逻辑式,当A为真时,find(A)=1,当A为假时,find(A)=0;k=1,...,|Rf|;f=1,2,...,|F|。
步骤206,判断g是否满足指定演化代数Gen,如果演化代数太少,优化不够充分,航班路径选择概率模型P不能收敛,所需优化的目标函数(航班拥挤度和航班总延误)也不能得到较优的解,故本发明Gen=1000;若是,优化结束;此时EPop中的非支配解构成最终的非支配解集,非支配解集中的个体即为优化的航班起飞降落时间和飞行路径;否则令g=g+1,返回步骤202。
所述的步骤203,流程如图3所示,具体包括如下步骤:
步骤2031,输入种群Pop(g);
步骤2032,计算Pop(g)中的每一个体的扇区负荷均衡指数;扇区负荷均衡指数定义如下:
首先,建立扇区—时间网络结构,其中每个节点表示在一定时间内的一个扇区,每条边表示相邻节点间有航班飞行;
其次,定义扇区—时间网络中节点介数为经过节点的航班数目;定义扇区负荷均衡指数为扇区—时间网络节点介数的标准方差值,计算如公式(4)所示:
SDB = 1 N T · N S Σ j = 1 N T Σ i = 1 N S ( b ij - b ‾ ) 2 - - - ( 4 )
其中,bij表示扇区—时间网络中扇区i时隙j节点的介数;b为扇区—时间网络中所有节点介数的平均值;NS和NT表示包括的扇区和时隙数目;
根据复杂网络理论,网络中节点介数差异越小,网络节点容量越均衡,故对于扇区—时间网络而言,扇区负荷均衡指数越小,扇区间负荷越均衡;
步骤2033,根据扇区负荷均衡指数值对Pop(g)中个体进行升序排序,选择前α·NP个个体,作为有效路径解,复制到种群SPop(g)中,其中α是值为0-1的系数,本发明中α=0.5;
步骤2034,输出种群SPop(g);因SPop(g)中个体均为Pop(g)中较优个体,从而可缩小RTA算法的搜索空间。
所述的步骤204,流程如图4所示,具体包括如下步骤:
步骤2041,搜索负荷峰值最大扇区,提取经过该扇区的航班,根据“削峰填谷”策略,判断航班起飞时间需提前还是延误,计算相应平移步长,执行时间平移操作,均衡各扇区负荷;
航班起飞时间提前或延误判断规则为:当满足公式(5)时航班起飞时间提前,否则延误。
|tin-tpeak|>|tout-tpeak|   (5)
其中,tpeak为扇区Sk中负荷峰值时刻,tin和tout分别为航班f进入和离开扇区Sk的时刻。时间平移的步长计算方法如公式(6)所示:
M = max { x | W S k ( t out - i ) - W S k ( t in - i - 1 ) > 0 , i &Element; [ 0 , x ] } , | t in - t peak | &GreaterEqual; | t out - t peak | max { x | W S k ( t in + i ) - W S k ( t out + i + 1 ) > 0 , i &Element; [ 0 , x ] } , | t in - t peak | < | t out - t peak | - - - ( 6 )
其中,M代表步长数,单位步长即单位时隙。
步骤2042,时间缩放操作是指航班在经过路径搜索和时间平移操作后,根据已优化的航班路径和起飞时间,分析航班最优到达时间,对航班飞行时间进行缩放操作,优化航班的降落时间,达到减少航班延误的目的;实际空中交通流量管理中,可通过航班调速或空中等待等流量调控措施改变航班降落时间;确定最优降落时间如公式(7)所示:
Figure BDA00003528954000082
其中,
Figure BDA00003528954000083
表示航班f经过时间缩放操作后更新的起飞时间;
Figure BDA00003528954000084
表示航班f时间缩放操作后更新的降落时间;
Figure BDA00003528954000085
表示是航班f在
Figure BDA00003528954000086
时刻起飞,以最大的飞行速度沿着选定的航路rf到达目的机场的降落时间;
Figure BDA00003528954000087
为没有航班延误时的降落时间。
图5为本发明时间平移操作示意图,图中所示为某一扇区内不同时刻的负荷值。假定航班f在第8个时隙进入(tin=8)该扇区,在第15个时隙离开(tout=15)该扇区,且在第13个时隙(tpeak=13)时达到负荷峰值。通过对航班f实施时间平移操作,即航班进入该扇区时间提前5个单位时隙,使得第11到第15个时隙间的负荷减少了,从而时间平移确可以减少空域拥挤。
图6为本发明时间优化示意图,图中所示为某一航班在更新飞行路径后,起飞降落时间优化示意图;图中
Figure BDA00003528954000088
Figure BDA00003528954000089
分别表示航班f当前起飞时间和降落时间;
Figure BDA000035289540000810
Figure BDA000035289540000811
分别表示航班f经过时间平移操作和时间缩放操作后更新的起飞时间;
Figure BDA000035289540000812
Figure BDA000035289540000813
表示航班f经过时间平移操作和时间缩放操作后更新的降落时间,并且有
本实施例提供了一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,建立了更加实用的MATNFO模型,考虑了地面等待、改航、空中调速和空中等待等ATFM策略,提出了一种更适用求解MATNFO问题的RTA算法。RTA运用时空分离的思想,分别通过路径搜索和时间优化模块对路径和起飞、降落时间进行优化,并通过预选择模块减小搜索空间,提高运行效率,满足实时性要求。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。

Claims (3)

1.一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,包括以下几个步骤:
步骤1、建立通用的双目标空中交通网络流优化数学模型,同时优化空域拥挤和航班延误;
本发明的双目标空中交通网络流优化数学模型建立基于如下假设:
(1)所有航班的起飞降落时间和飞行路径在一定的范围内可选;
(2)同一起飞机场—降落机场,简称OD,之间的航班可选路径集合相同;
本发明同时最小化空域拥挤和航班延误,即目标函数为:
(1)最小化空域拥挤,简称:AC,具体为:
Figure FDA00003528953900011
其中,
Figure FDA00003528953900012
表示扇区Sk在时段T内总工作负荷,
Figure FDA00003528953900013
表示扇区Sk在时段T内最拥挤时刻的工作负荷,即最大工作负荷;
Figure FDA00003528953900014
给扇区总工作负荷和最大工作负荷分配不同的权重;T表示时间集合,被划分成等长、连续的离散时隙nT表示时间段T被划分为nT个时隙,t1表示第1个时隙,
Figure FDA00003528953900015
表示第nT个时隙;
W S k ( t ) = W mo S k ( t ) + W co S k ( t ) + W cf S k ( t ) , 为时隙t监视负荷、协调负荷及冲突负荷之和;nS表示空域中扇区数目;
(2)最小化航班延误,简称TFD,具体为:
min TFD = &Sigma; f &Element; F | c f g &CenterDot; | &tau; f d - &tau; f 0 d | + c f a &CenterDot; ( &tau; f a - &tau; f d - t f min ) | - - - ( 2 )
其中,df表示航班f的起飞机场,af表示航班f的目的机场,
Figure FDA00003528953900018
表示航班f的起飞时间,
Figure FDA00003528953900019
表示航班f的降落时间,F表示航班集合,|F|表示总航班数目,rf0,
Figure FDA000035289539000110
分别表示初始飞行计划中的航班f的飞行路径、起飞时间和降落时间;
航班飞行延误成本由两部分构成:航班由于延迟起飞所产生的地面等待成本
Figure FDA000035289539000111
为单位时间地面等待成本;航班由于空中等待、调速或改航所产生的空中延误成本
Figure FDA000035289539000112
为单位时间空中飞行成本,为航班f在网络中沿着所选飞行路径到达目的机场的最短飞行时间;FPL表示航班飞行计划,
FPL={fplf|f=1,...,|F|},
Figure FDA00003528953900021
其中rf表示航班f的飞行路径;
模型中约束条件:
(1)rf∈Rf,其中
Figure FDA00003528953900022
表示航班f从起飞机场到降落机场的可选路径集合;|Rf|是可选路径集合的数目;
(2) &tau; f d &Element; [ T &OverBar; f d , T &OverBar; f d ] &tau; f a &Element; [ T &OverBar; f a , T &OverBar; f a ] , 其中
Figure FDA00003528953900025
为航班f的起飞、降落时间, [ T &OverBar; f d , T &OverBar; f d ]
Figure FDA00003528953900028
分别为航班有效的起飞、降落时间范围;
(3)
Figure FDA00003528953900029
其中vf表示航班f的平均飞行速度,L(rf)是航路rf的长度,[vmin,vmax]是飞机飞行速度的变化范围;
综上,本发明中通用的双目标空中交通网络流优化数学模型表示为:
目标函数:
Figure FDA000035289539000210
min TFD = &Sigma; f &Element; F | c f g &CenterDot; | &tau; f d - &tau; f 0 d | + c f a &CenterDot; ( &tau; f a - &tau; f d - t f min ) |
约束条件:
rf∈Rf,f∈F
&tau; f d &Element; [ T &OverBar; f d , T &OverBar; f d ] , f &Element; F
&tau; f a &Element; [ T &OverBar; f a , T &OverBar; f a ] , f &Element; F
v f = L ( r f ) &tau; f a - &tau; f d &Element; [ v min , v max ] , f &Element; F
步骤2、基于时空间分离的问题求解思想,提出一个路径时间分配系统算法,通过分别对航班进行路径搜索和时间优化操作,对双目标空中交通流量优化问题求最优解,同时,设计预选算子提取有效路径解,将预选算子嵌入路径时间分配系统算法中;
具体包括如下步骤:
步骤201,初始化航班路径选择概率模型P={P1,P2,...,Pf,...,P|F|}为均匀分布,根据飞行计划初始化种群Pop(g),初始化非支配解集EPop为空,设定演化代数g=0;其中|F|为总航班数,
Figure FDA00003528953900031
f=1,2,...,|F|表示航班f可选路径集Rf中各路径被航班f选择飞行的概率向量,可选路径集Rf中各路径即为
Figure FDA00003528953900032
|Rf|为航班f可选路径规模,且 &Sigma; i = 1 | R f | p i = 1 ;
步骤202,根据航班路径选择概率模型P生成航班当前最优飞行路径,更新Pop(g);
步骤203,令种群SPop(g)=Φ,设计预选算子从种群Pop(g)中提取有效路径解放入种群SPop(g)中;
步骤204,将种群SPop(g)中每一个体进行时间优化,即利用时间平移和时间缩放两个操作对每一个体中时间参数进行优化;
步骤205,更新种群SPop(g)为SPop(g)=SPop(g-1)∪EPop,计算种群SPop(g)中的每一个体目标函数值,并根据目标函数值对SPop(g)中个体进行非支配排序;清空非支配排序解集EPop,将种群SPop(g)中前50%的个体复制到EPop中,统计EPop中个体路径解情况,根据公式(3)对航班路径选择概率模型P进行更新;
P f ( r f = r f k ) = 1 NP R &Sigma; i = 1 NP R find ( r f = r f k ) - - - ( 3 )
其中,NPR为更新航班路径选择概率模型P的个体数目,即EPop中个体数目;find为逻辑式,当A为真时,find(A)=1,当A为假时,find(A)=0;k=1,...,|Rf|;f=1,2,...,|F|;
步骤206,判断g是否满足指定演化代数Gen,若是,优化结束;此时EPop中的非支配解构成最终的非支配解集,非支配解集中的个体即为优化的航班起飞降落时间和飞行路径;否则令g=g+1,返回步骤202。
2.根据权利要求1所述的一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,所述的步骤203,具体包括如下步骤:
步骤2031,输入种群Pop(g);
步骤2032,计算Pop(g)中的每一个体的扇区负荷均衡指数,扇区负荷均衡指数定义如下:
首先,建立扇区—时间网络结构,其中每个节点表示在一定时间内的一个扇区,每条边表示相邻节点间有航班飞行;
其次,定义扇区—时间网络中节点介数为经过节点的航班数目;定义扇区负荷均衡指数为扇区—时间网络节点介数的标准方差值,计算如公式(4)所示:
SDB = 1 N T &CenterDot; N S &Sigma; j = 1 N T &Sigma; i = 1 N S ( b ij - b &OverBar; ) 2 - - - ( 4 )
其中,bij表示扇区—时间网络中扇区i时隙j节点的介数;
Figure FDA00003528953900042
为扇区—时间网络中所有节点介数的平均值;NS和NT表示包括的扇区和时隙数目;
根据复杂网络理论,网络中节点介数差异越小,网络节点容量越均衡,故对于扇区—时间网络而言,扇区负荷均衡指数越小,扇区间负荷越均衡;
步骤2033,根据扇区负荷均衡指数值对Pop(g)中个体进行升序排序,选择前α·NP个个体,作为有效路径解,复制到种群SPop(g)中,其中α是值为0-1的系数;
步骤2034,输出种群SPop(g)。
3.根据权利要求1所述的一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法,所述的步骤204,具体包括如下步骤:
步骤2041,搜索负荷峰值最大扇区,提取经过该扇区的航班,根据“削峰填谷”策略,判断航班起飞时间需提前还是延误,计算相应平移步长,执行时间平移操作,均衡各扇区负荷;
航班起飞时间提前或延误判断规则为:当满足公式(5)时航班起飞时间提前,否则延误;
|tin-tpeak|>|tout-tpeak|   (5)
其中,tpeak为扇区Sk中负荷峰值时刻,tin和tout分别为航班f进入和离开扇区Sk的时刻;
时间平移的步长计算方法如公式(6)所示:
M = max { x | W S k ( t out - i ) - W S k ( t in - i - 1 ) > 0 , i &Element; [ 0 , x ] } , | t in - t peak | &GreaterEqual; | t out - t peak | max { x | W S k ( t in + i ) - W S k ( t out + i + 1 ) > 0 , i &Element; [ 0 , x ] } , | t in - t peak | < | t out - t peak | - - - ( 6 )
其中,M代表步长数,单位步长即单位时隙;
步骤2042,时间缩放操作是指航班在经过路径搜索和时间平移操作后,根据已优化的航班路径和起飞时间,分析航班最优到达时间,对航班飞行时间进行缩放操作,优化航班的降落时间,确定最优降落时间如公式(7)所示:
其中,
Figure FDA00003528953900052
表示航班f经过时间缩放操作后更新的起飞时间;表示航班f经过时间缩放操作后更新的降落时间;表示是航班f在
Figure FDA00003528953900055
时刻起飞,以最大的飞行速度沿着选定的航路rf到达目的机场的降落时间;为没有航班延误时的降落时间。
CN201310301325.9A 2013-07-18 2013-07-18 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法 Active CN103413462B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310301325.9A CN103413462B (zh) 2013-07-18 2013-07-18 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310301325.9A CN103413462B (zh) 2013-07-18 2013-07-18 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法

Publications (2)

Publication Number Publication Date
CN103413462A true CN103413462A (zh) 2013-11-27
CN103413462B CN103413462B (zh) 2016-01-20

Family

ID=49606467

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310301325.9A Active CN103413462B (zh) 2013-07-18 2013-07-18 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法

Country Status (1)

Country Link
CN (1) CN103413462B (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504198A (zh) * 2014-12-19 2015-04-08 北京航空航天大学 一种基于双层协同进化的航路网络拓扑设计方法
CN105117792A (zh) * 2015-08-06 2015-12-02 北京航空航天大学 一种考虑跑道口等待时长的航班机场场面运行优化方法
CN105390030A (zh) * 2014-09-03 2016-03-09 北京航空航天大学 飞行流量调控方法及系统
CN105489066A (zh) * 2014-09-16 2016-04-13 北京航空航天大学 空中交通流量调控方法
CN105513429A (zh) * 2014-09-24 2016-04-20 北京航空航天大学 飞行冲突解脱方法及装置
CN106297414A (zh) * 2015-06-05 2017-01-04 北京航空航天大学 航班流量的调控方法和装置
CN107103169A (zh) * 2017-06-26 2017-08-29 上海交通大学 一种用于满足出行续航要求的交通网络均衡计算方法
CN107591034A (zh) * 2017-08-09 2018-01-16 南京莱斯信息技术股份有限公司 一种进港航班排序的实现方法
CN107862910A (zh) * 2017-11-15 2018-03-30 南京航空航天大学 基于机场场面滑行多种冲突模式的滑行延误时间和离场延误时间预测方法
CN108960533A (zh) * 2018-08-06 2018-12-07 中国民航大学 一种空中航路网优化方法
CN109598985A (zh) * 2019-01-14 2019-04-09 南京航空航天大学 航路资源协同分配方法
CN109598984A (zh) * 2019-01-14 2019-04-09 南京航空航天大学 航路资源优化分配系统
CN109979245A (zh) * 2017-12-28 2019-07-05 北京航空航天大学 飞行流量调控方法和装置
CN110349444A (zh) * 2018-04-06 2019-10-18 杭州坚果壳科技开发有限公司 基于大数据的空中交通流量管理方法
CN110706520A (zh) * 2019-10-31 2020-01-17 中国电子科技集团公司第二十八研究所 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法
JP2021517290A (ja) * 2019-04-08 2021-07-15 中国電子科技集団公司第二十八研究所The 28Th Research Institute Of China Electronics Technology Group Corporation 対流気象条件での最適通抜け経路に基づく航路の渋滞度の評価方法
CN113140134A (zh) * 2021-03-12 2021-07-20 北京航空航天大学 一种智慧空管系统的体系架构
CN113379065A (zh) * 2021-05-17 2021-09-10 百融云创科技股份有限公司 一种基于多目标语法进化的自动机器学习方法
CN113408907A (zh) * 2021-06-22 2021-09-17 中国人民解放军空军工程大学 空中交通管制系统威胁影响分析方法、系统及设备
WO2022194310A1 (zh) * 2021-08-18 2022-09-22 中国电子科技集团公司第二十八研究所 一种航班飞行轨迹多目标动态规划方法
RU2788101C1 (ru) * 2022-09-07 2023-01-16 Акционерное общество "Ордена Трудового Красного Знамени Всероссийский научно-исследовательский институт радиоаппаратуры" (АО "ВНИИРА") Способ управления потоком взлетов и посадок летательных аппаратов и устройство для его осуществления
CN116246496A (zh) * 2023-05-10 2023-06-09 北京航空航天大学 一种基于航班与流量协同环境的飞行计划规划方法
CN116910465A (zh) * 2023-05-31 2023-10-20 南京航空航天大学 基于机器学习的空中交通流量管理延误分布预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000260000A (ja) * 1999-03-05 2000-09-22 Nec Corp 航空交通流調整システム
CN101533566A (zh) * 2009-02-23 2009-09-16 民航数据通信有限责任公司 预测空中交通流量的方法和装置
CN101540113A (zh) * 2009-03-16 2009-09-23 民航数据通信有限责任公司 空中交通航路点流量控制方法、服务器及系统
CN102222412A (zh) * 2011-05-26 2011-10-19 北京航空航天大学 一种引入空域容量的航路汇聚点布局优化方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000260000A (ja) * 1999-03-05 2000-09-22 Nec Corp 航空交通流調整システム
CN101533566A (zh) * 2009-02-23 2009-09-16 民航数据通信有限责任公司 预测空中交通流量的方法和装置
CN101540113A (zh) * 2009-03-16 2009-09-23 民航数据通信有限责任公司 空中交通航路点流量控制方法、服务器及系统
CN102222412A (zh) * 2011-05-26 2011-10-19 北京航空航天大学 一种引入空域容量的航路汇聚点布局优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张强,曹义华,陈勇: "民用机场及附近空域空中交通流量优化", 《飞行力学》, vol. 25, no. 1, 30 March 2007 (2007-03-30), pages 84 - 87 *

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105390030A (zh) * 2014-09-03 2016-03-09 北京航空航天大学 飞行流量调控方法及系统
CN105489066A (zh) * 2014-09-16 2016-04-13 北京航空航天大学 空中交通流量调控方法
CN105513429A (zh) * 2014-09-24 2016-04-20 北京航空航天大学 飞行冲突解脱方法及装置
CN104504198B (zh) * 2014-12-19 2017-10-10 北京航空航天大学 一种基于双层协同进化的航路网络拓扑设计方法
CN104504198A (zh) * 2014-12-19 2015-04-08 北京航空航天大学 一种基于双层协同进化的航路网络拓扑设计方法
CN106297414A (zh) * 2015-06-05 2017-01-04 北京航空航天大学 航班流量的调控方法和装置
CN106297414B (zh) * 2015-06-05 2019-03-05 北京航空航天大学 航班流量的调控方法和装置
CN105117792B (zh) * 2015-08-06 2018-08-07 北京航空航天大学 一种考虑跑道口等待时长的航班机场场面运行优化方法
CN105117792A (zh) * 2015-08-06 2015-12-02 北京航空航天大学 一种考虑跑道口等待时长的航班机场场面运行优化方法
CN107103169A (zh) * 2017-06-26 2017-08-29 上海交通大学 一种用于满足出行续航要求的交通网络均衡计算方法
CN107103169B (zh) * 2017-06-26 2020-03-24 上海交通大学 一种用于满足出行续航要求的交通网络均衡计算方法
CN107591034A (zh) * 2017-08-09 2018-01-16 南京莱斯信息技术股份有限公司 一种进港航班排序的实现方法
CN107862910A (zh) * 2017-11-15 2018-03-30 南京航空航天大学 基于机场场面滑行多种冲突模式的滑行延误时间和离场延误时间预测方法
CN107862910B (zh) * 2017-11-15 2019-12-06 南京航空航天大学 基于机场场面滑行多种冲突模式的滑行延误时间和离场延误时间预测方法
CN109979245B (zh) * 2017-12-28 2020-10-30 北京航空航天大学 飞行流量调控方法和装置
CN109979245A (zh) * 2017-12-28 2019-07-05 北京航空航天大学 飞行流量调控方法和装置
CN110349444B (zh) * 2018-04-06 2023-11-10 杭州坚果壳科技开发有限公司 基于大数据的空中交通流量管理方法
CN110349444A (zh) * 2018-04-06 2019-10-18 杭州坚果壳科技开发有限公司 基于大数据的空中交通流量管理方法
CN108960533A (zh) * 2018-08-06 2018-12-07 中国民航大学 一种空中航路网优化方法
CN108960533B (zh) * 2018-08-06 2021-08-17 中国民航大学 一种空中航路网优化方法
CN109598985B (zh) * 2019-01-14 2022-05-20 南京航空航天大学 航路资源协同分配方法
CN109598985A (zh) * 2019-01-14 2019-04-09 南京航空航天大学 航路资源协同分配方法
CN109598984A (zh) * 2019-01-14 2019-04-09 南京航空航天大学 航路资源优化分配系统
JP2021517290A (ja) * 2019-04-08 2021-07-15 中国電子科技集団公司第二十八研究所The 28Th Research Institute Of China Electronics Technology Group Corporation 対流気象条件での最適通抜け経路に基づく航路の渋滞度の評価方法
JP7077413B2 (ja) 2019-04-08 2022-05-30 中国電子科技集団公司第二十八研究所 対流気象条件での最適通抜け経路に基づく航路の渋滞度の評価方法
CN110706520A (zh) * 2019-10-31 2020-01-17 中国电子科技集团公司第二十八研究所 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法
CN110706520B (zh) * 2019-10-31 2021-08-31 中国电子科技集团公司第二十八研究所 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法
CN113140134B (zh) * 2021-03-12 2022-07-08 北京航空航天大学 一种面向智慧空管系统的航班延误智能预测框架
CN113140134A (zh) * 2021-03-12 2021-07-20 北京航空航天大学 一种智慧空管系统的体系架构
CN113379065A (zh) * 2021-05-17 2021-09-10 百融云创科技股份有限公司 一种基于多目标语法进化的自动机器学习方法
CN113408907B (zh) * 2021-06-22 2023-05-23 中国人民解放军空军工程大学 空中交通管制系统威胁影响分析方法、系统及设备
CN113408907A (zh) * 2021-06-22 2021-09-17 中国人民解放军空军工程大学 空中交通管制系统威胁影响分析方法、系统及设备
WO2022194310A1 (zh) * 2021-08-18 2022-09-22 中国电子科技集团公司第二十八研究所 一种航班飞行轨迹多目标动态规划方法
RU2788101C1 (ru) * 2022-09-07 2023-01-16 Акционерное общество "Ордена Трудового Красного Знамени Всероссийский научно-исследовательский институт радиоаппаратуры" (АО "ВНИИРА") Способ управления потоком взлетов и посадок летательных аппаратов и устройство для его осуществления
CN116246496A (zh) * 2023-05-10 2023-06-09 北京航空航天大学 一种基于航班与流量协同环境的飞行计划规划方法
CN116246496B (zh) * 2023-05-10 2023-09-08 北京航空航天大学 一种基于航班与流量协同环境的飞行计划规划方法
CN116910465A (zh) * 2023-05-31 2023-10-20 南京航空航天大学 基于机器学习的空中交通流量管理延误分布预测方法
CN116910465B (zh) * 2023-05-31 2024-06-11 南京航空航天大学 基于机器学习的空中交通流量管理延误分布预测方法

Also Published As

Publication number Publication date
CN103413462B (zh) 2016-01-20

Similar Documents

Publication Publication Date Title
CN103413462A (zh) 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法
CN102222412B (zh) 一种引入空域容量的航路汇聚点布局优化方法
CN102254453B (zh) 一种民航多机场终端区空域功能性扇区划分方法
Liang et al. Integrated sequencing and merging aircraft to parallel runways with automated conflict resolution and advanced avionics capabilities
CN104751683B (zh) 机场航班滑行调度指挥系统及方法
CN109191925B (zh) 一种面向四维航迹运行的多空域轨迹规划与协商方法
CN103226899A (zh) 基于空中交通特征的空域扇区动态划分方法
CN107591034A (zh) 一种进港航班排序的实现方法
García et al. Optimization of airport ground operations integrating genetic and dynamic flow management algorithms
Baldi et al. Simulation-based synthesis for approximately optimal urban traffic light management
CN107025806B (zh) 一种单阶段临时航迹鲁棒优化方法
CN103473955A (zh) 一种基于图论和谱聚类算法的终端扇区划分方法
CN110428665B (zh) 一种航路与机场时隙协同分配的随机双层规划方法
CN103942623A (zh) 一种基于需求与容量不确定性的机场拥挤风险预测方法
CN105117792A (zh) 一种考虑跑道口等待时长的航班机场场面运行优化方法
Daniel et al. Airspace congestion smoothing by multi-objective genetic algorithm
CN114446094A (zh) 一种航班排序信息的时空转换方法
CN114626201A (zh) 一种模型化开放式空管运行仿真系统与方法
Sergeeva et al. Dynamic airspace configurations generated by evolutionary algorithms
CN105390030B (zh) 飞行流量调控方法及系统
CN110033111B (zh) 机场场面运动规划方法和装置
Ghorpade Airspace configuration model using swarm intelligence based graph partitioning
CN105489066B (zh) 空中交通流量调控方法
Koeners et al. Improving taxi traffic flow by real-time runway sequence optimization using dynamic taxi route planning
CN110046748B (zh) 一种通航作业的运力与任务匹配的遗传算法优化方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant