Next Article in Journal
Modelling and Simulation/Optimization of Austria’s National Multi-Energy System with a High Degree of Spatial and Temporal Resolution
Previous Article in Journal
Uniform and Non-Uniform Pumping Effect on Ce:Nd:YAG Side-Pumped Solar Laser Output Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiobjective Optimization of a Hybrid PV/Wind/Battery/Diesel Generator System Integrated in Microgrid: A Case Study in Djelfa, Algeria

1
Laboratory of Applied Automation and Industrial Diagnostics (LAADI), Faculty of Science and Technology, Ziane Achour University, Djelfa 17000, Algeria
2
Electrical and Electronics Engineering Department, Nisantasi University, Istanbul 34398, Turkey
3
Institute for High-Power Converter Systems (HLU), Technical University of Munich (TUM), 80333 Munich, Germany
4
LACoSERE Laboratory, Department of Electrical Engineering, Faculty of Technology, University Amar Telidji of Laghouat, Laghouat 03000, Algeria
5
Department of Electrical Engineering, Assiut University, Assiut 71516, Egypt
*
Authors to whom correspondence should be addressed.
Submission received: 7 April 2022 / Revised: 7 May 2022 / Accepted: 10 May 2022 / Published: 13 May 2022

Abstract

:
Hybrid Renewable Energy Sources (HRES) integrated into a microgrid (MG) are a cost-effective and convenient solution to supply energy to off-grid and rural areas in developing countries. This research paper focuses on the optimization of an HRES connected to a stand-alone microgrid system consisting of photovoltaics (PV), wind turbines (WT), batteries (BT), diesel generators (DG), and inverters to meet the energy demand of fifteen residential housing units in the city of Djelfa, Algeria. In this context, the multiobjective salp swarm algorithm (MOSSA), which is among the latest nature-inspired metaheuristic algorithms recently introduced for hybrid microgrid system (HMS) optimization, has been proposed in this paper for solving the optimization of an isolated HRES. The proposed multiobjective optimization problem takes into account the cost of energy (COE) and loss of power supply probability (LPSP) as objective functions. The proposed approach is applied to determine three design variables, which are the nominal power of photovoltaic, the number of wind turbines, and the number of battery autonomy days considering higher reliability and minimum COE. In order to perform the optimum size of HMG, MOSSA is combined with a rule-based energy management strategy (EMS). The role of EMS is the coordination of the energy flow between different system components. The effectiveness of using MOSSA in addressing the optimization issue is investigated by comparing its performance with that of the multiobjective dragonfly algorithm (MODA), multiobjective grasshopper optimization algorithm (MOGOA), and multiobjective ant lion optimizer (MOALO). The MATLAB environment is used to simulate HMS. Simulation results confirm that MOSSA achieves the optimum system size as it contributed 0.255 USD/kW h of COE and LPSP of 27.079% compared to MODA, MOGOA, and MOALO. In addition, the optimization results obtained using the proposed method provided a set of design solutions for the HMS, which will help designers select the optimal solution for the HMS.

1. Introduction

Population increase, rising energy consumption, rising energy production costs, greenhouse gas emissions, and damage caused by fossil resources to the environment have led to a greater emphasis on renewable energy sources [1,2,3,4].
Although there have been considerable increases in renewable energy technology, electricity shortages remain a problem in rural communities and islands. According to the United Nations Development Program (UNDP), more than a quarter of the world’s population, particularly those living in rural regions, suffer from lack of access to electricity. Rural communities are typically isolated from the national grid and situated in challenging terrain, such as steep terrain or dense jungle, making transmission line extension prohibitively expensive or impossible. Renewable resources, such as wind, solar, and hydropower, help these areas save energy because they possess potentially ubiquitous, abundant, accessible, clean, and easy-to-obtain energy [5,6].
The most accessible and critical renewable resources, wind and solar energy, are combined to create a higher quality and more reliable hybrid energy system than individual resources [1,2,3,4]. In fact, in a hybrid energy system, wind and solar energy are the primary energy sources. Battery units, energy storage, fuel cells, and DGs can all be used in hybrid systems to improve efficiency and eliminate flaws. Indeed, when wind speed or solar radiation falls below a certain threshold or when a peak in demand arises, the presence of these storage devices becomes critical [4,7,8].
Renewable energy in the form of an HMS System is the most efficient, dependable, and cost-effective way to use localized renewable energy resources which combine renewable energy sources with a diesel generator as a backup system to provide centralized electricity generation at the local level. These systems range in size from a basic 5 kW single-phase system for a single residence to a massive three-phase network that serves as a primary power supply for the entire neighborhood. When demand grows or the community expands, they can quickly scale up and connect to the national grid [6,9].
Microgrid systems based on autonomous renewable energy sources (RES) are the most viable and cost-effective option for electrifying off-grid areas [10,11]. Therefore, from an economic and technical standpoint, such a system’s planning and optimal design are complex for various reasons. One of these reasons would be the inconsistency of RESs and their reliance on meteorological conditions. Microgrid systems are frequently either over- or undersized in order to meet energy demand. An oversized system will have a high operating cost and will create extra energy. Conversely, an undersized microgrid system will fail to supply electricity to the required loads. To fully realize the benefits of an RES-based microgrid, optimum microgrid sizing combined with a robust energy management strategy is required [11,12].

1.1. Literature State of the Art

A large number of works has been reported in the field of hybrid microgrid system sizing. Previous methodologies in the area of sizing can be grouped into different categories [13,14,15,16]. The first category—software tools such as HOMER, HOMER Pro, PVSYST, HOGA, IHOGA, and RAPSIM—was applied to optimize microgrid systems [17,18,19,20,21,22,23]. Despite being simple to use, this category has the drawback of users being unable to select appropriate system components intuitively. Furthermore, users have no access to or visibility of the calculations and algorithms [19,24]. The second category contains deterministic methods such as iterative, analytical, numerical, graphical construction, etc. [25,26,27,28,29,30,31]. Although simple, these methods require considerable simulation time, as all system configurations are analyzed [26]. The third category includes metaheuristic algorithms for solving microgrid problems, many of which have been utilized in the literature to tackle optimization problems where authors [32] developed a recent methodology based on social spider optimizers (SSO) to determine the optimal sizing of an HRES-integrated microgrid. This group comprises PV, WT, battery, DG, and inverter with COE as fitness function and was presented for sensitivity analysis of sizing different topologies of MG, including PV/battery/DG, WT/battery/DG, and PV/WT/battery/DG in Aljouf Region, Saudi Arabia. Another study [11] focused on the application of the grasshopper optimization algorithm (GOA) in the area of microgrid system sizing design problems in order to determine the optimal system configuration comprising PV, WTs, BSS, and DGs with COE as a single objective function in Yobe State, Nigeria. In [33], a novel bonobo optimizer (BO) technique was applied to find the optimal design for an off-grid HRES that contains a DG, PV, a WT, and batteries as a storage system in Saudi Arabia. It is based on annualized system cost (ASC) minimization and power system reliability enhancement. A hybrid power generation system consisting of DGs, PV systems, and battery energy storage (BES) has been optimized on the basis of BES dispatch (OBD) by considering two objectives in Indonesia: the levelized cost of energy (LCOE) and renewable energy (RE) [34].
A particle swarm optimization (MOPSO) method was proposed in [6], which dealt with the optimization problem consisting of COE and LPSP minimization for an HMS system and has been tested on three selected stations in Iran. Elsewhere, the authors of [35] proposed a multiobjective problem formulation to solve the optimization task with a self-adaptive differential evolution (DE) algorithm. This latter research has been used to analyze LPSP, COE, and RF for the city of Yanbu, Saudi Arabia for three cases of load profile. Another study presented a multiobjective evolutionary algorithm to optimally design an HMS considering load uncertainty [36]. A MOGOA was applied in [37] to optimize a microgrid with a rule-based energy management scheme. In [38], an HMS system comprising PV, wind, and battery with a diesel generator was introduced to supply load demand in Shlateen, Egypt with two planning scenarios. The first is PV/wind/battery and the second is PV/wind/battery/diesel. A multiobjective optimization problem has been solved by considering three objective functions—COE, LPSP, and RF—based on MODA. A hybrid photovoltaic, diesel, and battery nanogrid system installed in Saudi Arabia has been introduced and optimized by the parallel multiobjective PSO-based approach (PMOPSO) in [39].
Table 1 summarizes the details of the reported approaches.
These algorithms are flexible; they avoid falling into the trap of local optima and offer a better solution than other methods. In addition, these methods have various advantages that make them worthy of solving any kind of optimization problem because they simulate the problem-solving methods used by creatures and they offer promising results in solving such problems, but the same algorithm may show poor performance in another optimization problem [11]. However, based on the free lunch theory, it is always possible to develop new methods or algorithms that can solve the problem in a better way than the existing techniques [36]. To overcome these above-mentioned defects, MOSSA is chosen because of its simplicity, ease of use, and requirement of fewer control parameters, as it has shown impressive results in solving other engineering problems.

1.2. Article Contribution and Organization

To date, no study has been performed to show the effectiveness of MOSSA in multiobjective optimization of an HMS. MOSSA is a novel nature-inspired meta-heuristic optimization algorithm. The objective of this research paper is to design the optimal configuration of an HMS system. The main contributions of the paper can be summarized as follow:
  • Determination of the optimal sizing of PV, WT, BT, DG, and inverter integrated HMS based on a recent approach called MOSSA;
  • A rule-based EMS that manages the energy flow between different HRESs is proposed;
  • Analysis is performed using solar radiation, wind speed, and ambient temperature data obtained from the Djelfa region in Algeria;
  • The multiobjective optimization approach considered COE and LPSP as objective functions and renewable factor (RF) as a constraint.
The obtained results from MOSSA are compared with three algorithms: MODA, MOGOA, and MOALO. The comparison clearly shows that MOSSA performs better, with very fast convergence and balance between exploitation and exploration.
The rest of this paper is organized as follows: Section 2 is dedicated to the modeling of the hybrid microgrid system components, Section 3 presents the study site and the system specifications, Section 4 focusses on the proposed EMS. Whereas, Section 5 presents the obtained results based on the used optimization technique, and Section 6 contains discussion on the obtained results. This paper ends with a conclusion.

2. Modeling of Hybrid Microgrid System Components

The stand-alone microgrid of the HRES proposed in the present study is composed of two renewable energy sources (PV system, WT), an energy storage unit (battery), a backup power source (DG), and two types of loads. The PV system, WT, and battery are connected to the DC bus via bidirectional converters; in contrast, DG, load, and dump load are connected directly to the single-phase AC bus. Figure 1 shows the whole configuration of the studied stand-alone microgrid system. It is worthy to clarify that this stand-alone microgrid system is designed to be a low-voltage distribution network 220 V, 50 Hz that supplies power to a single-phase AC power system.

2.1. PV Array Modeling

The photovoltaic (PV) panels are obtained based on appropriate connections of photovoltaic cells. To provide the required output voltage in the PV panel, PV cells are connected in series to form a string. On the other side, to increase the output current capability of a panel, similar strings are connected in parallel, which forms a PV panel or a PV module whose output power rating is obtained by multiplying the voltage and the current at the output side [26,32,41].
In [42], several models were developed for the calculation of PV panel output power P p v _ o u t . However, in this study, a simple model has been used which takes into consideration two variables: ambient temperature and solar irradiance. The output power of this model can be obtained as follows [6,11,43,44]:
P p v o u t t = P N p v × G t t G t S T C × 1 + α t T C t T C S T C
where P p v _ o u t is the output power of the PV module (W); P N _ p v is the rated power (W) of the PV module at standard test condition (STC), which is usually stipulated by the manufacturer, G t is the solar irradiance (W/m2); G t _ S T C is the solar irradiance at STC ( G t _ S T C = 1000   W / m 2 ); α t is the temperature coefficient, its value is α t = 3.7 × 10 3 1 / ° C for the mono- and polycrystalline (Si) solar cells [45,46]; T C _ S T C is the cell temperature at STC ( T C _ S T C = 25   ° ); and T a m b is the ambient temperature (°C).
The cell temperature T c is calculated by the following expression [44]:
T C t = T a m b t + 0.0256 × G t t

2.2. Wind Turbine Modeling

The output power of the wind turbine generator varies following variation of wind speed, which is a function of the height at the same location. Thus, the measured wind speed at the anemometer height must be adjusted to the used hub height of the wind turbine [12]. In this context, two mathematical models are used for the determination of the wind speed vertical profile of a specific location such as the log-law and the power-law [47,48].
In this study, the determination of the wind profile is obtained by using the power-law model as follows [48]:
V 2 V 1 = h 2 h 1 α
where V 2 m / s is the wind speed at the hub height h 2 m , V 1 m / s is the wind speed at the reference height h 1 m , and α is the friction coefficient (also known as: Hellmann exponent, wind gradient, or power-law exponent). α is a function of varying parameters such as wind speed, roughness of terrain, height above ground, temperature, hour of the day, and time of the year [6,49,50,51]. The value of the friction coefficient can be considered to be 0.11 at extreme wind conditions and 0.2 at normal wind conditions as specified by Standard IEC 61400-1 [49,52,53]. However, a value of (α = 1/7) is commonly accepted [6].
The output power of wind turbine generator can be evaluated based on the following equation [54,55,56]:
P w t t =   0 V < V c u t i n   V 3 P r V r 3 V c u t i n 3 P r V c u t i n 3 V r 3 V c u t i n 3 V c u t i n V < V r a t e d   P r V r a t e d V V c u t o u t   0 V > V c u t o u t
where P r is the rated power (kW); V is the wind speed (m/s); V c u t _ i n , V r a t e d , V cut out   represent the cut-in, rated, and cut-out wind speed of the WT, respectively in m/s. These values are usually stipulated by the manufacturer. Figure 2 shows the output power of the wind turbine generator versus the wind speed at the height of the hub.
The rated power P r (W) of a wind turbine generator is expressed as a function of the area swept by the blades A w i n d , the maximum power coefficient C p , the air density ρ a i r , and the rated wind speed, as follows [55]:
P r = 1 2 C p × ρ a i r × A w i n d × V r 3

2.3. Battery Bank Modeling

The battery system is an important unit in the stand-alone microgrid and ensures supply of the load when the power generated from renewable energy sources is insufficient.
Battery capacity can be calculated as follows [57,58]:
C B a t t = A D × P l o a d η I n v × η B a t t × D O D
where P l o a d is the power demand of the load, η I n v is the inverter efficiency, η B a t t is the battery efficiency, D O D is the battery depth of discharge, and A D is the days of autonomy, which is defined as the number of days that the battery will be able to supply the required power demand of the load without deficiency.
It is obvious that the energy generated from renewable energy sources (PV and WT) depends on wind speed and solar radiation, which are intermittent in nature; therefore, days of autonomy is of great importance and must be taken into account while sizing the battery bank to overcome the problem of power deficit production from these sources. In case of energy production excess, the excess is used to charge the battery. The power produced from the battery bank can be expressed as follows [40]:
P B a t t t = P p v t + P w t t P l o a d t η I n v
where P p v t , P w t t , and P l o a d t represent the power produced from PV, WT, and load power demand, respectively, and η I n v is the inverter efficiency.
When P B a t t t < 0 , this indicates that there is an energy generation deficit. Otherwise, if P B a t t t > 0 , it indicates that the energy generation exceeds the power demand. In the rare case that P B a t t t = 0 , the generated power from the renewable sources is equal to the load power demand.
For checking the state of the battery bank, the state of charge (SOC) of the battery is an important parameter that affects battery performance and indicates its current capacity. Indeed, the SOC can be defined following both the states of charging and discharging as follows [59,60]:
  • Charging process, if; P p v t + P w t t > P l o a d t
S O C t = S O C t 1 × 1 σ + P p v t + P w t t P l o a d t η I n v × η B a t t
  • Discharging process, if; P p v t + P w t t < P l o a d t
S O C t = S O C t 1 × 1 σ + P l o a d t η I n v P p v t + P w t t × η B a t t

2.4. Diesel Generator Modeling

The diesel generator is used as a backup source if the renewable energy sources and battery bank are insufficient to meet the load power demand requirement. The following equation can be used to calculate the hourly fuel consumption of the diesel generator [61,62]:
F D G t = α P D G t + β P r
where F D G t is the generator fuel consumption L / h r , P D G t is the generated power k W , P r is the capacity of the generator k W , α is the fuel curve slope coefficient ( L / h r / k W o u t p u t ), and β is the fuel intercept coefficient ( L / h r / k W r a t e d ). α and β used in the present study are taken as α = 0.246   and   β = 0.08415 [61,62].
The efficiency of the diesel generator can be calculated as [63]:
η o v e r a l l = η b r a k e t h e r m a l × η g e n e r a t o r
where η o v e r a l l , η g e n e r a t o r , and η b r a k e t h e r m a l represent the overall efficiency, generator efficiency, and brake thermal efficiency of the diesel generator, respectively.

2.5. Inverter Modelling

The inverter converts the electrical energy from direct current (DC) to alternating current (AC). The inverter’s efficiency is defined as follows [6,64,65]:
η i n v = P P + P 0 + K P 2
where P , P 0 , and K can be calculated by the following equations [6,64,65]:
  • P = P o u t / P n
  • P 0 = 1 99 ( 1 / η 10 1 / η 100 9 ) 2
  • K = 1 / η 100 P 0 1
where P n is the inverter’s rated power and η 10 and η 100 are the efficiencies of the inverter at 10 and 100% of the rated power, respectively. Both η 10 and η 100 are stipulated by the manufacturer.

3. Definition of the Study Site and System Specifications

3.1. Location and Meteorological Conditions

The proposed stand-alone microgrid of the hybrid renewable energy system is supposed to be located in Aïn El Ibel, Djelfa in the north-central region of Algeria at 34.346° latitude and 3.163° longitude. It is situated in a transitional zone between the dry high plains in the north and the desert in the south. The area is characterized by very hot weather in the summer and very cold weather in winter, with high wind speeds in the winter. Figure 3 shows the geographical location of the study site, while Table 2 presents the related background information and the period of data measurement.
Wind speed, solar radiation, and ambient temperature used for the simulations in this study were obtained from the National Aeronautics and Space Administration (NASA) [66] at the coordinates of Aïn El Ibel (34.346° latitude and 3.163° longitude). The average solar radiation of the selected location during the mentioned year is 0.2357 KW/m2, the average wind speed is 4.3467 m/s, and the average ambient temperature is 290.4010 K. Figure 4, Figure 5 and Figure 6 show the hourly solar radiation profile, the wind speed at the height of 10 m above ground level, and the ambient temperature over a year (8784 h).

3.2. Load Assessment

To design a reliable and efficient stand-alone microgrid based on renewable energy sources system for a specific location, the load profile must be defined precisely. Indeed, it is obvious that the fluctuating nature of the load affects the system reliability, the number of components, and the price of the produced electrical energy. Therefore, in the present study, it is supposed that the studied stand-alone microgrid will power a specific load of a residential area composed of a group of homes with each home presented as a residential unit as detailed in Table 3.
The microgrid is analyzed for ten residential units in an off-grid community. The studied load profile is presented within the four seasons of the year—winter, spring, summer, and autumn—and has an average annual power consumption of 12.04545 kW. Figure 7 and Figure 8 present the hourly load profile of the residential units.

3.3. Specifications of Hybrid Microgrid System Components

The economic and technical parameters of hybrid microgrid system components used in this study are listed in Table 4:

4. Energy Management Strategy of Hybrid Microgrid System

The energy management strategy (EMS) is one of the main criteria to be considered when designing or sizing an autonomous microgrid which is intended to ensure the distribution and management of power flow among the various elements of the studied stand-alone microgrid system. The main targets of the proposed EMS can be summarized as follows:
  • System efficiency enhancement, thus achieving low cost and energy-saving benefits;
  • Maximization of utilization of the renewable energy sources (PV and WT);
  • Protection of the battery bank and minimization of its degradation;
  • Minimization of fuel consumption.
In this study, four modes of energy management strategy (EMS) have been used:
  • Mode 1: In this mode, the generated power from renewable energy sources (PV and WT) is sufficient to supply the load demand requirement. The extra energy is used to charge the battery bank system;
  • Mode 2: In this mode, the generated power from the renewable energy sources exceeds the load demand requirement while the battery is fully charged. In this case, the surplus of energy is consumed in a dump load;
  • Mode 3: In this mode, the generated power from the renewable energy sources is less than the load demand requirements. In this case the battery bank will cover the power generation deficiency to fit the load demand requirements;
  • Mode 4: In this mode, the power generated from the renewable energy sources is not sufficient to meet the load demand requirement and at the same time, the battery bank storage level is low. In this case, the diesel generator will operate to cover the gap in power generation to fit the load demand requirement and further ensure the battery bank’s charging
Figure 9 shows the flowchart of the proposed EMS algorithm used in the present study.

5. Optimization Problem Formulation of the Studied Stand-Alone Microgrid

5.1. Multiobjective Optimization

In this work, the optimization problem of the studied hybrid microgrid system (HMGS) is based on solving a multiobjective optimization problem which is formulated as follows:
Minimize : f x = f 1 x f 2 x                         f n x
Subjected   to H x = 0 G x 0
where f x = f 1 x , f 2 x , , f n x is the vector of the objective functions and n is the number of objective functions, f i x is the objective function number i , x = x 1 , x 2 , x 3 , , x n is the vector of the variables to be designed, H x is the set of equality constraints, and G x is the set of inequality constraints.

5.2. Objective Functions

In order to evaluate microgrid performances, the loss of power supply probability (LPSP) and the cost of energy (COE) are proposed as objective functions, where the main aim is to minimize these two functions in order to maintain the high reliability and the lowest cost of the studied HMGS.

5.2.1. Loss of Power Supply Probability

The reliability of the microgrid system is evaluated based on LPSP. It is a statistical term which indicates the probability of power supply failing to meet the load demand requirement due to technical problems or because the energy produced from renewable sources is insufficient to meet the load demand requirement. The LPSP can be calculated using the following equation [67,68]:
LPSP =   ( P l o a d P p v P w t + P S O C m i n + P d i e s e l )   P l o a d
where the value of LPSP is in the range between [0, 1]. If its value is equal to zero, the load energy demand is fully satisfied. On other hand, if the LPSP equal to unity, the load energy demand is not satisfied.
In this work, the evaluation of the system’s reliability is taken within the worst conditions as follows [6,35]:
P ( t ) l o a d > P ( t ) g e n e r a t e
where P ( t ) g e n e r a t e denotes the power generated. It should be noted that under this condition, the total load energy demand is greater than the entire energy generation from all the available sources.

5.2.2. Cost of Energy

The COE is one of the most important indicators of the economic profitability of HRES integrated in microgrid [6,69]. It is defined as the unit of cost per unit of energy production from the HRES (USD/kW h), and it can be calculated for each component of the studied microgrid based on the following equations [32,70]:
C t P V = N   P V C C P V + C O & M P V × ( 1 + i ) n 1 i ( 1 + i ) n
C t W T = N   W T C C W T + C O & M W T × ( 1 + i ) n 1 i ( 1 + i ) n
C t B a t t = C C B a t t + C O & M B a t t × ( ( 1 + i ) n 1 i ( 1 + i ) n ) + C R B a t t ×   j = 1 ( n n B a t t 1 ) ( 1 + 1 ( 1 + i ) j n B a t t )
C t D G = C C D G + C O & M D G × ( ( 1 + i ) n 1 i ( 1 + i ) n ) + C R D G ×   j = 1 ( n n B a t t 1 ) ( 1 + 1 ( 1 + i ) j n D G )
The resulting net cost of the four sources can be calculated as follows:
NPC = C t P V + C t W T + C t Batt + C t DG + C t Inv
where C t P V , C t W T , C t B a t t , C t D G , and C t I n v are the energy costs of PV system, WT system, battery bank, and diesel generator, respectively; C C P V , C C W T , C C B a t t , and C C D G represent the investment costs put into these sources of energy, respectively; C O & M P V , C O & M W T , C O & M B a t t , and C O & M D G are the operation and maintenance costs of PV, WT, battery bank, and diesel generator, respectively; C R B a t t and C R D G are the replacement costs of the battery bank and the diesel generator, respectively; i is the annual interest; n is the system life time; and n B a t t and n D G are the life time of the battery bank and the diesel generator, respectively. The COE can be calculated as follows [69]:
COE = NPC     h = 1 8784   P   l o a d × C R F
where P l o a d is the hourly power consumption and C R F is the capital recovery factor, which is defined as follows [71]:
C R F = i ( 1 + i ) n ( 1 + i ) n 1

5.3. Constraints

Renewable Factor

The renewable factor ( R F ) is a factor that determines the amount of power generated from renewable resources in comparison to non-renewable resources (diesel generator) used in the microgrid, and it is expressed as follows [36]:
R F ( % ) = ( 1   P d i e s e l   P p v +   P w t ) × 100
when the R F equals 100%, it means that the system is in an ideal state and is dependent only on power generated from renewable energy resources. When it equals zero percent, it means that the power coming from the diesel generator is equivalent to the power generated from renewable energy resources.

5.4. Design Variables

The design variables considered in this study are the nominal power of the photovoltaic system ( P V ), the number of wind turbines ( N W T ), and the number of autonomy days ( N A D ). The constraints on the proposed design variables are given as follows:
Design   variables   : 0   P V 80 0 N   W T 10 0 N   A D 3

5.5. Multiobjective Salp Swarm Algorithm

Salp swarm algorithm (SSA) is a new swarm intelligence algorithm that belongs to the family of metaheuristic optimization algorithms. It has been proposed based on the foraging behaviors of salps and their intelligence, which was proposed by Mirjalili et al. in 2017 [72]. After the food source is detected in salp lookup mechanisms, salps gather in chains to search for food sources using these chains. Each series will follow the lead salp (the leader), and the leader will direct and lead the other salps to catch the food source [73]. The shape of a salp is shown in Figure 10a, and the salp chain is illustrated in Figure 10b.

5.5.1. Basic Salp Swarm Algorithm

The salp leader moves toward the food source named ( F j ) in the search space, while all followers can move toward the rest of the salp. The position of the salps is specified in the n-dimensional search space, where n is the number of variables for a given problem. Therefore, the position of all the salps is stored in a two-dimensional matrix called X i with a size of N × d , as described in the flowing equation [75]:
X i = x 1 1 x 2 1 x d 1 x 1 2 x 2 2 x d 2 x 1 N x 2 N x d N
The update of the leader’s position is performed based on the following equation [72]:
x j 1 = F j + c 1 u b j l b j c 2 + l b j c 3 0 F j c 1 u b j l b j c 2 + l b j c 3 < 0
where x j 1 is the position vector of the first salp (the leader) in the j th dimension, F j is the position vector of the food source in the j th dimension. u b j and l b j represent the upper bound and the lower bound of the j th dimension, respectively. c 1 , c 2 , and c 3 are random numbers that uniformly generated with values between [0, 1]. As indicated in Equation (7), the leader will only update its position with respect to the food source.
The coefficient c 1 balances between the exploration and exploitation, so it is considered one of the most important parameters in salp swarm algorithms (SSA), and it is defined as follows [72]:
c 1 = 2 e 4 l L 2
where l indicates the current iteration, and L indicates to the maximum number of iterations.
The followers update their positions based on Newton’s law of motion as follows [72]:
x j i = 1 2 a t 2 + v 0 t
where i 2 , x j i is the position vector of the i   t h follower salp in the j   t h dimension, t is the time, v 0 is the initial speed, and a = v f i n a l v 0 where v = x x 0 t .
Because the time in optimization represents the iteration based on a specified sampling unit, the discrepancy between iterations is equal to 1, and by considering that v 0 = 0 , the aforementioned equation can be expressed as follows [72]:
x j i = 1 2 x j i + x j i 1
where i     2 , and x j i is the follower salps’ position vector in the j t h dimension. According to the mathematical simulation described in the above equations, the swarm behavior of the salp chains can be clearly understood and easily simulated.

5.5.2. Multiobjective Salp Swarm Algorithm

A set of solutions called the Pareto group presents the solution to the problem of multiple objective function optimization. The SSA algorithm is based on the movement of the salps towards the food source and updating their position over the course of iterations. This first problem is addressed by equipping the SSA algorithm with a food source repository. This repository holds the best nondominant solutions obtained so far during optimization. The multiobjective salp swarm algorithm (MOSSA) approach inherits the SSA operators due to the similar population division (leader and follower salps) and the position updating process. The MOSSA algorithm is logically capable of finding the perfect Pareto solutions with high distribution across all objective functions [75,76,77]. The MOSSA pseudo-code is represented in Algorithm 1.
Algorithm 1. Pseudo code of the MOSSA algorithm [72,78]
1Set the hyper-parameter:
2Max_iter:Maximum of iteration
3ArchiveMaxSize:Max capacity of archive (repository)
4Dim:The number of parameters on each salp
5 u b j and l b j :The upper bound and the lower bound of salp population
6Obj-no:The objective number to be estimated
7Initialize the salp population x i i = 1 , 2 , , n considering u b j and l b j
8Define the objective function (loss function): @ Ob-func
9while (end criterion is not met) do
10 Calculate the fitness of each search agent (salp) with Ob-func
11 Determine the non-dominated salps
12 Update the repository considering the obtained non-dominated salps
13 if (the repository becomes full) then
14 Call the repository maintenance procedure to remove one repository resident
15 Add the non-dominated salp to the repository
16 end
17 Choose a source of food from repository: F = SelectFood (repository)
18 Update c 1   by c 1 = 2 e 4 l L 2
19 for each salp x i : do
20 if i = = 1 then
21 Update the position of the leading salp by:
22 x j 1 = F j + c 1 u b j l b j c 2 + l b j c 3 0 F j c 1 u b j l b j c 2 + l b j c 3 < 0
23 else
24 Update the position of the follower salp by:
25 x j i = 1 2 x j i + x j i 1
26 end
27 end
28 Amend the salps based on the upper and lower bounds of variables
29end
30return repository
To solve the sizing problem, MOSSA is used to follow several procedures. The proposed procedure of actions is defined in detail in Figure 11. When the simulation begins, like any other optimization algorithm, MOSSA places random particles in the search landscape, whose bounds are set by the user. According to the algorithm’s governing equations, these particles move around the search landscape, optimizing the defined objective function. To validate the effectiveness and performance of the MOSSA, the developers have carried out a series of tests, in which the performance of the MOSSA is compared with other well-known optimization algorithms like the multiobjective particle swarm optimization (MOPSO) and the non-dominated sorting genetic algorithm (NSGA-II) [72], and they concluded that MOSSA balances between exploration and exploitation in a satisfactory and improved manner compared to its counterparts.
MOSSA has recently been used to solve engineering problems, such as short-term load forecasting presented by J. Wang et al. [78] and for the enhancement of electrical distribution system performances including DG and DSTATCOM simultaneously proposed by A. Lasmari et al. [75]. The authors of [75,76,77,78] have affirmed that the result of the MOSSA outperforms other approaches, namely the multiobjective grasshopper optimization algorithm (MOGOA) and the multiobjective ant lion optimizer (MOALO). These studies have verified the effectiveness of the MOSSA in solving complex optimizations problems better than the other techniques.

6. Results and Discussion

In this paper, the optimization design of a stand-alone microgrid based on hybrid renewable energy system consisting of PV/WT/battery bank with a diesel generator system as a backup source is presented. These sources’ main role is to fulfill the load energy demand, which is composed of 15 residential housing units. The whole studied system has been implemented in the MATLAB R2018a environment. All tests were performed on an Intel Core i7-10510U CPU 2.30 GHz 16 GB RAM, Windows 10 Pro Version 21H2 (64-bit) personal computer.
This study comprises the application of four optimization algorithms to solve the main problem of optimization design presented in this paper:
  • First Algorithm: MOSSA optimization algorithm;
  • Second Algorithm: multiobjective dragonfly algorithm (MODA) [79];
  • Third Algorithm: multiobjective grasshopper optimization algorithm (MOGOA) [80];
  • Fourth Algorithm: multiobjective ant lion optimizer (MOALO) [81].
In a previous study [6,35], the optimization problem of HMS was treated only as a single-target optimization problem where different objective functions have been grouped into one objective function as single-objective optimization problem. This approach has the disadvantage of determining only one optimal solution. In this paper, multiple choices are offered for the optimal design of the studied microgrid based on the approach of multiobjective optimization. This approach creates a set of optimal solutions known as the Pareto front.
MOSSA, MODA, MOGOA, and MOALO have been executed along 100 iterations. Table 5 shows the different control parameters, which have been used for each proposed approach.
The solutions obtained by applying the previously mentioned methods for stand-alone microgrid configurations present an impressive and consistent distribution. Figure 12a–d presents the Pareto front of the stand-alone microgrid system obtained for each algorithm based on LPSP and COE functions.
On the Pareto front, the results show not only an optimal solution, but a group of optimal solutions (non-dominated solutions), a variety of design decision possibilities. The next section presents the obtained results of MOSSA, MODA, MOGOA, and MOALO algorithms.
Table 6 shows 20 selected solutions from the Pareto fronts for the first algorithm as shown in Figure 12a. For more clarification, these solutions have been organized based on the COE. If solution #1 is selected by the designer, the power generated by the PV panels is 65.883 kW, the autonomy days is equal to three days, and the required number of wind turbines is ten, which ensures a power generation of 20 kW. This solution corresponds to a COE of 0.255 USD/kW h, an LPSP of 27.079%, and an RF of 90.46%. It can be noted that the COE and the RE have higher values and the LPSP has a lower value than the other available solutions.
In case when solution #10 is selected by the designer, the power generated by the PV panels is 41.726 kW, the autonomy days is equal to 2.711, and ten wind turbines are required, which generates a total power of 20 kW. This solution corresponds to a COE of 0.187 USD/kW h, an LPSP of 28.424%, and an RF of 85.90%.
If solution #20 is selected by the designer, the total power generated by the PV panels is 24.093 kW, one autonomy day is required, and nine wind turbines are required to generate a total power of 18 kW. This solution corresponds to a COE of 0.133 USD/kW h, an LPSP of 31.121%, and an RF of 75.24%.
The annual contributions of energy generated by each source in the studied power generation system (i.e., PV, WT, BT, and DG) for solutions #1, #10, and #20 of the first algorithm are represented in Figure 13a–c, respectively. The annual contribution of the energy produced by the PV panels is 121.80 MW, which presents 42% of the total generated energy in solution #1 of the first algorithm, shown in Figure 13a and Table 6. It can be clearly observed that this production is higher compared to solution #10 and solution #20. The WT and DG provide 38.87 MW and 15.33 MW, respectively, which represent 13 and 5%, respectively, of the total generated energy in solution #1. It can be clearly noted that is lower than the energy contribution of the WT and DG in solution #10 and solution #20. The BT contributes 118.33 MW, which represents 40% of the total generated energy in solution #1. In Figure 13b, the annual contributions of the energy provided by the PV, WT, BT, and DG are 77.14 MW, 38.87 MW, 102.74 MW, and 16.36 MW, representing 33, 16, 44, and 7% of the total generated energy in solution #10, respectively.
Figure 13c shows the annual contribution of the energy provided by each source—PV, WT, BT, and DG—as 44.54 MW, 36.58 MW, 35.83 MW, and 20.09 MW, which represent 32, 27, 26, and 15% of the total generated energy in solution #20, respectively.
The same analysis can also be made for the results of MODA, MOGOA, and MOALO when the algorithms are changed while maintaining the same system.

Comparison between MOSSA and the Other Techniques

In this section, the performance of the MOSSA algorithm has been compared to the MODA, MOGOA, and MOALO algorithms for solving the optimization problem investigated in this paper. From the previously obtained results in Figure 12a–d and Table 6, Table 7, Table 8 and Table 9, MOSSA provided the best solution (Solution #1, Table 6) among the optimal solutions, corresponding to a COE of 0.255 USD/kW h and an LPSP of 27.079%. MODA, MOGOA, and MOALO algorithms provided more expensive solutions than MOSSA with a COE of 0.286 USD/kW h and an LPSP of 27.149%; a COE of 0.307 USD/kW h and an LPSP of 27.021%; and a COE of 0.290 USD/kW h and an LPSP of 27.018%, respectively. In addition, MOSSS balanced the contributions of renewable energy sources in a way that corresponds to the objective function by a PV of 42%, WT of 13%, batteries of 40%, and DG of 5%. The other contributions of energy generation by PV/WT/battery and DG for MODA, MOGOA, and MOALO are presented in Figure 14, Figure 15 and Figure 16, respectively. On the other hand, MOGOA provided the best RE result among the previous algorithms, and at the same time, it had the highest COE among them, indicating that it did not succeed in balancing between COE and LPSP.
As mentioned previously, this work attaches importance to preserving the environment. In this sense, MOSSA provided the lowest fuel consumed by DGs of 15.33 MW annual contribution compared to others, where contribution in MODA was 16.24 MW, contribution in MOGOA was 16.26 MW, and contribution in MOALO was 15.38 MW.
Figure 17 summarizes the comparison results between the previous algorithms; it can be noted clearly that the MOSSA algorithm offers better convergence than others, as it has a broader and consistent distribution in the Pareto front. It can be said that the MOSSA has good performance in improving system cost and reliability, and it provides better COE and LPSP. The obtained results indicate that the proposed MOSSA possesses good performance in solving such types multiobjective problems.
After comparing the previous algorithms, the solution #1 corresponding to the first algorithm is chosen, which ensures the optimal design of the studied stand-alone microgrid system of HRES consisting of PV/WT/battery and DG system to fulfill the load energy demand as shown in Figure 1. The obtained optimal solution allows ensuring total annual power generations of 121.80 MW, 38.87 MW, 118.33 MW, and 15.33 MW by PV, WT, battery, and DG, respectively, as presented in Table 6.
The annual generated power by the PV and wind turbine according to the selected optimal system configuration are presented in Figure 18a,b, respectively. The contribution of the backup system (diesel generator) is shown in Figure 18c. At the same time, the SOC of the battery is shown in Figure 18d, while Figure 19 presents the total annual power generated by all sources—PV, WT, battery, and DG—and the load power demand.
For example, if the period [8300; 8200] in hours during the winter season is chosen in order to observe the behavior of the SOC of the battery, a zoom is taken from Figure 18d and shown in Figure 20, while Figure 21 shows the contribution of the diesel generator in the same period.

7. Conclusions

This research paper has proposed an approach for the optimal sizing of a stand-alone microgrid system which is powered by hybrid sources of energy such as PV, WT, BT, and DG and comprises two buses such as a DC bus and an AC bus that are bidirectionally connected with a power electronics conversion system. The main aim of the designed hybrid stand-alone microgrid system is to meet the load energy demand of a remote area in an off-grid community in Djelfa, Algeria. Indeed, this paper has proposed an optimization design based on a recent nature-inspired metaheuristic optimization algorithm called the multiobjective salp swarm algorithm in order to find the minimum cost of electricity COE and minimum potential for electrical loss LPSP at the same time. Moreover, the renewable energy factor has been introduced to ensure that the system works mainly on renewable energy resources. To prove the validity and the effectiveness of the proposed approach using MOSSA, simulations have been carried out based on real data collected from the study site. The obtained results based on MOSSA have been compared with those obtained using MODA, MOGOA, and MOALO and some previously obtained results, where it has been proved that the proposed approach based on MOSSA yields better results; it contributed 90.46% of RF, a COE of 0.255 USD/kW h, and an LPSP of 27.079%. This outperforms MODA, MOGOA, and MOALO, where the offered optimal solutions show clear superiority compared to other presented results due to the design goal of the hybrid stand-alone microgrid system being to obtain a set of solutions to be adopted following different scenarios which can be faced in real cases. Finally, it can be concluded that the proposed approach based on the MOSSA can be a promising tool for the designer to obtain optimal design solutions.
As part of future work, other renewable energy sources can be added and diversified in line with the capabilities of each region. The newly discovered artificial intelligence algorithms can also be applied. Another perspective is the study of optimal sizing using hybrid optimization techniques for further explorations in future research and facing the complexities and challenges of hybrid systems.

Author Contributions

Conceptualization, Z.B., A.K., B.T., L.M. and A.B.; methodology, Z.B., A.K. and B.T.; validation, Z.B., A.K. and B.T.; formal analysis, Z.B., A.K., B.T., L.M., A.B. and M.A.; investigation, A.K., L.M. and M.A.; resources, A.K., B.T., L.M. and A.B.; data curation, Z.B.; writing—original draft preparation, Z.B.; writing—review and editing, A.K., B.T., L.M., R.K. and M.A.; visualization, Z.B., B.T. and A.K.; supervision, A.K.; project administration, A.K.; funding acquisition, R.K. and M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This project was financially supported by the Directorate General for Scientific Research and Technological Development—Algerian Ministry of Higher Education and Scientific Research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Wind speed, solar radiation, and ambient temperature data presented in this study are openly available at http://www.soda-pro.com/web-services/meteodata, accessed on 6 October 2021. The other data presented in this study are available in the body of this paper.

Acknowledgments

This work was supported by the German Research Foundation (DFG) and the Technical University of Munich (TUM) in the framework of the Open Access Publishing Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dalton, G.J.; Lockington, D.A.; Baldock, T.E. Case study feasibility analysis of renewable energy supply options for small to medium-sized tourist accommodations. Renew. Energy 2009, 34, 1134–1144. [Google Scholar] [CrossRef]
  2. Kaldellis, J.K.; Zafirakis, D. The wind energy (r)evolution: A short review of a long history. Renew. Energy 2011, 36, 1887–1901. [Google Scholar] [CrossRef]
  3. Fathabadi, H. Novel grid-connected solar/wind powered electric vehicle charging station with vehicle-to-grid technology. Energy 2017, 132, 1–11. [Google Scholar] [CrossRef]
  4. Ghorbani, N.; Kasaeian, A.; Toopshekan, A.; Bahrami, L.; Maghami, A. Optimizing a hybrid wind-PV-battery system using GA-PSO and MOPSO for reducing cost and increasing reliability. Energy 2018, 154, 581–591. [Google Scholar] [CrossRef]
  5. United Nations Development Programme (UNDP). UNDP and Energy Access for the Poor: Energizing the Millennium Development Goals. In Environment & Energy, United Nations Development Programme; UNDP: New York, NY, USA, 2011. [Google Scholar]
  6. Borhanazad, H.; Mekhilef, S.; Gounder Ganapathy, V.; Modiri-Delshad, M.; Mirtaheri, A. Optimization of micro-grid system using MOPSO. Renew. Energy 2014, 71, 295–306. [Google Scholar] [CrossRef]
  7. Celik, A.N. Techno-economic analysis of autonomous PV-wind hybrid energy systems using different sizing methods. Energy Convers. Manag. 2003, 44, 1951–1968. [Google Scholar] [CrossRef]
  8. Dufo-López, R.; Bernal-Agustín, J.L. Design and control strategies of PV-diesel systems using genetic algorithms. Sol. Energy 2005, 79, 33–46. [Google Scholar] [CrossRef]
  9. Wiemann, M.; Rolland, S.; Glania, G. Hybrid Mini-Grids for Rural Electrification: Lessons Learned. In Proceedings of the Region 10 Conference (TENCON), Singapore, 22–25 November 2016; pp. 1270–1273. [Google Scholar]
  10. Bukar, A.L.; Modu, B.; Gwoma, Z.M.; Mustapha, M.; Buji, A.B.; Lawan, M.B.; Tijjani, I.; Benisheik, U.A.; Bukar, A.; Mai, K.B. Economic Assessment of a PV/Diesel/Battery Hybrid Energy System for a Non-Electrified Remote Village in Nigeria. Eur. J. Eng. Res. Sci. 2017, 2, 21. [Google Scholar] [CrossRef] [Green Version]
  11. Bukar, A.L.; Tan, C.W.; Lau, K.Y. Optimal sizing of an autonomous photovoltaic/wind/battery/diesel generator microgrid using grasshopper optimization algorithm. Sol. Energy 2019, 188, 685–696. [Google Scholar] [CrossRef]
  12. Nadjemi, O.; Nacer, T.; Hamidat, A.; Salhi, H. Optimal hybrid PV/wind energy system sizing: Application of cuckoo search algorithm for Algerian dairy farms. Renew. Sustain. Energy Rev. 2017, 70, 1352–1365. [Google Scholar] [CrossRef]
  13. Al-falahi, M.D.A.; Jayasinghe, S.D.G.; Enshaei, H. A review on recent size optimization methodologies for standalone solar and wind hybrid renewable energy system. Energy Convers. Manag. 2017, 143, 252–274. [Google Scholar] [CrossRef]
  14. Khatib, T.; Ibrahim, I.A.; Mohamed, A. A review on sizing methodologies of photovoltaic array and storage battery in a standalone photovoltaic system. Energy Convers. Manag. 2016, 120, 430–448. [Google Scholar] [CrossRef]
  15. Sinha, S.; Chandel, S.S. Review of recent trends in optimization techniques for solar photovoltaic-wind based hybrid energy systems. Renew. Sustain. Energy Rev. 2015, 50, 755–769. [Google Scholar] [CrossRef]
  16. Anoune, K.; Bouya, M.; Astito, A.; Abdellah, A. Ben Sizing methods and optimization techniques for PV-wind based hybrid renewable energy system: A review. Renew. Sustain. Energy Rev. 2018, 93, 652–673. [Google Scholar] [CrossRef]
  17. Bashir, N.; Modu, B. Techno-economic analysis of off-grid renewable energy systems for rural electrification in Northeastern Nigeria. Int. J. Renew. Energy Res. 2018, 8, 1217–1228. [Google Scholar]
  18. Ajlan, A.; Tan, C.W.; Abdilahi, A.M. Assessment of environmental and economic perspectives for renewable-based hybrid power system in Yemen. Renew. Sustain. Energy Rev. 2017, 75, 559–570. [Google Scholar] [CrossRef]
  19. Sinha, S.; Chandel, S.S. Review of software tools for hybrid renewable energy systems. Renew. Sustain. Energy Rev. 2014, 32, 192–205. [Google Scholar] [CrossRef]
  20. Bernal-Agustín, J.L.; Dufo-López, R. Simulation and optimization of stand-alone hybrid renewable energy systems. Renew. Sustain. Energy Rev. 2009, 13, 2111–2118. [Google Scholar] [CrossRef]
  21. Al Garni, H.Z.; Abubakar Mas’ud, A.; Wright, D. Design and economic assessment of alternative renewable energy systems using capital cost projections: A case study for Saudi Arabia. Sustain. Energy Technol. Assess. 2021, 48, 101675. [Google Scholar] [CrossRef]
  22. Seedahmed, M.M.A.; Ramli, M.A.M.; Bouchekara, H.R.E.H.; Milyani, A.H.; Rawa, M.; Nur Budiman, F.; Firmansyah Muktiadji, R.; Mahboob Ul Hassan, S. Optimal sizing of grid-connected photovoltaic system for a large commercial load in Saudi Arabia. Alex. Eng. J. 2022, 61, 6523–6540. [Google Scholar] [CrossRef]
  23. Thirunavukkarasu, M.; Sawle, Y. A Comparative Study of the Optimal Sizing and Management of Off-Grid Solar/Wind/Diesel and Battery Energy Systems for Remote Areas. Front. Energy Res. 2021, 9, 752043. [Google Scholar] [CrossRef]
  24. Tozzi, P.; Jo, J.H. A comparative analysis of renewable energy simulation tools: Performance simulation model vs. system optimization. Renew. Sustain. Energy Rev. 2017, 80, 390–398. [Google Scholar] [CrossRef]
  25. Kazem, H.A.; Khatib, T.; Sopian, K. Sizing of a standalone photovoltaic/battery system at minimum cost for remote housing electrification in Sohar, Oman. Energy Build. 2013, 61, 108–115. [Google Scholar] [CrossRef]
  26. Ayop, R.; Isa, N.M.; Tan, C.W. Components sizing of photovoltaic stand-alone system based on loss of power supply probability. Renew. Sustain. Energy Rev. 2018, 81, 2731–2743. [Google Scholar] [CrossRef]
  27. Abdul Aziz, N.I.; Sulaiman, S.I.; Shaari, S.; Musirin, I.; Sopian, K. Optimal sizing of stand-alone photovoltaic system by minimizing the loss of power supply probability. Sol. Energy 2017, 150, 220–228. [Google Scholar] [CrossRef]
  28. Sanajaoba, S.; Fernandez, E. Maiden application of Cuckoo Search algorithm for optimal sizing of a remote hybrid renewable energy System. Renew. Energy 2016, 96, 1–10. [Google Scholar] [CrossRef]
  29. Zhang, X.; Tan, S.C.; Li, G.; Li, J.; Feng, Z. Components sizing of hybrid energy systems via the optimization of power dispatch simulations. Energy 2013, 52, 165–172. [Google Scholar] [CrossRef]
  30. Khatod, D.K.; Pant, V.; Sharma, J. Analytical approach for well-being assessment of small autonomous power systems with solar and wind energy sources. IEEE Trans. Energy Convers. 2010, 25, 535–545. [Google Scholar] [CrossRef]
  31. Lujano-Rojas, J.M.; Dufo-López, R.; Bernal-Agustín, J.L. Probabilistic modelling and analysis of stand-alone hybrid power systems. Energy 2013, 63, 19–27. [Google Scholar] [CrossRef]
  32. Fathy, A.; Kaaniche, K.; Alanazi, T.M. Recent Approach Based Social Spider Optimizer for Optimal Sizing of Hybrid PV/Wind/Battery/Diesel Integrated Microgrid in Aljouf Region. IEEE Access 2020, 8, 57630–57645. [Google Scholar] [CrossRef]
  33. Farh, H.M.H.; Al-Shamma, A.A.; Al-shaalan, A.M.; Alkuhayli, A. Technical and Economic Evaluation for Off-Grid Hybrid Renewable Energy System Using Novel Bonobo Optimizer. Sustainability 2022, 14, 1533. [Google Scholar] [CrossRef]
  34. Jufri, F.H.; Aryani, D.R.; Garniwa, I.; Sudiarto, B. Optimal battery energy storage dispatch strategy for small-scale isolated hybrid renewable energy system with different load profile patterns. Energies 2021, 14, 3139. [Google Scholar] [CrossRef]
  35. Ramli, M.A.M.; Bouchekara, H.R.E.H.; Alghamdi, A.S. Optimal sizing of PV/wind/diesel hybrid microgrid system using multi-objective self-adaptive differential evolution algorithm. Renew. Energy 2018, 121, 400–411. [Google Scholar] [CrossRef]
  36. Bouchekara, H.R.E.H.; Javaid, M.S.; Shaaban, Y.A.; Shahriar, M.S.; Ramli, M.A.M.; Latreche, Y. Decomposition based multiobjective evolutionary algorithm for PV/Wind/Diesel Hybrid Microgrid System design considering load uncertainty. Energy Rep. 2021, 7, 52–69. [Google Scholar] [CrossRef]
  37. Bukar, A.L.; Tan, C.W.; Yiew, L.K.; Ayop, R.; Tan, W.S. A rule-based energy management scheme for long-term optimal capacity planning of grid-independent microgrid optimized by multi-objective grasshopper optimization algorithm. Energy Convers. Manag. 2020, 221, 113161. [Google Scholar] [CrossRef]
  38. Omar, A.S. Multi-Objective Optimization of a Stand-alone Hybrid PV/wind/battery/diesel Micro-grid. In Proceedings of the IEEE Conference on Power Electronics and Renewable Energy, CPERE 2019, Aswan, Egypt, 23–25 October 2019; pp. 391–396. [Google Scholar]
  39. Bouchekara, H.R.E.H.; Shahriar, M.S.; Irshad, U.B.; Sha’aban, Y.A.; Mahmud, M.A.P.; Javaid, M.S.; Ramli, M.A.M.; Farjana, S.H. Optimal sizing of hybrid photovoltaic/diesel/battery nanogrid using a parallel multiobjective PSO-based approach: Application to desert camping in Hafr Al-Batin city in Saudi Arabia. Energy Rep. 2021, 7, 4360–4375. [Google Scholar] [CrossRef]
  40. Zhu, W.; Guo, J.; Zhao, G.; Zeng, B. Optimal sizing of an island hybrid microgrid based on improved multi-objective grey wolf optimizer. Processes 2020, 8, 1581. [Google Scholar] [CrossRef]
  41. Fathy, A. Reliable and efficient approach for mitigating the shading effect onphotovoltaic module based on Modified Artificial Bee Colony algorithm. Renew. Energy 2015, 81, 78–88. [Google Scholar] [CrossRef]
  42. Ayop, R.; Tan, C.W. A comprehensive review on photovoltaic emulator. Renew. Sustain. Energy Rev. 2017, 80, 430–452. [Google Scholar] [CrossRef]
  43. Huld, T.; Gottschalg, R.; Beyer, H.G.; Topič, M. Mapping the performance of PV modules, effects of module type and data averaging. Sol. Energy 2010, 84, 324–338. [Google Scholar] [CrossRef]
  44. Daud, A.K.; Ismail, M.S. Design of isolated hybrid systems minimizing costs and pollutant emissions. Renew. Energy 2012, 44, 215–224. [Google Scholar] [CrossRef]
  45. Hamrouni, N.; Jraidi, M. Chérif, a Solar radiation and ambient temperature effects on the performances of a PV pumping system. Rev. Energ. Renouvelables 2008, 11, 95–106. [Google Scholar]
  46. Sukamongkol, Y.; Chungpaibulpatana, S.; Ongsakul, W. A simulation model for predicting the performance of a solar photovoltaic system with alternating current loads. Renew. Energy 2002, 27, 237–258. [Google Scholar] [CrossRef]
  47. Wang, L.; Tan, A.C.C.; Cholette, M.; Gu, Y. Comparison of the effectiveness of analytical wake models for wind farm with constant and variable hub heights. Energy Convers. Manag. 2016, 124, 189–202. [Google Scholar] [CrossRef] [Green Version]
  48. Justus, C.G. Wind energy statistics for large arrays of wind turbines (New England and Central U.S. Regions). Sol. Energy 1978, 20, 379–386. [Google Scholar] [CrossRef]
  49. Rehman, S.; Al-Abbadi, N.M. Wind shear coefficients and energy yield for Dhahran, Saudi Arabia. Renew. Energy 2007, 32, 738–749. [Google Scholar] [CrossRef]
  50. Jaramillo, O.A.; Borja, M.A. Wind speed analysis in La Ventosa, Mexico: A bimodal probability distribution case. Renew. Energy 2004, 29, 1613–1630. [Google Scholar] [CrossRef]
  51. Farrugia, R.N. The wind shear exponent in a Mediterranean island climate. Renew. Energy 2003, 28, 647–653. [Google Scholar] [CrossRef]
  52. IES 2005 Wind Turbines-Part 1: Design Requirements. In International Standard; International Electrotechnical Commission: 2005; Volume 2005, p. 88. Available online: https://www.iec.ch/dyn/www/f?p=103:7:508880294409522::::FSP_ORG_ID:1282 (accessed on 6 April 2022).
  53. Patel, M.R. Book Review: Wind and Solar Power Systems—Design, Analysis, and Operation; CRC Press: Boca Raton, FL, USA, 2006; Volume 30, ISBN 0849316057. [Google Scholar]
  54. Wang, L.; Singh, C. PSO-based multi-criteria optimum design of a grid-connected hybrid power system with multiple renewable sources of energy. In Proceedings of the 2007 IEEE Swarm Intelligence Symposium, SIS 2007, Honolulu, HI, USA, 1–5 April 2007; pp. 250–257. [Google Scholar] [CrossRef]
  55. Guangqian, D.; Bekhrad, K.; Azarikhah, P.; Maleki, A. A hybrid algorithm based optimization on modeling of grid independent biodiesel-based hybrid solar/wind systems. Renew. Energy 2018, 122, 551–560. [Google Scholar] [CrossRef]
  56. Kharrich, M.; Kamel, S.; Alghamdi, A.S.; Eid, A.; Mosaad, M.I.; Akherraz, M.; Abdel-Akher, M. Optimal design of an isolated hybrid microgrid for enhanced deployment of renewable energy sources in Saudi Arabia. Sustainability 2021, 13, 4708. [Google Scholar] [CrossRef]
  57. Mahmoud, M.M.; Ibrik, I.H. Techno-economic feasibility of energy supply of remote villages in Palestine by PV-systems, diesel generators and electric grid. Renew. Sustain. Energy Rev. 2006, 10, 128–138. [Google Scholar] [CrossRef]
  58. Jayachandran, M.; Ravi, G. Design and Optimization of Hybrid Micro-Grid System. Energy Procedia 2017, 117, 95–103. [Google Scholar] [CrossRef]
  59. Satish Kumar Ramoji, B.J.K. Optimal Economical Sizing of A PV-Wind Hybrid Energy System Using Genetic Algorithm and Teaching Learning Based Optimization. Int. J. Green Energy 2011, 8, 25–43. [Google Scholar]
  60. Hatata, A.Y.; Osman, G.; Aladl, M.M. An optimization method for sizing a solar/wind/battery hybrid power system based on the artificial immune system. Sustain. Energy Technol. Assess. 2018, 27, 83–93. [Google Scholar] [CrossRef]
  61. Skarstein, O.; Uhlen, K. Design considerations with respect to long-term diesel saving in wind/diesel plants. Wind Eng. 1989, 13, 72–87. [Google Scholar]
  62. Azoumah, Y.; Yamegueu, D.; Ginies, P.; Coulibaly, Y.; Girard, P. Sustainable electricity generation for rural and peri-urban populations of sub-Saharan Africa: The “flexy-energy” concept. Energy Policy 2011, 39, 131–141. [Google Scholar] [CrossRef]
  63. Deshmukh, M.K.; Deshmukh, S.S. Modeling of hybrid renewable energy systems. Renew. Sustain. Energy Rev. 2008, 12, 235–249. [Google Scholar] [CrossRef]
  64. Diaf, S.; Diaf, D.; Belhamel, M.; Haddadi, M.; Louche, A. A methodology for optimal sizing of autonomous hybrid PV/wind system. Energy Policy 2007, 35, 5708–5718. [Google Scholar] [CrossRef] [Green Version]
  65. Darras, C.; Sailler, S.; Thibault, C.; Muselli, M.; Poggi, P.; Hoguet, J.C.; Melscoet, S.; Pinton, E.; Grehant, S.; Gailly, F.; et al. Sizing of photovoltaic system coupled with hydrogen/oxygen storage based on the ORIENTE model. Int. J. Hydrogen Energy 2010, 35, 3322–3332. [Google Scholar] [CrossRef]
  66. MERRA-2 Re-Analysis. WEB SERVICES—SODA. Available online: http://www.soda-pro.com/web-services/meteodata (accessed on 6 October 2021).
  67. Yang, H.; Zhou, W.; Lu, L.; Fang, Z. Optimal sizing method for stand-alone hybrid solar-wind system with LPSP technology by using genetic algorithm. Sol. Energy 2008, 82, 354–367. [Google Scholar] [CrossRef]
  68. Rajkumar, R.K.; Ramachandaramurthy, V.K.; Yong, B.L.; Chia, D.B. Techno-economical optimization of hybrid pv/wind/battery system using Neuro-Fuzzy. Energy 2011, 36, 5148–5153. [Google Scholar] [CrossRef]
  69. Kaabeche, A.; Belhamel, M.; Ibtiouen, R. Techno-economic valuation and optimization of integrated photovoltaic/wind energy conversion system. Sol. Energy 2011, 85, 2407–2420. [Google Scholar] [CrossRef]
  70. Fathy, A. A reliable methodology based on mine blast optimization algorithm for optimal sizing of hybrid PV-wind-FC system for remote area in Egypt. Renew. Energy 2016, 95, 367–380. [Google Scholar] [CrossRef]
  71. Kashefi Kaviani, A.; Riahy, G.H.; Kouhsari, S.M. Optimal design of a reliable hydrogen-based stand-alone wind/PV generating system, considering component outages. Renew. Energy 2009, 34, 2380–2390. [Google Scholar] [CrossRef]
  72. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  73. Aljarah, I.; Habib, M.; Faris, H.; Al-Madi, N.; Heidari, A.A.; Mafarja, M.; Elaziz, M.A.; Mirjalili, S. A dynamic locality multi-objective salp swarm algorithm for feature selection. Comput. Ind. Eng. 2020, 147, 106628. [Google Scholar] [CrossRef]
  74. Ibrahim, A.; Mohammed, S.; Ali, H.A.; Hussein, S.E. Breast Cancer Segmentation from Thermal Images Based on Chaotic Salp Swarm Algorithm. IEEE Access 2020, 8, 122121–122134. [Google Scholar] [CrossRef]
  75. Lasmari, A.; Zellagui, M.; El-sehiemy, R.A.; Chenni, R. Multi-objective salp swarm algorithm for performance enhancement of electrical distribution system including DG and DSTATCOM simultaneously. In Proceedings of the 4th International Conference on Artificial Intelligence in Renewable Energetic Systems (IC-AIRES), Tipasa, Algeria, 22–24 November 2020. [Google Scholar]
  76. Shaheen, A.M.; El-Sehiemy, R.A. A Multiobjective Salp Optimization Algorithm for Techno-Economic-Based Performance Enhancement of Distribution Networks. IEEE Syst. J. 2021, 15, 1458–1466. [Google Scholar] [CrossRef]
  77. El Sehiemy, R.A.; Selim, F.; Bentouati, B.; Abido, M.A. A novel multi-objective hybrid particle swarm and salp optimization algorithm for technical-economical-environmental operation in power systems. Energy 2020, 193, 116817. [Google Scholar] [CrossRef]
  78. Wang, J.; Gao, Y.; Chen, X. A novel hybrid interval prediction approach based on modified lower upper bound estimation in combination with multi-objective salp swarm algorithm for short-term load forecasting. Energies 2018, 11, 1561. [Google Scholar] [CrossRef] [Green Version]
  79. Mirjalili, S. Dragonfly algorithm: A new meta-heuristic optimization technique for solving single-objective, discrete, and multi-objective problems. Neural Comput. Appl. 2016, 27, 1053–1073. [Google Scholar] [CrossRef]
  80. Mirjalili, S.Z.; Mirjalili, S.; Saremi, S.; Faris, H.; Aljarah, I. Grasshopper optimization algorithm for multi-objective optimization problems. Appl. Intell. 2018, 48, 805–820. [Google Scholar] [CrossRef]
  81. Mirjalili, S.; Jangir, P.; Saremi, S. Multi-objective ant lion optimizer: A multi-objective optimization algorithm for solving engineering problems. Appl. Intell. 2017, 46, 79–95. [Google Scholar] [CrossRef]
Figure 1. Stand-alone microgrid system configuration.
Figure 1. Stand-alone microgrid system configuration.
Energies 15 03579 g001
Figure 2. Variation of wind turbine generator power with wind speed.
Figure 2. Variation of wind turbine generator power with wind speed.
Energies 15 03579 g002
Figure 3. The geographical map of the study location using Google Earth. (a) Map of the selected location; (b) close look at the selected location.
Figure 3. The geographical map of the study location using Google Earth. (a) Map of the selected location; (b) close look at the selected location.
Energies 15 03579 g003
Figure 4. Annual solar radiation of the studied site.
Figure 4. Annual solar radiation of the studied site.
Energies 15 03579 g004
Figure 5. Annual wind speed of the studied site.
Figure 5. Annual wind speed of the studied site.
Energies 15 03579 g005
Figure 6. The annual ambient temperature of the studied site.
Figure 6. The annual ambient temperature of the studied site.
Energies 15 03579 g006
Figure 7. Load profile in the winter, spring, summer, and autumn seasons.
Figure 7. Load profile in the winter, spring, summer, and autumn seasons.
Energies 15 03579 g007
Figure 8. Hourly load profile for the residential units.
Figure 8. Hourly load profile for the residential units.
Energies 15 03579 g008
Figure 9. The proposed rule-based energy management strategy.
Figure 9. The proposed rule-based energy management strategy.
Energies 15 03579 g009
Figure 10. Swarm of salps. (a) Individual salp; (b) salp chain [74].
Figure 10. Swarm of salps. (a) Individual salp; (b) salp chain [74].
Energies 15 03579 g010
Figure 11. Flowchart of the proposed approach.
Figure 11. Flowchart of the proposed approach.
Energies 15 03579 g011
Figure 12. Pareto front of stand-alone microgrid system considering LPSP and COE. (a) MOSSA; (b) MODA; (c) MOGOA; (d) MOALO.
Figure 12. Pareto front of stand-alone microgrid system considering LPSP and COE. (a) MOSSA; (b) MODA; (c) MOGOA; (d) MOALO.
Energies 15 03579 g012
Figure 13. The annual contribution of energy generation by PV/WT/battery and DG for MOSSA. (a) Solution #1; (b) solution #10; (c) solution #20.
Figure 13. The annual contribution of energy generation by PV/WT/battery and DG for MOSSA. (a) Solution #1; (b) solution #10; (c) solution #20.
Energies 15 03579 g013
Figure 14. The annual contribution of energy generation by PV/WT/battery and DG for MODA. (a) Solution #1; (b) solution #10; (c) solution #20.
Figure 14. The annual contribution of energy generation by PV/WT/battery and DG for MODA. (a) Solution #1; (b) solution #10; (c) solution #20.
Energies 15 03579 g014
Figure 15. The annual contribution of energy generation by PV/WT/battery and DG for MOGOA. (a) Solution #1; (b) solution #10; (c) solution #20.
Figure 15. The annual contribution of energy generation by PV/WT/battery and DG for MOGOA. (a) Solution #1; (b) solution #10; (c) solution #20.
Energies 15 03579 g015
Figure 16. The annual contribution of energy generation by PV/WT/battery and DG for MOALO. (a) Solution #1; (b) solution #10; (c) solution #20.
Figure 16. The annual contribution of energy generation by PV/WT/battery and DG for MOALO. (a) Solution #1; (b) solution #10; (c) solution #20.
Energies 15 03579 g016
Figure 17. Pareto front of stand-alone microgrid system considering LPSP and COE obtained by MOSSA, MODA, MOGOA, and MOALO.
Figure 17. Pareto front of stand-alone microgrid system considering LPSP and COE obtained by MOSSA, MODA, MOGOA, and MOALO.
Energies 15 03579 g017
Figure 18. Power generated by the various components of the HMS over one year for MOSSA solution #1. (a) PV array; (b) WT; (c) DG. The annual profile of the battery SOC (d).
Figure 18. Power generated by the various components of the HMS over one year for MOSSA solution #1. (a) PV array; (b) WT; (c) DG. The annual profile of the battery SOC (d).
Energies 15 03579 g018
Figure 19. Power generation of the hybrid stand-alone microgrid system over one year.
Figure 19. Power generation of the hybrid stand-alone microgrid system over one year.
Energies 15 03579 g019
Figure 20. Hourly variation of the battery SOC during the winter season.
Figure 20. Hourly variation of the battery SOC during the winter season.
Energies 15 03579 g020
Figure 21. Power generated by DG during winter season.
Figure 21. Power generated by DG during winter season.
Energies 15 03579 g021
Table 1. Summary of reported methods in optimizing HRES.
Table 1. Summary of reported methods in optimizing HRES.
ReferenceYearLocationObjective FunctionAlgorithm
Bouchekara, H.R.E.H. [36]2021Saudi ArabiaMinimize the Annual LPSP/COEMultiobjective Evolutionary Algorithm
Bukar, A.L. [37]2020NigeriaMinimize the Annual COE/DPSPMultiobjective Grasshopper Optimization Algorithm
Bouchekara, H.R. [39]2021Saudi ArabiaMinimize the Annual COE/LPSPParallel Multiobjective PSO (PMOPSO)
Farh, H.M.H. [33]2022Saudi ArabiaMinimize the total Annualized System Cost (ASC)Bonobo Optimizer
Seedahmed, M.M. [22]2022Saudi ArabiaMinimize the Annual COEHOMER
Thirunavukkarasu, M. [23]2021IndiaMinimize the Annual NPC/COEHOMER
Fathy, A. [32]2020Saudi ArabiaMinimize the Annual COESocial Spider Optimizer
Omar, A.S. [38]2019EgyptMinimize the Annual LPSP/COE and Maximize RFMultiobjective Dragonfly Algorithm (MODA)
Bukar, A.L. [12]2019NigeriaMinimize the Annual COEGrasshopper Optimization Algorithm
Zhu, W. [40]2020ChinaMinimize the Annual CACS/DPSPMultiobjective Grey Wolf Optimizer
Table 2. Geographical coordinates of the study site.
Table 2. Geographical coordinates of the study site.
ParticularsDetails
CountryAlgeria
StateDjelfa
DistrictAïn El Ibel
MunicipalityAïn El Ibel
Latitude34.346°
Longitude3.163°
Altitude above sea level1098 m
Study siteCentral PV Aïn El Ibel (SKTM)
Period of measurement1 January 2020–31 December 2020
Table 3. Power-consuming appliances needed for the residential units.
Table 3. Power-consuming appliances needed for the residential units.
AppliancesPower (W)QuantityElectric Load (W)
Refrigerator2202440
Television1503450
Mobile Charger12672
Water Pump4502900
Radio12112
Lamps Bulb755375
Lamps CFL188144
Fluorescent Light405200
Laptop463138
Desktop computer1202240
Mixer4501450
Deep freezer2601260
Air conditioner4302860
Washing machine4201420
Microwave9001900
Table 4. Technical and economic specifications of the HMGS components [6,32,36].
Table 4. Technical and economic specifications of the HMGS components [6,32,36].
ComponentParameterValueUnit
PhotovoltaicRated power
PV regulator efficiency
Lifetime
PV regulator cost
Initial cost
7.3
95
24
1500
2150
kW
%
Year
$
$/kW
Wind TurbineModel
Rated power
Cut-in wind speed
Rated wind speed
Cut-out wind speed
Number of blades
Tower height
Efficiency
Lifetime
Wind turbine regulator cost
Price
Eolica 2 kW
2
2.0
9.0
20.0
3
20
95
24
1000
2000
kW
m/s
m/s
m/s
m
%
Year
$
$/kW
BatteryRated power
Efficiency
Lifetime
SOC_Min
SOC_Max
DOD
Initial cost
40
85
2
30
100
70
220
kW h
%
Year
%
%
%
USD/kW h
Diesel GeneratorRated power
Lifetime
Initial cost
4
24,000
1000
kW
hours
USD/kW h
InverterLifetime
Efficiency
Initial cost
24922500Year
%
$
Economic ParametersProject lifetime
Fuel inflation rate
O&M + Running cost
Real interest
Discount rate
24
5
20
13
8
Year
%
%
%
%
Table 5. Parameters for each approach.
Table 5. Parameters for each approach.
AlgorithmsParameters
MOSSAPopulation size: 100
Number of iterations: 100
Archive size: 200
The   coefficient   c 1 : Equation (28)
The   coefficient   c 2 : rand
The   coefficient   c 3 :   rand ,   c 3 < 0.5
MODAPopulation size: 100
Number of iterations: 100
Archive size: 200
w = 0.9 0.2 ,   s = 0.1
a = 0.1 ,   c = 0.7
f = 1 ,   e = 1
MOGOAPopulation size: 100
Number of iterations: 100
Archive size: 200
The   coefficient   c M a x   : 1
The   coefficient   c M i n   : 0.00004
MOALOPopulation size: 100
Number of iterations: 100
Archive size: 200
c t = c t I ,   d t = d t I
Table 6. Selected solutions from the Pareto front of MOSSA.
Table 6. Selected solutions from the Pareto front of MOSSA.
Solution #PV (KW)NADNWTCOE (USD/kW h)LPSP (%)RE (%)The AEC from PV (MW)The AEC from WT (MW)The AEC from BT (MW)The AEC from DG (MW)
Solution # 165.8833100.25527.07990.46121.838.87118.3315.33
Solution # 264.7153100.25227.11890.31119.6438.87118.1815.36
Solution # 364.4363100.25127.15390.27119.1238.87118.0915.36
Solution # 463.7342.953100.24927.21390.17117.8238.87116.2115.4
Solution # 563.7732.979100.24927.21590.17117.938.87117.1815.4
Solution # 651.4542.721100.21427.85988.1895.1238.87105.315.84
Solution # 746.6422.781100.20128.10587.286.2338.87106.6116.01
Solution # 844.6383100.19628.12686.882.5238.87114.3416.02
Solution # 942.8363100.19128.25286.3679.1938.87113.7416.11
Solution # 1041.7262.711100.18728.42485.977.1438.87102.7416.36
Solution # 1132.4142.697100.16129.35382.4959.9238.8798.0717.22
Solution # 1230.9782.844100.15829.38282.0457.2738.87102.3517.27
Solution # 1329.0813100.15229.5781.1253.7638.87106.3817.49
Solution # 1426.4932.766100.14530.02879.4148.9838.8796.3518.09
Solution # 1525.1322.733100.14130.29478.4646.4638.8794.2118.38
Solution # 1624.9352.711100.13930.4877.9446.138.8792.9518.55
Solution # 1723.832.887100.13630.71577.1744.0538.8797.618.7
Solution # 1823.592.697100.13530.96276.6143.6138.8790.9818.9
Solution # 1923.8431.169100.13431.02675.6344.0838.8741.519.87
Solution # 2024.093190.13331.12175.2444.5436.5835.8320.09
AEC: Annual Energy Contribution.
Table 7. Selected solutions from the Pareto front of MODA.
Table 7. Selected solutions from the Pareto front of MODA.
Solution #PV (KW)NADNWTCOE (USD/kW h)LPSP (%)RE (%)The AEC from PV (MW)The AEC from WT (MW)The AEC from BT (MW)The AEC from DG (MW)
Solution #178.0571.215100.28627.14990.99144.338.8749.2416.24
Solution #276.8111.119100.28327.15490.9914238.8745.4816.3
Solution #376.4991.137100.28227.15590.96141.4238.8746.1516.29
Solution #475.3771.478100.2827.15890.99139.3538.8759.5216.07
Solution #574.6791.478100.27827.17790.91138.0638.8759.4916.08
Solution #673.3991.437100.27427.21590.76135.6938.8757.816.13
Solution #773.3681.452100.27427.22390.77135.6338.8758.4116.11
Solution #874.2821.13590.27327.64290.47137.3236.6545.6916.53
Solution #954.7342.271100.22128.13288.41101.1938.8788.2516.07
Solution #1044.1951.926100.19328.42586.1281.738.8774.0216.74
Solution #1143.871.882100.19228.49586.0181.138.8772.3216.77
Solution #1234.0752.6100.16529.18483.1662.9938.8795.4717.07
Solution #1333.131100.1629.52381.5161.2538.8737.8518.38
Solution #1429.6551.34290.14930.25779.6254.8236.6548.9418.64
Solution #1528.7431.43690.14630.48278.9353.1436.6551.8618.8
Solution #1628.7171.490.14630.49678.953.0936.6550.5918.82
Solution #1728.6831.4790.14630.50278.8853.0336.655318.79
Solution #1828.8531.38390.14430.71178.4853.3436.6549.7718.98
Solution #1925.706180.13431.34775.447.5233.8635.820.08
Solution #2025.6371.00380.13431.41175.2547.3933.8635.8320.11
AEC: annual energy contribution.
Table 8. Selected solutions from the Pareto front of MOGOA.
Table 8. Selected solutions from the Pareto front of MOGOA.
Solution #PV (KW)NADNWTCOE (USD/kW h)LPSP (%)RE (%)The AEC from PV (MW)The AEC from WT (MW)The AEC from BT (MW)The AEC from DG (MW)
Solution #185.8251.014100.30727.02191.74158.6638.2641.4716.26
Solution #274.9891.007100.27727.33290.68138.6338.2640.9116.47
Solution #366.4132.695100.25427.44290.32122.7838.26106.1015.53
Solution #464.1882.098100.24727.82589.81118.6638.2682.6315.91
Solution #564.1882.098100.24727.83189.80118.6638.2682.6215.91
Solution #664.3092.72190.24428.11989.64118.8934.06105.8015.85
Solution #763.9192.68690.24128.35689.43118.1734.06104.1715.99
Solution #850.8413.00090.20828.45987.5693.9934.06114.3416.03
Solution #948.9442.07980.19829.27586.1490.4832.2079.0217.01
Solution #1048.4461.11780.19629.30485.5589.5632.2043.2317.68
Solution #1143.2471.13090.18429.35084.5779.9534.0643.3817.69
Solution #1243.3261.11590.18429.38784.4980.1034.0642.8417.75
Solution #1337.0801.89790.16729.81182.9268.5534.0670.5117.60
Solution #1426.7331.005100.14330.33177.8849.4238.2636.7919.39
Solution #1525.4871.11790.13431.37875.5347.1234.0639.6919.94
Solution #1620.1611.10290.11932.60070.4237.2734.0637.5121.14
Solution #1716.4901.13180.10534.57063.4630.4832.2036.3422.50
Solution #1812.9202.68290.10234.95061.8423.8934.0676.8922.47
Solution #1915.1441.95070.10135.19958.4128.0027.2156.9822.96
Solution #2012.5651.24260.08438.81045.9423.2323.3234.4425.17
AEC: annual energy contribution.
Table 9. Selected solutions from the Pareto front of MOALO.
Table 9. Selected solutions from the Pareto front of MOALO.
Solution #PV (KW)NADNWTCOE (USD/kW h)LPSP (%)RE (%)The AEC from PV (MW)The AEC from WT (MW)The AEC from BT (MW)The AEC from DG (MW)
Solution #178.8132.39100.2927.01891.67145.738.8795.5615.38
Solution #275.3171.953100.2827.22891.19139.2438.8778.0915.69
Solution #364.592.666100.25127.33990.21119.4138.87105.1215.49
Solution #462.961.842100.24527.5489.64116.3938.8772.9816.08
Solution #551.3872.629100.21427.88288.139538.87101.815.89
Solution #646.7351.629100.228.35886.5686.438.8763.1516.81
Solution #742.9911100.18928.60285.1979.4838.8739.1517.52
Solution #838.0521.555100.17429.26283.8570.3538.8758.8317.43
Solution #935.0451.644100.16629.41482.8664.7938.8761.5317.61
Solution #1033.1121.615100.16129.56182.0861.2138.8759.9317.81
Solution #1129.1212.008100.15129.82980.5353.8438.8772.3918.05
Solution #1228.1851100.14730.00879.0252.1138.8737.0719.09
Solution #1324.5261.532100.13531.07276.3745.3438.8753.6119.43
Solution #1420.0761.517100.12332.05372.4437.1138.8751.1320.44
Solution #1515.9221.535100.11233.34467.5629.4338.8749.3321.57
Solution #1616.9571.4980.10734.32664.9631.3531.5347.2422.03
Solution #1712.1171.39590.09735.74158.9322.434.5541.2923.39
Solution #1812.169380.09536.58956.1222.531.5381.1823.37
Solution #1912.474170.08937.44351.3123.0627.5229.3924.63
Solution #2013.9731.10360.08738.10649.5525.8322.9431.8924.61
AEC: annual energy contribution.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Belboul, Z.; Toual, B.; Kouzou, A.; Mokrani, L.; Bensalem, A.; Kennel, R.; Abdelrahem, M. Multiobjective Optimization of a Hybrid PV/Wind/Battery/Diesel Generator System Integrated in Microgrid: A Case Study in Djelfa, Algeria. Energies 2022, 15, 3579. https://0-doi-org.brum.beds.ac.uk/10.3390/en15103579

AMA Style

Belboul Z, Toual B, Kouzou A, Mokrani L, Bensalem A, Kennel R, Abdelrahem M. Multiobjective Optimization of a Hybrid PV/Wind/Battery/Diesel Generator System Integrated in Microgrid: A Case Study in Djelfa, Algeria. Energies. 2022; 15(10):3579. https://0-doi-org.brum.beds.ac.uk/10.3390/en15103579

Chicago/Turabian Style

Belboul, Zakaria, Belgacem Toual, Abdellah Kouzou, Lakhdar Mokrani, Abderrahman Bensalem, Ralph Kennel, and Mohamed Abdelrahem. 2022. "Multiobjective Optimization of a Hybrid PV/Wind/Battery/Diesel Generator System Integrated in Microgrid: A Case Study in Djelfa, Algeria" Energies 15, no. 10: 3579. https://0-doi-org.brum.beds.ac.uk/10.3390/en15103579

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop