Effect of Battery Degradation on the Probabilistic Optimal Operation of Renewable-Based Microgrids

: In order to maximize the use of renewable-based distributed generators (DGs), in addition to dealing with the effects of the inherent power management uncertainties of microgrids (MGs), applying storage devices is essential in the electrical system. The main goal of this paper is to minimize the total operation cost as well as the emissions of MG energy resources, alongside the better utilization of renewable energy sources (RES) and energy storage systems. The uncertainties of wind speed, solar irradiation, market price and electrical load demand are modeled using reduced unscented transformation (RUT) method. Simulation results reveal that, as expected, by increasing the battery efﬁciency, the achievable minimum daily operational cost of the system is reduced. For ex-ample, with 93% battery efﬁciency, the operational cost equals EUR 9200, while for an efﬁciency of 97%, the achievable minimum daily operational cost is EUR 8900. Moreover, the proper eco-nomic/environmental performance of the suggested approach, which contributes to the possibility of selecting a compromise solution for the MG operator in accordance with technical and economic constraints, is justiﬁed.


Introduction
In order to fulfill the increasing energy demand and to mitigate the environmental concerns related to energy generation from fossil fuels, the development of RESs has been expedited [1][2][3][4][5]. MGs, as an aggregation of DERs and loads with distribution networks, were presented to provide three main advantageous features, including reliability, sustainability and cost-effectiveness [1]. DERs include both renewable (such as PV systems and WTs) and non-renewable technologies (such as micro-turbines and diesel generators). RESs attracted much attention owing to their availability and sustainability. However, their inherent uncertainties raise new challenges [2]. These challenges, however, can be resolved by combining energy storage devices to RESs in MGs [2]. Energy storage devices, such as batteries, play an important role in the stability, reliability and performance of renewable-based MGs [1][2][3][4][5]. A review of different storage technologies and some major storage devices that are used in MGs were deliberated in [3,4], respectively.
The multi-objective probabilistic economic/emission optimal operation of renewablebased MGs is considered in this paper. A variety of methodologies were suggested in the literature to solve the multi-objective problems related to energy systems, including 1.
Meta-heuristic algorithms deal simultaneously with a set of feasible solutions; this allows different solutions to be found in the Pareto optimal front in just one execution of the algorithm. While in the mathematical programming approaches, a sequence of independent executions should be dealt with. 2.
Meta-heuristic algorithms are not sensitive to the continuity and formation of the Pareto front, which is one of the drawbacks of mathematical programming.
Research was conducted that dealt with the optimal operation management of MGs, considering different constraints and objectives. Some papers focused on conventional economic dispatch [6,7], while some studied the scheduling of energy storage systems [8]. Soares et al. [9] proposed a "signaled particle swarm optimization" for MG scheduling. A day-ahead MG scheduling problem was considered in [10], using a hybrid harmony search algorithm with differential evolution. The authors of [11] proposed a robust optimization-based method for solving optimal MG management problems. Elsied et al. [12] applied a binary PSO to optimize MG performance in real-time operations, where the objectives were to minimize the energy cost, carbon dioxide and pollutant emissions, while maximizing the power of the available RESs. In [13], a "real-time energy management system" was introduced to take the advantage of GA for minimizing energy costs and carbon dioxide emissions, while maximizing the power of available RESs. A multi-objective PSO algorithm was used in [14] for the short-term minimization of operating costs and emissions in a smart MG. In [15], a "multi-period artificial bee colony-based approach", combined with a Markov chain, was proposed for solving the problem of the optimal energy management of an islanded MG. A grey wolf optimization technique was applied in [16] to minimize the operation cost of MG. Marzband et al. scrutinized the optimal operation scheduling of MGs using a gravitational search algorithm, while fulfilling technical constraints [17]. A two-stage optimization method was presented in [18], which applied to PSO along with a deterministic technique based on mixed-integer linear programming for the daily operation of a smart grid. A multi-objective, model-based optimization method was expressed in [19] for the optimal sizing of all components and determination of the power electronic layout. In order to minimize the economic and environmental objectives, the problem was divided into three optimization problems.
Some papers were devoted to ODED problems [20][21][22]. The DED of a hybrid MG was solved in [20] considering cost minimization as the objective function, while a buildingbased virtual energy storage system was taken into account. The authors of [21] applied a chaotic GA to decrease operational costs including the economic and environmental costs of the power plant with WTs. In [22], a group search optimizer was proposed, in which transmission losses and ramp rate limits were the constraints. Azizipanah et al. [23] implemented a "θ-improved cuckoo optimization" algorithm in order to resolve the deficiencies of some previous methods.
The MG energy management problem was investigated in [5,[24][25][26][27]. In order to deal with grid-connected MG optimal energy management, the authors in [24] applied a GA algorithm, while in [6], a modified multi-objective bird-mating optimization algorithm was applied to minimize operational cost and emissions. The modified multi-objective grey wolf optimizer was utilized to study the multi-objective day-ahead scheduling of MGs [25]. In [26], a whale optimization algorithm was used. The authors of [27] employed a modified PSO algorithm for dealing with single-objective MG energy management. However, in [24,26], renewable-based DGs are considered as dispatchable units, which is an incorrect assumption from a technical viewpoint. Moreover, in [5,[24][25][26][27], the uncertainties related to RESs, market price and load demand were not investigated.
For properly dealing with multi-objective optimization problems, a meta-heuristic algorithm should be so robust that it cannot become trapped in local optima, while also converging to a suitable Pareto optimal front in a relatively short period of time [5]. In comparison with other meta-heuristic algorithms, "JAYA", introduced by "R. Venkata Rao" [28] manifested an approximately accurate search speed and convergence in solving optimization problems. The major convenience of JAYA is that there is no need to tune the algorithmic parameters. In each iteration performed by this algorithm, particles move towards the best solution and avoid the worst solution so that the global solution be optimally found [28].
In this paper, in order to investigate the MG's MOOM problem, the EMOJAYA algorithm is suggested. Different modifications are proposed such that the search space can become more immense, which leads to the more precise identification of the optimum global (i.e., best-compromised) solution. Most algorithms need an onerous tuning process for controlling their parameters, while this cumbersome process is not required in EMOJAYA, which makes it predominant when compared to all other algorithms. In short, a fast convergence and low computational time can be highlighted as the outstanding features of the proposed algorithm as well as its ability to independently control the algorithmic parameters.
On the other hand, the application of RESs increased in new MGs in order to achieve clean energy production. However, the operation and control of MGs have become complicated due to load demand and market price fluctuations, as well as the inherent uncertainties of RESs. Consequently, numerous approaches have been developed for the sake of a probabilistic analysis of MGs' performance [29]. A comprehensive comparison of different probabilistic approaches in solving MGs' optimal operation management problem was conducted in our previous work [29]. The multi-objective, probabilistic MG optimal energy management was mentioned as a possible and valuable future work in [29], which is addressed in this article.
In order to deal with the uncertainties related to the forecasted values of load demand, market price and the available outputs of RESs, the reduced unscented transformation (RUT) approach is applied. In the suggested RUT method, the number of samples is considerably reduced compared to other probabilistic methods. Consequently, applying RUT in solving problems with too many random variables leads to a significant reduction in the computational time of the process. In addition to being time-efficient, RUT is very accurate for dealing with problems with high data rates [29].
In this paper, an RUT-EMOJAYA algorithm is proposed to tackle the probabilistic economic/emission optimal operation of renewable-based MGs. The effectiveness of the proposed RUT-EMOJAYA algorithm in solving probabilistic MG's MOOM problem is verified by applying the method on the low-voltage, grid-connected MG of Figure 1, including five DG units, i.e., PV, WT, FC and MT, along with a battery as the storage device of the test system [30]. The MG is connected to the main utility grid through the PCC as it is able to exchange power. To control the battery charging current, avoid battery overcharging and satisfy the SOC limits, a charging controller is needed. Each DG unit has an LC to communicate with the MGCC for applying operational commands.
The main contributions of the paper are summarized as follows: 1.
The power constraints of the storage device (i.e., Li-ion battery in this study) as well as the degradation cost are considered in the MG's MOOM problem. In order to consider different battery characteristics, including the battery efficiency, the battery's initial charge, SOC, and different scenarios of the battery's degradation cost are studied in the probabilistic MG's MOOM problem.

2.
Modifications are added to the JAYA algorithm that make it more efficient in dealing with multi-objective problems. The efficiency of the suggested algorithm is examined by comparing its performance with some other well-known algorithms.

3.
The total cost of day-ahead market transactions and fuel costs, along with the emission of MG, are minimized through the introduced optimal scheduling approach. The suggested RUT-EMOJAYA reduces the MG's dependency on the main grid and the electricity market, while maximizing the utilization of RESs in the studied region.
The remainder of this paper is organized as follows: the economic model of the storage device is described in Section 2. In Section 3, the problem is formulated and the objective functions along with the constraints are presented. The RUT method is discussed in Section 4. Section 5 is devoted to the presentation of the applied algorithm. In Section 6, the applied methodology for solving MOOM in an ordinary MG is explained. Results are investigated in Section 7. Finally, the conclusions are highlighted in Section 8.

Economic Model of Battery Storage Devices
An economic model of storage devices is presented in this section. As dispatchable units, storage devices can be charged in low demand hours, while they discharge in peak load hours when there is a power shortage in the system [31]. In this paper, the degradation cost of the storage device is added to the operational cost of the system in order to improve the cost objective function. During the discharge hours, the degradation cost is a function of a number of factors, including the battery DoD at start and end of the discharging process, the investment cost and the usable energy of the battery, as well as battery- Figure 1. A typical MG test system [30].
The remainder of this paper is organized as follows: the economic model of the storage device is described in Section 2. In Section 3, the problem is formulated and the objective functions along with the constraints are presented. The RUT method is discussed in Section 4. Section 5 is devoted to the presentation of the applied algorithm. In Section 6, the applied methodology for solving MOOM in an ordinary MG is explained. Results are investigated in Section 7. Finally, the conclusions are highlighted in Section 8.

Economic Model of Battery Storage Devices
An economic model of storage devices is presented in this section. As dispatchable units, storage devices can be charged in low demand hours, while they discharge in peak load hours when there is a power shortage in the system [31]. In this paper, the degradation cost of the storage device is added to the operational cost of the system in order to improve the cost objective function. During the discharge hours, the degradation cost is a function of a number of factors, including the battery DoD at start and end of the discharging process, the investment cost and the usable energy of the battery, as well as battery-specific parameters [32]. Different factors can affect the N cycle of a battery, including, DoD, discharge rate, ambient temperature, charging regime, etc. N cycle is defined as the number of complete charge-discharge cycles that the battery can undergo before its nominal capacity falls below 80% of its initial rated capacity [31].
DoD is an absolute discharge relative to the rated battery capacity that is a function of the battery's state of charge (SOC). SOC (SOC(t) = Q(t) Q n ) is defined as the ratio of the battery's current capacity, Q(t) (kWh), to the battery's nominal capacity, Q n (kWh). Q n is representative of the maximum amount of charge that can be stored in the battery and is given by the manufacturer [33]. In this paper, the DoD of the battery is considered as follows [32]: where E Batt,max is the maximum energy stored in the battery and E t is the energy stored in the battery at a specific hour of the day. The cycle life of the battery is as follows: where a and b are battery-specific parameters considered equal to 1331 and −1.825, respectively, for Li-ion batteries [32]. Consequently, the battery degradation cost Cost t Deg is formulated as: where C Batt is the investment cost of the battery, and E Batt is the usable energy of the battery [32].

Objective Functions
The first objective is to minimize the total operation cost as presented in the following [33]: → P Batt,s = P 1 Batt,s , P 2 Batt,s , . . . , P t Batt,s , . . . , P T Batt,s ; s = 1, 2, . . . , N Batt (12) Electricity 2022, 3 58 → P RES r = P 1 RES r , P 2 RES r , . . . , P t RES r , . . . , P T RES r ; r = 1, 2, . . . , N RES (14) → P Grid = P 1 Grid , P 2 Grid , . . . , P t Grid , . . . , P T Grid (15) where DG demonstrates the dispatchable units including fuel cell (FC) and micro-turbine (MT), while Grid and Batt present the utility grid and the battery, respectively. The WT and photovoltaic (PV) are shown with RES as non-dispatchable units. Note that the utility has to buy all of the electrical power produced by RES units; consequently, the RESs' output powers (P t RES,r ) are not included in the design variables' vector; u t i is applied to imply the on/off states of the ith dispatchable DG during each hour of the day.
The minimization of environmental pollutants is considered as the second objective as follows [33]: where the extent of produced pollutants in one hour of the MG's operation is shown by Emission t [5]: In the equations above, CO t are the levels of CO 2 , SO 2 and NO x emissions from the ith DG source; CO t 2 Batt,s , SO t 2 Batt,s and NO t x Grid are the amounts of CO 2 , SO 2 and NO x emission from the sth battery unit and CO t 2 Grid , SO t 2 Grid and NO t x Grid are the amounts of CO 2 , SO 2 and NO x emission from the utility at hour t, respectively [5].

Power Balance Constraint
The balance of electricity demand and power supply is one of the most important requirements in MG energy management [33]; hence: where P t LD is the total MG load at hour t.

Battery Limits
The constraints of the battery SOC limits during each time interval are considered as follows [34]: where W t Batt s is the SOC of the sth battery at the end of one-hour interval, which is associated with the time t as follows: Batt and E t−1 Batt are the energies stored inside the battery at hours t and t − 1, respectively. The lower and upper limits of the stored energy inside the sth battery are E Batt s,min and E Batt s,max . Additionally, E t Batt should not violate the following limits:

Real Power Constraints
Power generations of each dispatchable DG are limited as follows: The power exchange with the utility grid is constrained as follows: The following shows the constraints of the battery rates of charge and discharge during each hour:

Reduced Unscented Transformation (RUT)
Wind power generation, solar irradiation, load demand and market price are the most important, uncertain parameters to be considered for solving the problem of MGs' optimal operation. In this paper, the RUT method is used to model uncertainties based on [29]. The uncertain variables follow a normal distribution function with mean values equal to those expressed in [29] and standard deviations equal to 10% of their mean values. The main steps followed by the RUT method are summarized in the following [29].
Step 1: Choose 0 ≤ W 0 ≤ 1: Step 2: Choose weight sequence according to the following criterion: Step 3: Initialize vector sequence: Step 4: Expand vector sequence (ξ j k ) for j = 2, . . . , m according to the following: If a random variable vector → Z has a mean → µ z and a covariance matrix P zz , the kth sigma point is: where √ P zz is a matrix square root of P zz . For a positive definite matrix P, the matrix square root means that a matrix A = √ P exists such that P = AA T , which should be calculated using numerically efficient and stable methods such as Cholesky decomposition [35]. If there is no correlation between uncertain variables, the elements of P zz 's main diameter will be the square of uncertain variables' standard deviations (σ). However, in correlated conditions, depending on which uncertain variables are correlated, the corresponding elements of the rows and columns of P zz may have negative or positive values. Accordingly, the P zz matrix can be obtained as follows: where γ α,β is the correlation coefficient between the αth and βth elements of the covariance matrix P zz.
Step 5: Input m + 2 sigma points to the function to find the output samples: Step 6: Calculate the mean → µ y and covariance matrix P yy of the output → Y as follows:

A Brief Overview of the Original JAYA
The JAYA algorithm, which was proposed by "R. Venkata Rao" [28], is enhanced and applied in this article. The inspiration of the original JAYA algorithm is derived from the Sanskrit word "JAYA" the meaning of which is "victory". In the JAYA algorithm, each individual moves toward the best solution while trying to move far away from the worst in order to achieve the best solution. As the methodology is demonstrated in the following, due to the simplicity and rapidity of the algorithm, and since there is no need to have any information regarding control parameters, JAYA becomes more advantageous compared to other meta-heuristic algorithms [28,29].
The goal is to minimize F( → X) as the objective function. Assume that the best solution, 'X best ', and the worst solution, 'X worst ', obtain the best and the worst values of F( → X) in the entire candidate solutions, respectively. Let → X new i,iter be the ith candidate during the iterth iteration. Then, to modify the decision variables, individuals could be changed according to the following [28]: where r 1 and r 2 are random variables in [0, 1], → X best is the best solution in each iteration, while → X worst is the worst solution in the population; → X new i,iter will be accepted if its objective function value is better than that of → X i,iter .

Multi-Objective JAYA (MOJAYA)
Since, in this paper, a multi objective problem is implemented, new definitions should be considered for 'X best ' and 'X worst '. For this purpose, the initial population should be organized according to the following [35]: where f 1 and f 2 are the cost and emission of the objective functions, respectively. The initial guess for ω 1 and ω 2 is equal to 0.5 [36]. The worst solution is the most dominant solution in each iteration and is then selected as 'X worst '. In order to select the best solution, in each iteration, after choosing the non-dominated solutions and controlling the size of the repository according to Algorithm 1, the best solution is selected from the repository as 'X best ' based on (33) and is used in the subsequent iterations.

Enhanced MOJAYA (EMOJAYA)
When the dimensions of the problem increase, the original JAYA may not find the global optimum solution easily. Consequently, in order to achieve a real and well-scattered Pareto optimal front, and for improving the performance, convergence ability and the accuracy of the algorithm, some modifications are necessary. In the suggested EMOJAYA algorithm, the main superiority of the original JAYA, which is independency to algorithmic parameters, is also preserved. The modification process of the algorithm is as follows: Step 1: Set i = 1 Step 2: Select the ith individual from the initial population.
Step 3: Choose five constants k 1 = k 2 = k 3 = k 4 = k 5 , unequal to i, randomly from the initial population, and select three constants k 1 = k 2 = k 3 from the repository.
Step 4: Calculate four new individuals as follows: Step 5: Calculate objective functions for the new individuals.
Step 6: Determine and save non-dominated solutions.
Step 7: Update X i from the new non-dominated solutions of Step 6 using (32).
Step 8: If all individuals are updated, stop; otherwise, set i = i + 1 and return to Step 2.

Application of the Suggested EMOJAYA Algorithm
In order to apply the proposed EMOJAYA algorithm on the considered MOOM problem, the following steps should be taken: Step 1: Initialize the population size, number of design variables and termination criterion. Problem information including MG properties, as well as bids and power information of DGs, storages and utility and emission coefficients shall be specified. The initial charge of the battery is also defined in this step.
Step 2: Hourly WT and PV power forecasts, demand and electricity price are determined through RUT.
Step 3: Since a mixed-integer problem is considered in this study, two types of variables, i.e., binary and continuous, are assumed. For states of generators as binary variables, u i s according to (5) are generated as follows: However, in order to consider the states of dispatchable units, U t i should satisfy the following condition for all hours: If the equation above is satisfied, a random population for continuous variables based on the achieved u i s and according to (5) must be generated as follows: while the power constraints in (12)- (14) should be satisfied.
Step 4: Since in the considered problem, some limitations, such as battery constraints, depend on previous and future hours, constraints change in different hours of the day as follows: For discharging mode: For charging mode: where for both discharging and charging modes: P t ch,max = max (P Batt,min , The limitation of the released energy in the discharging mode and the stored energy in charging mode are mentioned in (27) and (28), respectively, while (30) and (31) demonstrate the power discharged by the battery in the discharge mode and the power charged by the grid in the charging mode. Using (30) and (31), the maximum and minimum rates of charging and discharging modes can be calculated [16].
Step 5: Calculate the objective functions for the population.
Step 6: Determine the non-dominated solutions in the population and store them in the repository.
Step 7: Choose 'X best ' and 'X worst ' as described in Section 5.1.
Step 8: Update the initial population as Section 5.
Step 9: Determine the non-dominated solutions in the new population.
Step 10: Update the repository based on new and old non-dominated solutions.
Step 11: If the number of non-dominated solutions is more than a predefined value, N L , control the size of the repository according to Algorithm 1.
Step12: Control the termination criterion (iteration = iteration max ). If it is satisfied, terminate the algorithm; otherwise, set iteration = iteration + 1 and return to Step 7.
Step 13: Report the BCS that is achieved using the technique proposed in Section 5.1.

Simulation Results
In order to justify the efficient performance and competitiveness of the proposed EMO-JAYA algorithm, two subsections are investigated in this section. First, the performance of the suggested EMOJAYA is compared with the original JAYA and PSO algorithms by applying these algorithms (40 independent runs of each algorithm are performed, and the results are reported) to three test functions of [39]. In the second sub-section, three different scenarios are considered to inspect the effectiveness of the EMOJAYA algorithm in solving MG energy management problems in the presence of storage devices.

A Comparison between the Performance of Proposed EMOJAYA with Original JAYA and PSO Algorithms on Different Test Functions
In this section, three different test functions of [39] are considered and the proposed EMOJAYA algorithm is applied to solve the problems.
-Test Function 1: The first test objective functions of [39] are as follows: The Pareto front of the suggested EMOJAYA and those of the original JAYA and PSO are compared in Figure 2, and the prior performance of EMOJAYA in solving the first chosen test function is observed. Electricity 2022, 3, FOR PEER REVIEW 12 section, three different scenarios are considered to inspect the effectiveness of the EMO-JAYA algorithm in solving MG energy management problems in the presence of storage devices.

A Comparison between the Performance of Proposed EMOJAYA with Original JAYA and PSO Algorithms on Different Test Functions
In this section, three different test functions of [39] are considered and the proposed EMOJAYA algorithm is applied to solve the problems.
-Test Function 1: The first test objective functions of [39] are as follows: The Pareto front of the suggested EMOJAYA and those of the original JAYA and PSO are compared in Figure 2, and the prior performance of EMOJAYA in solving the first chosen test function is observed.
-Test Function 2: The second proposed benchmark in [39] is as follows: where:   -Test Function 2: The second proposed benchmark in [39] is as follows: where: and 0 ≤ x 1 ≤ 1, −30 ≤ x 2 ≤ 30 Figure 3 shows the performance and the resulting Pareto front of EMOJAYA, original JAYA and PSO algorithms for solving the second test case study. It is observed that the suggested EMOJAYA has a better performance compared to the two other algorithms.
Electricity 2022, 3, FOR PEER REVIEW 13 Figure 3 shows the performance and the resulting Pareto front of EMOJAYA, original JAYA and PSO algorithms for solving the second test case study. It is observed that the suggested EMOJAYA has a better performance compared to the two other algorithms.

-Test Function 3: The third test function is [39]:
and In Figure 4, the Pareto front of the proposed EMOJAYA and the results of the original JAYA and PSO are demonstrated. It can be concluded that EMOJAYA has a better performance compared to PSO and JAYA algorithms.
From Figures 2-4, and by comparing the Pareto optimal fronts of the investigated algorithms, it is clear that the proposed EMOJAYA has an efficient performance in solving multi-objective problems.  -Test Function 3: The third test function is [39]: and 0.1 ≤ x 1 ≤ 1, 0.1 ≤ x 2 ≤ 1 In Figure 4, the Pareto front of the proposed EMOJAYA and the results of the original JAYA and PSO are demonstrated. It can be concluded that EMOJAYA has a better performance compared to PSO and JAYA algorithms.  Figure 3 shows the performance and the resulting Pareto front of EMOJAYA, original JAYA and PSO algorithms for solving the second test case study. It is observed that the suggested EMOJAYA has a better performance compared to the two other algorithms.

-Test Function 3:
The third test function is [39]: and In Figure 4, the Pareto front of the proposed EMOJAYA and the results of the original JAYA and PSO are demonstrated. It can be concluded that EMOJAYA has a better performance compared to PSO and JAYA algorithms.
From Figures 2-4, and by comparing the Pareto optimal fronts of the investigated algorithms, it is clear that the proposed EMOJAYA has an efficient performance in solving multi-objective problems. From Figures 2-4, and by comparing the Pareto optimal fronts of the investigated algorithms, it is clear that the proposed EMOJAYA has an efficient performance in solving multi-objective problems.

Solving the MG Energy Managemnet Problem
In this section, three different scenarios are investigated to verify the effectiveness of the proposed RUT-EMOJAYA algorithm in solving the probabilistic MG energy management problem in the presence of storage devices. The MG of Figure 1 is considered as the test system. A 24-hour scheduling scheme is assumed for the analysis of the simulated system. For all DGs, the unity power factor is considered; thus, they only produce active power. The MGCC is responsible for controlling the power exchange between the MG and the utility, which is allowed at any hour in a day in order to more profitably exploit the market. The hourly active power of PV and WT, the forecasted load demand, and the utility power production bid are shown in Figures 5 and 6. It is assumed that all loads are electrical [40].

Solving the MG Energy Managemnet Problem
In this section, three different scenarios are investigated to verify the effectivene the proposed RUT-EMOJAYA algorithm in solving the probabilistic MG energy man ment problem in the presence of storage devices. The MG of Figure 1 is considered a test system. A 24-hour scheduling scheme is assumed for the analysis of the simu system. For all DGs, the unity power factor is considered; thus, they only produce a power. The MGCC is responsible for controlling the power exchange between the MG the utility, which is allowed at any hour in a day in order to more profitably explo market. The hourly active power of PV and WT, the forecasted load demand, and utility power production bid are shown in Figures 5 and 6. It is assumed that all load electrical [40].
In Table 1, the entire bid data for all DGs and utility grid, as well as the amou emission of DGs, are available [41]. Since PV and WT units do not consume any fuel w they produce electric power, the utility grid has to buy all electric power produce these units [5]. Since the two considered objective functions (cost and emission) are flicting, a Pareto-optimal set is achieved for each considered case.
The proposed method was implemented in MATLAB. The number of both pop tions and maximum iterations is considered to be 100.  The on\off states of dispatchable DGs (i.e., FC and MT) are taken into account i considered MG's energy management problem. Therefore, in the algorithm process,  In Table 1, the entire bid data for all DGs and utility grid, as well as the amount of emission of DGs, are available [41]. Since PV and WT units do not consume any fuel when they produce electric power, the utility grid has to buy all electric power produced by these units [5]. Since the two considered objective functions (cost and emission) are conflicting, a Pareto-optimal set is achieved for each considered case. Accordingly, the comparison of cost and emission objective functions in the BCS point for different efficiencies, including 93%, 95% and 97%, are shown in Figures 7-9. Moreover, the Pareto optimal front for different efficiencies, along with their BCS points are illustrated in Figure 10. It should be mentioned that, in Figure 10, results are compared for the case that E is equal to   The proposed method was implemented in MATLAB. The number of both populations and maximum iterations is considered to be 100.
The on\off states of dispatchable DGs (i.e., FC and MT) are taken into account in the considered MG's energy management problem. Therefore, in the algorithm process, if the solutions for dispatchable units are lower than their minimum limit powers, then the power will be set to zero.
Three scenarios are studied in this section to investigate the effect of storage devices on solving the MOOM problem. In the two first scenarios, it is assumed that the battery has a degradation cost, as in (3), while in the third scenario, a comparison between the situations with and without considering the battery degradation cost is carried out.
Scenario I: In this case, the battery efficiencies are considered to be 93%, 95% and 97%, taking the influence of the Li-ion battery efficiency into account. Additionally, the impact of different initial and final charges is contemplated in the three aforementioned efficiencies. In each case, the battery's initial and final charges (initial charge = final charge = E) differ from E Batt,max to 0.2E Batt,max so that this variation effect can be investigated. Accordingly, the comparison of cost and emission objective functions in the BCS point for different efficiencies, including 93%, 95% and 97%, are shown in Figures 7-9. Moreover, the Pareto optimal front for different efficiencies, along with their BCS points are illustrated in Figure 10. It should be mentioned that, in Figure 10, results are compared for the case that E is equal to E Batt,max . Electricity 2022, 3, FOR PEER REVIEW 1 Figure 6. Estimated values of (a) load demand and (b) utility hourly prices [40].

Scenario II:
In the second scenario, in order to obtain the best value of the battery initial and final charges (initial charge = final charge = E), they are considered as decision variables. Figures 11-13 show comparisons of the Pareto optimal fronts for different E values in different efficiencies (i.e., 93%, 95% and 97%). It can be observed from Figures  11-13 that when E is considered as a decision variable, a better Pareto optimal front is achieved in all efficiencies and, consequently, the obtained BCS will be more convenient.

Scenario II:
In the second scenario, in order to obtain the best value of the battery initial and final charges (initial charge = final charge = E), they are considered as decision variables. Figures 11-13 show comparisons of the Pareto optimal fronts for different E values in different efficiencies (i.e., 93%, 95% and 97%). It can be observed from Figures  11-13 that when E is considered as a decision variable, a better Pareto optimal front is achieved in all efficiencies and, consequently, the obtained BCS will be more convenient.

Scenario II:
In the second scenario, in order to obtain the best value of the battery initial and final charges (initial charge = final charge = E), they are considered as decision variables. Figures 11-13 show comparisons of the Pareto optimal fronts for different E values in different efficiencies (i.e., 93%, 95% and 97%). It can be observed from Figures 11-13 that when E is considered as a decision variable, a better Pareto optimal front is achieved in all efficiencies and, consequently, the obtained BCS will be more convenient. Electricity 2022, 3, FOR PEER REVIEW 18 Figure 11. Comparison of Pareto optimal fronts for different E values (93% efficiency).

Scenario III:
In the third scenario, the battery degradation cost is neglected and results are compared with the situation in which the battery degradation cost is taken into account as Equation (3). In this case, the battery efficiency is considered equal to 95%. The comparison of Pareto optimal fronts with and without considering degradation cost is given in Figure 14. According to this figure, as already expected, the Pareto optimal front reveals a drop in cost objective function when the degradation cost is ignored.

Scenario III:
In the third scenario, the battery degradation cost is neglected and results are compared with the situation in which the battery degradation cost is taken into account as Equation (3). In this case, the battery efficiency is considered equal to 95%. The comparison of Pareto optimal fronts with and without considering degradation cost is given in Figure 14. According to this figure, as already expected, the Pareto optimal front reveals a drop in cost objective function when the degradation cost is ignored.
The advantage and perfection of an optimization algorithm for solving multi-objective problems is decided based on its effective and fast convergence, resulting in a smooth and accurate Pareto front. According to the simulation results, it is confirmed that the presented EMOJAYA algorithm successfully satisfies these criteria.

Conclusions
The probabilistic economic/emission optimal operation of a typical MG was dealt with in this paper using the new RUT-EMOJAYA algorithm. In order to investigate the role of storage devices, as well as the effective performance of the suggested algorithm, three different scenarios were studied. In the first case, different efficiencies for the storage device were considered, while the battery's initial and final charges were not included in the decision variables. Simulation results showed that, by increasing the battery efficiency from 93% to 97%, the daily operational cost of the MG system was reduced from EUR 9200 to EUR 8900. In the second case, the battery's initial and final charges were considered as decision variables, and it was observed that the Pareto optimal front in this case was better distributed than in the previous case. In the third case, the degradation cost of the storage device was ignored in order to compare the results with those of other two cases where The advantage and perfection of an optimization algorithm for solving multi-objective problems is decided based on its effective and fast convergence, resulting in a smooth and accurate Pareto front. According to the simulation results, it is confirmed that the presented EMOJAYA algorithm successfully satisfies these criteria.

Conclusions
The probabilistic economic/emission optimal operation of a typical MG was dealt with in this paper using the new RUT-EMOJAYA algorithm. In order to investigate the role of storage devices, as well as the effective performance of the suggested algorithm, three different scenarios were studied. In the first case, different efficiencies for the storage device were considered, while the battery's initial and final charges were not included in the decision variables. Simulation results showed that, by increasing the battery efficiency from 93% to 97%, the daily operational cost of the MG system was reduced from EUR 9200 to EUR 8900. In the second case, the battery's initial and final charges were considered as decision variables, and it was observed that the Pareto optimal front in this case was better distributed than in the previous case. In the third case, the degradation cost of the storage device was ignored in order to compare the results with those of other two cases where this cost was taken into account. It was seen that the Pareto optimal front, without considering the degradation cost, demonstrated a decrease (USD 270) in the minimum achievable daily operation cost of the MG system, which is expected to be true. It was concluded that the proposed method for dealing with mixed-integer optimization problems in the presence of uncertainties is accurate. It should be pointed out that as the RUT-EMOJAYA algorithm is comparably simple and not dependent on control parameters, its application is recommended for clarifying other problems. Future studies can include the following: i.
Investigating elements of the future smart grids, including demand response and the influence of electric vehicles on the considered MG's energy management problem. ii. Inspecting reliability as an objective function in the MG's optimal operation management.
iii. Comparing different energy storage devices, as well as a variety of battery technologies to decide on the most optimal economic design of the system.  amount of pollutants emission for storage device at hour t (kg/MWh) E t Grid amount of pollutants emission for the utility grid at hour t (kg/MWh) W t Batt amounts of energy stored inside the battery at hour t P ch (P disch ) permitted rate of charge (discharge) of the battery η c (η d ) efficiency of the battery during charge (discharge) process W Batt s,min lower limit of amounts of energy storage inside the battery W Batt s,max upper limit of amounts of energy storage inside the battery P ch,max (P disch,max ) maximum rate of battery charge (discharge) during each time interval ∆t