CN106846812A - 一种十字路口的交通流量评估方法 - Google Patents

一种十字路口的交通流量评估方法 Download PDF

Info

Publication number
CN106846812A
CN106846812A CN201710149512.8A CN201710149512A CN106846812A CN 106846812 A CN106846812 A CN 106846812A CN 201710149512 A CN201710149512 A CN 201710149512A CN 106846812 A CN106846812 A CN 106846812A
Authority
CN
China
Prior art keywords
time
car
distance
acceleration
rear car
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
CN201710149512.8A
Other languages
English (en)
Other versions
CN106846812B (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 CN201710149512.8A priority Critical patent/CN106846812B/zh
Publication of CN106846812A publication Critical patent/CN106846812A/zh
Application granted granted Critical
Publication of CN106846812B publication Critical patent/CN106846812B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0133Traffic data processing for classifying traffic situation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种十字路口的交通流量评估方法,采用物理建模方案,引入安全跟驰距离与最大车间距,通过研究启动迟滞时间与车间距关系及车辆启动模式,估算一次绿灯时长内通过路口的车流量,建立了前后车跟驰速度、加速度关系及绿灯时间长度、车间距与车流量的关系,并与实际观测结果进行了对比。

Description

一种十字路口的交通流量评估方法
技术领域
本发明涉及一种十字路口的交通流量评估方法,属于交通技术领域。
背景技术
十字路口交通管理不仅是城市交通的瓶颈,更成为一个对经济、社会有重大影响的课题。国内外已经有一定的研究,虽然有一些经典的交通流理论,但大多以经验公式为基础,并没有形成普适的理论方法,难以准确的对十字路口的交通流量进行优化。例如格林希尔治(Greenshields)模型:其中u表示实际车流速度,uf为自由流车速,k表示车流密度,kj为阻塞密度。格林伯(Greenberg)模型是从实际数据出发进行归纳拟合得到的,对于小流量的模型并不适用。而《美国通行能力手册》给出的工程方法如下:当绿灯开启时,驾驶员做出反应和车辆启动都需要一个时间延迟,这个过程会影响到后几辆车(一般估计为五辆),而后面的车流符合饱和流的要求。这样就能够求出在红绿灯时间给定的情况下通过的车流量。这种方法建立在经验公式上,没有完备的理论支撑,适用性较差。
发明内容
本发明的目的是为了解决上述问题,提出一种十字路口的交通流量评估方法,采用物理建模方案,引入安全跟驰距离与最大车间距,通过研究启动迟滞时间与车间距关系及车辆启动模式,估算一次绿灯时长内通过路口的车流量,建立了前后车跟驰速度、加速度关系及绿灯时间长度、车间距与车流量的关系,并与实际观测结果进行了对比。
一种十字路口的交通流量评估方法,包括以下几个步骤:
步骤一:获取安全跟驰距离与前后车速度、加速度之间的关系;
步骤二:获取启动迟滞时间与车间距的关系;
步骤三:获取后车加速度随前车加速度以及前后两车速度的变化规律;
步骤四:依据上述理论方法,建立交通流量评估(计算)模型,获取单位时长内通过路口的车流量;
本发明的优点在于:
(1)以刚体运动学理论为基础,形成了较为严谨的交通流评估方法;
(2)以刚体运动学为基础,建立了安全跟驰距离的概念,并以这一概念作为交通流评估的安全判据;
(3)以刚体运动学为基础,建立了车间距、前后车速度和加速度的关系,从而形成了具有理论严谨性的跟驰模式;
(4)以上述跟驰模型为基础,考虑启动因素,建立了车辆启动后行驶距离的计算方法,从而构成了以严谨理论为基础的交通流量评估模型;
附图说明
图1是十字路口交通流量评估流程图;
图2是一次绿灯时长内通过的车流量计算流程图;
图3是车流量-红灯时长图;
图4是车流量-车间距图;
图5是实测车流量-红灯时间对比图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明的一种十字路口的交通流量评估方法,如图1所示,包括以下几个步骤:
步骤一:获取安全跟驰距离与前后车速度、加速度之间的关系;
所谓“安全跟驰距离”,是汽车行驶时,能够保证在出现意外紧急刹车时,前后辆车恰恰不发生追尾事故的前后车间距。设某时刻前车的速度为vi,后车的速度为vi+1,前车突然以加速度a0′刹车,后车在经过刹车迟滞时间Ts后以最大加速度a0急刹车,则从前车紧急刹车到两车都停车,前车行驶的距离为:
后车行驶的距离可以分为两部分,当前车开始刹车时,后车需要刹车迟滞时间Ts才能做出反应,在这段时间内,后车依然按照前面的加速度进行匀加速行驶,这段时间内行驶的距离为S1。经过刹车迟滞时间后,后车迅速以最大的加速度进行减速,这段时间内行驶的距离为S2,设加速与减速时的加速度数值相同,则
则从前车紧急刹车到两车都停车,后车行驶的距离一共为:
为了保证两车不发生碰撞,两车的安全跟驰距离Ls为:
步骤二:获取启动迟滞时间与车间距的关系
对于排队通过道路交叉口的车队,启动时间间隔随两车间距增大而减小。定义车间距为前车车尾到后车车头的距离。当红灯变为绿灯时,认为第一辆车启动时的启动迟滞时间为零。设当两车的间距为基准车间距L0时,后车启动迟滞时间为TSS=T0;当两车间距增大到最大车间距LM后,后车启动迟滞时间为TSS=0。即前车不动,后车便进行加速,经过刹车迟滞时间Ts后再以最大加速度进行减速,而前后车不致发生碰撞,即:
设启动迟滞时间与车间距成反比,可以得到后车的启动迟滞时间随着车间距LP变化的关系式:
步骤三:获取后车加速度随前车加速度以及前后两车速度的变化规律
排队通过十字道路交叉口的车辆,后车启动后,前后两车的相对速度和距离随时会发生变化,为保持前后两车间的安全跟驰距离,后车须主动调整车速。设某时刻前后两车的车速分别为vi和vi+1,在某个很短的时间间隔Δt内,前车车速变化为Δvi,经过Ts的刹车迟滞时间后,后车Δt内车速变化为Δvi+1。在这一变化过程中,两车由于速度变化行驶的距离分别为:
Δsi=viΔt,Δsi+1=vi+1Δt (7)
由于Δt很小,这里假设两车在该时间段内均以前一时刻的速度行驶。变速过程后,两车间的安全跟驰距离为:
其中分别为后车与前车行驶的距离。为保持前后两车的安全跟驰距离,两车间距变化为:
Ls′-Ls=Δsi-Δsi+1 (9)从而可以得到:
Δt是一个小量,所以在这段时间内速度的变化量Δvi也为一个小量,忽略二次小项并由可得:
可以据此求出后车的加速度ai+1随前车加速度ai以及前后两车速度的变化规律为:
其中:ai≤a0
步骤四:依据上述理论方法,建立交通流量评估(计算)模型,获取单位时长内通过路口的车流量;
具体流程如图2所示,包括以下几个环节:
a.设定刹车迟滞时间TS、最大启动迟滞时间T0、前车最大加速度a0′、后车最大加速度a0、基准车间距L0、车间距LP、车辆长度Lc、黄灯时间恒为ty、红灯时长与绿灯时长相等,分别为tr、tg
b.绿灯亮时,第一辆车立即以加速度a0′启动,第二辆车经过步骤二(6)式中所定义的启动迟滞时间TSS后启动,此时令i=2,i表示车辆编号。
c.设经过时间(i-1)TSS后,第i辆车开始启动,其加速度变化规律符合式(12),且速度达到v0后不再增加。
d.对第i辆车的加速度,即式(12)进行积分,得到第i辆车在启动后tg-(i-1)TSS时间内的行驶距离为Ltgi
e.判断其行驶距离是否满足条件Ltgi>(n-1)(Lc+Lp),若满足,则第i辆车可以通过红绿灯,令i=i+1,返回环节c;若不满足,则在Tg时长的绿灯期间通过十字道路交叉口的车辆数为n=i-1。
f.得到每小时通过十字道路交叉口的车流量为qm=3600n/(tg+tr+ty)
通过以上过程,可以在环节a中给定的初始条件下,最终得到通过十字道路交叉口的车流量,可用于十字道路交叉口的车流量评估,同时可以用于分析环节a中各初始条件对车流量的影响,以便提高十字道路交叉口的车流量。
实施例:
采用本发明所述技术方案,对十字道路交叉口车流量进行模拟,并与实际对比,分析影响因素。
使用MATLAB进行建模计算,设刹车迟滞时间TS=1.2s,a0=8m/s2,a0′=5m/s2;最大启动迟滞时间T0=1.5s,L0=0m;另外取车辆长度Lc=4.5m,限制最大车速为v0=12.5m/s,黄灯时间恒为ty=4s,红灯时长与绿灯时长相等。分别设置车间距为1m,1.6m,4.4m和6m,改变红灯时长,得到四种车间距对应的“车流量-红灯时长”图如附图3。图中横坐标为红灯时长,纵坐标为每小时车流量。可看出车间距LP为4.4m时的车流量要大于车间距为1.6m与1m的情况,同时还大于车间距为6m的情况,这说明随着车间距的增大,车流量会先增大后减小。同时从图还可以得到,随着红灯时间延长,车流量增长速度逐渐放缓,最终趋近于一个稳定值。因此当交通流量较大时,建议将红灯时间长度设置为55秒左右,这样既能获得一个较优的车流量值,又不会使得排队等待时间过长,避免造成另一个方向的堵塞。这个红灯时间也比较符合实际情况下的红灯时间的设置情况,说明现实中红灯时间长度较为合理。
其他条件不变,当设置红灯时长等于绿灯时长等于55秒,改变车间距LP时,得到模拟结果“车流量-车间距”图如附图4。图中横坐标为车间距,纵坐标为每小时车流量。可以看出,随着车间距的增大车流量先增加后减小,且车间距取4~6米时车流量最大,最大值比实际观测车间距为1.6m时的车流量要大8%左右。说明当十字道路交叉口车流量大时,如果等待信号灯的汽车车间距过小,将使得车辆无法同时高速启动,造成车流量较小。
实际对比情况:
实际道路较为拥堵时,等待信号灯的汽车车间距约为1.6m,附图5给出了作者在北京市海淀区学院桥交叉口实地观测所得车流量、美国交通研究委员会道路通行能力手册所给车流量经验公式评估结果和本发明计算结果。图中横坐标为信号灯红灯时长,纵坐标为每小时车流量。
对比经验模型可以看出,本发明对车流量的模拟比经验模型更贴近实际,同时由于本发明结果来自于理论模型的理想启动状况(无刹车减速等意外情况),所以得到的车流量也比实际的车流量略大,符合常理。
综上所述,本发明是一种兼具理论严谨性和工程适用性的十字道路交叉口交通流量评估方法,本方法评估结果与实际观测结果取得了较好的一致性,为提高十字道路交叉口交通流量提供了理论分析和优化建议,将产生良好的社会效益。

Claims (1)

1.一种十字路口的交通流量评估方法,包括以下几个步骤:
步骤一:获取安全跟驰距离与前后车速度、加速度之间的关系;
设某时刻前车的速度为vi,后车的速度为vi+1,前车以加速度a′0刹车,后车在经过刹车迟滞时间Ts后以最大加速度a0急刹车,则从前车紧急刹车到两车都停车,前车行驶的距离为:
s c i = v i 2 2 a 0 ′ - - - ( 1 )
后车行驶的距离分为两部分,当前车开始刹车时,后车需要刹车迟滞时间Ts做出反应,在这段时间内,后车依然按照前面的加速度进行匀加速行驶,这段时间内行驶的距离为S1;经过刹车迟滞时间后,后车迅速以最大的加速度进行减速,这段时间内行驶的距离为S2,设加速与减速时的加速度数值相同,则
S 1 = v i + 1 T s + 1 2 a 0 T s 2 , S 2 = ( v i + 1 + a 0 T s ) 2 2 a 0 - - - ( 2 )
则从前车紧急刹车到两车都停车,后车行驶的距离一共为:
s c i + 1 = S 1 + S 2 = v i + 1 T s + 1 2 a 0 T s 2 + ( v i + 1 + a 0 T s ) 2 2 a 0 - - - ( 3 )
则两车的安全跟驰距离Ls为:
L s = s c i + 1 - s c i = v i + 1 T s + 1 2 a 0 T s 2 + ( v i + 1 + a 0 T s ) 2 2 a 0 - v i 2 2 a 0 ′ - - - ( 4 )
步骤二:获取启动迟滞时间与车间距的关系;
设车间距为前车车尾到后车车头的距离,当红灯变为绿灯时,设第一辆车启动时的启动迟滞时间为零,当两车的间距为基准车间距L0时,后车启动迟滞时间为TSS=T0,当两车间距增大到最大车间距LM后,后车启动迟滞时间为TSS=0,即前车不动,后车便进行加速,经过刹车迟滞时间Ts后再以最大加速度进行减速,而前后车不致发生碰撞,即:
L M = 1 2 a 0 T s 2 + ( a 0 T s ) 2 2 a 0 - - - ( 5 )
设启动迟滞时间与车间距成反比,得到后车的启动迟滞时间随着车间距LP变化的关系式:
T S S = L M - L P L M - L 0 T 0 - - - ( 6 )
步骤三:获取后车加速度随前车加速度以及前后两车速度的变化规律;
排队通过十字道路交叉口的车辆,后车启动后,前后两车的相对速度和距离随时会发生变化,为保持前后两车间的安全跟驰距离,后车须主动调整车速,设某时刻前后两车的车速分别为vi和vi+1,在时间间隔Δt内,前车车速变化为Δvi,经过Ts的刹车迟滞时间后,后车Δt内车速变化为Δvi+1;在这一变化过程中,两车由于速度变化行驶的距离分别为:
Δsi=viΔt,Δsi+1=vi+1Δt (7)
假设两车在该时间段内均以前一时刻的速度行驶,变速过程后,两车间的安全跟驰距离为:
L s ′ = s c i + 1 - s c i = ( v i + 1 + Δv i + 1 ) T s + 1 2 a i + 1 T s 2 + ( v i + 1 + Δv i + 1 + a i + 1 T s ) 2 2 a 0 - ( v i + Δv i ) 2 2 a 0 ′ - - - ( 8 )
其中:分别为后车与前车行驶的距离,两车间距变化为:
L′s-Ls=Δsi-Δsi+1 (9)
得到:
Δv i + 1 T s = 2 v i + 1 Δv i + 1 + Δv i + 1 2 + 2 a i + 1 Δv i + 1 T s 2 a 0 - 2 v i Δv i + Δv i 2 2 a 0 ′ = ( v i - v i + 1 ) Δ t - - - ( 10 )
忽略二次小项并由可得:
a i + 1 T s + v i + 1 a i + 1 + a i + 1 2 T s a 0 - v i a i a 0 ′ = v i - v i + 1 - - - ( 11 )
获取后车的加速度ai+1随前车加速度ai以及前后两车速度的变化规律为:
a i + 1 = - ( T s + V i + 1 a 0 ) + ( T s - V i + 1 a 0 ) 2 + 4 T s a 0 ( V i + V i a i a 0 i ) 2 T s / a 0 - - - ( 12 )
其中:ai≤a0
步骤四:建立交通流量评估模型,获取单位时长内通过路口的车流量;
包括以下几个环节:
a.设定刹车迟滞时间TS、最大启动迟滞时间T0、前车最大加速度a′0、后车最大加速度a0、基准车间距L0、车间距LP、车辆长度Lc、黄灯时间恒为ty、红灯时长与绿灯时长相等,分别为tr、tg
b.绿灯亮时,第一辆车立即以加速度a′0启动,第二辆车经过启动迟滞时间TSS后启动,此时令i=2,i表示车辆编号;
c.设经过时间(i-1)TSS后,第i辆车开始启动,其加速度变化规律符合式(12),且速度达到v0后不再增加;
d.对第i辆车的加速度,即式(12)进行积分,得到第i辆车在启动后tg-(i-1)TSS时间内的行驶距离为Ltgi
e.判断其行驶距离是否满足条件Ltgi>(n-1)(Lc+Lp),若满足,则第i辆车可以通过红绿灯,令i=i+1,返回环节c;若不满足,则在Tg时长的绿灯期间通过十字道路交叉口的车辆数为n=i-1;
f.得到每小时通过十字道路交叉口的车流量为qm=3600n/(tg+tr+ty)。
CN201710149512.8A 2017-03-14 2017-03-14 一种十字路口的交通流量评估方法 Expired - Fee Related CN106846812B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710149512.8A CN106846812B (zh) 2017-03-14 2017-03-14 一种十字路口的交通流量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710149512.8A CN106846812B (zh) 2017-03-14 2017-03-14 一种十字路口的交通流量评估方法

Publications (2)

Publication Number Publication Date
CN106846812A true CN106846812A (zh) 2017-06-13
CN106846812B CN106846812B (zh) 2019-01-25

Family

ID=59144505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710149512.8A Expired - Fee Related CN106846812B (zh) 2017-03-14 2017-03-14 一种十字路口的交通流量评估方法

Country Status (1)

Country Link
CN (1) CN106846812B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311930A (zh) * 2018-12-12 2020-06-19 阿里巴巴集团控股有限公司 一种获取车流量的方法及其装置
CN111724603A (zh) * 2020-07-01 2020-09-29 中南大学 基于交通轨迹数据的cav状态判定方法、装置、设备及介质
CN113936453A (zh) * 2021-09-09 2022-01-14 上海宝康电子控制工程有限公司 一种基于车头时距的信息识别方法及系统
CN115019508A (zh) * 2022-06-13 2022-09-06 华南理工大学 基于机器学习的道路监控车流仿真方法、装置、设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005062944A (ja) * 2003-08-12 2005-03-10 Toshiba Corp 交通流ミクロシミュレーションにおける車線変更判定方法およびそれを適用した交通流ミクロシミュレーションシステム
CN102314771A (zh) * 2011-04-27 2012-01-11 天津市市政工程设计研究院 基于平顺性分析确定快速路出入口最小间距的方法
CN102819954A (zh) * 2012-08-28 2012-12-12 南京大学 交通区域动态地图监控预测系统
CN103395419A (zh) * 2013-08-22 2013-11-20 贵州大学 基于安全间距策略的车辆队列行驶控制系统及其控制方法
CN104978853A (zh) * 2014-04-01 2015-10-14 中国移动通信集团公司 一种道路交通安全评估方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005062944A (ja) * 2003-08-12 2005-03-10 Toshiba Corp 交通流ミクロシミュレーションにおける車線変更判定方法およびそれを適用した交通流ミクロシミュレーションシステム
CN102314771A (zh) * 2011-04-27 2012-01-11 天津市市政工程设计研究院 基于平顺性分析确定快速路出入口最小间距的方法
CN102819954A (zh) * 2012-08-28 2012-12-12 南京大学 交通区域动态地图监控预测系统
CN103395419A (zh) * 2013-08-22 2013-11-20 贵州大学 基于安全间距策略的车辆队列行驶控制系统及其控制方法
CN104978853A (zh) * 2014-04-01 2015-10-14 中国移动通信集团公司 一种道路交通安全评估方法及系统

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311930A (zh) * 2018-12-12 2020-06-19 阿里巴巴集团控股有限公司 一种获取车流量的方法及其装置
CN111724603A (zh) * 2020-07-01 2020-09-29 中南大学 基于交通轨迹数据的cav状态判定方法、装置、设备及介质
CN111724603B (zh) * 2020-07-01 2021-05-18 中南大学 基于交通轨迹数据的cav状态判定方法、装置、设备及介质
CN113936453A (zh) * 2021-09-09 2022-01-14 上海宝康电子控制工程有限公司 一种基于车头时距的信息识别方法及系统
CN115019508A (zh) * 2022-06-13 2022-09-06 华南理工大学 基于机器学习的道路监控车流仿真方法、装置、设备及介质
CN115019508B (zh) * 2022-06-13 2023-09-29 华南理工大学 基于机器学习的道路监控车流仿真方法、装置、设备及介质

Also Published As

Publication number Publication date
CN106846812B (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
Yu et al. An improved car-following model considering relative velocity fluctuation
CN106846812A (zh) 一种十字路口的交通流量评估方法
CN109808685B (zh) 一种基于危险评估的汽车预警自动避撞控制方法
Kamal et al. On board eco-driving system for varying road-traffic environments using model predictive control
CN104794265B (zh) 一种基于加速度信息的移动跟随模型设计方法
CN107031621B (zh) 一种汽车防撞方法
CN108583576A (zh) 一种车辆经济学车速前瞻优化方法
Suzuki et al. A new approach to green light optimal speed advisory (GLOSA) systems for high-density traffic flowe
CN105788363B (zh) 一种行车预警方法、装置及系统
CN105667506A (zh) 一种汽车动态紧急避撞控制方法
Guo et al. Driving simulation study on speed-change lanes of the multi-lane freeway interchange
CN115497315B (zh) 一种车路协同环境下的动态公交专用道节能优化控制方法
CN106530695A (zh) 基于车联网的城市干道车辆行程时间实时预测方法
Ro et al. A new safety distance calculation for rear-end collision avoidance
CN108257402A (zh) 一种速度控制方法、装置及计算机存储介质
Dong et al. A comparative study of energy-efficient driving strategy for connected internal combustion engine and electric vehicles at signalized intersections
Wu et al. A longitudinal minimum safety distance model based on driving intention and fuzzy reasoning
CN110176138B (zh) 一种路口级的主动交通诱导方法
CN112417601B (zh) 一种aebs测试时车辆与仿真行人协同控制方法
CN112542061B (zh) 基于车联网的借道超车控制方法、装置、系统及存储介质
Ying et al. Research on adaptive cruise control systems and performance analysis using Matlab and Carsim
Guo et al. Variable time headway autonomous emergency braking control algorithm based on model predictive control
Clement et al. Simple platoon advancement: a model of automated vehicle movement at signalised intersections
Huang et al. Study on vehicle cumulative energy consumption of signalized intersections
CN115817460A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190125

Termination date: 20200314

CF01 Termination of patent right due to non-payment of annual fee