CN113393688B - 一种基于队列长度预测的交叉口生态驾驶优化方法 - Google Patents
一种基于队列长度预测的交叉口生态驾驶优化方法 Download PDFInfo
- Publication number
- CN113393688B CN113393688B CN202110523552.0A CN202110523552A CN113393688B CN 113393688 B CN113393688 B CN 113393688B CN 202110523552 A CN202110523552 A CN 202110523552A CN 113393688 B CN113393688 B CN 113393688B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- intersection
- time
- optimization problem
- acceleration
- 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
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/09—Arrangements for giving variable traffic instructions
- G08G1/0962—Arrangements for giving variable traffic instructions having an indicator mounted inside the vehicle, e.g. giving voice messages
- G08G1/0967—Systems involving transmission of highway information, e.g. weather, speed limits
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/16—Anti-collision systems
- G08G1/166—Anti-collision systems for active traffic, e.g. moving vehicles, pedestrians, bikes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于队列长度预测的交叉口生态驾驶优化方法。首先,利用车路通信技术获得交通流和信号灯相位信息,基于运动冲击波模型预测队列排队长度;接下来,考虑排队长度,构建燃油最优车辆轨迹优化问题,并将其简化为停车等待和不停车优化问题,以获取实时解;然后,将获取的实时解作为智能汽车的期望加速度,并且为了保证车辆的安全,必要时,将根据智能驾驶员模型所得的安全加速度作为智能汽车的期望加速度,以避免碰撞。本发明方法使得智能汽车利用车路通信技术获得交通流与信号灯信息,以预测车队长度,所设计的生态驾驶算法避免了车辆在交叉口处的剧烈加减速,从而减少智能汽车通过交叉口时的油耗与行驶时间。
Description
技术领域
本发明属于网联环境下生态驾驶技术领域,特别涉及网联环境下交叉口处考虑排队长度的生态驾驶优化方法。
背景技术
近年来,汽车普及率的不断提高导致汽车能耗和排放显著增加,交通运输部门已成为温室气体排放和能源消耗的主要部门,约占全球石油消费量的一半。化石燃料在交通部门的广泛使用已经引起了公众对环境变化和能源可持续性的关注。如何有效降低运输带来的燃油消耗,提高燃油效率是目前迫切需要解决的问题。目前,已经有一些研究通过改进车辆技术和优化交通灯信号来降低油耗。研究表明,司机的驾驶习惯也会影响车辆油耗。生态驾驶技术可以帮助司机最大限度地减少十字路口的加速、减速和怠速,以降低油耗。
目前大多数交叉口的生态驾驶策略只考虑交通灯的信息,旨在减少加速、减速和停车等待产生的怠速时间。在实际的交通环境中,当车辆在红灯到达路口时,必须减速停车等待,会在交通流中产生向上游扩散的排队冲击波,阻碍车辆通行。受排队长度的影响,即使采用生态驾驶策略,也会为了避免碰撞而发生减速停车的情况。然而,很少有研究考虑车辆排队长度对主干道生态行驶的影响,可能导致速度轨迹次优解。其次,考虑到目前很少使用车-车通信,车路通信技术采用车辆和道路基础设施通信的方式进行传递信息目前的应用范围更广,可以在更大范围内优化车辆轨迹。
发明内容
本发明的目的是提供一种基于队列预测的交叉口生态驾驶优化方法,利用了车路通信技术,使得智能汽车能够获取交通信息以及实时队列长度的预测,从而在通过交叉口时能够最大限度的减小燃油消耗和通行时间。
以上目的通过以下的技术方案实现:
一种基于队列长度预测的交叉口生态驾驶优化方法,包括如下步骤:
步骤1、实时预测队列长度:
智能汽车进入算法预定优化区间时,根据智能汽车与交叉口基础设施通信后获得的交通流信息与信号灯时间,获得排队冲击波、消散冲击波和驶离冲击波速度,基于运动学冲击波模型给出队列长度预测算法,来实时预测交叉口车辆排队长度,进而预测车辆到达交叉口前的最大排队长度。所述的交通流信息包括每个周期车流量和车辆平均驶入速度。
步骤2、构建燃油最优车辆轨迹优化问题:
将燃油最优问题简化为不停车优化子问题和停车等待优化问题,先求解并采用不停车优化问题的解,若不停车优化问题无解,再求解并采用停车等待优化问题的解。
步骤3、考虑预测排队长度的生态驾驶算法:
将燃油最优车辆轨迹优化问题的解作为智能汽车的期望加速度,并且为了保证车辆的安全,当智能汽车与前车距离小于安全行驶距离时,将智能驾驶员模型所得的安全加速度作为智能汽车的期望加速度,从而避免发生碰撞。
所述的队列长度预测算法为:
其中,sp(t)为预测的实时队列长度,v0为车辆自由驾驶的平均速度,为第n个定时周期中红灯开始时间,为第n个定时周期中排队长度最大的时间,为第n个定时周期中队列完全消散的时间,为第n+1个定时周期中红灯开始时间。
步骤2中所述的燃油最优车辆轨迹优化问题为:
其中Fuel(v(t),a(t),θ(t))为燃油消耗模型:
其中α0,α1,α2由智能汽车类型决定的系数,P(t)是智能汽车输出功率,v(t),a(t),θ(t)分别为速度,加速度,道路坡度。
将该最优问题优化为两个问题:
1、停车等待优化问题:
智能汽车在时间t0以速度v0进入交叉口的控制路段后,先以恒定加速度a1减速,然后以恒定的速度行驶,最后以恒定加速度amin减速停止。当消散冲击波传播到该智能汽车处,智能汽车开始以恒定加速度as加速到饱和交通流速度vm,行驶通过交叉口。停车等待优化子问题表示为:
其中的约束为:
θ(t)=β
amin≤a1≤0
0≤as≤amax
其中,t1为以a1匀减速结束的时间,t2为开始以amin减速的时间,tq为车辆减速到0的时间,td为开始加速的时间,ta为加速到饱和流速的时间,tl为驶出交叉口的时间,amin为最大减速度,amax为最大加速度。
2、不停车优化问题:
其中的约束为:
θ(t)=β
amin≤a2≤0
amin≤a3≤amax
其中,t'1和t'a分别为以a2匀减速结束和以a3加速度至vm的时间。
步骤3中所述的考虑队列长度预测的生态驾驶算法如下:
1、如果智能汽车能够以其初始速度v0通过停车线,而不撞上车队或遇到红灯,则不会对其采取任何控制措施。
2、如果智能汽车在前车车队到达步骤一预测的排队最大长度之前到达交叉口,则优先采用不停车优化问题的解来控制车辆;若不停车优化问题无解,则采用停车等待优化问题的解来控制车辆。当智能车辆与前方车辆的距离小于安全距离时,将切换到智能驾驶员跟驰模型继续行驶。
控制车辆与前车的安全距离:
其中dmin为最小安全距离,Ht为车头时距,vp为前车速度。
本发明所产生的有益效果如下:
本发明提出了一种基于队列长度预测的交叉口生态驾驶优化方法。该生态驾驶优化方法要求,首先,利用车路通信技术获得交通流和信号灯相位信息,基于运动冲击波模型预测队列排队长度;接下来,考虑排队长度,构建燃油最优车辆轨迹优化问题,并将其简化为停车等待和不停车优化问题,以获取实时解;然后,将获取的实时解作为智能汽车的期望加速度,并且为了保证车辆的安全,必要时,将根据智能驾驶员模型所得的安全加速度作为智能汽车的期望加速度,以避免碰撞。
由此可知,该发明使得智能汽车利用车路通信技术获得交通流与信号灯信息,以预测车队长度,所设计的生态驾驶算法避免了车辆在交叉口处的剧烈加减速,从而减少智能汽车通过交叉口时的油耗与行驶时间。
附图说明
图1是交通冲击波传播示意图;
图2是本发明实施例中信号交叉口的交通动态图;
图3是本发明实施例中智能车辆的速度轨迹图;
图4是交叉口不同智能车辆速度分布图;
图5是不同市场渗透率下交通流的油耗和延时减少程度图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提出了一种基于队列长度预测的交叉口生态驾驶优化方法。该发明包括以下步骤:
图2是本发明实施例中信号交叉口的交通动态图;
步骤2,通过求解燃油最优车辆轨迹问题对原本需要停车等待的车辆进行求解优化。如图3所示,将燃油最优车辆轨迹问题简化为2个子问题进行求解,将该最优问题优化为两个问题:
考虑车辆停车通过交叉口:
智能汽车在时间t0以速度v0进入交叉口的控制路段后,先以恒定加速度a1减速,然后以恒定的速度行驶(如果有可能),最后以恒定加速度amin减速停止。当消散冲击波传播到该智能汽车处,智能汽车开始以恒定加速度as加速到饱和交通流速度vm,行驶通过交叉口。停车等待优化子问题表示为:
其中的约束为:
θ(t)=β
amin≤a1≤0
0≤as≤amax
其中,t1为以a1匀减速的时间,t2为开始以amin减速的时间,tq为车辆减速到0的时间,td为开始加速的时间,ta为加速到饱和流速的时间,tl为驶出交叉口的时间,amin为最大减速度,amax为最大加速度。
其中,t1,t2,tq,td,ta,tl能够按如下公式计算:
ta=td+τs
其中,D为优化区间,L为小车长度,sp(tq)为车辆停车前方的排队长度,dmin为最小安全距离,τs为车辆以加速度as加速到vm的时间。
考虑车辆不停车通过交叉口:
类似于停车等待优化问题,智能汽车在时间t0以速度v0进入交叉口控制路段后,先以恒定加速度a2减速,然后以恒定速度行驶至时间当检测到前方队列驶离,智能汽车则以恒定加速度a3加速到vm,通过交叉口。不停车优化问题表示为:
其中的约束为:
θ(t)=β
amin≤a2≤0
amin≤a3≤amax
其中,t'1和t'a分别为以a2匀减速结束和以a3加速度至vm的时间;t'1,t'a,tl能够按如下公式计算:
步骤3,将考虑预测排队长度的生态驾驶算法得出的速度轨迹实时更新给该智能车辆。
1、如果智能汽车能够以其初始速度v0通过停车线,而不撞上车队或遇到红灯,则不会对其采取任何控制措施。
2、如果智能汽车在前车车队到达步骤一预测的排队最大长度之前到达交叉口,则优先采用不停车优化问题的解来控制车辆;若不停车优化问题无解,则采用停车等待优化问题的解来控制车辆。当智能车辆与前方车辆的距离小于安全距离时,将切换到智能驾驶员跟驰模型继续行驶。
控制车辆与前车的安全距离:
其中dmin为最小安全距离,Ht为车头时距,vp为前车速度。
智能驾驶员模型:
利用交通仿真软件VISSIM与MATLAB联合仿真验证本发明的有效性:交叉口入口处上游道路长度设置为500米,车辆进入上游道路的初始速度设置为70公里/小时;采用固定定时方法控制交叉口的交通灯,红灯和绿灯的持续时间分别为40秒和40秒;在交叉口上游500米处设置一个检测器,以便交叉口能够实时采集交通流量信息。
一、生态驾驶算法对单辆智能车辆的控制效果
如图4所示,红色线表示不受生态驾驶算法控制的速度轨迹,a中蓝色曲线表示生态驾驶算法获得的智能车辆A的停车等待速度轨迹,b中蓝色曲线表示生态驾驶算法获得的智能车辆B的不停车通行速度轨迹。与不受生态驾驶算法控制相比,A的油耗和行程时间分别降低了9.4%和0.3%,B的油耗和行程时间分别降低了13.9%和3.7%。由此可见,生态驾驶算法对单车的控制非常有效,会对减少油耗和时延产生积极影响。
二、生态驾驶算法对交通流的影响
研究该生态驾驶算法在不同市场渗透率(MPRs)条件下对整个交通流的影响。通过车路通信技术预测队列长度,在不同市场渗透率条件下,采用生态驾驶算法获取的车辆速度轨迹如图5所示。在图5中,当MPRs分别为20%、50%、80%和100%时油耗可分别节省3%、4.5%、6.0%和6.3%。随着市场渗透率提升,车辆整体油耗在减少,整体延误也将略有改善,这说明生态驾驶算法能够对整个道路的交通流产生积极的影响。
应当指出,以上详细描述了本发明的实施过程,但本发明并不限于上述实施方式中的具体细节。在不脱离本发明原理的前提下,还可以对具体细节做出调整,比如加减速限速以及安全行驶距离等,这些细节也应该视为本发明的保护范围。
Claims (2)
1.一种基于队列长度预测的交叉口生态驾驶优化方法,其特征在于,包括如下步骤:
步骤1、实时预测队列长度:
智能汽车进入算法预定优化区间时,根据智能汽车与交叉口基础设施通信后获得的交通流信息与信号灯时间,获得排队冲击波、消散冲击波和驶离冲击波速度,基于运动学冲击波模型给出队列长度预测算法,来实时预测交叉口车辆排队长度,进而预测车辆到达交叉口前的最大排队长度;所述的交通流信息包括每个周期车流量和车辆平均驶入速度;
步骤2、构建燃油最优车辆轨迹优化问题:
将燃油最优问题简化为不停车优化子问题和停车等待优化问题,先求解并采用不停车优化问题的解,若不停车优化问题无解,再求解并采用停车等待优化问题的解;
步骤3、考虑预测排队长度的生态驾驶算法:
将燃油最优车辆轨迹优化问题的解作为智能汽车的期望加速度,并且为了保证车辆的安全,当智能汽车与前车距离小于安全行驶距离时,将智能驾驶员模型所得的安全加速度作为智能汽车的期望加速度,从而避免发生碰撞;
所述的队列长度预测算法为:
其中,sp(t)为预测的实时队列长度,v0为车辆自由驾驶的平均速度,为第n个定时周期中红灯开始时间,为第n个定时周期中排队长度最大的时间,为第n个定时周期中队列完全消散的时间,为第n+1个定时周期中红灯开始时间;
步骤2中所述的燃油最优车辆轨迹优化问题为:
其中Fuel(v(t),a(t),θ(t))为燃油消耗模型:
其中α0,α1,α2由智能汽车类型决定的系数,P(t)是智能汽车输出功率,v(t),a(t),θ(t)分别为速度,加速度,道路坡度;
将该最优问题优化为两个问题:
1、停车等待优化问题:
智能汽车在时间t0以速度v0进入交叉口的控制路段后,先以恒定加速度a1减速,然后以恒定的速度行驶,最后以恒定加速度amin减速停止;当消散冲击波传播到该智能汽车处,智能汽车开始以恒定加速度as加速到饱和交通流速度vm,行驶通过交叉口;停车等待优化子问题表示为:
其中的约束为:
θ(t)=β
amin≤a1≤0
0≤as≤amax
其中,t1为以a1匀减速结束的时间,t2为开始以amin减速的时间,tq为车辆减速到0的时间,td为开始加速的时间,ta为加速到饱和流速的时间,tl为驶出交叉口的时间,amin为最大减速度,amax为最大加速度;
2、不停车优化问题:
其中的约束为:
θ(t)=β
amin≤a2≤0
amin≤a3≤amax
其中,t′1和t′a分别为以a2匀减速结束和以a3加速度至vm的时间;
停车等待优化问题中的t1,t2,tq,td,ta,tl能够按如下公式计算:
v0+aI(t1-t0)=-amin(tq-t2)
ta=td+τs
其中,D为优化区间,L为小车长度,sp(tq)为车辆停车前方的排队长度,dmin为最小安全距离,τs为车辆以加速度as加速到vm的时间;
不停车优化问题中的t′1,t′a,tl能够按如下公式计算:
2.根据权利要求1所述的一种基于队列长度预测的交叉口生态驾驶优化方法,其特征在于,步骤3中所述的考虑队列长度预测的生态驾驶算法如下:
(1)、如果智能汽车能够以其初始速度v0通过停车线,而不撞上车队或遇到红灯,则不会对其采取任何控制措施;
(2)、如果智能汽车在前车车队到达步骤一预测的排队最大长度之前到达交叉口,则优先采用不停车优化问题的解来控制车辆;若不停车优化问题无解,则采用停车等待优化问题的解来控制车辆;当智能车辆与前方车辆的距离小于安全距离时,将切换到智能驾驶员跟驰模型继续行驶;
控制车辆与前车的安全距离:
其中dmin为最小安全距离,Ht为车头时距,vp为前车速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110523552.0A CN113393688B (zh) | 2021-05-13 | 2021-05-13 | 一种基于队列长度预测的交叉口生态驾驶优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110523552.0A CN113393688B (zh) | 2021-05-13 | 2021-05-13 | 一种基于队列长度预测的交叉口生态驾驶优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113393688A CN113393688A (zh) | 2021-09-14 |
CN113393688B true CN113393688B (zh) | 2022-07-19 |
Family
ID=77617164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110523552.0A Active CN113393688B (zh) | 2021-05-13 | 2021-05-13 | 一种基于队列长度预测的交叉口生态驾驶优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113393688B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115035724B (zh) * | 2022-07-13 | 2022-12-02 | 吉林大学 | 一种基于生态编队的物流车准时化运输车路协同管控方法 |
CN115497314B (zh) * | 2022-09-03 | 2023-10-24 | 河海大学 | 一种智能网联汽车不停车通过交叉口的生态驾驶方法 |
CN116580561B (zh) * | 2023-05-24 | 2024-02-06 | 交通运输部规划研究院 | 一种基于加速度模型和短期轨迹数据的交叉口流量释放预测方法 |
CN116524722B (zh) * | 2023-06-26 | 2023-10-13 | 浙江大学 | 面向混行交通流的车辆生态驾驶控制方法、电子设备 |
CN117133140B (zh) * | 2023-07-27 | 2024-07-19 | 北京工业大学 | 一种信号交叉口智能网联车辆双模式最优速度轨迹控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106530695A (zh) * | 2016-11-09 | 2017-03-22 | 宁波大学 | 基于车联网的城市干道车辆行程时间实时预测方法 |
CN111275987A (zh) * | 2020-01-21 | 2020-06-12 | 东南大学 | 一种考虑路口待行队列影响的汽车行驶车速优化方法 |
CN111341152A (zh) * | 2020-03-03 | 2020-06-26 | 东南大学 | 考虑待行队列影响与安全避撞的网联汽车绿色通行系统与方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190206247A1 (en) * | 2018-01-03 | 2019-07-04 | Xiaofeng Xie | Smart In-Vehicle Decision Support Systems and Methods with V2I Communications for Driving through Signalized Intersections |
CN111791887B (zh) * | 2020-07-03 | 2021-06-04 | 北京理工大学 | 一种基于分层式车速规划的车辆节能驾驶方法 |
-
2021
- 2021-05-13 CN CN202110523552.0A patent/CN113393688B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106530695A (zh) * | 2016-11-09 | 2017-03-22 | 宁波大学 | 基于车联网的城市干道车辆行程时间实时预测方法 |
CN111275987A (zh) * | 2020-01-21 | 2020-06-12 | 东南大学 | 一种考虑路口待行队列影响的汽车行驶车速优化方法 |
CN111341152A (zh) * | 2020-03-03 | 2020-06-26 | 东南大学 | 考虑待行队列影响与安全避撞的网联汽车绿色通行系统与方法 |
Non-Patent Citations (3)
Title |
---|
基于V2X的车辆编队与交叉口信号协同优化控制方法研究;章权;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20210215(第02期);正文第10-18,22-27页 * |
基于生态驾驶工况曲线函数的坡道和交叉口车辆轨迹优化方法;黎明;《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑》;20180615(第06期);正文第14,63-64页 * |
章权.基于V2X的车辆编队与交叉口信号协同优化控制方法研究.《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》.2021,(第02期), * |
Also Published As
Publication number | Publication date |
---|---|
CN113393688A (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113393688B (zh) | 一种基于队列长度预测的交叉口生态驾驶优化方法 | |
CN110085037B (zh) | 一种车路协同环境下交叉口信号控制及车速引导系统 | |
CN108831143B (zh) | 基于车路协同技术的信号交叉口车队速度引导方法 | |
CN106997690B (zh) | 一种车联网环境下高速公路车辆非强制换道控制方法 | |
CN110910646B (zh) | 用于交叉口的无人驾驶公交车辆的协同控制方法 | |
CN108765982A (zh) | 车路协同环境下信号控制交叉口车速引导系统及引导方法 | |
CN111383481A (zh) | 城市拥堵路口智能网联汽车绿色通行车速优化方法 | |
CN114973733B (zh) | 一种信号交叉口处混合流下网联自动车轨迹优化控制方法 | |
CN107730922A (zh) | 一种单向干线绿波协调控制自适应调整方法 | |
CN111341152B (zh) | 考虑待行队列与安全避撞的网联汽车绿色通行系统与方法 | |
CN112767715B (zh) | 一种交叉路口交通信号灯与智能网联汽车的协同控制方法 | |
CN113257007B (zh) | 交叉口共享自动驾驶车道的信号控制方法及系统 | |
CN114419903B (zh) | 一种智能网联汽车队列路口通行控制方法、装置及车辆 | |
CN108986509A (zh) | 一种基于车路协同的城市区域路径实时规划方法 | |
CN113570875B (zh) | 一种绿波车速计算方法、装置,设备和存储介质 | |
CN110176138B (zh) | 一种路口级的主动交通诱导方法 | |
CN111932910A (zh) | 一种智能车路协同环境下实时动态可变车道安全控制方法 | |
CN114120666B (zh) | 一种基于车路协同的智能车速诱导系统 | |
CN116189462B (zh) | 一种面向混合交通流的车辆轨迹与交通信号协同控制方法 | |
CN115273500A (zh) | 一种信号交叉口网联车辆轨迹优化引导方法及系统 | |
CN112820126B (zh) | 一种非侵入式导向运输车辆路权优先运行控制及仿真方法 | |
CN110085025A (zh) | 一种快速公交多模态运行速度优化方法 | |
CN115565390B (zh) | 智能网联汽车多车道队列通行控制方法、系统及计算机可读存储介质 | |
CN107215333A (zh) | 一种采用滑行模式的经济性协同自适应巡航策略 | |
CN115497315A (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 |