 Open Access
 Total Downloads : 53
 Authors : Ms. Priyanka Kothoke , Dr. Anupama Deshpande , Mr. Yogesh Chaudhari
 Paper ID : IJERTV8IS080182
 Volume & Issue : Volume 08, Issue 08 (August 2019)
 Published (First Online): 28082019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Analysis and Determination of Partial Discharge Type using Statistical Techniques and Back Propagation Method of Artificial Neural Network for PhaseResolved Data
Ms. Priyanka Kothoke
Research Scholar Shri J. J. T University Rajasthan, India
Dr. Anupama Deshpande
Research Guide Shri J. J. T University
Rajasthan, India
Mr. Yogesh Chaudhari Assistant Professor Navrachana University Vadodara, Gujarat
Abstract : Partial discharge (PD) patterns are significant tool for the diagnosis of high voltage (HV) insulation systems. Human experts can discover possible insulation defects in various representations of the PD data. One of the most widely used representations is phaseresolved PD (PRPD) patterns. In order to ensure reliable operation of HV equipment, it is crucial to relate the observable statistical characteristics of PDs to the properties of the defect and ultimately to determine the type of the defect. In present work, we have obtained and analysed combined use of PRPD patterns (q), (n) and (nq) using statistical parameters such as skewness and kurtosis for two patterns viz. (q) and (n) along with mean, standard deviation, variance, skewness and kurtosis for (nq) to detect type of PD and we have verified the obtained results by providing obtained statistical parameters as an input for training of artificial neural network (ANN) in MATLAB tools.
Keywords – Partial Discharge, Phaseresolved, Statistical parameters, Artificial Neural Network.

INTRODUCTION
PD is a localized electrical discharge that partially bridges the insulation between conductors and which may or may not occur adjacent to a conductor [1]. In general, PDs are concerned with dielectric materials used, and partially bridging the electrodes between which the voltage is applied. The insulation may consist of solid, liquid, or gaseous materials, or any combination of them. PD is the main reason for the electrical ageing and insulation breakdown of high voltage electrical apparatus. Different sources of PD give different effect on insulation performance. Therefore, PD classification is important in order to evaluate the harmfulness of the discharge [2].
PD classification aims at the recognition of discharges of unknown origin. For many years, the process was performed by investigating the pattern of the discharge using the well known ellipse on an oscilloscope screen, which was observed crudely by eye. Nowadays, there has been extensive published research to identify PD sources by using intelligent technique like artificial neural networks, fuzzy logic and acoustic emission [2].
The recent upsurge of research on PD phenomena has been driven in part by development of new fast digital and computerbased techniques that can process and analyze signals derived from PD measurements. There seems to be
an expectation that, with sufficiently sophisticated digital processing techniques, it should be possible not only to gain new insight into the physical and chemical basis of PD phenomena, but also to define PD patterns that can be used for identifying the characteristics of the insulation defects at which the observed PD occur [3]. One of the undoubted advantages of a computeraided measuring system is the ability to process a large amount of information and to transform this information into an understandable output [4]. There are many types of patterns that can be used for PD source identification. If these differences can be presented in terms of statistical parameters, identification of the defect type from the observed PD pattern may be possible. As each defect has its own particular degradation mechanism, it is important to know the correlation between discharge patterns and the kind of defect [4]. Therefore, progress in the recognition of internal discharge and their correlation with the kind of defect is becoming increasingly important in the quality control in insulating systems [5]. Research has been carried out in recognition of partial discharge sources using statistical techniques and neural network [6]. We present a method for the automated recognition of PRPD patterns using an ANN for the actual classification task [12]. In our study, we have tested various internal and external discharges like void, surface and corona using statistical parameters such as skewness and kurtosis for (q) and (n) and mean, standard deviation, variance, skewness and kurtosis for (nq) as a initial step known as preprocessing and a rough draft is made about the results of PD Type and verified the results by using back propagation method of
Artificial Neural Network in MATLAB software.

STATISTICAL PARAMETERS
The important parameters to characterize PDs are phase angle , PD charge magnitude q and PD number of pulses n. PD distribution patterns are composed of these three parameters. Statistical parameters are obtained for phase resolved patterns (q), (n) and (nq).

Processing of data
The data to be processed obtained from generator includes , q, n and voltage V. From this data, phase resolved patterns are obtained.

Analysis of PhaseResolved (q) and (n) using Statistical Techniques
Fig.1 (a). Block diagram of discharge analysis for (q)
Fig.1 (b). Block diagram of discharge analysis for (n)
PD pulses are grouped by their phase angle with respect to 50 (Â± 5) Hz sine wave. Consequently, the voltage cycle is divided into phase windows representing the phase angle axis (0 to 360). If the observations are made for several voltage cycles, the statistical distribution of individual PD events can be determined in each phase window. The mean values of these statistical distributions results in two dimensional patterns of the observed PD patterns throughout the whole phase angle axis [6]. A twodimensional (2D) distribution q and n represents PD charge magnitude q and PD number of pulses n as a function of the phase angle [3].
+ +
+ +
– –
– –
The mean pulse height distribution Hqn () is the average PD charge magnitude in each window as a function of the phase angle . The pulse count distribution Hn () is the number of PD pulses in each window as a function of phase angle . These two quantity are further divided into two separate distributions of the negative and positive half cycle resulting in four different distributions to appear: for the positive half of the voltage cycle Hqn () and Hn () and for the negative half of the voltage cycle Hqn () and Hn () [5]. For a single defect, PD quantities can be described by the normal distribution. The distribution profiles of Hqn () and Hn () have been modeled by the moments of the normal distribution: skewness and kurtosis.
N (xi Âµ)3f(xi)
Skewness (S ) = i=1 (1)
k 3 N f(x )
i=1 i
N (xi Âµ)4f(xi)
Kurtosis: (K ) = i=1 3 (2)
u 4 N f(x )
i=1 i
where,
f(x) = PD charge magnitude q, = average mean value of q, = variance of q.
Skewness and Kurtosis are evaluated with respect to a reference normal distribution. Skewness is a measure of
asymmetry or degree of tilt of the data with respect to normal distribution. If the distribution is symmetric, Sk=0; if it is asymmetric to the left, Sk>0; and if it is asymmetric to the right, Sk<0. Kurtosis is an indicator of sharpness of distribution. If the distribution has same sharpness as a normal distribution, then Ku=0. If it is sharper than normal, Ku>0, and if it is flattr, Ku<0 [3] [7].

Analysis of PhaseResolved (qn) using Statistical Techniques

Fig. 2. Block diagram of discharge analysis for (nq)
Where,
S.D = standard deviation Sk = skewness
Ku = kurtosis
Statistical analysis is applied for the computation of several statistical operators. The definitions of most of these statistical operators are described below. The profile of all these discrete distribution functions can be put in a general function, i.e., yi=f(xi). The statistical operators can be computed as follows:
N 1(xi) f(xi)
Mean Value: () = i= (3)
N f(xi)
i=1
N 2 f(xi)
Variance: (2) = i=1(xi) (4)
N f(xi)
i=1
Standard Deviation = (5)
where,
x = number of pulses n,
f(x) = PD charge magnitude q,
= average mean value of PD charge magnitude q, = variance of PD charge magnitude q
Skewness and Kurtosis are evaluated with respect to a reference normal distribution. Skewness is a measure of asymmetry or degree of tilt of the data with respect to normal distribution. If the distribution is symmetric, Sk=0; if it is asymmetric to the left, Sk>0; and if it is asymmetric to the right, Sk<0. Kurtosis is an indicator of sharpness of distribution. If the distribution has same sharpness as a normal distribution, then Ku=0. If it is sharper than normal, Ku>0, and if it is flatter, Ku<0 [3] [8].


RESULTS AND DISCUSSION OF STATISTICAL PARAMETERS
Analysis involves determining unknown PD patterns by comparing those with known PD patterns such as void, surface and corona. The comparison is done with respect to their statistical parameters [9] [10][11].

Analysis for (q)
The phase resolved patterns are divided into two types: (q) and (n). The phase resolved patterns (q) are obtained for three known PD patterns: void, surface and corona (as
discussed in 3.1.1) and three unknown PD patterns: data1, data2 and data3 (as discussed in 3.1.3) [9]

2D distribution of (q) for known PD patterns
Fig.3 (a), Fig.3 (b) and Fig.3 (c) are the phase vs. charge q plot for void, surface and corona discharges respectively.
Fig.3(a).Phase plot (
q) of void discharge
Fig.3(b).Phase plot (
q) of surface discharge
Fig.3(c). Phase plot (q) of corona discharge

Parameters of known PD patterns
TABLE I. PARAMETERS OF KNOWN PD PATTERNS
Parameter
Void
surface
corona
Skewness Hqn+ ()
1.0013
1.2134
0.3555
Skewness Hqn ()
0.9901
1.8219
1.3659
Kurtosis Hqn+ ()
2.9046
3.6064
2.4354
Kurtosis Hqn ()
2.7872
5.4506
7.5947

2D distribution of (q) for unknown PD patterns
Fig.4 (a), Fig.4 (b) and Fig.4 (c) are the phase vs. charge q plot for data1, data2 and data3 respectively.
Fig.4 (a) Phase plot (q) of data1
Fig.4 (b) Phase plot (q) of data2
Fig.4 (c) Phase plot (
q) of data3
From Fig.4 (a), it is seen that the following plot is similar to void and surface discharge. Fig.4 (b), is also similar to void and surface discharge and Fig.4(c), is similar to void discharge.

Parameters of unknown PD Patterns

TABLE II. PARAMETERS OF UNKNOWN PD PATTERNS
Parameter
data1
data2
data3
Skewness Hqn+ ()
0.8991
0.7456
1.0013
Skewness Hqn ()
1.1833
0.8509
0.9901
Kurtosis Hqn+ ()
2.5719
2.1814
2.9046
Kurtosis Hqn ()
3.7467
2.6512
2.7872

Analysis for (n)
The phase resolved (n) patterns consist of three known PD patterns: void, surface and corona (as discussed in 3.2.1) and three unknown PD patterns: data1, data2 and data3 (as discussed in 3.2.3) [9]. The plots are discussed below:
3.2.1 Phase resolved plot (n) of known PD patterns
Fig.5(a).Phase plot (n) of void discharge
Fig.5(b).Phase plot (n) of surface discharge
Fig.5(c).Phase plot (n) of corona discharge
Fig.5(a).Phase plot (n) of void discharge
Fig.5(b).Phase plot (n) of surface discharge
Fig.5(c).Phase plot (n) of corona discharge
Fig.5 (a), Fig.5 (b) and Fig.5 (c) are the phase vs. number of pulses n for void, surface and corona discharges respectively.

Parameters of unknown PD Patterns
TABLE III. PARAMETERS OF KNOWN PD PATTERNS
Parameter
Void
surface
corona
Skewness H + ()
n
0.4954
1.0082
1.3942
– ()
Skewness Hn
0.4329
2.3686
1.3798
Kurtosis H + ()
n
2.0535
2.871
4.8337
– ()
Kurtosis Hn
1.9137
8.4788
7.3215

Phase resolved plot (n) of unknown PD patterns
Fig.6(a).Phase plot (
n) of data1
Fig.6(b).Phase plot (n) of data2
Fig.6(c).Phase plot (
n) of data3
Fig.6(a).Phase plot (
n) of data1
Fig.6(b).Phase plot (n) of data2
Fig.6(c).Phase plot (
n) of data3
Fig.6 (a), Fig.6 (b) and Fig.6 (c) are the phase vs. number of pulses n plot for data1, data2 and data3 respectively.
From Fig.6 (a), it is seen that the following plot is similar to void and surface discharge. Fig.6 (b), is similar to void discharge and Fig.6 (c), is also similar to void discharge

Parameters of unknown PD patterns

TABLE IV. PARAMETERS OF UNKNOWN PD PATTERNS
Parameter
data1
data2
data3
+ ()
Skewness Hn
0.8016
0.574
0.4954
Skewness H – ()
n
1.0169
0.42
0.4329
Kurtosis H + ()
n
2.3724
2.1091
2.0535
– ()
Kurtosis Hn
3.2011
1.8003
1.9137

Analysis for (nq)
The phase resolved patterns nq are obtained for three known PD patterns: void, surface and corona (as discussed in 3.3.1) and three unknown PD patterns: data1, data2 and data3 (as discussed in 3.3.3) [10].

2D distrbution of nq for known PD patterns
Fig. 7(a), Fig. 7(b), Fig. 7(c), Fig. 7(d) and Fig. 7(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for void discharge respectively.
Fig.7(a).Mean plot (n
q) of void discharge
Fig.7(b).Standard deviation plot (nq) of void discharge
Fig.7(c).Variance plot (nq) of void discharge
Fig.7 (d) Skewness plot (n
q) of void discharge
Fig.7 (e) Kurtosis plot (nq) of void discharge
Fig.9 (d) Skewness plot (nq) of corona discharge
Fig.9 (e) Kurtosis plot (nq) of corona discharge
Fig.7 (d) Skewness plot (n
q) of void discharge
Fig.7 (e) Kurtosis plot (nq) of void discharge
Fig.9 (d) Skewness plot (nq) of corona discharge
Fig.9 (e) Kurtosis plot (nq) of corona discharge
Fig.9 (a) Mean plot (n
q) of corona discharge
Fig.9 (b) Standard deviation (nq) of corona discharge
Fig.9 (c) Variance plot (nq) of corona discharge
Fig.9 (a) Mean plot (n
q) of corona discharge
Fig.9 (b) Standard deviation (nq) of corona discharge
Fig.9 (c) Variance plot (nq) of corona discharge
Referring to Fig. 7 (a), Fig. 7 (b) and Fig. 7 (c) of void discharge, it can be seen there is a peak occurring somewhere after 1500 cycle, which is a void discharge and in Fig. 7 (d) and Fig. 7 (e) of skewness and kurtosis, the value decreases at that cycle where peak occurs.
Fig. 8(a), Fig. 8(b), Fig. 8(c), Fig. 8(d) and Fig. 8(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for surface discharge respectively.
Fig.8 (a) Mean plot (n
q) of surface discharge
Fig.8 (b) Standard deviation plot (nq) of surface discharge
Fig.8 (c) Variance plot (nq) of surface discharge
Fig.8 (d) Skewness plot (nq) of surface discharge
Fig.8 (e) Kurtosis plot (nq) of surface discharge
In surface discharge, charges are distributed uniformly over all cycles for mean, standard deviation, variance, skewness and kurtosis as shown in Fig. 8(a), Fig. 8(b), Fig. 8(c), Fig. 8(d) and Fig. 8(e).
Fig. 9(a), Fig. 9(b), Fig. 9(c), Fig. 9(d) and Fig. 9(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for corona discharge respectively.
Referring to Fig. 9(a), Fig. 9(b) and Fig. 9(c) of corona discharge, it can be seen the charges starts occurring after 500 cycle increasing somewhere upto 1200 cycle and then decreasing after 2000 cycle, and in Fig. 9(d) and Fig. 9(e) of skewness and kurtosis, the value decreases from 500 cycle till 2000 cycle.

Parameters of known PD patterns
TABLE V. PARAMETERS OF KNOWN PD PATTERNS
Parameters
Void
Surface
Corona
Mean
13320.32
145.706
1.426
Standard deviation
7553.716
126.009
1.139
Variance
1.64*108
17921.01
2.279
Skewness
3.66*1017
0.809
0.26
Kurtosis
0.04878
2.442
0.966

2D distribution of (nq) for unknown PD patterns
Fig. 10(a), Fig. 10(b), Fig. 10(c), Fig. 10(d) and Fig. 10(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for data1 respectively.
Fig.10 (a) Mean plot (nq) of data1
Fig.10 (b) Standard deviation plot (nq) of data1
Fig.10 (c) Variance plot (nq) of data1
Fig.10 (d) Skewness plot (nq) of data1
Fig.10 (e) Kurtosis plot (nq) of data1
In Fig. 10(a), Fig. 10(b), Fig. 10(c), Fig. 10(d) and Fig. 10(e), the charges are uniformly distributed similar to surface discharge. Hence, it can be concluded that data1 is having
surface discharge. Fig. 11(a), Fig.11(b), Fig. 11(c), Fig. 11(d) and Fig. 11(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for data2 respectively.
Fig.11 (a) Mean plot (n
q) of data2
Fig.11 (b) Standard deviation plot (nq) of data2
Fig.11 (c) Variance plot (nq) of data2
Fig.11 (d) Skewness plot (nq) of data2
Fig.11 (e) Kurtosis plot (nq) of data2
In Fig. 11(a), Fig. 11(b), Fig. 11(c), Fig. 11(d) and Fig. 11(e), the charges are uniformly distributed similar to surface discharge. Hence, it can be concluded that data2 is having surface discharge.



OBSERVATIONS FROM STATISTICAL METHODS
From the above results, following observations are made:
+
+
–
–

Fig. 13(a), Fig. 13(b) and Fig. 13(c) are the characteristics of skewness and kurtosis (Hqn () and Hqn ()) of data 1, data2 and data3 against void., surface and corona discharges respectively.

From Fig. 13(a), it is observed that data3 characteristics overlaps void discharge characteristics, it can be concluded that data3 is void discharge. Data2 characteristics approximately fits against void, it can be concluded that data2 is also void discharge

From Fig. 13(b), it is observed that data 1 characteristics is similar to surface discharge characteristics, it can be concluded that data 1 is surface discharge.

From Fig. 13(c), it is observed that none of the data characteristics is similar to corona discharge characteristics, it can be concluded that none of the data has corona discharge.
Fig.12 (a) Mean plot (n
q) of data3
Fig.12 (b) Standard deviation plot
+ –
+ –
(nq) of data3
Fig.12 (c) Variance plot (nq) of data3
Fig.12 (d) Skewness plot (nq) of data3
Fig.12(e).Kurtosis plot (nq) of data3
Fig.13 (a) Characteristics of skewness and kurtosis(H +() and H –
qn qn
()) of data1, data2, and data3 against void
Fig.13 (b) Characteristics of skewness and kurtosis(H +() and
qn
H ()) of data1,
qn
data2, and data3 against surface
Fig.13 (c) Characteristics of skewness and kurtosis(H +() and
qn
H ()) of data1,
qn
data2, and data3 against corona
Fig.13 (a) Characteristics of skewness and kurtosis(H +() and H –
qn qn
()) of data1, data2, and data3 against void
Fig.13 (b) Characteristics of skewness and kurtosis(H +() and
qn
H ()) of data1,
qn
data2, and data3 against surface
Fig.13 (c) Characteristics of skewness and kurtosis(H +() and
qn
H ()) of data1,
qn
data2, and data3 aainst corona
Fig. 14(a), Fig. 14(b) and Fig. 14(c) are the characteristics of skewness and kurtosis (Hn () and Hn ()) of data 1, data2 and data3 against void., surface and corona discharges respectively.
Fig. 12(a), Fig. 12(b), Fig. 12(c), Fig. 12(d) and Fig. 12(e) are the nq plot of mean, standard deviation, variance, skewness and kurtosis for data3 respectively. In Fig. 12(a), Fig. 12(b) and Fig. 12(c), there is a occurrence of peak after 1500 cycle and in Fig. 12(d) and Fig. 12(e), the skewness and kurtosis value decreases at that peak which is similar to void discharge. Hence, it can be concluded that data3 is void discharge.

Parameters of unknown PD patterns
Parameters
data1
data2
data3
Mean
105.119
553.93
13320.32
Standard deviation
97.966
698.3
7553.716
Variance
16714.23
4.94*105
1.64* 108
Skewness
0.692
1.939
3.7*1017
Kurtosis
2.004
6.31
0.04878
Parameters
data1
data2
data3
Mean
105.119
553.93
13320.32
Standard deviation
97.966
698.3
7553.716
Variance
16714.23
4.94*105
1.64* 108
Skewness
0.692
1.939
3.7*1017
Kurtosis
2.004
6.31
0.04878
TABLE VI. PARAMETERS OF UNKNOWN PD PATTERNS
Fig.14 (a) Characteristics of skewness and
+() and H –
kurtosis(Hn n
()) of data1, data2, and data3 against void
Fig.14 (b) Characteristics of skewness and
+() and
kurtosis(Hn
H ()) of data1,
n
data2, and data3 against surface
Fig.14 (c) Characteristics of skewness and
+() and H –
kurtosis(Hn n
()) of data1, data2, and data3 against corona
Fig.14 (a) Characteristics of skewness and
+() and H –
kurtosis(Hn n
()) of data1, data2, and data3 against void
Fig.14 (b) Characteristics of skewness and
+() and
kurtosis(Hn
H ()) of data1,
n
data2, and data3 against surface
Fig.14 (c) Characteristics of skewness and
+() and H –
kurtosis(Hn n
()) of data1, data2, and data3 against corona
From Fig. 14(a), it is observed that data3 characteristics overlaps void discharge characteristics, it can be concluded that data3 is void discharge. Data2 characteristics approximately fits against void, it can be concluded that data2 is also void discharge. From Fig. 14(b), it is observed that data 1 characteristics is close to surface discharge characteristics, it can be concluded that data 1 is surface discharge. From Fig. 14(c), it is observed that none of the
data characteristics is similar to corona discharge characteristics, it can be concluded that none of the data has corona discharge.
Fig.15 (a) Statistical Characteristics of data3 against void discharge
Fig.15 (b) Statistical characteristics of data1 against surface discharge
Fig.15 (c) Statistical characteristics of data2 against surface discharge
Fig.15 (a) Statistical Characteristics of data3 against void discharge
Fig.15 (b) Statistical characteristics of data1 against surface discharge
Fig.15 (c) Statistical characteristics of data2 against surface discharge
Fig. 15(a) is the statistical characteristics of mean, standard deviation, variance, skewness and kurtosis of void discharge against data3. Fig. 15(b), Fig. 15(c) are the statistical characteristics of mean, standard deviation, variance, skewness and kurtosis of surface discharge against data1 and data2 respectively.
Plotting statistical parameters of void discharge against data3 in Fig. 15(a) shows data3 characteristics overlaps void characteristics, it can be concluded that data3 is void discharge. Similarly, for surface discharge, data1 and data2 characteristics (Fig. 15(b) and Fig. 15(c)) approximately fit surface discharge characteristics, it can be concluded that data1 and data2 is surface discharge


DISCUSSION ON RESULTS OF STATISTICAL PARAMETERS
From all the above observations using all the three phase resolved patterns (q), (n) and (nq), it can be concluded that, data1 is surface discharge in all the three phaseresolved patterns, data2 is surface discharge in (nq) pattern and void discharge in (q) and (n), hence data2 is both void and surface discharge and data3 is void discharge in all the three phaseresolved patterns.
From statistical parameters, the PD source cannot be concluded accurately so it needs to be applied to other classification methods such as neural network, Fuzzy logic etc. as a preprocessing parameters for getting accurate PD source.
However it is felt that in data2, the type of discharge is slightly vague and in data1, data3 the type can be further ascertained by using latest techniques of Neural Networks. Hence, method of BPM of ANN was tried.

BACK PROPAGATION METHOD OF ARTIFICIAL NEURAL NETWORK
This Neural Network is named as Back Propagation method owing to its method of dealing with the error.. It comprises of an input layer, at least one hidden layer and the output layer. Each layer is completely connected to the succeeding layer. The number of processing elements (neurons) in the input layer is equal to the number of feature parameters selected for classification. The number of neurons in the output layer is equal to the number of PD types classification [13].

Importing Data
Firstly, preprocessing of data is carried out using Statistical parameters like mean, standard deviation, variance, skewness and kurtosis and a rough result is made about PD type. For accurate detection of partial discharge type, nntool built in Matlab is used which gives better results in the following GUI (Graphical User Interface) or nntool window.
Fig. 16. Network/Data Manager Window
Input data statistical parameters are imported in column INPUT DATA and targets which are decided by the user are imported in column TARGET DATA as follows:
Fig. 17. Input data and Target data Imported in Network/Data Manager
Where,

P1 = Known data which are void, surface and corona having parameters such as mean, standard deviation, variance, skewness and kurtosis.

pp 1 cycles = Unknown data1 having parameters such as mean, standard deviation, variance, skewness and kurtosis.

pp 2 cycles = Unknown data2 having parameters such as mean, standard deviation, variance, skewness and kurtosis.

pp 3 cycles = Unknown data3 having parameters such as mean, standard deviation, variance, skewness and kurtosis.

p1 cycles = Unknown data2007 having parameters such as mean, standard deviation, variance, skewness and kurtosis.

p2 cycles = Unknown data2011 having parameters such as mean, standard deviation, variance, skewness and kurtosis.
 T1 = Own targets are decided [50, 60, 70].

T2 = Own targets are decided [60, 70, 50].

T3 = Own targets are decided [70, 50, 60].


Creating a New Network
Larger networks require large training datasets. Over fitting is more likely when the model is large.
Then, a network is created using NEW command in nntool feature of Matlab. Following window appears where network properties are set.
Fig. 20. Training Window
Where,
Fig. 18. Creating Network or Data Window
In train menu bar of network1, inputs and targets are fed in training data.

Inputs: P1

Targets: T1
Network1 is trained by clicking on Train Network. Neural Network Training (nntraintool) appears which displays the training of the network.

Network type: Feedforward back propagation

Input data: P1

Target data: T1

Training function: TRAINSCG

Adaption learning function: LEARNGDM

Performance function: MSE

Number of layers: 2

Properties for Layer1 are:

Number of neurons: 20

Transfer function: TANSIG
Then, using Create, NETWORK 1 is created.
Fig. 19. Network created window
Then, by clicking on network1, following window appears:


Training
Once a network has been structured for a particular application, that network is ready to be trained. To start this process the initial weights are chosen randomly. Then, the training, or learning, begins.
There are two approaches to training – supervised and unsupervised.
Fig. 21. Neural Network Training (nntraintool) Window

Simulation
After the neural network classifier has been trained, the next step is to determine whether the resulting model of neural network is with the set of best weights and accurate enough for putting on the test data? Although measures such as sum of squared deviations or crossentropy error can be used, a reasonable measure of accuracy is simply the percentage of correct predictions made by the model.
Fig. 22. Simulating Window
In Simulate menubar, inputs are fed in simulation data and output is obtained in simulation results.
Fig. 23. Simulated Output Window

Output
Classification of PD using neural networks is an important, fascinating and complex topic. The random initialization of network weights prior to each execution of the neural network training algorithm may cause final classification results to vary from execution to execution in some cases, even when all other factors (e.g., training data, learning rate, momentum, network topology) are kept constant. Particularly when working with very limited training datasets, the variation in results can be large. Under such circumstances, it is better to expand training data on the basis of improved ground truth. If this is not possible, generation of optimum results can sometimes be made through combination of the results of multiple neural network classifications.
Fig.24. Output Window
The simulated output obtained 59.9192, which is approximately close to the target of surface discharge because we have given target as T1 for training so in that 50 target is decided for void, 60 for surface and 70 for corona. If we are selecting T2 then 60 is for void, 70 for surface and 50 is corona and so on. Hence, it can be concluded that data1 is
surface discharge. Similarly, the simulated results are obtained for other unknown PD data.
TABLE VII: Outputs
Sr. No.
Unknown data
Partial discharge Source
1
Data 1
Surface
2
Data 2
Void/surface
3
Data 3
Void


CONCLUSION
The use of backpropagation algorithm is better choice for training ANN as it performs optimization of the weights by first propagating the input signal forward, then propagating the calculated error backwards to make adjustment to the initial weights.
The ANN on execution produce the same results as obtained in Statistical Method for data1 and data3 resulting in accuracy of the results but again produce dual results for data2 indicating the need for more input samples in training data set so need to use other advanced software Python and other methods of Artificial Neural Network such as Random Forest method, Self organizing map and Support Vector Machine so as to further increase the accuracy in detecting type of PD. From the present work, we can make observe that PD type can be classified using the phase resolved pattern and can be given as the input. It is also possible to detect multiple sources of breakdown those exist in insulation along with determining exact percentage of that type using self organizing map method of ANN in software Python. In future work we aim at developing ANN trained using methods other than back propagation.
REFERENCES

MICAMAXXTM plus Partial Discharge Basics

M. G. Danikas, The Definitions Used for Partial Discharge Phenomena, IEEE Trans. Elec. Insul., Vol. 28, pp. 10751081, 1993.

N.C. Sahoo, M. M. A. Salama, R. Bartnikas, Trends in Partial Discharge Pattern Classification: A Survey, IEEE Transactions on Dielectrics and Electrical Insulation, Vol. 12, No. 2; April 2005..
 [4] E. Gulski, J. Smith, R. Brooks, Partial Discharge Databases for Diagnosis Support of HV Components, IEEE Symposium on Electrical Insulation, pp. 424427, 1998

E. Gulski and F. H. Kreuger, Computeraided recognition of Discharge Sources, IEEE Transactions on Electrical Insulation, Vol. 27 No. 1, February 1002.

E. Gulski and A. Krivda, Neural Networks as a Tool for Recognition of Partial Discharges, IEEE Transactions on Electrical Insulation, Vol. 28 No.8, December 1993.

F. H. Kreuger, E. Gulski and A. Krivda, Classification of Partial Discharges, IEEE Transactions on Electrical Insulation, Vol. 28 No. 6, December1993.

C. Chang and Q. Su, Statistical Characteristics of Partial Discharges from a RodPlane Arrangement

Namrata Bhosale, Priyanka Kothoke Amol Deshpande, Dr. Alice Cheeran, Analysis of Partial Discharge using PhaseResolved(q) and (n) Statistical Techniques, International Journal of Engineering Research and Technology, Vol. 2 (05), 2013,ISSN22780181.

Priyanka Kothoke, Namrata Bhosale, Amol Deshpande, Dr. Alice Cheeran, Analysis of Partial Discharge using PhaseResolved (nq) Statistical Techniques, International Journal of Engineering Research and Applications.

Yogesh R. Chaudhari, Namrata R. Bhosale, Priyanka M. Kothoke, Composite Analysis of Phase Resolved PD Patterns using Statistical Techniques, International Journal of Modern Engineering Research (IJMER) Vol. 3, Issue. 4, pp19471957, 2013.

Christian Cachin, Hans Jiirg Wiesmann, PD Recognition with Knowledgebased Preprocessing and Neural Networks, IEEE Transactions on Dielectrics and Electrical Insulation Vol. 2, Issue 4, August 1995.

Priyanka Kothoke, Satyendra Kumar, Yogesh Chaudhari, Analysis of Partial discharge source using Artificial Neural Network, International Journal of Innovative Research in Science, Engineering and Technology (IJIRSET), Vol. 6, Issue 6, June 2017