CN102270387B - 一种公路交通疏导方法 - Google Patents

一种公路交通疏导方法 Download PDF

Info

Publication number
CN102270387B
CN102270387B CN 201110084057 CN201110084057A CN102270387B CN 102270387 B CN102270387 B CN 102270387B CN 201110084057 CN201110084057 CN 201110084057 CN 201110084057 A CN201110084057 A CN 201110084057A CN 102270387 B CN102270387 B CN 102270387B
Authority
CN
China
Prior art keywords
traffic
scope
dredge
dredging
starting point
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.)
Expired - Fee Related
Application number
CN 201110084057
Other languages
English (en)
Other versions
CN102270387A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN 201110084057 priority Critical patent/CN102270387B/zh
Publication of CN102270387A publication Critical patent/CN102270387A/zh
Application granted granted Critical
Publication of CN102270387B publication Critical patent/CN102270387B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

一种公路交通疏导方法:当发现交通意外事件后,确定事发点N和事件引起的排队尾端位置M,与M直接相连的上游交叉口集合O1,与N直接相连的下游交叉口集合D1,交叉口集合Oi为与交叉口集合Oi-1中交叉口直接相连的上游交叉口集合,交叉口集合Dj为与交叉口集合Dj-1中交叉口直接相连的下游交叉口集合,疏导范围表示为
Figure DDA0000053807070000011
当n=1,m=1时,确定初始疏导范围(O1,D1);采用基于增量分配的交通疏导方法,在疏导范围内对交通进行疏导;对疏导结果合理性判断:若疏导范围内各疏导路段的交通流量未超过其实际通行能力,则疏导结果合理,生成疏导方案,否则,扩大疏导范围,然后采用基于增量分配的交通疏导方法,在扩大后的疏导范围内对交通进行疏导。

Description

一种公路交通疏导方法
技术领域
本发明涉及在公路上发生突发交通事件导致拥堵排队后,根据事件引发的拥堵排队长度、事发路段周边路网布局及交通流量情况,确定交通疏导范围以及该范围内的疏导方案,为公路突发交通事件的处理、交通处置措施的实施提供依据。属于公路交通控制领域。
背景技术
公路上发生突发交通事件后,将严重影响交通的正常运行,需要及时对未进入拥堵路段的交通量进行合理疏导,避免交通拥堵恶化,缓解事件造成的不利影响。为了保证疏导的合理性,需要确定一种实用的疏导方法,生成有效的疏导方案。国内外学者在突发交通事件交通疏导方面进行了一些研究,但现有的研究存在一些不足:(1)已有的疏导模型或疏导方法多数针对于城市突发交通事件,由于公路网络结构及流量特点不同于城市道路网,已有疏导方法不能适用于公路突发交通事件下的疏导;(2)已有的交通疏导方法,建立的疏导目标基本围绕于网络费用最小、总时间最小或者可达性较好,疏导方案产生速度慢、效率低,而公路突发交通事件的处理对时间要求比较紧迫,所以应以快速疏导为目标,确定疏导方法。
经发明人长期研究发现,根据事件影响程度以及周边路段流量情况,可以合理确定交通疏导范围,确定各疏导起点疏导交通量,并结合改进的增量分配法,实现对拥堵排队点上游的交通量实施疏导,避免更大程度的拥堵排队,保证公路网最大程度地畅通,减少突发交通事件对公路网造成的负面影响,保障事件的快速处置。
发明内容
技术问题:
本发明的主要目的是提供一种公路交通疏导方法,该方法与突发交通事件影响程度及公路交通特性相适应,能起到及时有效缓解拥堵道路交通压力,避免影响周边道路上交通正常运行的作用。
技术方案:为达到上述目的,本发明提出的一种公路交通疏导方法,通过路网数据的采集,确定初始疏导范围,在初始疏导范围内按照所述交通疏导方法进行交通疏导,若疏导结果合理,则生成疏导方案;否则扩大疏导范围,在扩大后的疏导范围内再次进行交通疏导;以此循环,直至生成疏导方案。
具体疏导步骤为:
步骤(1)初始化:当发现交通意外事件后,确定事件发生点N和事件引起的排队尾端位置M,与M直接相连的上游交叉口集合O1,与N直接相连的下游交叉口集合D1,定义M上游的第i个交叉口集合Oi为与M上游的第i-1个交叉口集合Oi-1中交叉口直接相连的上游交叉口集合,其中i≥2且为正整数,并将Oi定义为第i级疏导起点集合;定义N下游的第j个交叉口集合Dj为与N下游的第j-1个交叉口集合Dj-1中交叉口直接相连的下游交叉口集合,其中j≥2且为正整数,并将Dj定义为第j级疏导讫点集合,疏导范围则表示为
Figure GDA00003014380900021
其中n=1,2,3……,m=1,2,3……,当n=1,m=1时,确定初始疏导范围(O1,D1),表示初始疏导范围内包含第1级疏导起点、第1级疏导讫点及疏导起点与疏导讫点之间的路段,
步骤(2)采用基于增量分配的交通疏导方法,在疏导范围内对交通进行疏导,
步骤(3)疏导结果合理性分析:若疏导范围内各疏导路段的交通流量未超过其实际通行能力,则疏导结果合理,生成疏导方案;否则,转至步骤(4),
步骤(4)扩大疏导范围:在疏导范围中,若n=m,则寻找与On直接相连的上游交叉口集合On+1,扩大后的疏导范围为
Figure GDA00003014380900023
转至步骤(2);若n≠m,则寻找与Dm直接相连的下游交叉口集合Dm+1,扩大后的疏导范围为
Figure GDA00003014380900031
转至步骤(2),
所述的基于增量分配的交通疏导方法为:
步骤①确定疏导起点疏导量:起点的疏导总量q=Q-C',Q为事发路段M→N方向,在正常交通的情况下的交通量,
Figure GDA00003014380900032
其中Qk为正常交通情况下,疏导范围内各疏导起点流经事发路段的交通量,k=1,2,3……p,p为疏导范围内所有疏导起点的总数;C'为根据事发路段检测设备,实时采集的事发路段M→N方向的实际通行能力,事发后每一个疏导起点需要疏导的交通量 Δ Q k = Q k × q Q ,
步骤②疏导范围内各路段流量、路段行驶时间及路径行驶时间的初始化:将每一个疏导起点需要疏导的交通量ΔQk从疏导范围
Figure GDA00003014380900034
内的含有事发路段交通量的路段上的正常交通流量中剔除,确定疏导范围
Figure GDA00003014380900035
内的疏导前的各路段初始化流量qr,其中r为疏导范围内的任意一条路段;将每一个疏导起点需要疏导的交通量ΔQk平均分为10份,即
Figure GDA00003014380900036
根据疏导前经过初始化的各路段交通量qr、路段实际通行能力cr、路段自由流出行时间tr 0,计算每个路段初始化行驶时间并计算起讫点路径行驶时间
Figure GDA00003014380900038
其中,k表示第k个疏导起点,Rk表示第k个疏导起点的任一路径内所有路段的集合,
步骤③交通分配:以行驶时间最短为目标,将每一个疏导起点的交通量Δqk分配到行驶时间最短的路径上,并更新各路段交通量qr
步骤④更新各起讫点的路径行驶时间:根据更新后的各路段交通量qr、路段通行能力cr、路段自由流出行时间tr 0,更新每个路段的行驶时间 t r = t r 0 [ 1 + 0.15 ( q r c r ) 4 ] , 并重新计算起讫点路径行驶时间 T k = Σ r ∈ R k t r ,
⑤判定分配是否结束:如果交通分配次数已达10,结束分配;否则,转至步骤③。
有益效果:
本发明根据公路突发交通事件影响程度、事发路段周边路网流量状况,对拥堵排队点上游未到达拥堵路段交通量在周边路网上进行合理疏导,保证整个路网的畅通,同时为交通事件的处置提供更多的时间与空间支持。
1、提出了疏导范围扩大的方法,可以适应各种影响程度下的交通疏导,避免了受突发交通事件影响程度的制约,能够解决更大范围公路网的交通疏导问题,为疏导范围的确定提供了有效保障。
2、给出了各疏导起点具体疏导交通量的确定方法,为进行交通疏导提供了可靠的数据保障,避免了人为判断误差的出现;
3、将传统增量分配法步骤引入到疏导过程中,同时为与疏导时路段及流量参数的特殊性相适应,对其进行了修正,避免了增量分配法在疏导中的不适性。
附图说明
图1公路交通疏导方法流程图;
图2以增量分配法为核心进行交通疏导的流程图;
图3公路网示意图。
具体实施方式
步骤(1)初始化:当发现交通意外事件后,确定事件发生点N和事件引起的排队尾端位置M,与M直接相连的上游交叉口集合O1,与N直接相连的下游交叉口集合D1,定义M上游的第i个交叉口集合Oi为与M上游的第i-1个交叉口集合Oi-1中交叉口直接相连的上游交叉口集合,其中i≥2且为正整数,并将Oi定义为第i级疏导起点集合;定义N下游的第j个交叉口集合Dj为与N下游的第j-1个交叉口集合Dj-1中交叉口直接相连的下游交叉口集合,其中j≥2且为正整数,并将Dj定义为第j级疏导讫点集合,疏导范围则表示为
Figure GDA00003014380900051
其中n=1,2,3……,m=1,2,3……,当n=1,m=1时,确定初始疏导范围(O1,D1),表示初始疏导范围内包含第1级疏导起点、第1级疏导讫点及疏导起点与疏导讫点之间的路段,
步骤(2)采用基于增量分配的交通疏导方法,在疏导范围内对交通进行疏导,
步骤(3)疏导结果合理性分析:若疏导范围内各疏导路段的交通流量未超过其实际通行能力,则疏导结果合理,生成疏导方案;否则,转至步骤(4),其中公路路段实际通行能力可采用交通届公认的计算方法求解获得,参见徐吉谦主编的《交通工程总论》:75-76页,
步骤(4)扩大疏导范围:在疏导范围
Figure GDA00003014380900052
中,若n=m,则寻找与On直接相连的上游交叉口集合On+1,扩大后的疏导范围为
Figure GDA00003014380900053
转至步骤(2);若n≠m,则寻找与Dm直接相连的下游交叉口集合Dm+1,扩大后的疏导范围为
Figure GDA00003014380900054
转至步骤(2),
所述的基于增量分配的交通疏导方法为:
步骤①确定疏导起点疏导量:起点的疏导总量q=Q-C',Q为事发路段M→N方向,在正常交通的情况下的交通量,
Figure GDA00003014380900055
其中Qk为正常交通情况下,疏导范围内各疏导起点流经事发路段的交通量,k=1,2,3……p,p为疏导范围内所有疏导起点的总数;C'为根据事发路段检测设备,实时采集的事发路段M→N方向的实际通行能力,事发后每一个疏导起点需要疏导的交通量 Δ Q k = Q k × q Q ,
步骤②疏导范围内各路段流量、路段行驶时间及路径行驶时间的初始化:将每一个疏导起点需要疏导的交通量ΔQk从疏导范围
Figure GDA00003014380900061
内的含有事发路段交通量的路段上的正常交通流量中剔除,确定疏导范围内的疏导前的各路段初始化流量qr,其中r为疏导范围内的任意一条路段;将每一个疏导起点需要疏导的交通量ΔQk平均分为10份,即
Figure GDA00003014380900063
根据疏导前经过初始化的各路段交通量qr、路段实际通行能力cr、路段自由流出行时间tr 0,计算每个路段初始化行驶时间
Figure GDA00003014380900064
并计算起讫点路径行驶时间
Figure GDA00003014380900065
其中,k表示第k个疏导起点,Rk表示第k个疏导起点的任一路径内所有路段的集合,
步骤③交通分配:以行驶时间最短为目标,将每一个疏导起点的交通量Δqk分配到行驶时间最短的路径上,并更新各路段交通量qr
步骤④更新各起讫点的路径行驶时间:根据更新后的各路段交通量qr、路段通行能力cr、路段自由流出行时间tr 0,更新每个路段的行驶时间 t r = t r 0 [ 1 + 0.15 ( q r c r ) 4 ] , 并重新计算起讫点路径行驶时间 T k = Σ r ∈ R k t r ,
⑤判定分配是否结束:如果交通分配次数已达10,结束分配;否则,转至步骤③。
结合附图,对本发明做进一步说明:
(1)如附图3所示,当发现突发交通事件发生在o1d1间,影响交通流方向为o1→d1时,确定初始的疏导范围为上游交叉口集合{o1}、下游交叉口集合{d1}和{o1}与{d1}之间的路径o1a0d1、o1b0d1
(2)根据历史数据,采集事发路段o1→d1方向,在正常交通情况下的交通量为6660pcu/h,事发后的实际通行能力为4400pcu/h,计算需要疏导的总交通量q=6660pcu/h–4400pcu/h=2260pcu/h。因此,疏导起点o1需要疏导交通量为2260pcu/h,并将疏导起点o1需要疏导交通量平均分为10份,即Δqk=226pcu/h。
(3)采集疏导前各路段初始化流量qr、路段实际通行能力cr、路段自由流出行时间tr 0,如附图3所示,计算每个路段初始化行驶时间
Figure GDA00003014380900071
并计算两条路径o1a0d1、o1b0d1的行驶时间。以路径行驶时间最短为目标,将疏导起点o1的交通量Δqk分配到最短路径上。分配后,更新各路段交通量qr。并根据路段实际通行能力cr、路段自由流出行时间tr 0,更新每个路段的行驶时间
Figure GDA00003014380900072
计算两条路径o1a0d1、o1b0d1的行驶时间。
(4)如果交通分配次数等于10,结束分配,分配结果如见表1:
表1疏导范围下的疏导结果
Figure GDA00003014380900073
根据表1的疏导结果可知:参与疏导路段o1a0与b0d1的最终流量均大于两者的实际通行能力。所以疏导结果不合理,需要扩大疏导范围,重新进行疏导。
(5)扩大的疏导范围为上游交叉口集合{o1,o2,a1,b1,}、下游交叉口集合{d1}及其之间的14条路径,见表2。通过上述步骤在新的疏导范围内对交通流重新分配,分配结果如表2所示:
表2扩大的疏导范围内交通流疏导结果
Figure GDA00003014380900074
Figure GDA00003014380900081
根据表2的疏导结果可知:参与疏导路段的交通量均小于路段实际通行能力,疏导结果合理,生成疏导方案。

Claims (1)

1.一种公路交通疏导方法,其特征在于,所述方法的具体步骤为:
步骤(1)初始化:当发现交通意外事件后,确定事件发生点N和事件引起的排队尾端位置M,与M直接相连的上游交叉口集合O1,与N直接相连的下游交叉口集合D1,定义M上游的第i个交叉口集合Oi为与M上游的第i-1个交叉口集合Oi-1中交叉口直接相连的上游交叉口集合,其中i≥2且为正整数,并将Oi定义为第i级疏导起点集合;定义N下游的第j个交叉口集合Dj为与N下游的第j-1个交叉口集合Dj-1中交叉口直接相连的下游交叉口集合,其中j≥2且为正整数,并将Dj定义为第j级疏导讫点集合,疏导范围则表示为其中n=1,2,3......,m=1,2,3......,当n=1,m=1时,确定初始疏导范围(O1,D1),表示初始疏导范围内包含第1级疏导起点、第1级疏导讫点及疏导起点与疏导讫点之间的路段,
步骤(2)采用基于增量分配的交通疏导方法,在疏导范围内对交通进行疏导,
步骤(3)疏导结果合理性分析:若疏导范围内各疏导路段的交通流量未超过其实际通行能力,则疏导结果合理,生成疏导方案;否则,转至步骤(4),
步骤(4)扩大疏导范围:在疏导范围
Figure FDA0000053807040000012
中,若n=m,则寻找与On直接相连的上游交叉口集合On+1,扩大后的疏导范围为
Figure FDA0000053807040000013
转至步骤(2);若n≠m,则寻找与Dm直接相连的下游交叉口集合Dm+1,扩大后的疏导范围为
Figure FDA0000053807040000014
转至步骤(2),
所述的基于增量分配的交通疏导方法为:
步骤①确定疏导起点疏导量:起点的疏导总量q=Q-C′,Q为事发路段M→N方向,在正常交通的情况下的交通量,
Figure FDA0000053807040000015
其中Qk为正常交通情况下,疏导范围内各疏导起点流经事发路段的交通量,k=1,2,3……p,p为疏导范围内所有疏导起点的总数;C′为根据事发路段检测设备,实时采集的事发路段M→N方向的实际通行能力,事发后每一个疏导起点需要疏导的交通量
Figure FDA0000053807040000021
步骤②疏导范围内各路段流量、路段行驶时间及路径行驶时间的初始化:将每一个疏导起点需要疏导的交通量ΔQk从疏导范围
Figure FDA0000053807040000022
内的含有事发路段交通量的路段上的正常交通流量中剔除,确定疏导范围
Figure FDA0000053807040000023
内的疏导前的各路段初始化流量qr,其中r为疏导范围内的任意一条路段;将每一个疏导起点需要疏导的交通量ΔQk平均分为10份,即
Figure FDA0000053807040000024
根据疏导前经过初始化的各路段交通量qr、路段实际通行能力cr、路段自由流出行时间tr 0,计算每个路段初始化行驶时间
Figure FDA0000053807040000025
并计算起讫点路径行驶时间
Figure FDA0000053807040000026
其中,k表示第k个疏导起点,Rk表示第k个疏导起点的任一路径内所有路段的集合,
步骤③交通分配:以行驶时间最短为目标,将每一个疏导起点的交通量Δqk分配到行驶时间最短的路径上,并更新各路段交通量qr
步骤④更新各起讫点的路径行驶时间:根据更新后的各路段交通量qr、路段通行能力cr、路段自由流出行时间tr 0,更新每个路段的行驶时间
Figure FDA0000053807040000027
并重新计算起讫点路径行驶时间
Figure FDA0000053807040000028
⑤判定分配是否结束:如果交通分配次数已达10,结束分配;否则,转至步骤③。
CN 201110084057 2011-04-02 2011-04-02 一种公路交通疏导方法 Expired - Fee Related CN102270387B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110084057 CN102270387B (zh) 2011-04-02 2011-04-02 一种公路交通疏导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110084057 CN102270387B (zh) 2011-04-02 2011-04-02 一种公路交通疏导方法

Publications (2)

Publication Number Publication Date
CN102270387A CN102270387A (zh) 2011-12-07
CN102270387B true CN102270387B (zh) 2013-10-16

Family

ID=45052674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110084057 Expired - Fee Related CN102270387B (zh) 2011-04-02 2011-04-02 一种公路交通疏导方法

Country Status (1)

Country Link
CN (1) CN102270387B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105023424B (zh) * 2014-04-24 2017-08-08 深圳市赛格导航科技股份有限公司 一种道路交通拥堵控制方法及控制系统
CN106504521B (zh) * 2016-11-30 2018-10-30 东南大学 一种大型活动交通疏导方法及系统
CN108682147A (zh) * 2018-05-03 2018-10-19 广东省高速公路营运管理协会 一种高速公路交通拥堵疏导决策方法
CN111275957A (zh) * 2018-12-05 2020-06-12 杭州海康威视系统技术有限公司 一种交通事故信息采集方法、系统及摄像机
CN110276975B (zh) * 2019-07-10 2021-05-11 中国联合网络通信集团有限公司 一种自动驾驶车辆分流调度方法及系统
CN115359671A (zh) * 2022-07-11 2022-11-18 东风汽车集团股份有限公司 一种路口车辆协同控制方法及相关设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001023835A3 (en) * 1999-09-27 2001-11-15 Decell Inc Apparatus and methods for providing route guidance for vehicles
CN101159098A (zh) * 2007-11-20 2008-04-09 东南大学 城市交通意外事件的自动检测确认方法
WO2010077225A2 (en) * 2008-12-30 2010-07-08 Tele Atlas North America, Inc. A method and system for transmitting and/or receiving at least one location reference, enhanced by at least one focusing factor

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001023835A3 (en) * 1999-09-27 2001-11-15 Decell Inc Apparatus and methods for providing route guidance for vehicles
CN101159098A (zh) * 2007-11-20 2008-04-09 东南大学 城市交通意外事件的自动检测确认方法
WO2010077225A2 (en) * 2008-12-30 2010-07-08 Tele Atlas North America, Inc. A method and system for transmitting and/or receiving at least one location reference, enhanced by at least one focusing factor

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
李文华 等.高速公路突发交通事件疏导措施决策.《公路交通科技》.2010,第27卷(第11期),
胡晓健 等.道路交通意外事件影响范围确定方法.《东南大学学报》.2007,第37卷(第5期),
道路交通意外事件影响范围确定方法;胡晓健 等;《东南大学学报》;20070930;第37卷(第5期);第934-939页 *
高速公路突发交通事件疏导措施决策;李文华 等;《公路交通科技》;20101130;第27卷(第11期);第149-153页 *

Also Published As

Publication number Publication date
CN102270387A (zh) 2011-12-07

Similar Documents

Publication Publication Date Title
CN102270387B (zh) 一种公路交通疏导方法
CN102496264B (zh) 一种高速公路突发交通事件影响范围的确定方法
CN105279981A (zh) 一种潮汐车道行驶方向动态控制方法
CN102254425A (zh) 修正速度的宏观交通流离散模型稳定性建模方法
CN102592440A (zh) 道路网关键节点诊断技术
CN102254423A (zh) 宏观交通流离散模型稳定性建模方法
CN113159994A (zh) 多泥沙河流的船闸淤积处理方法及装置
CN103456181A (zh) 一种改进的multiband干线协调控制方法
CN105788238A (zh) 基于量子门和自适应控制的类弹簧车辆跟驰模型建立方法
Lee et al. Evaluation of semi-actuated signals and pedestrian push buttons using a microscopic traffic simulation model
CN102930727B (zh) 基于FPGA及改进Ross模型的在线交通瓶颈预测控制方法
CN102945609B (zh) 基于FPGA及改进Papageorgiou-E模型的在线交通瓶颈预测控制方法
CN106504521B (zh) 一种大型活动交通疏导方法及系统
CN205507885U (zh) 一种基于gprs网络的灌区用水计量、墒情监测及科学用水管理系统
CN102938209B (zh) 基于fpga及改进离散宏观d模型的在线交通瓶颈预测控制方法
CN102945611B (zh) 基于fpga及改进离散宏观p模型的在线交通瓶颈预测控制方法
CN102930728B (zh) 基于fpga及改进吴正模型的在线交通瓶颈预测控制方法
Ahmed et al. Modeling of delay induced by downstream traffic disturbances at signalized intersections
CN205276149U (zh) 一种用于高速公路短期保通便道的涵洞
Shen et al. A mesoscopic model for bicycle flow
Jiahui et al. Research on optimization of traffic management and control measures for expressway emergency based on VISSIM simulation
CN102945606B (zh) 基于FPGA及改进Kerner-Konhauser模型的在线交通瓶颈预测控制方法
CN102938207B (zh) 基于fpga及改进lwr模型的在线交通瓶颈预测控制方法
CN102938205B (zh) 基于FPGA及改进Xue-Dai模型的在线交通瓶颈预测控制方法
CN102945607B (zh) 基于FPGA及改进Aw-Rascle模型的在线交通瓶颈预测控制方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131016

Termination date: 20160402

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