Next Article in Journal
An Enhanced Second-Order Cone Programming-Based Evaluation Method on Maximum Hosting Capacity of Solar Energy in Distribution Systems with Integrated Energy
Previous Article in Journal
Mitigating Adverse Impacts of Increased Electric Vehicle Charging on Distribution Transformers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Optimal Scheduling of a Microgrid Using Oppositional Gradient-Based Grey Wolf Optimizer

by
Arul Rajagopalan
1,
Karthik Nagarajan
2,
Oscar Danilo Montoya
3,4,
Seshathiri Dhanasekaran
5,*,
Inayathullah Abdul Kareem
1,
Angalaeswari Sendraya Perumal
1,
Natrayan Lakshmaiya
6 and
Prabhu Paramasivam
7
1
School of Electrical Engineering, Vellore Institute of Technology, Chennai 600127, Tamil Nadu, India
2
Department of Electrical & Electronics Engineering, Hindustan Institute of Technology & Science, Chennai 601301, Tamil Nadu, India
3
Grupo de Compatibilidad e Interferencia Electromágnetica, Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 110231, Colombia
4
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
5
Department of Computer Science, UiT The Arctic University of Norway, 9037 Tromsø, Norway
6
Department of Mechanical Engineering, Saveetha School of Engineering, SIMATS, Chennai 602107, Tamilnadu, India
7
Department of Mechanical Engineering, College of Engineering and Technology, Mattu University, Mettu 318, Ethiopia
*
Author to whom correspondence should be addressed.
Energies 2022, 15(23), 9024; https://doi.org/10.3390/en15239024
Submission received: 27 October 2022 / Revised: 18 November 2022 / Accepted: 23 November 2022 / Published: 29 November 2022
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
Optimal energy management has become a challenging task to accomplish in today’s advanced energy systems. If energy is managed in the most optimal manner, tremendous societal benefits can be achieved such as improved economy and less environmental pollution. It is possible to operate the microgrids under grid-connected, as well as isolated modes. The authors presented a new optimization algorithm, i.e., Oppositional Gradient-based Grey Wolf Optimizer (OGGWO) in the current study to elucidate the optimal operation in microgrids that is loaded with sustainable, as well as unsustainable energy sources. With the integration of non-Renewable Energy Sources (RES) with microgrids, environmental pollution is reduced. The current study proposes this hybrid algorithm to avoid stagnation and achieve premature convergence. Having been strategized as a bi-objective optimization problem, the ultimate aim of this model’s optimal operation is to cut the costs incurred upon operations and reduce the emission of pollutants in a 24-h scheduling period. In the current study, the authors considered a Micro Turbine (MT) followed by a Wind Turbine (WT), a battery unit and a Fuel Cell (FC) as storage devices. The microgrid was assumed under the grid-connected mode. The authors validated the proposed algorithm upon three different scenarios to establish the former’s efficiency and efficacy. In addition to these, the optimization results attained from the proposed technique were also compared with that of the results from techniques implemented earlier. According to the outcomes, it can be inferred that the presented OGGWO approach outperformed other methods in terms of cost mitigation and pollution reduction.

1. Introduction

1.1. Literature Review

The demand for power has drastically increased in recent years, whereas fossil fuels are depleting at a fast rate. In this scenario, engineers, researchers, and scientists across the globe are looking for promising alternative energy resources that are not only renewable, but also environmentally friendly. Environmental pollution and the resultant global warming are two major concerns raised around the globe today. Fossil fuel-based traditional power plants that use coal, gas and other non-renewable energy sources are identified as the primary environmental pollutants. In this background, hybrid energy systems that combine wind, solar, hydro, and other renewable energy sources have gained attention in the past two decades, owing to its renewable nature and pollution-less energy production practices.
Further, these hybrid energy systems exhibit high efficiency and reliability. However, fluctuation in terms of both quality as well as quantity, as a result of environmental conditions, remains the primary disadvantage in renewable energy sources such as Photovoltaic solar cells, wind energy systems, etc. Since it is challenging to predict the fluctuations in PV and small Wind Turbines (WT) at the time of power generation, it becomes inevitable to integrate these sources with other reliable energy sources such as Fuel Cell (FC), Microturbine (MT), battery storage and so on, to ensure a constant supply of quality power. These hybrid systems can be relied upon as sustainable energy systems. In this background, the microgrid is one of the promising candidates that is not only energy-efficient, but also can integrate multiple number of pilot-level Distributed Energy Resources (DERs), such as microturbine (MT), Wind Turbine (WT), battery units, Fuel Cell (FC), photovoltaic (PV), and other energy sources. Though it is an advanced concept that was introduced in recent times, Microgrid (MG) has gained immense popularity owing to multiple advantages provided by the system. A microgrid can integrate different types of distributed RES to fulfil the local load requirements and connect or disconnect from the utility grid [1,2].
In literature, a PV-based 17-bus Low Voltage microgrid, connected with the grid, was analysed in the literature. This microgrid had a Wind Turbine as well as a Fuel Cell (FC) for the on-site production of hydrogen. The researchers optimally scheduled the storage of energy so as to provide energy at the right time and to ensure reserve provision [3]. In [3], the Harmony Search Algorithm was utilized to schedule the storage of energy during a 24-h horizon to mitigate the costs incurred upon operations. Further, the study also widely leveraged the hydrogen system under all operational scenarios, to balance the uncertainties that tend to arise from loads and renewable energy sources.
In [4], multi-objective PSO algorithm was proposed to elucidate multi-objective energy management problem in a grid connected microgrid. A simulation was carried out in a microgrid test system comprising a wind turbine, solar cell, battery, microturbine, and a diesel generator. Operating cost and pollutant emission were considered as objective functions. In order to investigate the performance of MOPSO, three different scenarios were tested in the simulation. In the study conducted earlier [5], NSGA II algorithm was applied to implement multi-objective optimal operation of a grid connected microgrid test system by considering operation cost and pollutant emission as objective functions. In [5], the uncertainty of the generation capacity of solar and wind power generating units and the uncertainty of load was considered to elucidate the multi-objective energy management problem. In literature [6], the researchers evaluated the bi-level optimal operation of a grid connected microgrid that contains loads, storage devices, and distributed energy resources. The upper level optimal model for distribution network dispatch was solved by using a multi-objective optimization approach by considering the microgrid’s power loss and voltage profile as objective functions. A self-adaptive genetic algorithm was proposed to implement the optimal operation of the microgrid. The lower-level optimal model was developed by applying non-linear programming technique to seek the optimal daily operating scheme of different DGs in the microgrid. The microgrid operating cost was taken as objective function. In the study conducted by Xinhui Lu, the authors developed a multi-objective optimal dispatch model to be applied in the microgrid in grid-connected mode. This model holistically considered the costs incurred upon operations and environmental protection through the microgrid system [7]. The distribution generators, considered in this microgrid system, included Electric Vehicles (EVs), diesel engines, WTs, PV arrays, microturbines, etc. An improved PSO algorithm was developed to overcome the multi-objective optimal dispatch problem. In literature [8], a two-stage stochastic p-robust optimal energy trading management system was proposed to be applied in the microgrid. This system included a microturbine, diesel engine, PV and WT systems. Based on the components involved, the authors implemented a multi-scenario tree approach to develop different scenarios for unsure parameters such as WTs, loads, PV, and market-clearing prices. The discretization of all the probability density functions was accomplished via definite intervals. Then, the authors made use of a differential evolution clustering algorithm to mitigate the number of scenario sets to accomplish the general scenario.
A hybrid algorithm, integrating Gradient Boosting Decision Tree (GBDT) and Sandpiper Optimization Algorithm (SOA), was proposed by Prakash Arumugam [9] to achieve optimal energy management in MG systems under the grid-connected mode. This grid-connected microgrid had MT, WT, battery, and photovoltaic (PV) systems. The objective of the proposed technique was to mitigate the cost incurred upon fuel and other areas during operations and maintenance of grid-connected microgrid system and to reduce the hourly variations in power generation of the microgrid. The researchers listed state-of-charge for storage elements, power demand and renewable energy sources as some of the constraints faced in this research work. In the study conducted earlier [10], the authors proposed an energy management system on the basis of Grasshopper Optimization Algorithm (GOA) to find the optimal power generated by Distributed Generators in microgrid. This is crucial information when it comes to cost mitigation from total power generation. In this study, the researcher applied the proposed unit upon microgrid that had five generating units: industrial, commercial, feeding, and residential loads. Further, the study also made use of the Muddy Soil Fish Optimization Algorithm (MSFOA), a novel algorithm developed on the basis of the foraging pattern of fishes. The main aim of this algorithm is to reduce the overhead incurred upon production and mitigate the costs spent during the import of energy from the grid. Further, the study also considered the system constraints [11]. The researchers implemented the presented energy management system in a microgrid framework that had energy storage units, solar PV, diesel generators, and wind units. Further, numerous case scenarios for instance, isolation, and grid-tied scenarios were investigated under different aspects to prove the potentials of the presented MSFOA approach. In the study conducted earlier [12], a promising energy dispatch strategy was devised by the researchers for both the individual microgrid and the grid-connected ones. These microgrids were also inclusive of battery Energy Storage Systems (BESS), MT, WT, PV, FC, and a diesel generator. In this study, the researchers formulated Microgrid Energy Management (MGEM) mechanism as a mixed-integer linear programming technique. To achieve energy management in the microgrid, a novel multi-objective solution was proposed along with a demand response program. In [13], the authors incorporated a multi-objective Particle Swarm Optimization (PSO) technique to ensure energy resource management and optimal distribution of the resources in the proposed microgrid system. In this study, emission and operation cost were considered as objective functions. The researchers calculated the overall operational cost incurred from the microgrid and pollutant-induced emissions under three different scenarios. The microgrid was analyzed in detail earlier [14] as an energy management system when connected with primary power system. The study took parameters such as WT, PV, and load demand into consideration under probabilistic and deterministic conditions. Further, an Equilibrium Optimizer (EO) was also utilized to overcome the energy management problem. In this study, the authors considered operation cost, voltage deviation, and voltage stability index as objective functions. In the study conducted earlier [15], an optimization framework was proposed with three objectives for microgrid energy management in the case of smart homes and Demand Response (DR) programs. The proposed model was tested using an 83-bus distribution system with 11 microgrids. This study considered load demand variations and uncertainties present in the output power of Renewable Energy Sources (RESs). Further, a max–min fuzzy method was utilized to model the objective function as bi-objective and tri-objective models. In addition to these, the study also considered Peak-to-Average Ratio (PAR), emissions and operations cost as objective functions. Recently several optimization techniques were proposed to eliminate energy management problems in microgrids [16,17]. In [18], ANN-based backtracking search algorithm (ANN-BBSA) and ANN-based binary practical swarm optimization (ANN-BPSO) algorithm was implemented to elucidate multi-objective energy management problem in a virtual power plant, which was realized as a modified IEEE 14-bus system with five identical microgrids. In the study conducted earlier by the current study authors, the Local Energy Management System (LEMS) was proposed and implemented under Generalized Power Prediction Model (GPPM) to overcome the performance uncertainties for varying distributed generations of a microgrid [19]. In [20], the interior search algorithm was applied to solve the economic load dispatch problem in a microgrid which comprises diesel generators, fuel cells, and wind turbine. In this study [20], total operating cost was taken as the objective function. Table 1 shows the contributions made by research investigations in microgrid energy management. This table covers different aspects such as optimization algorithm, test system, Renewable Energy Sources (RESs), and storage technology.

1.2. Research Gap and Objective of the Paper

The authors conducted an extensive review of literature from which knowledgeable insights have been gained with regards to energy management challenges in dynamic systems and different test systems and entities. However, most of the research works conducted earlier emphasized the inclusion of specific multi-objective optimization algorithm to achieve energy management in dynamic test systems. These research studies were conducted on the basis of the Pareto front, utilizing multi-objective techniques and attained energy management. From the review of literature and comparative analysis conducted among two or more than two multi-objective energy management methods, the authors found a research gap that needs to be fulfilled. In the literature [41], the microgrid test considered single and multi-objective optimization to be implemented for the first case study. However, for the second and third case study, only single objective optimization was done for the optimal scheduling of the microgrid. In this current research study, we implemented the single and multi-objective optimization using the proposed optimization algorithm for the optimal scheduling of the microgrid for all the case studies considered. A comparison was made between the optimization results obtained by the proposed method and four other different optimization algorithms coded by us to find the better-compromised solution for less pollutant emissions and minimum generation cost.
In this research paper, we have proposed the Oppositional Gradient-based Grey Wolf Optimizer (OGGWO) for elucidating the multi-objective optimal scheduling of a grid connected microgrid. Opposition based learning is incorporated with the proposed algorithm in order to enhance the solution quality and convergence speed. The incorporation of the opposition-based learning encompasses the instantaneous consideration of an estimate and its equivalent opposite estimate so as to accomplish an enhanced approximation for the present candidate solution. Our proposed optimization algorithm utilized the advantages of the gradient feature that offers valuable information about the solution search space. This incorporation of the gradient feature explores the solution search space more perceptively with the consideration of the gradient direction in the problem search process and aids in attaining the global optimal solution. In order to enhance the exploration and exploitation capabilities of the proposed OGGWO algorithm, Gaussian walk and Levy flight techniques were incorporated, and the results were compared with the results obtained using EDNSGA-II, NSGA-II, GGWO, GWO, CSA, and PSO to prove the effectiveness of the proposed algorithm. GGWO, GWO, CSA, and PSO algorithms were coded by us for comparison purposes. The results of EDNSGA-II and NSGA-II were taken from the literature.

1.3. Contributions

The following is the list of contributions made by current research work for the state-of-the-art energy management system proposed in this study.
i.
At first, OGGWO was implemented as an optimization tool to overcome the problem. Further, the efficiency and robustness of the proposed model were determined and contrasted against GGWO and conventional GWO.
ii.
The authors compared the outcomes of the proposed method under three different scenarios to identify the optimization algorithm that produces the best compromised solution between generation cost and the emission of pollutants.
iii.
Multi-objective optimization was carried out under all the case studies considered so as to find the best-compromised solution for the multi-objective energy management problem.

2. Problem Formulation

The primary objective behind the optimal operations of the microgrid, with Distributed Energy Sources (DERs) that include non-renewable and renewable energy sources, is to identify the optimal operating points of DERs to generate power during the scheduled horizon. The key aim of the optimal operations of the microgrid that is incorporated with DERs such as RESs and non-RESs is to find out the optimal operating points of DERs so that the power can be produced at the time of scheduled horizon. This objective should be met to reduce the cost incurred upon operations and reduce the emission rate in parallel. Thus, the current study aims at achieving two different and conflicting objectives together, i.e., emission reduction and mitigation of cost incurred upon operations.
In the current study, the authors propose a precise mathematical model to achieve short-term energy management so as to mitigate the costs spent upon operations and reduce the emission of pollutants in microgrid operations.

2.1. Objective Functions

As mentioned above, the current study considered the following objective functions such as operating costs and costs incurred upon reducing the pollution emissions.

2.1.1. Mitigation of Operating Cost

In Equation (1), the authors calculate the total operating costs incurred in a microgrid. In this equation, T corresponds to the total time taken, by the study, in terms of hours, whereas N g denote the generation of energy and N s correspond to storage units. Here, U i ( t )   corresponds to the status of i t h unit at time t (either it is switched on or off) P G i ( t ) denote the output power amount of i t h storage unit whereas in case of   j t h   at time, t. The price of the energy rendered for i t h storage unit is denoted by B G i ( t )   while it was B S j ( t )   for j t h storage unit at time, t. S G i ( t )   and S s j ( t )   denotes the costs involved during start-up or shut-down functions spent upon i t h storage unit and j t h storage unit respectively. Further, P G r i d ( t ) and B G r i d ( t )   denote the quantities of power that was exchanged with the offered market, at time, t [13].
Min   f 1 ( X ) = i = 1 T { i = 1 N g [ U i ( t ) P G i ( t ) B G i ( t ) + S G i | U i ( t ) U i ( t + 1 ) | ] j = 1 N s [ U j ( t ) P s j ( t ) B s j ( t ) + S s j | U j ( t ) U j ( t + 1 ) | ] ( P G r i d ( t ) B G r i d ( t ) ) }

2.1.2. Mitigation of Emission

For the next objective function, i.e., emission mitigation, a few components were considered such as carbon dioxide (CO2), nitrogen oxides (NOx), and sulphur dioxide (SO2), i.e., atmospheric pollutants. Equation (2) shows the mathematical formulation for overall pollutant emission in kg [13]:
Min   f 2 ( X ) = t = 1 T { i = 1 N g [ U i ( t ) P G i ( t ) E G i ( t ) ] + j = 1 N s [ U j ( t ) P s j ( t ) E s j ( t ) ] + ( P G r i d ( t ) E G r i d ( t ) ) }
E G i ( t ) denotes the volume of pollution emitted by ith generation unit, while for the jth storage unit, it was denoted by E s j ( t ) . In the case of volume of pollution at market, it was denoted by E G r i d ( t ) at time, t. All the values are denoted in kg/MWh, respectively.

2.1.3. Constraints and Limitations

(a)
Power Balance constraints
k = 1 N k P L K ( t ) = i = 1 N g [ P G i ( t ) ] + j = 1 N s [ P s j ( t ) ] + ( P G r i d ( t ) )  
Here, P L K corresponds to K amount at load level, whereas N k denotes the whole count of load levels available in the grid.
(b)
Ramp Rate Constraints
This constraint occurs as a result of fluctuations in output power. It can be either an increasing fluctuation or vice versa in DGs. The following is the equation used to describe the constraint [13]:
R d o w n i . Δ t P ( h ) i P ( h 1 ) i R u p i . Δ t
Here, R d o w n i denotes the ramp-down of i t h DG output power whereas R u p i correspond to ramp-up of i t h DG output power. Further, Δ t denotes the time step taken, in terms of hours.
(c)
Inequality constraints
All the units considered in this study have both upper, as well as lower thresholds for their power generation capacity. The units include DGs, storage units, and the market [13].
P G i , m i n ( t ) P G i ( t ) P G i , m a x ( t ) P s j , m i n ( t ) P s j ( t ) P s j , m a x ( t ) P G r i d , m i n ( t ) P G r i d ( t ) P G r i d , m a x ( t )
Equation (6) covers the limitations of charging and discharging rates of the storage unit.
  S O C s j ( t ) = S O C s j ( t 1 ) + P c h g / D c h g ( t ) 0 | P c h g / D c h g ( t ) | P C D S j , m a x
Here, S O C s j ( t ) denotes the charging amount of the storage unit at current time and S O C s j ( t 1 ) denotes the charging amount of the storage unit at previous time. P c h g / D c h g ( t ) corresponds to the charging (discharging) amount at t t h hour while the maximum charging (discharging) rate is denoted by P C D S j , m a x .

3. Oppositional Gradient-Based Grey Wolf Optimizer

Mirjalili et al. (2014) proposed GWO inspired from the hunting behavior of grey wolves [41]. In general, there are four types of wolves, such as the Alpha, followed by Beta, Delta, and finally Omega [42]. In a pack of wolves, Alpha is the dominant wolf and is the decision-maker for the whole pack, while the rest of the members in the pack must obey the alpha’s commands and decisions. Further, only the alphas breed in the pack. The second position is occupied by Beta wolves, which act as a mediator between the alpha wolf and the rest of the wolves. It helps the alpha wolf and becomes the first choice of nominee for alpha, in case an Alpha wolf dies or is too old for swarm management. Though beta wolves control other wolves of the swarm, it fulfils the order of Alpha wolf. Hierarchically, Omega is the least ranked grey wolf [43] and follows the rest of the high-ranking wolves. Those wolves that are not categorized as the other three, i.e., Alpha, Beta, or Omega fall under the category Delta Wolf. Deltas tend to manage the Omega wolves, while at the same time it also helps the Alpha wolf and Beta wolves.
It is important for an algorithm to ensure there is equilibrium between exploration and exploitation skills so that the exploration can be achieved. GWO continues to look for solution space through constant upgradation of the position of dominated wolves such as Alpha, Beta, and Delta. Two new features have been added to the traditional GWO algorithm in this study, which resulted in Oppositional Gradient-Based Grey Wolf Optimizer (OGGWO). The first feature is the adoption of a novel procedure that uses gradient data to improve the algorithm’s exploitation as well as exploration skills. In most of the optimization challenges, the gradient tends to render the critical information with regards to the shape of the solution space. This is done by measuring the steepest slope at every point in the solution space. Afterwards, the particles are transferred to a nearby local optima using gradient information, while at the same time, an appropriate exploration capability is also maintained. GWO is heavily empowered by such updating operators. It results in the efficient searching of solutions in solution space and enhances the exploration potentials of the algorithm to side-step local optima. Opposition Based Learning (OBL) is one of the powerful optimization tools to enhance the convergence speed of numerous heuristic optimization techniques [44]. The successful implementation of OBL involves the evaluation of population at both the opposite site and current site in the same generation so that a better candidate solution can be obtained for the given problem. The concept of OBL has successfully been applied in various meta-heuristics to enhance the convergence speed. The concept of opposite number needs to be defined to explain OBL [44]. A novel procedure is incorporated to update the formulations for Omega wolves in OGGWO algorithm as shown in (7) and (8).
Both Gaussian walks as well as Levy Flight (LF) are added as second features to the existing traditional GWO. These random walks enhance the random nature of the presented OGGWO model and empower it for exploration. Though self-similar clusters (or trajectories) are created by Levy Flight and Gaussian walks, significant differences are found between their structures [45]. Gaussian walk produces a minor, but dense cluster that contains numerous small steps and executes the same number of iterations, such as with clusters [46]. Having been randomly chosen, these methods improve the proposed OGGWO model’s exploration capability by helping the algorithm avoid local optima. In this case, OGGWO switches randomly between Gaussian walk, as well as LF and make use of both the methods [46]. Figure 1 shows the flowchart of the proposed Oppositional Gradient-based Grey Wolf Optimization algorithm. The following formulations are used in OGGWO to upgrade the Omega wolves’ position in the solution space at the time of closing every iteration as shown in (12) and (13).
Step 1: The position of grey wolves (search agents) is randomly initialized in the search space. This step further fixes the number of iterations, followed by population size (number of wolves).
Step 2: For each search agent, the fitness value is determined as it denotes the distance between prey and the wolf.
Step 3: The opposite points are initialized and used to generate the opposite population so as to calculate the fitness values of every individual population.
Step 4: Sorting is executed for both current and opposite populations (pop and opop respectively), in line with their fitness values
Step 5: The n P number of the fittest solution is selected out of a combination of current and the corresponding opposite population.
Step 6: Based on fitness values, three solutions are found such as best (a), second-best (b) and finally, the third best (d). These solutions correspond to alpha, beta and delta category wolves, respectively.
Step 7: Based on Equation (7), the position of the grey wolves is modified.
X W i ( t + 1 ) = γ f M i n X i f X w i ( t ) < γ f M a x X i   f o r   i = 1 , , m   a n d   w = 1 , , n ,
X W i ( t ) r a n d ( 0 , 1 ) λ i ( t ) ( f X w i ( t ) ) ,   O t h e r w i s e , f o r   i = 1 , , m   a n d   w = 1 , , n ,  
Here, i corresponds to the index of decision variables in an optimization problem. On the contrary, n denotes the count of grey wolves. f M a x X i and f M i n X i correspond to the highest positive and the least negative slopes, at every dimension, during every iteration of the algorithm. Further, γ corresponds to a continuous parameter that is calculated as (0,1]. As per the formulation given above, λ i   is updated using Equation (9) as given herewith.
λ i ( t ) = 0.1 [ U b i L b i ] max ( | f M i n X i | , | f M a x X i | )
Here, LB corresponds to Lower Bound, whereas UB denotes the Upper Bound of the problem. This results in the following equation straightforwardly.
| λ i ( t ) f X w i ( t ) | 10 ( U b i L b i )  
The gradient of the problem in few optimization problems may remain unknown. This is attributed to the non-differentiability of objective function or discrete features possessed by decision variables. In order to overcome these issues, the following equation is presented.
f x = [ f ( t ) f ( t 1 ) ] [ X ( t ) X ( t 1 ) ]  
Step 8: The fitness value is updated based on the modified position of grey wolves.
Step 9: Alpha, beta and delta values are updated
Step 10: The position of omega wolves is updated using the Equations (12) and (13) that employ Gaussian walk and Levy flight.
X W , n e w i = X W i + K G a u s s i a n ( | ϑ i | , σ ) ( ξ × ϑ i ξ × X W i ) ,   f o r   i = 1 , , m   a n d   w = 1 , , n ,
  X W , n e w i = X W i + X W i L e v y ( η ) , f o r   i = 1 , , m   a n d   w = 1 , , n ,
Here, the best solution is denoted by ϑ i   whereas |σ| correspond to the standard deviation of Gaussian distribution. The Gaussian parameter is changed by OGGWO as σ = | K × ( x i B P ) | , while it also reduces the length of steps at the time of iterations by fixing = log ( l ) l . Here, l corresponds to the number of iteration. Further, X W , n e w i corresponds to the new position of the wolves, whereas X W i   denotes the current position. In addition to these, ξ and ξ denote random numbers in (0,1]. Equation (14) denotes the calculation for Levy flight.
L e v y ( x ) = [ 0.01 × σ × r 1 ] | r 2 | 1 β  
Here, the random numbers are denoted by r 1 and r 2 between (0, 1]. β corresponds to constant, which is equal to 1.5. In Equation (13), σ is computed with the help of following equation.
σ = [ Γ ( 1 + β ) sin ( π β 2 ) / ( Γ ( 1 + β ) 2 β [ 2 ( β 1 ) 2 ] ) ] 1 / β  
Step 11: The opposite population is created out of the current population using the jumping rate.
Step 12: The n P number of the fittest solution is selected after integrating the current and the opposite populations.
Step 13: The steps 3–7 are repeated until the maximum number of iterations is achieved.
Step 14: The best solution is achieved as the output.
As per the steps discussed above, the leading framework remains the same in both GGWO and OGGWO. However, significant changes are made in the latter. For example, OGGWO leverages a mix of original and gradient-based operators to update the position of wolves more than the conventional GWO operators. Further, a new function is also added in OGGWO to determine the gradient’s objective function at every point in the solution space [47]. In addition to the above, Gaussian walk and Levy Flight have been incorporated in OGGWO to increase randomness during the closure of every iteration. This phenomenon significantly increases both exploration and exploitation capability, as long steps and short steps shift the particles in solution space.

4. Fuzzy Logic-Based Collection of the Finest Compromise Solution

In general, it is important to select the optimal compromise solution out of the available options prior to making any decisions. In this study, a fuzzy membership approach was used by the author to identify the optimal compromised solution [48]. In j t h objective function,   f j   of individual k is characterized by a membership function µ j   k , owing to the indefinite characteristic of the conclusion of the decision maker as denoted herewith.
  µ j k = { 1     f j   f j m i n   f j m a x f j   f j m a x f j m i n     f j m i n < f j < f j m a x   0   f j     f j m a x    
Here, f j m a x corresponds to the highest value of the j t h fitness function. In addition to this, its least value is denoted by f j m i n among the rest of the non-dominated solutions. The calculation is done for normalized membership function µ k that is applicable for each non-dominated solution k.
  µ k = j = 1 N µ j k k = 1 r j = 1 N µ j k
Here, r corresponds to the whole number of non-dominated solutions. The highest value, i.e., µ k is achieved by the best compromise solution.

5. Modeling of the Microgrid

In the case of applying DGs individually, it creates multiple problems. So, it is suggested to follow a systematic approach so that generation, as well as associated loads, is considered as either the microgrid or sub-system to avoid the challenges involved in it and leverage the true potentials of DGs. In parallel, the DGs in the microgrid can be combined together and analyzed about the exploitation of renewable energy sources in huge quantities regarding economy, technology, and environment within the target system. These insights help in making knowledge-based decisions and manage the operations in a better manner. In addition to this, Distribution Generation is a compendium of a vast range of prime mover technologies that include fuel cells, gas turbines, micro turbines, solar PV, IC engines, and wind power. Such RES-based technologies are cost-efficient and emit fewer pollutants, contradicting a huge economy [17]. For instance, Fuel Cells generate electricity out of oxygen and hydrogen, whereas water vapor is the only emission from this source. Having said so, natural gas and other fuel reformation emit NOx and CO2. Fuel cells are comparatively more efficient than microturbines, in terms of emission, which incurs heavy costs. In the current research work, the researchers considered a commonly-known LV microgrid based on the study conducted earlier that included various types of DGs such as PV, WT, MT, low-temperature Fuel Cell (PAFC), and storage devices such as lead-acid batteries [38]. The authors assumed that the active power, produced by whole DG sources, stands as unity power, while at the same time, it does not request or produce reactive power. Furthermore, a power exchange link is also used between utility (LV network) and the mentioned microgrid for the purpose of energy trading at different hours in a day, in line with the decisions taken by the MGCC (Microgrid Central Controller).

6. Results and Discussion

As per the literature [41], the current study considered a microgrid test system with a distributor and a few DGs such as Wind Turbine (WT), Fuel Cell (FC), Photovoltaic panel (PV), Microturbine (MT), and battery. The model, suggested in this study, has two objective functions, such as whole emission of the pollutants and microgrid costs, which include the costs incurred upon power generation and start-up/shut-down costs of units. The chosen problem was rectified through three different scenarios such as main case, secondary case, and third case. In the main issue, the whole set of units were dispatched based on their original drawbacks. In the second case, photovoltaic (PV), as well as Wind Turbine (WT) functions were assumed to be at its maximum output (Max-Ren). For the final case, the utility was considered as an unconstraint unit which can exchange energy with the microgrid without any limits (WLE). The microgrid was expected to meet the total load demands as predicted in [41]. One prime residential area, one feeder with light usage commercial consumers, along with one industrial feeder that supplies power to a small workshop in a typical day was considered. The total energy demand of the day considered was 1695 kWh.
In addition to these, the energy price difference for a real-time market was sourced from the literature [41] for the day considered. In order to ensure flawless microgrid operations, three DG units with WT, PV, and MT were used. The optimization algorithm has the ability to assign these units with on and off states, when there is a power dispatch problem, considering both objectives. Likewise, the on/off state was regarded as 1 for the utility in all cases, since the microgrid was operated in grid-connected mode. To analyze the overall impact of PAFC and battery upon grid operations and to leverage the optimal benefits out of the resources, the on state was selected with an intention for relevant units. From the study conducted earlier, both minimum as well as maximum generation limits were taken for DG sources [41]. Further, the same source was also used to provide the bid coefficients, in terms of Euro cents per kWh. In order to perform the data analysis in a rapid fashion, all the units in the current research work were considered to function under electricity mode only. The microgrid is preferred for installation since it has a salient feature, i.e., the ability to combine multiple sorts of renewable energy sources together.
Microgrid forecasts are highly important for actual operations since this information can be utilized to prepare the systems to be flexible and appropriate. Since it is impossible to execute the renewable energy operations in a conventional manner, it is easy to predict its behavior. Thus, microgrid system efficiency can be enhanced with the help of forecasted information. An expert prediction model was used in the current research work to evaluate the power outputs from PV and WT. From the literature [41], the daily maximum power output from PV and WT values were extracted. Further, the authors used the same study to source information on both daily load power, as well as energy market price in a typical microgrid.
In order to validate the efficacy of the proposed optimization algorithm, simulation was conducted for 50 trial runs to migrate the cost incurred upon operations. The current study considered population size (50) and 1000 maximum iterations as its controlling parameters. Table 1 and Table 2 show the characteristics of 5 DG sources which generate electricity in the microgrid. Through the common coupling point, the microgrid either shares surplus energy or demands additional energy from the utility. The whole set of units that include the macrogrid (utility), are able to function within their power limits, while it can meet the requirements in parallel.

6.1. Case I: Operation of Distributed Energy Sources with Specified Limits

The first test case, considered in this study, was the operation of all DGs and grids within the specified limits as represented in [41]. Table 2 shows the optimal generation schedule for 24 h to minimize both cost and emission levels. From Table 3, it can be understood that the majority of the load demands, during the first few hours of the day, was met by fuel cells present inside the grid and utility through the common coupling point. This is attributed to fewer bids of relevant units than the rest of the others for a specific period. Since the demand increased for high number of utility bids in the next few hours, the DG units enhanced their output power based on the priorities such as low cost and low emission level.
Further, the microgrid regulatory controller sequentially supplied the units and exported the energy from the microgrid, which in turn increased the revenue and reduced the net emission for a specific period under study. On the other hand, the battery was charged for the first few hours a day, especially when the low cost is charged during this period. When the load curve achieved the highest value in the later hours of the day, the discharge action was executed. When RESs such as solar and wind are deployed, it emits less pollutants, whereas the operating cost increases. In other terms, these energy sources are economically unviable due to which the energy generated from these sources needing to be confined to a limit based on emission/economic limitations.
Table 3 shows the analytical outcomes of optimization algorithms and compares the results achieved for the main case. Based on the results attained, in terms of best and worst solutions, it can be inferred that the proposed optimization algorithm accomplished an excellent outcome for cost and emission objectives. Further, it also attained faster convergence characteristics. In addition to this, the proposed algorithm proved to have another advantage for optimization, based on the statistical indices of average and standard deviation.
Table 3 shows the SD values attained by the proposed algorithm for cost and emission objectives, which were 2.55 and 4.72, respectively, which infers that the proposed algorithm is an excellent choice for optimization. In order to showcase the performance of the proposed OGGWO algorithm, Figure 2 and Figure 3 portray the convergence characteristics held by OGGWO, as well as the rest of the optimization algorithms for the best solution and for every single objective. The illustrations show that the proposed method achieved the minimum value for cost objective function after 705 iterations and remained unchanged thereafter. On the other hand, the GWO algorithm converged at 713 iterations. With regards to emission objective function, the proposed method reached the minimum value at 734 iterations, whereas GWO algorithm required 745 iterations for convergence.
By applying the fuzzy logic based approach as discussed in Section 4, the global best compromise solution was attained by the proposed algorithm to mitigate the operating cost and emission.
Figure 4 shows the pareto-fronts of respective trade-off objectives, obtained by various methods, used for comparison and the best-compromised solution. The figure infers that the proposed algorithm was effective in elucidating the non-linear multi-objective optimization problem since the Pareto optimal fronts of non-dominated solutions were distrusted in a good manner. In addition to this, the computational time consumed by the proposed algorithm to mitigate both emission, as well as operating cost, was lesser. This reveals the quality of the proposed algorithm. The optimization results amply corroborate the ability of the proposed algorithm to tackle the issues associated with equality and inequality that tend to happen in energy management challenges.

6.2. Case II: Operation of Renewable Energy Sources at Maximum Limits

Both PV, as well as WT, was considered to operate at maximum power in this case for every hour, whereas the rest of the sources considered such as grid, battery, MT, and FC were operated at permitted levels, as per the literature [41]. Table 4 shows the comparison of statistical performances of OGGWO and the rest of the optimization algorithms for each objective. When the results achieved by GWO were compared for each objective with the proposed OGGWO algorithm, the outcomes inferred that the GWO algorithm was deficient and its performance was inconsistent. For instance, Table 5 shows the performance of the proposed algorithm with regards to the multi-objective optimization and allocation of optimal power to units. The table implies that heavily-performing units must be leveraged to attain economic and emission objectives simultaneously. As per this inference, both the battery, as well as PAFC, is suggested to be extensively used, whereas the rest of the units are to be used as and when required. In the case of light load requirements, i.e., during the first few hours of the day, surplus energy can be exported to the macrogrid. However, if the load requirement is high, the reverse option is preferred. So, utility-supplied energy plays an important role in sustainable operations.
In the second scenario, an assumption was made that RESs such as WT and PV produce the highest amount of power for every hour in a day. On the other hand, the other generators such as a battery, PAFC, MT, and utility function in normal mode, similar to their behavior found in the primary case. In this scenario, the proposed algorithm was again incorporated and validated for the multi-objective optimization problem. The results were attained and inferred through a comparison with that of the other algorithms. Table 5 compares the statistical performances of the algorithms discussed earlier. In this comparison, every single objective underwent 50 trials.
From the table, it can be inferred that the multi-objective optimization problem was appropriately handled by the proposed algorithm while it also ensured a minimum diversity for every objective and trial, compared with other optimal solutions. The total operation cost of the microgrid was increased in the second scenario when compared with the main case. However, the net emission was reduced compared to the previous case. This is attributed to the fact that when RESs penetrate heavily into the grid network, it reduces the net emission. However, it incurs high costs upon the operations for a specific period under study.
The convergence characteristics for operating cost and emission, attained by OGGWO, GGWO, GWO, CSA, DE, and PSO are shown in Figure 5 and Figure 6, respectively. These figures infer that the proposed OGGWO algorithm possesses stable and quick convergence characteristics to find the optimal solution. Figure 7 showcases the distribution of non-dominant solutions, which infers the Pareto optimal fronts of OGGWO and GGWO algorithms. The figure implies that the proposed OGGWO algorithm is a promising candidate to overcome the multi-objective energy management problem compared to the rest of the optimization techniques. In addition to the above, the proposed OGGWO algorithm ensured the distribution of non-dominated solutions in the best manner. It also assured that the solutions are feasible for large-scale standard test systems.

6.3. Case III: Unlimited Power Exchange between LV and MV

The overall set of DG units function as per their power limits in the third scenario. In this case, the utility behaves as an unconstraint unit whereas energy exchange occurs between utility and microgrid irrespective of the limitations. The data required to overcome the multi-objective optimization issue such as real-time market prices, technical specifications of DGs, and load curve remains unaltered. Like the earlier cases, evolutionary-based optimization algorithms were utilized in this case as well to resolve operation management issues found in the microgrid. The simulation results were collected from the experimentation.
Table 6 shows the simulation results for the third case, which was attained to mitigate the operation cost and emission. As per the results, both operating costs, as well as the emission rate of the microgrid, were reduced more than in than the basic scenario, under a few conditions such as allowing unlimited power exchange and economic and environmental objectives. When the second and third scenarios were compared for optimization, in terms of costs and emission, the results infer that the total cost of operation was reduced in the third scenario, while the net emission rate was reduced in the second scenario. Table 7 compares the statistical results of OGGWO and other optimization algorithms, pertaining to every single objective. The comparative analysis results infer that GGWO was deficient and its results were more inconsistent than for the proposed OGGWO algorithm.
One of the significant observations in the third scenario is that during the first few hours of the day, the utility occupied the top position in meeting load demands inside the grid. On the contrary, the energy was purchased from the microgrid in bulk quantities at peak times. Regarding both objectives, when the grid experienced a shortage in power generation, both WT and PV were started; otherwise, both were triggered in the case of requirement for more energy to be exported to the macrogrid. Other DGs such as battery and FC produced the maximum electricity for the majority of the hours in a day under study, whereas MT generated the highest power between 9.00 to 17.00 h.
Figure 8 and Figure 9 show the convergence characteristics of OGGWO, GGWO, GWO, CSA, DE, and PSO for operating cost and emission. The proposed OGGWO achieved a fast and smooth convergence characteristic curve towards reaching the optimal solution. The proposed algorithm successfully overcame the optimization problem for the third time, whereas minor differences can be found in optimal solutions for both the objective functions. Figure 10 represents the trade-off characteristics between operating cost and emission for OGGWO and GGWO which revealed the efficiency of OGGWO in elucidating multi-objective energy management problems. Further, the comparative analysis conducted between the Pareto fronts of the proposed OGGWO algorithm and GGWO algorithm showcases that the former afforded the least and was well-distributed than the latter one. The simulation outcomes established the superiority of the proposed OGGWO in managing equality and inequality issues that occur in energy management problems.

7. Conclusions

Environmental pollution is a significant threat to humanity and other living beings on earth. Across the globe, there is a tremendous increase observed in environmental pollution due to unsustainable and environmentally-harmful human activities. The very existence of human beings on earth is being questioned these days, whereas the measures taken to curb the environmental pollution will only decide the fate of human beings. Traditional power plants, fueled by fossil fuels to produce energy remain the primary source of environmental pollution. Different corrective measures are being taken these days, such as the utilization of advanced technologies, the prevalent implementation of renewable energy sources into the power sector, and the appropriate utilization of energy. These initiatives may curb environmental pollution in the coming years. Microgrids with different RESs such as PV, WT, battery, and so on enact an important role in curbing the environmental pollution caused during conventional energy generation methods. In this background, a novel algorithm has been proposed in this study, i.e., the Oppositional Gradient-based Grey Wolf Optimizer (OGGWO). The proposed algorithm was validated under three different scenarios to establish its effectiveness. A typical microgrid was used as a test system loaded with a battery, FC, WT, MT, and PV. The microgrid was connected with the utility grid, as well as for power exchange as decided by MGCC. Further, the researchers proposed a multi-objective formulation encompassing both emission and operational costs under equality and inequality constraints. The application results were retrieved on a 24-h horizon, whereas the results were compared with other existing techniques. The results establish that the proposed technique outperformed other algorithms and produced excellent results. Further, the results of the proposed algorithm were also compared against GGWO and different optimization algorithms. Based on the comparison results, it can be inferred that the OGGWO algorithm is capable of avoiding premature convergence. The proposed optimization algorithm can also be applied to solve the optimal sizing and placement of FACTS devices in distribution systems.

Author Contributions

Conceptualization, A.R. and K.N.; Methodology, K.N. and S.D.; Software, O.D.M., S.D. and P.P.; Validation, S.D. and I.A.K.; Formal analysis, A.R.; Investigation, S.D.; Resources, A.S.P.; Data curation, A.R., K.N. and S.D.; Writing—original draft, A.R., S.D., N.L. and P.P.; Writing—review & editing, A.R., S.D., N.L. and P.P.; Visualization, O.D.M., N.L. and P.P.; Supervision, O.D.M., S.D. and I.A.K.; Project administration, K.N., O.D.M. and A.S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANNArtificial Neural Network
BARONBranch and Reduce Optimization Navigator
BBSABacktracking search algorithm
BESSBattery Energy Storage System
BPSOBinary Particle Search Algorithm
CSACuckoo Search Algorithm
DERDistributed Energy Resources
DGDistributed Generator
DRDemand Response
DSMDemand Side Management
EDNSGA-IIEconomic Dispatch based Non-Dominated Sorting Genetic Algorithm
EOEquilibrium Optimizer
EVElectric Vehicle
FCFuel Cell
GAMSGeneral Algebraic Modeling System
GBDTGradient Boosting Decision Tree
GGWOGradient based Grey Wolf Optimizer
GOAGrasshopper Optimization Algorithm
GPPMGeneralized Power Prediction Model
GRGGeneralized reduced-gradient
GWOGrey Wolf Optimizer
LEMSLocal Energy Management System
LFLevy Flight
LVLow Voltage
MBAModified Bat Algorithm
MGMicrogrid
MGCCMicrogrid Central Controller
MGEMMicrogrid Energy Management
MOPSOMulti-Objective Particle Search Algorithm
MSFOAMuddy Soil Fish Optimization Algorithm
MTMicroturbine
MVMedium Voltage
NANot Available
NSGA-IINon-Dominated Sorting Genetic Algorithm
OBLOpposition Based Learning
OGGWOOppositional Gradient based Grey Wolf Optimizer
PAFCPhosphoric Acid Fuel Cell
PARPeak-to-Average Ratio
PMORLPreference-based multi-objective reinforcement learning
PSOParticle Search Algorithm
PVPhotovoltaic
QPSOQuantum Particle Search Algorithm
RESRenewable Energy Sources
SOASandpiper Optimization Algorithm
WTWind Turbine

References

  1. Nagarajan, K.; Rajagopalan, A.; Angalaeswari, S.; Natrayan, L.; Mammo, W.D. Combined Economic Emission Dispatch of Microgrid with the Incorporation of Renewable Energy Sources Using Improved Mayfly Optimization Algorithm. Comput. Intell. Neurosci. 2022, 15, 1–22. [Google Scholar] [CrossRef] [PubMed]
  2. Karthik, N.; Parvathy, A.K.; Arul, R. A review of optimal operation of microgrids. Int. J. Electr. Comput. Eng. 2020, 10, 3. [Google Scholar] [CrossRef]
  3. Konstantinopoulos, S.A.; Anastasiadis, A.G.; Vokas, G.A.; Kondylis, G.P.; Polyzakis, A. Optimal management of hydrogen storage in stochastic smart microgrid operation. Int. J. Hydrogen Energy 2018, 43, 1. [Google Scholar] [CrossRef]
  4. Aghajani, G.; Yousefi, N. Multi-objective optimal operation in a micro-grid considering economic and environmental goals. Evol. Syst. 2019, 10, 239–248. [Google Scholar] [CrossRef]
  5. Kamarposhti, M.A.; Colak, I.; Shokouhandeh, H.; Iwendi, C.; Padmanaban, S.; Band, S.S. Optimum operation management of microgrids with cost and environment pollution reduction approach considering uncertainty using Multi-objective NSGAII algorithm. IET Renew. Power Gener. 2022, 16, 1–13. [Google Scholar] [CrossRef]
  6. Lv, T.; Ai, Q.; Zhao, Y. A bi-level multi-objective optimal operation of grid-connected microgrids. Electr. Power Syst. Res. 2016, 131, 60–70. [Google Scholar] [CrossRef]
  7. Lu, X.; Zhou, K.; Yang, S. Multi-objective optimal dispatch of microgrid containing electric vehicles. J. Clean. Prod. 2017, 165, 1572–1581. [Google Scholar] [CrossRef]
  8. Kim, H.J.; Kim, M.K.; Lee, J.W. A two-stage stochastic p-robust optimal energy trading management in microgrid operation considering uncertainty with hybrid demand response. Int. J. Electr. Power Energy Syst. 2021, 124, 106422. [Google Scholar] [CrossRef]
  9. Arumugam, P.; Kuppan, V. A GBDT-SOA approach for the system modelling of optimal energy management in grid connected micro-grid system. Int. J. Energy Res. 2020, 45, 6765–6783. [Google Scholar] [CrossRef]
  10. Gad, Y.; Diab, H.; Abdelsalam, M.; Galal, Y. Smart Energy Management System of Environmentally Friendly Microgrid Based on Grasshopper Optimization Technique. Energies 2020, 13, 5000. [Google Scholar] [CrossRef]
  11. Veluchamy, K.; Veluchamy, M. A new energy management technique for microgrid system using muddy soil fish optimization algorithm. Int. J. Energy Res. 2021, 45, 14824–14844. [Google Scholar] [CrossRef]
  12. Jain, D.K.; Tyagi, S.K.S.; Neelakandan, S.; Prakash, M.; Natrayan, L. Metaheuristic optimization-based resource allocation technique for cybertwin-driven 6G on IoE environment. IEEE Trans. Ind. Inform. 2021, 18, 4884–4892. [Google Scholar] [CrossRef]
  13. Aghajani, G.; Ghadimi, N. Multi-objective energy management in a micro-grid. Energy Rep. 2018, 4, 218–225. [Google Scholar] [CrossRef]
  14. Ahmed, D.; Ebeed, M.; Ali, A.; Alghamdi, A.S.; Kamel, S. Multi-Objective Energy Management of a Micro-Grid Considering Stochastic Nature of Load and Renewable Energy Resources. Electronics 2021, 10, 403. [Google Scholar] [CrossRef]
  15. Mansouri, S.A.; Ahmarinejad, A.; Nematbakhsh, E.; Javadi, M.S.; Jordehi, A.R.; Catalao, J.P. Energy management in microgrids including smart homes: A multi-objective approach. Sustain. Cities Soc. 2021, 69, 102852. [Google Scholar] [CrossRef]
  16. Tabar, V.S.; Jirdehi, M.A.; Hemmati, R. Energy management in microgrid based on the multi objective stochastic programming incorporating portable renewable energy resource as demand response option. Energy 2017, 118, 827–839. [Google Scholar] [CrossRef]
  17. Mandal, S.; Mandal, K.K. Optimal energy management of microgrids under environmental constraints using chaos enhanced differential evolution. Renew. Energy Focus 2020, 34, 129–141. [Google Scholar] [CrossRef]
  18. GM Abdolrasol, M.; Hannan, M.A.; Hussain, S.M.S.; Ustun, T.S.; Sarker, M.R.; Ker, P.J. Energy Management Scheduling for Microgrids in the Virtual Power Plant System Using Artificial Neural Networks. Energies 2021, 14, 6507. [Google Scholar] [CrossRef]
  19. Majumder, I.; Dash, P.K.; Dhar, S. Real-time Energy Management for PV–battery–wind based microgrid using on-line sequential Kernel Based Robust Random Vector Functional Link Network. Appl. Soft Comput. 2021, 101, 107059. [Google Scholar] [CrossRef]
  20. Karthik, N.; Parvathy, A.K.; Arul, R.; Jayapragash, R.; Narayanan, S. Economic load dispatch in a microgrid using Interior Search Algorithm. In 2019 Innovations in Power and Advanced Computing Technologies (i-PACT); IEEE: Piscataway, NJ, USA, 2019; pp. 1–6. [Google Scholar]
  21. Esapour, K.; Abbasian, M.; Saghafi, H. Intelligent energy management in hybrid microgrids considering tidal, wind, solar and battery. Int. J. Electr. Power Energy Syst. 2021, 127, 106615. [Google Scholar] [CrossRef]
  22. Aghajani, G.R.; Shayanfar, H.A.; Shayeghi, H. Demand side management in a smart micro-grid in the presence of renewable generation and demand response. Energy 2017, 126, 622–637. [Google Scholar] [CrossRef]
  23. Jasim, A.M.; Jasim, B.H.; Kraiem, H.; Flah, A. A Multi-Objective Demand/Generation Scheduling Model-Based Microgrid Energy Management System. Sustainability 2022, 14, 10158. [Google Scholar] [CrossRef]
  24. Yin, N.; Abbassi, R.; Jerbi, H.; Rezvani, A.; Müller, M. A day-ahead joint energy management and battery sizing framework based on θ-modified krill herd algorithm for a renewable energy-integrated microgrid. J. Clean. Prod. 2021, 282, 124435. [Google Scholar] [CrossRef]
  25. Han, Y.; Yang, H.; Li, Q.; Chen, W.; Zare, F.; Guerrero, J.M. Mode-triggered droop method for the decentralized energy management of an islanded hybrid PV/hydrogen/battery DC microgrid. Energy 2020, 199, 117441. [Google Scholar] [CrossRef]
  26. Han, Y.; Zhang, G.; Li, Q.; You, Z.; Chen, W.; Liu, H. Hierarchical energy management for PV/hydrogen/battery island DC microgrid. Int. J. Hydrogen Energy 2019, 44, 11. [Google Scholar] [CrossRef]
  27. Kafetzis, A.; Ziogou, C.; Panopoulos, K.D.; Papadopoulou, S.; Seferlis, P.; Voutetakis, S. Energy management strategies based on hybrid automata for islanded microgrids with renewable sources, batteries and hydrogen. Renew. Sustain. Energy Rev. 2020, 134, 110118. [Google Scholar] [CrossRef]
  28. Xu, J.; Li, K.; Abusara, M. Preference based multi-objective reinforcement learning for multi-microgrid system optimization problem in smart grid. Memetic Comp. 2022, 14, 225–235. [Google Scholar] [CrossRef]
  29. Jordehi, A.R.; Javadi, M.S.; Catalão, J.P. Energy management in microgrids with battery swap stations and var compensators. J. Clean. Prod. 2020, 272, 122943. [Google Scholar] [CrossRef]
  30. Sedighizadeh, M.; Esmaili, M.; Mohammadkhani, N. Stochastic multi-objective energy management in residential microgrids with combined cooling, heating, and power units considering battery energy storage systems and plug-in hybrid electric vehicles. J. Clean. Prod. 2018, 195, 301–317. [Google Scholar] [CrossRef]
  31. Yaghi, M.; Luo, F.; El Fouany, H.; Junfeng, L.; Jiajian, H.; Jun, Z. Multi-Objective optimization for Microgrid Considering Demand Side Management. In Proceedings of the 2019 Chinese Control Conference (CCC), Guangzhou, China, 27–30 July 2019; pp. 7398–7403. [Google Scholar] [CrossRef]
  32. Luo, L.; Abdulkareem, S.S.; Rezvani, A.; Miveh, M.R.; Samad, S.; Aljojo, N.; Pazhoohesh, M. Optimal scheduling of a renewable based microgrid considering phototaic system and battery energy storage under uncertainty. J. Energy Storage 2020, 28, 101306. [Google Scholar] [CrossRef]
  33. Farinis, G.K.; Kanellos, F.D. Integrated energy management system for Microgrids of building prosumers. Electr. Power Syst. Res. 2021, 198, 107357. [Google Scholar] [CrossRef]
  34. Kakran, S.; Chanana, S. Operation management of a renewable microgrid supplying to a residential community under the effect of incentive-based demand response program. Int. J. Energy Environ. Eng. 2019, 10, 121–135. [Google Scholar] [CrossRef] [Green Version]
  35. Mosa, M.A.; Ali, A.A. Energy management system of low tage dc microgrid using mixed-integer nonlinear programing and a global optimization technique. Electr. Power Syst. Res. 2021, 192, 106971. [Google Scholar] [CrossRef]
  36. Hasankhani, A.; Hakimi, S.M. Stochastic energy management of smart microgrid with intermittent renewable energy resources in electricity market. Energy 2021, 219, 119668. [Google Scholar] [CrossRef]
  37. De, M.; Das, G.; Mandal, K.K. An effective energy flow management in grid-connected solar–wind-microgrid system incorporating economic and environmental generation scheduling using a meta-dynamic approach-based multiobjective flower pollination algorithm. Energy Rep. 2021, 7, 2711–2726. [Google Scholar] [CrossRef]
  38. Dey, B.; Bhattacharyya, B.; Márquez, F.P.G. A hybrid optimization-based approach to solve environment constrained economic dispatch problem on microgrid system. J. Clean. Prod. 2021, 307, 127196. [Google Scholar] [CrossRef]
  39. Dey, B.; Bhattacharyya, B.; Srivastava, A.; Shivam, K. Solving energy management of renewable integrated microgrid systems using crow search algorithm. Soft Comput. 2020, 24, 10433–10454. [Google Scholar] [CrossRef]
  40. Paliwal, N.K.; Singh, A.K.; Singh, N.K. Energy scheduling optimisation of an islanded microgrid via artificial bee colony guided by global best, personal best and asynchronous scaling factors. Int. J. Sustain. Energy 2020, 39, 6. [Google Scholar] [CrossRef]
  41. Sarshar, J.; Moosapour, S.S.; Joorabian, M. Multi-objective energy management of a micro-grid considering uncertainty in wind power forecasting. Energy 2017, 139, 680–693. [Google Scholar] [CrossRef]
  42. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  43. Abdel-Basset, M.; El-Shahat, D.; El-henawy, I.; de Albuquerque, V.H.C.; Mirjalili, S. A new fusion of grey wolf optimizer algorithm with a two-phase mutation for feature selection. Expert Syst. Appl. 2020, 139, 112824. [Google Scholar] [CrossRef]
  44. Dhargupta, S.; Ghosh, M.; Mirjalili, S.; Sarkar, R. Selective opposition based grey wolf optimization. Expert Syst. Appl. 2020, 151, 113389. [Google Scholar] [CrossRef]
  45. Pradhan, M.; Roy, P.K.; Pal, T. Grey wolf optimization applied to economic load dispatch problems. Int. J. Electr. Power Energy Syst. 2016, 83, 325–334. [Google Scholar] [CrossRef]
  46. Khalilpourazari, S.; Doulabi, H.H.; Çiftçioğlu, A.Ö.; Weber, G.W. Gradient-based grey wolf optimizer with Gaussian walk: Application in modelling and prediction of the COVID-19 pandemic. Expert Syst. Appl. 2021, 177, 114920. [Google Scholar] [CrossRef] [PubMed]
  47. Pahnehkolaei, S.M.A.; Alfi, A.; Sadollah, A.; Kim, J.H. Gradient-based water cycle algorithm with evaporation rate applied to chaos suppression. Appl. Soft Comput. 2017, 53, 420–440. [Google Scholar] [CrossRef]
  48. Karthik, N.; Parvathy, A.K.; Arul, R.; Padmanathan, K. Multi-objective optimal power flow using a new heuristic optimization algorithm with the incorporation of renewable energy sources. Int. J. Energy Environ. Eng. 2021, 12, 641–678. [Google Scholar] [CrossRef]
Figure 1. A flowchart of the Oppositional Gradient-based Grey Wolf Optimizer.
Figure 1. A flowchart of the Oppositional Gradient-based Grey Wolf Optimizer.
Energies 15 09024 g001
Figure 2. The cost convergence characteristic for Case-I.
Figure 2. The cost convergence characteristic for Case-I.
Energies 15 09024 g002
Figure 3. The emission convergence characteristic for Case-I.
Figure 3. The emission convergence characteristic for Case-I.
Energies 15 09024 g003
Figure 4. The emission cost trade-off characteristic for Case-I.
Figure 4. The emission cost trade-off characteristic for Case-I.
Energies 15 09024 g004
Figure 5. The cost convergence characteristic for Case-II.
Figure 5. The cost convergence characteristic for Case-II.
Energies 15 09024 g005
Figure 6. The emission convergence characteristic for Case-II.
Figure 6. The emission convergence characteristic for Case-II.
Energies 15 09024 g006
Figure 7. The emission cost trade-off characteristic for Case-II.
Figure 7. The emission cost trade-off characteristic for Case-II.
Energies 15 09024 g007
Figure 8. The cost convergence characteristic for Case-III.
Figure 8. The cost convergence characteristic for Case-III.
Energies 15 09024 g008
Figure 9. The emission convergence characteristic for Case-III.
Figure 9. The emission convergence characteristic for Case-III.
Energies 15 09024 g009
Figure 10. The emission cost trade-off characteristic for Case-III.
Figure 10. The emission cost trade-off characteristic for Case-III.
Energies 15 09024 g010
Table 1. A review of optimization techniques in the energy management of a microgrid.
Table 1. A review of optimization techniques in the energy management of a microgrid.
Reference No.ObjectivesControl AlgorithmSystem DescriptionStorage TechnologyRESYear
[21]Operation costModified bat algorithmHybrid AC and DC microgridBatteryPV, WT, tidal2021
[22]Operation cost and EmissionMulti-Objective Particle Swarm Optimization (MOPSO)Grid connected microgridBatteryPV, WT2017
[23]Operation cost, Emission Binary Orientation Search AlgorithmGrid connected microgridBatteryPV, WT2022
[24]Operation costΘ-modified krill algorithmGrid connected microgridBatteryWT, PV2021
[25]Decentralized energy
management
HOMER SoftwareIslanded DC microgridHydrogen, BatteryPV2020
[26]Hierarchical energy management
Strategy
HOMER pro
Software
Islanded DC microgrid Hydrogen, BESSPV2019
[27]Energy Management SystemHybrid automata algorithmIslanded Microgrid Hydrogen, BESSPV, WT, biomass2020
[28]Multi-microgrids energy managementPreference-based multi-objective reinforcement
learning (PMORL) technique
Grid connected MicrogridBatteryPV, WT, tidal2022
[29]Optimal power flow with reactive power cost of MG as objective functionGeneralized reduced-gradient (GRG) algorithmGrid connected MicrogridBESSPV, WT2020
[30]Stochastic optimal energy management of MG with operation cost and emission as objectivesGAMS using CPLEX solverGrid connected Microgrid BESS, PHEV, TESPV2018
[31]Energy Management SystemNSGA-IIGrid connected Microgrid BatteryPV, WT2019
[32]Energy Management SystemModified bat algorithm (MBA)Grid connected MicrogridBESSPV, WT2020
[33]Energy Management SystemParticle swarm optimization (PSO)Grid conncted and autonomous Microgrid BESS, PEVPV, WT2021
[34]Energy Management in MicrogridMixed integer linear programing technique using CPLEX solverGrid connected Microgrid with PV, FC, MT, batteryBatteryPV, WT2019
[35]Energy Management System of DC microgridBranch and reduce optimization navigator (BARON) algorithmMicrogrid with PV, FC, MT, DE, batteryBatteryPV, WT2021
[36]Stochastic energy management of smart microgridQuantum Particle Swarm Optimization (QPSO)Grid connected Microgrid BatteryPV, WT2021
[37]Energy management in microgridmeta-dynamic-approach-based multiobjective flower pollination algorithmGrid connected MicrogridBattery PV, WT2021
[38]Economic Emission Dispatch in MicrogridHybrid Modified version of GWO3-unit
RES integrated low voltage microgrid system
-PV, WT2021
[39]Energy management System Crow Search Algorithm 3,5,6,7,11, 38-units microgrid test systemBESSPV, WT2019
[40]Energy management SystemModified version of artificial bee colony algorithmIslanded Microgrid with stationary and the dynamic energy BidBESSPV, WT2020
Table 2. An optimal generation schedule for Case-I.
Table 2. An optimal generation schedule for Case-I.
HourMTFCPVWTESSGrid
116.173008.84−3036.99
219.7426.98010.62−4.797.45
3622.0108.01−20.0544.03
46.7530012.44−12.3824.19
58.6823.02012.01−19.9843.27
625.280013.56−18.1354.29
725.0223.9012.03−25.5748.62
8030013.372.3444.29
927.3927.030.1715.19−27.2348.45
1029.99302.0419.3429.05−14.42
113029.518.0324.2626.18−24.98
1229.1828.039.7222.0311.03−11.99
133029.8611.0319.2−29.9825.89
140308.9124.0929.13−6.13
1530308.4424.9927.24−29.67
1622.4824.663.9919.98−20.2345.12
1730301.9823.47−29.9846.53
18300018.989.8146.21
193022.19019.011.0235.78
203030022.31−0.5922.28
2117.9919.79012.996.9435.29
221.0228.93021.0827.985.99
2322.050913.09−8.1642.02
2400020.1129.4117.48
Table 3. A statistical comparison of simulation results for case-I.
Table 3. A statistical comparison of simulation results for case-I.
Optimization AlgorithmParameterCost Euro Emission Kg
EDNSGA-II [41]Min Cost148.671230.66
Min Emission163.971175.99
BCS158.811212.7
Std Dev2.95.08
NSGA-II [41]Min Cost165.491236.5
Min Emission173.771199.08
BCS168.271216.54
Std Dev3.157.08
PSOMin Cost168.151241.29
Min Emission176.041203.12
BCS171.381220.07
Std Dev4.288.14
CSAMin Cost167.381239.86
Min Emission175.421202.18
BCS170.461219.14
Std Dev4.067.99
GWOMin Cost161.141234.08
Min Emission170.781195.64
BCS164.261214.31
Std Dev3.096.87
GGWOMin Cost147.391228.14
Min Emission162.041173.67
BCS157.131212.03
Std Dev2.714.88
OGGWOMin Cost146.441226.71
Min Emission161.281172.48
BCS156.311211.28
Std Dev2.554.72
Table 4. A statistical comparison of simulation results for case-II.
Table 4. A statistical comparison of simulation results for case-II.
Optimization AlgorithmParameterCost EuroEmission Kg
EDNSGA-II [41]Min Cost190.6NA
Min EmissionNA1162.4
BCSNANA
Std DevNANA
PSOMin Cost194.241205.29
Min Emission203.151168.12
BCS198.571184.89
Std Dev5.987.63
CSAMin Cost193.381203.14
Min Emission201.081166.58
BCS196.781183.28
Std Dev5.797.58
GWOMin Cost189.731200.23
Min Emission198.291160.25
BCS192.041179.27
Std Dev4.766.64
GGWOMin Cost187.131197.38
Min Emission202.871158.03
BCS190.481177.49
Std Dev4.434.63
OGGWOMin Cost186.271210.24
Min Emission201.891156.79
BCS188.231175.82
Std Dev4.264.46
Table 5. An optimal generation schedule for Case-II.
Table 5. An optimal generation schedule for Case-II.
HourMTFCPVWTBatteryGrid
129.98009.15−28.1351
217.2216.89010.84−16.9331.98
36.7521.56011.56−18.3138.44
423.153012.68−20.8943.06
56.2319.23013.29−29.9858.23
626.993.73013.41−25.1255.99
726.8624.91014.78−24.8942.34
86.9125.99013.95−14.9858.13
90300.4519.1224.3117.12
1030302.2919.2527.45−12.99
1130308.0224.1130−29.13
1229.98010.2322.4329.92−4.56
13021.5610.7518.79−5.9940.89
14303010.1824.1223.11−31.41
15009.7824.8719.2337.12
1630305.0924.92−29.9835.97
1724.3528.150.5623.12−26.9952.81
1811.9820.03020.5418.2234.23
1929.9530020.79−29.9857.24
2029.953.56022.59−8.3456.24
213030022.34−3040.66
2227.980023.8929.733.4
2319.9819.45022.57−19.2335.23
248.7820.12020.61−20.4837.97
Table 6. A statistical comparison of simulation results for case-III.
Table 6. A statistical comparison of simulation results for case-III.
Optimization Algorithm ParameterCost EuroEmission Kg
EDNSGA-II [41]Min Cost127.3NA
Min EmissionNA1506.9
BCSNANA
Std DevNANA
PSOMin Cost133.121552.48
Min Emission142.391512.83
BCS137.031530.25
Std Dev3.789.83
CSAMin Cost132.281549.28
Min Emission140.591511.35
BCS135.251528.94
Std Dev3.599.64
GWOMin Cost126.741545.85
Min Emission135.381505.37
BCS132.061526.12
Std Dev2.618.57
GGWOMin Cost125.131556.29
Min Emission140.561503.32
BCS130.271522.48
Std Dev2.346.69
OGGWOMin Cost124.461556.73
Min Emission139.591502.68
BCS129.061520.83
Std Dev2.186.52
Table 7. An optimal generation schedule for Case-III.
Table 7. An optimal generation schedule for Case-III.
HourMTFCPVWTBatteryGrid
16004.97−13.9965.02
212.233.1800.34−10.9955.24
3000.2710.98−29.2377.98
40005.65−13.9969.34
500011.98−3085.02
603013.21−28.2387.02
703.25012.99−25.1292.78
86008.28−17.2392.95
930300.1318.12−15.2928.04
1030301.6819.2329.76−14.67
113029.997.132330−27.12
1230309.8922.2330−34.12
1330011.0316.0227.031.92
149.569.322.997.949.8946.3
1530308.0724.2918.99−20.35
1611.2311.281.868.758.7653.12
1729.743.931.0322.27−3075.03
1800018.99−30116.01
196.9315.24012.96−20.1292.99
207.123015.38−21.1199.61
213030020.0329.89−16.92
2200022.34−29.9992.65
2303.2308.0312.0254.72
2400020.01−3076.99
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rajagopalan, A.; Nagarajan, K.; Montoya, O.D.; Dhanasekaran, S.; Kareem, I.A.; Perumal, A.S.; Lakshmaiya, N.; Paramasivam, P. Multi-Objective Optimal Scheduling of a Microgrid Using Oppositional Gradient-Based Grey Wolf Optimizer. Energies 2022, 15, 9024. https://doi.org/10.3390/en15239024

AMA Style

Rajagopalan A, Nagarajan K, Montoya OD, Dhanasekaran S, Kareem IA, Perumal AS, Lakshmaiya N, Paramasivam P. Multi-Objective Optimal Scheduling of a Microgrid Using Oppositional Gradient-Based Grey Wolf Optimizer. Energies. 2022; 15(23):9024. https://doi.org/10.3390/en15239024

Chicago/Turabian Style

Rajagopalan, Arul, Karthik Nagarajan, Oscar Danilo Montoya, Seshathiri Dhanasekaran, Inayathullah Abdul Kareem, Angalaeswari Sendraya Perumal, Natrayan Lakshmaiya, and Prabhu Paramasivam. 2022. "Multi-Objective Optimal Scheduling of a Microgrid Using Oppositional Gradient-Based Grey Wolf Optimizer" Energies 15, no. 23: 9024. https://doi.org/10.3390/en15239024

APA Style

Rajagopalan, A., Nagarajan, K., Montoya, O. D., Dhanasekaran, S., Kareem, I. A., Perumal, A. S., Lakshmaiya, N., & Paramasivam, P. (2022). Multi-Objective Optimal Scheduling of a Microgrid Using Oppositional Gradient-Based Grey Wolf Optimizer. Energies, 15(23), 9024. https://doi.org/10.3390/en15239024

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