Optimal Multi-Operation Energy Management in Smart Microgrids in the Presence of RESs Based on Multi-Objective Improved DE Algorithm: Cost-Emission Based Optimization

: Today, in various leading power utilities in developing countries, achieving optimal operational energy management and planning, taking into account the costs reduction of generation, transmission and distribution of electricity, and also reducing the emission of an environmental pollutant becomes more and more important. Optimal use of renewable energy sources (RESs) is an effective way to achieve these goals. In this regard, in this research article, an improved multiobjective differential evolutionary (IMODE) optimization algorithm is suggested and performed to dispatch electricity generations in a smart microgrid (MG) system, taking into account economy and emission as competitive issues. In this paper, a nonlinear equation of multi-objective optimization issue with various equality and inequality limitations is formulated in order to lower the total operational costs of the MG considering environmental pollution effects simultaneously. In order to address the issue of optimal operation of the MG in single-objective and multi-objective forms, an intelligent method according to the improved differential evolutionary (IDE) optimization is utilized and performed and the proposed algorithm is implemented on different problems. First, it is assumed that there is no limit to the exchange of power overhead, and secondly, the limitation of power exchange with the upstream grid is considered. In multi-objective mode, these two modes are also considered. In order to show the impact of renewable energy on the cost, in the third part of the simulations, the operation is solved with maximum participation of renewable energy sources. In the ﬁnal section, the sensitivity analysis on the number of populations in this problem is performed. The obtained results of the simulation are compared to differential evolutionary (DE) and particle swarm optimization (PSO) techniques. The effectiveness of the suggested multi-operational energy management method is conﬁrmed by applying a study case system.


Introduction
In recent years, novel methods have been conducted in power system utilities towards distributed generation (DG) in energy conversion systems (ECSs) which mean that they are located close to the electricity customers; besides, large power plant units have been substituted by smaller ones. In addition to this, for consumers, the potential of lower costs, higher reliability of services and power quality, increasing in power efficiency and power independency can be the reasons for higher interests in DGs. Using RESs and green powers, based utilizing metaheuristic algorithms have been suggested by Sadegheih [34]. Similarly, authors of the paper [35] have surveyed different methods like genetic algorithm (GA), in that parameters have been regulated with fuzzy logic (FL) concepts. Additionally, in order to solve an economic and environmental dispatch (EED) issue, an improved decision shuffle frog leaping algorithm has been used in reference [36]. In paper [37], also, to address the EED issue, a stochastic method taking into account uncertainties was proposed; and in [38], authors performed a detailed research regarding evolutionary programming algorithm i.e., ant colony optimization (ACO) and GA to solve EED problems have been implemented. For optimal dispatch and operating electric energy, various multi-objective evolutionary algorithms were also suggested in the reference [39]. In reference [40], authors focused on the benefits of using neural-fuzzy methods to find an acceptable way to address multi-objective generation and energy management issue. The impacts of electric vehicles, smart homes, and smart buildings, including smart hospitals and smart schools, have been examined, optimized, and managed in references [41,42].
Among the above-mentioned optimization techniques, differential evolutionary (DE) algorithm has been noticeably utilized in optimal operation and energy management issue, basically because of its population-based search ability and higher efficiency, more speed of convergence, and robustness; therefore, an improved multi-objective differential evolution (IMODE) technique is suggested in this research. To dominate the issues of falling into the local optimum and premature convergence, in our proposed IMODE algorithm, an immune operation technique is introduced to power the local search capability and improve optimization precision. The suggested method is successfully carried out on various case studies, and the efficiency measures consisting of convergence and divergence metrics are computed. Besides, the multi-objective optimization of a multi-operation energy management in SGs in presence of RESs procedure is performed using the presented IMODE algorithm. In order to show the impact of renewable energy on the cost, in the simulations, the operation is solved with maximum participation of RESs. In the final section, the sensitivity analysis on the number of populations in this problem is performed.

Materials and Methods
In this section, the problem of optimal operation, energy management, and power dispatch in MGs considering environmental and economic indicators will be implemented as a multi-objective optimization formulation. In addition, MG modeling and proposed improved differential evolutionary (IDE) optimization algorithm are introduced. Within the multi-objective optimization process, both opposing goals should be considered, including the total operating expense of the MG and the air emission of pollutants, which have to be minimized simultaneously while satisfying the limitations of the system. Hence, we should formulate these problems as the following functions.

Objective Functions (OFs)
The problem of optimizing energy and operation management can be defined as a multi-objective optimization problem by considering both environmental and economic criteria in a smart MG at the same time. In this regard, two different objectives, one being the operating costs of the network and the other being the emission of environmental pollutants, should be minimized while meeting the limitations of the system. Hence, the mathematical forms of this issue can be stated in the following subsections: The total expenses of operation in MG consist of the fuel costs of the power plant units and their costs of start-up and shut-down. The result of this part is a collection of optimal load flow from energy resources to the loads economically for a specific time period. Storage system options, in addition to distributed generations (DGs), can be utilized to compensate expensive power buying from the company and storing power in off-peak Appl. Sci. 2021, 11, 3661 4 of 28 hours to predict the increasing projected prices. In this way, we can formulate the first OF as follows: where B Gi (t) represents the bids from the DG resources, B sj (t) is the storing options at hour t, S Gi is the start-up expenses for ith DG and jth storage component, S sj gives the shut-down expenses for ith DG and jth storage component; P Grid (t) also represents the active power that is purchased/sold from/to the companies at time t; B Grid (t) represents the bid of company at t; X also gives the vector of state variables such as active powers of power plants and their related states which can be defined as the following equations: where in the Equation (2), T represents the total hours number in the tested time period; n gives the state variables number; N g also represents the DGs' total number, and N s gives the storage total number of power plant units; P g represents the all DG units' active power; U g also gives the state vector defining the position of ON/OFF states for all power plant units in every hour of a day. All of these variables are given in Equations (3) and (4). P G = P G1 , P G2 , . . . , P G,N g P Gi = [P Gi (1), P Gi (2), . . . , P Gi (t), . . . , P Gi (T)]; i = 1, 2, . . . , N g + 1 P s = [P s1 , P s2 , . . . , P s,N s ] P sj = P sj (1), P sj (2), . . . , P sj (t), . . . , P sj (T) ; j = 1, 2, . . . , N s (3) where in the Equation (3), P Gi (t) and P sj (t) represent the outputs power for ith generator and jth storage at time t.
where u k (t) gives the unit status k at hour t.

Objective 2: Minimizing the Total Amount of Air Emissions in MG
The environmental emissions impact of pollutants as the second objective should be considered in this study, where 3 of the most serious air pollutants in the optimization problem including carbon dioxide (CO 2 ), Sulphur dioxide (SO 2 ) and nitrogen oxides (NO x ) are involved [43]. In addition, logical functions which link air pollutants to the energy generation of various DGs are used to provide a model that illustrates such effects; similarly, it can be utilized for economic power dispatch. In this regard, the second OF will be given as Equation (5): where, E Gi (t), E Si (t), and E Grid (t)} represent the level of air emissions in kg/MWh for every distributed generation, storage module and the utility at hour t. All of variables are presented in the Equations (6) to (8): where, CO 2DG i is the amount of carbon dioxide, SO 2DG i (t) gives the amount of Sulphur dioxide, and NO xDG i (t) represents the amount of emissions for nitrogen oxides from ith DG unit during hour t.
where, CO 2ndorage j (t) is the amount of carbon dioxide, SO 2ndorage j (t) gives the amount of Sulphur dioxide, and NO x (t) represents the number of emissions for nitrogen oxides from jth storage unit at each hour of a day.
where, CO 2Grid (t) is the amount of carbon dioxide, SO 2Grid (t) gives the amount of Sulphur dioxide, and NO xGrid (t) represents the number of emissions for NO x from the utility or the MG at hour t.

Load Balance
In order to have an optimal energy management in energy system operational state, the principal limitation should be to satisfy loads with available power generations. Thus, the active power balancing condition for the MG should be formularized as follows: where P Lk gives the amount of kth load, also N k represents the load levels total number.

Active Power Limitations of the Units
The active power limitations of the power plant units are formulized as following equations: P Gi, min (t) ≤ P Gi (t) ≤ P Gi, max (t) P si, min (t) ≤ P si (t) ≤ P si, max (t) P grid, min (t) ≤ P Grid (t) ≤ P grid, max (t) (10) where in Equation (10), P G, min (t) represents the minimum active powers of ith DG, P s, min (t) is the minimum active powers of jth storage, and P grid, min (t) represents the minimum active powers of the utility, at time t. Similarly, P G, max (t), P s, max (t), and P grid, max (t) give the maximum energy generations of corresponding units during hour t.

Charge and Discharge Level Constrains of Storage System Modules
Owing to constrain of the charge and discharge level of storing systems at every time interval, the formulations and limitations need to be considered as follows: W ess,min ≤ W ess,t ≤ W ess,max P charge,t ≤ P charge, max ; P discharge, t ≤ P discharge, max (12)  where, W ess,t and W ess,t−1 represent the amounts of power in battery storage system at hour t and hour t − 1, respectively; also P charge and P discharge give the levels of permitted charge and discharge via a definite time period ∆t. η charge and η discharge indicate the efficiency of the charge/discharge in the battery storage system. W ess,min and W ess,max represent the minimum and maximum constraints of battery storage system while P charge, max and P discharge, max give the maximum levels of charge or discharge state in every interval ∆t.

MG Modeling
Since using DGs can have different effects on power system characteristics, application and modelling of DGs can cause various issues. One possible solution to maximize the benefits of using DGs is to have an effective method to optimally model the generation of relevant power and loads in the network as a subsystem in a form of MG. Additionally, by integrating DGs into a MG and using RESs in large quantities, issues related to the optimization of cost and environmental pollutants reduction can be investigated in the targeted system. Such strategies can make appropriate decisions to a better energy management and enhance the performance of operation.
In addition, using DGs in forms of RESs such as fossil fuel engines, gas and micro turbines (MTs), photovoltaics (PVs) and fuel cells (FCs), and wind turbine (WT) in MG needs higher and complicated technologies, which can significantly lower emissions and air pollution. For instance, FCs which generate power from oxygen and hydrogen are emitting only vapor. However, at the same time of modification of natural gas, various forms of NO x and CO 2 air pollutants are emitted. Another example could be the use of MTs with lower emissions; however, currently they are more expensive.
In this research work, a typical low-voltage (LV) smart-MG is considered which is shown in Figure 1 consisting of different DGs including MTs, FCs, PVs, WT, and storage system devices (batteries).

Principles of IDE Algorithm and Multi-Objective Optimization
In general, real-world optimization issues involve various goals that should be optimized simultaneously, and this set of goals might have conflicting and different goals and incomparable applications [44,45]. In order to simultaneously address these issues to- In this regard, it is presumed that all DG resources generate active power at the power factor one, without any requesting or generating reactive power. In addition, there is a link for power exchange from the presented MG and the grid utilized for power trading in various hours of a day according to decisions of MG control center (MGCC).

Principles of IDE Algorithm and Multi-Objective Optimization
In general, real-world optimization issues involve various goals that should be optimized simultaneously, and this set of goals might have conflicting and different goals and incomparable applications [44,45]. In order to simultaneously address these issues towards optimal aims, multi-objective optimization techniques can be used effectively to obtain a set of optimal answers. In this section, the principles of our proposed IDE algorithm to optimize energy management in MG are presented; hence, we need to present Equation (13): In Equation (13), F gives the objective functions vector, and X represents the optimization variables vector, f i (X) also represents the ith objective function, g i (X) is the equality constraint and h i (X) represents inequality constraints; m also gives the number of objective functions.

Differential Evolutionary Algorithm
DE algorithm, today, acts one of the highly efficient alternatives to optimize complicated multi-objective issues. DE algorithm dives a fast rate of convergence, robustness with a global optimization capability. To generate the vectors of trial parameter, the principle operation of this algorithm consists of the mutation, crossover, and selection operations. The principal operating process of the basic DE has been provided in reference [46]. Of all the operational above-mentioned items, mutation operation can act as the core of DE algorithm. Several and various mutations mechanisms have been presented in the literature review and the frequently utilized five different techniques which are listed in the following equations: DE/current to best/1: where, t = 1, 2, 3, . . . represents the generation (time), indices r i (i = 1, 2, . . . , 5) give the mutually exclusive integers which generated randomly through [1 ,NP], that differ from index i; x t bset represents the best vector that denotes the best amount of fitness of the population in generation t; also F gives the scale factor. The crossover operation also consists of the binomial crossover approach and the index cross approach. The most commonly utilized bi-nominal crossover operation has been used in this research work. Hence, we should define the following equation as: where sn ∈ [1, n] represents a random integer; u t i, j gives the jth part of the ith new trail vector, also CR denotes the rate of crossover through the range [0, 1]. In addition, if the experimental vector performs better than the target vector, in the next generation, the experimental vector will be replaced by the greedy selection vector in the selection operation with the target vector.
In the proposed strategy, the degree of solution limitation violation x on the jth limit can be obtained using Equation (21).
The general limitation degree of solution x can be computed using Equation (22).

Improved Self-Adaptive Differential Evolutionary Algorithm
In this section, the proposed improved IDE has been introduced and the adaption mechanism has been described.

Adaptive Regulation Mechanisms for Scale Factor and Crossover Rate
Selecting F and CR depends highly on the specified issue that was applied, and that is almost difficult. Overall, small values of F and CR lead to a faster convergence speed, whereas large values of F and CR contribute to rise the population diversity and can prevent premature convergence. In order to obtain a balance from the speed of convergence to premature convergence, having an adaptive regulation mechanism for F and CR will be essential. So, to develop the global searching ability with escaping the local optimal and preventing the premature convergence, an improved approach and technique is proposed to develop the total efficiency of the IMODE algorithm in this research. In the proposed algorithm, a heuristic searching optimization approach according to the basics of genetics and immune adaptive (IA) method is proposed. During the process of population search, IA maintains the potential solutions populations. As a result, in order to reach the parameters of the optimal scale factor, the principles of immune approach are used. The strategies of parameter F for updating are adjusted as the following steps:

1.
Computing values of affinity: IA utilizes the affinity as a differentiator between the qualities of solutions; besides, antibodies with a higher affinity are more likely to be opted, and at the higher probability, are more likely to be surviving for the next generation. Hence, we describe the affinity in Equation (23).
where, the coefficient r ∈ [0, 0.5]; i also by arranging them based on their performance values, shows the numerical order of every antibody in the population. In order to calculate the values of affinity for the parameters of scale factor, we define Equation (24) as follows: Based on their affinity, the parameters of scale factor should be arranged between the best and the worst, where the best twenty percent mutation parameters MF T+1 (new) have been stored in the memory based on MF.

2.
Updating the MF memory based on MF: 3.
Generating the scale factor parameters for the next generation. Now, the scale factor parameters in the next generation should be updated as follows: where the initial F T+1 j has been generated by a normal distribution N(0.5, 1) through [0, 1] 2 , P i represents a random number in the distance within [0, 1], and P F represents the regulating probability. Considering these items and based on the reference [47], the crossover rate (CR) values have been updated as: where, τ = 0.9 represents the adjusting probability in our proposed method.

Limitation-Handling Technique
Because optimization issues in the real-world are often connected to various limitations, the alpha (α)-finite approach has been utilized to address limitations. The α-finite technique treats several infeasible solutions as feasible approaches whenever such nonapplicable solutions are under a specific level. In order to control the level α, Equation (28) is presented.
where x θ represents the top θth solution at the initialization phase; cp denotes a parameter for controlling the speed of limitations reducing relaxation; θ has been suggested to be 0.2 * NP; also cp ∈ [2,10]. For having an efficiently using of the information performed by the infeasible solutions, the concept of pseudo feasible has been described. Generally, the pseudo feasible solution gives the overall solution whose general limit violation becomes less than α. A detailed study about the limited-dominate connection between individuals consisting of pseudo feasible and pseudo infeasible solutions has been presented in the reference [48]. After the phase of sorting operation based on the limited-dominate connection criteria, exactly NP solutions should be selected from combined population X t ∪ U t to new population X t+1 . By combining the self-adaptive evolutionary method and immune operation, the step by step pseudo code of the new improved self-adaptive evolutionary is presented as below: Step 1: Set the parameters' values and initialize the external archive.
Step 2: Set the generation number G = 0 randomly initialize a population of NP individuals uniformly distributed in the range [X min ,X max ] Step 3: Evaluate the fitness of initial solutions Step 5: Divide the combination of parent and offspring population into two subpopulations according to pseudo feasible concept, SubP1 (pseudo feasible solutions) and SubP2 (pseudo infeasible solutions) Step 5.1: Order the SubP1 by non-dominated sorting procedure Step 5.2: Order the SubP2 according to overall violation degree Step 6: Add all the non-dominated solutions into external archive and remove the inferior ones.
Step 7: Select solutions for the main population of the next generation Step 8: Update the parameters Step 9: Increment the generation count G = G + 1.

END WHILE
All simulations are performed in a typical system that has MT resources, PVs, WTs, FCs, ES units with having the capability of upstream grid power exchange.

Simulation Results
For demonstrating the efficiency of the IDE algorithm in solving the problem of optimal operation of the MG in both single-objective and multi-objective modes, this algorithm has been performed on different problems. First, it is assumed that there is no limit to the exchange of power overhead, and secondly, the limitation of power exchange with the upstream grid is considered. In multi-objective mode, these two modes are also considered. In order to illustrate the impact of RE generations on the cost, in the third part of the simulations, the operation is solved with maximum participation of renewable energy sources. In the final section, the sensitivity analysis on the number of populations in this problem is performed. The figures for gas emissions of resources, costs, and minimum and maximum production are available in Table 1. The rate of daily loading is available in Table 2. Additionally, Table 3 shows the market price for upstream power exchange. Forecasts for photovoltaic and wind turbine sources are available, and their maximum output is given in Table 4. The cost coefficients ($/kWh) and emissions (kg/MWh) assumed by the DG resources are shown in Table 1.  For solving the single-objective problem in the unlimited power exchange mode with the network, both cost and pollution (gas emission) targets have been studied separately. The results are compared with the PSO and DE methods and the convergence diagrams are presented in Figures 2 and 3. Both forms show better convergence of the IDE algorithm compared to the other methods. Table 5 to Table 6 illustrate how to distribute power to reduce costs and gas emission without considering the limitations of upstream network power exchange with two objectives. The last column of the tables illustrates the unlimited power exchange with the upstream utility grid. Tables 7 and 8 show the results of 10 consecutive runs of each method. The proposed algorithm is better than all other methods in all the best, worst, and average of the indices, and therefore, the superiority in both targets is quite clear.        In this section, the power exchange with the utility grid is considered as finite. The distributions of power for different purposes are provided in Tables           In this section, in order to solve the problem in a multi objective format, the limitation of power exchange with the upstream network is not considered. Figure 6 shows that the different approaches have come to close conclusions and the Pareto Front is very close to all three methods. Where NSIDE is non-dominated solution improved differential evolution, NSDE defines non-dominated solution differential evolution and NDPSO is non-dominated solution particle swarm optimization. Table 13 also shows the output of a sample response from the Pareto front of the NSIDE method. While the power generation capacity of all units is within the permissible range, it can be observed that at some times the network has largely purchased or sold electricity from the upstream network, which was due to the importance of cost and pollution in the optimization. When this problem was solved as a single objective, values of 618 and 136 were obtained for cost and pollution, respectively, whereas in the sample, these values are 753.97 and 394.74. ple response from the Pareto front of the NSIDE method5 While the power generation capacity of all units is within the permissible range, it can be observed that at some times the network has largely purchased or sold electricity from the upstream network, which was due to the importance of cost and pollution in the optimization5 When this problem was solved as a single objective, values of 618 and 136 were obtained for cost and pollution, respectively, whereas in the sample, these values are 753.97 and 394.745   In order to illustrate the effectiveness of the NSIDE approach in solving the issue, this algorithm is implemented in the problem and its Pareto Front is compared with NSDE and NSPSO methods in Figure 7. It is quite clear that NSPSO is quite weak in solving this issue state and is not capable of solving this problem state, and NSDE has been defeated in all NSIDE cases. However, it has been shown in the preceding section that, without limiting the power exchange of different methods, it leads to similar solutions. The highlight of Pareto Front here is that it allows the user to choose the right answer. For example, a sample response from the Pareto Front with a total cost of $ 748.82 and a pollution of 582.81 is shown in Table 14. By comparing the amount of gas emission and cost mentioned with the single-purpose mode ($661 and 449 kg, respectively); that is, although in the multi-objective mode, each of the answers alone has a higher value than the single-objective mode, but for the intermediate mode, the answer is good, because both the cost and the contaminant are considered.

Simulation in the Mode of Maximum Generation of Renewable Resources
This section is intended to illustrate the importance of maximization of using RE sources5 In this regard, the experiment was conducted in two parts: single objective and

Simulation in the Mode of Maximum Generation of Renewable Resources
This section is intended to illustrate the importance of maximization of using RE sources. In this regard, the experiment was conducted in two parts: single objective and multi objective. In the single-objective part, each of the cost and pollution targets are considered separately, and in the multi-objective part, the NSIDE algorithm is used to optimize the problem.

Single Objective Optimization
In this section, the objectives of reducing costs and reducing the amount of pollution are considered separately, each in two modes, with and without the limitation of power exchange with the network. The results of the total costs and gas emission obtained from this experiment are presented in Table 15 and compared with the normal operation mode of the network. Obviously, increasing the costs and reducing the pollution can be very effective if using renewable energy sources, and Table 15 also illustrates this fact. The distribution of power in units is presented in Tables 16-19 for different scenarios with and without network exchange constraints. Figures 8 and 9 also show the convergence of graphs with and without the limitations of upstream grid exchange. Obviously, taking into account the limitation will cost more to the users. In this part of the simulation, as the amount of WT and PV pollution is zero, and the amount of cost is higher than the other units, these two units are not considered in this part.

Multi Objective Optimization
When multiple objectives are also considered in optimizing the utilization of the MG by maximizing the use of renewables, it has significant effects on the network. These effects are presented in the form of the Pareto front in Figure 10. The distance between these two fronts indicates the positive impact of maximization of using RE sources, which confirms the ability of reduction in the costs for maximizing use of RE sources. In this case, a comparison is also made between the proposed method and the other methods in two cases with and without grid limitations in Figures 11 and 12, which illustrates the superiority of the proposed method over other methods.

Sensitivity Analysis
In order to show the sensitivity, the algorithm has been implemented to solve the issue of optimal number of populations in the problem of optimal operation of MG with

Sensitivity Analysis
In order to show the sensitivity, the algorithm has been implemented to solve the issue of optimal number of populations in the problem of optimal operation of MG with

Sensitivity Analysis
In order to show the sensitivity, the algorithm has been implemented to solve the issue of optimal number of populations in the problem of optimal operation of MG with different initial populations in two single-objective and multi-objective modes. Populations of 25, 50, 75, and 100 are intended to solve the problem. Table 20 shows the impact of population change on cost and Table 21 displays the impact of initial population on pollution. Increasing the initial population has improved the answers. On the other hand, the running time of the program also increases. Figures 13 and 14 show the effects of the initial population on convergence respectively. It is clear that increasing population size improves the convergence of graphs. Figure 15 shows the Pareto Fronts with different populations; an increase in population makes the Fronts more dominant. of population change on cost and Table 21 displays the impact of initial population on pollution. Increasing the initial population has improved the answers. On the other hand, the running time of the program also increases. Figures 13 and 14 show the effects of the initial population on convergence respectively. It is clear that increasing population size improves the convergence of graphs. Figure 15 shows the Pareto Fronts with different populations; an increase in population makes the Fronts more dominant.  Figure 13. The influence of initial population on cost convergence. Figure 13. The influence of initial population on cost convergence.     Figure 16 compares normal loading and survivability of loads during different hours a day. As can be seen from Figure 16, a significant proportion of loads in the MG has been supplied, which in turn shows the effectiveness of proposed method to enhance power system resilience after long-time outage.

Future Works
In this paper, the power required by consumers is provided through a MG in which various and renewable distributed generation sources are used. Additionally, by managing production in these resources and using simultaneous optimization of targets using intelligent Pareto algorithm, it has caused the reduction of costs and the amount of environmental pollution. The use of important technology such as MGs, consideration of important goals such as environmental and economic issues, as well as the use of Pareto method and system optimization with multiple goals and improving the IDE algorithm are among the important points. The behavior of renewable resources due to their unpredictable performance, in terms of consumer load pattern and market competition, is another case that has been studied in this study. In addition, the system studied in this paper can be compared with other new evolutionary and improved algorithms.
With the development of communication systems and systems intelligence, the risks of cyber and malicious attacks are increasing, in this regard, systems should be designed

Future Works
In this paper, the power required by consumers is provided through a MG in which various and renewable distributed generation sources are used. Additionally, by managing production in these resources and using simultaneous optimization of targets using intelligent Pareto algorithm, it has caused the reduction of costs and the amount of environmental pollution. The use of important technology such as MGs, consideration of important goals such as environmental and economic issues, as well as the use of Pareto method and system optimization with multiple goals and improving the IDE algorithm are among the important points. The behavior of renewable resources due to their unpredictable performance, in terms of consumer load pattern and market competition, is another case that has been studied in this study. In addition, the system studied in this paper can be compared with other new evolutionary and improved algorithms.
With the development of communication systems and systems intelligence, the risks of cyber and malicious attacks are increasing, in this regard, systems should be designed and improved to be more resilient and resistant to these malicious attacks. So, in future work, intelligent power systems can be optimized to withstand cyber and physical attacks such as storms and earthquakes, and power management can be optimized and studied with applicable and efficient strategies.

Conclusions
Having an optimal energy management and power dispatch in smart MGs considering different factors such as cost or air pollutant effects can be highly effective in modern developed society. In this research, an improved multi-objective DE algorithm was proposed and used in a multi-operation energy management form in a smart MG system. One of the main innovations of the suggested technique is that the modifying meta-heuristic parameters within this model, unlike conventional methods, have been allocated to the self-adaptive system.
In order to demonstrate the high efficiency of the IMODE algorithm in solving the problem of optimal operation of the MG in both single-objective and multi-objective modes, this algorithm has been carried out on different problems. In order to show the impact of renewable energy on the cost, in the simulations, the operation has been solved with maximum participation of renewable energy resources. Additionally, gas emission and air pollution have been considered in simulation. In addition to this, presenting an efficient concept of energy exchange between the MG and the main power network can be highly beneficial for both objectives in the grid-connected mode. Additionally, the sensitivity analysis on the number of populations in the meta-heuristic algorithms has been performed, and the results of simulation were compared and analyzed. The feasibility and effectiveness of proposed method has been confirmed by applying on a study case system. Besides, using improved self-adaptive DE in optimization issues not only can provide the best respondent to the designing and planning requirements, but can satisfy emission and economic aims with an admissible precision.