Modelling and Control of Multi-Terminal MVDC Distribution Network
Modelling and Control of Multi-Terminal MVDC Distribution Network
Distribution Network
Patrobers Simiyu†, IEEE Student Member; Ai Xin IEEE, Senior Member; Lawrence Bibaya; IEEE
Student; Vedaste Ndiyishimiye; George Adwek, Girmaw T. Bitew,
S 3
2 ^ eqiq j eqid ed iq `
ed id
(7) It is observed that DC voltage can be controlled the same
way as active power by regulating id to its set point i*d as
The PLL effects the transformation of the AC voltage eabc at shown in Fig. 6. In this case, the compensator process the DC
the PCC to edq. The voltage vector oriented control (VOC) is voltage error in order to determine the set point i*d of the inner
used to decouple the active and reactive power in which the it current controller via a feed-forward compensation. The inner
adjusts the rotational speed of the dq-axis frame such that eq = 0 current controller ultimately regulates the DC voltage of the
by aligning the ed vector with the d-axis at steady-state. Thus, plant. The DC voltage controller alongside the reactive power
the active and reactive power injected to or absorbed from the controller in Fig. 5(b) forms the Vdc-Q control mode of a VSC
AC system at the PCC is given in (8). terminal.
3 2 P*
P ed id o id* Equivalent Inner Current
2 3ed DC Voltage Controller Controller
idmax G(S)
(8)
3 2Q* + + id* + ic
Q ed iq o iq* PI id 2ed 1 Vdc
2 3ed V*dc 3Vdc Cs
- + -idmax -
Therefore, based on (8), the active and reactive power at the iL
PCC can be controlled by id and iq respectively which Vdc
2 Pdc
constitutes the PQ control mode of a VSC terminal. Fig. 5 Pdc 3ed
shows the active and reactive power controllers. The
controllers’ outputs have limiter functions for restricting the
Fig. 6. VSC Controllers and the Plant
reference currents within permissible values. The active and From Fig. 5 and 6, it can be observed that id is involved in
reactive power flow can be effectively and independently both active power and DC voltage controllers. When the two
controlled by fast reference tracking of their corresponding are combined, a DC voltage droop controller is derived as
references commands; idq≡ idqref. shown in Fig. 7. The controller regulates the DC bus voltage
iqmax linearly by changing the active power reference in proportion
idmax
to the DC voltage, regulating power to its reference value
+ + PI
PI id* Q* iq* simultaneously with DC voltage. The slope of the droop
P*
control characteristics is the droop constant, Rdroop and its
- -idmax - -iqmax
i inverse the droop gain, Kdroop. This parameters determine the
controller sensitivity to changes in DC voltage with variations
P Q (b) in power. With the use of the PI controller, the steady-state
(a)
error of the DC voltage controller equals zero as shown in
Fig.5 Active and Reactive Power Controllers (12)[17].
Active power can be regulated as an intermediate variable by
id for ultimately controlling DC voltage to its corresponding P* P Kdroop (V * V ) 0 (12)
reference in the DC grid. The power balance equation between
the AC input and the DC-side in Fig. 3 is given in (9). The
power losses of the VSC are ignored as the reactor reactance
has very small resistance.
P Pdc Pc 0
(9)
3
e i Pdc Vdc ic 0
2 d d
¯ Z dt ¿
°X d (i ) ½
°
X X L id ¾ ª¬et1q riq Vq º¼
q Inner Vdq_ref
® (26) Vdq_ref(pu) PI
Current
abc
PWM
°
¯ Z dt °
¿ Controller
dq
Vabc
Vdq PLL
C. AC Load and Load-Side VSC
A constant impedance-constant current-constant power (ZIP) abc Vabc
Idq dq Iabc
static AC load model is connected to the AC system’s bus
Fig. 10 Load-Side VSC Controller
which in turn is linked to the MVDC grid via a 3-phase VSC
NPC. It is a 3-phase block widely used for implementing a D. Battery and Buck-Boost Converter
3-phase balanced load as a series combination of RLC elements The lead-acid battery model is adopted for the study as it is
as shown in Fig. 9. The equations describing the AC load model the cheapest on cost per kW basis, mature and widely used
are described in [16] using a polynomial in (27) and (28), where battery energy source storage (BESS) technology. When
Pi, Qi and Vi are the real power, active power and voltage at the analyzing the dynamic behavior of the lead acid battery, the
ith bus respectively. The variables with a 0-subscript denote the electrical equivalent models are preferable to the mathematical
nominal values while those with 1,2 and 3 represents the and electro-chemical models. The ideal, linear and Thevenin are
fraction of the constant; Z, I and P of the AC load respectively, the common equivalent circuit models[31]. The linear model
hence, Pi1+Pi2+Pi3=1 and Qi1+Qi2+Qi3=1. shown in Fig. 11 is the modest choice representing the battery as
ª
º
2 a voltage source with a series internal resistance [32],[33].
V V
Pi Pi 0 « pi1 V i pi 2 V i pi3 » (27)
¬ i0 i0
¼ Ibat
+
ª
º
2
V V
Qi Qi 0 « qi1 V i qi 2 V i qi 3 » (28) EMF +
Rint
Vbat
¬ i0 i0
¼ -
Ud T1 Uo
- +
+
PI IL_ref PI Uo_ref Dref
Ud_ref 1/(Ud) PWM
Fig. 13 Buck-Boost DC-DC Converter
- -
Assuming a steady-state operation, the voltage across the +
P VI YES
ΔI/ΔV>(-I/V) ΔI>0
YES
dP d ( IV ) dI dI (43)
I | I V NO
N NO
dV dV dV dV
Increment Vref Decrement Vref Decrement Vreff Increment Vref
The conditions for the true MPPT are given in (44) [42].
dP
0
dV
dI
I V
I(t-Δt)=I(t)
0 (44) V(t-Δt)=V(t)
dV
dI I
Return
dV V Fig. 17 The Incremental Conductance Method Flow Chart [41]
Equation (44) indicates that the MPP can be obtained by
comparing the instantaneous conductance (I/V) to the III. RESULTS AND DISCUSSION
incremental conductance (dI/dV). Equations (42) and (44) can
A. MVDC Distribution Network Model in PSCAD/EMTDC
be used to derive the incremental conductance algorithm. The
relations before the MPP, at the MPP and after the MPP are The 10kV radial MVDC distribution network in Fig. 1 is
respectively given in (45) [41],[42]. modelled and simulated in PSCAD/EMTDC. The converters at
the MVAC I and II as well as that at the AC load are
dP dI I
! 0 or !0 galvanically separated from the AC side by an equivalent
dV dV V transformer reactance. In the PV system, the incremental
dP dI I
0 or 0 (45) conductance algorithm is used for MPPT. The key parameters
dV dV V
dP dI I of the network are given in Table 1.
0 or 0 The performance of the MVDC network is assessed for
dV dV V
The expressions in (45) define the three operation zones for master-slave, voltage margin and droop control. In the
the incremental inductance method and the flow chart is given investigations, different operation situations such as start-up
in Fig.17. In this method, the MPP is determined by comparing behavior including step changes in PV output at 2s, normal
the incremental conductance (IΔ/VΔ) to the PV array operation and operation under some contingency are
conductance (I/V). When (I/V = IΔ/VΔ), the output voltage at considered. At t=3s and t=6s, an LLL-G fault is applied at the
Vref is the VMPP and the controller maintains this voltage until AC load and MVAC I respectively for 0.1s whereas at t=9s, the
MVAC I is disconnected.
the solar irradiation changes hence change in I and the
In the master-slave control, the MVAC I is the master while
process is repeated. The MPPT regulates the PWM control
MVAC II the slave. The PV is operated at MPPT while the AC
signal of the boost converter until: (∂I/∂V) + (I/V) = 0 is
load controller regulates its terminal output voltage. At the
satisfied by adjusting its duty cycle. In this method, the peak
battery, the charge controller continually monitors the charging
power of the PV array lies above 98% of its incremental
and discharging regimes though the emphasis for the study is on
conductance[40].
the charging operation. A two-stage voltage margin controller
reviewed in [10] is applied in which the DC voltage control
priority is given to MVAC I where the lowest and highest DC
voltage values for the margin are 9.8kV and 10.2kV
respectively. The minimum, reference and maximum active
power are given as -2.8MW, -3.0MW and 3.2MW respectively
for MVAC II while the reference power for MVAC I is 6MW.
The PV, battery and AC load controller remain the same as in
8.0
6.0
possible for the DC voltage to be fast-tracked to the set value 4.0
2.0
than the other two. 0.0
There are also startup transients mainly within the first 1s in 20.0
AC Grid I(MW)
2.0
margin resulting in a 3% DC voltage increase in the MVDC 1.0
0.0
-1.0
grid. It can be observed that the step changes have less than 1% -2.0
-3.0
increase in DC voltage when droop control is applied. Besides, 0.0
AC Load (MW)
Vdc (kV)
10.0 10.0
8.0 8.0
6.0 6.0
4.0 4.0
2.0 2.0
0.0 0.0
AC Grid I (MW)
20.0 20.0
AC Grid I (MW)
0.0
0.0
10.0
AC Grid II (MW)
9.0
AC grid I (MW)
-10.0
-8.0
4.0
3.0 5.0
PV(MW)
2.0
PV (MW)
1.0
0.0
-1.0
-2.0
-3.0
-4.0
0.0
AC Load (MW)
1.0
AC Load (MW)
-16.0
1.00 -16.0
Battery (MW)
1.00
Battery (MW)
-4.00
0.0 2.0 4.0 6.0 8.0 10.0
Time (s) -4.00
(b) 0.0 2.0 4.0 6.0 8.0 10.0
14.0 Time (s)
12.0 (a)
Vdc (kV)
10.0
8.0 14.0
6.0 12.0
Vdc kV)
10.0
4.0
2.0 8.0
0.0 6.0
4.0
AC Grid I (MW)
2.0
20.0 0.0
AC Grid I (MW)
20.0
-5.0
8.0 0.0
6.0
AC Grid II (MW)
4.0
10.0
AC Grid I (MW)
2.0
0.0
-2.0
-4.0
-6.0
4.0 -10.0
3.0
PV (MW)
2.0 5.0
1.0
PV (MW)
0.0
-1.0
-2.0
-3.0
0.0 -3.0
AC Load (MW)
1.0
AC Load (MW)
-16.0
1.00 -16.0
Battery (MW)
1.00
Battery (MW)
-4.00
0.0 2.0 4.0 6.0 8.0 10.0
Time (s) -4.00
Fig. 18 Start-up and Steady-State Characteristics for: a). Master-Slave, b). (b)
Voltage Margin and c). Droop Control Strategies
8.0
6.0
4.0
2.0
guaranteeing better DC voltage control and active power
0.0
balance.
20.0
However, the DC voltage droop strategy has its main
AC Grid I (MW)
6.0
4.0
2.0 to simultaneously and effectively fulfil these constraints in
0.0
-2.0
-4.0
MTDC grids. For instance, Chaudhuri et al., proposed an
-6.0 adaptive/variable droop scheme for effective power sharing in
5.0
MTDC grids [16], [19]. A modified droop control for DC
PV (MW)