Note: Descriptions are shown in the official language in which they were submitted.
CA 02772713 2012-02-29
=
1
A METHOD AND APPARATUS FOR CALCULATING INSERTION INDECES
FOR A MODULAR MULTILEVEL CONVERTER
TECHNICAL FIELD
The invention relates to the calculation of insertion
indices comprising data representing the portion of
available submodules of a modular multilevel converter
that should be in voltage insert mode.
BACKGROUND
The concept Modular Multilevel Converter (M2C) denotes a
class of Voltage Source Converter (VSC). It has one or
several phase legs connected in parallel between two DC
bars, a positive DC+ and a negative DC-. Each phase leg
consists of two series-connected converter arms. The
connection point between the converter arms constitutes
an AC terminal for the leg.
Each arm consists of a number (N) of submodules. Each
submodule has two terminals. Using these terminals the
submodules in each arm are series-connected so that they
form a string. The end terminals of the string constitute
the connection terminals of the arm. By controlling
individual modules in each arm, a voltage corresponding
to the accumulation of insertion voltages can be provided
on the AC terminal.
Such a converter is known from DE10103031. In this
document, a method to equalize the voltages in the
submodule capacitors within the arm is briefly described.
For each arm, a modulator determines when the number of
inserted submodules shall change. The principle for
equalizing is that, at each instant when a change of the
number of inserted submodules is commanded, a selection
2
mechanism chooses the submodule to be inserted or
bypassed depending on the actual current direction in the
arm (charging or discharging) and the corresponding
available submodules in the arm (bypassed highest
voltage/bypassed lowest voltage/inserted highest
voltage/inserted lowest voltage). Such a selection
mechanism aims to achieve that the DC voltage across the
DC capacitors in the submodules are equal, uc,sm(t).
A problem with the prior art is the presence of a
circulation current going through the legs between the DC
terminals.
SUMMARY
An object of the invention is to reduce circulation
currents going through the legs between the DC terminals.
In accordance with an aspect of the present invention there
Is provided a method for calculating insertion indices for
a phase leg of a DC to AC modular muliilevel converter,
the converter via a control apparatus, said control
apparatus having at least one controller, comprising one
phase leg between upper and lower DC source common bars
for each phase, each phase leg comprising two serially
connected arms, wherein an AC output for each phase leg is
connected between the two serially connected arms, wherein
each arm comprises a number of submodules, wherein each
slibmodule can be in a bypass state or a voltage insert
mode, the insertion index comprising data representing a
portion of available submodules that should be in the
voltage insert mode for a particular arm, the method
comprising the steps of: calculating, in said controller,
a desired arm voltage ucu(t) for an upper arm connected to
the upper DC source common bar and a desired arm voltage
CA 2772713 2017-06-06
1.2ci..(t) for a lower arm connected to the lower DC source
common bar, obtaining, in said controller, values
representing actual total arm voltages ii140 and u;t.,,(0 in
the upper arm and lower arm, respectively, by measuring
the sum of capacitor voltages in each arm using sensors in
the submodules, and calculating, in said controller, the
insertion indices nu(t) and ni,(t) for the upper and lower
arm, respectively, using the respective desired arm
voltage and the respective value representing the total
actual arm voltage, as
ucu(t)u(t)
nu = njt) = _____ =
ucu(t) ucL,M'
wherein the step of calculating desired arm voltages for
a phase leg comprises calculating
upp(t) - up/2 - ev(t) - udffl(t)
for the upper arm, and calculating
UCL(t) =up/2 -r ev ( ) UcLff (t)
for the lower arm, where up(t) represents desired upper
arm voltage where un(t) represents desired lower arm
voltage, up represents a voltage between the upper and
lower DC source common bars, ev(t) represents a reference
inner AC output voltage and udiff(t) represents a control
voltage to control a current passing through the whole
phase leg, and calculating
udiff ) Udiffl (t) Udiff2 (t)
where UcLffl(t) represents a voltage obtained by summing
energy in the arms of the leg and udiff2(t) represents a
CA 2772713 2017-06-06
4
voltage obtained by calculating a difference in energy
between the arms of the leg.
The step of obtaining a value representing actual arm
voltage may comprise calculating
uaiff2 (t ) = Udiff2 Cos (u)at + )
where adif=2 represents an error between total upper arm,
energy and total lower arm energy, col represents the
angular velocity of the network frequency and *
represents the angle given by * = L(R jwiL) where R
IC represents the resistance of the converter arm and L
represents the inductance of the converter arm.
The step of obtaining values representing actual arm
voltages may comprise: calculating 11,3(t), actual voltage
for the upper arm, using Carm, capacitance for the arm,
15 id:_ffo, DC current passing through the two serially
connected arms of the phase leg, WiT(t), desired average
energy in the upper arm, ev, amplitude of reference for
the inner AC output voltage, iv, amplitude of AC output
current, p, a phase difference between iv(t) and ev(t),
20 DC current circulating through the two series-connected
arms, and calculating la,.(t), actual voltage for the lower
arm, using Caim, capacitance for the arm, idiffo, DC
current passing through the two serially connected arms
of the phase leg, W(t), desired average energy in the
25 lower arm, ev, amplitude of reference for inner AC output
voltage, iv, amplitude of AC output current p, a phase
difference between iv(t) and ev(t), ia.iffo, DC current
circulating through the two series-connected arms.
CA 2772713 2017-06-06
5
The step of obtaining a value representing actual arm
voltage may comprise calculating
evivcosT
idiff0= /
Up ¨ 4Ri.icoscp
where cp represents a phase difference between iv(t) and
ev(t), up represents a voltage between the upper and
lower DC source common bars and R represents the
resistance of the converter arm.
The step of obtaining a value representing actual arm
voltage may comprise calculating
4(t) 11/12t)
where W(t) represents instantaneous energy in the upper
arm and is calculated as follows:
Wri=u (t--) = sin co.t. +
ii,
¨ ) __ v sin (o.),t + (p)
_fC sin(26),t + T)
2 2w, 8,
and calculating
z
(t) =12WoL(t)
VCam
where Wõ(t) represents instantaneous energy in the lower
arm and is calculated as follows:
CA 2772713 2017-06-06
ev f0
W,t(t) = +
CLC sin co.t.
'ev i
(-2 ¨ R3-diffc sin (wit + (13) v sin(2wit + T)
2w, 8w1
where col represents the angular velocity of the network
frequency, up represents a voltage between the upper and
lower DC source common bars, R represents the resistance
of the converter arm.
The step of obtaining a value representing actual arm
voltage may comprise measuring voltages of the submodules
of the arm, and summing these measured voltages.
The insertion index may comprise data representing a
direction of the inserted voltage.
In accordance with another aspect of the present
invention there is provided A method for calculating
insertion indices for a phase leg of a DC to AC modular
multilevel converter, via a control apparatus, said
control apparatus having at least one controller, the
converter comprising one phase leg between upper and
lower DC source common bars for each phase, each phase
leg comprising two serially connected arms, wherein an AC
output for each phase leg is connected between the two
serially connected arms, wherein each arm comprises a
number of submodules, wherein each submodule can be in a
bypass state or a voltage insert mode, the insertion
index comprising data representing a portion of available
submodules that should be in the voltage insert mode for
a particular arm, the method comprising the steps of:
calculating, in said controller, a desired arm voltage
ucu(t) for an upper arm connected to the upper DC source
CA 2772713 2017-06-06
6a
common bar and a desired arm voltage ua(t) for a lower
arm connected to the lower DC source common bar,
obtaining, in said controller, values representing actual
total arm voltages //7,1(0 and /4-.,(0 in the upper arm and
lower arm, respectively, and calculating, in said
controller, the insertion indices nu(t) and nL(t) for
the upper and lower arm, respectively, using the
respective desired arm voltage and the respective value
representing the total actual arm voltage, as
uM cucl(t)
nu(t) =-- nr,(t) = ti
ucu(t) uci,(0 r
wherein the step of calculating desired arm voltages for
a phase leg comprises calculating
ucu(t) = up/2 - ev(t) - udiff(t)
for the upper arm, and calculating
ucL(t) = up/2 + ev(t) - udi.f(t)
for the lower arm, where ucu(t) represents desired upper
arm voltage where un(t) represents desired lower arm
voltage, up represents a voltage between the upper and
lower DC source common bars, ev(t) represents a reference
inner AC output voltage and udiff(t) represents a control
voltage to control a current passing through the whole
phase leg, wherein the step of obtaining values
representing actual arm voltages comprises: calculating
I
u(t), actual total voltage for the upper arm, using Calm,
capacitance for the arm, id;ffo, DC current passing through
the two serially connected arms of the phase leg, W(t),
desired average energy in the upper arm, Ov, amplitude of
reference for the inner AC output voltage, iv, amplitude
CA 2772713 2017-06-06
6b
of AC output current, cp, a phase difference between iv(t) and
ev(t), DC current circulating through the two series-connected
arms, and calculating ucz,(t), actual total voltage for the lower
arm, using Cam, capacitance for the arm, idiffOt DC current
passing through the two serially connected arms of the phase
W(t), desired average energy in the lower arm, ev, amplitude
of reference for inner AC output voltage, iv, amplitude of AC
output current o, a phase difference between iv(t) and ev(t),
idiff0, DC current circulating through the two series-connected
arms.
=
In a further aspect of the present invention there is provided
An apparatus for calculating insertion indices for a phase leg
of a DC to AC modular multilevel converter, the converter
comprising one phase Leg between upper and lower DC source common
bars for each phase, each phase leg comprising two serially
connected arms, wherein an AC output for each phase leg is
connected between the two serially connected arms, wherein each
arm comprises a number of submodules, wherein each submodule can
be in a bypass state or a voltage insert mode, the insertion
index comprising data representing a portion of available
submodules that should be in the voltage insert mode for a
particular arm, the apparatus comprises:
a controller arranged to calculate a desired arm voltage ucu(t)
for an upper arm connected to the upper DC source common bar and
a desired arm voltage ucL(t) for a lower arm connected to the
lower DC source common bar, to obtain values representing actual
total arm voltages ucE,u(t) and u1_(t) in the upper arm and lower
arm, respectively, by measuring the sum of capacitor voltages
in each arm using sensors in the submodules, and to
CA 2772713 2017-06-06
. õ
6c
calculate the insertion indices nu(t) and ni,(t) for the
upper and lower arm, respectively, using the respective
desired arm voltage and the respective value representing
the total actual arm voltage, as
2ucc(t) ua(t)
n = u(t) ni,(t)
ticu4) uõ(t)'
wherein the calculating of desired arm voltages for a phase
leg comprises calculating
uou(t) = up/2 - ev(t) Udiff ( t )
for the upper arm, and calculating
uoL(t) = up/2 ev(t) - udiff(t)
for the lower arm, where ucu(t) represents desired upper
arm voltage, ucL(t) represents desired lower arm voltage,
up represents a voltage between the upper and lower DC
source common bars, ev(t) represents a reference inner AC
output voltage and uaiff(t) represents a control voltage to
control a current passing through the whole phase leg, and
calculating
Ud4 ff ( t ) = Udiffi (t) Udiff2
where udif--_ (t) represents a voltage obtained by summing
energy in the arms of the leg and Udiff2(t) represents a
voltage obtained by calculating a difference in energy
between the arms of the leg.
In accordance with another aspect of the Present invention
there is provided An apparatus for calculating insertion
indices for a phase leg of a DC to AC modular multilevel
converter, the converter comprising one phase leg between
upper and lower DC source common bars for each phase, each
CA 2772713 2017-06-06
6d
phase leg comprising two serially connected arms, wherein
an AC output for each phase leg is connected between the
two serially connected arms, wherein each arm comprises a
number of submodules, wherein each submodule can be in a
bypass state or a voltage insert mode, the insertion index
comprising data represen-oing a portion of available
submodules that should be in the voltage insert mode for
a particular arm, the apparatus comprises: a controller
arranged to calculate a desired arm, voltage ucp(t) for an
upper arm connected to the upper DC source common bar and
a desired arm voltage u01(t) for a lower arm connected to
the lower DC source common bar, to obtain values
representing actual total arm voltages u(t) and uL(0 in
the upper arm and lower arm, respectively, and to calculate
insertion indices for the upper and lower arm,
respectively, using the respective desired arm voltage and
the respective value representing the total actual arm
voltage, as
ucu(t)ucL,(t)
nu(t) = ni,(t) = I ticL,( '
=
ucu(t)
0
wherein the calculating of desired arm voltages for a phase
leg comprises calculating
ucu(t) = up/2 - ev(t) - udffl(t)
for the upper arm, and calculating
ucL(t) = up/2 + ev(t) - udiff(t)
for the lower arm, where ucu(t) represents desired upper
arm voltage, ucL(t) represents desired lower arm voltage,
up represents a voltage between the upper and lower DC
source common bars, ev(t) represents a reference inner AC
CA 2772713 2017-06-06
6e
output voltage and Udiff (t ) represents a control voltage to
control a current passing through the whole phase leg,
wherein the step of obtaining values representing actual
arm voltages comprises: calculating u.(t) , actual total
voltage for the upper arm, using Carm, capacitance for the
arm, icLffOr DC current passing through the two serially
connected arms of the phase leg, ICIET(t), desired average
energy in the upper arm, ev, amplitude of reference for
the inner AC output voltage, iv, amplitude of AC output
current, cp, a phase difference between iv ( t ) and ev ( t), DC
current circulating through the two series-connected arms,
and calculating ucI,(t) , actual total voltage for the lower
arm, using Car," capacitance for the arm, 1cliff , DC current
passing- through the two serially connected arms of the
phase leg, Wc.:=.(t), desired average energy in the lower arm,
ev, amplitude of reference for inner AC output voltage, iv,
amplitude of AC output current cp, a phase difference
between iv (t) and ev (t) , id.4'f0, DC current circulating
through the two series-connected arms.
Generally, all terms used in the application are to be
interpreted according to their ordinary meaning in the
technical field, unless explicitly defined otherwise
herein. All references to "a/an/the element, apparatus,
component, means, step, etc. '1 are to be interpreted
CA 2772713 2017-06-06
CA 02772713 2012-02-29
7
openly as referring to at least one instance of the
element, apparatus, component, means, step, etc., unless
explicitly stated otherwise. The steps of any method
disclosed herein do not have to be performed in the exact
order disclosed, unless explicitly stated.
BRIEF DESCRIPTION OF THE DRAWINGS
The invention is now described, by way of example, with
reference to the accompanying drawings, in which:
Fig 1 is a schematic diagram of phase legs and arms,
Fig 2 is a schematic circuit model of a phase leg of Fig
1,
Fig 3 is a Nichols plot for an open loop transfer
function,
Fig 4 is a Nichols plot for an open loop transfer
function when a PID controller is used,
Fig 5 is a Nichols plot for an open loop transfer
function when a PID controller, time delay and notch
filter is used,
Fig 6 is a graph showing a simulation result at a step in
the reference for the total energy in the converter leg,
Fig 7 is a graph showing the simulation result when the
current changes from 0.1 pu to 1.0 pu in the converter
leg,
Fig 8 shows a Nichols plot for a balance controller
according to another embodiment.
, CA 02772713 2012-02-29
8
DETAILED DESCRIPTION
The invention will now be described more fully
hereinafter with reference to the accompanying drawings,
in which certain embodiments of the invention are shown.
This invention may, however, be embodied in many
different forms and should not be construed as limited to
the embodiments set forth herein; rather, these
embodiments are provided by way of example so that this
disclosure will be thorough and complete, and will fully
convey the scope of the invention to those skilled in the
art. Like numbers refer to like elements throughout the
description.
In the description in the following, continuous variables
are used, corresponding to the simplifying assumption
that the arms have infinite number of submodules that are
switched with infinite switching frequency.
Fig 1 shows an M2C (Modular Multilevel Converter) having
a phase leg 7 that comprising an upper arm 5 and a lower
arm 6. Each arm 5, 6 comprising a number of serially
connected submodules 9. Each submodule 9 comprises a
switchable capacitor. An AC output 8 is connected between
the upper and lower arms 5, 6. Although only one phase
leg 7 is shown here, the M2C comprises one phase leg 7
for each phase, i.e. three phase legs 7 for a three phase
system, where each phase leg 7 comprises upper and lower
arms 5, 6 comprising submodules 9.
An upper DC source common bar (in this case DC+) and a
lower DC source common bar (in this case DC-) for each
phase is provided. It is to be noted that the upper and
lower DC source common bars can switch polarity.
CA 02772713 2012-02-29
9
Ideally the capacitors keep a constant DC voltage and the
AC terminal voltage is controlled by varying the number
of inserted modules in the upper and lower arms. If the
voltage between the DC bars is constant this obviously
requires that, in average, the total number of inserted
modules in the two arms remain constant. The arm
inductors however will limit the rate of change of the
arm currents, making it possible to accept minor short
deviations from this condition.
Now once the context is presented, define the insertion
index, nx(t), for the arm x to be the ratio between the
inserted number of submodules and the total number of
available submodules in the arm. The arm voltage then
becomes
(1) ti((t)=rij,(t)ticff(t)
In a simple approach the number of inserted modules in
each arm can be generated by the modulator much in the
same way as in PWM modulation for conventional VSCs.
Then, in order to generate an inner AC voltage with
amplitude ev the insertion indices for the upper and the
lower arms become
u c(Jkt)=n 1,t)u 1, /IA= n i(t)u
¨ t'n cos u; t
n t =1 + In cos co 1
n(Jt)
2 2
ev
in=
(2) 2
When the simple modulation approached described above is
used and the converter is loaded on its AC side the
desired waveform will be distorted due to the ripple in
the capacitor voltages that will be created when the load
, CA 02772713 2012-02-29
current passes through the converter arms. Specifically a
strong second harmonic current will circulate through the
converter leg and the DC side and/or the neighbor phases.
This undesired second harmonic current increases the peak
5 of the arms currents and causes extra losses in the
converter arms.
The problem can be solved by generating the insertion
indices for the arms, nu and nL, in other ways. Such
methods would aim to:
10 = eliminate and/or control the harmonic current in the
converter arms
= for each arm control the total energy stored in all
capacitors in that arm which is equivalent to control the
total voltage of all capacitors in the arm
= thereby control the total energy stored in the phase
leg as well as the balance between the upper and the
lower arms in the phase leg
According to the invention the insertion indices nu(t)
and nL(t) for the converter arms are being derived in
real-time according to the following procedure
= the reference for the converter inner voltage relative
the midpoint of the DC link is given in the form ev(t)=ev
cos colt; this reference typically is delivered by an AC
side controller operating on AC quantities like output
voltage, current or flux; the converter circuit
parameters like arm resistance and inductance may be used
by the controller
= the desired arm voltages ucu(t) and ucL(t) are
calculated as
, CA 02772713 2012-02-29
11
uthr (11 u( L (t )=¨+e
(3) 2 , 2 - thjf
where up is the voltage between the DC rails and udiff(t)
is a control voltage that is created by the control
system that will be described later in the this memo
= the total capacitor voltages, and of all
capacitors in the upper and lower arms respectively, are
measured or estimated as will be described later in this
memo
= the insertion indices are calculated as
,
. ucu ul,t) =
nul,t)= õ n19)= ______________
(4) u((/It)" u(7 )
According to the invention there are two different ways
to create the variables udiff(t), ti,(0 and te,,(t).
First approach, feedback control
In this approach the sum of the capacitor voltages in
each arm , and 717,(0, are measured using sensors in
the submodules. If the voltage sharing between the
modules is assumed to be even the total energies in each
arm can be calculated as
s C arm E C arm = ,2
W r(ikt )= w." ' __ ua
(5) 2 l= " 2 '
where Carrn=Csubmod/N. Alternatively the energy in each arm
capacitor can be calculated individually and the total
energy for each arm then can be created by summing the
energies in all submodules in each arm. The voltage
reference component udiff(t) is created as the sum of the
CA 02772713 2012-03-01
PCT/EP 2010/062 923 - 09-02-2012
9. FEE. E12 11:20 AB AS GF_C / IP 'A 570 '. 2
12
output signals from two independent controllers
uoiff (t) =1.21iff: (t) +ucliff2(t ) .
The first controller has a reference for the total energy
in both arms of the phase leg. The response signal is the
measured total energy AlQ)+17,1(t), Which may be filtered
using e.g. a notch filter tuned to he frequency 2el (m:
is the network frequency) or any ot er filter suppressing
the same frequency. The error, i.e. the difference
between the reference and the respo se signals, is
connected to a controller (normally of type PIP) that has
the output signal udinm(t).
7he second controller has a referen e for the difference
between the energies in the arms in the phase leg. This
reference typically is zero, meanin that the energy in
the arms in the phase leg shall be .alanced. The response
signal is created as the measured ,Z,1,(E)-Fq(i), filtered
by a notch filter tuned to 01 or any other filter
suppressing the same frequency The error is brought to a
controller (typically of P type), which has an output
signal adi2f2- The contribution to the total voltage
reference udiff(t) ía obtained by multiplying tleit42 by a
sinusoidal time function cos(eitt)r), which is phase-
shifted relative the inner voltage reference by the nngle
* given by *-.... Lift t joiL), where J. and L are the
resistance and inductance respectivfly in the converter
arm. Thus: 1
1
() u46.(t)=-u47,(r)+cpcos(toit+op) i
The first approach to stabilisationof the converter
accOrding to the prOCedure describe! in this section is
described in more detail in Appendi 1.
i
I
1
1
i
1
1
____________________________________________________________________________ _
Ilor7G902.20121020-49-0902.201210.2324
Thispige26of23wascwOetedat0902.2012102314
;eceived at the EPO on Feb 09,2012 1012324. Page 26 of 28
AMENDED SHEET
CA 02772713 2012-03-01
PCT/EP 2010/062 923- 09-02-2012
9. FEE, 2)12 11:20 A;E( A3 ULC / iP ',O. 570 ). 27
13
Remark 1: The reference for the AC Lde inner voltage may
comprise a minor third harmonic vollage component, which
is used to increase the available o tput voltage level in
11
a 3-phase converter. This does not impact significantly
on the behaviour described,
Remark 2: A third reference compone t may be added to the
control voltage udiff(t) . This Campo ent has the purpose
of intentionally creating a second armonic current in
the arms in order to increase the a ailable output
voltage for loads with certain powe factors.
second approach, open-loop Control
In this approach the AC side curren ivml is measured.
Its fundamental frequency component is extracted with
amplitude and phase relative the reference inner voltage
evw for the converter. Thus the AC'side current can be
written as
'
( 7) s ;v t (tH Vcos(Eo I i-k-ril
Assuming that the converter shall operate ideally La
steady-state, 1.5. it shall produce undisturbed AC output
voltage and the upper and lOwer arm shall carry half the
AC output current each, it is possitle to ealculate the
ideal derivative of the energies in each arm. The result
is
dry' 7
(8) dr 4
1 . 1
¨='her [,,,,p,COSeilit¨t-T¨Rivu,Te05 tiy +9)¨-.-cos1,2 ct.),t +(p)
75 (q dt
-
(y.3 022012 10-20-49 - ,70 02 2012 10.n.24. This page 27 of 28 was cornpleted
at 09 02.2012 1023:19
eceived at the EPO on Feb 09, 2012 1023:24. Page 27 of 28
AMENDED SHEET
CA 02772713 2012-03-01
PCT/EP 2010/062 923 - 09-02-2012
9.FE,S.212 11:21 A iS GFC i" IP \O. '--,7()
14
where 1,00 is a pc our:ent circulating through the two
$eries-connected erns and the DC Su ply
=
7 Loh'''. li,Ircos(p
OUP- ="-. ___ I Z
(10) u0I-Nup-4Ror;7oosip
When there is only a DC eirCulatin current latro then also
5 the control voltage c4gf(t) becomes a DC voltage wlth the
value u(t)= so that (3) becomels
uctr(t)=1 2-ar(i)-Rlop
(11) ucL(t)=)'+ v(t)-Aloos
2 2 '
moreover, equations (8) and (9) ca be integrated, each
with a freely selected integration oonstant, so that
al'7(flPsintuiti-(44;-Rlay.)ir siaitti-LP)-"sin(2wis+w)
10 (12) 2ro, 8cui
Fit(t)=Pqa+aid wit -(5A -it 1 ein (215.0+9 )
(13) 2 2w, I 4w1
Thus the instantaneOuS energies in each arm can be
calculated in real-time knowing only the references for
the inner converter voltage and the actual AZ current.
15 The integration constants are the references for the
desired average energy in each arm to the phase leg.
But if the energies are known then lalso the total
capacitor voltage in the arms are know due to the
connection equations
C. Cw1õ
2C (14)
.tiorr 09 022012 10:20:49 - CO2.2012 10:2324. This page 23 of 28 wa5
completed:1109.02.2012 102324
7eceived at the EPO on Feb 09, 2012 10:23:24. Page 28 of 28
AMENDED SHEET
, CA 02772713 2012-02-29
Now the insertion indices valid for the desired steady-
state operation can be calculated using equation (4).
Given these insertion indices the energies in the upper
and lower arms converges to the reference values given as
5 free integration constants in (12) and (13). Normally
these values are selected equal for both arms so that
balanced operation is obtained. The value of the energy
reference is selected to give the desired total capacitor
voltage in each converter arm.
10 The second approach to stabilisation of the converter
according to the procedure presented in this section is
described in more detail in Appendix 2.
Remark 1: If a third harmonic voltage component (to
increase the available voltage in a 3-phase converter) is
15 added the formulas for the energies in the upper and
lower arm will change somewhat. However the principle
described in this paper still can be applied.
Remark 2: If even order harmonics are intentionally added
to the circulating current the formulas for the energies
in the upper and lower arms will change somewhat. The
principle described above however still applies.
APPENDIX 1: Description of feedback control system
A purpose of this appendix is to disclose in more details
the control mechanisms (modulation) that can be utilized.
The targets are
- to derive a control strategy that provides main circuit
stability
- to produce as high AC output voltage as possible with
very low harmonic distortion
, CA 02772713 2012-02-29
16
- to control the DC voltages of the capacitors in the
modules
Continuous model
It is of course possible to investigate the M2C converter
by simulation. This approach however seems to be quite
cumbersome in the sense that it involves detailed models
of the arms (with tens of semiconductor devices in each).
A lot of data will be generated making it more difficult
to extract useful results.
Another approach, which will be followed here, makes use
of a modulation principle that has been proposed by in
DE10103031, in which a selection mechanism is used to
determine which individual module that shall be inserted
or bypassed when the number of devices in an arm shall be
changed. The selection is made in dependence of the
direction of the arm current (or phase current) and a
comparison of the DC voltages in the congregation of
modules in each arm, from which the modules having the
highest and lowest voltages are identified.
Simulation has shown that this mechanism successfully
keeps the DC voltages of the module capacitors quite
close to each other, even for low number of modules in
each arm (say e.g. five per arm). This functionality
seems to remain even if the total switching frequency is
low (a few hundred switchings per second for each
semiconductor device).
Now it is assumed that this mechanism is in use and that
accordingly there is no need to look at the DC voltages
in the individual modules any more. The modulation
process then can be described in terms of the total
, CA 02772713 2012-02-29
17
collective energy in each arm. As the total switching
frequency (for all modules in each arm) becomes quite
high continuous modelling can be used. The continuous
model is a lot simpler to grasp than the detailed model
and it is an ideal base for understanding the principles
for the function of the M2C converter and to formulate
control laws for different control aspects.
Due to the assumptions static relations exist between the
total capacitor energy in the upper and lower arm, w</ and
/57( 1', and the corresponding total voltage of all capacitor
modules in the arm, luand 111. Namely, if it is assumed
that the energy is evenly shared between the modules,
this relation becomes
lc c I 2
W<F, =
2 N 2N ' ueõ = WC,71
)
1C7141
Pir r'; = N =C IE2 - "' 12N WE
2 N 2A "CL \I __ , CI
(Al) ,
where N is the number of modules per arm and C is the
capacitance per module. In the following we will use the
quantity 'arm capacitance' Carm defined as follows
=-17
(A2)
Then
= __
TA/
,
I C 41/171 (õE y
i2 ___________________________________
\2
W(1- 2 kuci / 14(/.=
I
(A3) <=>
CA 02772713 2012-02-29
18
Derivation of the continuous model
The electrical circuit representing the phase leg of the
M2C converter is depicted in Fig 2. The inserted
capacitor voltages, Ucu and UcL, are created from the
total capacitor voltages, 4(land 1141, respectively, by
applying the insertion indices, nu and nL, which are
controlled by the control system.
n, 0
n, 0 < n <1
(A4)
In the following, however, the main circuit model will be
formulated using the real voltages as variables. If the
total capacitor voltages, limand UrL are measured, the
corresponding insertion indices can always be obtained
from (A4).
The capacitor modules serve as controlled electromotive
forces in the circuit. Let the total energy in the
capacitors be "'and 141
:Lin the upper and lower arms
respectively. Inspection of the circuit model in Fig 2
immediately yields
d147(1,
iuum
dW
(A5) &
In order to gain some more insight it is helpful separate
the arms currents in two parts. One part emerges from the
AC current, which naturally separates into two halves,
one passing though the upper and one passing through the
lower arm. The deviation from this "ideal" condition is
CA 02772713 2012-02-29
19
described by a difference current idiff which passes
through the series-connected arms and the DC source.
Define
'd
2 Iff
== _ U L
= ¨
2
(A6) 2
The circuit in Fig 2 now gives the equations
________________ R(i + ldi L __ + - =0
2 2 2 dt dr ,
r =
( I di fi , did,
_uõ - R - L +uõ,-uõ =
(A7) 0
2 dt dt
Adding and subtracting the equations give the results
= ¨ R = L div
2 2 2 dt
di
L + Rid U11(A8) dt ff ¨ 2
2
These equations show that
- the AC voltage only depends on the AC current iv and
the difference between the arm voltages ucl, and ucu
- the arm voltage difference acts as an inner AC voltage
in the converter and the inductance L and resistance R
form a fix, passive inner impedance for the AC current
- the difference current "cliff only depends on the DC link
voltage and the sum of the arm voltages
- the difference current 'cliff can be controlled
independently of the AC side quantities by subtracting
the same voltage contributions to both arms
CA 02772713 2012-02-29
Define
U -uct,
e, -
U))1 = -e) -udg
=
-11(-1.-"(
kilt/
(A9) 2 <=> - 2 + -
where ev is the desired inner voltage in the AC voltage
source and udiff is a voltage that controls the difference
5 current idiff.
Then (8) becomes
R L
= er -- --
2 ' 2 dt
di
L " + Ri, =udar
(A10) di
Inserting equations (A6) and (A9) in (A5) yields
dw( 7i
di , + .d )(
, 2 ,-2 - e,.
e, )(A11) 2 "g, 2 -um
10 It makes sense to investigate the total energy stored in
all capacitor banks in the whole leg and to examine the
balance between the energy in the upper and the lower
arm. Define
Wz _____________________________________
+147
= c
2
wA _wz
(Al2) c = cu (f, 2
15 The result is
dr4/(T
di = õ 2u,õ)7.)idõ. -eriv
dW
/m.
____________ = -2edd,g +
(A13) di \2
CA 02772713 2012-02-29
21
Equation (A13) indicates that the total energy in both
arms as well as the energy balance between the upper and
lower arms can be controlled primarily by ithft, which is
in its turn controlled by udiff through (A10).
The term in the upper equation in (A13) is recognized as
the instantaneous power delivered to the AC side.
(A14) Pt --etiv
Steady state solution
We shall look at the special case where the AC emf and
current is given. Thus let
e = e cos cot
"="; COSIWI+0
(A15) v r ( )
Assume that there is a solution where the difference
current idlft is a pure dc component. Thus
(A16 id tt (0= lag
) '
Then, according to (A10)
(A17) u d'ff
The derivative of the total and difference energies the
according to (A13) become
cit _________ = kup 2 Rid, _)1thir 2 r
rosyo+cos(2o1+4
dWA u,,
____________ = 28,,i,41 cos cot + cos(ot + yo)
(A18) 2
From (A18) some observations can immediately be made
CA 02772713 2012-02-29
22
- the derivative of the total energy contains only a
constant and a component having double network frequency
- the derivative of the difference energy only contains
components having network frequency
Steady-state condition requires that the constant
component of the total energy derivative disappears so
that
(uõ ¨2 Ri es,21} cos = 0
i = ________________ cosc,
diff
(A19) uõ + Vu;) -4R cos
With this difference current the remaining term becomes
dW ev.
COSOCOt+0
(A20) dt 2
The steady-state average energy can be freely selected so
that the total energy in steady-state becomes
147 P171) eviv
((t)= ________________ sink2cot+0
(A21) 4w
The expression for the difference energy In (A18) can be
directly integrated, also with a freely selectable
integration constant (which normally shall be zero)
Ur) -
d)
i
+2
i 4
1I(1 = !I/ft __________ sina
_____________________________________ Mn(cto+0
(A22)
The investigation shows that
, CA 02772713 2012-02-29
23
- solutions of the desired type, i.e. with a difference
component having only a DC component, exist with freely
selectable energy levels in each arm
- the steady state solutions for the energy time
functions contain only a double frequency component in
the total energy and a fundamental frequency component in
the difference energy
Linearised model for control studies
Let us go back and linearise the equations (A13) around a
steady state point as described in the preceding section.
Assume that the DC link voltage is constant. The
dAlt
_____________ = õ _2R1 ) id,fl ¨ 2 ithif A d, AP,
dt
_____________ -Ai&if - i Au ¨2;thff Ae +
t &ff f
(A23) di
2
Further the differential equation that governs ithif as
function of udiff applies so that
dAt,
__________________________ + RAid,ff = Audiff
(A24) di
Stability requirements
When the AC side current is stiff (A23) reduces to
dAW,
di ___________ = ku 2Ri11)Ai,,1 -2i AU j+
dAW,
_____________ = 2e, Aid,õ -
(A25) di
The linearized equations (A25) show that any control
system, which makes the sum of the inserted voltages, uci,
and ucu, perfectly match the voltage up on the DC side,
CA 02772713 2012-02-29
21
u 0
i.e. makes d'ff , also makes the difference current
become zero causing the derivatives of the energies in
the arms to vanish. The main circuit in the converter
then is marginally stable. Thus is not sufficient to
select the inserted voltages in (A9) according to the
desired ev , but an udiff that creates stability must also
be provided.
Control law for the total capacitor energy
The equation for the total energy equation can be
formulated in the Laplace domain
(uõ -2Cõ)
Au/if(s)- AP' (s)=
s R+sl,
_____________________ , Au(s) (s)
'
(A25) (R + so s
Applying a proportional gain in an energy controller
(A27) K {M *'(s)- (s)- W(s)}
yields
[un - 4 Ri du, -
wt1 (s) = __________ A (s)¨
(R + sL) s + [tt, ¨ - 2Li111s ]K
R + sL
Ap,,(s)
(A28) (R + sL)s +[u õ - 4R4,õ- 2 L4,õ s IK
The poles in the above transfer functions are mainly
determined by
= j uDKp
(A29) s2'2L +ul)Kp =0 L
CA 02772713 2012-02-29
The control system is investigated for a converter leg in
an example converter with the main parameters given in
Table 1.
3-ph rated power 30 MVA
rated frequency 50 Hz
line-line voltage 13.8 kV rms
rated phase current 1255 A rms
arm capacitance 500 pF/arm
arm inductance 3 mH
arm resistance 100 mQ
Table 1: Example converter main data
5
Fig 3 shows the Nichols plot for the open loop transfer
function in (A26) with the proportional gain K=0.001
V/J. The curve is almost independent of the active load.
As expected is the phase margin at 90 rad/s quite small,
10 which means that the response will be quite oscillatory.
In Fig 4 the Nichols plot is shown when a PID controller
is used.
The selected transfer function is given by
F(s)=K1, 1+K ST
S l+sT
(A30) //
15 with KF=0.002 V/J, K1=20 s TD=10 ms, TF=2 ms
It has been shown that the total energy response signal
contains a frequency component with twice the network
, CA 02772713 2012-02-29
26
frequency. This component can be removed from the
controller response using a notch filter. Further it is
advisable to assume that a delay occurs in the measured
total energy (total capacitor voltage). Fig 5 shows the
corresponding Nicols's diagram where the notch filter and
the time delay have been included.
The transfer function in Fig 5 is
Fl(s)= K p( 1 + K + sTD 32+0.4
(A31) s 1 + + 24'a)
(2,) s + Oco,Y
with KF=0.002 V/J, KT=20 s-1, TF=10 ms, TF=2 ms, Tdel=1 ms,
=0.05
Fig 6 shows the simulation result at a step in the
reference for the total energy in the converter leg.
Equation (A28) shows that the energy control system
having only a proportional feedback will have a static
error
(A32) -R
Ap,
For the values in Table 1 together with K=0.002 V/J this
energy dependence becomes approximately 0.002 J/W. Each
leg of the converter handles about 10 MW causing the
energy drop to be about 20 kJ (out of about 312 kJ) per
leg.
Fig 7 shows the simulation result when the current
changes from 0.1 pu to 1.0 pu in the converter leg.
Control law for balancing the capacitor energies in the
arms
, CA 02772713 2012-02-29
27
The general differential equation governing the balance
between the energies in the upper and the lower arm was
derived in (A13)
dWA
(
= +u
(A33) & 2
It was linearised in (A23)
dAIV
= ¨2e, Aithfl ¨ 21,,frAe, + (141) iif) A
(A34) dt 2 d,
If we consider linearising around the steady state
solution defined by (A15) the linearised equation becomes
_____________ = -2 cos cot Aithrf -1, cos(cot + yo)Auro -2õAe, + R/õ,, Ai,
(A35) & ,2 2
Assume first that the AC side quantities are constant.
Then
dAT/V,A
_____________ = 2e,coscotAithrf-i, cos(col +y)Atia,fi
(A36) di
Further assume that the controller produces a fundamental
frequency sinusoidal signal with phase relative the
inner emf in the converter leg
(A37) 14,,,,t (t) = (t)cos(cot +
Using the quasi-stationary solution to (A10) yields
cliff (1) __ COS(COt ¨ ri)
(A38) + (coL)2
ij = arg(R + j coL)
Inserting in (A36) we get
, CA 02772713 2012-02-29
28
dAY'r2e,
c ______________________ cos(cot)cos(ng + - 0-1v cos(cot +c)cos(tot +)
An,
(A39) /
di V-R2 (AY
The products of the cosine functions in (A39) are DC
quantities and terms with the double network frequency.
These components are
a' (IL; V v
COS
\I R2 +coLY 2
¨ev =
v
a ¨ ______ cosi:.2wt ¨ cos(2 t (p)
(A40) VR2A-kuL)2 2
The relation between the two terms at various frequencies
has been found to show that the first term dominates
completely even for operating frequencies down to 5 Hz.
Therefore it is sufficient to consider the first term.
The maximum DC component then is obtained when
(A41) =r7=argR+ja,L)
With this selection of the argument for the inserted
difference voltage we get the simplified formula
dAff"
________________________ And,ff
di R2 +(coif
(A42)
A proportional controller is sufficient to control the
balance between the energies as the function is indeed
just an integrator. However, the measured difference
energy has a strong fundamental frequency component,
which should be eliminated in the response to the
regulator. Thus the transfer function in the balancing
controller is given by
, CA 02772713 2012-02-29
29
S F-1(S)=K ______________ I
PL' 2 2
(A43) s+2Cco,+(u,
Fig 8 shows the Nichols plot for the balance controller
with parameters according to Table 1 and with control
parameters KE,A=-0.005 V/J, tpel=1 ms, =0.1.
Fig 5 shows that the closed loop for the energy
controller has unity gain up to about 300 rad/s and that
it amplifies frequencies in the range 100 - 200 rad/s
with more than 3 dB. Therefore the gain in the balance
controller has been kept low for these frequencies in
order to avoid interaction between the two controllers.
Appendix 2: Description of open loop control system
The aim of the investigation is to describe an M2C system
where the modulation operates in open-loop mode. The
meaning of the name "open-loop" in this context is that
the modulation system does not measure the total voltage
of the capacitors in the phase leg arms. Rather these
total voltages are estimated in run-time using the
desired AC emf and the measured AC current. The reference
for the inserted arm voltages are obtained assuming that
the instantaneous AC emf and AC current are steady state
values. Further it is assumed that a voltage sharing
system is provided to distribute the total arm voltage in
each arm evenly between all modules that constitute the
arm.
Steady state analysis
The starting point is that the converter produces a
sinusoidal emf
CA 02772713 2012-02-29
(B1) ev=evcoswit
and is loaded with a sinusoidal phase current
(B2) ir=ircos(wit+q9)
Under ideal conditions the arm currents only contains a
5 DC component I ditto so that the arm currents become
/v
cosyuit-1-kp)+'idviov
( B3) ' z COSi, w it + IP) dip
l' 2
When the difference current is iduio the difference
d¨Ri -
voltage becomes u
¨ diffl) so that the arm voltages become
/ID ui)
ucu=-2¨evcoscuit¨ ditio
R; u ,= ¨2 +evcostel t ¨ R dujo
(B4)
10 The derivatives of the arm energies are
dW __________ = uC dW
i II11 _________________________ = u i
(B5) di di
Inserting the expressions in (B3) and (B4) yields
dW("1,1/ = = up A A evircosy9
dt=( 2 Ri )i
4
I-
.11p - = e
¨ e v thfiocos (.4.) t ¨ R dzio ¨v cos (a)! t +9))¨ cos (2 re, t + )
(B6) 2 ' ' 2 4
ivcos(,0 ,
R idv).)idifto
dt __________ = 2 4
i 'e i V '
+eµ v id(11000S COI t COS (WI t +(p)¨ cos (2 to, t +Li))
(B7) " 2 2 4
15 In steady state the DC term must be zero. This condition
allows us to determine the DC component to
, CA 02772713 2012-02-29
,
31
,. erzvcoscp
i ¨ _____________________________________
(LOP ¨
i_j 2 A __ A
(B8) up,- v ui, -r Re. v iv cos Lp
Thus, in steady state, the energy variations are
.
dW(21, , _ ,k v _i
V .
______________________ = e V i,WICOS IC , t ( Ri ) cosi to
t+ip11¨ cosi,2noit+cp)
(B9) dt . = 2 d00, 2 % 1 ,
4
,,,.õ- ,
1.4 rr (7 , UD A , i v , ' evir
,
______________________ ' =-P 'e i cosco t¨i Ri 1
(B10) ________________________________________________________ dt cos(wit-
H49) cos(2teitd-Lp)
v thito 1 ' 2 Off . 2 4
These formulas can immediately be integrated to obtain
the instantaneous energy variations. Note that a freely
selectable integration constant appears in each
expression. Thus
..,.. ,,= ...y 'et' ithijo
sin to t +1 , Ri
u ) - = i l' ,
Witõ(t)= Wcuo do) __ sin ( re, t + 49 )
sin (2w1 t+)
1 '
(B11) to, 2 .. 2 co,
8w1
, .
2. . .
e v .
6v
_________________________________ sin re, t¨ (u1-) ¨R,,,11)--sins ( rii, t +
ip ) ¨tv sin (2, t+ 49 )
(B12) ' Oa, 1 2w1 ' 4 a:, =
The total capacitor voltages now are given by
12. WI u 1'11 12 WI"' , (t)
, /42u/ (1)= c,
(313) C arm - GM
and they can be used to determine the insertion indices
in run-time according to
, u(.ul,t) /foist)
n = ________________________________ ti = __ "
(1I IL, , z , N 1,(t1 , + 2 = ,
(B14) ucti1,0 11 I 1, I
C Ls
Open-loop control
The idea of the open-loop control mode is to
- measure the AC terminal current
CA 02772713 2012-02-29
32
- extract the amplitude and phase relative the created
emf in the converter as in (B2)
- perform the calculation as described above
- utilize the so obtained insertion indices according to
(B14) in the converter
The invention has mainly been described above with
reference to a few embodiments. However, as is readily
appreciated by a person skilled in the art, other
embodiments than the ones disclosed above are equally
possible within the scope of the invention, as defined by
the appended patent claims.
