CN104793106A - Distribution network line fault section positioning method based on current break rate - Google Patents

Distribution network line fault section positioning method based on current break rate Download PDF

Info

Publication number
CN104793106A
CN104793106A CN201510210539.4A CN201510210539A CN104793106A CN 104793106 A CN104793106 A CN 104793106A CN 201510210539 A CN201510210539 A CN 201510210539A CN 104793106 A CN104793106 A CN 104793106A
Authority
CN
China
Prior art keywords
fault
current
phase
centerdot
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510210539.4A
Other languages
Chinese (zh)
Other versions
CN104793106B (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201510210539.4A priority Critical patent/CN104793106B/en
Publication of CN104793106A publication Critical patent/CN104793106A/en
Application granted granted Critical
Publication of CN104793106B publication Critical patent/CN104793106B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Locating Faults (AREA)

Abstract

The invention discloses a distribution network line fault section positioning method based on current break rate. The distribution network line fault section positioning method includes analyzing phase current characteristics of fault phase current before faults and before an arc suppression coil acts for ground faults, extracting fault characteristics from the phase current characteristics, and positioning by adopting full-process waveform current break rate. The distribution network line fault section positioning method only needs to measure phase current of line faults, is simple and high in applicability, and the problems of fault current weakness, poor reliability and sensitivity when a small-current grounding system has single-phase grounding faults commonly can be well solved. Meanwhile, interference cannot be introduced.

Description

Based on the distribution line fault section location method of current break rate
Technical field
The present invention relates to distribution network line fault diagnostic method, specifically a kind of distribution line fault section location method based on current break rate.
Background technology
According to statistics, electric system is in operational process, and the power outage caused by Distribution Network Failure accounts for more than 95% of total power outage, and wherein the accident of 70% is caused by singlephase earth fault or bus-bar fault.And distribution neutral point extensively adopts non-effective earthing (small current neutral grounding) mode both at home and abroad, to avoid causing power failure when there is singlephase earth fault.For the singlephase earth fault of distribution, because fault characteristic value is faint, lack reliable fault line selection and fault locating method always.Along with people are to the raising of distribution automation level requirement, more in the urgent need to fundamentally solving the fault-location problem of distribution.
The Fault Locating Method that current Chinese scholars proposes is roughly divided into two classes: one is injecting signal, and two is locate based on the section of fault characteristic value.Injecting signal comprises " S " injection method, the comprehensive injection method of alternating current-direct current and parallel medium resistance method, these class methods increase the interference to system, and can not detect instantaneity and intermittent grounding fault.Section location based on fault characteristic value comprises zero mould electric current relative method, section zero sequence admittance method, zero sequence reactive power direction method, the location based on difference of phase currents, residual method of addition, traveling wave method etc., electrical power distribution automatization system mainly utilizes main website to realize the time synchronized of FTU, is at least several millisecond to time error.In the case, the method such as to compare no longer valid for the amplitude of transient signal, polarity, Waveform Correlation.
Observe known, the data after current method all only utilizes fault to occur, and have ignored the utilization to information before fault.Meanwhile, most of localization method only considers zero-sequence current information (needing three phase information), high to transformer request, and acquisition of information is complicated, and these methods will lose efficacy when data phase shortage.
Summary of the invention
The object of the invention is to overcome above-mentioned the deficiencies in the prior art, propose a kind of distribution line fault section location method based on current break rate.
Principle of the present invention:
If fault phase is A phase, before fault, check point phase voltage is after fault, phase voltage is by Fig. 1, easily know that the two meets
U · ′ A = U · A - U · 0
1) adjacent check point 1,2 is positioned at trouble spot heteropleural (without branched line between check point), if check point 1 is positioned at upstream, trouble spot, check point 2 is positioned at downstream, trouble spot.
Before fault, the phase current of check point 1 is
I · 1 A = I · LA + I · C 1 A
In formula, for load current, for A phase is from check point 1 to the capacitance current of this line end, can be expressed as
I · C 1 A = jω C 1 A U · A
In formula, C 1Afor A phase is from check point 1 to the electric capacity of this line end.
Before and after assumed fault occurs, a cycle internal loading electric current is constant, and after fault, the faulted phase current of check point 1 is
I · ′ 1 A = I · LA + I · ′ C 1 A + I · f
In formula, be respectively A phase capacitance current and fault current.Wherein
I · ′ C 1 A = jω C 1 A U · ′ A = jω C 1 A ( U · A - U · 0 )
The inductive current that when arc suppression coil utilizes fault, neutral point voltage skew produces carrys out bucking-out system capacitance current, and earth point residual flow (fault current) can be expressed as
I · f = I · L - I · CΣ = U · 0 Z L - jω C Σ U · 0
In formula, C Σfor the whole network line mutual-ground capacitor.Fault occurs and in the arc suppression coil not operating period, preset arc suppression coil resistance is very large, and post-set arc suppression coils reactance is very large, generally speaking | and Z l| very large, in formula, Section 2 can be ignored.Then have
I · f ≈ - jω c Σ U · 0
In like manner, the phase current of check point 2 before fault
I · 2 A = I · LA + I · C 2 A = I · LA + jω C 2 A U · A
Check point 2 faulted phase current after fault
I · ′ 2 A = I · LA + I · ′ C 2 A = I · LA + jω C 2 A U · ′ A
According to above analysis, before fault, the difference of adjacent check point phase current
I · 1 A - I · 2 A = I · C 1 A - I · C 2 A = jω C 12 A U · A - - - ( 1 )
In formula, C 12Afor the relatively electric capacity between check point 1,2, its numerical value is very little.
After fault, the difference of adjacent check point faulted phase current
I · ′ 1 A - I · ′ 2 A = I · ′ C 1 A + I f - I · ′ C 2 A = jω C 12 A ( U · A - U · 0 ) - jω C Σ U · 0 = jω C 12 A U · A - jω ( C 12 A + C Σ ) U · 0 - - - ( 2 )
In formula, C 12A+ C Σ> > C 12A, size relevant with fault resistance size, be generally greater than 20% of phase voltage.Therefore Section 2 amplitude is much larger than Section 1, and amplitude and the direction of the difference of faulted phase current are determined by Section 2.
By comparing (1), (2) two formulas, after fault occurs, the difference of the faulted phase current of fault section two ends check point is before fault generation, and amplitude enlarges markedly.
2) adjacent check point 1,2 is positioned at trouble spot homonymy.
Similar to 1) in analysis, suppose that check point is all positioned at upstream, trouble spot, have:
Before fault
I · 1 A - I · 2 A = I · C 1 A - I · C 2 A = jω C 12 A U · A - - - ( 3 )
After fault
I · ′ 1 A = I · LA + I · ′ C 1 A + I · f
I · ′ 2 A = I · LA + I · ′ C 2 A + I · f
I · ′ 1 A - I · ′ 2 A = I · ′ C 1 A - I · ′ C 2 A = jω C 12 A U · ′ A - - - ( 4 )
After earth fault generation, faulted phase voltage reduces.Relatively (3), (4) two formulas are easily known, after fault occurs, the difference of the faulted phase current of non-faulting section two ends check point is before fault generation, and amplitude reduces.
If check point is all positioned at upstream, trouble spot, after having fault to occur equally, the spread value of the faulted phase current of non-faulting section two ends check point reduces.
Comprehensive 1), 2) analysis, can judge whether section exists fault according to the variation characteristic of the difference of section two ends faulted phase current.Wherein, the measurement of amplitude size provides by the effective value calculated in power frequency one cycle.
If i 1A(n), i 2An phase current sampling sequence that () is adjacent check point.Be designated as zero under making fault that the data point in moment occurs, define adjacent check point phase differential current:
Δi A(n)=i 1A(n)-i 2A(n)n∈[-N,N-1]
In formula, N=0.02f s, be a periodic sampling point.When n gets negative integer, before representing that fault occurs; N get positive integer represent fault occur after.Before fault occurs, difference between current effective value is
Δ I A , before = 1 N Σ n = - N - 1 [ Δ i A ( n ) ] 2
Δ I A , after = 1 N Σ n = 0 N - 1 [ Δ i A ( n ) ] 2
Definition section current break rate
A I=ΔI A,after/ΔI A,before
Analyze according to upper joint, easily know the A of non-faulting section ibe less than 1, the A of fault section i.be greater than 1.
Technical solution of the present invention is as follows:
A kind of distribution network line fault Section Location, its feature is, said method comprising the steps of:
Step S1, determines fault phase and fault moment t f:
When after generation earth fault being detected, according to phase voltage Changing Pattern, select fault phase, determine fault moment t according to phase voltage sudden change moment or arc-extinction device chugging moment f;
Step S2, chooses pick-up unit [t f-0.02s, t f+ 0.02s] interval interior 2N the faulted phase current Wave data put (before and after fault generation a cycle, i.e. positive and negative 0.02s, totally 2 power frequency periods, i.e. 2N data point) altogether, calculating current mutation rate A i, formula is as follows:
A I = 1 N Σ n = 0 N - 1 [ Δ i A ( n ) ] 2 / 1 N Σ n = - N - 1 [ Δ i A ( n ) ] 2
Δi A(n)=i 1A(n)-i 2A(n) n∈[-N,N-1]
In formula, i 1A(n), i 2An phase current sampling sequence that () is adjacent check point, it is that a periodic sampling is counted that fault is designated as zero, N under there is the data point in moment;
According to current break rate size, step S3, judges whether each section is fault section, namely works as A i> A iset, be judged as fault section, otherwise be then non-faulting section;
A isetfor the current break rate threshold value artificially arranged, be the positive number of about 1, consider certain nargin, desirable 1.2 ~ 1.5;
Step 4, according to check point distributing position on the line, to travel through successively, until find out fault section.
Compared with prior art, the invention has the beneficial effects as follows: to earth fault, by analysis of failure phase current before fault occurs and the phase current feature of fault generation to arc suppression coil action, therefrom extract fault characteristic value, and adopt the waveform-shaped current mutation rate of overall process to position.Therefore, only need the faulted phase current of measuring circuit, only consider the tradition of zero sequence (needing three phase information) before having broken, data acquisition is simple, applicability is strong; Synchronous from the angle of signal, total system adopts through the synchronous Monitoring Data of GPS, makes the difference between different check point sensitiveer.The problem that fault current faint, poor reliability when therefore can solve current ubiquitous single-phase grounded malfunction in grounded system of low current very well, sensitivity are low, can not introduce interference to system simultaneously.
Accompanying drawing explanation
Fig. 1 is single-phase earthing of small current earthing system schematic diagram
Fig. 2 is distributed fault section positioning system structure figure
Fig. 3 is 10kV system emulation figure
Fig. 4 is earth fault current change quantity waveform along the line
Embodiment
Fault waveform required for the present invention comes from distributed fault section positioning system, and system architecture as shown in Figure 2.Distribution network failure section positioning system is by master station, transformer station's (bus) measurement mechanism and be distributed in distribution line node failure locating device everywhere and form.Circuit is divided into some sections by localization of fault node in topology, each node installation three groups of measurement mechanisms, and real-time synchronization gathers circuit three-phase current and voltage.
According to the Fault Locating Method of invention, in 10kV power distribution network analogue system, dissimilar fault is set.System construction drawing as shown in Figure 3, be 1., 2., 3. sector number, 2. go up at section by fault verification.Sample frequency is 20kHz (each cycle data point N=400), and the moment occurs fault is 0.7s, and during earth fault, arc-extinction device is set to 0.04s actuation time, and the threshold value of current break rate is set to 1.2.
Embodiment citing for fault section judges:
Step S1, system is according to after residual voltage start detection to the generation of earth fault, and record A phase in bus three-phase voltage and reduce, B, C phase raises, and is defined as A phase fault; Determine that fault moment is 0.71s according to the phase voltage sudden change moment;
Step S2, chooses the faulted phase current Wave data (being designated as zero under establishing fault that the data point in moment occurs) of totally 800 points in each pick-up unit [0,69,0.73] interval, according to the definition of current break rate, calculates A i, result is as shown in table 1, and table 1 is singlephase earth fault simulation result;
Wherein, during 500 Ω stake resistance, fault phase current change quantity waveform along the line is as shown in Figure 4.
Step S3, section current break rate 2. meets A i>1.2, is judged as its fault section; Otherwise 1., 3. section is non-faulting section.
Table 1 singlephase earth fault simulation result

Claims (3)

1., based on a distribution network line fault Section Location for current break rate, it is characterized in that, said method comprising the steps of:
Step S1, determines fault phase and fault moment t f;
Step S2, chooses pick-up unit [t f-0.02s, t f+ 0.02s] interval interior 2N the faulted phase current Wave data put altogether, calculating current mutation rate A i, formula is as follows:
A I = 1 N Σ n = 0 N - 1 [ Δ i A ( n ) ] 2 / 1 N Σ n = - N - 1 [ Δ i A ( n ) ] 2
Δi A(n)=i 1A(n)-i 2A(n) n∈[-N,N-1]
In formula, i 1A(n), i 2An phase current sampling sequence that () is adjacent check point, it is that a periodic sampling is counted that fault is designated as zero, N under there is the data point in moment;
According to current break rate size, step S3, judges whether each section is fault section, namely works as A i> A i, set, be judged as fault section, otherwise be then non-faulting section;
Step 4, according to check point distributing position on the line, to travel through successively, until find out fault section.
2. distribution network line fault localization method according to claim 1, is characterized in that, described step S1 comprises:
After systems axiol-ogy to the generation of earth fault, according to phase voltage Changing Pattern, select fault phase, determine fault moment t according to phase voltage sudden change moment or arc-extinction device chugging moment f.
3. distribution network line fault localization method according to claim 1, is characterized in that, A in described step S3 isetfor the current break rate threshold value artificially arranged, get 1.2 ~ 1.5.
CN201510210539.4A 2015-04-28 2015-04-28 Distribution line fault section location method based on current break rate Active CN104793106B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510210539.4A CN104793106B (en) 2015-04-28 2015-04-28 Distribution line fault section location method based on current break rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510210539.4A CN104793106B (en) 2015-04-28 2015-04-28 Distribution line fault section location method based on current break rate

Publications (2)

Publication Number Publication Date
CN104793106A true CN104793106A (en) 2015-07-22
CN104793106B CN104793106B (en) 2017-08-25

Family

ID=53558076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510210539.4A Active CN104793106B (en) 2015-04-28 2015-04-28 Distribution line fault section location method based on current break rate

Country Status (1)

Country Link
CN (1) CN104793106B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158640A (en) * 2015-09-15 2015-12-16 江苏金智科技股份有限公司 Multi-power-source power supply system and fault positioning method based on GPS and current waveform
CN105403779A (en) * 2015-09-29 2016-03-16 昆明理工大学 Direct current line fault recognition method based on polar line current gradient sum
CN105914718A (en) * 2016-04-28 2016-08-31 长沙理工大学 Phase current sudden change based earth fault protection method
CN107290629A (en) * 2017-07-12 2017-10-24 国网青海省电力公司海南供电公司 A kind of 10KV low-voltage networks Earth design method
CN110632462A (en) * 2019-10-29 2019-12-31 深圳供电局有限公司 Small current grounding fault positioning method and system, computer equipment and medium
CN110912089A (en) * 2019-11-26 2020-03-24 华北电力大学 Flexible direct-current power distribution system protection method based on local current mutation quantity polarity
CN111308278A (en) * 2020-05-11 2020-06-19 国网江西省电力有限公司电力科学研究院 High-resistance fault direction detection method for resonant grounding system
CN111596170A (en) * 2020-05-27 2020-08-28 国网上海市电力公司 Fault diagnosis comprehensive positioning method for intelligent distribution network
CN111596171A (en) * 2020-05-27 2020-08-28 国网上海市电力公司 Intelligent distribution network fault diagnosis and positioning integrated system under artificial intelligence deep learning
CN112782528A (en) * 2020-12-31 2021-05-11 西安理工大学 Power distribution network fault section positioning method using PMU
CN114089123A (en) * 2021-12-06 2022-02-25 国网上海市电力公司 Fault positioning method, system and storage medium based on phase current time domain characteristics

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6454378A (en) * 1987-08-26 1989-03-01 Hitachi Cable Trouble section locating system for overhead power transmission line
US5629825A (en) * 1990-04-02 1997-05-13 Square D Company Apparatus and method for detecting a fault in a distributed line network
CN101907677A (en) * 2010-07-02 2010-12-08 华北电力大学 High voltage cable-overhead line hybrid line fault phase ranging method
CN102621449A (en) * 2012-03-16 2012-08-01 河南理工大学 Single phase ground fault section locating method in small current grounding system
CN103688434A (en) * 2011-07-21 2014-03-26 西门子公司 Fault identification and location in a power supply line which is fed from one side
CN104155582A (en) * 2014-08-29 2014-11-19 上海交通大学 Distribution network line fault section positioning method based on full-waveform information
CN104237738A (en) * 2014-08-29 2014-12-24 珠海威瀚科技发展有限公司 Distribution feeder single-phase grounding location system and location method
CN104237742A (en) * 2014-10-10 2014-12-24 华北电力大学 Method for improving fault location tolerance of power distribution network by utilizing fault current relative errors

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6454378A (en) * 1987-08-26 1989-03-01 Hitachi Cable Trouble section locating system for overhead power transmission line
US5629825A (en) * 1990-04-02 1997-05-13 Square D Company Apparatus and method for detecting a fault in a distributed line network
CN101907677A (en) * 2010-07-02 2010-12-08 华北电力大学 High voltage cable-overhead line hybrid line fault phase ranging method
CN103688434A (en) * 2011-07-21 2014-03-26 西门子公司 Fault identification and location in a power supply line which is fed from one side
US20140159740A1 (en) * 2011-07-21 2014-06-12 Siemens Aktiengesellschaft Fault identification and location in a power supply line which is fed from one side
CN102621449A (en) * 2012-03-16 2012-08-01 河南理工大学 Single phase ground fault section locating method in small current grounding system
CN104155582A (en) * 2014-08-29 2014-11-19 上海交通大学 Distribution network line fault section positioning method based on full-waveform information
CN104237738A (en) * 2014-08-29 2014-12-24 珠海威瀚科技发展有限公司 Distribution feeder single-phase grounding location system and location method
CN104237742A (en) * 2014-10-10 2014-12-24 华北电力大学 Method for improving fault location tolerance of power distribution network by utilizing fault current relative errors

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LAVON B.PAGE ET AL.: "Standard Deviation As an Alternative to Fuzziness in Fault Tree Models", 《IEEE TRANSACTIONS ON RELIABILITY》 *
郑顾平 等: "配网自动化系统中小电流接地故障区段定位方法", 《中国电机工程学报》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158640A (en) * 2015-09-15 2015-12-16 江苏金智科技股份有限公司 Multi-power-source power supply system and fault positioning method based on GPS and current waveform
CN105158640B (en) * 2015-09-15 2018-07-20 江苏金智科技股份有限公司 Multi-power supply system and Fault Locating Method based on GPS and current waveform
CN105403779B (en) * 2015-09-29 2019-01-04 昆明理工大学 A kind of DC line fault recognition methods based on polar curve current gradient sum
CN105403779A (en) * 2015-09-29 2016-03-16 昆明理工大学 Direct current line fault recognition method based on polar line current gradient sum
CN105914718A (en) * 2016-04-28 2016-08-31 长沙理工大学 Phase current sudden change based earth fault protection method
CN105914718B (en) * 2016-04-28 2018-04-10 长沙理工大学 A kind of earth-fault protection method based on difference of phase currents
CN107290629B (en) * 2017-07-12 2020-09-01 国网青海省电力公司海南供电公司 10KV low-voltage distribution network ground fault positioning method
CN107290629A (en) * 2017-07-12 2017-10-24 国网青海省电力公司海南供电公司 A kind of 10KV low-voltage networks Earth design method
CN110632462A (en) * 2019-10-29 2019-12-31 深圳供电局有限公司 Small current grounding fault positioning method and system, computer equipment and medium
CN110632462B (en) * 2019-10-29 2022-07-26 深圳供电局有限公司 Small current grounding fault positioning method and system, computer equipment and medium
CN110912089A (en) * 2019-11-26 2020-03-24 华北电力大学 Flexible direct-current power distribution system protection method based on local current mutation quantity polarity
CN111308278A (en) * 2020-05-11 2020-06-19 国网江西省电力有限公司电力科学研究院 High-resistance fault direction detection method for resonant grounding system
CN111596170A (en) * 2020-05-27 2020-08-28 国网上海市电力公司 Fault diagnosis comprehensive positioning method for intelligent distribution network
CN111596171A (en) * 2020-05-27 2020-08-28 国网上海市电力公司 Intelligent distribution network fault diagnosis and positioning integrated system under artificial intelligence deep learning
CN112782528A (en) * 2020-12-31 2021-05-11 西安理工大学 Power distribution network fault section positioning method using PMU
CN112782528B (en) * 2020-12-31 2023-07-18 西安理工大学 Power distribution network fault section positioning method by utilizing PMU
CN114089123A (en) * 2021-12-06 2022-02-25 国网上海市电力公司 Fault positioning method, system and storage medium based on phase current time domain characteristics

Also Published As

Publication number Publication date
CN104793106B (en) 2017-08-25

Similar Documents

Publication Publication Date Title
CN104793106A (en) Distribution network line fault section positioning method based on current break rate
CN106707084B (en) A kind of single-phase grounded malfunction in grounded system of low current determination method
CN103675605B (en) A kind of power distribution network earth fault line selection method based on the correlation analysis of fault-signal transient state
CN101201379B (en) Method for faulty indication and subsection of power system low current grounding
CN102944814A (en) Power distribution network single-phase earth fault locating method based on transient state
CN106771870A (en) A kind of distribution net work earthing fault localization method and system
CN102590703B (en) Single-phase grounding failure route selection method of resonance grounding system based on zero-sequence transient charge
CN106646139A (en) Method for positioning faults of power distribution network based on amplitude analysis of three-phase current
CN101232176A (en) Non-effective earthing distribution system fault locating method based on neutral point of transient traveling wave
CN103792465A (en) Power distribution network one-phase grounding fault location method based on zero sequence voltage
CN103018632B (en) Small current grounding system single-phase ground fault line selection method based on fisher information
CN104155582A (en) Distribution network line fault section positioning method based on full-waveform information
CN101975911B (en) Earth fault judging method for overhead line fault indicator
CN104166067A (en) Single-phase earth fault positioning detection method and device
CN101201380A (en) Method for faulty orientation and subsection of power system low current grounding
CN104181442A (en) Power distribution network single-phase earth fault section locating method based on correlation analysis
CN104764981B (en) Distribution line fault section location method based on normalized excursion degree
CN106093714A (en) The selection method of single-phase grounded malfunction in grounded system of low current circuit
CN105372616A (en) CVT state on-line test system and on-line evaluation method
CN104808112B (en) Distribution line fault section location method based on section instantaneous power
CN103558460A (en) Medium-voltage system arc fault detection device
CN104865496A (en) Distribution network line fault segment positioning method based on differential offset
CN104898017A (en) Euclidean distance-based distribution network line fault section locating method
CN108845233A (en) Electric distribution network overhead wire single-phase earthing double-end monitor localization method
CN105044550A (en) Distribution network line fault positioning method based on fault current discharge path

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant