gohokar2005
gohokar2005
Received 5 April 2004; received in revised form 2 November 2004; accepted 10 January 2005
Abstract
This paper presents new techniques for locating short circuit and open conductor faults on automated distribution system. This is a unique
method where short circuit and open conductor faults are treated simultaneously. The proposed fault location scheme is capable of accurately
identifying the zone in which fault occurs. Temporary faults are also detected which may not result in a blown fuse. Data is made available
to the substation remote terminal unit from various automated distribution system equipments. Decision is taken based upon the magnitude
of fundamental frequency component of voltage and current phasor. Testing of proposed fault location method has been demonstrated on the
laboratory model for its practical use. Analogues simulation results are obtained using electromagnetic transients program (EMTP).
© 2005 Published by Elsevier B.V.
Keywords: Distribution automation; Symbolic method; Phasor estimation; Short circuit faults; Open conductor fault
where Cf is the complex Fourier coefficient. Eq. (1) can be and 12–13. Let us assume that the fault occurs in between
written as sections 11 and 12. The status is determined from each section
N−1 CT on this feeder. Polling for magnitude of current is started
1 2πm 2πm from substation end CT. Results of this polling is 1, 1, 1, 0
Cf = xm cos − j sin (2)
N N N and 0 means fault is in between buses 11 and 12.
m=0
However, the problem of fault location in the feeder hav-
or
ing laterals is more challenging. Consider the feeder-1 having
N−1
1 1 2πm 2πm sections on main feeder as 1–2, 2–5, 5–6 and sections on lat-
(Cfr − jCfm ) = xm cos − j sin (3) eral as 2–3, 3–4. Let us assume that the fault occurs in between
2 N N N
m=0 sections 2 and 3. The CT placed on bus 2 is measuring sum
Therefore of load current at buses 3, 4, 5 and 6. The status detected by
above logic is 1 and 0 for sections 2–3 and 2–5 both. Actu-
N−1
2 2πm ally the fault is in between buses 2 and 3. This problem is
Cfr = xm cos (4)
N N solved by updating current and voltage vectors using static
m=0
impedance load model [5] as follows.
N−1 The data provided by a load flow program can be used
2 2πm
Cfm = xm sin (5) to account for laterals and voltage magnitude at a load bus
N N can be used to obtain equivalent load admittance. The series
m=0
impedance for each line section is calculated from knowledge
Eqs. (4) and (5) call full-cycle discrete Fourier transform
of construction configuration and is then used to determine
where Cfr and Cfm are the real and imaginary part of am-
the voltage drop of line section. The voltage and current pha-
plitude of fundamental frequency component of a phasor.
sor at each bus are updated using the equations summarized
So the amplitude C and phase angle φ of fundamental
below. This is illustrated as follows. Referring to sections 2–3
frequency component of a phasor are given by Eqs. (6) and
and 2–5, the load current at the bus may be described as
(7):
[Iloadi ] = [Yloadi ][Vi ] (8)
C = Cfr2 + Cfm 2 (6)
where Iloadi is the load current vector at bus I, Vi the voltage
Cfm vector at bus i, and Yloadi is the matrix representing the load
φ = tan−1 (7)
Cfr admittance at bus i.
As mentioned earlier analog current and voltage signals are From the network topology starting bus, ending bus,
made available at substation RTU. Based upon the sam- impedance of line section and real and reactive loading of
ples collected from CTs and PTs the rms value of voltage the buses is known. The status of voltage read by PTs shows
and current is determined from above technique. High fre- that the lateral ending bus voltage collapsed due to fault on
quency components are filtered out from the distorted wave- that section. Main feeder voltage is normal. After this first
form. Only fundamental component is considered for further guess for fault location, the voltage drop in both the section
decision-making process. is estimated by following set of equations. The load current at
ending bus does not change while starting end bus current is
updated by adding fault current. KCL at every node is given
3.1. Short circuit faults
by Eq. (9):
Comparing the magnitude of the current signals after fault [A]T [Ib ] = [Iload ] (9)
in all the sections, immediately the zone in which fault occurs
is identified. In the feeder having no laterals, it is very easy In order to accurately calculate the load current at each bus,
to locate the fault. The comparison of magnitude of current the updated bus voltage must be known. The updated voltage
in each section shows that the fault is in that section in which vectors are calculated as follows:
starting bus CT have fault current and ending bus CT have
normal current. This particular condition is read as symbols [Vi ] = [Vline ] − [Vd ] (10)
like digits 1 and 0. This status is monitored continuously for
where Vline is the voltage measured at the line terminal, and
each section CT. Under normal conditions, when there is no
Vd is the total voltage drop in the line section up to bus i.
fault on the system, the status is normal load current read as 0.
The equation for voltage drop in any zone may be stated
After fault the sections of faulted current towards substation
as
will have high fault current read as 1 and towards far end have
normal current read as 0. So the fault is detected in that section [Vd ] = [Zb ][Ib ], [Ib ] = [K][Iload ] (11)
where this status is 1 and 0. Rest of the section has status either
1 and 1 or 0 and 0. This is illustrated as follows. Consider the where Zb is the impedance of branch, and Ib is the current
feeder-3, shown in Fig. 2, having sections 1–10, 10–11, 11–12 flow in the branch.
54 V.N. Gohokar, M.K. Khedkar / Electric Power Systems Research 75 (2005) 51–55
Acknowledgement
References
[1] R. Patterson, Signature and software find high impedance fault, IEEE
Trans. Comput. Appl. Power Syst. (1995) 12–15.
[2] A.E. Enanuel, E.M. Gulachenski, High impedance fault occurring o
sandy soil in 15 kV distribution feeders: contribution to the evalua-
tion of low frequency spectrum, IEEE Trans. Power Delivery 5 (2)
(1990) 676–686.
[3] D.C.T. Wai, X. Yibin, A novel technique for high impedance fault
identification, IEEE Trans. Power Delivery 13 (3) (1998) 738–
744.
[4] T. Takagi, Y. Yamakoshio, J. Baba, K. Uemura, T. Sakaguchi, A
new algorithm of an accurate fault location for EHV/UHV transmis-
sion lines. Part II. Laplace transform method, IEEE Trans. Power
Apparatus Syst. PAS-101 (3) (1982) 564–573.
[5] D.J. Lawarence, L.Z. Cabeze, L.T. Hochberg, Development of an ad-
vance transmission line fault location system. Part II. Algorithm de-
velopment and simulation, IEEE Trans. Power Delivery 7 (4) (1992)
1972–1981.
[6] A.A. Girgis, M.B. Johns, A hybrid expert system for faulted sec-
tion identification, fault type classification and selection of fault lo-
cation algorithm, IEEE Trans. Power Delivery 4 (2) (1989) 978–
985.
[7] A.A. Girgis, C.M. Fallon, D.L. Lubkeman, A fault location technique
for rural distribution feeders, IEEE Trans. Ind. Appl. 29 (6) (1993)
Fig. 4. EMTP simulation results. 1170–1175.
[8] M.K. Khedkar, V.N. Gohokar, An integrated approach for automation
of distribution system, in: Proceedings of the 2002 IEEE Power
current flowing through zone towards extreme end of feeder
Engineering Society Transmission and Distribution Conference, vol.
is zero while rest of the zone has normal load current. 3, 2002, pp. 2106–2110.
[9] J.-C. Gu, S.-l. Yu, Removal of dc offset in current and voltage signals
using a novel filter algorithm, IEEE Trans. Power Delivery 15 (1)
5. Conclusion (2000) 73–79.
[10] K. Srinivasan, A.St. Jacques, A new fault location algorithm for
When a fault occurs in distribution network power supply radial transmission line with loads, IEEE Trans. Power Delivery 4
(3) (1989) 1676–1682.
to the whole feeder is interrupted. The proposed technique,