This article has been accepted for publication in a future issue of this journal, but has not been

fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 1

Load Frequency Control in Microgrids Based on a

Stochastic Non-Integer Controller

Mohammad-Hassan Khooban, Member, IEEE, Taher Niknam, Member, IEEE, Mokhtar Shasadeghi,
Tomislav Dragicevic, Senior Member, IEEE, and Frede Blaabjerg, Fellow, IEEE

 is quite significant as it can decrease power loss in

Abstract—In this paper, an adaptive multi-objective transmission lines and increase the overall reliability of the
Fractional-Order Fuzzy proportional–integral–derivative electrical power systems [4]. The concept of the microgrid
(MOFOFPID) controller is proposed for the load frequency (MG) has been proposed to locally co-optimize the operation
control (LFC) of islande d Microgrids (MGs), while benefiting
from the assets of electric vehicles (EVs) in this respect. Although
of a number of DGs and to further enhance the controllability,
the use of battery energy storage systems (BESS) can solve the reliability, flexibility and power quality of the overall
unbalance effects between the load and supply of an isolated MG, electrical systems, while postponing the construction of new
their high cost and tendency toward degradation are restrictive power stations at the same time [5].
factors, which call for the use of alternative power balancing From the utility point of view, MG is an independent unit
options. In recent years, the concept of utilizing the BESSs of which can be operated in both grid-connected and isolated
EVs, also known as vehicle-to-grid (V2G) concept, for frequency
support of MGs has attracted a lot of attention. In order to allow
modes [6]. In the grid connected mode, the main target is
the V2G controller operate optimally under a wide range of usually to perform auxiliary services to the overhead system
operation conditions caused by the intermittent behavior of and optimize the internal energy management at the same time
renewable energy resources (RESs), a new multi-objective [7]. On the other hand, the isolated operation of MG can be
fractional-order control strategy for the EVs in V2G scenarios is used for supplying remote applications (telecom, rural
proposed in this paper. Moreover, since the performance of the households, etc.) or for enhancing the power quality of weak
controller depends on its parameters, optimization of these
grids by manually disconnecting it from the overhead system.
parameters can play a significant role in promoting the output
performance of the LFC control; hence, a modified black hole In the latter mode, due to the low inertia of the MGs compared
optimization algorithm (MBHA) is utilized for the adaptive with the grid, it is a challenge to maintain the active and
tuning of the non-integer fuzzy PID controller coefficients. The reactive power balances between the supply and consumption,
performance of the proposed LFC is evaluated by using real especially in the presence of intermittent RES and frequent
world wind and solar radiation data. Finally, the extensive load variations [8, 9]. The most common solution to this
studies and hardware-in-the-loop (HIL) simulations are
problem in isolated MGs was the installation of backup
presented to prove that the proposed controller tracks frequency
with lower deviation and fluctuation and is more robust in systems such as micro turbines and energy storage systems
comparison with the prior-art controllers used in all the case (ESSs).
studies. Besides RES, electric vehicles (EVs) and plug-in hybrid
electric vehicle (PHEV) technologies are other key
Index Terms— Load Frequency Control (LFC), Modified technologies that will lead to a carbon-neutral society [10].
Black Hole Algorithm (MBHA), MicroGrid (MG), Fractional Since EVs are environmentally friendly and can significantly
Controller, Electric-Vehicle (EV).
reduce the use of gasoline, it is predicted that they will
I. INTRODUCTION dominate the future of automotive industry. In the PHEV
technology, the bidirectional flow of power between the grid
I n recent years, due to the decline of conventional fossil fuel
resources and the environmental pollutions they caused,
modern electrical power systems face increased deployment of
and the vehicle’s battery can be achieved, so that its remaining
capacity can be used for ancillary services such as peak
shaving or frequency support. However, the charging and
renewable energy sources (RESs) like photovoltaic (PV) and discharging of PHEV batteries can increase the complexity of
wind turbine generator (WTG) systems [1-3]. In this context, the overall frequency control of the system [11, 12].
the local harness of RES through distributed generators (DGs) When considering the autonomous MG, both RES and EVs
need to be coordinated together to limit the frequency
M.H. Khooban, T. Niknam and Mokhtar ShaSadeghi are with the fluctuation by compensating for the mismatch between
Department of Electrical and Electronic Engineering, Shiraz University of generation and demand. This functionality is often referred to
Technology, Shiraz, Iran (e-mail: khooban@sutech.ac.ir, and as load frequency control (LFC). The proper LFC can not
niknam@sutech.ac.ir).
T. Dragicevic and F. Blaabjerg are with the Electrical Engineering only guarantee the frequency stability of the MG but also
Department, University of Aalborg, Aalborg, Denmark, (e-mail: increase its efficiency (fuel saving). In order to improve the
fbj@et.aau.dk).

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 2

response of LFC, many controllers including conventional II. MODELING OF AN ISOLATED MICROGRID
PID control [13], intelligent control [14], adaptive control
A. The model of an isolated microgrid
[12], robust control [15] and MPC control [11] have been
applied to the DGs of islanded MGs. In [16], an MPC-based Fig. 1 depicts an isolated MG in which DGs such as PVs and
coordinated control of the blade pitch angles of the WTG and WTGs, and energy storage units such as BESS and FESS
supply the distributed loads [22].
PHEV has been proposed for the LFC. The MPC controller is
The power grid and the MG operation are controlled by the
applied in order to smooth wind power production and
distribution management system (DMS) and the MG dispatch
decrease the number of required PHEVs. In [17], a new system (MGDS), respectively. Also, the bidirectional
method is suggested for the LFC in MGs by using an information transfer can be achieved through communication
intelligent PI controller for improving the robustness of the links [11].
whole of system. Moreover, the integral square error (ISE) is
used for optimum tuning of PI’s parameters to increase the
performance of the proposed controller [17].
Since the operation conditions of the LFC can widely change,
the PI controller which is tuned for nominal conditions cannot
work properly in other conditions. So, in order to solve this
problem, the fuzzy logic (FL) which adjusts the control
parameters according to the operating point is proposed in
[11]. One drawback of the system is that its good performance
is achieved in only some specific member functions. In [18],
the application of robust H-infinity control in the LFC of an Fig. 1.General scheme of Micro-Grids.
isolated MG is studied. The control method, which has been B. Electric Vehicle Model
suggested in [18], is very complicated and hence not attractive Generally, there are different numbers of EVs in each EV
for real world applications. As a result, the small signal station. As a result, using an equivalent EV with a different
analysis of an isolated MG in the presence of energy storage inverter capacity can handle the modeling of the EVs'
unit has been proposed in [19]. Besides, in order to increase behavior. The equivalent EV model, proposed in [11, 23] for
the robustness of the isolated MG, the hierarchical control is the LFC study, is shown in Fig. 2. In this figure, Te and
presented in [20]. represent the time constant and the LFC command signal to
In order to overcome the aforementioned difficulties, this the EV, respectively. The interfacing inverter capacity limits
paper introduces a new adaptive approach by using a novel are shown by and which represent the ramp rate
optimization technique for the adaptive tuning of the most limits. The lower and upper limits of the EV’s controllable
common fractional-order fuzzy PID (FOFPID) controller for energy are represented by and , correspondingly.
the EV that supports the LFC in islanded MG(s). FOFPID The charging or discharging power of the EV is represented
parameters are tuned automatically according to the online by . The positive sign of indicates that the EV is in the
measurements, by applying MBHA. Unlike the classical discharging mode while the negative sign shows that EV is in
tuning methods which are not suitable for providing a useful the charging mode. expresses the idle operation mode
performance over a wide range of operation conditions, many of the EV.
advantages are offered by the proposed optimal tuning scheme
for the frequency control of an MG with many DGs and RESs.
Moreover, the proposed method is significantly less complex
in comparison to the above-mentioned approaches, making it
more convenient for practical applications. The simulation
study is performed on a complex MG, including different
loads and RESs to demonstrate the effectiveness of the
proposed control scheme; the superiority of the suggested Fig. 2.The electric vehicle model in Microgrid LFC.
controller over Multi-Objective Proportional–Integral–
Derivative (MOPID), Multi-Objective Fuzzy Proportional–
C. The Diesel power system model
Integral (MOFPI) and Multi-Objective Interval Type-2 Fuzzy
Logic (MOIT2FLC) controllers is demonstrated in Section VI. Due to the advantages of the diesel power system, including
To further corroborate the performance and robustness of the its fast starting speed, low maintenance and high efficiency, it
proposed control system, experimental validation using the has been a good backup option in isolated microgrids. The
Hardware-in-the-Loop (HiL) simulation is also given in this controllable DG can instantaneously track load demand
paper. changes through the power control mechanisms it employs
[17]. The fluctuation of uncontrollable DGs such as WTG, PV
and loads can be effectively compensated for by the diesel
power system.

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 3

However, power loss is inevitable in them due to the boundary

and external contact with the series resistance. Moreover, the
generated power of the PV is intermittent and depends on sun
irradiance and temperature; hence, in simulations a random
power source can be used to model the behavior of PV [13,
22]. The PV model which is considered as disturbance for the
LFC is depicted in Fig. 4.
Fig. 3.Diesel Power System Model

F. General scheme of microgrid with LFC controller

Fig. 3 shows the transfer function of the diesel power system
which represents the relationship between the LFC signal and The framework of the proposed LFC and an isolated MG
the output power of DG. As shown in Fig. 3, the model which consists of EVs, DGs such as PV, WT and diesel
generator, energy storage units like BESS and FESS and loads
consists of first order inertia models of the governor and diesel
is depicted in Fig. 5. As shown in this figure, the PV, Fuel Cell
generator and two saturation blocks.
(FC), BESS and FESS units are connected to the AC MG via
In this figure, and represent frequency deviation and DC/AC interfacing inverters. The FC is modeled by a third
the LFC command signal from the DG, respectively. and order transfer function [22, 24]. All small scale DGs and
denote the time constant of the governor and the diesel energy storage units are connected to the AC bus via a circuit
generator, correspondingly. shows the condition of the breaker. The spinning reserve for the secondary frequency
governor’s valve. The speed regulation coefficient of the DG control is provided by the diesel power system.
is shown by R in this Fig. and represent the power
increment and ramp rate limits, respectively. The output
power increment of the diesel power system is represented
by . means that demand and generation are in a
balance condition and there is not any need for changing the
power. means that the required power is higher than
the actual power while the represents the condition
that the actual power is less than the demand [23].

Fig. 5.Case Study model for the LFC in MG system.

The parameters of the MG system which is depicted in Fig. 4

are listed in Table.1.

Table 1: MG power system’s parameters

S ymbol and Values S ymbol and Abbreviation Values
Abbreviation
2.000 s 0.5 s
1.000 s 1.8 s
3.000 2s
0.010 1/300
0.025 1
0.950 1
0.800 1/500
D 0.012 0.6

2H 0.200 -1/100

0.100 s -1/300
Fig.4. The overall microgrid scheme including EV for LFC
0.100 s 1/100
D. Model of Wind Turbine 4.000 s ‘s’ means second
Because of the intermittent behavior of wind, the output power
of wind turbines changes according to the different speed and
directions wind takes [13, 22]. By controlling the DG units III. THE STRUCTURE OF FRACTIONAL ORDER FUZZY PID
CONTROLLER
and EVs, the uncontrolled WTG can be modeled as a
disturbance and a power fluctuation source in LFC. The A. Fractional-order calculus (FOC) in Control Systems
dynamic model of WTG is shown in Fig. 4.
Fractional calculus is one of the most important branches of
E. Photovoltaic (PV) Generation calculus in which the power of the differential and integration
PV cells which are made from semiconductor materials can operators can take a non–integer value. During the last few
directly convert the energy of photons to electrical energy. decades, the FOC has been applied in many fields of

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 4

automatic control systems [25]. Fractional Calculus-Based Moreover, the former parts of each rule will be made by using
Control Systems can be illustrated by where AND function (with hermeneutics of minimum). Besides, the
is the order of operation while t , are the upper and Mamdani fuzzy inference engine is also applied in this study
lower limits, respectively [26]. There exist a large number of [29]. In fuzzy control systems, changing the scaling factors of
definitions, such as Grünwald–Letnikov (GL), Cauchy integral the fuzzy output has a more positive influence on the
formula and Riemann–Liouville (RL), which are applied to performance of the FOFPID controller than does making
define the FOC. However, in automatic control systems, the alterations in the form of the MF(s). As a result, the tuning of
Cauchy integral formula is frequently applied for the purpose all coefficients in the FOFPID cannot be equally effective in
of realizing the fractional-order differentiations and affecting the total performance of the control system. The goal
integrations of the FO-Fuzzy-PID (FOFPID) controller [27]. of this study is to explore the result of tuning fractional-orders
∫ , ( and since keeping the rules set and the form of MFs
unaltered can noticeably increase the total closed loop
(1) performance of the FOFPID controller. Consequently, this
B. Design of fractional order proportional integral- method can be more effective in reducing the volume of
derivative (FOPID) controller calculations, and as a result, the output of the FOFPID
controller can be extremely influenced by tuning fractional-
The main advantage of the FOFPID controller which has been orders instead of tuning the MF(s) variables or other fuzzy
used in this study is that it encompasses the benefits of the inferencing parameters. So, the next section presents a new
combination of fuzzy logic with PID controller which is modified multi-objective optimization algorithm to turn the
comprehensively discussed in [25]. In the proposed method, parameters of proposed controller.
and are considered as the fuzzy logic input and, and
as fuzzy logic outputs. The advantages of this kind of
FOFPID structure over the Model Predictive Control (MPC) IV. OVERVIEW OF THE ORIGINAL BLACK HOLE ALGORITHM
and the conventional PID are presented in [28]. In the original (BHA)
FPID controller, the power of the input error derivative is
integer. In this study, in contrast, this power is fractional-order The idea of BHA was taken from the intelligent collective
( ). Moreover, the order of the integral in the output is behavior of the stars around a black hole. The movement
characteristic of the stars can be expressed as in [30] and [31]:
changed by the fractional order counterpart ( ). The general
scheme of the proposed method is shown in Fig. 6.
(3)
where and represent the position of the target
and the updated agent in iteration iter, respectively. A
comprehensive introduction with details of this algorithm is
provided in [31].
A. IBH Mechanism [31]:
A novel method is suggested for the collapsing process so that
the exploration of the properties of the BHA will be enhanced.
Fig. 6.The general scheme of Fractional-order PID (FOPID) Primarily, the updating mechanism of the BHA will be
controller. enhanced through (4):
By considering the general scheme of Fig. 6, the control law
can be written as: (4)
where r can be regarded as a random number within the range
(2) of [1, NPop] and . In the second step, the modification
The other conventional controllers such as the integer PID and of the event horizon of the BH ( ) based on the distribution
integer Fuzzy PID can be obtained by considering the values and collection of the stars should be done using (5) and (6).
of and as one. Table 2 displays the set of optimal FOFPID
rules. The fuzzy linguistic variables which are similar to the ‖ ‖ (5)
input and output variables are classified as (NS), (NM), (NL), ∑ (6)
(ZR), (PS), (PM), and (PL), referring correspondingly to
Negative Small, Negative Medium, Negative Large, Zero, Moreover, the new concept of Absorption Capacity (AC) is
Positive Small, Positive Medium and Positive Large. The proposed for the BH. The elimination of the hyper dispersion
categorization is based on triangular MF(s). of the solution and the control of the number of the stars
Table 2: The MOFOFPID controller Rules Set located in event horizon is achieved in this method. In BHA,
NL NM NS PS PM PL
after a star collapses in the BH, a new star is born randomly.
̇ The optimum usage of the data obtained by the members of
S NL NM NS PS PS PM the population is achieved by the new formulation expressed
M NL NL NM PS PM PM in (7).
L NL NL NL PM PM PM

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 5

minimized in the Interactive Fuzzy Satisfying Method can be

(7)
proposed in terms of the following equation:
where
V. MULTI-OBJECTIVE INTERACTIVE FUZZY SATISFYING (IFS) ∫ (11)
METHOD
∫ (12)
The value and the number of control signals are among the
In this equation, the first objective function is used for
most important factors to consider in designing and
providing the fast tracking of the desired set-point while the
implementing hardware. In other words, in an application with
accuracy of tracking the set point is achieved by the second
large values of control signals, a big size actuator is required,
objective function. A schematic representation of the online
which correspondingly increases the total cost of the system.
tuning of the FOFPID controller based on the MOMBHA
Consequently, effective control of a system requires a smooth
technique is depicted in Fig. 7.
control that can track the set point fast. Using the multi-
objective optimization principle can overcome complexity in
optimal control problem. Weighting the cost function of the
optimum problem with some indexes such as the Squared
Deviation of Controller Output (ISDCO), Integral of Absolute
Error (IAE), Integral of Time multiplied Squared Error (ITSE)
and the Integral of Time Multiplied Absolute Error (ITAE)
can prevent the high control signal. It should be mentioned
that the aforementioned objectives (IAE, ISDCO, ITSE and
ITAE) can act in contrast with each other. Consequently,
formulating the control system as a multi-objective stochastic
optimization problem can reduce the value of the control
effort.
Using the MG system can inherently emulate some of the
Fig.7.The general scheme of the suggested controller for LFC
uncertainties of the system pertaining, for instance to its
output variations, and storage batteries; hence, using integer
VI. SIMULATION RESULTS
and non-integer controllers can be a good choice for a system
in which the operation point is widely subject to change. To In this section, in a trial run to examine the proposed control
solve the optimization problem, in this paper a new online method, the MG which is depicted in Fig. 4 is simulated in
multi-objective optimization method for the optimum tuning MATLAB/Simulink software. The parameters of the isolated
of the common FOFPID controller is proposed. In this microgrid are listed in Table 1. The performance of the
method, using online measurement, the Multi-objective proposed controller is compared with that of the MOPID,
MBHA (MOMBHA) is applied for the optimum tuning of the MOFPI and MOIT2FLC controllers. Since the respond of the
FOFPID controller. The MOMBHA problem can be system depends on the parameters of these controllers, all of
formulated by (8). the parameters are optimized by the MOMBHA optimization
algorithm in different controllers. In order to evaluate the
performance of the proposed control method in the context of
(8) the MG depicted in Fig. 4, the Hardware-In-the Loop (HIL)
simulation approach is utilized. The real time HIL method is
In this paper, an IFS approach is used to resolve the used to emulate errors and delays that do not exist in the
aforementioned multi-objective problem [32]. In this method, classical off-line simulations. Fig. 8 sketches the HIL setup,
the sufficient degree of each objective can be defined by the consisting of: 1) OPAL-RT as a Real Time Simulator (RTS)
operator. Generally speaking, the non-inferior solution set in which simulates the MG depicted in Fig. 4; 2) a PC as the
the most optimal solution would be chosen in such a way that command station (programming host) in which the
the preferences of the operators are satisfied: Matlab/Simulink based code that is executed on the OPAL-RT
{ | |} (9) is generated; and 3) a router used as a connector of all the
setup devices in the same sub-network. The OPAL-RT is also
Here, μref is the satisfying degree of the related target,
connected to the DK60 board through Ethernet ports. More
determined by the operator in the range of [0, 1]. , on the details about the components of this setup can be found in [33,
other hand, denotes the MF value of the wth objective and is 34].
designed by the trapezoidal fuzzy membership as represented
below:

{ (10)

It should be noted that a heuristic algorithm like the MBHA

merely requires information about the fitness function. As a
consequence, in order for the online tuning of the FOFPID
controller, the objective functions vector that should be Fig.8.The Real time experimental setup.

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 6

A. Case 1: simulation. As illustrated in Fig. 12, the performance of the

In the first step, it is assumed that the load demand in the LFC in eliminating the load disturbance can be improved
isolated MG is constant i.e., . In other words, the when using the proposed controller; moreover, as shown in
power fluctuation of WTG ( ) and PVG ( ) is this figure, the settling time of the frequency response is
consistent in the LFC system. The wind power fluctuation data decreased significantly.
extracted from an offshore wind farm in Sweden [35] is
depicted in Fig.9.a, while Fig.9.b shows the solar radiation
data in Aberdeen (United Kingdom) [36], which was also used
in this case study. Fig.10 depicts the frequency response of the
simulated system.
Fig. 11. Step changes of the load in the interval of 90 seconds .

Fig. 12.Frequency response according to the power fluctuation of WPG,

Fig. 9.Power fluctuation, a- WPG, b- PV. PV and load disturbances.

C. Case 3
To evaluate the robustness of the proposed controller, some
parameters of the isolated MG are changed in the third case
study. The changing of the parameters is made in a scenario
which is displayed in Table 3.
Table 3. Uncertain parameters of the MG system
Parameters Variation Range
R +15%
D -25%
H +30%
Fig. 10. Frequency response according to the power fluctuation of WPG and -25%
PV +30%
As shown in Fig. 10, employing MOFOFPID reduces the peak -20%
value of frequency deviation and this reduction is achieved in +30%
the least amount of time as compared to the time when the As shown in this table, a more severe changing of the
MOPID, MOFPI and MOIT2FLC controllers are used. As a parameters is applied to evaluate the robustness of the
result, MOFOFPID can considerably enhance the life of both proposed control method. Fig. 13 depicts the response of the
the batteries and the equipment used in the DG by providing a LFC in the scenario by using the proposed as well as the
faster stabilization in the output frequency of the MG while MOPID, MOFPI and MOIT2FLC controllers.
requiring less frequency adjustment. .
Moreover, the simulation results in this case suggest that using
the MOFOFPID results in the LFC tracking the reference
frequency with fewer overshoots and much smaller settling
time in comparison to the other three controllers.

B. Case 2
In this section, multi-step load variation is applied to the LFC
as a disturbance. The load steps are depicted in Fig.11, while
Fig. 13.The frequency deviation of the micro-grid according to Scenario 3.
the frequency deviation responds of MOPID, MOFPID,
MOIT2FLC and MOFOFPID controls are depicted in Fig.12.
As revealed in Fig. 13, the proposed controller enhances the
As displayed in Fig. 12, the frequency deviation overshoot is
performance of the LFC more effectively than the other three
decreased and the LFC control can eliminate the effect of load
control methods, especially when the overshoots are
disturbance more effectively when the proposed method,
concerned. In other words, the results suggest that the
rather than the others, is employed. The performances of these
proposed controller is more robust against the changes in
controllers are further examined by creating a severe
parameters as compared to the other controllers. It is also
condition, applying a large load step at t=60s of the

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 7

revealed that the MOPID, MOFPI and MOIT2FLC controllers VII. CONCLUSION
do not have a satisfactory performance in this scenario where
there are severe parameter changes. In this paper, a new adaptive and time-varying controller was
presented for the load frequency control of an isolated
D. Case 4 microgrid. A stochastic multi-objective optimization algorithm
Finally, the last scenario examines the efficiency of the was used to optimize the parameter of the controller to track
proposed controller in a situation where there are limitations in the reference frequency in the presence of PV, wind
using EV(s) in the MG. To this end, the EV is disconnected generation, V2G electrical vehicles and load disturbance. To
from the grid at t=47s. At the same time, the same load enhance the robustness of the controller against the stated
disturbance and wind/PV power fluctuations considered in uncertainties of the microgrid, it was designed in two levels,
Cases A and B are applied to the system. The simulation namely a Fractional-Order Fuzzy and a conventional PI
results are presented in Figure 14. In addition, three common controller. Since the performance of the fuzzy systems
error measurement criteria are used to evaluate the efficacy of depends on their membership functions, the membership
the MOPI, MOFPI, MOIT2FLC and the suggested control function parameters were optimized by using
scheme in this scenario. These included Sum of the Squared modification/improvement of MOBHA optimization
Errors (SSE), Mean Absolute Error (MAE), and Mean Square algorithm. This novel approach enhances the performance of
Error (MSE). These approaches can lead to an optimal the LFC with low computation burden and complexity. In this
performance if the values of SSE, MAE, and MSE are close to study, both load disturbance and the output power of DG units
zero. Table 4 presents the evaluation results obtained for these were considered as ΔP in the model of microgrid; hence, the
controllers. proposed control approach was found to be adaptive enough to
give an appropriate performance in face of uncertain loads,
renewable energy sources and MGs of any typology. In order
to validate the performance and robustness of the controller,
Hardware-in-the-Loop (HiL) simulation was used in this
study. Furthermore, the performance of the controller was
compared with MOPID, MOFPI and MOIT2FLC controllers
which have been suggested in the most recent research studies
on LFC.

VIII. ACKNOWLEDGMENTS
Fig. 14.The frequency deviation of the micro-grid according to Scenario 4. Appreciation is offered to the Iran National Science
Foundation for a grant to collect data for this research. The
The frequency deviations presented in figure 14 demonstrate authors wish to thank the Iran National Science Foundation for
that the MOFOFPID controller can better improve the supporting helps to do this research and providing this
performance of the MG in comparison to the other three opportunity for the survey team. Moreover, the authors would
control structures especially as the overshoots are concerned. like to thank the reviewers for their very useful comments,
Indeed, as it can be seen, despite the EV battery restriction that which helped us to improve the quality of the paper. The
occurs in the MG in this scenario, the performance of the reviewers’ efforts are gratefully appreciated.
suggested controller is superior to that of the other three that
appear in figure 14. Further noticeable in this figure are the References
results pertaining to the robustness of the proposed controller 1. Bose, B.K., "Global Warming: Energy, Environmental Pollution,
against both parameter changes and the restriction that was and the Impact of Power Electronics," IEEE Industrial Electronics
imposed on the system by disconnecting the EV. As it can be Magazine, vol.4, no.1, pp.6-17, March 2010.
2. Dehghani M, Khooban MH, Niknam T, Rafiei SM. Time-Varying
noticed, none of the MOPID, MOFPI and MOIT2FLC Sliding Mode Control Strategy for Multibus Low-Voltage
controllers has such an agreeable performance so far as their Microgrids with Parallel Connected Renewable Power Sources in
frequency response against the severe changes of parameters Islanding Mode. Journal of Energy Engineering. 2016 Feb
(Scenario 3) and their control of the MG system are 29;142(4):05016002.
3. T.C. Green, M. Prodanović, "Control of inverter-based micro-
concerned. It is worth mentioning here that unlike the previous grids," Electric Power Systems Research Vol.77 pp.1204–1213,
research in this area, the current investigation was carried out 2007.
longitudinally to further contribute to the reliability of the 4. Chowdhury, S. and P. Crossley, "Microgrids and active
results concerning the satisfactory functioning of the proposed distribution networks". The Institution of Engineering and
Technology (IET), 2009:
controller. 5. Guerrero, J.M.; Chandorkar, M.; Lee, T.; Loh, P.C., "Advanced
Control Architectures for Intelligent Microgrids Part I:
Table 4.The performance analysis of controllers Decentralized and Hierarchical Control," IEEE Transactions on
Industrial Electronics, vol.60, no.4, pp.1254,1262, April 2013.
Criteria M OPID M OFPI M OIT2FLC M OFFOPID 6. H. Jiayi, J. Chuanwen, and X. Rong, “A review on distributed
SSE 10.8713 8.9211 1.8758 0.6919
energy resources and MicroGrid,” Renew. Sustain. Energy Rev.,
M AE 0.0799 0.0327 0.0055 0.0013 vol. 12, no. 9,pp. 2472–2483, Dec. 2008.
M SE 7.6311e-4 2.3527e-4 5.9094e-5 4.9006e-5 7. A. Hooshmand, B. Asghari and R. K. Sharma, "Experimental
Demonstration of a Tiered Power Management System for
Economic Operation of Grid-Tied Microgrids," in IEEE

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 8

Transactions on Sustainable Energy, vol. 5, no. 4, pp. 1319-1327, 28. Pan I, Das S. Enhancement of Fuzzy PID Controller with
Oct. 2014. Fractional Calculus. InIntelligent Fractional Order Systems and
8. R. H. Lasseter and P. Paigi, "Microgrid: a conceptual solution," Control 2013 (pp. 159-193). Springer Berlin Heidelberg.
Power Electronics Specialists Conference, 2004. PESC 04. 2004 29. Khooban MH, Niknam T. A new intelligent online fuzzy tuning
IEEE 35th Annual, 2004, pp. 4285-4290 Vol.6. approach for multi-area load frequency control: Self Adaptive
9. Yang, Jun, ZhiliZeng, Yufei Tang, Jun Yan, Haibo He, and Modified Bat Algorithm. International Journal of Electrical Power
Yunliang Wu. "Load Frequency Control in Isolated Micro-Grids & Energy Systems. 2015 Oct 31;71:254-61.
with Electrical Vehicles Based on Multivariable Generalized 30. Hatamlou A. Black hole: A new heuristic optimization approach
Predictive Theory." Energies 8, no. 3 (2015): 2145-2164. for data clustering. Information sciences. 2013 Feb 10;222:175-84.
10. Khooban MH, Niknam T, Blaabjerg F, Dehghani M. Free 31. Azizipanah-Abarghooee R, Terzija V, Golestaneh F, Roosta A.
chattering hybrid sliding mode control for a class of non-linear Multiobjective dynamic optimal power flow considering fuzzy -
systems: electric vehicles as a case study. IET Science, based smart utilization of mobile electric vehicles. IEEE
Measurement & Technology. 2016 Oct 1;10(7):776-85. Transactions on Industrial Informatics. 2016 Apr;12(2):503-14.
11. Yang, Jun, ZhiliZeng, Yufei Tang, Jun Yan, Haibo He, and 32. Khooban MH, Sha-Sadeghi M, Niknam T, Blaabjerg F. Analysis,
Yunliang Wu. "Load Frequency Control in Isolated Micro -Grids Control and Design of Speed Control of Electric Vehicles Delayed
with Electrical Vehicles Based on Multivariable Generalized Model: Multi-Objective Fuzzy Fractional-Order PIλDμ Controller.
Predictive Theory." Energies 8, no. 3 (2015): 2145-2164. IET Science, Measurement & Technology. 2016 Nov 18.DOI:
12. Zheng, S., Tang, X., Song, B., Lu, S. and Ye, B., 2013. Stable 10.1049/iet-smt.2016.0277.
adaptive PI control for permanent magnet synchronous motor drive 33. Zhang H, Zhang Y, Yin C. Hardware-in-the-Loop Simulation of
based on improved JITL technique. ISA transactions, 52(4), Robust Mode Transition Control for a Series –Parallel Hybrid
pp.539-549. Electric Vehicle. IEEE Transactions on Vehicular Technology.
13. Bevrani, H., F. Habibi, P. Babahajyani, M. Watanabe, and Y. 2016 Mar;65(3):1059-69.
Mitani. "Intelligent frequency control in an acmicrogrid: online 34. Matar M, Paradis D, Iravani R. Real-time simulation of modular
PSO-based fuzzy tuning approach." Smart Grid, IEEE multilevel converters for controller hardware-in-the-loop testing.
Transactions on 3, no. 4 (2012): 1935-1944. IET Power Electronics. 2016 Jan 20;9(1):42-50.
14. Yesil, Engin. "Interval type-2 fuzzy PID load frequency controller 35. www.winddata.com [Online; accessed 10.10.14].
using Big Bang–Big Crunch optimization." Applied Soft 36. www.solargis.info/doc/solar-and-pv-data [Online; accessed
Computing 15 (2014): 100-112. 10.10.14]
15. Zheng, Shiqi, Xiaoqi Tang, Bao Song, Shaowu Lu, and Bosheng
Ye. "Stable adaptive PI control for permanent magnet synchronous
motor drive based on improved JITL technique." ISA transactions Mohammad-Hassan Khooban (M’13) was born in
52, no. 4 (2013): 539-549. Shiraz, Iran, in 1988. He received the Ph.D. degree
16. Pahasa, Jonglak, and IssarachaiNgamroo. "Coordinated control of from Shiraz University of Technology, Shiraz, Iran, in
wind turbine blade pitch angle and PHEVs using MPCs for load 2017. He was a research assistant with the University of
frequency control of microgrid." IEEE Syst. J. (2014) , 99, 1–9.. Aalborg, Aalborg, Denmark from 2016 to 2017
17. Khooban MH, Niknam T, Blaabjerg F, Dragičević T. A new load conducting research on microgrids and Marine Systems.
His research interests include control theory and
frequency control strategy for micro-grids with considering
application, power electronics and its applications in
electrical vehicles. Electric Power Systems Research. 2017 Feb
power systems, industrial electronics, and renewable
28;143:585-98. energy systems.
18. Singh, Vijay P., Soumya R. Mohanty, NandKishor, and Prakash K.
Ray. "Robust H-infinity load frequency control in hybrid Taher Niknam (M’14) received the B.Sc. degree from
distributed generation system." International Journal of Electrical Shiraz University, Shiraz, Iran, in 1998, and the M.Sc.
Power & Energy Systems 46 (2013): 294-305. and Ph.D. degrees from Sharif University of
19. Lee, Dong-Jing, and Li Wang. "Small-signal stability analysis of Technology, Tehran, Iran, in 2000 and 2005,
an autonomous hybrid renewable energy power generation/energy respectively, all in power electrical engineering. He is a
storage system part I: time-domain simulations." Energy Faculty Member of the Department of Electrical
Conversion, IEEE Transactions on 23, no. 1 (2008): 311-320. Engineering, Shiraz University of Technology. His
20. Mojica-Nava, Eduardo, Carlos Andrés Macana, and research interests include power sys- tem restructuring,
NicanorQuijano. "Dynamic population games for optimal dispatch impacts of distributed generations on power systems,
on hierarchical microgrid control." Systems, Man, and optimization methods, and evolutionary algorithms.
Cybernetics: Systems, IEEE Transactions on 44, no. 3 (2014): 306-
317. Mokhtar Shasadeghi received the B.Sc. degree in
21. Simon, C., Robert, J.:"Geometric Type-1 and Type-2 Fuzzy Logic electronics engineering from Shiraz University,
Shiraz, Iran, in 1996, and the M.Sc. and Ph.D.
Systems", IEEE Transactions on Fuzzy Systems, 2007, 15, (1), pp.
degrees from Tarbiat Modares University, Tehran,
3-15.
Iran, in 2001 and 2007, respectively. He is currently
22. Khalghani, Mohammad Reza, Mohammad Hassan Khooban, an Associate Professor at Shiraz University of
EsmaeilMahboubi-Moghaddam, NavidVafamand, and Mohammad Technology, Shiraz. His research interests include
Goodarzi. "A self-tuning load frequency control strategy for robust control, fuzzy control, time-delay systems,
microgrids: Human brain emotional learning." International linear matrix inequalities, and sum- of-squares
Journal of Electrical Power & Energy Systems 75 (2016): 311-319. decomposition and optimization.
23. Khooban MH, Niknam T, Blaabjerg F, Davari P, Dragicevic T. A
robust adaptive load frequency control for micro-grids. ISA Tomislav Dragicevic (S’09-M’13-SM’17) received
transactions. 2016 Nov 30;65:220-9. the M.E.E. and the industrial Ph.D. degree from the
24. Bevrani H, Feizi MR, Ataee S. Robust Frequency Control in an Faculty of Electrical Engineering, Zagreb, Croatia, in
Islanded Microgrid: and-Synthesis Approaches. IEEE Transactions 2009 and 2013, respectively. From 2013 until 2016 he
on Smart Grid. 2016 Mar;7(2):706-17. has been a Postdoctoral research associate at Aalborg
25. Pan I, Das S. Intelligent Fractional Order Systems andControl. An University, Denmark. From March 2016 he is an
Introduction. Studies in Computational Intelligence.;438. Associate Professor at Aalborg University, Denmark.
26. Das S, Pan I. On the Mixed Loop-Shaping Tradeoffs in Fractional- His field of interest is overall system design of
autonomous and grid-connected DC and AC
Order Control of the AVR System. IEEE Transactions on
microgrids, and industrial application of advanced modeling, control and
Industrial Informatics. 2014 Nov;10(4):1982-91.
protection concepts to shipboard power systems, remote telecom stations,
27. Zhong J, Li L. Tuning Fractional-Order Controllers for a Solid- domestic and commercial facilities and electric vehicle charging stations.
Core Magnetic Bearing System. IEEE Transactions on Control
Systems Technology. 2015 Jul;23(4):1648-56.

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TSTE.2017.2763607, IEEE
Transactions on Sustainable Energy
IEEE TRANSACTIONS ON SUSTAINABLE ENERGY 9

Frede Blaabjerg (S’86-M’88-SM’97-F’03) was with

ABB-Scandia, Randers, Denmark, from 1987 to 1988.
From 1988 to 1992, he was a Ph.D. Student with
Aalborg University, Aalborg, Denmark. He became an
Assistant Professor in 1992, an Associate Professor in
1996, and a Full Professor of power electronics and
drives in 1998. His current research interests include
power electronics and its applications such as in wind
turbines, PV systems, reliability, harmonics and adjustable speed drives. He
has received 17 IEEE Prize Paper Awards, the IEEE PELS Distinguished
Service Award in 2009, the EPE-PEMC Council Award in 2010, the IEEE
William E. Newell Power Electronics Award 2014 and the VillumKann
Rasmussen Research Award 2014. He was an Editor-inChief of the IEEE
TRANSACTIONS ON POWER ELECTRONICS from 2006 to 2012. He has
been Distinguished Lecturer for the IEEE Power Electronics Society from
2005 to 2007 and for the IEEE Industry Applications Society from 2010 to
2011. He is nominated in 2014 and 2015 by Thomson Reuters to be between
the most 250 cited researchers in Engineering in the world.

1949-3029 (c) 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.