Ferroresonance in Power Systems
Ferroresonance in Power Systems
Ferroresonance in Power Systems
solutions
J.R. Marti, BEng, PhD, MEEPE
A.C. Soudack, PhD
/ /I
a
introduction
I1 d
Fig. 1
solution
oi A
"I i
soluilon
oi B
(unstable)
vc
soluilon
at C
Fig. 2
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
32 1
.. .
closed-form solution is then used to draw contour maps
of the system, showing the safe and ferroresonant regions
as a function of the circuits parameters.
To verify the validity of the analytical results, the
circuit was solved numerically with a Runge-Kutta integration routine. Excellent agreement was found between
the analytically-predicted steady-state values and the
numerical solution.
An important observation derived from this study is
the effect of the losses in determining the value of CcririCol
(proportional to the feeders length) beyond which ferroresonance can occur.
The methods of analysis presented are being extended
for the analytical prediction of harmonics and subharmonics in steady state solutions and for the characterisation of jump and chaotic phenomena in certain regions
of the parametric maps.
2
Problem description
The operation of a power transformer under ferroresonance can be illustrated with the diagrams of Fig. 2,
which show the relationship between the fundamental
frequency components of voltages and currents in the
circuit of Fig. 1 ignoring the losses R for simplicity. In the
development of the general case, losses are shown to be
an important consideration in the determination of the
critical feeder length. Upper case letters are used to indicate fundamental frequency phasor quantities.
In Fig. 2a, the intersection of the straight-line E, - V,
(representing the voltage applied to the transformer coil)
with the transformers magnetisation characteristic is
possible at three points: points A and C are stable operating points. Point B is unstable. The instability of point
B can be seen, for instance, by increasing the source
voltage by a small amount. In Fig. 2a, increasing E, displaces the line E, - V, upwards, parallel to itself. When
this happens, the current at intersection points A and C
also increases, as expected. For point 8, the current
decreases, which is not physically possible. Point A corresponds to normal operation in the linear region, with flux
and excitation current within the design limits. Point C
corresponds to the ferroresonant condition, characterised
by saturated flux and large excitation current. Phasor
diagrams for these operating conditions are shown in Fig.
2b.
Mathematically, both points A and C are equally valid
solutions to the circuit of Fig. 1. Which operating point
0.4
Transformer modelling
+ b4
(1)
0.4
I . . . . .
I. p.u.
Fig. 3
322
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
01
0.2
0.4
0.6
I,
Fig. 4
0.8
1.0
p.u.
+ 0)
(5)
sin w, t
+ mycos ost
(6)
where
CP=
~ CP
COS
e
(7)
+"
Losses
where
sin"(w, t
+ 0) = k, sin (U, t + 0)
+ k, sin [3(w,t + e)]
+ ... + k , sin [n(w, t + e)]
(8)
(9)
Substitutingeqn.6 into eqn. 3, then the steady-state solution (4 and 0 constant) is
Fundamental solutions
With the magnetisation curve of the transformer approximated by eqn. 1, with a general exponent 'n',
i = ~4 b @ , the differentia1 equation for the circuit in
Fig. 1 is
d24
dt2
1 d+
+- + -1 (a4 + b4")= w,E cos oat
RC dt C
1 d4
-- + CO:+
RC dt
(2)
(3)
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
323
= klw: V"-l(aX
sin ost + aycos w,t)
Substituting into eqn. 10, the resulting equation can be
separated into sine and cosine terms and the following
relationships can be established:
[-(of - W@
0
Qr = 0
+ k1~:W"-']@, - RC
[-(of - w@
+ k,w:
@"-']ay
+ ax= w,E
(1 1)
1, = P l c - Po
I2 = p 2 p +1)/2 - 5"
(12)
(21)
(22)
with
az= @:
+ a;
In what follows, eqns. 11-13, which define the fundamental component of the transformer flux in the circuit of
Fig. 1, are manipulated into a more convenient form for
the graphical location of the solution points. The effects
of the circuit parameters, e,, C and R, in the location of
the borders separating the normal and ferroresonant
regions are investigated.
Eqns. 11 and 12 can be combined by squaring and
adding. After some algebraic manipulations (Appendix),
the following single relationship is obtained:
5" - p z , p + l ) / Z + P l t - Po = 0
(14)
where
c = az
(15)
Fig. 5
in which
+ k,w:W"-'
A =(wf - wt)
B=O"
RC
324
o = +Jc;
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
G-12
551
mS
//G=lOmS
5-
4-
3-
E=O
0
1
--22 y
0.0
0.5
1.5
1.0
2.5
2.0
5
0
0.0
1.0
0.5
1.5
2.5
2.0
I
0
GI
Fig. 6
62
G2
G
E2
b General cax
e Gemral case
Fig. 7
E
b
a Teat system
b Gcncralcase
e
cicncraluae
(n + l)K10:(0f
[y](n
p$'
-(U:
= Io
where
- 03 (n - l)<("+')/'
- E' = 0 (24)
~@r("-~)/'
- n(K, o:) r
(26)
Changing the value of the capacitance C affects coeficients pi and p 2 in eqns. 21 and 22 and, therefore, both
curves I, and 1, in Fig. 5. A plot of q5 against C, similar
to the ones previously obtained for G and E (Figs. 6b and
7b), can be obtained more directly starting from the original form of the solution to the circuit (eqn. 14). After
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
325
some algebra
(a:[ - E2)C2- 2(a[
+ klb[("+')'Z)C
(28)
For a given value of [ (or of the transformer flux
4 = +,/<), eqn. 28 gives a second-order equation for C
of the form aCz + /?C + y = 0. Since C has to be real, the
discriminant of this equation, A = /?' - &y, has to be
positive. With this constraint considered, the graph of {
(or 4 = + J [ )against C has the form indicated by the
solid curves in Fig. 8. The dashed curves in this Figure
are the asymptotes of 4 against C for the case with no
damping (G = 0).
10
I
c:i
2
I
ccr
C
Effect of capacitance on circuit solution
damping
_ _ ~ _no damping
Fig. 8
(29)
In Fig. 8, the point of vertical tangency C,, which is
slightly different from the value of C at + b , separates the
possible regions of operation of the circuit. For values of
C less than C, no ferroresonance can occur. The value
of C, could be exactly determined from eqn. 28 with
dC/d[ = 0. The expression for C, is rather elaborate in
this case. The value of C at 4b is much easier to determine (eqn. 29) and is very close to C, . Point C, can also
be determined graphically from a 4 against C plot, as in
Fig. 14 for the test system.
In the case of the capacitance, there are only two possible regions of operation (as opposed to three regions for
the conductance G or the source voltage E). In region I
there is only one solution point. In region I1 ferroresonance may or may not occur. There is no region 111,
as for G (Fig. 66) or for E (Fig. 7 4 , where ferroresonance
will always occur.
When losses are ignored (dashed curves in Fig. 8), the
criticd capacitance value (Cg) occurs at a smaller value
of C than when losses are considered (Cm). As losses
increase from zero, the dashed curves in Fig. 8 'strangle'
the zone around Cm and at some value of G they split
326
"1
N
"2
r
b
Fig. 9
Test system
a Systemdiagram
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
0
x
Fig. 10
=-,
= 0.5 v, A
c, + 2c.
c = c, + ZC,
The magnetisation curve of the transformer is approximated by the 11th order two-term polynomial of eqn. 1.
The resulting approximation,shown in Fig. 4, is
i = 0.28 x 10-'1$
+ 0.72 x
10-2Q1'
(30)
-0.51
0.0
1.21
0.8-
ki = 0.451
o2= 85 rad/s (13.5 Hz)
2.5
1.0-
2.0
1.5
1.0
(31)
111
0.5
0.40.6
0.00
IO
15
20
25
30
G . pu.x IO-'
0.0
0.5
1.0
15
2.0
2.5
E , P.U.
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
321
Critical resistances
Critical source voltages
Critical capacitance
R, = n/a
E, = n/a
C,,
R , = 21.4 kfl
E, = 1.02 P.U.
= 509 nF
2.4~
-6
Y
1.2
1.0
0.0
0.5
0.01
1.0
1.5
2.0 2.5
c. p.u.
0.02
3.0
0.03
3.5
4.0
Fig. 15
328
Conclusions
An analytical technique for obtaining steady-state ferroresonant conditions in iron core transformers supplied
through capacitive coupling was presented. The analysis
proposed is general for the type of circuit considered and
includes the effect of the system losses.
The solution to the circuits nonlinear differential
equation is based on Ritzs method of harmonic balance.
To apply this technique, the transformer saturation curve
is approximated by a two-term 11th-order polynomial.
This 11th order polynomial provides a much better fit to
the saturation region than can be obtained with seventh
and lower order approximations. Differences of about
20% in the magnitude of the transformer flux were
observed in comparisons of the 11th order against
seventh order approximations.
Graphical solutions were used to map the boundaries
between normal and ferroresonant regions and to locate
the possible operating points as a function of the circuits
equivalent source voltage, series capacitance and shunt
losses.
A practical example of a system under ferroresonant
conditions was presented. The solution points as well as
parametric maps were calculated for this system. One
interesting observation from this study is the importance
of considering the transformer core losses in determining
the margin of safety before fenoresonance can occur.
This margin went from less than 2.6 km with no losses to
65 km when losses were considered. Another interesting
observation is that for certain combinations of circuit
parameters, ferroresonance can occur even for very small
values (E + 0) of the source voltage.
12
References
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
13
Appendix
and
H=o,E
(32)
+ bx = H
(33)
ay
a2(x2 y 2 ) + b2(x2 + y 2 ) = H 2
where
x=
(34)
= cPZ, and with the
+ k:o:@2"-2 - Yo: - 03
x klog W-' + O,'/(RC)~]@
-~U: E'
(k:o:)@'
- 0;)
+klo;w-l]
= [-CO,'
b = oJ(RC)
- 2(of
- o@k,w:@"+'
= @'
+ o . Z / ( R C ) ~ ] @-~0,'E'
=0
ir I
=0
Rearranging
Y=
a
+ @);
II
Authorized licensed use limited to: Universidad de chile. Downloaded on December 1, 2008 at 01:43 from IEEE Xplore. Restrictions apply.
329