CN102142195B - Method for acquiring driving condition information of urban bus rapid transit - Google Patents

Method for acquiring driving condition information of urban bus rapid transit Download PDF

Info

Publication number
CN102142195B
CN102142195B CN 201110064222 CN201110064222A CN102142195B CN 102142195 B CN102142195 B CN 102142195B CN 201110064222 CN201110064222 CN 201110064222 CN 201110064222 A CN201110064222 A CN 201110064222A CN 102142195 B CN102142195 B CN 102142195B
Authority
CN
China
Prior art keywords
time
dead time
vehicle
website
overall travel
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 201110064222
Other languages
Chinese (zh)
Other versions
CN102142195A (en
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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN 201110064222 priority Critical patent/CN102142195B/en
Publication of CN102142195A publication Critical patent/CN102142195A/en
Application granted granted Critical
Publication of CN102142195B publication Critical patent/CN102142195B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

The invention relates to a method for acquiring the driving condition information of an urban bus rapid transit, which comprises the following steps: firstly, selecting a BRT (Bus Rapid Transit) vehicle and carrying out information acquisition to the BRT vehicle through the sampling frequency of 1s, wherein the acquired initial information comprises the driving time ti, the longitude Xi, the latitude Yi and the instantaneous velocity Vi; secondly, defining the acceleration, deceleration, constant speed and idle speed driving conditions in the driving process of the vehicle; thirdly, converting the initial information; and fourthly, analyzing the vehicle information. The method has the beneficial effects that the working condition information of the vehicle can be accurately or precisely obtained, and more working condition information of the vehicle can be obtained to be convenient for understanding the working conditions of the vehicle.

Description

The acquisition methods of driving condition information of urban bus rapid transit
Technical field
The invention belongs to intelligent transportation field, particularly a kind of acquisition methods of driving condition information of urban bus rapid transit.
Background technology
In order effectively to alleviate Urban Traffic Jam Based and automotive emission pollution situation, energy-saving and emission-reduction, domestic increasing city is actively carried out the public traffic in priority strategy and is widelyd popularize fast public traffic system (Bus Rapid Transit is called for short BRT) in traffic policy.BRT is a kind of novel mass transit system between rail rapid transit and conventional public transport, it is a kind of large conveying quantity mode of transportation, usually also be called by the people " ground subway system ", it is to utilize modern public transport technology to cooperate intelligent transportation and operation management (such as the signal priority system), open up busway and build Novel bus station, the realization rail transportation operation is served, and reaches the urban transit system of a kind of uniqueness of light rail excellent service standard.
The domestic mainly total kilometrage by gathering at times BRT vehicle starting point and total working time, and then the operation conditions of BRT vehicle is grasped and estimated to the average running speed that obtains the different periods.Thereby there is a following problem:
1, can not accurately grasp the BRT vehicle in the stand-by period of each bus stop and signalized intersections, and each stand-by period is in total proportion relation working time;
2, lack the scientific basis that follow-up relevant intelligent transportation and operation management are further improved and optimize, such as the optimization of each intersection signal priority system etc.
Summary of the invention
For adopt existing method can't reflect the problem of BRT traveling state of vehicle comprehensively, the present invention proposes a kind of acquisition methods that can comprehensively reflect the driving condition information of urban bus rapid transit of BRT travel conditions.
Technical scheme of the present invention:
The acquisition methods of driving condition information of urban bus rapid transit is characterized in that: its step is as follows:
(1) a selected BRT vehicle carries out information acquisition with the sample frequency of 1s to the BRT vehicle, and the initial information of collection comprises that running time is t i, longitude is X i, latitude is Y i, instantaneous velocity is V i
(2) to the acceleration in the Vehicle Driving Cycle process, deceleration, at the uniform velocity reach the idling travel conditions and define; Accelerating mode refers in the process of moving acceleration a>0.1m/s of vehicle 2Decelerating mode, refer to vehicle in the process of moving acceleration a<-0.1m/s 2Operating mode at the uniform velocity, refer to vehicle in the process of moving-0.1m/s 2≤ a≤0.1m/s 2Idling operation refers to that vehicle motor is in running order but travel speed is 0;
(3) initial information is transformed, the instantaneous velocity after the definition process is v i=1000V i/ 3600, unit is m/s, acceleration a i=v i/ v I-1, unit is m/s 2
(4) information of vehicles is carried out analyzing and processing;
A, total kilometres are fixed values;
B, overall travel time are through collecting;
C, average running speed are drawn by average running speed=total kilometres/overall travel time formula, and unit is km/h;
D, acceleration time T aBe all stacks that meet the time point of accelerating mode definition in overall travel time, unit is s; The acceleration time ratio is by acceleration time ratio R a=acceleration time T a/ overall travel time formula draws;
E, deceleration time T jBe all stacks that meet the time point of decelerating mode definition in overall travel time, unit is s; Deceleration time, ratio was by ratio R deceleration time j=deceleration time T j/ overall travel time formula draws;
F, time T at the uniform velocity cBe that all meet the at the uniform velocity stack of the time point of operating mode definition in overall travel time, unit is s; At the uniform velocity time scale is by time scale R at the uniform velocity c=deceleration time T j/ overall travel time formula draws;
G, dead time T dFor all instantaneous velocity values in overall travel time are the stack of the time point of " 0 ", unit is s; The dead time ratio is by dead time ratio R d=dead time T d/ overall travel time formula draws;
H, average overall travel speed are drawn by average overall travel speed=total kilometres/(overall travel time-dead time) formula, and unit is km/h;
I, website P jThe on-board and off-board dead time
Figure BDA0000050580520000031
And dead time ratio
Figure BDA0000050580520000032
Website P jThe wait dead time that enters the station And dead time ratio
Figure BDA0000050580520000034
Crossing CR kDead time
Figure BDA0000050580520000035
And dead time ratio
Figure BDA0000050580520000036
The analyzing and processing process as follows:
Step 1: establishing website is P j, j=0,1,2...n, longitude and latitude is If the crossing is CR k, k=0,1,2...N, longitude and latitude is
Figure BDA0000050580520000038
Extract the idling point latitude and longitude coordinates longitude (X of BRT vehicle i, Y i), stop website with each The latitude and longitude coordinates of crossing
Figure BDA00000505805200000310
Find the solution respectively apart from d;
Step 2: if with idling point be crossing CR apart from the d minimum k, then the continuous dead time of this idling point formation is the BRT vehicle at crossing CR kThe dead time at place
Figure BDA00000505805200000311
If with idling point distance minimum for stopping website P j, then the continuous dead time of this idling point formation is the BRT vehicle and is stopping website P jThe dead time at place;
Step 3: extract next idling point, process by step 1 method, if with this idling point apart from the d minimum be the crossing, then the continuous dead time that forms of this idling point is this correspondence crossing dead time; If with this idling point apart from the d minimum for stopping website, and the idling point is judged to be when stopping website in the step 2, judge then whether the stop website in this stop website and the step 2 is same website, if same stop website, then the dead time in the step 2 belongs to the BRT vehicle at stop website P jThe wait dead time that enters the station
Figure BDA0000050580520000041
And the continuous dead time that this idling point forms is stop website P jThe on-board and off-board dead time If not same stop website, then the continuous dead time of this idling point formation is and stops website P J+1Dead time;
Step 4: after dead time is entered the station in the on-board and off-board dead time that obtains each website and wait, can determine respectively the dead time ratio With
Figure BDA0000050580520000044
After the dead time that obtains each crossing, can determine the dead time ratio of each crossing
Figure BDA0000050580520000045
Further, data acquisition in the step (1) is finished by acquisition system, described acquisition system comprises BRT vehicle, GPS terminal, signal receiving module, GPS information management system, satellite, and described GPS terminal and signal receiving module are installed in the front end of BRT vehicle; Described GPS terminal adopts the mode of high frequency sampling, and the sampling period is set as 1s; Described signal receiving module is finished real-time reception and the transmission of information; Described GPS information management system in real time to the BRT vehicle monitor, historical driving trace playback, electronic report forms statistical query; Described satellite gathers the longitude and latitude of BRT vehicle.
It is the GPS high frequency method of sampling in cycle that the present invention adopts 1s, BRT vehicle driving-cycle information is carried out data acquisition, proposed BRT vehicle each website on-board and off-board dead time, waited for the acquisition methods of the characteristic parameters such as dead time and each crossing dead time that enter the station.
Beneficial effect of the present invention: can be accurately or obtain accurately vehicle working condition information, and can obtain more vehicle working condition information, so that more understand the operating mode of vehicle.
Description of drawings
Fig. 1 is the structural drawing of acquisition system of the present invention.
Embodiment
The acquisition methods of driving condition information of urban bus rapid transit, its step is as follows:
(1) a selected BRT vehicle carries out information acquisition with the sample frequency of 1s to the BRT vehicle, and the initial information of collection comprises that running time is t i, longitude is X i, latitude is Y i, instantaneous velocity is V i
(2) to the acceleration in the Vehicle Driving Cycle process, deceleration, at the uniform velocity reach the idling travel conditions and define; Accelerating mode refers in the process of moving acceleration a>0.1m/s of vehicle 2Decelerating mode, refer to vehicle in the process of moving acceleration a<-0.1m/s 2Operating mode at the uniform velocity, refer to vehicle in the process of moving-0.1m/s 2≤ a≤0.1m/s 2Idling operation refers to that vehicle motor is in running order but travel speed is 0;
(3) initial information is transformed, the instantaneous velocity after the definition process is v i=1000V i/ 3600, unit is m/s, acceleration a i=v i/ v I-1, unit is m/s 2
(4) information of vehicles is carried out analyzing and processing;
A, total kilometres are fixed values;
B, overall travel time are through collecting;
C, average running speed are drawn by average running speed=total kilometres/overall travel time formula, and unit is km/h;
D, acceleration time T aBe all stacks that meet the time point of accelerating mode definition in overall travel time, unit is s; The acceleration time ratio is by acceleration time ratio R a=acceleration time T a/ overall travel time formula draws;
E, deceleration time T jBe all stacks that meet the time point of decelerating mode definition in overall travel time, unit is s; Deceleration time, ratio was by ratio R deceleration time j=deceleration time T j/ overall travel time formula draws;
F, time T at the uniform velocity cBe that all meet the at the uniform velocity stack of the time point of operating mode definition in overall travel time, unit is s; At the uniform velocity time scale is by time scale R at the uniform velocity c=deceleration time T j/ overall travel time formula draws;
G, dead time T dFor all instantaneous velocity values in overall travel time are the stack of the time point of " 0 ", unit is s; The dead time ratio is by dead time ratio R d=dead time T d/ overall travel time formula draws;
H, average overall travel speed are drawn by average overall travel speed=total kilometres/(overall travel time-dead time) formula, and unit is km/h;
I, website P jThe on-board and off-board dead time
Figure BDA0000050580520000061
And dead time ratio
Figure BDA0000050580520000062
Website P jThe wait dead time that enters the station
Figure BDA0000050580520000063
And dead time ratio
Figure BDA0000050580520000064
Crossing CR kDead time
Figure BDA0000050580520000065
And dead time ratio
Figure BDA0000050580520000066
The analyzing and processing process as follows:
Step 1: establishing website is P j, j=0,1,2...n, longitude and latitude is
Figure BDA0000050580520000067
If the crossing is CR k, k=0,1,2...N, longitude and latitude is
Figure BDA0000050580520000068
Extract the idling point latitude and longitude coordinates longitude (X of BRT vehicle i, Y i), stop website with each The latitude and longitude coordinates of crossing
Figure BDA00000505805200000610
Find the solution respectively apart from d;
Step 2: if with idling point be crossing CR apart from the d minimum k, then the continuous dead time of this idling point formation is the BRT vehicle at crossing CR kThe dead time at place
Figure BDA00000505805200000611
If with idling point distance minimum for stopping website P j, then the continuous dead time of this idling point formation is the BRT vehicle and is stopping website P jThe dead time at place;
Step 3: extract next idling point, process by step 1 method, if with this idling point apart from the d minimum be the crossing, then the continuous dead time that forms of this idling point is this correspondence crossing dead time; If with this idling point apart from the d minimum for stopping website, and the idling point is judged to be when stopping website in the step 2, judge then whether the stop website in this stop website and the step 2 is same website, if same stop website, then the dead time in the step 2 belongs to the BRT vehicle at stop website P jThe wait dead time that enters the station
Figure BDA0000050580520000071
And the continuous dead time that this idling point forms is stop website P jThe on-board and off-board dead time
Figure BDA0000050580520000072
If not same stop website, then the continuous dead time of this idling point formation is and stops website P J+1Dead time;
Step 4: after dead time is entered the station in the on-board and off-board dead time that obtains each website and wait, can determine respectively the dead time ratio
Figure BDA0000050580520000073
With
Figure BDA0000050580520000074
After the dead time that obtains each crossing, can determine the dead time ratio of each crossing
Figure BDA0000050580520000075
Data acquisition in the step (1) is finished by acquisition system, described acquisition system comprises BRT vehicle 1, GPS terminal 3, signal receiving module 2, GPS information management system 5, satellite 4, and described GPS terminal 3 and signal receiving module 2 are installed in the front end of BRT vehicle 1; Described GPS terminal 3 adopts the mode of high frequency sampling, and the sampling period is set as 1s; Described signal receiving module 2 is finished real-time reception and the transmission of information; Described GPS information management system 5 in real time to BRT vehicle 1 monitor, historical driving trace playback, electronic report forms statistical query; The longitude and latitude of 4 pairs of BRT vehicles 1 of described satellite gathers, and sees Fig. 1.
The described content of this instructions embodiment only is enumerating the way of realization of inventive concept; protection scope of the present invention should not be regarded as only limiting to the concrete form that embodiment states, protection scope of the present invention also reaches in those skilled in the art conceives the equivalent technologies means that can expect according to the present invention.

Claims (2)

1. the acquisition methods of driving condition information of urban bus rapid transit, it is characterized in that: its step is as follows:
(1) a selected BRT vehicle carries out information acquisition with the sample frequency of 1s to the BRT vehicle, and the initial information of collection comprises that running time is t i, longitude is X i, latitude is Y i, instantaneous velocity is V i
(2) to the acceleration in the Vehicle Driving Cycle process, deceleration, at the uniform velocity reach the idling travel conditions and define; Accelerating mode refers in the process of moving acceleration a of vehicle〉0.1m/s 2Decelerating mode, refer to vehicle in the process of moving acceleration a<-0.1m/s 2Operating mode at the uniform velocity, refer to vehicle in the process of moving-0.1m/s 2<=a<=0.1m/s 2Idling operation refers to that vehicle motor is in running order but travel speed is 0;
(3) the instantaneous velocity V of initial acquisition iUnit is km/h, and the instantaneous velocity after the definition process is v i=1000V i/ 3600, unit is m/s; So acceleration a i=v i-v I-1, unit is m/s 2, a wherein iBe t iAcceleration constantly, v iBe t iInstantaneous velocity constantly, v I-1Be t I-1Instantaneous velocity constantly;
(4) information of vehicles is carried out analyzing and processing;
A, total kilometres are fixed values;
B, overall travel time are through collecting;
C, average running speed are drawn by average running speed=total kilometres/overall travel time formula, and unit is km/h;
D, acceleration time T aBe all stacks that meet the time point of accelerating mode definition in overall travel time, unit is s; The acceleration time ratio is by acceleration time ratio R a=acceleration time T a/ overall travel time formula draws;
E, deceleration time T jBe all stacks that meet the time point of decelerating mode definition in overall travel time, unit is s; Deceleration time, ratio was by ratio R deceleration time j=deceleration time T j/ overall travel time formula draws;
F, time T at the uniform velocity cBe that all meet the at the uniform velocity stack of the time point of operating mode definition in overall travel time, unit is s; At the uniform velocity time scale is by time scale R at the uniform velocity c=time T at the uniform velocity c/ overall travel time formula draws;
G, dead time T dFor all instantaneous velocity values in overall travel time are the stack of the time point of " 0 ", unit is s; The dead time ratio is by dead time ratio R d=dead time T d/ overall travel time formula draws;
H, average overall travel speed are drawn by average overall travel speed=total kilometres/(overall travel time-dead time) formula, and unit is km/h;
I, website P jThe on-board and off-board dead time
Figure FDA00002083678100021
And dead time ratio
Figure FDA00002083678100022
Website P jThe wait dead time that enters the station
Figure FDA00002083678100023
And dead time ratio
Figure FDA00002083678100024
Crossing CR kDead time
Figure FDA00002083678100025
And dead time ratio
Figure FDA00002083678100026
The analyzing and processing process as follows:
Step 1: establishing website is P j, j=0,1,2...n, longitude and latitude is
Figure FDA00002083678100027
If the crossing is CR k, k=0,1,2...N, longitude and latitude is
Figure FDA00002083678100028
Extract the idling point latitude and longitude coordinates longitude (X of BRT vehicle i, Y i), stop website with each The latitude and longitude coordinates of crossing
Figure FDA000020836781000210
Find the solution respectively apart from d;
Step 2: if with idling point be crossing CR apart from the d minimum k, then the continuous dead time of this idling point formation is the BRT vehicle at crossing CR kThe dead time at place If with idling point distance minimum for stopping website P j, then the continuous dead time of this idling point formation is the BRT vehicle and is stopping website P jThe dead time at place;
Step 3: extract next idling point, process by step 1 method, if with this idling point apart from the d minimum be the crossing, then the continuous dead time that forms of this idling point is this correspondence crossing dead time; If with this idling point apart from the d minimum for stopping website, and the idling point is judged to be when stopping website in the step 2, judge then whether the stop website in this stop website and the step 2 is same website, if same stop website, then the dead time in the step 2 belongs to the BRT vehicle at stop website P jThe wait dead time that enters the station And the continuous dead time that this idling point forms is stop website P jThe on-board and off-board dead time
Figure FDA00002083678100032
If not same stop website, then the continuous dead time of this idling point formation is and stops website P J+1Dead time;
Step 4: after dead time is entered the station in the on-board and off-board dead time that obtains each website and wait, can determine respectively the dead time ratio
Figure FDA00002083678100033
With
Figure FDA00002083678100034
After the dead time that obtains each crossing, can determine the dead time ratio of each crossing
Figure FDA00002083678100035
2. the acquisition methods of driving condition information of urban bus rapid transit according to claim 1, it is characterized in that: the data acquisition in the step (1) is finished by acquisition system, described acquisition system comprises BRT vehicle, GPS terminal, signal receiving module, GPS information management system, satellite, and described GPS terminal and signal receiving module are installed in the front end of BRT vehicle; Described GPS terminal adopts the mode of high frequency sampling, and the sampling period is set as 1s; Described signal receiving module is finished real-time reception and the transmission of information; Described GPS information management system in real time to the BRT vehicle monitor, historical driving trace playback, electronic report forms statistical query; Described satellite gathers the longitude and latitude of BRT vehicle.
CN 201110064222 2011-03-16 2011-03-16 Method for acquiring driving condition information of urban bus rapid transit Expired - Fee Related CN102142195B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110064222 CN102142195B (en) 2011-03-16 2011-03-16 Method for acquiring driving condition information of urban bus rapid transit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110064222 CN102142195B (en) 2011-03-16 2011-03-16 Method for acquiring driving condition information of urban bus rapid transit

Publications (2)

Publication Number Publication Date
CN102142195A CN102142195A (en) 2011-08-03
CN102142195B true CN102142195B (en) 2013-01-16

Family

ID=44409673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110064222 Expired - Fee Related CN102142195B (en) 2011-03-16 2011-03-16 Method for acquiring driving condition information of urban bus rapid transit

Country Status (1)

Country Link
CN (1) CN102142195B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104875731B (en) * 2015-05-19 2017-05-17 交通运输部公路科学研究所 Method for identifying rapid acceleration or rapid deceleration of vehicle in real time by using satellite positioning data
CN107133374A (en) * 2017-01-23 2017-09-05 长安大学 A kind of construction method of mode operating mode
CN106997662B (en) * 2017-04-26 2020-04-10 长安大学 Urban bus working condition construction method
CN107527140A (en) * 2017-07-28 2017-12-29 西安理工大学 A kind of bullet train operating condition based on fuzzy membership determines method
CN113064933B (en) * 2021-03-22 2023-04-07 四川云智造科技有限公司 Vehicle weight prediction method and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101256717A (en) * 2008-04-09 2008-09-03 重庆市易博数字技术有限公司 Quick public transport vehicle positioning based on GPS and traffic signal control method and apparatus
CN101751788A (en) * 2009-12-31 2010-06-23 浙江工业大学 Public transit exclusive lane resource sharing method and system for guaranteeing bus rapid transit priority

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008046998A (en) * 2006-08-19 2008-02-28 Masahiro Watanabe Route bus operation management method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101256717A (en) * 2008-04-09 2008-09-03 重庆市易博数字技术有限公司 Quick public transport vehicle positioning based on GPS and traffic signal control method and apparatus
CN101751788A (en) * 2009-12-31 2010-06-23 浙江工业大学 Public transit exclusive lane resource sharing method and system for guaranteeing bus rapid transit priority

Also Published As

Publication number Publication date
CN102142195A (en) 2011-08-03

Similar Documents

Publication Publication Date Title
CN103258430B (en) Road traveling time calculating and traffic road condition judging method and road traveling time calculating and traffic road condition judging device
CN102142195B (en) Method for acquiring driving condition information of urban bus rapid transit
CN102184641B (en) Running information based road condition management method and system
CN108198425A (en) A kind of construction method of Electric Vehicles Driving Cycle
CN103106794B (en) Road traffic energy-saving method based on probe vehicle technology
CN102855760B (en) On-line queuing length detection method based on floating vehicle data
CN107204115B (en) Method for estimating parking rate of main line traffic flow
CN101577049B (en) Map-matching method based on filtering of parking point data of floating vehicles
CN103761876A (en) Vehicle-road cooperation-based method for traffic information collection and status evaluation
CN102708689B (en) Real-time traffic monitoring system
CN102722984B (en) Real-time road condition monitoring method
CN104575035A (en) Intersection self-adaptation control method based on car networking environment
CN104318770A (en) Method for detecting traffic jam state of expressway in real time based on mobile phone data
CN103198713A (en) Traffic accident reduction vehicle regulation and control method based on traffic data and weather data
CN202472944U (en) Urban traffic information acquisition and processing system based on integration of data of multiple sensors
CN103786733A (en) Environment-friendly driving behavior prompting method for automatic transmission automobile
CN105046983A (en) Traffic flow prediction system and method based on vehicle-road cooperation
CN105427638A (en) Man-machine interactive vehicle-speed control and early-warning system for expressway construction operation area
CN109067857A (en) A kind of car networking big data analysis intelligent optimization algorithm
CN105575140A (en) City traffic control method based on RFID electronic labels
CN105427604A (en) Realization method for manifesting public-transportation smooth pass index
CN104346953A (en) Dynamic parking system based on cloud computing
CN106683406B (en) A kind of current bottleneck detection method of the public transportation lane based on public transport vehicle-mounted GPS data
CN109816183A (en) A kind of optimization method and device of accurate bus passenger flow big data
CN103900553A (en) Regional map recording method and system

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: 20130116

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