CN111090929B - Heavy haul railway throughput computing method - Google Patents

Heavy haul railway throughput computing method Download PDF

Info

Publication number
CN111090929B
CN111090929B CN201911146733.5A CN201911146733A CN111090929B CN 111090929 B CN111090929 B CN 111090929B CN 201911146733 A CN201911146733 A CN 201911146733A CN 111090929 B CN111090929 B CN 111090929B
Authority
CN
China
Prior art keywords
time
skylight
train
railway
data
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
Application number
CN201911146733.5A
Other languages
Chinese (zh)
Other versions
CN111090929A (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.)
China Railway Design Corp
Original Assignee
China Railway Design Corp
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 China Railway Design Corp filed Critical China Railway Design Corp
Priority to CN201911146733.5A priority Critical patent/CN111090929B/en
Publication of CN111090929A publication Critical patent/CN111090929A/en
Application granted granted Critical
Publication of CN111090929B publication Critical patent/CN111090929B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/20Administration of product repair or maintenance
    • G06Q50/40

Abstract

The invention discloses a calculation method for the passing capacity of a heavy haul railway, which comprises the following steps: loading infrastructure data, solving basic parameters, solving skylight occupation time, loading skylight data and solving heavy haul railway passing capacity. The invention brings the constraint factors of railway line infrastructure, skylight type, skylight time and station to departure number into the category of throughput capability calculation; the method considers the unequal running time of the train section, the number of stations to departure lines and the time constraint of train tracking intervals, causes the difference of the loss time outside the skylight of the V-shaped skylight and the vertical skylight and the influence on the calculation of the passing capacity, provides a method for calculating the passing capacities of the V-shaped skylight and the vertical skylight of the heavy-load railway, and has high calculation accuracy.

Description

Heavy haul railway throughput computing method
Technical Field
The invention relates to the field of heavy haul railways, in particular to a method for calculating the passing capacity of a heavy haul railway.
Background
At present, the railway passing capability calculation method mainly comprises a deduction coefficient method, a minimum average interval method, a train operation diagram paving method through a computer and the like; the analysis of the factors influencing the passing ability mainly comprises an overrun scheme between passenger and cargo trains, a train stop scheme, a passenger and cargo train running line combination scheme, skylight arrangement and the like.
But currently there are the following problems:
(1) The heavy haul railway only runs heavy haul freight trains, the running speeds of the trains are basically the same, no over-running condition exists among the freight trains, the existing deduction coefficient method is difficult to solve the passing capacity of the heavy haul railway, and the passing capacity calculation method for the heavy haul railway needs to be researched according to the characteristics of the deduction coefficient method.
(2) Existing researches do not consider the unequal factors of the section length in the section and the running time of the train section into the calculation range of the passing capacity of the heavy haul railway, but in actual operation, the unequal factors cause larger loss time outside the skylight and have great influence on the passing capacity, so that the calculation result is inaccurate.
(3) The existing research does not consider the number of arrival/departure lines of the intermediate station for the train to stop into the calculation range of the passing capacity of the heavy haul railway, but in actual operation, when a vertical skylight is arranged, part of trains need to stop at the intermediate station, and when the distance between the stations is large, the influence of the number of arrival/departure lines of the intermediate station for the train to stop on the passing capacity is large, so that the calculation result is inaccurate.
(4) The existing researches do not consider the V-shaped skylight and the vertical skylight in the calculation range of the passing capacity of the heavy haul railway at the same time, but in actual operation, the V-shaped skylight and the vertical skylight are required to be opened for ensuring normal operation, so that the calculation result is inaccurate.
Therefore, it is necessary to provide a new method for calculating the throughput of the heavy haul railway, so as to improve the accuracy of calculating the throughput of the heavy haul railway.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for calculating the passing capability of a heavy-load railway, which comprises the steps of loading infrastructure data, solving basic parameters, solving the occupation time of a skylight, configuring skylight data, solving the passing capability of the heavy-load railway, and integrating factors such as infrastructure, different types of sky windows and the like into the category of the calculation of the passing capability, so that the calculation accuracy of the passing capability of the heavy-load railway is improved.
The invention is realized according to the following technical scheme.
A calculation method for the passing capacity of a heavy haul railway comprises the following steps:
s1, acquiring infrastructure data of a heavy haul railway, and classifying the infrastructure data according to a railway line interval and a power supply arm range;
s2, solving by adopting a train traction calculation method according to the loaded infrastructure data, wherein a solving result comprises interval running time, start-stop additional time and train tracking interval time;
s3, solving the skylight occupation time of each section and the power supply arm range by considering the skylight type, the number of stations reaching departure lines, train tracking intervals and the skylight occupation time model and constraint conditions of train section operation time factors;
s4, acquiring the opening types of the heavy-duty railway skylight and the opening time of different skylight types, and solving the proportion of the V-shaped skylight time and the vertical skylight time to the total skylight time;
s5, solving the passing capacity of the heavy haul railway according to the constraint conditions.
Further, the infrastructure data comprise railway line horizontal and vertical section data, station diagram and station to departure line data, tunnel gap data, power supply arm distribution data, maintenance skylight type data, train type and marshalling and rolling stock data;
the railway line flat profile data comprises line mileage, gradient, slope length, elevation, curve starting point mileage, curve length, curve radius and line speed limit;
the train type and consist and rolling stock data includes locomotive type and number, vehicle type and number, rolling stock performance data, train load.
Further, the basic parameter is solved by adopting a train traction calculation method according to the loading infrastructure data, and the solving result comprises interval running time
Figure GDA0004211376770000031
Additional time for starting and stopping>
Figure GDA0004211376770000032
Train tracking interval time (I);
the train tracking interval time is station departure interval time (I f ) Station inter-arrival time (I) d ) Station transit time interval (I) t ) And interval tracking interval time (I q ) The maximum of the four.
Further, the skylight occupation time refers to the time occupied by the loss of passing capacity caused by the fact that a train running line cannot be paved due to the fact that the skylight is opened;
the skylight occupation time (T w ) Comprising sunroof effective working time (T y ) Sunroof auxiliary operation time (T) f ) And sunroof outer loss time (T s );
The effective operation time of the skylight is the time for operation such as line maintenance, contact net maintenance and the like;
the auxiliary operation time of the skylight is the operation time required for handling confirmation and the like for maintaining the skylight;
the sunroof outer loss time includes a V-shaped sunroof outer loss time
Figure GDA0004211376770000033
And vertical skylight outer loss time->
Figure GDA0004211376770000034
The loss time outside the skylight is the ratio (B) of the loss time outside the V-shaped skylight and the time of the V-shaped skylight to the total skylight time V ) And the ratio of the vertical sunroof outer loss time to the vertical sunroof time to the total sunroof time (B H ) The sum of the products of (a) is given by:
Figure GDA0004211376770000035
furthermore, the effective operation time of the skylight is the actual operation time for providing equipment maintenance in the skylight, and is generally determined according to the actual skylight time of a railway operation department, wherein the actual operation time is related to the size of the adopted maintenance mechanical equipment and maintenance workload.
Furthermore, the auxiliary operation time of the skylight is the auxiliary operation time before and after the operation of the skylight is performed, and is the time required by the railway transportation department, the departments such as power supply and work and the inside of each department to perform a series of operations such as contact, confirmation, release, recitation commands and the like, and is usually determined according to the actual operation time of the railway operation department.
Further, the loss time outside the V-shaped skylight is a large value of the distance between stations in the range of a line section of the V-shaped skylight and the running time of a train in the range of the length of a power supply arm.
Further, the outside of the vertical skylight loses time
Figure GDA0004211376770000041
Train operation time in the line section for vertical sunroof>
Figure GDA0004211376770000042
Additional time for starting and stopping train>
Figure GDA0004211376770000043
And deducting the sum of all the train occupation time of each station of the tracking operation stop along the road which can be paved in the time range, wherein the formula is as follows:
Figure GDA0004211376770000044
wherein n is the sum of the number of departure lines available for the train to stop at all stations in the section.
Further, the calculation formula of the passing capacity of the heavy haul railway is as follows:
Figure GDA0004211376770000045
wherein I is I f 、I d 、I t And I q The maximum of the four.
Further, the constraints include railway line infrastructure constraints, sunroof type and sunroof time constraints, station to departure number constraints.
The invention has the following advantages and beneficial effects.
In the calculation and analysis process, the invention considers the restrictions of railway line infrastructure, skylight type, skylight time and station arrival and departure number, overcomes the limitation that the traditional throughput calculation cannot simultaneously consider skylight, station and line factors, improves the reliability of the calculation result of the throughput of the heavy-load railway, and has obvious popularization and application values.
Drawings
FIG. 1 is a flow chart of a heavy haul railway throughput calculation method of the present invention;
FIG. 2 is a schematic view of the occupancy time of the V-shaped skylight of the present invention;
fig. 3 is a schematic view of the vertical sunroof occupancy time of the present invention.
Wherein in FIG. 2, T w The time occupied by the skylight; t (T) y The effective operation time of the skylight; t (T) f 2 is 1/2 of the auxiliary operation time of the skylight; t (T) s Time is lost outside the skylight; a. b, c and d are respectively the names of the stations;
in FIG. 3, T y The effective operation time of the skylight; t (T) f 2 is 1/2 of the auxiliary operation time of the skylight;
Figure GDA0004211376770000051
run time for a-b interval; />
Figure GDA0004211376770000052
Additional time for a-b start-stop; i is train tracking interval time; a. b, c and d are respectively the names of the stations.
Detailed Description
In order to more clearly illustrate the present invention, the present invention will be further described with reference to the accompanying drawings and examples. It is to be understood by persons skilled in the art that the following detailed description is illustrative and not restrictive, and that this invention is not limited to the details given herein.
The invention discloses a calculation method of the passing capacity of a heavy haul railway, which is characterized in that under the constraint of unequal running time of a train section, the number of stations to departure lines and the tracking interval of trains, a skylight outer loss time model of a V-shaped skylight and a vertical skylight is established, and the calculation method of the passing capacity of the heavy haul railway is established according to the proportion of the opening time of the V-shaped skylight and the vertical skylight on the basis.
As shown in the flowchart of fig. 1, the method for calculating the passing capability of the heavy haul railway according to the present invention includes the steps of: loading infrastructure data (S1), solving basic parameters (S2), solving skylight occupation time (S3), loading skylight data (S4), and solving heavy haul railway passing capacity (S5).
S1, loading infrastructure data
The infrastructure data comprise railway line flat profile data, station diagram and station to departure line data, tunnel gap data, power supply arm distribution data, maintenance skylight type data, train type and marshalling and rolling stock data.
The railway line flat profile data comprises line mileage, gradient, slope length, elevation, curve starting point mileage, curve length, curve radius and line speed limit.
The train type and consist and rolling stock data includes locomotive type and number, vehicle type and number, rolling stock performance data, train load.
S2, solving basic parameters
The basic parameter solving is to adopt the section 1 of train traction calculation according to the loading infrastructure data: locomotive traction train (TB/T1407.1-2018) train traction calculation method is used for solving, and the solving result comprises interval running time
Figure GDA0004211376770000061
Additional time for starting and stopping>
Figure GDA0004211376770000062
Train tracking interval (I).
The train tracking interval time is station departure interval time (I f ) Station inter-arrival time (I) d ) Station transit time interval (I) t ) And interval tracking interval time (I q ) The maximum of the four.
S3, solving the skylight occupation time
The skylight occupation time refers to the time occupied by the loss of passing capacity caused by the fact that a train running line cannot be paved due to the fact that the skylight is opened.
The skylight occupation time (T w ) Comprising sunroof effective working time (T y ) Sunroof auxiliary operation time (T) f ) And sunroof outer loss time (T s )。
The effective operation time of the skylight is the actual operation time for line maintenance, contact net maintenance and the like, and the effective operation time of the skylight is valued according to the actual skylight time of a railway operation department.
The auxiliary operation time of the skylight is the operation time of the railway related departments for handling and maintaining the skylight, such as contact, confirmation, release, recitation and the like, and is valued according to the actual operation time of the railway operation department.
The sunroof outer loss time includes a V-shaped sunroof outer loss time
Figure GDA0004211376770000063
And vertical sunroof outer lost time
Figure GDA0004211376770000064
The loss time outside the skylight is the ratio (B) of the loss time outside the V-shaped skylight and the time of the V-shaped skylight to the total skylight time V ) And the ratio of the vertical sunroof outer loss time to the vertical sunroof time to the total sunroof time (B H ) And the sum of the products of (a) is provided.
The loss time outside the V-shaped skylight is a large value of the distance between stations in the range of the V-shaped skylight opening line section and the train running time in the range of the length of the power supply arm.
The outside loss time of the vertical skylight
Figure GDA0004211376770000071
Train operation time within the line section for vertical skylights
Figure GDA0004211376770000072
Additional time for starting and stopping train>
Figure GDA0004211376770000073
And deducting the sum of all the train occupation time of each station of the tracking operation stop along the road which can be paved in the time range, wherein the formula is as follows:
Figure GDA0004211376770000074
wherein n is the sum of the number of departure lines available for the train to stop at all stations in the section.
S4, loading skylight data
According to the time for opening the V-shaped skylight and the time for opening the vertical skylight of the heavy-duty railway, the total time for opening the skylight of the heavy-duty railway is obtained, so that the proportion of the time for opening the V-shaped skylight and the vertical skylight respectively accounting for the total skylight time, namely the proportion B of the time for loading the V-shaped skylight, is calculated V And vertical sunroof time ratio B H
The outside loss time of the skylight is obtained:
Figure GDA0004211376770000075
s5, solving the passing capacity of heavy haul railway
Solving the passing capacity of the heavy haul railway according to the obtained parameters:
Figure GDA0004211376770000076
wherein I is I f 、I d 、I t And I q The maximum of 4.
The foregoing examples of the present invention are merely illustrative of the present invention and are not intended to limit the embodiments of the present invention, and other variations or modifications of various forms may be made by those skilled in the art based on the foregoing description, and it is not intended to be exhaustive of all embodiments, and all obvious variations or modifications that come within the scope of the invention are defined by the following claims.

Claims (7)

1. The method for calculating the passing capacity of the heavy haul railway is characterized by comprising the following steps of:
s1, loading infrastructure data;
s2, solving basic parameters;
the basic parameter solving is to solve by adopting a train traction calculation method according to the loading infrastructure data, and the solving result comprises interval running time
Figure FDA0004211376760000011
Additional time for starting and stopping>
Figure FDA0004211376760000012
Train tracking interval time I;
the train tracking interval time is station departure interval time I f Station inter-arrival time I d Station passing interval time I t Interval tracking interval time I q The maximum of the four;
s3, solving the skylight occupation time;
the skylight occupation time refers to the time occupied by the loss of passing capacity caused by the fact that a train running line cannot be paved due to the fact that the skylight is opened;
the skylight occupation time T w Including sunroof effective working time T y Skylight auxiliary operation time T f And sunroof outer loss time T s
The effective operation time of the skylight is the operation time for line maintenance and contact net maintenance;
the auxiliary operation time of the skylight is the confirmation operation time required by handling and maintaining the skylight;
the sunroof outer loss time includes a V-shaped sunroof outer loss time
Figure FDA0004211376760000013
And vertical skylight outer loss time->
Figure FDA0004211376760000014
The loss time outside the skylight is the ratio B of the loss time outside the V-shaped skylight and the time of the V-shaped skylight to the total skylight time V And the ratio B of the vertical skylight time to the total skylight time H The sum of the products of (a) is given by:
Figure FDA0004211376760000015
s4, loading skylight data;
s5, solving the passing capacity of the heavy haul railway according to the constraint conditions;
the calculation formula of the passing capacity of the heavy haul railway is as follows:
Figure FDA0004211376760000021
wherein I is I f 、I d 、I t And I q The maximum of the four.
2. The method for calculating the passing capability of a heavy haul railway according to claim 1, wherein the infrastructure data comprises railway line flat profile data, station yard and departure to station data, tunnel gap data, power supply arm distribution data, maintenance skylight type data, train type and consist and rolling stock data;
the railway line flat profile data comprises line mileage, gradient, slope length, elevation, curve starting point mileage, curve length, curve radius and line speed limit;
the train type and consist and rolling stock data includes locomotive type and number, vehicle type and number, rolling stock performance data, train load.
3. The method for calculating the passing capacity of the heavy haul railway according to claim 1, wherein the effective working time of the skylight is the actual working time for providing equipment maintenance in the skylight, and is determined according to the actual skylight time of a railway operation department according to the adopted maintenance mechanical equipment and the maintenance workload.
4. The method for calculating the passing capacity of the heavy haul railway according to claim 1, wherein the auxiliary operation time of the skylight is the auxiliary operation time before and after the operation of the skylight is performed, and is the time required for the railway transportation department to perform a series of connection, confirmation, release and recitation commands with the power supply and the service department, and is determined according to the actual operation time of the railway operation department.
5. The method for calculating the passing capacity of a heavy haul railway according to claim 1, wherein the V-shaped skylight outer loss time is a large value of a station distance in a range of a V-shaped skylight opening line section and a train running time in a range of a power supply arm length.
6. The method for calculating the passing capacity of a heavy haul railway according to claim 1, wherein the vertical type sunroof has a lost time
Figure FDA0004211376760000031
Train operation time in the line section for vertical sunroof>
Figure FDA0004211376760000032
Additional time with train start-stop
Figure FDA0004211376760000033
And deducting the sum of all train occupation time of each station along the tracked operation stop paved in the time range, wherein the formula is as follows:
Figure FDA0004211376760000034
wherein n is the sum of the number of departure lines available for the train to stop at all stations in the section.
7. A method of calculating the throughput of heavy haul railway as defined in claim 1, wherein said constraints include railway line infrastructure constraints, sunroof type and sunroof time constraints, station-to-departure number constraints.
CN201911146733.5A 2019-11-21 2019-11-21 Heavy haul railway throughput computing method Active CN111090929B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911146733.5A CN111090929B (en) 2019-11-21 2019-11-21 Heavy haul railway throughput computing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911146733.5A CN111090929B (en) 2019-11-21 2019-11-21 Heavy haul railway throughput computing method

Publications (2)

Publication Number Publication Date
CN111090929A CN111090929A (en) 2020-05-01
CN111090929B true CN111090929B (en) 2023-06-23

Family

ID=70393477

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911146733.5A Active CN111090929B (en) 2019-11-21 2019-11-21 Heavy haul railway throughput computing method

Country Status (1)

Country Link
CN (1) CN111090929B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111874045B (en) * 2020-07-10 2022-03-08 交控科技股份有限公司 Method and device for determining main line passing capacity of urban rail transit
CN111680849B (en) * 2020-08-11 2020-12-29 北京交通大学 Method for calculating station passing capacity under abnormal event, storage medium and terminal
CN113415325B (en) * 2021-08-09 2022-01-28 中国铁道科学研究院集团有限公司通信信号研究所 Rail transit station arrival interval calculation system and method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107226094A (en) * 2017-06-12 2017-10-03 中铁第六勘察设计院集团有限公司 A kind of heavy haul railway parking lot " machine of 2 wire clamp 1 is walked " wiring harness structure design
CN108491950A (en) * 2018-01-25 2018-09-04 北京交通大学 A kind of high-speed railway handling capacity computational methods considering multiple resources constraint

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107226094A (en) * 2017-06-12 2017-10-03 中铁第六勘察设计院集团有限公司 A kind of heavy haul railway parking lot " machine of 2 wire clamp 1 is walked " wiring harness structure design
CN108491950A (en) * 2018-01-25 2018-09-04 北京交通大学 A kind of high-speed railway handling capacity computational methods considering multiple resources constraint

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
崔石军 ; .天窗开行方式对线路通过能力的影响分析.硅谷.2011,(05),149. *
牛会想 ; .重载铁路通过能力计算方法研究.铁道运输与经济.2009,(09),76-78. *

Also Published As

Publication number Publication date
CN111090929A (en) 2020-05-01

Similar Documents

Publication Publication Date Title
CN111090929B (en) Heavy haul railway throughput computing method
CN109657845B (en) Urban rail transit train timetable optimization system for time-varying passenger flow
CN107704950A (en) A kind of city rail train figure optimization method based on trip requirements and energy saving of system
US10902688B2 (en) System and method for determining grade errors of a route
CN107871179B (en) Railway freight train operation diagram compiling method based on arrival time limit
CN107016633B (en) Method for generating planned departure time of receiving, transporting and transporting buses of urban rail transit
Zubkov et al. Capacity and traffic management on a heavy-traffic railway line
Numonjonovich et al. Investigation of the Throughput of Railway and Road Transport of Angren-pop Through the Kamchik Pass and Ways of Its Further Improvement
CN114047765B (en) Virtual marshalling-oriented multi-train self-adaptive cruise control method and system
WO2015145900A1 (en) Service management device, train control method, and program
CN107555109A (en) A kind of tree-shaped track for tubular type logistics system
CN104527649A (en) Railway shunting-free marshalling station system
CN102060201B (en) Arrangement method of coal-unloading device of wagon tipper with wagon pulled in or out from transfer platform side
CN201459610U (en) D-type construction temporary beam of improved railway
Mlinarić et al. Energy efficiency of railway lines
CN103226666B (en) Sophisticated functions intercommunication road crossing system efficacy assessments method
Rickett Intermodal train loading methods and their effect on intermodal terminal operations
Butunov et al. Determination of factors affecting wagon turnover time
CN107392478A (en) A kind of preparation method of railway fast freight traffic production scheme
CN115796518A (en) Comprehensive evaluation method for traffic organization and transportation scheduling for intelligent construction site
CN109466572A (en) A kind of train is without sensor load calculation method
CN105292190B (en) Energy-saving train operation figure optimizes system and method
Belošević et al. Challenges in the railway yards layout designing regarding the implementation of intermodal technologies
CN110406568B (en) Method and system for balanced control of multiple trains in interval
CN109886467B (en) Urban ULS and road integrated cargo transportation network flow distribution system and method

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