 Open Access
 Authors : Sunil V Chandran , Sarath J , Jithinkumar I J
 Paper ID : IJERTV10IS050332
 Volume & Issue : Volume 10, Issue 05 (May 2021)
 Published (First Online): 31052021
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Fault Port Method for Fault Analysis
Sunil V Chandran Assistant Professor, EEEDepartment
GEC Palakkad Kerala, India
Sarath J Assistant Professor, EEEDepartment
GEC Palakkad Kerala, India
Jithin Kumar I J Assistant Professor, EEEDepartment
GEC Palakkad Kerala, India
Abstract Faults are the basic hurdle for any power system network. As the complexities of the system increased, the probability of the occurrence of the fault also increased. Hence it is always better to take precautions which involve study of system under various fault conditions. This paper describes the generalized method for solving any combination of simultaneous balanced and unbalanced faults and calculating symmetrical fault port currents. This method can be applicable irrespective of any combination of faults, any number of faults and their location.
KeywordsSimultaneous faults, fault port voltage, fault port current

INTRODUCTION
Load flow analysis requires the calculation of voltages and currents of all buses. These values assist in the control and planning of the system. Simultaneous faults are those cases in which there are one or more faults occurring simultaneously. The reason for the occurrence of such faults may be due to lightening stroke, human made accidents etc. and it is quite difficult to analyze more than one fault occurring at the same time. Faults mainly occur in sequentially rather than simultaneously. Keeping the aspect of protection in mind it is desirable to compute the voltage and current during and after the occurrence of first and second fault. For accomplishing this, system has to be monitored under all kinds of faulted situations. There are different types of faults which may be categorized as balanced and unbalanced faults. These simultaneous faults may be any combination of these balanced and unbalanced faults. Balanced faults can be of both symmetrical or unsymmetrical mode. [1].
An important method of simultaneous fault analysis
coupling is present between the lines but that is not considered in this paper.
The impedance matrix is to be modified whenever a fault occurs in a system. This is the mandatory requirement while doing fault analysis.. This modification can be anything ranging from addition or deletion of a bus as given in [8]. The sequence network interconnections are also another important point that needs consideration [7].

SEQUENCE NETWORK INTERCONNECTION The interconnection between the sequence networks
must be such that that all the limitations and the boundary conditions for that particular fault must be fulfilled. From table II it can be seen that for some faults the currents of all sequence networks at the fault point are equal and the sum of the voltages of all sequence network is zero from which we can conclude that the sequence networks are connected in series with each other. Similarly, for other types of faults voltages of all sequence networks are equal and the sum of all the sequence currents is equal to zero from which we can conclude that the sequence networks are connected in parallel.
In case of simultaneous faults it is not as simple as in case of single fault analysis. In very few combinations of faults these networks fulfill the connection constraints.
This difficulty mainly arises due to the occurrence of faults on different phases. To overcome this problem Z. X. Han has recommended using phase shifting transformers to shift the voltages and currents accordingly in order to meet the boundary conditions [4]. Ideal turns ration is given as:
comes from G. Gross and H.W. Hong in 1982 in the form of twostep compensation method [2], [3]. Papers [4], [5] and [6] deals with various other methods that could be
= 1 =
2
2 (1)
1
adopted for simultaneous fault analysis.
This paper is based on the method of multiple port theory and theory of symmetrical components given in [4]. Current and voltage constraints for various fault types with reference to symmetrical components are shown in table I [7].
Under the normal fault analysis condition analyzing occurrence of fault on any phase is not a difficult problem but in case of simultaneous faults occurrence of faults on different phases has to be considered. Usually, mutual
where I1 and V1 are the primary current and voltage respectively of ideal phase shifting transformer and I2 and V2 are secondary current and voltage of ideal phase shifting transformer.
TABLE I. CONSTRINTS BETWEEN SYMMETRICAL CURRENTS AND VOLTAGES FOR DIFFERENT FAULTS
S.No
Fault
Current Constraints
Voltage Constraints
1
Single Line to Ground Fault
I1=I2=I0 I1=a2 I2=aI0 I1=aI2=a2 I0
V1+V2+V0=0 V1+a2 V2+aV0=0 V1+aV2+a2 V0=0
2
Line to Line Fault
I1= I2, I0 = 0 I1= – a2 I2, I0 = 0
I1= – aI2, I0 = 0
V1 = V2 V1 = a2 V2 V1 = aV2
3
Double Line to Ground Fault
I1+I2+I0=0 I1+a2 I2+aI0=0 I1+aI2+a2 I0=0
V1=V2=V0 V1=a2 V2=aV0 V1=aV2=a2 V0
4
Single Phase Open Fault
I1+I2+I0=0 I1+a2 I2+aI0=0 I1+aI2+a2 I0=0
V s1=V s2=V s0
Vs1=a2Vs2=aVs0 Vs1=aVs2=a2 Vs0
5
Two Phase Open Fault
I1=I2=I0 I1=a2 I2=aI0 I1=aI2=a2 I0
Vs1+V s2+V s0=0
Vs1+a2Vs2+aVs0=0 Vs1+aVs2+a2Vs0=0
Z
C
Ic = 0
Va+IaZ = 0
Vb+IbZ = 0
Double Phase open
Z
A
Ib = Ic = 0 Va + IaZ = 0
B
Ic = Ia = 0 Vb + IbZ = 0
C
Ia = Ib = 0 Vc + IcZ = 0

BASIC NETWORK EQUATION
Basic equations for any n bus system may be given as [9]:
[] = [] (2)where Ybus is nÃ—n admittance matrix and Vbus and Ibus are nÃ—1 voltage and current vector respectively.
Vbus can be obtained by using equation (2)
Vbus = [Zbus] Ibus (3)
where [Zbus] = [Ybus]1 i.e. impedance matrix of size nÃ—n. In matrix form we can write the above equation as follows:
[1] = [ 1 12] [1 ] (4)2
21 22
2
TABLE II. BOUNDARY CONDITIONS FOR DIFFERENT TYPES
OF FAULTS
Type of fault
Specific phase
Boundary condition
Single line to ground fault
Zd
A
Ib = Ic = 0
Va+IaZd=0
B
Ic = Ia = 0
Vb+IbZd=0
C
Ia = Ib = 0
Vc+IcZd=0
Double line to ground fault
Z Z
Zd
A
Vb+Ib Z=Vc+Ic Z = (Ib+Ic)Zd Ia = 0
B
Vc+Ic Z=Va+Ia Z = (Ic+Ia)Zd Ib = 0
C
Va+Ia Z=Vb+Ib Z = (Ia+Ib)Zd Ic = 0
Single phase open
A
Ia = 0
Vb+IbZ = 0
Vc+IcZ = 0
B
Ib = 0
Vc+IcZ = 0
Va+IaZ = 0
where Z1 and Z2 are known as drivig point impedance and Z12 and Z21 are known as transfer impedance [4].
Hence, when fault occurs anywhere on the system impedance matrix changes accordingly. This paper uses fault port matrix to calculate parameters related to simultaneous faults as given in [4]. According to the multiple port theory each fault point is considered as a port. In a port the current coming out of one terminal should be equal to the current entering another terminal [1].

FAULT PORT MATRIX
In [2] although the short circuit fault was occurring but during simulating it the impedance matrix i.e. Zbus was used unchanged. But usually we cannot consider it like that as the fault may occur at any point on the network. Always considering a fault point as an addition of bus may make the calculations difficult. Hence Z. X. Han has devised the method of fault port matrix which is easy to understand and lesser the computation time as well.

Fault port matrix elements
Taking case of a general short circuit fault at any point k on any line as shown in fig. 2 and assuming the injected
current to be Ik and voltage at that point be Vk the equations may be:
= (5)
= (1 ) (6)
= (7) where
r is the distance between point i and k.
Fig. 1. Short Circuit fault at point k
In the formation of fault port matrix it is assumed that the unity current is injected at any bus m and the voltage at
Further these equations are manipulated using the method given in [9]. Considering the line ij is switched open and two sample buses s and t have been inserted in the circuit.
Method used in [9] demonstrates the modification of Zbus by addition and deletion of the bus. So here opening of line st will act as opening of the branch.
Transfer impedance between the fault point k and any bus m after opening of the branch is given as:
= ( 1 + )1 (13) driving point impedance is given as:
= ( 1 + )21 (14) transfer impedance between fault point k and terminal t is given as:
= ( 1 + )1( 1) (15) where
port k is equal to the impedance between m and k.
the ter s C
, C , C and P can be described as:
Therefore,
m m i k
= ,
x=m,k,i (16)
= = (1 ) + (8) where
=
(17)
Zmk is the transfer impedance.
For finding the driving point impedance the equation is:
= = (1 )2 + 2 + 2(1 ) +
(1 ) (9)
If two short circuit faults are occurring simultaneously on the same line as shown in fig 3 then the transfer impedance can be given as:
To understand these equations in a better way one may refer to reference [9].

Methods
The method given in this paper is the review of the method given by Z. X. Han in [4]. Considering N number of faults occurring out of which M are series faults and N M are the parallel faults.
The equation for these faults may be given as follows
= (1 )(1 )
+ (1 )
+ (1
[4]:) + + (1 ) (10)
() () () ()
[ ] = [ ] [ ] +(0)
] (18)
()
() ()
()
[ (0)Fig. 2. Two short circuit faults occurring simultaneously at k and k
If a short circuit fault and an open circuit fault occur simultaneously on the same line as shown in fig 4 then the driving impedance will be given as follows:
= + (1 ) (11)
S
S
(0)
(0)
where VS(i), IS(i) and V (0) are voltage, current and open circuit voltage of series type fault having size MÃ—1.
VP(i), IP(i) and VP are voltage, current and open circuit voltage of parallel type fault having size (NM)Ã—1.
ZSS(i), ZSP(i), ZPS(i) and ZPP(i) are the impedance matrices of sizes MÃ—M, MÃ—(NM), (NM)Ã—M and (NM)Ã—(NM)
respectively and i = 0,1,2 representing zero, positive and negative sequence network respectively.
Performing partial inversion on equation (18) we obtain the following equations:
transfer impedance for the same will be given as:
()
[] = [
() ()
] [
()
(0,0)
] + [
] (19)
= (1 ) (12)
()
() ()
()
(0)
Fig. 3. Short circuit fault and open circuit fault occurring simultaneously
where ASS(i) = ZSS(i) ZSP(i)ZPP(i)ZPS(i) ASP(i) = ZSP(i) ZPP(i)1
APS(i) = ZPP(i)1 ZPS(i) APP(i) = ZPP(i)1
V = V Z Z V
V = V Z Z V
(0,0) (0) 1 (0)
S S SP(1) PP(1) P
I = – Z V
I = – Z V
(0) 1 (0)
P PP(1) P
Since here simultaneous faults are considered therefore phase shifting transformers might be needed. Using equation (19)

will be equal to the voltages obtained from the load flow. Hence, it is one of the major parts of the flowchart.
The admittance matrices formed here are sparse in nature. The computation time is saved using admittance matrix in load flow analysis. Then these matrices are transformed into impedance matrices for further calculation as the fault analysis generally uses z bus values for its
[ ()] = [()
0 ()
] [
()
]
calculations..
()
0 ()
() ()
Once the type of the fault is specified we can proceed with
() 0 1
[ 0 ()]()
()
()
[
] + [
(1) (0,0)
(1) (0) ]
the formation of fault port matrix. Fault port matrix is formed based on the method given above. Bus voltages and currents of the buses affected by occurrence of the fault are
(20)
where NS(i) and NP(i) are turn ratio matrices of the ideal phase shifting transformer having sizes MÃ—1 and (NM)Ã—1 respectively.
According to the boundary conditions given in table 2 for series and parallel faults equation (20) can be further simplified as:
calculated. The value of currents and voltages are aligned with the boundary values.
Once a particular fault analysis is completed we can proceed with the other faults, if present. Else the algorithm stops and result is obtained.
Start
(0,0)
Enter the data required for load flow
0
()
(1)
analysis
[ ] = [] [
] + [
(0) ] (21)
0
where
()
(1)
Perform the base case load flow
()
()
= () ()1
=0,1,2
= ()
=0,1,2
()
1
()
Form Y bus matrices for all three sequence networks
()
()
= () ()1
Input types of faults
=0,1,2
= ()
=0,1,2
()
1
()
Form fault port matrices for all sequence network and open circuit voltage for positive sequence network is calculated
Solving equation (21) series fault port currents and parallel fault port voltages can be calculated as given below:
Solve for the boundary currents and voltages
()
() 0
1
1
(1) (0,0)
[ ] = [
] [ ] [
] (22)
()
0 ()
(1) (0)
Calculate bus voltages and currents on the selective basis
Substituting equation (22) in equation (19) series fault port voltages and parallel fault port currents can be calculated
This method has been computationally implemented and the results are shown in this paper using the algorithm.


FLOWCHART[4]
Below given is an algorithm implementing the above given method. First and foremost load flow is to be done to get the reference values to see the changes occurred after
Next fault
no
Stop
yes
faults. Also, the open circuit voltage mentioned in equation
Fig. 4. Flowchart

RESULT
Data for a four bus system was taken from [10] and the above method was implemented on that system. Giving different kinds of faults as input results are given below.
Fig. 5. Circuit Diagram
Fig. 5 shown above represents a four bus system whose line data and bus data are given below in table III and table IV.
Table V. THREE COMPONENTS OF FAULT CURRENTS
Bus No.
Positive Sequence
Negative Sequence
Zero Sequence
2
1.1672
1.1672
0.32661
4
1.1212
1.1212
0.2324
Fault port currents for a single line to ground fault on bus 2 and open circuit fault on bus 4 is shown in fig. 7:
TABLE III. LINE DATA
from bus
to bus
R
X
Y/2
R0
Xo
1
2
0.0101
0.05
0.05
0.003
0.013
1
3
0.0074
0.037
0.04
0.002
0.009
2
4
0.0074
0.037
0.04
0.002
0.009
3
4
0.0127
0.064
0.06
0.003
0.016
TABLE IV. BUS DATA
Bus
Pg
Qg
Pload
Qload
V(p.u)
Angle
type
1
0
0
50
30.99
1
0
1
2
0
0
170
105.4
1
0
2
3
0
0
200
123.9
1
0
2
4
318
0
80
49.58
1.02
0
3
For single line to ground fault occurring on bus 2 and 4 simultaneously as shown in fig. 6 fault port currents are given in table V:
Fig. 6. Single line to ground fault on bus 2 and 4 simultaneously
Fig. 7. Single line to ground fault on bus 2 and open circuit fault on bus 4
The fault port currents for the above mentioned condition is given in table VI.
TABLE VI. THREE COMPONENTS OF FAULT CURRENTS
Bus No.
Positive Sequence
Negative Sequence
Zero Sequence
2
0.064093
0.064093
0.064247
4
0.031969
0.031969
0.031969

CONCLUSION
The generalized method for doing simultaneous fault analysis have been discussed and results of calculated fault port current and voltages have been shown.
REFERENCES

P. M. Anderson, Analysis of Faulted Power Systems, The Iowa State University Press, 1973.

G.GROSS, H.W.HONG, A twostep compensation method for solving short circuit problems, IEEE Transaction on Power Apparatus and Systems, vol.PAS101, pp. 13221331, 1982.

G.GROSS, H.W.HONG, A twostep compensation method for solving short circuit problems, IEEE Transaction on Power Apparatus and Systems, vol.PAS101, pp. 13221331, 1982.

Z.X. Han, "Generalized Method of Analysis of Simultaneous Faults in Electric Power System". IEEE Transactions on Power Apparatus and System &, Vol PAS101, pp. 39333942, October 1982.

V.BRANDWAJN, W.F.TINNEY, Generalized method of fault analysis, IEEE
Transaction on Power Apparatus and Systems, vol.PAS104, pp. 13011306, 1985

R.BUALOTI, P.Pugliese, F.Torelli, M.Trovato,A Generalised Method For Simultaneous Fault Analysis, MELCON96,8th Mediterranean Electrotechnical Conference, pp 721725, 1996

J. R. Mortlock,The Evaluation of Simultaneous Faults on Three Phase Systems, Journal of Institution of Electrical Engineers
Part II : Power Engineering,, vol. 94, pp. 166190, 1947

T. E. DyLiacco and K. A. Ramarap, "Shortcircuit calculation for Multiline Switching and End Fault", IEEE Transaction on Power Apparatus and systems,vol. PAS89, No.6, pp.12261237, 1970.

M.EtezadiAmoli, Simultaneous Fault Analysis Using The Generalised Method of Fault Analysis, Proceedings of the TwentyFirst Annual Norht American Power Symposium, pp 188191, 1989

John J. Grainger and William D. Stevenson, Jr, Power System Analysis,McGrawHill Publication, Inc.