CN107016881B - 一种多跑道机场进场航班多效能优化排序方法 - Google Patents
一种多跑道机场进场航班多效能优化排序方法 Download PDFInfo
- Publication number
- CN107016881B CN107016881B CN201710329224.0A CN201710329224A CN107016881B CN 107016881 B CN107016881 B CN 107016881B CN 201710329224 A CN201710329224 A CN 201710329224A CN 107016881 B CN107016881 B CN 107016881B
- Authority
- CN
- China
- Prior art keywords
- flight
- runway
- time
- flights
- approach
- 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
Links
- 238000012163 sequencing technique Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 238000013459 approach Methods 0.000 claims abstract description 131
- 238000011156 evaluation Methods 0.000 claims abstract description 17
- 230000010006 flight Effects 0.000 claims description 89
- 238000000034 method Methods 0.000 claims description 50
- 238000012545 processing Methods 0.000 claims description 40
- 230000001133 acceleration Effects 0.000 claims description 36
- 230000008014 freezing Effects 0.000 claims description 10
- 238000007710 freezing Methods 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 9
- 230000003111 delayed effect Effects 0.000 claims description 7
- 230000001934 delay Effects 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/02—Automatic approach or landing aids, i.e. systems in which flight data of incoming planes are processed to provide landing data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0043—Traffic management of multiple aircrafts from the ground
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
- Road Paving Structures (AREA)
Abstract
本发明公开了一种新型、高效和动态的多跑道机场进场航班多效能优化排序方法,旨在为终端区进场航班提供安全、科学、可行的进场方案,包括通过定义多效能跑道分配评价函数来为航班分配跑道及确定进场航线,并根据航班综合优先级优化延误分配、确定航班的着陆时间,进而实现机场跑道流量的均衡分布,减少航班进场飞行时间和航班延误,提高空域资源利用效率,降低管制工作负荷,进一步提高多跑道机场进场航班的运行效率。
Description
技术领域
本发明涉及航班排序方法技术领域,尤其是一种多跑道机场进场航班多效能优化排序方法。
背景技术
随着空中交通需求的不断增长,空域拥挤、航班大面积延误以及航班正点率下降等问题不断涌现。为有效缓解空域供需矛盾,国内外不断加大力度进行空域和流量管理相关技术的研究,其中,航班排序问题一直是相关研究的热点和难点。当前研究大多是通过构建多目标优化模型,设计相应的智能算法进行求解,包括遗传算法、蚁群算法、鱼群算法等,该类方法虽然能够大幅度提升优化效果,但同时存在排序结果稳定性差、计算效率低下,模型过度抽象等问题。在实际工程应用中常采用经典的先到先服务(FCFS:first come firstservice)方法进行排序,该类方法计算效率较高,但因缺乏优化效果,在航班量较大情况下容易产生过量延误现象,造成空域资源浪费。
发明内容
本发明所要解决的技术问题在于,提供一种多跑道机场进场航班多效能优化排序方法,通过定义多效能跑道分配评价函数来为航班分配跑道及确定进场航线,并根据航班综合优先级优化延误分配过程、确定航班的着陆时间。
为解决上述技术问题,本发明提供一种多跑道机场进场航班多效能优化排序方法,包括如下步骤:
(1)搜集数据信息;搜集目标机场的跑道和进场航线数据,航班计划数据,航班预计着陆时间数据;进场航班每条进场航线的进场时间;假设一条跑道和一个进场点固定一条进场航线;
(2)设置进场航班排序时间区间W;按照扩大终端的概念,当前设置为W=90min;
(3)航班数据筛选与集合划分;确定本场在排序时间区间W内的所有进场航班集合F,所有进场航班数量为n;根据航班的跑道分配情况,将所有进场航班集合划分为已分配跑道航班集合和未分配跑道航班集合两部分,其中未分配跑道航班集合包含所有的未进入进场点的进场航班;根据航班是否已经经过排序,将所有进场航班集合划分为已排序进场航班集合和未排序进场航班集合两部分,已排序航班表示航班的可行着陆时间ui和跑道已经确定,ui表示航班经过排序后得出的航班可行的着陆时间;
(4)航班进场信息整理;根据航班当前的位置信息确定航班i,i∈F的可用进场路径集合和可用跑道集合以及航班i所有可能的进场飞行时间集合
(5)跑道初次分配;对未分配跑道的航班找出航班i的最小进场路径飞行时间min(FTf),所对应的初始路径和初始跑道并计算得到航班i的预计着陆时间跑道初次分配用于确定航班初始排序次序,并非正式分配跑道,因此经跑道初次分配处理的航班仍属于跑道未分配航班,后续需经跑道再次分配处理;为保证航班安全着陆和排序结果稳定性,定义距离着陆时间小于TD(含)的进场航班为冻结航班,对应的冻结航班集合为D;其中,TD为冻结航班划分时间临界值,可根据需要自行设置,当前设置为15min;该类航班所分配的跑道、着陆时间和顺序不再改变,则距离着陆时间大于TD以上的进场航班为非冻结航班,对应的非冻结航班集合为D',F=D∪D',另di为0-1变量,di=1表示航班为冻结航班,di=0表示航班为非冻结航班;
(6)确定当前排序航班;根据跑道初次分配所得航班预计着陆时间找出未排序进场航班集合中预计着陆时间最早的一架进场航班i,
(7)跑道再次分配;
(8)确定可行着陆时间;对航班i进行冲突探测与解脱,确定航班i在跑道ri的上可行着陆时间ui;此处的冲突探测包括跑道和终端区飞行过程中的冲突;
(9)优化航班延误;对比航班i预计着陆时间ti与可行着陆时间ui,并计算航班i的着陆时间偏差Δti=ui-ti,将其作为航班的延误时间,并对其进一步优化。
优选的,步骤(7)中,跑道再次分配具体为:若进场航班i已经分配跑道r',则不再重新分配跑道,ri=r',确定航班的进场航线p'和预计着陆时间t',则pi=p',ti=t';若进场航班i还未分配跑道则结合多效能跑道分配评价函数(式(1)),计算进场航班i的多效能跑道评价系数为航班i分配跑道ri和进场航线pi,计算预计着陆时间ti;并将航班i从未分配跑道航班集合中剔除,加入到已分配跑道航班集合
多效能跑道分配评价函数考虑进场飞行时间,跑道流量均衡性和航班延误三个方面,多效能跑道评价系数的计算如下:
式中,为航班i在跑道r上进场飞行时间评价系数;为航班i可用的最短飞行时间对应航线的影响系数:μ1和k为调节参数,用于调整的增长幅度;Tmin为机场最小的进场航线飞行时间,Ti min为航班i可行的最小进场航线飞行时间;Tmax为机场最长的进场航线飞行时间,Ti max为航班i可行的最长进场航线飞行时间;为航班i在跑道r着陆的进场航线飞行时间偏差系数;μ2为调节参数,用于调整的增长幅度;为航班i在跑道r着陆对应进场航线飞行时间Ti r的偏差比;为航班i所有可行进场航线飞行时间偏离最小进场航线飞行时间的标准差;为航班i所有可用跑道数量;
为航班i在跑道r上的流量评估系数;μ3>0用于衡量航班i所有可用跑道的流量均衡性;μ3为调节参数,用于的增长幅度;C为机场容量;为航班i可用跑道集合Ri的流量标准差, 为航班i在跑道r的流量偏差系数,μ4为调节参数,用于调整的增长幅度;为航班i在跑道r的流量偏差: 为跑道r的当前已排序航班数量;为航班i在所有可用跑道上的平均流量;为航班i在跑道r着陆的准点系数;xi∈[1,2]为延误重视度系数;为航班i的自身优先级,其取值范围为{1,2,3,4},该值由航班自身类型确定,1表示最低优先级,一般为货运航班,2表示普通优先级,一般为普通客运航班,3表示较高优先级,4表示最高优先级,也可以根据需求设置和调整;ω1为调节参数,用于调整的增长幅度;航班i在跑道r着陆的准点率;为航班i在跑道r着陆所需的延误占比;为航班i在跑道r上的着陆延误;ω2+ω3=1;→ω2=0.5,ω3=0.5;为航班i在跑道r着陆所能承受的最大正延误时间;M为普通优先级航班在终端区内的最大可接收延误阀值,可自行设置,当前已经起飞的航班设置为45分钟,未起飞的航班设置为2小时;ω2,ω3是的权值;为航班i的着陆趋势优先级,则:且IF: 为航班i在跑道r上的预计着陆时间;tsys为***当前时间;tB为终端区外界时间范围对应于进场航班排序时间区间W,tB=W;
为了把统一到一个参考系,需要进行归一化处理;另 分别为hk所能取到的最小值和最大值,则:
其中,g3'=(1-ω4)·(g1-1)+ω4,ω4∈[0.5,1),
优选的,步骤(9)中,优化航班延误,具体为:
(91)优化非冻结航班的延误;若航班i为非冻结航班di=0且着陆时间偏差Δti>0;
(911)前序航班加速;计算航班i综合优先级Li,判断跑道ri上先于航班i前一架进场航班j是否满足前序航班加速条件;
航班i在跑道r上的综合优先级定义如下:
其中,为航班i的自身优先级,定义参考步骤(7)为航班i的延误动态优先级,则且 定义为航班延误占比,可知当前排序航班i在跑道r上的着陆延误, 为航班i的着陆趋势优先级,定义参见步骤(7);设置了航班的可加速最低优先级Lmin,根据需求自行设置,当前设置为Lmin=2;
前序航班加速条件如下:
1)前序航班j综合优先级小于航班i综合优先级,Li>Lj;
2)航班i综合优先级大于最低可加速优先级Lmin,Li>Lmin;
3)前序航班j为非冻结航班,dj=0;
4)前序航班j未经过延误处理,Dj=0,其中,Dj表示航班j是否经过延误处理,0表示没有经过延误处理,1则表示经过延误处理;
5)前序航班j未经过加速处理,Aj=0,其中,Aj表示航班j是否经过加速处理,0表示没有经过加速处理,1则表示经过加速处理;
6)前序航班j与航班i之间存在冲突;
若前序航班j满足上述加速条件1)-6),则对航班i和j进行加速处理和冲突探测解脱,重新计算航班i和j的可行着陆时间、综合优先级和加速处理信息ui,uj,Li,Lj,Ai,Aj,以及着陆时间偏差Δti;若不满足以上加速条件,则直接转至步骤(912);
(912)判断是否满足延误上限要求;判断航班i的着陆时间偏差Δti是否满足航班i延误上限要求,对应于航班i的可接受延误阈值
(9121)查找可替换航班;若航班i的着陆时间偏差Δti不满足延误上限 则从跑道ri上所有前序航班中找出航班k,判断航班k是否满足以下前序航班可剔除条件:
1)前序航班k优先级小于航班i优先级,Li>Lk;
2)前序航班k为非冻结航班,dk=0;
3)航班i和前序航班k之间存在冲突;
4)航班i和前序航班k的使用不同的进场航线,pi≠pk;
5)前序航班k已参与剔除的次数Nk满足航班避让次数阈值K要求,Nk<K,为防止某架航班被多次剔除;其中,航班避让次数阈值K,可根据需求自行设置,当前设置为K=2;
6)为避让航班i而被剔除的前序航班数量PNi应满足航班剔除次数阈值PK要求,PNi<PK,为防止航班i的排序结果过分偏离其初始次序;其中,航班剔除次数阈值PK,可根据需求自行设置,当前设置为PK=2;
若存在前序航班k满足上述前序航班可剔除条件1)-6),则将前序航班k剔除,更新uk之后的已排序航班信息,包括可行着陆时间,着陆时间偏差,更新航班避让次数信息Nk,清空航班k可行着陆时间uk,将航班k从已排序航班集合中剔除,加入到未排序航班集合然后判断剔除航班k是否已经飞越进场点,若航班k没有飞越进场点,则清空航班k的分配跑道信息ri,进场路径信息pi,将航班从已分配跑道航班集合中剔除,并加入到未分配跑道航班集合转至步骤(8);若航班k已经飞越进场点,转至步骤(8);
若没有找到可剔除航班,则判断航班i是否已经飞越进场点,若航班i已经飞越进场点,则提高航班i的自身优先级转至步骤(8);若航班i没有飞越进场点,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7);
(9122)若航班i的着陆时间偏差Δti=ui-ti满足延误上限要求则将航班i从未排序航班集合剔除,并加入已排序航班集合转至步骤(6);
(92)冻结航班的延误判定;若航班i为冻结航班di=1且着陆时间偏差Δti>0,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7);
(93)挖掘可用空闲时隙;若航班i为非冻结航班di=0且着陆时间偏差Δti≤0,则判断航班i是否满足当前航班可加速条件;
当前航班可加速条件如下:
1)航班i后续航班压力大,Ei>Emin,Ei为根据航班初次分配的结果统计航班i着陆跑道ri上后续排序区间W内未排序航班的数量;Emin表示当前航班可加速的最低后续航班压力,可以自行设置,当前设置为Emin=20架次;
2)航班i优先级小于最低优先级标准(含),Li≤Lmin;
若航班i满足上述可加速的条件1)和2),则对航班i进行加速处理,更新航班i的优先级Li、加速处理信息Aj,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6);若航班i不满足加速条件,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6);
(94)冻结航班处理;若航班i为冻结航班且着陆时间偏差Δti≤0,则计算航班i的可行着陆时间ui,将航班i从未排序航班集合剔除,并加入已排序航班集合更新航班i的综合优先级加速处理信息Aj,转至步骤(6)。
本发明的有益效果为:本发明所述优化排序方法能够兼顾终端区内进场航班的延误、进场飞行时间、以及各跑道流量均衡状况,并进行优化,计算效率高且计算结果稳定,能够满足实际应用需求。
附图说明
图1是本发明的多跑道机场航班进场排序方法流程图。
图2是本发明的多跑道机场航班进场排序方法运行结果示意图。
图3是本发明所采用的航班飞行计划样本数据示意图。
图4是本发明基于管制习惯的排序方法运行结果示意图。
图5是本发明基于最小航线飞行时间的排序方法运行结果示意图。
图6是本发明基于最小延误的排序方法运行结果示意图。
图7是本发明的不同排序方法下航班平均延误对比示意图。
图8是本发明的不同排序方法下航班平均进场飞行时间结果对比示意图。
图9是本发明的不同排序方法下跑道流量结果对比示意图。
具体实施方式
如图1所示,一种多跑道机场进场航班多效能优化排序方法,包括如下步骤:(1)搜集数据信息。
搜集目标机场的跑道和进场航线数据,航班计划数据,航班预计着陆时间数据;进场航班每条进场航线的进场时间;假设一条跑道和一个进场点固定一条进场航线。
(2)设置进场航班排序时间区间W。
按照扩大终端的概念,当前设置为W=90min。
(3)航班数据筛选与集合划分。
确定本场在排序时间区间W内的所有进场航班集合F,所有进场航班数量为n;根据航班的跑道分配情况,将所有进场航班集合划分为已分配跑道航班集合和未分配跑道航班集合两部分,其中未分配跑道航班集合包含所有的未进入进场点的进场航班;根据航班是否已经经过排序,将所有进场航班集合划分为已排序进场航班集合和未排序进场航班集合两部分,已排序航班表示航班的可行着陆时间ui和跑道已经确定,ui表示航班经过排序后得出的航班可行的着陆时间。
(4)航班进场信息整理。
根据航班当前的位置信息确定航班i,i∈F的可用进场路径集合和可用跑道集合以及航班i所有可能的进场飞行时间集合
(5)跑道初次分配。
对未分配跑道的航班找出航班i的最小进场路径飞行时间min(FTf),所对应的初始路径和初始跑道并计算得到航班i的预计着陆时间跑道初次分配用于确定航班初始排序次序,并非正式分配跑道,因此经跑道初次分配处理的航班仍属于跑道未分配航班,后续需经跑道再次分配处理。
为保证航班安全着陆和排序结果稳定性,定义距离着陆时间小于TD(含)的进场航班为冻结航班,对应的冻结航班集合为D;其中,TD为冻结航班划分时间临界值,可根据需要自行设置,当前设置为15min;该类航班所分配的跑道、着陆时间和顺序不再改变,则距离着陆时间大于TD以上的进场航班为非冻结航班,对应的非冻结航班集合为D',F=D∪D',另di为0-1变量,di=1表示航班为冻结航班,di=0表示航班为非冻结航班。
(6)确定当前排序航班。
根据跑道初次分配所得航班预计着陆时间找出未排序进场航班集合中预计着陆时间最早的一架进场航班i,
(7)跑道再次分配。
若进场航班i已经分配跑道r',则不再重新分配跑道,ri=r',确定航班的进场航线p'和预计着陆时间t',则pi=p',ti=t';若进场航班i还未分配跑道则结合多效能跑道分配评价函数(式(1)),计算进场航班i的多效能跑道评价系数为航班i分配跑道ri和进场航线pi,计算预计着陆时间ti;并将航班i从未分配跑道航班集合中剔除,加入到已分配跑道航班集合
多效能跑道分配评价函数考虑进场飞行时间,跑道流量均衡性和航班延误三个方面,多效能跑道评价系数的计算如下:
式中,为航班i在跑道r上进场飞行时间评价系数;为航班i可用的最短飞行时间对应航线的影响系数:μ1和k为调节参数,用于调整的增长幅度;Tmin为机场最小的进场航线飞行时间,Ti min为航班i可行的最小进场航线飞行时间;Tmax为机场最长的进场航线飞行时间,Ti max为航班i可行的最长进场航线飞行时间;为航班i在跑道r着陆的进场航线飞行时间偏差系数;μ2为调节参数,用于调整的增长幅度;为航班i在跑道r着陆对应进场航线飞行时间的偏差比;为航班i所有可行进场航线飞行时间偏离最小进场航线飞行时间的标准差;为航班i所有可用跑道数量;
为航班i在跑道r上的流量评估系数;μ3>0用于衡量航班i所有可用跑道的流量均衡性;μ3为调节参数,用于的增长幅度;C为机场容量;为航班i可用跑道集合Ri的流量标准差, 为航班i在跑道r的流量偏差系数,μ4为调节参数,用于调整的增长幅度;为航班i在跑道r的流量偏差: 为跑道r的当前已排序航班数量;为航班i在所有可用跑道上的平均流量;为航班i在跑道r着陆的准点系数;xi∈[1,2]为延误重视度系数;为航班i的自身优先级,其取值范围为{1,2,3,4},该值由航班自身类型确定,1表示最低优先级,一般为货运航班,2表示普通优先级,一般为普通客运航班,3表示较高优先级,4表示最高优先级,也可以根据需求设置和调整;ω1为调节参数,用于调整的增长幅度;航班i在跑道r着陆的准点率;为航班i在跑道r着陆所需的延误占比;为航班i在跑道r上的着陆延误;为航班i在跑道r着陆所能承受的最大正延误时间;M为普通优先级航班在终端区内的最大可接收延误阀值,可自行设置,当前已经起飞的航班设置为45分钟,未起飞的航班设置为2小时;ω2,ω3是的权值;为航班i的着陆趋势优先级,则:且 为航班i在跑道r上的预计着陆时间;tsys为***当前时间;tB为终端区外界时间范围对应于进场航班排序时间区间W,tB=W。
为了把统一到一个参考系,需要进行归一化处理;另 分别为hk所能取到的最小值和最大值。则
其中,
g3'=(1-ω4)·(g1-1)+ω4,ω4∈[0.5,1),
(8)确定可行着陆时间。
对航班i进行冲突探测与解脱,确定航班i在跑道ri的上可行着陆时间ui;此处的冲突探测包括跑道和终端区飞行过程中的冲突。
(9)优化航班延误。
对比航班i预计着陆时间ti与可行着陆时间ui,并计算航班i的着陆时间偏差Δti=ui-ti,将其作为航班的延误时间,并对其进一步优化。
Step 9.1:优化非冻结航班的延误。
若航班i为非冻结航班di=0且着陆时间偏差Δti>0。
Step 9.1.1:前序航班加速。
计算航班i综合优先级Li(式(2)),判断跑道ri上先于航班i前一架进场航班j是否满足前序航班加速条件。
航班i在跑道r上的综合优先级定义如下
其中,为航班i的自身优先级,定义参考步骤(7)为航班i的延误动态优先级,则且 定义为航班延误占比(参见(7)计算方法),可知当前排序航班i在跑道r上的着陆延误, 为航班i的着陆趋势优先级,定义参见步骤(7);设置了航班的可加速最低优先级Lmin,根据需求自行设置,当前设置为Lmin=2;
前序航班加速条件如下:
1)前序航班j综合优先级小于航班i综合优先级,Li>Lj;
2)航班i综合优先级大于最低可加速优先级Lmin,Li>Lmin;
3)前序航班j为非冻结航班,dj=0;
4)前序航班j未经过延误处理,Dj=0,其中,Dj表示航班j是否经过延误处理,0表示没有经过延误处理,1则表示经过延误处理;
5)前序航班j未经过加速处理,Aj=0,其中,Aj表示航班j是否经过加速处理,0表示没有经过加速处理,1则表示经过加速处理;
6)前序航班j与航班i之间存在冲突;
若前序航班j满足上述加速条件1)-6),则对航班i和j进行加速处理和冲突探测解脱,重新计算航班i和j的可行着陆时间、综合优先级和加速处理信息ui,uj,Li,Lj,Ai,Aj,以及着陆时间偏差Δti;若不满足以上加速条件,则直接转至Step 9.1.2。
Step 9.1.2:判断是否满足延误上限要求。判断航班i的着陆时间偏差Δti是否满足航班i延误上限要求,对应于航班i的可接受延误阈值
Step 9.1.2.1:查找可替换航班。若航班i的着陆时间偏差Δti不满足延误上限 则从跑道ri上所有前序航班中找出航班k,判断航班k是否满足以下前序航班可剔除条件:
1)前序航班k优先级小于航班i优先级,Li>Lk;
2)前序航班k为非冻结航班,dk=0;
3)航班i和前序航班k之间存在冲突;
4)航班i和前序航班k的使用不同的进场航线,pi≠pk;
5)前序航班k已参与剔除的次数Nk满足航班避让次数阈值K要求,Nk<K,为防止某架航班被多次剔除;其中,航班避让次数阈值K,可根据需求自行设置,当前设置为K=2;
6)为避让航班i而被剔除的前序航班数量PNi应满足航班剔除次数阈值PK要求,PNi<PK,为防止航班i的排序结果过分偏离其初始次序;其中,航班剔除次数阈值PK,可根据需求自行设置,当前设置为PK=2;
若存在前序航班k满足以上前序航班可剔除条件1)-6),则将前序航班k剔除,更新uk之后的已排序航班信息,包括可行着陆时间,着陆时间偏差,更新航班避让次数信息Nk,清空航班k可行着陆时间uk,将航班k从已排序航班集合中剔除,加入到未排序航班集合然后判断剔除航班k是否已经飞越进场点,若航班k没有飞越进场点,则清空航班k的分配跑道信息ri,进场路径信息pi,将航班从已分配跑道航班集合中剔除,并加入到未分配跑道航班集合转至步骤(8);若航班k已经飞越进场点,转至步骤(8);
若没有找到可剔除航班,则判断航班i是否已经飞越进场点,若航班i已经飞越进场点,则提高航班i的自身优先级转至步骤(8);若航班i没有飞越进场点,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7)。
Step 9.1.2.2:若航班i的着陆时间偏差Δti=ui-ti满足延误上限要求则将航班i从未排序航班集合剔除,并加入已排序航班集合转至步骤(6)。
Step 9.2:冻结航班的延误判定;若航班i为冻结航班di=1且着陆时间偏差Δti>0,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7)。
Step 9.3:挖掘可用空闲时隙;若航班i为非冻结航班di=0且着陆时间偏差Δti≤0,则判断航班i是否满足当前航班可加速条件;
当前航班可加速条件如下:
1)航班i后续航班压力大,Ei>Emin,Ei为根据航班初次分配的结果统计航班i着陆跑道ri上后续排序区间W内未排序航班的数量;Emin表示当前航班可加速的最低后续航班压力,可以自行设置,当前设置为Emin=20架次;
2)航班i优先级小于最低优先级标准(含),Li≤Lmin;
若航班i满足可加速的条件1)和2),则对航班i进行加速处理,更新航班i的优先级Li、加速处理信息Aj,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6)。若航班i不满足加速条件,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6)。
Step 9.4:冻结航班处理;若航班i为冻结航班且着陆时间偏差Δti≤0,则计算航班i的可行着陆时间ui,将航班i从未排序航班集合剔除,并加入已排序航班集合更新航班i的综合优先级加速处理信息Aj,转至步骤(6)。
根据本发明设计思路,开发和实现了多跑道机场航班进场排序原型***,鉴于***截图清晰度问题,手绘原型***其运行界面示意图如图2所示,图中包含三根时间轴,代表机场三条不同跑道的进离场航班队列,每根时间轴左侧为进场降落航班队列,右侧为离场航班队列;图中每个方框代表一架航班,方框内包含航班号以及该航班经排序后的延误信息,时间轴顶部置灰航班为冻结航班,其余为非冻结航班;选取北京首都机场某典型日不同时刻的终端区航班运行数据进行实例验证,其中选取8:30、10:30、14:30、17:30、22:30五个时刻,对应航班架次为22,46,29,41,56,原型***运算时间分别为1秒,2秒,1秒,2秒,2秒。图3为航班计划数据和优化后的结果数据示例,从两图可以看出通过本发明方法,航班进场排序优化效果明显。
为说明本发明方法的优化效果,在***中还将所提优化排序方法(曲线一)与以下三种不同排序方法的结果进行了对比,分别为基于管制习惯的排序方法(曲线二)、基于最小航线飞行时间的排序方法(曲线三)和基于最小延误的排序方法(曲线四)。其中,基于管制习惯的排序方法是基于历史数据统计分析,而基于最小进场飞行时间的排序方法和基于最小延误的排序方法则基于贪婪算法,本发明所提方法简称为优化排序方法。三种排序方法的运行结果图4、图5和图6所示。
下面分别从平均进场飞行时间、平均延误时间、跑道流量均衡三个方面对上述四种排序方法进行对比。其中跑道流量均衡性为各跑道的流量均衡情况。
图7为不同排序方法下航班平均延误结果对比。由图可知,从整体延误水平上看,优化效果由从优到劣依次为基于最小延误排序方法、优化排序方法、基于管制习惯的排序方法、基于最小进场飞行时间排序方法,其中所提优化排序方法与基于最小延误排序方法的评价结果贴合度较高,在8:30、14:30两个时刻表现优于最小延误排序方法。此外,随着航班量的变化,以基于最小飞行延误排序方法结果为基准,优化排序方法结果偏差较小,优化效果较为稳定,明显优于其他两种排序方法,由此说明优化排序方法对航班延误进行了有效控制。
图8为不同排序方法下航班平均进场飞行时间结果对比,由图可知,从航班平均进场飞行时间上看,优化效果由从优到劣依次为基于最小进场飞行时间排序方法、优化排序方法、基于管制习惯的排序方法、基于最小延误排序方法。其中,以基于最小进场飞行时间排序方法结果为基准,优化排序方法结果偏差最小,优化效果稳定,而基于最小延误排序方法与基于管制习惯的排序方法结果偏差较大且波动明显。由此说明优化排序方法对航班进场飞行时间进行了有效控制。
图9为同排序方法下跑道流量结果对比,从跑道均衡性上来看,优化效果从优到劣依次为优化排序方法、基于最小延误排序方法、基于最小进场飞行时间排序方法、基于管制***稳;基于管制习惯的排序方法则导致跑道流量存在明显波动;由此说明优化排序方法对各跑道流量均衡状况进行了有效控制。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。
Claims (2)
1.一种多跑道机场进场航班多效能优化排序方法,其特征在于,包括如下步骤:
(1)搜集数据信息;搜集目标机场的跑道和进场航线数据,航班计划数据,航班预计着陆时间数据;进场航班每条进场航线的进场时间;假设一条跑道和一个进场点固定一条进场航线;
(2)设置进场航班排序时间区间W;按照扩大终端的概念,当前设置为W=90min;
(3)航班数据筛选与集合划分;确定本场在排序时间区间W内的所有进场航班集合F,所有进场航班数量为n;根据航班的跑道分配情况,将所有进场航班集合划分为已分配跑道航班集合和未分配跑道航班集合两部分,其中未分配跑道航班集合包含所有的未进入进场点的进场航班;根据航班是否已经经过排序,将所有进场航班集合划分为已排序进场航班集合和未排序进场航班集合两部分,已排序航班表示航班的可行着陆时间ui和跑道已经确定,ui表示航班经过排序后得出的航班可行的着陆时间;
(4)航班进场信息整理;根据航班当前的位置信息确定航班i,i∈F的可用进场路径集合和可用跑道集合以及航班i所有可能的进场飞行时间集合
(5)跑道初次分配;对未分配跑道的航班找出航班i的最小进场路径飞行时间min(FTf)所对应的初始路径和初始跑道ri 0,并计算得到航班i的预计着陆时间跑道初次分配用于确定航班初始排序次序,并非正式分配跑道,因此经跑道初次分配处理的航班仍属于跑道未分配航班,后续需经跑道再次分配处理;为保证航班安全着陆和排序结果稳定性,定义距离着陆时间小于或者等于TD的进场航班为冻结航班,对应的冻结航班集合为D;其中,TD为冻结航班划分时间临界值,可根据需要自行设置,当前设置为15min;该类航班所分配的跑道、着陆时间和顺序不再改变,则距离着陆时间大于TD以上的进场航班为非冻结航班,对应的非冻结航班集合为D',F=D∪D',另di为0-1变量,di=1表示航班为冻结航班,di=0表示航班为非冻结航班;
(6)确定当前排序航班;根据跑道初次分配所得航班预计着陆时间找出未排序进场航班集合中预计着陆时间最早的一架进场航班i,
(7)跑道再次分配;具体如下:若进场航班i已经分配跑道r',则不再重新分配跑道,ri=r',确定航班的进场航线p'和预计着陆时间t',则pi=p',ti=t';若进场航班i还未分配跑道则结合多效能跑道分配评价函数,计算进场航班i的多效能跑道评价系数为航班i分配跑道ri和进场航线pi,计算预计着陆时间ti;并将航班i从未分配跑道航班集合中剔除,加入到已分配跑道航班集合
多效能跑道分配评价函数考虑进场飞行时间,跑道流量均衡性和航班延误三个方面,多效能跑道评价系数的计算如下
式中,为航班i在跑道r上进场飞行时间评价系数;为航班i可用的最短飞行时间对应航线的影响系数:μ1>0,k≥0,μ1和k为调节参数,用于调整的增长幅度;Tmin为机场最小的进场航线飞行时间,Ti min为航班i可行的最小进场航线飞行时间;Tmax为机场最长的进场航线飞行时间,Ti max为航班i可行的最长进场航线飞行时间;为航班i在跑道r着陆的进场航线飞行时间偏差系数;μ2为调节参数,用于调整的增长幅度;ΔTi r≥0为航班i在跑道r着陆对应进场航线飞行时间Ti r的偏差比;为航班i所有可行进场航线飞行时间偏离最小进场航线飞行时间的标准差;为航班i所有可用跑道数量;
为航班i在跑道r上的流量评估系数;μ3>0用于衡量航班i所有可用跑道的流量均衡性;μ3为调节参数,用于的增长幅度;C为机场容量;为航班i可用跑道集合Ri的流量标准差, 为航班i在跑道r的流量偏差系数,μ4>0;μ4为调节参数,用于调整的增长幅度;为航班i在跑道r的流量偏差: 为跑道r的当前已排序航班数量;为航班i在所有可用跑道上的平均流量;为航班i在跑道r着陆的准点系数;xi∈[1,2]为延误重视度系数;为航班i的自身优先级,其取值范围为{1,2,3,4},该值由航班自身类型确定,1表示最低优先级,一般为货运航班,2表示普通优先级,一般为普通客运航班,3表示较高优先级,4表示最高优先级,也可以根据需求设置和调整;ω1为调节参数,用于调整的增长幅度;航班i在跑道r着陆的准点率;为航班i在跑道r着陆所需的延误占比;为航班i在跑道r上的着陆延误;ω2+ω3=1;→ω2=0.5,ω3=0.5;为航班i在跑道r着陆所能承受的最大正延误时间;M为普通优先级航班在终端区内的最大可接收延误阀值,可自行设置,当前已经起飞的航班设置为45分钟,未起飞的航班设置为2小时;ω2,ω3是的权值;为航班i的着陆趋势优先级,则: 且IF:IF: 为航班i在跑道r上的预计着陆时间;tsys为***当前时间;tB为终端区外界时间范围对应于进场航班排序时间区间W,tB=W;
为了把统一到一个参考系,需要进行归一化处理;另 k∈1,2,3,分别为hk所能取到的最小值和最大值,则:
其中,g3'=(1-ω4)·(g1-1)+ω4,ω4∈[0.5,1),
(8)确定可行着陆时间;对航班i进行冲突探测与解脱,确定航班i在跑道ri的上可行着陆时间ui;此处的冲突探测包括跑道和终端区飞行过程中的冲突;
(9)优化航班延误;对比航班i预计着陆时间ti与可行着陆时间ui,计算航班i的着陆时间偏差Δti=ui-ti,将其作为航班的延误时间,并对其进一步优化。
2.如权利要求1所述的多跑道机场进场航班多效能优化排序方法,其特征在于,步骤(9)中,优化航班延误方法,具体为:
(91)优化非冻结航班的延误;若航班i为非冻结航班di=0且着陆时间偏差Δti>0;
(911)前序航班加速;计算航班i综合优先级Li,判断跑道ri上先于航班i前一架进场航班j是否满足前序航班加速条件;
航班i在跑道r上的综合优先级定义如下
其中,为航班i的自身优先级,定义参考步骤(7);为航班i的延误动态优先级,则 且IF: 定义为航班延误占比,可知当前排序航班i在跑道r上的着陆延误, 为航班i的着陆趋势优先级,定义参见步骤(7);设置了航班的可加速最低优先级Lmin,根据需求自行设置,当前设置为Lmin=2;
前序航班加速条件如下:
1)前序航班j综合优先级小于航班i综合优先级,Li>Lj;
2)航班i综合优先级大于最低可加速优先级Lmin,Li>Lmin;
3)前序航班j为非冻结航班,dj=0;
4)前序航班j未经过延误处理,Dj=0,其中,Dj表示航班j是否经过延误处理,0表示没有经过延误处理,1则表示经过延误处理;
5)前序航班j未经过加速处理,Aj=0,其中,Aj表示航班j是否经过加速处理,0表示没有经过加速处理,1则表示经过加速处理;
6)前序航班j与航班i之间存在冲突;
若前序航班j满足上述加速条件1)-6),则对航班i和j进行加速处理和冲突探测解脱,重新计算航班i和j的可行着陆时间、综合优先级和加速处理信息ui,uj,Li,Lj,Ai,Aj,以及着陆时间偏差Δti;若不满足以上加速条件,则直接转至步骤(912);
(912)判断是否满足延误上限要求;判断航班i的着陆时间偏差Δti是否满足航班i延误上限要求,对应于航班i的可接受延误阈值
(9121)查找可替换航班;若航班i的着陆时间偏差Δti不满足延误上限 则从跑道ri上所有前序航班中找出航班k,判断航班k是否满足以下前序航班可剔除条件:
1)前序航班k优先级小于航班i优先级,Li>Lk;
2)前序航班k为非冻结航班,dk=0;
3)航班i和前序航班k之间存在冲突;
4)航班i和前序航班k的使用不同的进场航线,pi≠pk;
5)前序航班k已参与剔除的次数Nk满足航班避让次数阈值K要求,Nk<K,为防止某架航班被多次剔除;其中,航班避让次数阈值K,可根据需求自行设置,当前设置为K=2;
6)为避让航班i而被剔除的前序航班数量PNi应满足航班剔除次数阈值PK要求,PNi<PK,为防止航班i的排序结果过分偏离其初始次序;其中,航班剔除次数阈值PK,可根据需求自行设置,当前设置为PK=2;
若存在前序航班k满足上述前序航班可剔除条件1)-6),则将前序航班k剔除,更新uk之后的已排序航班信息,包括可行着陆时间,着陆时间偏差,更新航班避让次数信息Nk,清空航班k可行着陆时间uk,将航班k从已排序航班集合中剔除,加入到未排序航班集合然后判断剔除航班k是否已经飞越进场点,若航班k没有飞越进场点,则清空航班k的分配跑道信息ri,进场路径信息pi,将航班从已分配跑道航班集合中剔除,并加入到未分配跑道航班集合转至步骤(8);若航班k已经飞越进场点,转至步骤(8);
若没有找到可剔除航班,则判断航班i是否已经飞越进场点,若航班i已经飞越进场点,则提高航班i的自身优先级转至步骤(8);若航班i没有飞越进场点,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7);
(9122)若航班i的着陆时间偏差Δti=ui-ti满足延误上限要求则将航班i从未排序航班集合剔除,并加入已排序航班集合转至步骤(6);
(92)冻结航班的延误判定;若航班i为冻结航班di=1且着陆时间偏差Δti>0,则将跑道ri从航班i的可用跑道集合Ri中剔除,并更新航班i的可用路径集合Pi,清空航班i的分配跑道信息ri,进场路径信息pi,将航班i从已分配跑道航班集合中剔除,加入到未分配跑道航班集合转至步骤(7);
(93)挖掘可用空闲时隙;若航班i为非冻结航班di=0且着陆时间偏差Δti≤0,则判断航班i是否满足当前航班可加速条件;
当前航班可加速条件如下:
1)航班i后续航班压力大,Ei>Emin,Ei为根据航班初次分配的结果统计航班i着陆跑道ri上后续排序区间W内未排序航班的数量;Emin表示当前航班可加速的最低后续航班压力,可以自行设置,当前设置为Emin=20架次;
2)航班i优先级小于或者等于最低优先级标准,Li≤Lmin;
若航班i满足可加速的条件1)和2),则对航班i进行加速处理,更新航班i的优先级Li、加速处理信息Aj,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6);若航班i不满足加速条件,将航班i从未排序航班集合中剔除,加入到已排序航班集合转至步骤(6);
(94)冻结航班处理;若航班i为冻结航班且着陆时间偏差Δti≤0,则计算航班i的可行着陆时间ui,将航班i从未排序航班集合剔除,并加入已排序航班集合更新航班i的综合优先级加速处理信息Aj,转至步骤(6)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329224.0A CN107016881B (zh) | 2017-05-11 | 2017-05-11 | 一种多跑道机场进场航班多效能优化排序方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329224.0A CN107016881B (zh) | 2017-05-11 | 2017-05-11 | 一种多跑道机场进场航班多效能优化排序方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016881A CN107016881A (zh) | 2017-08-04 |
CN107016881B true CN107016881B (zh) | 2019-05-31 |
Family
ID=59450159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710329224.0A Active CN107016881B (zh) | 2017-05-11 | 2017-05-11 | 一种多跑道机场进场航班多效能优化排序方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016881B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107591034A (zh) * | 2017-08-09 | 2018-01-16 | 南京莱斯信息技术股份有限公司 | 一种进港航班排序的实现方法 |
CN107944625B (zh) * | 2017-11-23 | 2021-09-07 | 南京航空航天大学 | 基于历史运行数据驱动的单机场航班换季时隙优化方法 |
CN108230757B (zh) * | 2018-01-29 | 2020-12-15 | 上海空地信息技术有限公司 | 一种基于水平航迹的终端区进场效率测度方法 |
CN110349444B (zh) * | 2018-04-06 | 2023-11-10 | 杭州坚果壳科技开发有限公司 | 基于大数据的空中交通流量管理方法 |
CN109583627B (zh) * | 2018-10-31 | 2020-09-29 | 北京航空航天大学 | 飞机着陆排队优化方法及装置 |
CN110322161A (zh) * | 2019-07-10 | 2019-10-11 | 中国民航信息网络股份有限公司 | 航班生效批次的生效调整方法及装置 |
CN110706520B (zh) * | 2019-10-31 | 2021-08-31 | 中国电子科技集团公司第二十八研究所 | 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法 |
CN110909946B (zh) * | 2019-11-28 | 2020-11-24 | 北京航空航天大学 | 一种基于公路换乘的航班计划优化方法 |
CN111899571B (zh) * | 2020-07-29 | 2021-11-02 | 成都民航空管科技发展有限公司 | 基于aman***的自动调整航班进近程序的方法及*** |
CN112447067A (zh) * | 2020-10-09 | 2021-03-05 | 广州市中南民航空管通信网络科技有限公司 | 多机场多跑道进港航班的动态调度方法、终端、存储装置 |
CN112447068A (zh) * | 2020-10-09 | 2021-03-05 | 广州市中南民航空管通信网络科技有限公司 | 多跑道机场的跑道分配方法、终端、存储装置 |
CN113643571B (zh) * | 2021-10-18 | 2022-02-08 | 中国电子科技集团公司第二十八研究所 | 一种基于航班正常性目标的空域网络优化方法 |
CN115310731B (zh) * | 2022-10-12 | 2022-12-30 | 珠海翔翼航空技术有限公司 | 航班进离场调度方法、***、装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1987002493A1 (en) * | 1985-10-17 | 1987-04-23 | Toyo Communication Equipment Co., Ltd. | Alarm restraining system in an apparatus for preventing collision of aircraft |
JP3913623B2 (ja) * | 2002-07-05 | 2007-05-09 | 株式会社エヌ・ティ・ティ・データ | 航空交通管制制御装置 |
US20070010921A1 (en) * | 2005-07-05 | 2007-01-11 | Honeywell International Inc. | Method, apparatus, and database products for automated runway selection |
CN101465064B (zh) * | 2009-01-15 | 2011-03-30 | 北京航空航天大学 | 终端区飞行冲突解脱方法和*** |
CN105679102B (zh) * | 2016-03-03 | 2018-03-27 | 南京航空航天大学 | 一种全国飞行流量时空分布预测推演***及方法 |
CN105741614B (zh) * | 2016-03-14 | 2018-03-13 | 中国民航大学 | 一种能够提高繁忙机场跑道飞机降落利用率的***化方法 |
-
2017
- 2017-05-11 CN CN201710329224.0A patent/CN107016881B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107016881A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016881B (zh) | 一种多跑道机场进场航班多效能优化排序方法 | |
CN107591034A (zh) | 一种进港航班排序的实现方法 | |
CN110602633B (zh) | 一种面向爆发性流量的移动边缘计算无人机群辅助通信方法 | |
CN103413462B (zh) | 一种综合考虑空域拥堵和航班延误的空中交通网络流量优化方法 | |
CN111160770A (zh) | 一种进场航班动态协同排序方法 | |
CN112396871B (zh) | 基于航迹预测的进场延误分配吸纳方法 | |
CN110428665B (zh) | 一种航路与机场时隙协同分配的随机双层规划方法 | |
CN109741638B (zh) | 一种进离场管理***一体化协同运行方法 | |
CN105825717B (zh) | 一种基于不确定到达时刻的空域时隙资源优化分配方法 | |
CN116580601A (zh) | 基于空域约束和进场流量的进港排序方法 | |
WO2022213756A1 (zh) | 一种航路时隙与高度层协同优化分配方法 | |
CN105813218A (zh) | 无线自组织网络的业务调度方法 | |
CN116307542A (zh) | 基于多计量点约束的进港排序方法及电子终端 | |
CN109948844A (zh) | 一种停机位分配鲁棒性的优化方法、装置、设备及介质 | |
CN110706520B (zh) | 一种基于概率容量的航路与机场时隙鲁棒分配双层规划方法 | |
CN116029534B (zh) | 一种机场停机位的分配方法、装置、电子设备及存储介质 | |
CN102890876B (zh) | 基于流量比率的区域尾随间隔限制值的计算方法 | |
EP1616313B1 (en) | Method for sequencing landing aircrafts | |
Isaacson et al. | Knowledge-based runway assignment for arrival aircraft in the terminal area | |
CN110751857A (zh) | 一种机会约束空中交通流量管理方法 | |
Soomer et al. | Fairness in the aircraft landing problem | |
CN108960649A (zh) | 一种外卖配送调度方法 | |
CN112201082B (zh) | 一种航班放行时隙的n点置换方法 | |
CN113986531A (zh) | 一种任务调度以及资源分配方法 | |
CN114550505A (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 |