Control Engineering Practice: Youssef Krim, Dhaker Abbes, Saber Krim, Mohamed Faouzi Mimouni
Control Engineering Practice: Youssef Krim, Dhaker Abbes, Saber Krim, Mohamed Faouzi Mimouni
1. Introduction                                                                                   Renewable Distributed Generators (RDGs) (Jia, Mu, & Qi, 2014). The
                                                                                                  wind technology application stresses the Batteries (BT) storage system,
     The significant increase in the population and areas of using electri-                       because it uses a large fraction of the energy stored in it (Matthieu,
cal energy expected in the last few years must cope with the growing                              Arnaud, Karl, Moe, Marc, & Richard, 2017). The BT lifespan is limited
energy consumption associated with it. To do this, this evolution                                 by the charge/discharge cycle. Compared to BT, the Super-Capacitors
must result in intelligent management of an electrification network                               (SC) storage system life is much longer and has much higher power
(Petronela et al., 2016). Therefore, the power grid is becoming more                              density. Moreover, SC can provide a fast and effective energy output
and more interconnected and meshed. One of the envisaged solutions
                                                                                                  because of its high power density and high efficiency (Wang, Liu,
concerns the reinforcement of the electricity supply network by the
                                                                                                  Pan, & Chen, 2017). Thus, the adoption of SC in HESS is an effective
integration of innovative energy storage devices, new renewable and
                                                                                                  solution to prolong the lifespan of BT in renewable energy production
relocated production means, as well as an energy optimization of the
                                                                                                  applications. The HESS technology can be passive, semi-active or purely
grid architecture (Glasnovic & Margeta, 2011). The wind technology
has become a favored form of renewable energy technologies because                                active. Passive HESS is the simplest configuration, such as SC and
it is seen as clean and sustainable (Islam, Djamel, Abdel-Moumen, Bilal,                          BT directly connected to a DC bus without any power converter, and
& Ikram. El, 2018). The availability of wind energy depends on the                                is characterized by its low cost. However, the performance of this
climatic and geographical contexts of the installation region, which                              configuration is limited because SC cannot be used effectively (Ziyou,
represents the major problem of wind turbines. Therefore, it is naturally                         Jun, Heath, Jianqiu, & Minggao, 2017). Semi-active HESS uses a single
very unlikely to have a concomitance between production and demand.                               power converter. This topology combines good performance and low
In order to guarantee production, the aggregation of substation storage                           system cost (Ziyou et al., 2017). The purely active HESS configuration
units will make it possible to overcome this problem. Hybrid Energy                               uses two DC/DC converters. It gives the possibility to control SC and
Storage Systems (HESSs) have become more and more important in                                    BT currents simultaneously. Therefore, this configuration combines
    ∗ Corresponding author.
      E-mail address: krim_enim@hotmail.com (Y. Krim).
https://doi.org/10.1016/j.conengprac.2018.09.013
Received 21 January 2018; Received in revised form 8 July 2018; Accepted 12 September 2018
Available online xxxx
0967-0661/© 2018 Elsevier Ltd. All rights reserved.
Y. Krim et al.                                                                                                      Control Engineering Practice 81 (2018) 215–230
efficiency and good performance. For this reason, the purely active                3. Renewable distributed generator model and control
configuration is focused on in this study.
    The key problem with HESS is its protection against over-discharge             3.1. Model of wind generator
and overcharge, as well as the BT protection against the rapid fluctu-
ation of charge/discharge cycles using SC. For this reason, the major                 The aerodynamic power produced by the turbine can be expressed
task is to optimize the power flow management between BT and SC,                   by (Youssef, Saber, & Mohamed, 2017):
which manages the reference current for each storage device. The most
common of Energy Management Strategies (EMSs) are rule-based EMSs                  𝑃𝑎𝑒𝑟 = 0.5𝜌𝜋𝑅2 𝑣3 𝐶𝑝 (𝜆, 𝛽)                                                (1)
(Song et al., 2014). However, focusing on renewable energy, if RDG has
                                                                                   where R is the radius of the rotation of each blade, 𝜌 ≈ 1, 2 kg/m3
to operate on different modes, rule development will become complex
                                                                                   denotes the air density, v is the wind speed, and 𝐶𝑝 can then be
due to several variables and case studies. In this framework, fuzzy logic
                                                                                   represented as a polynomial function of 𝜆 and 𝛽. In our case, this
appears as a promising tool owing to the possibility of avoiding most
                                                                                   polynomial is defined by (Youssef et al., 2017):
of the mathematical stiffness and complexity in problem formulation                                 (                                 )   (       )
and representing it based on human reasoning (Victor et al., 2016).                                     151                                  18.4
                                                                                   𝐶𝑝 (𝜆, 𝛽) = 0.53         − 0.58𝛽 − 0.002𝛽 2.14 − 10 exp −        (2)
EMS based on fuzzy logic is used in several applications, such as the                                    𝜆𝑖                                   𝜆𝑖
energy management of HESS used in a tramway (Victor et al., 2016) and                            1
                                                                                   𝜆𝑖 =                                                             (3)
EMS for standalone RDGs (Erdinc & Uzunoglu, 2011). In this study, EMS                       1
                                                                                                  − 0.003
                                                                                         𝜆−0.02𝛽    𝛽 3 +1
based on a Fuzzy Logic Supervisor (FLS) is proposed to monitor RDG in
standalone and grid-connected modes. The purpose of suggested FLS is                   We deduce the expression of the generated aerodynamic torque as
to ensure a balance between production and consumption, to guarantee               follows:
a continuous load supply by maintaining the State Of Charge (SOC)                            𝑃𝑎𝑒𝑟   0.5𝜌𝜋𝑅2 𝑣3 𝐶𝑝 (𝜆, 𝛽)
of SC (𝑈𝑠𝑐 ) and BT (SOCbat ) at acceptable levels, and to assist RDG to           𝑇𝑎𝑒𝑟 =         =                                                           (4)
                                                                                             𝛺𝑚           𝛺𝑚
contribute to the improvement of electrical network performance.
    On the other hand, the studied RDG is associated with a load to                  This torque makes it possible to turn the rotor of PMSG, which is
form an Active Generator (AG). AG can operate according to the grid                modeled in the park frame by the following expression (Masmoudi,
stability in two main operation modes. First, the grid connection mode             Abdelkafi, & Krichen, 2011):
participates in system services by adjusting the frequency and amplitude                    𝑑𝑖𝑠𝑑
of the grid voltage. Second, the standalone mode ensures a continuous              𝑉𝑠𝑑 = 𝐿𝑠      + 𝑅𝑠 𝑖𝑠𝑑 − 𝑝𝛺𝑚 𝑖𝑠𝑞
                                                                                             𝑑𝑡
power supply of the load in case of a grid fault. There is an intermediate                  𝑑𝑖𝑠𝑞                                                              (5)
mode to safety reconnecting AG to the utility grid. Thus, some research            𝑉𝑠𝑞 = 𝐿𝑠      + 𝑅𝑠 𝑖𝑠𝑞 + 𝑝𝛺𝑚 𝑖𝑠𝑑 + 𝑝𝛺𝑚 𝜙𝑚
                                                                                             𝑑𝑡
work has processed the regulation of the frequency and voltage in a                𝑇𝑒𝑚 = 𝑝𝜙𝑚 𝑖𝑠𝑑
connected or standalone mode. Table 1 show that the droop control
method is commonly employed in RDG, which has more than one                            Fig. 2 illustrates a vector control of PMSG in the PARK frame. This
operation mode.                                                                    strategy consists in keeping axis d constantly aligned with the flux vector
    Consequently, this study is focused on a novel control technique               of the magnet. The reference of direct current 𝑖𝑑 is kept at zero. The
based on EMS using FLS and adaptive FLDC, capable of monitoring AG                 reference for quadratic current 𝑖𝑞 is determined by the electromagnetic
in different operation modes and injecting active and reactive powers              torque deduced by the Maximum Power Point Tracking (MPPT) strategy
into the grid with high precision ensuring its stability.                          as follows (Krim, Abbes, Krim, & Mimouni, 2017):
    This paper is structured as follows. Section 2 provides an overview                                                                       5 3
                                                                                            𝑇𝑒𝑚−𝑀𝑃 𝑃 𝑇                  𝑃         1 𝜌𝐶𝑝 max 𝑅 𝛺 𝑚
of the wind hybrid generator configuration. Section 3 presents the                 𝑖∗𝑠𝑞 =              𝑤𝑖𝑡ℎ 𝑇𝑒𝑚−𝑀𝑃 𝑃 𝑇 = 𝑀𝑃 𝑃 𝑇 =                             (6)
                                                                                               𝑝𝜙𝑚                        𝛺𝑚      2           3
modeling and control of RDG. Section 4 discusses the working principle                                                                   𝜆𝑜𝑝𝑡
of EMS based on FLS. Section 5 focuses on the FLDC technique. The
simulation results and the discussion are given in Section 6. Section 7            3.2. BT model and control
concludes this report.
                                                                                      We choose for a main storage unit a BT bank, which must be
2. Active generator configuration                                                  connected to the DC bus through a bidirectional converter (DC/DC
                                                                                   converter 3). We have retained the recent technology of the CIEMAT
    Before undertaking modeling, we need to define in more detail the              (Research Center for Energy, Environment and Technology, Espagne)
architectures of the different parts that make up RDG, detailed in Fig. 1:         BT model, which has very high energy density. An equivalent model for
a wind generator, SC, BT, and a DC bus link. RDG is associated with                BT is depicted in Fig. 3. It is composed by voltage source 𝐸𝑏 in series
a balanced AC load to form AG. The latter can operate in a grid-                   with internal resistance 𝑅𝑖 .
connected mode or an islanded one, according to the stability status                  The expressions of BT quantities are expressed in the following
of the main power grid. In addition, we need to choose the most                    (Cabrane, Ouassaid, & Maaroufi, 2017):
appropriate converter type to control each AG component, which will                   − The general expression of the BT voltage is as follows:
ensure the adaptation of this one to the continuous DC bus. The three-
bladed wind generator absorbs mechanical power 𝑃𝑚 , captured by wind               𝑉𝑏𝑎𝑡 = 𝑛𝑏 𝐸𝑏 + 𝑛𝑏 𝑅𝑖 𝑖𝑏𝑎𝑡                                                  (7)
speed v passing through surface S covered by its blades. The turbine
                                                                                       − The expression of SOC is as follows:
then generates torque 𝑇𝑚 , which will drive at angular velocity 𝜔𝑚 the
rotor of a rotating machine. This transforms the absorbed mechanical                              𝑄𝑑
                                                                                   𝑆𝑂𝐶 = 1 −                                                                  (8)
power into exploitable electrical power. We choose to use a Permanent                             𝐶𝑏𝑎𝑡
Magnets Synchronous Generator (PMSG), able to operate at different                     Capacity 𝐶𝑏𝑎𝑡 of BT is expressed as a function of charging and
speeds until it stops without stalling, which require little maintenance.          discharging current 𝑖𝑏𝑎𝑡 .
To complete the power lack or excess, BT/SC HESS having two types
                                                                                   𝐶𝑏𝑎𝑡             1.67
of storage units is suggested. The role of this AG is double: maintaining               =             ( )0.9 (1 + 0.005𝛥𝑇 )                                   (9)
the frequency and voltage of the main grid in the desired stability range          𝐶10                  𝑖
                                                                                             1 + 0.67 𝑖𝑏𝑎𝑡
and ensuring the continuous supply of the load in case of a grid failure.                                10
Indeed, a new control strategy based on FLS and FLDC will be applied               where charge/discharge current 𝑖10 corresponds to rated capacity 𝐶10 .
to make the studied AG intelligent and able to meet the aforementioned                Fig. 5(a) represents the control of DC/DC converter 3 and the BT
objectives.                                                                        current:
                                                                             216
Y. Krim et al.                                                                                 Control Engineering Practice 81 (2018) 215–230
                                                       217
Y. Krim et al.                                                                                                                    Control Engineering Practice 81 (2018) 215–230
                 Table 1
                 Summary of different droop control methods. (See Furtado et al. (2008), Guerrero et al. (2006), Hou et al. (2015), Il-Yop et al. (2010), Juan et al.
                 (2009), Lee et al. (2009), Matas et al. (2010), Sahyoun et al. (2015), Shuai et al. (2016), Tayaba et al. (2017), Tolani and Sensarma (2012).).
                                                                                                        𝑈𝑚𝑏𝑎𝑡
                                                                                               𝑚𝑏𝑎𝑡 =                                                                      (11)
                                                                                                        𝑈𝐷𝐶
                                                                                                   Faced with the slow dynamics and the risk of the premature wear of
                                                                                               BT, we add a secondary storage unit to the system, intended to absorb
                                                                                               or provide power peaks. The technology of SC is chosen. The model of
                 Fig. 3. Equivalent model of nb BT elements in series.
                                                                                               SC used in this work is presented in Fig. 4(a). An SC module composed
                                                                                               of multiple in-series and parallel SC is thus modeled by capacitance 𝐶𝑠𝑐
                                                                                               in series with resistance 𝑅𝑠𝑐 (Fig. 4(b)):
    − A PI controller is used to control the charge and discharge current
                                                                                                            𝑁𝑠            𝑁𝑝
of BT and make it equivalent to its reference ‘‘𝑖𝑏𝑎𝑡−𝑟𝑒𝑓 ’’:                                   𝑅𝑠𝑐 = 𝑅𝑠𝑠𝑐      ; 𝐶 = 𝐶𝑠𝑠𝑐                                                  (12)
                                                                                                            𝑁𝑝 𝑠𝑐         𝑁𝑠
                  (                )                                                           where 𝐶𝑠𝑠𝑐 is the nominal capacitance and 𝑅𝑠𝑠𝑐 is the equivalent in-series
𝑈𝑚𝑏𝑎𝑡 = 𝑉𝑏𝑎𝑡 − 𝑃 𝐼 𝑖𝑏𝑎𝑡−𝑟𝑒𝑓 − 𝑖𝑏𝑎𝑡                                                (10)         resistance.
                                                                                         218
Y. Krim et al.                                                                                                        Control Engineering Practice 81 (2018) 215–230
   Fig. 5(b) illustrates the control of DC/DC converter 2 and the SC                       If 𝑃𝑠𝑡𝑜 is positive and 𝑆𝑂𝐶 𝑏𝑎𝑡 ≥90% then 𝑃𝑠𝑐−𝐶𝑆3 = 𝑃𝑠𝑡𝑜 else 𝑃𝑠𝑐−𝐶𝑆3 =
current:                                                                               𝑃𝑠𝑐−𝑠𝑡𝑜
   − Current control: A PI regulator is used to maintain the charging                      If 𝑃𝑠𝑡𝑜 is negative and 𝑆𝑂𝐶 𝑏𝑎𝑡 ≤ 30% then 𝑃𝑠𝑐−𝐶𝑆3 = 𝑃𝑠𝑡𝑜 else
and discharging current of SC equivalent to its reference value ‘‘𝑖𝑠𝑐−𝑟𝑒𝑓 ’’.          𝑃𝑠𝑐−𝐶𝑆3 = 𝑃𝑠𝑐−𝑠𝑡𝑜
This PI controller is expressed as follows:                                                To elaborate the reference powers of BT and SC, the SOC of each
                 (             )                                                       storage system must be taken in account. Switches 2 and 4 are used
𝑈𝑚𝑠𝑐 = 𝑈𝑠𝑐 − 𝑃 𝐼 𝑖𝑠𝑐−𝑟𝑒𝑓 − 𝑖𝑠𝑐                                        (13)
                                                                                       to extract the exact reference powers of SC and BT by maintaining the
   − Converter control: DC/DC converter 2 permits adapting the SC                      SOC of every storage system into acceptable levels, to protect them from
output voltage to the adequate inverter input voltage. The control of                  overcharge and discharge excess.
DC/DC converter 2 is given by the duty ratio as follows:                                   ∙ Switch 2 enables choosing between 0 and 𝑃𝑏𝑎𝑡−𝐶𝑆1 :
        𝑈𝑚𝑠𝑐                                                                               If 𝑃𝑏𝑎𝑡−𝐶𝑆1 is positive and 𝑆𝑂𝐶 𝑏𝑎𝑡 ≥𝑆𝑂𝐶 𝑏𝑎𝑡−𝑚𝑎𝑥 then 𝑃𝑏𝑎𝑡−𝑟𝑒𝑓 = 0 else
𝑚𝑠𝑐 =                                                                     (14)         𝑃𝑏𝑎𝑡−𝑟𝑒𝑓 = 𝑃𝑏𝑎𝑡−𝐶𝑆1 .
        𝑈𝐷𝐶
                                                                                           If 𝑃𝑏𝑎𝑡−𝐶𝑆1 is negative and 𝑆𝑂𝐶 𝑏𝑎𝑡 ≤ SOC 𝑏𝑎𝑡−𝑚𝑖𝑛 then 𝑃𝑏𝑎𝑡−𝑟𝑒𝑓 = 0 else
    Finally, the DC bus voltage is controlled according to the principle               𝑃𝑏𝑎𝑡−𝑟𝑒𝑓 = 𝑃𝑏𝑎𝑡−𝐶𝑆1 .
described in Fig. 5(c). The PI regulator calculates the reference of                       ∙ Switch 4 allows selecting between 0 and 𝑃𝑠𝑐−𝐶𝑆3 :
DC bus current 𝑖∗𝐷𝐶 to maintain the DC bus voltage equivalent to its                       If 𝑃𝑠𝑐−𝐶𝑆3 is positive and 𝑈𝑠𝑐 ≥𝑈𝑠𝑐𝑚𝑎𝑥 then 𝑃𝑠𝑐−𝑟𝑒𝑓 = 0 else 𝑃𝑠𝑐−𝑟𝑒𝑓 =
reference 𝑈𝐷𝐶∗ . The DC bus voltage can be determined by the following
                                                                                       𝑃𝑠𝑐−𝐶𝑆3 .
relationship:                                                                              If 𝑃𝑠𝑐−𝐶𝑆3 is negative and 𝑈𝑠𝑐 ≤ 𝑈 𝑠𝑐𝑚𝑖𝑛 then 𝑃𝑠𝑐−𝑟𝑒𝑓 = 0 else 𝑃𝑠𝑐−𝑟𝑒𝑓 =
  𝑑𝑈𝐷𝐶                                                                                 𝑃𝑠𝑐−𝐶𝑆3 .
𝐶        = 𝑖𝑚𝑔 − 𝑖𝑚𝑠𝑐 − 𝑖𝑚𝑏𝑎𝑡 − 𝑖min 𝑣                                    (15)
    𝑑𝑡                                                                                     Indeed, the low rate of the BT charge is 𝑆𝑂𝐶 𝑏𝑎𝑡−𝑚𝑖𝑛 =30% and the
where 𝑖𝑚𝑔 , 𝑖𝑚𝑠𝑐 , 𝑖𝑚𝑏𝑎𝑡 , and 𝑖𝑚𝑖𝑛𝑣 represent the currents modeled by the             high rate of the BT charge is 𝑆𝑂𝐶 𝑏𝑎𝑡−𝑚𝑎𝑥 =90%. Also, the low rate of
wind generator, SC, BT, and the load and the grid, respectively. C is the              the SC charge is 𝑈 𝑠𝑐𝑚𝑖𝑛 =58 V and the high rate of the SC charge is
capacity of the DC bus.                                                                𝑈 𝑠𝑐𝑚𝑎𝑥 =98 V.
    The HESS can also be used to regulate the DC bus voltage. In this
case, BT and SC power references ‘‘𝑃𝑠𝑐−𝑟𝑒𝑓 ’’ and ‘‘𝑃𝑏𝑎𝑡−𝑟𝑒𝑓 ’’ are calculated         4.2. FLS strategy
by considering the ‘‘𝑃𝑔 ’’ fluctuation wind power and the required power
for the DC bus voltage regulation (𝑃DC   ∗ ).                                              EMS based on FLS has two main objectives. The first one is to control
    The reference currents of BT and SC (respectively 𝑖𝑠𝑐−𝑟𝑒𝑓 and 𝑖𝑏𝑎𝑡−𝑟𝑒𝑓 )           the power flow and minimize the number of charging/discharging cycles
are delivered by FLS.                                                                  of BT to improve their lifespan. The second one is to select the exact
                                                                                       reference powers of BT and SC by controlling the SOC of each storage
4. Fuzzy logic supervisor                                                              system and make them in acceptable margins.
                                                                                           ∙ Fuzzy logic for power flow control
4.1. FLS structure                                                                         The power flow control by fuzzy logic (Fuzzy toolbox ‘‘Power flow
                                                                                       control’’) includes three inputs and two outputs, as shown in Fig. 6.
    The principle of determining the reference currents of BT and SC is                The inputs are the difference between the wind generated power and
detailed in Fig. 6. A Low-Pass Filter (LPF) is used to filter the power                the required power (𝑃𝑠𝑡𝑜 ), the SOC of BT ‘‘𝑆𝑂𝐶 𝑏𝑎𝑡 ’’ and the SOC of SC
quantity to be stored (𝑃𝑠𝑡𝑜 ). The purpose of this LPF is to construct                 ‘‘𝑈𝑠𝑐 ’’. The outputs are the command signal of switch 1 ‘‘CS1’’ and the
the ‘‘ 𝑃𝑏𝑎𝑡−𝑠𝑡𝑜 ’’ power of BT and to divert the rapid fluctuations in ‘‘              command signal of switch 3 ‘‘CS3’’.
𝑃𝑠𝑐−𝑠𝑡𝑜 ’’ power into SC. However, it reduces the peak power demand                        The determination of the membership functions for the fuzzification
and the charging/discharging cycle on BT. After that, the SC power is                  of the input and output variables of the energy supervisor are an
determined by the difference between 𝑃𝑠𝑡𝑜 and 𝑃𝑏𝑎𝑡−𝑠𝑡𝑜 .                               important phase of the fuzzy algorithm.
    To select the exact reference powers of BT and SC with the considered                  In this work, the membership functions have been chosen in an
SOC of each storage system, four switches are used. They are controlled                empiric manner relying on our expertise on the system.
by fuzzy logic as a function of CS1, CS2, CS3 and CS4.                                     The membership functions for the three input variables and for the
    Switches 1 and 3 are responsible for improving the system perfor-                  two output variables must be defined (Fig. 7). Since the number of
mance, in terms of dynamic behavior of BT and their lifespan.                          fuzzy rules depends on the number of fuzzy sets of inputs, we will only
    ∙ Switch 1 allows selecting between 𝑃𝑠𝑡𝑜 and 𝑃𝑏𝑎𝑡−𝑠𝑡𝑜 :                            consider the sets relevant to the case study.
    If 𝑃𝑠𝑡𝑜 is positive and 𝑈𝑠𝑐 ≥𝑈𝑠𝑐𝑚𝑎𝑥 then 𝑃𝑏𝑎𝑡−𝐶𝑆1 = 𝑃𝑠𝑡𝑜 else 𝑃𝑏𝑎𝑡−𝐶𝑆1 =               For SOC (𝑆𝑂𝐶 𝑏𝑎𝑡 and 𝑈𝑠𝑐 ), the membership functions consist of three
𝑃𝑏𝑎𝑡−𝑠𝑡𝑜                                                                               levels (‘‘L’’, ‘‘M’’, ‘‘H’’) which correspond to the three operating modes.
    If 𝑃𝑠𝑡𝑜 is negative and 𝑈𝑠𝑐 ≤ 𝑈 𝑠𝑐𝑚𝑖𝑛 then 𝑃𝑏𝑎𝑡−𝐶𝑆1 = 𝑃𝑠𝑡𝑜 else 𝑃𝑏𝑎𝑡−𝐶𝑆1 =         The ‘‘L’’ and ‘‘H’’ sets ensure the storage availability while avoiding
𝑃𝑏𝑎𝑡−𝑠𝑡𝑜                                                                               low and high saturations. The ‘‘M’’ set is used to compensate for the
    ∙ Switch 3 permits choosing between 𝑃𝑠𝑡𝑜 and 𝑃𝑠𝑐−𝑠𝑡𝑜 :                             oversubscribed power and to store the surplus of renewable production.
                                                                                 219
Y. Krim et al.                                                                                                 Control Engineering Practice 81 (2018) 215–230
This set also makes it possible to ensure the adjustment of the long-term            For the membership functions of the command signals of switches 1
and short-term storage instructions.                                              and 3 (CS1 and CS3), the choice of sets is such possible values of the
                                                                                  output variable are in the interval [−1, 1]. These functions are based on
    The membership functions of the power deviation 𝑃𝑠𝑡𝑜 can be ‘‘N’’ or          two levels: N stands for negative and P stands for positive.
‘‘P’’, where N and P respectively represent the discharging and charging             The used inference matrix is described by Tables 2 and 3 to deter-
of the storage system.                                                            mine the rules of control in order to associate the fuzzy inputs and
                                                                            220
Y. Krim et al.                                                                                                             Control Engineering Practice 81 (2018) 215–230
Table 2                                                                                  Table 4
CS1 rules.                                                                               CS2 rules.
CS1                                 𝑈𝑠𝑐                                                  CS2                                      𝑆𝑂𝐶𝑏𝑎𝑡
                                    L                      M                   H                                                  L                  M              H
 𝑃𝑠𝑡𝑜              P                N                      N                   P          𝑃𝑏𝑎𝑡−𝐶𝑆1                 P              N                  N              P
                   N                P                      N                   N                                   N              P                  N              N
Table 3                                                                                  Table 5
CS3 rules.                                                                               CS4 rules.
CS3                                𝑆𝑂𝐶𝑏𝑎𝑡                                                CS4                                       𝑈𝑠𝑐
                                   L                   M                   H                                                       L                  M                 H
 𝑃𝑠𝑡𝑜              P               N                   N                   P              𝑃𝑠𝑐−𝐶𝑆3                  P               N                  N                 P
                   N               P                   N                   N                                       N               P                  N                 N
outputs. In fact, Table 2 presents the rules between 𝑃𝑠𝑡𝑜 and 𝑈𝑠𝑐 to                     5. Fuzzy logic droop control technique
determine the control signal of switch 1, and Table 3 presents the rules
between 𝑃𝑠𝑡𝑜 and 𝑆𝑂𝐶 𝑏𝑎𝑡 to determine the control signal of switch 3.
    ∙ Selection of the exact reference power of BT and SC                                   In this section, based on the power management supervisor outputs,
    A fuzzy logic used to select the exact reference powers of BT and                    adaptive FLDC is proposed, which is able to inject or absorb active and
SC (Fuzzy toolbox ‘‘BT and SC control’’) includes four inputs and two                    reactive power ‘‘ 𝑃𝑓 and 𝑄𝑣 ’’ into the grid with high accuracy and ensure
outputs, as shown in Fig. 6. The inputs (Fig. 8) are the BT power                        a continuous supply of load.
(𝑃𝑏𝑎𝑡−𝐶𝑆1 ), the SOC of BT (𝑆𝑂𝐶 𝑏𝑎𝑡 ), the SOC of SC (𝑈𝑠𝑐 ), and the                        In the inverter output, the active and reactive powers transferred to
SC power (𝑃𝑠𝑐−𝐶𝑆3 ). The outputs (Fig. 8) are the command signals of
                                                                                         the load and the grid can be expressed as follows (Vasquez, Guerrero,
switches 2 and 4 (CS2 and CS4) to determine the exact reference power
                                                                                         Luna, Rodriguez, & Teodorescu, 2009):
of SC and BT while respecting these limits of charging and discharging.
The purpose of this fuzzy-logic toolbox is to maintain the SOC of SC and                              𝑉𝑐2    ( (                  )               )
BT at acceptable levels.                                                                  𝑃 =                 𝑅𝑔 𝑉𝑐 − 𝑉𝑝𝑐𝑐 cos (𝛿) + 𝑋𝑉𝑝𝑐𝑐 sin (𝛿)
                                                                                                 𝑅2𝑔 + 𝑋 2
    The membership functions of 𝑃𝑏𝑎𝑡-CS1 and 𝑃𝑠𝑐-CS3 are based on two                                                                                               (16)
levels, N and P, to accommodate the needs of the proposed strategy,                                   𝑉𝑐2    ( (                 )               )
                                                                                          𝑄=                  𝑋 𝑉𝑐 − 𝑉𝑝𝑐𝑐 cos (𝛿) − 𝑅𝑉𝑝𝑐𝑐 sin (𝛿)
where N and P respectively represent the discharge and charge of each                            𝑅2𝑔 + 𝑋 2
storage system. The output membership functions present the command
of switch 2 by the CS2 command signal and the command of switch 4                            The powers passing through the line depend on the reactance of the
by the CS4 command signal. These functions are based on two levels: N                    line, the voltage levels and the voltage angles. After using the orthogonal
stands for negative and P stands for positive.                                           linear rotational T transformation matrix, the active and reactive powers
    The used inference matrix is described by Tables 4 and 5 to deter-                   become:
mine the rules of control to associate the fuzzy inputs and outputs.                     ( ′)        ( )
                                                                                           𝑃           𝑃
Actually, Table 4 presents the rules between 𝑃𝑏𝑎𝑡-CS1 and 𝑆𝑂𝐶 𝑏𝑎𝑡 to                         ′   =𝑇                                                             (17)
                                                                                           𝑄           𝑄
determine the control signal of switch 2, and Table 5 presents the rules
between 𝑃𝑠𝑐-CS3 and 𝑈𝑠𝑐 to determine the control signal of switch 4.                     with:
    In this paper, the Mandani method is chosen for defuzzification using
the center of gravity (Arun & Mohan, 2017; Nikita & Garg, 2017). It is                     ⎛𝑋         −
                                                                                                       𝑅⎞
                                                                                           ⎜           𝑍⎟
                                                                                         𝑇 ⎜𝑍
                                                                                                      𝑋 ⎟⎟
intuitive, it has widespread acceptance and it is well suited to human
inputs.                                                                                    ⎜𝑅
                                                                                           ⎝𝑍         𝑍 ⎠
                                                                                   221
Y. Krim et al.                                                                                                         Control Engineering Practice 81 (2018) 215–230
                                                                               222
Y. Krim et al.                                                                                                      Control Engineering Practice 81 (2018) 215–230
Table 6                                                                           Table 7
Rules of frequency control.                                                       Rules of voltage control.
𝛥𝜔𝑡              𝑒𝑝                                                               𝛥V               𝑒𝑞
                 NBB      NB   NM   NS   Z      PS    PM      PB     PBB                           NBB        NB   NM    NS     Z      PS     PM     PB     PBB
         N       PS       PS   PS   PS   PS     PS    PS      PS     PS                     N      PS         PS   PS    PS     PS     PS     PS     PS     PS
 𝑑𝑒𝑝     M       PBB      PB   PM   PS   Z      NS    NM      NB     NBB           𝑑𝑒𝑞      M      PBB        PB   PM    PS     Z      NS     NM     NB     NBB
         P       NS       NS   NS   NS   NS     NS    NS      NS     NS                     P      NS         NS   NS    NS     NS     NS     NS     NS     NS
NBB: Negative Big Big, NB: Negative Big, NM: Negative Medium, NS:                 smooth transition between different operation modes using automatic
Negative Small, Z: Zero, PBB: Positive Big Big, PB: Positive Big, PM:             switches (S5, S6, S7 and S8). They are controlled by fuzzy logic.
Positive Medium, and PS: Positive Small.                                             ∙ Switches S5 and S6 permit selecting between 𝑉𝑛 and 𝑉𝑝𝑐𝑐 , and 𝑓𝑛
    For the same reason, the forms of the membership functions of the             and 𝑓𝑝𝑐𝑐 , respectively.
output variables are also symmetrical. However, we introduce nine                    If 𝐴𝑖𝑠𝑙 is positive then S5 and S6 are in states 𝑉𝑝𝑐𝑐 and 𝑓𝑝𝑐𝑐 ,
fuzzy subsets.                                                                    respectively. Else, S5 and S6 are in states 𝑉𝑛 and 𝑓𝑛 , respectively.
    The membership functions of the error derivatives are chosen to                  ∙ Switch S7 enables choosing between 0 and 𝑅2 .
specify the system response. Nevertheless, we introduce three fuzzy                  If 𝐴𝑣 is positive then 𝑄𝑣 = R2 ; else, 𝑄𝑣 =0.
subsets: P: Positive, M: Medium, and N: Negative.                                    ∙ Switch S8 allows choosing between 0 and 𝑅1 .
                                                                                     If 𝐴𝑓 is positive then 𝑃𝑓 = R1 ; else, 𝑃𝑓 =0.
    The next step is the development of a droop controller rule base. By
                                                                                     These switches are controlled by control signals ‘‘ 𝐴𝑖𝑠𝑙 , 𝐴𝑣 and 𝐴𝑓 ’’
describing step by step the behavior of the process and the action of
                                                                                  generated by a fuzzy logic toolbox named ‘‘Fuzzy Logic Islanding
the variation in the command to apply, we deduce tables (basic fuzzy
                                                                                  Detection (FLID)’’. The input and output membership function of this
controller tables) which actually correspond to rules. The responses of
                                                                                  FLID is to ensure the transition between different operation modes, as
the frequency and voltage control are summarized in Tables 6 and 7.               provided in Fig. 12.
    This set of rules groups together all the possible situations of the             The inputs are grid voltage 𝑉𝑝𝑐𝑐 and grid frequency 𝑓𝑝𝑐𝑐 . The outputs
system evaluated by the different values attributed to ‘‘ 𝑒𝑖 ’’ and its           are the control signals of automatic switches ‘‘ 𝐴𝑖𝑠𝑙 , 𝐴𝑣 and 𝐴𝑓 ’’. We find
variation ‘‘ 𝑑𝑒𝑖 ’’ and all the corresponding values of the command               that the membership functions of the inputs have a symmetrical form.
variation, where 𝑖 = 𝑝, 𝑞.                                                        However, we introduce five fuzzy sets for each input: LB: Low Big, L:
    In this work, FLDC is used to exchange active and reactive powers             Low, M: Medium, H: High, and HB: High Big.
from the three-phase voltage source inverter to the grid by automatically            For the same reason, the forms of the membership functions of the
adjusting the amplitude and frequency of the output voltage in both               output variables are also symmetrical. Yet, we introduce two fuzzy sets:
standalone and grid-connected modes. Furthermore, it allows having a              N stands for negative and P stands for positive.
                                                                            223
Y. Krim et al.                                                                                                   Control Engineering Practice 81 (2018) 215–230
   The fuzzy logic rules used for the islanding detection are obtained            associate the fuzzy inputs to the fuzzy outputs are made up of three
from the analysis of the grid fluctuations. The control roles can help
AG to improve the performances of grid stability. The control rules that          parts, which are as follows:
                                                                            224
Y. Krim et al.                                                                                                     Control Engineering Practice 81 (2018) 215–230
Table 8                                                                              - In the grid connected mode, automatic switches S5 and S6 are respec-
𝐴𝑖𝑠𝑙 rules.                                                                                 tively in states 𝑉𝑝𝑐𝑐 and 𝜔𝑝𝑐𝑐 t, and S7 and S8 are respectively in
𝐴𝑖𝑠𝑙                        𝑉𝑝𝑐𝑐                                                            states 𝑅2 and 𝑅1 . The adaptive FLDC system reduces the voltage
                            HB         H          M            L          LB                and frequency fluctuations by adopting the power supplied by
                 HB         N          N          N            N          N                 the inverter to total powers ‘‘ 𝑃𝐿 +𝑃 𝑓 ’’ and ‘‘ 𝑄𝐿 +𝑄𝑣 ’’ consumed
                 H          N          P          P            P          N                 by the load and the grid. It controls exported or imported powers
 𝑓𝑝𝑐𝑐            M          N          P          P            P          N                 (𝑃𝑓 and 𝑄𝑣 ) with the grid ensuring its stability.
                 L          N          P          P            P          N
                 LB         N          N          N            N          N
                                                                                     - In the standalone mode, once FLID detects an isolating condition, RDG
                                                                                            will pass to the islanded mode and the additional injection of ‘‘
Table 9                                                                                     𝑃𝑓 and 𝑄𝑣 ’’ will have no sense, so they will be zero in this mode.
𝐴𝑣 rules.                                                                                   Automatic switches S7 and S8 are changed to state 0, and S5 and
                  𝑉𝑝𝑐𝑐                                                                      S6 are respectively changed to states 𝑉𝑛 and 𝜔𝑛 t.
                  HB               H         M             L              LB
 𝐴𝑖𝑣              N                P         N             P              N
                                                                                     6. Simulations and interpretations
                                                                               225
Y. Krim et al.                                                                                                                    Control Engineering Practice 81 (2018) 215–230
Fig. 13. Simulation conditions: (a) Wind speed, (b) Active power demanded by load, (c) Grid voltage, (d) Grid frequency.
Table 11                                                                                         demanded ones (Fig. 16). This simulation test is executed with the
Parameters of studied system.                                                                    classic control in MATLAB/SIMULINK and with FLS. Fig. 15(a) and
                        Parameters                                 Value                         15(b) represent the reference BT storage power deduced with classical
 PMSG                   Stator resistance 𝑅𝑠                       0.82 Ω                        control and fuzzy logic, respectively. Similarly for SC, Fig. 16(a) and
                        Stator inductance 𝐿𝑠                       0.0151 mH                     16(b) present the reference SC storage power deduced with classical
                        Magnet flux 𝛷𝑚                             0.4832 Wb
                                                                                                 control and fuzzy logic, respectively.
                        Nominal power                              3.5 kW
                        Number of pole pairs p                     4                                 The low rate of the charge of BT is 30% and the high rate of the
                        Inertia J                                  99*10−4 kg m2                 charge of BT is 90%. Also, the low rate of the charge of SC is 58 V and
                        Friction f                                 10−3 N m s rad−1              the high rate of the charge of SC is 98 V. Figs. 15(c) and 16(c) show the
 Wind turbine           Nominal power                              3.5 kW                        BT and SC currents, respectively. BT reacts more slowly to the needs.
                        Blade radius R                             2m
                                                                                                 Contrariwise, SC provides the transient currents. When SC reaches its
                        Optimal tip speed ratio 𝜆𝑜𝑝𝑡               8.15
                        Maximum power coefficient 𝐶𝑝𝑚𝑎𝑥            0.4794                        low rate of charge (𝑈𝑠𝑐𝑚𝑖𝑛 ), the totality of lack will be provided by BT.
                        Density of air 𝜌                           1.225 Kg m3                   Similarly, when BT reaches its low rate of charge, the totality of lack
 Battery                𝑆𝑂𝐶 𝑚𝑎𝑥                                    0.9                           will be provided by SC, in order to protect them from discharge excess.
                        𝑆𝑂𝐶 𝑚𝑖𝑛                                    0.3                           When SC reaches 98 V, its charge will be stopped and the production
                        Internal resistance 𝑅𝑖                     0.15 Ω
                                                                                                 excess must be absorbed by BT. Likewise, when BT reaches 90%, its
                        Nominal voltage 𝐸𝑏𝑎𝑡0                      60 V
                        Battery inductor filter 𝐿𝑏𝑎𝑡               10−6 H                        charge will be stopped and the production excess must be absorbed
 Supercapacitor         Capacitance 𝐶𝑠𝑐𝑐                           94 F+20%–0%                   by SC, in order to protect them from overcharge. Thus, the simulation
                        Rated voltage                              78 V                          results demonstrate the efficiency of proposed FLS by maintaining 𝑈𝑠𝑐
                        DC maximum current                         50 A                          (Fig. 15(d)) and 𝑆𝑂𝐶 𝑏𝑎𝑡 (Fig. 16(d)) at acceptable levels.
                        SC inductor filter 𝐿𝑠𝑐                     10−6 H
                                                                                                     By comparing the classical and fuzzy-logic simulation tests, we
                        Equivalent series resistance 𝑅𝑠𝑐𝑐          12.5 mΩ
                        [𝑈𝑠𝑐𝑚𝑖𝑛 , 𝑈𝑠𝑐𝑚𝑎𝑥 ]                         [58 V, 98 V]
                                                                                                 notice that EMS, using FLS, gives excellent results.
                        Leakage current                            0.15 A, 75 h, 25 ◦ C.             Consequently, the proposed fuzzy logic energy management supervi-
                        Operating temperature                      −40 ◦ C to +65 ◦ C            sor represents a reliable and efficient energy management. However, the
 DC bus                 Capacitance C                              2200 μF                       simulation results prove the effectiveness of the suggested strategy by
                        DC voltage                                 400 V
                                                                                                 ensuring the balance between production and consumption (Fig. 17(a))
 RLC filter             Filter resistance 𝑅𝑓                       0.2 Ω
                        Filter inductance 𝐿𝑓                       20 mH
                                                                                                 to keep a DC bus voltage stable at 400 V (Fig. 17(b)). In the proposed
                        Filter capacitance 𝐶𝑓                      10 μF                         simulations, BT and SC are able to respond to demands of the grid and
                                                                                                 the load. During the transitions of the wind power, the load and the grid,
                                                                                                 SC replies directly to the need of the load and the grid by providing or
                                                                                                 absorbing peak currents. BT reacts more slowly to the needs knowing
                                                                                                 that SC provides the transient currents (Fig. 17(c)).
                                                                                                     However, for the overproduction of a wind generator, and if the
                                                                                                 HESS reaches its high rate of charge, the wind turbine will operate in
                                                                                                 a limited mode to generate only the power quantity demanded by the
                                                                                                 load and the grid. When HESS reaches its low rate of charge and the
                                                                                                 wind power is unable to meet the demands, a no-priority load will be
                                                                                                 disconnected to ensure the balance between produced and consumed
                                                                                                 powers. Then the real, produced, stored and demanded powers are
                                                                                                 depicted in Fig. 17(a).
                                 Fig. 14. Stored power.
charge/discharge cycles (Fig. 15). Whereas, SC responds to the rapid                                According to the grid frequency and the voltage state, the active
changes in the difference between the produced wind power and the                                generator can operate in grid-connected and standalone modes (FLID).
                                                                                           226
Y. Krim et al.                                                                                                              Control Engineering Practice 81 (2018) 215–230
Fig. 15. BT storage system, (a) BT reference power with FLS, (b) BT reference power without FLS, (c) BT current, (d) BT SOC.
Fig. 16. SC storage system, (a) SC reference power with FLS, (b) SC reference power without FLS, (c) SC current, (d) SC SOC.
Fig. 17. (a) Simulated power assessment, (b) DC-bus voltage, (c) Simulation of 𝑃sc * and 𝑃bat * with variable demanded power.
Fig. 18(c) shows the periods in which AG operates in a separated                         to make them at more demanding margins than those mentioned in the
‘‘𝐴𝑖𝑠𝑙 = −0.5 ’’ or connected ‘‘ 𝐴𝑖𝑠𝑙 = 0.5 ’’ state. In a grid-connected
                                                                                         standards. These quantities are depicted in Fig. 18(a) and Fig. 18(b),
mode, the quantities of active and reactive powers are exchanged with
the grid to reduce the fluctuations in the voltage and the frequency, and                respectively.
                                                                                   227
Y. Krim et al.                                                                                                               Control Engineering Practice 81 (2018) 215–230
                                                                                               7. Conclusion
Fig. 18. (a) Active power exchanged with grid, (b) Reactive power exchanged with grid,
(a) Signal of operation mode.
                                                                                                   In this paper, an intelligent control strategy based on fuzzy logic
                                                                                               has been proposed, which facilitates the improvement of electrical
                                                                                               network stability by means of an active generator. The main purpose
    In this part, we present a comparative study to facilitate choosing                        of the intelligent control is to make an active generator participate
between GDC and FLDC.                                                                          in system services by adjusting the grid voltage and frequency and
    Fig. 19 presents the grid frequency before and after regulation with                       ensure a continuous supply load in case of grid failure. This control
fuzzy and generalized droop control. According to the latter, the impact                       has three control parts. The first one is the fuzzy logic supervisor,
of active power ‘‘ 𝑃𝑓 ’’ exchanged with the grid on the reduction in                           which manages the power flow to keep the DC bus voltage constant
the frequency fluctuation is remarkable. These variations are well mini-                       in order to deal with the load, the wind-power variations and the grid
mized, which demonstrates the importance of the renewable generator.                           fluctuations by maintaining BT and SC at admissible intervals of their
Similarly, with a grid voltage, the impact of reactive power ‘‘ 𝑄𝑣 ’’                          SOC. The second one is FLDC used to control the transferred active
exchanged with the grid on the reduction in the voltage fluctuation is                         and reactive powers with the grid ensuring its stability by ameliorating
remarkable (Fig. 20). With elaborate FLDC, it is remarkable that the                           the frequency and amplitude of its output voltage. The third part is
voltage and frequency levels are well minimized and become in optimal                          FLID, which ensures the transition between the grid-connected and
areas ([50.1, 49.9] and [235, 225], respectively).                                             separated modes and calculates the necessary quantity of active and
    The frequency and voltage errors are given in Figs. 21 and 22. We                          reactive powers to be injected or absorbed by the grid to guarantee its
find that with FLDC, the voltage and frequency levels are well minimized                       stability. Based on the comparative study between generalized and fuzzy
                                                                                         228
Y. Krim et al.                                                                                                                               Control Engineering Practice 81 (2018) 215–230
                  Table 12
                  Comparative study between proposed control approach and conventional control.
                    Comparison criteria         Grid frequency f 𝑝𝑐𝑐 stability                                      Grid voltage V 𝑝𝑐𝑐 stability
                                                At t = 45 s                        At t = 60 s                      At t = 45 s                     At t = 60 s
                                                GDC                FLDC            GDC               FLDC           GDC                FLDC         GDC               FLDC
                    Percentage of overrun       6%                 0.2%            5.8%              0.4%           5.21%              1.7%         10.43%            1.6%
                    Stabilization time (s)      0.17               0.04            0.2               0.06           0.16               0.035        0.172             0.06
                    Rising time (s)             0.2                0.11            0.2               0.15           0.2                0.11         0.2               0.15
                    Control capability                      Yes                                Yes                              Yes                             Yes
                    Control performance         High peak          Good            High peak         Good           High peak          Good         High peak         Good
                    Robustness                  Poor               Robust          Poor              Robust         Poor               Robust       Poor              Robust
logic control, the enhanced voltage and frequency regulation and the                                 Guerrero, J. M., Matas, J., de Vicuna, L. G., Berbel, N., & Sosa, J. (2006). Wireless-control
power management performance have been observed with the suggested                                      strategy for parallel operation of distributed-generation inverters. IEEE Transactions
                                                                                                        on Industrial Electronics, 53(5), 1461–1470.
fuzzy logic. The simulation results have demonstrated the validity of the
                                                                                                     Hou, G., Xing, F., Yang, Y., & Zhang, J. (2015). Virtual negative impedance droop method
proposed intelligent control strategy.
                                                                                                        for parallel inverters in microgrids. In Proceedings of the 10th conference on industrial
   We will start the concept and experimental implementation of the                                     electronics and applications (pp. 1009–1013). IEEE.
generalized control and the proposed control method for the studied                                  Il-Yop, C., Wenxin, L., David, A. C., Emmanuel, G. C., & Seung-Il, M. (2010). Control
renewable distributed generator, in order to make a real-time com-                                        methods of inverter-interfaced distributed generators in a microgrid system. IEEE
parative study between these two control techniques. In addition, in                                      Transactions on Industry Applications, 46, 1078–1088.
future work, we will develop an intelligent control based on fuzzy                                   Islam, Z., Djamel, B., Abdel-Moumen, D., Bilal, A., & Ikram. El, A. (2018). Hierarchical
logic by utilizing nonlinear and asymmetrical membership functions to                                    control for flexible microgrid based on three-phase voltage source inverters operated
                                                                                                         in parallel. Electrical Power and Energy Systems, 95, 188–201.
ensure more power grid stability. An optimization technique will be put
                                                                                                     Jia, H., Mu, Y., & Qi, Y. (2014). A statistical model to determine the capacity of battery–
forward to improve the choice of the membership functions by using
                                                                                                          supercapacitor hybrid energy storage system in autonomous microgrid. Electrical
neural networks or genetic algorithms for instance.                                                       Power and Energy Systems, 54, 516–524.
                                                                                                     Juan, C. V., Josep, M. G., Alvaro, L., Pedro, R., & Remus, T. (2009). Adaptive droop control
References                                                                                               applied to voltage-source inverters operating in grid-connected and islanded modes.
                                                                                                         IEEE Transactions on Industrial Electronics, 56, 4088–4096.
Arun, N. K., & Mohan, B. M. (2017). Modeling, stability analysis, and computational                  Krim, Y., Abbes, D., Krim, S., & Mimouni, M. F. (2017). Classical vector, first-order sliding
    aspects of some simplest nonlinear fuzzy two-term controllers derived via center of                  mode and high-order sliding-mode control for a grid-connected variable speed wind
    area/gravity defuzzification. ISA Transactions, 70, 16–29.                                           energy conversion system: a comparative study. Wind Engineering, 42(1), 16–37.
Cabrane, Z., Ouassaid, M., & Maaroufi, M. (2017). Battery and supercapacitor for photo-              Krim, Y., Abbes, D., Krim, S., & Mimouni, M. F. (2018). Power management and second-
    voltaic energy storage: a fuzzy logic management. IET Renewable Power Generation,                    order sliding mode control for standalone hybrid wind energy with battery energy
    11(8), 1157–1165.                                                                                    storage system. Proceedings of the Institution of Mechanical Engineers, Part I (Journal of
Edmund, C. C. C. (2004). Field measurement and experimental study of wind speed profile                  Systems and Control Engineering), 1–23.
    during thunderstorms. Journal of Wind Engineering and Industrial Aerodynamics, 92(3),
                                                                                                     Lee, C., Chuang, C., Chu, C., & Cheng, P. (2009). Control strategies for distributed energy
    275–290.
                                                                                                         resources interface converters in the low voltage Microgrid. In 2009 IEEE energy
Erdinc, O., & Uzunoglu, M. (2011). The importance of detailed data utilization on
                                                                                                         conversion congress and exposition, (p. 20-22–9).
    the performance evaluation of a grid-independent hybrid renewable energy system.
    International Journal of Hydrogen Energy, 36, 12664–12677.                                       Masmoudi, A., Abdelkafi, A., & Krichen, L. (2011). Electric power generation based on
Furtado, E. C., Aguirre, L. A., & Torres, L. A. B. (2008). UPS parallel balanced operation              variable speed wind turbine under load disturbance. Energy, 36, 5016–5026.
    without explicit estimation of reactive power-a simpler scheme. IEEE Transactions on             Matas, J., Castilla, M., Vicuña, L. G. D., Miret, J., & Vasquez, J. C. (2010). Virtual
    Circuits and Systems II: Express Briefs, 55(10), 1061–1065.                                         impedance loop for droop-controlled single-phase parallel inverters using a second-
Glasnovic, Z., & Margeta, J. (2011). Vision of total renewable electricity scenario.                    order general-integrator scheme. IEEE Transactions on Power Electronics, 25, 2993–
    Renewable Sustainable Energy Review, 15(4), 1873–1884.                                              3002.
                                                                                               229
Y. Krim et al.                                                                                                                             Control Engineering Practice 81 (2018) 215–230
Matthieu, D., Arnaud, D., Karl, S., Moe, T., Marc, M., & Richard, R. (2017). Battery energy           Standard EN 50160. (2004). Voltage characteristics in public distribution systems.
    storage system battery durability and reliability under electric utility grid operations:         Tayaba, U. B., Roslana, M. A. B., Hwaia, L. J., & Kashifa, M. (2017). A review of droop
    analysis of 3 years of real usage. Journal of Power Sources, 338, 65–73.                              control techniques for microgrida review of droop control techniques for microgrid.
Mouna, R., Achraf, A., & Lotfi, K. (2015). A micro-grid ensuring multi-objective control                  Renewable and Sustainable Energy Reviews, 76, 717–727.
    strategy of a power electrical system for quality improvement. Energy, 88, 351–363.               Tolani, S., & Sensarma, P. (2012). An improved droop controller for parallel operation of
Nikita, G., & Garg, Rachna (2017). Tuning of asymmetrical fuzzy logic control algorithm                   single-phase inverters using R-C output impedance. In Proceedings of the 2012 IEEE
    for SPV system connected to grid. International Journal of Hydrocarbon Engineering, 42,               international conference on power electronics, drives and energy systems (pp. 1–6).
    16375–16385.                                                                                      Vasquez, J. C., Guerrero, J. M., Luna, A., Rodriguez, P., & Teodorescu, R. (2009). Adaptive
Petronela, P., Dhaker, A., Christophe, S., Stephane, B., Julien, P., & Benoit, R. (2016).                 droop control applied to voltage-source inverters operating in grid-connected and
    Multi-criteria fuzzy-logic optimized supervision for hybrid railway power substations.                islanded modes. IEEE Transactions on Industrial Electronics, 56(10), 4088–4096.
    Mathematics and Computers in Simulation, 130, 236–250.                                            Victor, H., Aitor, M., Haizea, G., Ion, E-O., Igor, V., & Haritza, C. (2016). Adaptive energy
Randa, K., Ghada, B., & Lotfi, K. (2015). Demand side management of household                             management strategy and optimal sizing applied on a battery-supercapacitor based
    appliances in stand-alone hybrid photovoltaic system. Renewable Energy, 81, 123–135.                  tramway. Applied Energy, 169, 831–845.
Sahyoun, S., Djouadi, S. M., & Shankar, M. (2015). Optimal control of droop controlled                Wang, Y., Liu, C., Pan, R., & Chen, Z. (2017). Modeling and state-of-charge prediction
    inverters in islanded microgrids. IFAC-Papers, 48, 363–368.                                           of lithium-ion battery and ultracapacitor hybrids with a co-estimator. Energy, 121,
Saleh, A., Shoresh, S., & Hassan, B. (2015). A fuzzy logic-based droop control for                        739–750.
    simultaneous voltage and frequency regulation in an AC microgrid. Electrical Power                Youssef, K., Saber, K., & Mohamed, F. M. (2017). Control and management strategy
    and Energy Systems, 64, 148–155.                                                                      of a generation system made by variable-speed wind turbine and super - capacitor
Shuai, Z., Shanglin, M., Jun, W., John, S., Wei, T., & Yan, F. (2016). Droop control method               energy storage system for standalone Applications. International Journal on Electrical
    for load share and voltage regulation in high voltage microgrids. Journal of Modern                   Engineering and Informatics, 9(3), 567–584.
    Power System Clean Energy, 4, 76–86.                                                              Ziyou, S., Jun, H., Heath, H., Jianqiu, L., & Minggao, O. (2017). Sliding-mode and
Song, Z., Hofmann, H., Li, J., Hou, J., Han, X., & Ouyang, M. (2014). Energy management                   lyapunov function-based control for battery/supercapacitor hybrid energy storage
    strategies comparison for electric vehicles with hybrid energy storage system. Applied                system used in electric vehicles. Energy, 122, 601–612.
    Energy, 134, 321–331.
230