On Voltage and Power Stability in Ac/Dc Systems: Working Group 14.05
On Voltage and Power Stability in Ac/Dc Systems: Working Group 14.05
Working group
14.05
April 2003
ON VOLTAGE AND POWER STABILITY
IN AC/DC SYSTEMS
Working Group
14.05
T. Adhikari (India), J.B. Davies (Canada), M.A Eitzmann (USA), A. Golé (Canada), R. Hirvonen
(Finland), J.Makino (Japan), V. Khoudiakov (Russia), D. Kumar (India), P.C.S. Krishnayya (Canada),
H.A.D. Lee (Malaysia), W. Long (USA), R.S. Moni (India), C. Parker (Australia), C.A.O. Peixoto
(Brazil), J. Reeve (Canada), C.V. Thio (Canada)
Copyright © 2002
Tout détenteur d'une publication CIGRE sur support papier ou électronique n'en possède
qu'un droit d'usage.Sont interdites,sauf accord express du CIGRE, la reproduction totale ou
partielle autre qu'à usage personnel et privé, et toute mise a disposition de tiers, dont la
diffusion sur un réseau intranet ou un réseau d'entreprise.
Copyright © 2002
Ownership of a CIGRE publication, whether in paper form or on electronic support only
infers right of use for personal purposes..Are prohibited, except if explicity agreed by CIGRE,
total or partial reproduction of the publication for use other than personal and transfer to a
third party; hence circulation on any intranet or other company network is forbidden.
To the memory of Clarence V. Thio
Contents
Preface vii
Summary ix
v
vi Contents
References 77
1 This section is identical with section 1.2, but since it is essential for the subsequent
This report is the result of work done in Cigre WG14.05 over the years 1996 -
1999. By various reasons all the various contributions have not been put together
in a single document until now. When this work started, it was not clear to the
convenor, nor to the working group members, the amount of work required to
cover the topic in a systematic way. To do this substantial additional work would
be needed. Still, the report contains significant new results, e.g. concerning the
definitions of voltage and power stability an AC/DC systems, section 1.5. This
definition is adopted to comply with the definitions introduced in power systems
without any HVDC systems. But also other chapters contain new material.
As agreed at the SC14 meeting in Paris 2000, the material available has been
put together without any excessive editing. Therefore, it might occur repetitions
and the coherency might not be complete. Despite this, the report contains a
lot of material that is deemed to be useful for the power industry.
Since the work started in 1996 new types of converters have been introduced
for high power application, i.e voltage source converters, which not are discussed
or analyzed in this report. Only line commutated current source converters and
capacitor commutated converters are covered. However, the analysis methods
and tools in the report could also be applied to systems using voltage source
converters or other types of converters.
The main contributors of the chapters in this report are:
Chapter 1: G. Andersson
Chapter 2: H.A.D. Lee
Chapter 3: A. Golé
Chapter 4: J. Reeve
Chapter 5: H.A.D. Lee
Chapter 6: C.V. Thio and J.B. Davies
During the course of this work the working group suffered a significant loss
when Clarence Thio passed away. Clarence was one of the most active members
of the working group and his contributions were essential for the progress of the
work documented in this report. His insights and experience together with his
sincere interest to understand the problems dealt with inspired all members of
the working group and had a decisive influence of the work pursued here. We
will all remember Clarence as a prominent engineer and a great friend.
vii
viii Preface
Göran Andersson
Convenor WG14.05
Summary
Voltage and power instability has been a limiting factor for many HVDC projects
during the last decades. Pioneering and systematic work to understand and
analyse the interactions leading to this instability started in the early eight-
ies, and joint Cigré and IEEE working group presented the first comprehensive
overview of the topic in the early nineties, ref. [4] of this report. After that
report was published substantial work has been done, both in the academia and
the industry, to further develop analytical and simulation tools in order to im-
prove the understanding and the computational tools. These later developments
are reported in this report. A brief overview of the report is given below.
In the late eighties and early nineties extensive work was done within Cigré
and IEEE to better understand and analyse voltage stability in (pure) AC sys-
tems. Among other things definitions regarding voltage stability were formu-
lated. These definitions are reviewed in chapter 1 and definitions applicable for
AC/DC systems that are consistent with those of AC systems are formulated.
The concepts of Maximum Power Curves (MPCs) and Voltage Sensitivity
Factors (VSFs) were early introduced for the analysis of simple AC/DC systems.
In chapter 2 the theoretical foundations of these tools are further developed and
extensions to more complex systems are carried through. Systems with loads,
parallel AC lines, and multi-infeed HVDC systems are included. In chapter 5
further considerations for multi-infeed systems are discussed.
The methods dealt with in chapter 2 are based on quasi steady-state consid-
erations. The virtue is that very powerful computational tools can be developed,
which enable very fast analysis of different system conditions as well as analysis
of the influence of different design and system parameters. However, the dy-
namics of the system is not explicitly modelled and the analysis based on quasi
steady-state methods must be extended with time simulations for a complete
analysis of system. In chapter 3 comparisons between the static analytical meth-
ods and full detailed time simulations are given. The conclusion is that static
methods, i.e. the ones based on MPC and VSF, can be used for determining the
limits of voltage/power stability, but for determining the dynamic properties a
full dynamic model of the system is required.
In chapter 4 the influence of loads, parallel AC lines, etc, on the stability
limits is studied. It is also demonstrated that the stability limits of the Capacitor
Commutated Converters (CCCs) are significantly higher as compared with those
of the standard line commutated converters. (VSC converters were not included
in this analysis, see below.)
In chapter 6 a discussion how the developed tools can be applied to real
planning and design problems is given. Particularly, the interplay between the
use of the static and dynamic analysis tools are exemplified. The Nelson River
HVDC scheme is the system mostly considered in this chapter.
The analysis methods described in this report can be used for all converter
types. However, in this report the methods have only been applied to the con-
ventional line commutated converters, and Capacitor Commutated Converters
(CCCs), but not to systems with Voltage Source Converters (VSCs). The ex-
tension to this latter kind of HVDC systems should be straightforward.
ix
1
Basic Theory and Definitions
1.1 Introduction
Modern power systems are becoming increasingly complex and large. To be able
to predict the behaviour of these systems under different operating conditions
engineers are more and more dependent on sophisticated software. Since the
models of the power system and the software used in analyses are getting more
and more complex the requirements on the engineers involved in these studies
get higher. Even if the analytical tools are much more powerful today, it is still
of utmost importance that the engineers have a thorough understanding of the
basic physical processes and interactions in the system. Without this knowledge
it is not possible to understand phenomena in complex systems, and to design
appropriate control and protection systems that can alleviate possible adverse
interactions.
Another important prerequisite for the establishment of a good understand-
ing of a phenomenon is that commonly accepted and appropriate definitions
exist. This will facilitate the discussion between experts and contribute to the
development of the knowledge and understanding. It is not uncommon that
disagreements between experts on an issue are due to lack of accepted defini-
tions. The importance of clear definitions, that are commonly accepted, should
be obvious.
In this chapter a basic explanation of the power and voltage instability phe-
nomena is given. This is done by introducing a simple model of an AC/DC
system, which exhibits the basic interactions that are dominating also in real
size power systems. Definitions of power and voltage stability are also given.
These definitions are based on definitions of voltage stability for AC systems.
The chapter is organised as follows. First a brief physical explanation of the
power and voltage instability in AC/DC systems is given. Then the definitions
of voltage stability in AC systems are reviewed and discussed. A simple AC/DC
system model is then introduced and definitions of voltage and power stability
in AC/DC systems are proposed.
1
2 1. Basic Theory and Definitions
As seen these definitions are quite general and refer to how power system volt-
ages react to disturbances in the system.
Ze U l ∠φ
Infinite bus
U N ∠0 Load
P l + jQ l
This set of definitions are obviously more restrictive than those of the previous
subsection. This will be illustrated in subsection 1.3.4.
The rationale for this definition is obvious since voltage control is mostly ac-
complished by controlling reactive power from controllable sources, such as syn-
chronous machines, reactive shunt elements, SVCs, etc. The authors of ref. [3]
suggest that in addition to the requirements given in the above definition, the
system needs also to be voltage regular in order to be voltage stable.
A
1
0.8
B
Load Voltage
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Active Load Power
Figure 1.2. Nose curve of simple AC power system. Power factor of load
=1. Impedance of line = 0.4 p.u., impedance angle = 90◦ . Voltage of
infinite bus = 1.0 p.u.
If the voltage of the infinite bus is kept constant, both magnitude and phase,
and the active load demand Pl is varied, the so-called nose curve is obtained.
This curve shows how the load voltage Ul varies as Pl is changed. It is assumed
that for a given active load there corresponds a unique reactive load, for instance
by assuming a constant power factor. A typical example of a nose curve is shown
in Figure 1.2.
The curve in Figure 1.2 shows that for each load power level below the power
corresponding to the point B, there exist two solutions: One with “high” voltage
and one with “low” voltage. At the point B these two solutions coalesce, and
above this power level no solutions exist with the given data of the system.
The point corresponding to B is often referred to as the Point of Maximum
Loadability (PML).
Physically the curve in Figure 1.2 can be constructed in the following way.
At the load bus a shunt admittance (impedance) is connected, and the value
of this admittance is increased from zero, corresponding to point A, to infinity,
corresponding to point C. The ratio between the active and reactive parts of the
admittance should be kept constant to a value determined by the power factor
of the load. For any given value of the load admittance, the load voltage can
be calculated and the load power determined, which means that to any given
value of the admittance there exists a unique point on the nose curve.
kp Ustab
≥1 0
0.75 0.45
0.5 0.58
0 0.71 (“B”)
Table 1.1. Regions of voltage regularity for different active load voltage
dependencies, kp Operating points on the nose curve of Figure 1.2 with
voltages greater than Ustab are voltage regular, i.e. V SF > 0, those with
voltages less than Ustab are not voltage regular. The transition point for
constant power loads, i.e. loads with kp = 0 corresponds to the point B
in Figure 1.2.
∆U
V SF = >0 (1.1)
∆Q
As seen this means that the system is voltage regular according to the definition
of ref. [3], and a necessary condition for stability is thus voltage regularity. It is
straightforward to calculate the VSF for the simple system above. To simplify
the calculations it is assumed that the load is purely active and it has voltage
dependence according to
U kp
P = P0 (1.2)
U0
With these assumptions it can be shown that regions of the nose curve where
VSF is positive depends on the exponent kp . It turns out that the nose curve
is divided into two regions, one where VSF is positive and one where VSF is
negative. At the transition point VSF goes to infinity. If it is assumed that the
line impedance Ze is lossless and that the voltage of the infinite bus is kept at
1 p.u., the stability regions for some exponents kp are given in Table 1.1 can be
summarized as:
For constant power loads, i.e. kp = 0, only operating points on the upper
part of the nose curve are voltage regular and the transition takes place at the
point B of the nose curve. If the voltage dependence is as constant current load
or with a higher exponent, i.e. kp ≥ 1, all points on the nose curve are voltage
regular. For values of kp in between, the border between voltage regularity and
irregularity lies on the lower part of the nose curve as indicated in Table 1.1.
Negative values of kp are deemed unphysical and not considered.
It should be noted that in principle a voltage control strategy based on that
reactive power injection decreases the voltage could be used in the regions found
voltage irregular above. However, such a voltage control scheme had to detect
6 1. Basic Theory and Definitions
-Pd + jQd
Id
Us 0
Z θ Ut δ Ud
p
τ:1
Qc
the sign of VSF before acting, and such a controller is deemed to be highly
hypothetical and not to be considered in a practical power system.
Furthermore, one could here imagine that the load characteristics are such
that the initial voltage dependence ensures stability, but slower control actions
are such that the operating point moves away from the initial one, leading
to either too high voltages in the system or to a voltage collapse. Therefore,
operation on the lower part of the nose curve might be tolerable for short period
of times, but to guarantee stability, normal (steady state) operation must be on
the upper part of the nose curve.
Conclusions
Two main conclusions can be drawn from the simple example analysed above:
• According to the IEEE definition, operation on the lower part of the nose
curve is unstable.
These two conclusions will be very important when we discuss stability for
AC/DC systems in the next section.
• Short Circuit Ratio, SCR, is the per unit admittance of the impedance
Z∠θ in Figure 1.3. The per unit value is based on the nominal MW rating
of the HVDC converters and the nominal AC voltage of the converter AC
bus.
From the definitions of SCR and ESCR it is clear the these quantities are com-
plex numbers. However, in most cases the angle θ is close to 90◦ . Furthermore,
as shown in ref. [4], system performance is rather insensitive to the angle θ as
long as its value is close to 90◦ . Therefore, if not explicitly stated, it is assumed
that θ = 90◦ , and only the magnitude of SCR is given. In this case ESCR can
be written as ESCR = SCR − Qc .
The discussion in the following will be based on the model system of Fig-
ure 1.3, and the assumptions briefly reviewed above are adopted.
being that when the DC current is varied that the controls uniquely determine the states of
the system so that the load flow equations could be solved. Other control modes that are
conceivable are constant firing angle control, constant DC voltage control, etc. see section 1.6
8 1. Basic Theory and Definitions
0.8
D.C.−Power
0.6
0.4
SCR = 3.5
0.2
0
0 0.5 1 1.5 2 2.5
DIRECT CURRENT
Figure 1.4. Maximum Power Curve (MPC) of the system in Figure 1.3.
SCR = 3.5, γ = 17◦ . Initial operating point: Id = 1 p.u. and Ut = 1 p.u.
1.8
1.6
1.4
DC POWER & AC−VOLTAGE
1.2
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5
DIRECT CURRENT
Figure 1.5. Maximum Power Curve (MPC) of the system in Figure 1.3.
SCR = 1.5, γ = 17◦ . Initial operating point: Id = 1 p.u. and Ut = 1 p.u.
1.2
1.1
0.9
MAP
0.8
0.7
0.6
0.5
0.4
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Power
Figure 1.6. Nose curve for the AC/DC system in Figure 1.3, SCR = 3.5.
As indicated above the nose curves of the pure AC system and the MPCs
exhibit great similarities. To verify this, curves showing the variation of AC
voltage as function of the DC power are plotted in Figure 1.6 and Figure 1.7
for the MPCs in Figure 1.4 and Figure 1.5, respectively.
1.8
1.6
1.4
AC−Voltage
1.2 MAP
0.8
0.6
0.4
0.2
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Power
Figure 1.7. Nose curve for the AC/DC system in Figure 1.3, SCR = 1.5.
Table 1.2. Sign of Voltage Sensitivity Factors of the systems in Figure 1.4
and Figure 1.5 at the nominal operating point. C.C. = Constant DC
Current Control. C.P. = Constant DC Power Control.
but when calculating the V SF other quantities than the DC current could be
varied. When calculating V SF , the control mode of the HVDC is crucial to
the result. For the systems of Figures 1.4 and 1.5 the control modes constant
power control and constant DC current control will be considered. Other con-
trol modes can also be analysed, but for the discussion here it is sufficient to
consider these two modes. A simple analysis gives the results of Table 1.2.
As seen from Table 1.2 all cases are voltage regular except the last one, which
corresponds to constant DC power control in the weak AC system case. It is
thus seen that for the weak system (SCR = 1.5), constant power control implies
a negative V SF , while constant DC current control implies positive V SF .
systems. From the discussion above the following two definitions are proposed
for an AC/DC system:
1.6 Discussion
It is clear that operation to the left of MAP is stable according to the definitions
above.
From the proposed definitions it is obvious that operation to the right of the
MAP will not be power stable. When it comes to voltage stability it is clear
that in constant DC power control operation to the right of MAP is unstable4 .
A control action to increase the voltage magnitude, e.g. by connecting a shunt
capacitor, will result in a decrease in voltage, see Table 1.2. (It is here assumed
that the constant power control loop can control the DC power to the desired
value.) This will in a real system lead to an unstable situation.
In constant DC current control and to the right of MAP, the system is
initially voltage stable since the V SF is positive. However, the system will
not transmit the desired amount of power and any attempt to adjust the DC
current with a standard DC power controller will lead to unstable conditions.
Therefore, the system can operate transiently to the right of MAP, but in order
to use the HVDC link for power control it must be brought back to the left side
of the MAP.
The definition of voltage stability involves arbitrary disturbances, such as
earth faults and trippings of lines. For such cases analyses based on MPCs are
not sufficient, since the large signal dynamics of the system plays an important
role. The restart strategies of the HVDC after faults are very important when
analysing (large disturbance) stability of the interconnected AC/DC system.
In the examples above it has been assumed that the inverter is in constant
γ control. However, as indicated above MPCs can be computed for other types
of HVDC controls such as constant DC voltage control or other types of γ
4 It is here assumed that the DC power controller has a positive gain in its main loop, i.e.
1.2
0.8
DC POWER
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
DIRECT CURRENT
Figure 1.8. MPC for an AC/DC system with constant DC voltage control.
control. The MPC for such a control mode will look different, but there will be
a M AP point for such modes also. The M AP will occur for the point where γ
has reached its minimum permissible value, γmin . The MPC to the right of this
point would coincide with an MPC with constant γ = γmin . The same reasoning
for these types of control as for the constant γ control applies, which means that
only steady state operation to the left of M AP is stable in a practical system.
An example of an MPC with constant DC voltage control is given in Figure 1.8.
To distinguish between MPCs plotted for different control modes, an index
indicating which DC controls that are used when calculating can be used. The
standard MPC with constant γ control should thus be denoted MPCγ and the
MPC in Figure 1.8 should be MPCUd .
2
Analytical Tools and Methods
2.1 Introduction
E ∠ψ U ∠δ P d, – Q d
τ:1 Ud
z ∠θ
Xc γ, µ
Id
bc
13
14 2. Analytical Tools and Methods
dPd
=0 (2.1)
dId
where
1
ESCR = − bc
z
2Ud
r= b = 2 Qd − Pd tan(γ + µ)
cos γ + cos(γ + µ)
4 2 2 ∂Qd 2 ∂Pd
det JI = U ESCR − Qd + U U ESCR + Qd − Pd Pd − U
∂U Id ∂U Id
and ∂Pd /∂UI |Id and ∂Qd /∂U |Id are the voltage dependence of the converter
active and reactive power, respectively, for the constant current control mode
(see eq. (2.9)). From eq. (2.2), the M AP condition thus occurs when
U 4 ESCR2 + b U 2 ESCR + Qd − Q2d − Pd2 = 0 (2.3)
If the system operating point corresponds with the M AP condition under nom-
inal conditions, then the prevailing ESCR is said to be critical, CESCR.
Eq. (2.3) then solves to
1 b b 2 2
CESCR = 2 − + − Qd + Pd (2.4)
U 2 2
2.2. Power Stability 15
ESCR=
2.1
1.8
Pd [p.u] 1.48
1.2
IMAP
Id [p.u]
1.2
1.48
1.8
2.1
U [p.u] ESCR=
Id [p.u]
ESCR=1.4775
Pd
dPd
dI d
(equation 2.2)
Id [pu]
induces a small and same directional change in the converter AC bus voltage.
This characteristic is consistent with stable system behaviour. An increasing
V SF corresponds to a decreasing voltage stability margin and when the V SF
becomes infinite, the transition to an unstable system occurs. When the V SF
is negative, a small change of the reactive power induces an opposite change
in the converter AC bus voltage. This characteristic corresponds with unstable
system behaviour.
The V SF can be derived mathematically from the same power flow equa-
tions used to formulate the gradient of the MPC in eq. (2.2). By the same
linearization of these equations with respect to U and δ, and eliminating δ, the
V SF is given by
∆U/U
V SF = =
∆Q
U 2 ESCR + Qd (2.7)
∂Qd 2 ∂Pd
U 4 ESCR2 − Q2d + U U ESCR + Qd + Pd Pd − U
∂U ∂U
where ∂Pd /∂U and ∂Qd /∂U are the voltage dependence of the converter active
and reactive power, respectively, and are given by the control mode as follows
CEA/Constant Power Control
∂Pd ∂Qd 2
=0; = Qd − Pd tan[γ + µ) (2.8)
∂U ∂U U
Example 2 Using the same parameter values as in Example 1 for the single-
infeed configuration of Figure 2.1, the V SF as a function of ESCR for different
HVDC control modes are shown in Figure 2.5. It is seen that HVDC systems
with constant DC voltage control have small V SF even for low ESCR, implying
that this control mode is inherently stable. In contrast, constant power control
is inherently unstable at low ESCR as shown by the negative V SF below the
CESCR.
a
b
VSF
c
a. const. power
a b. const. current
c. const. DC voltage
1.4775
ESCR
γ
Kp
H
+
- +
H + α
KI
s L
γref
L
The power flow equations comprise AC power balance equations at the converter
AC bus, and converter steady-state DC power, current, and voltage equations.
This set of six equations can be found in [8].
The operating state vector comprises system states which will vary dur-
ing dynamic operation of the system. Its composition depends on the control
mode under consideration but includes the controlled and controlling system
states. For the system model of Figure 2.1, the operating state vector vxy is
(α, φ, δ, Ut , Ud , Id , Pd ) for constant power (constant γ), (α, γ, φ, δ, Ut , Ud , Pd ) for
CEA (constant current) and constant DC voltage (constant current), control.
α and φ are the inverter firing angle and power factor at inverter AC bus,
respectively, and other symbols are self-explanatory from Figure 2.1.
Thus the Jacobian, say Jxy , computed at the operating point (power flow
solution) is given by
Jxy · ∆vxy = 0 (2.11)
Rewriting eq. (2.11) gives
dP d
d Id
ESCR
P , –Q
E 1 ∠ψ 1 d1 d1
U 1 ∠δ 1 U d1
z 1 ∠θ 1 d x1 γ 1, µ 1
1 Id1
P12, Q12 b c1
z ∠θ
12 12 P d2, – Q d2
z ∠θ U d2
2 2
2 d x2 γ 2, µ2 Id2
E 2 ∠ψ 2 U 2 ∠δ 2 b c2
U 1 ∠δ 1 1 P12, Q12 2 U 2 ∠δ 2
z 12 ∠θ 12
P d1, – Q d1 P d2, Q d2 E 2 ∠ψ 2
E 1 ∠ψ 1 Id
U U d2
z 1 ∠θ 1 γ 1, µ 1, α 1d1 µ 2, α 2 z 2 ∠θ 2
b c2 b c1
Similar to the single-infeed case, the power stability methods for these higher-
order system configurations are derived from their linearized power flow equa-
tions. These equations also assume constant Thevenin voltages, thus the source
2.5. Extension to Higher-Order Models 21
where;
∆Pd is a vector of incremental DC power at the converter DC buses.
∆Id , ∆δ, and ∆U/U are vectors of incremental change in the system states
comprising the DC current, converter AC bus voltage angle and relative
magnitude, respectively.
The Jacobian submatrices are the partial derivatives of the power flow equa-
tions with respect to the system states. The subscript labels d, ac denote the
quantities related to the DC power, AC active and reactive power, respectively.
Note that Id is additionally chosen as a system state so that the Jacobian is an
extended form of the well known Newton-Raphson model. The expressions for
these Jacobian submatrices are given in [17].
For no active or reactive power injections at the converter AC buses, ∆Pac
,∆Qac may be assumed to be zero. Thus eq. (2.14) may be reduced to
where
−1
JMP C = JdI − JdU JR1 JR2 , JR1 = JQU − JQδ JP−1
δ JP U
and
JR2 = JQI − JQδ JP−1
δ JP I
where JMP Cii is the i-th diagonal element of JMP C . When JMP Cii vanishes,
eq. (2.16) is akin to eq. (2.1) for the single-infeed case. Thus JMP Cii is used to
determine the M AP condition for these higher- order configurations. Similar
to the single-infeed MPC, the SMPC for the higher-order configurations may
also be derived. This is done practically by numerically computing the DC
22 2. Analytical Tools and Methods
Side
P1 Front
Pd1 [pu.]
Plan
P2
where
2.5. Extension to Higher-Order Models 23
P1 1.6 pu.
ESCR1=1, P2
ESCR2=3
z12=0.6 pu.
Pd1 [pu.]
IdN2=0.2 pu.
Id1 [pu.]
JMPCii=0 P1
Id1 [pu.]
unstable dP d1
= 0
d I d1
stable
P2
IdN2 [pu.]
where ∆pd = η∆Pd and ∆id = η∆Id are the vector of incremental modal DC
power and DC current, respectively. For the i-th eigenmode, eq. (2.18) can be
written as
∆pdi
= λi (2.19)
∆idi
where
nt
∆pdi = ηij ∆Pdj , i = 1, . . . , nt (2.20)
j=1
24 2. Analytical Tools and Methods
nt
∆idi = ηij ∆Idj , i = 1, . . . , nt (2.21)
j=1
∆pdi
=0 (2.22)
∆idi
This is akin to eq. (2.1) for the single-infeed case. However, this is a modal
relationship, thus corresponding to the MPC for the single-infeed configuration,
the Modal Maximum Power Curve (MMPC) may be defined using eq. (2.19).
The system operating conditions resulting in eq. (2.22) would thus represent the
Modal Maximum Available Power of the system. Alternatively, eq. (2.19) may
be expressed as
∆idi 1
= (2.23)
∆pdi λi
Again, positive eigenvalues cause positive modal DC current changes in response
to positive modal DC power order variations. This corresponds to stable DC
power control behaviour. For a zero minimum eigenvalue, an infinitesimally
small change in the modal DC power order causes an infinite change in the
modal DC current, hence instability of the DC power control. Therefore, for
stable DC power behaviour eq. (2.19) and eq. (2.23) imply
Eq. (2.19) may be used to derive graphical plots in parameter space to define
the MMPC stability boundary. In such case, the plots are values of parameters
for which the system minimum eigenvalue, λmin , becomes zero. Thus, the plots
are the MMPC stability boundaries that demarcate the parameter space into
stable and unstable regions.
Example 5 Consider the multi-infeed configuration of Figure 2.8 with the fol-
lowing system parameter values: ESCR1 = 1, ESCR2 = 3, z12 = 0.6 p.u. The
MMPC stability boundary in the PdN1 − PdN2 parameter space (the subscript
label N denotes rated quantities) is shown by the plot in Figure 2.13. In Fig-
ure 2.13, the region above and below the curve are thus unstable and stable
regions, respectively. The curve comprises values of the parameters for which
λmin becomes zero, thus it is the system MMPC stability boundary.
unstable
PdN1 [pu.] λmin=0
stable
PdN2 [pu.]
where
∆P is a vector of incremental change in converter AC bus active powers.
∆Q is a vector of incremental change in converter AC bus reactive powers.
∆δ is a vector of incremental converter AC bus voltage angles.
∆U/U is a vector of incremental converter AC bus relative voltage magnitudes.
JP δ , JP U , JQδ , JQU are the partial derivatives of the power flow equations with
respect to the voltage angles and magnitudes.
If the active power is held constant at the converter AC buses and only
incremental reactive power is applied, then ∆P may be assumed to be zero, and
eq. (2.25) thus reduces to
∆U
∆Q = JR (2.26)
U
where JR = JQU − JQδ JP−1 δ JP U is the QU reduced Jacobian matrix. Alterna-
tively, eq. (2.26) can be expressed as
∆U −1
= JR ∆Q (2.27)
U
Though eq. (2.27) appears similar to the V SF in eq. (2.6), there is no direct
relationship for the converter AC bus in question, due to the coupling from
−1
the other non-zero off-diagonal elements of JR . To obtain a decoupled form,
an eigenvalue decomposition technique similar to the MMPC approach in sub-
−1
section 2.5.1 may be used. Thus applying the same decomposition to JR in
eq. (2.27) gives
∆qi
∆ui = (2.28)
λi
26 2. Analytical Tools and Methods
where
nt
∆Uj
nt
∆ui = ηij ∆qi = ηij ∆Qj i = 1, . . . , nt
j=1
Uj j=1
Example 6 Consider the multi-infeed configuration of Figure 2.8 with the fol-
lowing system parameters; PdN 1 = 1 p.u., PdN 2 = 0.5 p.u., z12 = 0.6 p.u.
Figure 2.14 shows the M V SF stability boundary in the ESCR1 − ESCR2
parameter space. For all ESCR values on the boundary, λmin is zero.
Participation Factors
For the modal approach, not only is the minimum eigenvalue useful in estimating
the proximity to voltage instability but the eigenvectors also contain information
concerning how critical system locations are. The eigenvectors may be used to
define a participation factor [11] given by
where
pij is the participation factor for the converter AC bus i in the j-th voltage
variation mode
∆ηji is the i-th element of the j-th row left eigenvector matrix.
2.6. Incorporation of Local Devices 27
stable
ESCR1
λmin=0
unstable
ESCR2
∆Ψij is the i-th element of the j-th column right eigenvector matrix.
nt is the number of converter AC buses.
Physically, Ψij is a measure of the activity of converter AC bus i in the j-th
voltage variation mode, ηji is the weighting of the contribution of this activity,
their product pij is a measure of the net participation of bus i in the j-th voltage
variation mode.
In the context of multi-infeed configurations, bus participation factors com-
puted from the eigenvectors associated with the minimum eigenvalue provide
information on the critical converter AC bus. The bus with the largest par-
ticipation factor is the critical bus. Consequently it is also the most effective
system location for implementation of remedial measures [13].
2.6.1 Loads
In this subsection, the basic equations associated with power stability are pre-
sented for the various system configurations with static loads connected at the
converter AC buses as shown in Figures 2.15 - 2.17. It is assumed that the load
models have static load characteristics given by
P
= KP C + KP I U + KP Z U 2 + KP U U nP U (2.31)
PL
Q
= KQC + KQI U + KQZ U 2 + KQU U nQU (2.32)
QL
where
28 2. Analytical Tools and Methods
E ∠ψ U ∠δ P , –Q
d d
τ:1 U
d
z ∠θ
d
x
γ, µ
PL I
d
QL b
c
U 1 ∠δ 1 1 2 U 2 ∠δ 2
z ∠θ
12 12
P d1, – Q d1 P d2, Q d2 E 2 ∠ψ 2
E 1 ∠ψ 1 Id
U U d2
z 1 ∠θ 1 γ 1, µ 1, α 1d1 µ 2, α 2 z 2 ∠θ 2
b c2 b c1
PL1 , QL1 PL2 , QL2
KP I = 1 − (KP C + KP Z + KP U )
KP C , KQC is the per unit of constant MVA active and reactive load, respec-
tively.
KP I , KQI is the per unit of constant current active and reactive load, respec-
tively.
KP Z , KQZ is the per unit of constant impedance active and reactive load,
respectively.
KP U , KQU is the per unit of voltage dependent active and reactive load,
respectively.
nP U , nQU is the voltage sensitivity exponent for active and reactive load,
respectively.
Using the same approach as in subsections 2.2.1 and 2.5.1, the power flow
equations incorporating the load models of eqs . (2.31) and (2.32) are similarly
linearized and reduced to give the following.
For the single-infeed configuration, Figure 2.15, the following holds
dPd rU cos(γ + µ) U 4 ESCR2 + (b + KQ) U 2 ESCR + Q 2 − Pd P d
d − Q
d
=
dId det JI
(2.33)
2.6. Incorporation of Local Devices 29
P , –Q
E ∠ψ U ∠δ d1 d1
1 1 1 1
z ∠θ U
1 1 d1
PL1 d γ ,µ
1 x1 1 1 I
QL1 b d1
c1
z ∠θ
12 12
E ∠ψ P , –Q
2 2 d2 d2
z ∠θ U
2 2 d2
U 2 ∠δ 2 2
d γ ,µ
x2 2 2 I
d2
PL2 b c2
QL2
where
Pd = Pd − PL , P d = Pd + KP and Q
d = Qd + QL
∂PL
U = [KP I U + 2KP Z U 2 + nP U KP U U nP U ]PL
∂U
The basic equations for the multi-infeed and single-infeed with a parallel AC
line configurations are identical in form, given by
1 b + KQ1 b + KQ1 2
CESCR1 = 2 − +
− Qd1 + (Pd1 P d1 + P0 + Q0 )
U 2 2
(2.35)
where the quantities involved are given in chapter 8 of ref. [17].
as
dPd 4 2 ∂Qd 2
= rU cos(γ + µ) U ESCR + b + U U ESCR + Qd
dId ∂U
(2.36)
2 ∂Pd
− Qd − Pd Pd − U / det JI
∂U
where r, b, and det JI are as previously defined. Here, ∂Pd /∂U , ∂Qd /∂U are
the voltage dependence of the active and reactive power output, respectively, of
the locally connected device. For example, an SVC whose power out given by
PSV C = 0 QSV C = KSV C U (2.37)
within its linear control range would have voltage dependence given by
∂PSV C ∂QSV C
=0 = KSV C (2.38)
∂U ∂U
From eq. (2.36), the M AP condition for the single-infeed configuration of Fig-
ure 2.1 with a general device locally connected may be similarly specified and
the CESCR derived.
CESCR=1.4775
bc=0.54 pu, γ=17o
dx=0.08 pu.
VSF
dP d
d Id
VSF
ESCR
stable
λVSF=0 unstable
λMPC=0
z12 [pu.]
Example 7 These relationships are also illustrated in Figures 2.21 - 2.23 for
the various system configurations of Figures 2.15 - 2.17 having the following
sample parameters: PL1 = PL2 = 0.1 p.u. (constant impedance), QL1 = QL2 =
0.075 p.u. (constant MVA), z12 = 0.6 p.u., γ1 = γ2 = 17◦ , bc1 = bc2 = 0.54 p.u.,
dx1 = dx2 = 0.08 p.u. Initial conditions are P12 = 0.2 p.u. (for multi-infeed)
and - 0.2 p.u. (for single-infeed with parallel AC line), and Ud1 , Ud2 , Id1 , Id2 ,
U1 , U2 are all 1 p.u. For the multi-infeed and single-infeed with parallel AC line
configurations, the constituent converter of interest is the referenced system and
32 2. Analytical Tools and Methods
Single-Infeed MPC V SF
Single-Infeed
with parallel MPC M SV F
AC line (constant power)
M V SF
Multi-Infeed SMPC (constant power,
constant current)
M V SF
Multi-Infeed MMPC (constant power,
constant power)
the neighbouring constituent converter is then the influencing system. For the
single-infeed configuration, the parameter subscript indices are omitted and only
the relevant parameters apply. Note that the relationship for the multi-infeed
configuration without incorporating load models is not shown in Figure 2.22,
but the correspondence is similar to Figure 2.19.
2.7. Relationship Between Power and Voltage Stability Methods 33
dPd
1 d Id
-----------
VSF
ESCR
Figure 2.22. Stability boundaries for multi-infeed configuration with load models.
34 2. Analytical Tools and Methods
MVSF(circle) ,MPC(solid)
Several indices have been developed for the stability analysis of HVDC Systems.
Ainsworth [5] introduced the concept of Maximum Available Power (M AP ).
This index defines the stability limit for an HVDC System with the rectifier
controlling the DC power and the inverter operating at constant extinction
angle. Hammad [7] has introduced the Voltage Stability Factor (V SF ) in which
the change in sign of the ratio (∆Q/∆V ); or incremental reactive power to
incremental voltage determines the onset of instability. Similar indices have
been introduced by Franken and Andersson [10], and by Nayak et al [8]. These
indices are typically based on the steady-state DC equations and thus do not
represent the full dynamic behavior of the system.
More recently, several investigators have attempted to determine the accu-
racy of these indices. Franken et al [10] have used stability programs, as Pilotto
[18], Nayak [8] have attempted to use electromagnetic transient (EMTP-type)
programs. These analyses generally confirm the usefulness of the steady state
indices. However, recent studies by Reeve et al [9] that consider AC and DC
lines in parallel and a more detailed load representation, show variations from
results obtained from the straightforward calculation of these indices.
Similar discrepancies may arise when synchronous compensators are present
at the converter bus. In this case the effective short circuit ratio (ESCR) is
calculated assuming the transient impedance (Xd ) of the synchronous compen-
sator to be lumped with the impedances of the external AC network and filter.
This ESCR value is later used in the calculation of the MPC. The details of
the excitation control system are also usually ignored. In this section the MPCs
obtained with the simplified model are compared with those obtained using a
detailed electromagnetic transient simulation model.
The study system is derived from the first CIGRE HVDC benchmark [19]
(60 Hz). The following changes are made to the benchmark: Inclusion of a local
load at the inverter bus of 300 MVA, 0.95 power factor local load at the inverter
bus. Inclusion of a synchronous compensator rated -165/+300 Mvar on the bus.
The resultant system is shown in Figure 3.1.
37
38 3. Comparison of MPC with Detailed Simulations
Filters/Ac
Local Load
Capacitors
2.5 Ω, 0.5 H 2.5 Ω, 0.5 H (300 MVA,
0.95 pf)
21.67 µF
Ac Filters (300 Mvar)
Syn. Con.
System Equiv.
ESCR 2.5 , 75 System Equiv.
deg. ESCR 2.5 , 80
deg.
Z2 (local load)
Zf (filter)
Vd= 500
kV
Zs
ESCR=1.554
Theoretical MAP with Xd’ Theoretical Vt
Simulation MAP with Xd’ Simulation Vt
1.25
Pdc (pu), Vt (pu)
1.125
1
0.875
0.75
0.6 0.7 0.8 0.9 1 1.1 1.2
Idc (pu)
current (see solid curve in Figure 3.3). A system in power control is stable when
its operating point lies to the left of the peak of the resultant inverted parabola
shaped curve. This is because the power controller effects an increase in power
by requesting an increase in DC current, consistent with the rising curve to the
left of the peak. Obviously, to the right of the peak, a current increase causes a
power decrease resulting in the instability of the power control mode.
Figure 3.3. shows the theoretically calculated curves for DC power and AC
bus voltage obtained from steady-state equations (see Appendix, section 3.4)
3.2. Inclusion of a Synchronous Compensator at the Converter Bus 39
This investigation was carried out in order to check whether the MPC approach
can be used to determine the stability limits of inverters with synchronous com-
pensators on the AC bus bar. In this particular study, the exciter response was
not modelled and it was assumed that the field voltage of the machine remained
at a fixed value. As the flatness of the MPC makes it difficult to determine the
exact point at which the peak occurs, another equivalent approach is to plot
the index CSI = dP/dId as a function of the parameter being varied. The
point at which this quantity crosses the horizontal axis is thus the transition
point from stable to unstable operation. This index is also referred to as the
controller sensitivity index for the constant power mode of operation [8], see
also section 2.4.
It can be argued that lumping the synchronous compensator with the net-
work equivalent as in Model 1 (Figure 3.1) is an oversimplification. This is
because the actual synchronous compensator does not absorb or deliver any
real power in the steady state. No such assumption is made in the traditional
M AP which means that real power may be supplied or delivered by the syn-
chronous machine. The equations in the Appendix, section 3.4, can be modified
to include an additional constraint for the synchronous compensator so that it
produces zero power. A new set of M AP (or equivalently CSI) curves now
results as shown in Figure 3.4. The method for obtaining the simulated curves
is as described in section 3.1.
In simplifying the machine to an equivalent, two assumptions are possible:
Assume a synchronous reactance of Xd or Xd for the machine. As can be seen,
the added assumption makes a significant difference to the CSI and hence to
the instability point on the M AP characteristic. In comparison with a detailed
model the curves generated with the Psyn = 0 summation match the observed
response more closely.
40 3. Comparison of MPC with Detailed Simulations
100
0
-100
1.4 1.6 1.8 2 2.2 2.4
ESCR(Xd")
Z2 (local load)
Zf (filter)
Vd= 500
kV
M Deatiled Model of
Sync machine
with fixed Efield
Zs
ESCR=1.597
Theoretical with Xd’ Theoretical with Xd"
Simulation with H=10pu
1.05
Pdc (pu) 0.98
0.91
0.84
0.77
0.7
0.6 0.7 0.8 0.9 1 1.1 1.2
Idc (pu)
1.25
1.14
1.03
0.92
0.81
0.7
0.6 0.7 0.8 0.9 1 1.1 1.2
Idc (pu)
assumed reactance values (Xd and Xd ) are superposed on the simulated char-
acteristic. The simulated characteristic (obtained as described in section 3.1)
appears to follow the Xd based analytical result more closely at least at lower
DC current levels. After the peak is reached, the simulated curve deviates from
the Xd based characteristic and appears to be closer to the Xd based character-
istic. However it can be seen that the peak of the actual curve is more closely
approximated by the Xd based curve.
The oscillations in the simulated waveform, particularly in the area where
M AP is exceeded are due to the electrical and mechanical transients excited in
the model by the slow ramp-up procedure.
Z2 (local load)
Zf (filter)
Vd= 500
kV
Exciter
Vref
M _+
Detailed
M/c model
Zs
ESCR=1.597
Pd without Vt reg. Pd with Vt reg.
Vt without reg. Ideal Vt
3
Pdc (pu), Vt (pu)
2.5
2
1.5
1
0.5
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Idc (pu)
Figure 3.8. MPC and bus voltage characteristics with and without ideal excitation.
M AP now has a very large value and occurs at a very large DC current as shown
in Figure 3.8. The MPC and voltage characteristics for the earlier case (as in
Figure 3.2) are also superposed for comparison.
However, when a detailed simulation is carried out, this expected result does
not occur. The system follows the theoretical curve as seen in Figure 3.9 (solid
and dashed traces). Note that only a portion of the x-axis range of Figure
3.8 is plotted. The system becomes unstable at Id = 1.08 p.u. This is due to
instability caused by the excitation system dynamics, which were ignored in the
theoretical calculation. The exact onset level is also seen to vary depending on
the controller gains. This shows that the results obtained from a simple MPC
analysis could be excessively optimistic when excitation control is ignored.
If instead of modeling the machine in detail, a simplified exciter-machine
model as shown in Figure 3.10 is used, the results are in agreement with the
curve derived in Figure 3.8. In this approach the exciter-machine system is
modeled as a Proportional-Integral controller whose output is the Thevenin
voltage behind the machine impedance Xd . The resulting curves are shown as
the dotted lines in Figure 3.9. In this case the instability does not occur as
in the detailed model, but the curve flattens out after a current of 1.1 p.u. is
reached due to the attainment of a limit in the DC current controller. This
3.3. Conclusions 43
Thoretical Vt Simulation Vt
Simulation Vt (mach.&exciter)
Bus voltage, Vt(pu)
1.1
1
0.9
0.8
0.7
0.6
0.5
0.6 0.71 0.82 0.93 1.04 1.15
Idc(pu)
Zs Vref
Kp+Ki/s _+
Vpu
proves that the instability observed with the detailed machine model is indeed
due to dynamics which are ignored in the theoretical calculation of MPC.
3.3 Conclusions
The following observations can be made based on the above results:
impedance) do not follow the simulated MPC over the entire range. How-
ever, a reasonable figure for the M AP can be arrived at using Xd in the
theoretical derivation.
• When Modelling Synchronous Compensators in MPC analysis, the addi-
tional condition of P = 0 gives more accurate results.
• Inclusion of an exciter results in an actual response that tracks the M AP
result, but only up to a point. The margin of stability is not improved.
3.4. Appendix: Equations used in MPC Calculations 45
1. Test System without the inclusion of load models and with voltage sources
generating and absorbing the transmitted power for validation with a
Thevenin source AC model used in basic SCR (ESCR) calculation;
2. Constant P Q load model, as often used in steady state analyses but as-
sumed to be overly pessimistic in voltage stability studies;
47
48 4. Influence of System Operating Conditions and Converter Design
The approach has been to normalize the loads to provide the same P and Q at
the nominal operating point of the DC system while permitting the interaction
away from nominal to respond to above mentioned load models. The curves in
Figure 4.2 indicate that M AP and IMAP are sensitive to the assumption of load
model and the operation above rated DC current shows a spread in the MPCs.
For the Test System of ESCR of 1.26, the ESCR is less than the critical
value for which nominal operation coincides with maximum power. Figure 4.3
shows the MPCs for different load models, provided that the nominal operating
point for all load models used is the same.
1 2
3
8 11 12
5 6
10 9 13 7
1.2
1
Pdi & Vac (pu)
0.8
0.6
0.4
−. voltage dep. loads
* impedance loads
− no loads
0.2 − − PQ loads
0
0 0.2 0.4 0.6 0.8 1 1.2
Idi (pu)
1.2
0.8
Pdi (pu)
0.6
− no load
0.4 −. voltage dep. loads
− − impedance loads
* PQ loads
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
Idi (pu)
1.6
1.5
1.4
1.3
1.2
Vac (pu)
1.1
− no loads
1
−. voltage dep. loads
0.9 * impedance loads
− − PQ loads
0.8
0.7
0.6
0 0.2 0.4 0.6 0.8 1
Pdi (pu)
Figure 4.4. Voltage variations for different load models (ESCR=1.26) [9].
The discrepancies between the curves become more marked as the ESCR of
the inverter is reduced, as shown in Figure 4.3. A comparison of Figure 4.2 and
Figure 4.3 suggests that load representation is not a critical need in high SCR
situations but is progressively significant as the SCR reduces.
With constant current control rather than constant power control at the
rectifier, the load flow was able to converge in all four cases for a current larger
than and a common nominal operating point above was achievable. It is notable
from Figure 4.3 that the maximum available powers for the P Q and voltage-
dependent loads of 1.086 and 1.037, respectively, are higher than for the case
with no load included when the M AP was 1.01.
The four curves in Figure 4.4, so-called nose curves as commonly referred to
in voltage stability studies, show AC voltage variations at the inverter AC bus for
50 4. Influence of System Operating Conditions and Converter Design
1.6
1.4 AC voltage
1.2
MAP curve
1
Pdi & Vac (pu)
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
Idi (pu)
Figure 4.5. MPC and AC voltage for test system with no parallel AC
line (ESCR=1.69) [9].
different DC powers. The turning point between stable and unstable operation
(between dUac /dPd ≤ 0 on the upper branch and dUac /dPd ≥ 0 on the lower
branch) corresponds to the different M AP values in Figure 4.3. The spread
between AC voltages is particularly noticeable below the nominal operating
point of 1.0 pu for DC current and power. This suggests that provisions for
voltage control during start up, and scheduled changes in DC power, need to
take into consideration different load models in assessing stability issues.
1.4
1.2 o
C
o
B
1 o
A
0.8
Pdi (pu)
0.6
− Gamma=18 (TAP=1.23)
−. Gamma=18 (TAP=1.17)
... Gamma=22 (TAP=1.17)
0.4 − − Constant Vdc
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Idi (pu)
1.6
1.4
1.2
1
Vac (pu)
0.8
− Gamma=18 (TAP=1.23)
0.6
−. Gamma=18 (TAP=1.17)
... Gamma=22 (TAP=1.17)
0.4
* Constant Vdc (TAP=1.17)
− − Constant Vdc (TAP=1.23)
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
Pdi (pu)
Table 4.1. Data and for the systems which MPCs are given in Figures 4.8
to 4.11. It should be noted that the UT OV values given refer to complete
load rejection without taking any non-linear effects, e.g. transformer sat-
uration and arrester action, into account.
4.5 Conclusions
The study has shown the necessity of plotting maximum power curve and ob-
taining the maximum available power for the particular AC/DC system and
particular control mode at the inverter. It has also been shown that the basic
SCR calculation, while still very helpful in the first stage of design, is insuffi-
cient in itself for evaluating the power transfer capability and power stability
of the DC system. The presence and characteristics of AC loads as well as
the presence of any parallel AC lines influence the maximum available power
transfer and make the nose curves more diverse. Furthermore, the power trans-
fer capability is affected by the adopted strategy for tap changing. It was also
shown that CCC has superior performance when it comes to voltage and power
stability.
4.5. Conclusions 53
1.6
1.4
DC Power (p.u.) & AC−voltage (p.u.)
1.2
AC Voltage
0.8
0.6
MPC
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Current (p.u.)
Figure 4.8. Maximum Power Curve (MPC) of a system with CCC and SCR = 2.0.
1.6
AC Voltage
1.4
DC Power (p.u.) & AC Voltage (p.u.)
1.2
0.8
0.6
MPC
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Current (p.u.)
Figure 4.9. Maximum Power Curve (MPC) of a system with LCC and SCR = 2.0.
54 4. Influence of System Operating Conditions and Converter Design
1.4
DC Power (p.u.) & AC−voltage (p.u.)
1.2
0.8
0.6
0.4
MPC
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Current (p.u.)
Figure 4.10. Maximum Power Curve (MPC) of a system with CCC and
SCR = 0.93.
AC Voltage
1.6
1.4
DC Power (p.u.) & AC Voltage (p.u.)
1.2
0.8
MPC
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
DC Current (p.u.)
Figure 4.11. Maximum Power Curve (MPC) of a system with LCC and
SCR = 0.93.
5
Multi-Infeed HVDC Systems
55
56 5. Multi-Infeed HVDC Systems
sight into various potential problems arising from multi-infeed HVDC systems
was given. These include small-signal instability due to control interactions
among the constituent HVDC links, voltage instability and collapse, increased
commutation failures in one constituent HVDC link due to AC faults occurring
in the vicinity of a neighboring one, transient AC voltage depression due to
simultaneous recovery of constituent HVDC links after an AC fault. Similar to
[29], control coordination among the constituent HVDC links or modification of
the HVDC control modes were identified as potential solutions to overcome the
associated problems. In [32] a digital simulator study was performed to investi-
gate the transient response and recovery of multi-infeed HVDC systems based
on different control strategies of the constituent HVDC links. Transient recovery
of the studied system was found to be highly system dependent although certain
control strategies could provide robust system response. In [33] synchronizing
and damping torque modulation applied to the controls of constituent HVDC
links of a multi-infeed HVDC system were investigated. Cooperative coordi-
nation of controls among the constituent HVDC links was shown to provide
effective damping and synchronizing torque in a de-coupled manner, to achieve
overall system integrity with respect to electromechanical stability.
As seen, a broad spectrum of system phenomena could potentially arise in
multi-infeed HVDC systems [29], [30], [32], [33]. Conceivably these system phe-
nomena are related to those arising from single-infeed HVDC systems, since
multi-infeed system configurations historically evolved from the latter. More-
over, voltage/power instability problems clearly hold special interest since these
had been among the most important technical concerns for single-infeed HVDC
systems. In this respect it is appropriate to know how close these system phe-
nomena correspond, under what conditions they differ, and whether the ana-
lytical techniques developed for the single-infeed situation could be similarly
applied to the multi-infeed case. A host of other fundamental questions arise
and there is clearly a need to address them. These questions essentially provide
the main motivation for the work in [12] - [17], and [35] - [37]. Specifically these
works aim to achieve the following objectives;
1st
1st
2nd
N-th 2nd
N-th
3rd
system model suitable for systematic investigations should have the following
characteristics;
• A system topology that derives its roots from the classical single-infeed
HVDC system configuration. This strong analogy between the single- and
multi-infeed HVDC system configuration is essential to allow analytical
comparisons between them.
• A system model that is simple but representative of the multi-infeed
HVDC system configuration. This allows simulation studies based on
a realistic system model, at the same time theoretical analysis on it can
be done without being distracted by intractable mathematical rigour.
These desired characteristics motivated much use of the system configuration
as shown in Figure 5.1 for the work in [12] - [17], and [35] - [37]. Figure 5.1
shows two variants of this system configuration, namely the chain and ring
multi-infeed HVDC system configuration, but their essential characteristics are
58 5. Multi-Infeed HVDC Systems
that they comprise point-to-point HVDC links and all the converters are in-
verters. This is perceived to be the most common and onerous situation with
respect to voltage and power instability. Other system configurations and con-
verter type combinations are of course possible, such as shown in Figure 5.2
where the multi-infeed HVDC system configuration is constituted by multi-
terminal HVDC links. Ultimately a most general multi-infeed HVDC system
configuration would conceivably be as shown in Figure 5.3, being constituted by
point-to-point as well as multi-terminal HVDC links. This system configuration
could potentially come about as a result of the Baltic Ring, East-West Europe,
and southern Scandinavia HVDC projects.
As for converter type combinations, in general there are a number of possi-
bilities as categorized in [29], described below and schematically shown in Figure
5.4;
• Type 1A: converters of the same type connected to different HVDC sys-
tems.
• Type 1B: converters of the same type connected to a common HVDC
system.
• Type 2A: converters of mixed type connected to different HVDC systems.
5.3. Analysis of Voltage/Power Stability of Multi-Infeed HVDC Systems 59
Static Approach
Sensitivity Technique The underlying assumption for the static approach to
study voltage/power instability in HVDC systems is that the phenomenon oc-
curs so fast that the AC/DC system dynamics do not appreciably affect it. Thus
the DC line and converter pole control dynamics, and voltage control dynamics
in the AC system are neglected. The stability of the AC/DC system is thus es-
sentially determined by the algebraic system states that govern the power flow
in the system, and sensitivity of system states to small changes in controlling
system quantities is used as a measure of system stability.
Two main concepts based on the sensitivity technique had been proposed
in the past to analyze the voltage/power stability of weak single-infeed HVDC
systems. One approach known as the Maximum Power Curve (MPC) method is
based on the concept of Maximum Available Power (MAP ) and was introduced
in [5] and [6]. The other approach, known as the Voltage Stability Factor (V SF )
method, is based on the concept of voltage sensitivity and was first presented
in [7].
Taking cue from these concepts for the single-infeed HVDC system, they
were extended to analyze the voltage and power stability of multi-infeed HVDC
systems in [12], [17], [37]. This entails use of the system power flow model
and modal analysis technique. Briefly, this mathematical technique decomposes
the multi-dimensional AC power flow Jacobian into its modal equivalent. The
60 5. Multi-Infeed HVDC Systems
STATIC DYNAMIC
APPROACH APPROACH
Figure 5.5. Frameworks for voltage and power stability analysis of multi-
infeed HVDC systems.
resulting incremental modal voltages thus relate with the modal reactive power
injections in a de-coupled form and through an eigenvalue of the corresponding
mode. This is intuitively akin to the VSF of the single-infeed case and can be
interpreted as the modal voltage sensitivity factor MVSF for the multi-infeed
HVDC system.
Based on the above technique, a general and comprehensive method for volt-
age stability analysis of multi-infeed HVDC systems was developed in [12], [17],
[37]. The technique was further developed in [13] through use of participation
factors of the critical modes to determine the critical system location in a multi-
infeed HVDC system where the most severe system phenomena are expected to
occur, and for evaluating the effectiveness of remedial actions implemented.
Similarly the modal analysis technique incorporating the concept of max-
imum available power was proposed in [12] to investigate the power stability
of multi-infeed HVDC systems. It was shown in [12], [17] that the DC power
flow Jacobian could be decomposed into its modal equivalent and the resulting
incremental modal DC power and modal DC currents are related in a man-
ner similar to that for the single-infeed case. A limiting condition, akin to the
Maximum Available Power (MAP ) condition for the single-infeed case, is also
reached when the eigenvalue of the corresponding mode vanishes. Thus the
concept of Modal Maximum Available Power (MMAP ) for multi-infeed HVDC
systems could be similarly proposed.
Similar to the relationship between the power and voltage analysis methods
for the single-infeed case shown in [10], such relationships were also shown in
[12], [17] to exist for the multi-infeed HVDC case. It was also shown in [12],
[17] the inter-relationship between the methods for the single-infeed and multi-
5.3. Analysis of Voltage/Power Stability of Multi-Infeed HVDC Systems 61
infeed HVDC systems. From such relationship structure it was established that
the modal power/voltage analysis methods for multi-infeed HVDC systems are
the multi-dimensional equivalent of the corresponding ones for the single-infeed
case. In [16] it was further shown that these relationships are preserved when
load models are incorporated into the system models.
In [9] the influence of static load characteristics on the MPC was shown
for a HVDC link with a parallel AC line interconnection. Similar to this, the
influence of static load characteristics on the voltage/power stability of multi-
infeed HVDC systems was shown in [15], [16], [17]. In [15], [16], [17] the basic
voltage and power stability equations incorporating static load characteristics
were derived and used to show the effect of load and system parameters on
the voltage/power stability margins of multi-infeed HVDC systems. Conditions
for ascertaining the most unfavorable static load characteristics with respect to
degrading the voltage and power stability margins were derived in [16]. Indices
to characterize and estimate a bound for the effect of static load characteristics
on the system stability margins were also defined and studied in [16].
Dynamic Approach
The dynamic approach to voltage and power stability analysis of HVDC systems
considers the full or at least the most essential dynamics of the system. The
AC/DC system is thus represented by a nonlinear differential-algebraic system
model and solution techniques are applied according to the problem of interest,
that is whether the instability problem is due to a large or small-signal system
disturbance.
namely the saddle-node and Hopf bifurcation, were studied for the single-infeed
and multi-infeed HVDC systems. It was also shown in [17], [36] that the quasi-
static and dynamic conditions for voltage collapse considered via saddle-node
bifurcation in the HVDC system models considered were equivalent.
6.1 Background
In the actual application and installation of HVDC schemes in the last few
years, there is an impression that not so much attention has been paid to the
power/voltage stability phenomenon. This is probably because most, if not all,
recent HVDC schemes have been applied into systems where the inverter AC
system is relatively strong, that is, with a high Equivalent Short Circuit Ratio
(ESCR). In these cases, it is unlikely that power instability would be a problem
unless unusual operating conditions were imposed. However, this phenomenon
is still very real and a major issue whenever a conventional, line-commutated
thyristor-controlled HVDC inverter is applied into an AC system with a low
ESCR.
In the past, a number of DC schemes experienced the full effects of this
power instability, including in some cases the resulting consequences of system
collapse. In these systems, the problem has been mitigated by invoking special
supplementary controls on the DC or by a combination of controls and special
compensation considerations [4], [42], [43].
In recent years, new developments in HVDC technology provide the potential
to significantly reduce any manifestation of the power stability problem in most
systems. Developments such as capacitive-commutated converters (CCC) and
voltage-source inverters (VSI) now make it theoretically possible to apply DC
links into much weaker AC systems without these problems. Both the controls
and the parameters incorporated into these new technologies are such that they
present positive impacts towards the prevention or mitigation of instability.
The application development of these technologies is still progressing and a full
treatment of their theory, operation and system effects are beyond the scope
of this document. However, it is important to note that in these cases the
basic principles of the stability phenomenon do not change. But the modes of
operation and limits of DC application can change.
63
64 6. Planning and Analysis of Systems - Principles and Concepts
terms, the basic mechanism of the instability occurs as follows: Suppose the AC
voltage drops just a small amount; then the DC voltage at the inverter drops
proportionately (they are directly related). But the main function of the DC
controls in many DC schemes is to hold constant power order (Pd = Ud · Id ).
If the DC voltage (Ud ) goes down, the only way the controls can maintain
constant power is to increase the DC current (Id ). However, as soon as it does
this, the reactive power consumption of the DC increases because this is directly
related to the DC current. Unless the additional reactive power is immediately
(very important - this will be explained later) available, the AC voltage will dip
further and start the cycle all over again. In other words, an unstable cycle will
be started and perpetuated which, unless it is stopped in some way, can lead to
a runaway condition and a total collapse of the entire AC/DC system.
system stays constant; that is, before any dynamic compensation or volt-
age regulators have a chance to respond. From the AC system viewpoint,
the critical value of ESCR and hence the maximum available DC power
should then be dependent on the Thevenin voltage behind reactance and
the equivalent system reactance in the classical case of AC system steady
state stability. Real systems, however, do not act as fixed Thevenin equiv-
alents over long times beyond a few hundred milliseconds. The Thevenin
voltage in this case is not fixed. The general effect can be instability oc-
curring in the form of relaxation oscillations at a low frequency, generally
less than 5 Hz, depending on the response lags of machine voltage regu-
lators correcting AC voltages, and on communication and master control
lags in the DC The criterion for stability, however, will be practically the
same as for a fixed Thevenin equivalent. If the DC control is for constant
frequency, then the relaxation oscillation frequency will be lower because
of system inertia effects, but the stability criterion will again be similar.
The criteria for power/voltage instability addressed above does not say anything
about electro-mechanical stability because reactive power flow is the critical
quantity, while in electro-mechanical stability, phase-angle shifts and frequency
are predominately governed by active power flow. In addition, it is only valid for
fundamental frequency. Steady state stability of the AC voltage is a prerequisite
for normal system operation. However, the problem can become worse if the AC
system is weak in the inertia sense, that is, if the total system inertia constant
is so low as to cause substantial frequency drop during a fault or other event.
With everything fixed in time and condition (including controls), any at-
tempt to operate that fixed system above its MAP, or, conversely, any attempt
to weaken the system when it is operating at its MAP, implicitly implies the
onset of voltage collapse in that system. The effects that the real, non-fixed,
system controls have on stability will be treated in subsequent sections below.
6.5. Relationship Between DC Stability and Recovery from AC Faults 67
While a Udhold control can be a good back-up protection, there are other
quasi-static or single action reactive controls, that are more effective in arresting
a potential collapse, but still within the DC controls. This is by a SUVC (System
Undervoltage Control) or by a PRBC (Power Runback Control). In principle,
once the AC voltage falls to a certain level, these controls reduce or ramp-down
the DC power by a set amount. This not only stops the collapse but it assists
the AC voltage to recover by releasing reactive power back into the AC system,
thereby returning the system back to a satisfactory operating state (it closes
the door again).
The controls described above can only have limited capability to protect the
system, and this is by design. In SUVC or PRBC controls, if the DC power
were to be reduced by a large amount whenever an AC undervoltage occurs,
then instability or collapse could almost always be prevented. But this would
be generally unsatisfactory for system operation and for system requirements.
In any system, the critical question will revert to determining how severe a con-
tingency condition, reflected on the AC bus, that these controls must cater for
and protect against. As examples, in many systems a determining contingency
may be the loss of a large filter, capacitor bank or other large reactive power
compensation device, including possible loss of nearby generating machines.
where a DC system has older machines with slow exciters in close proximity,
then this is a different matter.
ESCR(X’d)
4
ESCR
ESCR(Xd)
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Relative Gain of Voltage Regulator
Figure 6.1. Curves showing the onset of voltage instability for various
excitation gains on the synchronous compensators connected to the in-
verter AC bus (Nelson River HVDC system), as a comparison against the
critical ESCR based on the classic steady state stability equation.
collapse problem. The first is the System Undervoltage Control (SUVC), which
performs a relatively minor (about 300 MW) reduction when a collapsing AC
voltage is detected. The second is a Udhold control, which substitutes a fixed
DC voltage into the constant-power-division circuit upon low AC voltage. While
Udhold should prevent a total collapse initiated by compensators on a limit, the
violation of the criterion can be severe enough to cause a large change in AC
voltage.
The main benefit of the SUVC is not in the small increase in ESCR from
the reduction in DC power, but rather in the return of controllability to the
compensators. Thereafter, the contribution of the compensators to the fault
level in the ESCR calculation is based on Xd instead of Xd . This is a dramatic
increase in ESCR, which now makes the system stable. The SUVC is thus
a dominant control in preventing voltage instability on this system for many
possible conditions.
action before more restrictive limits are applied. To ensure a fast-acting control
system, the effective time constant of the voltage measurement circuit of the
AVR should be low, perhaps 25 milliseconds or less. For modern AVR designs,
this should not be a problem. The field-forcing positive-limit may have some
affect on stability for more severe disturbances by quickly returning the AC
voltage to a safe level. Finally, the contribution of the machine to the fault
level should be maximized from the perspective of stability, meaning the lowest
realistic total of subtransient reactance plus transformer reactance. Of course,
the down side of this is the possible implication on the fault rating of system
circuit breakers and other equipment.
From all the foregoing discussion, the impression has no doubt been made
that the worst single problem associated with incorporating an HVDC link into
low ESCR weak-connections or weak AC systems may well involve the interre-
lated issues of DC power stability and reactive compensation at the inverter.
Also, these issues are major problems from both technical and economic per-
spectives. All this applies to conventional line-commutated converters.
For relatively large power DC links, any requirement for sophisticated com-
pensation equipment such as synchronous compensators or SVCs will represent
major capital investment - a significant proportion of the total link cost. This
increases the difficulty for conventional HVDC to compete against evolving new
technologies in AC transmission systems and other supply alternatives. In some
areas, it will also increase the difficulty to compete against alternative energy
solutions that may also be aided by the restructuring and deregulation in the
electrical industry.
If new DC technology, such as capacitive-commutated-converters or voltage-
sourced converters, can eliminate the stability/compensation problems and be-
come more economic, this will greatly enhance the possibility of new, large
power HVDC links being applied in relatively weak or moderate AC systems.
References
[1] Cigré Task Force 38-02-12, “Criteria and Countermeasures for Voltage Col-
lapse”, Paris 1995.
[3] Hill, D.J., Hiskens, I.A., and Popovic, D., “Load Recovery in Voltage Sta-
bility Analysis and Control”, Bulk Power System Phenomena III, Voltage
Stability, Security and Control, ECC, Davos, 1994, pp. 579 - 595.
[4] Cigré Working Group 14.07, “Guide for Planning DC Links Terminating
at AC Systems Locations Having Low Short-Circuit Capacities, Part I:
AC/DC Interaction Phenomena”, Cigré Guide no. 95, Paris 1992.
[5] Ainsworth, J.D., Gavrilovic, A., Thanawala, H.L., ”Static and Synchronous
Compensations for HVDC Transmission Converters Connected to Weak AC
Systems”, CIGRE General Session, Paris, France, Paper 31-01, 1980.
[7] Hammad, A.E., Sadek, K., Kauferle, J., “A New Approach for the Analysis
of and Solution of AC Voltage Stability Problems at HVDC Terminals”,
Proceedings of International Conference on DC Power Transmission, Mon-
treal, Canada, pp.164-170, Jun.1984.
[8] Nayak, O.B., et.al, “Control Sensitivity Indices for Stability Analysis of
HVDC Systems”, IEEE Transactions on Power Delivery, Vol.10 No.4,
pp.2054-2060, Oct.1995.
[9] Reeve, J., and Uzunovic, E., “Study of Power Transfer Capability of DC
Systems Incorporating AC Loads and a Parallel AC Line,” IEEE Trans. on
Power Delivery, Vol. 12, No. 1, pp. 426 - 432, January 1997.
77
78 REFERENCES
[13] Lee, H.A.D., Andersson, G., “Use of Participation Factors in Modal Voltage
Stability Analysis of Multi-Infeed HVDC Systems”, IEEE Transactions on
Power Delivery, vol.1 no..13, pp. 203-211, Jan. 1998.
[17] Lee, H.A.D., “Voltage and Power Stability of HVDC Systems”, Ph.D The-
sis, Royal Institute of Technology, Sweden, TRITA-EES-9801, ISSN 1100-
1607, March 1998.
[19] Szechtman, M., et al.,(1991): “First Benchmark Model for HVDC Control
Studies” Electra 135. 1991. 54-73
[20] Reeve, J., “Short Circuit Ratio and Related Factors Characterizing AC
Systems for DC Transmission,” International Conference on Large High
Voltage Electric Systems - CIGRE, Colloquium “Interaction Between AC
and DC Systems”, England, September 1985.
[22] Lefebre, S., Gole, A. M., Reeve, J., Pilotto, L., Martins, N., and Bhat-
tacharya, S., “Report on Test Systems for AC/DC Interaction Studies,”
IEEE Trans. on Power Delivery, Vol. 10, No. 4, pp. 2027 - 2034, October
1995.
REFERENCES 79
[23] Jonsson, T., and Björklund, P., “Capacitor Commutated Converters for
HVDC”, Proceedings of Stockholm Power Tech, Part: Power Electronics,
pp. 44 -51, Stockholm, 1995.
[24] Diago, J.R. et al, “Garabi the Argentina - Brazil 1000 MW Interconnec-
tion”, Proceedings of VIII ERLAC, Cidade del Peste, Paraguay, 1999.
[27] Krishnayya, P.C.S., Piwko, R.J., Weaver, T.L., Bahrman, M.P., Hammad,
A.E., “DC Transmission Terminating at Low Short-Circuit Ratio Loca-
tions”, IEEE Transactions on Power Delivery, Vol.PWRD-1 No.3, pp.308-
318, July 1986.
[28] Thallam, R. S., “Review of the Design and Performance Features of HVDC
Systems Connected to Low Short-Circuit Ration AC Systems”, IEEE
Transactions on Power Delivery, Vol.7 No.4, pp.2065-2073, October 1992.
[29] Bui, L. X., Sood, V. K., Laurin, S., “Dynamic Interactions between HVDC
Systems Connected to AC Buses in Close Proximity”, , IEEE Transactions
on Power Delivery, Vol.6 No.1, pp.223-230, January 1991.
[30] Szechtman, M., et. al, “The Behavior of Several HVDC Links Terminating
in the Same Load Area”, Cigre General Session, Paris, France, Paper 14-
201, 1992.
[31] Lips, H. P., “Aspects of Multiple Infeed of HVDC Inverter Stations into a
Common AC System”, IEEE Transactions on Power Apparatus and Sys-
tems, Vol.92 No.2, 1973.
[32] Reeve, J., Lane-Smith, S. P., “Multi-Infeed HVDC Transient Response and
Recovery Strategies”, “, IEEE Transactions on Power Delivery, Vol.8 No.4,
pp.1995-2001, October 1993.
[33] Pilotto, L. A. S., et. al, “Synchronizing and Damping Torque Modulation
Controllers for Multi-Infeed HVDC Systems”, IEEE Transactions on Power
Delivery, Vol.10 No.3, pp.1505-1513, July 1995.
[34] Ärnlöv, B., et. al, “Aspects of Cooperation of Three HVDC Links between
Southern Denmark. Southern Sweden and Northern Germany/Jutland”,
Cigre Symposium, Tokyo, Paper 220-02, 1995.
[35] Lee, H.A.D., Andersson, G., “Impact of Dynamic System Modeling on the
Power Stability of HVDC Systems”, IEEE Transactions on Power Delivery,
Vol. 14, No 4 (1999), pp. 1427 - 1437
80 REFERENCES