1. Introduction
A robust power grid is a vital element that ensures the sustained supply of electrical power to all customers connected to the grid, with low cost, better quality, and minimal environmental harm [
1]. Establishing such a grid is very challenging in a centralized system because of long-distance transmission lines, which increase the outage threat. Moreover, factors like high carbon emission levels, high transmission losses, and the challenge of entertaining remote consumers have called reliance on an extensive grid system into question. The distributed generation (DG) concept has provided an alternative to such a situation; it deregulates highly interconnected power systems and takes care of the issues described above. This approach could significantly reduce transmission loss and increase system efficiency [
2,
3]. The microgrid concept is one of the DG implementation outcomes, which aims to ensure the optimal power supply to connected loads [
4]. Integration of renewable energy sources like solar and wind in microgrids ensures the mitigation of environmental effects. It also helps control the power flow and eases connection with the main grid [
5]. A nanogrid is essentially a scaled-down version of a microgrid, and has received significant attention from researchers for its promising prospects in power system engineering.
The nanogrid is a single electrical power distribution domain in terms of capacity, voltage profile, expense, control strategy, and overall administration. This distribution system is usually used for one consumer: a single-unit house, a small-scale load or a small building [
6]. The factors which make a nanogrid divergent from other minigrids are: consumer size, power rating, load size, complexity, hardware configuration, and control strategy. The margin of these factors between a nanogrid and other grids is still debated and nicely presented by Burmester et al. [
6]. In a nanogrid system, there must be a local unit of power production. It can be a renewable source or a fossil fuel-dependent one. A hybrid nanogrid (HN) is a single unit of nanogrid in which different power production sources are present simultaneously. Such a unit is often connected to other neighboring grids by a gateway channel. Any sort of storage device is used in the nanogrid to ensure the stability of the system. Important technoeconomical parameters—cost of energy, storage cost, annual life cycle cost etc.—directly influence the storage devices chosen [
7]. A controller, also known as the brain of a nanogrid, must be present, in order to coordinate the necessary operational strategies [
8].
The classification of a nanogrid can be made via several factors. One of them is the type of transferred energy in the network—DC or AC. DC distribution offers better efficiency than AC despite some limitations [
9]. Among the sources of energy, both renewable (solar, wind, biogas) and traditional (diesel, fuel cell) forms can be used [
10]. However, to ensure the desired energy flow, suitable converters are widely used in both AC and DC grid systems.
Several methods are presented and discussed in the literature to optimize the size and operation of the nanogrid. A mixed-integer linear programming (MILP)-based algorithm was proposed by Atia et al. [
11], which successfully optimized a residential load-connected hybrid nanogrid. This approach has a significant advantage: easy calculation. A hybrid grid system was optimized with a MILP-based multi-objective optimization problem [
12]. Both the renewable and traditional sources were presented, and the fuzzy model searched for the optimal solution. Moreover, this paper’s problem formulation included the cost-emission factor to reduce energy cost and CO
2 emission. For a small power system to support a ship, a similar approach was addressed by Lan et al. [
13], including the optimal navigation of route tracking. Lokeshgupta et al. [
14] offered an optimal energy management system, where a multi-objective optimization strategy simultaneously reduced the energy bill and the peak load. A decision support mechanism was proposed by Li et al. [
15] to optimize the nanogrid, where the problem formulation was made by an augmented–constrained method. The characteristics and behavior of every unit of a hybrid network were well discussed by Hosseinalizadeh et al. [
16]. The formulation complexity and optimization time were significantly reduced by linear programming formulation [
17] while dealing with many variables. In a recent article [
18], the branch-and-cut approach was proposed to search the Pareto front (PF) solution of the multi-objective problem. The paper addressed the issue of real-time variation of energy demand in a nanogrid.
A multi-objective approach was proposed by Brandoni et al. [
18], where different load demand scenarios are considered in the load demand. The matter of climate classification was used to optimize a grid-connected system by Shivam et al. [
19], and the authors investigated the issue in four different Taiwanese regions. A study was carried out in Savannah, Georgia, USA in which the impact of grid-connected nanogrids in tropical climates was discussed [
20]. The article [
21] presented a geographical map for optimally installing the standalone and grid-connected systems. It summarized that Cambridge Bay’s location (in Canada) was the most suitable place from an economic point of view, whereas Toamasina was the optimal location in terms of energy efficiency. Voltage quality improvements [
22] in a nanogrid system and cybersecurity insurance [
23] were also addressed in some recent articles. The review article [
24] gathered all the relevant works on nanogrid optimization and demonstrated a critical comparison among the attempted approaches.
Camping is a part of tradition and culture in all the gulf countries; the Kingdom of Saudi Arabia (KSA) is not an exception. Camps are established in the desert, as a large portion of the countryside is covered by bare desert land. Al-Falahi et al. and Hilden et al. [
24,
25] demonstrated life in a desert tent, along with the life of tribal citizens of KSA. Some tents are built for a limited period, mostly during the winter season. In contrast, some are built for an extended period, mainly to graze domestic livestock. The temporary tents are primarily established for recreation, where local people spend time with their family and friends in cold weather. This paper used Hafr Al-Batin city, located in the Eastern Province of Saudi Arabia, to study the effect of nanogrids in a desert camp. Many tent-establishing shops in the city enjoy patronage in winter, selling and setting up hundreds of tents in the desert. However, current practices utilize diesel generators in these tents, which has evident shortcomings, e.g., high cost and negative environmental impact. These tents can be considered nanogrids. Upgrading such energy-consuming units by integrating a renewable energy source has vast potential. This upgrade is what we propose in this paper. This approach is suitable for any city or country, where desert camps are present. This paper presents an optimization strategy for a single nanogrid and proposes combining multiple nanogrids to form a microgrid—and operate them optimally. Even though much research has already been conducted on grid optimization, this paper will fill the gap by addressing the practical application in desert camps, where both the concepts of nanogrids and microgrids are concurrently implemented.
A multi-objective evolutionary algorithm is used in this paper to solve the complex multi-objective optimization of the nanogrid model. This algorithm offers an acceptable balance between diversity and convergence [
26]. Several optimization techniques are noted in the literature, with a shortcoming: successfully dealing with multi-objective problems is more difficult than accounting for only two or three objectives. However, this paper’s adopted technique is free from such drawbacks and effectively combines the properties of dominance and decomposition. Separately, dominance [
27] and decomposition-based [
28,
29] evolutionary techniques have been used in the literature to deal with similar problems of microgrid characteristic optimization.
The main contribution of this paper is the optimal design of a smart grid, which can be modeled as a nanogrid or microgrid, depending upon the size of the desert camp load. This was achieved as follows: first, we formulated the nanogrid/microgrid design as a multi-objective problem. Second, we solved this problem using a performant multi-objective evolutionary algorithm based on dominance and decomposition—which has not been used before to solve a similar problem. Another contribution of this paper is the adoption of a variety of models for each equipment. The combination of these equipment models will lead to better results than the traditional use of a single model for each equipment.
The remaining paper is organized as follows: the subsequent section presents the fundamentals of the problem formulation. The description of the case study is discussed in
Section 4. The detailed solution strategy of the addressed problem is explained in
Section 5.
Section 6 discusses the attained results with the necessary explanation. Finally,
Section 7 concludes the paper.
5. Problem Solution Approach
The approach proposed in this paper optimized the sizing of hybrid photovoltaic/diesel/battery nanogrids using a multi-objective evolutionary algorithm based on both dominance and decomposition (MOEA/DD) [
26]. To solve the multi-objective problem presented in (11), MOEA/DD decomposed the problem into several scalar subproblems using classic multi-objective optimization approaches, e.g., the weighted sum, weighted Tchebycheff, and boundary intersection methods [
29]. The penalty-based boundary intersection (PBI) approach, known for good performance on a large class of optimization problems, was adopted in this work [
48,
49].
The algorithm for MOEA/DD is presented in Algorithm 2 [
26]. The MOEA/DD algorithm was initialized by generating
initial solutions and their corresponding weights. The parent population was then updated using elite-preserving mechanisms—offspring generated from parents through a mating procedure. In general, MOEA/DD is made up of the initialization, reproduction, and update procedures.
The initialization procedure (Algorithm 3) began with randomly sampling for the initial parent population
from
, uniformly distributed. This step was followed by the identification of the nondomination level structure of
, after which a set of weight vectors was generated before the assignment of the neighborhood. Weight generation methods discussed in [
46,
47] were associated with the explosion of computational complexity and diversity reduction, respectively. As such, the two-layer weight vector generation, proposed by Li et al. [
26] was adopted. First, the sets of weight vectors in the boundary (
and inside (
layers of a simplex were generated such that
. A coordinate transformation was then adopted to reduce the coordinates of weight vectors in the inside layer. This allowed for the evaluation of the jth component of
according to:
where
represents the shrinkage factor and
The set
) was formed by combining B and
The weight vectors,
respectively defined unique subregions in the objective space,
:
where
. is the acute angle between
and
. The neighborhood of each of the weight vectors
was made up of the
closest weight vectors in a Euclidean sense. The fast nondominated sorting method [
50] was used to divide the solutions in
into nondomination levels
. Finally, each solution in
was initially randomly associated with a unique subregion.
The primary role of the reproduction procedure, presented in Algorithm 3 [
26], was to update the parent population after generat1ing offspring solutions. This update was achieved using two main steps: mating selection and variation operation. In the mating selection, some parents are selected for offspring generation from a neighborhood. Each solution is associated with a uniquely weight-specified subregion based on Euclidean distance in the method used. This association allowed for consideration of neighboring solutions from neighboring subregions before randomly selecting mating parents from the whole population when there were no associated solutions. As proposed by Li et al. in [
26], this work adopted the binary crossover (SBX) [
51] and polynomial mutation [
52] for the variation operation. Any other genetic operator could also be used.
The generated offspring was used to update the parent population according to Algorithm 4 [
26]. The algorithm involved identifying the subregion of the offspring solution
and combing with the parent population
to generate a hybrid population
. The nondominated level structure of
was then determined using the method presented by Li et al. [
53]. Other considerations for updating were adopted from [
26]. The optimal values of hypermeters were determined by trial and error aided by experience working with evolutionary algorithms.
Algorithm 2: General Framework of multi-objective evolutionary algorithm based on both dominance and decomposition (MOEA/DD) |
| Output: |
1 | INITIALIZATION;→ // —Parent population, —weight vector set, —neighborhood index set. |
2 | while stopping criterion not fulfilled do |
3 | fortodo |
4 | MATING_SELECTION |
5 | VARIATION |
6 | foreach// is an offspring |
7 | UPDATE_POPULATION |
8 | end |
| end |
| end |
| return. |
Algorithm 3: Mating Selection. |
| Input: | |
| Output: | |
1 | ifthen |
2 | Randomly choose indices from |
3 | If no solution in the selected subregion then |
4 | Randomly choose solutions from to form |
5 | else |
6 | Randomly choose solutions from the selected subregions to form |
7 | end |
8 | else |
9 | Randomly choose solutions from to form |
12 | end |
13 | return |
Algorithm 4. Update Procedure. |
| Input: | |
| Output: | |
1 | Find the subregion associated with according to (6) |
2 | |
3 | Update the nondomination level structure of |
4 | Ifthen |
5 | OCATE_WORST |
6 | |
7 | else |
8 | If then |
| Solution |
9 | If // is the associated subregion of |
10 | |
11 | else // |
12 | |
13 | |
14 | end |
15 | else |
16 | Identify the most crowded subregion associated with those solutions in |
17 | ifthen |
18 | Find the worst solution |
19 | |
20 | else // |
21 | LOCATE_WORST |
22 | |
23 | end |
24 | end |
25 | end |
26 | Update the nondominated level structure of ; |
27 | return P |
6. Application, Results and Discussion
The details of the models used for PV panels, batteries and inverters are tabulated in
Table A1,
Table A2 and
Table A3 of the
Appendix A, respectively. The design variables were constrained as follows:
,
,
,
,
and
.
The proposed approach was applied to the following two case studies using a population size of 100 for 200 iterations:
CASE 1: Nanogrid operation mode
CASE 2: Microgrid operation mode
The following subsections describe the investigated cases and discuss the obtained results with some analysis and recommendations.
6.1. CASE 1: Nanogrid Operation Mode
This case investigated the operation of each tent separately (i.e., as a nanogrid). The MOEA/DD was run to solve the formulated problem for this case, and the obtained Pareto front was plotted in
Figure 7. In this figure, the MOEA/DD generated 23 solutions; all of which occured between (COE = 0.3396
$/kWh, LPSP = 0.3015%) and (COE = 0.5544
$/kWh, LPSP = 0.1946%). It can be seen that the designer/operator can use a variety of solutions for this case.
In addition to that, the obtained PF solutions are given in
Table 2. For convenience, these solutions are sorted in ascending order based on their COEs. If the designer/operator selects solution # 1, the nanogrid will be made of 2 PV panels and two diesel generators and will have an autonomy of 24 hours. The models for PV panel, battery and inverter are model # 7, model # 1 and model # 3, respectively, for this case. Finally, the objective functions obtained for this solution are COE = 0.3396
$/kWh and LPSP = 0.3015%.
A second example can be the selection of solution # 8. The designed nanogrid has an autonomy of more than 26 h. It comprises 13 PV panels (Model #1), two diesel generators, the PV panel, battery, and inverter model are model # 1, mode #5 and model # 3, respectively. This solution yields a COE of 0.4182 $/kWh and an LPSP of 0.2383%.
A third example can be illustrated by solution #23. The sized nanogrid, composed of 30 PV panels and two diesel generators, has more than 25 h of autonomy for this solution. The optimized component models are model #1, model #3 and model #3 for the PV panel, battery, and inverter, respectively. The solution has a COE of 0.5544
$/kWh and an LPSP of 0.1946%. The COE obtained for this solution is the most expensive due to the high number of PV panels used. Furthermore,
Figure 8 exposes each component’s contribution (i.e., PV panels, diesel generator, and battery) to the nanogrid over 50 h for solution #23.
6.2. CASE 2: Microgrid Operation Mode
In this second case study, four tents were connected to form one microgrid. The proposed approach based on MOEA/DD was applied to this case, and the obtained PF is plotted in
Figure 9. The number of obtained solutions was 37, and they are well spread between the extreme points with the following objectives (COE = 0.1933
$/kWh, LPSP = 0.4541%) and (COE = 0.278
$/kWh, LPSP = 0.3028%). This solution offers the designer/operator many options to operate the tents together as a microgrid, instead of 4 separate nanogrids.
Furthermore,
Table 3 tabulates the final set of PF solutions sorted in ascending COE-based order. Many options can be selected for this case. One option could be solution #1, where the resulting microgrid will be composed of 60 PV panels based on model # 1, and 5 diesel generators. The batteries have autonomy for almost one and a half days. The battery model #16 and inverter model #1 are the best options for this solution. The objective functions for this solution are COE = 0.1933
$/kWh and LPSP = 0.4541%.
Similarly, if solution #12 is selected, the microgrid will be composed of 77 PV panels, five diesel generators and the system will have an autonomy of fewer than 30 hours. For this solution scheme, the PV panel, battery, and inverter models are model #1, model #16, and model #1, respectively. The objective functions for this solution are COE = 0.1998 $/kWh and LPSP = 0.4272%.
Another solution could be solution #29. This solution is composed of 98 PV panels of model 1 and 5 diesel generators. The models for battery and inverter are model #16 and model #1, respectively.
It can also be seen from
Table 3 that for almost all the cases the best models for PV panels are model #1, the best model among batteries is model #16, and the best model for the inverter is model #1.
Furthermore, solutions from 30 to 37 are not very acceptable. Although they are part of the PF numerically, they represent nonfeasible solutions.
6.3. Comparative Study
A comparative study with other multi-objective optimization algorithms was conducted to assess the quality of the obtained results. Three well-known algorithms were selected for this study: (1) multi-objective evolutionary algorithm based on decomposition (MOEA/D), (2) novel multi-objective particle swarm optimization (NMPSO) and (3) speed-constrained multi-objective particle swarm optimization (SMPSO).
Figure 10 and
Figure 11 show the comparison of the PFs obtained using the four algorithms for CASE 1 and CASE 2, respectively.
For CASE 1, MOEA/DD, MOEA/D, NMPSO and SMPSO generated a set of 23 solutions, 15 solutions, 11 solutions and 11 solutions, respectively. When these results were compared to each other to determine which solutions were nondominated, it was found that:
Among the 23 solutions generated using MOEA/DD, 16 were nondominated by any other solutions found using the remaining algorithms.
Among the 15 solutions generated using MOEA/D, 5 were nondominated by any other solutions found using the remaining algorithms.
Among the 11 solutions generated using NMPSO, all were dominated by the solutions found using the remaining algorithms.
All 11 solutions generated using SMPSO were nondominated.
For CASE 2, the four used algorithms—MOEA/DD, MOEA/D, NMPSO and SMPSO—generated 37 solutions, 28 solutions, 4 solutions, and 4 solutions, respectively. These results were compared to each other to determine which solutions were nondominated. In this case, (i.e., they are part of the PF), the 37 solutions generated using MOEA/DD were nondominated by any other solutions found using the remaining algorithms. In contrast, all solutions generated using MOEA/D, NMPSO and SMPSO were dominated by those obtained using the MOEA/DD.
Therefore, it can be concluded that the proposed approach, using MOEA/DD, is the most performant one for both operation modes. This conclusion was verified in the microgrid operation mode case (i.e., CASE 2.).
7. Conclusions
This paper presents a dominance and decomposition-based multi-objective evolutionary algorithm to optimize a hybrid nanogrid/microgrid’s size and operation. Hafr Al-Batin, a city of the Eastern Province of Saudi Arabia, was selected for implementation. The reason behind choosing that particular city is the common practice of using desert camps, which can easily be considered nanogrid units. The nanogrid consisted of a solar system, storage batteries, diesel generators, invertor, and load components. Two different modes of operation were investigated in this paper: one with the single nanogrid and the second a combination of nanogrids, coordinated to form a microgrid. In the formulation of the optimization problem, the reliability and the cost of the system were considered. A set of diverse, acceptable, and widespread solutions were obtained from the developed program for both the case studies, which will guide the designer/operator of a nanogrid/microgrid to model such a system. The obtained results offered competitive and practical solutions to the addressed problem using a multi-objective optimization strategy.
However, the addressed issue of desert camping is not a regional problem; rather, it can be implemented in any part of the world. The designed nanogrid/microgrid concept applies to desert camps, but is suitable for any islanded load of a similar kind. Moreover, using such desert camps is widely seen in all the world’s desert areas, particularly in the Gulf countries. Therefore, the paper’s aspect is global, even though it addresses a particular region of the globe.
The presented work can be further continued to investigate possible recovery strategies, subject to the failure of a source like a PV panel or diesel generator. The effect of ageing on associated equipment (PV system, diesel generator, batteries, inverter, etc.) or of equipment modification (due to the connection of a new tent, for example) of load size can be considered future work of this research.