Combined Diagnosis of PD Based on the Multidimensional Parameters

This paper presents a comprehensive multiparameter diagnosis method based on multiple partial discharge (PD) signals which include high-frequency current (HFC), ultrasound, and ultrahigh frequency (UHF). The HFC, ultrasound, and UHF PD are calculated under different types of faults. Therefor the characteristic values, as nine basic characteristic parameters, eight phase characteristic parameters, and the like are calculated. Diagnose signals are found with the method based on information fusion and semisupervised learning for HFC PD, adaptive mutation parameters of particle entropy for ultrasonic signals, and IIA-ART2A neural network for UHF signals. In addition, integrate the diagnostic results, which are the probability of fault of various defects and matrix, of different PD diagnosis signals, and analysis with Sugeno fuzzy integral to get the final diagnosis.


Introduction
Only part of the power equipment insulation generates a discharge.The effect that a fixed discharge channel is not formed in the discharge area is called PD.The PD can reflect the presence of insulation faults in the power equipment.The development of PD will speed up the insulation deterioration of the internal power equipment, leading to internal electrical equipment insulation fault occurrence, resulting in equipment fault and reducing the reliability of power supply.So it is vital to monitor the existence of PD in electrical equipment, as well as the type of discharge and early detection of equipment insulation faults, to deal with it in time, and to avoid the failure of electrical equipment and reduce the serious economic losses caused by power equipment insulation fault.
Kim et al. [1] examined an intelligent spacer built into the internal type UHF PD sensors.Three-dimensional electromagnetic simulations were performed to analyze electricfield distribution of the single-phase GIS and three-phase GIS.After considering the spacer's specifications, sensor structures were designed and analyzed using the 3D EM simulator.
Kaneko et al. [2] used a method for diagnosing the insulation of a gas-insulated switch gear which is to detect PDs using signals in the UHF band.
Schwarz et al. [3] presented that the aim was the detection of beginning destruction in the electrical insulation as a result of electrical stress.A multiplicity of different PD sources and their appearances show different physical and electrical characteristics [4].
Sahoo et al. [5] presented PD detection, measurement, and classification which constitute an important tool for quality assessment of insulation systems utilized in HV power apparatus and cables.The patterns obtained with PD detectors contain characteristic features of the source/class of the respective PD process involved.Gao et al. [6] used the UHF methods for detecting PD owing to its sensitivity and antiinterference capability.
They modeled five types of typical PD sources to imitate the defects in GIS online.
Wu et al. [7] presented a method based on statistical parameters to optimize the algorithm of multi-PD.They used this method for estimation of residual breakdown voltage of generator bars and the optimized multiparameter diagnosis algorithm.Wu et al. [8] used a method to quantitatively determine the reliability of multiparameter diagnosis and to optimize the algorithm of multiparameter diagnosis.Gulski and Kreuger [9] used a computer-aided discharge analyzer, a combination of statistical and discharge parameters to discriminate between different discharge sources.

Modelling and Simulation in Engineering
Saha [10] presented cellulosic paper and oil insulation in a transformer degrade at higher operating temperature.
In this paper, the available PD charging detection methods include ultrasonic method, UHF method, and HFC method.The combined diagnosis of PD based on the multidimensional characteristic parameters is brought by Sugeno fuzzy integral to process the collected signals.A large amount of experimental data shows that the method for PD signal has well diagnostic accuracy and is able to judge the internal insulation of electrical equipment.(

PD Extraction
(B) Average PD Quantity avg .The total discharge quantity within cycle/discharging times is (C) Average Discharge Current .The discharge quantity per unit time (one cycle) is (D) Mean Square Rate.The mean square value of PD pulse discharge quantity per unit time (one cycle) is (E) NQN.This parameter reflects the number of PD pulses and the discharge intensity.
In (5),   , , and  represent the number of pulses per second in the quantization window , the number of amplitude windows, and the gain of PD detector, respectively.FS is the maximum number of quantization window in the unit mV gain.
(F) max 95.It is the maximum of the discharge quantity with 5% of the maximum PD quantity of discharge pulse removed.
(G) avg 95.It is the average of the discharge quantity with 5% of the maximum PD quantity of discharge pulse removed.
(H) max 10, max 20, max 50, max 100, max .These are the minimum amplitude of xx maximum PD within each cycle (in mV).
(K) avg max 10, avg max 20, avg max 50, avg max 100, avg max .These are the average amplitude of xx maximum PD within each cycle (in mV).
The parameters of rank (A) to (E) are values obtained which were assumed in the absence of interference.Actual testing will inevitably introduce severe interference.Therefore, based on a large number of experiments and statistical data, four very valuable parameters of rank (F) to (K) were proposed.A lot of statistics show that, in the field detection with interference, the credibility of the insulation condition reflected by parameters above is greatly improved compared to conventional parameters [11].Table 1 shows the extraction of partial basic characteristic parameters of PD.

Phase Characteristic Parameters.
When a PD abnormal signal was found, the type of PD needs to be determined.Judging the type of PD is generally based on PD phase distribution spectrums, which mainly include phase distribution of the maximum discharge quantity  max (Φ); phase distribution of the average discharge quantity   (Φ); phase distribution of discharge times   (Φ); amplitude distribution of PD () ; energy distribution of PD () .There are totally five two-dimensional discharging spectrums and a threedimensional discharging spectrum   (Φ, ).Suppose the measured number of fundamental frequency cycle is , and each frequency cycle is divided into  phase windows, with each phase window represented by Φ  ; then it is clear that (Φ  ) is the sum of the discharge numbers in the phase window Φ  , and (  ) and (  ) are the number of discharge times and energy where the discharge quantities are   and   , respectively.In (7),  is the discharge number of the phase window Φ  .In this system,  is the maximum number of fundamental frequency cycle intercepted from actual data.Discharge quantity is divided into 25 intervals.According to the above five characteristic spectrums based on phase analysis, PD can be derived from the phase characteristic parameters as shown in Table 2. Also (a) Sk + is defined as below: in (8), , , , and  are discrete value, probability of , average value of distribution, and the standard deviation of the distribution, respectively, (b) the definition of Sk − is in the following: (d) in the following equation, the definition of Ku − has been presented: (e, f) the number of local peaks Pe + and the number of local peaks Pe − are the ratio of the pulses' number and the number of phase intervals in positive and negative half cycle of atlas, respectively, (g) discharge asymmetry  is the ratio of discharge corresponding parameters' average value in positive half cycle and negative half cycle: +  and (Φ) are the sum of corresponding parameter values distributed in the positive half cycle and negative half cycle. + and  − are the number of nonzero value distributed in positive half cycle and negative half cycle of the voltage, (h) correlation factor (cc) is used to assess the difference on shape and distribution between  + (Φ) and  − (Φ) [11]: and  are discharge parameters in positive and negative half-cycle phase window.Also,  is the number of phase window of each half cycle [12].PD generated by different types of insulation defects have a significant difference in the phase distribution atlas, and the fingerprint database established by the simulation test in laboratory can be taken as a reference in detecting and identifying a defect at the scene.Table 3 shows the PD data phase characteristic parameters.
Use the information gain feature fusion method for dimension reduction

Results
Enter the HFC PD HFC PD characteristic parameters Diagnosis using soft computing method

High-Frequency Current PD Data's Diagnostic Method Based on Information Fusion and Semisupervised Learning.
In this paper, for the diagnosis of high-frequency current PD data, a method based on information fusion and semisupervised learning is used.Figure 1 shows the flow chart of this method.
The implementation steps of the algorithm can be described as follows.
Table 4 shows the four types of discharge model and the recognition result of semisupervised learning.From Table 4, it is clear that the accuracy rate of spike discharge is 79%, particulate discharge is 88%, suspended discharge is 85%, and air discharge is 83%.

Ultrasonic Signals Diagnostic Method Based on Adaptive
Mutation Parameters of Particle Entropy.In this paper, the parameters are optimized with a method based on particle entropy parameter adaptive mutation and then used for PD ultrasound data' diagnosis.Figure 2 shows the flow chart of this method.Process of parameters' adaptive mutation algorithm based on particle entropy is as follows: (1) Randomly determine the initial position and velocity of each particle in the data field, and set  0 as the stable threshold value of particle entropy.
(2) Set the particles'  best to the current optimal position, and set   as the global optimal position of initial population.
(3) Update the position and velocity of the particles.
(4) If the particle fitness is better than the fitness of   , update   with the current position.
(5) Calculate the collection of particle entropy and determine whether the entropy of each particle is less than the predetermined threshold value  0 , and perform (6) when it is true, otherwise skip.

UHF Signal Diagnostic Method Based on IIA-ART2A.
Parameters of ART2 neural network (NN) are optimized with a method based on immune algorithm with the improved variation of natural cycle and then used for the UHF PD data diagnostic.Figure 3 shows the flow chart of this method based on NN.IIA-ART2A NN is unsupervised learning.If we select the correct input vector, only a few samples, then we can get high recognition accuracy [14,15].One thousand new sample data set are selected in each type of defect, characteristic parameters are extracted, and pattern recognition with IIA-ART2A NN is done [16].Table 5 shows the recognition results, where the fractal dimension numerical lists only the 4th node of 4floor wavelet packet decomposition, and matching nodes 1-4 represent spike, particulate, suspended, and air, respectively.It is clear that the recognition rate of spike discharge is 91%, PD is 85%, suspended discharge is 84%, and air discharge is 77%.Combining a plurality of different classifiers to obtain high accuracy is an important research topic.There is interaction between the classifiers, rather than being independent.Fuzzy integral is nonlinear integral based on fuzzy measurement, and fuzzy measurement is a nonnegative nonadditive set function, while the nonadditive characteristic of fuzzy measurement can precisely describe interactions between classifiers.Therefore, exploring fusion technology of multiple classifiers based on fuzzy integral is the content of this article.Let  = { 1 ,  2 , . . .,   } be a collection of  target categories and  = { 1 ,  2 , . . .,   } is a collection of  classifiers.  is the th identified object.After sample   is identified by each classifier, we can define a matrix referred to as PD(  ), which is the cross section in decision-making model [12].
Each column vector ℎ  = (ℎ  1 , ℎ  2 , . . ., ℎ   ),  = 1, 2, . . ., , represents recognition results of each classifier on the sample   in category   .When the sample is fixed, ℎ  can be thought of as a function ℎ  :  → [0, 1] (if the output of the classifier is not in the interval [0, 1], it can be normalized to meet the conditions), which maps the classifier   to the corresponding components ℎ  of fusion vector   .Each intersection ℎ   of the output vector   and the fusion vector   represents the degree of determination on the fact that the class   assigned sample   to category   , also known as the objective estimate of the determination that   belongs to   of the classifier   .When ℎ   = 1, the classifier   determines that   belongs to the category   ; on the contrary, when ℎ   = 0,   confirms that   does not belongs to category   .Let  be the Sugeno fuzzy measurement of the defined power set () of , fuzzy measurement on a single set of points (the fuzzy density).
∀ ∈ (), () represents degree of reliability of local decision made by a subset  of .Fuzzy integral fuses the objective evaluation of the determination that   belongs to   (fusion vectors of   ) and the degree of reliability of classifiers credibility; the integral value is the total objective evaluation of the system that the sample belongs to category   .Thus, the system has integrated value for each category, taking the corresponding category of maximum integrated value as the system's determined category of sample   .
Therefore, the basic steps of fusion with fuzzy integral can be summarized as follows: (1) Determine fuzzy density.
(2) Construct fuzzy measurements.Parameters are determined through the fuzzy density, and then fuzzy measurements are also determined.
(3) Calculate fuzzy integral, the integrated extent of determination of the fused object to be identified belonging to each category.
(4) Size comparison of the integral value is carried out, to determine the category of the object to be identified.
In order to facilitate the integration of fuzzy arithmetic, we must first set the output of the probability of failure, provided   ( = 1, 2, . . ., ) the corresponding characteristic parameters' values of sample data and the th fault; the occurrence probability of th fault is defined as And 0 <   ≤ 1 and   has a clear sense; it reflects the difference between the binding energy of sample and the th fault and the binding energy of sample and the fault second to th.When the two values are equal   = 0.5, critical state   > 0.5, it can be classified as  type, and the greater the   (1) Category 1: for one thousand training samples, C1 classifier is constructed in accordance with the proposed HFC PD data's diagnostic method based on information fusion and semisupervised learning.
(2) Category 2: for one thousand training samples, C2 classifier is constructed in accordance with the proposed ultrasonic signals diagnostic method based on adaptive mutation parameters of particle entropy.
(3) Category 3: for one thousand training samples, C3 classifier is constructed in accordance with the proposed UHF signal diagnostic method based on IIA-ART2A NN.
Table 6 shows the accuracy rates obtained when identifying four types of samples.
We choose a test sample as an example so the sample is identified by the C1, C2, and C3 classifiers and then we calculate the probability of its belonging to A, B, C, and D categories; the results are shown in Tables 7, 8, and 9 and then directly classified in Table 7; it is divided into class D by C1 and C3 and into A category by C2.Table 8 shows the results of fuzzy measure and fuzzy integral, within which ℎ() is measurable function; the results of the descending corresponding item are in Table 8; (⋅) is the fuzzy measurement calculated according to the rank of fuzzy density and measurable function ℎ() and  V express the Sugeno fuzzy integral value of ℎ() and (⋅), wherein the fuzzy density is calculated from the classification accuracy.Make judgments based on the value of  V , and then the test sample belongs to class D. As can be seen, the fuzzy integral gets the "nonlinear" mean-value of ℎ().Since fuzzy measurements are from the classification accuracy of different categories, taking them as references when calculating the mean-value of multiclassification can reduce the uncertainty of multiclassifier [13].
Take PD as blurred density, to make fuzzy integral decision for various types of PD samples; recognition results are shown in Table 9; compared with Tables 4-8, it can be

Field Application of Combined Diagnosis of PD Based on the Multidimensional Characteristic Parameters
Detect PD on a 1000 kV high-voltage reactor, field test with the integrated application of UHF, and ultrasonic and highfrequency method, and then implement the multiple classifier fusion analysis based on Sugeno fuzzy integral.Test results show that there are two distinct PD phenomena in the high resistance; positioning discharge source 1 is located in the upper clip part of reactor column X (the body of 1) and discharge source 2 is located at the bottom of the reactor winding.The accuracy of the above PD diagnosis and localization is verified by the disintegrated inspection [14].
(1) With the detecting method of high-frequency PD pulse current, it is detected that, at the clip of body 1 of the high resistance iron core clamps grounding line, the amplitude of high-frequency PD signal is the highest, reaching 500 mV; at the iron core of body 1, the amplitude reaches 300 mV; at the clip and the iron core of body 2, the amplitude reaches 250 mV.PD phase spectra shows the characteristic of symmetric distribution in first and third quadrants, the high-frequency PD pulse current testing results are shown in Figure 5.
(2) Detecting with the UHF method, the PD detecting results are shown in Figure 6.In PRPD spectra, it is obvious that there are two discharge sources, which are shown in the red circle in the first and third quadrants and green circle in the twice and fourth quadrants, respectively.
(3) Locate the source of discharge.The discharging source is located at a position 136 cm away from the east side of the reactor's wall of box, 200 cm away from the north side of the box's wall, and at a height of 391 cm.The discharge source is located at a middle and high position of the box, close to the east and west box's wall, and the variation trend keeps consistent with UHF detection with normalized amplitude.
(4) Significant discharge signals are detected on the reactor by the ultrasonic sensor.
(5) Two discharging points have been found in disintegrated inspection of the high resistance, and the results are consistent with the PD detection, as follows.
Discharging Point 1. Iron core's lead insulation of high resistance column A of A phase collides with the upper clamps (Figure 7); carbon footprint is found at the side of the clip (Figure 8); there are totally six corresponding lead insulation crepe papers, wherein the outer three layers of crepe paper have disruptive discharge phenomena (Figures 9 and 10).
Traces of discharge at the position of iron core's lead insulation colliding with the upper clamps may be the oil gap's PD under a high electric field (the upper end of column A is in the operation of 500 kV voltage level, which will form a strong electric field on the outer surface of the iron core's lead insulation; and by testing the standby phase, it can be calculated that potential difference of the iron core's lead and clamp's lead of column A is significantly greater than column X under the rated working voltage).
Discharging Point 2. The shield cap of fixing bolt is not installed in place and collides with the fixing bolts at the bottom of the column X clamps.There exists carbon footprint (Figure 11).Because the shield cap is not installed in place, leading to being too close to clip bolts, and contact discharge happened under running vibration.The discharge here and the ultrasonic testing result are in agreement and also in agreement with the high-frequency detection result of clip and iron core ground loop anastomosis.

Conclusion
This article proposes nine PD basic characteristic parameters to reflect insulation of device and 8 phase characteristic parameters to reflect the type of discharge.The collected signals are processed with the method of high frequency, UHF, and ultrasonic.With the multiclassifier fusion method based on Sugeno fuzzy integral, integrate all kinds of defects  fault probability matrix, namely, the diagnosis results of high-frequency current PD data's diagnostic method based on information fusion and semisupervised learning ultrasonic signals diagnostic method based on adaptive mutation parameters of particle entropy and UHF signal diagnostic method based on IIA-ART2A NN.The results show that this method is significantly better than single detection method when diagnosing defects of spike, particulate, suspended, and air.

2. 1 .
Basic Characteristic Parameters.The basic parameters of PD are as follows.(A) Maximum PD Quantity max .The maximum discharge quantity within one sample data is max  = max {  } .

Figure 1 :
Figure 1: The flow chart of diagnosis method.

Table 3 :
Partial PD phase characteristic parameters.

Figure 3 :Figure 4 :
Figure 3: The flow chart of UHF signal diagnostic method.

Fuzzy
density is defined based on the accuracy of the method, and then the fuzzy density matrix is obtained as follows:

Figure 10 :
Figure 10: After poking iron core's lead insulation of column A.
The position of contact Traces of carbon wiped off

Figure 11 :
Figure 11: The lower lead supporting shield fixing bolts cap of high resistant's column X.

Table 1 :
Extraction of parameters of PD.

Table 2 :
List of characteristics parameters.
[13]ne the initial  ×  matrix .Steps of iterative solver are as follows[13]:(1) Construct undirected weighted graph  = (, );  is a collection of samples and  is a set of edges of figure .

Table 4 :
The recognition result.

Table 5 :
The recognition result of IIA-ART2A NN.

Table 6 :
Accuracy of the classifiers.

Table 7 :
The conversion of probability on sample characteristics among different classifiers.