Errors analysis in distance relay readings with presence of facts devices

Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

Errors Analysis in Distance Relay Readings with Presence of
FACTS Devices
Gaber El-Saady1*Rashad M.Kamel2 Essam M.Ali3
Electrical Engineering Department, Assiut University, Assiut, Egypt, gaber1@yahoo.com
Abstract
This paper presents studying the performance of distance impedance relay with the presence of FACTS
devices .Also the measured impedance at the relaying point in the presence of series FACTS device SSSC ,
and shunt FACTS device called STATCOM are obtained. A detailed model of both SSSC and STATCOM is
introduced and then the faulty system is studied analytically , where the errors in the measured impedance of
distance impedance relay are introduced as a result of the presence of series and shunt FACTS devices. The
analysis results show the effect of impacting the FACTS devices location, the values of inserted SSSC and
STATCOM voltages (operational conditions) and also the fault resistance values.
Keywords: Distance Relay, FACTS, STATCOM, SSSC.
1. Introduction
The measured impedance at the relaying point is the basis of distance protection operation. In addition to the power
system parameters, the fault resistance could greatly influence the measured impedance in the case of single phase to
ground, double phase to ground and three phase to ground faults, in such a way that for zero fault resistance the
power system parameters do not affect the measured impedance (Kazemi & Jamali 2009). FACTS devices can
be broadly classified into three types, (a) Shunt (b) Series and (c) composite series and shunt (Tarlochan et al. 2005).
The use of power electronics devices to improve power transfer capability of long transmission lines forms
the basis of the concept of FACTS. The Static Synchronous Series compensator (SSSC) is a series device of the
FACTS family using power electronics to control power flow and improve power oscillation damping on power
grids (Khederzadeh 2009).
More than 70% of transmission line faults are single phase to earth faults, and phase to phase faults are the next
common fault type . Double phase to earth and three phase faults are less common in power transmission
systems (Jamali & Shateri,2011). The operating principle of distance protection is based on the fact that , from
any measuring point in power system the line impedance to a fault in that system can be determined by measuring
the voltage and current at the measuring point (Khederzadeh 2009).
Some researchers have been done to evaluate the performance of a distance relay for transmission line in the
presence of FACTS controllers. The work in (Kazemi & Jamali 2009) presented a steady state analysis on
impacting STATCOM and SSSC and compared between SSSC and STATCOM impacting effects on the measured
impedance by distance relay. The work in (Khederzadeh 2009) presented analytical and simulation results of the
application of distance relays for the protection of transmission employing Static Synchronous Series Compensator
(SSSC).
This paper studied the effects of the installation of SSSC and STATCOM on a doubly fed circuit transmission line
in the case of single line to ground fault by means of presenting the errors in measured impedance at the relaying
point and the distance relay ideal tripping impedance values. The measured impedance is presented for two cases of
SSSC and STATCOM exclusion and inclusion in the fault loop. For three installation points of SSSC, i.e. at near
end, mid-point, and far end, the ideal tripping characteristic is investigated. Also the operational conditions of
FACTS that may be affected by loading condition are considered in this work.

2-FACTS devices models
The Facts devices under study have different types; the model of each type is described in the following sections.
2.1. SSSC model
The SSSC is placed in a group of series connected FACTS devices. As shown in figure (1) SSSC consists of a
voltage source connected in series with the transmission line through voltage source inverter (VSI) and coupling
transformer.

112

Vol.4, No.14, 2013

www.iiste.org

Figure (1) (a) Basic construction

(b): an equivalent circuit of SSSC

2.2. STATCOM model
A STATCOM is a controlled reactive power source and provides voltage support, which is its primary duty, by
generating or absorbing reactive power at the point of common coupling without the need of large external
reactors or capacitor banks (Prof. Detlef 2007).
The STATCOM basically consists of a step down transformer Tr, reactance LF, a three phase GTO or IGBT based
multi- pulse VSC and a DC capacitor as shown in figure (2a). STATCOM can be modelled as a shunt branch
consisting of an impedance Zsh due to the coupling transformer, and a voltage source Esh which is in phase with the
voltage of its connection point (Kazemi & Jamali 2009).

Figure (2) (a) STATCOM model construction

(b) STATCOM representation

2.3. The system model
The system model under studying here is a single circuit transmission line fed doubly from two generators EA and
EB where EB=EAhe-jδ. Our FACTS devices are inserted between bus A and bus B to control the transferred power,
voltage and current of line, etc. The distance relay is placed at bus A, so its function is protecting the transmission
line till point B see figure (3).

Figure (3): the model of the system consists of doubly fed transmission line
3. Fault analysis
Operating times of protective relays are typically 1.0–3 cycles and the tripping times of breakers are 1.5 – 3
cycles, giving a total fault tripping time of 2.5– 6 cycles (50 –120 ms at 50 Hz, 42–100 ms at 60 Hz) . The relays
are installed on each feeder at a substation, to which the secondary voltages and currents of the VT (voltage
transformers) and CT (current transformers) of the associated feeder line are applied as shown in Fig4. The signal
voltage V (t) and current I (t) from the secondary terminals of the associated PT and CT are applied to the relay, and
the new quantities u (t), ν’ (t) compose the detected part of the relay circuit as shown in the figure (4b). Then the
phase angular difference between u (t) and ν’ (t) is continuously compared in the phase comparator part of the relay
circuit (Xia et al.1994) and (Yi Hu et al. 2002).
In the absence of FACTS devices with zero fault resistance, the measured impedance by a distance relay only
depends on the length of the line section between the fault location and the relaying points. In Figure (3) this

113

Vol.4, No.14, 2013

www.iiste.org

impedance is equal to pZL1, where p is per unit length of the line section between the fault and the relaying points,
and ZL1 is the line positive sequence impedance in ohms. In the case of a non-zero fault resistance, The measured
impedance is not equal to the impedance of the line section between the relaying and the fault points. In this
case, the structural and operational conditions of the power system affect the measured impedance (Kazemi &
Jamali 2009).
The operational conditions prior to the fault instance can be represented by the load angle of the line δ and the ratio
of the voltage magnitude at the line ends h or EB / EA = h e-jδ. The power system before the fault at point F can be
drawn as a symmetrical circuit in the 1–2–0 domain as shown in figure (5). The fault calculations for impedance
readings will derive now from the first principles case of SLG fault and it is given in appendix A .

Figure (4): (a) principle of fault detection

(b) MHO relay characteristics

Figure (5): the 1-2-0 domain of the model case of fault with no FACTS connection.
3.1. Fault impedance calculation case of no FACTS
See appendix A for more details.
From the equations given in the appendix A the predicted measured impedances with changing fault location can be
given by the following equations:
ZA1 = ZSA1 + p ZL1

(1)

ZB1 = ZSB1 + (1− p) ZL1

(2)

ZA0 = ZSA0 + p ZL0

(3)

ZB0 = ZSB0 + (1-p) ZL0

(4)

Z1=Z2
each positive sequence impedance equal negative sequence one. Where ZA1, ZA2, ZA0 are positive, negative
and zero sequence impedances looking into the left-hand side at point F.
ZB1, ZB2, ZB0 are positive, negative and zero sequence impedances looking into the right-hand side at point F.
Now the two measured quantity Vrelay and Irelay at point A
The relaying current is as follows:
I relay = IA + KOL IA0 = (Cd +2 C1+C0 (1 + 3K0L)) If (1)
The relaying voltage

114

(5)

Vol.4, No.14, 2013

www.iiste.org

Vrelay = VA = VF + IAf (p ZL1) = 3 If (1) RF + Irelay (p ZL1)

(6)

From equations (5) and (6) Where Zrelay = Vrelay / Irelay, we found that:
Zrelay = pZL1 +

(7)

C1 = ZB1 / (ZA1 + ZB1)

(8)

C0 = ZB0 / (ZA0 + ZB0)

(9)

Where

C1=C2

because

Z1=Z2
Ze = ZA1+ZB1

(10)

K0L = (ZL0-ZL1) / (3ZL1)

(11)

Kd =1- h e-jᵟ

(12)
-j

D = ZB1+ZA1h e ᵟ

(13)

Cd = (Ze +3 RF) Kd / D

(14)

It can be seen when the fault resistance is equal to zero, the measured impedance at the relaying point is equal to
the impedance of the line section between the relaying and the fault points, and otherwise the measured
impedance deviates from its actual value depending on the power system conditions. Note K0L is Zero sequence
compensated factor.
3.2. Relay impedance with presence of SSSC
Once a FACTS devise is installed on the transmission line, the measured impedance at the relaying point is affected.
The device is installed at the length of m per unit from the relaying point. In the following, the measured impedance
in the presence of SSSC and STATCOM is presented. Inserting SSSC depending on its exclusion or inclusion in the
fault loop, the measured impedance would change. And the following equations are derived with the same previous
sequence. Firstly, some factors are defined:
CSC1 = Zse / ZL1

(15)

CSC0 = Zse / ZL0

(16)

ZAm1 = ZSA1 + mZL1

(17)

ZAm0 = ZSA0 + mZL0

(18)

Inserting SSSC may be at any point on the line and also the fault, so the measured impedance will be affected by
both the FACTS and fault location. Let us check impedance calculations at two locations (1) nearer to far end
than fault point (i.e. SSSC out of fault loop) (2) nearer to near end than fault point (i.e. SSSC inside the fault loop).
3.2.1. SSSC nearer to far end (out of fault loop)
Some equations are modified according to variation of VSSSC and Zse. Fault distance p is smaller than inserting
location m (m > p).
ZBm1 = ZSB1 + (1 – m + CSC1) ZL1

(19)

ZBm0 = ZSB0 + (1 - m + CSC0) ZL0

(20)

ZAf1 = ZSA1 + pZL1

(21)

ZmF1 = |m – p| ZL1

(22)

ZmF0 = |m – p| ZL0

(23)

ZA1 = ZAf1

(24)

ZB1 = Zmf1 + ZBm1

(25)

ZAf0 = ZSA0 + pZL0

(26)

115

Vol.4, No.14, 2013

www.iiste.org

ZAm0 = (m - p) ZL0

(27)

ZA0 = ZAf0

(28)

ZB0 = Zmf0 + ZBm0

(29)

Dn = ZBn1+ZA1 [h e-jᵟ - re jϒ]

(30)

jᵧ

Kdn = 1+ r e - h e-jδ

(31)

Using the same equation (10) of measured impedance here, the measured impedance in the case of zero fault
resistance is equal to the impedance of the line section between the relaying and the fault points as will be shown
later.
3.2.2. SSSC at near end (inside of fault loop)
In this case the fault location p is larger than FACTS inserting point m (m < p). The equations in previous case will
be modified as follow:
Zmf1 = (p – m + CSC1) ZL1

(32)

ZA1 = ZAm1+ Zmf1

(33)

ZB1 = ZBf1

(34)

Zmf0 = (p – m + CSC0) ZL0

(35)

ZA0 = ZAm0 + Zmf0

(36)

ZB0 = ZBf0

(37)
jᵧ

jᵧ

Dn = [ZAm1 (1 + re ) + Zmf1] h e-jδ + ZBf1 (1 + re )

(38)

CZ = (CSC0 - CSC1) C0 (1 + 3KL0) ZL1

(39)

jᵧ

Cv = - (Ze + 3RF) (ZAn1+ZB1) re / Dn

(40)

Zrelay = (p + CSC1) ZL1 +

(41)

It can be seen obviously that when (RF = 0) no fault resistance , the measured impedance will not be the same
as pZL1 because of the effect FACTS presence.
3.3. Relay impedance with presence of STATCOM
Inserting STATCOM in the line depending on its position with respect to fault location (inside or outside fault
location) affect the measured impedance, so some equations are modified with STATCOM with STASTCOM
voltage = Esh and an impedance ZSh as follows :
ZAm1 = ZSA1 + mZL1

(42)

ZBm1 = ZSB1 + (1 - p) ZL1

(43)

ZAm0 = ZSA0 + m ZL0

(44)

ZBm0 = ZSB0 + (1 - p) ZL0

(45)

Zmf1 = |m - p| ZL1

(46)

Zmf0 = |m - p| ZL0

(47)

These equations specify the behaviour of no fault impedances independently of inserting FACTS inside or outside
fault loop.
3.3.1. STATCOM nearer to far end (out of fault loop)
The following equations are modified also according to insertion point location.
ZA1 = ZAf1

(48)

ZB1 = Zmf1 +

(49)

ZA0 = ZAf0

(50)

116

Vol.4, No.14, 2013

www.iiste.org

ZB0 = Zmf0 +

(51)

Dn =

(52)

Kdn = 1- esh (1 – CB1) – CB1 h e

-jδ

(53)

And using equation (10) it can be observed that in case of zero fault resistance, the relay measured impedance
is equal to pZL1. Otherwise the relay impedance will have impedance different from pZL1.
3.3.2. STATCOM nearer to near end (inside fault loop)
Some modifications are performed according to point m and p locations with respect to each other.
Dn =ZAn1 + ZB1 – (CA1 + eSh (1 – CA1) – h e -jδ) (Cm ZAm1 + ZAn1)
Kdn= 1 – hn e

(54)

-jδ

(55)

Cdn = (Ze + 3 RF) Kdn / Dn

(56)

Csh = Zmf1 (- CV + (1-CA1) (Cdn + 2C1) + C0 (1- CA0) (1 + 3K0L))

(57)

Zrelay = pZL1 +

(58)

It is obvious also that when RF = 0, the relay impedance is not equal to pZL1 as will be discussed later.
4. Case Study and Results
We choose in our study three important points of impacting of FACTS through transmission lines (1) at near end (2)
at far end (3) at midpoint impacting. Power system calculations gave these impedance values of the system shown in
fig .3.
ZL1 = 0.01133 + j 0.3037

Ώ/KM

ZL0 = 0.1535 + j 1.1478

Ώ/KM

ZSA1 = 1.3945 + j 15.9391

Ώ

ZSA0 = 7.454 + j 27.8187

Ώ

ZSB1 = 0.6972 + 7.9696

Ώ

ZSB0 = 3.7270 + j 13.9093

Ώ

h= 0.96
δ=16˚
The ideal tripping impedance readings without FACTS are shown in figure (6a) and the errors in readings due to
presence of fault resistance as in figure (6b).
tripping characteristics

the reading errors in resistance

100

350
Rf = 0
Rf = 100
Rf= 200

300
R errors in ohm

90
80
70

250
200
150
100
50

60
0

0.1

0.2

0.3

0.4

0.5
percentage distance

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

the reading errors in reactance
10

40

5

30
X errors in ohm

X relay

0

50

20
10

0
-5
-10
-15

0

0

50

100

150
200
R relay

250

300

350

-20

0

0.1

0.2

Figure (6) (a) tripping characteristics

0.3

0.4

0.5
percentage distance

(b) errors due to RF presence

117

Vol.4, No.14, 2013

www.iiste.org

From the above figure it is observed that when a bolt fault occurred, the angle of measured impedances were
constant and the magnitudes only changed or say that both the measured resistance and reactance values
changed with the same percentage with changing the fault location. Unlike the previous case, case of fault
occurrence through resistance Rf≠0 both the angles and magnitudes changed. As shown in figure 6 the resistance
increased but the reactance decreases compared with case of bolt fault.
4.1. FACTS at far end (outside fault loop)
4.1.1. SSSC impacting
From equation (10) , there is no error at any value of VSSSC when a bolt fault occurred at any value of SSSC
injected voltage. The results in fig.7 shows the errors occurred in impedance (resistance and reactance)
readings case of fault occurrence through resistance Rf = 50 ohm . With changing the fault location, the
errors in resistance were negative what means that resistance decreased compared with original case in fig6,
also the reactance decreased with small changes, but at the end of the line fault occurrence the error was larger.
With increasing the values of injected voltage we observe that readings error increased negatively in both resistance
and reactance.
errors in Rrelay with Rf = 50.0 ohm
20

R errors in ohm

0

-20
r=0.0
r= 0.1
r = 0.2

-40

-60

-80

0

10

20

30

40

50
percentage distance p

60

70

80

90

100

60

70

80

90

100

errors in Xrelay with Rf = 50.0 ohm
5

X errors in ohm

0
-5
-10
-15
-20
-25
-30

0

10

20

30

40

50

Figure (7) errors in readings with SSSC insertion out of loop with RF =50 Ώ
With increasing the fault resistance Rf, Comparing fig8 with fig.7 we can note that readings error increased
comparing leading SSSC with inactive one . In some details, the errors in resistances increased unlike reactance
errors. Resistances are lower values compared with original case.
50

R errors in ohm

0

-50
r=0.0
r= 0.1
r = 0.2

-100

-150

0

10

20

30

40

50

60

70

80

90

100

70

80

90

100

10

X errors in ohm

0

-10

-20

-30

-40

0

10

20

30

40

50

60

Figure (8) errors in readings with SSSC insertion out of loop with R F =200 Ώ

118

Vol.4, No.14, 2013

www.iiste.org

4.1.2 STATCOM impacting
From equation (10), it is obvious that case of a bolt fault the errors in readings disappeared at any value of
STATCOM injected current (or voltage). The results in figure (9) shows the errors in impedance relay readings
case of fault occurrence through resistance Rf = 50 Ώ. The errors in resistances are positive what shows that the
resistance readings are higher than original case. but the reactance errors are negative values although the
errors are smaller compared with SSSC impacting .
errors in Rrelay with Rf = 50.0
6
Ef=1.0
Ef=1.1
Ef=0.9

R errors in ohm

5
4
3
2
1
0

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

errors in Xrelay with Rf = 50.0
1

X errors in ohm

0

-1

-2

-3

0

0.1

0.2

0.3

0.4

0.5

Figure (9) errors in readings with STATCOM insertion out of loop with R F =50 Ώ
Giving a careful consideration to figure 10 that shows the errors with changing fault location and
STATCOM voltage (Esh) but at higher fault resistance (Rf = 200 Ώ). We can observe that impedance
readings diverge slowly from the original case at inactive STATCOM. Comparing leading STATCOM
with inactive one we can say that the errors are within a convergent values but errors in reactance are
higher than resistance one.
5
Ef=1.0
Ef=1.1
Ef=0.9

R errors in ohm

4

3

2

1

0

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

4

X errors in ohm

2

0

-2

-4

-6

0

0.1

0.2

0.3

0.4

0.5

Figure (10) errors in readings with STATCOM insertion out of loop with R F =200 Ώ
4.2. FACTS at near end
Back to equation (43) when a belt fault occurred, there are errors due to the presence of SSSC inside fault
loop. By scrutiny in figure (11) the errors case of an inactive SSSC are very small in both resistance and reactance

119

Vol.4, No.14, 2013

www.iiste.org

readings but these values differed obviously in leading SSSC devices by increasing the injected voltage as
shown.
Increasing the injected voltage of SSSC resulted in observable values in impedance readings errors as seen in
figure (11).
20
r=0.0
r= 0.1
r = 0.2

R errors in ohm

15

10

5

0

-5

0

10

20

30

40

50

60

70

80

90

100

60

70

80

90

100

4

X errors in ohm

2

0

-2

-4

-6

0

10

20

30

40

50

Figure (11) errors in readings with SSSC insertion in loop with R F =0.0 Ώ
By increasing the fault resistance the situation has became worst and the errors had larger values as obvious by
comparing the results in figs 11 and 12. The errors commensurate with voltage injection were in a negative manner.
The errors in impedance readings are at unacceptable levels in reactance and resistance especially where a fault
occurs at far end side.
20

R errors in ohm

0
-20
-40
r=0.0
r= 0.1
r = 0.2

-60
-80
-100
-120

0

10

20

30

40

50

60

70

80

90

100

70

80

90

100

10

X errors in ohm

0

-10

-20

-30

-40

0

10

20

30

40

50

60

Figure (12) errors in readings with SSSC insertion in loop with R F =200.0 Ώ
4.2.2. STATCOM impacting
During a bolt fault in the presemce of STATCOM the errors in resistance can be neglected but in reactance; errors
are considered as shown in fig13.

120

Vol.4, No.14, 2013

www.iiste.org

6
Ef=1.0
Ef=1.1
Ef=0.9

R errors in ohm

5
4
3
2
1
0

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

errors in X
relay with Rf = 00.0
14

X errors in ohm

12
10
8
6
4
2
0

0

0.1

0.2

0.3

0.4

0.5

Figure (13) errors in readings with STATCOM insertion in loop with RF =0.0 Ώ
By increasing the fault resistance the resistance readings decreased but reactance readings increased slowly see
figure 14.
30

R errors in ohm

20
10
0
Ef=1.0
Ef=1.1
Ef=0.9

-10
-20
-30

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

40

X errors in ohm

20

0

-20

-40

-60

0

0.1

0.2

0.3

0.4

0.5

Figure (14) errors in readings with STATCOM insertion in loop with R F =200.0 Ώ
4.3. FACTS at midpoint
In this case, both SSSC and STATCOM are out of the fault loop for faults on the near half of the line, while
they are introduced in the fault loop for faults on the far half one.
In the presence of SSSC at the midpoint the impedance errors are split into two parts. The right part (from p = 0.5
1 ) is for the faults on the near half of the line, while the left part (p=0  0.5 ) is corresponding to the faults
on the far half.
Figure (15) shows that at a bolt fault in the right part; increasing the injected voltage has no effect on impedance
relay readings. On the other hand in the right side the injected voltage magnitude resulted in decreasing the
measured reactance and increasing the measured resistance.

121

Vol.4, No.14, 2013

www.iiste.org

20

R errors in ohm

15
r=0.0
r= 0.1
r = 0.2

10

5

0

-5

0

20

40

60

80

100

120

140

160

180

200

120

140

160

180

200

10

X errors in ohm

0
-10
-20
-30
-40
-50

0

20

40

60

80

100

Figure (15) errors in readings with SSSC insertion at midpoint with R F =0.0 Ώ
Comparing fig 15 and 16, We can note that case of inactive SSSC the reactance and resistance readings in the left
part are unnoticeable unlike in the right side where errors are more obvious . With increasing the injected
voltage the errors are in negative values in both resistance and reactance readings change obviously.
100
r=0.0
r= 0.1
r = 0.2

R errors in ohm

50

0

-50

-100

-150

0

20

40

60

80

100

120

140

160

180

200

120

140

160

180

200

errors in Xrelay with Rf =200.0
10

X errors in ohm

0
-10
-20
-30
-40
-50

0

20

40

60

80

100

Figure (16) errors in readings with SSSC insertion at midpoint with R F = 200.0 Ώ
4.3.2. STATCOM Impacting
Figure (17) shows case of a bolt fault occurrence during impacting STATCOM at midpoint; in the left part
increasing the injected voltage has no effect on impedance relay readings. On the other hand in the right part the
injected voltage magnitude nearly has no effect on the errors in measured reactance and very small effect on the
measured resistance.

R errors in ohm

15

10
Ef=1.0
Ef=1.1
Ef=0.9
5

0

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

35

X errors in ohm

30
25
20
15
10
5
0

0

0.1

0.2

0.3

0.4

0.5

122

Vol.4, No.14, 2013

www.iiste.org

Figure (17) errors in readings with STATCOM insertion at midpoint with RF = 0.0 Ώ
By increasing the fault resistance; in the left part of fig18 the measured resistance errors increased slightly with
increasing the injected voltage unlike reactance readings that decreased slightly. On the other hand at the
right part with increasing the injected voltage (Esh), very small effects appeared on resistance also on
reactance but with higher manner than resistances. Fortunately the errors in STATCOM are lower than that of
SSSC.
50
Ef=1.0
Ef=1.1
Ef=0.9

R errors in ohm

0

-50

-100

-150

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.6

0.7

0.8

0.9

1

errors in X
relay with Rf = 200.0
20

X errors in ohm

15
10
5
0
-5
-10
-15

0

0.1

0.2

0.3

0.4

0.5

Figure (18) errors in readings with STATCOM insertion at midpoint with RF = 200.0 Ώ
5. Conclusion
The measured impedance has been studied for two cases of FACTS in and out of the fault loop for the installation
points at near end , mid-point , and far end of a transmission line.
Comparing the errors in the measured impedance from its original value, the impedance of the line section
between the relaying and the fault points in the presence of SSSC and STATCOM at the line ends, It can be
concluded as follows :
1- In case of FACTS installed at the far end of the line
a- The measured impedance at absence of fault resistance is the original value. When there was a fault resistance
some errors appeared that could be neglected for STATCOM but can’t for SSSC.
b- When a fault occurred through a resistance RF; increasing the injected voltage of SSSC affected negatively
errors more than increasing STATCOM voltage with the same injected reactive power.
2- In case of FACTS devices at near end of the line
a- There were some errors in impedance readings in SSSC that disappeared with STATCOM case of a bolt fault
occurrence.
b- By increasing the fault resistance the effects of FACTS are more obvious but SSSC effect was higher.
From the above two notes it can be said that the location of impacting FACTS devices from being inside or
outside fault loop is an important uncontrollable factor that affect impedance readings. The operational conditions
(changing Esh and VSSSC) of FACTS also play an effective role in appearance of impedance readings error.

3- In case of FACTS devices at midpoint of the line
When FACTS devices were inserted at midpoint the impedance readings changed obviously between two nearby
points because one of them located inside fault loop and the other one was outside. This study may be considered a
contribution to the protection engineering to easy modify the relay setting according to FACTS impacting
locations and operational condition and that will be our future work (Automatic modification of Distance Relay
Setting Avoiding the Effects of FACTS Devices Impacting).
Appendix (A)
Derivations of relay impedance ZA formulas
1. without facts
Fault current can be calculated as follows:

123

Vol.4, No.14, 2013

www.iiste.org

Ifa (1) = Ifa (2) =Ifa (2)

(A1)

Where Vf0: pre-fault voltage
Z1, Z2, and Z0: are positive, negative and zero sequence
Where RF: fault impedance
ZA1, 2 = ZSA1 + p ZL1, 2

(A2)

ZA0 = ZSA0 + p ZL0

(A3)

ZB1, 2 = ZSB1 + (1- p) ZL1

(A4)

ZB0 = ZSB0 + (1 – p) ZL0

(A5)

Z1, 2 =

(A6)

Z0 =

(A7)

Ze = Z1 + Z2 + Z0 = 2

+

(A8)

To complete the fault current equation, we need now to calculate the pre-fault voltage Vfo at point f. from kirchauf’s
voltage law (KVL) in figure (2), we can deduce these equations.
Vfo = EA – IApf ZA1
(A9)

IApf =

=

Kd = 1 – h e

=

-jδ

(A10)
(A11)

Using Equation (A11) in Equation (A9), Vfo is derived as follows:
Vf0 = EA -

ZA1 =

=

=
Vf0 =

(A12)
(A13)

D = ZB1 + ZA1 h e-jδ

(A14)

Ifa (1) =

(A15)

(1)

Standing on previous equations; we got Ifa .

As been discussed previously, the work of distance relays depends on two measured quantities; current and voltage
measurable values. We now will deduce these two equations of relay voltage and current at bus A. Using
Kirchhoff’s current law (KCL) in figure (5), we can directly get both negative and zero sequence current values as
follows.
IAf (2) =

Ifa (1) = C1 Ifa (1)

(A16)

IAf (0) =

Ifa (1) = C0 Ifa (1)

(A17)

C1, 0 =

(A18)

KCL is applied only on zero and negative sequence currents only, but positive one has voltage sources. So KVL is
used to get Ifa (1) as shown in figure (A3).
EA – EB = IAf (1) ZA1 - IAf (1) ZB1 = IAf (1) ZA1 - (Ifa (1) - IAf (1)) ZB1
= IAf (1) [ZA1 + ZB1] - Ifa (1) ZB1

124

(A19)

Vol.4, No.14, 2013

www.iiste.org

IAf (1) = C1 Ifa (1) +

(A20)

FigureA3. Positive sequence circuit

=

=

x

= Cd Ifa (1)

(A21)

IAf (1) = (C1 + Cd) Ifa (1)

(A22)

Cd =

(A23)

The relay current is:
IAf = IAf (1) + IAf (2) + IAf (0)

(A24)

For single line to ground fault, the measured voltage will be as follows.
VA= IAf (1) Z1+ IAf (2) Z2+ IAf (0) Z0 =Z1 (IAf (1) + IAf (2)) + Z0 (IAf (0))
Z0 IAf

(0)

= Z1 IAf

(0)

- Z1 IAf

= Z1 IAf

(0)

(0)

+ Z0 IAf

(A25)

(0)

+ [Z1 – Z0] IAf (0)

(A26)

Back to Equation (A25), we can get.
V = Z1 (IAf (1) + IAf (2) + IAf (0)) + [Z1 – Z0] IAf (0)
= Z1 (IAf (1) + IAf (2) + IAf (0) (1+3K0L))

(A27)

Where K0L is called zero sequence compensated factor. And it is used to modify zero sequence impedance to be
positive sequence to treat with relay impedance as positive sequence Z1 only.
K0L = 3

(A28)

IAf = (Cd+2C1+ C0 (1 + 3 K0L)) Ifa (1)

(A29)

After finishing deducing the first quantity (current one) formula, we now want to complete the relay measurable
feedings.
From KVL in figure (2) the relay voltage can be calculated as follows.
VA = VF + IAf p ZL1 =3 Ifa (1) RF + (Cd+2C1+ C0 (1 + 3 K0L)) Ifa (1)

(A30)

So after deducting the two quantity formulas, we can get predictable relay impedance at bus A.
Zrelay =

=

(

)

(

)

Zrelay = PZL1 +

(A31)

A.2. with the presence of SSSC
Regardless the location of SSSC the following equations are valid
CSC1, 0 = ZSC / ZL1, 0

(A32)

ZAf1, 0 = ZSA1, 0 + p ZL1, 0

(A33)

ZBf1, 0 = ZSB1, 0 + (1 - p) ZL1, 0

(A34)

A.2.1. Outside fault loop

125

Vol.4, No.14, 2013

www.iiste.org

According to figA:
Zmf 1, 0= (m - p) ZL1, 0

(A35)

ZAm1, 0 = ZSA1 + m ZL1

(A36)

ZAm1, 0 = ZSA1, 0 + m ZL1, 0

(A37)

ZA1 = ZAf1

(A38)

ZA0 = ZAf0

(A39)

Calculating the pre-fault current
IApf =
=

(A40)

Pre-fault voltage
Vf0 = EA – IApf ZA1 = EA - ZA1 x
[

=

]

=

(A41)

And with the same previous algorithm in (1), we found that
Dn = ZBn1+ZA1 [h e-jᵟ - re jϒ]
-jϒ

(A42)

-j

Kdn = 1+ re - h e ᵟ

(A43)

A.2.2. Inside fault loop
Zmf1, 0 = (p – m + CSC1, 0) ZL1, 0

(A44)

Zmf1, 0 = (p - m) ZL1, 0

(A45)

ZA1, 0=ZAm1, 0+Zmf1, 0

(A46)

ZB1, 0 = ZBf1, 0

(A47)

Now calculating the pre-fault voltage but noting that the pre-fault current didn’t change because the series voltage
location doesn’t affect the line current.
Vf0 = EA + VSSSC – IApf ZAn1

(A48)

Vf0 = EA + VSSSC [

Vf0 =

ZAn1

(A49)

]

(A50)

Dn = ZB1 (1 + r e jγ) + ZAn1 he-jδ

(A51)

Ifa (1) =

(A52)

Calculating fault current

Calculating the fault relay current. Firstly; using KVL in positive sequence circuit, the positive sequence is derived
as follows:
EA + VSSSC – IAf (1) ZAn1 = EB – IBf (1) ZB1
But IBf

(1)

=Ifa

(1)

– IAf

(A53)

(1)

IAf (1) =
+ C1 Ifa (1) =

=
=

+ C1 Ifa (1)
x

Ending to eq (A14); the relay fault voltage can be calculated as follows:

126

= Cd Ifa (1)

(A54)

Vol.4, No.14, 2013

www.iiste.org

VA = VF – VSSSC + IAfC (p + CSC1) ZL1

(A55)

IAfC: is the compensated relay current due to using ZL1.
VA = 3RF Ifa (1) - VSSSC + (Cd + 2C1 + C0 (1 + 3K0L)) Ifa (1) (p + CSC1) ZL1
- (CSC1 - CSC0) C0 (1 + 3K0L) Ifa (1)

(A56)
(A57)

VSSSC / Ifa (1) =
= (Ze + 3 RF) (ZAn1 + ZB1) re jϒ / Dn

(A59)

CZ = - (CSC1 - CSC0) C0 (1 + 3K0L) ZL1
CV = - (Ze + 3 RF) (ZAn1 + ZB1) re

(A58)

jϒ

/ Dn

(A60)
(A61)

A.3. with the presence of STATCOM
CSh1, 0 = ZSC / ZL1, 0

(A61)

Zmf1, 0 = |p – m| ZL1, 0

(A62)

ZAm1, 0= ZSA1, 0 + m ZL1, 0

(A63)

ZBm1, 0 = ZSB1, 0 + (1- m) ZL1, 0

(A64)

A.3.1. STATCOM Outside fault loop
From figure (A2a), we derived the fault current
calculations are as follows:

using curtailment the Esh and EB sources to one source eq, the fault

CB1, 0 = ZSh / (ZSh + ZBm1, 0)

(A65)

ZB1, 0 = Zmf1, 0 + CB1, 0 ZBm1, 0

(A66)

EBn =

(A67)

ZBm1, 0n =

(A68)

hn =

(A69)

δn = ∟EA - ∟EBn
ZA1 , 0= ZSA1, 0 + p ZL1, 0
ZB1, 0 =

+ ZmF1, 0 = ZBm1, 0n + ZmF1, 0

127

(A70)
(A71)
(A72)

Vol.4, No.14, 2013

www.iiste.org

Figure (A2). The system circuit during fault (a) outside fault loops

(b) inside fault loop

The pre-fault current IApf can be calculated as follows:
–

IApf =

=
=

=

Kdn = 1- esh (1 – CB1) – CB1 h e
Eq (A73) is used to get the pre-fault voltage at point p.
VF0 = EA – IApf ZA1 =EA -

-jδ

(A74)

ZA1

= EA =

(A73)

ZA1

(

)
(

=

)

=

(A75)

Dn =

(A76)

CB1 =

(A77)

Ifa (1) =

(A78)

IAF (2) =
IAF

(0)

Ifa (1) = C1n Ifa (1)

=

Ifa

(1)

= C0n Ifa

(A79)

(1)

(A80)

KCL can be applied only on zero and negative sequence circuits only, but positive one has voltage sources. So we
will use KVL to get Ifa (1) in positive sequence circuit.
EA – EBn = IAf (1) ZA1 - IAf (1) ZBn1 = IAf (1) ZA1 - (Ifa (1) - IAf (1)) ZBn1
= IAf (1) [ZA1 + ZBn1] - Ifa (1) ZBn1
(A80)
IAf (1) = C1n Ifa (1) +

(A81)

=
=

x

IAf (1) = (Cdn + C1) Ifa (1)

128

= Cdn Ifa (1)

(A82)
(A83)

Vol.4, No.14, 2013

www.iiste.org

Cdn =

(A84)

The compensated relay current is:
Irelay = (Cdn + 2C1n + C0n (1 + 3K0L)) Ifa (1)
Now we want to calculate the second quantity (voltage quantity) with the same technique used previously.
VA = VF + IAF (1) P ZL1 + IAF (2) P ZL1 + IAF (0) P ZL0

(A85)

= 3 RF Ifa (1) + P ZL1 [IAF (1) + IAF (2) + (1+3K0L) IAF (0)]
= 3 RF Ifa (1) + P ZL1 (Cdn +2C1n + (1 + 3K0L) C0n) Ifa (1)
Now we can get the relay measured impedance.

(A86)

Zrelay = PZL1 +

(A87)

CA1, 0 = ZSh / (ZSh + ZAm1, 0)

(A88)

ZA1, 0 = Zmf1, 0 + CA1, 0 ZAm1, 0

(A89)

A. 3.2. STATCOM Inside fault loop:

Following the same sequence in calculating fault current
CV =

(A90)

IAm (1) = CV Ifa (1) + CA1 Imf (1) = CV +CA1 (Cdn + C1) Ifa (1)
IAmc = [CV +CA1 (Cdn + 2C1) + CA0 C0 (1 + 3K0L)] Ifa

(1)

(A91)
(A92)

The relay voltage VA is driven as follows:
VA = VF +Imf Zmf1 + IAf mZL1
= 3 Ifa

(1)

RF + Zmf1 (Cdn + 2C1 +C0 (1 + 3K0L)) Ifa (1)
+ mZL1 (CV +CA1 (Cdn + 2C1) + C0 CA0 (1 + 3K0L) Ifa (1)

(A93)

Zmf1 (Cdn + 2C1 +C0 (1 + 3K0L)) =
Zmf1 (CV +CA1 (Cdn + 2C1) + C0 CA0 (1 + 3K0L))
+ Zmf1 (- CV + (1-CA1) (Cdn + 2C1) + C0 (1- CA0) (1 + 3K0L))

(A94)

VA = 3 Ifa (1) RF + pZL1 (CV +CA1 (Cdn + 2C1)
+ C0 CA0 (1 + 3K0L)) Ifa (1) + Csh Ifa (1)

(A95)

Csh = Zmf1 (- CV + (1-CA1) (Cdn + 2C1) + C0 (1- CA0) (1 + 3K0L) (A96)
From equations (A96) and (A93) by compensating the bus A current, we can get the relay measured impedance as
follows:
Zrelay = pZL1 +

(A97)

Where
Dn =ZAn1 + ZB1 – (CA1 + eSh (1 – CA1) – h e -jδ) (Cm ZAm1 + ZAn1)

(A98)

Cm =

(A99)

References
A. Kazemi, S. Jamali, (Sep 2009) “Comparing Impacts of SSSC and STATCOM on Measured Impedance at
Relaying Point” IEEE transactions, power delivery.
Tarlochan Singh Sidhu, Rajiv K. Varma, Pradeep Kumar Gangadharan and Fadhel Abbas Albasri, (JULY 2005),
“Performance of Distance Relays on Shunt FACTS Compensated Transmission Lines “IEEE transactions, power
delivery, VOL. 20, NO.3.

129

Vol.4, No.14, 2013

www.iiste.org

M. Khederzadeh, sep2009, “Impact of SSSC on the Digital Distance Relaying”, IEEE transactions, power delivery.
S Jamali and H Shateri, (2011), “Locus of apparent impedance of distance protection in the Presence of SSSC
“EUROPEAN TRANSACTIONS ON ELECTRICAL POWER Euro. Trans. Electr. Power.
Prof. Detlef Schulz, (2007), “Improvement in the Quality of Delivery of Electrical Energy using Power Electronics
Systems”, Helmut Schmidt University, pp27-39.
Yoshihide Has, (2007), “Handbook of Power System Engineering”, pp315-325.
Y. Q. Xia, K. K. Li and A. K. David , (Jan 1994) ,“Adaptive Relay Setting for Stand-alone Digital Distance
Protection”, IEEE Transactions on Power Delivery, Vol. 9, No. 1,.
Yi Hu, M, Damir Novosel and Murari Mohan Saha, (Jan 2002), “An Adaptive Scheme for Parallel-Line Distance
Protection “, IEEE transactions, power delivery, VOL. 17, NO. 1.

130

This academic article was published by The International Institute for Science,
Technology and Education (IISTE). The IISTE is a pioneer in the Open Access
Publishing service based in the U.S. and Europe. The aim of the institute is
Accelerating Global Knowledge Sharing.
More information about the publisher can be found in the IISTE’s homepage:
http://www.iiste.org
CALL FOR JOURNAL PAPERS
The IISTE is currently hosting more than 30 peer-reviewed academic journals and
collaborating with academic institutions around the world. There’s no deadline for
submission. Prospective authors of IISTE journals can find the submission
instruction on the following page: http://www.iiste.org/journals/
The IISTE
editorial team promises to the review and publish all the qualified submissions in a
fast manner. All the journals articles are available online to the readers all over the
world without financial, legal, or technical barriers other than those inseparable from
gaining access to the internet itself. Printed version of the journals is also available
upon request of readers and authors.
MORE RESOURCES
Book publication information: http://www.iiste.org/book/
Recent conferences: http://www.iiste.org/conference/
IISTE Knowledge Sharing Partners
EBSCO, Index Copernicus, Ulrich's Periodicals Directory, JournalTOCS, PKP Open
Archives Harvester, Bielefeld Academic Search Engine, Elektronische
Zeitschriftenbibliothek EZB, Open J-Gate, OCLC WorldCat, Universe Digtial
Library , NewJour, Google Scholar

Errors analysis in distance relay readings with presence of facts devices

More Related Content

What's hot (20)

Viewers also liked (8)

Similar to Errors analysis in distance relay readings with presence of facts devices (20)

More from Alexander Decker (20)

Recently uploaded (20)

Errors analysis in distance relay readings with presence of facts devices