T.L Fault Detection
T.L Fault Detection
net/publication/251998942
CITATIONS READS
28 11,343
3 authors:
R.P. Maheshwari
Indian Institute of Technology Roorkee
105 PUBLICATIONS 2,021 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by R.P. Maheshwari on 24 December 2013.
Abstract—Transmission line protection is an important issue in resistance, mutual-coupling from adjacent parallel lines, etc.
power system engineering because 85-87% of power system faults They do not have the ability to adapt dynamically to the system
are occurring in transmission lines. This paper presents a operating conditions, and to make correct decisions if the
technique to detect and classify the different shunt faults on a signals are uncertain.
transmission lines for quick and reliable operation of protection
schemes. Discrimination among different types of faults on the Fault detection and classification is a very challenging task.
transmission lines is achieved by application of evolutionary Different attempts have been made for fault classification
programming tools. including approaches based on traveling waves [1-2], adaptive
PSCAD/EMTDC software is used to simulate different Kalman filtering [3], fuzzy logic, neural networks[4], and the
operating and fault conditions on high voltage transmission line, fusion of different artificial intelligence techniques. Several
namely single phase to ground fault, line to line fault, double line researchers have proposed different techniques for fault
to ground and three phase short circuit. The discrete wavelet classification of transmission lines using different types of
transform (DWT) is applied for decomposition of fault transients, neural networks and their combination with different
because of its ability to extract information from the transient transforms, such as wavelet and hyperbolic-s [5]. Although the
signal, simultaneously both in time and frequency domain. The neural-network based approaches have been quite successful in
data sets which are obtained from the DWT are used for training determining the correct fault type, the main disadvantage of
and testing the SVM architecture. After extracting useful neural-network is that it requires a considerable amount of
features from the measured signals, a decision of fault or no fault
training effort for good performance, especially under a wide
on any phase or multiple phases of a transmission line is carried
out using three SVM classifiers. The ground detection task is
variation of operating conditions (such as system loading level,
carried out by a proposed ground index. Gaussian radial basis fault resistance, source impedance, etc.). Moreover, another
kernel function (RBF) has been used, and performances of disadvantage of neural-network-based algorithms is that the
classifiers have been evaluated based on fault classification training may not converge in some cases, as the starting point is
accuracy. In order to determine the optimal parametric settings chosen at random and can end up in a local minimum [6-7].
of an SVM classifier (such as the type of kernel function, its This paper presents a new approach for fault detection and
associated parameter, and the regularization parameter c), five-
classification on transmission lines using SVMs for training
fold cross-validation has been applied to the training set. It is
with an online wavelet-based preprocessing-stage. The
observed that an SVM with an RBF kernel provides better fault
classification accuracy than that of an SVM with polynomial
algorithm is based mainly on calculating the RMS value of
kernel. It has been found that the proposed scheme is very fast transient energy of pre-fault and post-fault signals of the three
and accurate and it proved to be a robust classifier for digital line currents and three line voltages. These preprocessed
distance protection. signals are trained with SVM. Fault classification accuracy is
Keywords—Discrete wavelet transforms, Support vector estimated for different loading and for source impedance
machine (SVMs), Feature extraction, Fault classification accuracy
and Ground index.
Table I
Transmission line Parameters
Positive sequence resistance R1 Ω/KM 0.02056
Table II
Source Parameters
Parameters Source 1 Source 2
Positive sequence 5∟85 10∟85
impendence
16
sufficient condition for wavelets is that it must be oscillatory, data. When formulated under SVMs the optimal hyper-plane
must decay quick1y to zero and must have an average value of is the decision boundary that attains the maximum margin of
zero. In addition, for the discrete Wavelet transform separation between the two classes. Optimal hyper-plane
considered here, the Wavelets are orthogonal to each other under these conditions can be uniquely determined by solving
[11]. Wavelets can provide multiple resolutions in both time a constrained optimization problem whose solution has an
and frequency. The signals can be accurately reproduced with expansion
the Wavelet analysis using relatively small number of w = ∑iai xi
components. The analyzing wavelets are called the “mother
wavelets” and it’s dilated and translated versions are called the which is the weighted sum of the training patterns that lie on
“daughter wavelets”. It has a digitally implementable the margin [24]. The training vectors which lie on the margin
counterpart called the discrete Wavelet transform (DWT). The are referred to as the support vectors. Formally, given a set of
generated waveforms are analyzed with Wavelet multi-
features vectors xi and corresponding class labels yi we form
resolution analysis (MRA) to extract sub-band information
from the simulated transients. Daubechies Four (db-4) wavelet the index pair {( xi , yi )} where xi ∈ R and yi ∈ {+1,−1} .
is used in this work for the analysis as it closely matches the
signal to be processed which is of utmost importance in We denote those input features xi with label +1 as belonging
wavelet applications. Wavelet co-efficient of the signal are to Class 1 and those feature vectors with label -1 as belonging
obtained by the decomposition of a discrete fault current and to Class 2. In the linearly separable case we can find a hyper-
voltage signals using Mallat’s algorithm [12]. plane which is able to separate the training data. This
Due to the unique feature of providing multiple separating hyper-plane will take on the following form,
resolutions in both time and frequency by wavelets, the sub- f ( x) = wx + b
band information can be extracted from the original signal.
When applied to faults, this sub-band information of faulted Where wi ∈ R and is the normal vector to the separating
power system is seen to provide useful signatures of faults, so hyper-plane and b ∈ R is the scalar bias. Hence, if x belongs
that fault classification can be done elegantly. By randomly to Class 1(+1 label) then f ( x ) ≥ 0 and conversely if x
shifting the point of fault on the transmission line, a number of
simulations are carried out employing PSCAD/ EMTDC. The belongs to Class 2 (-1 label) we have f ( x ) ≤ 0 .
generated time domain signal for each case is analyzed using Consequently, we have the following constraints,
Wavelet transform. From the different decomposed levels [13- yi ( wi xi + b) ≥ 0 for i= 1…n
17]. Only first level output is considered for the analysis. As
for a weight vector w and bias b we define the margin
the absolute values of the summation of first level output for
separation ρ as the separation between the hyper-plane and
all the inception angles considered in the analysis are found
higher as compared to that of other level outputs. It indicates the closest feature vector. In particular an SVM finds the
that the total area under the characteristics of first level outputs separating hyper-plane which maximizes the margin of
is more than that of other level outputs. Another reason for the separation ρ . The idea behind maximizing the margin of
first level output to be selected as the parameter for fault separation is rooted in concepts of structural risk
Classification is due to discrimination in transient energy minimization.
content of faulty and healthy signals is higher as compare to Roughly, when the margin of separation is maximized the
other levels. Moreover, as the sampling time considered in the generalization error on unseen test data tends to its minimum
analysis is 0.0002 s which corresponds to a sampling value [21-25].
frequency of 5 kHz and the total number of wavelet levels The SVM classifier is also defined as
considered is 6, therefore, 6th level wavelet output corresponds f ( x ) = wΦ ( x ) + b (1)
to a frequency band of 70 – 160 Hz. Level 1 corresponds to
frequency band of 1.25 – 5 kHz, which are rich in higher order where Φ ( x ) is a mapping function to map the input pattern
harmonics, predominant in case of faults. The RMS value of x into higher dimensional space . This classifier f ( x) is
transient energy of current and voltage signals is utilized as a linear in terms of the transformed data, but non-linear in terms
feature for detection and discrimination among the faults. of the original data
However, the fault is to be discriminated from other types of Following non-linear transformation, for the parameters of
transient disturbances before using the classification the decision function f ( x ) the minimization of following cost
algorithm.
function is
IV. SUPPORT VECTOR MACHINE 1 2
j ( w, b) = w + c ∑i ξ i
Support vector machines (SVMs) were first introduced 2
by Vapnik [19-20] in the late 1960’s as a part of his Subjected to constraints
development of statistical learning theory. Support vector
machines are a type of hyper-plane classifier which attempts
yi ( w.φ ( xi ) + b ≥ 1 − ξ For i = 1…n (2)
to find an optimal separating hyper-plane for the input training The solution to minimize the above cost function subject to the
constraints in Eq. (2) is given from following dual:
17
Maximize and voltage are exported to support vector machine model for
1 training. The output of each SVM indicates weather a particular
Ld = ∑ α i − ∑ α iα j yi y jφ ( xi )φ ( x j ) phase is involved in fault or not [34].
i 2 i, j
Ground fault detection is carried through detection of
Zero-Sequence component. Phases which do not involve the
1
Ld = ∑ α i − ∑ α i α j y i y j k ( xi , x j ) (3)
ground has almost zero, zero-sequence component, where as
2 i, j for phases involving ground have observable amount of zero-
i
sequence component.
Where k ( xi , x j ) known as kernel, is a non-linear
function and is defined as
k ( xi , x j ) = φ ( xi , x j ) (4)
Subjected to constraint
0 ≤ αi ≤ c
∑α y
i
i i =0 (5)
Solution is given as
Ns
w = ∑ α i yiφ ( xi )
i
k ( x, y ) = e
(−γ x− y ) 2
Where γ =
1
(6) A. Data Generation
2σ 2 For training and testing of SVM, the fault simulation has
been carried out on transmission system as shown in Fig.1.
V. PROPOSED FAULT CLASSIFICATION SCHENE PSCAD software is used to simulate the faults on the
transmission line. Table III shows the different conditions that
A 400kv, 50 Hz power system shown in Fig.1, consisting of
are simulated on transmission line for data generation used for
two sources representing two areas connected by 128 km long
training the SVMs and Table IV shows condition
transmission line is used for simulation studies. Faults take
corresponding to the data simulated for testing the trained
place at different position of transmission line. All the data is
models.
collected from sending end of the system. Simulations of faults
which occur at 0.6 sec are shown in Fig. 2. . Fault signals Table III
(voltages and currents) are collected for a cycle at 50 Hz. These Line and system parameters used for generation of training patterns
signals are taken in discrete form, 100 samples of current and Fault type a-g, b-g, c-g, a-b, b-c, c-a, a-b-g, b-c-g, c-a-g, a-
b-c and a-b-c-g
100 samples of voltage are taken for all three phases for all 11 Fault location (KM) 25.6, 51.2, 76.8 and 102.4 (128 KM line length)
types of faults [25]. Fault resistance (Ω) 0 , 25 , 50 , 75 , 100 , 125 and 150
The proposed scheme makes use of SVM classifier for Fault inception angle 0, 45, 90, 135, 180, 225, 270, 315 and 360
(0’)
detection and classification of fault. Fig. 3 shows the flow chart Load angle variation 10% and 20% ( source 1)
of proposed algorithm. The fault classifier has three classifiers (0’)
for each phase. The extracted features are fed to these Table IV
classifiers and their output indicates the type of fault. The Line and system parameters used for generation of testing patterns
voltages and currents signal are sampled at a frequency of 5 Fault type a-g, b-g, c-g, a-b, b-c, c-a, a-b-g, b-c-g, c-a-g, a-b-c
kHz. The RMS value of transient energy of voltages and and a- b-c-g
Fault location (KM) 19.2 and 57.6
currents are extracted using db4 as mother wavelet. The Fault resistance ( ) 5 Ω, 30 Ω, 55 Ω, 80 Ω, 105 Ω, 135 Ω and155Ω
decomposition is performed at level 1, which contains rich. Fault inception angle
harmonics of order 25-50 Hz. From the wavelet decomposition (0’) 0, 60, 120, 180, 240, 300 AND 360
of simulated signals (voltage and current), it is clearly
Load angle variation 5% and 15% ( source 1)
detectable that RMS values of transient energy of current (0’)
signals for faulty phases is much more as compare to healthy
phases. Then these RMS value of transient energy of current
18
B. Selection of mother wavelet for feature Extractions
The Daubechies family is one of the most suitable wavelet
families in analyzing power-system transients although there
are no definite criteria for the selection of wavelets; the best
choice is a wavelet that most strikingly exhibits the
phenomena to be studied. It has been found that db4, with its
characteristic as the shortest wavelet in its family, is an
excellent choice to detect abrupt change in RMS value of
transient energy of the pre-fault and post-fault signals. Even
though other families of wavelet-transform also give change in
transient energy of the pre-fault and post-fault signals, but this
change is more in case of db4 mother wavelet. Hence db4 is Fig.5. Wavelet coefficients of phase “a” for a-g
used as mother wavelet. The Fig.4.shows the variation of
wavelet coefficients decomposed at level 1 using db4 as
mother wavelet. On the top of Fig. 4 comparison of phase ‘a”
fault current for single line to ground fault vs. normal load
current is shown, middle corresponds to phase ‘b’ current and
bottom to phase ‘c’ current for the same type of fault.
Fig.6.Transient energy of current and voltage for phase “a” for a-g fault
19
Table V
Selection of parameters for RBF kernel function
C Accuracy variation with various value of sigma (ߪ)
20
Table VII
Fault classification result for different training samples and testing samples for different fault location, including source variation
Line Training Testing TP Ph TN Ph TP Ph TN Ph TP Ph TN Ph Accuracy Ph Accuracy Ph Accuracy Ph
length samples samples A A B B C C A (%) B (%) C (%)
P 8008 1313 1078 253 1202 111 1118 195 82.1021 91.5461 85.1485
21
[20] C. Cortesand and V. Vapnik, “Support vector network,” [39] Thomas Dalstein and Bernd Kulicke, “Neural Network
Machine Learning, Vol. 20, pp. 273-297, 1995. Approach to Fault Classification for High Speed Protective
[21] C. J. C. Burges, “A tutorial on support vector machine for Relaying,” IEEE Transactions on Power Delivery, Vol.10, No.
pattern recognition,” Data Mining and knowledge Discovery, 2, pp. 1002-1011, April 1995.
Vol.2, No.2, pp. 121-167, 1998. [40] C. H. Kim and R. Aggarwal, “Wavelet transform in power
[22] Haykin Simon “Neural networks: A comprehensive system,” IET Power Engineering Journal, Vol.15, pp. 193-200,
foundation”, 2nd Edition, Prentice-Hall, Inc., 1999. August 2001.
[23] Chih-Wei Hsu and Chih-Jen Lin, A, “Comparison of methods [41] B. Ravikumar, H. P. Khincha, D. Thukaram, “An intelligent
for multiclass support vector machines,” IEEE Transactions on approach using support vector machines for monitoring and
Neural Networks, Vol.13, No.2, March 2002. identification of faults on transmission systems,” IEEE Power
[24] Chih-Wei Hsu and Chih-Jen Lin, “LIBSVM: A library for India Conference, pp. 184-190, 05, June, 2006.
support vector machines,” 2001 [42] P. K. Dash, A. K. Pardhan, and G. Panda, “Application of
[25] Joe-Air Jiang, ching-Shah chen, and Chih-Wen Liu, “A new minimal radial basis function neural networks to distance
protection scheme for fault detection, direction discrimination, protection,” IEEE Transactions on Power Delivery, Vol.16, No.
classification, and location in transmission lines,” IEEE 1, pp. 68-74, January 2001.
Transactions on Power Delivery, Vol.18, No.1, pp. 1253-1261, [43] R. N. Mahanty, P. B. and Dutta, “Application of radial basis
January 2003. function neural networks to fault classification and location in
[26] U. Parikh, R. P. Maheshwari, B. Das and B. Bhalja, “A new transmission lines,” IEE Proc. Gener. Transm. Distrib., Vol.151,
fault classification technique for protection of series No. 2, March 2004.
compensated transmission lines,” International Conference on [44] Jeo-Air Jiang, Ching-shan Chen, and chin-Wen Liu, “A new
Power System Protection, CPRI, Bangalore, India, 20–21 protection technique for fault detection, direction
February 2007. discrimination, classification, and location in transmission
[27] R. K. Aggarwal, Q. Y. Xuan, A. T. Johns, R. W. Dunn and lines,” IEEE Transactions on Power Delivery, Vol.18, No. 1,
A.Bennett, “A novel fault classification technique for double- pp. 34-42, January 2003.
circuit lines based on a combined unsupervised/supervised [45] A. S. Omar Youssef, “New algorithm to phase selection based
neural network,” IEEE Transactions on Power Delivery, Vol. on wavelet transform,” IEEE Transactions on Power Delivery,
14, No. 4, pp. 1250–1256, October 1999. Vol.17, No.4, pp. 908-914, October 2002.
[28] Bhavesh, Bhalja and R. P. Maheshwari, “Trends in adaptive
distance protection of multi – terminal and double-circuit lines,”
Int. J. Elect. Power Components Syst., Vol. 34, pp. 603–617,
June 2006.
[29] Bhavesh, Bhalja and R. P. Maheshwari, “Protection Of Double
Circuit Line Using Wavelet Transform,” IE Journal-EL, Vol.
87, pp. 67-70 September 2006.
[30] Bhavesh, Bhalja and R. P. Maheshwari, “Philosophy of
protection for multi-terminal and double-circuit lines,” Int.
Conf. on Computer Applications in Electrical Engineering
Recent Advances, Roorkee, India, Sept. 28-October 1, 2005, pp.
582-588.
[31] Yi, Hu, Damir, Novosel, Murari Mohan Saha, and Volker
Leitloff, “An Adaptive Scheme For Parallel Line Distance
Protection,” IEEE Transactions on Power Delivery, Vol. 17, No.
1, pp. 105-110 January 2002.
[32] Xiangning Lin, Hanil Weng, and Bin Wang, “Identification of
cross country fault of power transformer for fast unblocking of
differential protection,” IEEE Transactions on Delivery, Vol.
24, No. 3, pp. 1079-1086 July 2009.
[33] Bhavesh, Bhalja and R. P. Maheshwari, “New differential
protection scheme for tapped transmission line,” IET Gener.
Transm. Distr., Vol. 2, pp. 271–279, 2008
[34] Bhavesh, Bhalja and R. P. Maheshwari, “Wavelet-based fault
classification scheme for a transmission line using a support
vector machine,” Electric Power Components and Systems,
Vol.36, Pp. 1017-1030, 2008.
[35] V. S. Kale, S. R. Bhide, R. P. Bedekar and G. V. K. Mohan,
“Detection and classification of faults on parallel transmission
lines using wavelet transform and neural network,” International
Journal of Electric and Electronics Engineering Vol. 1, pp.364-
368, 2008.
[36] P. S. Bhowmik, P. Purkait, K. Bhattacharya, “A novel wavelet
transform aided neural networks based transmission line fault
analysis method,” Electrical Power & Energy System, Vol.31,
No. 5, pp. 213-219, 2009.
[37] Robert Salat and Stanislaw Osowski, “Accurate fault location in
the power transmission line using support vector machine
approach,” IEEE Transaction on Power Sysyem, Vol.19, No.2,
May 2004.
[38] S.M. De Souza, A.P. Alves da Suva, and A.C.S. Lima,
“Attribute selection for fault location recognition in
transmission lines,” IEEE Power Engineering Society General
Meeting, pp. 6, October, 2006.
22