Power Quality Enhanced Operation and Control of A Microgrid Based Custom Power Park
Power Quality Enhanced Operation and Control of A Microgrid Based Custom Power Park
Power Quality Enhanced Operation and Control of A Microgrid Based Custom Power Park
Ghosh, Arindam and Majumder, Ritwik and Ledwich, Gerard F. and Zare, Firuz (2009)
Power quality enhanced operation and control of a microgrid based custom power park.
In: 7th IEEE International Conference on Control & Automation (ICCA'09), December 9‐
11, 2009, New Zealand, Christchurch. (In Press)
© Copyright 2009 IEEE
Abstract − This paper describes the operation of a microgrid ics and are somewhat unreliable [7]. However, a power elec-
that contains a custom power park (CPP). The park may contain tronic converter can mitigate harmonic and unbalanced load or
an unbalanced and/or nonlinear load and the microgrid may source problems. In [8] a compensator is added with each in-
contain many distributed generators (DGs). One of the DGs in dividual DG system, which is controlled to achieve an en-
the microgrid is used as a compensator to achieve load compensa-
hancement of both the quality of power within the microgrid
tion. A new method is proposed for current reference generation
for load compensation, which takes into account the real and and the utility grid. The compensator is connected both in se-
reactive power to be supplied by the DG connected to the com- ries and shunt. The series element can compensate for all the
pensator. The real and reactive power from the DGs and the util- unwanted positive-, negative-, and zero-sequence voltage dur-
ity source is tightly regulated assuming that dedicated communi- ing any utility grid voltage unbalance, while the shunt element
cation channels are available. Therefore this scheme is most suit- is controlled to ensure balanced voltages within the microgrid
able in cases where the loads in CPP and DGs are physically lo- and to regulate power sharing among the parallel-connected
cated close to each other. The proposal is validated through ex- DG systems.
tensive simulation studies using EMTDC/PSCAD software pack- All customers in a custom power park benefit from high-
age (version 4.2).
quality power supply which is superior to the normal power
Index Terms: Distributed generation, microgrid, power quali- supply from a utility [9]. The electrical power to the park is
ty, load compensation and power sharing. supplied through two feeders from two independent substa-
tions. Both these feeders are coupled together through a solid
I. INTRODUCTION state transfer switch that can perform a seamless sub-cycle
transfer from the preferred feeder to the alternate feeder. In
MICROGRID contains a cluster of loads and distributed
A generators operating as a single controllable system to
provide power to its local area. Many of the DGs are inter-
addition, a power-converter-based device is used to eliminate
harmonic current and/or to correct voltage sag/swell [10-11].
In this paper, the concept of microgrid and custom power
faced to the microgrid through power electronic converters, are combined together to come up with a microgrid that share
which are controlled to deliver the desired active and reactive their local load power requirement with the utility. Moreover,
power to the system. The most common method used for this one the DGs is used as a compensator to cancel the effects the
purpose is the droop control for load sharing without any harmonics and unbalance in the local load such that the PCC
communication, which can be achieved by controlling two voltage is balanced and distortion-free. The other DGs in the
independent quantities – the power angle, and the fundamental system are used a voltage sources providing part of the power
voltage magnitude [1-4]. Local signals are used as feedback to requirement. The compensating DG is interfaced through a
control to the DGs. voltage source converter (VSC), while the remaining DG can
The Consortium for Electric Reliability Technology be either inertial or inertia-less. The real and reactive power
Solutions (CERTS) microgrid concept assumes a relative through the DGs and the utility source is tightly regulated as-
small system in a localized site, where the microgrid can suming that dedicated communication channels are available.
seamlessly separate from the utility grid during contingencies The efficacy of the proposed scheme is validated through ex-
where the microgrid is separated from the utility and continues tensive simulation studies using EMTDC/PSCAD software
to supply its own local loads [1, 5]. The CERTS microgrid package (version 4.2).
assumes a single point of common coupling (PCC) with the
utility. In this, the DGs are controlled locally by power elec- II. SYSTEM STRUCTURE
tronic devices incorporating the droop characteristics. An in-
teresting feature of the CERTS microgrid is that a DG can be The basic power system model with two DG sources is
inserted in a plug-and-play fashion and not all the DGs are shown in Fig. 1 in which the real and reactive power
required at any given time making the microgrid completely drawn/supplied are denoted by P and Q respectively. The mi-
flexible. A review of various microgrid projects is given in crogrid is connected to the utility grid at the point of common
[6]. coupling (PCC). In Fig. 1, DG-1 is connected through the
The power quality issues are important for DG interconnec- compensator and is connected directly to the microgrid
tions as the power electronic converters can generate harmon- through circuit breaker CB-3, while injecting P1 and Q1. DG-2
is connected to the microgrid through CB-4 and an inductor
A. Ghosh, R. Majumder, G. Ledwich and F. Zare are with the School of L2. It injects P2 and Q2. The microgrid supplies an unbalanced
Engineering Systems, Queensland University of Technology, Brisbane, Qld and nonlinear load that draws PL and QL. The utility supply is
4001, Australia.
Email: a.ghosh@qut.edu.au denoted by vs and the feeder resistance and inductance are
denoted respectively by Rs and Ls. The utility supplies Pg and III. REFERENCE GENERATION AND COMPENSATOR CONTROL
Qg to the microgrid and the balance Ps − Pg and Qs − Qg are
supplied to the utility load, which is shown as motor in Fig. 1. In this section, the reference generation and for both the
The breakers CB-1 and CB-2 can isolate the microgrid and the compensator and DG-2 is presented. Also the control issues of
utility load respectively from the supply. the compensator are presented.
The structure of the compensator is shown in Fig. 2. In this i1a + i1b + i1c = iLa + iLb + iLc (3)
DG-1 is assumed to be an ideal dc voltage source supplying a
voltage of Vdc1 to the compensator. The compensator contains Since ig and i2 balanced due to the action of the compensa-
three H-bridge converters. The outputs of the H-bridges are tor, the voltage vp will also become balanced. Hence the in-
connected to three single-phase transformers that are con- stantaneous real powers Pg and P2 will be equal to their aver-
nected in wye for required isolation and voltage boosting [12]. age components. Therefore we can write
The resistance Rf represents the switching and transformer ( ) ( ) ( )
v pa i ga + i2 a + v pb igb + i2b + v pc igc + i2c = Pg + P2 (4)
losses, while the inductance Lf represents the leakage reac-
tance of the transformers. The filter capacitor Cf is connected From the KCL of (2), (4) can be written as
to the output of the transformers to bypass switching harmon-
ics. v pa (iLa − i1a ) + v pb (iLb − i1b ) + v pc (iLc − i1c ) = Pg + P2 (5)
Equations (3), (5) and (7) form the basis of the algorithm.
From these three, the following can be written
⎡i1a ⎤ ⎡iLa ⎤ ⎡ 0 ⎤
⎢ ⎥
A⎢⎢i1b ⎥⎥ = ⎢ ⎥ (
A⎢iLb ⎥ + ⎢ − Pg + P2 ) ⎥ (8)
Fig. 2. DG-1 with the compensator structure. ⎢⎣i1c ⎥⎦
⎣ (
⎢⎣iLc ⎥⎦ ⎢− 3 Qg + Q2 )
⎥
⎦
where regulator (LQR) design. From uc(k), the switching function is
generated as
⎡ 1 1 1 ⎤
⎢ ⎥ If uc (k ) > h then u = +1
A = ⎢ v pa v pb v pc ⎥ (15)
elseif uc (k ) < −h then u = −1
⎢ v pb − v pc v pc − v pa v pa − v pb ⎥⎦
⎣
where h is a small number.
The determinant of the matrix A is given by With respect to Fig. 3, the state vector is chosen as
( ) (
A = v pa v pc + v pb − 2v pa + v pb v pa + v pc − 2v pb ) [
xT = i1 icf vp ] (16)
(
+ v pc v pb + v pa − 2v pc ) (9)
It is stipulated that the compensator does not regulate the PCC
If vp is balanced, then the following is true voltage. Therefore the positive sequence of the fundamental of
the PCC voltage is extracted as in [13] and it is used as the
v pa + v pb + v pc = 0 (10) reference for vp. Also note that the PCC voltage vp is the same
as the voltage across the filter capacitor vcf. Hence the refer-
Substituting (10) in (9), the determinant of A is given as ence for the current icf can easily be computed by phase shift-
(
A = −3 v 2pa + v 2pb + v 2pc ) (11)
ing the reference for vp by 90° and from the knowledge of the
value of the filter capacitor Cf.
The reference for the current i1 can be obtained from (12).
In similar manner, the solution of (8) is computed as
However we stipulate that DG-1 supplies a fraction of the real
⎡i1a ⎤ ⎡iLa ⎤ ( ) (
⎡3 Pg + P2 v pa + 3 Qg + Q2 v pb − v pc ⎤
⎢ ⎥
)( ) and reactive power demanded by the load. Then noting that
⎢i ⎥ = ⎢i ⎥ + 1
( )
⎢ 1b ⎥ ⎢ Lb ⎥ A ⎢3 Pg + P2 v pb + 3 Qg + Q2 v pc − v pa ⎥( )( ) Pg + P2 = PL − P1
( ) ( )( )
⎢3 P + P v + 3 Q + Q v − v ⎥ (17)
⎢⎣i1c ⎥⎦ ⎢⎣iLc ⎥⎦ Qg + Q2 = QL − Q1
⎢⎣ g 2 pc g 2 pa pb ⎥
⎦
(12) Also note that Pg + P2 is a dc quantity, while PL will contain a
double frequency component due to the unbalance in the load.
B. Compensator Control It will also be distorted due to the nonlinear nature of the load.
Hence it is evident from (17) that P1 should cancel out both
The equivalent circuit of one phase of the converter is these components. Similarly, Q1 should also cancel out both
shown in Fig. 5. In this, u⋅Vdc1 represents the converter output the double frequency and nonlinear component of QL.
voltage, where u = ± 1. The main aim of the converter control Now let us stipulate that P1 supplies λ1P times the average
is to generate u. From the circuit of Fig. 1, a state space de- load power PLav and Q1 supplies λ1Q times the average load
scription of the system can be given as reactive power QLav. Then we can rewrite (17) as
x& = Ax + Buc (
Pg + P2 = PL − λ1 p × PL = PLav 1 − λ1 p = PL1)
where uc is the continuous time version of switching function (
Qg + Q2 = QL − λ1Q × QL = QLav 1 − λ1Q = QL1 ) (18)
⎡i1a ⎤ ⎡iLa ⎤ (
⎡3PL1v pa + 3QL1 v pb − v pc
⎢
)⎤⎥
⎢ i ⎥ = ⎢i ⎥ + 1
(
⎢ 1b ⎥ ⎢ Lb ⎥ A ⎢3PL1v pb + 3QL1 v pc − v pa )⎥ (19)
⎢⎣i1c ⎥⎦ ⎢⎣iLc ⎥⎦ ⎢3P v + 3Q v − v
⎢⎣ L1 pc L1 pa ( pb )⎥⎥⎦
C. Reference Generation for DG-2
2
where K is a gain matrix and xref is the reference vector. The VP V2 cos δ − V p
gain matrix, in this paper, is obtained through linear quadratic Q2 = (21)
X2
where δ = δ2 − δp. Let us assume that P2 is equal to λ2P times the load, i.e., λ1Q = 1 and λ2Q = 0. Furthermore the real power
the average load power PLav and Q2 is equal to λ2Q times the sharing between the two DGs are set as λ1P = 0.2 and λ2P =
average load reactive power QLav. Then we have 0.3.
The load currents are shown in Fig. 4 (a). In this figure and
P2 = λ2 p × PLav all the subsequent figures showing three-phase quantities, the
(22)
Q2 = λ2Q × QLav phases a, b and c are plotted with solid, dashed and dotted
Equations (20) and (21) can be written as lines respectively. It can be seen that they are both unbalanced
and distorted. The PCC voltages are shown in Fig. 4 (b). They
λ2 p PLav X 2 = VP V2 sin δ attain the steady state within 0.1 s (i.e., 5 cycles). The real and
2
reactive powers are shown in Fig. 5 (a) and (b) respectively. It
λ2Q QLav X 2 + V p = VP V2 cos δ can be seen that Pg and P2 attain constant values in the steady
state in accordance to their power share, i.e., 50% and 30%
From the above equation δ is given by respectively. The compensator supplies 20% power as well as
the nonlinear, unbalanced component of the load. The reactive
⎛ ⎞
−1 ⎜ λ2 p PLav X 2
⎟ powers supplied by the source and DG-2 are zero and hence
δ = tan ⎜ 2 ⎟
(23) the compensator supplies the entire reactive power demand of
⎜ λ2Q QLav X 2 + V p ⎟ the load, as evident from Fig. 5 (b).
⎝ ⎠
Then the magnitude of the back emf is calculated from (21) as
2
λ2Q QLav X 2 + V p
V2 = (24)
VP cos δ
Fig. 7. (a) RMS PCC voltage and (b) instantaneous grid power during island-
ing and resynchronization.