CN113933749A - Method, device, equipment and storage medium for detecting high-resistance earth fault of power distribution network - Google Patents
Method, device, equipment and storage medium for detecting high-resistance earth fault of power distribution network Download PDFInfo
- Publication number
- CN113933749A CN113933749A CN202111210767.3A CN202111210767A CN113933749A CN 113933749 A CN113933749 A CN 113933749A CN 202111210767 A CN202111210767 A CN 202111210767A CN 113933749 A CN113933749 A CN 113933749A
- Authority
- CN
- China
- Prior art keywords
- distribution network
- frequency data
- power distribution
- data
- phase
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 230000009466 transformation Effects 0.000 claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 38
- 238000004364 calculation method Methods 0.000 claims description 30
- 238000000819 phase cycle Methods 0.000 claims description 9
- 238000009825 accumulation Methods 0.000 description 10
- 239000003990 capacitor Substances 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000004020 conductor Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/52—Testing for short-circuits, leakage current or ground faults
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
The application discloses a method, a device, equipment and a storage medium for detecting a high-resistance ground fault of a power distribution network, wherein the method comprises the following steps: performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data; acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data; when the energy spectrum data are larger than a first threshold and the energy coefficient is smaller than a second threshold, acquiring a plurality of second voltage data of the power distribution network to be detected; performing discrete wavelet transform on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data; and judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not, and if so, judging that the distribution network to be detected has a high-resistance grounding fault. The method solves the technical problems that misjudgment and misjudgment are easy to occur during detection, and further the power supply reliability is influenced in the existing detection method for the high-resistance ground fault of the power distribution network.
Description
Technical Field
The application relates to the technical field of power distribution networks, in particular to a method, a device, equipment and a storage medium for detecting a high-resistance ground fault of a power distribution network.
Background
Under the influence of factors such as natural environment, overhead distance of lines and the like, single-phase high-resistance grounding faults of non-ideal conductors occur frequently in the power distribution network, such as the fact that leads fall on grasslands, roads, sand lands, ponds and the like. The high resistance ground accounts for approximately 10% of the total number of ground faults. If the high-resistance grounding fault exists for a long time, the system can generate a new grounding point due to overvoltage generated by the fault, so that the accident is further expanded, even a fire disaster is caused by electric arcs accompanied by the fault, the personal safety and the property safety are threatened, and therefore the high-resistance grounding fault needs to be identified in time.
Although the existing detection for the high-resistance ground fault can be carried out to a certain degree, the normal operation and the high-resistance ground fault cannot be effectively distinguished during the detection, and misjudgment are easy to occur, so that the power supply reliability is influenced.
Disclosure of Invention
In view of this, the application provides a method, an apparatus, a device and a storage medium for detecting a high-resistance ground fault of a power distribution network, which solve the technical problem that the existing method for detecting a high-resistance ground fault of a power distribution network is prone to misjudgment and misjudgment during detection, and further affects power supply reliability.
The application provides a method for detecting a high-resistance ground fault of a power distribution network in a first aspect, which comprises the following steps:
performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data;
acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data;
when the energy spectrum data are larger than the first threshold value and the energy coefficient is smaller than the second threshold value, acquiring a plurality of second voltage data of the power distribution network to be detected;
performing discrete wavelet transform on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data;
and judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not, if so, judging that the distribution network to be detected has a high-resistance ground fault, wherein the preset accumulated sum is the accumulated sum of the high-frequency data with the preset number of layers in each second high-frequency data.
Optionally, the obtaining of the energy spectrum data and the first threshold corresponding to the first high-frequency data, and the energy coefficient and the second threshold corresponding to the first low-frequency data specifically includes:
calculating corresponding energy spectrum data based on the first high-frequency data, and calculating corresponding energy coefficients based on the first low-frequency data;
and calculating a corresponding first threshold according to the energy spectrum data, and calculating a second threshold according to the energy coefficient.
Optionally, a calculation formula corresponding to the first threshold is as follows:
ρ＝(1+α).max(es)；
where ρ is a first threshold, α is an expansion coefficient, and es is energy spectrum data.
Optionally, a calculation formula corresponding to the second threshold is as follows:
τ＝(1-α).min(ec)；
where τ is the second threshold, α is the expansion coefficient, and ec is the energy coefficient.
Optionally, when it is determined that the power distribution network to be detected has a high-resistance ground fault, the method further includes:
acquiring three-phase voltage and three-phase current corresponding to any two adjacent feeder line nodes in a power distribution network to be detected;
constructing a left positioning matrix and a right positioning matrix for fault positioning based on the three-phase voltage and the three-phase current;
acquiring a fault type corresponding to the power distribution network to be detected and a fault positioning distance calculation formula corresponding to the fault type;
and calculating the fault positioning distance corresponding to the power distribution network to be detected by taking one of the two feeder line nodes as a reference point based on the fault positioning distance calculation formula and according to the left positioning matrix and the right positioning matrix.
Optionally, the calculation formula of the left positioning matrix is:
wherein ,
for left location matrix LFtThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, x represents a phase sequence of a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
the t-th current in the x-phase for the feed node p,
the t current in phase x for feed node q.
Optionally, the calculation formula of the right positioning matrix is:
wherein ,
for right positioning matrix RItThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, L tableShowing the line length between the feeder nodes p and q, wherein x represents the phase sequence and is a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
for the t-th voltage at phase x at the feed node p,
for the t-th voltage at phase x at feed node q,
the t-th current in phase a for the feed node p,
the t-th current in phase b for the feed node p,
the t-th current in phase c for the feed node p.
This application second aspect provides a detection device of distribution network high resistance ground fault, includes:
the first transformation unit is used for performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data;
a first obtaining unit, configured to obtain energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data;
the second obtaining unit is used for obtaining a plurality of second voltage data of the power distribution network to be detected when the energy spectrum data are larger than the first threshold value and the energy coefficient is smaller than the second threshold value;
the second transformation unit is used for carrying out discrete wavelet transformation on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data;
and the judging unit is used for judging whether the preset accumulated sum of the second high-frequency data is greater than a preset judgment index or not, and if so, judging that the distribution network to be detected has a high-resistance ground fault, wherein the preset accumulated sum is the accumulated sum of the high-frequency data with the preset number of layers in each second high-frequency data.
The third aspect of the application provides a detection device for a high-resistance ground fault of a power distribution network, which comprises a processor and a memory;
the memory is used for storing program codes and transmitting the program codes to the processor;
the processor is configured to execute any one of the methods for detecting a high impedance earth fault in a power distribution network according to the first aspect according to instructions in the program code.
A fourth aspect of the present application provides a storage medium, which is used to store program codes, where the program codes are used to execute any one of the methods for detecting a high impedance ground fault of a power distribution network according to the first aspect.
According to the technical scheme, the method has the following advantages:
the application provides a method for detecting a high-resistance ground fault of a power distribution network, which comprises the steps of firstly carrying out discrete wavelet transformation on acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data, then acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data, then when the energy spectrum data is larger than a first threshold value and the energy coefficient is smaller than a second threshold value, a plurality of second voltage data of the distribution network to be detected are obtained, then discrete wavelet transform is carried out on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data, whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not is judged, if yes, the distribution network to be detected is judged to have high-resistance grounding fault, the preset accumulation sum is the accumulation sum of the high-frequency data with the preset number of layers in each second high-frequency data.
The inventor discovers that the misjudgment and the missed judgment of the existing detection method are caused by the fact that high-frequency data are generated in normal operations such as switching operation, capacitor access and heavy load switching, and the existing detection method only uses the high-frequency data to detect the high-resistance ground fault, so that the normal operation and the high-resistance ground fault cannot be effectively distinguished in detection. Therefore, the technical problem that misjudgment and misjudgment are easy to occur during detection so as to influence the power supply reliability in the conventional method for detecting the high-resistance ground fault of the power distribution network is solved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise.
Fig. 1 is a schematic flowchart of a first embodiment of a method for detecting a high-resistance ground fault of a power distribution network according to an embodiment of the present application;
fig. 2 is a schematic flowchart of a method for detecting a high-resistance ground fault of a power distribution network according to a second embodiment of the present disclosure;
fig. 3 is a schematic structural diagram of a detection apparatus for a high-resistance ground fault of a power distribution network in an embodiment of the present application.
Detailed Description
The embodiment of the application provides a method, a device, equipment and a storage medium for detecting a high-resistance ground fault of a power distribution network, and solves the technical problem that misjudgment and misjudgment easily occur during detection so as to influence power supply reliability in the existing method for detecting the high-resistance ground fault of the power distribution network.
In order to make the technical solutions of the present application better understood, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The first aspect of the embodiments of the present application provides an embodiment of a method for detecting a high-resistance ground fault of a power distribution network.
Referring to fig. 1, a flowchart of a method for detecting a high-resistance ground fault of a power distribution network according to a first embodiment of the present application is shown.
The method for detecting the high-resistance ground fault of the power distribution network in the embodiment comprises the following steps:
It is understood that the first voltage data in the present embodiment is 1 cycle voltage signal. When discrete wavelet transform is carried out on the first voltage data, m layers of discrete wavelet transform processing are carried out on the first voltage data, and m first high-frequency data h are obtained respectivelyiAnd m first low frequency data li. Specifically, when m-layer discrete wavelet transform processing is performed, for each layer of data, a high-pass filter is used for acquiring high-frequency data in each layer, and a low-pass filter is used for acquiring low-frequency data.
In this embodiment, both the low frequency data and the high frequency data are used as the judgment basis, and therefore, after the first high frequency data and the first low frequency data are obtained, the energy spectrum data and the first threshold corresponding to the first high frequency data, and the energy coefficient and the second threshold corresponding to the first low frequency data are calculated.
And 103, when the energy spectrum data are larger than the first threshold and the energy coefficient is smaller than the second threshold, acquiring a plurality of second voltage data of the power distribution network to be detected.
And when the energy spectrum data is larger than the first threshold and the energy coefficient is smaller than the second threshold, indicating that a potential high-resistance ground fault exists, otherwise, outputting a low-resistance ground fault and finishing the judgment.
It can be understood that one second voltage data is 1 cycle voltage signal, but in this embodiment, after the high-resistance ground fault possibly existing in the power distribution network to be detected is obtained based on the first voltage data, in order to detect accuracy, the second voltage data is used for judging secondary detection.
In order to effectively identify that the high-resistance grounding fault is a long-period fault generated by intermittent arc light, when the second voltage data is acquired, the acquisition of the second voltage data is carried out after a period of time is delayed after the acquisition of the first voltage data. Specifically, the delay time may be set to 150ms, or may be set to another value, which may be selected by those skilled in the art as needed.
The number n of the second voltage data may be set to 10, 20, etc., which is not specifically limited in this embodiment.
And 104, performing discrete wavelet transform on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data.
After the n second voltage data are obtained, discrete wavelet transform is performed on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data.
It is understood that the discrete wavelet transform of the second voltage data may also adopt the processing manner of the first voltage data, so that after the discrete wavelet transform is performed on the second voltage data, each second voltage data corresponds to m second high frequency data
And 105, judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index, and if so, judging that the distribution network to be detected has a high-resistance grounding fault.
In this embodiment, the preset accumulated sum is an accumulated sum of the high frequency data of the preset number of layers in each second high frequency data. I.e., the sum of the w-th layer high frequency data of the m second high frequency data. It can be understood that, in this embodiment, the value of w is 4, and the corresponding calculation formula of the preset cumulative sum is:
wherein hif is a predetermined sum; h isj4J is the high frequency data of the 4 th layer in the jth second high frequency data, and the value of j is 1, 2, …, m.
And judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index, if so, indicating that the distribution network to be detected has a high-resistance ground fault, otherwise, indicating that the distribution network to be detected has a low-resistance ground fault.
In the embodiment, firstly, discrete wavelet transform is performed on the acquired first voltage data of the distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data, then acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data, then when the energy spectrum data is larger than a first threshold value and the energy coefficient is smaller than a second threshold value, a plurality of second voltage data of the distribution network to be detected are obtained, then discrete wavelet transform is carried out on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data, whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not is judged, if yes, the distribution network to be detected is judged to have high-resistance grounding fault, the preset accumulation sum is the accumulation sum of the high-frequency data with the preset number of layers in each second high-frequency data. In the embodiment, the low-frequency data and the high-frequency data are simultaneously used as judgment bases, so that normal operations such as voltage drop, switching operation, capacitor access, heavy load switching and the like and high-resistance grounding faults can be effectively distinguished, and misjudgment and missed judgment of the high-resistance grounding faults are avoided. Therefore, the technical problem that misjudgment and misjudgment are easy to occur during detection so as to influence the power supply reliability in the conventional method for detecting the high-resistance ground fault of the power distribution network is solved.
The above is a first embodiment of a method for detecting a high-resistance ground fault of a power distribution network provided by the embodiment of the present application, and the following is a second embodiment of a method for detecting a high-resistance ground fault of a power distribution network provided by the embodiment of the present application.
Referring to fig. 2, a flowchart of a method for detecting a high-resistance ground fault of a power distribution network according to a second embodiment of the present application is shown.
The method for detecting the high-resistance ground fault of the power distribution network in the embodiment comprises the following steps:
It is understood that the description of step 201 is the same as that of step 101 in the first embodiment, and reference may be specifically made to the description of step 101, which is not described herein again.
It can be understood that the energy spectrum data es is calculated by the formula:
the calculation formula of the energy coefficient ec is as follows:
The calculation formula corresponding to the first threshold is as follows:
ρ＝(1+α).max(es)；
where ρ is a first threshold, α is an expansion coefficient, es is energy spectrum data, and max () represents a maximum value before the determination time, that is, the first threshold is an adaptive value and is applicable to different failure environments.
The calculation formula corresponding to the second threshold is as follows:
τ＝(1-α).min(ec)；
in the formula, τ is a second threshold, α is an expansion coefficient, ec is an energy coefficient, and min () represents a minimum value before the determination time, that is, the second threshold is an adaptive value and is suitable for different fault environments.
And 204, when the energy spectrum data are larger than a first threshold and the energy coefficient is smaller than a second threshold, acquiring a plurality of second voltage data of the power distribution network to be detected.
It is understood that the description of step 204 is the same as that of step 103 in the first embodiment, and reference may be specifically made to the description of step 103, which is not described herein again.
And 205, performing discrete wavelet transform on each second voltage data to obtain a plurality of corresponding second high-frequency data.
It is understood that the description of step 205 is the same as that of step 104 in the first embodiment, and reference may be specifically made to the description of step 104, which is not described herein again.
And step 206, judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index, and if so, judging that the distribution network to be detected has a high-resistance grounding fault.
The preset accumulation sum is the accumulation sum of the high-frequency data with the preset number of layers in each second high-frequency data.
It is understood that the description of step 206 is the same as that of step 105 in the first embodiment, and reference may be specifically made to the description of step 105, which is not described herein again.
And step 207, when it is determined that the high-resistance ground fault occurs to the power distribution network to be detected, acquiring three-phase voltage and three-phase current corresponding to any two adjacent feeder line nodes in the power distribution network to be detected.
The feeder nodes (tie switch, branch switch, section switch) of the power distribution network are provided with measuring terminals for collecting three-phase voltage and three-phase current of any two adjacent feeder nodes (p, q), respectively
x represents a phase sequence and is a, b and c; ip and Vp are the current and voltage of the P node; iq and Vq are current and voltage of a q node; t represents a sampling number, the value is 1-T, and T is the total number of sampling points.
And 208, constructing a left positioning matrix and a right positioning matrix for fault positioning based on the three-phase voltage and the three-phase current.
It is understood that the formula for the left positioning matrix is:
wherein ,
for left location matrix LFtThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, x represents a phase sequence of a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
the t-th current in the x-phase for the feed node p,
the t current in phase x for feed node q.
The calculation formula of the right positioning matrix is as follows:
wherein ,
for right positioning matrix RItThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, L represents the line length between feeder nodes p and q, and x represents the phase sequence and is a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
for the t-th voltage at phase x at the feed node p,
for the t-th voltage at phase x at feed node q,
the t-th current in phase a for the feed node p,
the t-th current in phase b for the feed node p,
the t-th current in phase c for the feed node p.
And 209, acquiring a fault type corresponding to the power distribution network to be detected and a fault positioning distance calculation formula corresponding to the fault type.
For obtaining the fault type of the power distribution network to be detected, the single-phase and two-phase ground faults can be obtained by judging the maximum voltage phase
And zero sequence voltage thereof
The phase relationship therebetween. If it is
Lags behind
Single-phase earth fault occurs with lag behind
The one phase of (a) is a ground phase;
ahead of
A two-phase ground fault occurs and,
the two phases are non-failure phases and the remaining two phases are failure phases. And for the three-phase short-circuit fault, detecting whether the three-phase current is increased at the same time and exceeds a preset value, and if so, determining that the three-phase short-circuit fault occurs.
It will be appreciated that the general formula for calculating the fault location distance is:
wherein x1, x2 and x3 represent fault phases, can be a, b and c, x1-x2-x3 represents fault combination, and lox1-x2-x3The corresponding fault location distance is combined for the fault,
for right positioning matrix RItThe t-th element of the medium x1 phase,
left side is decidedBit matrix LFtThe t-th element of the medium x1 phase,
for right positioning matrix RItThe t-th element of the medium x2 phase,
for left location matrix LFtThe t-th element of the medium x2 phase,
for right positioning matrix RItThe t-th element of the medium x3 phase,
for left location matrix LFtThe t-th element of the medium x3 phase.
And for different fault types (namely fault combinations), corresponding to different fault positioning distance calculation formulas. For a three-phase power distribution network, the following fault types are provided:
(1) a phase is grounded in a single phase; x1 ═ a, x2 ═ 0, and x3 ═ 0 (i.e., where 1 is a and the other two are 0).
(2) b, grounding the phase single phase; x1 ═ b, x2 ═ 0, and x3 ═ 0 (i.e., where 1 is b and the other two are 0).
(3) c, grounding the phase single phase; x1 ═ c, x2 ═ 0, and x3 ═ 0 (i.e., where 1 is c and the other two are 0).
(4) a-b are short-circuited to ground; x1 ═ a, x2 ═ b, and x3 ═ 0 (i.e., where 1 is 0 and the other two are a, b).
(5) a-c are short-circuited to ground; x1 ═ a, x2 ═ c, and x3 ═ 0 (i.e., where 1 is 0 and the other two are a, c).
(6) c-b short circuit to ground; x1 ═ c, x2 ═ b, and x3 ═ 0 (i.e., where 1 is 0 and the other two are c, b).
(7) a-b-c three-phase short circuit; x1 ═ a, x2 ═ b, and x3 ═ c (i.e., 3 are a, b, and c, respectively).
It can be understood that when the values of x1, x2 and x3 are 0, the corresponding values are 0
And
also 0.
For example: when a-phase single-phase grounding occurs, the fault location distance x1 ═ a is as follows:
when bc-phase short-circuit fault occurs, x1 is b, x2 is c, and the fault location distance is as follows:
and step 210, calculating the fault location distance corresponding to the power distribution network to be detected by taking one of the two feeder line nodes as a reference point according to the left location matrix and the right location matrix based on a fault location distance calculation formula.
It can be understood that the calculation of the fault location distance in the present embodiment is not only applicable to the calculation of the high resistance ground fault, but also applicable to the calculation of the corresponding fault location distance in the case of the low resistance ground fault.
Compared with the prior art, the embodiment has the following advantages:
1) meanwhile, the low-frequency component and the high-frequency component are used as judgment features, daily operations such as voltage drop, switching operation, capacitor access, heavy load switching and the like can be effectively distinguished, and misjudgment of high-resistance grounding faults is avoided.
2) After the potential high-resistance grounding fault is judged, secondary judgment is carried out after a period of time is delayed, and long-period faults of the high-resistance grounding fault due to intermittent arc light can be effectively identified
3) The positioning algorithm (namely calculation of fault positioning distance) is not only suitable for low-resistance earth faults and high-resistance earth faults, but also suitable for different fault types, is not influenced by fault resistance and fault phases, and has universality.
In the embodiment, firstly, discrete wavelet transform is performed on the acquired first voltage data of the distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data, then acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data, then when the energy spectrum data is larger than a first threshold value and the energy coefficient is smaller than a second threshold value, a plurality of second voltage data of the distribution network to be detected are obtained, then discrete wavelet transform is carried out on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data, whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not is judged, if yes, the distribution network to be detected is judged to have high-resistance grounding fault, the preset accumulation sum is the accumulation sum of the high-frequency data with the preset number of layers in each second high-frequency data. In the embodiment, the low-frequency data and the high-frequency data are simultaneously used as judgment bases, so that normal operations such as voltage drop, switching operation, capacitor access, heavy load switching and the like and high-resistance grounding faults can be effectively distinguished, and misjudgment and missed judgment of the high-resistance grounding faults are avoided. Therefore, the technical problem that misjudgment and misjudgment are easy to occur during detection so as to influence the power supply reliability in the conventional method for detecting the high-resistance ground fault of the power distribution network is solved.
The second aspect of the embodiments of the present application provides an embodiment of a device for detecting a high-resistance ground fault of a power distribution network.
Referring to fig. 3, a schematic structural diagram of a detection apparatus for detecting a high-resistance ground fault of a power distribution network in an embodiment of the present application is shown.
In this embodiment, a detection apparatus for a distribution network high resistance ground fault includes:
the first transformation unit is used for performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data;
the first acquisition unit is used for acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data;
the second acquisition unit is used for acquiring a plurality of second voltage data of the power distribution network to be detected when the energy spectrum data is larger than the first threshold and the energy coefficient is smaller than the second threshold;
the second transformation unit is used for carrying out discrete wavelet transformation on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data;
and the judging unit is used for judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judging index or not, and if so, judging that the high-resistance ground fault occurs to the power distribution network to be detected, wherein the preset accumulated sum is the accumulated sum of the high-frequency data with the preset number of layers in each second high-frequency data.
Further, the first obtaining unit specifically includes:
the first calculating subunit is used for calculating corresponding energy spectrum data and corresponding energy coefficients based on the first high-frequency data and the first low-frequency data;
and the second calculating subunit is used for calculating a corresponding first threshold and a corresponding second threshold according to the energy spectrum data.
Further, the calculation formula corresponding to the first threshold is as follows:
ρ＝(1+α).max(es)；
where ρ is a first threshold, α is an expansion coefficient, and es is energy spectrum data.
Specifically, the calculation formula corresponding to the second threshold is as follows:
τ＝(1-α).min(ec)；
where τ is the second threshold, α is the expansion coefficient, and ec is the energy coefficient.
Further, the power distribution network high-resistance ground fault detection device in this embodiment further includes:
the third acquisition unit is used for acquiring three-phase voltage and three-phase current corresponding to any two adjacent feeder line nodes in the power distribution network to be detected;
the building unit is used for building a left positioning matrix and a right positioning matrix for fault positioning based on the three-phase voltage and the three-phase current;
the fourth acquisition unit is used for acquiring the fault type corresponding to the power distribution network to be detected and the fault positioning distance calculation formula corresponding to the fault type;
and the calculating unit is used for calculating the fault positioning distance corresponding to the power distribution network to be detected by taking one of the two feeder line nodes as a reference point according to the left positioning matrix and the right positioning matrix based on a fault positioning distance calculation formula.
Optionally, the calculation formula of the left positioning matrix is:
wherein ,
for left location matrix LFtThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, x represents a phase sequence of a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
the t-th current in the x-phase for the feed node p,
the t current in phase x for feed node q.
Further, the calculation formula of the right positioning matrix is as follows:
wherein ,
for right positioning matrix RItThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, L represents the line length between feeder nodes p and q, and x represents the phase sequence and is a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
for the t-th voltage at phase x at the feed node p,
for the t-th voltage at phase x at feed node q,
the t-th current in phase a for the feed node p,
the t-th current in phase b for the feed node p,
the t-th current in phase c for the feed node p.
In the embodiment, firstly, discrete wavelet transform is performed on the acquired first voltage data of the distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data, then acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data, then when the energy spectrum data is larger than a first threshold value and the energy coefficient is smaller than a second threshold value, a plurality of second voltage data of the distribution network to be detected are obtained, then discrete wavelet transform is carried out on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data, whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not is judged, if yes, the distribution network to be detected is judged to have high-resistance grounding fault, the preset accumulation sum is the accumulation sum of the high-frequency data with the preset number of layers in each second high-frequency data. In the embodiment, the low-frequency data and the high-frequency data are simultaneously used as judgment bases, so that normal operations such as voltage drop, switching operation, capacitor access, heavy load switching and the like and high-resistance grounding faults can be effectively distinguished, and misjudgment and missed judgment of the high-resistance grounding faults are avoided. Therefore, the technical problem that misjudgment and misjudgment are easy to occur during detection so as to influence the power supply reliability in the conventional method for detecting the high-resistance ground fault of the power distribution network is solved.
A third aspect of embodiments of the present application provides an embodiment of a device for detecting a high-resistance ground fault of a power distribution network.
A detection device for a high-resistance ground fault of a power distribution network comprises a processor and a memory; the memory is used for storing the program codes and transmitting the program codes to the processor; the processor is configured to execute the method for detecting a high impedance earth fault of a power distribution network according to the first aspect.
A fourth aspect of embodiments of the present application provides an embodiment of a storage medium.
A storage medium for storing program code for performing the method for detecting a high impedance earth fault of a power distribution network of the first aspect.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of a unit is only one logical functional division, and there may be other divisions when the actual implementation is performed, for example, a plurality of units or components may be combined or may be integrated into another grid network to be installed, or some features may be omitted or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.
Claims (10)
1. A detection method for a high-resistance ground fault of a power distribution network is characterized by comprising the following steps:
performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data;
acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data;
when the energy spectrum data are larger than the first threshold value and the energy coefficient is smaller than the second threshold value, acquiring a plurality of second voltage data of the power distribution network to be detected;
performing discrete wavelet transform on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data;
and judging whether the preset accumulated sum of the second high-frequency data is larger than a preset judgment index or not, if so, judging that the distribution network to be detected has a high-resistance ground fault, wherein the preset accumulated sum is the accumulated sum of the high-frequency data with the preset number of layers in each second high-frequency data.
2. The method for detecting the high-resistance ground fault of the power distribution network according to claim 1, wherein acquiring energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data specifically comprises:
calculating corresponding energy spectrum data based on the first high-frequency data, and calculating corresponding energy coefficients based on the first low-frequency data;
and calculating a corresponding first threshold according to the energy spectrum data, and calculating a second threshold according to the energy coefficient.
3. The method for detecting the high-resistance ground fault of the power distribution network according to claim 2, wherein the calculation formula corresponding to the first threshold is as follows:
ρ＝(1+α).max(es)；
where ρ is a first threshold, α is an expansion coefficient, and es is energy spectrum data.
4. The method for detecting the high-resistance ground fault of the power distribution network according to claim 2, wherein the calculation formula corresponding to the second threshold is as follows:
τ＝(1-α).min(ec)；
where τ is the second threshold, α is the expansion coefficient, and ec is the energy coefficient.
5. The method for detecting the high-resistance ground fault of the power distribution network according to claim 1, wherein when it is determined that the power distribution network to be detected has the high-resistance ground fault, the method further comprises the following steps:
acquiring three-phase voltage and three-phase current corresponding to any two adjacent feeder line nodes in a power distribution network to be detected;
constructing a left positioning matrix and a right positioning matrix for fault positioning based on the three-phase voltage and the three-phase current;
acquiring a fault type corresponding to the power distribution network to be detected and a fault positioning distance calculation formula corresponding to the fault type;
and calculating the fault positioning distance corresponding to the power distribution network to be detected by taking one of the two feeder line nodes as a reference point based on the fault positioning distance calculation formula and according to the left positioning matrix and the right positioning matrix.
6. The method for detecting the high-resistance ground fault of the power distribution network according to claim 5, wherein the calculation formula of the left positioning matrix is as follows:
wherein ,
for left location matrix LFtThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, x represents a phase sequence of a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
the t-th current in the x-phase for the feed node p,
the t current in phase x for feed node q.
7. The method for detecting the high-resistance ground fault of the power distribution network according to claim 5, wherein the calculation formula of the right positioning matrix is as follows:
wherein ,
for right positioning matrix RItThe t-th element of the phase of (1), wherein
Comprises the following steps:
in the formula, L represents the line length between feeder nodes p and q, and x represents the phase sequence and is a, b and c; z is a radical ofxa，zxb，zxcRepresenting the impedances between the x-phase and a, b, c, respectively, t representing the sample number,
for the t-th voltage at phase x at the feed node p,
for the t-th voltage at phase x at feed node q,
the t-th current in phase a for the feed node p,
the t-th current in phase b for the feed node p,
the t-th current in phase c for the feed node p.
8. The utility model provides a detection apparatus for distribution network high resistance ground fault which characterized in that includes:
the first transformation unit is used for performing discrete wavelet transformation on the acquired first voltage data of the power distribution network to be detected to obtain a plurality of first high-frequency data and first low-frequency data;
a first obtaining unit, configured to obtain energy spectrum data and a first threshold corresponding to the first high-frequency data, and an energy coefficient and a second threshold corresponding to the first low-frequency data;
the second obtaining unit is used for obtaining a plurality of second voltage data of the power distribution network to be detected when the energy spectrum data are larger than the first threshold value and the energy coefficient is smaller than the second threshold value;
the second transformation unit is used for carrying out discrete wavelet transformation on each second voltage data to obtain a plurality of second high-frequency data corresponding to each second voltage data;
and the judging unit is used for judging whether the preset accumulated sum of the second high-frequency data is greater than a preset judgment index or not, and if so, judging that the distribution network to be detected has a high-resistance ground fault, wherein the preset accumulated sum is the accumulated sum of the high-frequency data with the preset number of layers in each second high-frequency data.
9. The detection equipment for the high-resistance ground fault of the power distribution network is characterized by comprising a processor and a memory;
the memory is used for storing program codes and transmitting the program codes to the processor;
the processor is configured to execute the method for detecting a high impedance earth fault of a power distribution network according to any one of claims 1 to 7 according to instructions in the program code.
10. A storage medium for storing a program code for executing the method for detecting a high impedance earth fault of a power distribution network according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111210767.3A CN113933749B (en) | 2021-10-18 | 2021-10-18 | Method, device, equipment and storage medium for detecting high-resistance ground fault of power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111210767.3A CN113933749B (en) | 2021-10-18 | 2021-10-18 | Method, device, equipment and storage medium for detecting high-resistance ground fault of power distribution network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113933749A true CN113933749A (en) | 2022-01-14 |
| CN113933749B CN113933749B (en) | 2023-06-16 |
Family
ID=79279871
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111210767.3A Active CN113933749B (en) | 2021-10-18 | 2021-10-18 | Method, device, equipment and storage medium for detecting high-resistance ground fault of power distribution network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113933749B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU4547593A (en) * | 1992-07-10 | 1994-01-31 | RMIT Innovation Ltd | High impedance fault detector |
| CN101187687A (en) * | 2007-12-21 | 2008-05-28 | 清华大学 | High resistance earthing fault detection method based on transient traveling wave |
| CN101673930A (en) * | 2008-07-24 | 2010-03-17 | 通用电气公司 | Arc detection using discrete wavelet transforms |
| CN104297636A (en) * | 2014-10-23 | 2015-01-21 | 广东电网有限责任公司电力科学研究院 | Pole-mounted distribution substation traveling wave detection method |
| US20190128944A1 (en) * | 2017-10-27 | 2019-05-02 | Siemens Aktiengesellschaft | Method and detection device for detecting a high-impedance ground fault in an electrical energy supply network with a grounded neutral point |
| CN111487507A (en) * | 2020-06-01 | 2020-08-04 | 南方电网调峰调频发电有限公司西部检修试验分公司 | Power distribution network high-resistance fault identification method based on wavelet packet energy ratio |
| CN112255492A (en) * | 2020-09-07 | 2021-01-22 | 西安理工大学 | Power distribution network single-phase grounding high-resistance fault identification method under strong noise background |
| CN112505474A (en) * | 2020-10-30 | 2021-03-16 | 国网山东省电力公司聊城供电公司 | Fault positioning method and device based on phase current increment |
-
2021
- 2021-10-18 CN CN202111210767.3A patent/CN113933749B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU4547593A (en) * | 1992-07-10 | 1994-01-31 | RMIT Innovation Ltd | High impedance fault detector |
| CN101187687A (en) * | 2007-12-21 | 2008-05-28 | 清华大学 | High resistance earthing fault detection method based on transient traveling wave |
| CN101673930A (en) * | 2008-07-24 | 2010-03-17 | 通用电气公司 | Arc detection using discrete wavelet transforms |
| CN104297636A (en) * | 2014-10-23 | 2015-01-21 | 广东电网有限责任公司电力科学研究院 | Pole-mounted distribution substation traveling wave detection method |
| US20190128944A1 (en) * | 2017-10-27 | 2019-05-02 | Siemens Aktiengesellschaft | Method and detection device for detecting a high-impedance ground fault in an electrical energy supply network with a grounded neutral point |
| CN111487507A (en) * | 2020-06-01 | 2020-08-04 | 南方电网调峰调频发电有限公司西部检修试验分公司 | Power distribution network high-resistance fault identification method based on wavelet packet energy ratio |
| CN112255492A (en) * | 2020-09-07 | 2021-01-22 | 西安理工大学 | Power distribution network single-phase grounding high-resistance fault identification method under strong noise background |
| CN112505474A (en) * | 2020-10-30 | 2021-03-16 | 国网山东省电力公司聊城供电公司 | Fault positioning method and device based on phase current increment |
Non-Patent Citations (2)
| Title |
|---|
| 曾志辉等: "基于变异系数与高阶累积量的小电流接地故障选线", 《电力***保护与控制》, vol. 48, no. 13, pages 99 - 109 * |
| 朱晓娟等: "基于小波能量矩的高阻接地故障检测方法", 《电力自动化设备》, vol. 36, no. 12, pages 161 - 168 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113933749B (en) | 2023-06-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2814661C (en) | Apparatus and method for fault detection and location determination | |
| Jafarian et al. | High-speed superimposed-based protection of series-compensated transmission lines | |
| Musa et al. | A covariance indices based method for fault detection and classification in a power transmission system during power swing | |
| EP3482472A1 (en) | A method and system for locating a fault in a mixed power transmission line | |
| CN111983510B (en) | Single-phase ground fault phase selection method and system based on phase voltage and current abrupt change | |
| Zubić et al. | Speed and security improvements of distance protection based on Discrete Wavelet and Hilbert transform | |
| CN111537837A (en) | Method and system for positioning small current ground fault of power distribution network | |
| Akorede et al. | Wavelet transform based algorithm for high-impedance faults detection in distribution feeders | |
| CN105914718B (en) | A kind of earth-fault protection method based on difference of phase currents | |
| Pandakov et al. | Misoperation analysis of steady-state and transient methods on earth fault locating in compensated distribution networks | |
| Kujur et al. | Detection of high impedance fault in distribution system considering distributed generation | |
| CN111969575A (en) | Arc extinction method, device and equipment for single-phase earth fault of power distribution network and storage medium | |
| de Alvarenga Ferreira et al. | A novel high impedance arcing fault detection based on the discrete wavelet transform for smart distribution grids | |
| CN104062555B (en) | The discrimination method of distribution line high resistance earthing fault characteristic harmonics | |
| CN104090211B (en) | A kind of online test method of distribution line high resistance earthing fault | |
| Chatterjee et al. | Dual use line relays to improve power swing deblocking function | |
| Dhend et al. | Efficient fault diagnosis in smart grid using non conventional mother wavelet function | |
| CN106950465B (en) | Electric leakage detection and selection method and system for mine electric power network | |
| Radhakrishnan et al. | Location of high impedance faults using smart meters in distribution systems with DGs, power electronic loads and Electric Arc furnaces | |
| Patel | Superimposed components of Lissajous pattern based feature extraction for classification and localization of transmission line faults | |
| CN113933749B (en) | Method, device, equipment and storage medium for detecting high-resistance ground fault of power distribution network | |
| CN110221182B (en) | Distributed small-current grounding line selection method and line selection system | |
| CN107632237A (en) | A kind of fault line selection method for single-phase-to-ground fault based on wavelet character amount correlation after failure | |
| US20220252644A1 (en) | Fault detection in an electric power system | |
| CN108254650B (en) | Quick judgment method for single-phase earth fault of substation bus |
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 |