02 Whole
02 Whole
by
Liying WANG
Doctor of Philosophy
in
School of Electrical and Electronic Engineering
Faculty of Engineering,
Computer and Mathematical Sciences
The University of Adelaide, Australia
August 2013
Copyright 2013
Liying WANG
All rights Reserved
THE UNIVERSITY
of ADELAIDE
Typeset in Word2010
Liying WANG
ABSTRACT
Abstract
The primary aim of this thesis is to investigate the small-signal dynamic performance
of high voltage direct current (HVDC) transmission links based on voltage source
converter (VSC) technology operating in parallel with the existing longitudinal Australian
power system. This thesis presents the principle design methodology to achieve robust
controllers for VSCs including inner current controller, outer power and voltage controllers
as well as the supplementary damping controllers for enhancing the small-disturbance
rotor-angle stability of a weak multi-machine power system with embedded VSC-HVDC
links.
Three types of linear current controller schemes (proportional-integral, proportionalresonant and Dead-Beat schemes) are investigated and discussed in detail to identify the
most suitable control method. Due to its wider bandwidth and superior performance under
unbalanced operating conditions, the Dead Beat current controller is set as the inner current
controller that has not been analysed in detail in the literature.
A new methodology for the selection and optimization of the parameters of the
proportional-integral compensators in the various control loops of a VSC-HVDC
transmission system using a decoupled control strategy is also proposed in this thesis. It
was found that the new methodology is effective in a relatively strong system. However,
since the method did not take various operating conditions and system disturbances into
account, it will not be effective in a relatively weak system. The analysis shows that the
ABSTRACT
design of robust outer loop controllers is challenging due to the limited bandwidth of the
inner current controller in a weak AC system. Therefore, the second primary objective of
the project was to develop a simple fixed parameter controller, which can perform well
over a wide range of operating points within the active/reactive power (PQ) capability
chart of the VSCs. To achieve this second objective, various grid conditions including
various Short Circuit Ratios (SCRs), different X/R ratios and PQ capabilities of the VSC
system were studied.
To support the primary objectives, a detailed higher order small-signal model of the
DB controlled VSC is developed and systematically verified. As an original contribution,
the study developed a new methodology to linearize the modulator/demodulator blocks
which are used to develop the small signal models for several key components such as the
sampling block, the delay block and the DB inner current controller.
The initial values of the PI/PID compensator parameters are obtained by applying the
classical frequency response design methods to a set of detailed linear models of the openloop transfer functions of the VSC-HVDC control system. It was concluded that an
iterative process may be required after examining the co-operation performance of these
controllers designed.
In the final chapter of this thesis, the small-signal rotor-angle stability of a model of
the Australian power system with embedded VSC based HVDC links was examined. For
the analytical purposes of this thesis a simplified model of the Australian power system is
used to connect the high capacity, but as yet undeveloped, geothermal resource in the
region of Innamincka in northern South Australia via a 1,100 km HVDC link to Armidale
in northern New South Wales. It is observed that the introduction of the new source of
geothermal power generation has an adverse impact on the damping performance of the
system. Therefore, two forms of stabilization are examined: (i) generator power system
stabilisers (PSS) fitted to the synchronous machines which are used to convert geothermal
energy to electrical power; and (ii) power oscillation damping controllers (PODs) fitted to
the VSC-HVDC link. In the case of the PODs two types of stabilizing input signals are
considered: (i) local signals such as power flow in adjacent AC lines and (ii) wide-area
signals such as bus voltage angles at key nodes in the various regions of the system. It was
concluded that the small-signal rotor-angle stability of the interconnected AC/DC system
has been greatly enhanced by employing the designed damping controllers.
ii
STATEMENT OF ORIGINALITY
Statement of Originality
This work contains no material which has been accepted for the award of any other
degree or diploma in any university or other tertiary institution and, to the best of my
knowledge and belief, contains no material previously published or written by another
person, except where due reference has been made in the text.
I give consent to this copy of my thesis when deposited in the University Library, being
made available for loan and photocopying, subject to the provisions of the Copyright Act
1968.
I also give permission for the digital version of my thesis to be made available on the
web, via the Universitys digital research repository, the Library catalogue and also
through web search engines, unless permission has been granted by the University to
restrict access for a period of time.
______________________________
______________________________
Signed
Date
iii
ACKNOWLEDGEMENTS
Acknowledgements
Firstly, I would like to express my sincere appreciation to my supervisor, Associate
Professor Nesimi Ertugrul for all your help, support, guidance, time, and, most of all, the
lessons you have taught me. Special thanks to Mr David Vowles, my co-supervisor, for his
excellent guidance, insightful conversations and endless encouragement throughout the
duration of the research. Working side-by-side with him was an honour and a privilege. I
also extend my gratitude to Professor Boon-Teck Ooi, Mr Jian Hu and Dr Lianxiang Tang
for their assistance in understanding and implementing the Dead-Beat control algorithm.
I am also grateful to the China Scholarship Council (CSC) and the University of
Adelaide (UA) for their financial support of this work. In addition, my thanks go also
towards all the members in the school of Electrical and Electronic Engineering of the
University of Adelaide, in particular Dr Nicolangelo Iannella, Ahmed Abdolkhalig and
Qiming Zhang for their great help with thesis writing, interesting discussions and
providing data analysis tool for producing the technical eigenvector compass plots.
Finally, I would like to thank my parents for the invaluable support they provided
during the period of research. Special thanks to my beloved husband and daughter, Xike.
Thank you so much for your endless love, support and understanding. Thank you, Yaoyi,
for the joy and encouragement you have brought to my life.
iv
TABLE OF CONTENTS
Table of Contents
Abstract ................................................................................................................................. i
Statement of Originality...................................................................................................... iii
Acknowledgements ............................................................................................................. iv
Table of Contents ................................................................................................................. v
List of Figures ................................................................................................................... viii
List of Tables .................................................................................................................... xix
List of Publications ........................................................................................................... xxi
Symbols ........................................................................................................................... xxii
Acronyms ........................................................................................................................ xxiv
Chapter 1:
Introduction ........................................................................................................................ 1
1.1
Background ................................................................................................................ 1
1.2
VSC-HVDC System................................................................................................... 2
1.2.1
1.2.2
1.2.3
1.3
1.3.1
1.3.2
1.3.3
1.4
Thesis Overview....................................................................................................... 15
1.4.1
TABLE OF CONTENTS
1.4.2
Chapter 2:
2.1
Introduction .............................................................................................................. 18
2.2
2.2.1
2.2.2
2.3
2.3.1
2.3.2
2.3.3
2.4
Conclusion ................................................................................................................ 56
Chapter 3:
3.1
3.2
3.2.1
3.2.2
3.2.3
3.2.4
3.2.5
3.2.6
3.2.7
3.2.8
3.2.9
3.3
Comparison of the Small Signal Linear Model and PSCAD Simulation ............... 109
3.4
3.4.1
3.4.2
3.4.3
3.5
Chapter 4:
vi
4.1
4.2
4.2.1
4.2.2
TABLE OF CONTENTS
4.3
4.3.1
4.3.2
4.3.3
4.3.4
4.3.5
4.4
Chapter 5:
Stability of VSC-HVDC Links Embedded with the Weak Australian Grid .................. 165
5.1
5.2
5.2.1
Admittance Matrix Representation of the Integrated Grid and Filter ............ 168
5.2.2
5.2.3
5.3
5.3.3
5.4
Chapter 6:
Conclusion...................................................................................................................... 212
6.1
6.2
6.2.1
6.2.2
6.3
vii
LIST OF FIGURES
List of Figures
Figure 1-1 A typical VSC topology connected to the AC grid (a) simplified diagram,
(b) the details of the circuit including three phase two-level VSC using IGBTs................... 3
Figure 2-1 Current control classifications .................................................................. 19
Figure 2-2 Representation of rotating vector in converter including stationary
reference frame and abc natural reference frame. .......................................................... 21
Figure 2-3 Representation of rotating vector in dq reference frame and eference
frame .................................................................................................................................... 22
Figure 2-4 Single line diagram representation of VSC-HVDC .................................. 22
Figure 2-5 (a) The diagram of the inner and outer controllers; (b) Inner current
control loop. ......................................................................................................................... 25
Figure 2-6 Flow chart of the process for determining the PI compensator parameters
for VSC-HVDC control system ........................................................................................... 27
Figure 2-7 Detailed inner current control loop in pu system ...................................... 28
Figure 2-8 Block diagram of DC voltage control scheme in pu system..................... 29
Figure 2-9 (a) Block diagram of active power control scheme in pu system; (b) block
diagram of reactive power control scheme in pu system. .................................................... 30
Figure 2-10 VSC-HVDC system diagram .................................................................. 32
viii
LIST OF FIGURES
Figure 2-11 Open-loop Bode plots of the current controller transfer function (i)
without PI control (blue solid line); (ii) with PI compensation using initial parameters
(green dashed line). .............................................................................................................. 33
Figure 2-12 Step responses of current controller transfer function (i) without PI
control (blue solid line); (ii) with PI compensation using initial parameters (green dashed
line) ...................................................................................................................................... 34
Figure 2-13 Open-loop Bode plots of the outer DC controller transfer-function (i)
without PI compensation (solid blue line); and (ii) with PI compensation using the initial
parameters (dashed green line). ........................................................................................... 35
Figure 2-14 Step responses of the outer DC controller transfer-function (i) without PI
compensation (solid blue line) and (ii) with PI compensation using the initial parameters
(dashed green line). .............................................................................................................. 35
Figure 2-15 Open-loop Bode plots of the outer active/reactive power controller
transfer-function (i) without PI compensation (blue solid line) and (ii) with PI
compensation using the initial parameters (dashed green line) ........................................... 36
Figure 2-16 Responses of VSC-HVDC following a reversal of the power order for (i)
the initial set of PI compensator parameters (red line); (ii) the parameters obtained with the
optimization function (2-34) (blue line) and (iii) the parameters obtained with the refined
optimization function (2-35) (cyan line). The reference values of the variables are shown in
green line.............................................................................................................................. 37
Figure 2-17 Structure of the paralleled harmonic compensators ............................... 40
Figure 2-18 Diagram of reference frame mathematical model of VSC system .... 41
Figure 2-19 The simplified model of PR current controlled VSC system in
reference frame .................................................................................................................... 41
Figure 2-20 The categories of DB current controllers ............................................... 42
Figure 2-21 Digital DB current controlled VSC system ............................................ 43
Figure 2-22 The closed-loop DB current Control ...................................................... 44
Figure 2-23 Pole-zero map of the one sample delay system: red cross: one sample
delay DB without considering the computation delay time; blue cross: represents one
sample delay DB considering the computation delay time; green cross: represents reducing
the proportional gain; ........................................................................................................... 45
Figure 2-24 One sample delay DB control with Smith Predictor ............................ 46
ix
LIST OF FIGURES
Figure 2-25 Internal model control design for DB implementation block ................. 47
Figure 2-26 The block diagram of solving feedback transfer function DB current
control .................................................................................................................................. 48
Figure 2-27 The structure of the solving feedback transfer function DB current
controller .............................................................................................................................. 48
Figure 2-28 Frequency response and step response of various DB current controllers
.............................................................................................................................................. 49
Figure 2-29 PSCAD step response simulation results of the DB current controller .. 54
Figure 3-1 Small signal model of VSC with the discrete DB current controller........ 60
Figure 3-2 Structure of phase locked loop.................................................................. 62
Figure 3-3 Non-linear model for PLL ........................................................................ 65
Figure 3-4 Block diagram of PLL linearized model................................................... 66
Figure 3-5 Frequency response of PLL in the case study........................................... 68
Figure 3-6 Linearized model verification, 1% step change in phase.......................... 68
Figure 3-7 50% step (180o) change in phase .............................................................. 69
Figure 3-8 1% step change in frequency .................................................................... 70
Figure 3-9 PLL responses to 90% magnitude step change of input voltage .............. 70
Figure 3-10 PLL output phase angle in comparison with input phase angle when
being subject to a sudden 90% magnitude step change of input voltage ............................. 70
Figure 3-11 PLL output frequency in comparison with input frequency when being
subject to a sudden 90% magnitude step change of input voltage ....................................... 71
Figure 3-12 PLL output phase angle when being subject to an A-phase to ground
fault ...................................................................................................................................... 71
Figure 3-13 PLL output frequency when being subject to an A-phase to ground fault
.............................................................................................................................................. 72
Figure 3-14 Relationship between converter reference frame and grid RI reference
frame .................................................................................................................................... 73
Figure 3-15 Controlled voltage source representation with an inductance ................ 74
Figure 3-16 Test results for grid RI reference frame transformation ......................... 74
Figure 3-17 The test circuit for small signal mathematical equation ......................... 75
Figure 3-18 Small signal test results for the grid frame transformation ..................... 75
Figure 3-19 Current compensation block verification ................................................ 77
LIST OF FIGURES
Figure 3-40 Small signal model for DB current controller in equivalent dq reference
frame .................................................................................................................................... 98
Figure 3-41 Compare Yd/Yq outputs from (i) the modulator/de-modulator system; (ii)
the small-signal equivalent transfer function and iii) the small-signal equivalent state space.
(a)UdSTEP = 0, UqSTEP = -0.02 and (b) UdSTEP = -0.02, UqSTEP = 0; ....................................... 99
Figure 3-42 Difference between Yd/Yq outputs from (i) the modulator/de-modulator
system versus the small-signal equivalent transfer function; (ii) the modulator/de-
xi
LIST OF FIGURES
modulator system versus the small-signal equivalent state space. (a) UdSTEP = 0, UqSTEP = 0.02 and (b) UdSTEP = -0.02, UqSTEP = 0; .............................................................................. 99
Figure 3-43 Compare Yd/Yq outputs from the (i) modulator/de-modulator system; (ii)
the small-signal equivalent transfer function and iii) the small-signal equivalent state space
............................................................................................................................................ 102
Figure 3-44 Test circuit for the detailed VSC model ............................................... 106
Figure 3-45 The equivalent circuit diagram of grid ................................................. 107
Figure 3-46 Small signal model of the grid .............................................................. 108
Figure 3-47 (a) Grid side current and (b) filter bus voltage responses of i) large signal
model ii) small signal model following a input voltage step change order on VcqC = 0.56kV. ............................................................................................................................... 109
Figure 3-48 Comparison of the converter output current ICcd , ICcq using the (i) large
signal model simulating with VSC simplified as a controlled voltage source; (ii) the smallsignal model for AC system with different SCRs. ............................................................. 111
Figure 3-49 Comparison of the converter output currents ICcd , ICcq from the (i) large
signal model using the detailed VSC model, and (ii) the small-signal model for the weak
AC system .......................................................................................................................... 112
Figure 3-50 The model DC transmission link .......................................................... 113
Figure 3-51 Test circuit for DC link test .................................................................. 114
Figure 3-52 The converter currents behind the DC capacitor from both of the rectifier
side and inverter side (circle: Large signal model; rectangular: Small signal model) ....... 114
Figure 3-53 Demonstration of sub-modular interconnection ................................... 115
Figure 3-54 The base model for the converter side controller ................................. 117
Figure 3-55 Comparison of the converter output currents ICcd , ICcq using the (i)
large-signal model simulated in the detailed VSC model including the DC link, and (ii) the
small-signal model for a weak AC system together with the DC link ............................... 118
Figure 4-1 Single phase equivalent circuit ............................................................... 122
Figure 4-2 Simplified DB current controlled VSC ................................................... 123
Figure 4-3 Zeros/poles in z-domain for (a) weak system (b) strong system of a DB
current controlled VSC system .......................................................................................... 124
xii
LIST OF FIGURES
Figure 4-4 (a) Bode plots and (b) step responses of a DB current controlled VSC
system ................................................................................................................................ 124
Figure 4-5 Flow chart of the controller design methodology ................................... 125
Figure 4-6 PQ capability chart (a) demonstration chart; (b) defined operating points
........................................................................................................................................... 126
Figure 4-7 Single-line diagram of AC side converter including a high-pass filter .. 129
Figure 4-8 Influence of the DC link model: (a) pole-zero map of Iq(s)/Iqref(s); (b)
Bode plots of Iq(s)/I qref(s) .................................................................................................. 129
Figure 4-9 Basic structure for outer power controller design .................................. 130
Figure 4-10 (a) In the DB current controlled VSC with P=Pmax (a) eigenvalue maps
with SCR reduced from 7.5 to 0.5 at a step of -1 (i.e. A to G); (b) Bode plots of the TF
(blue: SCR=7.5, red: SCR=4.5, cyan: SCR=2). ................................................................ 133
Figure 4-11 The time-domain simulation results: PSCAD SCR=7.5 (solid cyan line);
Matlab SCR=7.5 (dash-dot blue line); PSCAD SCR=2 (solid black line); Matlab
SCR=2(dash-dot red line). ................................................................................................. 134
Figure 4-12 Four main low frequency modes for A: rectifier mode with a 90o angle;
B: rectifier mode with a 75o angle; C: rectifier mode with a 60o angle; a: inverter mode
under 90o angle; b: inverter mode under 75o angle; c: inverter mode under 60o angle; .... 135
Figure 4-13 The Bode plots of the DB current controlled VSC. blue: rectifier mode
under 90o angle; cyan: Rectifier mode under 75o angle; red: rectifier mode under 60o angle;
magenta: inverter mode under 90o angle; black: inverter mode under 75o angle; green:
inverter mode under 60o angle; .......................................................................................... 135
Figure 4-14 Step responses of DB current controlled VSC for the rectifier (left) and
the inverter (right). rectifier mode at angle 90o (blue); rectifier mode at angle 75o (cyan);
rectifier mode at angle 60o (red); inverter mode at angle 90o (magenta); inverter mode at
angle 75o (black); inverter mode at angle 60o (green); ...................................................... 136
Figure 4-15 For DB controlled VSC with SCR=2 (a) the eigenvalue map with P
equals to be 1pu (A), 0.5pu (B) and 0pu (C); (b) Bode plots of the TF with different
loading conditions. blue: P=1pu, cyan: P=0.5pu, red: P=0pu). ......................................... 137
Figure 4-16 (a) The eigen-value map and (b) Bode plot of the transfer function for
the DB current controlled VSC at SCR=2 and different power factors (Q=Qmax, Q=0 and
xiii
LIST OF FIGURES
Q=Qmin) under full load condition. (dark blue: lagging power factor, light blue: unity
power factor, red: leading power factor). ........................................................................... 138
Figure 4-17 The analysis of the transfer functions using Bode plots and at 36
operating conditions in total (a) SCR=2 at 90o (b) SCR=2 at 75o; (c) SCR=7.5 at 90o; (d)
SCR=7.5 at 75o; ................................................................................................................. 139
Figure 4-18 Right half-plane zeros for the transfer function of the open loop power
controller ............................................................................................................................ 141
Figure 4-19 Bode plots of the selected cases for the power controller design ......... 142
Figure 4-20 The Bode plots with (a) low-pass filter; (b) DB current controlled VSC
with a low-pass filter; (c) designed PI compensator; (d) open loop transfer function of
P(s)/Pref(s). .......................................................................................................................... 143
Figure 4-21 Step responses of P(s)/Pref(s) ................................................................ 143
Figure 4-22 Control block for the power controller ................................................. 144
Figure 4-23 A typical MIMO system, in which the hidden feedback loop is shown in
red lines .............................................................................................................................. 144
Figure 4-24 Bode plots of selected cases in the AC voltage controller design ........ 146
Figure 4-25 The Bode plots of (a) low-pass filter; (b) Uf(s)/Idref (s) of DB current
controlled VSC with a low-pass filter; (c) the designed PI compensator and (d) open loop
transfer function of Uf(s)/Ufref(s) at the rectifier side. ........................................................ 147
Figure 4-26 Step responses of Uf(s)/Ufref(s) ............................................................. 148
Figure 4-27 Control block for AC voltage controller ............................................... 148
Figure 4-28 The frequency responses of (a) Udc(s)/Iqref(s); (b) Udc(s)/Idref(s); (c)
Uf_inv(s)/ Iq_ref_inv(s) and (d) Uf_inv(s)/Id_ref_inv(s)
variation of the rectifier side converter (36 scenarios), but with a constant operating point
at the inverter side converter that is PmaxQ0, SCR=2 and with angle equals to 90o. .......... 149
Figure 4-29 (a) The Bode plots of the selected cases of the DC voltage controller
design Udc(s)/Iqref_inv(s); (b) designed PID compensator (Kp_dc+Ki_dc/s+Kd_dcs/(1+TfUdcs)=0.024-0.088/s+0.0009s/(1+0.06s)); (c) frequency responses of open loop transfer
functions of Udc(s)/ Udc_ref (s) and (d) step responses of the closed-loop Udc(s)/ Udc_ref(s).
............................................................................................................................................ 151
Figure 4-30 The control block diagram of the DC voltage controller ...................... 151
xiv
LIST OF FIGURES
Figure 4-31(a) Bode plots of selected cases for inverter end AC voltage controller
design Uf_inv(s)/Idref_inv(s); (b) designed filter+PI compensator; (c) frequency responses of
open loop transfer functions of Uf_inv(s)/Uf_inv_ref(s); (d) step responses of closed-loop
Uf_inv(s)/ Uf_inv_ref(s)............................................................................................................ 153
Figure 4-32 1% step of AC voltage reference at the rectifier end; (a) responses of
power at the rectifier end; (b) responses of the filter bus voltage Uf_rec at the rectifier end;
(c) responses of DC voltage Udc_inv at the inverter end; (d) responses of the filter bus
voltage Uf_inv at the inverter end. ....................................................................................... 155
Figure 4-33 1% step of power reference at the rectifier end; (a) responses of the
power at the rectifier end ; (b) responses of the filter bus voltage Uf_rec at the rectifier end;
(c) responses of the DC voltage Udc_inv at the inverter end; (d) responses of the filter bus
voltage Uf_inv at the inverter end. ....................................................................................... 156
Figure 4-34 1% Step of the AC voltage reference at the inverter end; (a) response of
the power at the rectifier end; (b) response of filter bus voltage Uf_rec at the rectifier end; (c)
response of the DC voltage Udc_inv at the inverter end and (d) response of the filter bus
voltage Uf_inv at the inverter end. ....................................................................................... 157
Figure 4-35 1% step of the DC reference voltage at the rectifier side; (a) the
responses of power at the rectifier; (b) the responses of the filter bus voltage Uf_rec at the
rectifier; (c) the responses of the DC voltage Udc_inv at the inverter end; (d) the responses of
the filter bus voltage Uf_inv at the inverter end. .................................................................. 158
Figure 4-36 The step responses of Ufinv_ref (a) and Udcinv_ref (b) and their induced
performances in their cross-coupling loops as the system working at an operating point
PmaxQmax, SCR=2 with =90o. red: Matlab 1%; cyan: PSCAD 1%; black: PSCAD 2%/2;
magenta: PSCAD 4%/4; green: PSCAD 8%/8. ................................................................. 161
Figure 4-37 The step responses of Udcinv_ref (a) and Ufinv_ref (b) and their induced
performances in their cross-coupling loops as the system working at an operating point
P0Qmin, SCR=7.5 with =75o. red: Matlab 1%; cyan: PSCAD 1%; black: PSCAD 2%/2;
magenta: PSCAD 4%/4; green: PSCAD 8%/8. ................................................................. 162
Figure 5-1 The diagram of the extended Simplified South-East Australian power grid
with VSC-HVDC links. ..................................................................................................... 167
Figure 5-2 Grid model as an admittance .................................................................. 168
xv
LIST OF FIGURES
Figure 5-3 1% DC voltage step responses of the dynamic and grid admittance models
(a) for rectifier side; (b) for inverter side. .......................................................................... 170
Figure 5-4 The grid and filter model as an admittance............................................. 171
Figure 5-5 1% DC voltage step response comparison of the admittance representation
of grid model only, and the gird and filter models adopting admittance representation: (a)
the rectifier side; (b) the inverter side. ............................................................................... 172
Figure 5-6 System parameter scaling scheme ......................................................... 173
Figure 5-7 Scaling of inner current controllers ........................................................ 174
Figure 5-8 Scaling of PI controller with additional filter ......................................... 176
Figure 5-9 Scaling of the PID controller .................................................................. 176
Figure 5-10 Flow chart of the verification methodology of modeling ..................... 177
Figure 5-11 Eigenvalue map of the original and scaled systems ............................. 178
Figure 5-12 The comparison of the system eigenvalues between the new updated DC
link system and the original system ................................................................................... 179
Figure 5-13 The frequency responses for the rectifier side (a) and for inverter side (b)
for the higher order grid impedance models (I to V). ........................................................ 181
Figure 5-14 The open loop frequency responses, P_rec/P_ref ..................................... 182
Figure 5-15 The step responses test of the power reference at the rectifier side ...... 183
Figure 5-16 The open loop frequency responses, Uf_ref_rec/Uf_rec .............................. 183
Figure 5-17 The step response test of the AC voltage reference at the rectifier side
............................................................................................................................................ 184
Figure 5-18 The open loop frequency responses, U_dc_ref_rec/U_dc _rec ....................... 184
Figure 5-19 The step response test of the DC voltage reference at the inverter side
............................................................................................................................................ 185
Figure 5-20 The open loop frequency responses, Uf_ref_inv/Uf_inv ............................ 185
Figure 5-21 The step response of the AC voltage reference at the inverter side ...... 186
Figure 5-22 Eigenvalue map of the interconnected system ...................................... 187
Figure 5-23 Right eigenvector prototype of mode I40 ............................................. 188
Figure 5-24 Participation factor for mode I40 .......................................................... 189
Figure 5-25 Right eigenvector prototype of mode I35 ............................................. 190
Figure 5-26 Participation factor for mode I35 .......................................................... 190
Figure 5-27 Right eigenvector prototype of mode I25 ............................................. 191
xvi
LIST OF FIGURES
xvii
LIST OF FIGURES
Figure 5-42 The performance evaluation of the VSC controllers with the
supplementary controllers in service following a step change of 0.01 pu on power reference
of the rectifier side converter. ............................................................................................ 206
Figure 5-43 Frequency response of lead compensator Q(s) of the WAPOD ........... 207
Figure 5-44 The eigenvalue map of system fitted with WAPOD under high loading
condition, where the letters represent the gains of the WAPOD, as KSS is increased from
3.64210-4 to 0.12 with a step size of 0.02. ........................................................................ 207
Figure 5-45 The eigenvalue map of system fitted with WAPOD under light loading
condition, where the letters represent the gains of the WAPOD, while KSS is increased
from 3.64210-4 to 0.12 with a step size of 0.02. ................................................................ 208
Figure 5-46 Eigenvalue analysis of WAPOD with PSS under high loading condition,
where KSS is increased from 3.642 10-4 to 0.12 with a step size of 0.02 .......................... 209
Figure 5-47 Power outputs of Innamincka generator #1 following a step change of
0.01 pu on voltage reference of Innamincka generators. ................................................... 209
Figure 5-48 The performance evaluation of the VSC controllers with PSS and
WAPOD in service following a step change of 0.01 pu on power reference of the rectifier
side converter. .................................................................................................................... 210
xviii
LIST OF TABLES
List of Tables
Table 1-1 Categories of SCR........................................................................................ 8
Table 2-1 PI parameters comparison .......................................................................... 38
Table 2-2 Open loop transfer function of DB current controllers .............................. 49
Table 2-3 Sensitivity to plant parameters ................................................................... 50
Table 2-4 Simulation parameter ................................................................................. 51
Table 2-5 Summary of specific data obtained from Figure 2-29 ............................... 54
Table 2-6 Mean derivation of simulation results........................................................ 56
Table 3-1 PLL parameters .......................................................................................... 67
Table 3-2 Characteristics of the linearized model ...................................................... 69
Table 3-3 Characteristics of the step response corresponding to the pure one sample
delay ..................................................................................................................................... 81
Table 3-4 Characteristics of the step response corresponding to ZOH ...................... 85
Table 3-5 Summary of transfer functions for DB current controller, sampling block
and delay block .................................................................................................................... 97
Table 3-6 Equivalent dq reference frame transfer function for the lead-lag block .... 98
Table 3-7 State space equations for DB controller in equivalent dq reference frame
........................................................................................................................................... 100
Table 3-8 Simulation results for VSC model verification ........................................ 107
xix
LIST OF TABLES
Table 3-9 Summary of the large signal model of the AC grid in abc reference frame
and small signal model in dq reference frame ................................................................... 108
Table 3-10 Large and small signal for DC link ........................................................ 113
Table 3-11 Parameters for the DC link test .............................................................. 114
Table 4-1 Equivalent Thevenin impedance of the AC grid ...................................... 127
Table 4-2 Basic power flow equation ....................................................................... 127
Table 4-3 Equations for the power flow study ......................................................... 128
Table 4-4 Modes of the inner current controller (SCR=2) ....................................... 131
Table 4-5 Participation factors of dominant state variables for selected modes of the
inner DB current controller ................................................................................................ 132
Table 5-1 The units derivation for the PID coefficients ........................................... 177
Table 5-2 Inter-area modes of the integrating system .............................................. 186
Table 5-3 Approximate improvements on system damping: comparison of the results
obtained from adding with equivalent damping torque (Figure 5-31) and the analysis for
participation factor. ............................................................................................................ 194
Table 5-4 Residue analysis of the system ................................................................. 202
xx
LIST OF PUBLICATIONS
List of Publications
Published
[1] W. Liying and N. Ertugrul, "Selection of PI compensator parameters for VSC-HVDC
system using decoupled control strategy," in Universities Power Engineering Conference
(AUPEC), 2010 20th Australasian, pp. 1-7.
[2] W. Liying, N. Ertugrul, and M. Kolhe, "Evaluation of dead beat current controllers for
grid connected converters," in Innovative Smart Grid Technologies - Asia (ISGT Asia),
2012 IEEE, 2012, pp. 1-7.
Papers in Preparation
[1] W. Liying, David J. Vowles and N. Ertugrul, "Reference Frame Transformation
Approach for Small Signal Modeling of VSC with Stationary Frame Controllers"
[2] W. Liying, David J. Vowles and N. Ertugrul, " Generalized small signal modeling of
DB controlled VSC"
[3] W. Liying, David J. Vowles and N. Ertugrul, "Robust controllers design for DB
controlled VSC linked to a weak AC system"
[4] W. Liying, David J. Vowles and N. Ertugrul, "Damping controller design of DB
controlled VSC operating in parallel with a Weak Multi-machine AC Power System "
xxi
SYMBOLS
Symbols
xxii
system matrix
active power
ac
base value
capacitance
current
De
damping torque
Kd
derivative gain
dc
feedthrough matrix
inductance
input matrix
SYMBOLS
Ki
integral gain
Laplace factor
max
maximum value
min
minimum value
output matrix
peak
peak value
phase angle
abc
phase quantities
PLL
PLL quantity
Kp
proportional gain
reactive power
ref
reference value
Ts
db
Kv
voltage
xxiii
ACRONYMS
Acronyms
xxiv
AC
Alternating Current
AVM
AVR
DB
Dead Beat
DC
Direct Current
DG
Distributed Generation
EMT
Electro-Magnetic Transient
FIR
GM
Gain Margin
GPS
GTOs
H Infinity
HSV
IGBT
Im
Imaginary
IMC
ITAE
ACRONYMS
KCL
KVL
LCC
Line-Commutated Control
LL
Lead-Lag
LMI
LTI
Max
Maximum
MIMO
Multi-Input Multi-Output
Min
Minimum
MTDC
PCC
PM
Phase Margin
PI
Proportional-Integral
PLL
POD
PQ
Active-/Reactive- Power
PR
Proportional-Resonant
PSS
pu
Per Unit
PWM
Re
Real
RHP
SCR
SISO
Single-In Single-Out
SVCs
TF
Transfer Function
VCO
VSC-HVDC
WAPOD
ZOH
xxv
Chapter 1: Introduction
1.1 Background
As is known, the limited availability and environmental concern of fossil fuels, as well
as the continuous growing demand of electricity, have caused renewable energy to become
commercially attractive. However, the integration of large scale renewable energy sources
into the power grid changes the characteristics of the supply in terms of the composition of
energy generation, the transmission network, the technology and the economics of the
electricity industry. This will be very different from the current situation in which fossilfuel generation occupies the dominant position [1].
To allow the large scale integration of renewable energy sources into the grid, VSCHVDC is a commonly developed power electronic converter. This converter topology is
particularly suitable for long distance power transmission, such as in offshore wind farms.
Since VSCs do not require commutating voltage from the connected AC grid, they are also
highly effective in supplying power to isolated and remote loads and convenient when
interconnecting distributed sources.
CHAPTER 1: INTRODUCTION
The existing transmission networks and power stations will not be able to meet the
growing challenges of renewable energy transmission. The increasing penetration of power
electronic converters into existing power systems have also caused several power system
stability issues [2, 3]. The modern VSC-HVDC technology can help to solve the adverse
impact on system stability and power quality issues. Furthermore, it can also significantly
increase the transmission capacity and can provide flexible control of power flow.
Therefore, it is envisaged that there is a strong need to investigate the design, operation and
control of VSC-HVDC transmission links. This includes the links integration into the
existing power network which may be notably weak as in the Australian case which is
investigated in this thesis.
Constant Voltage
T
Constant Current
L
C
VSC
S2
Filter
S1
(a)
Pdc
U c c
P,Q
Ra
Usa
n
T1
T3
idc +
T5
La
Usb
Rb
Usc
Rc
ia
Lb
Uca
ib
U dc
Ucb
ic
Lc
Ucc
T4
U s s
DC Line
T6
T2
C
DC Line-
System Control
Converter
Control
Firing Control
Inner Current
Control
Outer Control
Other Control
(b)
Figure 1-1 A typical VSC topology connected to the AC grid (a) simplified diagram, (b)
the details of the circuit including three phase two-level VSC using IGBTs.
U sU c
sin
XL
U (U U c cos )
Q s s
XL
(1-1)
CHAPTER 1: INTRODUCTION
outputs, which consequently regulate P and Q. This allows us to operate the VSC in four
quadrants.
The VSC-HVDC technology can also be applied in a wide range of applications due to
the above mentioned characteristics [4, 5, 10-13]. For example, it can
The VSC-HVDC projects around the world are summarized in [8] and [14]. Among
these the Hellsjn-Sweden prototype project 3MW, 10kV was the first VSC-HVDC
transmission project, which was commissioned in 1997. The worlds first HVDC project
employing VSCs in a modular multilevel converter (MMC) topology was the Trans Bay
Cable Project (TBC) of USA. It has a transmission capacity 400MW and 200kV of DC
voltage rating which was the highest rating for this type of technology until 2010 [15].
After that, the Caprivi link interconnector was commissioned in Oct 2010. This was the
first ABBs HVDC-Light transmission system (350kV, 300MW) to employ overhead line
which is 950km long. The latest milestone is the Skagerrak 4 link from Denmark to
Norway, for which the voltage rating is 500kV, the highest to date [16]. Note that it will be
used as a reference for determining the DC voltage rating in this project.
CHAPTER 1: INTRODUCTION
control irrespective of the direction of power flow. The reactive power at either end of the
link is controlled separately by the respective converters. The control of reactive power is
used to control reactive power directly or indirectly as a means of controlling the powerfactor or AC voltage at a designated bus. Due to the simplicity and robustness, double
closed-loop vector oriented PI controllers have been utilized in compensating the system to
achieve the desired performance. The inner-loop controllers employ feed-forward
decoupled control to make the active- and reactive-current track the reference values
produced by the outer-loop controllers [25].
1.3.1.2 Controller Design for VSC Operating Under Strong AC System Condition
The selection of the parameters of the PI compensators of the various VSC control
systems is a key issue to ensure an adequate dynamic performance (including a fast
response, and a sufficient stability margin for the VSC-HVDC transmission system). Note
that the dynamic performance criteria typically results in conflicting requirements when
tuning the PI controller parameters. Consequently, a compromise between the speed of
response and stability for small disturbances is needed. Further compromises may also be
needed to achieve adequate performance in response to large disturbances. The need to
coordinate the tuning of several PI compensators is particularly challenging. Typically, the
PI controller parameters of the various compensators are obtained by trial-and-error which
relies heavily on the experiences and skills of the design engineers.
Very few publications were identified that quantify the selection and optimization of PI
parameters associated with VSC control systems. The approach adopted in [26] applies a
trial and error method with the objective of ensuring the transient stability of the system
without using any theoretical evaluation criteria. In [27] frequency response analysis is
used to obtain an envelope of PI compensator parameter values which satisfy specified
stability criteria. Then an objective function is evaluated for each set of PI parameters in
the above envelope using a detailed electro-magnetic transients (EMT) model. The set of
PI parameters which minimizes the objective function is selected. The approach proposed
in [27] was found to be promising but requires a significant amount of computation time
because there is no mathematically derived procedure for selecting the next set of PI
parameters on the basis of accumulated experience gained from the previous trials [28]. In
addition, it does not provide any insight on how to adjust several PI parameters
simultaneously. Although the papers in [29] and [30] apply optimization techniques based
CHAPTER 1: INTRODUCTION
on the Simplex Algorithm to different VSC-HVDC control systems, the selection of the
initial set of parameters for input to the optimization algorithm is ad-hoc. In addition, the
method tends to find a local optimum rather than the global optimum set of PI compensator
parameters. References [31] and [32] propose approaches for determining the PI
compensator parameters, but these methods do not include the simultaneous adjustment of
several PI parameters. Therefore, it is desirable to develop a new methodology to
determine these parameters.
1.3.1.3 Controller Design for VSC Operating Under Weak AC System Condition
It is noted that there are only a few research papers that consider the impacts of VSCHVDC transmission systems connected to a weak grid. It is identified that the outer power
loop, the inner current loop, the synchronization method and the input-output impedance
can be studied in such grids that accommodate the VSC-HVDC system. The strength of the
AC system to which the link is connected is determined by the short circuit ratio (SCR).
The SCR is a ratio of the AC-system short-circuit capacity divided by the rated power of
the HVDC link. The AC system strength is classified in [33] and given below in Table 1-1.
SCR3
Weak System
3 > SCR 2
SCR<2
is considered. In a weak system, the voltage is sensitive to converter current, whereas the
sensitivity can be ignored when the grid is strong. Moreover, the dynamics of the PLL and
the bus filter also become significant in the dynamic performance of the system.
The literature survey has also considered how to improve the design and performance
characteristics of the various outer loop controllers. It was reported in [35] that the flexible
AC voltage control can be realized through optimum power and reactive power
management to offer additional voltage support to reduce the impact of fast varying filter
bus voltage. A non-linear eight state model of the VSC was developed in [36] by taking
into account the dynamics of the PLL and the bus filter to accommodate the connection to
weak AC grids. Based on this multi-input multi-output (MIMO) model, a robust H infinity
(H) controller over a range of operating points is designed by employing linear matrix
inequality (LMI) techniques [36, 37]. Two genetic algorithms based on direct search
methods are also applied in [38] to optimize this type of controller.
B. Inner Loop Control
A challenging feature for the design of controllers associated with distributed
generation (DG) systems is high grid impedance and resistance. Hence, the advanced
techniques developed in the DG research area provide useful insights to controller design
for VSC-HVDC system [39-42]. As discussed below it was suggested in previous studies
that the instability problems associated with the connection of VSC-HVDC systems to
weak AC grids can be divided into two categories: i) high frequency instability which can
be yielded due to grid-converter resonance (e.g. LCL filter resonance) and interaction
between the inner current control and the active damping controller [41, 43]; ii) low
frequency instability which is highly influenced by the inner current control loop.
a) High Frequency Instability
Mitigation of the high frequency instability can be achieved by a number of alternative
active damping control schemes which include: i) the utilization of poorly damped
complex poles [43-45], ii) split capacitor [46, 47], iii) full state feedback control [48] and
iv) virtual resistance [39-41, 49]. With the employment of active damping control, the open
loop gain and bandwidth of the inner current controller are consequently increased and
current harmonic distortion is reduced. As a result, the increased bandwidth facilitates the
design and tuning of the outer-loop controller parameters.
CHAPTER 1: INTRODUCTION
10
11
CHAPTER 1: INTRODUCTION
12
13
CHAPTER 1: INTRODUCTION
The effectiveness of the POD design depends on a number of factors including the
location of controlled plant (relates to controllability); the modal content of the selected
stabilizing signal (relates to observability); and the measurement accuracy of the selected
signal. In [96], to damp out the inter-area mode in the Hydro-Quebec network, two
different approaches on feedback signal selection were investigated. They were i) the
geometric approach (also known as controllability/observability index) and ii) the residue
approach. It was concluded that the geometric approach is more reliable and useful in this
application, rather than the residues approach which is widely used. Moreover, it was also
concluded that the global signal input is more effective than the local signal. However, it
should be noted that, this is not always true. Differences in performance between local
POD and WAPOD vary depending on several factors including: i) the actual tuning of
those PODs, ii) the availability of measurements with good modal content, iii) the power
system configuration and disturbances, iv) the location of VSC and v) the specific
operating point [89].
Later, in addition to the aforementioned controllability/observability and residue
methods, the Hankel singular value (HSV) approach is added to complete the comparison
in [97]. It was concluded that based on a small two-area system, the residue and Hankel
singular value methods perform in a similar way. However, in a medium scale system, the
system damping performance was better when signal selection was based on the HSV
method. There are also a number of studies on investigations of the design approaches of
PSS and POD, which are documented in [98-103].
Recently, several investigations on enhancing system stability through POD are
published based on the Hydro-Qubec network [96], the Nordic power system [89], the
reduced Great Britain model [104], the China Southern Power Grid (CSG) [105], and the
New England test system integrated with the New York power system [94]. In fact,
Australia also has an important track record in utilizing PODs for damping inter-area
modes of oscillation. References [98] and [100] formed the basis of PODs fitted to SVCs
in Brisbane that were required to stabilize the interconnection between NSW and
Queensland [106]. Furthermore, improved damping controls fitted to the SVCs in the
Brisbane area were investigated in [107] in the mid to late 2000s. The new controls were
based on local bus frequency measurements.
14
A simple and effective methodology for the controller design and optimization of
VSC-HVDC systems connected to weak grids is not available.
Stability studies involving the DB controlled VSC-HVDC links are only based on
time-domain analysis in the previous studies. The small signal analysis, which is
essential to characterize the system characteristics, of such discrete power system
has not yet been reported due to the difficulty that the time varying nature prevents
the direct application of small signal studies.
15
CHAPTER 1: INTRODUCTION
The possible interactions between VSC-HVDC links and the extended weak
Australian grid have not been reported.
Finally, the compensation methodologies (PSS, POD and WAPOD) for a secure
operation of VSC-HVDC links embedded in an extended weak Australian grid also
require further investigation.
16
Chapter 6: The chapter summarizes the conclusions drawn in each chapter and also
reiterates the original contributions.
17
2.1 Introduction
In a typical VSC-HVDC system, the hierarchical control structure is conventionally
adopted including system control, converter control and firing control from top to the
bottom, as shown in Figure 1-1b. Here, only two levels of the converter control that the
outer power-, voltage- control and the inner current control are discussed. The main
function of the outer loop controller is to control voltage or power of the three-phase
converter. Yet, the inner current control is mainly responsible for current control according
to current references that are generated by the outer loop controllers. For example, to
18
2.1. INTRODUCTION
realize sinusoidal current wave tracking with a unity power factor by using inner current
control. Since the inner current loop control is much faster than the outer loop control,
hence, the inner current control loop plays a significant role in improving the overall
performance of the control system.
Usually the adopted outer loop control contains a) constant DC voltage control; b)
constant real power control; c) constant reactive power control; d) constant frequency
control and e) constant AC voltage control. Generally speaking, a, b and d are a group of
incompatible control objects, which means a VSC station can only work in either of the
three control modes. In order to maintain the real power balance and DC voltage stability,
there must be one VSC station adopts constant DC voltage control, and others can work in
b or d mode. For the same reason, the c and e are also two incompatible control modes.
Note that one VSC station can only realize one control mode.
Converter control is the core strategy of the VSC-HVDC system, and it is therefore
being set as the research emphasis of this chapter.
In accordance with the linearity properties, the control method can be classified into
two categories: i) linear control scheme and ii) non-linear control scheme [61]. The linear
controllers can also be classified as the conventional decoupled PI-controller in the
synchronously rotating reference frame, the PR controller in the stationary reference frame,
and the predictive controller involving feed-forward and DB control scheme, as shown in
Figure 2-1. All these linear controllers are the focus of this chapter and will be explained
further in the following section.
PI
Linear
Predictive
Proportional
Resonant
Current Control
Method
Feed-Forward
Dead-Beat(DB)
Fuzzy
Non-Linear
Passivity
Hysterisis
19
equation (2-1), where up and PLL denote the peak value and initial phase angle of the
rotating vector respectively.
ua (t ) u p sin[0 (t )t PLL (t )]
2
ub (t ) u p sin[0 (t )t PLL (t ) ]
3
2
uc (t ) u p sin[0 (t )t PLL (t ) ]
3
(2-1)
The adopted Clark Transformation matrix for transformation from the converter
stationary abc natural reference frame to
[85].
1
1
1 2 2
ua
u
3
3
u 2 0
ub
3
2
2
uc
u0
1
1
1
2 2
2
(2-2)
u
sin[0 (t )t PLL (t )]
u u p
cos[0 (t )t PLL (t )]
(2-3)
Therefore, the positive sequence voltage for the balanced three phase supply voltage
can be represented using the stationary reference frame and abc natural reference frame
as given in equation (2-4) and shown in Figure 2-2 [110], where the scaling of the factor K
20
can be equal to 2/3 to ensure power invariance, alternatively it can also be chosen to
2/3
j
j
U u (t ) ju (t ) K (ua (t ) ub (t ) e 3 uc (t ) e 3 )
(2-4)
The inverse voltage transformation equations are also given below in equation (2-5),
0
1
ua
3
u 1
b 2
2
uc
3
1
2 2
1
2
u
1
u
2
u0
1
2
0 t
PLL
(2-5)
'
reference
Reference Frame
Reference to the relative motion theory, it can be noted that the dq reference frame
rotates at the same frequency 0 with the rotating vector . Therefore observing from
vector it will appear stationary, and will contain a set of DC terms which will facilitate
the controller design and modeling of VSC system. As shown in Figure 2-3, the q-axis is
chosen aligning with the opposite direction of the rotating vector, and d-axis is orthogonal
to the q-axis, lagging it by 90o. This transformation can be given in a vector form as in
equation (2-6), and in matrix form as in equation (2-7). Equation (2-8) and equation (2-9)
provide the reverse transformation from dq reference frame to
21
udq (t ) u (t )e j (0 (t )t PLL (t ))
(2-6)
u (t ) udq (t )e j (0 (t )t PLL (t ))
(2-7)
(2-8)
(2-9)
d
0 t
PLL
0 t
'
q
Figure 2-3 Representation of rotating vector in dq reference frame and reference frame
iL
idc
L A
B R
Us
is
Udc
Uc
-
22
three-phase AC voltage at the VSC terminals (i.e. point A in Figure 2-4) uc_abc are able to
be given by equation (2-10) [112-114].
disa
L
Risa
sa
ca
dt
disa
Risb
U sb U cb L
dt
disc
U sc U cc L dt Risc
(2-10)
Due to the DC link capacitor C, the converter DC voltage (Udc) and currents (idc and iL)
are related, which can be given by equation (2-11),
C
dU dc
idc iL
dt
(2-11)
U s U c L s Ris
dt
U U L dis Ri
s
c
s
dt
(2-12)
Since the fact that the relationship between and dq reference frames is determined
by equation (2-8) .Then, the voltage and current equations in reference frame becomes,
U s U sdq e jt
jt
U c U cdq e
i i e jt
s sdq
(2-13)
Substituting equation (2-13) into equation (2-12) yields to the following set of equations,
U s dq e
jt
U cdq e
jt
d (isdq e jt )
U s dq e jt U cdq e jt e jt L
U s dq e jt U cdq e jt e jt L
dt
disdq
dt
disdq
Risdq e jt
Lisdq
d (e jt )
Risdq e jt
dt
(2-14)
j Lisdq e jt Risdq e jt
dt
If we extend these equations to the matrix form and apply the Laplace Transformation,
then VSC systems mathematical model can be given in the dq reference frame by equation
(2-15).
usd sL R L isd ucd
u
sL R isq ucq
sq L
(2-15)
23
R
1
L
L
si R i i 1 (u u )
sq
sq
sd
sq
cq
L
L
(2-16)
Let us assume that the VSC output voltage is determined by the following feed-forward
decoupled PI controller [113],
u ( K KiI )(i* i ) Li u
cq
ip
sq
sq
sd
sq
(2-17)
where K ip and K iI are proportional and integral gains, i*sd and i*sq are the current
references of the d,q axes respectively, and isd and isq are the measured converter output
currents in terms of the d-, q- axis respectively.
Substituting equation (2-17) into equation (2-16),
24
[ R ( Ki p iI )] / L
0
*
isd 1
sisd
KiI isd
s
)
( Kip
si
i
K
L
s isq*
sq
0
[ R ( Ki p iI )] / L sq
s
(2-18)
From the above equation, it can be easily seen that, d- and q- axis current id and iq are
successfully decoupled. Hence, the control strategy structure using the above discussions
can be illustrated as in Figure 2-5b.
Converter
is_abc
P,Q
R
icap
L
u s_abc
PI
PI
U dc_ref /Pref
abc-dq
Q
+-
isd_ref
isq_ref
isd
PWM
isq
u c_abc _ref
cd _ref
Inner
Current
dq-abc
Control u cq _ref
Loop
u sq
u sd
U dc /P
dc
u c_abc
is_abc
Q ref
iL
idc
abc-dq
us _ abc
PLL
(a)
u sd
isd_ref
+
+
isd
isq
+
isq_ref
PI
u cd_ref
u cq_ref
PI
+
u sq
(b)
Figure 2-5 (a) The diagram of the inner and outer controllers; (b) Inner current control loop.
25
Figure 2-5a shows a three-phase, two-level VSC adopting a voltage oriented vector
control structure implemented in the synchronously rotating dq-reference frame, where the
three-phase AC voltage at the network node of the VSC and the three-phase current
flowing from the network into the VSC are denoted by us_abc and is_abc respectively, C is the
capacitance of the DC filter, and the entire control system including the outer- power and
voltage loop, the inner current control loop together with phase locked loop (PLL) control
are also given. The error of Pref and P is the input to the sending end controller, while the
error of Udc_ref and Udc is the input to the receiving end controller. Note that, only the
positive sequence current is the control target considered. The reference value of the
converter terminal voltage uc_abc_ref is generated from the control system then fed into the
VSC-HVDC system through the firing control that PWM control. Lastly, it should be
noted that the reference phase angle of control system is extracted from the filter bus
voltage by the PLL control.
2.3.1.1 Controller Design and Optimization for PI Controlled VSC-HVDC System
As is well known, the performance of a system highly depends on the performance of
its controllers, which was emphasised here about the importance of controller design. This
section proposes an approach to the selection and optimization of the parameters of the PI
compensators in the various control loops of a VSC-HVDC transmission system using a
decoupled control scheme.
In this study, an optimization algorithm based on the simplex method is adopted. The
main objective of this method is to simultaneously minimize the weighted sum of the
integral of the time absolute-error products (ITAE) of the active power, the reactivepower, the DC voltage and the inner current controllers of the respective VSCs. The initial
values of the PI compensator parameters for input to the optimization algorithm are
obtained by the application of classical frequency response design methods to simplified
linear models of the open-loop transfer functions (TFs) of VSC-HVDC control system.
The optimization process is applied to a detailed electromagnetic transient (EMT) model of
the VSC-HVDC system to which a large disturbance is applied.
The effectiveness of the optimized PI compensator parameter settings are assessed in
terms of rise-time, overshoot, and settling times for a range of disturbances applied to the
detailed EMT model. On the basis of these assessments, modifications to the weightings of
26
Convergence ?
Y
End
Figure 2-6 Flow chart of the process for determining the PI compensator parameters for
VSC-HVDC control system
In the below paragraphs the design parameters of the controllers are explained
systematically.
A. Transfer Function of Control Loops
a) Inner Current Control loop
The inner current loop control shown in Figure 2-7 is composed of 4 parts: i) PI
compensator
(Hc1(s)=(Kp1+Ki1/s));
ii)
equivalent
simplified
filter
model
27
F1 (s) H c1 s Gp1 s
(2-19)
K PWM
1
1
1
X ( pu )
1 Ts s 1 Ts s R( pu ) s L
2
(2-20)
Here, Ts is the sampling interval of the inner current control loop, KPWM is the
equivalent gain of the PWM, which is normally defined by KPWM= 2 uac/Udc, where uac is
the peak value of the terminal line-to-ground voltage of the converter.
isd.max/sq.max
i dref/qref
-+
PI Controller
PWM
u dis
PI (1,2)
K PWM
T
1+s s
2
-+
-isd.max/sq.max
Load
i d/q
1
R(pu)+s
X L (pu)
1
1+sTs
equivalent
inner
current
control
loop
(G 4 (s)=1/(1+Teq s)) , where Teq is the equivalent time constant of first order approximation
of current control loop [31, 32]; (iii) line circuit (G5 (s)=1/(s Cpu )) ; (iv) current
transmission relationship according to power balance between AC side and DC side
(G6 (s)=3usq.pu /2Udc.pu ) and (v) measurement circuit (G7 (s)=1/(1+Ts s)) [115]. The
compensated open loop transfer function of the outer DC voltage control loop can then be
written as
F2 (s) H c 2 s Gp 2 s
where the plant transfer function G p2 (s) is
28
(2-21)
G p 2 ( s) G4 ( s) G5 ( s) G6 ( s) G7 ( s)
2 U dc.pu
isq.max
PI controller
Vdcref
PI(3)
+-
-isq.max
usq. pu
1
1
1
Teq s 1 U dc. pu s C pu Ts s 1
+-
i qref
1
Teq s+1
isq.pu
(2-22)
I L.pu
+
s Cpu
2 U dc.pu
Vdc
1
1+sTs
(2-23)
According to Figure 2-3, the q-axis is aligned with the opposite direction of Us phasor.
Therefore, usd equals to be 0. The power equation can be simplified as,
Ps 3 / 2 usq isq
Qs 3 / 2 usq isd
(2-24)
The above equation shows that the active power and reactive power can be controlled
independently by controlling the q axis current isq and the d axis current isd.
The outer- active and reactive power control loops are shown in Figure 2-9. For initial
design purposes the network voltage can be assumed to be fixed with a value of one perunit, i.e. u sq =1.0 .
The outer- active/reactive power control loops is composed of 4 parts: i) PI
compensator
equivalent
inner
current
control
loop
(G8 (s)=1/(1+Teq s)) ; (iii) power balance transmission relationship (G9 (s)=3/2u sq (pu) and
(iv) measurement circuit (G10 (s)=1/(1+Ts s)) . The compensated open loop transfer function
of the outer active/reactive power control loop is shown in equation (2-25),
F3 (s) H c3 s Gp3 s
(2-25)
29
isq.max
isqref (pu)
Pref (pu)
+-
PI(4)
-isq.max
isq (pu)
1
1+Teq s
P(pu)
3
u sq (pu)
2
1
1+s Ts
(a)
isd.max
Q ref (pu)
+-
isdref.pu
PI(5)
-isd.max
1
isd.pu 3
u sq (pu)
1+Teq s
2
Q pu
1
1+s Ts
(b)
Figure 2-9 (a) Block diagram of active power control scheme in pu system; (b) block
diagram of reactive power control scheme in pu system.
where the plant transfer function G p3 s is
1
3
1
usq ( pu)
Teq s 1 2
Ts s 1
(2-26)
F3 (s) ( K p Ki / s)
1
3
usq ( pu)
4Ts s 1 2
(2-27)
30
the rise time must be less than 30 ms for underdamped system or less than 70
ms for a damped system;
J ITAE t edt
0
(2-28)
Where t is the time since the disturbance is applied, |e| is the absolute value of the
control system error and T is a finite time chosen so that the integral approaches to a
steady-state value and is usually chosen as the settling time Tf.
For a VSC-HVDC transmission system which includes multiple control loops, the
objective is to minimize the weighted sum of ITAE indices Of (X) associated with each of
the control systems [28].
31
O f ( X ) mi fi ( X )
(2-29)
fi ( X ) j t e j (t ) dt
(2-30)
AC power grid 1
Converter 1
Pc1 Qc1
Pc1 Qc1
U s11
U c10
R1
Converter 2
DC Cable
I dc1 I L1
I L 2 I dc 2
L1
AC power grid 2
Pc 2 Qc 2
Qc 2 Pc 2
U c 2 0
U s 2 2
L2
R2
I cap 2
I cap1
Phase Reactor
Phase Reactor
AC filter 1
U dc1
U dc 2
32
F1 ( s) ( K p1 Ki1 / s)
0.78509
1
1
4
4
1 3.7 10 s 1 7.4 10 s 9.5493 104 s 0.0015
0.78509K p1 s 0.78509Ki1
(2-31)
The frequency response of the plant transfer function Gp1(s) in equation (2-20) is
calculated and displayed in Bode plot as shown in Figure 2-11 (solid blue line). According
to the requirements listed in section 2.3.1.1 B, it can be seen that the phase margin of 39.4
deg is smaller than the minimum requirement of 40 deg. For the uncompensated system the
gain margin is 11.4dB, the overshoot of the step response is 30%, and the settling time is
0.01s, all of which are acceptable (see solid blue line in Figure 2-12). PI compensation is
applied (Kp1=0.62, Ki1= 53) to reduce the gain cross-over frequency thereby increasing the
phase margin to 50.4 deg and the gain margin to 17.1dB (dashed green line in Figure
2-11). The overshoot of 23.6% (dashed green line in Figure 2-12) is lower than the original
30%. Although the settling time has been increased to 0.025s, it is still acceptable.
Magnitude(dB)
200
without PI
with PI
100
0
-100
-200
-4
10
10
-2
10
10
10
10
Phase(deg)
0
-100
-200
-300
-4
10
10
-2
10
10
10
10
Frequency (rad/sec)
Figure 2-11 Open-loop Bode plots of the current controller transfer function (i) without PI
control (blue solid line); (ii) with PI compensation using initial parameters (green dashed
line).
33
without PI
with PI
1.2
0.8
0.6
0.4
0.2
0.01
0.02
0.03
time (sceonds)
0.04
0.05
Figure 2-12 Step responses of current controller transfer function (i) without PI control
(blue solid line); (ii) with PI compensation using initial parameters (green dashed line)
PI compensation is applied (Kp1=0.62, Ki1= 53) to reduce the gain cross-over frequency
thereby increasing the phase margin to 50.4 deg and the gain margin to 17.1dB (dashed
green line in Figure 2-11). The overshoot of 23.6% (dashed green line in Figure 2-12) is
lower than the original 30%. Although the settling time has been increased to 0.025s, it is
still acceptable.
b) Outer Constant DC Voltage Control
Substituting numerical values of the system parameters in equation (2-22) yields the
following compensated open-loop transfer-function for the DC voltage controller,
F2 ( s)
Ki 3 K p 3 s
3.33592 104 s3 0.1127 s 2
(2-32)
The frequency response of the plant transfer function Gp2(s) in equation (2-22) is
calculated and displayed in Bode plot as shown by the solid blue line in Figure 2-13. The
phase margin of 88.5 deg is significantly higher than the specified maximum of 60 deg,
and the settling time of 0.432 (solid blue line in Figure 2-14) is considered excessive. The
bandwidth of the controller is increased with the PI compensation (KP=20, Ki=1000). The
compensated open-loop transfer-function (dashed green line in Figure 2-13) has a phase
margin of 47.1 deg, decreased settling time 0.046s (dashed green line in Figure 2-14). The
overshoot of 30.2% (dashed green line in Figure 2-14) is relatively high, but it will be
revised in the next optimization procedure.
34
Magnitude(dB)
100
without PI
with PI
50
0
-50
-100
0
10
10
10
10
10
Phase(deg)
-50
-100
-150
-200
0
10
10
10
Frequency (rad/sec)
10
10
Figure 2-13 Open-loop Bode plots of the outer DC controller transfer-function (i) without
PI compensation (solid blue line); and (ii) with PI compensation using the initial
parameters (dashed green line).
Step Response Based on Simplified Linear Model
1.4
without PI
with PI
1.2
DC voltage (pu)
1
0.8
0.6
0.4
0.2
0
0.2
0.4
0.6
time(sceonds)
0.8
Figure 2-14 Step responses of the outer DC controller transfer-function (i) without PI
compensation (solid blue line) and (ii) with PI compensation using the initial parameters
(dashed green line).
c) Outer Active/Reactive Power Control
Similarly, inserting the numerical values into equation (2-25), the following
compensated open-loop transfer-function of the active/reactive power controller is
obtained.
Ki 4 K p 4 s
3
F3 ( s)
2 29.6 104 s 3 s 2
(2-33)
35
Magnitude(dB)
100
0
-100
-200
-1
10
10
10
10
10
10
Phase(deg)
-50
-100
-150
-200
-1
10
10
10
10
10
10
Frequency (rad/sec)
Figure 2-15 Open-loop Bode plots of the outer active/reactive power controller transferfunction (i) without PI compensation (blue solid line) and (ii) with PI compensation using
the initial parameters (dashed green line)
t eid _inner _ rec dt t eiq _inner _ rec dt t eid _inner _ inv dt t eiq _inner _ inv dt
(2-34)
Where ej(t) is the difference between the jth controlled variable and its reference value.
The value of j is set to unity indicating that all of the controlled variables have the same
importance.
The vector X= (Kp1 Ki1 Kp2 Ki2 Kp3 Ki3 Kp4 Ki4 Kp5 Ki5 Kp6 Ki6) represents the
variables.
36
pu
pu
Time (seconds)
Active Power at Rectifier end
Time (seconds)
(b) Reactive Power at Rectifier end
pu
pu
(a)
Time (seconds)
(c) DC voltage at inverter end
Time (seconds)
(d)Reactive Power at Inverter end
37
settling times and the relatively high overshoot of the other reactive power controllers and
DC voltage controller are considered unsatisfactory (see Figure 2-16c). Therefore the
definition of the objective function is refined according to the controller design
requirement as illustrated in Figure 2-16. As an example, an increased weighting of 4 is
applied to the DC voltage controller and the active power controller which results in the
following objective function as equation (2-35).
Table 2-1 PI parameters comparison
Optimized Values Initial Values of PI
Of(X)
of PI Parameters
Parameters
3.2939 (3s)
3.4139 (3s)
Refined Values of
PI Parameters
4.5148 (4s)
Kp
Ti(1/Ki)
4.4717 (4s)
Kp
Ti(1/Ki)
Kp
Ti(1/Ki)
Q(rec)
0.4809
3.6540
0.0436
0.0512
0.6200
3.6500
0.0189
0.0500
0.4850
3.6577
0.0484
0.0508
P(rec)
3.6467
0.0145
3.6500
0.0500
3.6500
0.0139
Q(inv)
3.6569
0.0436
3.6500
0.0500
3.6569
0.0427
Udc(inv)
20.0061
0.0032
20.0000
0.0010
20.0107
0.0019
Inner Current(inv)
0.4765
0.0312
0.6200
0.0189
0.4688
0.0346
t eid _inner _ rec dt t eiq _inner _ rec dt t eid _inner _ inv dt t eiq _inner _ inv dt
(2-35)
It can be seen from Figure 2-16 that the performance of the active and reactive power
controllers at the sending end are barely altered as a result of the refined optimization
function. However, the dynamic performance of the DC voltage controller and the
receiving end reactive power controller are improved significantly with refined objective
function. The values of the new objective function have been reduced from 4.51 to 4.4729
within duration of 4s after optimization. To conclude, the dynamic behaviour of the VSCHVDC system is improved with the refined set of PI parameters.
38
39
branches or all of them (which is the exact PR controller). Results showed that it is this
elimination that caused the increase of the frequency sensitiveness. This study offers a
systematic method to understand this type of control in detail, and explains the
discrepancies from the frequency analysis and experimental verification point of views.
A more complicated small signal model of PR controlled VSC was first proposed in
paper [79] and thesis [80], which provides researchers a valuable tool to understand how
such controller can affect the stability of an entire system.
One of the most preferable characteristics of PR controllers is it is easy to compensate
on the low-order harmonics. This is done simply by positioning the generalized integrators
that tuned to resonate at the harmonic frequencies to be compensated, in parallel with the
main controller without affecting the dynamics of the entire closed-loop system as shown
in Figure 2-17 [54, 120].
is_ref
Kp
sK i3
s +302
sK i5
s +502
sK i7
s +702
u*
is
GPR ( s) K p Ki
s
s 02
2
(2-36)
s
s 2 (h0 )2
(2-37)
This is because the PR controller works at a very narrow band around the resonant
frequency 0. However the bandwidth of these compensators has to be smaller, lower
enough than the bandwidth of the system to prevent triggering the system low frequency
40
( R sL) is us uc
u c
u s
(2-38)
is
1
R+sL
K p1
is
sK i1
s 2 +02
+
+
2 m
U dc
Figure 2-19 The simplified model of PR current controlled VSC system in reference
frame
41
Then the closed-loop transfer function of the system yields to equation (2-39),
GIc ( s)
I s
I s _ ref
K p1 ( s 2 02 ) sK s1
(2-39)
( R sL)( s 2 02 ) K p1 ( s 2 02 ) sK s1
The proportional gain of the controller Kp1 can be tuned by root locus theory provided
that the integrator gain K i1 is set to be zero first. Then the integrator of the controller Ks1
can be obtained by frequency response analysis of the open loop transfer function. As
suggested in [123], the controller parameters for PR scheme should be kept the same as the
parameters of PI controllers.
Smith Predictor
DB
Two samples delay DB with compensation
(Achieve current command in two sampling periods )
IMC
Solving Feedback TF
42
signal) is composed of two parts: the output of voltage prediction block Uv_abc_2 and the
output of the DB current control block Uv_abc_1.
Uf_abc
Usabc
-
Uv_abc
Rac
Lac
Ia
Rv
Lv
Rac
Lac
Ib
Rv
Lv
Rac
Lac
Ic
Rv
Lv
C
-
U dc
Triangular
Carrier
SPWM
Iabc
Uv_abc_ref
Id_ref
Iq_ref
Iabc_ref2
Dead- beat
Control
d-q-0 to a-b-c
Transformation
Uf_a
Uf_b
Uf_c
Uv_abc1
Uv_abc2
U line_
Uf_abc
Phase
Forward
U line_
a-b-c to d-q-0
Transformation
Voltage
Prediction
Phase Lock
Loop
43
interval [53, 54, 124]. The DB current control block in Figure 2-21 generates a precalculated additional voltage that is equal and opposite to the voltage reduction caused by
the Rv and Lv. This allows the current in the inductance is not affected by the PCC filter
bus voltage. The PLL is used for synchronizing the PCC voltage. A
to abc reference
frame transformation block is also needed in the control system. The phase forward loop
block in the control system is employed for time delay compensation.
2.3.3.2 DB Current Control with One Sample Time Delay
The closed-loop DB current control block without considering the computation time
delay is illustrated in Figure 2-22. The internal model control for design of DB control can
be expressed by using the following set of equations [51],
G 1 ( z ) z 1
GDB ( z )
(
Lv Rv ) / ( z 1)
z 1
Ts
(2-40)
Lv
L
( Rv v ) z 1 ) / (1 z 1 )
Ts
Ts
(2-41)
Lv
L
I (kTs ) ( Rv v )I ((k 1)Ts )
Ts
Ts
(2-42)
where GDB (z) is the transfer function of the current control, uva_1 is the first part of
VSC input voltage generated by DB current controller, and kTs represents the kth time
interval, uva is the total input voltage for VSC, ufa((k+1)Ts) is the predicted filter bus
voltage at (k+1)Ts time interval, and I represents the current change order. Hence the
control algorithm can be implemented by using equation (2-42).
Ia_ref
G DB (z)
GZOH(z)
G(z)
Ia
Ia
Figure 2-22 The closed-loop DB current Control
The closed-loop transfer-function of the one sample time delay DB current control with
the presence of the one sample computation time delay can be given by,
44
GDB _ CL ( z )
z 1 / ( z 1)
z 1
1
1 z 1 / ( z 1) z 1 z 1 z 2 z 1
(2-43)
Therefore, one sample delay DB current control does work at the stability margin with
computational delay).
A. Reducing the Proportional Gain Method
This method used in papers [54, 124] does move the poles inside the unit circle, but at
the cost of changing the closed-loop gain, which increases the phase difference between
the reference and output current. For example, by changing the proportional gain part Lv/Ts
of the controller to Lv/(2Ts) and Rv to Rv/2, the system characteristic equation becomes
z2 -z+1/2=0 . This results in two new poles as z1,2 = 1/2 (1j), which are inside the unity
Pole-Zero Map
1
0.6/T
0.8
0.5/T
0.10.3/T
0.7/T
0.6
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.8/T
Imaginary Axis
0.4
0.9/T
0.2
0
-0.2
0.4/T
0.2/T
0.1/T
1/T
1/T
0.9/T
0.1/T
-0.4
-0.6
-0.8
0.8/T
0.2/T
0.7/T
0.3/T
0.6/T
-1
-1
-0.5
0.5/T
0
Real Axis
0.4/T
0.5
Figure 2-23 Pole-zero map of the one sample delay system: red cross: one sample delay
DB without considering the computation delay time; blue cross: represents one sample
delay DB considering the computation delay time; green cross: represents reducing the
proportional gain;
45
Ia_ref
G DB (z)
VDB
z -1
Ia
G(z)
Ia
G(z)(1-z -1 )
G DB (z)
VDB
z -1
G(z)
Ia
Ia
G(z)(1-z -1 )
GDB ( z )
G ( z)
DB 1
1
1 GDB G ( z )(1 z ) 1 z
G 1 ( z )
( z 1)(1 z 1 )
(2-44)
46
2
s 2 2s
G '1 ( s)
(2-45)
where, represents the bandwidth of the low pass filter adopted for damping purpose.
After discretizing the equation (2-45), the following controller is obtained,
GDB 2 (z)=
G '1 (z)
1
Ts
(2-46)
G (z)
(1 z 1 ) (z 1)
Lv
L
I ((k 1)Ts ) ( Rv v )I (kTs ) uva _1 ((k 1)Ts )
Ts
Ts
(2-47)
The developed control algorithm is very similar to the one used in [54], which clearly
demonstrate the derivation of the control algorithm from a different point of view (IMC).
The algorithm can be implemented in PSCAD according to Figure 2-25.
Sample&Hold
u va_ 2
Ia
Ia_ref
Sample&Hold
K1
+
++
-K 2
u va_1
u va_ref
Z-1
Gp(z)
Ia
e-sTs
e-2sTs
47
Ia_ref
u va_1
Ia_ref (k)
Execution
-1
delay z
G db (z)
Ia
Gp(z)
GZOH(z)
T
GF(Z)
Figure 2-26 The block diagram of solving feedback transfer function DB current control
In the above figure, the block Gdb(z) and the feedback loop transfer function GF(z) are
given by,
Gdb ( z) N1 / (1 N2 z 1 )
(2-48)
GF ( z ) N 2 2
(2-49)
where, N1 Rv / (1 N2 ) , N2 e( Rv / Lv )Ts
The control algorithm can be implemented in PSCAD using the following equation
(2-50) and is also illustrated in Figure 2-27.
(2-50)
G F (s)
Ia_ref (t)
Sample&Hold
u va_ 2
Ia
u va_1
N1
e-sTs
++
u va_ref
Ia
z-1
Gp(z)
N2
Figure 2-27 The structure of the solving feedback transfer function DB current controller
48
49
As they are shown and observed in Figure 2-28 that the Smith Predictor method, the
IMC method and the solving feedback TF method have nearly the same bandwidth, the
same cut-off frequency and the same settling time. The performance of the reducing gain
method is observed slightly poor specifically with a relatively narrower bandwidth, higher
overshoot and longer settling time. Note that, the results in Figure 2-28 are obatained based
on the simplified linear models under the assumption that the voltage applied across the
AC side L/R circuit is unlimited. It is done for analytical purpose only, and which is not
the authentic case in practice.
2.3.3.5 Sensitivity to Plant Parameters
The expressions of the sensitivity to plant parameters for each of the four controllers
are shown in the Table 2-3 which are derived according to the sensitivity formula defined
as the percentage change in overall transfer function divided by percentage change in the
plant transfer function [12], The sensitivity formula is depicted by,
S (s)
G ( s ) / Go ( s)
G ( s ) GPo ( s)
GP ( s ) / GPo ( s) GP ( s) Go ( s)
GP ( s )
Go ( s )
(2-51)
1
1 GF ( s)Gdb ( s) G p ( s)
(2-52)
50
Sensitivity Expressions
Sensitivity Results
Reduce Gain
1 z 1
2 2 z 1 z 2
0.2239
1 z 2
0.6173
Solving Feedback
1 e2 aTs z 2
0.6152
75 MW
AC bus voltage
62.5kV
Inductance, Lv
0.15 pu
Resistance, Rv
0.0015 pu
0.001s
(pu)
(pu)
Time (seconds)
51
Time (seconds)
(pu)
(pu)
a)
Time (seconds)
Time (seconds)
52
(pu)
(pu)
Time (seconds)
Time (seconds)
(pu)
(pu)
c)
Time (seconds)
53
Time (seconds)
Figure 2-29 PSCAD step response simulation results of the DB current controller
It can be observed from Figure 2-29a that the reducing gain method has a frequently
altered current as seen from dq synchronous reference frame plots due to the existence of
phase difference caused by reducing the gain. The settling time of this controller is fast, but
its overshoot in reactive current is also the highest among these four controllers, which is
not desirable. The Smith Predictor method has approximately the same performance with
IMC based DB when being subject to the step response tests. However, the overshoot in
reactive current is observed slightly higher than IMC based DB controller. Compared with
the IMC based method, the solving feedback TF method has a much shorter settling time,
but with a slightly higher reactive current overshoot. The summary of the characteristics of
the simulation results are given in Table 2-5.
Table 2-5 Summary of specific data obtained from Figure 2-29
Overshoot
Controller
Step Test
Settling
Overshoot
Reactive
Type
Type
Time (ms)
(%)
Current (%)
Error
(pu)
Large Step
14.9
16
84.6
0.479
Small Step
2.5
0.2
Smith
Large Step
26
25
48.8
0.349
predictor
Small Step
20
5.46
0.110
IMC
Large Step
25
17.2
34
0.358
Designed
Small Step
21
4.9
0.067
Solving
Large Step
14.5
20
44.78
0.3348
Feedback
Small Step
14.9
16
84.6
0.479
Reduce Gain
54
Largest
Deviation
( I d I d 0 )2 ( I q I q 0 )2
I d 02 I q 02
(2-53)
where, Id0 and Iq0 represent the resulted current under original plant parameter.
The simulation results shown in the Table 2-6 demonstrate that the resistance has little
effect on the sensitivity, which is mainly determined by the inductance. The sensitivity of
solving the feedback transfer function method is somewhat lower than the IMC based DB
controller in this particular case, and it is the same with the analysis given in section
2.3.3.4. However, the simulation result is significantly smaller than the results given in
section 2.3.3.4, which is likely due to the inaccuracy introduced during simplification. The
simulation results show that the output of control system have an approximate 3%
variation when with the 10% change in an inductance value.
55
Original Plant
Iq
Id
Iq
Id
Iq
Id
-1
-1
-1
-1
R L
0.8382
0.2557
0.5
0.203
0.354
0.358
R L
0.559
0.487
0.65
0.254
0.703
R L
2.761
2.716
2.75
3.198
3.147
2.974
R L
3.25
2.795
3.238
3.758
3.522
R L
3.509
3.27
2.761
3.198
3.041
2.974
R L
3.25
2.795
3.5
3.238
3.487
3.48
2.4 Conclusion
This chapter has investigated several linear inner current controllers from operational
principles introduction and their controller design methodologies. Firstly, the operation
principle of the PI control is introduced and discussed. Furthermore, a methodology for the
design and optimization of the various loops of controllers of the PI controlled VSCHVDC transmission system is developed. This is adequately effective for a relatively
strong system. However, for the weak AC system application, the inaccurate feed-forward
compensation caused by low-order harmonics distortion in the case of PI control structures
stands for the major disadvantage [61, 120]. In addition, the cross-coupling terms also adds
to these shortcomings further. The theory and the simple control design method for the PR
control strategy are also introduced in this chapter. It was reported that this method has a
serious frequency sensitiveness problem, which makes it hard to implement when
interconnected with a weak AC grid characterized by large grid impedance.
As suggested in [54], the DB controller is the best-tailored scheme in a weak AC
application in terms of robustness against parameter variation and system disturbances.
Therefore, this chapter also studied four different types of DB controllers by including one
sample delay with the i) reducing gain and ii) Smith Predictor, two sample delays based on
iii) IMC control design and iv) solving feedback transfer function. It was also found that
56
2.4. CONCLUSION
the two sample delays DB current controller based on IMC derived and introduced in this
chapter demonstrates similar performance as the one employed by A. Timbus, M. Liserre,
R. Teodorescu, P. Rodriguez, and F. Blaabjerg in [54].
The analysis of the results showed similar step responses for both of two sample delays
DB current controllers. However, simulation results indicated that solving feedback
transfer function method proves to be superior to the other methods specifically during the
large step response with a fast settling time and with an acceptable reactive current
overshoot. In addition, it was observed that simulation results of the same method presents
smaller variations under parameter sensitivity tests compared with the IMC based DB
current controller. Therefore, the solving feedback TF based DB control method is chosen
as the definite inner current control scheme which is utilized in the remaining chapters of
this thesis.
57
58
synchronous reference frame was proposed as an extension to the Charles Saos and Peter
W Lehns works in [78-80].
This newly proposed reference frame transformation approach, which facilitates the
understanding the DB controlled VSC, is one of the original contributions of this thesis.
The study aims at developing a detailed higher order small signal model for the discrete
DB current controlled VSC, which can facilitate further linear controller design, system
parameterizations when implementing such type of VSC based applications, and
establishing the foundation for the system small-signal stability analysis that Modal
analysis (e.g. Eigen values, mode shapes, participation factor, etc.).
and
are pre-set based on the specified active power Q and reactive power P, where Ic stands for
converter output current, the superscript C represents the converter frame, the subscript dq
implies the variable is in the dq synchronous reference frame, the subscript RI denotes
Real-Imaginary, and the subscript ref indicates it is a reference value. Note that the
similar format is used in the rest of the thesis. Then the input current references are
compensated with a two sample period ahead values, followed by a sampling block to
obtain the discrete form of the input current references
and
. In the
meantime, the measured converter output currents pass through a ZOH block and feedback
gain GF(s), which arrive at
and
respectively.
Based on the reference values and measured values of the input currents, the current
errors
and
these current errors, the voltage drop across the controlled plant can be predicted by the
DB control algorithm, including
and
voltage should also be accessed and compensated by a 1.5 sample period ahead values.
Then the second part of the converter input voltage references
and
,
and
respectively.
59
PLL
c
cq _ q
I (k )
c
cd
(k )
(k )
Vcdc _ ref 1 (k )
I cdc _ e (k )
Dead- beat
Control
+
I cqc _ e (k )
Vcqc _ ref 1 (k )
I (k )
GF ( s )
+
c
cq
Vcdc _ ref (k )
Vcqc _ ref (k )
c
cd
d
Computation
Vcqc _ ref (k )'
q
Delay
VcR _ ref
R
I
VSC
Sampling
V
c
cq _ ref 2
Vcdc _ ref 2
Vcqc _ ref 2
(k )
I
I cqc _ ref (k )' Sampling
U dc E g
PLL
c
'
cd _ ref
I cIg
u gfR
u gfI
EIgg
Rg
c
fd
u cfq
d
q
R
I
PLL
u cfd '
cos(1.50Ts ) sin(1.50Ts )
sin(1.5 T ) cos(1.5 T )
c '
0 s
0 s u fq
Voltage Prediction
Grid
g
cI
Filter
Vcdc _ ref 2 (k )
I cRg
VcRg
VcI _ ref
PLL
R
I
d
q
I cqc
Sampling
Filter
+
+
c
cd _ d
cos(20Ts ) sin(20Ts )
sin(2 T ) cos(2 T )
0 s
0 s
I cdc _ ref
I cqc _ ref
Figure 3-1 Small signal model of VSC with the discrete DB current controller
The computation delay needs to be taken into account as well which is approximately
one sample period (Z-1). Before conducting the calculation in grid reference frame, the
VSC input references (
be transformed into the grid reference frame, dq to RI block, as shown in Figure 3-1.
Similarly, the grid output current and voltage should also be back transformed from the
grid frame to the converter frame, RI to dq block, before conducting the converter frame
calculation.
By ignoring the dynamics of the DC voltage, the VSC is supposed to be ideal. This
means that the VSC ideally outputs voltage as required. Utilizing this voltage, the
converter reactance and the filter bus voltage, the converter output current can be easily
obtained.
60
The procedure on how to develop the small signal model for each individual
component needs to be carefully addressed. Therefore, the block diagram and its
corresponding non-linear equations are established first, which is followed by a
linearization process to obtain the small signal model and the small signal state-space
equations. Finally, simulation results based upon tests setting on both of Simulink/Matlab
and PSCAD platforms are utilized to assess the accuracy of the linearized models. The
following sub-sections of this thesis will provide the details of these studies. It should be
clarified that in the remnent of the thesis the large signal model can also be termed as nonlinear model which is detailed model transient simulation (PSCAD). Similarly, the small
signal model can also be named as linear model which is simplified s-domain simulation
(Matlab).
61
blocks sin () and cos () of the VCO block, which are magnified by a gain KL, then output
uv, uv. Finally the close loop PLL is formed.
u v
u s
u s
1
u sb
1
u sb
K p_PLL +
K i_PLL
uf
Kv
1
s
Central frequency 0
u v
u sb =voltage base
cos
KL
Phase Detector
Low Pass
Filter
KL
sin
Voltage Controlled
Oscillator
d us uv us uv
usp sin(1t 1 ) K L cos usp cos(1t 1 ) K L sin
(3-1)
Since the PLL is assumed to be locked to its input signal, then the output of PLL can
be given by,
2t 2
(3-2)
(3-3)
d usp K L sin[et (1 2 )]
d usp K L (sin et cos(1 2 ) cos et sin(1 2 ))
(3-4)
Note that if the et is small enough, then sin e t e t and cos e t 1 . However, this
is not strictly true since the grid frequency is a time dependent term. Therefore,
it is worth noting that under some condition, small variation in grid frequency may results
in large phase difference. Hence,
62
(3-5)
Assuming the variations in grid frequency and initial phase angle is small enough,
which means 1=2 and 1= 2,
d usp K L sin(1 2 )
(3-6)
d usp K L (1 2 )
Therefore, the double frequency components 21t are fully eliminated by this
arrangement. However, under unbalanced conditions, phase can be locked to the positive
sequence component of its input signal.
3.2.1.1 Derivation of Linearized Model of PLL
Linearized model of the PLL will be derived in this sub-section considering dynamic of
each individual component of PLL independently.
For facilitating the analysis, the nominal system frequency 0 of the VCO which is
assumed to be constant is involved linking the input phase angle 1 and the output phase 2.
Using Figure 3-2, the input reference frame voltages to PLL block can be given by
us
sin(0t 1 )
u usp
cos(0t 1 )
s
(3-7)
sin(0t 2 )
v
and
(3-8)
can be given by
1 1 (t )t 1 (t ) 0t
(3-9)
2 (t ) 0t
(3-10)
A. Phase Detector
Note that using equation (3-9), one can derive that 1 (t)=1 +0 t-1 (t)t . If 2 =1 ,
using equation (3-2) and (3-10), 2 (t)=2 +0 t-1t can be obtained. If 1(t) and 2(t) are
substituted into equation (3-6), the controlled error can be obtained as
d usp K L sin(1 2 )
d usp K L sin(0t 1 1t (0t 2 1t ))
d usp K L sin(1 2 )
(3-11)
63
d usp K L (1 2 )
(3-12)
d / (1 2 ) usp K L Kd
(3-13)
u f ( s)
d ( s)
K p _ PLL
Ki _ PLL
s
(3-14)
Ki _ PLL
s
)d
(3-15)
(t ) [0 ]d 0t KV _ PLLu f ( )d
t
2 (t ) (t ) 0t 2 (t ) KV _ PLLu f ( )d
2 ( s)
u f ( s)
(3-16)
(3-17)
K v _ PLL
s
(3-18)
K v _ PLL
s
u f
(3-19)
d us cos us sin
2
3
1
3
1
3
3
3
usb
usc ) sin
3
3
usa
2
2
2
d cos( ) cos( ) cos( ) usb
3
3
3
usc
64
(3-20)
d usd
This means
d usd usp K L (1 2 )
Therefore,
(3-21)
d usd usp K L (1 2 )
(3-22)
u sd
K p_PLL
K i_PLL
KV
1/s
usd l
dt 2 KV K p _ PLL
KV 0
and,
0 KV K p _ PLLusd Kvl
(3-23)
(3-24)
When the above state space equations are linearized, the small signal model for PLL
can be obtained as,
0
d l Ki _ PLL
usd l
dt 2 KV K p _ PLL
KV
(3-25)
2 (s)/1 (s) , where 2 (s) denotes the output of PLL and 1 (s) represents the input to
the PLL control block. The block diagram of the linearized model of the PLL is illustrated
below,
65
1 (s)
K d =u sp K L
d (s)=u sd
u f (s)
PI
Kv
s
2 (s)
2n s n 2
2 ( s)
2
1 ( s) closeloop s 2n s n 2
K
n K d KV Ki _ PLL , n p _ PLL
2 Ki _ PLL
(3-26)
(3-27)
(3-28)
L 2n Kd KV K p _ PLL
(3-29)
n L / 2
(3-30)
(3-31)
where, SMargin is the stability margin, max and min are the maximum and minimum
frequency that the PLL can lock to. In practice, the VCO control signal is usually limited to
a range smaller than the VCO supply voltage, which is mostly +5kV [52]. For example, the
66
values of ufmax and ufmin, with +5kV supply voltage, can be chosen as 90% 5kV and
10% 5kV respectively.
C. Let Kd = 1, then Determine Kp_PLL and Finally Determine Ki_PLL
It should be noted that the lock-in range and speed response to step change of PLL is
proportional to the characteristic frequency, but inversely proportional to the capability of
harmonics rejection. Decision on the value of n has to be made by trading off the
response speed and noise rejection capability of PLL according to the requirements of
applications.
3.2.1.4 Case study
Using the PLL parameters shown in Table 3-1, the transfer function of PLL can be
defined as equation (3-32).
K d K v ( Ki _ PLL K p _ PLL s )
2 (s)
2
1 ( s) closeloop s K d K v K p _ PLL s K d K v Ki _ PLL
(3-32)
Us_phase_peak
51.031
kV
KL
Kd
51.031
KV
90
PI Parameter
Nature
Kp_PLL
0.163
Ki_PLL
31.25
213.74
Damping Ratio
1.751
Lock in Range
Frequency
748.61
Figure 3-5 shows the frequency response of the PLL, where the gain margin is infinite
and the phase margin is 84.5o. This indicates that the designed PLL is stable.
67
80
Magnitude (dB)
60
40
20
0
Phase (deg)
-20
-90
-135
Pm=85.4 deg (at 2.336e+003 rad/sec)
-180
1
10
10
10
Frequency (rad/s)
10
Rad
0.04
0.02
0
-0.02
0
0.005
0.01
time(s)
PSCAD
MATLAB
0.015
0.02
68
Rise
Time
0.00094
0.9553
1.059
5.925
1.059
0.0025
Rad
2.5
2
1.5
1
0.5
0
0.01
0.02
0.03
PSCAD
Matlab
0.04
t(s)
69
Time (seconds)
Time (seconds)
Figure 3-9 PLL responses to 90% magnitude step change of input voltage simulated in
PSCAD
Time (seconds)
Figure 3-10 PLL output phase angle in comparison with input phase angle when being
subject to a sudden 90% magnitude step change of input voltage simulated in PSCAD
70
Time (seconds)
Figure 3-11 PLL output frequency in comparison with input frequency when being subject
to a sudden 90% magnitude step change of input voltage simulated in PSCAD
The results given in Figure 3-9 to Figure 3-11 show that the sudden magnitude
reduction has no effect on the PLL outputs.
E. Single Phase to Ground Fault
In this simulation study, while the initial phase is zero at a constant frequency of 50Hz
and at a constant magnitude, a ground fault is created at A-phase at t=1s.
Figure 3-12 shows the phase response of PLL under a phase to ground fault. As is
shown in the figure, the result of twice the fundamental frequency, which can be expressed
below,
y 0.00346cos(200 t )
(3-33)
Time (seconds)
Figure 3-12 PLL output phase angle when being subject to an A-phase to ground fault
simulated in PSCAD
Figure 3-13 shows the frequency response of PLL under the same fault, which can be
expressed by,
71
y 50.3475cos(200 t )
(3-34)
Time (seconds)
Figure 3-13 PLL output frequency when being subject to an A-phase to ground fault
simulated in PSCAD
3.2.1.6 Summary of PLL
The developed linearized model and non-linear model of PLL both have the fast
response and wide frequency range acquisition abilities. In addition, the linear model and
non-linear model match each other very well. Its performance for a large reduction in
voltage magnitude and for a single phase to ground fault are also examined in this section.
u cfdq e
72
j ( _ PLL /2)
c
'
Vcdq
_ ref (k )
j ( _ PLL /2)
j ( _ PLL /2)
(3-35)
g
I cRI
(3-36)
u gfRI
(3-37)
qc
0 t
_ PLL
Rg
dc
0 t
s t
Ig
U f
Figure 3-14 Relationship between converter reference frame and grid RI reference frame
(3-38)
VcIg _ ref (k )' cos PLL Vcdc _ ref (k )' sin PLL Vcqc _ ref (k )'
The phase angle PLL in the above set of equation is the phase difference between the
filter bus voltage and grid phase reference point, which is a constant value by excluding the
time varying component 0t.
The linearized form of equation (3-38) can be given as
VcRg _ ref (k )' sin PLL 0 Vcdc _ ref (k )' cos PLL 0 Vcqc _ ref (k )' VcIg _ ref (k )'0 PLL
VcIg _ ref (k )' cos PLL 0 Vcdc _ ref (k )' sin PLL 0 Vcqc _ ref (k )' VcRg _ ref (k )'0 PLL
(3-39)
; ii) measuring the phase of the filter bus voltage Us, and the phase
difference PLL between the grid source voltage and the filter bus voltage. Note that using
the angle , the representation of Us in the converter reference frame (
) can be
73
computed, which is then followed by a transformation to grid frame using equation (3-35).
This calculation utilizes PLL to obtain
U G G
LG U s
Time (seconds)
varied at t=0.5s on q-axis of the filter bus voltage. This is followed by a secondary 0.99 pu
C
large signal u sq0 that varied at t=0.5s imposed in another simulation before transforming to
the grid RI reference frame. The errors in the output quantities of both simulations are
G
G
given by u sR_Large
and u sI_Large
. In addition, a 0.01pu small step change on q-axis of the
filter bus voltage is applied simultaneously to the small signal model at t=0.5s. Then the
G
G
outputs u sR_small
and u sI_small
are compared with the results obtained from the non-linear
74
0pu
1.0pu
0pu
usdC 0
u
usqC 0
usqC 0
usdC
+0.99pu
G
dq to RI
usR
0
sin PLL cos PLL
G
usI 0
cos
sin PLL
PLL
C
sd 0
usqC +-
usdC
sin PLL
cos
PLL
usqC
usdC
sin PLL
cos
PLL
usqC
Non-Linear
g
usR
_ L arg e
-+
-
usIg _ L arg e
G
usR
cos PLL
G
sin PLL usI
g
usR
_ small
usIG 0
+-
cos PLL
sin PLL
++
G
usR
0
PLL
PLL
g
usR
_ small
Linear
Figure 3-17 The test circuit for small signal mathematical equation
The simulation results are given in Figure 3-18. Note that the exponentially slow rising
trend in the figure is actually caused by the combination of the voltage ramp up time which
is set as 0.02s in this test and the dynamics introduced by PLL utilized in the
transformation process of voltage from dq frame to RI frame. These results successfully
confirm the accuracy of the developed small signal model for the grid frame
transformation.
Time (seconds)
Figure 3-18 Small signal test results for the grid frame transformation
The large signal and the small signal algebraic equations for converter terminal current
Ic transforming from the grid RI reference frame to the converter dq reference frame are
given in equations (3-40) and (3-41).
(3-40)
75
g
I cdc sin PLL 0 I cR
cos PLL 0 I cIg I cqc 0 PLL
g
I cqc cos PLL 0 I cR
sin PLL 0 I cIg I cdc 0 PLL
(3-41)
Note that this transformation is similar to the VSC input voltage reference, which is
simply an opposite transforming direction. Therefore, no verification studies are needed.
The large signal and the small signal models for the filter bus voltage u Cf transforming
from grid RI reference frame to converter dq reference frame can be given by
u Cfd sin PLL u gfR cos PLL u gfI
u Cfq cos PLL u gfR sin PLL u gfI
u Cfd sin PLL 0 u gfR cos PLL 0 u gfI u Cfq 0 PLL
u Cfq cos PLL 0 u gfR sin PLL 0 u gfI u Cfd 0 PLL
(3-42)
(3-43)
u cfd '
u
c
fd
u cfq '
u
c
fq
1
1 Tf s
(3-44)
u fd
1
u cfd '
dt
Tf
Tf
du cfq '
dt
(3-45)
1 c ' u fq
u fq
Tf
Tf
(3-46)
After linearizing, the small-signal state-space model for the grid voltage signal filter
can be given by,
d u cfd '
dt
d u cfq '
dt
u fd
1
u cfd '
Tf
Tf
c
u fq
1
u cfq '
Tf
Tf
(3-47)
(3-48)
Note that there is no need to compare the large signal model and the small signal model
with respect to this component, since they do have the same formula. The only difference
is the step size to the input signals, which will naturally output proportionally results.
76
c
I c ref sin( _ PLL 0 20Ts ) cos( _ PLL 0 20Ts ) I cqref
(3-49)
c
'
I cdref
cos _ PLL 0 sin _ PLL 0 I cc ref
c '
c
I cqref sin _ PLL 0 cos _ PLL 0 I c ref
(3-50)
(3-51)
The below set of equation denotes the linearized form of the above equations,
c
'
c
c
I cdref
cos(20Ts ) I cdref
sin(20Ts ) I cqref
c
'
c
c
I cqref
sin(20Ts ) I cdref
cos(20Ts ) I cqref
(3-52)
To verify these equations, a 0.01pu (1pu-0.99pu) step change is applied to the inputs of
both the large signal model and the small signal model, and the tests results are given in
Figure 3-19, which shows a good agreement.
-3
10
x 10
8
6
Icref(pu)
4
2
0
Icdref-matlab
-2
Icqref-matlab
-4
Icdref-PSCAD
-6
0
Icqref-PSCAD
0.005
0.01
Time(seconds)
0.015
0.02
77
For the voltage prediction block in Figure 3-1, the relationship between the input
voltage u Cfd ' , u Cfq ' and the output predicted voltages in reference frame vCc_ref2 , vCc_ref2 can
be given by,
vcc _ ref 2 cos( PLL 0 1.50Ts ) sin( PLL 0 1.50Ts ) u cfd '
c
c '
vc _ ref 2 sin( PLL 0 1.50Ts ) cos( PLL 0 1.50Ts ) u fq
(3-53)
Similarly, the back transformation from the reference frame to the converter dq
reference frame for the predicted voltages can be written as,
c
'
vcq _ ref 2 sin PLL 0 cos PLL 0 vc _ ref 2
(3-54)
Substituting equation (3-54) into equation (3-53), the direct link between the input
voltage and output predicted voltage are obtained as equation (3-55),
vcdc _ ref 2 cos(1.50Ts ) u cfd ' sin(1.50Ts ) u cfq '
vcqc _ ref 2 sin(1.50Ts ) u cfd ' cos(1.50Ts ) u cfq '
(3-55)
(3-56)
The verification process for the accuracy of the linearization of this component is
similar to the current compensation block. Hence, it will not be repeated.
(3-57)
I cdc d (k ) GF I cdc (k )
I cqc d (k ) GF I cqc (k )
78
(3-58)
I cdc d (k ) GF I cdc (k )
I cqc d (k ) GF I cqc (k )
(3-59)
Note that, there is also no need to verify the linearization of this feedback loop since it
is only a coefficient gain.
The large signal model and the small signal model for generation of the DB current
input errors can be given by
c
I cdc _ e (k ) I cdref
(k )' I cdc d (k )
c
I cqc _ e (k ) I cqref
(k )' I cqc d (k )
c
I cdc _ e (k ) I cdref
(k )' I cdc d (k )
c
I cqc _ e (k ) I cqref
(k )' I cqc d (k )
(3-60)
(3-61)
79
F4 (s)=(1-sTs /6)3 /(1+sTs /6)3 and iii) a comparably crude approximation given by a first
order lag F5 (s)=1/(1+sTs ) are investigated. The frequency responses of the delay are
compared with that of the above three approximations. Furthermore, the time-responses of
the exact delay and the approximations are compared for (i) a step signal; and (ii) a
sinusoidal input signal.
a) Frequency Response
Figure 3-20 shows the frequency responses of the pure one sample delay in s- and zdomain and the three s-domain approximations. It can be seen from the figure, the z -1 and
e-st are overlaid with each other as expected. The third order Pad approximation is the
most accurate equivalence among the three approximations. However, for the powersystem modeling purposes outlined in section 3.2.6.1, the first order PadApproximation
is more than adequate (where the sampling time Ts equals to 1ms) since it matches the
exact delay at least up to 389rad/s (62Hz), and our interested frequency range is within this
limit.
Magnitude(dB)
10
0
-10
-20
-30
0.1/Ts
e-sTs
First-Order Pade
Third-Order Pade
First Order Lag
Z-1
1/Ts
10/Ts
1/Ts
Frequency (rad/sec)
10/Ts
Phase(deg)
200
100
0
-100
-200
0.1/Ts
80
b) Step Response
The unit step responses of the one sample delay and its approximations are shown
below in Figure 3-21.
1.5
1
pu
0.5
Z-1
0
-0.5
-1
0
Ts
2Ts
e-sTs
First-Order Pade
Third-Order Pade
First Order Lag
3Ts
4Ts
5Ts
Time in units of Ts
Over
Under
Peak
Peak
Time
Time
Min
Max
Shoot
Shoot
e-st
8Ts
9.8Ts
100 Ts
z-1
Ts
Ts
1.1Ts
1.3Ts
0.8026
1.0138
1.3755
100
1.0138
1.7Ts
1.1Ts
2 Ts
0.8057
100
2.2Ts
3.9 Ts
0.9029
10.5 Ts
Third order
Pad
First order
Pad
First order lag
Time
Table 3-3 presents the characteristics of the step response. As it can be observed in
Figure 3-21, the one sample delay e-st and Z-1 performed very well as expected (exactly as
one sample delay). Among the three approximation methods, the first-order lag has the
longest settling time and rise time. The third order PadApproximation and the first-order
PadApproximation have the same rise time which is 1.1Ts. It was found that the third
order Pad Approximation has a comparably shorter settling time which is preferred.
However this approximation also has a higher overshoot which is undesirable. Note that,
both approximation methods have a common defect that imposes a positive zero (RHP
81
zero) to the system. Reflected on the unit step response performance, the plots start arising
from -1. Taking all the factors into consideration, the first-order Pad Approximation is
found to be the most suitable and sufficiently accurate to approximate the pure one sample
delay in this study.
B. Zero Order Holding Block
a) Transfer Function of Zero Order Holding Block in S Domain
The mathematical representation of zero order holding block is
(3-62)
If we input an ideal unit impulse signal to the ZOH block, a rectangular impulse with
magnitude 1 and the lasting time Ts will be obtained as shown in Figure 3-22, which can be
divided into a unit step followed by a delayed negative unit step with a time constant T s.
Note that, the ZOH output signal has two typical characteristics in Figure 3-22: i) contain
high frequency component (single frequency input) and ii) has the fundamental frequency
shifted by Ts/2.
Time (seconds)
82
r (t )
Ts
D( z )
r[k ] r[k Ts ]
ZOH
y[k ]
y (t )
(3-63)
0, t 0
u (t ) 0.5, t 0
1.0, t 0
(3-64)
y
y(t)
y[k]
0
(k 1)Ts
kTs
-y[k]
1 skTs
s k 1 T
e
e s y k
s
1 esTs
y k e skTs
s
sTs
1 e
y(Ts )
s
y(t)
(3-65)
GZOH ( s)
1 e sTs
s
(3-66)
83
Ts
1 e sTs 1 (1 sTs / 2) / (1 sTs / 2)
( s)
s
s
1 sTs / 2
(3-67)
Magnitude(dB)
-50
-100
0.1/Ts
1/Ts
10/Ts
1/Ts
10/Ts
Phase(/deg)
0
-50
-100
-150
-200
0.1/Ts
Exact
First-Order Pade
Frequency (rad/sec)
84
1.4
1.2
1
pu
0.8
0.6
0.4
0.2
0
0
Ts
Matlab Exact
Matlab Pade
4Ts
5Ts
2Ts
3Ts
Time in units of Ts
0.4605
0.5641
262.5
0.4605
0.5641
0.5756
Verification by Matlab
According to the Bode plot of ZOH, the magnitude passing through the exact transfer
function should be 99.6% of the input magnitude and 98.8% for the approximate function.
The phase shift is 9 degree which is equal to half sampling time. Figure 3-27 shows the
Matlab simulation results, in which the amplitudes for the exact transfer function and the
first-order PadApproximation are 98.8% and 99.6% respectively. These results match the
desirable outputs in the frequency analysis.
85
Exact
First-Order Pade
pu
0.5
-0.5
-1
0.95
0.955
0.96
0.965
0.97
0.975
0.98
Time (s) (seconds)
0.985
0.99
0.995
Figure 3-27 Verification of responses to the sinusoid input for ZOH by Matlab
ii. Verification by PSCAD
PSCAD simulation studies are presented in Figure 3-28. The magnitude of the signal
passing through sampling block built in PSCAD is 1pu, which is ideal. The peak amplitude
of the output from the exact ZOH is 99.3% of the input signal and the output from the first
order PadApproximation is 98.8% of the input signal. The output amplitudes also match
(pu)
well with the attenuation value observed in the frequency response analysis.
Time (seconds)
86
Matlab Pade
Matlab Exact
PSCAD Sin Input
PSCAD Pade
PSCAD Exact
PSCAD Sample
0.8
0.6
0.4
pu
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
0.95
0.955
0.96
0.965
0.97
0.975
0.98
Time (s) (seconds)
0.985
0.99
0.995
-0.8
pu
-0.85
-0.9
-0.95
-1
0.973
0.974
0.975
0.976
Time (s) (seconds)
0.977
Matlab Pade
Matlab Exact
PSCAD Sin Input
PSCAD Pade
PSCAD Exact
PSCAD Sample
0.978
87
(3-68)
( a b s ) / (c d s )
a= N1 ; b= N1
where,
T
Ts
; c= (1 N 2 ) ; d= s (1 N 2 )
2
2
Magnitude(dB)
here.
80
60
Exact
First-Order Pade
First Order Lag
40
20
0.1/Ts
1/Ts
10/Ts
1/Ts
Time (s)
10/Ts
Phase(deg)
200
-200
0.1/Ts
Figure 3-30 Frequency response of inner current controller and its approximations
As it can be observed in Figure 3-30, the first-order lag approximation performs poorly
compared to the first-order Pad Approximation. Therefore, the first-order Pad
Approximation can still be considered as the most appropriate method. The attenuation
and phase shift for both the inner current DB controller and its first-order Pad
Approximation are found to be identical at 50 Hz.
88
Gdb ( z )
N1
N
z
z
A
B
1z
N1 z (
)
1
1 N 2 z z 1 N 2 ( z 1/ N 2 )( z 1)
z N2 z 1
A
where,
z
( z 1)
Gdb ( z ) N1 A
Hence,
B
z N 2
z
z N2
(3-69)
(3-70)
z 1
1
1
N1B
1
1 N2 z
1 z 1
(3-71)
Therefore, the unit step response of the z domain controller can be easily obtained in
time domain by using inverse of equation (3-71).
In this study, by substituting the parameter listed in Appendix A, the coefficient N1 and
N2 of the DB current controller are able to be obtained which are -24.904 and 0.9997
respectively. A and B are both 1/2. Therefore, the step response of the z domain DB
current controller is,
y (t ) 12.4695 12.4695(0.9997)
t
Ts
(3-72)
Ts 0.001s
0
-20
pu
-40
-60
-80
-100
-120
0
0.02
0.04
0.06
Figure 3-31 Step response of inner DB current controller and its approximation
89
Hence, the step response of the inner current controller should be a pulse train (high
frequencies oscillation) with an exponentially reduced magnitude to which the time
constant is long, as shown in Figure 3-31.
d) Response to the Sinusoid Input for Inner Current Controller together with ZOH
The aim of this section is to compare the actual discrete system with the continuous
Pad Approximation model and also to examine the feasibility of replacing the discrete
inner current DB controller with the continuous approximation to be utilized in power
system simulation. Figure 3-32 shows the continuous and discrete implementation schemes
of the inner current control strategy.
Input
GZOH _ Pade ( s )
GDB ( s )
GZOH ( z )
Sampler
GDB ( z )
N1
N2
est
90
10
pu
-5
-10
-15
-20
0.95
0.955
0.96
0.965
0.97
0.975
0.98
0.985
PSCAD Discrete
Matlab Pade Approximation
Matlab Discrete
0.99
0.995
1
Figure 3-33: Comparison of the responses of (i) the discrete DB current controller and (ii)
its first order PadApproximation to a sinusoidal input signal.
3.2.6.2 Transformation of DB Controller from Reference Frame to dq Reference
Frame
In the previous discussion, the discrete controller has been successfully transformed to
continuous domain. However, the DB current controller considered was defined in abc
natural reference frame as shown in Figure 3-34, which contains time-varying terms or
sinusoidal components which are inappropriate for developing the small signal model.
Therefore, it is necessary to transform the controller to dq reference frame. This will be
done in four steps as will be explained below.
C
c _ abc _ ref
C
c _ abc
C
c _ abc _ ref
( s)
I cC_ abc ( s )
S
+_
a bs
c ds
e sTs
GF
Figure 3-34 The abc natural reference frame DB current control block
91
Step 1: The model is simplified by reducing all the sampling blocks (circled area in Figure
3-34), moving to behind the computation block, and leaving the rest of the model remained
unchanged, as illustrated Figure 3-35.
a bs
c ds
+_
I cC_ abc
C
VcC_ abc _ ref 1 ( s) Vc _ abc _ ref
e sTs
GF
Figure 3-35 The simplified abc natural reference frame DB current control block
Step 2: The abc natural reference frame controller is transformed to reference frame by
using Clark Transformation (see equation (2-2)) as shown in Figure 3-36.
VcC_ _ ref 2
I cC_ _ ref
a bs
c ds
+_
I cC_
VcC_ _ ref ( s)
VcC_ _ ref _ d ( s )
e sTs
GF
Demodulator
VcC_ _ ref 2
C
c _ dq _ ref 2
I cC_ dq _ ref
I cC_ dq
e j0t
I cC_ _ ref
e j 0 t
+_
C
c _
a bs
c ds
C
c _ _ ref
VcC_ _ ref ( s)
VcC_ _ ref 1 ( s) V
C
c _ _ ref _ d
e sTs
(s)
VcC_ dq _ ref _ d ( s)
GF
Figure 3-37 Simplified dq reference frame DB current control block including the timevarying terms (modulator/de-modulator)
92
Step 4: The modulator sectors are eliminated and transformed to dq reference frame which
do not contain time varying terms (modulator/de-modulator).
Note that after these transformations, each variable in reference frame going
through the modulator system is transformed to dq reference frame which is a linear time
invariant (LTI) system. The next step is to shift the rotary factor ejt to the right hand
through three blocks: i) DB block; ii) sampling block and iii) delay block step by step until
reaching to e-jt. During this process both of the terms are nullified without containing the
time-varying term since ejte-jt =1. Therefore, the controller then can be transformed to
dq synchronous reference frame.
A. Methodology Development for Cancelation of Modulator/De-modulator
The primary question is how to eliminate the modulator/de-modulator to achieve the
equivalent linearized small-signal model system. Figure 3-38 illustrates a generic
modulator/de-modulator system. The objective here is to derive the transfer function or
state-space equation for G(s) = Ydq(s)/Udq(s) of the system. It should be noted that the
derivation process with respect to the elimination of the system modulator/de-modulator is
one of the original contributions of this thesis. It should be emphasized that two
approaches are developed from different point of views which result in the same
conclusion. The first methodology accommodates the transfer matrix state-space equation,
which avoids the presence of the imaginary components during the derivation process but
is comparably more complicated, which will be discussed in the rest of this section. The
second approach is derived from the rotary viewpoint which is also based on the state
space equations. Due to the space consideration in this thesis, this approach is not included
in the main text but included in Appendix B.
Ud
Uq
e j0t
cos
sin
sin
cos
Modulator
X A X B U
Y C X D U
e j0t
cos
sin
sin
cos
Demodulator
Yd
Yq
93
It should be noted that the transmission path between the modulator and de-modulator
is represented by a linear set of time dependent state equations. Therefore, to facilitate the
development of the linearized system of equation, they can be defined as
(3-73)
(3-74)
where, Zdq [Zd1 , Zd 2 , Zdn , Zq1 , Zq 2 , Zqn ]T , udq [ud ,uq ] ydq [ yd , yq ]
T
X i jX i (Zdi jZqi ) e j (t )
X i jX i ( Z di jZ qi )(cos (t ) j sin (t ))
X i jX i ( Z di cos Z qi sin ) j ( Z di sin Z qi cos )
(3-75)
In
the
below
paragraphs
the
complete
set
(3-76)
of
the
state
variables,
0
0
0
Xn 0
0 cos
0
sin
0
0
cos
0
X 1
sin 0
0
cos
X 0
2
0
0
0
0 sin
0
X 0
n
Z
d1
Zd 2
0
0
Z
0 sin dn
0 Z
q1
0
Zq2
0
0
cos Z
qn
0
(3-77)
(3-78)
Zd
cos In sin In d X d
d Zd
Rn
, Rn Rn
dt Z q
sin In cos In dt X dt Z q
d cos
d d sin
d
sin
;
cos
dt
dt
dt
dt
(3-79)
We can obtain
d X sin In cos In d
dt X cos In sin In dt
Z d cos In sin In d Z d
Z
q sin In cos In dt Z q
(3-80)
U jU (U d jU q )(cos (t ) j sin (t ))
u cos
u
sin
sin ud
cos uq
(a)
ud cos
u
q sin
sin u
cos u
(b)
(3-81)
X A 0 X B 0 u
X 0 A X 0 B u
(3-82)
X A 0 X B 0 cos
X 0 A X 0 B sin
sin ud
cos uq
(3-83)
sin yd
cos yq
Z dq Adq Z dq Bdqudq
ydq Cdq Z dq Ddqudq
(3-84)
(3-85)
(3-86)
95
where,
In0
B 0
C 0
D 0
, Bdq
, Cdq
, Ddq
A
0 B
0 C
0 D
A
Adq
In0
Gdd ( s) C ( sI A)[( sI A) 2 02 ]1 B D
Gdq ( s) C0 [( sI A) 2 02 ]1 B
(3-87)
Gqd ( s) Gdq ( s)
Gqq ( s) Gdd ( s)
ud
uq
cos 0t sin 0t
sin t cos t
0
0
ud ( s )
G ( s )
u
G ( s )
yd
cos 0t sin 0t
sin t cos t y
q
0
0
yd ( s )
Gdd ( s )
Gdq ( s )
uq ( s )
Gqd ( s )
yq ( s )
Gqq ( s )
(3-88)
reference frame
to dq reference frame can be done simply by adding a j term to the Laplace factor.
Although the above methodology is given for state-apace equation transforming from
reference frame to dq reference frame, similar relationship is also valid for abc natural
reference frame to dq reference frame transformation since the abc natural reference frame
96
and
constant matrix (see equation(2-2)). However, the dq reference frame rotates at a frequency
corresponding to the and abc reference frame variables. It can be noted that this
approach is valid for any linear component with transfer function derived in abc stationary
natural frame such as a grid filter, or linear controllers implemented in abc reference frame
or even in a physical grid model.
B. Small Signal Model for DB Current Controller, Sampling and Delay Blocks
As is discussed previously there is a common characteristic for the equivalent
continuous domain models references to the DB controller, sampling block and delay
block, since all of them can be written either in the lead-lag format or in the lead-lag-gain
format which are summarized in Table 3-5. Note that the third column in Table 3-5 is a
case example obtained by substituting the parameter listed in Appendix G and which will
be used for simulation verifications later in this chapter.
Table 3-5 Summary of transfer functions for DB current controller, sampling block and
delay block
Lead-Lag
a bs
c ds
DB
Block
Sampling
Block
Delay
Block
Lead-Lag-Gain
Case Study
1 Ta s
1 Tb s
T
1 s s
N1
2
(
)
1 N 2 1 Ts 1 N 2 s
2 1 N2
-0.01247s-24.94
0.000001566s 1.997
1
1 sTs / 2
1
1 sTs / 2
1
1 0.0005s
1 sTs / 2
1 sTs / 2
1 sTs / 2
1 sTs / 2
-0.0005s +1
0.0005s +1
N1 N1
(1 N 2 )
Ts
s
2
Ts
(1 N 2 ) s
2
97
functions for the target blocks can be obtained as shown in Figure 3-40. It should be noted
that so far, the transfer functions of the DB controllers are successfully transformed to dq
reference frame.
VcC_ d _ ref 2
VcC_ d _ ref 1 ( s )
I cC_ d _ ref
+_
C
c_d
GDB1
VcC_ d _ ref
VcC_ d _ ref ( s )
GSample1 ( s )
GF
GDB 2
GDB 2
I cC_ q _ ref
+_
I cC_ q
GDB1
C
c _ q _ ref 1
GF
GSample 2 ( s )
GDelay 2 ( s )
GSample 2 ( s )
GDelay 2 ( s )
GSample1 ( s )
( s ) V
C
c _ q _ ref
GDelay1 ( s )
VcC_ q _ ref ( s )
VcC_ q _ ref 2
Figure 3-40 Small signal model for DB current controller in equivalent dq reference frame
Table 3-6 Equivalent dq reference frame transfer function for the lead-lag block
reference frame
transfer function
Y ( s)
a b( s j )
G1 ( s) j G2 ( s)
X ( s) dq c d ( s j )
Lead lag
block
Y ( s)
a bs
X ( s) c ds
Main loop
Crosscoupling
G1 ( s)
G2 ( s)
ad0 cb0
;
d s 2dcs c 2 d 202
2 2
C. Simulation Verification
To test the accuracy of the DB controller in terms of equivalent transfer function and
state space equations in the dq reference frame, two step tests are performed: i) -0.02pu
step change on Uq and ii) -0.02pu step change on Ud imposed on the three target blocks in
series (which included DB control block, sampling block and computation delay block).
98
(a)
(b)
Figure 3-41 Compare Yd/Yq outputs from (i) the modulator/de-modulator system; (ii) the
small-signal equivalent transfer function and iii) the small-signal equivalent state space.
(a)UdSTEP = 0, UqSTEP = -0.02 and (b) UdSTEP = -0.02, UqSTEP = 0;
(a)
(b)
Figure 3-42 Difference between Yd/Yq outputs from (i) the modulator/de-modulator
system versus the small-signal equivalent transfer function; (ii) the modulator/demodulator system versus the small-signal equivalent state space. (a) UdSTEP = 0, UqSTEP = 0.02 and (b) UdSTEP = -0.02, UqSTEP = 0;
The results presented in Figure 3-41 indicate that there is a good agreement on the test
results, since the difference between Yd/Yq outputs from the modulator/de-modulator
system compared with the small-signal equivalent transfer function is small. In addition,
the calculated errors in the modulator/de-modulator system compared to the small-signal
99
equivalent state-space model further confirm the conclusion drawn above, as shown in
Figure 3-42 .
Equivalent State-Space Equation in dq Frame
To obtain the equivalent state-space equations in dq reference frame, the state-space
equation for the reference frame transfer function should be defined first
which is written in the lead-lag-gain format (see Appendix D),
sz Az Bx
y Cz Dx
Where A
(3-89)
T
1
k k Ta
;B
; C 1; D k a
Tb
Tb Tb Tb
Tb
Table 3-7 State space equations for DB controller in equivalent dq reference frame
Large signal state space model
dwd _ db
Dead-Beat
dt
Vcdc _ ref 1 Cdb wd _ db Ddb I cdc _ e
Block
dwq _ db
dt
dZ d _ sp
dt
Sampling
dZ q _ sp
dt
d wd _ db
dt
block
d Z d _ sp
dt
Vcdc _ ref (k ) Csp Z d _ sp Dsp Vcdc _ ref
d Z q _ sp
dt
Delay block
dM q _ dl
dt
c
cq _ ref
100
(k ) Cdl M q _ dl D V
'
c
dl cq _ ref
(k )
d M q _ dl
dt
By adding with the cross-coupling term j , we can obtain the state space equation in
dq reference frame for lead-lag block as shown in equation (3-90). The state-space
equations for DB control block, sampling block and delay block in equivalent dq reference
frame are obtained reference to Table 3-5, and equations (3-89) and (3-90), which are
summarized in Table 3-7 (in which the second column stands for the large signal models,
while the third column denotes the small signal models).
d zd
ALL zd zq BLL xd
dt
yd CLL zd DLL xd
d zq
wd ALL zq BLL xq
dt
yq CLL zq DLL xq
(a)
(3-90)
(b)
(a)
(b)
101
(c)
(d)
Figure 3-43 Compare Yd/Yq outputs from the (i) modulator/de-modulator system; (ii) the
small-signal equivalent transfer function and iii) the small-signal equivalent state space. (a)
Yd computes from UdSTEP = 0, UqSTEP = -0.02, (=0+2Ku(t), K=0.01, t=0.01s); (b) Yq
computes from UdSTEP = 0, UqSTEP = -0.02, (=0+2Ku(t), K=0.01, t=0.01s); c)Yd
computes from UdSTEP = 0, UqSTEP = -0.02, (=0+Ksin(t), K=0.01, =10Hz, t=0.01s);
(d) Yq computes from UdSTEP = 0, UqSTEP = -0.02, (=0+Ksin(t), K=0.01, =10Hz,
t=0.01s);
From the tests results shown in Figure 3-43a and b, it can be seen that the value of the
output is smaller than the exact one, which results in an under compensation on the
generated voltage reference. The problem can be solved by using the PLL output frequency
as a feedback to the TF coefficient calculation, or by increasing the controllers gain to
enhance the compensation. Therefore, it can be concluded that the increase or decrease in
frequency can cause the outputs to be under or over compensated in the equivalent dq
reference frame. Moreover, the results given in Figure 3-43c,d for sinusoidal input
frequency disturbance test reveal that the magnitude of the consequential oscillation
increases with the simulation time. Finally, it can be concluded that it is better to design the
linear controller in dq reference frame but implemented in abc natural reference frame as it
can offer better capability on frequency sensitiveness resistance.
102
(3-91)
(3-92)
Vdc (t )
mabc (t )
2
(3-93)
where, the Vc_abc (t) is the vector of VSC terminal voltages, Vdc (t) represents the dc
voltage, and the mabc (t) is the modulation signal, which can be given by,
ma (t ) M cos( )
mb (t ) M cos( 2 / 3 )
(3-94)
mc (t ) M cos( 2 / 3 )
Note that in the above equation M controls the magnitude of the AC output voltage,
is the angle of the bus voltage at PCC point that is provided by a PLL and the derivative of
103
(d/dt) is the angular frequency of the AC supply. In addition, is the relative angle
of VSC terminal voltage Vc_abc which is controlled by the control system of the converter.
Hence, the dq reference frame averaged model of three-phase two-level VSC can be
given by [2, 72],
Vdc
V
M q dc
2
2
Vdc
Vdc
Vcd M sin
Md
2
2
Vcq M cos
(3-95)
Vdc
) and the angle is equal to
2
tan -1 (Vcd /Vcq ) which is known as the modulation ratio or VSCs terminal voltage angle.
If we substitute M and into the above equation, we will obtain,
c
c
Vcd
Vcd
_ ref Vdc / Vdc 0
c
c
Vcq
Vcq
_ ref Vdc / Vdc 0
(3-96)
(3-97)
If we linearize the above equations, we can obtain the small signal model that provides
the link between the AC and the DC voltages of VSC.
g
g
g
'
VcR
VcR
_ ref (k ) VcR _ ref 0 Vdc / Vdc 0
(3-98)
Since the efficiency of the converters is high in such power applications, the
conduction losses can often be ignored. Therefore, the DC-side and AC-side terminal
quantities of the two-level VSC can be interconnected by the power balance equation as
Pdc(t)=Pac(t) , which can be written using terminal quantities as [72, 77].
3
U DC (t )iDC (t ) [Vcq (t )icq (t ) Vcd (t )icd (t )]
2
(3-99)
If we substitute equation (3-95) into the above equation, we can obtain the DC link
current as,
3
iDC (t ) [ M q (t )iq (t ) M d (t )id (t )]
4
104
(3-100)
Similarly, linearized version of the above equation can provide the small signal model
as [56],
3
iDC (icq M q M q icq icd M d M d icd )
4
(3-101)
The advantage of this model is that it does not include the complex switching states
and ignore the influence of the high frequency harmonics. Therefore, such model can run
faster than the detailed model simply acting as a voltage filtered model. However, ignoring
the high frequency harmonics will also ignore the impact of high harmonics on the grid
and the control. It should be noted that such model is only valid for frequencies up to about
one-third of the VSC switching frequency and the switched model operates in the linear
modulation region without over modulation, which is Vcd2 +Vcq2 Vdc /2 [131].
It should be reported that the above assumption will not represent the real operation.
However, a more accurate model is a very challenging task, since it depends on specific
operating point as well as the loading conditions. As reported in paper [130], a parametric
AVM model can be used for an accurate model. However, the establishment of the
parametric AVM model requires simulations of a wide range of loading conditions based
on a detailed model (closed loop simulation involving the control functions) to restore the
coefficients ((), (),
coefficient between the DC voltage and AC terminal voltage of the loading conditions; ()
represents the algebraic function coefficient between the AC terminal current and the DC
current of the loading conditions, and
expressed in terms of the dq components of the AC side voltage and the current as,
i
V
() tan 1 ( cd ) tan 1 ( cd )
icq
Vcq
(3-102)
If we replace the original constant value with these new algebraic coefficient functions
as in [130], we can obtain a more accurate AVM model suitable for various loading
conditions and operating points. However, it should be emphasized that these coefficients
which vary with the operating points and load conditions can only be evaluated for a
specific case, which cannot be formulated to represent all the operating conditions.
105
and
is an approximation.
Figure 3-44 illustrates the test circuit for the detailed VSC model. In this method a
small step change
input reference
(from
to
is set to be
Then the reference voltage is transformed to abc reference frame as an input to the
converter. To ensure the linearity of the VSC, the on-state resistance of the transistor and
diode in the converter has to be small enough, which is set as 0.001 in the case study. The
converter terminal output voltage
transformed to dq reference frame before connecting to the first order low-pass filter to
eliminate the high switching frequency harmonics. In the first stage, dq reference frame
converter terminal voltage
are obtained.
e jt
Vcq0_ ref
T0
Vcq_ref
d
q
Vca
Vca_ref
a
T b Vcb_ref
c Vcc_ref
Vcb
VSC
Vcc
a b c
T
d
Vcd
F
I
L
T
E
R
jt
Gf _ dq (s)
SIGNAL
FILTER
1
1 0.01s
Vcq
106
whole system are closed. But it should be noted that the AVM model causes an inaccurate
loop gain which should be considered in the controller design.
Test for no Load Current Icd=0 Icq=0(Vcqref Step) G_qq & G_dq
_ / _
Vcqref(Vcqref=-8.11)
Vcq
Vcd
-49.6566
-49.7538
-8.1407
0.864
-0.02
-50.5129
-50.501
-8.1169
0.8173
0.0178
-50.67
-50.6294
-8.1197
--------
--------
-50.8271
-50.7472
-8.1183
0.7498
-0.0089
-51.6834
-51.7274
-8.1319
1.0835
0.012
Test for no Load Current Icd=0 Icq=0 (Vcdref Step) G_dd & G_qd
Test
Vcdref (Vcqref=-50.67)
Vcd
Vcq
-7.9478
-7.9607
-50.6363
0.9803
-0.0425
-8.0849
-8.0799
-50.6218
1.5857
0.302
-8.11
-8.1197
-50.6294
--------
--------
-8.1351
-8.1465
-50.6271
1.0678
-0.0916
-8.2722
-8.2757
-50.6166
0.9618
-0.079
U sg
Rg
Lg
u gf
Lc
Cf
Rc
icg
VCg
Rf
107
Using the circuit analysis (KVL, KCL) for the gird circuit given above, a set of
equations can be obtained as given in the second column of Table 3-9. As concluded in
section 3.2.6.2, the method used to obtain the equivalent linearized small-signal model of
modulator/de-modulator system is also suitable for the grid model. This can be done by
replacing the A matrix with a new Adq matrix which involves the cross coupling
components j, which is listed in the third column of Table 3-9. Figure 3-46 illustrates the
small signal model of the grid.
U sRg
VcRg
-+
+
1
Rv sLv
g
icR
-+
+
Lv
1
Rv sLv
1
Rg sLg
-+
+
Cf
Lv
VcIg
1
sC f
g
fR
Lg
Lg
Cf
icIg
1
sC f
g
igR
u gfI
1
Rg sLg
g
igI
U sIg
Converter Reactor
Table 3-9 Summary of the large signal model of the AC grid in abc reference frame and
small signal model in dq reference frame
abc reference frame large signal model
RI reference frame small signal model
dic _ abc
dt
Rc
1
ic _ abc (u f _ abc Vc _ abc )
Lc
Lc
Grid Model
Grid Filter
108
dFftabc
dt
dig _ abc
dt
1
(ig _ abc ic _ abc )
Cf
Rg
1
(U s _ abc u f _ abc ) ig _ abc
Lg
Lg
d icg_ R
R
1
c icg_ R 0 icg_ I (u g f _ R V g c _ R )
dt
L
L
c
c
g
d
i
Rc g
1
c_I
g
g
g
dt 0 ic _ R L ic _ I L (u f _ I V c _ I )
c
c
u g f _ R FftR (i g g _ R i g c _ R ) R f
g
g
g
u f _ I FftI (i g _ I i c _ I ) R f
d FftR
(1/ C f )i g g _ R (1/ C f )i g c _ R 0 FftI
dt
(U sg_ R u g f _ R )
i g _ R 0 i g g _ I
dt
L
L
g
g
g
Rg g
1
d i g _ I
(U sg_ I u g f _ I )
i g _ I 0 i g g _ R
dt
L
L
g
g
3.3. COMPARISON OF THE SMALL SIGNAL LINEAR MODEL AND PSCAD SIMULATION
Figure 3-47 shows the test results of the i) grid side current and ii) filter bus voltage
responses of the VSC model, which are obtained by applying a small step change -0.56kV
G
on the input of the grid model VcR
. Note that there is a good agreement on the results of
(a)
(b)
Figure 3-47 (a) Grid side current and (b) filter bus voltage responses of i) large signal
model ii) small signal model following a input voltage step change order on VcqC = 0.56kV.
Finally, the characterization and validation tasks with regards to all the individual
component of the VSC small signal model are complete and prepared so far for the whole
system validation.
109
strong system with SCR equals to 20 following a -0.0095pu step change on Icq_ref, while
the results for a weak AC system with SCR equals to 2 following a 0.095pu step change on
Icq_ref are denoted by Figure 3-48c and d. Note that, for Figure 3-48a, the green line
denotes the reference value, the cyan line represents the output from the large signal
model, the purple line represents the and q-axis converter terminal current in the converter
reference frame
outputs from the small signal model; for Figure 3-48b, the blue line
denotes the reference value, the red line stands for the output from the large signal model,
the black line represents the d-axis converter terminal current in the converter reference
frame
outputs from the small signal model; for Figure 3-48c, the d- and q- axes
and
for an AC system with SCR=2; and a zoomed in figure for the initial section of Figure
3-48c is presented in Figure 3-48d.
The results presented in Figure 3-48a and b (where the VSC is simplified as a
controlled voltage source in PSCAD) for a strong AC system indicate that the method has
a limitation as it displays an extra error that is introduced when using the pad
Approximation. Note that the d-axis cost nearly 9 times sampling period (where the
sampling frequency is 4000Hz), while the q-axis settles within about 7 sampling periods.
However, the discrete controller itself only cost 2 sampling time to settle down as deigned.
Despite these limitations, the dynamic simulation results reveal a good agreement on the
system performance. However, in a weak AC system, Figure 3-48 (c) and (d), the dynamic
simulation results reveal a better agreement on the system performance. Both of the d and
q axes cost 4 times sampling periods to closely track the current output from the PSCAD
large signal model. In addition, the DB controller presents a slower current tracking ability
in the weak AC system compared to the strong AC system.
(a)
110
(b)
3.3. COMPARISON OF THE SMALL SIGNAL LINEAR MODEL AND PSCAD SIMULATION
(d)
(c)
Figure 3-48 Comparison of the converter output current ICcd , ICcq using the (i) large signal
model simulating with VSC simplified as a controlled voltage source; (ii) the small-signal
model for AC system with different SCRs.
It can be easily seen from the results shown in Figure 3-49 more harmonics with
switching frequency appeared in the plots (due to switched converter modeling) compared
with that of the results shown in Figure 3-48, However, it can be clearly concluded that the
linear small signal model for the DB controlled VSC demonstrates an accurate
approximation of the detailed model developed in PSCAD within the frequency range of
interest.
Iccd-PSCAD(S)
0.1
Iccq-PSCAD(S)
0.08
Iccd-Matlab
0.06
Iccq-Matlab
0.04
Iccd-PSCAD(D)
0.02
Iccq-PSCAD(D)
Iccd-PSCAD(D-FIR)
-0.02
-0.04
0
Iccq-PSCAD(D-FIR)
0.02
0.04
Time (s)
0.06
0.08
(a)
111
0.12
Iccd-PSCAD(S)
0.1
Iccq-PSCAD(S)
0.08
Iccd-Matlab
0.06
Iccq-Matlab
0.04
Iccd-PSCAD(D)
0.02
Iccq-PSCAD(D)
0
-0.02
-0.04
0
0.002
0.004
0.006
Time (s)
0.008
0.01
112
I cdc _ rec
Rdc _ rec
I cdc _ inv
Ldc _ rec
Rdc _ inv
Ldc _ inv
I dc _ rec
I dc _ inv
Vdc _ rec
Cdc _ rec
Vcdc
Vdc _ inv
Cdc _ line
Cdc _ inv
dVdc _ rec
dt
dI cdc _ rec
dt
Vdc _ rec
Ldc _ rec
I dc _ rec
Cdc _ rec
I cdc _ rec
d Vdc _ rec
Cdc _ rec
dt
Ldc _ rec
Ldc _ rec
d I cdc _ rec
dt
dt
dVdc _ inv
dt
I cdc _ inv
Cdc _ inv
Vdc _ rec
Ldc _ rec
I dc _ rec
Cdc _ rec
I cdc _ rec
Cdc _ rec
Ldc _ rec
Ldc _ rec
dt
Cdc _ line
Cdc _ line
dt
Cdc _ line Cdc _ line
dI cdc _ inv
d I cdc _ inv
dt
I dc _ inv
d Vdc _ inv
Cdc _ inv
dt
I cdc _ inv
Cdc _ inv
I dc _ inv
Cdc _ inv
Note that, although the established DC link model appears to represent a mono-polar
link, it is also applicable to represent bi-polar links. But we need to apply the following
two rules to specify the line, cable and shunt capacitors parameters when representing a
bipolar link [134],
sum the respective quantities for each line and capacitance to specify the total
line, cable resistance, inductance and capacitance.
specify the shunt capacitance by summing the capacitance between the positive
polar to neutral point and the capacitance between the negative polar to neutral
point.
113
Figure 3-51 illustrates the block diagram of the small signal model for the DC link,
which will be verified against the PSCAD large signal model. In the verification study, the
parameters given Table 3-11 are used.
500.0 [uF]
7[ohm]
0.5968 [H]
26.0 [uF]
7[ohm]
0.5968 [H]
500.0 [uF]
I cdc _ rec
I dc _ rec
Vcdc
Vdc _ rec
1
sCdc _ rec
1
sLdc _ rec Rdc _ rec
I cdc _ rec
I cdc _ inv
1
sCdc _ line
Vcdc
Vdc _ inv
I cdc _ inv
1
sLdc _ inv Rdc _ inv
I dc _ inv
1
sCdc _ inv
Time (seconds)
Figure 3-52 The converter currents behind the DC capacitor from both of the rectifier side
and inverter side (circle: Large signal model; rectangular: Small signal model)
114
U dc _ inv
The results shown in Figure 3-52 are obtained by applying a step change 0.577kA to
the input of the DC grid model Idc_rec, which indicate the accuracies of the small signal
model since the large signal model and the small signal model currents are very similar.
U iEXT
U j
X i J XXi
0 J ZXi
Y j X J XXj
j
0 J ZXj
J XZi X i J XUi
J
U i XUEi U iEXT
J ZZi Z i J ZUi
J ZUEi
Yi
U i
J XZj X j J XUj
J XUEj
U
J
j
jEXT
J ZZj Z j J ZUj
ZUEj
U jEXT
0 J ZXi
X j J XXj
0 J ZXj
J XZi X i J XUi
J XUEi
i
J
U iEXT
J ZZi Zi J ZUi
ZUEi
J XZj X j J XUj
J XUEj
J
U jEXT
j
J ZZj Z j J ZUj
ZUEj
(3-103)
(3-104)
115
(3-105)
U j Yi ;
Y j JYXj X j JYZj Z j JYUj U j ;
(3-106)
U i Y j ;
0 J ZXi
J XZi X i J XUi
J
U i XUEi U iEXT
J ZZi Zi J ZUi
J ZUEi
(3-107)
X
JYZi i JYUi U i Yi
Zi
(3-108)
J XZj X j J XUj
J XUEj
U j
U jEXT
J ZZj Z j J ZUj
J ZUEj
(3-109)
0 JYXi
X j J XXj
0 J ZXj
0 JYXj
X j
JYZj
JYUj U j Y j
Z j
(3-110)
0 U j Yi
(3-111)
0 Ui Yj
(3-112)
0
xi J XXi
x 0
J XXj
j
0 J ZXi
0
0 0 J ZXj
0 J
0
YXi
J YXj
0 0
0 0
0
0 0
0
J ZX (C )
J XZ (C )
J XUE (C )
X i J XUei
0 J XZj 0 0 0
J XUj X j 0
J ZZi
0
0
0 J ZUi
0 Z i J ZUei
0 J ZZj 0
0
0
J ZUj Z j 0
J YZi
0 0 J YUi
0 Yi J YUei
0 J YZj
0 0
JYUj Y j 0
0
0 0
0
U i JUUei
0
0
0
0 U j 0
J XZi
0 0 J XUi
0
J XUej
0
J ZUej uiext
0 0
J YUei
0
JUUei
0
(3-113)
u jext
J ZUE (C )
J ZZ (C )
and
T
X (C ) X i
X j ;
Z (C ) Zi
Z j
Yi
Y j
U i
U j ;
Note that this allows us to obtain a single set of equation from two separate sets of
equations. Hence, using this block interconnection approach, all the sub-components of the
116
system can be mathematically linked together, which completes the interconnection task of
the subsystems.
Ic
u *f
RECTIFIER
C *f
INNER
CURRENT
CONTROLLER
AC Source
I cd _ ref
DC Link
DC Source
I cq _ ref
Figure 3-54 The base model for the converter side controller
Similar to the results given in Figure 3-49, Figure 3-55 also reveals a good agreement
on the system dynamic that includes the dynamics of the DC link, which is needed for the
outer loop controller design. The results are obtained by applying a 0.095pu small step
change on the q-axis of the input current reference of the converter.
0.14
Iccd-PSCAD(D)
0.12
Iccq-PSCAD(D)
0.1
Iccd-Matlab
0.08
Iccq-Matlab
0.06
Iccd-PSCAD(D-FIR)
0.04
Iccq-PSCAD(D-FIR)
0.02
0
-0.02
-0.04
-0.06
0
0.05
0.1
Time (s)
0.15
0.2
(a)
117
0.15
Iccd-PSCAD(D)
Iccq-PSCAD(D)
0.1
Iccd-Matlab
Iccq-Matlab
0.05
-0.05
0
0.01
0.02
Time (s)
0.03
0.04
(b)
Figure 3-55 Comparison of the converter output currents ICcd , ICcq using the (i) largesignal model simulated in the detailed VSC model including the DC link, and (ii) the
small-signal model for a weak AC system together with the DC link. The representations
of the abbreviation in terms of D and D-FIR are the same as in Figure 3-49.
3.5 Conclusion
This chapter developed a complete small signal model for the DB controlled VSC
including the DC link. In order to obtain the small signal model for the three phase discrete
controller, four steps are involved. In the first step, the PadApproximation is adopted to
transfer the discrete controller to the continuous domain controller. In the second step, the
abc three phase continuous domain controller is transformed to the
reference frame
reference frame
reference frame to dq
reference frame, a frame transformation approach to eliminate the modulator and demodulator block was proposed in this chapter which has also been verified. It was noticed
that, although the implementation of the PadApproximation introduces a degree of error,
it can still capture the main properties of the system, which can be a valid tool for the
controller design. The small signal model development is presented systematically in the
chapter including the derivation of the entire system equations which are then verified
against the PSCAD simulations. In the simulation studies both of the large signal model
118
3.5. CONCLUSION
and small signal model are verified for every individual system components, as well as the
entire integrated VSC system. It is demonstrated that the contribution of this study can be
used as a solid foundation work that will be used in the remaining chapters of this thesis.
119
4.1 Introduction
Integrating a VSC-HVDC into a weak AC system which is characterized by large grid
impedance and its corresponding control system design were a challenging task for the
researchers. Traditionally, a cascade control scheme involving outer loop control and inner
current loop control is utilized. In such an approach, the performance of the system
predominantly depends on the performance of the inner current control, since a fast inner
120
current controller can provide sufficient frequency space for the outer power loop
controller design.
It can be noted that the behavior of the inner current controller highly depends on the
interaction of the current control, PLL control and filter bus capacitor. Specifically, under
weak AC condition, the dynamics of PLL and filter bus capacitor play a critically
important role. This is due to the success of a control system is based on the accuracy and
speed response of PLL. However, the control of PLL is also highly governed by the
voltage that is being locked to. In such a case, however, the filter bus voltage fluctuates
dramatically, which is extremely sensitive to the converter current. This posed significant
challenges to the PLL controller design. Therefore, it can be concluded that it is a highly
complex control issue. Although two papers [54, 68] addressed this problem using
simulation studies, it has not been studied systematically and quantitatively so far.
One of the motivations of this thesis is to perform a complete analysis for VSCs
operating in weak AC grids by using a detailed accurate mathematical model from the
small signal stability point of view. Furthermore, another motivation was to develop a
fixed parameter controller with a simple structure (e.g. PI or PID). Although this is
constrained by the limited bandwidth of the inner current controller it can still perform
sufficiently well for a wide range of operating points.
change. In such a system, the plant to be controlled is only the converter reactor which is a
first-order system. However, the controlled plant becomes a third-order system since the
dynamics of the grid impedance and the filter capacitance are taken into account, such as in
a weak AC grid. The explicit circuit analysis in a weak AC grid is given in the following
section.
121
Egg
*
g
ugf
Lc
Cf
g
c
Vcg
1
sL*g
Egg
*1 L*g Lg Ltr
1
sC f
Rf
ugf
Rf
icg
1
sLc
Vcg
i gg
g
z12 ig
z22 icg
1
*
sLg sC R f
z11 z12
f
1
21 z22
Rf
sC f
Hence, the currents can be calculated as
igg Y11
g
ic Y21
1
sLc
Rf
sC f
1
Rf
sC f
Y12 Eg z11
Y22 Vcg z21
g
(4-1)
z12
z22
Egg
g
Vc
(4-2)
(4-3)
If we expand the above equation, the individual current functions can be given by
g
g
g
ig Y11 Eg Y12Vc
g
g
g
ic Y21 Vc Y22Vc
122
(4-4)
1 sC f R f s 2C f L*g
icg
Y22 3 *
Vcg
s Lg Lc C f s 2C f R f ( L*g Lc ) s( L*g Lc )
(4-5)
The above equation shows that the plant to be controlled is a third-order system by
taking into account the grid impedance L*g and filter capacitance Cf. If it is applied to a
strong AC grid, these two factors can be ignored and equation (4-5) is simply reduced to
1/sLc only, which is a first-order system.
Ia_ref
Ia_ref (k)
u va_1
G db (z)
Execution
-1
delay z
GZOH(z)
Ia
Gp(z)
T
GF(Z)
123
(a)
(b)
Figure 4-3 Zeros/poles in z-domain for (a) weak system (b) strong system of a DB current
controlled VSC system
Frequency response of Ia(s)/Iaref(s)-open loop
(a)
(b)
Figure 4-4 (a) Bode plots and (b) step responses of a DB current controlled VSC system
124
of GM >= 6dB and PM >= 40o respectively. Figure 4-5 shows the flow chart of the process
for designing such controllers.
Determine OP points
(PQ,SCR,R/X Ratio)
Newton Raphson
(Initial Condition Calculation)
System Matrices
(A,B,C,D)
i>=36?
Y
A sets of systems
U
V ( s) U f _ INV ( s)
P( s )
f _ REC ( s )
(
;
; dc
;
)
I cq _ REC ( s) I cd _ REC ( s) I cq _ INV ( s) I cd _ INV ( s)
PI/PID Compensation
(GM>=6dB;PM>=40deg)
EMTDC Verification
End
125
obtained with different radiuses. This can help us to understand how the VSC terminal
voltage will affect the PQ capability chart of VSC system.
The Automatic Access Standard in Australia states that [138]:
Any level of active power and any voltage at the connection point within the limits
established under cause S5.1a.4 without a contingency event, must be capable of supplying
and absorbing continuously at its connection point an amount of reactive power of at least
the amount equal to the product of the rated active power of the generating system and
0.395.
To meet the above requirement, an additional vertical purple line in Figure 4-6a is
added. In this study, the DC voltage is chosen reasonably high, which is not a limiting
factor. Therefore, only two limits are imposed on the system as shown in Figure 4-6b. The
operating scenarios in terms of VSC PQ capability considered in this section are shown in
Figure 4-6b which involves different loading conditions (Pmax, Phalf and P0) and various
level of power factors (leading, unity and lagging). In addition, the grid conditions that
SCR and X/R ratios are also considered as operating scenarios which are listed in Table
4-1.
For weak AC system, the SCR is chosen as SCR=2, while for a strong system it is
chosen as SCR=7.5. For X/R ratios, the pure impedance condition (at 90o) and line resistor
Rg with an impedance angle 75o are considered. The combination of these conditions
allows us to study 36 different operating scenarios, which covers a wide range of operating
points in such grid-connected VSC system.
DC Voltage Limit
(a)
(b)
Figure 4-6 PQ capability chart (a) demonstration chart; (b) defined operating points
126
7.5
: (|Z|cos)( )
: (|Z|sin/2 f) (H)
|Z|
Angle of |Z|
26.042
90o
0.083
75o
6.740
0.080
90o
0.022
75o
1.797
0.021
6.944
X n 1 X n
n 0,1,2...
F(Xn)
F'(Xn)
Ic
Ig
U f E g 2 Pgc jQgc
Ig
Zg
3
U f
Ig Ig I f
jB U f
3
3
U f
Vc
(4-6)
Figure 4-7 shows a single line of the AC side of the converter that includes a high-pass
filter. Using the Kirchhoffs voltage law and the power balance relationship, the converter
current, the filter current and the grid side current can easily be obtained as shown in Table
127
4-2. If we expand these equations and apply linearization, the F(X) and its Jacobian Matrix
equations F(X) can be defined which are summarized in Table 4-3, which form the basis
for steady-state operating point calculation.
No.
Pc
3
(VCR I CR VCI I CI )
2
3
Pc (VCR ICR _ 0 ICR VCR _ 0 VCI I CI _ 0 I CI VCI _ 0 )
2
Qc
3
(VCI I CR VCR I CI )
2
3
Qc (VCI ICR _ 0 I CR VCI _ 0 VCR I CI _ 0 I CI VCR _ 0 )
2
I CR
Rc
(U _ fR V_ cR ) ...
R X c2
2
c
ICR
Xc
(U _ fI V_ cI )
R X c2
Rc
X
(U _ fR V_ cR ) 2 c 2 (U _ fI V_ cI )
2
R Xc
Rc X c
2
c
2
c
I CI
Xc
(U _ fR V_ cR ) ...
Rc2 X c2
ICI
Rc
(U _ fI V_ cI )
R X c2
Xc
R
(U _ fR V_ cR ) 2 c 2 (U _ fI V_ cI )
2
R Xc
Rc X c
2
c
2
c
Pgc
3
(U fR I gR U fI I gI )
2
3
Pgc (U fR I gR _ 0 I gR U fR _ 0 U fI I gI _ 0 I gI U fI _ 0 )
2
Qgc
3
(U fI I gR U fR I gI )
2
3
Qgc (U FI I gR _ 0 I gR U fI _ 0 U fR I gI _ 0 I gI U fR _ 0 )
2
I gR
Rg
R X g2
2
g
Xg
R X g2
2
g
I gI
( E_ gR U _ fR ) ...
( E_ gI U _ fI )
Xg
Rg2 X g2
Rg
Rg2 X g2
I gR
( E_ gR U _ fR ) ...
( E_ gI U _ fI )
I gI
Rg
R X
2
g
2
g
Xg
R X
2
g
2
g
(E_ gR U _ fR )
(E_ gR U _ fR )
Xg
R X g2
2
g
Rg
R X g2
2
g
(E_ gI U _ fI )
(E_ gI U _ fI )
I gR I cR B U fI
I gR I cR B U fI
10
I gI I cR B U fR
I gI I cR B U fR
11
U fR U f 0 cos( _ PLL )
12
U fI U f 0 sin( _ PLL )
128
Pgc jQgc
Eg 0o
Z g
Pfc jQ fc
U f _ PLL
Pc jQc
Vc
Zc
AC network Side
VSC Side
Ig
If
Ic
jB_ REC
Figure 4-7 Single-line diagram of AC side converter including a high-pass filter
(a)
(b)
Figure 4-8 Influence of the DC link model: (a) pole-zero map of Iq(s)/Iqref(s); (b) Bode
plots of Iq(s)/I qref(s)
Figure 4-8 shows the pole-zero map and Bode plots of the inner DB current controller
with and without the DC link model. As it can be seen in the Figure 4-8a that there are two
more pole-zero pairs presented in the model with DC link than the model without DC link
at very low frequency, which are nearly cancelled with each other. However, at frequency
129
around 30rad/s, a pair of cancellable pole-zero pair exists in the model without DC link,
while in the model with DC link there is only one pole pair which is un-cancellable (see
the circled area in the Bode plot). In addition, small magnitude difference was observed in
the Bode plots of two systems. The results reveal that it is possible to develop more
accurate and realistic controllers when the dynamics of DC link are being considered.
4.3.3.1 Controller Design Circuit for Rectifier
The model adopted for the outer power controller design at the rectifier side is shown
in Figure 4-9. The DC side transmission line is modelled as a T-section ended with a
constant DC voltage source.
u *f
RECTIFIER
PREC
C *f
INNER
CURRENT
CONTROLLER
AC Source
DC Link
DC Source
I d _ ref
Pref
isq _ max
I q _ ref
PI
isq _ max
Note that the control scheme used here is a cascade control strategy. The DB control
scheme is adopted as the inner current control for both of the rectifier and inverter sides,
while for the outer loops, the rectifier is responsible for controlling the active power and
the constant filter bus AC voltage, and the inverter controls the DC voltage and the filter
bus AC voltage to be constant.
The above circuit diagram is the basis for rectifier side controllers design. Due to the
fact that the objective of the controller design is to control the DC voltage at the inverter
end, the inverter can be simplified as a constant DC voltage source. Although the
assumptions made here that the constant DC voltage controller is implemented ideally and
the response speed of the DC voltage controller is sufficiently fast does not 100 precent
hold it is still adequate to be used as a start point in the controller design.
130
Eigenvalue
1,2
-2.04107314j
50
3,4
-7.91103307j
48.834
5,6
-1.6251033.03103j
0.473
482.18
7,8
-1.621032.48103j
0.547
394.78
9,10
-31033.04102j
0.994
48.345
11,12 -2.791026.64102j
0.387
105.69
13,14
-78.7263.6j
0.286
41.955
15,16
-6.094182j
0.033
28.949
17,18
-25.37918.384j
0.81
2.926
19,20
-12.2127.37j
0.407
4.3558
As shown in Table 4-5, the above mentioned modes 15, 16, 19 and 20 are also greatly
influenced by the DC link states. Therefore, any attempts to improve these modes must
131
take into accounts these DC link states. PLL is mainly associated with the modes 17 and 18.
While, the time constant of the AC voltage signal filter and the AC grid states are the main
contributors to modes 11-14. The participation factor table also suggests that the states
with respect to the circuit in adjacent area of VSC terminal (converter-reactance and grid
voltage filter) and the DB controller have dominant impacts on modes 5-10. Lastly, the
modes 1-4 are solely determined by the DB control block. In the following section, the
eigen-sensitivity analysis in terms of these states will be explicitly discussed.
Table 4-5 Participation factors of dominant state variables for selected modes of the inner
DB current controller
Sub-System State
Grid Controller
DB control
AC network
DC grid
Mode
PLL
Tf
DeadBeat
Sampling
Delay
Con-Reactance
Cf
AC grid
C_dc_rec
Ldc_rec
Cdc_line
Ldc_inv
1
2
100
3
4
5
6
7
8
9
10
11
12
0.892
19.1
3.079
2.772
5.352
10.87
11.55
46.39
13
15
14
16
6.24
36.2
1.54
2.58
5.1
5.32
12.3
29.3
1.49 1.455
26.33
48.58
23.63
17
19
18
20
86.4 12.75
1.48 0.609
0.1
0.19
0.34
2.68
3.35
2.68
0.09
2.67
1.003
39.9
21.28
1.087
23.37
132
down as the consequence of the high frequency modes evolving to the low frequency
region. This makes the design of outer loop controller design even more complex and
challenging.
(b)
(a)
Figure 4-10 (a) In the DB current controlled VSC with P=Pmax (a) eigenvalue maps with
SCR reduced from 7.5 to 0.5 at a step of -1 (i.e. A to G); (b) Bode plots of the TF (blue:
SCR=7.5, red: SCR=4.5, cyan: SCR=2).
Figure 4-10b shows the Bode plots of the transfer function for the DB current
controlled VSC. The upper left figure plot shows the magnitude plot for power. Note that,
only the magnitude plot is shown here for simplicity purpose. It can be observed that the
magnitude of the strong system has a higher gain compared with that of the weak AC
system. However, it should be noted that this is not a regular pattern for all loading
conditions. The test results of the zero loading condition (static synchronous compensator
(STATCOM) mode) show a reverse pattern. Nevertheless, the Bode plots of the filter bus
AC voltage (the lower plots of Figure 4-10b) confirm that it is more sensitive to the power
variation in a weak AC grid. Furthermore, the upper right and lower left plots also confirm
the fact that the cross-coupling effects increase as SCR reduces.
Figure 4-11 shows the time-domain simulation results in MATLAB and
PSCAD/EMTDC following a 0.01pu step change on iqref at t=0s. The PSCAD simulation
confirms the accuracy of the small signal model again. In addition, the results also suggest
that the time domain simulation agrees with the conclusion drawn from the frequency
response analysis. For example, the speed of the current tracking response of the weak AC
grid is slower compared with that of the strong AC grid, since the weak AC system has a
133
narrower bandwidth. As it is can be seen in the middle power plot of Figure 4-11, the
system with a higher SCR has more reduction since the un-shown phase angle plot for the
power curve starts from -180o. Furthermore, the third plot in Figure 4-11 shows that a
much higher filter bus voltage (black and red plots) increase for weak AC system is
Icq(kA)
observed compared with that of the strong system (see blue and cyan plots).
0.02
0
-0.02
0
0.01
0.02
0.03
Time(Seconds)
0.04
0.05
0.01
0.02
0.03
Time(Seconds)
0.04
0.05
0.01
0.02
0.03
Time(Seconds)
0.04
0.05
P(MW)
1
0
-1
0
Uf(kV)
0.2
0
-0.2
0
Figure 4-11 The time-domain simulation results: PSCAD SCR=7.5 (solid cyan line);
Matlab SCR=7.5 (dash-dot blue line); PSCAD SCR=2 (solid black line); Matlab
SCR=2(dash-dot red line).
134
I1
I3
I4
I2
I3
I4
Figure 4-12 Four main low frequency modes for A: rectifier mode with a 90o angle; B:
rectifier mode with a 75o angle; C: rectifier mode with a 60o angle; a: inverter mode under
90o angle; b: inverter mode under 75o angle; c: inverter mode under 60o angle;
40
|G PId(jw)|(dB)
|G PIq(jw)|(dB)
40
30
20
10 0
10
10
30
20
10 0
10
(rad/s)
(rad/s)
40
|G UfId(jw)|(dB)
30
|G UfIq(jw)|(dB)
10
20
10
0 0
10
10
(rad/s)
30
20
10 0
10
10
(rad/s)
Figure 4-13 The Bode plots of the DB current controlled VSC. blue: rectifier mode under
90o angle; cyan: Rectifier mode under 75o angle; red: rectifier mode under 60o angle;
magenta: inverter mode under 90o angle; black: inverter mode under 75o angle; green:
inverter mode under 60o angle;
The Bode plots of the DB current controlled VSC are shown in Figure 4-13. As it can
be observed in the top-left figure, the rectifier mode operation, the system with higher
135
impedance has the highest gain in the active power loop. At the inverter side the system
with a higher resistance has the highest gain. It should also be emphasized here that a main
trend can be observed the inverter has a relatively higher gain in comparison with that of
the rectifier. Alternatively, for a rectifier, it can be interpreted as to transmit the same
amounts of power the system with a higher resistance experiences larger current variation
which results in the current level reaching its limit. Therefore, it can be concluded that the
minimum SCR required in such system is higher than the more inductive system, while
such trend reverses at the inverter reference to the Bode plots analysis. Moreover, the
inverter is more robust when embedded in a weak AC system since it generally has a much
higher gain than that of the rectifier side as shown in Figure 4-13a. The same conclusions
have also been reported in ref [140].
It was also observed that the rate of change in P/Id or Uf/Id has a similar trend as
in P/Iq. However, under above two conditions, the change in the relative gain is found to
be not as big as that of in the power loop. Similarly, the rate of change in Uf/Iq has an
inverse trend. The more resistive system has the highest gain reference to the rectifier side,
vice versa for the inverter side. In conclusion, the rectifier mode has a higher gain, but the
inverter mode has a wider relative gain change. The time domain simulation results shown
in Figure 4-14 also confirm the above conclusions.
Figure 4-14 Step responses of DB current controlled VSC for the rectifier (left) and the
inverter (right). rectifier mode at angle 90o (blue); rectifier mode at angle 75o (cyan);
rectifier mode at angle 60o (red); inverter mode at angle 90o (magenta); inverter mode at
angle 75o (black); inverter mode at angle 60o (green);
136
I1
(a)
(b)
Figure 4-15 For DB controlled VSC with SCR=2 (a) the eigenvalue map with P equals to
be 1pu (A), 0.5pu (B) and 0pu (C); (b) Bode plots of the TF with different loading
conditions. blue: P=1pu, cyan: P=0.5pu, red: P=0pu).
It can also be observed that the magnitude of the transfer function of the cross-coupling
behaviors increased dramatically with loading conditions (see upper right and lower left
plots in Figure 4-15b). This indicates that more interaction between active and reactive
power loops could be observed once the outer power loops are closed.
d) Operation under Various Power Factors at SCR=2
As shown in Figure 4-16 modes I1 and I2 move to the opposite direction when varying
the operating modes of VSC from receiving reactive power from grid to sending reactive
power to grid. However, Mode I1 is the main limitation of the bandwidth of the GPIq(s).
Note that the system with a leading power factor (case C in Figure 4-16a) has the
narrowest bandwidth. Moreover, the impacts of changing power factor are more significant
on the active power related loop (left column of Figure 4-16b), whereas the impacts on the
reactive power loop are comparably small (right column of Figure 4-16).
137
I1
I2
(a)
(b)
Figure 4-16 (a) The eigen-value map and (b) Bode plot of the transfer function for the DB
current controlled VSC at SCR=2 and different power factors (Q=Qmax, Q=0 and Q=Qmin)
under full load condition. (dark blue: lagging power factor, light blue: unity power factor,
red: leading power factor).
It can be concluded from these results that SCR, X/R ratios, loading conditions and
various power factors all have significant impacts on system small signal stabilities.
Therefore, it is expected that these analyses can provide a better understanding on the
characteristics of the plant to be controlled and facilitate ongoing controller design.
4.3.3.3 Power Controller Design for Rectifier
As it was mentioned at the beginning of this chapter, the primary aim of this chapter is
to design a linear controller with a set of fixed parameter that can work under a wide range
of operating conditions. However, it is important to note that the inner DB current
controller of the weak AC grid connected VSC has a limited bandwidth as mentioned
before.
A. Frequency Response Analysis of all the Covering Operating Points
It is important to examine the characteristics of the transfer functions for 36 scenarios
mentioned earlier before conducting the controller design.
As shown in Figure 4-17a, b, c and d, the gain varies from 31.97dB to 36.79dB, which
corresponds to a difference of 4.82dB (factor of 1.7418 which is acceptable under steady
state). Under the weak AC system condition, the gain varies from 31.97dB to 36.77dB,
which is much wider than the case of a strong system (in which the gain varies from 36.79
dB to 36.3 dB).
138
(a)
(b)
(c)
(d)
Figure 4-17 The analysis of the transfer functions using Bode plots and at 36 operating
conditions in total (a) SCR=2 at 90o (b) SCR=2 at 75o; (c) SCR=7.5 at 90o; (d) SCR=7.5 at
75o;
It can be seen in Figure 4-17a that the operating points under the same loading levels
belong to the same group. In other words, the system delivering the same amounts of
power have similar gains. For example, the operating points 7-9 which transfer zero power
to the DC side have the highest set of gains. Similarly, the operating points 4-6 (half
loading condition) have a set of moderate gains in the middle group, and the remaining 3
operating points which are the full-loading modes have the lowest set of gains. Moreover,
the operating modes which absorb the maximum reactive power from grid has a
comparably higher gains compared with the modes where the maximum reactive power is
139
sent to grid. Figure 4-17b has a similar trend with that of Figure 4-17a, while the rest of the
results c and d of Figure 4-17 have reverse trends compared with the pattern observed in
Figure 4-17a, b.
B. RHP Zero Limitation on Outer Loop Bandwidth
It is well known that the poles of a system have significant impact on system stability.
However, there are limited studies on how the zeros of a system limit the controller design
of a VSC-HVDC system. It was reported in [63] that the bandwidth of the closed-loop
system should be chosen at least lower than half position of the RHP zero. A more accurate
analysis was given in [141] indicated that the achievable bandwidth b of the outer loop
controller is significantly limited by the RHP zeros, which are
Re( z ) Im( z )
| z | /4
(4-7)
where, z denotes the complex form of the system zero which is the closest to the origin.
It can be summarized that to have RHP zeros close to origin is inappropriate, and it is even
worse to have RHP zeros closer to the real axis rather than the imaginary axis [141]. It was
also found in [141] that the root causes resulting in the RHP zeros is the competing effects
among sub-components of a system.
As given in the figure below, there are 4 RHP zeros in Figure 4-18a and 1 RHP zero in
Figure 4-18b involved in case C and F locating in different frequency range. This caused
1080o change crossing low frequency to high frequency (see Figure 4-17a). It should be
noted that the cases A-I given in Figure 4-18 correspond to the cases 1-9 given earlier in
Figure 4-6b. Note that, these cases impose limits to the achievable bandwidth of the outer
loop controller. The case F has the closest zero to the origin, which is located at 44.28
Np/s+264.4 rad/s. This indicates that the designed outer loop controller should not exceed
the value |Z|/2 which corresponds to 134.04 rad/s. In other words, a controller with an
extremely fast responding speed is hard to achieve. Similarly, 3 RHP zeros included in the
cases A, B, D, E, and G, induce 720o change varying from low frequency to high
frequency. Similarly, only 1 RHP zero in the cases H and I produce 360o change in the
phase plot.
140
(a)
(b)
Figure 4-18 Right half-plane zeros for the transfer function of the open loop power
controller
C. Power Controller Design and its Small Signal Model
As indicated previously, the main aim of this sub-section is to design a power
controller and then to develop its small signal model. Therefore, the design methodology is
proposed first, which is not just limited to the power controller design, followed by step
tests to evaluate the performance of the designed power controller.
a) Design Methodology
The proposed design methodology is composed of four steps:
Step 1: Observe the frequency responses of P(s)/Iqref(s), and then choose the cases with
the lowest bandwidth, the highest bandwidth and the bandwidth in the middle of those two.
Step 2: Add a low-pass filter to ensure the open-loop gain is substantially less than DC
gain at all frequencies (above the gain cross-over frequency, if necessary).
Step 3: Apply the classic loop-shaping technique to the system obtained in Step 2, and
obtain the controller parameters Kp, Ki or Kp, Ki, Kd depending on the type of controller
designed.
Step 4: Evaluate the controller performance thoroughly against the detailed model
established on PSCAD/EMTDC simulator considering different size steps, operating points
and SCRs.
In the first step, the frequency response plots presented in Figure 4-17 are analysed. In
this design, cases 4, 6 and 9 are chosen as a reference for the design as given in Figure
4-19.
141
Phase(deg)
Magnitude(dB)
10
10
10
10
500
0
-500 -1
10
PhalfQmax
PhalfQmin
P0Qmin
0
10
10
10
10
Frequency (rad/sec)
Figure 4-19 Bode plots of the selected cases for the power controller design
In step 2, a first-order low-pass filter (see Figure 4-20a) is added to the system. This
indicates that an additional attenuation (where 20dB/decay rolling off when exceeding
60rad/s) is added to the gain plot (see Figure 4-20b). To meet the stability criteria
mentioned in section 4.3 (in which the gain margin (GM) and the phase margin (PM) of
the compensated open-loop systems are required to be within the ranges GM >= 6dB and
PM >=
designed in Figure 4-20c, these transfer functions yield a phase margin of better than about
62o at a gain cross-over frequency of about 7.46 rad/s (a little bit lower than the first dip in
the magnitude plot that is actually a zero caused by DC link) and a gain margin better than
about 27 dB (see Figure 4-20d). This suggests that the corner frequency may exceed the
frequency where zero occurs, which is caused by the dynamics of the DC link, and may
cause an inadequate performance or even instability due to the ignored dynamics, therefore
the dynamics of DC link should be considered in the controller design.
(a)
142
(b)
(d)
(c)
Figure 4-20 The Bode plots with (a) low-pass filter; (b) DB current controlled VSC with a
low-pass filter; (c) designed PI compensator; (d) open loop transfer function of P(s)/Pref(s).
The closed loop responses of P(s)/Pref(s), including all 36 loading scenarios being
subject to a 1pu step change in the power reference are employed to evaluate the
performance of the designed power controller which is given in Figure 4-21.
1.4
1.2
Case 9: Overshoot(%):4.35
At time (seconds):0.219
P(pu)
1
Case 2:
Settling time (seconds):0.506
0.8
0.6
0.4
0.2
0
0
0.1
0.2
0.3
0.4
0.5
Time (s) (seconds)
0.6
0.7
0.8
143
Not that the responses are satisfactory as they settle to within 2% of the final value at
about 0.506s (the maximum settling time for the entire transfer functions). The overshoot
(height of peak relative to the final value) is about 4.35% which is within the acceptable
range (< 10%) suitable for the reference tracking.
b) Small Signal Model of Power Controller
Figure 4-22 shows the control block diagram of the power controller, where Pm
denotes the measured value of power and PX represents the power state. Tf_P, Kp_P and
Ki_P denote the control parameters derived from the previous sub-section. Hence, the
corresponding small signal model is given by the below equations.
d Pf
dt
Pref Pm Pf
d PX
Ki _ P Pf
dt
(4-9)
iqref K p _ P Pf PX
(4-10)
Pm
Pref
+-
(4-8)
Tf _ P
1/ (1 T f _ P s)
Pf
K p_P
iqref ( s )
+
Ki _ P / s
PX
4.3.3.4 AC Voltage Controller Design for the Rectifier and its Small Signal Model
Note that the power control loop should be taken into account when designing the AC
voltage controller to improve the accuracy.
Ysp1 - E1
Gc1
+
U1
G p11
Y1
G p12
Ysp 2
E2
Gc 2
U2
G p 21
G p 22
Y2
Figure 4-23 A typical MIMO system, in which the hidden feedback loop is shown in red
lines
144
Figure 4-23 shows the typical multi-in multi-out (MIMO) system, where Gp11, Gp22 are
two main inner loops, while Gp12, Gp21 denote the transfer functions of two cross-coupling
loops, and Gc1, Gc2 represents the transfer functions of two outer loops. In addition, Ysp1,
Ysp2 are inputs of the outer-loops, while U1(s) and U2(s) are the outputs of the outer-loops
and the inputs to the inner loops and Y1, Y2 are the outputs of the MIMO system.
The closed-loop transfer function of the output Y1(s) to the input Ysp1(s) can be given
by,
Y1 ( s)
G p11
U1 ( s )
(4-11)
For simplicity purpose, all the (s) are ignored in the following analysis. Hence, the
following equations can be given.
Y1 U1 Gp11 U 2 Gp12
(4-12)
U 2 (Ysp 2 Y2 ) Gc 2
(4-13)
Y2 U1 G p 21 U 2 G p 22
Y2 U1 G p 21 (Ysp 2 Y2 ) Gc 2 G p 22
(4-14)
Y2 (1 Gc 2 Gp 22 ) U1 Gp 21 Ysp 2 Gc 2 Gp 22 ,
(4-15)
Y2
U1 G p 21
(4-16)
1 Gc 2 G p 22
U 2 Y2 Gc 2
Gc 2 G p 21
1 Gc 2 G p 22
U1
(4-17)
Y1 U1 G p11
Y1 (G p11
Gc 2 G p 21 G p12
1 Gc 2 G p 22
Gc 2 G p 21 G p12
1 Gc 2 G p 22
U1
) U1
(4-18)
(4-19)
Without considering the feedback loop in red in Figure 4-23, which means leaving the
outer loop appended, then the transfer function becomes,
145
Y1 Gp11 U1
(4-20)
By comparing equation (4-19) and equation (4-20), we can come to the conclusion that
that to involve the power controller in the model can achieve a more precise basement for
AC voltage controller design.
A. AC Voltage Controller Design
Using the methodology proposed in section 4.3.3.3 for power controller design, the
first step is to choose the baselines for controller design. Scenarios 3, 7, 9 and 36 (i.e.
PmaxQmax, P0Qmax, P0Qmin for SCR=2 with angle
Magnitude(dB)
10
10
10
Phase(deg)
-200
-400 0
10
Case3
Case7
Case9
Case36
.
1
10
10
Frequency (rad/sec)
10
Figure 4-24 Bode plots of selected cases in the AC voltage controller design
40
Magnitude(dB)
Magnitude (dB)
0
-10
-20
Filter: 1/(1+0.0075*s)
-30
Phase(deg)
Phase (deg)
-20
1
10
10
10
-45
10
10
10
Frequency (rad/s)
(a)
146
-40 0
10
-40
0
-90
0
10
20
10
-200
-400
-600 0
10
Case3
Case7
Case9
Case36
1
10
10
Frequency (rad/sec)
(b)
10
(c)
(d)
Figure 4-25 The Bode plots of (a) low-pass filter; (b) Uf(s)/Idref (s) of DB current controlled
VSC with a low-pass filter; (c) the designed PI compensator and (d) open loop transfer
function of Uf(s)/Ufref(s) at the rectifier side.
Similar to the power controller design, a first-order low-pass filter (see Figure 4-25a) is
added to the system to provide an additional attenuation that 20dB/decay which rolls off
when exceeding 133.33rad/s (see Figure 4-25b). Then a set of PI parameters Kp=
0.0314252 and Ki= 1.575 are chosen, which yield a phase margin of better than about 60.9
degree at a gain cross-over frequency of about 41.7 rad/s and a gain margin better than
about 9.6 dB (see Figure 4-25d). The closed-loop responses of Uf(s)/Ufref(s) (including all
36 scenarios) being subject to a 1pu step change in the AC voltage reference are employed
to evaluate the performance of the designed AC voltage controller as shown in Figure 4-26.
1.2
Case 3 Oershoot(%):9.91
At time(seconds):0.0761
Case 36
Settling time(seconds):0.336
Uf(pu)
0.8
0.6
0.4
0.2
0
-0.2
0
0.1
0.2
0.3
0.4
Time (s) (seconds)
0.5
0.6
147
U f m
U f ref +
U f_ f
1/(1+TfUf_rec s)
K P_Uf
+
K i_Uf /s
i dref (s)
U fX
U fref U fm U f _ f
d U fX
dt
T fU f _ rec
Ki _ U f U f _ f
idref K P _U f U f _ f U fX
(4-21)
(4-22)
(4-23)
148
(a)
(b)
(d)
(c)
Figure 4-28 The frequency responses of (a) Udc(s)/Iqref(s); (b) Udc(s)/Idref(s); (c) Uf_inv(s)/
Iq_ref_inv(s) and (d) Uf_inv(s)/Id_ref_inv(s) being subject to the operating point variation of the
rectifier side converter (36 scenarios), but with a constant operating point at the inverter
side converter that is PmaxQ0, SCR=2 and with angle equals to 90o.
Note that 36 operating points are considered in the controller design for the rectifier.
However, it is not feasible to consider 36 operating points for both of the inverter and
rectifier operating points together (362) when designing an inverter side controller.
Therefore, to reduce the number of the total operating points considered it is necessary to
examine the sensitivity of the transfer functions for a given frequent operating point of the
inverter first. In such an approach, if the sensitivity is less, any operating point can be
chosen as a base operating point for further controller design, normally the most frequent
operating point is chosen as the base scenario. In such a case, only 36 inverter side
operating points are considered, which significantly reduce the initial total number 362.
Figure 4-28 shows the frequency responses of the transfer functions Udc(s)/Iqref(s),
Udc(s)/Idref(s), Uf_inv(s)/ Iq_ref_inv(s) and Uf_inv(s)/ Id_ref_inv(s) which correspond to the 36
149
operating point variations of the rectifier. In this figure, a constant operating point of the
inverter (PmaxQ0, SCR=2, 90o angle) is assumed. Note that no significant discrepancies are
observed in the results of the two main loops up to 100rad/s (Figure 4-28a, d). However,
results are found to be different in the cross-coupling loops (Figure 4-28b, c). This
indicates that the outer loop controllers for the main loops need a careful design approach
to minimize the cross-coupling effects. Anyhow, it can be concluded that the sensitivity of
the transfer functions of operating points is low, in fact, of which the impact can be ignored
for simplification purpose. Therefore, based on this assumption, only a single frequent
operating point (PmaxQ0, SCR=2, 90o angle) is sufficient for the base scenario of the
rectifier.
It should be emphasized that the VSC losses are not considered in the small signal
model. However, the VSC losses cause inaccuracies in the steady state calculation. This
results in a minor power imbalance further causing minor inaccuracies in the filter bus
voltage. Consequently, the DC voltage also becomes inaccurate, since the calculated DC
voltage is utilized rather than the real DC transmission line feedback data in modeling of
the control system. This approach avoids the potential problems caused by the feedback
dynamics of the transmission line which can complicate the calculation process. Although
the above approach reduced the accuracy of the model, the deviation in the final results is
found small once the loops of the entire system are closed.
4.3.4.1 DC Voltage Controller Design for the Inverter and its Small Signal Model
Figure 4-29a shows the frequency responses of the selected cases in the DC voltage
controller design. By employing the compensator shown in Figure 4-29b, it can be
predicted that a satisfactory performance can be achieved when being subject to a step
change in the DC voltage reference (see Figure 4-29d) based on the analysis of the Bode
plots of the compensated system (see Figure 4-29c).
By testing the frequency response for all 36 scenarios, it is found that the system using
the designed controller yields a phase margin of better than 85.9 degree at a gain crossover frequency of about 9.71 rad/s and a gain margin of better than 24.3 dB (see Figure
4-29c). Furthermore, the results of the step response tests reveal that a fairly low maximum
overshoot 1.89% (at t=0.212s) is obtained in the scenario that PminQmin, SCR=2 and angle
=90o. In addition, a satisfactory maximum settling time at t=0.532 is obtained at the
scenario that PminQmax, SCR=7.5 and angle =90o.
150
(a)
(b)
(c)
(d)
Figure 4-29 (a) The Bode plots of the selected cases of the DC voltage controller design
Udc(s)/Iqref_inv(s); (b) designed PID compensator (Kp_dc+Ki_dc/s+Kd_dcs/(1+TfUdcs)=-0.0240.088/s+0.0009s/(1+0.06s)); (c) frequency responses of open loop transfer functions of
Udc(s)/ Udc_ref (s) and (d) step responses of the closed-loop Udc(s)/ Udc_ref(s).
The control block diagram of the DC voltage controller is shown in Figure 4-30 and its
small signal model is given by equation(4-24).
U dcm
U dcref
+-
K P _ U dc
K i _ U dc / s
K d _ U dc s
U dcX 1
iqref ( s )
U dcX 2
1 T fU dc s
151
GPID ( s) K p Ki / s K d s / (1 T f s )
GPID ( s)
( K p T f K d ) s 2 ( K p Ki T f ) s Ki
(4-24)
Tf s2 s
For facilitating the derivation of the state space equations, it can be written as a secondorder transfer function as given in equation (4-24).
If we substitute the values A= K p Tf +K d ; B= K p +Ki Tf ; C=Ki; D=Tf; E=1and F=0
into equation (4-24), the general form of the second order form transfer function can be
given by
GPID ( s)
A s2 B s C
D s2 E s F
(4-25)
The state space equation for a general form second order transfer function can be given
by
F
AF
0
C
w
w
d d1
D d1
D
ud
E wd 2 AE
dt wd 2
1
B
D
1 wd 1 A
yd 0
ud
D wd 2 D
(4-26)
Hence, if by i) substituting two states of the general form of the second order transfer
function w d1 and w d2 with UdcX1 and UdcX2 ; ii) relacing u d with Udcref -Udcm and iii)
substituting y d with iqref (s) , the small signal for the DC-PID controller can be obtained as
F
AF
AF
0
C
C
d U dcX 1
D U dcX 1
D
D U dcref
U dcm
E U dcX 2 AE
AE
dt U dcX 2
1
B
B
D
D
1 U dcX 1 A
A U dcref
iqref 0
D U dcX 2 D
D U dcm
(4-27)
4.3.4.2 AC Voltage Controller Design for Inverter and its Small Signal Model
Note that the model used in the AC voltage controller design is a MIMO model by
placing the designed rectifier side power controller, AC voltage controller and the inverter
side DC voltage controller in service. This approach ensures the design is more accurate
152
since the impacts of the dynamics of the co-operating controllers are also being taken into
account.
(a)
(b)
(c)
(d)
Figure 4-31(a) Bode plots of selected cases for inverter end AC voltage controller design
Uf_inv(s)/Idref_inv(s); (b) designed filter+PI compensator; (c) frequency responses of open
loop transfer functions of Uf_inv(s)/Uf_inv_ref(s); (d) step responses of closed-loop Uf_inv(s)/
Uf_inv_ref(s).
The same methodology is used in this controller design as proposed earlier in section
4.3.3.3 and shown in Figure 4-31a. Using the loop-shaping technique, a set of PID
parameters (Kp_Uf_inv= 0.00842, Ki_Uf_inv=0.752, Kd_Uf_inv=-8.910-5 and Tf_Uf_inv=0.018) are
selected for the PID compensator. The frequency responses of the designed AC voltage
153
controller is shown in Figure 4-31b, and the compensated system is shown in Figure 4-31c.
Note that, the AC voltage loop yields a phase margin of better than 63.2o at a gain crossover frequency of 27.3 rad/s and a gain margin better than 17.8 dB.
The small signal model with respect to the AC voltage is similar to the method
explained earlier in section 4.3.4.1. Figure 4-31d shows the step test responses of the
system operating under 36 different scenarios following a 1pu step change in the reference
value of the filter bus voltage. A satisfactory performance is obtained with a maximum rise
time of 0.404s (reaching to 90% of the final value). In addition, a settling time of better
than 0.73s is obtained in the STATCOM mode in which VSC is embedded in a strong
system which is represented as an inductance in series with a resistance. Moreover, a
maximum overshoot of 7.02% is also obtained within the satisfactory range.
154
which half of the DC voltage step change occurred within a short time 0.5s. It can be
concluded that the cross-coupling effects between the power and the AC voltage induced
by changing the power reference are low. However, the cross-coupling effects between the
rectifier side AC voltage and the power or DC voltage caused by changing the voltage
reference are not insignificant, which is noticeable.
(a)
(b)
(c)
(d)
Figure 4-32 1% step of AC voltage reference at the rectifier end; (a) responses of power at
the rectifier end; (b) responses of the filter bus voltage Uf_rec at the rectifier end; (c)
responses of DC voltage Udc_inv at the inverter end; (d) responses of the filter bus voltage
Uf_inv at the inverter end.
155
(a)
(b)
(c)
(d)
Figure 4-33 1% step of power reference at the rectifier end; (a) responses of the power at
the rectifier end ; (b) responses of the filter bus voltage Uf_rec at the rectifier end; (c)
responses of the DC voltage Udc_inv at the inverter end; (d) responses of the filter bus
voltage Uf_inv at the inverter end.
156
coupling. Although it is not the case in a DC meshed power grid, this is beyond the scope
of this research. The peak of the DC voltage transient is found to be about 5 times less than
the AC voltage change, and the rising time is less than 0.05s. In addition, it settles down to
half of the change in less than 0.5s, which are all satisfactory.
(a)
(c)
(b)
(d)
Figure 4-34 1% Step of the AC voltage reference at the inverter end; (a) response of the
power at the rectifier end; (b) response of filter bus voltage Uf_rec at the rectifier end; (c)
response of the DC voltage Udc_inv at the inverter end and (d) response of the filter bus
voltage Uf_inv at the inverter end.
157
(a)
(b)
(c)
(d)
Figure 4-35 1% step of the DC reference voltage at the rectifier side; (a) the responses of
power at the rectifier; (b) the responses of the filter bus voltage Uf_rec at the rectifier; (c) the
responses of the DC voltage Udc_inv at the inverter end; (d) the responses of the filter bus
voltage Uf_inv at the inverter end.
By observing Figure 4-35d, it is found that the 1% step change in the DC voltage
reference results in about 0.196% change in the AC voltage which is not quite satisfactory
but acceptable. However, the cross-coupling effects can be improved by further slowing
down the DC voltage controller through multiplying the controller coefficient with a gain
less than 1 (e.g. 0.8), but at cost of reducing the response speed of the DC voltage
controller. Moreover, further reducing the speed of the AC voltage controller can also
provide another alternative solution to reduce the cross-coupling effects.
From the above analysis, it can be concluded that the overall performance of the
system controllers behaviour is a trade-off between the primary controller performance
158
and the cross-coupling effects. In other words, a better cross-coupling effect can be
achieved by slowing down the primary controller, meanwhile the performance of the
primary controllers have to be ensured as well. This requires further fine tuning of the
controller parameters by iterating the controller design procedure to get the final desirable
control system. Also, the root locus analysis and PSCAD simulation verifications can also
assist in tuning such controllers.
4.3.5.3 Analytical Model Verification against Detailed Nonlinear Simulation
It is well known that the small signal models are only accurate around a certain
operating points. Therefore the analytical model verification against detailed nonlinear
simulation (such as PSCAD) is always required.
Since 36 operating scenarios are involved in this research, the bunch file simulations
are required. However, PSCAD is only applicable for a single case simulation. The Multirun block of PSCAD is also not applicable for this particular application. Hence, a selfdefined multi-study component needs to be developed to successively read the input data
from a text file. The procedures for the PSCAD verification is given below,
Compute grid operating points and generate a text file containing these data points
in a form that is suitable as an input to PSCAD.
Create a multi-study component that will read a file from an input file, and read
the operating-point values of the system successively from the above file, and
modify grid impedance, and set control input references accordingly.
Proceed the simulation results in Matlab for comparison of PSCAD results with the
small signal model at various step sizes, operating points and SCRs.
Performing the step tests for 36 operating points and with different step sizes, it is
found that the developed small single model is accurate enough to capture the main
characteristics of the system model developed in PSCAD. It was also observed that the
developed small signal model can be used as a powerful tool for controller design as well
as in the small signal stability analysis of the system. However, due to the space limitations,
for verification purpose, the results are only given for two operating points of the inverter: i)
PmaxQmax, SCR=2 with =90o and ii) P0Qmin, SCR=7.5 with =75o for various step size
changes of 1%, 2%, 4% and 8%.
159
Figure 4-36 shows the closed-loop responses of the complete control system, under the
step size changes of 1%, 2%, 4% and 8% of the reference value of the inverter side filter
bus voltage (Figure 4-36 a) and the DC voltage (Figure 4-36b). As it can be seen in the
figure the inverter side filter bus voltage results obtained from both of the small signal and
the non-linear model match very well with each other. Although some minor errors are
observed in the cross-coupling loops at the rectifier side, the small signal model does
capture the main characteristic of the system. It can be deduced that the errors may be due
to the inaccuracy introduced during the calculation process of the steady-state value, since
the losses on both side of the converters are not included. However, the level of error is
acceptable and cannot significantly affect the controller design.
Udcinv(%)
0.5
-0.02
0
0.5
I cdrec,I cqrec(%)
0.5
0.5
0
-0.5
0
0.5
Time(Seconds)
1
(a)
160
-0.5
20
Ufinv(%)
Ufrec(%)
-0.2
0
0.02
I cdinv,I cqinv(%)
P(%)
0.2
0.5
0.5
0.5
Time(Seconds)
0
-2
0
2
0
-2
0
0.5
0
-0.1
0
1
0.5
0
-1
0
0
-2
0
0.5
Ufinv(%)
Ufrec(%)
Udcinv(%)
0
-0.5
0
0.1
I cdrec,I cqrec(%)
I cdinv,I cqinv(%)
P(%)
0.5
0.5
1
Time(Seconds)
0.5
0.5
0.5
Time(Seconds)
0
-0.5
0
5
0
-5
0
(b)
Figure 4-36 The step responses of Ufinv_ref (a) and Udcinv_ref (b) and their induced
performances in their cross-coupling loops as the system working at an operating point
PmaxQmax, SCR=2 with =90o. red: Matlab 1%; cyan: PSCAD 1%; black: PSCAD 2%/2;
magenta: PSCAD 4%/4; green: PSCAD 8%/8.
Figure 4-37 is another step test performed in the STATCOM mode. Firstly, a 1% step
change is applied to the small signal model of the inverter side DC voltage reference,
which is then compared with four results from the detailed PSCAD non-linear model. The
testes set on PSCAD include 1%, 2%, 4% and 8% step changes on both of the AC filter
bus reference voltage and the DC reference voltage. While, the final output results are
manipulated by dividing 1, 2, 4 and 8 respectively. To do in such a way is to examine the
linearity of the PSCAD model. The results especially forcefully confirm the accuracy of
the developed small signal model.
161
Udcinv(%)
P(%)
0.5
0
-0.5
0
0.5
0
-2
0
Ufinv(%)
0
-0.02
0
I cdrec,I cqrec(%)
0.5
0.5
Time(Seconds)
0.1
0.5
0
-1
0
0.5
Time(Seconds)
-0.1
0
I cdinv,I cqinv(%)
Ufrec(%)
0.02
0.5
2
0
-2
0
(a)
Udcinv(%)
x 10
0.5
0.5
0.5
0
-1
0
0.05
0
-0.05
0
-0.05
0
-3
0
-1
0
I cdrec,I cqrec(%)
0.5
Ufinv(%)
Ufrec(%)
-0.01
0
0.05
I cdinv,I cqinv(%)
P(%)
0.01
0.5
1
Time(Seconds)
10
0
-10
0
0.5
1
Time(Seconds)
(b)
Figure 4-37 The step responses of Udcinv_ref (a) and Ufinv_ref (b) and their induced
performances in their cross-coupling loops as the system working at an operating point
P0Qmin, SCR=7.5 with =75o. red: Matlab 1%; cyan: PSCAD 1%; black: PSCAD 2%/2;
magenta: PSCAD 4%/4; green: PSCAD 8%/8.
162
163
controller with the previous designed controller in service, since the more detailed model
employed will result in a more accurate design for the later controller. In addition, such
approach also allows the cross-coupling effects to be considered.
Finally, a self-defined block is developed in PSCAD to enable the successive PSCAD
simulations for verification purpose. The co-operating performance of the various designed
controllers are examined against the time domain verification. It is found that the designed
power and voltages controllers work very well within all the considered 36 operating
points.
164
5.1 Introduction
As is known, the limited availability and environmental concern of fossil fuels, as well
as the continuous growing demand of electricity, have caused renewable energy to become
commercially attractive. Meanwhile, as a result of increasing greenhouse gas emissions,
the Australian government has announced an emissions trading scheme in 2010 termed as
the Carbon Pollution Reduction Scheme [142] and setup new green energy targets to
165
the
longest interconnected
power
Douglas
of
Queensland and Port Lincoln of South Australia with an end-to-end distance of more than
4000 km as shown in Appendix E [144]. In addition, the potential Paralana Geothermal
source is located very remote from the existing shared transmission network. There are
several investigations [145-147] on exploring the potential network extension solution
upon integrating the staged Innamincka generation of 500MW, 2000MW and 5000MW
[145, 146].
This chapter examines the small signal stability performance of this potential extended
Australian grid using an alternative HVDC transmission solution assuming that a 2000MW
capacity is available at Innamincka. It is assumed that the generation 200MW is consumed
locally, and the remaining 1200MW is delivered via a 500kV, 1100km, bi-pole VSCHVDC links to Armadale converter substation located at the half way between Queensland
and New South Wales. It is considered that this arrangement can make full utilization of
the existing transmission corridor, consequently benefiting both States. The remaining 600
MW are transferred to Roxby Down from Innamincka, using the 490km 275kV
transmission line, and then being integrated to the main grid at Port Augusta by a 290km
transmission line. The proposed diagram of the extended South-East Australian power
system with VSC-HVDC links is illustrated in Figure 5-1, which is obtained based on an
IEEE 14 generator 59 bus test system. The test system represents a simplified model of the
Southern and Eastern Australian network which is composed of five areas in which 14
large generators and 5 static var compensators (SVCs) are involved. The data for this
system is available in [148]. Comparing with the other possible solutions, it is considered
that this approach is a realistic option and will be investigated here.
166
406
SVC
SPS_4
GPS_4
N_SVC
AREA4:QSL
407
403
404
405
CPS_4
408
402
409
TPS_4
606
601
602
603
410
AREA6:Innamincka
Brisbane
401
VAI Innamincka
604
VAB
411
413
201
414
412
415
ASVC_2
PL34
SVC
SVC BSVC_4
416
VPS_2
BPS_2
205
201
207
203
PL12
206
609
208
EPS_2
202
NPS_3
609 608
209
211
MPS_2
501
TPS_3
204
504
210
Sydney
502
VAS
215
212
505
213
PSVC_5
AREA2:NSW
214
SVC
AREA5:SA
507
Adelaide
216
508 VAA
506
217
503
PPS_3
509
SVC
102
SSVC_5
305
315
311
310
AREA1:Canberra
309
LPS_3
304
308
HPS_1
307
301
303
306
AREA3:VIC
YPS_3
RSVC_3
314
302
SVC
313 VAM Melbourne
312
Figure 5-1 The diagram of the extended Simplified South-East Australian power grid with
VSC-HVDC links.
167
u gfI
u gfR
VcRg
g
I cR
REACTANCE
(Dynamic)
VSC
VcIg
g
I gR
FILTER
(Dynamic)
I cIg
I gIg
Grid
(Adimittance)
168
di g g _ R
(1/ Lg ) E g g _ R (1/ Lg )u g f _ R ( Rg / Lg )i g g _ R 0i g g _ I
dt
g
di g _ I
(1/ Lg ) E g g _ I (1/ Lg )u g f _ I ( Rg / Lg )i g g _ I 0i g g _ R
dt
(5-1)
which can be replaced with an grid admittance model and resulting in the following set
of equations,
i gg
i gg
( E g g _ R jE g g _ I ) (u g f _ R ju g f _ I )
Rg jX g
( E g g _ R u g f _ R ) Rg ( E g g _ I u g f _ I ) X g
Rg2 X g2
( E g g _ R u g f _ R ) X g ( E g g _ I u g f _ I ) Rg
(5-2)
Rg2 X g2
In the above equations (5-2), the real part of the grid current i ggR and the imaginary part
of the grid current i ggI can be expressed by the following set of equations,
i
g
gR
i
g
gI
( E g g _ R u g f _ R ) Rg ( E g g _ I u g f _ I ) X g
Rg2 X g2
( E g g _ R u g f _ R ) X g ( E g g _ I u g f _ I ) Rg
(5-3)
Rg2 X g2
After linearizing the above equations (5-3), the small signal model for the admittance
representation of the grid model can be given by
i g gR
Rg
R X
i g gI
2
g
2
g
(E g g _ R u g f _ R )
Xg
Rg2 X g2
Xg
R X
(E g g _ R u g f _ R )
2
g
2
g
(E g g _ I u g f _ I )
Rg
Rg2 X g2
(5-4)
(E g g _ I u g f _ I )
Figure 5-3 shows the simulation results of the dynamic and the admittance
representation of the grid model following a 1% step change on the inverter side input DC
voltage reference. As it can be seen in the figure, the responses of the DC voltage tracking
and the cross-coupled rectifier side power responses are similar. However, the crosscoupling voltage responses in the rectifier and the inverter side have some minor
differences which can be ignored, as highlighted by the circled areas in the graphs.
169
AC vltage-rec(%)
Power-rec(%)
0.5
Grid Dynamic
Grid Admittance only
0
-0.5
0
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
0.7
0.05
-0.05
0
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
0.5
(%)
0.9
Grid Dynamic
Grid Admittance only
Ic
dq
0.8
-0.5
0
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
0.7
0.8
0.9
AC vltage(%)-inv
DC voltage-inv(%)
(a)
2
Grid Dynamic
Grid Admittance only
0
-2
0
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
0.7
0.2
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
Ic
dq
(%)
-2
0
0.9
Grid Dynamic
Grid Admittance only
0
-0.2
0
0.8
0.1
0.2
0.3
0.4
0.5
0.6
Time(Seconds)
0.7
0.8
0.9
(b)
Figure 5-3 1% DC voltage step responses of the dynamic and grid admittance models (a)
for rectifier side; (b) for inverter side.
170
+
+
g
I gR
g
EgR
EgIg
g
gI
Grid
(Admittance)
u gfR
+
+ g
I fI
g
fR
Filter
(Admittance)
I cRg
I cIg
VSC
VcRg
VcIg
u gfI
Finally, obtain the injected currents to the grid source I ggR and I ggI according to the second
equation of the set of equation (5-5).
u fg Egg Z g I gg
g
g
g
I g I f Ic
I g ug B j
f
f
f
(5-5)
If this set of equations (5-5) is expended and linearized, the small signal model of the
filter and grid admittance representation of the system model can be obtained as,
u g f _ R E g g _ R Rg I g g _ R X g I g g _ I
g
g
g
g
u f _ I E g _ I Rg I g _ I X g I g _ R
(5-6)
g
g
g
I g _ R I f _ R I c _ R
g
g
g
I g _ I I f _ I I c _ I
(5-7)
g
g
I f _ R u f _ I B f
g
g
I f _ I u f _ R B f
(5-8)
171
Power-rec(%)
0.4
0.2
0
-0.2
-0.4
0.2
0.4
0.6
Time(Seconds)
0.8
AC vltage-rec(%)
0.06
Grid Admittance only
Grid,Filter Admittance
0.04
0.02
0
-0.02
0.2
0.4
0.4
0.6
Time(Seconds)
0.8
Ic (%)
0.2
dq
0
Grid Admittance Only Icd
Grid Admittance Only Icq
Grid and filter Admittance Icd
Grid and filter Admittance Icq
-0.2
-0.4
0.2
0.4
0.6
Time(Seconds)
0.8
DC voltage-inv(%)
(a)
1.5
1
0.5
0
-0.5
0.2
0.4
0.6
Time(Seconds)
0.8
AC vltage(%)-inv
0.15
0.1
0.05
0
Grid Admittance only
Grid,Filter Admittance
-0.05
-0.1
0.2
0.5
0.4
0.6
Time(Seconds)
0.8
dq
Ic (%)
0
Grid Admittance Only Icd-inv
Grid Admittance Only Icq-inv
Grid and filter Admittance Icd-inv
Grid and filter Admittance Icq-inv
-0.5
-1
-1.5
0.2
0.4
0.6
Time(Seconds)
0.8
(b)
Figure 5-5 1% DC voltage step response comparison of the admittance representation of
grid model only, and the gird and filter models adopting admittance representation: (a) the
rectifier side; (b) the inverter side.
172
The test results presented in Figure 5-5 reveal that the further simplification of the filter
model does not weaken the accuracy of the results in this application. Hence, the
admittance representation of both the grid and the filter models will be set as the standard
modelling methodology in the following analysis.
x ( new)
x ( original ) ( new)
xbase
( original )
xbase
Model Parameters
(Original)
(5-9)
New Rating
Parameter Conversion
Model Parameters
(New in SI units)
Model Equations
in SI units
173
I ( SI )
I (P)
V (P)
N1( p )
I base
Vbase
V ( SI )
I base1
N1( sys1) / Zbase1
Vbase1
(5-10)
Ibase 2
N1( sys 2) / Zbase 2
Vbase 2
(5-11)
Zbase 2
Zbase1
(5-12)
(1)
R (1)
N1 1 N (1) ;
(2)
N (2) R
;
1
1 N 2(2)
Set T
(1)
Set T
(2)
R (1)
( (1) );
L
R (2)
( (2) );
L
(5-13)
Then, calculate the per unit value of coefficients R(p), L(p) and T (p) of the systems.
174
( P ) R (1)
L(1)
; L( P )
R
Rb1
Lb1
(1)
(1)
T ( P ) R / Rb1 R Lb1
(5-14)
Then calculate the updated coefficients R(2) and L(2) of the new system
(2)
( P)
(1) Rb 2
R R Rb 2 R R ;
b1
L(2) L( P ) L L(1) Lb 2 ;
b2
Lb1
(5-15)
Finally, It can be obtained that the coefficient N2 for both systems are the same by
calculating the time constant T of N2 and further to calculate N2,
Rb 2
Rb1
R
R (1)
( (2) )
( (1) ) T (1)
L
L
L
L(1) b 2
Lb1
R (1)
(2)
T (2)
(5-16)
(1)
Therefore, it can be derived that N(2)
2 =N 2 . Then using the calculated value of N2, we
can obtain.
N1(2)
(2)
R
1 N 2(2)
N1(2) N1(1)
Rb 2
)
Rb1
;
1 N 2(2)
R (1) (
Rb 2
Z
N1( sys1) base 2
Rb1
Zbase1
(5-17)
(5-18)
From the above analysis, it can be concluded that the DB controller coefficient directly
calculated using the equation (2-48) should be the same as the scaled value using the
original system.
5.2.2.2 PLL Controller Coefficient
The PLL controller coefficient can be modified by multiplying by Vb, as given in the
below equation (5-19).
Y /U
Y /U
SI
pu
/ VSI
/ b
VSI / Vb
Vb / VSI
(5-19)
175
K p(2) K p(1) (
Pbase1 I base 2
)(
)
Pbase 2
I base1
(5-20)
Ki(2) Ki(1) (
Pbase1 I base 2
)(
)
Pbase 2
I base1
(5-21)
Kp
Pe( MW ) / U (f kV )
1
1 Tf s
+
+
Ki
s
( kA )
( kA )
I qref
/ I dref
Kp
U dc( kV )
Ki
s
Kd s
1 Tf s
( kA )
I qref
176
Integral part
Udc Ki/s=Iqref
Derivative part :
Iqref/ dt)
Iqref
Therefore, according to Table 5-1, the new PID coefficients for the DC voltage loop
can be given by
K p(2) K p(1) (
Vdc _ base1
Vdc _ base 2
Ki(2) Ki(1) (
K d(2) K d(1) (
Vdc _ base1
Vdc _ base 2
Vdc _ base1
Vdc _ base 2
)(
I base 2
)
I base1
(5-22)
)(
I base 2
)
I base1
(5-23)
)(
I base 2
)
I base1
(5-24)
177
250
200
Im
150
100
50
0
-300
-250
-200
-150
Re
-100
-50
178
in
[16],
the
highest
pole
voltage
in
recorded
operation
was
(Skagerrak 4, Denmark Norway). Therefore, the selection of the DC link parameters in this
research is considered realistic.
Figure 5-12 shows the eigenvalue map of the system with a new series of updated DC
link parameters. It can be seen from Figure 5-12 that there is significant changes in a pair
of modes: 41.83 rad/s moves to 84.91rad/s and 184 rad/s moves towards 569.2 rad/s.
Therefore, it is important to determine whether this change will affect the efficiency of the
designed controllers.
To investigate this issue, a review on performance evaluation in terms of all the power
and voltage controllers with a new series of DC link parameters are also undertaken, and
the simulation results are given in Appendix F, which are then being compared with the
results in Chapter 4. It was concluded that no significant effects are observed on the
controllers within the frequency range considered. As it was also observed in the power
and the DC voltage control loops, the pattern of the frequency responses is much more
consistent due to the changes in the eigenvalues. It can be noted that this change decreases
the difficulties of controller design for these VSC-HVDC links embedded in a weak AC
system. Hence it can be concluded that the robustness and sensitivity to the DC link
parameters of the designed controllers are acceptable.
600
400
Im
200
-200
-400
-600
-100
-80
-60
-40
-20
Re
Figure 5-12 The comparison of the system eigenvalues between the new updated DC link
system and the original system
179
Vr(s)/Ir(s), Vi(s)/Ir(s), Vr(s)/Ii(s) and Vi(s)/Ii(s) for both of the rectifier side at Innamincka
and the inverter side at Armadale. Note that if the source is represented by a constant
voltage behind source impedance (Zs=Rs+jXs), then the above transfer functions will have
the following four forms: Vr(s)/Ir(s) =Rs, Vi(s)/Ir(s) =Xs, Vr(s)/Ii(s) =-Xs and Vi(s)/Ii(s) =Rs.
This confirms the results given in Figure 5-13a (Rs) to be the high frequency responses of
the system analyzed. Further analysis of Figure 5-13a also demonstrates the resistance and
inductance values at the inverter side,
(5-25)
base 100 MW
(5-26)
base 100 MW
(5-27)
base 100 MW
(5-28)
base 100 MW
However, at very low and at intermediate frequencies, the source cannot be represented
by a voltage behind impedance. For example, at very low frequencies at the inverter side,
the results of the transfer functions are,
180
base 100 MW
base 100 MW
(5-29)
(5-30)
base 100 MW
base 100 MW
(5-31)
(5-32)
As shown above Vr (s)/Ir (s) is not equal to Vi (s)/Ii (s) , and Vi (s)/Ir (s) is not equal to
-Vr (s)/Ii (s) . This indicates that the VSC controllers adopted in chapter 4 are in fact to
(a)
(b)
Figure 5-13 The frequency responses for the rectifier side (a) and for inverter side (b) for
the higher order grid impedance models (I to V).
181
Magnitude(dB)
200
100
0
-100
-200
-300 -4
10
-2
10
10
10
10
10
Phase(deg)
-100
-200
-300
-400
-500 -4
10
Grid simplified
Grid detailed
-2
10
10
10
Frequency (rad/sec)
10
10
182
0.15
Inverter side DC voltage(%)
Simplified
Detailed
0.5
-0.5
0
0.01
0
-0.01
-0.02
-0.03
0
2
3
Time(Seconds)
Simplified
Detailed
4
5
1.5
Simplified
Detailed
0.1
0.05
-0.05
0
0.01
0
-0.01
-0.02
-0.03
-0.04
0
2
3
Time(Seconds)
Simplified
Detailed
4
5
Figure 5-15 The step responses test of the power reference at the rectifier side
Magnitude(dB)
200
0
-200
-400 -4
10
-2
10
10
10
10
10
Phase(deg)
0
-200
-400
-600 -4
10
Grid simplified
Grid detailed
-2
10
10
10
Frequency (rad/sec)
10
10
183
Simplified
Detailed
0.6
0.4
0.2
0
-0.2
0
1.5
5
Inverter side AC voltage(%)
0.8
Simplified
Detailed
1
0.5
0
-0.5
0
2
3
Time(Seconds)
0.15
Simplified
Detailed
0.1
0.05
0
-0.05
0
0.04
Simplified
Detailed
0.02
0
-0.02
-0.04
0
2
3
Time(Seconds)
Figure 5-17 The step response test of the AC voltage reference at the rectifier side
Magnitude(dB)
200
100
0
-100
-200 -4
10
-2
10
10
10
10
10
Phase(deg)
0
-100
-200
-300
-400 -4
10
Grid simplified
Grid detailed
-2
10
10
10
Frequency (rad/sec)
10
10
184
Simplified
Detailed
0.1
0.05
0
-0.05
0
0.5
0.01
Simplified
Detailed
0.005
0
-0.005
-0.01
0
0.5
Time(Seconds)
0.05
Simplified
Detailed
0
-0.05
-0.1
-0.15
0
0.15
0.5
1.5
Simplified
Detailed
1
0.5
0
-0.5
0
0.5
Time(Seconds)
Figure 5-19 The step response test of the DC voltage reference at the inverter side
Magnitude(dB)
200
100
0
-100
-200
-300 -4
10
-2
10
10
10
10
10
Phase(deg)
-100
-200
-300
-400
-500 -4
10
Grid simplified
Grid detailed
-2
10
10
10
Frequency (rad/sec)
10
10
185
0.1
0.05
-0.05
0
0.01
Rectifier side AC voltage (%)
0.05
Simplified
Detailed
0.005
-0.005
2
3
Time(Seconds)
-0.05
-0.1
-0.15
0
1.5
Simplified
Detailed
-0.01
0
Simplified
Detailed
0.15
Simplified
Detailed
0.5
-0.5
0
2
3
Time(Seconds)
Figure 5-21 The step response of the AC voltage reference at the inverter side
Real
Imag
Damping Ratio
Mode Description
I 40
-0.286
3.907
0.073
I 35
-0.290
3.417
0.085
Innamincka VS VIC& SA
I 25
-0.242
2.273
0.106
I 15
-0.154
1.703
0.090
The eigenvalue map of the main electro mechanical modes of the interconnected
system is given in Figure 5-22, which is obtained by using Matlab m-file. The main inter-
186
area modes of I40, I35, I25, and I15 are found to be lightly damped in the figure. The
damping criterion used here is chosen according to the Australian National Electricity
Rules which require the halving time of the least damped electro-mechanical mode of
oscillation is not more than five seconds (1/T ln(2)/5) [144], which corresponds to a
damping constant with a magnitude greater than -0.14 Nepers/Second (Np/s).
11
Damping Criterion
10
-0.14Np/s
9
8
Im
7
6
5
4
I40
I35
3
I25
2
1
-1.5
I15
-1
-0.5
Re
187
VAR
W
.ING3
W
.ING1
W
.ING2
0.472 -164.51 W
0.432 -162.87 W
0.128 -147.06 W
.YPS_3
.LPS_3
.HPS_1
CASE
(A)
(A)
(A)
CASE
(A)
(A)
(A)
VAR
W
.NPS_5
W
.TPS_5
W
.PPS_5
Mag
Ph()
0.835
9.34
0.815
4.18
1.000
0.00
(A)
(A)
(A)
(A)-0.286 +3.907 i
p
2H
(5-33)
where, DPSS represents damping gain of the PSS on the machine rating, p is the
magnitude of the participation of the machine rotor speed in the selected electromechanical mode and H is the inertia constant of the machine on the machine rating.
188
0.0553
) 0.0259
2 3.2
(5-34)
State Variable
0.02
0.04
0.06 0.08
0.1
0.12
Participation factor
0.14
0.16
0.056
) -0.026
2 3.2
(5-35)
189
VAR
W
.LPS_3
W
.YPS_3
W
.HPS_1
0.467
0.404
0.468
0.122
0.172
0.166
W
W
W
W
W
W
-138.14
-136.55
-123.08
-110.61
-100.43
-95.23
.PPS_5
.TPS_5
.NPS_5
.MPS_2
.BPS_2
.VPS_2
CASE
(A)
(A)
(A)
CASE
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
VAR
W
.ING2
W
.ING1
W
.ING3
(A) W
.EPS_2
Mag
Ph()
1.000
0.00
1.000
0.00
1.000
0.00
0.169
-89.07
(A)-0.290 +3.417 i
State Variable
0.05
0.1
0.15
0.2
0.25
Participation factor
0.3
190
0.35
0.334
0.300
0.339
0.370
0.341
-179.40
-178.95
-174.12
-166.95
-164.17
W
W
W
W
W
.MPS_2
.VPS_2
.HPS_1
.YPS_3
.LPS_3
CASE
(A)
(A)
CASE VAR
(A) W
.ING2
(A) W
.ING3
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
(A)
W
W
W
W
W
W
W
W
.ING1
.NPS_5
.CPS_4
.TPS_5
.PPS_5
.TPS_4
.GPS_4
.SPS_4
Mag
Ph()
1.000
0.00
1.000
0.00
1.000
0.617
0.628
0.515
0.569
0.438
0.658
0.662
-0.00
-17.57
-19.39
-20.59
-20.96
-22.98
-23.52
-24.44
(A)-0.242 +2.273 i
0.0205
) -0.0096
2 3.2
(5-36)
191
State Variable
0.02
0.04
0.06 0.08
0.1
0.12
Participation factor
0.14
0.16
3 (
192
0.0671
) -0.031453
2 3.2
(5-37)
VAR
W
W
W
W
.SPS_4
.GPS_4
.TPS_4
.CPS_4
CASE
(A)
(A)
(A)
(A)
CASE VAR
(A) W
.ING3
(A)
(A)
(A)
(A)
(A)
(A)
(A)
W
W
W
W
W
W
W
.ING1
.ING2
.NPS_5
.PPS_5
.TPS_5
.LPS_3
.YPS_3
Mag
Ph()
1.000
0.00
1.000
1.000
0.777
0.776
0.729
0.224
0.221
-0.00
-0.00
-6.24
-8.26
-8.37
-27.44
-31.87
(A)-0.154 +1.703 i
State Variable
0.02
0.04
0.06
0.08
Participation factor
0.1
0.12
193
0p.u. A
5p.u. B
10p.u. C
15p.u. D
20p.u. E
25p.u. F
30p.u. G
10
9
8
Im
7
6
5
4
I40
I35
I25
I15
2
1
-1.5
-1
-0.5
Re
Figure 5-31 Eigenvalue evolution for the PSS damping gain on each generator is increased
from zero (no PSS in service) to 30 pu on machine base (750MVA) in 5 pu steps.
Table 5-3 Approximate improvements on system damping: comparison of the results
obtained from adding with equivalent damping torque (Figure 5-31) and the analysis for
participation factor.
Participation factor Analysis
Mode
Figure 5-31
(Approximation based on 30pu)
I40
-0.3699-(-0.2862)= -0.0837
I35
-0.3771-(-0.2899)= -0.0872
I25
-0.2683-(-0.2425)= -0.0158
I15
-0.2638-(-0.1542)= -0.1096
Table 5-3 gives the comparison results of the improvements on the system damping. In
the table, the third column data is calculated by the participation factor, and the data in the
second column is computed by using the system model by adding an equivalent damping
194
gain to the torque. In this study, the base capacity of the system (Sbase) is 100MVA, while
the machine base (Mbase) is set to 750MVA. It is shown that both of these cases have
similar results despite an error of 0.02 occurred in each mode. This also implies that the
approximation based on the participation factor is an effective way to evaluate the damping
performance of PSS.
5.3.3.1 PSS Design and Assessment
The approach for PSS design employed in this study is explained in detail in [99, 152].
As stated in [99],
The objective of PSS design is to induce on the shaft of the generator a torque of
electromagnetic origin which is in-phase with rotor-speed perturbations of the generator
Here, the procedure to design the PSS is given below:
Step 1: Determine the torque coefficients induced by PSS. In other words, specify the
DAMP(s) coefficient, add wash out filter and low pass filter to the coefficient making the
transfer function proper to obtain the transfer function of the stabilizer DAMP(s),
Step 2: Obtain the gain and phase characteristic of the excitation system, the generator
and the power system (i.e. the PVr characteristic),
Step 3: Apply curve fitting to the selected PVr characteristic,
Step 4: Derive the PSS transfer function,
Step5: Assess the performance of PSS by using eigenvalue analysis and by comparing
the torque coefficients induced by the PSS specified with the specified value.
A. Specifying the DAMP(s)
It was also stated in [99] that The wash-out and low-pass filter time constant are
chosen so that, over the frequency range of interest, DAMP(s) has near constant gain
(
) and phase between 0 and -30 degrees (depend on the damping requirements of
where
DAMP(s) De B(s)
(5-38)
sTw
1
B( s ) (
)(
)N
1 sTw 1 sTp
(5-39)
195
Bode Diagram
Magnitude (dB)
30
20
10
0
Phase (deg)
-10
90
45
0
-45
-90
-2
10
-1
10
10
10
Frequency (rad/s)
10
10
Figure 5-32 Bode plot for DAMP(s) with DPSS=20pu on Mbase 750 MVA, Tw=3s and
Tp=0.05s.
PVr (s)
18.6s 3268
0.0056916s 1
54.805
s 15.25s 59.63
0.01677s 2 0.25574s 1
2
(5-40)
Then converting PVr(s) from pu on Sbase (100 MVA) to pu on Mbase (750 MVA) and
specify the DAMP(s) function,
PVr _ M ( s) K sm G( s) K sm
1 0.0056916s
0.01677s 2 0.25574s 1
196
(5-41)
Magnitude(dB)
40
20
0
-20 -1
10
10
10
10
Phase(deg)
0
-100
-200
-300 -1
10
10
10
10
Frequency (rad/sec)
Figure 5-33 PVr characteristic and curve fitting for Innamincka generator #1
(5-42)
197
Figure 5-34 The comparison of the rotor modes of the system with the equivalent damping
torque (blue star) and the system fitted with PSS (red dot) under the high loading condition
and light loading condition. The damping gain in Innamincka generators are increased
from zero (no PSS in service) to 30 pu in 5 pu steps.
0.3
without PSS(high)
with PSS(high)
without PSS(Light)
with PSS(Light)
P.ING1(p.u.)
0.2
0.1
0
-0.1
-0.2
0
10
Time(Seconds)
15
20
Figure 5-35 Comparison on the power output of generator #1 from Innamincka under the
condition with and without PSS in service by applying a small disturbance 0.01pu on the
reference voltage of Innamincka generator #1
198
(a)
(b)
Figure 5-36 (a) Damping torque coefficients introduced by the PSS in the Innamincka
generator #1 (red curve) compared with the specified damping torque coefficient (blue
curve) (b) induced synchronizing torque (De equals to be 20 pu on Mbase)
5.3.3.2 POD Design and Assessment
The auxiliary device fitted to the FACTS devices for enhancing power system stability
are called the power oscillation damping controllers, which is designed to enhance the
damping performance of certain inter-area modes. Based on the source of the feedback
signals, two types of POD, local POD and wide area POD (WAPOD) are designed
systematically based on the residue method. Note that the residue analysis gives an
indication on the selection of the best suitable feedback signal which is the one with the
largest residue, since such signal is defined to be capable of improving the targeting
oscillation mode with less efforts [153].
Figure 5-37 shows the general structure of the AC/DC hybrid network with POD
damping controllers. The POD controllers are composed of an amplification block, a wash-
199
G( s) k
sTw 1 sTlead mc
1
[
]
1 sTLP 1 sTw 1 sTlag
(5-44)
where, TLP denotes the time constant of the low-pass filter; similarly, Tw stands for the
time constant of wash-out filter; Tlead and Tlag imply time constants for the lead-lag block;
mc is the number if compensation stages (usually mc = 2).
1
1 sTLP
1 sTlead
1 sTlag
sTw
1 sTw
PL 34 , BI
Pmax
Pmin
POWER
Rdc _ rec
SYSTEM
Ldc _ rec
Ldc _ inv
Rdc _ inv
I dc _ inv
I dc _ rec
Cdc _ rec
U dc _ rec
Rdc _ rec
+-
Ki
s
Ldc _ rec
Ldc _ inv
Rdc _ inv
Cdc _ inv
Power Plant
I cq
Kp
1
1 Tf s
U dc _ inv
U cdc
+
+
+-
I cqref
DB
Vcqref
, where Rh is
the open-loop residue from the selected input to the output of the MIMO system for the
particular hth complex mode. Furthermore, assuming that the POD controller G(s) is in
place, the closed-loop transfer function for the POD loop and the balanced equation for h
via the eigenvalue calculation can be obtained. Furthermore, by applying the small signal
technique and first order Taylor series expansion, it can be obtained the left-shift of a
target mode caused by the POD controllers satisfies the below equation,
200
h k Rh G(h )
(5-45)
1
G ( s)
Rh (
)
s s h
(5-46)
(5-47)
Tlead
1 sin(
1 sin(
i c
comp
mc
comp
mc
(5-48)
)
; Tlag c Tlag
(5-49)
A. Residue Analysis
Table 5-4 presented in the following shows the residue analysis for the interconnected
system. For the local POD, the power flow in adjacent AC lines PL34 is selected as the
input signal since it has relatively higher residue [98]. For the same reason, the bus voltage
angle difference between Brisbane and Innamincka BI is chosen as the input signal for
the WAPOD.
Note that in Table 5-4, 17 different modes are considered. In this table, two types of
stabilizing input signals are considered: (i) local signals such as power flow in adjacent AC
lines (PL12, PL34) and (ii) wide-area signals such as bus voltage angles at key nodes
(VAB,VAS,VAM,VAA,VAI) in the various regions of the system. For explanations of the
symbols, refer to Figure 5-1. In this table, and denote the bus voltage angle and
frequency respectively, while the subscript B, S, M, A, I denote the 5 areas Brisbane,
Sydney, Melbourne, Adelaide and Innamincka respectively.
201
10
11
12
13
14
15
16
17
202
57.841
70.108
-71.047 62.521
57.854
-67.306 30.107
-0.7921 46.431
58.184
7.3503
53.306
53.037
33.911
8.6735
1.06E-05 6.53E-05
-61.049 -87.127
3.12E-07 2.02E-06
8.6169
49.935
61.9
-0.6005 -32.847
-43.598 88.446
70.661
53.367
89.679
6.28E-10 1.02E-08 3.68E-08 2.68E-08 2.93E-12 2.73E-10 2.02E-09 1.34E-08 9.96E-09 2.55E-12 3.23E-09 1.06E-10 7.83E-06
-48.311 -32.77
2.9391
0.005093
-54.016 84.466
-50.364 67.417
69.314
-14.135
6.00E-12 7.60E-11 4.03E-10 3.39E-10 3.52E-13 3.28E-11 5.14E-11 3.43E-10 2.54E-10 6.43E-14 8.78E-11 2.93E-12 2.33E-07
-61.132 -87.197
2.76E-05
-3.3829 -29.084
33.923
-39.292 44.134
3.12E-07 2.03E-06
8.5585
4.4388
-88.415 0.1195
62.615
-23.45
54.818
1.07E-05 6.53E-05
47.711
0.0036817 6.01E-10 1.24E-09 1.32E-08 1.71E-08 9.24E-23 1.10E-20 1.93E-09 1.36E-08 4.12E-08 2.55E-12 3.24E-09 1.08E-10 1.93E-05
78.327
-64.707 10.01
-62.937 5.6807
83.933
-24.17
63.324
5.96E-12 9.05E-12 1.45E-10 2.14E-10 3.08E-24 3.69E-22 4.84E-11 3.50E-10 1.02E-09 6.45E-14 8.79E-11 3.00E-12 5.54E-07
-86.186 72.057
2.01E-05
32.683
34.932
-60.123 -87.256
43.344
81.031
3.05E-07 2.02E-06
7.2523
30.911
17.075
49.151
1.04E-05 6.53E-05
43.926
-38.049 44.862
0.0014634 2.75E-10 1.21E-08 2.04E-08 1.01E-08 1.09E-22 9.96E-21 1.79E-09 2.65E-08 3.69E-09 2.36E-12 3.23E-09 1.56E-10 1.36E-07
-67.376 41.226
-66.381 52.442
4.4266
-55.874
18.144
2.72E-12 8.86E-11 2.24E-10 1.26E-10 3.65E-24 3.34E-22 4.49E-11 6.83E-10 9.13E-11 5.97E-14 8.76E-11 4.34E-12 3.91E-09
26.285
7.99E-06
-88.562
-88.66
-3.715
3.76E-07 2.16E-06
21.912
76.37
-46.966 -41.833
67.228
-2.3052
-30.137 49.683
8.3962
1.28E-05 6.96E-05
-50.128 19.393
-62.319 28.663
0.0006068 7.50E-11 1.19E-08 6.81E-09 2.14E-09 1.47E-22 1.39E-20 2.57E-09 2.35E-08 5.26E-10 2.52E-12 3.40E-09 2.51E-10 1.62E-07
-45.802 21.757
-60.697
2.31E-06
26.162
7.96E-08
-59.069
2.16E-06
25.2
7.31E-08
-35.928
1.43E-06
48.34
4.82E-08
-23.082
1.00E-04
61.186
3.40E-06
Local Signal
-54.671 40.21
8.9947
75.545
-89.091
17.186
7.43E-13 8.67E-11 7.46E-11 2.67E-11 4.90E-24 4.66E-22 6.45E-11 6.05E-10 1.30E-11 6.37E-14 9.23E-11 7.00E-12 4.64E-09
3.31E-06
-86.268 57.894
63.213
-75.286 52.865
73.425
89.684
7.977
7.11E-07
18.854
0.0066914 2.54E-09 4.01E-08 2.04E-07 1.92E-07 1.45E-22 5.76E-20 1.94E-09 6.46E-09 2.80E-07 9.22E-14 4.28E-12 4.51E-11 0.0002839 1.59E-07 2.73E-07
25.881
Mode
78.691
9.8428
-3.4094
0.0045244 5.79E-10 9.15E-08 5.56E-08 1.83E-08 8.53E-21 1.20E-18 2.20E-08 1.97E-07 1.03E-08 2.46E-11 3.18E-08 2.71E-09 8.74E-06
Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang Mag Ang
(5-50)
Figure 5-38 shows the effectiveness of the designed POD control. The result meets the
compensation requirement, since sum of the compensated angle of the POD (92.3o) and the
residue angle of the target mode (89.684o) is approximate 180o, which ensures the
horizontally left-shift of the mode.
Bode Diagram
Magnitude (dB)
40
20
0
-20
-40
Frequency(rad/s:1.7)
Phase(deg):92.3
Phase (deg)
-60
135
90
45
0
-45
-90
-2
10
10
10
Frequency (rad/s)
10
Figure 5-38 Frequency response of lead compensator Q(s) for the local POD
203
Damping Criterion
Damping Criterion
-0.14Np/s
I50
-0.14Np/s
I40
I40
I35
I20
I30
I20
I25
I15
I25
Figure 5-39 The eigenvalue map of system fitted with POD, where the letters represent the
gains of the POD and KSS is increased from 3.642 10-3 to 7.03642 10-2 with a step size of
0.01.
204
12
Damping Criterion
10
-0.14Np/s
Damping Criterion
Im
I40
I35
I25
I15
I20
0
-4
-3
-2
Re
-1
-0.14Np/s
I50
I40
I30
I25
2
1
0
0 -1.5
-1
-0.5
Re
0.3
0.2
P.ING1(p.u.)
Im
0.1
0
-0.1
-0.2
-0.3
0
4
6
Time(Seconds)
10
Figure 5-41 Power output of Innamincka generator #1 following a 0.01 pu step change on
the voltage reference of Innamincka generators.
205
DC-INV (p.u.)
P-REC (p.u.)
0.015
0.01
0.005
10
15
20
10
x 10
-5
Time(Seconds)
-1
-2
15
20
-6
x 10
AC-INV (p.u.)
AC-REC (p.u.)
-5
10
Time(Seconds)
10
15
20
Time(Seconds)
x 10
0
-5
-10
-15
10
15
20
Time(Seconds)
Figure 5-42 The performance evaluation of the VSC controllers with the supplementary
controllers in service following a step change of 0.01 pu on power reference of the rectifier
side converter.
to
of the lightly damped mode 17 is 70.661o. Similarly, for the WAPOD, the results of
the calculated Tlead and Tlag are 1.8435s and 0.18697s respectively. The transfer function
for the WAPOD can be obtained by
10s
1
1+1.8435s 2
POD( s) KSS Q( s) KSS (
)(
)(
)
1 10s 1 0.01s 1+0.18697s
(5-51)
Figure 5-43 shows the frequency response of the designed lead compensator Q(s) for
WAPOD, which meet the design requirements.
206
Bode Diagram
Magnitude (dB)
40
20
0
-20
-40
Phase (deg)
-60
135
90
45
0
-45
-90
-2
10
10
10
Frequency (rad/s)
10
Damping Criterion
Damping Criterion
-0.14Np/s
-0.14Np/s
F*
Figure 5-44 The eigenvalue map of system fitted with WAPOD under high loading
condition, where the letters represent the gains of the WAPOD, as KSS is increased from
3.64210-4 to 0.12 with a step size of 0.02.
207
Damping Criterion
-0.14Np/s
Damping Criterion
-0.14Np/s
Figure 5-45 The eigenvalue map of system fitted with WAPOD under light loading
condition, where the letters represent the gains of the WAPOD, while KSS is increased
from 3.64210-4 to 0.12 with a step size of 0.02.
The simulation results suggest that the involvement of the WAPOD does not
deteriorate the light loading scenario, but not significantly improve it as well (see Figure
5-45).
b) Eigenvalue Analysis on WAPOD in cooperation with PSS
Figure 5-46 shows the eigenvalue analysis of the system operating under high loading
condition by employing the WAPOD damping controller and the PSS controller. It can be
observed from the figure that the damping performance of the system has been improved
significantly with employing the PSS only. This is further confirmed by the time response
plot shown in Figure 5-47, which demonstrate that the adequate damping performance can
be solely obtained by the WAPOD controller in service. In addition, it is shown that the
performance of WAPOD together with PSS is better than that of local POD with PSS
controller (red curve and magenta curve in Figure 5-47). Figure 5-48 illustrates the
performance of the designed VSC-controllers is still satisfactory as they operate well even
with the additional WAPOD controller.
208
Damping Criterion
-0.14Np/s
Damping Criterion
-0.14Np/s
Figure 5-46 Eigenvalue analysis of WAPOD with PSS under high loading condition,
where KSS is increased from 3.642 10-4 to 0.12 with a step size of 0.02
0.4
0.3
P.ING1(p.u.)
0.2
0.1
0
-0.1
-0.2
-0.3
0
4
6
Time(Seconds)
10
Figure 5-47 Power outputs of Innamincka generator #1 following a step change of 0.01 pu
on voltage reference of Innamincka generators.
209
-5
DC-INV (p.u.)
0.01
0.005
AC-REC (p.u.)
10
15
20
-5
10
15
Time(Seconds)
-6
x 10
-1
Time(Seconds)
-2
10
-5
AC-INV (p.u.)
P-REC (p.u.)
0.015
x 10
10
15
20
Time(Seconds)
20
x 10
0
-5
-10
-15
10
15
20
Time(Seconds)
Figure 5-48 The performance evaluation of the VSC controllers with PSS and WAPOD in
service following a step change of 0.01 pu on power reference of the rectifier side
converter.
210
assumption that designing the control system of VSC-HVDC based on the first order grid
model is sufficiently accurate.
Furthermore, for analytical purposes, a simplified model of the Australian power
system is used to connect the high capacity, but as yet undeveloped, geothermal resource in
the region of Innamincka in northern South Australia via a 1,100 km HVDC link to
Armadale in northern New South Wales. By diagnosing the lightly damped inter-area
modes of the hybrid AC/DC grid, it is found that the introduction of the new source of
geothermal power generation has an adverse impact on the damping performance of the
system. Hence, the supplementary damping controllers (i) generator power system
stabilizers (PSS) fitted to the synchronous machines which are used to convert geothermal
energy to electrical power and (ii) power oscillation damping controllers fitted to the VSCHVDC link are well designed to provide adequate damping in rotor modes of oscillation.
In the design of POD, two types of stabilizing input signals are considered: (i) local signals
such as power flow in adjacent AC lines and (ii) wide-area signals such as bus voltage
angles at key nodes in the various regions of the system. It was concluded that the smallsignal rotor-angle stability of the interconnected AC/DC system was greatly enhanced by
employing the designed damping controllers. Lastly, it should be emphasized that the
analysis tool utilized in this chapter is MATLAB m-file.
211
CHAPTER 6: CONCLUSION
Chapter 6: Conclusion
6.1 Summary
Due to the increasing demand on producing energy from renewable resources, VSCHVDC has become a common and effective solution for the grid integration of such energy
resources, which also promoted the ongoing developments in the field. As it is a widely
employed and proven VSC-HVDC technique, it can not only increase the transmission
capability but also enhance the damping ability and transient ability of the system by
employing an additional damping controller.
Significant levels of renewable energy resources (such as wind, solar and geothermal)
are available worldwide and economically harvested locally. However, the locations of
such resources are usually far from the main load centres and cities. Therefore, utilization
of these energies is very challenging. For example, the power grid in Australia is one of the
worlds longest interconnected power systems starting from Port Douglas in Queensland
and ending in Port Lincoln in South Australia with a total distance more than 4000
kilometres.
212
6.1. SUMMARY
The motivation of this research study has been to design and develop a detailed
mathematical understanding of the control system of a VSC-HVDC that can be integrated
into such a long and weak AC system.
As presented in the thesis, the outer power loop, the inner current loop, synchronization
method, and the input-output impedance of a weak AC interconnected VSC-HVDC system
have been partially studied in the literature. However, most of these earlier studies
typically employed a single-machine infinite-bus test system and did not study a wide
range of operating conditions that is likely to be encountered in practice. Therefore, one of
the primary aims of this thesis was to design a robust control system with a simple
structure and study a wide range of loading conditions that may occur in practice.
To achieve this main aim, this research can be divided into two main sections. The first
section (Chapter 2-4) of the research has investigated three potential linear inner current
control schemes to determine the most suitable type. Then an analytical model of the DB
controlled VSC was developed as the suitable controller. The discussions in the thesis
include a set of linearized representation of system components which prepared the ground
work. This resulted in the development of a new method for eliminating the
modulator/demodulator blocks. Using this new model, the classical frequency response
technique was also applied to a set of linear models using a number of operating points to
develop a robust outer loop controller.
In the second section of the research (Chapter 5), the impact of an admittance
representation of the grid to the controller design in the VSC-HVDC system and a scaling
approach were investigated. The investigations of this section also include the integration
of the designed VSC-HVDC links into the Extended Simplified South-East Australian
power grid, which aimed to form a large AC/DC integrated power system. In addition,
extensive eigen-sensitivity studies were performed to determine the characteristics of the
interaction and to identify the critical parameters. Finally, the power damping controllers
as PSS, local POD and WAPOD were designed to enhance the power system stability.
213
CHAPTER 6: CONCLUSION
A new method for the selection and optimization of the parameters of the PI
compensators in the various control loops using a decoupled control strategy was
proposed. In this study, the initial value of the PI compensator parameters for
input to the optimization algorithm were obtained using the classical frequency
response design approach to simplified linear models of the open-loop transfer
functions of VSC-HVDC control system. An optimization algorithm based on the
simplex method was adopted. The objective was to simultaneously minimize the
weighted sum of the integral of the time absolute-error products (ITAE) of the
active power, the reactive-power, the DC voltage and the inner current controllers
of the respective VSCs. It was concluded that if necessary the weightings of the
error signals may need to be modified in this stage.
214
Four different DB controllers were also studied in the thesis by including one
sample delay with the reduced gain and Smith Predictor, two sample delays
based on IMC control design and solving feedback transfer function. It was
concluded from the simulation results that the solving feedback transfer function
method proves to be more superior to the other methods specifically during the
large step response with a fast settling time and with an acceptable reactive
current overshoot
2. A new small signal model for the digital DB controlled VSC was developed
(Chapter 3).
Three approximation methods were investigated and compared: the first order
Pad Approximation, the third order Pad Approximation and first order lag
approximation since it was essential to model the ZOH, the DB control block and
the delay block. In this investigation the first-order Pad Approximation was
found to be accurate enough to approximate the pure one sample delay.
PLL was also investigated in detail. The principle of PLL, and the parameter
design were studied and verified using the small signal and the large signal
system disturbances.
In this part of the study, a new approach to eliminate the modulator and
demodulator block is proposed and realized the freely transformation among abc
natural, and dq synchronous reference frames. In addition, the limitation of
this methodology was also examined in this thesis, based on which it is suggested
that the best strategy for controller design is to design in dq synchronous
reference frame but implemented in abc natural reference frame.
The small signal model for VSC itself adopted here was the average value model
(AVM) which was verified in this project. Although it was found that to
represent it as an ideal converter is not so accurate due to its non-linearity nature,
the deviation of final results contributed by these differences was small once the
loops of the whole system are closed.
215
CHAPTER 6: CONCLUSION
Every component of the control system was characterized and verified step by
step by comparing the small signal and the large signal models. This provided a
solid foundation to verify the accuracy of the complete model. Furthermore, the
state-space models of the system components were also developed to obtain the
model of a large system.
The limitation imposed by the RHP zero on the bandwidth of the outer loop
controllers was also discussed.
216
Finally, the small signal model development for each controller such as PID, PI
with low pass filter was presented.
4. To avoid the adverse impacts of the introduction of the new source of geothermal
power generation where the power are delivered via VSC-HVDC links fed into a
weak AC grid (Australian), three types of damping controllers (PSS, POD and
WAPOD) were designed to enhance the small-signal rotor-angle stability of the
power system. To achieve this, the following tasks were completed (Chapter 5).
A scaling technique was introduced which made the model applicable to the
systems with different ratings.
The parameter sensitivity, in terms of the length of the DC link, was also
examined in the thesis, which did not weaken the performance of the designed
controller, but enhanced the system stability to some extent.
Back examination of the basement for the controller design against a higher order
detailed Australian grid confirmed the accuracy of the assumption that design the
control system of VSC-HVDC based on first order grid model.
It was observed that the introduction of the new source of geothermal power
generation has an adverse impact on the damping the inter-area modes by using
the small signal stability analysis.
Three types of damping controllers (PSS, POD and WAPOD) were designed and
verified systematically to enhance the power system stability.
217
CHAPTER 6: CONCLUSION
The dynamics of the grid frequency can be introduced for the small signal
model of DB controlled VSC, which may offer a better precision for the model.
The control strategy employed in this thesis can be improved to allow a much
weaker system to be studied to increase the operating limits of the conventional DB
controlled VSC.
Finally, the DB current controller can be applied to a grid integrated energy storage
system, and the system stability can be studied. Such study can provide in-depth
understanding on how the energy storage system can enhance the power system
security.
218
219
APPENDICES
Rated Power
rec
inv
Rated ac Voltage
Rated Frequency
Switching Frequency
Sampling
ratio
commuting transformer
equivalent
inductance
L(pu)
R(pu)
Equivalent Filter
L(pu)
rec
inv
rec
0.15 pu
inv
rec
inv
rec
inv
Rated dc Voltage
DC-link Capacitor
0.0015 pu
0.15 pu
130 kV
C=500F
rec
inv
Pn=75MW
Vn,LL(RMS)=13.8kV
Vn,LL(RMS)=115kV
fn=50Hz
fsw=1350Hz
fsamp=1350Hz
Tsamp=7.410-4
13.8 kV / 62.5 kV
62.5 kV / 115 kV
Ub ( rec ) 2 / 3 U n, LL ( rec )
Ub (inv ) 2 / 3 U n, LL (inv )
rec
2 S
Ib ( rec ) b
3 U b ( rec )
inv
2 S
I b (inv ) b
3 U b (inv )
Ub ( dc ) 130kV
Ib( dc ) Sb / Ub ( dc )
C pu CUb2( dc ) / Sb [116]
220
Z A0 Z B0u
y C0 Z D0u
(B.1)
u
y
u ; y
u
y
G( s)
y
u
C0 ( sI A0 )1 B0 D0
(B.2)
221
APPENDICES
Zt At Zt Bt udq
ydq Ct Zt Dt udq
(B.3)
ud
yd
udq ; ydq
yq
uq
u e j (0t )udq
j (0t )
ydq
y e
(B.4)
Z AZ Be j (0t )udq
j (0t )
udq
y CZ De
(B.5)
cos(0t ) sin(0t )
Z AZ BRudq
Set R
sin(0t ) cos(0t )
y CZ DRudq
(B.6)
In which,
A
A 0
0
0
B
;B 0
A0
0
0
C
;C 0
B0
0
0
D
;D 0
C0
0
0
D0
Z AZ BRudq
(B.7)
R1 (Z AZ ) R1BRxdq
(B.8)
0 1
B0
R 1 ( Z AZ ) B0udq
R 1Z R 1 AZ B0udq
222
(B.10)
R 1 Z ( R 1 ) Z - ( R 1 ) Z R 1 AZ B0udq
( R 1 Z ) ( R 1 ) Z R 1 AZ B0udq
R 1 e j0t ;( R 1 ) j0 e j0t j0 R 1
( R 1 Z ) j0 R 1 Z R 1 AZ B0udq
cos(0t ) sin(0t ) A0 0 cos(0t ) A0 sin(0t ) A0
R 1 A
( R 1 Z ) j0 ( R 1 Z ) A( R 1Z ) B0udq
s ( R 1 Z ) j0 ( R 1 Z ) A( R 1Z ) B0udq
R B0 udq
B0
R 1 Z
Z
udq
[( s j ) A]
[( s j ) A]
(B.11)
(B.12)
G( s)
ydq
udq
(B.13)
G ( s) C0 [( s j ) A]1 B0 D0
Therefore, the transfer function transferred from
y CZ Du
y CZ Du
(B.14)
( s j ) Z dq Adq Z dq Bdqudq
sZ dq ( Adq j ) Z dq Bdqudq
(B.15)
223
APPENDICES
Z d Adq Z d Bdq 0 ud
s
Z q Adq Z q 0 Bdq uq
0 ud
yd Cdq 0 Z d Ddq
y
0 Cdq Z q 0 Ddq uq
q
(B.16)
G ( s) dq (C
B D) jC
B
2
2
( sI A)
( sI A) 2 2
sI A
Yd jYq [(C
B D ) jC
B](U d jU q )
2
2
( sI A)
( sI A) 2 2
sI A
Yd (C
B D)U d C
BU q
2
2
( sI A)
( sI A) 2 2
sI A
Yq C
BU d (C
B D)U q
2
2
( sI A)
( sI A) 2 2
Yq
Yq
Yd
Y
sI A
sI A
C
B D; d C
B ; C
B; C
BD
2
2
2
2
2
2
Ud
( sI A)
Uq
( sI A) U d
( sI A) U q
( sI A) 2 2
G ( s ) dq C
224
q
B
0
(C.1)
sin ud
cos uq
cos In sin In
sin In cos In
Aside: Compute
Zd
Z
q
(C.2)
According to the Woodbury matrix identity theory, this matrix has the form that
1
A1 A1BPCA1 A1BP E F
A B
C D
PCA1
P G H
(C.3)
225
APPENDICES
where, P ( D CA1B)1
A cos In
B sin In
C sin In
D cos In
(C.4)
1
In
cos
(C.5)
A1
P [cos In sin In(
1
) In( sin In)]1
cos
cos 2 sin 2 1
P [
] In cos In
cos
E ( A1 A1 BPCA1 )
(C.6)
from(C.3)
1
1
1
In (
In)( sin In) (cos In)(sin In) (
In)
cos
cos
cos
1
sin 2
In
In
cos
cos
1
(
In)(1 sin 2 ) In
cos
cos In
F A1 BP
1
In)( sin In)(cos In)
cos
F sin In from(C.3)
F (
G PCA
1
) In
cos
from(C.3)
Thus, equation (C.2) can be written as equation (C.7) according to equation (C.3).
1
cos In sin In
E
sin In cos In G
F cos In sin In
H sin In cos In
226
(C.7)
Z q sin In cos In 0 A sin In cos In Z q
cos In sin In sin In cos In d
Zd
Z
q
cos In sin In B 0 cos sin ud
sin In cos In 0 B sin cos uq
cos A sin A cos In sin In Z d
sin A cos A sin In cos In Z q
(cos sin sin cos ) In
d
(cos 2 sin 2 ) In
2
2
(sin cos ) In
( sin cos cos sin ) In dt
sin B cos B sin In cos In uq
Zd
Z
q
Zd
(cos 2 sin 2 ) A
( cos sin sin cos ) A Z d 0 In d Z d
Z
(sin 2 cos 2 ) A
q In 0 dt Z q
Z q ( sin cos cos sin ) A
(cos 2 sin 2 ) B
( cos sin sin cos ) B ud
u
(sin 2 cos 2 ) B
( sin cos cos sin ) B
q
d
Now assume that, which is significantly important,
0
dt
Zd A
Z q In0
In0 Z d B 0 ud
A Z q 0 B uq
(C.8)
yd cos In sin In C
y
q sin In cos In 0
cos C
sin C
yd C 0 Z d D
y
q 0 C Zq 0
0 cos
D sin
sin ud
cos uq
sin ud
cos uq
sin ud
cos uq
0 ud
D uq
(C.9)
227
APPENDICES
In0 Z d B 0 ud
A Z q 0 B uq
Zd A
Z q In0
(C.10)
yd C 0 Z d D 0 u d
Z 0 D u
y
0
C
q
q
q
These two equations can be combined as equations (3-86).
B. Compute the Frequency Response Matrix
sZ dq Adq Z dq Bdqudq
( sI Adq ) Z dq Bdqudq
Z dq ( sI Adq ) 1 Bdqudq
ydq
udq
(sIn Adq )
sIn 0 In Adq 0 In
(C.11)
1
(C.12)
a sIn Adq ;
b 0 In;
(C.13)
c 0 In;
d sIn Adq ;
P (d ca 1b) 1
a ca 1b
a 0 Ina 1 (0 In)
(a 02 a 1 ) 1
228
(C.14)
E a 1 a 1bpca 1
a 1 a 1b(a 02 a 1 ) 1 ca 1
a 1 02 a 1 (a 02 a 1 ) 1 a 1
a 1 ( In 02 (a 02 a 1 ) 1 a 1 )
(C.15)
a 1 ( In 02 (a (a 02 a 1 )) 1 )
a 1 ( In 02 (a 2 02 In) 1 )
a(a 2 02 In) 1
F a 1bp
a 10 (a 02 a 1 ) 1
0 a 1 (a 02 a 1 ) 1
(C.16)
0 ((a 02 a 1 )a) 1
0 (a 2 02 In) 1
G pca 1
(a 02 a 1 ) 10 Ina 1
0 ((a 02 a 1 ) 1 a 1 )
(C.17)
0 (a(a 02 a 1 )) 1
0 (a 2 02 In) 1
Hp
(a 02 a 1 )1
The Woodbury identity is
Then:
a 1 ( In 02 (a(a 02 a 1 )) 1 )
a 1 ( In 02 (a 2 02 In) 1 )
a(a 2 02 In) 1
H=E
(C.18)
Thus
E ( s) F ( s) a(a 2 02 In)1 0 (a 2 02 In)1
( sI Adq )1
2
2
1
a(a 2 02 In)1 (C.19)
F ( s) E ( s) 0 (a 0 In)
229
APPENDICES
2
2 1
( sI A)[( sI A)2 02 ]
0 [( sI A) 0 ]
1
0 C F (s) E (s) 0 B 0
CE ( s ) CF ( s) B 0 D 0
CF ( s ) CE ( s) 0 B 0 D
CF ( s ) B
CE ( s ) B D
CF ( s ) B CE ( s ) B D
Gdd ( s ) Gdq ( s )
Gqd ( s ) Gqq ( s )
(C.20)
(C.21)
0
D
(C.22)
Gdd ( s ) CE ( s) B D
C ( sI A)[( sI A) 2 02 ]1 B D
Gdq ( s) CF ( s) B
C 0 [( sI A) 2 02 ]1 B
(C.23)
Gqd ( s ) CF ( s ) B
C 0 [( sI A) 2 02 ]1 B
Gdq ( s )
230
x
k
1 Ta s
1 Tb s
x k (1 Ta s ) y (1 Tb s )
kx kxTa s y yTb s
kx y s (Tb y KTa x)
(D.1)
T
kx y
s( y k a x)
Tb Tb
Tb
set y k
sz
Ta
T
x z y zk a x
Tb
Tb
T
k
1
x ( z k a x)
Tb
Tb
Tb
(D.2)
(D.3)
231
APPENDICES
1
k k Ta
sz T z ( T T T ) x
b
b
b b
y zk a x
Tb
T
1
k k Ta
Set A ; B
; C 1; D k a
Tb
Tb Tb Tb
Tb
232
(D.4)
Appendix E: Regional Boundaries for the National Electricity Market & Committed
Developments
233
APPENDICES
234
Appendix F: Performance Evaluation of all the Power and Voltage Controllers with a New
set of DC Link Parameters
Magnitude(dB)
Magnitude(dB)
40
20
0
-20 -1
10
10
10
10
10
10
10
10
10
Phase(deg)
Phase(deg)
-500 -1
10
-50
-100 -1
10
500
Phalf-Qmax
Phalf-Qmin
P -Qmin
0
10
10
Frequency (rad/sec)
(a)
10
10
-500
-1000 -1
10
Phalf-Qmax
Phalf-Qmin
P -Qmin
0
10
10
Frequency (rad/sec)
10
10
(b)
Figure F-1 Open loop frequency response (a) without and (b) with power controller under
the updated DC link condition
235
APPENDICES
1.4
System: untitled9
Peak amplitude: 1.04
Overshoot (%): 4.35
At time (seconds): 0.219
1.2
1
System: untitled12
Settling time (seconds): 0.506
pu
0.8
0.6
0.4
0.2
0
0
0.2
0.4
0.6
0.8
Time (s) (seconds)
1.2
Figure F-2 Step response on power reference with updated DC link parameters
B.
50
Magnitude(dB)
Magnitude(dB)
40
20
0
-20 0
10
10
10
10
10
10
10
Phase(deg)
Phase(deg)
-400 0
10
-50
-100 0
10
-200
case 3
case 7
case 9
case 36
1
10
10
Frequency (rad/sec)
10
(a)
-200
-400
-600 0
10
case 3
case 7
case 9
case 36
1
10
10
Frequency (rad/sec)
10
(b)
Figure F-3 Open loop frequency response (a) without and (b) with rectifier side AC voltage
controller under the updated DC link condition
236
Appendix F: Performance Evaluation of all the Power and Voltage Controllers with a New
set of DC Link Parameters
1.2
System: untitled12
Peak amplitude: 1.09
System: untitled36
Overshoot (%): 9.02
At time (seconds): 0.123 Settling time (seconds): 0.337
0.8
pu
0.6
0.4
0.2
0
-0.2
0
0.1
0.2
0.3
0.4
Time (s) (seconds)
0.5
0.6
Figure F-4 Step response on rectifier side AC voltage reference with updated DC link
parameters
C.
Magnitude(dB)
Magnitude(dB)
50
0
-50
-100
10
10
10
-100
-200 0
10
200
10
10
10
10
Phase(deg)
Phase(deg)
100
0
-100
-200
10
10
Frequency (rad/sec)
(a)
10
-500
-1000 0
10
10
10
Frequency (rad/sec)
10
10
(b)
Figure F-5 Open loop frequency response (a) without and (b) with inverter side DC voltage
controller under the updated DC link condition
237
APPENDICES
1.2
System: untitled1
Peak amplitude: 1.02
Overshoot (%): 1.97
At time (seconds): 0.426
1
0.8
System: untitled18
Settling time (seconds): 0.483
pu
0.6
0.4
0.2
0
-0.2
0
0.5
1
Time (s) (seconds)
1.5
Figure F-6 Step response on inverter side DC voltage reference with updated DC link
parameters
D.
A.
(a)
(b)
Figure F-7 Open loop frequency response (a) without and (b) with inverter side DC voltage
controller under the updated DC link condition
238
Appendix F: Performance Evaluation of all the Power and Voltage Controllers with a New
set of DC Link Parameters
1.2
1
System: untitled3
Peak amplitude: 1.06
Overshoot (%): 5.82
At time (seconds): 0.135
0.8
System: untitled36
Settling time (seconds): 0.452
pu
0.6
0.4
0.2
0
-0.2
0
0.1
0.2
0.3
0.4
Time (s) (seconds)
0.5
0.6
0.7
Figure F-8 Step response on inverter side AC voltage reference with updated DC link
parameters
E.
Step Tests for System with All Controllers in service under Updated DC
Link Condition
239
APPENDICES
(c)
240
Appendix F: Performance Evaluation of all the Power and Voltage Controllers with a New
set of DC Link Parameters
c) Step test on inverter side DC voltage controller
(c)
241
APPENDICES
242
75MW
=51.03kV
=0.9798kA
0.0015pu=0.0015*
0.15pu=0.15*
Base frequency
2 50Hz
0.0781
=0.0249H
*The base for per unit transformation is chosen as to achieve a power invariant transformation,
so that the ac and dc side power is the same.
Basement for DC system
243
APPENDICES
Nominal DC voltage
130kV
Nominal DC current
0.5769kA
Nominal DC impedance
225.33
Time constant
244
Reference
Reference
[1]
[2]
[3]
[4]
[5]
[6]
[7]
W. Pan, Y. Chang, and H. Chen, "Hybrid Multi-terminal HVDC system for large
scale wind power," in Power Systems Conference and Exposition, 2006. PSCE'06.
2006 IEEE PES, 2006, pp. 755-759.
245
APPENDICES
[8]
[9]
J. Liang, O. Gomis-Bellmunt, J. Ekanayake, and N. Jenkins, "Control of multiterminal VSC-HVDC transmission for offshore wind power," in Power Electronics
and Applications, 2009. EPE'09. 13th European Conference on, 2009, pp. 1-10.
[10]
[11]
[12]
[13]
[14]
(12June).
List
of
HVDC
http://en.wikipedia.org/wiki/List_of_HVDC_projects
[15]
[16]
[17]
B. T. Ooi and X. Wang, "Voltage angle lock loop control of the boost type PWM
converter for HVDC application," Power Electronics, IEEE Transactions on, vol.
5, pp. 229-235, 1990.
[18]
X. Wang and B. T. Ooi, "High voltage direct current transmission system based on
voltage source converters," in Power Electronics Specialists Conference, 1990.
PESC '90 Record., 21st Annual IEEE, 1990, pp. 325-332.
[19]
[20]
246
projects.
Available:
Reference
[21]
[22]
[23]
[24]
[25]
[26]
[27]
[28]
[29]
Z. Chengyong, L. Xiangdong, and L. Guangkai, "Parameters Optimization of VSCHVDC Control System Based on Simplex Algorithm," in Power Engineering
Society General Meeting, 2007. IEEE, 2007, pp. 1-7.
[30]
[31]
[32]
[33]
247
APPENDICES
[34]
W. Liying and N. Ertugrul, "Selection of PI compensator parameters for VSCHVDC system using decoupled control strategy," in Universities Power
Engineering Conference (AUPEC), 2010 20th Australasian, 2009, pp. 1-7.
[35]
[36]
[37]
[38]
[39]
[40]
[41]
[42]
[43]
[44]
[45]
248
Reference
[47]
[48]
J. Dannehl, F. W. Fuchs, Th, x00F, and P. B. gersen, "PI State Space Current
Control of Grid-Connected PWM Converters With LCL Filters," Power
Electronics, IEEE Transactions on, vol. 25, pp. 2320-2330, 2010.
[49]
[50]
[51]
[52]
J.Hu, "Deatbeat Controlled PWM converter," Master, Dept. Elect. Eng., McGill
Univ., Montreal,QC, Canada, 1999.
[53]
[54]
[55]
[56]
D. Jovcic, N. Pahalawaththa, and M. Zavahir, "Analytical modelling of HVDCHVAC systems," Power Delivery, IEEE Transactions on, vol. 14, pp. 506-511,
1999.
249
APPENDICES
[57]
[58]
[59]
[60]
[61]
[62]
Z. Lidong, L. Harnefors, and H. P. Nee, "Power-Synchronization Control of GridConnected Voltage-Source Converters," Power Systems, IEEE Transactions on,
vol. 25, pp. 809-820, 2010.
[63]
[64]
[65]
[66]
[67]
[68]
250
Reference
[69]
[70]
[71]
[72]
[73]
[74]
[75]
[76]
[77]
[78]
[79]
[80]
251
APPENDICES
[81]
[82]
[83]
[84]
[85]
[86]
[87]
[88]
[89]
[90]
[91]
[92]
252
Reference
[93]
H. F. Latorre and M. Ghandhari, "Improvement of voltage stability by using VSCHVdc," in Transmission & Distribution Conference & Exposition: Asia and
Pacific, 2009, 2009, pp. 1-4.
[94]
[95]
[96]
[97]
N. Modi, M. Lloyd, and T. K. Saha, "Wide-area signal selection for power system
damping controller," in Universities Power Engineering Conference (AUPEC),
2011 21st Australasian, 2011, pp. 1-6.
[98]
[99]
253
APPENDICES
254
Reference
255
APPENDICES
Analysis of Power Systems," Power Delivery, IEEE Transactions on, vol. 25, pp.
2655-2669, 2010.
[130] J. Jatskevich, S. D. Pekarek, and A. Davoudi, "Parametric average-value model of
synchronous machine-rectifier systems," Energy Conversion, IEEE Transactions
on, vol. 21, pp. 9-18, 2006.
[131] A. Yazdani, "Electromagnetic transients of grid-tied photovoltaic systems based on
detailed and averaged models of the voltage-sourced converter," in Power and
Energy Society General Meeting, 2011 IEEE, 2011, pp. 1-8.
[132] C. Karawita and U. D. Annakkage, "Multi-Infeed HVDC Interaction Studies Using
Small-Signal Stability Assessment," Power Delivery, IEEE Transactions on, vol.
24, pp. 910-918, 2009.
[133] S. Cole, J. Beerten, and R. Belmans, "Generalized Dynamic VSC MTDC Model for
Power System Stability Studies," Power Systems, IEEE Transactions on, vol. 25,
pp. 1655-1662, 2010.
[134] D. J. Vowles and M. J. Gibbard, "Mudpack Software Package for the SmallSignal
Analysis of Power Systems (Version 10R-19)," The University of Adelaide April
2008.
[135] P. Rault, F. Colas, X. Guillaud, and S. Nguefeu, "Method for small signal stability
analysis of VSC-MTDC grids," in Power and Energy Society General Meeting,
2012 IEEE, 2012, pp. 1-7.
[136] R. Teodorescu, F. Blaabjerg, M. Liserre, and A. Dell'Aquila, "A stable three-phase
LCL-filter based active rectifier without damping," in Industry Applications
Conference, 2003. 38th IAS Annual Meeting. Conference Record of the, 2003, pp.
1552-1557 vol.3.
[137] M. Liserre, F. Blaabjerg, and S. Hansen, "Design and control of an LCL-filterbased three-phase active rectifier," Industry Applications, IEEE Transactions on,
vol. 41, pp. 1281-1291, 2005.
[138] (24-May-13). National Electricity Rules: Current Rules (55 ed.). Available:
http://www.aemc.gov.au/electricity/national-electricity-rules/current-rules.html
[139] Wikipedia. (24 May ). Available: https://en.wikipedia.org/wiki/Newton's_method
[140] J. Z. A. M. Gole, "VSC transmission limitations imposed by AC system strength
and AC impedance characteristics," presented at the The 10th International
Conference on AC and DC Power Transmission ACDC 2012 4 - 6 Birmingham,
UK, December,2012.
[141] S. Skogestad and I. Postlethwaite, Multivariable Feedback Control: Analysis and
Design, 2nd Edition. West Sussex, England: John Wiley & Sons Ltd., 2005.
256
Reference
[142] Wikipedia, the free encyclopedia. (15 October 2009, 30th Oct). Mandatory
renewable
energy
targets.
Available:
http://en.wikipedia.org/wiki/Mandatory_renewable_energy_targets
[143] Electricity Supply Industry Planning Council "ANNUAL PLANNING REPORT ",
ed, June 2009.
[144] AEMO. "AEMO
www.aemo.com.au
Australian
Energy
Market
Operator".
Available:
257