CN112114214A - Transformer fault diagnosis method - Google Patents

Transformer fault diagnosis method Download PDF

Info

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
Application number
CN202010934013.1A
Other languages
Chinese (zh)
Other versions
CN112114214B (en
Inventor
张霖
李诗勇
姜龙
杨超
薛静
张丽娟
杨俊秋
靳斌
李俊文
申峻
郝越峰
施艳
汪德军
丁健
谢荣斌
王瑞果
赵世钦
刘波
刘栖
陈湫林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guizhou Power Grid Co Ltd
Original Assignee
Guizhou Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guizhou Power Grid Co Ltd filed Critical Guizhou Power Grid Co Ltd
Priority to CN202010934013.1A priority Critical patent/CN112114214B/en
Publication of CN112114214A publication Critical patent/CN112114214A/en
Application granted granted Critical
Publication of CN112114214B publication Critical patent/CN112114214B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0027General constructional details of gas analysers, e.g. portable test equipment concerning the detector
    • G01N33/0036General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0027General constructional details of gas analysers, e.g. portable test equipment concerning the detector
    • G01N33/0036General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
    • G01N33/0047Organic compounds
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0027General constructional details of gas analysers, e.g. portable test equipment concerning the detector
    • G01N33/0036General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
    • G01N33/005H2
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/26Oils; Viscous liquids; Paints; Inks
    • G01N33/28Oils, i.e. hydrocarbon liquids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/26Measuring 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/2688Measuring quality factor or dielectric loss, e.g. loss angle, or power factor
    • G01R27/2694Measuring dielectric loss, e.g. loss angle, loss factor or power factor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/12Testing 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/1227Testing 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/1263Testing 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/1281Testing 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/62Testing of transformers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411Classification 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial 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

Transformer fault diagnosis method
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,
Figure BDA0002671263060000021
wherein R isiAnd SiIs the value grade of the observed value i,
Figure BDA0002671263060000022
and
Figure BDA0002671263060000023
is 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,
Figure BDA0002671263060000024
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,
Figure BDA0002671263060000031
Figure BDA0002671263060000032
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:
Figure BDA0002671263060000061
wherein R isiAnd SiIs the value grade of the observed value i,
Figure BDA0002671263060000062
and
Figure BDA0002671263060000063
is the average grade of the two variables, N is the total number of observations;
the Kendall correlation coefficients include,
Figure BDA0002671263060000064
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;
Figure BDA0002671263060000065
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:
Figure BDA0002671263060000071
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:
Figure BDA0002671263060000072
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:
Figure BDA0002671263060000073
Figure BDA0002671263060000074
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:
Figure BDA0002671263060000081
wherein, ai、biRespectively is the minimum value and the maximum value of the ith parameter random walk,
Figure BDA0002671263060000082
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):
Figure BDA0002671263060000083
Figure BDA0002671263060000084
wherein,
Figure BDA0002671263060000085
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:
Figure BDA0002671263060000086
Figure BDA0002671263060000087
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:
Figure BDA0002671263060000091
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:
Figure BDA0002671263060000092
wherein,
Figure BDA0002671263060000093
indicating the ant positions obtained by roulette in the t-th iteration;
Figure BDA0002671263060000094
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:
Figure BDA0002671263060000095
Figure BDA0002671263060000101
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.
Figure BDA0002671263060000102
Figure BDA0002671263060000111
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.
5. The transformer fault diagnosis method according to claim 4, characterized in that: the Spearman correlation coefficients include,
Figure FDA0002671263050000011
wherein R isiAnd SiIs the value grade of the observed value i,
Figure FDA0002671263050000012
and
Figure FDA0002671263050000013
is the average grade of the two variables and N is the total number of observations.
6. The transformer fault diagnosis method according to claim 5, characterized in that: the Kendall correlation coefficients include,
Figure FDA0002671263050000014
wherein, P is the numerical logarithm with consistent arrangement relation, and n is the number of statistical vectors.
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,
Figure FDA0002671263050000021
Figure FDA0002671263050000022
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.
CN202010934013.1A 2020-09-08 2020-09-08 Transformer fault diagnosis method Active CN112114214B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (18)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
章彬: "基于数理统计的变压器油色谱联合预警策略", 《高压电器》 *
简春菲: "基于ALO-SVM的模拟电路故障诊断研究", 《自动化技术与应用》 *
董翔: "电力变压器故障类型与关键状态量关联规则分析", 《浙江电力》 *

Cited By (11)

* Cited by examiner, † Cited by third party
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