Watson Thesis
Watson Thesis
System
Equival nts fo nic Studi s
an
Transie Convertor Simulatio
A thesis
presented for the degree of
Doctor of Philosophy in Electrical Engineering
in the
University of. Canterbury,
New Zealand
by
1987
i.
TABLE OF CONTENTS
Page
List of Figures v
List of Tables x
List of Principal Symbols xi
Abbreviations xii
Abstract xiv
Acknowledgements xv
CHAPTER 1 INTRODUCTION
CHAPTER 2 OVERVIEW
3.1 Introduction 13
3.2 Representation of Components 13
3.2.1 Static Convertor 13
3.2.1.1 Modes of operation 14
3.2.1.2 Convertor control systems 14
3.2.1.3 Types of control schemes 15
3.2.2 Synchronous Machines 16
3.2.3 Transmission Lines 17
3.2.4 Transformers 17
3.2.5 Static Shunt Elements 18
3.2.6 Static Series Elements 18
3.2.7 A.C. System 19
ii.
Page
Page
Page
REFERENCES 164
LIST OF FIGURES
Figure Page
Figure Page
Figure Page
Figure Page
Figure Page
LIST OF TABLES
Table Page
6.1 Phase "A" Current Harmonics for simple test system 135
6.2 Phase "A" Voltage Harmonics for simple test system 135
6.3 Phase "A" Current Harmonics for modified test system 131
6.4 Phase "A" Voltage Harmonics for modified test system 131
6.5 IHA/TCS comparison with Thevenin a.c. system representation
used in TCS
(a) Phase "A" harmonic currents 138
(b) Harmonic Voltages 138
(c) Comparison of d.c. Current Harmonics 138
6.6 IHA/TCS comparison with frequency-matched a.c. system
equivalent used in TCS.
(Implicit mutual coupling representation)
(a) Phase "A" harmonic currents 139
(b) Harmonic Voltages 139
(c) Comparison of d.c. Current Harmonics 139
6.1 Comparison of Computational Efficiency 141
6.8 Harmonic Voltages throughout the Lower South Island of
New Zealand's primary transmission network 143
6.9 Harmonic branch currents at the sending end throughout the
Lower South Island (N.Z.) test system 144
SYMBOLS
Subscripts:
a - alpha node
B - beta node
y - gamma node
1· 1.th e 1 emen t 0f vec t or
ij element in ith row and jth column of matrix
~ - inductive branch
r - resistive branch
c - capacitive branch
ec - Equivalent circuit
cb - correction branch
rr - Required response
- denotes a vector
Superscripts :
Abbreviations
ABSTRACT
Acknowledgements
Finally thanks are due to my family and friends for their support
and understanding.
1.
CHAPTER 1 INTRODUCTION.
ratings increase the relative a,c. system strength diminishes and the
problems of dynamic response and harmonics increase greatly.
Load flow studies determine the voltage, current, power and power-
factor or reactive power at various points in a power system. This is a
3.
Under fault conditions the a.c. system parameters vary due to the
finite time required for the flux to change in components such as
synchronous machines. A means of obtaining accurate time varying network
equivalents was developed by Heffernan (1980). This involved interfacing
the dynamic simulation program with a multi-machine a.c./d.c. transient
stability program. The network equivalents obtained from the latter
program provided a more accurate indication of the time response of the
a.c. system by virtue of the synchronous machine models used in the
transient stability program.
It has been known for a long time that the frequency response of
the power system components affects the transient response due to the
multitude of frequencies present in transient waveforms. However,
incorporating such frequency-dependence in dynamic simulation programs
has proved difficult, The present work develops methods of obtaining a
practical and computationally efficient equivalent circuit that
accurately represents the frequency-dependence of the actual system being
represented.
5,
Chapter 2 OVERVIEW.
representation is a problem that has not been solved with this method.
this stage four different a.c. system models were catered for. A more
realistic d.c. fault model was introduced based on experiments performed
by Kohler (1967). At this stage the two main deficiencies were apparent;
the transformer model, and the representation of the frequency-dependence
of an a.c. system. The development of a transformer model that represents
both saturation and hysteresis effects was accomplished by Joosten
(1987), while this thesis reports on the later problem.
EMTP TCS
3.1 Introduction.
The Graetz bridge is well accepted as the basic building block for
HVDC schemes and its operation well documented (Adamson and Hingorani
1960, Cory 1965, Kimbark 1971, Uhlmann 1975, Arrillaga 1983). Both 6 and
12 pulse convertor models are available, with the 12 pulse consisting of
two 6 pulse bridges connected in series on the d.c. side and via a
star-star and star-delta convertor transformers respectively. This
configuration of convertor transformers gives the 30 degree phase shift
required for 12 pulse operation.
1 4.
A Non-commutating mode
Two valves on different sides and arms of the bridge are ON.
B Normal commutation
Three valves in ON state ; one on each arm and at least one
conducting on each side of bridge.
C Non-commutating arm short-circuit
Both valves on one of the arms in ON state with no other valves ON.
D Commutating arm short-circuit
Both valves on one of the arms in ON state with one or both valves
on another arm being ON. There is a single commutation process if
one valve on the second arm is conducting, otherwise there are
commutations taking place on both sides of the bridge
simultaneously.
E a.c. short-circuit
Four or more valves conducting, involving all three arms of the
bridge. This has the effect of short-circuiting the convertor
transformer's three secondary terminals.
- Invertor
E.A.C.
C.C.C.
Individual Phase-Control (IPC) was the norm for early HVDC schemes.
With this scheme valve firing instances are determined independently for
each valve based on the relevant a.c. voltage zero crossing. IPC reflect
any unbalance or distortion in the a.c. supply voltages by producing
variations in valve conduction periods. These in turn produce further
distortion and unbalance, with the result of uncharacteristic harmonic
production in the a.c. system, and harmonic instability if it reinforces
the original distortion (Ainsworth 1968). These phenomena are common with
weak a.c. systems. Filtering the a.c. system waveform entering the
16.
vg peL
g
. I )
g
+ Rg . I g
+ I + R • I (3.1)
g g g
L is a 6x6 matrix and its elements are given in appendix A1. The
g
stator values can accommodate fourth harmonic terms if available, as
these can be of significance in convertor-generator units (Campos-Barros
1976) .
17.
3.2.4 Transformers.
LI2
II RI R2 I 2
'" L21
~
t
VI LI I L22
f
V2
I I
0
Neglecting the iron losses, the matrix equation for the coupled coil
model is ;
(3.2)
Both a.c. and d.c. harmonic filters and surge capacitors greatly
influence the convertor waveforms following a disturbance and even in the
steady state, and hence need to be accurately represented. Fig. 3.3
displays the various shunt elements that can be used to represent plant
components and faults. Three phase a.c. components are assumed to be star
connected with the neutral point earthed.
I I
Figure 3.3 The various static shunt element models.
These models are essentially the same as for the static shunt
elements except that the components are connected between phases. This
allows circuit breakers, delta loads and other such components to be
modelled. Care must be taken to ensure that the topology constraints, as
discussed in section 3.3.6, are not contravened when using these
elements.
19.
The size and complexity of a.c. power systems necessitates the use
of less complex, yet sufficiently accurate a.c. system models. Explicit
busbar information will generally not be available at all busbars and
hence system components of particular interest must be represented in
detail. A number of a.c. system equivalents have already been developed
and are displayed in figure 3.4. The major feature of the a.c. system
models developed by Heffernan (1980) was their time varying nature. The
present work involves the incorporation of the frequency response of the
system, as expressed by impedance loci, into an a.c. system model. The
details of obtaining the impedance loci and the synthesis of frequency
matched a.c. system models is presented in chapter 5.
E L
s 5
(l-a_
_ IL 5
a
(l-aIL
_ _ 5
a
where
u represents the input voltages and currents.
x is the vector of state variables.
¥ represents the output voltages and currents.
z represents a set of dependent variables.
For each time step the change in a state variable is equal to the
integral of the area under its derivative and the trapezoidal integration
21.
(3.8)
(3.9)
(iv) Steps (ii) and (iii) are performed iteratively until convergence
is reached. Convergence is deemed to have occurred when all the
state variables satisfy
(3.10)
(3.11)
lowered so that the simulation step falls exactly on the required time
the change is to take place. As the state space coefficient matrices are
not functions of step size, changing it does not impose a dramatic
computational burden.
x A.sin(wt+<jl) (3.12)
wA.cos(wt+<jl) (3.13)
where
lk is the inductance
Lk the inductive reactance
c is the capacitance
k
C the capacitive susceptance
k
Wo the base angular frequency
(-i
if node i is the sending end of branch b
t
Kbi if node i is the receiving end of branch b
if branch b is not connected to node i.
Kt = [ Kt ,K t S ,0 J <3.16)
rn ra r
Kt = [ Kt 0 ,0 J <3.17)
cn ca'
ec I
s
Resistive branches.
0.18)
Inductive branches.
E (L I ) - R I + Kt V + Kt V Kt V
~-p ~ ~ ~ ~ ~a a ~8 8+ ~y y o 0.19)
Capacitive branches.
t V )
ec P(K Ca I 0.21 )
a c
27.
KaJl.IJI.+ K I = 0
as s
or KaJl.pIJI.+ K pI = 0
as s I C3. 23)
K I + K~Jl.IJI.+ K I = 0 C3. 24)
~r r ~s s
K I + K I + K I = 0 ,I C3. 25)
'(c c '(r r K,(Jl.IJI.+ '(s s
t
Premultiplying equation (3.21) by K and substituting into
ac
equation (3.25) yields:
pQ = - K I - K I - K I C3. 26)
a aJi. JI. ar r as s
where
28.
where
State-Related Variables.
-1
I~= L~ 'I'~ (3.29)
V = C- 1Q (3.30 )
a a a
Dependent Variables.
1 t
V = - R (K I + Koolo+ K R- K V ) (3.31 )
8 8 8s S ~N N 8r r ra a
1 t
I = R- (K V + KrtoVo) (3.32 )
r r ra a ~ ~
(3.33 )
State Equations.
P'l' ~ = E~- pL ~ I ~- R~ I ~ + Kt~a Va + Kt~8 V8 + Kt~y Vy (3.34 )
pQ = - KNo1o- K I - K I (3.35)
a ~N N ar r as s
where
29.
3.4 Extensions.
removal of this constraint from the TCS formulation would require major
modifications to the program. A possible solution is the addition of a
very small capacitance between earth and each phase in order to
convert the coupling nodes into alpha nodes; however, this would result
in small integration step lengths and thus deterioration of computational
efficiency.
~e [~p,u'J
-R
wL vp.u. 0.36)
p.u. [ /j Zb
wL
p.u.
where V is the per unit voltage across the RLC branch.
Figures 3.7 and 3.8 provide some insight into possible iterative
algorithms. The internal iterative loop estimates the extra network in
figure 3.7 and the standard TCS circuit in figure 3.8. An alternative
approach, and the method adopted, consists of embedding the extra state
variables as part of the standard variable iterative loop as shown in
figure 3.9. The estimates for the new state variables and the state
variable derivatives are evaluated for the standard and extra variables
at the same stage and the convergence criteria reflects the convergence
of both sets of variables. Compared to the use of separate iterative
loops for the standard and extra variables this solution permits
considerable computational saving. However the separate iterative loop
method requires far less modification to the original TCS algorithm.
Has
convergence No
been reached for the >---~~---------- __~
extra network
variables
Has
No convergence
been reached for the
standard TCS
variables
?
Yes
Has
convergence No
been reached for the)---~~ ____________~
standard TeS
variables
?
Yes
Has
No convergence
been reached for the
extra network
variables
?
Yes
convergence
No
been reached all
variables (both
&
The spectral analysis facility has been used for the analysis of
earthquake records as well as the time waveforms derived from TCS, The
results may be displayed either in tabular form, as a histogram or as a
two-dimensional plot,
The usual practice in the past has been to have separate programs
to perform each distinct major step in the analysis. This therefore
required read and write operations from disc at each step. The speed of
operation can be greatly increased by minimizing the number of read and
write operations. This is achievable by lumping several of the required
steps into one program, as has been done in ACREP. The results of one
step are stored in the program variables (virtual memory), ready to be
used by the next step, rather than writing out a data file and then
reading it in. Virtual memory is the area set aside by the computer for
storing the program variables while the program is executing. When the
program stops execution the contents are lost and this memory is
allocated to store the variables of any another image that starts
execution.
This use of program vectors rather than data files also leads to
increased accuracy and reduced build-up of error. The writing out of a
formatted file windows the data in the program variables by the format of
the write statement, hence causing a loss of accuracy. Where there are
several programs in a chain, each introduces a loss of significance by
writing a data file and the total error builds up. By keeping data in
program vectors maximum accuracy for the given data type is obtained.
The disc 1/0 subroutines are devoted to the reading and writing of
data files. They also perform the additional function of transferring
data between the working set and data base regions of memory.
Local
Program
Vectors
Plotter
DISC =8 Computatioo
Debug
Software --<--
~ ;/
Tenninal
~--
------
Main
Program
Vectors
I;
:>
.
0..
. 0.4
v
W
0
:>
I- 0m 2
H
Z
(j
-<
l:
0.0 I I
IZI 5 11Z1 15 2121 25
FREQUENCY (HARMONIC ORDER)
" . 1. 2
::J
0..
. 0. 8
v
w 121. 4
o
::J
I-
H 121. 121
Z
(j
-< -IZI. 4
l:
I-
Z -IZI. 8
W
!r
!r
:>-1.2
U IZI 6121 12121 18121 24121 3121121 360
T I ME (DEGREES)
121.8
W
0 121.6
:J
I-
H
Z
!.j 0. 4
-<
2:
0. 2
o. 121 -!I.1.Wl1JJ.W.LLW.Iw.w.
1. 121 ~----~---r--~------~------r-----~------~------~~
0. 8
w
0 0.6
:J
I-
H
Z
!.j 0. 4
«
2:
0. 2
50 1 00 1 50
FREQUENCY (HARMONIC ORDER)
histogram as illustrated in figures 4.4 and 4.5, the latter being the
continuous version of the former. In addition to the graphical output the
information can also be written to a data file or terminal screen, in the
form of a table.
w
0 5
~
~
H 4
~
Z rl
~ 3
~
C
~
2
2
W Lru
~
~ 1
~
~
~ 0
0 1 2 3 4 5
>
ORDER OF HARMONIC
Figure 4.6 Isometric Plot
w
o
::J
f-
H~
Z+l
·rl
19 C
<1::J
L
W L
19 ill
<1~
f-
.J
o
>
ORDER OF HARMONIC
Figure 4.8 Busbar voltage harmonics throughout the Lower South Island
of New Zealand when one per unit of each harmonic is injected
w
o 1.0
:J
I-
H~
Z.j.J
.r!
l'J C
i::J 0.5
W L
l'J UJ
<t:-~
I-
-.J 0.0
o
>
ORDER OF HARMONIC
0.5
.........
r 0.4
H
Z
:J
a: 0.3
W
n.
0.2
W
(!)
<t:
r 0.1
-.J
0
>
0.0
1 2 3 4 5 6 7 B 9
BUSBAR NUMBER
0.5
t-
H 0.4
Z
:J
..... .. '
.. ' .. ,
a: 0.3
W .
.....I 1'...,.
[l. . I ..~ I
>
0.0
1 2 3 4 5 6 7 8 9
8US8AR NUMBER
W
0 O.B A PHASE
:J
I-
H~
Z.j.J
('Jrl 0.6
<t C
L:J
W L 0.4
('J I]J
<{~
I-
~
0.2
0
>
0.0
1 5 9 13 17 21 25
ORDER OF HARMONIC
0- o ...,. f\) PJ o ~ 0 ~ 0
C!.l 'V
CD ::s
.0 AOXElUAGH-220
Pl ROxeURGH-Z20
'"rJ C Ul
TO TO
1-'- CD RQXeURGH10lil
CD ROXeURGH:1 0 11
()q ::l
C 0
'J CD ;p
=
ROXeUAGH-O 11 AOXeUFlGH-O 11
CD TO TO
0 AQX15URGH-220
= ROXBUAGH-220
I
.!=" 0 I
a l> »
'D ]] :n
I-IANAPOURI220 MANAPQUR.I220
W 0 3: TO -
3:
TO a
::l MANAPOUR101lA a MANAPOUA" 0 14
CD z Z
'V ...,. H ...,. H
::l
0 cT ....,j [1 ---J [1
1-'- MANAPOUAX:ii:ZO MANAPQUAI2ze
Ul TO rt TO rt
::l MANAPOUR201A :r [1 OJ MANAPOUA201A "J [1
cT OJ
]] c ]] C
'V l> I ]] l> I ]]
i-' Z !NVERCAAG033 l> ]] Z If'oIVERC.a..RG033 l> ]]
[1 TO [1 TO - ]] rn
0 INVERCARGZ20
]] III INVERCAAG220
cT I 3: Z I 3: Z
III a -j III a -j
0 (JJ (JJ Z
....., INVERCARG220
Z :INVERCARG220
H 3:
H 3: TO
TO [1 l>
MANAPOURI220 [1 l> MANAPOURI220
0 (j) (j)
C Z Z
'J H H
'J INVERCARG220 INVEACAF'lC220
-j -j
CD TO TO -
M.4.NAPOURI220 C MANAPOURI220 C
::l 0
cT 0
III III
::r: INVERCA.RG220 INVEI=ICARG220
PJ TO TO -
ROXBUFIGH-220 AOXBUROH-220
'J
a0
::l INVERCARG220 INVERCAAG220
1-'- TO TO -
AOXeUAGH-220 FlOXBU~GH-220
0
Ul
INVERCARG220 INVEF\CA~G220
TO TO -
TII-IA.I----220 TIWAI----220
INVERCARG220 ! g INVERC:AAG220
TO TO -
TlwAI----220 ~ T:IWAI----220
.
"
! "~
MANAPOURI220 MANAPOUAl:220
TO TO
TIWAI----220 TlwAI----220
w
0 3
~
I-
H~
Z-iJ
·rl 2
~ C
<t::J
L
I- L 1
Z QJ
w9::;
II
II
~ 5 9 13 17 21 25
U
ORDER OF HARMONIC
500
400
300
200
100
o
800
-100 R
5
-200
-300
-400
"(j)
E 6121121
..c
0
v
W 4121121
U
Z
-<
0
W 2121121
0..
l:
I-t
121
5121 25121 45121 65121 85121 11215121 125121
FREQUENCY (H",,)
F F F F
T
51.
~
I I
Figure 4.19 Frequently used models
The series branches are handled in the same manner. The more
frequently used filter branches (and subnetworks) have their own separate
model and hence data requirements, thereby allowing a simpler and more
economic data input. Figure 4.19 displays these more frequently used
models.
The program can also plot the response of an RLC network together
with the required response, stored in a file, thereby allowing easy
comparison. This shows the accuracy of a frequency-matched a.c. system
equivalent.
(1) (2 ) (4 ) (5 )
R3 R4
Ll L2 C
L3 L4 s
c C C C C
1 2 3 4 s
I-- -- I I
- -
Figure 4.20 RLC circuit
BOO
600
ADO
200
o ~----~rT--~----~-----+~----~--+--J
400 1200
R
-:200
-400
-600
-800
-1000
Figure 4.21 Impedance loci of RLC circuit
53.
U)
L
I 800
o
UJ
U
z
« 400
o
w
D..
L
H
o
50 450 850 1250
FREQUENCY (Hz)
(a) Impedance versus frequency plot
1200
---UJ
L
I
o 800
UJ
U
z
«
r 400
UJ
H
UJ
W
II
o .J \. .J \... J \...
50 450 850 1250
FREQUENCY (Hz)
800
UJ
L 400
I
o
W 0
U
z
«
r
u -400
«
W
II
-800
50 450 850 1250
FREQUENCY (Hz)
(c) Reactance versus frequency plot
0.50
0.40
lJJ
0 0.30
:J
I-
H
Z 0.20
t!J
«
:r
0.10
1 2 3 4 56789 10 11 12 13
X Variable
Hingorani and Burbery (1970) were the first to propose the use of
an RLC network to represent the frequency-dependence of an a.c. system.
Their paper described the derivation of the network parameters so that
the modelled frequency response matched the required response and
discussed the numerical problems that were encountered. Figure 5.1 shows
the equivalent circuit proposed by Hingorani and Burbery (1970). The
, - - - - - . - - - -..................... ---.-----1
I I
Figure 5.1 A.C. System Equivalent proposed by
Hingorani and Burbery
57.
stages of the circuit synthesis were ; first the required Land C values
are determined based on a lossless model (ie neglecting the resistance in
the model of Figure 5.1) ; resistances are then added which correspond
with the real part of the impedance at the appropriate resonant
frequencies; finally Land C values are scaled so as to give the correct
impedance level at the peaks while maintaining the correct resonant
frequencies,
(i) the system must be left in an operating condition and hence the
numerous outage conditions which are liable to give rise to
resonances at harmonic frequencies cannot be studied,
(ii) the need to make measurements while the system is energized
demands the use of a high power source of harmonics,
(iii) the measurements are restricted to a few harmonic frequencies and
do not give a satisfactory picture of the resonances at
intermediate frequencies,
(iv) the accuracy of measurements made on the actual system are prone
to significant experimental error.
59.
6
5 1
Power-Flow 7
4 1
Transient Convertor
Program Simulator
(PF) (TCS)
8 1
Graphics and Data
manipulation Program
(JAUGPS)
Figure 5.3 Data Flow and sequence for a dynamic study using
a frequency-matched a.c. system equivalent.
0'1
o
61.
14 kV kV
~lanapouri Roxburgh
220 kV 220 kV
90 MW 54 MVAR
220 kV
220 kV
Figure 5.4 The Lower South Island of New Zealand test system
1 axa a
axa a 1
62.
150 15
R(S< )
-50 -5
15 HIO
-5 -50
-250 -250
15 150
-5 -50
Za Zb Zc
250 250 250
RW) RW)
0 0 0 --'---'
and extracting the diagonal terms. This method is valid for systems with
little phase current asymmetry, as is normally the case with convertor
plant. Applying this method to the above test system results in the three
loci displayed in figure 5.6.
o 0 a
50 450 850 1250 50 450 8:50 1250 SO 450 650 1250
I II
(n) (nl (nl
sao 600 600
300
.,0
100
o
00 .4:50 650 12~O
300
.,0
100
a
vV
50 450 SESO
\d'00
300
100
12:50:S0
a
450 B:50 12:10
o a a
!so 450 650 1250 50 450 850 1250 50 450 850 1250
Figure 5.7 Impedance versus Frequency Matrix derived from figure 5.5
0'\
-I::"
65.
I
a +
r-------------~--~ 2:aa~-----<
Zba1a + Zbc1c
Ib
+
2:bb
Z I
ca a + Zcb1b
I
c +
"-
2: cc
1 t
Use Hingorani's topology Decide equivalent
and direct synthesis to circuit topology and
derive a multi-phase optimize the circuit
equivalent parameters
~ ~
Is the equivalent .....
optimize
circuit acceptable
No Yes No
~
·1 Matching Completed I
Figure 5.9 Flow Diagram for the Direct and Optimization Algorithms
(i) The type of study. For steady state harmonic penetration studies
minimal error at the harmonic frequencies is required while the errors at
intermediate frequencies are unimportant. This is clearly unacceptable
for transient convertor simulations as some intermediate frequencies
could be excited.
(ii) The feasibility of altering automatically the equivalent circuits
to reduce the error to a prescribed accuracy.
(iii) The availability of an economically viable alternative to model
the actual system explicitly. This relates to the amount of CPU time
required to provide the simulation and accuracy required.
67.
5 •• 1 Direct Method.
The first step is to obtain values for the inductors and capacitors
of the equivalent circuit based on a lossless model. The resistance of
the branches are ignored as they do not affect the resonant frequencies.
The admittance of an n branch network is
n n
s I {( 1 /L ) II (s 2 + w2 )
n l
s k=l k l=l .. k
yes) I }
n
(5.1)
k=l 2
L (s + W )
2
(s 2 + wi)
2
k k ill1
and s = jw
n
II o (5.2)
i =1
n n
s I {(l/L) II o
k=l l=l;tk
68.
Let F1 ,F
, ... F be the frequencies at which the admittance is
2 n
zero. Then equation (5.3) can be written as :
s o (5.4)
where M.=21fF.
1 1
n
I (1/L )
k
(5.5)
k=1
n n n-1
2 M2
I { (1/L )
k I w.
1
I k
(5.6)
k=1 i=1 k=1
.,k
n n n n-1 n-1
2 2
I ( 1 ILk) I I wiw j } I I M~M~
1 J
(5.7)
k=1 i=1 j =1 i=1 j =1
.,i .,i
Coefficient for S1
n n n
2
I {( 1/Lk ) IT wi} IT M2
k
(5.8)
k=1 j=1 k=1
.,k
Thus n equations are formed with 1/L ,1/L , ...• 1/L as the
1 2 n
unknowns to be solved for. i.e.
[A] x b (5.9)
where
t
x [1 IL1 ' 1/L 2 , .... 1 ILn]
n n n-l
2 2 2
I w.
1
I w.
1
I w.
1
i=2 i=l i=l
*2
A
n n n-l
2 2 2
II w. II w. II w.
1 1 1
i=2 i=l i=l
*2
(5.10)
where k=l .2 •... n
The admittance of the tuned network at any particular frequency is
given by:
n n R - j (Lk/w)
k
(i- 1 I (LkC ) )
k
Y Y I I 2
(5.11)
k=l k k=l R2+ [(Lk/w)(w - 1 I(LkC )) J2
k k
n
w
y - j I (5.12)
k=l
Feature Extraction.
Figure 5.10 displays the phase "B" loci based on equal current for
an equivalent circuit formed using the intersection of real axis feature
extraction (c.f. figure 5.6). The first loop between 260 to 290 Hz is
300
19
150
480 600
R
-150 5
-300
300
19
150
480 600
R
-150
-300
completely missing while the loop consisting of the loci from 500 to 860
Hz does not correspond well with the required loop (displayed in Figure
5.6), as the latter only clips the real axis. Using the new feature
extraction routine the loci of figure 5.11 is obtained. Although the
first loop is too large it can easily be improved by adjusting the
resistance values.
0
50 450 850 1250
FREQUENCY (Hz)
Figure 5.12 Equivalent cct based on new feature extraction
600
A 500
W
L
I 400
av
W 300
U
Z
< 200
0
W
~
L 100
H
0
50 450 850 1250
FREQUENCY (Hz)
trough too low. By increasing the resistance on the first RLC the
trough's minimum value increases and the parallel resonance peaks on
either side diminish, hence the resistance of the RL branch needs to be
decreased to maintain the level of the first peak. Figure 5.13
illustrates the response after these two resistances have been altered.
This gives a better match across the frequency range, but may cause
slight errors in the fundamental frequency (50 Hz) parameters. This
discrepancy, due to the compromises made to get the troughs and peaks of
the correct size, can easily be corrected by the inclusion of a
correction filter, as will be discussed later.
74.
All three performed equally well, although unlike the first two,
the Success/Failure algorithm has the advantage that it does not need two
initial guesses. A discussion of these three optimization techniques is
given in appendix A4. Although the optimization technique determines the
efficiency of solution, it is the objective function's definition that
determines the optimum point and hence ultimately the extent of success
of the optimization technique. Of the objective functions tried the least
squares function has been found to be the most satisfactory. Applying the
least square technique to the total frequency range is undesirable. A
better fit may be achieved by accepting a large initial error in a
particular region and using a correction branch at a later stage. To stop
regions that are going to be corrected by subsequent use of correction
filters detrimentally affecting the scale factor, the frequency range
can be partitioned into up to five regions, with different weighting
factors for each region. Different weightings can also be given depending
on the sign of the error.
75.
300
19
150
480 600
-150
6fZH21
W 3121121
U
Z
< 2121121
0
W
11.
L 1121121
1-1
121
5121 45121 85121 125121
FREQUENCY (Hz)
Z
0
H 81ZlfZllZllZllZllZl
I-
U
Z
::J 61ZllZlfZllZllZllZl
LL
W
>
H
41Zl1Zl1Zl1Zl1Zl1Zl
I-
U
.,
W 21Zl1Zl1Zl1Zl1Zl1Zl
CD
0
IZl
51Zl lfZ1lZl 151Zl 21Zl1Zl 251Zl
SCALE FACTOR
Figure 5.15 Objective function in conjunction with intersection
with x axis feature extraction
z
8 81Zl1Zl1ZllZlfZllZl
I-
U
Z
:J 61ZlfZllZllZlfZllZl
LL
W
>
H
41ZlfZllZllZlfZllZl
I-
U
~ 21ZlfZllZl1Zl01Zl
rn
o
1ZlL---~--~~~==~==~--~--~--~
51Zl llZllZl 151Zl 21Zl1Zl 251Zl
SCALE FACTOR
A 5121121
Ul
2:
I 4121121
0
v
W 3121121
U
Z
-<
0 2121121
W
lL
2: 1121121
H
.. .
'
121
50 45121 85121 1250
FREQUENCY (Hz)
(a) Equivalent based on intersection with x axis feature
extraction
A 5121121
Ul
2:
I 4121121
o
v
W 3121121
U
Z
-<
o 21210
W
lL
2: 11210
H
o
5121 450 850 1250
FREQUENCY (Hz)
(b) Equivalent based on new feature extraction
r. 51Z11Z1
U1
2:
I 41Z11Z1
a
v
W 31Z10
U
Z
«
0 21Z10
W
!l
2: l1Z10
H
..........
0
50 451Z1 851Z1 1250
FREQUENCY (Hz)
(a) Equivalent based on intersection with x axis feature
extraction and scale factor optimized
r. 51Z10
U1
2:
I
av 41Z10
l.JJ 31Z11Z1
U
Z
« 21Z11Z1
0
W
!l
2: llZl0
H
...........
IZI
50 451Z1 851Z1 1251Z1
FREQUENCY (Hz)
r.. 51210
U1
L
I 400
0
v
W 31210
U
Z
-< 21210
0
W
(L
L 100
H
.......
0
50 450 850 125121
FREQUENCY (Hz)
Figure 5.19 Frequency reponse of equivalent when the scale factor
is optimized from 50 to 1250 Hz
60121 r----r---.----~--_r--_,----~--_r--~~.~.--T---~---,--~
·.
r.. 500
U1
L
I 41210
ov ·.
·· ..
W 30121
U
Z
-<
o 2121121
W
(L
L 10121
H
... , .
. . '
121 ~---r----~--~--~----~--~----~--~----~~~--~--~
5121 450 85121 125121
FREQUENCY (Hz)
scale factor. Figure 5.20 shows the frequency response when this new
optimum scale factor is applied.
5121121
'"
(f)
L
I 4121121
0
v
W 3121121
U
Z
«
0 2121121
W
0..
L 1121121
H
121
5121 450 850 125121
FREQUENCY (Hz)
Figure 5.21 Impedance versus frequency
180
'"
(f)
W 120
W
0:: .............. .
(j
W 6121
0 ......
v
W 121
-.l
(j
Z -60 .' .
« .... , ..
W
UJ- 1 20
«
I
0..-180
50 450 850 1250
FREQUENCY (Hz)
Figure 5.22 Phase Angle versus frequency
r, 500
lD
L
I
ov 400
W 300
U
Z
«
I- 200
lD
H
~ 100
(t
0~~__~~~__~__~~==~~~ __ __
~ ~~~'~"~'"
50 450 850 1250
FREQUENCY (Hz)
Figure 5.23 Resistance versus frequency
r,
lD
L 200
I
0
v
W 0
U
z
«
f-
U -200
«
W
(t
-400 ~--~----~--~--~----~--~----~---r----r----r----r---~
50 450 850 1250
FREQUENCY (Hz)
10: v
50
~
450 850
.\j 10: k!
1250 50
~
450 850
\ld
1250
10:
50 450 550 1250
(n) (Il) (Il)
500 500 500
400
::;00
~
f l · ·
Zba 400
::;00·
J 400
::;00
1 1\ Zbc
..
: : lJ \2JJ : :
a
50 450 850
....
1250
0
:: a
50 450 850 1250
: : lJ \2JJ : :
0
50 450 850
. ..
1250
a
50 450 850 1250
Correction Filters.
Y
-cb
Y
-rr
Y
-ec (5.13)
where
Y Required admittance response.
-rr
Y Equivalent circuits admittance.
-ec
Y Admittance required by correction branch.
-cb
2
IZl 12 = R2 + (w1L - 1/(w 1C» (5.14)
2
IZ21 2 = R2 + (w 2L - 1/(w 2 C» (5.15)
2
Iz31 2 = R2 + (w 3 L - 1/(w 3 C» (5.16)
84.
(5.17)
21
1 Z2 - 1Z321 = (2
W 2- 2) L2
W3 + 1 w22 - 1 I"'~)
(I w, IC 2 (5.18)
(5.20)
(5.21 )
Hence
2
IZ212- I Z 31
2 2
(1/w 2 - 1/( 3 )
(5.22)
Illustrative Example.
121
. . ....... .
---'.-:".~
'
A 5121121
U1
2:
I 4121121
0
v
W 3121121
U
Z
« '.
0 2121121 ....
W
(l
2: 11210
H
----:-:-....... ..
121
50 450 850 125121
FREQUENCY (Hz)
Figure 5.26(b) Impedance match of (1,2) mutual element when one
correction branch has been added
61210 r---~----~--~--~----~---r----r---~--~----'---~----,
A 500
U1
2:
I 41210
0
v
W 31210
U
Z
«
0 2121121
W
(l
2: 10121
H
......
121
5121 45121 85121 125121
FREQUENCY (Hz)
Figure 5.26(c) Impedance match of (1,2) mutual element when two
correction branches are used
87.
W 3121121
U
Z
-<
0 2121121
W
0..
2: 1121121
H
0
50 450 850 1250
(a) Phase A FREQUENCY (Hz)
6121121
A 5121121
(J)
2:
I 4121121
0
v
W 3121121
U
Z
-<
0 2121121
W
0..
2: 11210
H
121
50 450 85121 1250
(b) Phase B FREQUENCY (Hz)
600
A 5f210
(f)
2:
I 4f21f21
0
v
W 300
U
Z
-<
0 2f210
W
0..
2: 100
H
0
50 450 850 1250
(c) Phase C
FREQUENCY (Hz)
512l12l
""
Ul
2:
I 412l12l
a
v
W 312l12l
U
Z
« 212l12l
0
W
[L
2: 100
H
........... . ....... . .......
..... . ........... .
0
50 75 112l1Zl 125 1512l
( a) Phase A FREQUENCY (Hz)
"" 51Zl1Zl
Ul
2:
av 41Zl1Zl
I
W 31Zl1Zl
U
Z
« 212.112.1
o
W
[L
2: 112.1!Zl
H ............................. ~
.........................
...... . ......... .
12.1 ~-r--~-r--.--.--.--r--~~--~~--~~--~~--,--,--~~r-~
5!Zl 75 1!Zl0 125 1512.1
(b) Phase B FREQUENCY (Hz)
"" 512.10
Ul
2:
I
a 40!Zl
v
IlJ 3!Zl12.1
U
Z
-< 212.1!Zl
o
W
[L
2: 112.112.1
H
................
..... . ................ .
!Zl +--r--.--r--~~--.-~--~~--,--,--~-,--,-~--~~r--.--r-~
5!Zl 75 101Zl 125 150
(c) Phase C FREQUENCY (Hz)
W 0
.J
(j
Z -60
«
W
(J)-120
<
I
~-180 +--r--r--r--r--r~~'-~--~-.--~-r--r--r~r-~~--~~~
50 75 100 125 150
(a) Phase A FREQUENCY (Hz)
180
1"\
UJ
W 120
W
(t
(j
W 60
0 ' .......................................... , .................................. .
v
W 0
.J
(j
Z -60
«
W
UJ- 1 20
«
I
~-180 +--r--r--r--r-~~--~-.--~-.--~-r--~~~r-~~--~~~
50 75 100 125 150
(b) Phase B FREQUENCY (Hz)
180
1"\
UJ
W 120
W
(t
(j
W 60
0 ..............................................................................
v
W 0
.J
(j
Z -60
«
W
UJ- 120
«
I
~ - 1 80 +---.---,----,---.---.--r---r--r----r---,-----,---,---,.---,----,----,---,--,---,--!
50 75 100 125 150
(c) Phase C FREQUENCY (Hz)
Figure 5.29 Phase Angle Match with a corrective branch for the
50 Hz parameters
91.
(5.23)
(5.24)
Least Squares is the most widely used error criterion. The general
form is :
m 2
Minimize e L (w.G.(x)) }
1 1 -
(5.25)
i=l
92.
where w ,w ' .... w are termed weights or penalties and have the
1 2 m
effect of emphasizing G.(x) in regions of importance.
1 -
Minimize 8 = f W.G.(x)
11-
(5.28)
k k
Sex) - L (A.C.(x))
11 -
+ (1/2) L (cr.
11 -
C~(x)) (5.30)
'1
1= '1
1=
W 300
U
Z
<
0 200
W
[L
2: 100
H
0
50 450 850 1250
FREQUENCY (Hz)
(a) Least Squares
6121121
51Zl1Zl
"
Ul
::L
I 400
0
v
W 300
U
Z
<
0 200
W
[L
::L 11210
H
......
121
. ' . .. , .
50 45121 85121 1250
FREQUENCY (Hz)
(b) Mini-max
6121121
A 50121
UJ
2:
I 40121
0
v
W 3121121
U
Z
<
0 2121121
W
[L
2: 1121121
H
121
5121 45121 85121 125121
FREQUENCY (Hz)
(c) Mini-average
95.
6lZllZl
A 5lZllZl
ill
2:
I 4lZl121
0
v
W 3lZl121
U
Z
«
0 2lZl121
W
[L
2: llZll2l
H
.. " .........
121
5121 45121 85lZl 125121
FREQUENCY (Hz)
(d) Mini-area
squares error criterion, which accounts for its popularity. The various
CPU requirements to reach these solutions from a nominal starting point
are displayed in Table 5.1 . In general the Gradient method takes in the
order of 1 to 4 minutes of CPU time to optimize a typical equivalent
circuit consisting of 11 variables. This compares well with the 20
minutes to 3 hours required for the same optimization problem using the
pattern search.
Estimated
I
I,
I
I I :
I~
h
J ~h~
h - Optimization Algorithm's Step Size
A.C.
system
(5.34)
where I 9. is the current in the RL branch and Ib the phase "B" current.
r
, - - - - - - , - - -..................... - - . , - - - - l
I I
Figure 5.33 Extension of Hingorani and Burbery's equivalent by
incorporating the mutual voltage as part of the
fundamental frequency source
Zb a I a+Zb I
ce
VNlvmSITY OF
100.
r-------,----..........................- - - - - , r - - - l
Figure 5.35 The self circuit for the frequency-matched a.c. system
as implemented in TCS
nuplicit Representation
of Mutual Coupling.
r
11-1
1
II~~
1 1
~I ~ 'V~ ~ T
o
Figure 5.36 Various a.c. system equivalent circuits [\)
103.
The solution process for the Self <-) Mutual circuit interface
is depicted in figure 5.37.
Has No
Convergence
occurred
?
Yes
Standard Torn
TCS Branches
Circuit
Standard Mutual
TCS E
m Circuits
Circuit
(mutual circuits) are used to calculate the mutual voltage at each time
step and do not represent actual branches in the network. The early
attempt used tearing to remove branches and the interface consisted of a
voltage source input for the torn branches and an equivalent current
source input for the main part of the network, as depicted in figure
5.38(a). In the frequency-dependent a.c. system equivalent as
implemented, the mutual voltage is the input to the main part of the
network while the current source the input to the mutual networks, as
shown in figure 5.38(b). In this case the voltage source models the
voltage induced due to a current in the other two phases, while the
mutual network calculates the voltage to be applied, based on the phase
currents and frequency-dependent mutual impedance.
only one frequency at a time was injected. Due to the time considerations
the testing of the complete model by this method, at best, can be
performed only at a few selected frequencies.
50
(1 ,1) A 47.031 46.85
B 47.093 46.81
(2,1) A 15. 122 81.12
B 15.144 81 .09
( 3 , 1) A 12.482 79.21
B 12.488 79.19
150
( 1 , 1) A 160.503 63.10
B 160.941 63.01
(2,1) A 57.962 85.14
B 58.120 85.12
(3,1) A 46.952 84.41
B 47.018 84.40
275
(1 , 1 ) A 474.444 3.86
B 475.135 3.78
(2,1 ) A 499.388 -2.89
B 498.435 -3.50
( 3 , 1) A 397 .566 13.17
B 397.542 12.93
450
(1 , 1 ) A 37.714 -57.18
B 37.575 -57.09
(2,1) A 25.656 -49.13
B 25.470 -48.81
(3, 1) A 21.430 -45.43
B 21.295 -45.13
750
(1 ,1) A 163.500 66.79
B 164.089 66.68
(2,1) A 64.642 45.29
B 64.890 45.32
(3,1) A 49.613 63.65
B 50.183 63.48
1000 Hz
( 1, 1) A 146.361 -6.29
B 145.495 -7.01
(2,1) A 7.202 169.25
B 7.031 162.38
(3,1) A 30.873 -11 .86
B 30.195 -11.15
50
(1 ,2) A 15.122 81. 12
8 15.144 81 .09
(2,2) A 47.023 46.69
B 47.085 46.66
(3,2) A 15 . 122 81 .12
8 15.144 81 .09
150
(1 ,2) A 57.962 85.14
B 58. 119 85.10
(2,2) A 164.500 61 .82
8 164.994 61 .71
(3,2) A 57.962 85.14
8 58.119 85.10
275
(1 ,2) A 499.388 -2.89
8 498.435 -3.50
(2,2) A 507.092 5.89
B 507.685 5.79
(3,2) A 499.388 -2.89
8 498.435 -3.50
450
( 1 ,2) A 25.656 -49.13
8 25.468 -48.80
(2,2) A 41 .059 -59.29
8 41.013 -59.23
(3,2) A 25.656 -49.13
8 25.468 -48.80
750
(1 ,2) A 64.642 45.29
8 64.886 45.33
(2,2) A 152.786 65.97
8 153.404 65.70
(3,2) A 64.642 45.29
8 64.886 45.33
1000 Hz
( 1 ,2) A 7.202 169.25
B 7.032 162.35
(2,2) A 155.870 -4.00
B 152.301 -5.29
(3,2) A 7.202 169.25
8 7.032 162.35
50
(1 ,3) A 12.482 79.21
B 12,489 79.18
(2,3) A 15. 122 81 .12
B 15.144 81 .08
C3 ,3) A 47.031 46.85
B 47.097 46.80
150
( 1 ,3) A 46.952 84.41
B 47.018 84.39
(2,3) A 57.962 85.14
B 58.120 85.11
(3,3) A 160.503 63.10
B 160.942 63.01
275
( 1 ,3) A 397.566 13.17
B 397.541 12.93
(2,3) A 499.388 -2.89
B 498.435 -3.50
C3, 3) A 474.444 3.86
B 475.135 3.78
450
( 1 ,3) A 21.430 -45.43
B 21 .294 -45.13
(2,3) A 25.656 -49.13
B 25.468 -48.81
C3,3) A 37.714 -57.18
B 37.572 -57.09
750
(1 ,3) A 49.613 63.65
B 50.181 63.49
(2,3) A 64.642 45.29
B 64.886 45.33
(3,3) A 163.500 66.79
B 164.089 66.69
1000 Hz
( 1 ,3) A 30.873 -11 .86
B 30.184 -11.14
(2,3) A 7.202 169.25
B 7.032 162.35
(3,3) A 146.361 -6.29
B 1 45.460 -7.01
Consider now the case when two convertors are connected to the
same a.c. system as shown in figure 5.39. This system requires two
points of observation, i.e. busbars A and B. The convertor connected to
busbar A will affect the impedances of the a.c. system when observed from
busbar B and vice versa. This inter-dependence increases the number of
mutual terms as each phase has a self and five mutual terms.
D.C.
System
A.C.
System
D.C.
System
2
(5.31 )
w
I·DCLINK .,
,
~
Initialize frequency Repeat until all section C 1
Repeat for each frequency /
Range and number of data has been entered
sections
Prompt for section Form the data C Repeat for each section Convert from
number, type and file for the ABCD to [y]
relevant data section and reduce to
a 2x2 matrix
I
Read the admittance Form the [y] Write the
and impedance data matrix, convert admittance
from data file to ABCD form and data into a
chain with data file
previous section
-!=
115.
This gives
(5.32)
(5.33)
r.
REQUIRED RESPONSE
SYNTHESIZED CIRCUIT'S RESPONSE
~ 91Z11Z1
I
ov
W 61Z1lZl
U
Z
«
o
W 31Z1lZl
(L
L
H
lZl +----r----.---,---~----~--_r----r_--_r--~----~--~----~
lZl 4lZllZl BIZIIZI 121Z11Z1
FREQUENCY <Hz)
I I
Figure 5.43 Partitioned equivalent circuit
118.
(5.34)
( R22 + X2)
2
Zl Z2 (R l + jX I )(R 2 - jX 2 )
Zt
Zl + Z2) (R l + R2 ) + j (Xl - X2 )
(R 1 R2 + XI X2 ) + j(X 1 R2 - X2Rl )
(5.35)
(R l + R2 ) + j (Xl - X2 )
The impedance magnitude and phase angle are both very important in
harmonic assessment studies. This is due to the time waveforms being
dependent on the magnitude as well as phase angle of the harmonics.
Therefore after the circuit has been synthesized the match in both needs
to be inspected. Both the direct and optimization methods can correct for
regions of inadequate matching of phase angle. The direct method can use
correction filter branches while the optimization method can incorporate
both impedance and phase angle as part of the objective function. The
effect of impedance phase angle on the dynamic behaviour depends on the
cause of the transient and the system and is difficult to predict in
advance.
[ZJ p . =
[A] [B]]
rlm. [ [C] [0]
3.0
a.s
-3.0
6.1 Introduction.
(HARMAC)
I I Calculate
I
Input current Solve
1
Calculate
Input frequency Read shunt Read line
data for all system injection [Ih)=[Y ] [V ] and output
range and capacitors, h h
frequencies harmonic busbar and for all the line
busbars at transformers,
and include impedance for three phase frequencies current
which harmonic filters and
in harmonic a reduced injection to obtain the flows for
impedances are unbalanced
loads at admittance system data. three phase all
required.
matrices. voltages. frequencies.
fundamental
frequency and Output the
matrix. voltages
[\)
U1
126,
(6,1 )
When HVDC convertors are the sources of the harmonic currents then
determining the injected harmonic currents is a complex problem as they
cannot be considered fixed due to the strong interaction which exists
between the system and a large power convertor, HARMAC requires the
harmonic current injections to be specified in order to perform a
harmonic penetration study, Therefore a program external to HARMAC is
required for the simulation of this interaction and thus determination of
the injected current harmonics, Section 6,4 deals with two methods while
section 6,5 compares their results,
h pq ± (6.2)
h pq (6.3)
Since the six pulse Graetz bridge is the basic building block of
an HVDC convertor the pulse number will normally be 6 or 12 or perhaps
some higher multiple of 6. The characteristic harmonics produced by the
three phase bridge are well known and were reported as early as 1945 by
Reid (1945). With the increase in the number and power of convertors the
assumptions used are becoming invalid. The production of non-
characteristic harmonics due to non-ideal conditions have become a
problem. The reasons for non-ideal conditions are :
(a) firing errors.
(b) a.c. voltage unbalance and/or distortion.
(c) Ripple or modulation of d.c. current.
(d) Unbalance commutation inductance.
-1
I
Read the
1
Repeat Until Converged
harmonic
information
for the a.c.
and d.c. For each convertor The a.c. system The convetor The zero crossings
systems. solve for the a.c. voltages at the voltage and d.c. and firing
current injections convertor system model are instants are
Read the and the d.c. terminal are used to calculate updated because of
initial voltage. found for each the d.c. current the a.c. voltage
power flow harmonic. waveform. distortion
conditions
for the
convertor.
w
o
131.
Figure 6,3 shows the data flow and programs used for a complete
harmonic penetration study based on TCS derived assessment of the
convertor interaction. ACREP has been split to show the main functions it
performs and to aid clarity, A corresponding diagram showing data flow
and programs used when the iterative harmonic algorithm is used is given
in figure 6.4,
1
Transmission Line Interacti ve I Harmonic 11 Graphical Display
Parameter Program Data Program Penetration Program of the Harmonics
(TL) (INTER) , "" (HARMAC) (ACREP)
5 Synthesis of
frequency-matched
equivalent (ACREP)
4 10
3
6 71
4
--...;:;....
I Power-Flow
Program
5
;;..
Convertor
Interaction Program
(PF) (IHA)
w
w
134.
E
c
power system and the iterative algorithm uses the distorted voltage V as
t
the commutating voltage. The delay angle is measured from the zero
crossings of the Ec phase-to-phase voltages and 20 degrees was selected
for this study. The other parameters were; Vc= 13 p.u.,
X =X =0.1
1 2
p.u. and I =1.0 p.u . . The solution of equation 6.4 yields
d
19.07 degrees as the commutation angle. This was indeed verified by
measurement of the TCS waveforms.
where X
c
A certain amount of data manipulation was required in order to
compare the TCS and iterative harmonic algorithm phase angle results. The
TCS algorithm uses a cosine wave on phase "A" as its reference while the
iterative algorithm uses a sine wave as its reference. The iterative
algorithm also treats a.c. current harmonics as injections while TCS
treats them as a load current. Therefore to compare results the TCS
waveforms were advanced by 1/4 of a cycle and an FFT performed. The
iterative algorithm expresses the harmonics as the coefficients of a sine
series while the TCS results obtained via ACREP are the coefficients for
an exponential or sine series. As, sin(e+90) = cos(e), 90 degrees is
added to the phase angle of the TCS results to obtain the sine series
form.
The results for the harmonic current injected into the a.c. system
for phase "A" are tabulated in table 6.1. The analytic solution assuming
ideal conditions shows a vast disparity with the other results,
135.
TCS IHA
Harmonic
Order
Magn. Phase Magn. Phase
especially for the higher order harmonics. When the waveshape of the
commutation current is considered in the analytic solution the harmonic
levels obtained are far closer to the simulation results of TCS and the
iterative harmonic algorithm.
E
c
TCS IHA
Harmonic
.order
Magn. Phase Magn. Phase
TCS IHA
Harmonic
.order
Magn. Phase Magn. Phase
Next the Lower South Island (N.Z.) system was used as the
frequency-dependent a.c. test system. Two TCS simulations were performed,
the first using a Thevenin equivalent, the second the frequency-matched
equivalent circuit with implicit incorporation of mutual effects. The
iterative algorithm was then used to assess the convertors interaction
with the a.c. system and compared with the TCS results. The results are
summarized in Tables 6.5 (a),(b)&(c) and tables 6.6 (a),(b)&(c) for a
Thevenin and frequency-matched a.c. equivalent respectively.
138.
Harmoni c Phase A
Order
TCS IHA
1 1.5136 1.5094
3 0.0005 0.0014
5 0.2951 0.2742
7 0.1784 0.1760
9 0.0001 0.0011
11 0.0971 0.0825
13 0.0699 0.0555
15 0.0001 0.0006
17 0.0352 0.0260
19 0.0248 0.0187
21 0.0000 0.0003
23 0.0084 0.0143
25 0.0056 0.0130
0 1.94473 1.94473
2 0.00073 0.00102
4 0.00039 0.00012
6 0.02321 0.02053
8 0.00015 0.00013
10 0.00011 0.00007
12 0.00237 0.00302
14 0.00005 0.00005
16 0.00005 0.00004
18 0.00175 0.00199
20 0.00004 0.00004
22 0.00003 0.00003
24 0.00171 0.00102
Harmonic Phase A
Order
TCS IHA
1 1 .4794 1.4782
3 0.0026 0.0031
5 0.2851 0.2849
7 0.1730 0.1708
9 0.0004 0.0009
11 0.0916 0.0900
13 0.0653 0.0632
15 0.0004 0.0006
17 0.0312 0.029 11
19 0.0209 0.0197
21 0.0003 0.0006
23 0.0068 0.0079
25 0.0054 0.0069
Harmonic DC Current
Order
TCS IHA
0 1.8991 1.8991
2 0.0043 0.0044
4 0.0003 0.0005
6 0.0211 0.0223
8 0.0002 0.0004
10 0.0001 0.0002
12 0.0024 0.0025
14 0.0001 0.0001
16 0.0002 0.0001
18 0.0016 0.0019
20 0.0002 0.0001
22 0.0000 0.0001
24 0.0014 0.0014
Thevenin 106
283 Clast cycle) 11 31 148
Frequency-dependent 183
439 Clas t cycl e) 1903 241
w
o 0.03
::J
f- TIWAI----220
H~
z+l
·rl O. 02
c.9 [
<t::J
2:
wL 0.01
c.9 OJ
<tf!:,
f-
..J 0 . 0 0 ..JLLb.LL.'L.L.}.L.~e:.L.LI...<,;.L..L-L..J,L-.L..'-""""""""'~
o 1 5 9 13 17 21
> BUSNAMES
ORDER OF HARMONIC
Busbar
Harmonic
INVER CARGILL MANAPOURI ROXBURGH
Order
33 kV 220 kV 14 kV a 14 kV b 220 kV 11 kV 1 11 kV g 220 kV
Jlarmonic BRANCHES
Order
A B C D E F G H I
Branch Designation
WINAPOUR1014 TO HlINI\POURI220 A
HIINAPOUR2014 '1'0 /-IIINAPOURI220 B
~IANAPOURI220 TO TlI'iAI----220 C
HANAPOURI220 TO INVERCARG220 D
ItlVERCARG220 'ro ROXBURGII-220a E
INVERCARG220 TO ROXBURGH-220b f'
HNERClIRG220 TO INVERCIIRG033 G
ROXDURGlIl0l1 TO ROXDURGII-220 II
ROXBURGII-220 '1'0 HOXBURGII-Ol1 I
Table 6.9 Harmonic Branch Currents at sending end. throughout the Lower South Island
of New Zealand test system. Based on TCS assessed harmonic currents using
implicit accounting for frequency-dependent mutual coupling.
145.
w
o .08
:::::J
I-
H-:;.06
Z .r!
l'J c
i:::::J .04
I- L
Z ~. 02
W~
0::
0:: . 0 0 .JL:.':;;"";:;'.L!...)L..I.-'-<:::L....:"""-"""'-L....:::.LL""",-,L.:>.L:"::"""'-"'-",,,-,,-
:::::J 1 5 9 13
U
ORDER OF HARMONIC
6.8 Conclusions.
In general, increasing the number of sample points used for the FFT
improves the accuracy. However beyond a certain point the accuracy
deteriorates as the round-off errors and errors associated with
obtaining equally spaced data points, build up. Although the results of
the previous sections were all obtained using 2048 sample points,
subsequent tests have shown that the use of 1024 sample points provides
better results.
146,
7.1 Introduction.
The lower South Island of New Zealand is again used as the test
system and comparisons are made between a Thevenin equivalent of this
system and a frequency-matched equivalent. The filter data used is the
monopolar equivalent of the New Zealand hvdc link filters (Robinson
1966). A simple Thevenin representation is used for the a.c. system at
the receiving end of the d.c. link (system 2) in sections 7.2 and 7.3.
Finally, in section 7.4, the detailed a.c. representation and the a.C.
fault location are placed at the receiving end (system 2) in order to
assess the effect of the disturbance at the invertor end.
Transmission Line
III,
I JlI
A.C.
System
Figures 7.2 (a)&(b) show the rectifier current and d.c. voltage on
the convertor side of the rectifier's smoothing reactor, when a d.c. line
fault is applied midway between the two convertors. The control strategy
on the rectifier is such that on detection of the fault the convertor is
7.5
;J.~~
7.121
6.5
~!.~
--._- FREQUENCY-MATCHED EQUIVALENT No. 3
5.121 ----- FREQUENCY-MATCHED EQUIVALENT No. 4
4.121
3.5
~
~\
3.121
2.5
2.121 .~
1.5
1.121
~
121.5
121. 121
~
-121.5
51214121.121 522121.121 54121121. 121 558121.121 576121.121
T I ME (DEGREES)
4.5r-------.-------~--------r_------,_------~--------r_------._------_,
1\
3. 121
-;0;~12I-.~12I---5~2~2~12I~.~12I~~5~4~12I~12I~.-12I--~5~5~8~12I-.~12I---5~7~6~12I-.-I2I----5-9-4~12I-.-I21----6-1-2LI2I-.I21----6-3-12I~12I-.-I2I----6-4~8121.121
T I ME (DEGREES)
A.C.
System
The d.c. voltage waveforms show (figure 7.2 (b)) greater variation
between the various alternative models. This is due to its sensitivity to
the trapped energy in the d.c. line once blocking of the convertors has
occurred. The d.c. voltage is primarily determined by the inductor
voltage and hence dependent on the rate of change of current.
150.
A.C. system
CUrrent
A.C.
System
10
-5
-10
-15
-20
-25
-30L-__________
5040.0
~~~~--------~~~~--------~~~----------~~
5220.0 5400.0 5580.0 5760.0
T I ME (DEGREES)
(a) A.C. System Current
1.4
1.2
1.0
0. 6
0.6
0. 4
0.2
.J
0.0
! \ '\
-0.2
~ ,
v\ .
-0. 4
, ~\(\ / .\ I
-0.6
1/ . J ~ .
-0.8
'.,..1\
..- \1\1 I'j .
\/-V
~
-1. 0
-1. 2
-1.4L-__________~~~~--------~~~~--------~~~----------~~[
5040.0 5220.0 5400.0 5580.0 5760.0
T I ME (DEGREES)
a.c. system equivalent also models the coupling between phases and thus
any distortion present in the two unfaulted phases is reflected into the
faulted phase. The fault effectively shorts out the harmonic filters on
the faulted phase and thus the dynamic response is that of the a.c.
system following the sudden application of a short circuit.
301r----------------,----------------1I----------------r---------------~
25
20
15
10
\//
-1121 ."-..,.. X.... _____ / /
-15
-2121
-25
-301L-__________~~~~--------~~~~--------~~~~--------~~
512140.121 522121.121 54121121.121 558121.121 576121.121
T I ME <DEGREES)
~
1.2
1.121
121.8
121.6
\
121.4 \
121.2
\
121. 121
-121.
\
\
-121.8 \ .~.
/'
-1. 121 -V
-1. 2
-1.4L-____________~~~----------~~~----------~~~~--------~~
512140. IZI 522121.121 54121121.121 SSB0.121 S76121.121
TIME <DEGREES)
25
2121
15
1121
,
5 \
\
-5
'\ \)(.-r--~
-}121 '----..../ '--- - /
-15
-2121
-25
-3121L-____________~~~~----------~==~----------~~=_=_----------~~
51214121.121 522121.121 54121121.121 558121.121 576121.121
T I ME (DEGREES)
1.4
1.2
1.121
121. 8
121.6
121. 4
121.2
J
121.121
-121.2
I I
I
-121. 4
-121.6
\
'V"'" ,../
I
I \ ~. )
-121.8
-1. 121
" \."-J.\-1
-I. 2
-1.4L-______________~______________=_~-----------------L-------------~~
51214121.121 522121.121 54121121.121 558121.121 576121.121
T I ME (DEGREES)
30'~-----------r------------,------------r----------~
25
20
15
I /
''I' /
, ~--
-15
-20
-25
-30~------------~~~~----------~5~4~0~0~.~0~----------!5~5518~0a..~0;------------;577660.0
5040.0 5220.0 TIME <DEGREES)
/ \
-0.8
\
--..,./
/
/
\ /
-1.0
\ Ii
-I. 2 \ .r
-1.4L-----------------L-------------~5~4~0~0~.~0~----------~5~5~8~0;.~0;------------5~776~0.0
5040.0 5220.0 TIME <DEGREES)
A.C. System
Current
A.C.
System
2
1.2
- - - Phase A
1. '" _
"'.9
\
-"'.4 \\ /
-0.6
-121.8
\. ).~
-1. 121
\. J .
\.
-1.4L-____________
51214121_ t2I
~~L--------------~~~------------~~~------------~
S22121~
121 5400. 121 5580. 0 576121. 12)
TIME <DEGREES)
1.2
1.121
(\
" ,
J \
/ /\ \
. 1'"'/
I
\
.
.~.
\ \
~ l\, I \
\
-1. 2
\)
-1.4L-______________-L____________~~~~----------~~~=_----------~~
5040. I2l 5220.0 5400.0 5580.121 576121.121
T I ME (DEGREES)
1.121
1'.
121.8
.'\
/ .\ / /\ ,....F' \
.
\
-0. 2
-121. 4
-121.6
-0.8
-1. 0
-1. 2
-1.4L-______________ ____________
~ ~~~~------------~--------------~_:".
16
12
.r--" ____... ~'\
/
8
'
/I/\,,_V\,
-8
-12
-16
-20
-24L-________________L-________________ ~ ______________~~~------------~=
5040.0 5220.0 5400.0 5580.0 5760.0
TIME (DEGREES)
(a) Thevenin Representation
24~----------------.-----------------._----------------.-----------------,
20
16
12
-16
-20
20
35
30
25
20
15
10
0
-5
-10
-15
-20
-25
-30
-35
-40
5040.0 5220.0
T I ME (DEGREES)
response of the circuit that models the frequency-response and hence the
resonances of the a.c. system. The remaining two phases experience
enormous distortion due to the upset in the valve conduction pattern.
When a frequency-matched e·quivalent with mutual coupling is used (figure
7.11 (c)) this distortion is reflected into phase "A" causing the faulted
phase current to be grossly distorted.
CHAPTER 8 CONCLUSIONS
With the rapid growth in both size and number of HVDC schemes there
is a greater need for the accurate simulation of transient phenomena,
particularly in cases where the d.c. convertor station feeds into a
relatively weak a.c. system. Although explicit modelling of essential
a.c. components (e.g. convertor transformer, harmonic filters) has been
retained, due to the size and complexity of modern power systems some
form of equivalent circuit must be used for the remainder of the a.c.
system. Towards this end the main aim of this thesis has been to present,
and show the effect of, a three phase a.c. system model that accurately
represents the a.c. system at harmonic frequencies. Both the
frequency-dependence of the self and mutual terms associated with each
phase have been accurately modelled. Several frequency-matched a.c.
system models have been presented with any inherent assumptions outlined.
Both a.c. and d.c. system disturbances have been presented and the
response of the proposed models show considerable difference with the
Thevenin representation, particularly in the case of an asymmetrical
fault.
The use of a Thevenin model based on the S.C.R. of the a.c. system
only represents the system accurately at fundamental frequency and
therefore gives inaccurate results at harmonic frequencies. The
resonances between the a.c. system and harmonic filters as well as
resonances within the a.c. system must be modelled for accurate dynamic
simulation. It is the combined impedance response of the a.c. system in
parallel with harmonic filters that needs to be accurately modelled. The
resonances in the parallel combination may be excited by disturbances,
161.
With reference to the New Zealand HVDC link the initially installed
th th th
filters were tuned to the 5 , 7 , 11th & 13 . The N.Z. Electricity
Department did this as some HVDC authorities recommended a "wait and see"
policy about harmonic problems, for there is no justification in
providing expensive filter branches to guard against trouble that might
not occur. However, early operating experience showed a need for a
th
high-pass filter and 9 harmonic filters to reduce the wide spread
th
telephone interference that occured. The 9 harmonic problem was a
result of a resonance between the filters and a.c. system (Robinson
1966). The development of the ACREP program described in this thesis
allows the a.c. system, filter and combined reponse to be inspected and
any resonance problems forecasted. This allows remedies to be found
before the troubles occur. Also the combined response of "a.c. system
equivalent" and filters can be compared to the combined response of
"actual" a.c. system and filters in order to assess the sophistication
needed in the a.c. system equivalent for accurate steady state and
dynamic studies.
Future Work.
Although the modelling of d.c. line faults has been carried out in
detail, further work is needed for the a.c. fault Simulation, in
particular the arc characteristics of the a.c. fault and more flexible
circuit breaker models are required.
the Lower South Island of New Zealand where Tiwai and Benmore busbars
both have large a.c./d.c. convertors connected.
REFERENCES
Al-Khashali (1976)
"Generalised Dynamic Modelling of High Voltage AC-DC Transmission
Systems"
Ph.D. Thesis, University of Manchester Institute of Science and
Technology, Great Britian.
Ametani A. (1976)
"A highly efficient method for calculating transmission line
transients"
IEEE Trans., Vol. PAS-95, No.5, pp 1545-1551
Arrillaga J., Arnold C.P., Heffernan M.D. and Campos-Barros J.G. (1977a)
"Stability of Isolated Generator-HVDC Convertor Units"
IEEE/PES Summer Power Meeting, Mexico City, Paper A77-508-5
Arrillaga J., Bradley D.A., Gellen S.J. and Densem T.J. (1985b)
"Computer Graphics as a Tool for Understanding the Harmonic
Behaviour of Power Systems"
Journal of Electrical and Electronic Engineering, Australia,
Vol. 5. No.1, pp 1-7
166.
Arrillaga J., Watson N.R., Eggleston J.F. and Callaghan C.D. (1987)
"Comparison of steady-state and dynamic models for the calculation
of ACIDC System Harmonics"
lEE Proc., Vol. 134, Pt. C, pp 31-37
Budner A. (1970)
"Introduction of Frequency-Dependent Line Parameters into an
Electromagnetic Transients Program"
IEEE Trans., Vol. PAS-80, No.1, pp 88-95
Breuer G.D., Chow J.H., Gentile T.J., Lindh C.B., Numrich F.H.,
Lasseter R.H., Addis G. and Vithayathil J.J. (1982)
"HVDC - AC Harmonic Interaction.
Part II - AC System Harmonic Model with Comparison of Calculated
and Measured Data tl
IEEE Trans., Vol. PAS-101, No.3, pp 709-718
Concordia C. (1956)
"The Transient Network Analyzer for Electric Power System
Problems"
Supplement to CIGRE Committee No. 13 Report for 1956.
Cory (1965)
"High Voltage Direct Current Convertors and Systems"
Macdonald & Co. Ltd.
DEC (1984)
"Programming in VAX FORTRAN" (Version 4.0)
Digital Equipment Corporation, Maynard, Massachusetts.
Gear C. W. (1 971 )
"Numerical Initial Value Problems in Ordinary Differential
Equations"
Prentice-Hall, Inc., New Jersey.
Greenwood A. (1971)
"Electric Transients in Power Systems"
J. Wiley & Sons, Inc.
170.
Gunn H. R. (1966)
"Commissioning and Early Operating Experience with the Benmore-
Haywards H.V.D.C. Transmission Scheme"
New Zealand Engineering, Vol. 21. No.3. pp 93-101
Harwell (1986)
"Harwell Subroutine Specification, VF04"
Computer Science and Systems Division.
Atomic Energy Research Establishment,
Harwell, Oxfordshire, England.
IEEE (1983)
"Power System Harmonics : An Overview"
IEEE Working Group on Power System Harmonics,
IEEE Trans., Vol. PAS-102, No.8, pp 2455-2460
172.
IEEE (1984)
"Bibliography of Power System Harmonics", Parts I and II,
IEEE Trans., Vol. PAS-103, No.9, pp 2460-2479
Joosten A. P. B. (1987)
"Convertor Transformer Simulation for T.C.S. Studies"
Ph.D. Thesis to be submitted, University of Canterbury,
New Zealand.
Kohler A. (1967)
"Earth Fault Clearing on an HVDC Transmission Line, with Special
Consideration of the Properties of the DC Arc in Free Air"
IEEE Trans., Vol. PAS-86, No.3, pp 298-304
Kron G. (1939)
"Tensor Analysis of Networks"
J. Wiley
Kron G. (1959)
"Tensors for Circuits" (2nd Edition)
Dover Publications.
Kron G. (1963)
"Diakoptics - the Piecewise Solution of Large Scale Systems"
Macdonald, London.
Marti L. (1982)
"Accurate Modelling of Frequency-dependent Transmission Lines in
Electromagnetic Transient Simulations"
IEEE Trans., Vol. PAS-101, No.1, pp 147-155
Marti L. (1983)
"Low-Order Approximation of Transmission Line Parameters for
Frequency-Dependent Models"
IEEE Trans., Vol. PAS-102, No. 11, pp 3582-3589
Pender J. T. (1 969)
"Fast Transients in Electric Power Systems"
Int. J. elect. Engng Educ, Vol. 7, pp 419-429
Reeve J. (1 970)
Discussion of Hingorani and Burbery (1970) paper.
IEEE Trans., Vol. PAS-89, No. 5/6, pp 827
Rumpf E. (1 980 )
"The Operational Performance of HVDC Systems throughout the world
during 1975~1978"
Uhlmann E. (1975)
"Power Transmission by Direct Current"
Springer-Verlag Berlin/ Heidelberg
UMIST (1981)
"Harmonics in Power Systems"
An International Conference on Harmonics in Power Systems,
UMIST, Manchester.
APPENDIX A1
L
g [
~ss ~srJ
rs rr
where
-L
aa -L aa 200S(28» -L ao -L ao 2008(2(8+150)
-L OC -L 00 200S(2(8+120»
LakdCOS(e) -L a kq sineS) ]
L Lbkdcos(S-120) -L sin(S-120)
sr bkq
Lckdcos(S+120) -L c k q sin(S+120)
L Lt
rs sr
Lfd Lfkd 0
L Lfkd Lkd 0
rr
0 0 L
kq
[
d L
d L
Ci8g
aEi ss a-e
aL ]
sr
d L aL
a-e rs a-e rr
a Lrr a L d L }t
where
aS
0 and
a-e rs ae sr
180.
APPENDIX A2
(S1/km)
where
R .. - a.c. resistance of conductor i (S1/km)
11
GMR.- geometric mean radius of conductor i (S1/km)
1
f - frequency in Hz,
w = 21ff
~Rii' ~Xii - are Carson's correction terms for earth return effects
Note for the self impedance terms that the geometric mean radius
(GMR) instead of the actual radius r is used to account for the
contribution which the internal inductance makes to the total inductance.
The series impedance matrix is also symmetric as Zik= Zki
L Lab L L
aa ac ae
Lba Lbb Lbc Lbe
[L]
L LCb L L
ca cc ce
L Leb L L
ea ec ee
182.
a at
Lab (
Lbb
b b'
Lee
c et
~ f..-
C
Lee s
C
1
Appendix A3
R L C C
s
(rt) ( H) ( llF) (llF)
Appendix A4
Let x and x be two x values that bracket the actual minima. Then
1 2
the method that successively halves the interval known to contain the
minima by placing two test points close together at the centre of the
interval, is known as a Dichotomous Search. The dichotomous search is
illustrated in figure A4.1 . A pair of test pOints are required at each
step in order to be able to distinguish which half of the interval being
divided contains the minima. The small difference ~x between the
evaluations of each test point is limited by the accuracy in which the
objective function can be evaluated and sets a limit on the resolution of
the optimization.
--
lk ----------~
I
I
I L
~I~~-----------Ik+l
I
R I
I~~~--------+--Ik+l ----------~)~I
: L
I ~ Ik+2 I >- !
I! I
I I R '
,
1-1.«.«---, Ik+2 ---.;>--1
I ,
I I I
I L ' :
~Ik+3~
I R
r--:
,,
Ik+3~
I
,I ,,
I
'L I I ,R >-
x xa,k I xb,k x
k , I k
L R
x x b ,k+l xk+l
,xk+l I
a,k+l
I I
I I
L R
x k +2 x a,k+2 x ,k+2
b xk+2
(A4 . 1 )
Note that ~t the final step the interval In+l is created equal in
length to I n but is not used. However, it is needed in equation (A4.1) to
determine I n- l' The ratio of the kth interval length to the final
interval length then forms the well known Fibonacci sequence.
I n+l I
n
I 21
n-l n
I 31
n-2 n
I 51
n-3 n
I = 81 etc ...
n-4 n
The problem with the previous methods is the need to specify two
points that bracket the true minima. The Step Success/Failure Search,
like the dichotomous search, uses two closely positioned test points.
However, only one initial starting point is required. The algorithm then
steps in the direction of the minima and evaluates a new pair of test
points. Based on these results, the algorithm either: steps again in the
same direction; halves the step length and steps in the new direction
indicated; or remains at the same point and reduce the distance between
the two test pOints. A flow diagram for the Step Success/Failure search
is displayed in Figure A4.3.
187.
Do
these two tes t Ves
poInts IndIcate the0-____________________-,
minIma lies
at higher
vaiue
?
Do Was
the last No
step taken In this
dIrectIon
?
Was No
the last
step taken In this step In direction
direction of Minima
?
lIave
too many Yes
halvlngs of test
point spacing
occurred
?
Appendix A5
PUblished Papers.
189.
394
Given the size and complexity of a power As an illustration, Fig. 1 shows the New
system it is not practical to model each of Zealand South Island transmission system in
its components individually. Equivalent the vicinity of an aluminium smelter (at
circuit models are usually derived to suit Tiwai) and the three phases of the
the fundamental frequency, the simplest being corresponding impedance locus is illustrated
the short circuit impedances used in fault in Figure 2.
studies. Such models, however, are
inadequate to represent the system behaviour Derivation of Harmonic Equivalents
at harmonic frequencies.
In an earlier model (1) the selected minima
AC networks are generally inductive for the and maxima frequencies used for the
lower harmonics. Since the convertor filters derivation of the equivalent circuit are
are often capacitive for such frequencies, a those where the impedance locus crosses the
parallel resonant condition may exist, real axis (refer to Figure 2). However the
typically at around the third or fourth impedance versus frequency plots derived from
harmonic, or even at higher frequencies the harmonic impedance locus provide more
depending on the a.c. system strength. explicit information of the maxima and minima
~Ioreover a complex system/f i 1 ter combination harmonic impedances.
will produce several oscillatory harmonic
frequencies of relatively low orders. The equivalent circuit will consist of a
Disturbances may excite these natural number of single-tuned parallel branches as
frequencies and produce considerable waveform shown in figure 3 and the first step is to
distortion and overvoltages which will in obtain values for the inductors and
turn interact with the convertor control. capacitors of the equivalent circuit that
will produce the same minima and maxima
Symmetrical impedance loci are currently resonance frequencies. The resistances of
measured or calculated for the complete the branches are ignored as they have very
harmonic spectrum at various operating little influence on the magnitude of the
conditions. These can be used to derive admittance at frequencies that are
accurate harmonic equivalents (1)(2) for the sufficiently higher or lower than the series
a.C. system. However it has been shown (3) resonant frequency of the branch.
that there is considerable diversity between
the resonant frequencies seen by different The admittance of an n branch network is;
phases and thus in general the use of
three-phase loci and corresponding n n
independent phase equivalent circuits will be
n
s L .l.. Ti(s'+w')
more appropriate. I s
k=l Lk J;i k J
Y (s) (1)
k=l Lk(s'+w~) n
This paper describes the derivation of TI(s'+w')
three-phase harmonic equivalent circuits to 1=1 I
represent a.c. networks and illustrates their
effect on the dynamic behaviour of hvdc where:
convertors.
1
w'I LIC 1=1,2 ... etc and s = jw
HARMONIC MODELLING OF THE A.C. NETWORK r
For an LC network the impedance is zero at
System Harmonic Impedances minima and infinite at maxima frequencies.
395
n 1 n n n
s L -
k=l Lk
n
J~
(s '+w ') = 0
J
(3 ) Y
jk!lL k (W 2 =r/L k C )
k
-j \"
L L
k=l k
(",2'" _",2 )
k
(11 )
Let 1'1' 1'2' ... F n be the frequencies at which As Lk appears in the denominator of equation
the admittance is zero, then equation (3) can (11), to increase the admittance contribution
be written as: of each branch, so as to reduce the
impedance, the inductors are divided by the
scaling factor. In order to keep the same
s n
n-l
k=l
(s'+M k) where Mk = 21lFk (4 ) resonance frequencies the capacitor values
are multiplied by the scaling factor.
By equating the coefficients of s in (3) and These equivalent circuits, developed in the
(4) the following equations result: frequency domain, can also be used in the
time domain according to standard circuit
For s2n-l laws. When using transfer functions in the
frequency domain, multiplication by the
n Laplace transform of the input and
1
L 1 (5 ) application of the inverse Laplace transform
k=l Lk obtains the time solution.
396
waveform information from the TCS programme 2. Van Que Do and Gavrilovic, M.M. 'An
without altering the basic algorithm. iterative pole-removal method for
synthesis of power system equivalent
Illustrative Example networks', IEEE Trans. Power apparatus
and systems, Vol. PAS-103, No.8, August
Figure 6 illustrates the test system used to 1984, pp. 2065-70.
demonstrate the effect of alternative a.c.
system representations during a d.c. short- 3. Densem, T., Badger, P.S. and Arrillaga,
circuit. J. 'Three-phase transmission system
modelling for harmonic pentration
The a.c. system on the rectifier side is studies', Trans. IEEE, PAS-103, 1984, pp.
first represented as a Thevenin equivalent to 310-17.
give the specified Short Circuit Ratio. The
dynamic response prior to, during, and after 4. Arrillaga, J., Densem, T. and Bodger,
fault clearance is illustrated in Figures 7 P. S. 'Zero sequence harmonic current
and 8 (continuous line). generation in transmission lines
connected to large convertor plant',
The dynamic study is then repeated with the Trans. IEEE, PAS-I02, 1983, pp. 2357-63.
a.c. system represented by the tuned
equivalent circuit developed above with all 5. Heffernan, M.D. 'Analysis of AC-DC
other parameters, controls and specifications system disturbances', Ph.D. Thesis,
remaining unchanged. The results are plotted University of Canterbury, School of
in the same figures (dash and dot). Engineering, Electrical Engineering
Department, 1980.
It is not the purpose of this paper to
discuss the fault behaviour in any detail,
but rather to show that the predicted
behaviour can be inaccurate in the absence of
a frequency-dependent a.c. system model.
CONCLUSIONS
REFERENCES
397
--r-r-~~r-ROXBURGH
nWAI 90 MW 54 MVAR
X(rI)
25121
2121121
150
1121121
5121
121
5121121
-5121
/
-11210 / /
/ (/
-150
-2121121
-25121
.
----
~
--._._.--- ./
//
.
398
R (12)
500
400
300
200
100
....
a
50 250 450 650 850 1050 1250
(f)Hz
X(I1)
300
(ii )
200
.....
100
f(Hz)
a
50 250
-100
-200
-300
399
ld (per unit)
10
o
-1
6 6.25 6.5 6.75 7.25 7.5
cycles
Vd (per unit)
7
Figure 8 Variation of the voltage at the conVertor d.c. terminals following a line short-circuit
(i) Improved equivalent model
(ii) Thevenin equivalent
195.
87 WM 171-2
FREQUENCY-DEPENDENT A.C. SYSTEH EQUIVALENTS FOR HARHONIC STUDIES AND TRANSIENT CONVERTOR SIHULATION
1
© 1987 IEEE
196.
150 15
RIll )
-50 -5
-250
15 150
-5 -50
150
-50
Za Zb Zc
250 250 250
Rln) R (Q)
0 0 0
WI (1/ I (1/1
aoo
Zac
.00
000
'00
12f10 00
Zbb
'00
eoo
'00 eoo
Zca
'250
Self circuits
+
functions of frequency).
and h is the Error Criterion Function
Among the various types of Error Criterion, the e"",
Least Squares method has been found to be the most
suitable in our problem. Its general form is: 5"'..
minimize 0 = Sum{(wi·g(~))2) ..."
(a)
where wI,---w m are weighting factors to emphasize 3"'"
9i{~) in regions of importance.
2 ....
The optimization method is more flexible in the
selection of equivalent circuit topology; the only ",,.,
limitation in this respect is the possibility of
evaluating its frequency response. However, the m f
circuit topology must be decided before the 50 250 ."0 6"01 SS0 112150 1250
optimization process.
The CPU time requirements of the multi-variable
optimization techniques is several orders of magnitude aD ..
greater than those of the direct approach but they can
be substantially reduced if the circuit derived by the 5""
direct method is used as a first approximation in the
optimization approach. ..." (b)
EFFEC~'IVENESS OF THE IILGORITHH 3""
2 ....
Figure 8(a) displays the frequency response of
one of the mutual elements of the 'actual' test system
1 ....
of Figure 1, together with that of Hingorani's [21
harmonic equivalent. Although the results are in good
agreement at low frequencies the second and third "5111 2"0 .5" 65111 85" lUSt!! 12513
peaks (at 845 and 876 Hz) show substantial disparity.
The use of a correction filter in the region of
largest discrepancy (862 to 906 Hz) modifies the
harmonic response as shown in Figure 8(b). The new 6""
branch reduces considerably the error in the regions
where the impedances of the equivalent circuit were 5""
too high, but no marked improvement is noticeable
where the equivalent circuit impedances are too low. 4""
As explained earlier these are best corrected by a (e)
combination of scaling factor and subsequent use of 3""
correction filters.
The region of greatest discrepancy of the new 2""
equivalent circuit is now between 750-810 Hz. The
addition of a second correction filter in this region '"''
produces a far better match as shown in Figure 8(c).
With reference to multi-variable optimization "5"
techniques, the effectiveness of the synthesised
equivalent circuit is very dependent on the choice of
error criteria. If the optimization criterion is the
minimization of the maximum deviation, the response, Fig. 8. Direct Synthesis of a Mutual Element of the
illustrated in Figure 9(a), is unsatisfactory, i.e. a Impedance/Frequency Hatrix
large discrepancy is introduced in the first peak. In Actual system response (continuous line)
this case it is necessary to use more elaborate Initial equivalent (dotted line (a))
weighting functions which are a function of the With one additional correction filter
discrepancy as well as the frequency. On the other (dotted line (b))
hand the use of the least squares as the error With two additional correctioll filters
criterion, illustrated in Figure 9(b), shows a marked (dotted line (e))
improvement. Although the sum of the errors squared
is minimized in this case, there are larger In order to illustrate the effectiveness of the
discrepancies at Some individual frequencies when frequency-matched equivalents described above, the
compared with the direct solution. test system of Figure 1 Was modified to include an
Thus an interactive approach based on the direct h.v.d.c. convertor at the TIWAI Bus, Figure 10 shows
method enhanced by the use of correction filters the d.c. voltage waveform, derived from Transient
generally has les& error over most of the frequency Convertor Simulation, when a short-circuit occurs on
range as compared to an optimization derived the d.c. line. The continuous line shows the response
equivalent. However, each correction filter increases with a freqUency-matched a.c. system equivalent and
the complexity of the equivalent circuit and hence the the dotted line that of a Thevenin equivalent derived
computational burden when carrying out time domain from the short-circuit capacity of the a.c. system.
studies. On the other hand the CPU time required for There is a significant difference between the two
a multi-variable optimization study is several orders responses which shows the importance of using accurate
of magnitude greater than the direct approach. Hence frequency-dependent equivalents for the a.c. system.
if a relatively small number of time domain The time domain simulation algorithm includes
simulations are required the direct. approach is accurate modelling of the non-linear characteristics
preferable while for many, the optimization technique of the convertor plant, such as convertor transformer
is probably best as a comparable accuracy equivalent saturation, surge arresters, etc.
can be obtained with fewer circuit elements.
199.
---1 optimize
I
I
No
branch to
Fig. 6. Two-convertor System circuit
4
200.
5"
REFERENCES
8,m
[ 1) Kitchin, R.II. "New mothod for digital comlluter
5"''' evaluation of convertor harmonics in power
4<111J
systems using state variable analysis n , Proc.
lEE, Pt. C, Vol. 128, No.4, pp. 196-207.
3111111 (b) [ 2) }lingorani, N.G. & Burhery, M.P. "Simulation of
A.C. system impcdnnce in IIVDC system studies",
2111" IEEF. Trans. Vol. PAS-89, No. 5/6, May/Juno 1970,
pp. 820-826. .
1"111 [ 3) watson, N.R., Arrillaga, J. & Joosten, A.P.B.
"A.C. system equivalents for the dynamic
I1J simulation of HVDC convertors", IER Conference
5" Publication 255 on AC and DC Power Transmission,
Sept. 1985, London, pp. 366-371.
Fig. 9. Optim~zed Mutual Element of the Impedance/ [ 4) Arrillaga, J., Arnold, C.P. llarker, B.J.
Frequency Hatrix "Computer Modelling of Electrical Power
Actual system response (continuous line) Systems", J. Wiley (1983).
Using maximum deviation criterion (dotted [ 5) Adby, P.R. and Dempster, M.A.H. "Introduction
line (a)) to optimization Methods", Chapman and Hall
Using least squares criterion (dotted line (1974).
(b) ) [6) Fletcher, R. Practical Methods of Optimization,
Vol. 1 Unconstrained Optimization, John Wiley &
Sons (1980). Vol. 2 Constrained Optimization,
John Wiley & Sons (1981).
APPENDICES
requirements.
It should be noted that the equivalent circuit
does not give explicit harmonic information at the Impedance maxima will occur when the numerator
6
201.
n n
y )'
k=1
-w
Lk (wLI/L C )
k k
-j r
k=l
W
Lk(w2-w~)
( 11)
L
k=1 Lk J;ik
.l.. n (s2+wJ) 0 (3 )
As Lk appears in the denominator of equation
(11) ,to increase the admittance contribution of each
Let FI' F2, ... F n be the frequencies at which the branch, so as to reduce the impedance, the inductors
admittance is zero, then equation ()) can be written are divided by the scaling factor. In order to keep
as: the same resonance frequencies the capacitor values
n-l are multiplied by the scaling factor.
5 f! (5 2HI~) where Mk (4 )
k=l B. Derivation of Correction Circuit Parameters
2
By equating the coefficients of 5 in (3) and (4) 1 zl1
2 R2 + (wl L - _1_) (12 )
wJC
the following equations result:
n n n-l
)' .l.. IT w~ iT (8 ) ( 17)
k=l Lk J=k k=l
Thus n equations are formed which are then substituting equation (17) into equation (15)
solved for the n unknown inductances (Ll,L2, ... Ln l. gives)
The appropriate capacitor values are found using:
n n Rk - ::x
w
(w 2
1
- LC)
k k
y = I. Yk I. (10)
k=l k=l
R~ + [.~ (w 2 - LC)]1 2
w
k k
Hence
7
202.
the frequency components of the phase-current wave- As a result of initial mismatches between the real and
forms. In conventional harmonic analysis these are the assumed voltages and currents throughout the system,
harmonic currents used for the derivation of harmonic the dynamic simulation must be run for a number of
voltages throughout the (otherwise linear) power system. cycles until the solution reaches a steady state.
Each of the harmonic currents is injected in turn into The structure of the TCS algorithm is described in
the AC system to derive the individual levels of harmonic Appendix 12.2 and a more detailed formulation and the
voltage distortion. numerical solution can be found in Reference 1.
Finally the steady-state time-domain current wave-
3 Iterative harmonic analysis (IHA) forms are subject to a fast Fourier transform to obtain
their harmonic components.
The direct solution described in the previous section
assumes that the convertor-current waveform derived 5 Verification of the algorithms
from the single-frequency power-flow solution is not
affected by voltage distortion. However, the approx- The IHA and TCS programmes have been developed
imations used in the power-flow solution can introduce independently using different conventions, per-unit
considerable error in the derivation of the harmonic cur- systems etc. It is thus essential to validate their per-
rents. formance before attempting to compare their relative
On completion of the three-phase AC/DC power-flow merits. For a reliable verification, the results of IHA and
solution, which provides reasonable initial information TCS must be compared with an exact solution.
for the harmonic analysis, the approximations used can The AC/DC power-flow algorithm provides an exact
then be removed and an iterative solution carried out solution when the system under consideration satisfies
involving the following steps: the approximations set out in Section 2.1. The idealised
(i) Input information is read consisting of AC and DC circuit diagram of Fig. 1 represents such a system which
system harmonic impedances as well as (initial) convertor includes a six-pulse bridge rectifier.
operating conditions derived from a three-phase AC/DC
"~~L
load flow.
(ii) The alternating (harmonic) current injections and
the direct (harmonic) voltage are calculated.
(iii) For each harmonic, the AC system voltages at the
convertor terminals are found.
(iv) The convertor direct voltage and DC system har-
monic impedances are used to calculate the direct-current 1.73p.u.
waveform.
Fig. 1 Basic equivalent AC/DC sysrem
(v) The voltage zero crossings and firing instants are
updated with the latest information of AC convertor- The AC system voltage VG (the commutating voltage)
voltage waveforms. behind inductance X, is assumed perfectly sinusoidal. On
Steps (ii) to (v) are repeated until convergence is achieved. the DC side the smoothing inductance is very large so
In the calculation of the alternating-current injections that the current is perfect DC. The actual values of the
from the convertor (step (ii» the effects of waveform test-system components are indicated (in p.u.) in the
asymmetry and distortion are included, both during and Figure. Details of the per-unit system are given in Appen-
outside the commutation intervals. The basic formulation dix 12.3.
is given in Appendix 12.1. For a specified DC current of say 1 p.u., the power-
The AC and DC system models used in the derivation flow solution provides all the information required to
of the alternating-voltage and direct-current waveforms derive the phase-current waveform (Le. IX = 20°, Jl. = 19°)
(steps (iii) and (iv» are frequency dependent, i.e. the trans- and with it the current harmonics. Under these condi-
mission lines are represented by their equivalent-It cir- tions the conventional (direct) solution can be used to
cuits. calculate the voltage harmonic distortion and the results
are shown in Table 1.
4 Transient convertor simulation (TCS)
5.1 Verification of IHA
State-space and nodal analysis can be used to derive a set Instead of the voltage source, the intermediate node VT ,
of dynamic equations which allows a systematic solution shown in Fig. 1, is used as the commutating voltage in
of the network in the time domain. the iterative solution.
At each step of the solution the instantaneous values The fundamental component of the phase current,
of the network voltages and currents are checked against derived in the previous section, is used to calculate the
specified logical constraints, to decide whether a topo- corresponding voltage component at VT and this value,
logical change is required. together with those of the direct current Id and firing
angle ct. constitute the initial conditions for the iterative parallel resonance between the filters' capacitance and
process described in Section 3. At each iteration the latest the system inductance
value of VT is used as the commutating voltage and the (ii) the relative strength of the AC system, normally
process is repeated until convergence. represented by the SCR (short circuit ratio).
Preliminary results obtained differed markedly from
those of Section 5. The deviations were traced to different The sensitivity of the IHA algorithm to these two factors
interpretations of the specified conditions, and in particu- is investigated in Sections 6.1 and 6.2.
lar the implementation of the zero crossing points used as
a reference for the delay angle. Once the above anomalies 6.1 Symmetrical system
were corrected, the results, illustrated in Table 1, show Fig. 4 represents a six-pulse rectifier connected to a sym-
very good correspondence with the exact solution. metrical system (represented by its Thevenin equivalent)
La= 00
10 (ii)
new results, illustrated by the dotted line in Fig. 5, show lion does not converge. This conclusion is reached by
that the algorithm converges in two or three iterations .observing that the convertor alternating-voltage wave-
irrespective of the SCR. forms settle down to a steady state. For instance, Fig. 6
shows two cycles of the waveforms predicted by TCS
6.2 Asymmetrical system conditions (with an SCR of 5.3) after ten cycles of simulation.
The system of Fig. 7 represents the New Zealand South
Island network in the vicinity of a large power convertor
(at TlW AI) which, for simplicity, is assumed to be a six-
pulse rectifier. The convertor is fed from a long-distance
untransposed transmission line which introduces some
asymmetry into the test system, and therefore produces
triplen harmonic currents.
Fig. 8 illustrates the three-phase impedance loci of the
test system, seen from the convertor busbar; these loci
have been derived by injecting 1 p.u. harmonic currents
at that busbar (with the convertor disconnected).
With this test system the IHA algorithm failed to con-
verge when the short circuit ratio was reduced to 5.3.
Earlier studies carried out with iterative algorithms had cycles
concluded that computational divergence reflects a har- Fig. 6 TCS results Jar a case olltside t"e converging region oj tlte
monic instability in the real system or at least a weakness IHC algoritltm (SCR = 5.3)
which may result in unreliable convertor operation (i) Phase a voltage at the filters
(ii) Phase b voltage at the filters
[5, 7]. So far, however, no definite proof has been offered (iii) Phase c vohage at the filters
to substantiate such claims. Should they be correct, the
use of dynamic simulation should also fail to settle into a 7 Limitation of TeS in a practical system
steady state.
Therefore the system of Fig. 7, with a six-pulse conver- The complexity of a power system cannot be accurately
tor at TIW AI, was analysed in the time domain using the represented by means of a single Thevenin equivalent and
frequency-dependent model described in Section 7.1. therefore the algorithmic verification of TCS carried out
The TCS results show that it is possible to operate the in Section 5 gives no indication about its accuracy in a
system with SCR levels ill the region where the IHA solu- real system. Even with a more comprehensive equivalent
circuit the state-space model makes no provision for the
frequency dependence of the individual parameters
involved.
To assess the error incurred by ignoring such fre-
quency dependence, the test system of Fig. 7 was used
instead.
In this case the results of an accurate IHA model, with
the transmission lines represented as equivalent-n
matrices were used as a reference. These results are
compared in Table 2 with a dynamic simulation using
a Th6venin equivalent of the AC network and they
Table 2: TCS/IHA comparison for the system of Fig. 7 with
Fig. 7 Practical AC system witlt a 6-plllse bridge rectijier at THV AI ThBvenin AC system equivalent
(a) Harmonic current injections into AC system (phase a)
350
TCS IHA
300
1 1.5136 1.5098
250 5 0.2951 0.2822
7 0.1784 0.1661
11 0.0971 0.0791
13 0.0698 0.0522
17 0.0351 0.0223
19 0.0248 0.Q155
23 0.0083 0.0131
25 0.0056 0.0119
clearly show that, in the absence of detailed three-phase Table 3: TCS/IHA comparison for the system of Fig. 7 with
and frequency-dependence representation, the TCS solu- frequency matched AC system equivalant
tion cannot provide accurate information. (a) Harmonic current injection into the AC system (phase A)
TCS IHA
7.1 Modified res to include the effect of frequency
dependence 1 1.4794 1.4782
5 0.2850 0,2849
It is not practical to represent the frequency dependence 7 0.1729 0.1707
of each individual system component in the dynamic 11 0,0916 0,0900
simulation model. A more realistic solution is to derive 13 0,0652 0,0631
the simplest AC system equivalent that can provide a 17 0,0311 0,0294
19 0,0208 0.0196
reasonable match at the main resonant frequencies. 23 0,0068 0.0078
A simple algorithm based on the impedance loci of 25 0,0054 0,0068
Fig. 8 has recently been developed by the authors [6] to
derive a three-phase dynamic equivalent with reasonably (b) Voltage harmonics
accurate harmonic-frequency dependence. The applica- Phase a Phase b Phase c
tion of this algorithm to the system of Fig. 7 results in the
equivalent system of Fig. 9 which provides a greatly res IHA Tes IHA res IHA
improved frequency match, as shown in Fig. 10. 1 1,0467 1,0467 1,0476 1.0476 1,0297 1,0297
5 0.0096 0.0098 0,0097 0.0098 0,0100 0.0101
7 0,0033 0,0059 0,0062 0.0057 0.0059 0.0055
11 0.0021 0.0015 0.0014 0.0015 0,0014 0.0016
13 0.0018 0.0014 0,0016 0,0014 0.0014 0,0014
17 0.0119 0,0079 0,0072 0.0076 0.0093 0.0089
19 0,0169 0,0070 0.0055 0.0060 0.0047 0.0050
23 0.0006 0,0009 0.0010 0,0009 0.0006 0.0010
25 0,0001 0,0001 0,0003 0,0002 0,0001 0,0002