CN113483773B - 一种对流天气下的航路使用方法 - Google Patents
一种对流天气下的航路使用方法 Download PDFInfo
- Publication number
- CN113483773B CN113483773B CN202110657978.5A CN202110657978A CN113483773B CN 113483773 B CN113483773 B CN 113483773B CN 202110657978 A CN202110657978 A CN 202110657978A CN 113483773 B CN113483773 B CN 113483773B
- Authority
- CN
- China
- Prior art keywords
- flight
- weather
- path
- navigation
- route
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000010006 flight Effects 0.000 claims abstract description 10
- 238000012216 screening Methods 0.000 claims abstract description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/34—Route searching; Route guidance
- G01C21/3453—Special cost functions, i.e. other than distance or default speed limit of road segments
- G01C21/3492—Special cost functions, i.e. other than distance or default speed limit of road segments employing speed data or traffic data, e.g. real-time or historical
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种对流天气下的航路使用方法,包括如下步骤:首先指定航路,其次获取并筛选历史航班数据,之后利用历史当时的气象条件和改航情况生成天气避让场,然后根据当前天气条件计算对流天气的边界,参考场景分类,最后生成航路的使用策略。本发明使得空管人员能够在对流天气条件下有效利用航路,提高航路的飞行安全水平,使得航班尽可能按时飞行,减少空中延误和地面堵塞,降低对流天气下的经济损失。
Description
技术领域
本发明涉及航路规划设计,特别涉及一种对流天气下的航路使用方法。
背景技术
随着最近几年机场的快速建设和航班量的增长,对流天气下航路的使用和通行安全的重要性逐步增加。对流天气条件下航路使用策略的生成,可以提高航路利用效率,在保障航班安全性的同时得到最大的经济效益,实现安全性和经济性之间的最佳平衡。
目前对流天气下航路是否使用,取决于管制员和飞行员的历史经验,多数情况下会关闭航路,致使航班不能按计划飞行,造成大量空中延误和地面延误,造成严重的经济损失。
发明内容
发明目的:本发明目的是提供一种对流天下的航路使用策略。
技术方案:本发明提供一种对流天气下的航路使用方法,包括如下步骤:
步骤1:指定航路(或航段);
步骤2:获取并筛选历史航班数据;
步骤3:根据历史数据的雷达轨迹和飞行计划路径判断是否改航,结合当时的气象情况,构建天气避让场;
步骤4:根据当前的气象条件,对照天气避让场,确定对流天气避让范围;
步骤5:给出避让原则;
进一步地,所述步骤2方法:获取航班历史数据,筛选在航路T范围内飞行的航班历史数据,包括历史雷达轨迹数据、飞行计划数据和气象指标数据,所述历史雷达轨迹数据包括起飞时间、飞行时间、着陆时间、呼号、地面速度、飞行高度、经纬度;所述飞行计划数据包括航班号、起飞时间、着陆时间、起飞机场、着陆机场、飞行计划路径;所述气象指标数据包括CR、ET和VIL。
进一步地,所述步骤3的方法:判断是否改航的依据是计算实际雷达轨迹和飞行计划路径之间的距离,进行航班聚类,寻找聚类中心线,聚类中心线到飞行计划路径的距离即为阈值,以此阈值来判断航班是否改航,针对改航航班,统计气象情况,根据历史轨迹改航的航班数,计算改航概率,改航概率=改航航班数/在此气象条件下的航班总数,生成以VIL为横轴,横坐标为0-0.05、0.05-0.14、0.14-0.7、0.7-3.5、3.5-6.9、6.9-12、12+,单位为kg/m3;deltaZ为纵轴,纵坐标为-10km到10km,每隔2.5km为一个间隔,每一格代表对应改航概率的天气避让概率。
进一步地,所述步骤4方法:当VIL介于3.5kg/m3和6.9kg/m3,deltaZ介于-7.5km和0时,设改航概率为0.4,由黄色表示;当VIL大于6.9kg/m3,deltaZ介于-10和-7.5km时,改航概率设为0.7,由红色表示,因此规定对流天气避让范围由外及内分别为黄色和红色,分别代表改航概率为0.4-0.7和0.7-1,确定避让对流范围的具体算法是:首先使用Canny边缘检测器检测出天气避让场的边界,然后用DBSCAN算法将边缘像素进行聚类,最后将每一种聚类结果包含在凸包中,形成避让范围,如果两个凸包的距离小于8km,则合并为一个。
进一步地,所述步骤5方法:避让原则由对流天气的位置以及个数决定,令沿飞行计划路径方向最先与黄色天气避让区相交的点为飞入交汇点Ri,另一交汇点为飞出交汇点Ci,选取距离Ri和Ci最近的计划路径上的点分别作为改航起始点Si和改航结束点Ei(i=1,2,3,......,为飞行计划路径上依次遇到的对流天气区域编号),当飞行计划路径上只有单个的对流天气时,策略一:若对流天气存在于飞行计划路径中间,则先沿飞行计划路径飞行,然后找到改航点ri,将ri和改航起始点Si以及改航结束点Ei相连,作为改航的航段,以此进行改航,然后继续沿着飞行计划路径飞行;策略二:若该对流天气存在于飞行计划路径起点,则一开始就进行改航,找到改航点ri,并做平行于计划路径的航段qiri,连接ri和改航结束点Ei,然后沿着飞行计划路径继续飞行;策略三:若对流天气存在于飞行计划路径终点,则先沿飞行计划路径飞行,找到改航点ri,连接ri和改航起始点Si,再以ri为起点作平行于飞行计划路径的航段;当飞行计划路径上存在多个对流天气时,策略四:计算前一个飞出交汇点Ci-1和后一个改航起始点Si之间的距离L,若L<7.4km,则进行连续改航,首先找到每个对流天气所对应的改航点,连接Si和距离最近的改航点,再连接两个改航点;策略五:若L≥7.4km,则进行组合改航,策略六:若改航后的路径总距离≥1.5飞行计划路径总距离,则直接改航到其他航路。
有益效果:本发明通过该对流天气下的航路使用方法,使得空管人员能够在对流天气条件下有效利用航路,保障航路安全合理地使用,充分利用航路资源。能够保证对流天气下的航班飞行安全,提高对流天气条件下航路飞行的安全水平。能够减少空中延误和地面堵塞,使得航班尽可能按时飞行,减少对流天气下的经济损失。
附图说明
图1为本发明的方法流程图;
图2为航路图;
图3为VIL图;
图4为ET图;
图5为天气避让图;
图6为改航点的确定;
图7为单个对流天气无影响;
图8为单个对流天气存在于飞行计划路径中间;
图9为单个对流天气存在于飞行计划路径起点;
图10为单个对流天气存在于飞行计划路径终点;
图11为多个对流天气连续改航;
图12为多个对流天气组合改航。
具体实施方式
以下将结合附图,对本发明的技术方案及有益效果进行详细说明。
如图1所示,本发明公开了一种对流天气下的航路使用方法,具体步骤如下:
步骤1:获取航路;选定航路T,如图2所示。
步骤2:获取并筛选历史航班数据
获取2018年8月航班历史数据,筛选在航路T范围内飞行的航班历史数据,包括历史雷达轨迹数据、飞行计划数据和气象指标数据。所述历史雷达轨迹数据包括起飞时间、飞行时间、着陆时间、呼号、地面速度、飞行高度、经纬度;所述飞行计划数据包括航班号、起飞时间、着陆时间、起飞机场、着陆机场、飞行计划路径;所述气象指标数据包括CR、ET和VIL。
步骤3:根据历史数据中的航班的雷达轨迹和飞行计划路径判断是否改航,统计改航航班气象情况,计算改航概率,构建天气避让概率表;
判断是否改航的依据是计算实际雷达轨迹和飞行计划路径之间的距离,进行航班聚类,寻找聚类中心线,聚类中心线到飞行计划路径的距离即为阈值,以此阈值来判断航班是否改航。针对改航航班,统计气象情况,根据历史轨迹改航的航班数,计算改航概率,改航概率=改航航班数/在此气象条件下的航班总数。生成以VIL为横轴,横坐标为0-0.05、0.05-0.14、0.14-0.7、0.7-3.5、3.5-6.9、6.9-12、12+,单位为kg/m3;deltaZ(飞行高度-90%ET高度)为纵轴,纵坐标为-10km到10km,每隔2.5km为一个间隔,每一格代表对应改航概率的天气避让概率。
步骤4:根据当前的气象条件,对照天气避让概率表,确定对流天气避让图;
当VIL介于3.5kg/m3和6.9kg/m3,deltaZ介于-7.5km和0时,由于航班总数据量小,航班改航和未改航的数据量比较接近,具有高度不确定性,因此设改航概率为0.4,由黄色表示;当VIL大于6.9kg/m3,deltaZ介于-10和-7.5km时,改航概率设为0.7,由红色表示。因此规定对流天气避让范围由外及内分别为黄色和红色,分别代表改航概率为0.4-0.7和0.7-1,确定避让对流范围的具体算法是:首先使用Canny边缘检测器检测出天气避让场的边界,然后用DBSCAN算法将边缘像素进行聚类,最后将每一种聚类结果包含在凸包中,形成避让范围。如果两个凸包的距离小于8km,则合并为一个。
步骤5:给出避让原则;
在飞机性能上:转弯半径不能大于90°;航段距离不能小于7.4km。
在改航方法上:首先确定改航起始点和改航结束点,然后根据不同情况实行不同的航路使用策略。
当红色避让区域与黄色避让区域在边界有重合时,可将红色区域继续外扩10km,然后当成黄色区域处理,与原先的黄色区域再一起形成新的避让边界。
具体方法:令沿飞行计划路径方向最先与黄色区域相交的点为飞入交汇点Ri,另一交汇点为飞出交汇点Ci,选取距离Ri和Ci最近的计划路径上的点分别作为改航起始点Si和改航结束点Ei(i=1,2,3,......,为飞行计划路径上依次遇到的对流天气区域编号)。
定义RiCi的中点Zi,RiCi将避让范围分为两部分,将分属两部分的顶点分别放入集合Pr和Pl中,判断方法为:将各项点带入RiCi直线的一般方程,如果大于0则顶点属于Pr,反之属于Pl。在图1中,Pr={P2,P3,P4},Pl={P1,P5,P6},在Pr与Pl中分别搜索距点Zi最远的点(图1中为P4和P5),并分别将此距离记为dr、dl,取Rmin=min(dr,dl)。过Zi,Rmin所在一侧(P5所在一侧),作垂于RiCi的射线l。以Zi为圆心,Rmin为半径作弧,与射线l的交点ri为所求改航点。如图6所示。
策略一:当航路上只有单个的对流天气时,判断该对流天气是否与飞行计划路径相交。若没有相交,则按飞行计划路径继续飞行,如图7所示;若相交,则考虑对流天气位置,详见策略二、三、四。
策略二:若该单个对流天气存在飞行计划路径中间,先找改航点ri,若ri仍在航路内,则先沿飞行计划路径飞行,然后进行改航,再沿着飞行计划路径飞行,如图8所示;若ri不在航路内,则联系管制员,征询改航后是否回到原航路。若不回,则在其他航路内飞行。
策略三:若该单个对流天气存在于飞行计划路径起点,则一开始就进行改航,找到改航点ri,若ri仍在航路内,则做平行于计划路径的航段qiri,连接Eiri,如图9所示;若ri不在航路内,则联系管制员,征询改航后是否回到原航路。若不回,则在其他航路内飞行。
策略四:若该单个对流天气存在于飞行计划路径终点,先找到改航点ri,若改航点仍在航路内,则先沿飞行计划路径飞行,连接Siri,再以ri为起点作平行于飞行计划路径的航段,如图10所示;若ri不在航路内,则联系管制员,征询改航后是否回到原航路。若不回,则在其他航路内飞行。
策略五:当航路上存在多个对流天气时,计算前一个飞出交汇点Ci-1和后一个改航起始点Si之间的距离L,若L<7.4km,则进行连续改航。首先找到每个对流天气所对应的改航点,连接Si和距离最近的改航点,再连接两个改航点。若改航点都在航路内,则如图11所示;若存在改航点不在航路内,则联系管制员,征询改航后是否回到原航路。若不回,则在其他航路内飞行。
策略六:当航路上存在多个对流天气时,计算前一个飞出交汇点Ci-1和后一个改航起始点Si之间的距离L,若L≥7.4km,且改航点都在航路内,则进行组合改航,如图12所示。若改航之后的路径总距离≥1.5飞行计划路径的距离,则直接选择其他航路进行改航。若存在改航点不在航路内,则联系管制员,征询改航后是否回到原航路。若不回,则在其他航路内飞行。
注意:改航路径不能进入颜色区域。
Claims (4)
1.一种对流天气下的航路使用方法,其特征在于:包括如下步骤:
步骤1:指定航路T;
步骤2:获取并筛选历史航班数据;
步骤3:根据历史数据中的航路上的雷达轨迹和飞行计划路径判断是否改航,结合当时的气象情况,构建天气避让场;
步骤4:根据当前的气象条件,对照天气避让场,确定对流天气避让范围;
步骤5:给出避让原则,避让原则由对流天气的位置以及个数决定,令沿飞行计划路径方向最先与黄色天气避让区相交的点为飞入交汇点Ri,另一交汇点为飞出交汇点Ci,选取距离Ri和Ci最近的计划路径上的点分别作为改航起始点Si和改航结束点Ei,i=1,2,3,……,为飞行计划路径上依次遇到的对流天气区域编号,当飞行计划路径上只有单个的对流天气时,策略一:若对流天气存在于飞行计划路径中间,则先沿飞行计划路径飞行,然后找到改航点ri,将ri和改航起始点Si以及改航结束点Ei相连,作为改航的航段,以此进行改航,然后继续沿着飞行计划路径飞行;策略二:若该对流天气存在于飞行计划路径起点,则一开始就进行改航,找到改航点ri,并做平行于计划路径的航段qiri,连接ri和改航结束点Ei,然后沿着飞行计划路径继续飞行;策略三:若对流天气存在于飞行计划路径终点,则先沿飞行计划路径飞行,找到改航点ri,连接ri和改航起始点Si,再以ri为起点作平行于飞行计划路径的航段;当飞行计划路径上存在多个对流天气时,策略四:计算前一个飞出交汇点Ci-1和后一个改航起始点Si之间的距离L,若L<7.4km,则进行连续改航,首先找到每个对流天气所对应的改航点,连接Si和距离最近的改航点,再连接两个改航点;策略五:若L≥7.4km,则进行组合改航,策略六:若改航后的路径总距离≥1.5飞行计划路径总距离,则直接改航到其他航路。
2.根据权利要求1所述的对流天气下的航路使用方法,其特征在于:所述步骤2方法:获取航班历史数据,筛选在航路T范围内飞行的航班历史数据,包括历史雷达轨迹数据、飞行计划数据和气象指标数据,所述历史雷达轨迹数据包括起飞时间、飞行时间、着陆时间、呼号、地面速度、飞行高度、经纬度;所述飞行计划数据包括航班号、起飞时间、着陆时间、起飞机场、着陆机场、飞行计划路径;所述气象指标数据包括CR、ET和VIL。
3.根据权利要求2所述的对流天气下的航路使用方法,其特征在于:所述步骤3的方法:判断是否改航的依据是计算实际雷达轨迹和飞行计划路径之间的距离,进行航班聚类,寻找聚类中心线,聚类中心线到飞行计划路径的距离即为阈值,以此阈值来判断航班是否改航,针对改航航班,统计气象情况,根据历史轨迹改航的航班数,计算改航概率,改航概率=改航航班数/在此气象条件下的航班总数,生成以VIL为横轴,横坐标为0-0.05、0.05-0.14、0.14-0.7、0.7-3.5、3.5-6.9、6.9-12、12+,单位为kg/m3;deltaZ为纵轴,纵坐标为-10km到10km,每隔2.5km为一个间隔,每一格代表对应改航概率的天气避让概率。
4.根据权利要求3所述的对流天气下的航路使用方法,其特征在于:所述步骤4方法:当VIL介于3.5kg/m3和6.9kg/m3,deltaZ介于-7.5km和0时,设改航概率为0.4,由黄色表示;当VIL大于6.9kg/m3,deltaZ介于-10和-7.5km时,改航概率设为0.7,由红色表示,因此规定对流天气避让范围由外及内分别为黄色和红色,分别代表改航概率为0.4-0.7和0.7-1,确定避让对流范围的具体算法是:首先使用Canny边缘检测器检测出天气避让场的边界,然后用DBSCAN算法将边缘像素进行聚类,最后将每一种聚类结果包含在凸包中,形成避让范围,如果两个凸包的距离小于8km,则合并为一个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110657978.5A CN113483773B (zh) | 2021-06-11 | 2021-06-11 | 一种对流天气下的航路使用方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110657978.5A CN113483773B (zh) | 2021-06-11 | 2021-06-11 | 一种对流天气下的航路使用方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113483773A CN113483773A (zh) | 2021-10-08 |
CN113483773B true CN113483773B (zh) | 2024-02-13 |
Family
ID=77934904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110657978.5A Active CN113483773B (zh) | 2021-06-11 | 2021-06-11 | 一种对流天气下的航路使用方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113483773B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106323295A (zh) * | 2016-08-29 | 2017-01-11 | 中国船舶重工集团公司第七0九研究所 | 一种基于气象雷达数据的危险气象条件下的航空器改航方法 |
CN111127285A (zh) * | 2019-12-30 | 2020-05-08 | 南京航空航天大学 | 一种对流天气下两机场间航路通行能力获取的方法 |
CN111508279A (zh) * | 2020-03-18 | 2020-08-07 | 南京航空航天大学 | 一种天气避让场划设的方法 |
CN111984878A (zh) * | 2020-07-22 | 2020-11-24 | 南京航空航天大学 | 一种对流天气下机场终端区各进离港点流量控制方法 |
CN112071123A (zh) * | 2020-07-29 | 2020-12-11 | 南京航空航天大学 | 一种计算航路阻塞度的方法 |
-
2021
- 2021-06-11 CN CN202110657978.5A patent/CN113483773B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106323295A (zh) * | 2016-08-29 | 2017-01-11 | 中国船舶重工集团公司第七0九研究所 | 一种基于气象雷达数据的危险气象条件下的航空器改航方法 |
CN111127285A (zh) * | 2019-12-30 | 2020-05-08 | 南京航空航天大学 | 一种对流天气下两机场间航路通行能力获取的方法 |
CN111508279A (zh) * | 2020-03-18 | 2020-08-07 | 南京航空航天大学 | 一种天气避让场划设的方法 |
CN111984878A (zh) * | 2020-07-22 | 2020-11-24 | 南京航空航天大学 | 一种对流天气下机场终端区各进离港点流量控制方法 |
CN112071123A (zh) * | 2020-07-29 | 2020-12-11 | 南京航空航天大学 | 一种计算航路阻塞度的方法 |
Non-Patent Citations (3)
Title |
---|
Comprehensive evaluation for airspace classification;Wang Yinan;IEEE;全文 * |
危险天气下的战术改航方法;毛利民;彭瑛;贾子君;;华东交通大学学报(02);全文 * |
基于飞行受限区划设的航班改航研究;蒋昕;胡明华;张颖;田文;;华东交通大学学报(03);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113483773A (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106781707B (zh) | 一种针对低空中长距离转场飞行的航迹规划方法 | |
CN108196575B (zh) | 一种无人机任务分配及路线规划方法 | |
CN110111608B (zh) | 基于雷达轨迹构建机坪场面运动目标运行意图识别的方法 | |
Patron et al. | Flight trajectories optimization under the influence of winds using genetic algorithms | |
CN109493644A (zh) | 一种基于历史航迹数据挖掘的四维航迹推测方法 | |
Liang et al. | Integrated sequencing and merging aircraft to parallel runways with automated conflict resolution and advanced avionics capabilities | |
CN110930770A (zh) | 一种基于管制意图和飞机性能模型的四维航迹预测方法 | |
CN109542876B (zh) | 基于Hadoop数据挖掘航空器经验轨迹模型关键因子的提取方法 | |
CN103473955B (zh) | 一种基于图论和谱聚类算法的终端扇区划分方法 | |
CN104991895A (zh) | 一种基于三维空域网格的低空救援航空器航迹规划方法 | |
CN108803656B (zh) | 一种基于复杂低空的飞行控制方法及系统 | |
CN111583724A (zh) | 一种面向四维航迹运行的预战术阶段间隔管理方法 | |
CN112102650B (zh) | 改航航路生成方法、装置及存储介质 | |
CN110060513A (zh) | 基于历史轨迹数据的空中交通管制员工作负荷评估方法 | |
CN104615849A (zh) | 一种适用于通用航空的飞行计划评估系统及实现方法 | |
CN112880684A (zh) | 一种城市空间无人机安全航路规划方法 | |
CN111192481B (zh) | 一种基于碰撞风险的进离场程序无人机管控区边界确定方法 | |
CN106875755A (zh) | 一种基于复杂度的空中交通冲突管理方法和装置 | |
CN107067822A (zh) | 一种终端管制区的进离场航线动态管理方法和装置 | |
CN109918818A (zh) | 一种基于性能导航的pbn单发程序保护区划设方法 | |
CN110531782A (zh) | 用于社区配送的无人机航迹路径规划方法 | |
CN105868448A (zh) | 一种基于点融合系统的进近运行方法 | |
CN112071123A (zh) | 一种计算航路阻塞度的方法 | |
Meng et al. | Flight conflict resolution for civil aviation based on ant colony optimization | |
CN112083734B (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 |