Control of MMC in HVDC Applications Masters Thesis WPS4 1054
Control of MMC in HVDC Applications Masters Thesis WPS4 1054
Control of MMC in
HVDC Applications
Master Thesis 30/05/2013
Artjoms Timofejevs
Daniel Gamboa
Title:
Semester:
9 10th
Semester theme:
Masters Thesis
Project period:
01/10/12 to 30/05/13
ECTS:
50
Supervisor:
Remus Teodorescu
Marco Liserre
Sanjay K. Chaudhary
Project group:
WPS4-1054
_____________________________________
[Artjoms Timofejevs]
_____________________________________
[Daniel Gamboa]
Copies:
Pages, total:
Appendix:
Supplements:
[4]
[90]
[1]
[ 1 CD ]
SYNOPSIS:
The Modular multilevel converter
(MMC) is the latest converter
topology suitable for transformerless
applications in HVDC transmission.
HVDC are required to remain
connected during grid fault, provide
grid support and completely decouple
the healthy side from the faulty one.
Due to its complex structure, the
inner dynamics of the MMC converter
are challenged by this particular
condition and by these demands.
This thesis demonstrates the effect of
negative and zero sequence current
control in MMC-HVDC during
asymmetric grid faults. A modified
circulating
current
suppression
controller is proposed to eliminate
the voltage ripples in the DC-link.
A current limitation strategy for MMC
is derived and its impact on the
performance
of
the
HVDC
transmission system is verified
through simulations.
By signing this document, each member of the group confirms that all group members have participated
in the project work, and thereby all members are collectively liable for the contents of the report.
Furthermore, all group members confirm that the report does not include plagiarism.
Acknowledgements
This thesis would not have been possible without the help and support of our supervisors. We would like to
express our gratitude to Prof. Remus Teodorescu for motivating us to choose this topic and guiding us
through the learning process of this master thesis. We are thankful to Prof. Marco Liserre for his valuable
discussions and constructive feedbacks, particularly during the second part of this project. The good advice,
continuous encouragement and always available support of Dr. Sanjay K. Chaudhary, has been invaluable
on both an academic and a personal level, for which we are sincerely grateful.
Amongst our fellow master students in the Department of Energy Technology, we thank Csaba Kopacz and
Lorand Bede for their support and discussions during the complete master program. We look forward to
calling you friends for many years to come.
I thank Vestas Wind Systems A/S and Aalborg University for awarding me the Vestas Scholarship, providing
me with the financial means to complete the master studies Artjoms Timofejevs.
I would acknowledge the CONACYT Mexico and Aalborg University for the award of a scholarship that
provided me the necessary financial support during the master program Daniel Gamboa.
Table of Contents
Introduction ............................................................................................................................................... 1
1.1
1.2
Objectives .......................................................................................................................................... 4
1.3
Limitations ......................................................................................................................................... 4
1.4
Thesis outline..................................................................................................................................... 4
Background ................................................................................................................................................ 6
2.1
2.1.1
2.1.2
2.1.3
2.1.4
2.2
2.2.1
2.2.2
Multi-carrier PWM................................................................................................................... 12
2.2.3
2.3
2.3.1
2.3.2
2.4
2.4.1
2.4.2
2.4.3
3.2
3.2.1
3.2.2
3.2.3
3.3
3.3.1
3.3.2
3.4
3.5
3.5.1
3.6
4
4.1.1
4.1.2
4.1.3
4.2
4.2.1
4.2.2
4.2.3
4.2.4
4.2.5
4.2.6
4.3
5
5.2
5.3
5.4
5.4.1
5.4.2
5.5
5.5.1
5.6
5.6.1
5.6.2
5.6.3
5.7
5.7.1
5.7.2
5.8
6
Conclusions .............................................................................................................................................. 74
6.1
Future Work..................................................................................................................................... 75
References ....................................................................................................................................................... 76
Appendix .......................................................................................................................................................... 80
1 Introduction
Thanks to the global energy consciousness and the government support, more and more renewable energy
sources (RES) are being installed in order to cope with the increasing energy demand [1]. In 2009, the EU
Renewable Energy Directive set the goal of producing the 20% of the overall energy mix from RES by 2020.
[2]
Aiming to meet this target, the share of renewables in the total new power installations has grown from
20.7% in the year 2000 to 70% in 2012, as shown in Figure 1.1 [3]. The installation of Photovoltaic (PV)
generation in the EU has grown at remarkable rates over the last 5 years, moving from 2.0 MW installed in
2007, to 17.2 MW in 2012; covering the 37% of the total installed capacity of that year [4].
The renewable resources in the world are unequally distributed. In Europe, while most of the wind energy
resources are located in the north, the solar resources are located in the south, in countries like Italy and
Spain [5]. In addition, the availability of RES has a strong daily and seasonal pattern. [6]
The integration of more renewable generation into the energy mix combined with a rapid growth of the
energy consumption are evidence that the actual AC grid will not be suitable for the transmission and
distribution of the new power generation.
In order to utilize efficiently the vast potential of RES, the selection of the transmission technology is of
critical importance [7]. The high losses associated with long AC lines make HVDC an attractive alternative
for the transport of bulk power over long distances. Additionally, HVDC interconnections allow power
transmission between unsynchronized AC systems and prevent cascading failures to propagate through
Introduction
wide transmission grids; increasing the reliability of the system. Due to its stability benefits and reduced
losses, many transmission system operators (TSO) have considered a wider use of HVDC technology. [8]
Before the emergence of voltage source converters (VSC), the implementation of HVDC was limited to very
large installations located at a long distance or to under-sea cable transmission; due to the elevated costs
of the converter infrastructure, the reactive power demand and limited controllability [9]. In Europe HVDC
transmission has been used since 1954, mostly for submarine cable transmission, as shown in Figure 1.2.
Existing connection
Under construction
Future plan
The appearance of VSC-HVDC, having a smaller converter size, reduced filters, dynamic reactive power
support and fast power reversal, has facilitated the implementation of HVDC technology in a wider range of
applications [9]. Since the first successful application of VSC-HVDC in 1997 in Gotland, Sweden; many other
installations provided asynchronous interconnections of AC grids all-over the world. Due to its smaller
footprint and ability to act as a virtual synchronous generator, VSC is especially suitable for connecting
offshore wind power plants.
In 2002, Marquardt and Lesnicar proposed the modular multilevel converter (MMC) topology [11]. Due to
the fact that it can achieve high power and high voltage levels using proven semiconductor technology,
MMC has been widely accepted in the industry [12]. The HVDC-Plus from Siemens, HVDC-Light from ABB
and HVDC-MaxSine from Alstom are examples of the implementation of the MMC concept in applications
for VSC-HVDC transmission [13][14][15].
With advantages such as modularity, increased efficiency and reliability, the MMCs aims to substitute the
two-level converters in VSC-HVDC applications, becoming a backbone for the future HVDC transmission
systems [16].
1.1
The origins of the MMC circuits go to the 1970s when synthesising waveform converter topologies were
patented [17]. In the last decade, after the design and control of the MMC was proposed [11], an intense
research has been done on all aspects of the system.
During steady state operation, a circulating current is observed flowing through the converter phase-leg. In
[18] the voltage ripples in the capacitors are studied as a cause for this circulating current component. Later
in [19], the circulating current is analyzed and found to have double line-frequency and negative sequence.
Having a large number of sub-modules increases the complexity of the control of the MMC. To assure
accurate and stable operation, the sub-modules must share the voltage equally. Different control strategies
have been proposed aiming to regulate the equal share of charge within the capacitors in the arm.
Akagi et al. propose in [20] a voltage balancing algorithm for the capacitors per phase-leg. The control is
performed by adding a balancing component to the modulation index of each sub-module. Later in [21] a
modification of this method is proposed to take into consideration the voltage balance between arms.
Another method to eliminate the circulating current and balance the converter arm voltages is proposed in
[22]. The method is based on the control of the stored energy in the converter, and two control loops are
added to ensure stable operation. A simplification of this method using open-loop approach for the
estimation of the arm energy is proposed in [23], proving to have a better stability.
As alternative, a straightforward method to suppress unwanted current component is derived in [24]. The
method is based on the vector control of the measured current. Later in [25] a combination of the energy
approach presented in [22] with the method presented in [24] is implemented, showing an improved
performance for grid connected applications.
The MMC topology has the possibility of synthesising high voltage levels by increasing the number of series
connected sub-modules, therefore is especially suitable for HVDC systems. In [16] the MMC is highlighted
as the most promising VSC topology to be used in the future HVDC grids. In 2010 Siemens announced the
commercial installation of the first MMC-HVDC transmission [13]. Furthermore, several projects using
MMC-HVDC for the interconnection of large offshore wind farms are announced to be commissioned in the
next few years [14].
Lately, a lot of research is focused on evaluation of MMC performance in HVDC applications under different
grid conditions. In [26] the converter operation under unbalanced grid conditions has been analyzed, using
a -Y transformer to remove the zero sequence current from the converter terminals. The presence of DClink voltage ripple due to negative sequence components is described and a controller to compensate for
the disturbance is proposed. The transformerless grid connection of the MMC is studied in [27], observing
that under unbalanced conditions; the zero sequence current is looped between two AC grids. In [28] an
Introduction
additional zero sequence current control is proposed for transformerless connections under unbalanced
conditions, proving stable HVDC system operation and enhanced fault ride-through capability of the
converter.
Technical literature and the growing amount of installations show the success and market acceptance of
the MMC topology. However, there are very little publications discussing the implementation of MMC in
HVDC applications with transformerless connection.
Motivated by the fast penetration of MMC into the HVDC market and the lack of research published on that
topic, the performance of transformerless MMC-HVDC under unbalanced conditions has been studied here.
The analysis is focused on the fault propagation over HVDC transmission and its impact on the converter
inner dynamics.
1.2 Objectives
The objective of the thesis is modelling and analysis of a VSC-HVDC transmission system based on Modular
Multilevel Converters. The main goals of the project are:
Good understanding of operating principles and inner controls of MMC and VSC-HVDC
transmission systems;
Modelling of the MMC-HVDC system and implementation of the control strategies in
PSCAD/EMTDC software;
Study of the fault propagation in transformerless MMC-HVDC transmission systems;
Analysis of the system behaviour as response to unbalanced grid conditions;
1.3 Limitations
The main limitations in the project are:
In Chapter 3 the inner dynamics of the MMC are studied in order to have a better understanding of the
converter response for the subsequent chapters. An average model of the converter is derived and
validated through simulations. The arm currents present in the converter and the sub-module dynamics are
also analyzed.
Chapter 4 describes the development of the MMC-HVDC model under study. The inner control strategies
for MMC are discussed. The control loops employed in VSC-HVDC applications are presented with the
considerations required for MMC converters.
The operation of the system under asymmetric grid conditions is assessed in Chapter 5. An analysis of the
propagation of unbalanced faults in a MMC-HVDC system with a transformerless connection is performed.
A control structure is proposed to eliminate fault reflections in the DC-link considering the three sequence
components. An evaluation of the power controllability of the converter is presented. The performance of
the MMC with current limitation is also evaluated.
The conclusions about the project are drawn in Chapter 6.
Background
2 Background
In this chapter the most common multilevel converter topologies are reviewed. Particular concentration is
addressed in the Modular Multilevel Converter, central for this project. Several modulation strategies
applicable for MMC are also reviewed. Then, an overview of HVDC systems is presented. Finally, the grid
codes applicable for HVDC interconnections are discussed.
2.1.1
The NPC voltage-source converter was initially proposed as a three-level inverter. It is a modification of the
two-level converter topology having two extra power semiconductor switches per phase leg as shown in
Figure 2.1. The midpoint of the switches is connected to the neutral point of the converter through
clamping diodes, enabling the generation of the zero voltage level. By this means, for the same DC-link
voltage, the voltage level that the devices have to withstand is reduced to half comparing to the two-level
topology.
Following the same philosophy, the NPC converter can be extended to more than three voltage levels.
However this topology has several drawbacks. Under certain operating conditions the NPC may experience
capacitor voltage unbalances, creating a potential between the neutral point and ground and causing
distorted output waveforms. This implies the necessity of neutral point or a capacitor balancing control
which is a challenging task when the number of output voltage levels is above three. Moreover there is a
quadratic relation between the required voltage blocking rate of the clamping diodes and the number of
converter levels, which hinders the implementation with a high number of levels. [12]
2.1.2
The topology of the FC converter is presented in Figure 2.2. Each capacitor in the phase is charged to a
different voltage level, therefore by changing the states of the switches, various output voltage levels can
be obtained. [30]
Vdc /2
Vca1
Vdc
Vca2
Vcb1
Vcb2
Vcc1
Vcc2
Vdc /2
Va
Vb
Vc
This topology can have phase redundant switching states that can be used for capacitor voltage regulation,
showing an advantage over the NPC topology. Thanks to the energy storage in the capacitors, the converter
can ride through short duration outages and deep voltage sags. As a drawback, the pre-charge of capacitors
before the start-up, also known as initialisation, is required. Also this topology presents unequal duty
distribution between the switches. Even though the FC topology can be extended to an arbitrary number of
cells, the addition of capacitors leads to an increase in cost and footprint; therefore the number of levels is
usually limited to fours. [12][31]
2.1.3
The CHB topology is based on the series connection of single-phase full-bridge inverter cells with isolated
DC supplies (Figure 2.3). The absence of clamping diodes or flying capacitors, as in case of NPC and FC,
results in the use of minimum components to produce the desired voltage levels. Each inverter cell can
Background
generate three voltages, i.e. both polarities of the DC supply voltage and zero. The output phase voltage is
the result of superimposing the voltages generated by all the cells in the leg. Because the DC sources are
usually supplied by multi-pulse secondary windings of the input transformer, there are no balancing or
initialisation problems as with NPC or FC. [12][32]
The main advantages of the CHB over the NPC and FC are its modular structure and the possibility to have
an independent control over the zero-sequence component in the current. In case of rectifier applications,
the need of many isolated DC sources in series limits the number of cells in the leg, keeping this topology
unfavourable for bidirectional power applications [12][33]. However, a proposal for CHB in HVDC
applications using a reinjection circuit can be found in [31].
2.1.4
The MMC topology is based on a series connection of identical elements, called sub-modules or cells. Each
sub-module represents the basic component of the MMC, shown in Figure 2.4-a. The series connection of
sub-modules in one phase is known as leg. The leg is divided into upper and lower arms such that the
number of the sub-modules in each arm is equal. The AC voltage terminal is the common connection point
between both arms. Since the leg capacitors share a common DC-link voltage there is no need of bulky DClink capacitors, as in case of two-level, NPC or FC topologies. Inductors (
) are placed in the arms to
limit transient currents. [34]
Different sub-module topologies can be applicable to the MMC depending on the application (STATCOM,
HVDC, BTB) [34][35][36]. The difference in the cell structure results in different possible voltage levels at
the terminals of the sub-module. However, with the increase of elements, the capacitor balancing becomes
more complicated. According to the experimental studies performed in [35] evaluating the capacitor
balance and switching losses, the half-bridge topology is the most favourable topology to be implemented
in the sub-modules when bidirectional power conversion is required.
In this project, the term sub-module refers to a half-bridge formed by two bidirectional switches with antiparallel diodes and a DC capacitor, as shown in Figure 2.4-b. The capacitor acts as an energy buffer and a
voltage source. The switches execute the insertion of the sub-module into the arm circuit while the antiparallel diodes ensure uninterruptable current flow.
Since all the sub-modules are identical, the operation principle of MMC can be resumed to the cell level
operation. Each sub-module has two states depending on the switch positions. When the switch S1 in
Figure 2.4-b is ON and the switch S2 is OFF, the sub-module is inserted into the circuit. The voltage between
the terminals
is equal to the capacitor voltage . When the lower switch is ON and the upper is OFF
the sub-module is bypassed and the terminal voltage is zero. As it can be derived from the sub-module
topology, the switches have to operate in complementary way in order not to short circuit the capacitor. By
controlling the number of the sub-modules inserted and bypassed, a staircase output voltage can be
obtained at the AC terminals of the converter.
The direction of the arm current affects the capacitor voltage profile. In Figure 2.5 is shown the current
flow in the sub-module for different states. The assumed positive direction of the arm current
is
represented with red colour and the negative with blue. When the sub-module is inserted, the positive
current will charge the capacitor, passing through the upper diode (a) whereas the negative current will
discharge the capacitor (b). When the sub-module is bypassed the capacitor voltage remains constant.
Background
ON
OFF
ix
ON
ix
OFF
ON
a)
OFF
ix
OFF
ix
ON
c)
b)
d)
Figure 2.5 Positive and negative current flow in a sub-module with different switching states
In Table 2.1 the sub-module terminal voltages and capacitor status depending on the switching states and
the direction of the arm current are summarized. The condition when both switches are off can be used for
the initial charging of the cell capacitors.
Table 2.1 - Switch States of a Sub-module
Switch state
SM terminal voltage
Arm current
polarity
Status of the
capacitor
(+)
Charging
(+)
By-passed
(-)
Discharging
(-)
By-passed
S1
S2
Capacitor Shorted
Open Circuit
0
0
Modularity. The converter can be easily scaled in terms of power or voltage ratings;
Increased output quality. Because the converter can be easily scaled to a large number of the submodules, nearly sinusoidal output can be obtained;
Reliability. Use of redundant sub-modules in the case of a cell failure;
Increased efficiency. Due to low switching frequency of each sub-module;
Reduced footprint. Due to significant reduction or even elimination of the AC filters and no use of
bulky DC-link capacitors.
Due to its possibility to be scaled to high voltage levels, achievable efficiency, ease of implementation, low
harmonics output and reliability, the MMC proves to be the most suitable topology for modern HVDC
applications. [16]
10
11
multilevel modulation methods can be split into two main categories: Space Vector Modulation (SVM) and
Voltage level Based Modulation; i.e. Carrier PWM (CPWM) and Nearest Level Modulation (NLM) [12].
2.2.1
The Space Vector Modulation theory is well established nowadays. Due to its advantages, such as easy
digital implementation and the possibility of optimizing the switching sequences, it is an attractive
modulation technique for multilevel converters. The principle applied for the calculation of the voltage
vectors in two or three level converters can be extended to multilevel converters. However, the complexity
of the algorithms for the calculation of the state vectors and computational costs increase with the number
of levels. Recent publications have presented strategies where simpler algorithms are used; accordingly the
computational efforts are significantly reduced, comparing with conventional SVM techniques. [39][40][41]
The space vector plane for a N-level converter is shown in Figure 2.6. Each connection point on the plane
represents a specific state of the three-phase voltages of the converter. The point (2, 1, 0), for example,
means that with respect to ground, phase A is at 2 , phase B is at 1 , and phase C is at 0; where is the
voltage of the DC capacitor in one sub-module.
The phase voltages of the converter can be represented in matrix form by the switching states of the sub) and the voltages of the DC capacitors. For converters with N+1 voltage levels, the following
modules (
equations applies,
_
,
=
_
(2.1)
,
(2.2)
(2.3)
Background
(2.4)
Because some output voltages can be generated by several switching combinations, redundant switching
states are possible. The number of possible switching combinations is equal to the cube of the converter
voltage levels (N3) while the number of unique states can be calculated by [42]:
!"#
= $ % '$ 1(%
(2.5)
Following a similar approach as with two-level converters, the reference vector for N-level converter is
obtained from the three nearest stationary vectors, which form the vertices of a triangle in which the
reference vector lies. The dwell times of the stationary vectors should satisfy the equation:
)
*+
,- = *+ , . , + *+ , . , + *+%, .%, 0
,- = . , + . , + .%,
(2.6)
In [40] is derived a transformation matrix for the calculation of the dwell times of a N-level converter,
having the times for a two-level. Further improvements regarding multilevel SVM can be found in [43].
2.2.2
Multi-carrier PWM
12
13
Figure 2.7 - Level shifted PMW carriers. (a) Phase Disposition (PD) (b) phase opposition disposition (POD) (c)
alternate phase opposition disposition (APOD)
The LS-PWM methods produce an unequal duty and power distribution among the sub-modules since the
vertical shifts relate each carrier and output level to a particular cell [12]. These can be corrected by
implementing carrier rotation and signal distribution techniques [44].
The Carrier phase shifted method (PS-PWM) has N-1 carrier signals with the same amplitude and
frequency. To achieve a staircase multilevel output waveform, the phase shift between the carriers is
calculated as = 3605 /'N 1( [12]. The multicarrier PS-PWM process is shown in Figure 2.8, the
This approach provides equal duty and power distribution between the cells and, by selecting an adequate
carrier frequency, capacitor voltage balancing can be achieved.
A comprehensive analysis of the Multicarrier PWM techniques was performed in [44], where the
mentioned methods were extended and analysed particularly for MMC applications. The harmonic
distortion of the generated waveforms and the possibility of sub-module capacitor voltage balancing were
the main assessment criteria. It was concluded that the Carrier phase shifted PWM is more suitable for
control of MMC.
Background
2.2.3
The Nearest Level Modulation strategy has been proposed for converters with an arbitrary number of
voltage levels [45][46]. The main idea lies in deciding the number of cells to be inserted and bypassed
'.( with the voltage steps that represent idealised
based on the comparison of the modulating signal
cell capacitor voltages.
For MMC, assuming that the cell voltages are constant, '.( =
one of the $ + 1 discrete voltage levels (0, 8 /$, 2 8 /$,,
inserted and bypassed can be calculated as
9 ,!
9 ,C
'.(
@A ,
1
= :;< = >$ ?
2
1
= :;< = >$ ? +
2
'.(
@A ,
9BB,!
9BB,C
=$
=$
9 ,!
9 ,C
(2.7)
By inserting or bypassing the cells according to (2.7), the average of the generated output voltage matches
the reference voltage, as shown in Figure 2.9. This modulation strategy is suitable for converters with a
large numbers of cells due to the small voltage steps and fundamental switching frequency.
Having a low number of voltage levels, the harmonic generation can be improved by modulating one cell in
each arm. When the modulating signal lies between two adjacent voltage levels D '.( <
'.( <
'D + 1( '.( with D = 0 '$ 1(; D cells have to be selected as ON to provide the base voltage and
one cell should be pulsewidth-modulated to generate the voltage remaining. The number of cells ON to
provide the base voltage can be calculated as:
9 ,!
9 ,C
1
= FG;;: >$ ?
2
1
= FG;;: >$ ? +
2
'.(
@A ,
'.(
@A ,
9BB,!
9BB,C
=$'
=$'
9 ,!
9 ,C
+ 1(
+ 1(
(2.8)
14
15
The modulating signal to be used for the PWM cell can be obtained as,
H
,!
,C
1
= $?
2
1
= $? +
2
'.(
@
'.(
@
9 ,!
9 ,C
(2.9)
In Figure 2.10 the generated arm voltage waveform is presented. In black is marked the reference voltage
for the arm [
'.(], in red is shown the base voltage calculated by (2.8), the PWM voltage is marked
_
blue.
Figure 2.10 Nearest Level Modulation, arm voltage waveform with SM modulation
Traditional CSCs with mercury-arc valves were used since 1950s, until they were substituted by thyristors in
the mid-1970s. Thanks to the rapid development of self-commutated devices and micro controllers, an
Background
alternative as VSC became economically feasible; resulting in the first VSC-HVDC project installed in 1997.
Both converter technologies have different operational principles as well as advantages and drawbacks.
The decision of which option to select depends on the requirements of particular project. [47][49]
2.3.1
Depending on functional aspects, three main HVDC configurations shown in Figure 2.11 are used. [47][48]
Figure 2.11 - HVDC system configurations. (a) Monopolar. (b) Bipolar. (c) Back-to-back
2.3.2
Monopolar configuration (a) - interconnects two converter stations via a single line, with the
possibility to operate at both DC polarities. Ground, sea or metallic conductor can be used for
return path.
Bipolar configuration (b) - involves two conductors, operating at opposite polarities. This results in
two independent DC circuits, rated at half capacity each. During outages of one pole, a monopolar
operation can be used. This is the most common configuration for modern HVDC transmission.
In Back-to-Back configuration (c) - the DC sides of two converters are directly connected, having no
DC transmission line. This arrangement is used for the interconnection of asynchronous AC
systems.
VSC-HVDC transmission
Even though traditional CSC-HVDC transmission is well established for high power and voltage ratings
(typically up to several GW and 800 kV), it is predicted, that from now on the VSCs will be dominant in the
future high power HVDC interconnections due to numerous advantages in economic and technical features
[47][49]. The main advantages of VSC-HVDC over CSC-HVDC are summarized below [48][49][50]:
Independent reactive power control at the both terminals, possibility of four quadrant operation
(Figure 2.12). The elimination of reactive power compensation devices results in significant
footprint reduction;
Dynamic support of the AC grid voltage. Operation as STATCOM increases transfer capability and
stability of the AC grid;
Possibility of connection to the weak and passive grids. Low short-circuit capacity requirements of
the AC grid. Since a VSC can be considered as a virtual synchronous generator, it can be used for
forming offshore AC collector systems for wind power parks;
Possibility of safe fault ride-through and black start capability;
16
17
The typical configuration of modern VSC-HVDC transmission system is shown in Figure 2.13. Two DC
conductors of opposite polarity interconnect two converter stations. The polarity of the DC-link voltage
remains the same while the DC current is reversed when the direction of the power transfer has to be
changed. The DC side capacitors ensure support and filtering of the DC voltage. The converter AC terminals
are connected with phase reactors and harmonic filters. The phase reactors ensure control of power
exchange between the converter and AC system, the limitation of fault currents and blocking of current
harmonics appearing due to PWM. The AC filters reduce harmonics content on the AC bus voltage. Power
transformers are used to interface the AC system, adapting converter and AC system voltages as well as
participate in power regulation by means of tap changers. [51][52]
Theoretically all the multilevel topologies presented in this chapter can be used in VSC-HVDC
configurations. However, due to the complex structure, voltage balancing issues and economical
considerations, most of the real life applications of VSC-HVDC systems rely on the proven two-level and
three-level NPC converter technologies. [16][49]
Background
With the introduction of MMC, the application areas of VSC-HVDC transmission can be broadened
significantly. Due to the numerous advantages such as modularity, increased efficiency and reliability that
MMC presents, it aims to substitute the existing VSC-HVDC topologies in the nearest future. [16]
The configuration of MMC-HVDC transmission system is shown in Figure 2.14. Compared to the topology
presented in Figure 2.13, it can be noticed, that depending on the number of voltage levels and quality
requirements of output voltage; AC filters can be significantly reduced or eliminated. Transformers become
also optional, since the converter can be scaled to meet the voltage levels of the transmission systems. Due
to distributed energy storage in the leg sub-modules, the DC capacitors are also eliminated. [37]
The mentioned technical aspects result in reduced complexity and footprint of the converter station,
making it especially suitable for offshore platforms, where size and reliability are the major selection
criteria. [17]
Operation requirements
Since the AC grid is a dynamic system, its main parameters such as voltage at different network points and
frequency are subjected to variation due to change in power balance and structure of the grid. These
variations to some extent should not affect the operation of the HVDC transmission. The minimum
operation capabilities are defined in [53].
18
19
Operate continuously at constant active power output at transmission system frequencies in the
range 49.5Hz to 50.5Hz;
Operate and remain connected to the transmission system at frequencies within the range 47.5Hz
to 52.0Hz;
Remain connected to the transmission system at frequencies within the range 47.0Hz to 47.5Hz for
a duration of 30s required each time the frequency is below 47.5Hz;
Remain synchronised to the transmission system during rate of frequency change up to 1 Hz/s;
Remain connected providing constant active power output at transmission system voltage
variations within the ranges of 10% of nominal.
The relation of the converter active power and the grid frequency for both inverter and rectifier modes is
shown in Figure 2.15. It can be observed, that in inverter mode, when grid frequency is decreasing, active
power output can be reduced only by 5% of rated. In contrary, active power input has to be reduced up to
40% when rectifying. This restriction provides frequency support when unbalanced generation-load
conditions in the AC system.
2.4.2
Remain connected to the transmission system during a negative phase sequence load unbalance;
Be capable of reversing the power at a rated capacity within 5s when emergency;
Provide ramp-up and ramp-down capability not less than the greater of 10% of the rated capacity
per minute or 50 MW per minute.
V/Q Control requirements
The VSC-HVDC interconnection is required to have full control over the reactive power at any point
between the 0.95 lagging and 0.95 leading power factor. The reactive power limits in accordance to the
active power output are specified in Figure 2.16-a. At steady-state operation, the converter should also
provide continuous voltage control at the PCC. Set points and slope characteristics are detailed in Figure
2.16-b.
Background
Figure 2.16 (a) Reactive power requirements. (b) Voltage control requirements. [53]
The reactive power response should start within 0.2s of the application of the step. Full reactive power
capability should be reached within 1s.
2.4.3
Fault ride-through capability is defined as the ability of the power converter to withstand different types of
faults. According to [53], the converter must remain stable and connected to the system without tripping
for a close-up solid three-phase or any unbalanced short circuit fault for up to 140 ms as shown in Figure
2.17. Each point on the line represents grid voltage level and the associated time duration which converter
must stay connected. When voltage level is below the line, converter disconnection is allowed.
Upon clearance of the fault and within 0.5 seconds of the restoration of the voltage to the 90%, active
power output should be restored to at least 90% of the pre-fault level. The fault will affect the level of
transferred active power, therefore a load reduction or rejection in the other side of HVDC transmission is
acceptable. With active power reduction, the converter should generate maximum reactive current without
exceeding the transient rating limits of the device, thus supporting the grid voltage during the fault.
Moreover, presence of the fault should not be indicated in the other side of the HVDC transmission. [53]
20
21
H<
=N
<
JLFF <
J<
= FF
HG
MN
JLFF G
JG
Under these considerations it is possible to represent the insertion of the arm capacitors with a continuous
value (H!/C ); going from 0, when all the sub-modules are bypassed, to 1, when all the sub-modules are
inserted. If the sum of the voltages in the arm capacitors is
represented as:
'.( = H '.(
(3.1)
J
$ H '.(
(3.2)
When the arm current '.( is flowing through the effective capacitance, the total capacitor voltage
dynamics are described by:
=
=.
'.(
(3.3)
Having defined the direction of the arm currents according to Figure 3.1, the output phase current
calculated by Kirchhoff's current law (KCL) as:
!
(3.4)
8#
8#
8#
+
2
(3.5)
(3.6)
(3.7)
The expression of capacitor voltage dynamics from (3.2) and (3.3) can be expanded for the upper and lower
arms:
=
=.
=
=.
O
!
O
C
$H! !
J
$HC C
=
J
(3.8)
(3.9)
Analyzing the circuit given in Figure 3.1 it is possible to derive expression for the generated AC voltage (LR )
as:
LR =
LR =
+I
=!
H!
=.
=C
+ HC
=.
O
!
O
C
(3.10)
(3.11)
22
23
Subtracting (3.10) from (3.11) the expression for the difference current can be obtained as follows:
0=
2I
8#
8#
=.
HC
O
C
+ H!
O
!
(3.12)
From (3.10), (3.11) and substituting for the currents from (3.7) in (3.12) the phase leg of the converter can
be described by the following system of differential equations:
I
W
V
8#
=
V $H!
S O! T = V
J
=.
O
V
C
V $HC
U J
H!
0
0
HC
Z
Y
Y
YS
Y
Y
X
W
V 2
8#
V$H!
O
! T + V 2J
O
V $H
C
C
V
2J
U
8
Z
Y
Y
Y
Y
Y
X
(3.13)
8#
time, 8# only depends on the DC-link voltage, 8 , and the total inserted arm voltage, H [ . Therefore,
the difference current, 8# , can be influenced without disturbing the AC side quantities by adjusting the
lower and upper arm insertion index in the same amount [22].
It can be noticed from (3.10) and (3.11) that the output voltage does not depend on
. At the same
The selection of the cell capacitance is a trade-off between the voltage requirements of the sub-module
and the capacitor size. In [57] the total cell capacitance is proposed to be 30-40 kJ per MVA of the
converter, resulting in 10% voltage ripple. For the simulation model cell capacitance is calculated as 220 F
for 40 kJ/MVA. However, in [11] is presented an analytical expression that can be used for the capacitance
calculation based on the desired ripple factor.
3.2.2
Inductors are placed in the converter arms to suppress transients in the circulating and fault currents. For
the simulation model the arm inductance is selected to be 0.15 p.u. on converter base which is a
conventional value used for HVDC projects [55]. Nevertheless, several methods have been proposed to
calculate the inductance based on the desired circulating current amplitude or the limits of fault currents,
as presented in [19].
3.2.3
Model Parameters
A switching model of MMC is implemented in PSCAD/EMTDC software with the circuit parameters
summarised in Table 3.1. The semiconductors have 0.01 ON-state resistance.
Table 3.1 Circuit parameters used for MMC model 8-SM
Description
Rated apparent power
Rated cell voltage
DC link voltage
Arm resistance
Arm inductance
Cell capacitance
Number of sub-modules per arm
Abbreviation
]
Value
850 MVA
80 kV
320 kV
0.1
70 mH (0.15 p.u.)
220 F (40 kJ/MVA)
8
Additionally, the PS-PWM method with carrier frequency F = 301 ^ is implemented in the switching
model (the switching frequency is selected to have a natural balancing of the sub-module capacitors). Both
models are subjected to direct modulation with modulation index 0.9. The load at the AC terminal is
I_ = '180 + a0.05( c, corresponding to the nominal load of the converter system.
3.3.1
In Figure 3.2 the response of the switching and average models is shown with black and red lines
respectively. The step is applied at the instant t=1.4s. It can be noticed that the capacitor voltages of both
models are close during steady operation having error of 3%. After the transient the DC value of the sum
capacitor voltages is changed from 640 kV to 512 kV which represents the 20% reduction in DC-side
voltage. The relative error between both models is increased to 8% during the transient. The difference
currents have an average error within 5%.
24
800
kV
600
400
800
kV
600
400
(a)
4.0
2.0
kA
0.0
-2.0
-4.0
1.350
1.400
1.450
1.500
1.550
3.3.2
1.600
[
!,
(b)
Figure 3.2 Average vs Switching model, response to DC step. (a)
[
C.
1.650
In Figure 3.3 the response of both models to the AC reference change is shown. The step is applied at the
instant t=1.4s. The sum of the capacitor voltages of both models are close during steady and transient
having relative error within 4%. After the step is applied the ripple is reduced as a consequence of change
in the power output. The difference currents have an average error within 5%.
kV
800
600
kV
800
600
(a)
2.0
0.0
kA
25
-2.0
-4.0
1.350
1.400
1.450
1.500
1.550
(b)
Figure 3.3 - Average vs Switching model, response to AC step. (a)
[
!,
1.600
[
C.
1.650
It can be concluded that the average model represents the dynamics of the MMC system acceptably,
having mean error of 4% in the simulations performed, with a maximum error of 8% during transients. The
difference between the switching model and the average model is proportional to the number of submodules and the switching frequency, thus in simulations with a higher number of sub-modules the
accuracy should be significantly improved. The average model can be used for preliminary design of the
outer control loops and steady-state analysis of the MMC system. However, with the average model, it is
not possible to verify control strategies implemented at cell level.
LR =
+ LR =
LR N;f 'g.(
(3.14)
+ LR N;f 'g.(
(3.15)
where LR is the peak value of the generated emf. Considering a pure DC difference current, in a 1-phase
system, it is possible to substitute 8# in (3.5) and (3.6) by the DC current 8 as follows,
!
C
h
N;f 'g. + i(
2
h
N;f 'g. + i(
2
+
(3.16)
(3.17)
where h is the peak value of the AC current. Using (3.14) (3.17), the power in the upper and lower arms
is calculated as,
j! =
! d!
jC = C dC
8
8
h LR
N;f'i(
4
h LR
N;f'i( +
4
8
8
LR N; f'g.( +
LR N; f'g.(
N; f'g. + i(
N; f'g. + i(
h LR
N;f'2g. + i(
4
h LR
N;f'2g. + i(
4
(3.18)
(3.19)
Adding the power in the upper and lower arms to calculate the power in the phase-leg results in,
jC
= j! + jC =
h LR
h LR
N;f'i(
N;f'2g. + i(
2
2
(3.20)
where the first two terms represent the DC power and the third term is a double line-frequency AC
component present in the arm. Integrating (3.20), the AC component of the energy stored in the phase-leg
is given by,
mC
l,
= n jC l =. =
h LR
f '2g. + i(
4g
(3.21)
26
27
The energy stored in the phase is located in the sub-module capacitors. The energy stored in the capacitors
has a quadratic relation with the capacitor voltage, thus from (3.21) can be derived that the voltage across
the phase-leg must contain a double line-frequency AC component and a DC component.
In a three-phase system the DC components of the leg voltages are equal, and the double frequency
components are shifted by 2o3 in a-c-b sequence [19]. Including the double frequency component into
the arm voltage equations presented previously (3.14)(3.15),
d
f '2g. + i(
2
2
d
8
dC =
+ LR N; f'g.(
f '2g. + i(
2
2
d! =
LR N; f'g.(
(3.22)
(3.23)
The double line-frequency voltage excites the leg impedance, creating a double line-frequency component
on the difference current also called circulating current # ,
=
d
4q
f '2g. + i M:r'q
(( = h
8#
8#
+ h
+ h
N;f '2g. + s(
h
N;f 'g. + i(
2
h
N; f'2g. + s(
N;f 'g. + i(
2
N; f'2g. + s( +
(3.24)
.
(3.25)
(3.26)
Analyzing (3.25) and (3.26) can be concluded that the circulating current has no effect on the AC current.
This current flows through the converter phase-leg, increasing the RMS value of the arm currents, resulting
in higher converter losses. From (3.20) can be noticed that the DC component of the difference current is
responsible for the DC/AC power transfer.
S2
J]t
Q
]t
S1
The arm current interacts with the capacitors of the sub-modules by means of their switching function (]# ),
defined as follows,
1, the sub module is inserted 'S1 = ON, S2 = OFF( 0
]# '.( = u
0, the sub module is bypassed 'S1 = OFF, S2 = ON(
(3.27)
= ]# '.( '.(
(3.28)
During normal operation, due to the power transfer the capacitors get charged and discharged, imposing a
ripple component (d _# ) on the capacitor DC voltage ( _8 ). In general, the z-th harmonic of the capacitor
voltage ripple (d
'(
_# )
'(
_# '.(
ag J
(3.29)
The volume of the capacitor is inversely proportional to the amplitude of the voltage ripple; therefore the
selection of the capacitors is a trade off between the acceptable ripple and the capacitor size. The voltage
in the sub-module capacitor can be generally expressed as:
_# '.(
_8
+ d
'(
_# '.(
(3.30)
The capacitor voltage is reflected in the terminal voltage of the sub-module by means of its switching
function as follows:
_# '.(
= ]# '.(
_8
+ d
'(
_# '.(
(3.31)
28
29
_# '.(
= ]# '.( d
'(
_# '.(
(3.32)
Subsequently, the voltages generated by the upper and lower arms of the converter can be expressed as:
! '.(
C '.(
= ]!# '.(
#
= ]C# '.(
_8
+ d
#
(3.33)
_C_# '.(
(3.34)
(3.35)
_8
+ d
_!_# '.(
The sum of the sub-module switching functions in the arm is called insertion index H '.(, presented in
Section 3.1. In ideal conditions, at any time instant the sum of the inserted sub-modules in the upper and
lower arms should be equal to the number of the sub-modules in the arm ($). That implies the following
relation between the switching functions of the arms:
Since each leg of the converter is composed of two arms, the voltage across the leg can be obtained by
adding the voltages generated by the arms. Adding (3.33) to (3.34) and using relation from (3.35) the
expression for the leg voltage can be obtained,
where dC
C l '.(
=$
_8
+ dC l '.(
(3.36)
'(
_!_# '.(
+ d
#
'(
_C_# '.(
(3.37)
Having the DC-link voltage of the converter formed by C l and the voltage drop across the arm impedance;
from (3.36) and (3.37) can be concluded that, by reducing the ripples in the capacitors, the DC-link voltage
ripples can be reduced.
3.5.1
Assuming a large number of sub-modules or high switching frequency, the higher order components can be
eliminated and the upper and lower arm insertion index of one phase become continuous, given as:
1
1 1
H'.( = tN;f'g.(
2
2 2
1
1 1
HC = + H'.( = + tN;f'g.(
2
2 2
H! =
(3.38)
(3.39)
!_
1
1
1
1
1
th N;f'i( t 8 N;f'g.( + h N;f'g. + i( t # N;f'g. + s(
8
2
4
4
2 8
R '.( =
F< =MHL .MG N;H; L .
=N N;H; L .
1
1
1
th N;f'2g. + i( +
t # N;f'3g. + s(
# N;f'2g. + s(
8
2
4
2"8 M:H; N N;H; L .
C_
(3.40)
3 8 M:H; N N;H; L .
1
1
1
1
1
th N;f'i(+ t 8 N;f'g.( + h N;f'g. + i( + t # N;f'g. + s(
8
2
4
4
2 8
R '.( =
F< =MHL .MG N;H; L .
=N N;H; L .
1
1
1
N;f'2g. + s( + t # N;f'3g. + s(
th N;f'2g. + i( +
#
2
4
8
2"8 M:H; N N;H; L .
(3.41)
3 8 M:H; N N;H; L .
Multiplying the fundamental term of the arm currents with the capacitor reactance, the fundamental
voltage ripple for both arm capacitors can be obtained:
d
d
' (
! '.(
' (
C '.(
=
=+
t8
2gJ
N;f'g.( +
t8
2gJ
N;f'g.( +
h
4gJ
h
4gJ
N;f'g. + i(
N;f'g. + i( +
t #
4gJ
t #
4gJ
N;f'g. + s(
N;f'g. + s(
(3.42)
(3.43)
In similar way the 2nd and 3rd harmonic components of the arm capacitor ripples can be derived:
d
d
' (
! '.(
' (
C '.(
=
=
th
16gJ
th
16gJ
'%(
! '.(
N;f'2g. + i( +
N;f'2g. + i( +
'%(
C '.(
t #
12gJ
t #
12gJ
4gJ
#
4gJ
N;f'2g. + s(
N;f'2g. + s(
N;f'3g. + s(
N;f'3g. + s(
(3.44)
(3.45)
(3.46)
(3.47)
From (3.44)-(3.47) the relation between upper and lower arm capacitor ripples can be noticed. The
corresponding components are equal in magnitude but have different signs as shown in (3.48).
30
31
(3.48)
Having determined the arm capacitor voltage ripples as the result of the average arm currents, the voltages
across the upper and lower arms can be expressed from (3.33), (3.34) and using arm switching functions
from (3.38) and (3.39):
where d
! '.(
= H!
C '.( = HC
_8
_8
+ H! d ! '.(0
+ HC d C '.(
(3.49)
and d C represent the total voltage ripple on the arm capacitors and are described as:
' (
' (
'%(
d ! '.( = d ! '.( + d ! '.( + d ! '.(0
' (
' (
'%(
d C '.( = d C '.( + d C '.( + d C '.(
(3.50)
By substitution of (3.50) into (3.49) using the relationship given in (3.48) the arm voltages can be
represented as follows:
! '.(
C '.(
_8
_8
H'.(
+ H'.(
_8
_8
d8 '.( + d
+ d8 '.( + d
d8 '.( =
'.( =
'.(
'.(
(3.51)
(3.52)
$
' (
' (
'%(
d C '.( + d C '.( tN;f'g.( + d C '.(
2
$
' (
' (
'%(
d C '.( tN;f'g.( + d C '.( + d C '.( tN;f'g.(
2
(3.53)
(3.54)
From (3.51) and (3.52) it becomes clear, that converter arm voltages during steady operation are composed
of several frequency components. The second harmonic of d generate the circulating current observed
in the previous section. Since the third harmonic in all three phases is equal, the corresponding ripple of
d generates a zero sequence current that flows to the DC-link. The differential mode voltage ripple,
d8 , produces a current component that flows to the AC side.
In this chapter the average model of the MMC was derived and compared with a 9-level switching model.
Analyzing their response to a step change, the average model demonstrated an accurate representation of
the converter dynamics.
The converter arm currents were studied, demonstrating the presence of a double frequency AC
component product of the DC/AC power transfer. It was observed that this component does not contribute
to the power transfer; therefore it only increases the RMS value of the arm current, increasing the
converter losses.
The effect of the sub-module dynamics on the converter was mathematically analyzed, showing the
influence of the capacitor current harmonics on the ripples of the converter arm voltages.
32
33
At the end of this section, a control strategy or combination of strategies is selected for further
implementation into the HVDC transmission system.
4.1.1
Energy Control
This method was first proposed by Antonopoulos in [22], where the arm capacitor voltages are kept to a
reference through the control of the total stored energy m [ in the phase leg and the difference between
the energy stored in the upper and lower arms m .
Following the closed loop approach proposed in [22]; aiming to compensate for the capacitor ripples, the
ideal insertion indices can be calculated as,
H! =
HC =
2 LR
2 + LR
[
!
[
C
=I
d8#
d8#
8#
(4.1)
(4.2)
8#
8
when using this selection of insertion indexes two additional control loops must be implemented in order
to gain stability.
An open loop approach using the estimation of the stored energy is proposed in [23] with the intention of
increasing the stability of the system and avoiding the need of a continuous measurement of the capacitors
voltage to calculate the converter stored energy.
In this approach the measurement of the sum of the voltages in the arm capacitors (
an estimated value ( [!,C ). The insertion indices are calculated as follows:
H! =
HC =
2 LR d8#
[!
2 + LR d8#
[C
[
!,C )
is replaced with
(4.3)
(4.4)
The estimated values are obtained by integrating (3.18) and (3.19), and considering a constant difference
current 8# and difference voltage d8# .
The obtained insertion indexes compensate for the ripples in the capacitor voltages, as the result
circulating current is suppressed.
A block diagram of this control strategy is presented in Figure 4.1. The nearest level modulation is used
(Section 2.2.3) and the sub-module balance is done by means of sorting algorithm.
34
=N /2
=N
MN
N<
m
Ld
:LF
d= FF
=N /2
NG
m
N<
NG
H<
HG
Gating pulses
d= FF
Gating pulses
kV
90
80
70
(a)
100
kV
90
80
70
(b)
1.50
1.00
0.50
kA
35
0.00
-0.50
70
0.90
1.00
1.10
1.20
1.30
1.40
1.50
(c)
Figure 4.2 Open-Loop Energy Control, simulation results. (a) Upper arm capacitor voltages. (b) Lower arm
capacitor voltages. (c) Difference current
The method presented good response to a change in the arm energy reference and effective elimination of
the circulating current. The capacitor voltages remain well balanced in steady state and transients, which is
achieved by the sorting algorithm.
However, implementing this control strategy, it was observed that distortions appear in the converter
output currents. In [61] a method is proposed to eliminate this problem, but due to time constrains it was
not further investigated.
4.1.2
Distributed control
In this method, proposed by Akagi in [20], the cell capacitor voltages are controlled independently. The
control is implemented in two parts:
-
In Figure 4.3 the block diagram of the distributed control method is presented. As it can be observed, the
averaging control is implemented in two loops, outer voltage loop and inner current loop. The voltage loop
is responsible of controlling the mean value of the capacitor voltages in the leg by influencing each cell
individually. The error signal is processed in the controller, resulting in the reference signal for the
difference current loop. Under the balanced conditions, the DC component of the difference current is
equal to 1/3 of the DC-link current, therefore a feed-forward term is added to increase the response of the
controller as highlighted in Figure 4.3.
If the average voltage , R is lower than the desired value,
, a positive current reference is obtained.
The current reference is subtracted from the measured value, reducing the control command. By this
means, the DC component of the difference current is increased, rising the charge in the capacitors.
:LF
= FF
J,Md
:LF
J
=N /3
= FF
M:H
Ld
:LF
The average charge of the capacitors depends on the DC component of the difference current. If only an
integral compensator is used, the DC value of the difference current is controlled, making in no effect on
36
the circulating current. The compensator in the current loop acts on the AC component of the difference
current. The Balancing control is implemented in each sub-module individually. The control signal is
generated based on the capacitor voltage and the direction of the corresponding arm current.
The final sub-module voltage reference is obtained by adding both averaging and balancing control signals
to the voltage reference.
4.1.2.1 Simulation results
Figure 4.4 depicts the voltages of the phase-leg capacitors with only balancing control. The controls is
activated at t=0.5s, it can be observed that after 0.1s the capacitor voltages are effectively balanced.
90
kV
80
70
-1.00
0.450
0.500
0.550
0.600
0.650
0.700
0.750
In Figure 4.5 the arm capacitor voltages (a,b) and difference current (c) of converter phase A are shown.
The averaging control is activated at t=1s, after the stabilisation the amplitude of the circulating current is
reduced by 35% and the ripple factor in the capacitor voltages is reduced from 9% to 5%. At time t=1.2s a
10% step change in arm voltage reference is done. The new voltage value of 88 kV is reached within 2
fundamental cycles.
100
kV
90
80
70
(a)
100
kV
90
80
70
(b)
1.50
1.00
0.50
kA
37
0.00
70
-0.50
s
0.90
1.00
1.10
1.20
1.30
1.40
1.50
(c)
Figure 4.5 - Distributed Control, simulation results. (a) Upper arm capacitor voltages. (b) Lower arm capacitor
voltages. (c) Difference current
The method has shown a fast response to changes in the voltage reference. However, the circulating
current amplitude was just reduced by 35%, when the objective is to suppress it. The capacitor voltages
were effectively balanced in steady state and during transients.
4.1.3
This method first proposed in [24] aims to eliminate the double line-frequency circulating current. The
circulating current suppressing controller (CCSC) is implemented as vector current control, as shown in
Figure 4.6.
The three phase difference currents are calculated and transformed to a double line-frequency negativesequence rotational reference frame. The current references for the control loops are set to zero since the
objective is to eliminate the circulating current. The output signal is transformed back to three phase
signals, which are subtracted from the upper and lower arm voltage commands of the corresponding
phase. In this manner, the arm currents can be influenced without affecting the output voltage.
<= FF ,=
N :N ,MN
N :N ,=
N :N ,
:LF
<= FF
:LF
:LF
<= FF ,
<= FF ,M
:LF
Ld
:LF
38
kV
90
80
70
(a)
kV
90
80
70
(b)
1.50
1.00
0.50
kA
39
0.00
-0.50
0.900
0.950
1.000
1.050
1.100
1.150
1.200
1.250
(c)
Figure 4.7 - Ciculating Current Suppression Control, simulation results. (a) Upper arm capacitor voltages. (b) Lower
arm capacitor voltages. (c) Difference current
A modification of this method is presented in [62], implementing quasi-PR controllers on each phase to
suppress higher order harmonics of the circulating current. By this means, the harmonic content of the
output voltage can be reduced even more. Because of the limited advantages presented in this method,
adding a more complex structure, the modified CCSC was not implemented and tested.
Due to the satisfactory results shown in simulations and simple implementation, the Direct suppression
strategy together with the Distributed control were selected for further integration into an HVDC
transmission system.
The current loop is responsible for fast tracking of references generated in the power controller, DC or AC
voltage controllers. When operating in inverter mode, the converter controls the DC-link voltage at
predefined value. To achieve this, the DC voltage controller adjusts the active current reference in such a
way, that the net imbalance of power exchange between the DC and AC systems is kept to zero. In rectifier
mode, the converter tracks active power references directly. The reactive power at both sides can be
controlled independently. It can be regulated to track a reference, thus regulating power factor at the PCC;
or to control of the AC grid voltage at the PCC [33]. A phase locked loop (PLL) is used for the
synchronisation with the grid voltage. The PLL mechanism is able to detect phase angle and the magnitude
of the grid voltage, to be later used in the controls. The grid frequency can also be obtained from PLL.
In the following sections, the application of the control loops presented in Figure 4.8 is investigated for
MMC-HVDC system.
4.2.1
The synchronous reference frame PLL is the most extended technique used for the synchronisation with
the three-phase grid. In Figure 4.9 is presented the structure of PLL where the three-phase voltages are
transformed into vectors in the dq rotating frame by means of Parks transformation (Appendix). The feedforward loop includes an integrator, which is resettable at g = 2o, to calculate the grid angle and a PI
controller, which adjusts the angle to the objective'.( g. + 5 , thus making grid voltage component
- in the steady-state equal to zero. In that case -8 represents the amplitude of the balanced three-phase
voltage. The obtained angle and angular grid frequency g are used for abc/dq voltage and current
transformations and are fed into the control loops as part of de-coupling terms. [33][31]
f,MN
f=
g0
40
41
f +
,#
f + f +
,#
(4.5)
The structure of 'f( is similar to the general second-order transfer function, therefore PI parameters can
be obtained as follows:
9.2
..-
,# =
2.3
=
(4.6)
(4.7)
Where .- is the desired settling time and is the damping factor. [63]
4.2.2
The dynamics of the converter AC side in dq reference frame are given by [33]:
=
=.
=
=.
= 8g
8I
-8
(4.8)
(4.9)
is
Based on (4.8) and (4.9) the dq current control structure is designed, as shown in Figure 4.10. Since the
tracking signals are DC values, a proportional-integral (PI) controller is used. A voltage feed-forward and
current de-coupling terms are added in order to improve the performance of the controller. [33]
:LF
=
MN
f=
:LF
:LF
.=
:LF
.
:LF
.,=
Since both d and q loops have the same dynamics, the tuning of PI controllers was realised for one loop.
The current control loop is designed to achieve fast response. The tuning is based on the modulus optimum
principle, due to cascaded control structure [64]. Neglecting the coupling and feed-forward terms, the
reduced block diagram of current control loop is shown in Figure 4.11.
:LF
= , +
where ,
PI Controller
, 1 + f, ,
f, ,
System Transfer
Function
Converter delay
:LF
.=,
1
1 + f,Md
.=,
1
I'1 + f(
=,
is the average time delay, which equals to half of the switching cycle.
,C
,# 1 + f,#,#
1
f,#,#
1 + f,
1
I'1 + f(
(4.10)
where = /I is the time constant of the plant. By choosing the controller zero as ,#,# = , the
cancellation of slow system pole is achieved. The open loop transfer function becomes:
,C
,#
1
fI '1 + f, R (
(4.11)
1+
,C
Resulting in:
,# =
,C
=1
I
2, R
(4.12)
(4.13)
, C
1
2, R f + 2, R f + 1
(4.14)
Substituting system values into derived equations, controller parameters are obtained as ,# = 194,
,#,# = 0.81. These values were used as starting point for, to be adjusted in PSCAD simulations. The final
42
kA
43
-1.0
s
0.10
0.20
0.30
0.40
0.50
0.60
4.2.3
0.70
, (green)
0.80
DC voltage Control
The DC voltage control loop adjusts active current reference to keep the net power exchange between the
DC and AC systems to zero. The tuning method is based on the symmetrical optimum criterion [64][65]. The
simplified block diagram of the DC voltage control is shown in Figure 4.13. The inner current controller is
substituted by an equivalent first order approximation, having , = 2, R [64].
:LF
=N
,d 1 + f, ,d
f, ,d
:LF
=
1
1 + f,L
3
2
=N
1
fJ
=N
The DC system transfer function is obtained from the DC-link capacitor current equation, linearised around
the steady-state voltage reference. The open loop transfer function of the system is:
RC,C =
,R 1 + f,#,R
1
3
f,#,R
1 + f, 2
1
fJ
(4.15)
Applying symmetric optimum conditions [64], the controller parameters can be obtained as:
,#,R = M ,
,R =
2 8 J
3 8,
where M is the symmetrical tuning parameter, usually taken within the range of 2-4 [65].
(4.16)
(4.17)
In the MMC topology there is no capacitor between the DC terminals; the DC-link voltage is formed by the
inserted sub-modules in the three phases and the voltage drops across arm inductances. Assuming that at
any time instant the number of inserted sub-modules in one leg is equal to $, the equivalent DC-link
capacitance is:
J8
=3
J
$
(4.18)
where coefficient 3 comes from the number of converter legs. Using the equivalent capacitance and system
parameters, the controller gains are obtained as ,R = 0.025, ,#,R = 0.015. In PSCAD implementation the
gains are = 0.025, is ,# = ,#,R /,R = 0.62. The response of the DC voltage control to a 10% step
change is presented in Figure 4.14. It can be observed that the voltage settles around the reference after
0.1s, having an overshoot of 21%.
650
kV
600
550
500
s
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
4.2.4
The instantaneous power delivered to the AC grid in dq reference frame can be calculated by [33]:
3
j-,8 '.( =
2
-,8 '.( =
-8 '.( 8 '.(
3
2
-8 '.( '.(
(4.19)
(4.20)
where -8, are the grid voltage components and cannot be controlled. The references for the currents can
be obtained from (4.19) and (4.20) as follows [33]:
8
'.( =
'.( =
-8
j-
-8
'.(
- '.(
(4.21)
(4.22)
To ensure precise tracking of power references, the power control is implemented in a closed-loop manner
as shown in Figure 4.15. This implies additional pair of controllers that provide error reference signals for
currents.
44
45
:LF
=
jf
:LF
jf
:LF
:LF
:LF
The controller gains are calculated based on symmetrical optimum criterion and adjusted in the simulation
model to achieve desired response. In PSCAD implementation the gains are = 0.0006, ,# = 0.2.
4.2.5
AC voltage control
Alternatively to the power controller, the reactive power can be controlled indirectly, through the AC
voltage control at the PCC. The AC voltage control loop provides reactive current reference based on the
measured grid voltage as seen in Figure 4.16.
MN
:LF
MN
:LF
The control signal is determined comparing the grid voltage RMS with the reference vale. If the grid voltage
level at the PCC is lowered, the reactive power output of the converter increases, boosting the AC voltage.
The described control loops were implemented in the simulation of HVDC transmission system, shown in
Figure 4.17. The system consists of two MMC stations, interconnected with 100 km long DC cable
(parameters in Appendix). The three-phase MMC with parameters from Table 3.1 are transformerlessly
connected to the AC grids. The grid stiffness is indicated by the short-circuit ratio of 10. The AC line-to-line
RMS voltage is 350 kV resulting in converter operation at 0.9 modulation index, which ensures good
controllability margins.
Both MMC stations have similar control loops as shown in Figure 4.18 having the possibility of changing the
operation mode, i.e. power transfer reversal. The output of the MMC inner controls is added to voltage
reference from the current control loop and sent to the modulator.
:LF
=N
=N
:LF
:LF
=
jf
f=
f=
f,MN
MN
f=
jf:LF
:LF
:LF
=
:LF
<]t
:LF
:LF
.
MN
:LF
MN
:LF
.=
MW
0
-400
-800
(a)
400
MVAr
200
0
s
0.40
0.50
0.60
0.70
0.80
0.90
1.00
1.10
(b)
Figure 4.19 HVDC system response to P,Q ramps. (a) blue - MMC-1 Active power, green - MMC-2
Active power. (b) blue - MMC-1 Reactive power, green - MMC-2 Reactive power.
46
Voltage profile at both sides of the DC transmission is shown in Figure 4.20. Before the power transfer is
activated, DC voltage level at the converter stations is maintained to the reference of 640 kV. After the
active power transmission starts from MMC-2 to the DC-link, the voltage rises by 4% at the rectifier side
(green) and 1.5% at the inverter (blue) during the power ramp. Since DC voltage control is active at MMC-1,
the reference voltage is reached within 0.05s after active power stabilisation. The differences between the
voltages at the two terminals are a consequence of the DC-cable impedance. The noise in the DC voltage
reflects switching of the converter sub-modules.
680
kV
660
640
620
s
0.40
0.50
0.60
0.70
0.80
0.90
1.00
1.10
Figure 4.20 HVDC system response to ramp, (blue) DC-link voltage in MMC-1 terminal, (green) DClink voltage in MMC-2 terminal.
Since operational conditions of both converter stations are similar, only the capacitor voltages of MMC-1
are presented (Figure 4.21). The capacitor voltages present ripple factor of 5% at rated active power, after
t=0.8s the ripple factor is increased to 9% as a result of the reactive power exchange between the
converter and the grid. During the power changes, the capacitor voltages remain balanced and follow the
reference of 80 kV.
P ramp
Q ramp
85.0
kV
80.0
75.0
(c)
85.0
80.0
kV
47
75.0
0.40
0.50
0.60
0.70
0.80
0.90
1.00
1.10
(d)
Figure 4.21 - HVDC system response to ramp. (a) Upper arm capacitor voltages. (b) Lower arm capacitor
voltages.
MW. At time t=1.5s a ramp command for full active power reversal is activated. When power reference is
crossing zero at t=4s, the task of DC voltage control automatically is assigned to MMC-2. During the whole
simulation time the DC-link voltage at both terminals closely follows the reference (b).
0.8k
MW, MVAr
0.4k
0.0
-0.4k
-0.8k
(a)
680
660
kV
640
620
600
72.0
1.0
2.0
3.0
4.0
5.0
6.0
7.0
(b)
Figure 4.22 - HVDC system power reverse. (a) blue Active power MMC-1, green Active power
MMC-2. (b) blue DC voltage MMC-1, green DC voltage MMC-2.
48
49
d+
dM+
dM
dM0
d0
dN0
dN
To extract the sequence components, the three phase voltages are first converted to stationary
reference frame using Clarks transformation (Appendix). To obtain the direct and in-quadrature
components, the Decoupled double synchronous reference frame (DDSRF) and the second-order
generalized integrator (SOGI) methods were investigated. The SOGI-QSG was selected due to its improved
filtering characteristics. The sequence components are calculated as follows, [69]
d
= , , d
d
1 1
> A =
d
2
= , , ,
d
1 d
(5.1)
d
= , , d
d
1 1
d =
= , , ,
d
1 d
(5.2)
Where = L is a phase lagging operator by 90 used to obtain the quadrature components of the input
signal. From (5.1) and (5.2) can be noticed that the positive component d is in phase with d while d is
lagging d by 90. Also the negative sequence component d is in phase with d while d leads d by 90.
The implementation diagram of the SOGI-QSG, together with the positive and negative sequence
calculation is shown in Figure 5.2.
d
dMN
d
d
n
n
d+
d+
+
d
d
d
d
The transfer functions of the SOGI-QSG for the in-phase and quadrature components are [69],
d
2Dg f
(f)
=
d
f + 2Dg f + g
d
2Dg
(f)
=
d
f + 2Dg f + g
(5.3)
(5.4)
Where D determines the bandwidth, chosen to be 2 to have a damping factor = 12, and the centre
frequency obtained from the PLL is g g = 2o 50 rad/s.
The same procedure can be applied to obtain the sequence components of the current. In Figure 5.3 are
presented the positive and negative components of the grid voltage. At t=0.5s the voltage in phase A is
reduce to 20%, as shown in (a). It can be observed that, when the unbalance is imposed, the negative
component of the voltage appears and the amplitude of the positive component is reduced. Although the
calculated settling time of SOGI-QSG is 20.7 ms, because of the PLL dynamics, the time obtained in
simulations is approximately 40 ms.
50
400
kV
0
-400
(a)
400
kV
0
-400
(b)
400
kV
0
-400
(c)
200
kV
100
0
-100
-200
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(d)
Figure 5.3 Extraction of sequence components from unbalanced three phase voltage using DSOGI-QSG. (a)
The angle of the positive and negative voltage components, and , is obtained by implementing the
PLL structure presented in Section 4.2.1 for the positive and negative voltage components respectively. In
Figure 5.4 are shown and when the unbalance is introduced at t=0.5s. It can be observed that the
positive angle is not affected by the unbalance, and the calculation of the negative angle settles 0.025ms
after the fault appears. Since the q-axis component is used in the PI controller to get the frequency and the
phase angle; a phase shift of 90 appears between the calculated angle and the grid.
6.0
rad
4.0
2.0
0.0
(a)
0.0
-2.0
-4.0
rad
51
-6.0
x
0.480
0.500
0.520
0.540
0.560
0.580
0.600
0.620
(b)
Figure 5.4 DSOGI-PLL, simulation results. (a) blue Positive sequence grid angle , green Positive
PLL angle . (b) blue Negative sequence grid angle , green Negative PLL angle
The positive and negative sequence PLL is implemented to align the d-axis with the grid voltage, thus the qaxis components under steady state and unbalanced conditions are zero ( - = - = 0). By this means, the
instantaneous power equations used in the subsequent sections can be simplified.
Description
Rated apparent power
Rated cell voltage
DC link voltage
Arm resistance
Arm inductance
Cell capacitance
Number of sub-modules per arm
Abbreviation
Value
850 MVA
25 kV
400 kV
0.1
90 mH (0.15 pu)
570 F (40 kJ/MVA)
32
The Distributed control method together with the CCSC presented in Chapter 4 was implemented as inner
control strategies.
To assess the performance of the MMC based HVDC transmission system under unbalanced conditions; the
system is subjected to a single-line-to-ground (SLG) fault on the bus PCC-1 (Figure 5.5), dropping the
voltage of phase A in the AC system-1 to 20%. The MMC-1 station operates in inverter mode, controlling
DC-link voltage while the MMC-2 station has active power reference set to 800 MW. Reactive power
references for both stations are set to 0 MVAr. The fault is imposed at 0.5s and cleared after 0.2s.
MMC-1
MMC-2
52
kA
6.0
4.0
2.0
0.0
-2.0
-4.0
-6.0
(a)
kA
4.0
3.0
2.0
1.0
0.0
-1.0
(b)
kA
4.0
3.0
2.0
1.0
0.0
-1.0
(c)
Io
4.0
2.0
0.0
kA
53
-2.0
-4.0
20.0
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(d)
Figure 5.6 Fault propagation, currents in MMC-1. (a)Three phase currents. (b) blue blue - 8 , green - . (d) 5 .
8,
0.850
green -
(c)
In Figure 5.7 the voltages in MMC-1 are shown. In the absence of a grounded midpoint in the DC-link, the
zero sequence current from the grid is transmitted to the DC side, causing undesired oscillations at
fundamental frequency between the DC-link neutral point and ground (a). Moreover, as a result of the
power oscillations caused by the presence of negative sequence components in the grid, double linefrequency ripples with amplitude of 5% appear in the DC-link (b). The voltage profile of one capacitor in
upper and lower arms of phase A, B and C are shown in (c,d,e) respectively.
It can be observed that capacitor voltage waveforms are affected by the corresponding phase currents in
the converter arms. As demonstrated in Section 0; the voltage ripples of the capacitors in the faulted phase
reflect the increment of the phase current, having ripple factor of 28% (c). Phase B and C capacitors have
change in voltage shape due to current distortions.
200
100
kV
0
-100
-200
(a)
840
kV
800
760
(b)
35.0
30.0
kV
25.0
20.0
(c)
35.0
30.0
kV
25.0
20.0
(d)
35.0
30.0
kV
25.0
20.0
20.0
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(e)
Figure 5.7 - Fault propagation, voltages in MMC-1. (a) Neutral point voltage. (b) DC-link voltage. (c,d,e)
Phase A, B, C capacitor voltages respectively.
The presence of negative and zero-sequence current components affects the performance of the whole
HVDC transmission system. The currents of MMC-2 are shown in Figure 5.8, where a strong presence of the
zero-sequence in the converter output currents can be noticed. This could force the overcurrent
protections to trip as well as have negative impact on the AC System-2. Furthermore, the double linefrequency ripples present in the DC-link voltage are reflected in the three-phase difference currents of the
MMC-2 (b), stressing the converter components.
54
4.0
2.0
kA
0.0
-2.0
-4.0
IdiffA2
IdiffB2
(a)
IdiffC2
0.20
0.00
-0.20
-0.40
-0.60
-0.80
-1.00
20.0
kA
55
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(b)
Figure 5.8 - Fault propagation, currents in MMC-2. (a)Three phase AC currents. (b) Difference currents
In transformerless MMC-HVDC transmission systems, without the proper control, faults are propagated
through DC-link and reflected in the healthy converter station affecting the AC grid. Moreover, the fault
currents presented could damage the converter components or require unrealistic over ratings. The next
section focuses on reducing the fault reflection in the HVDC system and mitigating the impact of the
unbalance in the faulty side.
+ :LF
=
+
=
+ :LF
.=
+ :LF
.
+ :LF
+
f
f=
:LF
=
MN
:LF
.=
:LF
MN
-g
:LF
.
:LF
:LF
0 =0
f0
:LF
To overcome the problems mentioned previously, a current control method that involves three sequence
components is implemented. In Figure 5.9 is presented current control structure. The measured AC
currents are transformed into positive and negative sequence components in dq reference frame, using the
appropriate transformation angles. The control of the zero-sequence current, due to its sinusoidal
behaviour, is implemented in stationary reference frame with a PR controller, as proposed in [71].
5.4.1
In order to reduce the fault reflection in the other side of the HVDC system, one of the control targets for
the current control is the elimination of the zero sequence current, i.e.
In Figure 5.10 the simulation waveforms of the MMC-1 system with only positive and zero sequence control
are shown. The control is active during both balanced and unbalanced conditions. The presence of the
negative sequence currents is reflected in the distorted AC currents (a). According to the control objective,
the zero sequence current is effectively suppressed (d); as a result, the oscillations between the DC-link
neutral point and ground are eliminated.
Iabc
Ia_rms
Ib_rms
-2
Ic_rms
Iabc2
kA
6.0
4.0
2.0
0.0
-2.0
-4.0
-6.0
(a)
kA
4.0
3.0
2.0
1.0
0.0
-1.0
(b)
kA
4.0
3.0
2.0
1.0
0.0
-1.0
(c)
Io
4.0
2.0
kA
0.0
-2.0
20.0
-4.0
s
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(d)
Figure 5.10 Zero sequence control, results in MMC-1. (a)Three phase currents. (b) blue (c) blue - 8 , green - . (d) 5 .
0.850
8,
green -
In Figure 5.11 can be observed that the DC-link voltage ripples of 2% are still present (a). These ripples are
product of the power oscillations caused by the negative sequence components in the AC system-1.
56
The elimination of zero sequence current implies lower converter phase currents, therefore the amplitude
of the capacitor ripple is reduced. However, in the faulted phase voltage ripple it is still 22% (b). Phase B
and C capacitors have change in voltage shape due to unbalanced currents in the grid, having voltage ripple
factor of 14%.
840
kV
800
760
(a)
35.0
30.0
kV
25.0
20.0
(b)
35.0
30.0
kV
25.0
20.0
(c)
35.0
30.0
kV
25.0
20.0
20.0
s
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(d)
Figure 5.11 Zero sequence control, results in MMC-1. (a) DC-link voltage (b,c,d) Phase A, B, C
capacitor voltages respectively.
Having eliminated the zero sequence current, the fault is not reflected in the output currents of the MMC-2
as shown in Figure 5.12-a. Though, the ripples in the DC-link are transmitted to the three phases of MMC-2
causing double line-frequency zero sequence ripples in the difference currents (b).
4.0
2.0
kA
0.0
-2.0
-4.0
(a)
-0.20
kA
57
-0.40
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(b)
Figure 5.12 - Zero sequence control, currents in MMC-2. (a) AC currents. (b) Difference currents
5.4.2
Aiming to suppress the ripples in the DC-link, the references for the negative sequence currents are set to 0
(
= 0.
In Figure 5.13 the simulation waveforms of the MMC-1 system with negative sequence control are shown.
According to the control objective, the negative sequence current (a) is effectively suppressed during the
unbalanced conditions; subsequently the ripples in the DC-link voltage are reduced to 1% (b). The
amplitude of the capacitor voltage ripples comparing to the previous case is decreased, resulting in 10% for
the faulted phase (c) and 8% for B, C converter phases (d,e). Since only positive sequence in the AC currents
is present, the voltage profile of B, C phase capacitors is similar to normal operation conditions.
1.0
kA
0.0
-1.0
Vdc_ref
Vdc
(a)
Vdc2
840
kV
800
760
(b)
35.0
30.0
kV
25.0
20.0
(c)
35.0
30.0
kV
25.0
20.0
(d)
35.0
30.0
kV
25.0
20.0
20.0
s
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(e)
Figure 5.13 Negative sequence control, results in MMC-1. (a) blue - 8 , green voltage. (c,d,e) Phase A, B, C capacitor voltages respectively.
0.850
(b) DC-link
58
The presence of the double line-frequency ripples in the DC-link voltage, similar to the previous case, is
reflected in the difference currents of the MMC-2 station, as shown in Figure 5.13. However, the amplitude
of the ripples is decreased, comparing to Figure 5.12-b.
-0.20
kA
59
-0.40
20.0
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
It can be observed that suppressing the negative sequence current results in the reduction of the DC-link
ripples; yet, the presence of the unbalance is still reflected on the difference currents of the unfaulted
converter station.
The CCSC implemented in Chapter 4, acts on the negative sequence double line-frequency component of
the difference currents, transformed into dq signals. Under asymmetric grid conditions, the converter
difference currents get distorted, appearing a positive sequence component. Having unbalanced difference
currents, the PI controller in the dq frame is not able to track the reference. To overcome this problem,
implementation of the CCSC in reference frame using PR controllers is proposed. The PR controllers are
tuned at double line-frequency, acting on both sequences, therefore ensuring elimination of the desired
components.
Moreover, the presence of double line-frequency zero sequence difference current, observed in the
previous section may be produced by the current control strategy used for unbalanced conditions, as
proven in [26]. Aiming to compensate the unbalances in the grid, the converter generates an unbalanced
emf, which produces a zero sequence power component at the converter arms. This component appears as
a double line-frequency current circulating through the converter arms and transmitted over the DC-link.
To eliminate this current, a supplementary zero-sequence control loop with a PR controller tuned at double
line-frequency was added to the CCSC.
The control diagram of the modified circulating current suppression controller is shown in Figure 5.15.
= FF ,
= FF
= FF ,
= FF ,0
:LF
= FF , =0
<= FF ,MN
:LF
:LF
= FF , =0
:LF
= FF ,0 =0
In Figure 5.16 the performance of standard and modified CCSC is compared. As seen in (a), the control
implemented in dq frame does not ensure effective suppression of the circulating current during
unbalanced conditions. Comparing both results, it can be observed, that the modified control eliminates
the oscillatory terms in the difference currents of unfaulted phases and reduces significantly the ripples in
the faulted phase.
0.60
0.40
kA
0.20
0.00
-0.20
(a)
0.60
0.40
kA
0.20
0.00
-0.20
20.0
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(b)
Figure 5.16 Difference currents in MMC-1 with (a) conventional CCSC. (b) modified CCSC
The effective reduction of the zero-sequence oscillations in the difference currents is reflected in the DClink voltage ripples, as shown in Figure 5.17-a. This, results in the suppression of the ripples in MMC-2
60
difference currents (b), which remain as pure DC values except for the short transients at the fault
occurrence and clearance.
840
kV
800
760
(a)
-0.20
kA
61
-0.40
20.0
s
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
(b)
Figure 5.17 Results on MMC-2 with modified CCSC. (a) DC-link voltage. (b) Difference current
The presented simulation cases have shown the effectiveness of using negative and zero sequence controls
under unbalanced conditions to avoid the fault propagation through the DC-link.
(5.5)
Where j, , j5 are the average powers and j , j- , , - , j5 represent the magnitude of the
oscillating terms. The sub-index 0 is used to denote the active power delivered by the zero sequence. If the
zero sequence current is eliminated by appropriate control, j5 term can be neglected. Transforming the
three-phase currents and voltages into synchronous reference frame and splitting into positive and
negative, the power components can be represented as [70]:
j
W
j
V
V jV
V
U-
W -8
Z
V -8
Y
3V
Y = V -
Y 2 V -
V
Y
V -
X
U -8
-8
-8
-8
-
-8
-8
-
-8
-
Z
- Y
Y W
-8
YV
V
-8 Y V
Y U
-8
Y
- X
8Z
Y
Y
8Y
(5.6)
The given relation of sequence components and instantaneous power can be used for calculating the
current references depending on the grid conditions and application requirements. It can be noticed in
(5.6) that the converter has four degrees of freedom to control the AC currents. Selecting the average
active and reactive powers to follow the reference,
j = j = -
(5.7)
two degrees of freedom are left, which means that two control targets can be selected. Three simulation
cases are presented with the following objectives:
5.6.1
The elimination of the negative sequence currents results in balanced converter AC currents under
unbalanced grid voltages. This feature is important for Grid Connected Inverters due to the restrictions
provided by Grid Operators regarding the injection of unbalanced currents into the grid. The control
objectives are defined as:
8
=0
=0
(5.8)
As a result, having - = - = 0 ensured by positive and negative sequence PLL, the references for the
positive sequence current are:
8
j-8
-8
(5.9)
The response of MMC-1 is shown in Figure 5.18. During the fault, the positive component of the d-axis
current 8 is increased by 38% (c), in order to maintain the average output power when phase A voltage
drops by 80%. The negative sequence components are effectively controlled to zero, resulting in balanced
AC converter currents with the increased amplitude. The active and reactive power experience double linefrequency oscillations of 300 MW/MVAr, with the average value maintained at the reference level (e).
62
Cc
Bb
Aa
Vabc_t
Vabc_g
Vabc_t2
400
kV
0
-400
Iabc
(a)s
Ib_rm
Ia_rms
Ic_rms
Iabc2
4.0
2.0
kA
0.0
-2.0
-4.0
(b)
3.0
2.0
kA
1.0
0.0
-1.0
(c)
1.0
kA
0.0
-1.0
(d)
1.5k
1.0k
0.5k
0.0
-0.5k
MW/MVAr
63
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(e)
Figure 5.18 Case 1, results in MMC-1. (a) Grid three phase voltage. (b) Three phase AC currents. (c) blue
8 , green - . (d) blue - 8 , green - . (e) blue Active power, green Reactive power.
In Figure 5.19 are shown the inner parameters of the MMC-1 together with the DC-link voltage. Due to
change in the active power distribution, the difference currents of phase B and C have an increment of 38%
(a). The reduction in the A-phase difference current corresponds to the 74% decrease in active power
supplied by the phase.
The voltage ripples in the sub-module capacitors are affected by the unbalanced conditions. In (b), the
voltage profile of one capacitor in upper and lower arm of phase A is shown. During the fault, the capacitor
voltage ripple contains mainly fundamental component. This is the reflection of the almost pure AC arm
current, since loading of the phase is drastically reduced. The ripple amplitude is increased from 6% to 12%,
having transient peak of 22%. The ripple factor in the capacitor voltages of phase B and C (c,d) is raised up
to 10% due to the increased loading of the mentioned phases. The phase angle of the ripple components in
all the phases remains unchanged due to balanced currents provided by the converter. The DC-link voltage
replicates transients in the sub-modules, with 2.5% peak increase in magnitude. After 0.05s all system
parameters are reaching steady values, showing stable operation under unbalanced grid conditions. When
the fault is cleared, pre-fault characteristics of the system are reached within 3 fundamental cycles.
IdiffA
IdiffB
IdiffC
Idiff_sum
kA
0.80
0.60
0.40
0.20
0.00
-0.20
(a)
30.0
kV
25.0
20.0
Vc_u
Vc_l
Vc_uB
(b)
Vc_lB
Vc_uC
Vc_lC
30.0
kV
25.0
20.0
(c)
30.0
kV
25.0
20.0
Vdc_ref
Vdc
(d)
Vdc2
kV
840
820
800
780
760
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(e)
Figure 5.19 Case 1, results in MMC-1. (a) Difference currents. (b,c,d) Phase A, B, C capacitor voltages
respectively. (e) DC-link voltage
5.6.2
In a HVDC system rated in hundreds of MW, the active power oscillations between the converter station
and the grid might affect the AC system. This could trip the grid relay protections, worsen the fault
conditions. To compensate for the active power oscillations, the control objectives are defined as:
j =0
j- = 0
Using (5.6), the positive and negative sequence current references are obtained:
(5.10)
64
-8 j-
(
' -8
2 -8
=
(
(
3' -8 + ' -8
2 -8
j=
(
(
3' -8 ' -8
2 -8 =
(
(
3' -8
+ ' -8
3'
-8 (
(5.11)
In Figure 5.20 are shown currents and powers of MMC-1. During the fault, the positive component of the
current 8 is increased by 56%. The faulted phase has an increment in the amplitude of 212% and the
unfaulted phases of 38%. The objective of the control strategy is met, having suppressed oscillations in the
active power (d). However, the presence of the 8 increases reactive power oscillations to 650 MVAr.
Iabc
Ia_rms
Ib_rms
Ic_rms
Iabc2
4.0
2.0
kA
0.0
-2.0
-4.0
(a)
3.0
2.0
kA
1.0
0.0
-1.0
(b)
1.0
kA
0.0
-1.0
(c)
1.5k
1.0k
0.5k
0.0
-0.5k
MW/MVAr
65
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
(d)
Figure 5.20 - Case 2, results in MMC-1. (a) Three phase AC currents. (b) blue - 8 , green green - . (d) blue Active power, green Reactive power
0.800
(c) blue -
8,
In Figure 5.21 are presented the inner parameters of the MMC-1 and the DC-link voltage. It can be
observed in (a) a 30% increase in the difference currents of the two unfaulted phases, in order to
compensate the 60% reduction in the active power delivered by phase A. The sub-module capacitors have
ripples of 20% in phase A and 10% in B and C respectively. The phase angle of the voltage ripples in the
phase-legs is changed according to the AC currents, which is slightly reflected in the DC-link voltage ripple.
When the fault is imposed, the DC-link voltage presents 2.5% ripples, reaching steady-state after 0.05s
under unbalanced grid conditions. The pre-fault characteristics of the system are reached within 3
fundamental cycles after the unbalance is cleared.
kA
0.80
0.60
0.40
0.20
0.00
-0.20
(a)
30.0
kV
25.0
20.0
(b)
30.0
kV
25.0
20.0
(c)
30.0
kV
25.0
20.0
(d)
kV
840
820
800
780
760
s
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(e)
Figure 5.21 - Case 2, results in MMC-1. (a) Difference currents. (b,c,d) Phase A, B, C capacitor voltages
respectively. (e) DC-link voltage
5.6.3
Reactive power oscillations between the converter station and the grid have negative impact on the AC
voltage profile, which could worsen the fault conditions. To compensate for the mentioned oscillations, the
control objectives are defined as:
=0
- = 0
(5.12)
66
-8 j-
(
' -8
2 -8
=
(
(
3' -8 ' -8
2 -8
j=
(
(
3' -8 ' -8
2 -8 =
(
(
3' -8
' -8
3'
-8 (
(5.13)
During the fault, the positive component of the current 8 is increased by 23%, as shown in Figure 5.22.
Grid currents are unbalanced, having 40% decrease in the magnitude of the faulted phase and 50% increase
in the unfaulted phases respectively (a). The objective of the control strategy is fulfilled having suppressed
oscillations in the reactive power after 0.05s. The presence of 8 increases active power oscillations which
are 500 MW (66% higher than in Case-1). The average value of power is following the reference.
4.0
2.0
kA
0.0
-2.0
-4.0
(d)
3.0
2.0
kA
1.0
0.0
-1.0
(e)
1.0
kA
0.0
-1.0
(f)
1.5k
1.0k
0.5k
0.0
-0.5k
MW/MVAr
67
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
(g)
Figure 5.22 - Case 3, results in MMC-1. (a) Three phase AC currents. (b) blue - 8 , green green - . (d) blue Active power, green Reactive power
0.800
(c) blue -
8,
In Figure 5.23 are presented the inner dynamics of the MMC-1 and the DC-side. Since the power delivered
by phase A is reduced by 75% during the fault, the difference currents of the two unfaulted phases is
increased by 40% to avoid a drop in the average active power produced (a). It can be noticed that, phases B
and C share active power loading unequally due to unbalance in the AC currents. In phases B and C the 12%
raise in the ripple amplitude corresponds to the increase in the loading (c,d). The phase angle of the
capacitor voltage ripples is changed according to the present unbalanced currents. The DC-link voltage has
transient peaks of 4.5%. After 0.08s converter reaches stable operation with unbalanced grid conditions.
The pre-fault characteristics of the system are reached within 4 fundamental cycles after the fault is
cleared.
kA
0.80
0.60
0.40
0.20
0.00
-0.20
(a)
30.0
kV
25.0
20.0
(b)
30.0
kV
25.0
20.0
(c)
30.0
kV
25.0
20.0
(d)
kV
840
820
800
780
760
0.400
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
(e)
Figure 5.23 - Case 3, results in MMC-1. (a) Difference currents. (b,c,d) Phase A, B, C capacitor voltages
respectively. (e) DC-link voltage
68
69
Control objective
8
=0
j = j- = 0
= - = 0
Amplitude of AC
currents (%)
A
Amplitude of
Power
oscillations
(MW/MVAr)
+38
Amplitude of
capacitor
voltage ripples
(%)
DC-link
voltage
peak (%)
Comments
on AC
currents
300
300
12
10
10
2.5
balanced
+212
+38
+38
650
20
10
10
2.5
unbalanced
-40
+50
+50
500
12
12
4.5
unbalanced
Analyzing the results it can be noticed that elimination of the negative sequence currents results in lower
phase currents. Between the two power compensation methods, the elimination of reactive power
oscillations implies less stresses to the converter in terms of currents and sub-module voltage ripples.
However, it injects unbalanced currents to the grid, which might have limitations defined by the Grid
Operator. To have flexible control over the power oscillations, a modification to the calculation of injected
current sequence components is proposed in [69].
During simulations it was observed that the ripples in the DC-link, when the fault occurs and when is
cleared, were reflected in the difference current of MMC-2. The highest peak current ripple was of 6.3%
and appeared in Case 3. As analyzed in Chapter 3, the ripples in the difference current are not reflected in
the output currents; therefore AC-system 2 is not affected.
All the control strategies ensured stable operation of the converter under unbalanced conditions. The inner
controls of the MMC maintained leg voltages at the reference and balancing effectively the cell voltages
during transients (detailed plots in Appendix). The pre-fault parameters of the system were reached within
3-4 fundamental cycles after the fault was cleared, depending on the control strategy. This complies with
the Grid Code requirements presented in Chapter 2.
]"
(5.14)
Knowing the permissible overcurrent, it is possible to obtain the maximum AC currents allowed. However,
the currents in the converter arms do not represent the AC currents directly. As it was presented in Chapter
3, the converter arm currents are composed of the DC, fundamental AC and double line-frequency
circulating components. Assuming effective suppression of the circulating current, the arm currents can be
expressed as shown:
!
C
(5.15)
(5.16)
Where the DC component represents the active power transfer and the AC component is half of the grid
current, due to its equal share between the arms. Neglecting the losses in the converter, the DC component
of the arm current in a three-phase converter under balanced conditions can be calculated as:
8
j-
(5.17)
Knowing the desired maximum active power and the allowable currents in semiconductors (
limit for the AC grid current amplitude can be obtained as:
h
= 2
,C#
8 ,
), the
(5.18)
Then, having the limits obtained, the grid currents can be effectively controlled by estimating the
instantaneous active or reactive power references from (5.14), depending on the grid voltage:
j-
5.7.1
,C#
,C#
j-
(5.19)
(5.20)
For the converter with rated apparent power of 850 MVA connected to the 400 kV grid, the peak rated
current is:
h
]"
2 850
3 400
= 1.74 D
(5.21)
,"
h
800
1.74
=
+
= 1.2 D
2
3 800
2
(5.22)
70
,C#
= 2h
8 ,
= 2 1.25 1.2
800
= 2 D
2 800
(5.23)
In order to highlight the results obtained with the current limitation strategy the unbalance in Phase A is
imposed at t=0.5s, dropping the voltage to zero. At t=0.6s the current limitation is activated. Two cases are
simulated; in the first case the active power reference is changed to comply with the current limitations. In
the second case the active power is reduced to 0 while full reactive power injection is activated to provide
grid voltage support. The objective of the current-control strategy is the elimination of negative sequence
currents (Section 5.4.2).
5.7.2.1
In Figure 5.24 is presented the response of the MMC-1 when the grid unbalance is imposed. The amplitude
of AC currents is increased by 50% with peak value of 2.8 kA, exceeding the calculated limit of 2 kA.
400
kV
0
-400
Iabc
Ia_rms
Ib_rms
(a)
-2
Ic_rms
Iabc2
2.0
kA
0.0
-2.0
(b)
2.0
kA
0.0
-2.0
(c)
1.0k
MW/MVAr
71
0.5k
0.0
18.0
0.450
0.500
0.550
0.600
0.650
0.700
(d)
Figure 5.24 Current limitation, results in MMC-1 with maximum active power injection. (a) Grid
voltages (b) Three phase AC currents. (c) Upper arm currents in the three phases. (d) blue Active
power, green Reactive power
When the current limitation is activated, a new power reference is calculated from (5.19) in order to
produce maximum active power without exceeding the current limits (d). The new power reference is 630
MW, which represents the reduction in 21% of active power delivered to the grid. The upper arm currents
of the three-phases are shown in (c). After the limitation is activated, the arm currents are effectively kept
under de limit of 1.5 kA. It can be noticed that the arm current of phase A (blue) is purely sinusoidal due to
no active power transfer while B and C contain a DC component responsible of the active power transfer.
5.7.2.2
It can be observed in Figure 5.25 that during the fault, the amplitude of AC currents increases to 2.8 kA
while the arm currents rise to 2 kA, exceeding the limits by 0.5 kA. At time t=0.6s the current limitation is
activated and the active power is ramped to zero. The new AC current limit is calculated using (5.23) as 3 kA
(b). The reactive power reference is calculated according to (5.20) having a value of 1030 MVAr (d).
Cc
Bb
Aa
Vabc_t
Vabc_g
Vabc_t2
400
kV
0
-400
(a)
4.0
2.0
kA
0.0
-2.0
-4.0
(b)
2.0
kA
0.0
-2.0
(c)
Q2
P2
Qo
Pdc
MW/MVAr
1.0k
0.5k
0.0
(d)
350
kV
300
250
200
18.0
0.450
0.500
0.550
0.600
0.650
0.700
(e)
Figure 5.25 - Current limitation, results in MMC-1 with maximum reactive power injection. (a) Grid
voltages (b) Three phase AC currents. (c) Upper arm currents in the three phases. (d) blue Active
power, green Reactive power. (e) Positive sequence component of the Grid voltage
72
73
The arm currents fall into the imposed limits within 3 fundamental cycles, because of the ramped change of
power references. However, after stabilisation, the limits for the arm and AC currents are not exceeded.
With the injection of reactive power, the grid voltage is raised by 8% (e), thus the converter provides grid
voltage support with maximum allowed reactive current injection.
Conclusions
6 Conclusions
The simulation model of the transformerless MMC-HVDC transmission system implemented in
PSCAD/EMTDC was presented. For the inner controls of the converter, a combination of Distributed control
and Direct suppression of circulating current was selected. The HVDC system controls were based on fast
inner current loop, implemented in dq reference frame, and outer power and voltage loops. Grid
synchronisation was achieved by PLL. The step response of the current and DC voltage control was assessed
independently in the simulation. The performance of the MMC-HVDC transmission system was verified for
ramp changes in active and reactive power as well as full active power reverse. The effectiveness of the
inner converter controls was proven, having balanced capacitor voltages and effectively suppressed
circulating current throughout the simulations.
The operation of the MMC-HVDC system under unbalanced grid conditions was tested. First, the
propagation of unbalanced grid faults through the DC-link was investigated. It was concluded, that zero
sequence currents flow through the DC-link, causing voltage oscillations between the neutral point and
ground. These oscillations are directly transmitted to the three phase output currents of the healthy AC
system, which is not acceptable by the Grid Codes. For that reason, an additional loop with zero sequence
current control was implemented. It was observed that by suppressing the zero sequence current, it was
possible to eliminate the propagation of the fault to the healthy AC system. This enhances the feature of
HVDC transmission to act as isolation between the two AC grids. Moreover, the voltage oscillations in the
DC-link neutral point were removed and the capacitor voltage ripples were reduced, minimizing the
stresses on the faulted converter station. However, the presence of the unbalance was still reflected as a
ripple in the DC-link voltage, disturbing the inner currents of the unfaulted converter station. To solve this
problem a negative sequence current control was implemented. It was observed that, although the
negative sequence currents got eliminated, the ripples in the DC-link were still present. A modification to
the inner current control of the MMC was proposed to eliminate the zero sequence component of the
difference current and, as a consequence, the voltage ripples in the DC-link. The disturbance to the
unfaulted converter was effectively reduced to short transients in the arm currents.
Three strategies to calculate the current references under unbalanced grid conditions were studied. The
possibility of compensating the power oscillations by means of negative sequence current injection was
verified in simulations. The effect of different strategies on the inner parameters of the converter and the
DC-link was analyzed. It was concluded that by maintaining balanced AC currents, less stresses on the
converter are applied due to smaller voltage ripples in the capacitors. The inner controls of the converter
presented good performance under unbalanced conditions.
Although the mentioned control strategies presented a satisfactory performance, the current in the
converter during the fault increased to unacceptable values that could trip an overcurrent protection or
damage the converter components. For this reason, a strategy to limit the converter currents was
implemented. The relation between the grid and converter arm currents was derived, making current
limitation approach straightforward. Two cases with current limitation were simulated. First, the active
power of the converter was reduced based on the grid voltage and the converter current limits. This
approach showed a fast response without significant transients. In the second case, the converter active
74
75
power was reduced to zero with full reactive power injection. The injection of reactive power boosted the
grid voltage, thus supporting the grid during the fault. The effectiveness of the converter current limitation
was verified, demonstrating the fault ride-through capability of the converter. Both studied cases required
partial or full reduction in active power transmission through the DC-link. This might have a negative
impact, since the other converter station will be force to reduce or stop power transmission.
The results obtained in the simulations demonstrate the great potential that MMC has for HVDC
applications, having the possibility of being implemented in transformerless connection with a good fault
ride-through capabilities.
Analysis of the converter performance under different types of balanced and unbalanced faults;
Evaluation of delays on the controllability of the converter;
Implementation of different control architectures using advanced control techniques (model
predictive control, deadbeat control, multivariable control);
Investigation of MMC-HVDC performance in meshed DC transmission grids;
Implementation of the control strategies in small-scale converter prototype.
References
References
[1]
[2]
[3]
[4]
[5]
[6]
Blaabjerg, F.; Chen, Z.; Teodorescu, R.; Iov, F., "Power Electronics in Wind Turbine Systems," Power Electronics
and Motion Control Conference, 2006. IPEMC 2006. CES/IEEE 5th International , vol.1, no., pp.1,11, 14-16 Aug.
2006
[7]
O. Heyman, L. Weimers, M. Bohl, HVDC - A key solution in future transmission systems, WEC Montreal 2010
[8]
[9]
Morton, A.B.; Cowdroy, S.; Hill, J. R A; Halliday, M.; Nicholson, G.D., "AC or DC? economics of grid connection
design for offshore wind farms," AC and DC Power Transmission, 2006. ACDC 2006. The 8th IEE International
Conference on , vol., no., pp.236,240, 28-31 March 2006
[10] http://en.wikipedia.org/wiki/List_of_HVDC_projects
[11] R. Marquardt, A. Lesnicar, J. Hildinger, Modulares Stromrichterkonzept fr Netzkupplungsanwendungen bei
hohen Spannungen, ETG-Fachtagung, Bad Nauheim, VDE 2002
[12] Franquelo, L.G.; Rodriguez, J.; Leon, J.I.; Kouro, S.; Portillo, R.; Prats, M.A.M., "The age of multilevel converters
arrives," Industrial Electronics Magazine, IEEE , vol.2, no.2, pp.28,39, June 2008
[13] Siemens, HVDC plus Basics and Principle of Operation, online:
http://www.energy.siemens.com/hq/pool/hq/powertransmission/HVDC/HVDC_Plus_Basic%2520and%2520Princ
ipals.pdf
[14] N. Mahimkar, G. Persson, C. Westerlind, HVDC Technology for Large Scale Offshore Wind Connections,
Smartelec, Vadodara 2013
[15] Alstom, HVDC-VSC: transmission technology of the future, online:
http://www.alstom.com/Global/Grid/Resources/Documents/Smart%20Grid/Think-Grid-08-%20EN.pdf
[16] Ahmed, Noman; Norrga, Staffan; Nee, H-P; Haider, A.; Van Hertem, D.; Lidong Zhang; Harnefors, Lennart, "HVDC
SuperGrids with modular multilevel converters The power transmission backbone of the future,"Systems,
Signals and Devices (SSD), 2012 9th International Multi-Conference on , vol., no., pp.1,7, 20-23 March 2012
[17] C. Oates, Future Topologies for VSC HVDC, Vestas symposium, Aalborg, 2012
[18] Minyuan Guan; Zheng Xu; Huijie Li, "Analysis of DC voltage ripples in modular multilevel converters," Power
System Technology (POWERCON), 2010 International Conference on , vol., no., pp.1,6, 24-28 Oct. 2010
[19] Qingrui Tu; Zheng Xu; Hongyang Huang; Jing Zhang, "Parameter design principle of the arm inductor in modular
multilevel converter based HVDC," Power System Technology (POWERCON), 2010 International Conference on ,
vol., no., pp.1,6, 24-28 Oct. 2010
[20] Hagiwara, M.; Akagi, H., "Control and Experiment of Pulsewidth-Modulated Modular Multilevel Converters,"
Power Electronics, IEEE Transactions on , vol.24, no.7, pp.1737,1746, July 2009
[21] Hagiwara, M.; Maeda, R.; Akagi, H., "Control and Analysis of the Modular Multilevel Cascade Converter Based on
Double-Star Chopper-Cells (MMCC-DSCC)," Power Electronics, IEEE Transactions on , vol.26, no.6, pp.1649,1658,
June 2011
[22] Antonopoulos, A.; Angquist, Lennart; Nee, H-P, "On dynamics and voltage control of the Modular Multilevel
Converter," Power Electronics and Applications, 2009. EPE '09. 13th European Conference on , vol., no., pp.1,10,
8-10 Sept. 2009
76
77
[23] Angquist, Lennart; Antonopoulos, A.; Siemaszko, D.; Ilves, K.; Vasiladiotis, M.; Nee, H-P, "Inner control of Modular
Multilevel Converters - An approach using open-loop estimation of stored energy," Power Electronics Conference
(IPEC), 2010 International , vol., no., pp.1579,1585, 21-24 June 2010
[24] Qingrui Tu; Zheng Xu; Lie Xu, "Reduced Switching-Frequency Modulation and Circulating Current Suppression for
Modular Multilevel Converters," Power Delivery, IEEE Transactions on , vol.26, no.3, pp.2009,2017, July 2011
[25] Bergna, G.; Berne, E.; Egrot, P.; Lefranc, P.; Arzande, A.; Vannier, J.-C.; Molinas, M., "An Energy-Based Controller
for HVDC Modular Multilevel Converter in Decoupled Double Synchronous Reference Frame for Voltage
Oscillation Reduction," Industrial Electronics, IEEE Transactions on , vol.60, no.6, pp.2360,2371, June 2013
[26] Qingrui Tu; Zheng Xu; Yong Chang; Li Guan, "Suppressing DC Voltage Ripples of MMC-HVDC Under Unbalanced
Grid Conditions," Power Delivery, IEEE Transactions on , vol.27, no.3, pp.1332,1338, July 2012
[27] Allebrod, S.; Hamerski, R.; Marquardt, R., "New transformerless, scalable Modular Multilevel Converters for
HVDCtransmission," Power Electronics Specialists Conference, 2008. PESC 2008. IEEE , vol., no., pp.174,179, 15-19
June 2008
[28] Minyuan Guan; Zheng Xu, "Modeling and Control of a Modular Multilevel Converter-Based HVDC System Under
Unbalanced Grid Conditions," Power Electronics, IEEE Transactions on , vol.27, no.12, pp.4858,4867, Dec. 2012
[29] Florin Iov, Anca Daniela Hansen, Mapping of grid faults and grid codes, Ris National Laboratory, DTU, 2007
[30] Jing Huang; Corzine, K., "Extended operation of flying capacitor multilevel inverters," Industry Applications
Conference, 2004. 39th IAS Annual Meeting. Conference Record of the 2004 IEEE , vol.2, no., pp.813,819 vol.2, 3-7
Oct. 2004
[31] Liu, Y.H.; Arrillaga, J.; Watson, N.R., "Cascaded H-Bridge Voltage ReinjectionPart II: Application to HVDC
Transmission," Power Delivery, IEEE Transactions on , vol.23, no.2, pp.1200,1206, April 2008
[32] S. Khomfoi and L. M. Tolbert, Chapter 31. Multilevel Power Converters, Transform, pp. 1-50, 1975.
[33] Amirnaser Yazdani, Reza Iravani, Voltage-Sourced Converters in Power Systems, Wiley-IEEE press, 2010
[34] Marquardt, R., "Modular Multilevel Converter: An universal concept for HVDC-Networks and extended DC-Busapplications," Power Electronics Conference (IPEC), 2010 International , vol., no.,pp.502,507, 21-24 June 2010
[35] Solas, E.; Abad, G.; Barrena, J.A.; Aurtenetxea, S.; Carcar, A.; Zajac, L., "Modular Multilevel Converter With
Different
Submodule
ConceptsPart
II:
Experimental
Validation
and
Comparison
for
HVDC Applications,"Industrial Electronics, IEEE Transactions on , vol.60, no.10, pp.4536,4545, Oct. 2013
[36] Akagi, H., "Classification, Terminology, and Application of the Modular Multilevel Cascade Converter
MMCC)," Power Electronics, IEEE Transactions on , vol.26, no.11, pp.3119,3130, Nov. 2011
[37] Lesnicar, A.; Marquardt, R., "An innovative modular multilevel converter topology suitable for a wide power
range," Power Tech Conference Proceedings, 2003 IEEE Bologna , vol.3, no., pp.6 pp. Vol.3,, 23-26 June 2003
[38] Holtz, J., "Pulsewidth modulation-a survey," Industrial Electronics, IEEE Transactions on , vol.39, no.5, pp.410,420,
Oct 1992
[39] Bakhshai, A.R.; Saligheh Rad, H.R.; Joos, G., "Space vector modulation based on classification method in threephase multi-level voltage source inverters," Industry Applications Conference, 2001. Thirty-Sixth IAS Annual
Meeting. Conference Record of the 2001 IEEE , vol.1, no., pp.597,602 vol.1, Sept. 30 2001-Oct. 4 2001
[40] Weidong Jiang; Qunjing Wang; Quan Chen; Cungang Hu, "A simplified SVPWM control algorithm for multilevel
VSI based two-level SVPWM control algorithm," Electrical Machines and Systems, 2007. ICEMS. International
Conference on , vol., no., pp.43,49, 8-11 Oct. 2007
[41] Dorin O. Neacsu, Space vector modulation an introduction, 2001
[42] M. Vasiladiotis, Analysis, Implementation and Experimental Evaluation of Control Systems for Modular
Multilevel Converter. KTH Electrical Engineering, 2009
[43] Rongfeng Yang; Gaoling Wang; Yong Yu; Dianguo Xu; Zhuang Xu, "A fast generalized multilevel voltage SVPWM
algorithm based on voltage decomposition," Industrial Electronics and Applications (ICIEA), 2010 the 5th
IEEE Conference on , vol., no., pp.1220,1224, 15-17 June 2010
[44] Konstantinou, G.S.; Agelidis, V.G., "Performance evaluation of half-bridge cascaded multilevel converters
operated with multicarrier sinusoidal PWM techniques," Industrial Electronics and Applications, 2009. ICIEA
2009. 4th IEEE Conference on , vol., no., pp.3399,3404, 25-27 May 2009
References
[45] Rohner, S.; Bernet, S.; Hiller, M.; Sommer R., "Modulation, Losses, and Semiconductor Requirements of Modular
Multilevel Converters," Industrial Electronics, IEEE Transactions on , vol.57, no.8, pp.2633,2642, Aug. 2010
[46] Zixin Li; Ping Wang; Haibin Zhu; Zunfang Chu; Yaohua Li, "An Improved Pulse Width Modulation Method for
Chopper-Cell-Based Modular Multilevel Converters," Power Electronics, IEEE Transactions on , vol.27, no.8,
pp.3472,3481, Aug. 2012
[47] Vijay K. Sood, HVDC and FACTS Controllers, Boston: Kluwer Academic Publishers, 2004
[48] Bahrman, M.P., "HVDC transmission overview," Transmission and Distribution Conference and Exposition, 2008.
T&D. IEEE/PES , vol., no., pp.1,7, 21-24 April 2008
[49] Flourentzou, N.; Agelidis, V.G.; Demetriades, G.D., "VSC-Based HVDC Power Transmission Systems: An
Overview," Power Electronics, IEEE Transactions on , vol.24, no.3, pp.592,602, March 2009
[50] S. Chaudhary, Control and Protection of Wind Power Plants with VSC-HVDC Connection, PhD thesis, Aalborg
University, Institute of Energy Technology, Denmark, 2011.
[51] ABB, Its time to connect - Technical description of HVDC Light technology.
[52] Andersen, B.R.; Xu, L.; Horton, P. J.; Cartwright, P., "Topologies for VSC transmission," Power Engineering
Journal , vol.16, no.3, pp.142,150, June 2002.
[53] The Grid Code, Issue 5 Rev.1, 2012, National Grid Electricity Transmission.
[54] EWEA, European offshore statistics 2012, online:
http://www.ewea.org/fileadmin/files/library/publications/statistics/European_offshore_statistics_2012.pdf
[55] Peralta, J.; Saad, H.; Dennetiere, S.; Mahseredjian, J.; Nguefeu, S., "Detailed and Averaged Models for a 401-Level
MMCHVDC System," Power Delivery, IEEE Transactions on, vol.27, no.3, pp.1501,1508, July 2012
[56] ABB, DolWin2 Project, online:
http://www05.abb.com/global/scot/scot245.nsf/veritydisplay/391e0bfbf3e6b11f48257a390046bf5e/$file/dolwi
n2%20flyer%201jnb002900d0006%20rev-.pdf
[57] Bjorn Jacobson, et al., VSC-HVDC Transmission with Cascaded Two-Level Converters, CIGRE 2010
[58] Qiang Song; Wenhua Liu; Xiaoqian Li; Hong Rao; Shukai Xu; Licheng Li, "A Steady-State Analysis Method for a
Modular Multilevel Converter," Power Electronics, IEEE Transactions on , vol.28, no.8, pp.3702,3713, Aug. 2013
[59] Harnefors, L.; Antonopoulos, A.; Norrga, S.; Angquist, L.; Nee, H.-P., "Dynamic Analysis of Modular Multilevel
Converters," Industrial Electronics, IEEE Transactions on , vol.60, no.7, pp.2526,2537, July 2013
[60] Harnefors, L.; Norrga, S.; Antonopoulos, A.; Nee, H.-P., "Dynamic modelling of modular multilevel
converters," Power Electronics and Applications (EPE 2011), Proceedings of the 2011-14th European Conference
on , vol., no., pp.1,10, Aug. 30 2011-Sept. 1 2011
[61] Ilves, K.; Antonopoulos, A.; Norrga, S.; Angquist, L.; Nee, H.-P., "Controlling the ac-side voltage waveform in a
modular multilevel converter with low energy-storage capability," Power Electronics and Applications (EPE 2011),
Proceedings of the 2011-14th European Conference on , vol., no., pp.1,8, Aug. 30 2011-Sept. 1 2011
[62] Xu She; Huang, A.; Xijun Ni; Burgos, R., "AC circulating currents suppression in modular multilevel
converter," IECON 2012 - 38th Annual Conference on IEEE Industrial Electronics Society , vol., no., pp.191,196, 2528 Oct. 2012
[63] A. Timbus, Grid Monitoring and Andvanced Control of Distributed Power Generation Systems. PhD thesis,
Aalborg University, Institute of Energy Technology, Denmark, May 2007.
[64] Chandra Bajracharya, Control of VSC-HVDC for wind power, MSc thesis, Norwegian University of Science and
Technology, Department of Electrical Power Engineering, Norway 2008
[65] Chandra Bajracharya, Marta Molinas, Understanding of tuning techniques of converter controllers for VSCHVDC, Nordic Workshop on Power and Industrial Electronics, June 9-11, 2008
[66] EirGrid plc, HVDC Interconnectors - Proposed Modifications to the Grid Code, 2009
[67] Lopez, O.; Teodorescu, R.; Doval-Gandoy, J., "Multilevel transformerless topologies for single-phase gridconnected converters," IEEE Industrial Electronics, IECON 2006 - 32nd Annual Conference on , vol., no.,
pp.5191,5196, 6-10 Nov. 2006
[68] R. E. Fehr, A Novel Approach for Understanding Symmetrical Components and Sequence Networks of ThreePhase Power Systems, 2006
78
79
[69] R. Teodorescu, M. Liserre, P. Rodriguez, Grid Converters for Photovoltaic and Wind Power Systems, Wiley-IEEE
press, 2011
[70] Gnanarathna, U.N.; Gole, A.M.; Jayasinghe, R.P., "Efficient Modeling of Modular Multilevel HVDC Converters
(MMC) on Electromagnetic Transient Simulation Programs," Power Delivery, IEEE Transactions on , vol.26, no.1,
pp.316,324, Jan. 2011
[71] K. Ma, M. Liserre, and F. Blaabjerg, Power Controllability of Threephase Converter with Unbalanced AC Source.
[72] Saeedifard, M.; Iravani, R., "Dynamic Performance of a Modular Multilevel Back-to-Back HVDC System," Power
Delivery, IEEE Transactions on , vol.25, no.4, pp.2903,2912, Oct. 2010
Appendix
Appendix
Parks Transformation
Three phase variables are transformed into rotating reference frame dq.
d8
d
d5
W N;f
2V
= Vf
3V
V 1
U 2
2o
4o
N;f Z
3
3 Y d
2o
4o
f Y d
3
3 Y d
1
1
Y
X
2
2
N;f
Clarkes Transformation
1
1
W1
Z
2
2Y d
V
2
3
3
= V0
Y d
3V
2
2Y d
1 Y
V1 1
U2 2
2 X
DC cable parameters
DC cable is implemented in frequency-dependent PSCAD model.
80
Case-1. Elimination of negative sequence currents. Sub-module capacitor voltages of MMC-1 (left) and MMC-2 (right). Fault at MMC-1 at t=0.5-0.7s
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(f)
(g)
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
Vc_lB
Vc_uB
(h)
Vc_u
kV
25.0
25.0
20.0
20.0
(j)
18.0
s
0.450
(i)
Vc_l
30.0
30.0
kV
0.500
0.550
0.600
(k)
18.0
0.650
0.700
0.750
0.800
0.850
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
Appendix
Case-2. Elimination of active power oscillations. Sub-module capacitor voltages of MMC-1 (left) and MMC-2 (right). Fault at MMC-1 at t=0.5-0.7s
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(a)
(b)
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(c)
Vc_u
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(e)
18.0
s
0.450
(d)
Vc_l
0.500
0.550
0.600
(f)
18.0
0.650
0.700
0.750
0.800
0.850
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
Case-3. Elimination of reactive power oscillations. Sub-module capacitor voltages of MMC-1 (left) and MMC-2 (right). Fault at MMC-1 at t=0.5-0.7s
Vc_l
Vc_u
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(a)
(b)
30.0
kV
kV
30.0
25.0
25.0
20.0
20.0
(c)
Vc_u
30.0
25.0
25.0
20.0
20.0
(e)
18.0
s
0.450
(d)
Vc_l
30.0
kV
kV
0.500
0.550
0.600
(f)
18.0
0.650
0.700
0.750
0.800
0.850
0.450
0.500
0.550
0.600
0.650
0.700
0.750
0.800
0.850
Appendix