CN111144431B - Transformer fault diagnosis method based on CBBO-SVM - Google Patents
Transformer fault diagnosis method based on CBBO-SVM Download PDFInfo
- Publication number
- CN111144431B CN111144431B CN201811312408.7A CN201811312408A CN111144431B CN 111144431 B CN111144431 B CN 111144431B CN 201811312408 A CN201811312408 A CN 201811312408A CN 111144431 B CN111144431 B CN 111144431B
- Authority
- CN
- China
- Prior art keywords
- habitat
- cbbo
- svm
- fault diagnosis
- species
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000003745 diagnosis Methods 0.000 title claims abstract description 31
- 230000000739 chaotic effect Effects 0.000 claims abstract description 26
- 238000012706 support-vector machine Methods 0.000 claims abstract description 21
- 238000005457 optimization Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 230000000694 effects Effects 0.000 claims abstract description 4
- 238000007781 pre-processing Methods 0.000 claims abstract description 3
- 238000013508 migration Methods 0.000 claims description 26
- 230000005012 migration Effects 0.000 claims description 26
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 15
- 238000013507 mapping Methods 0.000 claims description 12
- 239000007789 gas Substances 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 8
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 238000005291 chaos (dynamical) Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 claims description 3
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 claims description 3
- 239000005977 Ethylene Substances 0.000 claims description 3
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 claims description 3
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 3
- 239000001569 carbon dioxide Substances 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000010891 electric arc Methods 0.000 claims description 3
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 claims description 3
- 229910052739 hydrogen Inorganic materials 0.000 claims description 3
- 239000001257 hydrogen Substances 0.000 claims description 3
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 238000013021 overheating Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 10
- 241000894007 species Species 0.000 description 18
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 241000254158 Lampyridae Species 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/08—Computing arrangements based on specific mathematical models using chaos models or non-linear system models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Nonlinear Science (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Housings And Mounting Of Transformers (AREA)
Abstract
The invention discloses a CBBO-SVM-based transformer fault diagnosis method, which comprises the following specific steps of step 1, collecting a sample set S = { (x) of an oil-immersed transformer with class labels 1 ,x 2 ,x 3 ,...x n ),(y 1 ,y 2 ,y 3 ,...,x m ) }; step 2, preprocessing and normalizing the acquired data: normalizing both input and output quantities to [0,1]To achieve the purpose of convenient calculation; step 3, utilizing parameter c of CBBO to support vector machine i And σ i Optimizing; step 4, utilizing the parameter c obtained after the optimization in the step 3 i And σ i And establishing a multi-stage SVM model, and classifying by using the sample set data so as to achieve the effect of fault diagnosis of the oil-immersed transformer. The method can optimize parameters of the SVM (support vector machine) algorithm by using the CBBO (chaotic biophysics) algorithm, and effectively improve the accuracy of classification.
Description
Technical Field
The invention belongs to the technical field of transformer fault online monitoring, and particularly relates to a CBBO-SVM-based transformer fault diagnosis method.
Background
Nowadays, the demand for electricity in the whole country is rapidly increased, the power system is gradually advanced to the linkage direction of a large power grid with high voltage grade and high automation, and the stable operation of the power transformer is not only beneficial to promoting the extra-high voltage construction, but also is closely related to the life of the power system, medical treatment, industry and even consumers. Due to the complex structure of the transformer and the interactive influence of various links such as manufacturing, installation, operation and maintenance, maintenance and the like, the transformer is easy to fail in the operation process. Oil-immersed transformers represent a significant proportion of all transformer classes. Therefore, it is necessary to diagnose a fault in the oil-filled transformer.
The traditional diagnosis method of the oil-immersed transformer mainly relates to the following steps: the three-ratio method, the great defense triangle, the gas graph method, the characteristic gas method, etc., which are well established, but are continuously developed today, the diagnostic methods gradually expose many disadvantages, such as: low accuracy, failure to judge the coexistence of multiple faults, incomplete ratio code combination and the like. People add intelligent algorithms such as BP neural network, bayesian algorithm, firefly algorithm and the like to the traditional diagnosis method, and although the accuracy is improved to a certain extent, new problems occur, such as: the required fault data is more, the iteration times are high, the calculation time is long, and the like.
Therefore, the invention provides a transformer fault diagnosis method based on CBBO-SVM (chaotic biophysical optimization support vector machine), which can obtain optimal parameters on the premise of lower iteration number and greatly improve the diagnosis accuracy.
Disclosure of Invention
The invention aims to provide a CBBO-SVM-based transformer fault diagnosis method, which can optimize parameters of an SVM (support vector machine) algorithm by using a CBBO (chaotic biophysics) algorithm and effectively improve the classification accuracy.
The technical scheme adopted by the invention is that the transformer fault diagnosis method based on the CBBO-SVM is implemented according to the following steps:
step 1, collected sample set S = { (x) of oil-immersed transformer with class label 1 ,x 2 ,x 3 ,...x n ),(y 1 ,y 2 ,y 3 ,...,x m ) Of which 50% are training samples, 50% are testing samples, x i Representing the properties of the sample, including 7 properties of hydrogen, methane, ethane, ethylene, acetylene, carbon monoxide and carbon dioxide, y i The category labels 1, 2, 3, 4, 5 and 6 are represented and respectively correspond to 6 states of normal state, medium and low temperature overheating, high energy discharge, low energy discharge and arc discharge;
step 2, carrying out pretreatment and normalization treatment on the collected data: normalizing the input quantity and the output quantity to [0,1] by the following formula (1) so as to achieve the purpose of convenient calculation;
in the above formula, X max 、X min Respectively is collected H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 Maximum and minimum values of the concentration of the gas, X i Is represented by H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 The concentration of the gas;
step 3, utilizing parameter c of CBBO to support vector machine i And σ i Optimizing;
step 4, utilizing the parameter c obtained after the optimization in the step 3 i And σ i And establishing a multi-stage SVM model, and classifying by using the sample set data so as to achieve the effect of fault diagnosis of the oil-immersed transformer.
The present invention is also characterized in that,
step 3 is specifically implemented according to the following steps:
step 3.1, SVM parameter c i And σ i After the initialization operation is performed, c is i And σ i As each habitat P i Is vector of fitness X i Then generating an initial population W i 1 ;
Step 3.2, calculating initial population W i 1 Population habitat suitability index f i 1 Obtaining f from the formula (2) i 1 ;
In the above formula, i is the number of the habitat variable; n is the total number of habitat variables; j represents each habitat variable P i A level of (d); c. C j As a habitat variable P i The total number of stages; beta is a ij Is a multiple regression coefficient;
ω ij is a weighting factor derived from the frequency distribution of habitat variables; a is a constant coefficient;
step 3.3, initial population W i 1 Carrying out migration operation, and calculating the migration rate lambda and the migration rate mu according to the formula (3);
in the above formula, I is the maximum mobility; e is the maximum migration rate; lambda is the mobility rate; mu is an migration rate; s. the max The maximum number of species that can be accommodated by the habitat; s is the number of species in the habitat;
after the migration operation, a new species W is formed according to the migration in and out i 2 Then calculating the suitability index f of the population habitat i 2 ;
Step 3.4, performing variation operation on the population after the migration operation, wherein the variation rate m (x) of the migration habitat is along with the variation rate P of the species s Is increased, i.e.:
in the above formula: m (x) is habitat variation rate; m is max The maximum rate of variation; p s The probability of the species corresponding to the species number S in the habitat is shown; p max Is P s Maximum value of (d);
P s relationship to the mobility λ and the mobility μ:
randomly generating new species W in random habitats based on mutation operations i 3 Instead of the original species W i 2 Calculating the suitability index f of the population habitat i 3 ;
Step 3.5, adding the C chaos theory and utilizing the subsection LThe g-istic chaotic mapping equation (6) updates the population to generate a new initialization fitness vector f i 4 :
Wherein x is i Is a chaotic variable; k is the number of chaotic iterations; in step 3.5, when μ =4, the mapping equation (6) is in a chaotic state.
In order to ensure the diversity of species, chaotic search is carried out on the basis of Logistic chaotic mapping, and the optimal solution in the new fitness solution is recorded as J best To J with best Chaotic search is performed on the basis of:
J n =J best +ωx i (7)
wherein x is i Is the solution of the chaotic mapping equation; omega is an adjustment coefficient;
for the obtained optimal solution J n And J best Comparing, taking the optimum of the two, repeating step 3.2, and updating the output parameter c according to the requirement i And σ i 。
Step 4 is specifically implemented according to the following steps:
step 4.1, firstly, kernel function selection is carried out;
step 4.2, taking training samples accounting for 50% of the total samples as model input for verification;
step 4.3, using c obtained in step 3.5 i And σ i Introducing and obtaining a new objective function f (omega, sigma) i ):
Wherein c is a penalty factor; sigma i Is a relaxation factor;
step 4.4, solving by utilizing a Lagrange multiplier method, and combining an objective function and a constraint condition formula (10)
y i (ωx i +b)≥1-σ i (10)
solving omega and b to obtain a final decision function;
and 4.5, obtaining 6 decision functions through 6 layers of training, and establishing a 6-level support vector machine transformer fault diagnosis model.
The kernel function in step 4.4 is a Gaussian kernel function K (X, X) i )。
Gaussian kernel function K (X, X) i ) The method specifically comprises the following steps:
where ξ is a parameter that controls the height of the gaussian kernel function.
The invention has the beneficial effects that:
(1) According to the transformer fault diagnosis method based on the CBBO-SVM, the BBO algorithm is optimized by using the C chaos theory, and the formed CBBO algorithm greatly improves the parameter optimization capability;
(2) According to the transformer fault diagnosis method based on CBBO-SVM, CBBO is used for optimizing SVM parameters, the principle is simple, and classification performance of SVM can be improved on the premise of ensuring less iteration times;
(3) According to the transformer fault diagnosis method based on the CBBO-SVM, the CBBO-SVM algorithm is combined with gas detection in transformer oil, so that the accuracy of fault diagnosis of the oil-immersed transformer is greatly improved.
Drawings
FIG. 1 is a schematic diagram of a cosine mobility model involved in a CBBO-SVM-based transformer fault diagnosis method of the invention;
FIG. 2 is a schematic diagram of a multi-stage support vector machine involved in the CBBO-SVM based transformer fault diagnosis method of the invention;
FIG. 3 is a flow chart of the transformer fault diagnosis method based on CBBO-SVM of the invention.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The invention relates to a CBBO-SVM-based transformer fault diagnosis method, the flow of which is shown in figure 1 and is implemented according to the following steps:
step 1, collected sample set S = { (x) with class label of oil-immersed transformer 1 ,x 2 ,x 3 ,...x n ),(y 1 ,y 2 ,y 3 ,...,x m ) Wherein, the training sample accounts for 50%, the testing sample accounts for 50%, x i Representing the properties of the sample, including 7 properties of hydrogen, methane, ethane, ethylene, acetylene, carbon monoxide and carbon dioxide, y i The category labels 1, 2, 3, 4, 5 and 6 are represented and respectively correspond to 6 states of normal state, medium and low temperature overheating, high energy discharge, low energy discharge and arc discharge;
step 2, preprocessing and normalizing the acquired data: normalizing both the input quantity and the output quantity to [0,1] by the following formula (1) so as to achieve the purpose of convenient calculation;
in the above formula, X max 、X min Respectively collected H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 Maximum and minimum values of the concentration of the gas, X i Represents H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 The gas concentration;
step 3, utilizing parameter c of CBBO to support vector machine i And σ i Optimizing;
step 3 is specifically implemented according to the following steps:
step 3.1, SVM parameter c i And σ i After the initialization operation is performed, c is i And σ i As each habitat P i Is (SIV) X i Then generating an initial population W i 1 ;
Step 3.2, calculating initial population W i 1 Population Habitat Suitability Index (HSI) f i 1 Obtaining f from the formula (2) i 1 ;
In the above formula, i is the number of the habitat variable; n is the total number of habitat variables; j represents each habitat variable P i A level of (d); c. C j As a habitat variable P i The total number of stages; beta is a ij Is a multiple regression coefficient;
ω ij is a weighting factor derived from the frequency distribution of habitat variables; a is a constant coefficient;
as shown in fig. 2, the cosine mobility model is specifically established as follows:
step 3.3, for the initial population W i 1 Carrying out migration operation, wherein in four common mobility models, namely a linear mobility model, a secondary mobility model, an index mobility model and a cosine mobility model, in order to make result parameters more optimal, a cosine mobility model which more accords with a natural rule is selected and established, and the migration rate lambda and the migration rate mu are calculated according to a formula (3);
in the above formula, I is the maximum mobility; e is the maximum migration rate; lambda is the mobility; mu is the migration rate; s max The maximum number of species that can be accommodated by the habitat; s is the number of species in the habitat;
after the migration operation, according to the migrationThe new species W is formed by emigration i 2 Then calculating the suitability index f of the population habitat i 2 ;
Step 3.4, performing variation operation on the population after the migration operation, wherein the variation rate m (x) of the migration habitat is along with the variation rate P of the species s Is increased, i.e.:
in the above formula: m (x) is habitat variation rate; m is max The maximum rate of variation; p s The probability of the species corresponding to the species number S in the habitat is shown; p max Is P s Maximum value of (d);
P s relationship to the mobility λ and the mobility μ:
randomly generating new species W in random habitats based on mutation operations i 3 Instead of the original species W i 2 Calculating the suitability index f of the population habitat i 3 ;
Step 3.5, adding a C chaos theory, updating the population by using a segmented Logistic chaos mapping equation (6) and generating a new initialized fitness vector f i 4 :
Wherein x is i Is a chaotic variable; k is the number of chaotic iterations; in step 3.5, when μ =4, the mapping equation (6) is in a chaotic state.
In order to ensure the diversity of species, chaotic search is carried out on the basis of Logistic chaotic mapping, and the optimal solution in the new fitness solution is recorded as J best In addition to J best Chaotic search basedNamely:
J n =J best +ωx i (7)
wherein x is i Is the solution of the chaotic mapping equation; omega is an adjustment coefficient;
for the obtained optimal solution J n And J best Comparing, taking the optimum of the two, repeating step 3.2, and updating the output parameter c according to the requirement i And σ i ;
Step 4, utilizing the parameter c obtained after the optimization in the step 3 i And σ i Establishing a multi-stage SVM model, and classifying by using sample set data so as to achieve the effect of fault diagnosis of the oil-immersed transformer;
as shown in fig. 3, the specific establishment method of the multi-stage support vector machine fault diagnosis model is as follows:
step 4 is specifically implemented according to the following steps:
step 4.1, firstly, kernel function selection is carried out;
selecting the most common Gaussian (RBF) kernel function from the linear kernel function, the polynomial kernel function, the sigmoid kernel function and the Gaussian (RBF) kernel function:
where ξ is a parameter controlling the height of a Gaussian (RBF) kernel;
step 4.2, taking training samples accounting for 50% of the total samples as model input for verification;
step 4.3, using c obtained in step 3.5 i And σ i Introducing and obtaining a new objective function f (omega, sigma) i ):
Where C is a penalty factor (meaning that the classification is strict when C tends to be large; meaning that greater error tolerance is possible when C tends to be small), a constant that needs to be specified,
parameters that need to be optimized; sigma i Is a relaxation factor.
Is also the invention
Step 4.4, solving by utilizing a Lagrange multiplier method, and combining an objective function and a constraint condition formula (10)
y i (ωx i +b)≥1-σ i (10)
solving omega and b to obtain a final decision function;
and 4.5, obtaining 6 decision functions through 6 layers of training, and establishing a 6-level support vector machine transformer fault diagnosis model.
Claims (5)
1. A transformer fault diagnosis method based on CBBO-SVM is characterized by comprising the following steps:
step 1, collected sample set S = { (x) with class label of oil-immersed transformer 1 ,x 2 ,x 3 ,...x n ),(y 1 ,y 2 ,y 3 ,...,x m ) Wherein, the training sample accounts for 50%, the testing sample accounts for 50%, x i Representing the properties of the sample, including 7 properties of hydrogen, methane, ethane, ethylene, acetylene, carbon monoxide and carbon dioxide, y i The category labels 1, 2, 3, 4, 5 and 6 are represented and respectively correspond to 6 states of normal state, medium and low temperature overheating, high energy discharge, low energy discharge and arc discharge;
step 2, preprocessing and normalizing the acquired data: normalizing both the input quantity and the output quantity to [0,1] by the following formula (1) so as to achieve the purpose of convenient calculation;
in the above formula, X max 、X min Respectively collected H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 Maximum and minimum values of the concentration of the gas, X i Represents H 2 、CH 4 、C 2 H 2 、C 2 H 6 、C 2 H 4 CO and CO 2 The gas concentration;
step 3, utilizing parameter c of CBBO to support vector machine i And σ i Optimizing;
step 4, utilizing the parameter c obtained after the optimization in the step 3 i And σ i And establishing a multi-stage SVM model, and classifying by using the sample set data so as to achieve the effect of fault diagnosis of the oil-immersed transformer.
2. The CBBO-SVM based transformer fault diagnosis method as claimed in claim 1, wherein the step 3 is implemented by the following steps:
step 3.1, SVM parameter c i And σ i After the initialization operation is performed, c is i And σ i As each habitat P i Is vector of fitness X i Then generating an initial population W i 1 ;
Step 3.2, calculating initial population W i 1 Population habitat suitability index f i 1 Obtaining f from the formula (2) i 1 ;
In the above formula, i is the number of the habitat variable; n is the total number of habitat variables; j represents each habitat variable P i A level of (d); c. C j As a habitat variable P i The total number of stages; beta is a ij Is a multiple regression coefficient; omega ij Is a weighting factor derived from the frequency distribution of habitat variables; a is a constant systemCounting;
step 3.3, for the initial population W i 1 Carrying out migration operation, and calculating the migration rate lambda and the migration rate mu according to the formula (3);
in the above formula, I is the maximum mobility; e is the maximum migration rate; lambda is the mobility; mu is the migration rate; s max The maximum number of species that can be accommodated by the habitat; s is the number of species in the habitat;
after the migration operation, a new species W is formed according to the migration in and out i 2 Then calculating the suitability index f of the population habitat i 2 ;
Step 3.4, carrying out variation operation on the population after the migration operation, wherein the variation rate m (x) of the migration habitat is along with the variation rate P of the species s Is increased, i.e.:
in the above formula: m (x) is habitat variability; m is a unit of max The maximum rate of variation; p s The probability of the species corresponding to the species with the number of S in the habitat is set; p max Is P s Maximum value of (d);
P s relationship to the mobility λ and the mobility μ:
randomly generating new species W in random habitats based on mutation operations i 3 Replacing the original species W i 2 Calculating the suitability index f of the population habitat i 3 ;
Step 3.5, adding a C chaos theory, and utilizing a segmented Logistic chaos mapping equation (6) to perform population matchingUpdating to generate new initialized fitness vector f i 4 :
Wherein x is i Is a chaotic variable; k is the number of chaotic iterations; in the step 3.5, when μ =4, the mapping equation (6) is in a chaotic state;
in order to ensure the diversity of species, chaotic search is carried out on the basis of Logistic chaotic mapping, and the optimal solution in the new fitness solution is recorded as J best In addition to J best Chaotic search is carried out on the basis of the chaotic search, namely:
J n =J best +ωx i (7)
wherein x is i Is the solution of the chaotic mapping equation; omega is an adjustment coefficient;
for the obtained optimal solution J n And J best Comparing, selecting the optimal one, repeating step 3.2, and updating output parameter c according to requirement i And σ i 。
3. The CBBO-SVM based transformer fault diagnosis method as claimed in claim 1, wherein the step 4 is implemented by the following steps:
step 4.1, firstly, kernel function selection is carried out;
step 4.2, taking training samples accounting for 50% of the total samples as model input for verification;
step 4.3, c obtained in step 3.5 i And σ i Introducing and obtaining a new objective function f (omega, sigma) i ):
Wherein c is a penalty factor; sigma i Is a relaxation factor;
step 4.4, solving by utilizing a Lagrange multiplier method, and combining an objective function and a constraint condition formula (10)
y i (ωx i +b)≥1-σ i (10)
solving omega and b to obtain a final decision function;
and 4.5, obtaining 6 decision functions through 6 layers of training, and establishing a 6-level support vector machine transformer fault diagnosis model.
4. The CBBO-SVM based transformer fault diagnosis method as claimed in claim 3, wherein the kernel function in step 4.4 is a Gaussian kernel function K (X, X) i )。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312408.7A CN111144431B (en) | 2018-11-06 | 2018-11-06 | Transformer fault diagnosis method based on CBBO-SVM |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312408.7A CN111144431B (en) | 2018-11-06 | 2018-11-06 | Transformer fault diagnosis method based on CBBO-SVM |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111144431A CN111144431A (en) | 2020-05-12 |
CN111144431B true CN111144431B (en) | 2023-03-24 |
Family
ID=70515224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811312408.7A Active CN111144431B (en) | 2018-11-06 | 2018-11-06 | Transformer fault diagnosis method based on CBBO-SVM |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111144431B (en) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701940B (en) * | 2009-10-26 | 2012-05-30 | 南京航空航天大学 | On-line transformer fault diagnosis method based on SVM and DGA |
CN103077402B (en) * | 2012-12-28 | 2016-05-11 | 国家电网公司 | Partial discharge of transformer mode identification method based on singular value decomposition algorithm |
CN107644127A (en) * | 2017-09-08 | 2018-01-30 | 西安工程大学 | Diagnosis Method of Transformer Faults based on IMABC Support Vector Machines Optimizeds |
-
2018
- 2018-11-06 CN CN201811312408.7A patent/CN111144431B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111144431A (en) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108596212B (en) | Transformer fault diagnosis method based on improved cuckoo search optimization neural network | |
Dai et al. | Dissolved gas analysis of insulating oil for power transformer fault diagnosis with deep belief network | |
Zhang et al. | A fault diagnosis model of power transformers based on dissolved gas analysis features selection and improved krill herd algorithm optimized support vector machine | |
CN105930901B (en) | A kind of Diagnosis Method of Transformer Faults based on RBPNN | |
CN110175386B (en) | Method for predicting temperature of electrical equipment of transformer substation | |
Neo et al. | Forecasting of photovoltaic power using deep belief network | |
US20210334658A1 (en) | Method for performing clustering on power system operation modes based on sparse autoencoder | |
CN115563563A (en) | Fault diagnosis method and device based on transformer oil chromatographic analysis | |
Fan et al. | Transformer fault diagnosis method based on improved whale optimization algorithm to optimize support vector machine | |
CN111104972A (en) | Method for identifying low-voltage risk of distribution room based on genetic algorithm optimization support vector machine multi-classifier | |
CN116842337A (en) | Transformer fault diagnosis method based on LightGBM (gallium nitride based) optimal characteristics and COA-CNN (chip on board) model | |
CN109919178A (en) | Failure prediction method based on characteristic quantity preferably with Wavelet Kernel Function LSSVM | |
CN116562114A (en) | Power transformer fault diagnosis method based on graph convolution neural network | |
CN111737907A (en) | Transformer fault diagnosis method and device based on deep learning and DGA | |
Huang et al. | Resilient power network structure for stable operation of energy systems: A transfer learning approach | |
CN112990593A (en) | Transformer fault diagnosis and state prediction method based on CSO-ANN-EL algorithm | |
CN111275204B (en) | Transformer state identification method based on hybrid sampling and ensemble learning | |
Wang et al. | Prediction method of wind farm power generation capacity based on feature clustering and correlation analysis | |
CN116960978A (en) | Offshore wind power prediction method based on wind speed-power combination decomposition reconstruction | |
CN111697560A (en) | Method and system for predicting load of power system based on LSTM | |
CN114881106A (en) | Transformer fault diagnosis method and device based on MPA-SVM | |
Nanfak et al. | Hybrid DGA method for power transformer faults diagnosis based on evolutionary k-means clustering and dissolved gas subsets analysis | |
CN111144431B (en) | Transformer fault diagnosis method based on CBBO-SVM | |
CN116883057A (en) | XGBoost-based high-precision power customer marketing channel preference prediction system | |
Tran et al. | A new grid search algorithm based on XGBoost model for load forecasting |
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 |