Abstract

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.

1. 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 electric-field 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 anti-interference 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.

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.

2. PD Extraction

2.1. Basic Characteristic Parameters

The basic parameters of PD are as follows.

(A) Maximum PD Quantity . The maximum discharge quantity within one sample data is

(B) Average PD Quantity . 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) . It is the maximum of the discharge quantity with 5% of the maximum PD quantity of discharge pulse removed.

(G) . It is the average of the discharge quantity with 5% of the maximum PD quantity of discharge pulse removed.

(H) , , , , . These are the minimum amplitude of xx maximum PD within each cycle (in mV).

(K) , , , , . 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.

2.2. 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 ; phase distribution of the average discharge quantity ; phase distribution of discharge times ; amplitude distribution of ; energy distribution of . There are totally five two-dimensional discharging spectrums and a three-dimensional 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) 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 is in the following: (c) is as below:(d)in the following equation, the definition of has been presented: (e, f)the number of local peaks and the number of local peaks 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.

3. Single PD Diagnosis Algorithm

3.1. 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.

Given a data collection and a tag collection , the first l ones () are already labeled samples, labeled as , and the last ones () are unlabeled samples. The targets are the labels of predicted unlabeled samples. Define an matrix . () is the energy vector of and is tagged as . Define 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 .(2)Construct similar matrix .(3)Matrix with is calculated and is (4)The problem is solved by Iterative method.(5)When (), mark each as .

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%.

3.2. 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 as the stable threshold value of particle entropy.(2)Set the particles’ 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 , and perform (6) when it is true, otherwise skip.(6)Vary to , calculate the variation value of , and then continue to update the iteration.(7)Determine whether the convergence criteria of algorithm are met; if true, implement (8); otherwise, continue iterating.(8)Output and end the algorithm.

3.3. 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 4-floor 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%.

4. Multiple Classifier’s Fusion Analysis Based on Sugeno Fuzzy Integral

In this work, a multiple classifier fusion based on Sugeno fuzzy integral is used, to integrate the final diagnosis matrix of different PD diagnosis signals (HFC, ultrasound, and UHF) and calculate the final diagnosis result. Figure 4 shows the flow chart of this strategy.

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 be a collection of target categories and 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 , which is the cross section in decision-making model [12].And each row vector , , is the recognition results of a classifier on the sample in various categories, and we call it output vector of classifier .

Each column vector , , represents recognition results of each classifier on the sample in category . When the sample is fixed, can be thought of as a function (if the output of the classifier is not in the interval , 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 , the classifier determines that belongs to the category ; on the contrary, when , 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).

, , indicates the credibility of the decisions made by classification .

, 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 () the corresponding characteristic parameters’ values of sample data and the th fault; the occurrence probability of th fault is defined asAnd 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 , critical state , it can be classified as type, and the greater the value of the stronger the deterministic of classification. PD fault data of four states including spikes on the surface of high-voltage conductor, free metal particulates, suspended electrode, and a solid insulating gap take 1000 samples from each state, within which training samples change depending on the conditions, 500 test samples, and construct classifiers of three categories:(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.

Fuzzy density is defined based on the accuracy of the method, and then the fuzzy density matrix is ​​obtained as follows: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 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 , 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 48, it can be seen that the recognition accuracy of error-prone categories’ sample (such as C, D) has significantly improved with fuzzy integral fusion decision. Because differences among original characteristic vectors of error-prone subsample are small, the randomness and uncertainty of recognition process make it difficult to correctly classify. Multiple classifier fusion analysis based on fuzzy integral can avoid this kind of “uncertainty” and improve the recognition accuracy to some extent.

5. 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 high-frequency 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.

6. 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.

Competing Interests

The author declares that they have no competing interests.