CN112114214A - Transformer fault diagnosis method - Google Patents
Transformer fault diagnosis method Download PDFInfo
- Publication number
- CN112114214A CN112114214A CN202010934013.1A CN202010934013A CN112114214A CN 112114214 A CN112114214 A CN 112114214A CN 202010934013 A CN202010934013 A CN 202010934013A CN 112114214 A CN112114214 A CN 112114214A
- Authority
- CN
- China
- Prior art keywords
- data
- transformer
- fault diagnosis
- fault
- support vector
- 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 67
- 238000003745 diagnosis Methods 0.000 title claims abstract description 46
- 238000012360 testing method Methods 0.000 claims abstract description 39
- 238000012706 support-vector machine Methods 0.000 claims abstract description 35
- 241001206881 Myrmeleon inconspicuus Species 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 230000003449 preventive effect Effects 0.000 claims abstract description 13
- 238000012549 training Methods 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000010219 correlation analysis Methods 0.000 claims abstract description 4
- 241000257303 Hymenoptera Species 0.000 claims description 21
- 241000258923 Neuroptera Species 0.000 claims description 21
- 238000012544 monitoring process Methods 0.000 claims description 8
- 238000013021 overheating Methods 0.000 claims description 6
- 238000005295 random walk Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 241000282320 Panthera leo Species 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 239000002994 raw material Substances 0.000 claims 3
- 239000000284 extract Substances 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 description 1
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 description 1
- 239000005977 Ethylene Substances 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 229910002091 carbon monoxide Inorganic materials 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 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
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0027—General constructional details of gas analysers, e.g. portable test equipment concerning the detector
- G01N33/0036—General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0027—General constructional details of gas analysers, e.g. portable test equipment concerning the detector
- G01N33/0036—General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
- G01N33/0047—Organic compounds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0027—General constructional details of gas analysers, e.g. portable test equipment concerning the detector
- G01N33/0036—General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
- G01N33/005—H2
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/26—Oils; Viscous liquids; Paints; Inks
- G01N33/28—Oils, i.e. hydrocarbon liquids
-
- 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
-
- 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/2688—Measuring quality factor or dielectric loss, e.g. loss angle, or power factor
- G01R27/2694—Measuring dielectric loss, e.g. loss angle, loss factor or power factor
-
- 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/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
- G01R31/1263—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation
- G01R31/1281—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation of liquids or gases
-
- 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/62—Testing of transformers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Immunology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Biochemistry (AREA)
- Data Mining & Analysis (AREA)
- Combustion & Propulsion (AREA)
- Power Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Biomedical Technology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
The invention discloses a transformer fault diagnosis method, which comprises the steps of collecting historical preventive test data and fault data of a transformer in a power system and preprocessing the historical preventive test data and the fault data; analyzing the correlation among the preprocessed data items by utilizing a correlation analysis strategy to obtain a characteristic quantity set; optimizing the support vector machine by combining an ant lion optimization strategy to obtain an optimized support vector machine model; and carrying out fault diagnosis on the transformer by using the support vector machine model. The method combines a large amount of preventive test data and historical fault data of the transformer, analyzes the relevance among the data, extracts fault information in the transformer pre-test data, establishes a training data set and a test data set, and performs transformer fault diagnosis by using the support vector machine model optimized by the ant lion, so that the method has the advantages of high accuracy and short training time, and has the practical application value of engineering.
Description
Technical Field
The invention relates to the technical field of power systems, in particular to a transformer fault diagnosis method.
Background
Whether the power system can safely and stably operate and avoid accidents as much as possible depends on power equipment in the system to a great extent, and once the equipment fails, the safe operation state of the power grid is possibly damaged. Therefore, the safe operation of the power equipment is the first line of defense for avoiding major accidents of the power grid, the power transformer is used as the junction equipment of the power system, the cost is high, the maintenance cost is high, whether the operation condition of the transformer is stable, whether the transformer breaks down or not and whether the transformer is in a healthy state or not directly influences whether the power grid can safely and stably operate or not.
The transformer preventive test detects and evaluates each characteristic quantity aiming at the operation condition of the transformer, brings massive pre-test data, represents the instant condition of each characteristic parameter generated in the operation process of the transformer, can reflect the current overall operation condition of the transformer, and particularly reflects various fault conditions of the transformer under the current operation condition. The fault type can be judged through preventive test data of the transformer. Most of the current common fault diagnosis methods are intelligent algorithms based on operation data, and the problems that the data application amount is insufficient and the associated data cannot be selected efficiently exist. How to efficiently apply the transformer pre-test data and adopt the improved intelligent algorithm to realize the rapid and efficient specific diagnosis of the transformer fault is also the key problem of the current research.
Disclosure of Invention
This section is for the purpose of summarizing some aspects of embodiments of the invention and to briefly introduce some preferred embodiments. In this section, as well as in the abstract and the title of the invention of this application, simplifications or omissions may be made to avoid obscuring the purpose of the section, the abstract and the title, and such simplifications or omissions are not intended to limit the scope of the invention.
The present invention has been made in view of the above-mentioned conventional problems.
Therefore, the transformer fault diagnosis method provided by the invention can solve the problem of insufficient application of a large amount of preventive test data of the power enterprise operation information platform.
In order to solve the technical problems, the invention provides the following technical scheme: the method comprises the steps of collecting historical preventive test data and fault data of a transformer in a power system and preprocessing the historical preventive test data and the fault data; analyzing the correlation among the preprocessed data items by utilizing a correlation analysis strategy to obtain a characteristic quantity set; optimizing the support vector machine by combining an ant lion optimization strategy to obtain an optimized support vector machine model; and carrying out fault diagnosis on the transformer by using the support vector machine model.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the historical preventative test data comprises nameplate information, oil monitoring data and chromatographic monitoring data; the fault data comprises equipment numbers, fault time and fault types.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: preprocessing the data includes removing erroneous data and supplementing partially missing data by interpolation.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: respectively calculating the correlation between the data items by using a Spearman correlation coefficient and a Kendall correlation coefficient; counting, analyzing and calculating results and dividing various characteristic quantities into five types; and taking the occurrence frequency of the correlation as the standard of the characteristic quantity classification, and combining two correlation results to comprehensively analyze to obtain the characteristic quantity set for fault diagnosis.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the Spearman correlation coefficients include,
wherein R isiAnd SiIs the value grade of the observed value i,andis the average grade of the two variables and N is the total number of observations.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the Kendall correlation coefficients include,
wherein, P is the numerical logarithm with consistent arrangement relation, and n is the number of statistical vectors.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: screening and obtaining fault result data from each transformer pre-test data center and classifying according to the fault types; coding the fault result data, and dividing the data into a training set and a test set according to the proportion of 10: 1; initializing parameters of the ant lion optimization strategy; iteratively updating the positions of the ants and the ant lions based on the ant lions optimization strategy, calculating the fitness value of each ant and the ant lions by using a fitness function, and finally outputting the optimal (C, g) parameter value according to the end condition; and extracting the (C, g) parameter values as parameters of the support vector machine, constructing the support vector machine model by combining a support vector machine strategy and the training set, analyzing and calculating the test set by using the optimized support vector machine model, and outputting a fault diagnosis result.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the support vector machine model includes a model of a support vector machine,
K(xi,xj)=exp(-g||xi-xj||2)
wherein, omega represents hyperplane normal vector, C represents punishment factor, punishment degree of control error division sample, n represents sample number, xi represents relaxation factor, indicates allowable error division rate under linear inseparable condition, yiRepresents the sample output, and yi∈{-1,1},xiRepresenting the sample input, b representing the threshold, and g being a gaussian radial basis function parameter.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the fault types include normal operation, medium and low temperature overheating, high temperature overheating, low energy discharge and high energy discharge.
As a preferable aspect of the transformer fault diagnosis method of the present invention, wherein: the ant lion optimization strategy includes that ants seek food in a search space through random walk behavior, and ant lions capture the ants with traps, as follows,
xi=[0,cumsum(2r(t)-1)]
wherein x isiIs the ant i wandering position coordinate, cumsum represents the calculation cumulative sum, t represents the iteration, and r (t) represents the random value.
The invention has the beneficial effects that: the method is based on an electric power company integrated data platform and an operation and maintenance data platform, the relevance between each item of data and a fault is obtained through data mining analysis, a characteristic set is obtained, an improved intelligent algorithm is adopted to diagnose the fault of the transformer, and the method has important significance for the operation and maintenance of the transformer; on the other hand, the method combines a large amount of preventive test data and historical fault data of the transformer, analyzes the relevance among the data, extracts fault information in the transformer pre-test data, establishes a training data set and a test data set, and utilizes the ant lion optimized support vector machine model to diagnose the fault of the transformer, so that the method has the advantages of high accuracy and short training time, and has the practical engineering application value.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise. Wherein:
fig. 1 is a schematic flow chart of a transformer fault diagnosis method according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a ant lion optimized support vector machine model of the transformer fault diagnosis method according to an embodiment of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, specific embodiments accompanied with figures are described in detail below, and it is apparent that the described embodiments are a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making creative efforts based on the embodiments of the present invention, shall fall within the protection scope of the present invention.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those specifically described and will be readily apparent to those of ordinary skill in the art without departing from the spirit of the present invention, and therefore the present invention is not limited to the specific embodiments disclosed below.
Furthermore, reference herein to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one implementation of the invention. The appearances of the phrase "in one embodiment" in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments.
The present invention will be described in detail with reference to the drawings, wherein the cross-sectional views illustrating the structure of the device are not enlarged partially in general scale for convenience of illustration, and the drawings are only exemplary and should not be construed as limiting the scope of the present invention. In addition, the three-dimensional dimensions of length, width and depth should be included in the actual fabrication.
Meanwhile, in the description of the present invention, it should be noted that the terms "upper, lower, inner and outer" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation and operate, and thus, cannot be construed as limiting the present invention. Furthermore, the terms first, second, or third are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
The terms "mounted, connected and connected" in the present invention are to be understood broadly, unless otherwise explicitly specified or limited, for example: can be fixedly connected, detachably connected or integrally connected; they may be mechanically, electrically, or directly connected, or indirectly connected through intervening media, or may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Example 1
Referring to fig. 1 and 2, a transformer fault diagnosis method is provided for a first embodiment of the present invention, and includes:
s1: historical preventive test data and fault data of a transformer in a power system are collected and preprocessed. Wherein, it is required to be noted that:
historical preventative test data includes nameplate information, oil monitoring data, chromatographic monitoring data;
the fault data comprises equipment number, fault time and fault type;
the fault types comprise normal operation, medium and low temperature overheating, high temperature overheating, low energy discharge and high energy discharge;
the nameplate information comprises equipment type, equipment name, voltage, model and delivery date;
the oil monitoring information mainly reflects the electrical, chemical and physical properties of the insulating oil, including test date, appearance (oil), acid value, pH value, flash point (closed), breakdown voltage resistance, volume resistivity and dielectric loss factor;
the chromatographic monitoring data is mainly obtained by analyzing the concentration and the composition of gas dissolved in the insulating oil of the oil-immersed transformer and comprises test date, gas components and moisture;
the gas components include methane, ethane, ethylene, acetylene, total hydrocarbons, hydrogen, carbon monoxide, carbon dioxide.
Specifically, preprocessing the data includes:
and eliminating wrong data, and supplementing partial missing data by using an interpolation method.
S2: and analyzing the correlation among the preprocessed data items by using a correlation analysis strategy to obtain a characteristic quantity set. Referring to fig. 2, the steps to be described are:
respectively calculating the correlation between the data items by using a Spearman correlation coefficient and a Kendall correlation coefficient;
counting, analyzing and calculating results and dividing various characteristic quantities into five types;
and (3) taking the occurrence frequency of the correlation as the standard of characteristic quantity classification, and combining two correlation results to comprehensively analyze to obtain a characteristic quantity set for fault diagnosis.
Specifically, Spearman correlation coefficients include:
wherein R isiAnd SiIs the value grade of the observed value i,andis the average grade of the two variables, N is the total number of observations;
the Kendall correlation coefficients include,
wherein, P is the numerical logarithm with consistent arrangement relation, and n is the number of statistical vectors.
S3: and optimizing the support vector machine by combining the ant lion optimization strategy to obtain an optimized support vector machine model. It should be noted that the ant lion optimization strategy includes:
ants seek food in the search space by random walk behavior, while ant lions capture ants with traps, as follows,
xi=[0,cumsum(2r(t)-1)]
wherein x isiIs the coordinates of the walking position of the ant i, cumsum represents the calculation accumulated sum, t represents iteration, and r (t) represents a random value;
wherein rand is a random number over a uniform interval [0,1 ];
in the foraging process of ants, the position of each ant is stored in the matrix M under different statesAntThe method comprises the following steps:
wherein A isi,jRepresenting the position of the ith ant in the j dimension, n representing the number of the ants, and d representing the dimension of the parameter;
in the ALO algorithm, the position of each ant represents an attempted solution to the corresponding problem, and the matrix M is consideredAntThe positions of all ants generated in the optimization problem, i.e. all trial solutions corresponding to the problem, may be saved, and in the optimization process, the position of each ant (i.e. the trial solution) is evaluated using a fitness function, and the fitness value is saved in a matrix MOAIn (1), the following:
wherein f () represents a fitness function;
just as the ant positions correspond to the trial solution of the problem, in the ALO algorithm, the position of each lion represents a locally optimal solution of the problem, similar to the ant positions, using the matrix MAntLionAnd MOALThe position and fitness value of the ant lion are respectively stored:
wherein M isAntLionFor preserving the position of each ant lion, MOALFor storing the adaptation value, AL, corresponding to each ant lion positioni,jRepresenting the position of the ith ant lion in the j dimension, n representing the number of ant lions, d representing the dimension of a parameter, and f () representing a fitness function;
the random movement of each ant is limited by the search boundary of the target problem, once the random movement exceeds the search boundary to form an out-of-range problem, the obtained optimal solution loses practical significance, and in order to prevent the out-of-range problem, the random walking behavior of the ants is limited by adopting a maximum and minimum normalization mode:
wherein, ai、biRespectively is the minimum value and the maximum value of the ith parameter random walk,the evaluation formulas are respectively the minimum value and the maximum value of the ith parameter random walk in the tth iteration, and the evaluation formulas are shown as formulas (5-8) and (5-9):
wherein,indicates the position of the selected i-th ant lion in the t-th iteration, ct、dtThe minimum value and the maximum value of all variables of the ith dimension ant are expressed;
to obtain higher accuracy in the search process, the definition is:
wherein, T represents the maximum iteration times, w is a constant defined based on the current iteration condition, and the accuracy of the search can be adjusted;
in the process of random walking of ants, the ants randomly enter traps arranged on the ant lions, when the ants enter the traps and gradually reach the bottom of the pit to be caught by the ant lions, the process is the last stage of hunting, the ant lions enter the sand pit and slide to the ant lions to catch, in order to simulate the process, when the ants are more suitable than the ant lions catching the ant lions, the ants are caught by the ant lions, and then the positions of the ant lions are updated to the positions of the ants catching the ants for the last time; here, the update of the local optimal solution is realized, and is expressed by a logic language as follows:
for each ant, selecting one ant lion by Russian roulette, and recording as REIn the optimization process, the optimal solution obtained at each step is called elite lion, which has the best value in the problem, denoted as RAIn an iterative process, each ant will surround R simultaneouslyEAnd RARandom walk is performed, and the process is described as follows:
wherein,indicating the ant positions obtained by roulette in the t-th iteration;representing the random position of the elite lion in the t-th iteration.
When the adaptability of ants is stronger than that of ant lion REWhen it comes, the ant lion REThe position of (A) is replaced by an ant position, namely 'predation';
after searching all the ant lions, marking the ant lion with the best current fitness as the elite ant lion RA;
If the algorithm satisfies the stop condition, the current Elite lion RAThe position is the global optimal solution, and the output is the optimal solution of the problem;
if the algorithm does not meet the stop condition, a roulette is required.
S4: and carrying out fault diagnosis on the transformer by using a support vector machine model. It should be further noted that the step includes:
screening and obtaining fault result data from each transformer pre-test data center and classifying according to fault types;
coding fault result data, and dividing the data into a training set and a test set according to the proportion of 10: 1;
initializing parameters of the ant lion optimization strategy;
iteratively updating the positions of the ants and the ant lions based on an ant lions optimization strategy, calculating the fitness value of each ant and each ant lions by using a fitness function, and finally outputting the optimal (C, g) parameter value according to the end condition;
extracting (C, g) parameter values as parameters of a support vector machine, and constructing a support vector machine model by combining a strategy of the support vector machine and a training set;
and analyzing and calculating the test set by using the optimized support vector machine model, and outputting a fault diagnosis result.
Further, the support vector machine model includes:
K(xi,xj)=exp(-g||xi-xj||2)
wherein, omega represents hyperplane normal vector, C represents punishment factor, punishment degree of control error sample, n represents sample quantity, xi represents relaxation factor, which refers to allowance under linear inseparable conditionError rate of allowance, yiRepresents the sample output, and yi∈{-1,1},xiRepresenting the sample input, b representing the threshold, and g being a gaussian radial basis function parameter.
Preferably, in this embodiment, it should be further explained that the conventional fault diagnosis method for the transformer uses a self-contained spss decision tree for analysis, and the original method is not changed; the existing transformer fault diagnosis method based on the intelligent integrated algorithm integrates the traditional ratio method and the machine learning algorithm, and then carries out fault classification on each algorithm weight; the method is an improved support vector machine method, optimizes the parameters of the support vector machine and then classifies the faults, and the method has a difference in model selection, namely, optimizes the parameters of the support vector machine method to improve the accuracy of fault diagnosis.
Example 2
In order to better verify and explain the technical effects adopted in the method, the method selects the traditional transformer fault diagnosis method based on the intelligent integration algorithm and the method for carrying out comparison test, compares the test results by means of scientific demonstration, and verifies the real effect of the method.
The fault diagnosis accuracy of the traditional transformer fault diagnosis method based on the intelligent integrated algorithm is low, and in order to verify that the method has higher accuracy compared with the traditional method, the traditional method and the method are adopted to respectively carry out real-time measurement and comparison on transformers of the same type and operated by the same type of equipment in the embodiment.
And (3) testing environment: (1) selecting historical preventive test data and fault data of transformers in 100 groups of power systems, and respectively performing parameter processing by using two methods;
(2) MATLB software simulation, inputting the preprocessed 5 groups of parameters for testing;
the test results were as follows:
table 1: accuracy vs. data table.
Referring to table 1, it can be seen visually that under the test condition facing the same device, the fault identification accuracy of the conventional method is much lower than that of the method of the present invention, and the real effect of the method of the present invention is verified.
It should be noted that the above-mentioned embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, which should be covered by the claims of the present invention.
Claims (10)
1. A transformer fault diagnosis method is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
collecting historical preventive test data and fault data of a transformer in a power system and preprocessing the historical preventive test data and the fault data;
analyzing the correlation among the preprocessed data items by utilizing a correlation analysis strategy to obtain a characteristic quantity set;
optimizing the support vector machine by combining an ant lion optimization strategy to obtain an optimized support vector machine model;
and carrying out fault diagnosis on the transformer by using the support vector machine model.
2. The transformer fault diagnosis method according to claim 1, characterized in that: the historical preventative test data comprises nameplate information, oil monitoring data and chromatographic monitoring data;
the fault data comprises equipment numbers, fault time and fault types.
3. The transformer fault diagnosis method according to claim 2, characterized in that: preprocessing the data includes removing erroneous data and supplementing partially missing data by interpolation.
4. The transformer fault diagnosis method according to claim 3, characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
respectively calculating the correlation between the data items by using a Spearman correlation coefficient and a Kendall correlation coefficient;
counting, analyzing and calculating results and dividing various characteristic quantities into five types;
and taking the occurrence frequency of the correlation as the standard of the characteristic quantity classification, and combining two correlation results to comprehensively analyze to obtain the characteristic quantity set for fault diagnosis.
7. The transformer fault diagnosis method according to any one of claims 1 to 6, characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
screening and obtaining fault result data from each transformer pre-test data center and classifying according to the fault types;
coding the fault result data, and dividing the data into a training set and a test set according to the proportion of 10: 1;
initializing parameters of the ant lion optimization strategy;
iteratively updating the positions of the ants and the ant lions based on the ant lions optimization strategy, calculating the fitness value of each ant and the ant lions by using a fitness function, and finally outputting the optimal (C, g) parameter value according to the end condition;
extracting the (C, g) parameter values as parameters of the support vector machine, and combining a support vector machine strategy and the training set to construct the support vector machine model;
and analyzing and calculating the test set by using the optimized support vector machine model, and outputting a fault diagnosis result.
8. The transformer fault diagnosis method according to claim 7, characterized in that: the support vector machine model includes a model of a support vector machine,
K(xi,xj)=exp(-g||xi-xj||2)
wherein, omega represents a hyperplane normal vector, C represents a penalty factor, controls the penalty degree of the wrong divided samples, and n represents the number of the samplesQuantity xi denotes the relaxation factor, meaning the allowable error rate in the case of linear inseparability, yiRepresents the sample output, and yi∈{-1,1},xiRepresenting the sample input, b representing the threshold, and g being a gaussian radial basis function parameter.
9. The transformer fault diagnosis method according to claim 8, characterized in that: the fault types include normal operation, medium and low temperature overheating, high temperature overheating, low energy discharge and high energy discharge.
10. The transformer fault diagnosis method according to claim 9, characterized in that: the ant lion optimization strategy comprises the steps of,
ants seek food in the search space by random walk behavior, while lion captures the ants with traps, as follows,
xi=[0,cumsum(2r(t)-1)]
wherein x isiIs the ant i wandering position coordinate, cumsum represents the calculation cumulative sum, t represents the iteration, and r (t) represents the random value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010934013.1A CN112114214B (en) | 2020-09-08 | 2020-09-08 | Transformer fault diagnosis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010934013.1A CN112114214B (en) | 2020-09-08 | 2020-09-08 | Transformer fault diagnosis method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112114214A true CN112114214A (en) | 2020-12-22 |
CN112114214B CN112114214B (en) | 2022-12-27 |
Family
ID=73803403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010934013.1A Active CN112114214B (en) | 2020-09-08 | 2020-09-08 | Transformer fault diagnosis method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112114214B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112613584A (en) * | 2021-01-07 | 2021-04-06 | 国网上海市电力公司 | Fault diagnosis method, device, equipment and storage medium |
CN113253709A (en) * | 2021-06-07 | 2021-08-13 | 江苏中车数字科技有限公司 | Health diagnosis method and device suitable for rail transit vehicle |
CN114019298A (en) * | 2021-09-28 | 2022-02-08 | 中电华创(苏州)电力技术研究有限公司 | PCC-SVM-based generator rotor turn-to-turn short circuit online monitoring method |
CN114152825A (en) * | 2021-11-16 | 2022-03-08 | 国网北京市电力公司 | Fault diagnosis method and device of transformer and fault diagnosis system of transformer |
CN116990625A (en) * | 2023-09-27 | 2023-11-03 | 国网江苏省电力有限公司电力科学研究院 | Function switching system and method of intelligent quick-checking device of distribution transformer |
CN117233540A (en) * | 2023-11-15 | 2023-12-15 | 广东电网有限责任公司江门供电局 | Metering pipeline fault detection method and system based on deep learning |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101709995A (en) * | 2009-12-24 | 2010-05-19 | 浙江大学 | Methods for vibration online monitoring and fault diagnosis of power transformer |
CN105203876A (en) * | 2015-09-15 | 2015-12-30 | 云南电网有限责任公司电力科学研究院 | Transformer on-line monitoring state assessment method utilizing support vector machine and correlation analysis |
CN106501693A (en) * | 2016-12-08 | 2017-03-15 | 贵州电网有限责任公司电力科学研究院 | A kind of Diagnosis Method of Transformer Faults based on fuzzy Boltzmann machine |
WO2017146877A1 (en) * | 2016-02-03 | 2017-08-31 | General Electric Company | Systems and methods for monitoring and diagnosing transformer health |
CN107169514A (en) * | 2017-05-05 | 2017-09-15 | 清华大学 | The method for building up of diagnosing fault of power transformer model |
CN108663582A (en) * | 2017-11-30 | 2018-10-16 | 全球能源互联网研究院有限公司 | A kind of fault diagnosis method and system of transformer |
CN108680814A (en) * | 2018-08-29 | 2018-10-19 | 国网河北省电力有限公司电力科学研究院 | A kind of various dimensions running state of transformer evaluation method |
CN108872771A (en) * | 2018-07-11 | 2018-11-23 | 云南电网有限责任公司电力科学研究院 | A kind of transformer preventive trial equivalent method |
CN109061453A (en) * | 2018-08-02 | 2018-12-21 | 国网福建省电力有限公司 | Substation's disconnecting link secondary circuit failure prediction technique of meter and related coefficient |
CN109245943A (en) * | 2018-10-19 | 2019-01-18 | 北京奇艺世纪科技有限公司 | A kind of Fault Locating Method and device |
CN109767036A (en) * | 2018-12-28 | 2019-05-17 | 北京航空航天大学 | Support vector machines failure prediction method based on the optimization of adaptive ant lion |
CN110006649A (en) * | 2018-12-24 | 2019-07-12 | 湖南科技大学 | A kind of Method for Bearing Fault Diagnosis based on improvement ant lion algorithm and support vector machines |
CN110129832A (en) * | 2019-06-21 | 2019-08-16 | 广西大学 | A kind of Multipurpose Optimal Method of aluminium electrolysis process tank voltage |
CN110766175A (en) * | 2019-10-25 | 2020-02-07 | 长沙理工大学 | Pitch system fault detection method and device based on optimal interval distribution machine |
CN110969262A (en) * | 2019-12-03 | 2020-04-07 | 广东电网有限责任公司 | Transformer fault diagnosis method |
CN111237134A (en) * | 2020-01-14 | 2020-06-05 | 上海电力大学 | Offshore double-fed wind driven generator fault diagnosis method based on GRA-LSTM-stacking model |
US20200202181A1 (en) * | 2018-12-19 | 2020-06-25 | Netskope, Inc. | Multi-label classification of text documents |
CN111537836A (en) * | 2020-05-15 | 2020-08-14 | 国网山东省电力公司济宁供电公司 | Automatic power distribution network fault diagnosis method and system based on wave recording data |
-
2020
- 2020-09-08 CN CN202010934013.1A patent/CN112114214B/en active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101709995A (en) * | 2009-12-24 | 2010-05-19 | 浙江大学 | Methods for vibration online monitoring and fault diagnosis of power transformer |
CN105203876A (en) * | 2015-09-15 | 2015-12-30 | 云南电网有限责任公司电力科学研究院 | Transformer on-line monitoring state assessment method utilizing support vector machine and correlation analysis |
WO2017146877A1 (en) * | 2016-02-03 | 2017-08-31 | General Electric Company | Systems and methods for monitoring and diagnosing transformer health |
CN106501693A (en) * | 2016-12-08 | 2017-03-15 | 贵州电网有限责任公司电力科学研究院 | A kind of Diagnosis Method of Transformer Faults based on fuzzy Boltzmann machine |
CN107169514A (en) * | 2017-05-05 | 2017-09-15 | 清华大学 | The method for building up of diagnosing fault of power transformer model |
CN108663582A (en) * | 2017-11-30 | 2018-10-16 | 全球能源互联网研究院有限公司 | A kind of fault diagnosis method and system of transformer |
CN108872771A (en) * | 2018-07-11 | 2018-11-23 | 云南电网有限责任公司电力科学研究院 | A kind of transformer preventive trial equivalent method |
CN109061453A (en) * | 2018-08-02 | 2018-12-21 | 国网福建省电力有限公司 | Substation's disconnecting link secondary circuit failure prediction technique of meter and related coefficient |
CN108680814A (en) * | 2018-08-29 | 2018-10-19 | 国网河北省电力有限公司电力科学研究院 | A kind of various dimensions running state of transformer evaluation method |
CN109245943A (en) * | 2018-10-19 | 2019-01-18 | 北京奇艺世纪科技有限公司 | A kind of Fault Locating Method and device |
US20200202181A1 (en) * | 2018-12-19 | 2020-06-25 | Netskope, Inc. | Multi-label classification of text documents |
CN110006649A (en) * | 2018-12-24 | 2019-07-12 | 湖南科技大学 | A kind of Method for Bearing Fault Diagnosis based on improvement ant lion algorithm and support vector machines |
CN109767036A (en) * | 2018-12-28 | 2019-05-17 | 北京航空航天大学 | Support vector machines failure prediction method based on the optimization of adaptive ant lion |
CN110129832A (en) * | 2019-06-21 | 2019-08-16 | 广西大学 | A kind of Multipurpose Optimal Method of aluminium electrolysis process tank voltage |
CN110766175A (en) * | 2019-10-25 | 2020-02-07 | 长沙理工大学 | Pitch system fault detection method and device based on optimal interval distribution machine |
CN110969262A (en) * | 2019-12-03 | 2020-04-07 | 广东电网有限责任公司 | Transformer fault diagnosis method |
CN111237134A (en) * | 2020-01-14 | 2020-06-05 | 上海电力大学 | Offshore double-fed wind driven generator fault diagnosis method based on GRA-LSTM-stacking model |
CN111537836A (en) * | 2020-05-15 | 2020-08-14 | 国网山东省电力公司济宁供电公司 | Automatic power distribution network fault diagnosis method and system based on wave recording data |
Non-Patent Citations (3)
Title |
---|
章彬: "基于数理统计的变压器油色谱联合预警策略", 《高压电器》 * |
简春菲: "基于ALO-SVM的模拟电路故障诊断研究", 《自动化技术与应用》 * |
董翔: "电力变压器故障类型与关键状态量关联规则分析", 《浙江电力》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112613584A (en) * | 2021-01-07 | 2021-04-06 | 国网上海市电力公司 | Fault diagnosis method, device, equipment and storage medium |
CN113253709A (en) * | 2021-06-07 | 2021-08-13 | 江苏中车数字科技有限公司 | Health diagnosis method and device suitable for rail transit vehicle |
CN113253709B (en) * | 2021-06-07 | 2021-09-21 | 江苏中车数字科技有限公司 | Health diagnosis method and device suitable for rail transit vehicle |
CN114019298A (en) * | 2021-09-28 | 2022-02-08 | 中电华创(苏州)电力技术研究有限公司 | PCC-SVM-based generator rotor turn-to-turn short circuit online monitoring method |
CN114019298B (en) * | 2021-09-28 | 2023-12-05 | 中电华创(苏州)电力技术研究有限公司 | On-line monitoring method for turn-to-turn short circuit of generator rotor based on PCC-SVM |
CN114152825A (en) * | 2021-11-16 | 2022-03-08 | 国网北京市电力公司 | Fault diagnosis method and device of transformer and fault diagnosis system of transformer |
CN114152825B (en) * | 2021-11-16 | 2023-11-14 | 国网北京市电力公司 | Transformer fault diagnosis method and device and transformer fault diagnosis system |
CN116990625A (en) * | 2023-09-27 | 2023-11-03 | 国网江苏省电力有限公司电力科学研究院 | Function switching system and method of intelligent quick-checking device of distribution transformer |
CN116990625B (en) * | 2023-09-27 | 2023-12-08 | 国网江苏省电力有限公司电力科学研究院 | Function switching system and method of intelligent quick-checking device of distribution transformer |
CN117233540A (en) * | 2023-11-15 | 2023-12-15 | 广东电网有限责任公司江门供电局 | Metering pipeline fault detection method and system based on deep learning |
CN117233540B (en) * | 2023-11-15 | 2024-02-20 | 广东电网有限责任公司江门供电局 | Metering pipeline fault detection method and system based on deep learning |
Also Published As
Publication number | Publication date |
---|---|
CN112114214B (en) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112114214B (en) | Transformer fault diagnosis method | |
Huang | A new data mining approach to dissolved gas analysis of oil-insulated power apparatus | |
Huang et al. | Evolving wavelet networks for power transformer condition monitoring | |
Malik et al. | Fuzzy reinforcement learning based intelligent classifier for power transformer faults | |
CN106896219B (en) | The identification of transformer sub-health state and average remaining lifetime estimation method based on Gases Dissolved in Transformer Oil data | |
CN117323578B (en) | Intelligent control method and system for infrared therapeutic apparatus | |
Das et al. | Accurate identification of transformer faults from dissolved gas data using recursive feature elimination method | |
CN116010884A (en) | Fault diagnosis method of SSA-LightGBM oil-immersed transformer based on principal component analysis | |
CN111060652A (en) | Method for predicting concentration of dissolved gas in transformer oil based on long-term and short-term memory network | |
CN115018512A (en) | Electricity stealing detection method and device based on Transformer neural network | |
CN114358116A (en) | Oil-immersed transformer fault diagnosis method and system and readable storage medium | |
CN116842337A (en) | Transformer fault diagnosis method based on LightGBM (gallium nitride based) optimal characteristics and COA-CNN (chip on board) model | |
CN117290404A (en) | Method and system for rapidly searching and practical main distribution network fault processing method | |
CN113721000A (en) | Method and system for detecting abnormity of dissolved gas in transformer oil | |
CN112817954A (en) | Missing value interpolation method based on multi-method ensemble learning | |
Leauprasert et al. | Intelligent machine learning techniques for condition assessment of power transformers | |
CN118446118A (en) | Highway electromechanical system fault prediction method and device based on neural network, electronic equipment and storage medium | |
CN117786516A (en) | Oil immersed transformer fault diagnosis method based on multi-source data fusion | |
CN117454113A (en) | Power transformer fault early warning method based on CNN-BiLSTM and RF model | |
CN113112067A (en) | Method for establishing TFRI weight calculation model | |
Kim et al. | Anomaly pattern detection in streaming data based on the transformation to multiple binary-valued data streams | |
CN116150604B (en) | Transformer fault diagnosis method and device and electronic equipment | |
CN116662862A (en) | Transformer fault early warning method based on prediction and diagnosis of content of dissolved gas in oil | |
Fosić et al. | Network traffic verification based on a public dataset for IDS systems and machine learning classification algorithms | |
CN110426612B (en) | Two-stage transformer oiled paper insulation time domain dielectric response characteristic quantity optimization method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |