CN114089046B - ZPW-2000 series track circuit compensation capacitance parameter estimation method and device - Google Patents
ZPW-2000 series track circuit compensation capacitance parameter estimation method and device Download PDFInfo
- Publication number
- CN114089046B CN114089046B CN202111397610.6A CN202111397610A CN114089046B CN 114089046 B CN114089046 B CN 114089046B CN 202111397610 A CN202111397610 A CN 202111397610A CN 114089046 B CN114089046 B CN 114089046B
- Authority
- CN
- China
- Prior art keywords
- induced voltage
- data
- zpw
- track circuit
- voltage amplitude
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000003137 locomotive effect Effects 0.000 claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 241000282461 Canis lupus Species 0.000 claims abstract description 11
- 239000003990 capacitor Substances 0.000 claims description 34
- 230000006870 function Effects 0.000 claims description 24
- 238000001514 detection method Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 238000010292 electrical insulation Methods 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/26—Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
- G01R27/2605—Measuring capacitance
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Train Traffic Observation, Control, And Security (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
The invention discloses a ZPW-2000 series track circuit compensation capacitance parameter estimation method and device, and relates to the technical field of railway signal equipment fault processing. Comprising the following steps: when the actual measurement data of the induced voltage of the locomotive is utilized to estimate the compensation capacitance parameter, firstly, sample data are preprocessed, the capacitance parameter of the compensation capacitance is taken as a decision variable, the difference between the minimized actual data of the envelope of the induced voltage amplitude and the theoretical data of the envelope of the induced voltage amplitude calculated based on a mathematical model is taken as a target, an adaptability function is constructed, and according to the characteristics that a gray-wolf optimization algorithm is insensitive to an initial solution value, has high optimization efficiency, good global optimizing performance and the like, the optimal compensation capacitance parameter value is searched for by adopting the gray-wolf optimization algorithm in an iterative manner, and the estimation of the compensation capacitance parameter of the ZPW-2000 series track circuit can be realized, and the adaptability and the accuracy are high.
Description
Technical Field
The invention relates to the technical field of railway signal equipment fault processing, in particular to a ZPW-2000 series track circuit compensation capacitance parameter estimation method and device.
Background
The track circuit uses the steel rail as a conductor, is used for automatically and continuously detecting whether the circuit is occupied by the vehicle, and simultaneously transmits control information to the train. The ZPW-2000 series non-insulated track circuit commonly adopted in high-speed railways in China leads the track transmission characteristic to tend to be resistive by adding a compensation capacitor on the track. Factors such as surge voltage, lightning and track current skin effect of the traction section can cause faults of a compensation capacitor, the transmission characteristic of a track circuit can be deteriorated due to the faults of the compensation capacitor, the effective transmission distance of a track signal is shortened, and the faults of a red light band are caused. Therefore, the method for diagnosing the faults has important practical significance for maintenance personnel.
At present, a railway signal dynamic detection system is widely applied, and the system is arranged on a comprehensive detection vehicle and generates a corresponding locomotive induction voltage signal in an electromagnetic induction mode by utilizing a track signal receiving antenna. The comprehensive detection vehicle is adopted by the railway electric service department to detect the ZPW-2000 track circuit for 2 times each month, and the detection data provides data support for a fault detection algorithm.
Disclosure of Invention
The invention aims to overcome the defect that misdiagnosis is easy to occur when railway faults are diagnosed in the prior art, and provides a ZPW-2000 series track circuit compensation capacitance parameter estimation method and device.
In order to achieve the above object, the present invention provides the following technical solutions:
a ZPW-2000 series track circuit compensation capacitance parameter estimation method comprises the following steps:
s1: acquiring induction voltage data of a locomotive signal sample acquisition point in a dynamic detection system in real time;
s2: preprocessing the acquired induced voltage data;
s3: envelope extraction is carried out on the preprocessed induced voltage data, and actual data A Icv (x) of an induced voltage amplitude envelope (Locomotive Signal Amplitude Envelope, LSAE) is obtained;
s4: for a track section, calculating induced voltage amplitude envelope theoretical data corresponding to the track section in a mathematical modeling mode, and establishing a basic database according to the induced voltage amplitude envelope theoretical data;
The induced voltage amplitude envelope theory data is expressed as:
ARcv(x)＝F(X，C)
Wherein, C= [ C 1,c2,...cN ] is the compensation capacitance, N is the number of the compensation capacitance; x= [ X 1,x2,...,xK ] is the position corresponding to the compensation capacitor, K is the number of locomotive signal sample acquisition points from the transmitting end to the train shunt point X K, and a Rcv (X) is the induced voltage amplitude and comprises theoretical data;
s5: constructing a difference function according to the induced voltage amplitude envelope theoretical data and the induced voltage amplitude envelope actual data, and taking the minimization of the difference function as a target construction objective function;
The difference function is:
The objective function is:
Wherein A Icv(xi) is the actual data of the induced voltage amplitude envelope of the position corresponding to the ith compensation capacitor, C= [ C 1,c2,...cN ] is the compensation capacitor, and N is the number of the compensation capacitors; x= [ X 1,x2,...,xK ] is the position corresponding to the compensation capacitor, and K is the number of locomotive signal sample acquisition points from the transmitting end to the train shunt point X K;
s6: according to the standard range of the ZPW-2000 series track circuit compensation capacitor, taking the compensation capacitor as a decision variable to construct constraint conditions;
S7: and (3) solving the objective function in the step (S5) by using a gray wolf optimization algorithm (Grey Wolf Optimizer, GWO) based on the objective function and the constraint condition, and estimating the parameter value of the ZPW-2000 series track circuit compensation capacitor to obtain an optimal compensation capacitor calculation value.
According to the technical scheme, when the actual measurement data of the induced voltage of the locomotive is utilized to estimate the compensation capacitance parameter, firstly, sample data are preprocessed, the compensation capacitance parameter is taken as a decision variable, the difference between the minimized actual data of the envelope of the induced voltage amplitude and the theoretical data of the envelope of the induced voltage amplitude calculated based on a mathematical model is taken as a target, an adaptability function is constructed, and according to the characteristics of insensitivity to an initial solution value, higher optimization efficiency, good global optimizing performance and the like of a gray-wolf optimization algorithm (Grey Wolf Optimizer, GWO), the optimal compensation capacitance parameter value is iteratively searched by adopting the gray-wolf optimization algorithm (Grey Wolf Optimizer, GWO), so that the estimation of the compensation capacitance parameter of the ZPW-2000 series track circuit can be realized, and the adaptability and the accuracy are higher.
As a preferred embodiment of the present invention, in the step S2, the preprocessing includes filtering out harmonics, reducing noise, and normalizing.
As a preferred embodiment of the present invention, in the step S4, the mathematical modeling is performed by using a uniform transmission line theory.
In a preferred embodiment of the present invention, in the step S6, the standard range is between [ C t-5％Ct,Ct+ 5％Ct ], where C t is a standard value of the compensation capacitor, and generally 40 μf is taken.
As a preferred embodiment of the present invention, in the step S6, the constraint condition is expressed as:
0≤cj≤Ct+5％Ct
Wherein c j is the jth compensation capacitor.
On the other hand, the fault diagnosis device of the electric insulation joint of the track circuit is arranged on the comprehensive detection vehicle, and the real-time interface of the railway signal dynamic detection system comprises at least one processor and a memory which is in communication connection with the at least one processor; the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of any one of the above.
Compared with the prior art, the invention has the beneficial effects that: when the actual measurement data of the induced voltage of the locomotive is utilized to estimate the compensation capacitance parameter, firstly, sample data are preprocessed, the capacitance parameter of the compensation capacitance is taken as a decision variable, the minimized difference between the actual data of the envelope of the induced voltage amplitude and the theoretical data of the envelope of the induced voltage amplitude calculated based on a mathematical model is taken as a target, an adaptability function is constructed, and according to the characteristics of insensitivity to an initial solution value, higher optimization efficiency, good global optimizing performance and the like of a gray-wolf optimization algorithm (Grey Wolf Optimizer, GWO), the optimal compensation capacitance parameter value is searched for by adopting the gray-wolf optimization algorithm (Grey Wolf Optimizer, GWO) in an iterative manner, so that the estimation of the compensation capacitance parameter of a ZPW-2000 series track circuit can be realized, and the adaptability and the accuracy are higher.
Drawings
FIG. 1 is a flow chart of a method for estimating the compensation capacitance parameters of ZPW-2000 series track circuits according to embodiment 1 of the present invention;
fig. 2 is a block diagram of a ZPW-2000 series track circuit compensation capacitance parameter estimation device according to embodiment 3 of the present invention.
Detailed Description
The present invention will be described in further detail with reference to test examples and specific embodiments. It should not be construed that the scope of the above subject matter of the present invention is limited to the following embodiments, and all techniques realized based on the present invention are within the scope of the present invention.
Example 1
A ZPW-2000 series track circuit compensation capacitance parameter estimation method, as shown in figure 1, comprises the following steps:
S1: data acquisition, namely acquiring induced voltage data of a locomotive signal sample acquisition point in a dynamic detection system in real time;
S2: preprocessing data, namely preprocessing the acquired induced voltage data, wherein the preprocessing comprises harmonic wave filtering, noise reduction and normalization;
S3: data conversion, namely carrying out envelope extraction on the preprocessed induced voltage data to obtain actual data A Icv (x) of an induced voltage amplitude envelope (Locomotive Signal Amplitude Envelope, LSAE);
s4: establishing a basic database, calculating induced voltage amplitude envelope theoretical data A Rcv (x) corresponding to the region by using a mathematical modeling mode for the track section, and establishing the basic database by using the induced voltage amplitude envelope theoretical data A Rcv (x), wherein the induced voltage amplitude envelope theoretical data A Rcv (x) is expressed as:
ARcv(x)＝F(X，C)
wherein, C= [ C 1,c2,...cN ] is the compensation capacitance, N is the number of the compensation capacitance; x= [ X 1,x2,...,xK ] is the position corresponding to the compensation capacitor, K is the number of locomotive signal sample acquisition points from the transmitting end to the train shunt point X K, and a Rcv (X) is the induced voltage amplitude and comprises theoretical data;
s5: constructing a difference function according to the induced voltage amplitude envelope theoretical data and the induced voltage amplitude envelope actual data, and taking the minimization of the difference function as a target construction objective function, wherein the difference function is as follows:
The objective function is:
wherein, A Icv(xi) is the actual data of the induced voltage amplitude envelope of the position corresponding to the ith compensation capacitor;
S6: constructing constraint conditions, namely constructing constraint conditions by taking the compensation capacitor as a decision variable according to the standard range of the ZPW-2000 series track circuit compensation capacitor;
S7: and (3) estimating parameters, namely solving the objective function in the step (S5) by using a gray-wolf optimization algorithm (Grey WolfOptimizer, GWO) based on the objective function and the constraint condition, and estimating the parameter values of the ZPW-2000 series track circuit compensation capacitors to obtain optimal compensation capacitor calculation values.
In the step S2, the preprocessing includes filtering out harmonics, reducing noise, and normalizing.
In the step S4, the mathematical modeling method adopts a uniform transmission line theory to perform modeling.
In the step S6, the standard range is between [ C t-5％Ct,Ct+5％Ct ], where C t is a standard value of the compensation capacitor, and generally 40 μf is taken.
In the step S6, the constraint condition is expressed as:
0≤cj≤Ct+5％Ct
Wherein c j is the jth compensation capacitor.
Example 2
The ZPW-2000 series track circuit compensation capacitance parameter estimation method is different from embodiment 1 in that:
The method provided by the invention is verified by selecting a group of data acquired from a signal dynamic detection system on the comprehensive detection vehicle and using the detection result of the ZPW-2000 series track circuit as a comparison object, and the field detection result shows that: the capacitance of the 7 th compensation capacitor is reduced to 10 mu F;
The method provided by the invention is used for estimating the compensation capacitance parameter of the track circuit, the obtained result is shown in table 1, the estimated result of C 7 is 10.2867 mu F, the estimated result is reduced to be less than one fourth of the standard value, and field maintenance personnel can replace the capacitance in advance according to the result.
Table 1 verification results and errors
Example 3
The embodiment also provides a fault diagnosis device for the track circuit electric insulation joint, which is arranged on the comprehensive detection vehicle, and a real-time interface of the railway signal dynamic detection system, as shown in figure 2, and comprises at least one processor and a memory in communication connection with the at least one processor; the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of any one of the above.
According to the technical scheme, when the actual measurement data of the induced voltage of the locomotive is utilized to estimate the compensation capacitance parameter, firstly, sample data are preprocessed, the compensation capacitance parameter is taken as a decision variable, the difference between the minimized actual data of the envelope of the induced voltage amplitude and the theoretical data of the envelope of the induced voltage amplitude calculated based on a mathematical model is taken as a target, an adaptability function is constructed, and according to the characteristics that the gray-wolf optimization algorithm (GreyWolfOptimizer, GWO) is insensitive to an initial solution value, the optimization efficiency is high, the global optimizing performance is good and the like, the gray-wolf optimization algorithm (Grey WolfOptimizer, GWO) is adopted to iteratively search the optimal compensation capacitance parameter value, so that the estimation of the compensation capacitance parameter of the ZPW-2000 series track circuit can be realized, and the adaptability and the accuracy are higher.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (6)
1. A ZPW-2000 series track circuit compensation capacitance parameter estimation method is characterized by comprising the following steps:
s1: acquiring induction voltage data of a locomotive signal sample acquisition point in a dynamic detection system in real time;
s2: preprocessing the acquired induced voltage data;
s3: envelope extraction is carried out on the preprocessed induced voltage data, and actual induced voltage amplitude envelope data are obtained;
s4: for a track section, calculating induced voltage amplitude envelope theoretical data corresponding to the track section in a mathematical modeling mode, and establishing a basic database according to the induced voltage amplitude envelope theoretical data;
The induced voltage amplitude envelope theory data is expressed as:
ARcv(x)＝F(X，C)
wherein, C= [ C 1,c2,...cN ] is the compensation capacitance, N is the number of the compensation capacitance; x= [ X 1,x2,...,xK ] is the position corresponding to the compensation capacitor, K is the number of locomotive signal sample acquisition points from the transmitting end to the train shunt point X K, and a Rcv (X) is the induced voltage amplitude and comprises theoretical data;
s5: constructing a difference function according to the induced voltage amplitude envelope theoretical data and the induced voltage amplitude envelope actual data, and taking the minimization of the difference function as a target construction objective function;
The difference function is:
The objective function is:
Wherein A Icv (xi) is the actual data of the induced voltage amplitude envelope of the position corresponding to the ith compensation capacitor, C= [ C 1,c2,...cN ] is the compensation capacitor, and N is the number of the compensation capacitors; x= [ X 1,x2,...,xK ] is the position corresponding to the compensation capacitor, and K is the number of locomotive signal sample acquisition points from the transmitting end to the train shunt point X K;
s6: according to the standard range of the ZPW-2000 series track circuit compensation capacitor, taking the compensation capacitor as a decision variable to construct constraint conditions;
s7: and based on the objective function and the constraint condition, solving the objective function in the step S5 by using a gray wolf optimization algorithm, and estimating the parameter value of the ZPW-2000 series track circuit compensation capacitor to obtain an optimal compensation capacitor calculated value.
2. A ZPW-2000 series track circuit compensating capacitance parameter estimation method as in claim 1, wherein the preprocessing includes filtering out harmonics, reducing noise, and normalizing in step S2.
3. The ZPW-2000 series track circuit compensating capacitance parameter estimation method as recited in claim 1, wherein in step S4, the mathematical modeling is performed using a uniform transmission line theory.
4. The ZPW-2000 series track circuit compensating capacitance parameter estimation method as claimed in claim 1, wherein in the step S6, the standard range is between [ C t-5％Ct,Ct+5％Ct ], wherein C t is a standard value of the compensating capacitance.
5. The ZPW-2000 series track circuit compensating capacitance parameter estimation method as recited in claim 1, wherein in the step S6, the constraint condition is expressed as:
0≤cj≤Ct+5％Ct
Wherein C j is the j-th compensation capacitor, and C t is the standard value of the compensation capacitor.
6. A rail circuit electrical insulation joint fault diagnosis device, comprising at least one processor and a memory in communication with the at least one processor; the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111397610.6A CN114089046B (en) | 2021-11-23 | 2021-11-23 | ZPW-2000 series track circuit compensation capacitance parameter estimation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111397610.6A CN114089046B (en) | 2021-11-23 | 2021-11-23 | ZPW-2000 series track circuit compensation capacitance parameter estimation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114089046A CN114089046A (en) | 2022-02-25 |
| CN114089046B true CN114089046B (en) | 2024-04-26 |
Family
ID=80303535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111397610.6A Active CN114089046B (en) | 2021-11-23 | 2021-11-23 | ZPW-2000 series track circuit compensation capacitance parameter estimation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114089046B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102556118A (en) * | 2012-01-06 | 2012-07-11 | 北京交通大学 | Fault online diagnosis method of uninsulated track circuit tuning zone equipment |
| WO2014190759A1 (en) * | 2013-05-27 | 2014-12-04 | 国家电网公司 | Svc compensation strategy optimization method |
| CN108872715A (en) * | 2018-08-28 | 2018-11-23 | 北京交大思诺科技股份有限公司 | Track circuit compensation capacitance real-time detecting system |
| CN109768573A (en) * | 2019-01-29 | 2019-05-17 | 三峡大学 | Reactive power optimization method for distribution network based on multi-objective differential gray wolf algorithm |
-
2021
- 2021-11-23 CN CN202111397610.6A patent/CN114089046B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102556118A (en) * | 2012-01-06 | 2012-07-11 | 北京交通大学 | Fault online diagnosis method of uninsulated track circuit tuning zone equipment |
| WO2014190759A1 (en) * | 2013-05-27 | 2014-12-04 | 国家电网公司 | Svc compensation strategy optimization method |
| CN108872715A (en) * | 2018-08-28 | 2018-11-23 | 北京交大思诺科技股份有限公司 | Track circuit compensation capacitance real-time detecting system |
| CN109768573A (en) * | 2019-01-29 | 2019-05-17 | 三峡大学 | Reactive power optimization method for distribution network based on multi-objective differential gray wolf algorithm |
Non-Patent Citations (1)
| Title |
|---|
| 基于遗传算法的无绝缘轨道电路故障综合诊断方法;赵林海;冉义奎;穆建成;;中国铁道科学;20100515(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114089046A (en) | 2022-02-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105044524B (en) | A kind of arrester monitoring method and system suitable for intelligent substation | |
| AU2011313840B2 (en) | Apparatus and method for fault detection and location determination | |
| CN105738764B (en) | Fault Section Location of Distribution Network based on transient information Whole frequency band | |
| CN101907437A (en) | A Cable Fault Location Method Based on Wavelet Difference Algorithm | |
| CN111766470B (en) | Fault location method and system for high-voltage direct current transmission line and direct current transmission line | |
| CN108562820A (en) | Single-phase ground fault line selecting method of small-electric current grounding system and device | |
| CN114636896A (en) | Single-phase grounding high-resistance fault traveling wave positioning method for power distribution network by utilizing kurtosis | |
| CN113945806B (en) | Single-ended fault location method and system for flexible direct current transmission line | |
| CN113671314A (en) | Method for positioning and ranging single-phase earth fault section of ring network of power distribution network | |
| CN118731548A (en) | New energy vehicle charger automatic test system | |
| CN108646136A (en) | A kind of DC line method of single end distance measurement and system | |
| CN111157848A (en) | Traveling wave fault positioning method and system for high-speed railway traction power supply system | |
| CN114089046B (en) | ZPW-2000 series track circuit compensation capacitance parameter estimation method and device | |
| CN107797018B (en) | Grounding detection circuit, device and method for diesel locomotive | |
| CN119310493B (en) | Centralized station terminal with data real-time acquisition and analysis function | |
| CN116125196A (en) | High-voltage cable fault traveling wave ranging system and method | |
| CN114301175A (en) | Power distribution station area user transformation relation identification method and device based on injection signals | |
| CN119093599A (en) | An intelligent analysis method for the operation status of distribution network transformers | |
| CN117825876A (en) | Cable fault feature extraction method and device combining feature fragments and wavelet packets | |
| CN114779028B (en) | Generator partial discharge online monitoring device and monitoring method | |
| CN114200246B (en) | Fault location method and system for direct-current transmission line | |
| CN106680632B (en) | Iron tower lightning protection performance detection method and system based on electromagnetic radiation abnormity judgment | |
| CN111257629B (en) | New energy station power characteristic detection method, device and system | |
| CN115511111A (en) | Power station operation and maintenance state monitoring method based on time sequence data multi-scale features | |
| CN114118164A (en) | Fault diagnosis method and device for electric insulation joint of track circuit |
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 |