energies

Article
Distributed Generation Islanding Effect on
Distribution Networks and End User Loads Using
the Load Sharing Islanding Method
Maen Z. Kreishan 1 , George P. Fotis 2 , Vasiliki Vita 1 and Lambros Ekonomou 1, *
1 Department of Electrical and Electronic Engineering, City, University of London, London EC1V 0HB, UK;
[email protected] (M.Z.K.); [email protected] (V.V.)
2 Department of Electrical and Electronic Engineering Educators, School of Pedagogical and Technological
Education, N. Heraklion, 141 21 Athens, Greece; [email protected]
* Correspondence: [email protected]; Tel.: +44-20-7040-8111

Academic Editor: Josep M. Guerrero

Received: 28 July 2016; Accepted: 8 November 2016; Published: 16 November 2016

Abstract: In this paper a realistic medium voltage (MV) network with four different distributed
generation technologies (diesel, gas, hydro and wind) along with their excitation and governor
control systems is modelled and simulated. Moreover, an exponential model was used to represent
the loads in the network. The dynamic and steady state behavior of the four distributed generation
technologies was investigated during grid-connected operation and two transition modes to the
islanding situation, planned and unplanned. This study aims to address the feasibility of planned
islanding operation and to investigate the effect of unplanned islanding. The load sharing islanding
method has been used for controlling the distributed generation units during grid-connected and
islanding operation. The simulation results were validated through various case studies and have
shown that properly planned islanding transition could provide support to critical loads at the event
of utility outages. However, a reliable protection scheme would be required to mitigate the adverse
effect of unplanned islanding as all unplanned sub-cases returned severe negative results.

Keywords: distributed generation; islanding; load sharing; distribution networks; Neplan

simulation tool

1. Introduction
The increasing energy demand, the decay in finite traditional energy resources and the growing
concerns from environmental pollution brought considerable challenges on energy and electricity
supply industry [1]. Many measures have been taken so far to overcome these challenges or at
least mitigate their effect on electricity supply industry with the deregulation of electricity markets
across developed countries, the investments in renewable energy deployments, the developments in
electricity distribution networks and the legislation of new energy policies to be at the core of those
measures [2].
The developments in distribution networks had major impact on the diffusion of distributed
generation (DG). The traditional principle of distribution networks where electricity used to travel
long distances from distant centralized power generation plants through radial passive networks to
supply the end users has changed. It became active networks that support more dynamic energy flows
with real-time management. Information technology and power electronics have facilitated the control
of these new active networks in order to withstand the increasing number of DG units connected at
the low and medium voltage level of distribution networks [3].
The size of system support that DG can bring to electricity networks’ operators is huge, as it
contributes in: network loss reduction, improved voltage support, transmission and distribution

Energies 2016, 9, 956; doi:10.3390/en9110956 www.mdpi.com/journal/energies

Energies 2016, 9, 956 2 of 24

congestion relief, deferment in capital investments for transmission and distribution networks
infrastructure, and the overall improved system security, efficiency, and reliability [4].
Despite the promising benefits DG brought for distribution networks much work has been done
to reveal the adverse effect resulting from DG penetration on utility grid. Those operational effects can
be summarized into the following: increased thermal constraints on utility’s assets when DG capacity
exceeds the utility upstream firm capacity, voltage level control issues as DG may interfere with the
standard utility voltage control causing abnormal voltage levels at downstream loads, reversed power
flows since DG have changed the traditional power flow from unidirectional flow to bidirectional flow,
and network protection issues which form a significant challenge against further DG deployment [5–8].
Protection challenges arise from the fact that DG penetration may pose serious issues on the
normal operation of utility protection schemes, such challenges include: false tripping of feeders,
blinding of protection, increased or decreased fault levels, unintentional islanding, unsynchronized
reclosing, and prevention of automatic reclosing [9].
There is a close relationship between islanding and reclosing problems, as unwanted islands may
lead eventually to unsynchronized reclosing. Due to the unpredicted behavior of power systems with
unintentional islanding many technical standards state that any DG unit installed at a distribution
feeder must be automatically tripped should the feeder become disconnected from the network,
this situation is called: “Anti-Islanding Protection” [9,10].
Anti-islanding may reduce DG units’ efficiency as well as network reliability if DG units were
prevented operation despite their supply support benefits [3]. New developments in DG technologies
such as gas turbines, diesel engines, hydro turbines, and wind energy conversion system (WECS)
promise more efficient and greener power production [11]. Many studies have been conducted on
distributed generation in terms of generation technology, generation resources availability, supply
quality and reliability, as well as protection. While previous work on islanding phenomenon in power
systems focused on many different aspects those in particular to this study were related to intentional
(planned) and unintentional islanding, anti-islanding detection, and protection from islanding.
Those studies examined the dynamic and steady state behavior of selected systems during
grid-connected and islanding operation, hence certain conclusions about islanding whether or not it
was a viable option will be made depending on the proposed control strategies as well as the existing
DG capacity. While some studies focused on allowing islanding, others aimed to address the problems
related to islanding protection and tried to suggest new methods for islanding detection [12,13].
The results revealed that although wind turbines accounted for small portion of the installed
DG capacity, most wind turbines needed to be disconnected in order to preserve system frequency
within the permissible range if islanding operation was considered, this is related to the insufficient
reserve and inertia issues of intermittent energy sources [13]. The study also suggested a new
approach to enable successful islanding with the presence of wind generation and results were
promising. This approach was based on demand as frequency controlled reserve (the DFR technology).
This technology provides support to islanding frequency at the presence of intermittent generation by
performing demand side management on many end-users loads. Since it was found that disconnecting
large number of end-users loads such as electric heating, freezers, and refrigerators for short period
of time not affecting those loads performance but it returns considerable benefit in terms of load
reduction to preserve frequency [13].
The unintentional islanding was investigated as a part of an islanding study proposed a power
sharing approach to improve the dynamic behavior of an islanded micro-grid system [14]. The dynamic
behavior of the islanded part was studied by means of PSCAD/EMTDC software package. The control
strategy capability in providing frequency and voltage control to the islanded system was observed.
Due to the major contribution from the primary frequency control units, the gas unit added a smoother
regulation to the island frequency by injecting less real power. Furthermore, voltage regulation was
successful in bringing the bus voltage to the permissible range after unplanned islanding has been
detected [14].
Energies 2016, 9, 956 3 of 24

A new anti-islanding detection method was proposed [15] as part of new requirement for utilities
in order to provide protection from islanding negative impacts such as: out-of-phase reclosing, electric
shocks from supposedly de-energized islands, and power quality degradation. The new proposed
anti-islanding method detects islanding by calculating frequency slip and acceleration between two
systems using time synchronized measurements. This method was referred to as islanding detection
scheme based on wide-area measurement (IDS_WA) [15]. Results showed that islanding detection
scheme based on local area measurement method (IDS_LA) detects islanding successfully if the existing
power mismatch between utility grid and DG connected network is large. Meanwhile, IDS_LA failed
to establish reliable detection within certain period of time when the power mismatch was small.
IDS_WA method managed to detect islanding regardless of the existing power mismatch as simulation
results revealed [15].
The current paper aims to contribute the current studies on the distributed generation islanding
effect on distribution networks and end user loads. The paper studies the feasibility of planned
islanding operation to provide support at the event of utility outages using the load sharing islanding
method for controlling the distributed generation units during grid-connected and islanding operation.
Furthermore the paper investigates the effect an unplanned islanding on a realistic MV network.
The obtained simulation results were validated through various case studies and have shown that
properly planned islanding transition could provide support to critical loads at the event of utility
outages. The presented work contributes the conducted worldwide research on islanding phenomenon
aiming to increase the DG units and the power distribution networks’ efficiency.

2. Islanding Phenomenon in Power Systems

2.1. Islanding Definition

Islanding phenomenon (or power islands) is an electrical condition exists in power distribution
networks containing distributed generation when a zone or an area of distribution network supplied
by both power grid and the DG gets isolated from grid supply due to any reason, thus DG units
continue to energize some or the entire load presented within the isolated section. This implies that
the distribution utility loses control over supply in the island and cannot guarantee supply quality
neither its stability [16,17].
Traditionally, distribution networks used to be passive networks and it did not have any power
sources injecting power into grid, therefore any fault used to occur in the network was dealt with
by utilities’ protection schemes on that basis. However the introduction of DG at distribution level
made this assumption void and further investigation and research were needed to address the issue of
islanding in power systems [18].
In other words, when the amount of power generated by multiple DG units connected at the LV
grid matches the absorbed amount of power by loads connected to the same network; utility protection
units may fail to detect this situation as no power is flowing through utility grid and that may indicate
islanding situation, bearing in mind that DG units may self-disconnect following islanding occurrence
without any consideration to the loads presented in the island [17].

2.2. Unintentional Islanding

Islanding situation may pose threat to the distribution system especially if it occurs accidentally
to part of the system that was not designed to operate in such conditions. Below are some of the most
potential hazards unintentional islanding may bring to a system if it occurs [16,19]:

(1) Power quality degradation: Voltage and frequency levels in an unintentional island may severely
vary and usually will exceed the permissible limits which will put end user loads, DG units,
and utility assets at risk.
(2) Out of phase re-closure: Utility protection equipment are designed usually to auto-reclose after
fault clearance in order to reconnect service, however unintentional islands will cause voltage
Energies 2016, 9, 956 4 of 24

on both sides being out of phase if reclose occurred. This will pose serious risk on rotational
machines’ prime movers due to currents and high mechanical torques produced as result of out
of phase re-closure.
(3) Loss of grounding: Power systems usually have the earthling point at the utility side, thus
unintentional islanding will likely leave the existing DG units and loads within the island to be
improperly grounded.
(4) Safety concerns: Unwanted islanding will pose health risk on utility’s maintenance teams,
since they work on a network being thought de-energized when it is actually still energized.

2.3. Islanding Detection

Concept of islanding detection is based on monitoring the DG units’ output parameters and
sometimes the grid parameters to notice any change in them, and then conclude whether or not
islanding has occurred. On that basis, islanding detection techniques can be categorized into remote
and local techniques, while the latter can be divided into passive, active, and hybrid technique [18].
Islanding detection techniques face three main challenges. The first challenge is related to the
range in power difference between the DG output power and that absorbed by the load in which
islanding condition may not be detected; this is referred to as the non-detection zone (NDZ). The second
challenge comes from the amount of degradation the detection technique may pose upon DG output
power quality and stability; this is related mainly to active techniques that aim to inject certain
disturbances to detect islanding [20]. Finally, the type of loads presented in power islands play
an important role in islanding detection techniques success, since certain types of loads that can be
modelled as parallel RLC circuits form some obstacles against detection methods, while constant
power loads or current harmonics generating nonlinear loads make it easier for detection techniques
to detect islanding [20].
A brief description and some examples for each one of the islanding detection techniques
mentioned previously are as follow:

â Local techniques: the main principal here works by measuring the main electrical parameters
(voltage, current, frequency, etc.) at the DG side [20], local techniques are sub-categorized into:

• Passive techniques: those techniques specifically monitor the electric parameters at DG

side and compare them with predefined values; despite their ease of implementation and
non-harmful effect on DG stability those techniques suffer from wide NDZ. Such techniques
include: Over/under-voltage and over/under-frequency, phase jump detection (PJD),
detection of voltage and current harmonics [18,20].
• Active techniques: active techniques work by injecting certain disturbances at the output of
the DG unit and depending on the amount of change in the electric parameters, islanding
might be detected if that amount was considerable. The NDZ of active techniques is
very small or almost negligible compared to passive techniques; however the amount
of degradation they may produce is significant. Such techniques include: impedance
measurement, active frequency drift, sandia frequency shift [18,20].
• Hybrid techniques: these methods employ a combination of both active and passive
techniques, where active techniques are used when passive techniques suspect islanding
situation. Those techniques include: positive feedback and voltage imbalance, voltage and
reactive power shift [18,20].

â Remote techniques: these techniques use certain communication between utility substation and
DG. Remote techniques detect islanding with high degree of reliability and accuracy, however
their implementation costs are quite expensive compared to the local techniques. These techniques
include: impedance insertion and power line carrier communications [18,20].
Energies 2016, 9, 956 5 of 24

3. Methodology

3.1. Conditions for Islanding and Grid-Connected Operation

For the sake of system control and protection, DG units are not usually set to participate in
regulating the frequency of the grid neither the voltage during a grid-connected mode. They are set
to operate in a constant power mode and that is to provide fixed level of real power at high power
factor (PF), while the grid frequency and voltage levels are determined by the utility control systems.
The governor system attached with those DG units is responsible for providing that constant level of
real power, while the excitation and automatic voltage regulation system ensures the constant reactive
power output through automatic voltage regulation reference [14,21].
However, when there is a transition from the grid-connected mode to the islanding mode whether
it was for intentional or unintentional reasons, these external power control systems are bypassed and
new control systems take over depending on the islanding method agreed upon by the utility and
the DG units’ operators. If this is not the case, then existing governor and excitation control systems
remain operating at a predefined reference values chosen according to international standards for safe
and stable operation whilst in the islanding mode [14,21].
Moreover, operating part of the distribution networks in islanding mode needs to take into
consideration certain limits for frequency, voltage and harmonics. These limits must be derived from
international standards such as those listed in Table 1, and that ensure a safe, reliable, and standard
quality supply to all end users loads existing within the island [22].

Table 1. Standard Power Systems Voltage and Frequency Limitations.

Voltage Window Frequency Window Reference Standard

0.9 to 1.1 UN 0.99 to 1.01 fN IEEE P1547
0.8 to 1.15 UN 0.94 to 1.02 fN VDE 0126, OVE E2750 (partly)
0.8 to 1.2 UN 0.994 to 1.006 fN CEI 0-21, CEI 0-16

According to the IEEE 1547-2003 standard, distributed generation units are required to cease
energizing the part of the network they supply during certain clearing times should the voltage across
the system falls within the correspondent voltage ranges illustrated in Table 2. Clearing time is defined
as the time taken by the DG unit to de-energize its supplied part of the network completely after
an abnormal condition has occurred [23,24].

Table 2. IEEE Abnormal Voltage Range Clearing Times.

Voltage Range (Percentage of Nominal Voltage) Clearing Time in Seconds

V < 50% 0.16
50% ≤ V < 88% 2
88% ≤ V < 110% Normal Operation
110% ≤ V < 120% 1
V ≥ 120% 0.16

As for the system frequency, according to the same IEEE standard and IEEE 1547.2-2008, DG units
are required to cease energizing their supplied part of the network within a pre-defined clearing times
should the system frequency lies within the specified frequency ranges of Table 3. Depending on the
DG unit capacity the frequency limits and their correspondent clearing times may be fixed or field
adjustable. For instance, DG units with capacity over 30 kW those limits need to be field adjustable,
while for capacity less than or equal to 30 kW those limits can be either field adjustable or fixed [23].
If those abnormal operating situations listed above were avoided via properly implemented and
secured protection scheme, then there is no reason why islanding should not be utilized. Since it helps
Energies 2016, 9, 956 6 of 24

increasing the reliability of the distribution networks besides the many other benefits of the distributed
generation mentioned earlier. In order to make the islanding mode a feasible option, two key states
need to be addressed, studied, and planned fully and completely. Those states are: (1) The transition
from grid-connected into islanding mode; and (2) island operation after disconnection from utility
control [25].

Table 3. IEEE Abnormal Frequency Range Clearing Times.

DG Capacity Voltage Range (Percentage of Nominal Voltage) Clearing Time in Seconds

f > 50.4 0.16
DG ≤ 30 kW 49.4 ≤ f ≤ 50.4 Normal Operation
f < 49.4 0.16
f > 50.4 0.16
49.8 ≤ f ≤50.4 Normal Operation
DG > 30 kW
47.5 ≤ f < 49.8 Adjustable 0.16–300
f < 47.5 0.16

3.2. Load Sharing Islanding Method

Under this islanding strategy all DG units within an island are compelled to participate in the
islanding transients as well as to serve the existing load by sharing it in an active way [21,26,27].
This strategy can be achieved by forcing all DG units’ control systems following the islanding
occurrence to switch from fixed power control mode (P-Q) to fixed frequency and voltage regulation
control mode (V-f), these changes are based upon voltage following and frequency adjustments using
pre-defined voltage and frequency characteristics. This implies that as the island frequency increases
then the total real power output from all of the DG units decrease or vice versa, while the change in
the voltage level at any DG unit’s terminal causes the change in its reactive power control in order to
regulate the voltage across the islanded network.

4. Distribution System Description

The test network used for this study, shown in Figure 1, is part of a realistic MV network
with a nominal voltage of 20 kV [28]; it consists of the utility main feeder NETZ rated at 110 kV,
with (2000–1000) MVA for max–min initial symmetrical short circuit apparent power. The HV/MV part
of the network consists of two transformers each rated at 110/20 kV, and 30 MVA rated apparent power.
The MV network consists of five feeders, while the investigated part of that network Feeder#2
is connected to the MV network at the point of common coupling (PCC); this feeder consists of four
generators which is the scope of interest for the case studies. The remaining two generators are
connected to Feeders#3 and 4 respectively. Technical data of the six distributed generators can be
found in Table 4.
Bearing in mind that all generators were ungrounded in order to prevent stator damage caused
by fault currents from utility grid, moreover high resistive grounding or ungrounded systems are used
where allowed when the continuity of supply to loads during single line earth fault is desired. All of
the feeders’ lines are made of copper cable lines with a cross-sectional dimension of 3 × (1 × 95) mm2 ,
the corresponding lengths for all the lines in the MV network are illustrated in Table 5 [28].
The remaining part of the MV network, the loads, is as follows: Loads L1-1 to L1-5, L3-1 to L3-3,
L4-1 to L4-6, and L5-1 to L5-5 are all with 0.8 MVA apparent powers with 0.95 lagging PF [28]. As for
loads on the studied Feeder#2, they are as follows: L2-1, L2-2, L2-5, L2-8, L2-9, L2-11, and L2-12 are
rated at 0.4 MVA, while loads L2-3, L2-4, L2-6, L2-7, and L2-10 are rated at 0.8 MVA; all loads are with
0.95 lagging PF [28].
Energies 2016, 9, 956 7 of 24
Energies 2016, 9, 956 7 of 24

Figure
Figure 1.
1. Studied
Studied network
network single
single line
line diagram.
diagram.

Table 4. DG
Table 4. DG Units
Units Characteristics.
Characteristics.
Rated Nominal Rated Apparent
Generator
Generator Name Name
Rated Nominal Voltage Neutral Point
Voltage (kV) (kV) Rated
PowerApparent
(MVA) Power (MVA) Neutral Point
GEN1 GEN1 20 20 3 3 Ungrounded
Ungrounded
GEN2 GEN2 20 20 3 3 Ungrounded
Ungrounded
GEN3 GEN3 20 20 3 3 Ungrounded
Ungrounded
GEN4 GEN4 20 20 3 3 Ungrounded
Ungrounded
GEN5 GEN5 20 20 3 3 Ungrounded
Ungrounded
GEN6 20 3 Ungrounded
GEN6 20 3 Ungrounded

Table 5.
Table Test Network
5. Test Network Lines
Lines Lengths
Lengths in
in km.
km.

Line
Line Name Length
Name Length(km)
(km) Line Name
Line Name Length
Length(km)
(km) Line Name)
Line Name Length
Length(km)
(km)
L1-1 1.00 L2-7 0.50 L4-3 1.00
L1-1 1.00 L2-7 0.50 L4-3 1.00
L1-2 0.50 L2-8 0.70 L4-4 0.80
L1-2 0.50 L2-8 0.70 L4-4 0.80
L1-3
L1-3 0.80
0.80 L2-9
L2-9 0.70
0.70 L4-5
L4-5 0.50
0.50
L1-4
L1-4 0.80
0.80 L2-10
L2-10 0.70
0.70 L4-6
L4-6 0.50
0.50
L1-5
L1-5 1.00
1.00 L2-11
L2-11 1.00
1.00 L5-1
L5-1 0.50
0.50
L2-1
L2-1 0.50
0.50 L3-1 1.00
1.00 L5-2
L5-2 0.80
0.80
L2-2
L2-2 0.50
0.50 L3-2 1.00
1.00 L5-3
L5-3 0.80
0.80
L2-3
L2-3 0.80
0.80 L3-3 1.00
1.00 L5-4
L5-4 1.00
1.00
L2-4
L2-4 1.00
1.00 L3-4
L3-4 1.00
1.00 L5-5
L5-5 0.70
0.70
L2-5
L2-5 1.00
1.00 L4-1
L4-1 1.50
1.50
L2-6 1.00 L4-2 0.80
L2-6 1.00 L4-2 0.80
Energies 2016, 9, 956 8 of 24

5. Dynamic Models of the Network Components

5.1. The Load Model

The load model selected for the purpose of the dynamic simulation of this study is a static
exponential load model. Naturally, loads consist of static and dynamic components which are derived
from the nature of the load whether it was residential, commercial, or industrial. Static load models
describe the dependency of load active and reactive power as a function of the system voltage and/or
frequency at one specific instant of time. While on the other hand dynamic load models describe this
relation at any instant of time, hence static load models are not time dependent [29].
However, static load models have been sufficient enough to be used for a very long time to
study load behavior for both the static and dynamic components [29], therefore and for the sake of
simplicity the static exponential load model has been used throughout this study to describe the static
components for all the loads presented in the island and to approximate their dynamic components
based on the facts explained earlier.
An equation that can describe the load power dependency as a function of system voltage and
frequency under the static exponential load model is as follows [30,31]:

U nP ∆f
P = P0 ( ) (1 + Fp ) (1)
U0 f0

U nQ ∆f
Q = Q0 ( ) (1 + Fq ) (2)
U0 f0
where:
P: is the ongoing active power of the load
Q: is the ongoing reactive power of the load
P0 : is the nominal active power of the load
Q0 : is the nominal reactive power of the load
U: is the magnitude of ongoing node voltage
U0 : is the nominal system voltage
f0 : is the rated system frequency
∆f: is frequency difference from rated frequency
nP: is the dependence of load active power on voltage
nQ: is the dependence of load reactive power on voltage
Fp : is the dependence of load reactive power on frequency
Fq : is the dependence of load reactive power on frequency
Typical values for the loads used for all four case studies dynamic simulation were [29]: nP = 1,
nQ = 2, = 0 and = −1.

5.2. DG Units
The diesel, gas and hydraulic DG units are modeled as a single-mass synchronous machine.
The machine electrical system is represented in the d-q-0 frame with two rotor windings on each
axis. The excitation and governor systems of the machine are also included in the simulation model.
The parameters for the diesel, gas and hydraulic units dynamic model used in this study are shown in
Table 6.
The wind DG unit is modeled as a direct drive synchronous generator. Rotor and generator shafts
are mounted to the same shaft without gear-box. The generator is a permanent magnet synchronous
generator designed for low speed. The permanent magnet machine dynamic model used for this
study is the one embedded in Neplan simulation tool and has the parameters of Table 6. Permanent
magnet synchronous generator (PMSG) output is injected to the utility grid through back-to-back
Energies 2016, 9, 956 9 of 24

Energies 2016, 9, 956 9 of 24

pulse width modulator (PWM) converter. The converter system ensures that both generator and grid
currents are sinusoidal and it can enable for variable speed control. To achieve regulation of the current
for both d-axis and q-axis of the current components in a very fast inner control loop. Meanwhile, the
components for the grid side PWM, the regulator used ensures regulation for both d-axis and q-axis of
slower outer control loop that regulates the reactive and active powers is used to define the current
the current components in a very fast inner control loop. Meanwhile, the slower outer control loop
references [32].
that regulates the reactive and active powers is used to define the current references [32].
Table 6. Synchronous Machine Dynamic Model Parameters for the four different types of DG units.
Table 6. Synchronous Machine Dynamic Model Parameters for the four different types of DG units.
Synchronous Synchronous Stator Mechanical
Transient Data Sub-Transient Data
Machine Reactances
Synchronous Resistance
Stator Data
Mechanical
Synchronous Transient Data Sub-Transient Data
Dynamic
Machine Reactances Resistance Data
H (MW*s/
Model Model
Dynamic 0 0 0 T0 Xd00 Xq00 00 00
Tq0 R
(pu) Xd (pu) Xq (pu)Xi (pu)Xd (pu)Xq (sec) Td0
(sec)q0 (pu) (pu) T
(s)
d0 (s) R
H (MW*s/
MVA)
Parameters
Parameters (pu) (pu) (pu) (pu) (pu) (sec) (sec) (pu) (pu) (s) (s) MVA)
Diesel DGDG
Diesel unitunit 1.6 1.6 0.950.950.1050.105 0.330.33 0.70 0.70 4.04.0 1.2 1.2 0.24
0.24 0.30
0.30 0.05
0.05 0.05
0.05 -- 1.761.76
Gas Gas
DGDG unitunit 1.5 1.5 0.750.75 0.1 0.1 0.256 0.256 0 0 0.53 0.53 1.20
1.20 0.168
0.168 0.184
0.184 0.03
0.03 0.03
0.03 1 1
Hydro
Hydro DGDGunitunit1.2361.2360.750.750.155
0.1550345 0345 0.70 0.70 4.17
4.17 1.20
1.20 0.264 0.211
0.264 0.211 0.03
0.03 0.19
0.19 3.123.12
Wind DG unit 1.05 0.75 0 - - - - - - - - 0.042 5
Wind DG unit 1.05 0.75 0 - - - - - - - - 0.042 5

5.3. Voltage and

5.3. Voltage and Power
Power Control
Control
The
The automated
automated voltage
voltage regulator (AVR) and
regulator (AVR) and excitation
excitation system
system used
used for
for diesel,
diesel, gas
gas and
and hydro
hydro
units is modeled using the block diagram shown in Figure 2 [33]. For simplicity this model is selected
units is modeled using the block diagram shown in Figure 2 [33]. For simplicity this model is selected
from
from the
the pre-defined
pre-defined exciter
exciter models
models in
in Neplan
Neplan under the name
under the name Simple
Simple Exciter
Exciter Model
Model with
with Limits.
Limits.
Moreover,
Moreover, the
the excitation
excitation and
and AVR
AVRmodel
modelparameters
parametersare
areprovided
providedininTable
Table7.7.

Figure
Figure 2.
2. AVR and excitation
AVR and excitation model
model for
for diesel,
diesel, gas
gas and
and hydro
hydro generators.
generators.

Table 7. AVR and excitation model parameters for diesel, gas and hydro generators.
Table 7. AVR and excitation model parameters for diesel, gas and hydro generators.
Parameter Description Value
Parameter
TR (s) Description
Voltage Transducer time constant Value
0.01
KA (pu)
TR (s) Voltage
Voltage regulatortime
Transducer constant
constant 200
0.01
KATA(pu)
(s) Voltage regulator
Voltage time
regulator constant
constant 0.015
200
TA
TB(s)
(s) VoltageLag
regulator time constant
time constant 0.015
10
TB
TC(s)
(s) Lag time
Lead time constant
constant 110
TC (s) (pu)
EFDMAX Lead time
Maximum constant
field voltage 51
EFDMAX
EFDMIN(pu) (pu) Maximum field voltage
Minimum field voltage −55
EFDMIN (pu) Minimum field voltage −5

Figure 3 shows a simplified schematic of the proposed wind energy conversion system (WECS)
Figure
that has been3 shows
used fora simplified
the purpose schematic
of wind DGof the proposed
units modelingwind energy
[34]. conversion
Permanent magnet system (WECS)
synchronous
that has been
generator used output
(PMSG) for the purpose of to
is injected wind
the DG units
utility modeling
grid through[34]. Permanent
back-to-back magnet
pulse widthsynchronous
modulator
generator (PMSG) The
(PWM) converter. output is injected
converter to the
system utilitythat
ensures gridboth
through back-to-back
generator and gridpulse width
currents aremodulator
sinusoidal
(PWM) converter.
and it can enable forThe converter
variable system
speed ensures
control. that both system
This converter generator and topology
circuit grid currents are sinusoidal
as shown in Figure
and it canofenable
4 consists for variable
two PWM speed
converters control.together
connected This converter systemDC
via a common circuit
link. topology
Informationas shown in
about the
Figure 4 consists
principle of two
and design PWM converters
of back-to back PWM connected
used in WECS together viafound
can be a common DC link. Information
in [34,35].
about the principle and design of back-to back PWM used in WECS can be found in [34,35].
Energies 2016, 9, 956 10 of 24
Energies 2016,
Energies 2016, 9,
9, 956
956 10 of
10 of 24
24
Energies 2016, 9, 956 10 of 24

Figure 3.
Figure 3. Direct
Direct drive permanent
Direct drive permanent magnet
magnet WECS.
WECS.
Figure 3. Direct drive permanent magnet WECS.

Figure 4.
Figure 4. Back-to-back
Back-to-back PWM
PWM converter.
converter.
Figure 4.
Figure Back-to-back PWM
4. Back-to-back PWM converter.
converter.

The converter
The converter to to thethe PMSG
PMSG side side acts
acts as as aa PWM
PWM rectifier
rectifier while
while the
the converter
converter connected
connected to to the
the grid
grid
side The as
acts converter
PWM to the PMSG
inverter. PWM side actsconverts
rectifier as a PWM the rectifier
rectifier
terminal while the of
voltage converter
PMSG connected
into DC to thewhile
voltage, grid
side acts as PWM inverter. PWM rectifier converts the terminal voltage of PMSG into DC voltage, while
sideLC
an acts
acts as
asPWM
filter usedinverter.
isPWM inverter.
at the DC DC PWMPWM
link in rectifier
order converts
rectifier
in order converts
to achieve
achieve the terminal
the
smooth DCvoltage
terminal ofThe
voltage
output. PMSG
of
two into
PMSG DC
intovoltage,
converters DC while
voltage,
are isolated
isolated
an LC filter is used at the link to smooth DC output. The two converters are
an LC
while
by filter
an LCisfilter
used isatused
the DC atthe link
the DCinlink
order
link intoorder
achievetoas smooth
achieve DC output.
smooth The[36].
DC output.twoTo converters
The are
two converters isolated
are
by aa capacitor
capacitor connected
connected at
at the DC DC which
link which acts
acts as aa stiff
stiff voltage
voltage source
source [36]. To achieve
achieve aa settlement
settlement
by a capacitor
isolated
between bothby a
both the connected
capacitor
the generator at
generator and the
connected DC
and the link
at the which
the converterDC acts
link
converter ratings as
whicha
ratings as stiff
actsvoltage
as well as
well as a source
stiff
as to [36].
voltage
to minify To achieve
source
minify saturation, [36].
saturation, phase aTosettlement
phase angleachieve
angle ofof
between
between
a settlementboth
generator current the
betweengenerator
current needs
needs to both
to be and
the
be adjusted the
adjusted to converter
generatorto ensureand
ensure thatratings
the as
converter
that stator well
stator terminal as to
ratings minify
terminal voltage as well
voltage and saturation,
as to
and internal phase
minify
internal induced angle
saturation,
induced voltage of
voltage
generator
generator
phase angle
are equal
equal in current
of
in magnitudeneeds
generator
magnitude [34]. to be
current
[34]. To adjusted
needs
To achieve totoensure
be
achieve regulation
regulation ofthat
adjusted stator
to
of the terminal
ensure
the current that voltage
stator
current components and internal
terminal
components for for the induced
voltage
the grid
grid side and
side PWM, voltage
internal
PWM, the the
are
are equalvoltage
induced
regulator in Figure
of magnitude
are5 equal
is [34].inTo
used. It achieveregulation
magnitude
ensures regulation
[34]. To for of both
the current
achieve components
regulation
d-axis and of theof
q-axis for
current
the the grid components
side PWM,
components
current for the
in
regulator of Figure 5 is used. It ensures regulation for both d-axis and q-axis of the current components in
regulator
grid side of
PWM,Figurethe 5 is used.
regulator It ensures
of Figure regulation
5 is used. for
It both
ensures d-axis and
regulation q-axis
for of
boththe current
d-axis andcomponents
q-axis of in
the
aa very
very fast
fast inner
inner control
control loop.loop. Meanwhile,
Meanwhile, the the slower
slower outer
outer control
control loop
loop that
that regulates
regulates thethe reactive
reactive and
and
a very
current fast
active powers inner
components
powers is control
is used
used to in aloop.
very
to define
define the Meanwhile,
fast inner
the current the
control
current references slower
loop.
references [32]. outer control
Meanwhile,
[32]. loop
the that
slower regulates
outer the
controlreactive
loop and
that
active
active powers
regulates is used to
the reactive anddefine
active thepowers
currentis references [32]. the current references [32].
used to define

Figure 5.
Figure 5. Grid
Grid side
side PWM
PWM regulator.
regulator.
Figure
Figure 5.
5. Grid
Grid side
side PWM
PWM regulator.
regulator.
The generator
The generator side
side PWM
PWM regulation
regulation isis achieved
achieved through
through the
the regulator
regulator ofof Figure
Figure 6.
6. This
This regulator
regulator is
is
The
fittedThe generator
withgenerator side
fast current
current PWM regulation
controllers as the is
the grid
grid achieved
side through
PWMthrough the
regulatortheregulator
to ensure
ensure theof
theof Figure 6.
regulation This
of regulator
both DC and
andis
fitted with fast side PWM regulation
controllers as isside
achieved
PWM regulator to regulator Figure 6.of
regulation This
bothregulator
DC
fitted
AC with fast
voltages; current
however controllers as
DCcontrollers
and AC the
AC voltage
voltagegrid side PWM
regulators regulator
are used
used to ensure
to define
define the regulation
the current
current of both
referencesof DC
[32]. and
is fitted
AC with however
voltages; fast current
DC and as theregulators
grid side PWM
are regulator
to to ensure
the the regulation
references both DC
[32].
AC voltages; however DC and AC voltage regulators are used to define the current references [32].
and AC voltages; however DC and AC voltage regulators are used to define the current references [32].
Energies 2016, 9, 956 11 of 24
Energies 2016, 9, 956 11 of 24

Figure
Figure 6.
6. Generator
Generator side
side PWM
PWM regulator.
regulator.

The parameters used with the WECS components: wind turbine model, grid side PWM
The parameters used with the WECS components: wind turbine model, grid side PWM regulator,
regulator, and generator side PWM regulator are provided in Table 8.
and generator side PWM regulator are provided in Table 8.
Table 8. WECS components’ model parameters.
Table 8. WECS components’ model parameters.
Wind Turbine Grid Side PWM Regulator Generator Side PWM Regulator
Parameter
Wind Turbine Value Parameter
Grid Value
Side PWM Regulator Parameter Value
Generator Side PWM Regulator
ParameterTSP Value 0.05 TMI
Parameter 0.0001
Value TMI
Parameter 0.0001
Value
KPP 3.00 TMP 0.001 TMV 0.0001
TSP 0.05 TMI 0.0001 TMI 0.0001
KIP 0.60 KPQ 1 KPVT 0.1
KPP 3.00 TMP 0.001 TMV 0.0001
PMIN 0.10 KIQ 10 KIVT 1
KIP 0.60 KPQ 1 KPVT 0.1
PMIN PMAX 0.10 1.12 KPID
KIQ 0.02
10 KPID
KIVT 10 1
PMAX RPMIN 1.12−0.45 KPIDKIID 100
0.02 KIID
KPID 10 10
RPMINRPMAX −0.45 0.45 MIN_IDREF
KIID −1.0
100 MIN_IDREF
KIID −1.010
RPMAXRTBETA 0.4510.0 MIN_IDREF
MAX_IDREF −
1.01.0 MAX_IDREF
MIN_IDREF 1.0−1.0
RTBETA KPPT 10.0150.0 MAX_IDREF
MIN_IQREF 1.0
−1.0 MAX_IDREF
MIN_MQ −1.01.0
KPPTKIPT 150.025.0 MIN_IQREF
MAX_IQREF −
1.01.0 MIN_MQ
MAX_MQ 1.0−1.0
KIPTKPPC 25.03.00 MAX_IQREF KPP 1.0
1 MAX_MQ
KPVDC 5 1.0
KPPCKIPC 3.0030.00 KPPKIP 101 KPVDC
KIVDC 1 5
KIPC 30.00 KIP 10 KIVDC 1
TP 0.30 KPIQ 0.015 KPIQ 1.0
TP 0.30 KPIQ 0.015 KPIQ 1.0
BETAMIN
BETAMIN 0.000.00 KIIQ
KIIQ 100
100 KIIQ
KIIQ 10 10
BETAMAX
BETAMAX 27.0027.00 MIN_MD
MIN_MD −1.0
−1.0 MIN_IQREF
MIN_IQREF −1.0−1.0
BETA0BETA0 10.510.5 MAX_MD
MAX_MD 1.0
1.0 MAX_IQREF
MAX_IQREF 1.0 1.0
NOM_POWER 2.0 2.0
NOM_POWER MIN_MQ
MIN_MQ −1.0
−1.0 MIN_MD
MIN_MD −1.0−1.0
NOM_TURBSPEED 18 18
NOM_TURBSPEED MAX_MQ
MAX_MQ 1.0
1.0 MAX_MD
MAX_MD 1.0 1.0
BLADELENGTH
BLADELENGTH 50 50 MMAX
MMAX 11 MMAX
MMAX 1 1
AIRDENS
AIRDENS 1.2251.225 KPPLL
KPPLL 30
30 KPPLL
KPPLL 30 30
KIPLL
KIPLL 0.1
0.1 KIPLL
KIPLL 0.1 0.1
PLLMIN −0.1 PLLMIN −0.1
PLLMIN −0.1 PLLMIN −0.1
PLLMAX 0.1 PLLMAX 0.1
PLLMAX 0.1 PLLMAX 0.1

6. CaseStudies
6. Case Studies
The
The distribution network was
distribution network was built
built using
using Neplan
Neplan Desktop
Desktop power
power systems
systems simulation
simulation tool
tool and
and
according
according toto each
each case
case study
study the
the correspondent
correspondent distributed
distributed generation technology is
generation technology is injected to the
injected to the
network on Feeder#2 to enable the study of the dynamic behavior of each DG technology
network on Feeder#2 to enable the study of the dynamic behavior of each DG technology following following
an
an islanding
islanding occurrence
occurrence underunder the
the same
same operational
operational values
values for
for the
the rest
rest of
of the
the network
network components.
components.
The current IEEE
The current IEEE standards
standards prevent
prevent the
the operation
operation of
of distribution
distribution system
system in in islanded
islanded conditions,
conditions,
where
where the DG unit is still supplying part of the total load within an island at the event of outage.
the DG unit is still supplying part of the total load within an island at the event of utility utility
This is related
outage. This is to system
related to control, protection,
system control, and personnel
protection, safety issues
and personnel safety[25].
issuesThe[25].
current control
The current
control scheme in the network normally trips the DG units within 2 s as indicated earlier in the
standard voltage and frequency clearing times (Tables 2 and 3). However it is assumed that those
Energies 2016, 9, 956 12 of 24

scheme in the network normally trips the DG units within 2 s as indicated earlier in the standard
voltage and frequency clearing times (Tables 2 and 3). However it is assumed that those clearing times
were adjusted and relaxed to enable the planned and unplanned islanding operation to run smoothly
without interruption by the utility’s protection units during the tests for all of the four case studies.
The total rated system load distributed throughout Feeder#2 for the following four case studies
is set to be (6.46 + j2.125) MVA, however depending on the dynamic response of the load this rated
value may undergo slight changes according to the type of units presented in the island as well as the
islanding method used.
For the purpose of this study, it is assumed that the DG units are considered for autonomous
operation at the planning stage. This implies that all the used DG units in this study have the capability
by working together for each case study within the range of their rated power to supply the total
load in the islanded part of the network without affecting the island desirable operational voltage
and frequency limits. This level of power redundancy ensures that no loads in the islanded part of
the network have to be disconnected during the islanding mode to re-establish the generation with
demand balance in real time [36,37].
The dynamic response of the system is studied in the four case studies by observing the changes in
frequency, voltage, real power, and reactive power across the islanded section following the transition
from utility connected mode to islanded mode using the load sharing islanding method.
Two transition conditions were used throughout the first three case studies using the disturbances
of the dynamic transient stability module inside the used power simulation tool Neplan, these are:
CB opening (in Neplan by opening a logical switch on the upper or lower side of line LIN2-1) and
that demonstrated the planned islanding transition condition, while the other transition condition,
the unplanned islanding, was achieved by triggering a three phase short circuit fault at line (LIN2-1) at
a certain distance on the line (again by using the disturbances of dynamic transient stability module
in Neplan) and that caused the line to be tripped from both sides after some deci-seconds due to the
assumed protection scheme used by the local utility. As for the forth case study, only one transition
condition was used in order to less complicate the case, hence only the planned islanding transition
was used.

6.1. Case Study A: Islanding of Diesel Only Units

In this case study all four DG units (GEN1, GEN2, GEN3, GEN4) are diesel based units connected
to buses N2_4, N2_7, N2_9, and N2_11 respectively as illustrated in Figure 1. Tests were carried out
and results were obtained according to the following sub-case scenarios:

6.1.1. Over-Generation Mode (11.45% Power Mismatch) (Case A1)

This case was with 60% of their rated capacity at almost 0.98 PF (i.e., (1.8 + j0.4) MVA), this situation
made the islanded section of the network to operate in the over-generation level with (7.2/6.46 MW)
Generation to load ratio with around 11.455% real power mismatch, that residual power was exported
to the utility grid during grid-connected mode. There are two sub-case scenarios, the planned islanding
of the diesel units (case A1.1) and the unplanned islanding of the diesel units (case A1.2), whose results
are shown in Figure 7a,b, respectively.
Energies 2016, 9, 956 13 of 24
Energies 2016, 9, 956 13 of 24

(a)

(b)
Figure 7.
Figure 7. Simulation
Simulation results
results for
for islanding
islanding of
of diesel
diesel only
only units
units in
in over-generation
over-generation mode:
mode: (a)
(a) Planned
Planned
(case A1.1); and (b) Unplanned (case A1.2).
(case A1.1); and (b) Unplanned (case A1.2).

6.1.2. Under-Generation Mode (13.31% Power Mismatch) (Case A2)

6.1.2. Under-Generation Mode (13.31% Power Mismatch) (Case A2)
In this loading case the diesel units were with 46.66% of their rated capacity at approximately
In this loading case the diesel units were with 46.66% of their rated capacity at approximately
0.96 PF (i.e., (1.4 + j0.4) MVA), this loading level ensured that the islanded section operated in the
0.96 PF (i.e., (1.4 + j0.4) MVA), this loading level ensured that the islanded section operated in the
under-generation level with (5.6/6.46 MW) generation to load ratio. This situation enforced Feeder#2
under-generation level with (5.6/6.46 MW) generation to load ratio. This situation enforced Feeder#2 to
to import the remaining gap in real power from the utility grid whilst still operating in the grid-
import the remaining gap in real power from the utility grid whilst still operating in the grid-connected
connected mode since there was a 13.31% real power mismatch between generation and load. There
mode since there was a 13.31% real power mismatch between generation and load. There were two
were two sub-case scenarios, the planned islanding of the diesel units (case A2.1) and the unplanned
sub-case scenarios, the planned islanding of the diesel units (case A2.1) and the unplanned islanding
islanding of the diesel units (case A2.2), whose results are shown in Figure 8a,b, respectively.
of the diesel units (case A2.2), whose results are shown in Figure 8a,b, respectively.
Energies 2016, 9, 956 14 of 24
Energies 2016, 9, 956 14 of 24

(a)

(b)
Figure 8.
Figure 8. Simulation
Simulation results
results for
for islanding
islanding of
of diesel
diesel only
only units
units in
in under-generation
under-generation mode:
mode: (a)
(a) Planned
Planned
(case A2.1); and (b) Unplanned (case A2.2).
(case A2.1); and (b) Unplanned (case A2.2).

6.2. Case Study B: Islanding of Gas Only Units

6.2. Case Study B: Islanding of Gas Only Units
In this case study all four DG units (GEN1, GEN2, GEN3, GEN4) are gas type units connected
In this case study all four DG units (GEN1, GEN2, GEN3, GEN4) are gas type units connected to
to buses N2_4, N2_7, N2_9, and N2_11 respectively as illustrated in Figure 1. Tests were carried out
buses N2_4, N2_7, N2_9, and N2_11 respectively as illustrated in Figure 1. Tests were carried out and
and results were obtained according to the following sub-case scenarios.
results were obtained according to the following sub-case scenarios.
6.2.1. Over-Generation
6.2.1. Mode (11.45%
Over-Generation Mode (11.45% Power
Power Mismatch)
Mismatch) (Case
(Case B1)
B1)
Similar to
Similar to the
thecase
caseA1A1there
therewas a generation
was to load
a generation ratio ratio
to load with around 11.455%
with around real power
11.455% real
mismatch. There were two sub-case scenarios, the planned islanding of the diesel
power mismatch. There were two sub-case scenarios, the planned islanding of the diesel units units (case B1.1)
and the
(case unplanned
B1.1) islanding ofislanding
and the unplanned the dieselofunits (case B1.2),
the diesel whoseB1.2),
units (case results are shown
whose resultsin
areFigure
shown 9a,b
in
respectively.
Figure 9a,b respectively.
Energies 2016, 9, 956 15 of 24
Energies 2016, 9, 956 15 of 24

(a)

(b)
Figure 9. Simulation results for islanding of gas only units in over-generation mode: (a) Planned (case
Figure 9. Simulation results for islanding of gas only units in over-generation mode: (a) Planned
B1.1);B1.1);
(case and (b) Unplanned
and (case (case
(b) Unplanned B1.2).B1.2).

6.2.2. Under-Generation Mode (13.31% Power Mismatch) (Case B2)

6.2.2. Under-Generation Mode (13.31% Power Mismatch) (Case B2)
Similar to the A2 case there was a 13.31% real power mismatch between generation and load.
Similar to the A2 case there was a 13.31% real power mismatch between generation and load. There
There were two sub-case scenarios, the planned islanding of the diesel units (case B2.1) and the
were two sub-case scenarios, the planned islanding of the diesel units (case B2.1) and the unplanned
unplanned islanding of the diesel units (case B2.2), whose results are shown in Figure 10a,b,
islanding of the diesel units (case B2.2), whose results are shown in Figure 10a,b, respectively.
respectively.
Energies 2016, 9, 956 16 of 24
Energies 2016, 9, 956 16 of 24

(a)

(b)
Figure 10. Simulation
Figure 10. Simulation results for islanding
results for islanding of
of gas
gas only
only units
units in
in under-generation
under-generation mode:
mode: (a)
(a) Planned
Planned
(case B2.1); and (b) Unplanned (case B2.2).
(case B2.1); and (b) Unplanned (case B2.2).

6.3. Case Study C: Islanding of Hydro Only Units

6.3. Case Study C: Islanding of Hydro Only Units
In this case study all four DG units (GEN1, GEN2, GEN3, GEN4) are hydroelectric units
In this case study all four DG units (GEN1, GEN2, GEN3, GEN4) are hydroelectric units connected
connected to buses N2_4, N2_7, N2_9, and N2_11 respectively as illustrated in Figure 1. Tests were
to buses N2_4, N2_7, N2_9, and N2_11 respectively as illustrated in Figure 1. Tests were carried out
carried out and results were obtained according to the following sub-case scenarios:
and results were obtained according to the following sub-case scenarios:

6.3.1.
6.3.1. Over-Generation
Over-Generation Mode (11.45% Power
Mode (11.45% Power Mismatch)
Mismatch) (Case
(Case C1)
C1)
Similar
Similar to
to the
thecase
caseA1A1there
therewaswas
a generation
a generationto load ratio ratio
to load with around 11.455%
with around real power
11.455% real
mismatch.
power There were
mismatch. twowere
There sub-case scenarios,scenarios,
two sub-case the planned theislanding
planned of the diesel
islanding of units (case C1.1)
the diesel units
and the
(case unplanned
C1.1) and the islanding
unplanned ofislanding
the dieselof
units
the (case
dieselC1.2),
unitswhose resultswhose
(case C1.2), are shown in are
results Figure 11a,b,
shown in
respectively.
Figure 11a,b, respectively.
Energies 2016, 9, 956 17 of 24
Energies 2016, 9, 956 17 of 24

(a)

(b)
Figure
Figure 11.
11. Simulation
Simulation results
results for
for islanding
islanding of
of hydro
hydro only
only units
units in
in over-generation
over-generation mode:
mode: (a)
(a) Planned
Planned
(case C1.1); and (b) Unplanned (case C1.2).
(case C1.1); and (b) Unplanned (case C1.2).

6.3.2. Under-Generation Mode (13.31% Power Mismatch) (Case C2)

6.3.2. Under-Generation Mode (13.31% Power Mismatch) (Case C2)
Similar to the A2 case there was a 13.31% real power mismatch between generation and load.
Similar to the A2 case there was a 13.31% real power mismatch between generation and load. There
There were two sub-case scenarios, the planned islanding of the diesel units (case C2.1) and the
were two sub-case scenarios, the planned islanding of the diesel units (case C2.1) and the unplanned
unplanned islanding of the diesel units (case C2.2), whose results are shown in Figure 12a,b,
islanding of the diesel units (case C2.2), whose results are shown in Figure 12a,b, respectively.
respectively.
Energies 2016, 9, 956 18 of 24
Energies 2016, 9, 956 18 of 24

(a)

(b)
Figure 12.
Figure 12. Simulation
Simulation results
results for
for islanding
islanding ofof hydro
hydro only
only units
units in
in over-generation
over-generation mode:
mode: (a)
(a) Planned
Planned
(case C1.1); and (b) Unplanned (case C1.2).
(case C1.1); and (b) Unplanned (case C1.2).

6.4. Case Study D: Islanding of Mixed Units

6.4. Case Study D: Islanding of Mixed Units
In this case study two WECS were modelled via Neplan. The single line diagram that relates to
In this case study two WECS were modelled via Neplan. The single line diagram that relates to
each WECS is shown in Figure 13a,b. Each wind DG unit has a rated capacity of 2 MVA and 3.3 kV
each WECS is shown in Figure 13a,b. Each wind DG unit has a rated capacity of 2 MVA and 3.3 kV
rated voltage, and were connected to the investigated feeder via a (20/3.3) kV, 2 MVA transformer at
rated voltage, and were connected to the investigated feeder via a (20/3.3) kV, 2 MVA transformer
the correspondent bus. Each wind DG unit has a rated capacity of 2 MVA and 3.3 kV rated voltage,
at the correspondent bus. Each wind DG unit has a rated capacity of 2 MVA and 3.3 kV rated
and were connected to the investigated feeder via a (20/3.3) kV, 2 MVA transformer at the
voltage, and were connected to the investigated feeder via a (20/3.3) kV, 2 MVA transformer at the
correspondent bus as shown in Figure 13a,b. WECS single line diagrams were built under different
correspondent bus as shown in Figure 13a,b. WECS single line diagrams were built under different
network lawyer in Neplan and were shortened at the main network diagram as Gen_W1 and
network lawyer in Neplan and were shortened at the main network diagram as Gen_W1 and Gen_W2
Gen_W2 to avoid congestion on the main diagram. Both WECS’s were connected to Feeder#2 at buses
to avoid congestion on the main diagram. Both WECS’s were connected to Feeder#2 at buses N2-4 and
N2-4 and N2-7 respectively, this WECS layout and point of connection were the same for the
N2-7 respectively, this WECS layout and point of connection were the same for the following three
following three sub-case studies. Generators PMSM1 and PMSM2 were configured to operate in (P-
Q) mode before and after islanding, and were loaded at 75% of their rated capacity with a unity PF
(i.e., (1.5 + j0) MVA).
Energies 2016, 9, 956 19 of 24

sub-case studies. Generators PMSM1 and PMSM2 were configured to operate in (P-Q) mode before
and after
Energies islanding,
2016, 9, 956 and were loaded at 75% of their rated capacity with a unity PF (i.e., (1.5 + j0) MVA).
19 of 24

(a)

(b)
Figure 13.
Figure 13. WECS
WECS single
single line
line diagram:
diagram: (a)
(a) Unit
Unit 1;
1; and
and (b)
(b) Unit
Unit 2.
2.

Meanwhile, islanding methods were applied to the remaining other DG technologies (diesel, gas,
Meanwhile, islanding methods were applied to the remaining other DG technologies (diesel, gas,
and hydro), however no changes occurred to those units in terms of unit ratings. Therefore, each unit
and hydro), however no changes occurred to those units in terms of unit ratings. Therefore, each unit
had 3 MVA rated capacity and a 20 kV rated voltage.
had 3 MVA rated capacity and a 20 kV rated voltage.
The situation mentioned above regarding the four DG technologies remained unchanged during
The situation mentioned above regarding the four DG technologies remained unchanged during
the simulation of the following three sub-case studies, moreover Table 9 demonstrates DG units
the simulation of the following three sub-case studies, moreover Table 9 demonstrates DG units
arrangements for the purpose of the mixed case.
arrangements for the purpose of the mixed case.
Table 9. DG units’ arrangements for the purpose of the mixed case.
Table 9. DG units’ arrangements for the purpose of the mixed case.
Correspondent Gen. Correspondent Gen. Correspondent Gen. Correspondent
Generator
Name (Wind & Gen.Name
Correspondent (Wind & GasGen. Name
Correspondent (Wind & Gen.
Correspondent Node of
Correspondent
Name
Generator NameDiesel
Name (Wind & Diesel Sub-Case)
Sub-Case) Name (Wind & Gas HydroNameSub-Case)
(Wind & Hydro Connection
Node of
GEN Sub-Case)
Gen_W1 Sub-Case)
Gen_W1 Sub-Case)
Gen_W1 Connection
N2-4
GENGEN Gen_W2Gen_W1 Gen_W2
Gen_W1 Gen_W2
Gen_W1 N2-7
N2-4
GENGEN GEN_D1Gen_W2 Gen_W2
GEN_G1 Gen_W2
GEN_H1 N2-7
N2-9
GENGEN GEN_D2GEN_D1 GEN_G1
GEN_G2 GEN_H1
GEN_H2 N2-9
N2-11
GEN GEN_D2 GEN_G2 GEN_H2 N2-11
GEN GEN5 GEN5 GEN5 N3-2
GEN GEN5 GEN5 GEN5 N3-2
GENGEN
GEN6 GEN6
GEN6GEN6
GEN6GEN6
N4-4
N4-4

6.4.1. Planned Islanding of Two Diesel & Two Wind Units (Case D1)
6.4.1. Planned Islanding of Two Diesel & Two Wind Units (Case D1)
In this sub-case two DG units (Gen_W1, Gen_W2) are two wind based units connected to buses
In this sub-case two DG units (Gen_W1, Gen_W2) are two wind based units connected to buses
N2-4 and N2-7 respectively, while the other two DG diesel based units (Gen_D1 and Gen_D2) are
N2-4 and N2-7 respectively, while the other two DG diesel based units (Gen_D1 and Gen_D2) are
connected to buses N2-9 and N2-11, respectively of the scheme in Figure 1. For simplicity only one
connected to buses N2-9 and N2-11, respectively of the scheme in Figure 1. For simplicity only one
generation situation were investigated (over-generation) during this sub-case and the following two
load sharing sub-cases. Prior to the islanding mode, the two diesel units Gen_D1 and Gen_D2 were
loaded at (1.8 + j0.4) MVA, while wind units were loaded at (1.5 + j0) MVA. Hence, the existing
generation to load ratio was (6.6/6.46) MW prior to islanding, this implies a 2.16% over-generation
power mismatch. The results are shown in Figure 14.
Energies 2016, 9, 956 20 of 24

generation situation were investigated (over-generation) during this sub-case and the following two
load sharing sub-cases. Prior to the islanding mode, the two diesel units Gen_D1 and Gen_D2 were
loaded at (1.8 + j0.4) MVA, while wind units were loaded at (1.5 + j0) MVA. Hence, the existing
generation to load ratio was (6.6/6.46) MW prior to islanding, this implies a 2.16% over-generation
power
Energiesmismatch.
2016, 9, 956 The results are shown in Figure 14. 20 of 24
Energies 2016, 9, 956 20 of 24

Figure 14. Simulation results for planned islanding of two diesel and two wind units (case D1).
Figure 14. Simulation results for planned islanding of two diesel and two wind units (case D1).
Figure 14. Simulation results for planned islanding of two diesel and two wind units (case D1).
6.4.2. Planned Islanding of Two Gas & Two Wind Units (Case D2)
6.4.2.
6.4.2.Planned
PlannedIslanding
Islandingofof Two
TwoGasGas&&Two
TwoWind
WindUnits
Units(Case
(Case D2)
D2)
In this sub-case, loading of the gas and wind units followed sub-case D1 exactly. The results are
shownInInthis
this
in sub-case,
15. loading
sub-case,
Figure loading of
of the
the gas
gas and
and wind units followed
wind units followed sub-case
sub-caseD1
D1exactly.
exactly.The
Theresults
resultsare
are
shown in Figure
shown in Figure 15.15.

Figure 15. Simulation results for planned islanding of two gas and two wind units (case D2).
Figure 15. Simulation results for planned islanding of two gas and two wind units (case D2).
Figure 15. Simulation results for planned islanding of two gas and two wind units (case D2).
6.4.3. Planned Islanding of Two Hydro & Two Wind Units (Case D3)
6.4.3. Planned Islanding of Two Hydro & Two Wind Units (Case D3)
In this sub-case, loading of the gas and wind units followed sub-case D1 exactly. The results are
shown In in
this sub-case,
Figure 16. loading of the gas and wind units followed sub-case D1 exactly. The results are
shown in Figure 16.
Energies 2016, 9, 956 21 of 24

6.4.3. Planned Islanding of Two Hydro & Two Wind Units (Case D3)
In this sub-case, loading of the gas and wind units followed sub-case D1 exactly. The results are
shown
Energiesin Figure
2016, 9, 956 16. 21 of 24

Figure 16. Simulation results for planned islanding of two diesel two wind units (case D3).
Figure 16. Simulation results for planned islanding of two diesel two wind units (case D3).

7. Discussion of the Results

7. Discussion of the Results
Table 10 summarizes the main points and outcomes obtained from this work, hence they give a
Table 10 summarizes the main points and outcomes obtained from this work, hence they give
brief picture of the problem and the key index to understand it. In this table there are comments about
athe
brief picture of the problem and the key index to understand it. In this table there are comments
reasons that the islanding attempt has failed. Voltage and frequency is compared to the limits that
about the reasons
are introduced that the
in Table 1. islanding attempt has failed. Voltage and frequency is compared to the
limits that are introduced in Table 1.
Table 10. Results summary.
Table 10. Results summary.
DG Sub-Case Results Islanding
Islanding Attempt Comment
Technology
DG Number
Sub-Case Figure
Results Scenario
Islanding
Islanding Attempt Comment
Technology Number Figure
A1.1 Figure
7a Scenario
Planned Succeeded -
A1.2A1.1 Figure 7b
Figure 7a Unplanned
Planned Failed
Succeeded Voltage and frequency - limits violation
Diesel
A2.1A1.2 Figure 8a 7b
Figure Planned
Unplanned Succeeded
Failed -
Voltage and frequency limits violation
Diesel
A2.2A2.1 Figure
Figure 8b 8a Planned
Unplanned Succeeded
Failed Voltage and frequency - limits violation
B1.1A2.2 Figure 8b
Figure 9a Unplanned
Planned Failed
Succeeded Voltage and frequency limits violation
-
B1.2B1.1 Figure
Figure 9b 9a Planned
Unplanned Succeeded
Failed Voltage and frequency - limits violation
Gas
Gas B2.1B1.2 Figure
Figure
10a9b Unplanned
Planned Failed
Conditional Success Voltage andlimit
Frequency frequency
adjustedlimits
to tviolation
> 27 s
Frequency limit adjusted to t > 27 s
B2.2B2.1 Figure
Figure 10a
10b Planned
Unplanned Conditional Success
Failed Voltage and frequency limits violation
B2.2 Figure 10b Unplanned Failed Voltage and frequency limits violation
C1.1 Figure 11a Planned Failed Upper & lower frequency limits violation
C1.2C1.1 Figure
Figure 11b11a Planned
Unplanned Failed
Failed Voltage&and
Upper lower frequency
frequency limits
limits violation
violation
Hydro C1.2 Figure 11b Unplanned Failed Voltage and frequency limits violation
Hydro C2.1 Figure 12a Planned Conditional Success Frequency limit adjusted to t > 4.4 s
C2.1 Figure 12a Planned Conditional Success Frequency limit adjusted to t > 4.4 s
C2.2C2.2 Figure 12b12b
Figure Unplanned
Unplanned Failed
Failed Voltage
Voltageand
andfrequency
frequencylimits
limitsviolation
violation
Wind & Islanding may not be recommended if
D1 Figure 14 Planned Succeeded Islanding may not be recommended if
Diesel
Wind & Diesel D1 Figure 14 Planned Succeeded voltage critical loads exist
voltage critical loads exist
Islanding may not be recommended if
Wind & Gas D2 Figure 15 Planned Succeeded Islanding may not be recommended if
Wind & Gas D2 Figure 15 Planned Succeeded voltage critical loads exist
voltage critical loads exist
Wind & Islanding may not be recommended if
D3 Figure 16 Planned Succeeded Islanding may not beloads
recommended if
Hydro
Wind & Hydro D3 Figure 16 Planned Succeeded voltage critical exist
voltage critical loads exist

The planned islanding proved to be feasible option for most of the cases investigated, however
someThe planned
cases’ islanding
islanding proved
feasibility to becertain
needed feasible option forbetween
agreements most of the cases
utility andinvestigated, however
DG units’ operators.
some cases’
Besides theislanding
unstable feasibility
behavior of needed certain
the system agreements
following betweenislanding
unplanned utility and
forDG
all units’
cases operators.
tested, it
could also lead to DG units’ rotor damage, according to the power variations recorded from the
individual unplanned sub-case studies.
There is indeed a pivotal role for the anti-islanding protection schemes embedded within
distribution networks, since the unintentional islanding effect brings severe outcomes on end users’
loads, DG units, and utility assets. This was clear from all of the examined unplanned islanding sub-
Energies 2016, 9, 956 22 of 24

Besides the unstable behavior of the system following unplanned islanding for all cases tested, it could
also lead to DG units’ rotor damage, according to the power variations recorded from the individual
unplanned sub-case studies.
There is indeed a pivotal role for the anti-islanding protection schemes embedded within
distribution networks, since the unintentional islanding effect brings severe outcomes on end users’
loads, DG units, and utility assets. This was clear from all of the examined unplanned islanding
sub-cases, where DG units were given the permission to operate freely without interruption. That was
only to demonstrate the size of negative impact it could leave on real distribution systems with DG
units, whenever unintentional islanding occurs without having the proper anti-islanding detection
and protection schemes.
Frequency variations were related directly to the existing active power mismatch prior to islanding,
and implicitly to the islanding method used, the type of DG technology, and the governor system
associated with them. While voltage variations were function of the existing reactive power mismatch
prior to islanding, islanding method used, and the type DG unit technology; however only one
excitation system was used therefore it was not enough to determine the effect of different excitation
systems on voltage behavior.
Since among the nine investigated planned sub-cases five out of six over-generation sub-cases
and one out of three under-generation sub-cases returned reliable islanding results. While two
under-generation sub-cases required certain agreement between utility and DG units’ operators and
only one over-generation sub-case failed to satisfy the standard islanding operating conditions.
Among the four DG technologies, diesel units alone had the best planned islanding response,
while gas units had the longest response time, and hydro units had the most variation levels. Moreover,
loads had an impact on the frequency response of the system prior islanding and that was obvious in
the three mixed units load sharing sub-cases, while voltage variations had the direct impact on load
dynamic behavior.
The droop characteristic of the gas turbine in gas sub-cases had the advantage of preserving the
system stability during the load sharing method. Meanwhile, a slight improvement on frequency droop
response was observed after using wind with gas and with hydro during the mixed case; however the
stable diesel units’ response were reduced after mixing with wind. This was not the case for voltage
since all mixed cases had stronger voltage dips compared to individual cases.

8. Conclusions
The dynamic behavior of a realistic MV network with existing distributed generation following
a planned and unplanned islanding has been investigated. Four different distributed generation
technologies (diesel, gas, hydro and wind) along with their excitation and governor control systems
have modelled and simulated. Furthermore an exponential model was used to represent the loads in
the network, while the load sharing islanding method has been used for controlling the distributed
generation units during grid-connected and islanding operation. The obtained simulation results were
validated through various case studies and have shown that properly planned islanding transition
could provide support to critical loads at the event of utility outages. The current work contributes the
conducting research on intentional and unintentional islanding in an effort to increase the DG units
and the power distribution networks’ efficiency.

Acknowledgments: Funds have received from City, University of London, Systems & Control Center for covering
the costs for open access publishing.
Author Contributions: Maen Z. Kreishan and Lambros Ekonomou conceived the idea and led the whole work;
Maen Z. Kreishan performed the simulations; George P. Fotis, Maen Z. Kreishan and Vasiliki Vita analyzed the
data; Lambros Ekonomou contributed required material; George P. Fotis and Vasiliki Vita wrote the paper.
Conflicts of Interest: The authors declare no conflict of interest.
Energies 2016, 9, 956 23 of 24

References
1. Vahl, F.P.; Rüther, R.; Fihlo, N.C. The influence of distributed generation penetration levels on energy markets.
Energy Policy 2013, 62, 226–235. [CrossRef]
2. Newbery, D.M. Towards a green energy economy? The EU Energy Union’s transition to a low-carbon zero
subsidy electricity system—Lessons from the UK’s Electricity Market Reform. Appl. Energy 2016, in press.
[CrossRef]
3. Conti, S.; Nicolosi, R.; Rizzo, S.A.; Zeineldin, H.H. Optimal dispatching of distributed generators and storage
systems for MV islanded microgrids. IEEE Trans. Power Deliv. 2012, 27, 1243–1251. [CrossRef]
4. Karimi, M.; Mokhlis, H.; Naidu, K.; Uddin, S.; Bakar, A.H.A. Photovoltaic penetration issues and impacts in
distribution network—A review. Renew. Sustain. Energy Rev. 2016, 53, 594–605. [CrossRef]
5. Franz, P.; Talavera, I.; Hanson, J.; Sgoff, I. Optimized regulation of dispersed generation units for
minimization of reactive power consumption. In Proceedings of the IEEE Innovative Smart Grid Technologies
Conference (ISGT), Washington, DC, USA, 18–20 February 2015; pp. 1–5.
6. Prakash, P.; Khatod, D.K. Optimal sizing and siting techniques for distributed generation in distribution
systems: A review. Renew. Sustain. Energy Rev. 2016, 57, 111–130. [CrossRef]
7. Karatepe, E.; Ugranlı, F.; Hiyama, T. Comparison of single- and multiple-distributed generation concepts in
terms of power loss, voltage profile, and line flows under uncertain scenarios. Renew. Sustain. Energy Rev.
2015, 48, 317–327. [CrossRef]
8. Vita, V.; Ekonomou, L.; Christodoulou, C.A. The impact of distributed generation to the lightning protection
of modern distribution lines. Energy Syst. 2016, 7, 357–364. [CrossRef]
9. Chowdhury, S.P.; Chowdhury, S.; Ten, C.F.; Crossley, P.A. Operation and control of DG based power island
in smart grid environment. In Proceedings of the 20th International Conference and Exhibition on Electricity
Distribution—Part 1, CIRED, Prague, Czech Republic, 8–11 June 2009; pp. 1–5.
10. National Grid. Working Group Report: Small Embedded Generation Issue 1.0; National Grid: Warwick, UK, 2010.
11. Massucco, S.; Pitto, A.; Silvestro, F.A. Gas turbine model for studies on distributed generation penetration
into distribution networks. IEEE Trans. Power Syst. 2011, 26, 992–999. [CrossRef]
12. Balaguer, I.; Lei, Q.; Yang, S.; Supatti, U.; Peng, F.Z. Control for grid-connected and intentional islanding
operations of distributed power generation. IEEE Trans. Ind. Electron. 2011, 58, 147–157. [CrossRef]
13. Chen, Y.; Xu, Z.; Ostergaard, J. Frequency analysis for planned islanding operation in the Danish distribution
system—Bornholm. In Proceedings of the 43rd International Universities Power Engineering Conference,
UPEC, Padova, Italy, 1–4 September 2008; pp. 1–5.
14. Shahabi, M.; Haghifam, M.R.; Mohamadian, M.; Nabavi-Niaki, S. Dynamic behavior improvement in
a microgrid with multiple DG units using a power sharing approach. In Proceedings of the IEEE Bucharest
PowerTech, Bucharest, Romania, 28 June–2 July 2009; pp. 1–8.
15. Mulhausen, J.; Schaefer, J.; Mynam, M.; Guzman, A.; Donolo, M. Anti-islanding today, successful islanding
in the future. In Proceedings of the 63rd Annual Conference for Protective Relay Engineers, College Station,
TX, USA, 29 March–1 April 2010; pp. 1–8.
16. Mohamad, H.; Mokhlis, H.; Bakar, A.H.A.; Ping, H.W. A review on islanding operation and control for
distribution network connected with small hydro power plant. Renew. Sustain. Energy Rev. 2011, 15,
3952–3962. [CrossRef]
17. Trujillo, C.; Velasco, D.; Figueres, E.; Garcera, G. Local and remote techniques for islanding detection in
distributed generators. In Distributed Generation; InTech Publication: Rijeka, Croatia, 2010; Chapter 6.
18. Mahat, P.; Chen, Z.; Bak-Jensen, B. Review of islanding detection methods for distributed generation.
In Proceedings of the 3rd International Conference on Electric Utility Deregulation and Restructuring and
Power Technologies, DRPT, Nanjing, China, 6–9 April 2008; pp. 2743–2748.
19. Maki, K.; Kulmala, A.; Repo, S.; Jarventausta, P. Problems related to islanding protection of distributed
generation in distribution network. In Proceedings of the IEEE PowerTech, Lausanne, Switzerland,
1–5 July 2007; pp. 467–472.
20. Velasco, D.; Trujillo, C.; Garcerá, G.; Figueres, E. Review of anti-islanding techniques in distributed generators.
Renew. Sustain. Energy Rev. 2010, 14, 1608–1614. [CrossRef]
Energies 2016, 9, 956 24 of 24

21. Katiraei, F.; Abbey, C.; Tang, S.; Gauthier, M. Planned islanding on rural feeders-utility perspective.
In Proceedings of the IEEE Power and Energy Society General Meeting—Conversion and Delivery of
Electrical Energy in the 21st Century, Pittsburgh, PA, USA, 20–24 July 2008; pp. 1–6.
22. Bründlinger, R.; Bletterie, B. Unintentional islanding in distribution grids with a high penetration of
inverter-based DG: Probability for islanding and protection methods. In Proceedings of the IEEE Russia
Power Tech, St. Petersburg, Russia, 27–30 June 2005; pp. 1–7.
23. IEEE Standard for Interconnecting Distributed Resources with Electric Power Systems; IEEE: New York, NY,
USA, 2003.
24. Basso, T.S.; DeBlasio, R. IEEE 1547 Series of Standards: Interconnection issues. IEEE Trans. Power Electron.
2004, 19, 1159–1162. [CrossRef]
25. Ozeer, B.A.; Hernandez-Gonzalez, G.; El-Fouly, T.H. Investigation of Planned Islanding Performance of Rotating
Machine-Based DG Technologies, 2nd ed.; CanmetENERGY: Ottawa, ON, Canada, 2012.
26. He, J.; Li, Y.W. An accurate reactive power sharing control strategy for DG units in a microgrid. In Proceedings
of the IEEE 8th International Conference on Power Electronics and ECCE Asia (ICPE & ECCE), Seogwipo,
Korea, 30 May–3 June 2011; pp. 551–556.
27. He, J.; Li, Y.W.; Guerrero, J.; Blaabjerg, F.; Vasquez, J. An islanding microgrid power sharing approach using
enhanced virtual impedance control scheme. IEEE Trans. Power Electron. 2013, 28, 5272–5282. [CrossRef]
28. Conti, S.; Greco, A.; Messina, N.; Vagliasindi, U. Intentional islanding of MV microgrids: Discussion
of a case study and analysis of simulation results. In Proceedings of the International Symposium on
Power Electronics, Electrical Drives, Automation and Motion, SPEEDAM, Catania, Italy, 11–13 June 2008;
pp. 422–427.
29. Navarro, I.R. Dynamic Load Models for Power Systems-Estimation of Time-Varying Parameters during
Normal Operation. Ph.D. Thesis, Lund University, Lund, Sweden, September 2002.
30. Stojanović, D.P.; Korunović, L.M.; Milanović, J. Dynamic load modelling based on measurements in medium
voltage distribution network. Electr. Power Syst. Res. 2008, 78, 228–238. [CrossRef]
31. Li, Y.; Chiang, H.D.; Choi, B.K.; Chen, Y.T.; Huang, D.H.; Lauby, M.G. Load models for modeling dynamic
behaviors of reactive loads: Evaluation and comparison. Int. J. Electr. Power Energy Syst. 2008, 30, 497–503.
[CrossRef]
32. Achilles, S.; Poller, M. Direct drive synchronous machine models for stability assessment of wind farms.
In Proceedings of the 4th International Workshop on Large Scale Integration of Wind Power and Transission
Networks for Offshore Wind Farms, Billund, Denmark, 20–21 October 2003; pp. 1–9.
33. Reddy, P.B.; Hiskens, I.A. Limit-induced stable limit cycles in power systems. In Proceedings of the IEEE
Russia Power Tech, St. Petersburg, Russia, 27–30 June 2005.
34. Li, H.; Chen, Z. Design optimization and comparison of large direct-drive permanent magnet wind generator
systems. In Proceedings of the International Conference on Electrical Machines and Systems, Seoul, Korea,
8–11 October 2007.
35. Babu, B.; Mohanty, K. Doubly-fed induction generator for variable speed wind energy conversion
systems-modeling & simulation. Int. J. Comput. Electr. Eng. 2010, 2, 1793–8163.
36. Karimi, H.; Davison, E.; Iravani, R. Multivariable servomechanism controller for autonomous operation of
a distributed generation unit: Design and performance evaluation. IEEE Trans. Power Syst. 2010, 25, 853–865.
[CrossRef]
37. Bahrani, B.; Karimi, H.; Iravani, R. Stability analysis and experimental validation of a control strategy for
autonomous operation of distributed generation units. In Proceedings of the International Power Electronics
Conference (IPEC), Sapporo, Japan, 21–24 June 2010; pp. 464–471.

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC-BY) license (http://creativecommons.org/licenses/by/4.0/).