CN104155585A - GIS partial discharge type identification method based on GK fuzzy clustering - Google Patents
GIS partial discharge type identification method based on GK fuzzy clustering Download PDFInfo
- Publication number
- CN104155585A CN104155585A CN201410394763.9A CN201410394763A CN104155585A CN 104155585 A CN104155585 A CN 104155585A CN 201410394763 A CN201410394763 A CN 201410394763A CN 104155585 A CN104155585 A CN 104155585A
- Authority
- CN
- China
- Prior art keywords
- gis
- shelf depreciation
- fuzzy clustering
- partial discharge
- method based
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 27
- 230000007547 defect Effects 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000003909 pattern recognition Methods 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 11
- 239000012212 insulator Substances 0.000 claims description 10
- 238000005192 partition Methods 0.000 claims description 8
- 239000002245 particle Substances 0.000 claims description 7
- 239000000725 suspension Substances 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 6
- 238000002955 isolation Methods 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 4
- 241000208199 Buxus sempervirens Species 0.000 claims description 3
- 238000012706 support-vector machine Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 3
- 238000007635 classification algorithm Methods 0.000 abstract 2
- 238000013528 artificial neural network Methods 0.000 description 9
- 238000012360 testing method Methods 0.000 description 9
- 238000009413 insulation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 5
- 238000012549 training Methods 0.000 description 5
- 238000012423 maintenance Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 239000002923 metal particle Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 235000000981 Brassica parachinensis Nutrition 0.000 description 1
- 241001604868 Brassica rapa var. parachinensis Species 0.000 description 1
- 239000004698 Polyethylene Substances 0.000 description 1
- 229910052946 acanthite Inorganic materials 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011859 microparticle Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000012567 pattern recognition method Methods 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- XUARKZBEFFVFRG-UHFFFAOYSA-N silver sulfide Chemical compound [S-2].[Ag+].[Ag+] XUARKZBEFFVFRG-UHFFFAOYSA-N 0.000 description 1
- 229940056910 silver sulfide Drugs 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Landscapes
- Testing Relating To Insulation (AREA)
Abstract
The present invention discloses a GIS partial discharge type identification method based on GK fuzzy clustering. The method is characterized by comprising the steps of S01, constructing a GIS partial discharge gray scale image according to the acquired data; S02, extracting the fractal features of the GIS partial discharge gray scale image, namely a box dimension and an information dimension; S03, processing the fractal feature data further by a GK fuzzy clustering algorithm, isolating a GIS site interference signal; S04, designing a GIS partial discharge mode recognizer based on a least squares support vector machine classification algorithm; S05, identifying the GIS partial discharge type. By a GK fuzzy clustering method, the GIS site interference signal is isolated, thereby improving the accuracy of extracting the fractal features of the partial discharge signal. Meanwhile, by identifying the GIS partial discharge type by the least squares support vector machine classification algorithm, the correctness and rapidity of the discharge type identification are improved.
Description
Technical field
The present invention relates to a kind of GIS shelf depreciation kind identification method based on GK fuzzy clustering.
Background technology
Gas insulated combined electric appliance equipment (gas insulated switchgear, GIS, also claims SF6 fully closed combined electric unit) has the features such as floor area is little, operational reliability is high, easy to maintenance, in electric system, is widely applied.Yet, GIS because of manufacture or installation process in the fault that causes of defect under residual happen occasionally, its internal fault type is main mainly with insulativity fault, mainly comprises Four types: corona discharge, suspension electric discharge, the electric discharge of insulator surface defect and the electric discharge of free particle.When fault occurs, by defect distortion electric field, produce shelf depreciation.Therefore, by Partial Discharge Detection, judge that the potentiality fault of GIS is to grasping the insulation status of GIS and instructing its maintenance all to have great importance, and GIS local discharge signal is carried out to type identification, become the effective method of judgement GIS inherent vice type characteristic.
The feature extraction of GIS local discharge signal and pattern classification are two key issues in PD Pattern Recognition research.At present, the main method of extracting discharge signal characteristic parameter both at home and abroad comprises: statistical nature, wavelet character, waveform character, moment characteristics, assemblage characteristic and fractal characteristic.Wherein, fractal characteristic is due to its separating capacity and pattern description ability is superior, characteristic parameter is few and can reflect very complicated local discharge signal, and has obtained increasingly extensive application.But in current document, in the extraction of GIS local discharge signal fractal characteristic, do not consider the on-the-spot impact of disturbing fractal characteristic of GIS, if do not consider the on-the-spot impact of disturbing discharge signal, likely can, using undesired signal as discharge signal, cause thus fractal characteristic to extract the accurate problem of owing.
Aspect PD Pattern Recognition, document will have the sorter of the artificial neural network of good adaptive ability, robust performance and non-linear mapping capability as GIS shelf depreciation pattern at present, yet, artificial neural network, owing to having the determining of network structure, local minimum and owing study and cross the problems such as study, may bring certain impact to the correctness of recognition result.
Summary of the invention
For the problems referred to above, the invention provides a kind of GIS shelf depreciation kind identification method based on GK fuzzy clustering, adopt GK fuzzy clustering method, by the on-the-spot undesired signal isolation of GIS, greatly improve the accuracy that local discharge signal fractal characteristic extracts, meanwhile, adopt least square method supporting vector machine sorting algorithm (Least Squares Support Vector Machine, LS-SVM) identification GIS shelf depreciation type, improves correctness and the rapidity of electric discharge type identification.
For realizing above-mentioned technical purpose, reach above-mentioned technique effect, the present invention is achieved through the following technical solutions:
GIS shelf depreciation kind identification method based on GK fuzzy clustering, is characterized in that, comprises the steps:
S01: according to the data construct GIS shelf depreciation gray level image gathering;
S02: extract the fractal characteristic of GIS shelf depreciation gray level image, i.e. box counting dimension and information dimension;
S03: adopt GK fuzzy clustering algorithm further to process fractal characteristic data, the on-the-spot undesired signal of isolation GIS;
S04: based on least square method supporting vector machine design of algorithm GIS PD Pattern Recognition device;
S05: identification GIS shelf depreciation type.
The shelf depreciation sample data of the method based under corona discharge, suspension electric discharge, four kinds of different electric pressures of defect of the electric discharge of insulator surface defect and free particle electric discharge, extract the fractal property of shelf depreciation gray level image as recognition feature amount, consider the on-the-spot impact of disturbing local discharge signal of GIS simultaneously, utilize GK fuzzy clustering algorithm further to process fractal characteristic amount, to extract the fractal characteristic amount after isolation is disturbed.Finally, based on least square method supporting vector machine design of algorithm GIS PD Pattern Recognition device, can effectively identify GIS shelf depreciation type, than Artificial Neural Network, have advantages of that discrimination is high, good stability simultaneously.
The invention has the beneficial effects as follows:
1) the method takes into full account the impact of the on-the-spot undesired signal of GIS on local discharge signal fractal characteristic, adopts GK fuzzy clustering algorithm that undesired signal is isolated, and can greatly improve the accuracy that GIS local discharge signal fractal characteristic extracts.
2) the method can be avoided the determining of network structure in Artificial Neural Network, local minimum and owe study and cross the problems such as study, can effectively improve the correctness that electric discharge type is identified.
3) the GIS shelf depreciation type that the method identifies can be used as the direct basis of judgement GIS inherent vice type, and to grasping the insulation status of GIS and instructing its maintenance all to have great importance.
Accompanying drawing explanation
Fig. 1 is the process flow diagram that the present invention is based on the GIS shelf depreciation kind identification method of GK fuzzy clustering;
Fig. 2 is the fractal characteristic schematic diagram of the untreated four kinds of front GIS shelf depreciation defects of the present invention;
Fig. 3 is the fractal characteristic schematic diagram of four kinds of GIS shelf depreciation defects after GK fuzzy clustering of the present invention is processed;
Fig. 4 is the schematic diagram that GIS of the present invention inside manually arranges four kinds of defects;
Fig. 5 is the schematic diagram of GIS testing ground undesired signal in the embodiment of the present invention;
Fig. 6 is recognition result (σ=9.8463 before LS-SVM parameter optimization of the present invention; γ=94730.521);
Fig. 7 is recognition result (σ=13.5741528 after LS-SVM parameter optimization of the present invention; γ=23478.3341).
Embodiment
Below in conjunction with accompanying drawing and specific embodiment, technical solution of the present invention is described in further detail, so that those skilled in the art can better understand the present invention also, can be implemented, but illustrated embodiment is not as a limitation of the invention.
The invention discloses a kind of gas insulated combined electrical equipment (Gas Insulated Switchgear based on Gustafson-Kessel (GK) fuzzy clustering, GIS) shelf depreciation kind identification method, belongs to fault detect and the state estimation field of high voltage electric equipment.
Fig. 1 is the process flow diagram of the GIS shelf depreciation kind identification method based on GK fuzzy clustering, comprises the steps:
S01: according to the data construct GIS shelf depreciation gray level image gathering;
S02: extract the fractal characteristic of GIS shelf depreciation gray level image, i.e. box counting dimension and information dimension;
S03: adopt GK fuzzy clustering algorithm further to process fractal characteristic data, the on-the-spot undesired signal of isolation GIS;
S04: based on least square method supporting vector machine design of algorithm GIS PD Pattern Recognition device;
S05: identification GIS shelf depreciation type.
Be specially: in GIS local discharge signal feature extraction, first utilize the data construct GIS shelf depreciation gray level image gathering, then extract the fractal characteristic of GIS shelf depreciation gray level image, that is: the information dimension of gray-scale map negative half period
the information dimension of positive half cycle
the box counting dimension of negative half period
and the box counting dimension of positive half cycle
and using this characteristic parameter as GIS shelf depreciation type identification.Wherein, concrete methods of realizing can be referring to [Sun Caixin, permitted peak, Tang Ju, Deng. take the GIS Partial Discharge Pattern Recognition Method [J] that box counting dimension and information dimension be recognition feature amount. Proceedings of the CSEE, 2005, 25 (3): 100-104.] institute's extracting method, the result of the gray level image fractal property that the GIS Partial Discharge Data based on gathering builds as shown in Figure 2, wherein " ☆ " represents pollution severity of insulators electric discharge (electric discharge of insulator surface defect), "●" represents metal tip electric discharge (corona discharge), "+" represents suspension electrode electric discharge (suspension electric discharge), " * " represents metal particle electric discharge (free particle electric discharge).
As can be seen from Figure 2, the fractal property of 4 kinds of discharge defects (box counting dimension, information dimension) is all gathered in a certain region, and fractal characteristic has stronger separating capacity.Yet testing ground exists larger interference, inevitably can have undesired signal in the local discharge signal gathering, even stronger undesired signal can covering internal discharge signal.But in local discharge signal fractal characteristic extracts, prior art is not all considered the impact of undesired signal on normal discharge signal, also be in extracted fractal characteristic data as shown in Figure 2, to comprise most probably the fractal characteristic data of undesired signal, if do not processed disturbing a minute graphic data, and the directly input using the fractal characteristic extracting as electric discharge type mode discriminator, this will bring certain impact to the precision of recognition result.Therefore,, for eliminating the impact of disturbing, this method, on the basis of extracted fractal characteristic data, adopts GK fuzzy clustering algorithm further to process fractal characteristic data, to obtain pure electric discharge type recognition feature amount.
GK fuzzy clustering algorithm is to utilize the self-adaptation distance of cluster covariance matrix to measure, the center vector V=(v of its fuzzy clustering
1, v
2..., v
c)
tdegree of membership matrix (or fuzzy partition matrix) U=[μ with data set
ij]
c * Nby minimizing objective function, try to achieve the number that wherein N is sample, the number that c is cluster, μ
ijfor the degree of membership of data point with respect to cluster centre, and meet:
Suppose that data sequence to be identified is that the data sequence that step S02 obtains is X=(x
1, x
2..., x
n), it minimize objective function and can be expressed as:
Wherein, and m ∈ [1, ∞), be a weighted index that characterizes classification fog-level;
for data point x in data sequence
jto cluster centre v
idistance, and be a square of inner product apart from norm, be expressed as:
Wherein, M
ifor following matrix:
Wherein, F
iit is the covariance matrix of i cluster centre.
Utilize Lagrangian multiplication, to minimizing target function type (2), be optimized, make objective function obtain minimizing two necessary conditions and be respectively:
The concrete execution step of GK fuzzy clustering algorithm is as follows:
(I) data-oriented sequence X=(x
1, x
2..., x
n), and select clusters number 1<c<N, Weighting exponent m>=1 and termination allowable error ε >0.
(II) arranges fuzzy partition matrix U=[μ
ij]
c * Ninitial value, and make it meet constraint equation (1), set iteration number of times l=0,1 ... (iterations l is a variable, and algorithm repeats how many times), upgrades cluster centre v according to formula (6)
i.
(III) calculates the covariance matrix F of i cluster centre
i:
(IV) utilizes the covariance matrix F of step (III)
i, according to formula (4) compute matrix M
i, then according to formula (3), ask for data point x in data sequence
jto cluster centre v
idistance norm
(V) is according to the distance norm of step (IV) gained
utilize formula (5) to upgrade fuzzy partition matrix U.
For any one, stop allowable error positive number ε > 0, if meet fuzzy partition matrix || U
(l+1)-U
(l)|| < ε, computing stops, otherwise with the fuzzy partition matrix U iteration number of times l=0 after upgrading, 1 ..., according to formula (6), upgrade cluster centre v
i, and perform step successively (III) to (V), until satisfy condition.Namely when not satisfying condition, whole algorithm continues to carry out, i.e. iteration, until satisfy condition.
The fractal characteristic data that extract are further processed, in the situation that taking into full account on-the-spot undesired signal, extract the fractal characteristic data of undesired signal, and utilize above-mentioned GK fuzzy clustering algorithm to carry out cluster to normal shelf depreciation fractal characteristic data, reject the fractal characteristic data of mixing undesired signal therein.
The parameter of GK fuzzy clustering algorithm is selected as follows: c=23, m=1, ε=0.001.Fractal characteristic in Fig. 2 is carried out to fractal characteristic after GK fuzzy clustering processing as shown in Figure 3, and wherein circle represents the fractal characteristic of undesired signal after GK fuzzy clustering is processed.
In step S04, based on least square method supporting vector machine design of algorithm GIS PD Pattern Recognition device, specific implementation process is:
Suppose to have λ sample
training set, z wherein
i∈ R
nbe i input data; y
i{ 1 ,+1} is i output data to ∈.Constructing a sorter with following form is category support vector machines simulated target:
So that sample z can correctly be classified by f (z).
The Function Estimation problem of least square method supporting vector machine sorting algorithm can be described as the following problem that solves:
Meet equality constraint
Introduce Lagrangian function:
α wherein
i, (i=1,2 ..., λ) be Lagrange multiplier,
represent nuclear space mapping function, e
i∈ R represents error variance, w ∈ R
nhfor weight vector, b is departure, and γ is error penalty factor, and J (w, e) is the optimal function shown in formula (9).
The necessary condition existing according to extreme value, has following equation to form vertical:
Formula (12) is arranged and can be obtained:
Wherein
e=[e
1, e
2..., e
λ]
t, y=[y
1, y
2..., y
λ]
t,
α=[α
1, α
2..., α
λ]
t.
According to Mercer condition:
Thus, the system of equations obtaining is relevant with α and b.
Therefore, system of equations (13) can be converted into:
By formula (15), can be obtained:
The sorter that can obtain supporting vector machine model according to the first formula of formula (8), formula (12), formula (14) and formula (16) is:
ψ (x, x wherein
i) be any symmetric function and meet Mercer condition.Based on RBF type function, have advantages of, the kernel function of selecting is herein following radial basis function:
ψ(z
i,z
j)=exp(-||z
i-z
j||
2/(2σ
2)) (18)
In GIS partial discharge test, the artificial defect arranging has 4 kinds, therefore need the shelf depreciation mode type of identification also to have accordingly 4 classes, therefore when solving PD Pattern Recognition, LS-SVM recognizer need to be expanded to multicategory classification problem, yet multicategory classification problem is compared with two class classification problems, need to consider a plurality of samples, decision boundary is complicated, be difficult to training and precision not high.For the identification problem of 4 kinds of electric discharge types, this method adopts the mode of a plurality of 2 classification LS-SVM recognizer combinations to realize when structure LS-SVM pattern recognition classifier device, and this mode has simply, is easy to the advantages such as realization.
Based on above-mentioned analysis, for better identifying 4 kinds of GIS shelf depreciation types, need to design 62 classification LS-SVM recognizers, and it is defined as respectively to SVM
11, SVM
12, SVM
13, SVM
21, SVM
22, SVM
23so, by 6 LS-SVM recognizers, take the output of ballot method can judge the electric discharge type of GIS shelf depreciation sample, its voting rule is as shown in table 1.
Table 1 ballot method discharge mode recognition rule
In addition, LS-SVM recognizer exists kernel function ψ (z
i, z
j) the selection problem of width cs and error penalty factor γ.It is whether suitable that these parameters are selected, and will produce material impact to accuracy of identification.For obtaining optimum parameter, this method (specifically can be referring to document SUYKENS J.A.K by the method for cross validation (CV), VANDEWALLE J.Least squares support vector machine classifiers[J] .Neural Processing letters, 1999,9 (3): 293-30.) to LS-SVM kernel function ψ (z
i, z
j) width cs and error penalty factor γ choose.
For the GIS shelf depreciation kind identification method based on GK fuzzy clustering described in clear explanation this method, the GIS intensive care monitoring system of Jiangsu Province's DianKeYuan of take is experimental study platform, and using the local discharge signal gathering on its platform as fundamental research data.
In the collection of each test figure sample, corona discharge (bus spine and barrel spine), suspension electric discharge, the electric discharge of insulator surface defect and 4 kinds of artificial defects of free particle electric discharge are set in GIS inside respectively, the setting up procedure of each defect is as follows:
1) corona discharge: a long draw point that is 1mm for 53mm, its most advanced and sophisticated place equivalent redius is installed on GIS cavity bus surface and is carried out the most advanced and sophisticated defect of simulated high-pressure conductor metal, equally also same draw point can be fixed on GIS cavity inner wall, as shown in Fig. 4 (a).
2) floating potential defect: by polyethylene insulation screw rod, life is belonged to aluminium flake and be fixed on high-pressure conductor in GIS cavity, both are at a distance of 1.5mm, insulation blanket stud can be for regulating the clearance distance of suspension aluminium flake and guide rod, its length is enough to avoid edge flashing, with this, simulate floating potential defect, as shown in Fig. 4 (b).
3) pollution severity of insulators: coat at insulator inside surface-layer has the silver sulfide powder simulation surface filth of semiconduction character, produces along face electric field and concentrates, and then simulation creeping discharge, as shown in Fig. 4 (c).
4) free particulate electric discharge: due in GIS, electrically conductive microparticle is maximum to the harm of insulation, therefore the harm that it has been generally acknowledged that non-conductive particulate is much smaller, and is that the circular spheroid of stainless steel of 1mm is simulated free particle defects, by sprinkle diameter on cavity as shown in Fig. 4 (d).
Above-mentioned 4 kinds of insulation defects carry out respectively the partial discharge test under 47kV, 50kV, 55kV and 60kV voltage, during test, ultra-high frequency signal adopts the external mode of sensor, by the high-speed oscilloscope of 4GHz/s, sample respectively and store the Partial Discharge in continuous 100 power frequency periods under every kind of defect situation, under each trial voltage, gathering 30 groups of data.
In GIS partial discharge test, on-the-spot electromagnetic interference (EMI) and background interference have larger impact to local discharge signal.Obtaining Accurate local discharge signal from strong interference how is to realize that GIS real-time online detects and the key point of electric discharge type pattern-recognition.For the impact of Analysis interference on shelf depreciation, when carrying out above-mentioned 4 kinds of tests, also on-the-spot undesired signal is gathered, it disturbs as shown in Figure 5, from on-the-spot undesired signal, can find out, there is stronger interference source in testing ground, moreover power frequency correlativity is strong.Therefore, be necessary to take appropriate measures, to suppress the impact of undesired signal on normal discharge signal and electric discharge type pattern-recognition.
The present embodiment selects 30 groups of data that 4 kinds of defect local discharge signals processed through GK fuzzy clustering algorithms as training sample, and every group of sample data comprises: the information dimension of negative half period
, positive half cycle information dimension
the box counting dimension of negative half period
and the box counting dimension of positive half cycle
For checking LS-SVM kernel function ψ (x
i, x
j) width cs and error penalty factor γ choose the impact on accuracy of identification, take pollution severity of insulators defect as example, by the recognition result before parameter optimization with utilize cross validation method to compare to the identification of the electric discharge type after parameter optimization, its comparing result as shown in Figure 6,7.
By comparing result Fig. 6 and Fig. 7, can find out that the accuracy of identification after parameter optimization obviously improves, and illustrates kernel function ψ (x
i, x
j) selection of width cs and error penalty factor γ has a direct impact accuracy of identification.
For further checking the performance of institute's extracting method herein, build respectively the partial discharges fault recognizer based on BP neural network and RBF neural network, and with contrasting based on GK fuzzy clustering recognition device of carrying herein, by square error MSE and accuracy, weigh the quality of recognition performance.
Square error MSE computing formula is as follows:
Accuracy is described as:
Wherein, M is the right number of training data, y
kfor actual value,
identification Output rusults for LS-SVM.
The PD Pattern Recognition result of 4 kinds of defects based on BP neural network, RBF neural network and GK fuzzy clustering recognition method is as shown in table 2-5.
Corona discharge recognition result under three kinds of algorithms of table 2
Corona discharge | The inventive method | BP | RBF |
MSE | 0.001 | 2.283 | 1.896 |
Accuracy | 95.6% | 81.6% | 87.13% |
Time | 0.911 | 63.36 | 36.84 |
Floating potential discharge recognition result under three kinds of algorithms of table 3
Corona discharge | The inventive method | BP | RBF |
MSE | 0.00135 | 2.239 | 1.811 |
Accuracy | 91.1% | 85.6% | 88.4% |
Time | 0.911 | 63.36 | 36.86 |
Pollution severity of insulators electric discharge recognition result under three kinds of algorithms of table 4
Corona discharge | The inventive method | BP | RBF |
MSE | 0.00087 | 1.233 | 0.935 |
Accuracy | 97.9% | 94.6% | 95.4% |
Time | 0.911 | 63.36 | 36.86 |
Free particle electric discharge recognition result contrast under three kinds of algorithms of table 5
Corona discharge | The inventive method | BP | RBF |
MSE | 0.0013 | 2.201 | 1.813 |
Accuracy | 92.6% | 84.5% | 87.4% |
Time | 0.911 | 63.36 | 36.86 |
From the recognition result contrast of table 2-table 5, can find out, adopt the GIS partial discharges fault recognizer of institute of the present invention extracting method all to there is higher discrimination to 4 kinds of discharge faults, especially corona discharge and the electric discharge of free metal particle are had than BP network and the better recognition effect of RBF network, simultaneously the training time of algorithm the shortest (in table, the unit of Time is second), real-time performance is good.
The invention has the beneficial effects as follows:
1) the method takes into full account the impact of the on-the-spot undesired signal of GIS on local discharge signal fractal characteristic, adopts GK fuzzy clustering algorithm that undesired signal is isolated, and can greatly improve the accuracy that GIS local discharge signal fractal characteristic extracts.
2) the method can be avoided the determining of network structure in Artificial Neural Network, local minimum and owe study and cross the problems such as study, can effectively improve the correctness that electric discharge type is identified.
3) the GIS shelf depreciation type that the method identifies can be used as the direct basis of judgement GIS inherent vice type, and to grasping the insulation status of GIS and instructing its maintenance all to have great importance.
These are only the preferred embodiments of the present invention; not thereby limit the scope of the claims of the present invention; every equivalent structure or conversion of equivalent flow process that utilizes instructions of the present invention and accompanying drawing content to do; or be directly or indirectly used in other relevant technical fields, be all in like manner included in scope of patent protection of the present invention.
Claims (6)
1. the GIS shelf depreciation kind identification method based on GK fuzzy clustering, is characterized in that, comprises the steps:
S01: according to the data construct GIS shelf depreciation gray level image gathering;
S02: extract the fractal characteristic of GIS shelf depreciation gray level image, i.e. box counting dimension and information dimension;
S03: adopt GK fuzzy clustering algorithm further to process fractal characteristic data, the on-the-spot undesired signal of isolation GIS;
S04: based on least square method supporting vector machine design of algorithm GIS PD Pattern Recognition device;
S05: identification GIS shelf depreciation type.
2. the GIS shelf depreciation kind identification method based on GK fuzzy clustering according to claim 1, is characterized in that, step S03 specifically comprises the steps:
A1: suppose the data sequence X=(x that step S02 obtains
1, x
2..., x
n), select clusters number 1<c<N, Weighting exponent m>=1 and termination allowable error ε >0;
A2: it is degree of membership matrix U=[μ that fuzzy partition matrix is set
ij]
c * Ninitial value, and make it meet constraint equation (1):
A3: set iteration number of times l=0,1 ..., according to formula (6):
Upgrade cluster centre v
i;
A4: the covariance matrix F that calculates i cluster centre
i:
A5: the covariance matrix F that utilizes step a4
i, according to formula (4)
Compute matrix M
i, then according to formula (3)
Ask for data point x in data sequence
jto cluster centre v
idistance norm
A6: according to the distance norm of step a5 gained
utilize formula (5)
Upgrade fuzzy partition matrix U;
A7: if meet fuzzy partition matrix || U
(l+1)-U
(l)|| < ε, computing stops, otherwise increases the iterations of step a3, calculates successively, until satisfy condition according to step a3-a7;
In above formula: the number that N is sample, the number that c is cluster, μ
ijfor the degree of membership of data point with respect to cluster centre,
for data point x in data sequence
jto cluster centre v
idistance, F
ibe the covariance matrix of i cluster centre, and m ∈ [1, ∞), be a weighted index that characterizes classification fog-level.
3. the GIS shelf depreciation kind identification method based on GK fuzzy clustering according to claim 1, is characterized in that, in step S04, adopts several 2 classification LS-SVM recognizer, and adopts the type of ballot method judgement GIS local discharge signal.
4. the GIS shelf depreciation kind identification method based on GK fuzzy clustering according to claim 3, is characterized in that, adopts 62 classification LS-SVM recognizers.
5. the GIS shelf depreciation kind identification method based on GK fuzzy clustering according to claim 4, is characterized in that, the kernel function width cs of LS-SVM and error penalty factor parameter γ adopt cross validation method to determine.
6. according to the GIS shelf depreciation kind identification method based on GK fuzzy clustering described in claim 1-5 any one, it is characterized in that, GIS shelf depreciation type comprises corona discharge, suspension electric discharge, the electric discharge of insulator surface defect and the electric discharge of free particle Four types altogether.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410394763.9A CN104155585A (en) | 2014-08-12 | 2014-08-12 | GIS partial discharge type identification method based on GK fuzzy clustering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410394763.9A CN104155585A (en) | 2014-08-12 | 2014-08-12 | GIS partial discharge type identification method based on GK fuzzy clustering |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104155585A true CN104155585A (en) | 2014-11-19 |
Family
ID=51881137
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410394763.9A Pending CN104155585A (en) | 2014-08-12 | 2014-08-12 | GIS partial discharge type identification method based on GK fuzzy clustering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104155585A (en) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388402A (en) * | 2015-10-30 | 2016-03-09 | 山东电力研究院 | Cross-linked cable partial discharge mode recognition method based on parameter optimization SVM (Support Vector Machine) algorithm |
CN105785236A (en) * | 2016-03-02 | 2016-07-20 | 国网江西省电力科学研究院 | GIS local discharge detection external interference signal elimination method |
CN106646165A (en) * | 2016-12-23 | 2017-05-10 | 西安交通大学 | Method and system for classifying and positioning internal insulation defects of GIS |
CN106778855A (en) * | 2016-12-08 | 2017-05-31 | 华中科技大学 | A kind of method that utilization visible images colourity treatment divides the corona discharge stage |
CN107192686A (en) * | 2017-04-11 | 2017-09-22 | 江苏大学 | A kind of Possibility Fuzzy Clustering local tea variety discrimination method of fuzzy covariance matrix |
CN107329053A (en) * | 2017-05-26 | 2017-11-07 | 华南理工大学 | A kind of method that utilization high-speed camera measures suspended conductor induction discharge frequency |
CN109670536A (en) * | 2018-11-30 | 2019-04-23 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | A kind of local discharge signal clustering method in the case of multi-source electric discharge and interference superposition |
CN110470964A (en) * | 2019-08-13 | 2019-11-19 | 国网天津市电力公司电力科学研究院 | GIS point discharge stage judgment method and judgment means based on maintenance decision purpose |
JP2020046202A (en) * | 2018-09-14 | 2020-03-26 | 株式会社東芝 | Partial discharge detection device, partial discharge detection method, partial discharge detection system, and computer program product |
CN111537850A (en) * | 2020-05-21 | 2020-08-14 | 北京传动联普科技有限公司 | Big data identification management system for intelligent separation and classification diagnosis of partial discharge signals |
CN111610418A (en) * | 2020-05-28 | 2020-09-01 | 华乘电气科技股份有限公司 | GIS partial discharge positioning method based on intelligent ultrahigh frequency sensor |
CN111610417A (en) * | 2020-05-28 | 2020-09-01 | 华乘电气科技股份有限公司 | Discharge signal source separation method based on community discovery |
CN112147465A (en) * | 2020-08-13 | 2020-12-29 | 国网浙江省电力有限公司电力科学研究院 | GIS partial discharge optical diagnosis method based on multi-fractal and extreme learning machine |
CN112305381A (en) * | 2020-09-21 | 2021-02-02 | 国网山东省电力公司临沂供电公司 | Method and system for monitoring and positioning online partial discharge of distribution cable |
CN112816838A (en) * | 2021-01-06 | 2021-05-18 | 国网重庆市电力公司电力科学研究院 | GIS equipment defect diagnosis device and method based on FFT, VMD and LS-SVM |
WO2021213142A1 (en) * | 2020-04-23 | 2021-10-28 | 中车株洲电力机车研究所有限公司 | Traction motor fault diagnosis method and apparatus |
CN114037021A (en) * | 2021-12-02 | 2022-02-11 | 国网上海市电力公司 | Construction method of alternating current/direct current partial discharge defect type multi-classifier |
CN115166453A (en) * | 2022-09-08 | 2022-10-11 | 国网智能电网研究院有限公司 | Partial discharge continuous monitoring method and device based on edge real-time radio frequency pulse classification |
-
2014
- 2014-08-12 CN CN201410394763.9A patent/CN104155585A/en active Pending
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388402A (en) * | 2015-10-30 | 2016-03-09 | 山东电力研究院 | Cross-linked cable partial discharge mode recognition method based on parameter optimization SVM (Support Vector Machine) algorithm |
CN105785236B (en) * | 2016-03-02 | 2019-02-26 | 国网江西省电力科学研究院 | A kind of GIS partial discharge detection external interference signals method for removing |
CN105785236A (en) * | 2016-03-02 | 2016-07-20 | 国网江西省电力科学研究院 | GIS local discharge detection external interference signal elimination method |
CN106778855A (en) * | 2016-12-08 | 2017-05-31 | 华中科技大学 | A kind of method that utilization visible images colourity treatment divides the corona discharge stage |
CN106778855B (en) * | 2016-12-08 | 2019-05-10 | 华中科技大学 | A method of it is handled using visible images coloration and divides the corona discharge stage |
CN106646165A (en) * | 2016-12-23 | 2017-05-10 | 西安交通大学 | Method and system for classifying and positioning internal insulation defects of GIS |
CN107192686A (en) * | 2017-04-11 | 2017-09-22 | 江苏大学 | A kind of Possibility Fuzzy Clustering local tea variety discrimination method of fuzzy covariance matrix |
CN107192686B (en) * | 2017-04-11 | 2020-08-28 | 江苏大学 | Method for identifying possible fuzzy clustering tea varieties by fuzzy covariance matrix |
CN107329053A (en) * | 2017-05-26 | 2017-11-07 | 华南理工大学 | A kind of method that utilization high-speed camera measures suspended conductor induction discharge frequency |
CN107329053B (en) * | 2017-05-26 | 2019-11-15 | 华南理工大学 | A method of suspended conductor induction discharge frequency is measured using high-speed camera |
JP2020046202A (en) * | 2018-09-14 | 2020-03-26 | 株式会社東芝 | Partial discharge detection device, partial discharge detection method, partial discharge detection system, and computer program product |
JP7242225B2 (en) | 2018-09-14 | 2023-03-20 | 株式会社東芝 | Partial discharge detection device, partial discharge detection method, partial discharge detection system, and computer program |
CN109670536A (en) * | 2018-11-30 | 2019-04-23 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | A kind of local discharge signal clustering method in the case of multi-source electric discharge and interference superposition |
CN109670536B (en) * | 2018-11-30 | 2023-05-30 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Partial discharge signal clustering method under multi-source discharge and interference superposition condition |
CN110470964A (en) * | 2019-08-13 | 2019-11-19 | 国网天津市电力公司电力科学研究院 | GIS point discharge stage judgment method and judgment means based on maintenance decision purpose |
WO2021213142A1 (en) * | 2020-04-23 | 2021-10-28 | 中车株洲电力机车研究所有限公司 | Traction motor fault diagnosis method and apparatus |
CN111537850A (en) * | 2020-05-21 | 2020-08-14 | 北京传动联普科技有限公司 | Big data identification management system for intelligent separation and classification diagnosis of partial discharge signals |
CN111537850B (en) * | 2020-05-21 | 2022-05-27 | 北京传动联普科技有限公司 | Big data identification management system for intelligent separation and classification diagnosis of partial discharge signals |
CN111610418A (en) * | 2020-05-28 | 2020-09-01 | 华乘电气科技股份有限公司 | GIS partial discharge positioning method based on intelligent ultrahigh frequency sensor |
CN111610417B (en) * | 2020-05-28 | 2022-03-15 | 华乘电气科技股份有限公司 | Discharge signal source separation method based on community discovery |
CN111610418B (en) * | 2020-05-28 | 2022-11-22 | 华乘电气科技股份有限公司 | GIS partial discharge positioning method based on intelligent ultrahigh frequency sensor |
CN111610417A (en) * | 2020-05-28 | 2020-09-01 | 华乘电气科技股份有限公司 | Discharge signal source separation method based on community discovery |
CN112147465B (en) * | 2020-08-13 | 2022-06-03 | 国网浙江省电力有限公司电力科学研究院 | GIS partial discharge optical diagnosis method based on multi-fractal and extreme learning machine |
CN112147465A (en) * | 2020-08-13 | 2020-12-29 | 国网浙江省电力有限公司电力科学研究院 | GIS partial discharge optical diagnosis method based on multi-fractal and extreme learning machine |
CN112305381A (en) * | 2020-09-21 | 2021-02-02 | 国网山东省电力公司临沂供电公司 | Method and system for monitoring and positioning online partial discharge of distribution cable |
CN112816838A (en) * | 2021-01-06 | 2021-05-18 | 国网重庆市电力公司电力科学研究院 | GIS equipment defect diagnosis device and method based on FFT, VMD and LS-SVM |
CN114037021A (en) * | 2021-12-02 | 2022-02-11 | 国网上海市电力公司 | Construction method of alternating current/direct current partial discharge defect type multi-classifier |
CN115166453A (en) * | 2022-09-08 | 2022-10-11 | 国网智能电网研究院有限公司 | Partial discharge continuous monitoring method and device based on edge real-time radio frequency pulse classification |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104155585A (en) | GIS partial discharge type identification method based on GK fuzzy clustering | |
Peng et al. | Random forest based optimal feature selection for partial discharge pattern recognition in HV cables | |
CN103076547B (en) | Method for identifying GIS (Gas Insulated Switchgear) local discharge fault type mode based on support vector machines | |
CN102809718B (en) | Ultra-high-frequency partial discharge signal identification method for gas insulated switchgear (GIS) | |
CN108957251B (en) | Cable joint partial discharge mode identification method | |
Wu et al. | Defect recognition and condition assessment of epoxy insulators in gas insulated switchgear based on multi-information fusion | |
CN109633368A (en) | The method of duration power quality disturbances containing distributed power distribution network based on VMD and DFA | |
Dai et al. | Feature extraction of GIS partial discharge signal based on S‐transform and singular value decomposition | |
WO2019015311A1 (en) | Vibration signal support vector machine-based gil fault online monitoring system | |
Chan et al. | Time-frequency sparsity map on automatic partial discharge sources separation for power transformer condition assessment | |
CN102645621B (en) | Multisource partial discharge detection method and device of transformer substation based on space characteristic parameters | |
CN102645620A (en) | Multisource partial discharge detection method and device of transformer substation based on time-frequency characteristic parameters | |
CN103558519A (en) | GIS partial discharge ultrasonic signal identification method | |
CN105785236B (en) | A kind of GIS partial discharge detection external interference signals method for removing | |
CN108663202A (en) | GIS mechanical failure diagnostic methods based on chaos cuckoo algorithm and system | |
Zhao et al. | Interpretation of transformer winding deformation fault by the spectral clustering of FRA signature | |
CN115600088A (en) | Distribution transformer fault diagnosis method based on vibration signals | |
CN112180221B (en) | GIS unknown category partial discharge identification method based on double-measurement supervision rule | |
CN108919067A (en) | A kind of recognition methods for GIS partial discharge mode | |
CN110738232A (en) | grid voltage out-of-limit cause diagnosis method based on data mining technology | |
CN108304567A (en) | High-tension transformer regime mode identifies and data classification method and system | |
CN110703078A (en) | GIS fault diagnosis method based on spectral energy analysis and self-organizing competition algorithm | |
CN109376626A (en) | A kind of GIS switch defect diagnostic method based on radiated electric field characteristic parameter support vector machines | |
CN105403777A (en) | Classification and discrimination method of aging state of composite insulator | |
Zhu et al. | Faulty line identification method based on bayesian optimization for distribution network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141119 |
|
RJ01 | Rejection of invention patent application after publication |