Power Distribution Optimization Based on Demand Respond with Improved Multi-Objective Algorithm in Power System Planning

In this article, a novel dynamic economic load dispatch with emission based on a multi-objective model (MODEED) considering demand side management (DSM) is presented. Moreover, the investigation and evaluation of impacts of DSM for the next day are considered. In other words, the aim of economical load dispatch is the suitable and optimized planning for all power units considering different linear and non-linear constrains for power system and generators. In this model, different constrains such as losses of transformation network, impacts of valve-point, ramp-up and ramp-down, the balance of production and demand, the prohibited areas, and the limitations of production are considered as an optimization problem. The proposed model is solved by a novel modified multi-objective artificial bee colony algorithm (MOABC). In order to analyze the effects of DSM on the supply side, the proposed MODEED is evaluated on different scenarios with or without DSM. Indeed, the proposed MOABC algorithm tries to find an optimal solution for the existence function by assistance of crowding distance and Pareto theory. Crowding distance is a suitable criterion to estimate Pareto solutions. The proposed model is carried out on a six-unit test system, and the obtained numerical analyses are compared with the obtained results of other optimization methods. The obtained results of simulations that have been provided in the last section demonstrate the higher efficiency of the proposed optimization algorithm based on Pareto criterion. The main benefits of this algorithm are its fast convergence and searching based on circle movement. In addition, it is obvious from the obtained results that the proposed MODEED with DSM can present benefits for all consumers and generation companies.


Introduction
Electrical energy has particular downsides such as high cost, being time consuming, etc., and it is gradually is being replaced by other energies as a result. In recent years, many causes such as increasing electricity generation cost, the uncertainty growth of load consumption, the global increase of fuel cost, restrictions in the environment, determining the place of powerhouses, etc., have drawn important attention toward the logical growth of electricity consumption [1].
Therefore, the demand side management (DSM) of electricity consumption in recent years is considered one of the best and vital methods of maintaining common profits between the electricity industry and electrical energy consumers. DSM as a fundamental subject tries to change load consumption at peak-time to the form of long-term in dedicating resources regarding its numerous profits for both sides of the electricity industry and consumers [2,3]. Hereby, the power markets have paid significant attention to the programs of DSM due to their benefits in financial and load profile [4]. DSM not only helps the electricity industry and consumers, furthermore, it has effective benefits for generation programs of DSM due to their benefits in financial and load profile [4]. DSM not only helps the electricity industry and consumers, furthermore, it has effective benefits for generation companies. The implementation of DSM in electricity networks requires updated information and technologies. These technologies make a two-way connection between the electricity's suppliers and consumers. Therefore, the smart controlling and the measuring tools of central energy are used for their dynamic interaction. Dynamic bidding policies include time of use (ToU), division of the critical peak time, on-time bidding, the least off-peak time, and bidding for the next days of smart bidding tools in the process of DSM implementation [5][6][7][8]. According to recent papers, three types of DSM commonly are used: environmental, network, and market. DSM with the environmental option commonly concentrates the standards on environmental and social areas such as decrease of greenhouse gasses' emissions. DSM by using the network option aims to maintain the certainty of the network, and the aim of DSM is to trigger the financial frugality of public facilities and consumers [9]. The overall view of the suggested demand side management (DSM) is shown in Figure 1. In [10], a new integrated method of multi-objective dynamic economic and emission dispatch (MODEED) is presented by implementing demand side management (DSM) in order to consider the advantages of presented DSM on the group side. The presented method investigated the upcoming day on the basis of load changing and the DSM method.
In [11], the authors implemented genetic algorithm to solve the problem of DSM and dynamic economic dispatch (DED). The authors had deliberated the issues of DSM and DED in two interrelated levels during the optimization procedure. In fact, it had used load profile for each individual consumer. Thus, they used Genetic Algorithm (GA) in order to find out the enhanced problem to introduce electrical tools for consumers.
Lorenthiran et al. [12] implemented the DSM technique of load profile based on the next day in the environment of smart network with the assistance of a discovering algo- In [10], a new integrated method of multi-objective dynamic economic and emission dispatch (MODEED) is presented by implementing demand side management (DSM) in order to consider the advantages of presented DSM on the group side. The presented method investigated the upcoming day on the basis of load changing and the DSM method.
In [11], the authors implemented genetic algorithm to solve the problem of DSM and dynamic economic dispatch (DED). The authors had deliberated the issues of DSM and DED in two interrelated levels during the optimization procedure. In fact, it had used load profile for each individual consumer. Thus, they used Genetic Algorithm (GA) in order to find out the enhanced problem to introduce electrical tools for consumers.
Lorenthiran et al. [12] implemented the DSM technique of load profile based on the next day in the environment of smart network with the assistance of a discovering algorithm. The energy control strategy is developed for load demand management at peak times [13]. Another algorithm for energy management is suggested by Ma et al. [14] to access the bidding strategies and consumers' functional states. Dynamic emission and economical dispatch (DEED) are considered in the operation and control of an electricity network. The DEED in online production presents the predicted load demand period, which is done by minimizing the cost and emission simultaneously. Altinoz et al. [15] investigated a type of financial or emission load model based on the transmitting issue. The suggested model considered non-synchronous and reserved producers to varied demanding zones, which it contradicted with previous transmitting issues. Among them, only a few of the producers are considered for the demanding zones. Consequently, during the process, each producer has an effect over the other producers. In addition, the suggested model is enhanced by generated energy that could be purchased by other close locals. From the point that most of the dimensions, restrictions, and aims are presented within the provided model, it declines the frequency of aims. Alham et al. [16] recommended a novel DEED with wind power regarding the feature of intermittency and ambiguity. They used an energy storage system (ESS) and DSM in order to investigate the effect of charge, discharge, and the use of wind power.
In order to solve the EED, the GAMS software was used. The obtained results demonstrated the significance of applying ESS and DSM in declining the charge and discharge rate and raising the wind power consumption. Dehnavi et al. [17] evaluated TOU on the demand side by combining the problem of DED, which considers the supplier side. Zaman et al. [18] combined an algorithm of DE and GA for solving the issues of DEED by considering different generator units such as thermal, solar-thermal, and wind. The problem of DEED based on the developed PSO and wind's uncertainty is presented by Aghaei et al. [19].
In [20], new combined framework-based electric vehicles and wind farms are suggested for solving the dynamic economic emission dispatch problem. In [21], the authors presented the combination of renewable energy sources to a multi-objective DEED problem in a day through minimizing the energy cost and emission as well. In [22], the DEED problem is solved as a non-linear and complex objective based on a combined Flower Pollination Algorithm (FPA) with Sequential Quadratic Programming (SQP) algorithm. The authors in [23] provide MODEED for renewable anergy sources i.e., wind-solar-hydro power under tradable green certificates.
To cover the gap of economic load dispatch with demand side management, this article suggests a method of MODEED and DSM program based on next day for residential loads. In this suggested model, we had considered the benefit of consumers by presenting different time delays for all manageable demand with related lifestyle. The suggested model significantly focused on advantages of DSM programs in order to accommodate the supplier side with different residential load levels. In addition, in this paper, an enhanced multi-objective artificial bee colony (MOABC) based on entropy, arrangement, and fuzzy selection is proposed to minimize the MODEED simultaneously with DSM.
The rest of this paper is organized as follows. Section 2 explains the formulation of the DEED problem and DSM. The suggested integrated method is presented in Section 3. Section 4 presented the suggested improved multi-objective artificial bee colony algorithm. In Section 5, some different test systems are presented. The numerical analysis and their comparisons are presented in Section 6. At the last section, the conclusion is presented.

Proposed Multi-Objective Problem
This multi-objective problem is the vital optimization challenge in the electricity network operation and control for enhancing social and economic dimensions. In the DEED problem, the total fuel cost and emission cost should be minimized simultaneously [24]: (2) F 1 and F 2 are the total costs of fuel and the emission of the gth generator unit at time τ. C g,t and E g,t are the price and pollution functions through time t and P g,t is the output of the gth generator at unit time t. The fuel cost function is modeled by the effect of the valve point, and the total rate of gases such as NO x and SO x are modeled as the sum of quadratic expressions: C g,t P g,t = a g P 2 g,t + b g P g,t + c g + e g sin e g P min g − P g,t E g,t P g,t = γ g P 2 g,t + β g P g,t + α g + ξ g exp λ g P g,t where c g , b g , and a g are the of fuel cost coefficients of the gth generator unit and e g and d g are the valve point effect coefficients. γ g , β g , α g , and λ g are emission coefficients of the gth generator unit.

Constrains
The problem of DEED optimization has different equality and inequality constrains as follows: P L,t and P D,t are the transmission power loss and total load demand power in the tth time, respectively. P L,t could be calculated by: where B(g,j) denotes the coefficient of transmission power loss. P g,t and P j,t are the real electricity outputs of the gth and jth generator unit at the tth time, respectively.

(ii) Generation boundary
This inequality constrain can be modeled by: where P max g and P min g describe the lower and higher values for power output in the gth generator unit.

(iii) Ramp-up and Ramp-down
This inequality constrain can be modeled by: Therefore, the above equations can be rewritten as: max P min g , P g,t−1 − DR g ≤ P g,t ≤ min P max g , P g,t−1 + UR g where UR g and DR g are the upper and lower ramp bound rates of the gth of generation unit. They show that increasing or decreasing the output powers for thermal production units is a time-consuming action.

Problem of Demand Side Management
The strategies of DSM mostly are needed for public facilities in order to grow their economics and develop a demand shape. Typically, companies of public industry have interesting DSM methods for encouraging the higher incorporation of consumers. Therefore, the public industry companies could easily succeed in reaching their aims. According to recent papers that address the DSM problem, the clipping of peak, filling the path, load shifting, strategic consideration, growth of strategic demand, and supple load curve are considered as the six main approaches [7]. The first three approaches are basic levels and the others are innovative methods that monitor the total load curve by either limiting or extending.
The load shifting method is the sum of two restricting peak and filling path methods, and it is considered the best method among the all demand management approaches. The demand-shifting model could be applied by controllable demands in the user side. In this mode, the controllable demands shift from peak time to off-peak, excluding any change in consuming power. The six methods of DSM are illustrated in Figure 2.
where URg and DRg are the upper and lower ramp bound rates of the gth of generation unit. They show that increasing or decreasing the output powers for thermal production units is a time-consuming action.

Problem of Demand Side Management
The strategies of DSM mostly are needed for public facilities in order to grow their economics and develop a demand shape. Typically, companies of public industry have interesting DSM methods for encouraging the higher incorporation of consumers. Therefore, the public industry companies could easily succeed in reaching their aims. According to recent papers that address the DSM problem, the clipping of peak, filling the path, load shifting, strategic consideration, growth of strategic demand, and supple load curve are considered as the six main approaches [7]. The first three approaches are basic levels and the others are innovative methods that monitor the total load curve by either limiting or extending.
The load shifting method is the sum of two restricting peak and filling path methods, and it is considered the best method among the all demand management approaches. The demand-shifting model could be applied by controllable demands in the user side. In this mode, the controllable demands shift from peak time to off-peak, excluding any change in consuming power. The six methods of DSM are illustrated in Figure 2.  The DSM technique of load shifting based on the next day is used in this paper in order to investigate its impacts on the production side. This paper considers minimizing the bill of public energy as a basic aim for implementing DSM. Therefore, the prediction curve of load demand is estimated based on previous data. Consequently, the price signal is dedicated for each unique load, and the public company creates an aim load shape that is contrary with separate values. During the previous day, the DSM central controller receives an aim load shape as an input and computes the controlling measurements of ideal demand consumptions. Regarding the obtained results, the controlling measurements were Energies 2021, 14, 2961 6 of 18 conducted in ON time periods. During ON time operations, when the consumer sends the request connection of a device via bottom ON over the domestic devices, DSM controls or permits the attachment or provides another attachment time. The total ON time procedure is conducted efficiently by using mutual information and connection technologies. The total sell price of the electricity market of each individual break is dedicated by the ToU tariff approach, which is one of the methods of pricing tools. In the ToU tariff model, the pricing of the critical peak and lower price for off-peak periods are evaluated. According to theses pricing tools, DSM members want to place their plans in off-peak periods. The DSM can be controlled by electricity consumption patterns. In this paper, the aim of DSM with the suggested algorithm is discovering an optimized solution.

Utility Energy Modeling of DSM
The critical purpose of public industry power bill function is decreasing the space of the predicted load demand curves and aimed load demand trend through the controllable loads. The proposed function is formulated as follows: where ξ(t) and ψ(t) are the quantities of aimed load demands and the load after timing DSM at period t. τ denotes the numbers of existence total time in load demand profile in a day. The suggested flowchart is demonstrated in Figure 3.
the bill of public energy as a basic aim for implementing DSM. Therefore, the prediction curve of load demand is estimated based on previous data. Consequently, the price signal is dedicated for each unique load, and the public company creates an aim load shape that is contrary with separate values. During the previous day, the DSM central controller receives an aim load shape as an input and computes the controlling measurements of ideal demand consumptions. Regarding the obtained results, the controlling measurements were conducted in ON time periods. During ON time operations, when the consumer sends the request connection of a device via bottom ON over the domestic devices, DSM controls or permits the attachment or provides another attachment time. The total ON time procedure is conducted efficiently by using mutual information and connection technologies. The total sell price of the electricity market of each individual break is dedicated by the ToU tariff approach, which is one of the methods of pricing tools. In the ToU tariff model, the pricing of the critical peak and lower price for off-peak periods are evaluated. According to theses pricing tools, DSM members want to place their plans in off-peak periods. The DSM can be controlled by electricity consumption patterns. In this paper, the aim of DSM with the suggested algorithm is discovering an optimized solution.

Utility Energy Modeling of DSM
The critical purpose of public industry power bill function is decreasing the space of the predicted load demand curves and aimed load demand trend through the controllable loads. The proposed function is formulated as follows: where ( )   According to above equation, the ψ(t) can be computed as follows: where FL(t) is the quantity of the predicted load at period t, CL(t) and DL(t) are the quantities of attached load demands and shut-off during demand shifting at period t. It is calculated as follows: X lit is the amount of l manageable domestic devices that is changed from i to t time, N is the total type of controllable domestic devices, P 1l and P (t−k+1)l are the electricity load consumption due to the change device of l type and time steps (t -k + 1). Here, the (t − k + 1) is the time step of electricity consumption with device type 1 in t time, which is changed from j to k.
where X ltq is the number of controllable domestic devices from type 1 that had changed from t time to q period and in addition, the delayed maximized steps of m have the ability to manage the controllable domestic devices. For the load-shifting process, the following constrains should be considered for load management: where the Ctrd (i) is the total number of controllable domestic devices from type 1 in time t.
The other overall conditions are provided as follows: where D depicts the whole time period. Equation (18) presents that the DSM approach only has the delayed feature and it is not of advanced type. Equation (19) is the indicator of maximized m delayed for all domestic devices. In this paper, we had provided solutions to the multi-objective optimization model with the non-dominated sort and fuzzy tools. At this point, a linear membership function u i,k (kth objective function, ith solution) is proposed for each of the objective functions F i as follows: where O and S denote the number of objective functions and non-dominated solutions, respectively. F min k and F max k are the minimum and maximum of Fk(X) of the kth objective function, while u is the quantity of the kth objective function. The maximized quantity of u makes a decision about the ultimate solution of the multi-objective problem.

Improved Multi-Objective Artificial Bee Colony Algorithm
The artificial bee colony (ABC) has been proposed by Mr. Karaboka [25]. To improve the classic version of ABC, a chaotic version of this model is considered in this paper. The ergodicity and randomness specifications of chaos can understand a local deep search. An optimal local optimized ability is provided by looking for the space near superior individuals. The concept of a chaotic system can strengthen the local search ability. The restrictions of this region are determined by system parameters. Based on other features of chaos systems, they are looking at a particular area of primary conditions. In this paper, the mentioned improvement is presented by the following steps. First, the riotous coefficient appropriate with the determined region and the number of variables are generated. Then, the chaotic variables of the next iteration are calculated and then converted into decision variables. Finally, if the new results are better than the initial one/reach max iteration, they replace it as the final answer; otherwise, go to step 2. This paper had used the features of a riotous system to find ideal weights. The logistic function is a non-linear function and has one control parameter: where x(k) shows the system condition and A is a constant value. This function intensively depends on the parameter A. By changing A, this function shows behaviors such as (from periodic). By putting A = 4 in the equation, the logistic function shows chaos behavior. From the point of formulation, the riotous coefficient appropriate with a determined region and the number of variables (Ng) are: So, we achieve: where X 0 cls is the obtained primary place for the chaos variable. X j,min and X j,max are the lower and upper values of the chaos variable. N chaos is the number of chaos variables, respectively. On the other hand, the study and investigations of artificial bee colony show that while the particles are leading by g best , they will have better total searching. In the artificial bee colony algorithm, the investigation is done on the basis of worker bees, and it is not utilized from the obtained best solution in each stage. Based on this principle, the standard equation of bee colony algorithm (25) is changed to Equation (26): where φ ij and λ ij are the random numbers between (−1,1). In comparison to the other optimization methods, the proposed methods have the following significant benefits: The proposed algorithm can be followed in the below steps: 1.
Generate the initial solutions X ij 2.
Compute the initial solutions in objective function 3.
Supply new solutions based on Equation (26)  5.
Select the best resource or best answer between X ij and V ij 6.
Calculate the possibility rate for X ij solutions based on: Indeed, to achieve appropriate solutions, one gets: The solutions of P i are between (−1, 1).

7.
Generate new solutions V i based on ABC operator from solutions Xi and determine their possibility rate P i . 8.
Select the best solutions between the answers of X ij and V ij . 9.
Determine the distorted resources and substitute the random resources with created random resources by scout bee Xi by: 10. Save the best solution (the most qualified resource) that is obtained in this stage. 11. Cycle = cycle + 1. 12. Repeat the previous stages until arriving at the stop criterion.
For modeling the suggested method, let X be an n-dimension key of the search space through the decision-making vector X = x 1 , x 2 , x 3 , . . . , x n and X* be the optimized result of the fitness f (x) = [ f 1 (x), f 2 (x), . . . f k (x)] for which the non-linear condition of g i (x) ≥ 0, i = 1, 2, . . . , m and the P linear provision of h i (x) = 0, i = 1, 2, . . . , p are set. In the procedure of the suggested approach, all functions based on Pareto criteria are connected together. It can be noted that the problem of optimization could be converted to a maximized problem of optimization and vice versa: Essentially, the above formula could be answered in two methods. One is the combination of fitness as one objective function. Therefore, the objective functions with appropriate and precise weight coefficients should be integrated together. The improved method of this article or the 2nd approach are based on sets of optimized Pareto. Assume that the C is a probable region answer (sets of results of x ∈ X that cover all linear and non-linear provisions) and I = {1, 2, . . . , k}. The f x * ∈ C will be an optimized Pareto if ∀ x∈C is the ∀ i∈I f i (x * ) = f i (x) or ∃ i∈I f i (x * ) < f i (x). U = u 1 , u 2 , u 3 , . . . , u k and V = v 1 , v 2 , v 3 , . . . , v k are the vectors in the considered functions and U will dominate for V if and only if the ∀ i∈I , u i ≤ v i and ∃i ∈ I when U i < V i . The optimized Pareto (P*) is a set of possible solutions if P * = {x ∈ C| ∃x ∈ C, f( x ) ≤ f(x)} . In this work, the PF = {f(x) |x ∈ P * } as a Pareto front. The dominant count in a set of answers is equivalent to the numbers of solutions in the whole population that are dominant. The initial non-dominated front is equivalent to the sets of solutions whose domination count equals 0. The 2nd non-dominated front is the sets of results whose dominant count equals zero and whose 1st non-dominated front is uncountable. Meanwhile, the total ith of the non-dominant front ψ i is a set of solutions whose dominant count equals zero when the solutions are uncountable in front ψ 1 , ψ 2 , . . . , ψ i−1 . Now, all the obtained solutions throughout the 1st front include the non-dominant front value 1. Naturally, the solutions throughout the second front consist of non-dominant value 2. The mutuality approach is evaluated in guarantee the discovery of optimized solutions. In the method, the quantity of evaluated suitability for every aim function in level 1 is shared. Then, this quantity of suitability that is the little shared quantity in past functions is signed as a next Pareto front. Similarly, the method of sharing is done for the second level. Overall, this strategy will be continued until the quantity of suitability is joint in all levels. The possible equation for the shared approach (Share (d ij )) is defined as follows: where p 1 is the number of decisions made, x s is the quantity of the sth decision-making variable, and i and j are counters. µ share is the highest acceptable distance between two variables, and its signing as the member niche indicates that this counter could be computed for the total (N) population as follows: Additionally, to enhance the convergence procedure between the non-dominated sets, the crowding distance is evaluated for evaluating the set of results [26]. This method is based on the computation of average distance, and an appropriate estimation of the density rate of better results between populations is achieved. In other words, this technique leads to the concentration of Pareto solutions. The entropy distribution for the ith objective function could be expressed as follows: where du ij and dl ij are the current solutions with determined upper and limited restrictions of the ith objective function. As expressed before, the optimized Pareto scale is the indicator that vector x ∈ C is not present, which could modify the optimized Pareto vector. The fuzzy logic is implemented as a strategy in order to evaluate the appropriate cluster of their non-linear ones. The fuzzy decision-making function, which is presented by the membership function, should be able to substitute the accurate quantity of parameters.

Simulation Test Systems
In this article, the six thermal units considered as the case study are considered in different scenarios. The fuel cost coefficient of a generator, emission coefficient, the bound of the ramp rate, the predicted load demand profile for 24, and other corresponding system data can be found in [27][28][29], and for the reader's convenience, they are presented in Table 1 and Figure 4. As presented in this figure, the first unit production is higher than the other units. The matrix of transmission loss constant of the six studied system units is shown in Equation (35). In total, the demand of each network is a combination of domestic, commercial, and industrial demands. The residential area has several controllable domestic devices such as a dish washing machine, washing machine, kettle, oven, and other domestic devices that could be simply managed during the demand shifting approach in comparison with other industrial and commercials demands. Hence, only the demand shift technique of DSM based on the next day is implemented for residential loads. During the implementation process of DSM, four cooperating levels with 5%, 10%, 15%, and 20% of the whole habitable demand are evaluated. The electrical prices of the market are dedicated for total public facilities in each hour with the ToU tariff. In the mentioned test, the utilization of total controllable domestic devices is considered 2 kW, but their pattern of electricity consumption is different. In the mentioned model, the maximized late steps of every controllable device based on their consumption lifestyle are dedicated. For example, two or three maximized delayed steps for necessary devices such as a kettle, oven, and hair dryer are evaluated in this work. In same approach, 12 maximum of delayed steps for less necessary devices such as a cloth dryer, iron, and vacuum cleaner are considered.  Figure 4. The load demand profile of 24 h from six test system units.

The Results of Simulation and Discussion
In this article, the effects of DSM are studied over the public and production in try companies. The issue optimization of DSM and DEED is calculated in MAT software. In order to find the best answer, 30 different tests are conducted.

The Results of Simulation and Discussion
In this article, the effects of DSM are studied over the public and production industry companies. The issue optimization of DSM and DEED is calculated in MATLAB software. In order to find the best answer, 30 different tests are conducted.

The Effects of DSM on the Public Facility Side
The main aim of this paper is analyzing the effects of DSM evaluation on the production side. Thus, the DSM is evaluated by minimization of the public electricity bill as the basic goal. Figure 5 presents the predicted load shape and the load curves after implementation of DSM with different levels of residential cooperation. It is obvious that an improved demand curve for higher cooperator levels is achieved. The public industry company has achieved interesting advantages, and this is because of the DSM evaluation that is presented in Figure 6a,b. By the obtained results of the DSM problem, the daily frugality of energy bills of the public industry company for participation levels are 5%, 10%, 15%, and 20% from those whole inhabitant demand are $363, $731, $1950, and $1457, respectively. Similarly, the peak load demand also for different levels of participation is 12.32, 24.32, 36.12, and 49.01 MW, respectively.

The Results of Simulation and Discussion
In this article, the effects of DSM are studied over the public and production industry companies. The issue optimization of DSM and DEED is calculated in MATLAB software. In order to find the best answer, 30 different tests are conducted.

The Effects of DSM on the Public Facility Side
The main aim of this paper is analyzing the effects of DSM evaluation on the production side. Thus, the DSM is evaluated by minimization of the public electricity bill as the basic goal. Figure 5 presents the predicted load shape and the load curves after implementation of DSM with different levels of residential cooperation. It is obvious that an improved demand curve for higher cooperator levels is achieved. The public industry company has achieved interesting advantages, and this is because of the DSM evaluation that is presented in Figure 6a,b. By the obtained results of the DSM problem, the daily frugality of energy bills of the public industry company for participation levels are 5%, 10%, 15%, and 20% from those whole inhabitant demand are $363, $731, $1950, and $1457, respectively. Similarly, the peak load demand also for different levels of participation is 12.32, 24.32, 36.12, and 49.01 MW, respectively.  In order to analyze the DSM effects in production, three test vases by various participation levels of DSM are considered. In the first and second test cases, the problems of dynamic economical and emission dispatch are individually minimized by the DSM approach. Case 3 elaborates individually the suggested approach MODEED by the DSM model. In order to analyze the DSM effects in production, three test vases by various participation levels of DSM are considered. In the first and second test cases, the problems of dy-

Case 1: DSM Bundled with Dynamic Economic Dispatch
The single dynamic economic dispatch is considered here as the first test case by considering the DSM effects over a fuel cell. For qualitative study, four different levels of participation are presented in DSM. In this approach, dynamic economic dispatch and DSM have been minimized simultaneity based on the suggested approach. The obtained numerical results are presented in Table 2 by including/excluding the DSM. Figure 7 demonstrates the convergence issues of DEED optimization including/excluding DSM. Results demonstrate the useful effects of DSM on electricity production through decreasing their fuel costs. In addition, it shows that the quantity of fuel production cost is minimized in different participation levels of DSM. For instance, DSM with the participation level of 20% has about a 7.5% drop in daily fuel cost in compared to the condition without DSM. Production companies can save a daily fuel cost equal to $91.5, and this frugality is the result of DSM implementation.  In order to analyze the DSM effects in production, three test vases by various participation levels of DSM are considered. In the first and second test cases, the problems of dynamic economical and emission dispatch are individually minimized by the DSM approach. Case 3 elaborates individually the suggested approach MODEED by the DSM model.

Case 1: DSM Bundled with Dynamic Economic Dispatch
The single dynamic economic dispatch is considered here as the first test case by considering the DSM effects over a fuel cell. For qualitative study, four different levels of participation are presented in DSM. In this approach, dynamic economic dispatch and DSM have been minimized simultaneity based on the suggested approach. The obtained numerical results are presented in Table 2 by including/excluding the DSM. Figure 7 demonstrates the convergence issues of DEED optimization including/excluding DSM. Results demonstrate the useful effects of DSM on electricity production through decreasing their fuel costs. In addition, it shows that the quantity of fuel production cost is minimized in different participation levels of DSM. For instance, DSM with the participation level of 20% has about a 7.5% drop in daily fuel cost in compared to the condition without DSM. Production companies can save a daily fuel cost equal to $91.5, and this frugality is the result of DSM implementation.

Case 2: DSM Bundled with Dynamic Emission Distribution
The single dynamic emission dispatch objective is considered in these test cases including the DSM. Obtained numerical analyses are presented in Table 3 by considering the dynamic emission distribution and including/excluding the DSM. Additionally, in this table, the proposed algorithm is compared with two other optimization models i.e., multi-objective genetic algorithm (NSGA) and multi-objective particle swarm optimization (MOPSO). Obtained results prove the validity and superiority of the proposed algorithm in minimum cost and emission as well. Figure 8

Case 2: DSM Bundled with Dynamic Emission Distribution
The single dynamic emission dispatch objective is considered in these test cases including the DSM. Obtained numerical analyses are presented in Table 3 by considering the dynamic emission distribution and including/excluding the DSM. Additionally, in this table, the proposed algorithm is compared with two other optimization models i.e., multi-objective genetic algorithm (NSGA) and multi-objective particle swarm optimization (MOPSO). Obtained results prove the validity and superiority of the proposed algorithm in minimum cost and emission as well. Figure 8

Case 3: DSM Bundled with Economical and Emission Dispatch
The DEED model by considering the DSM is considered in this test case through the MOABS algorithm. The obtained numerical analysis is presented in Figure 9 based on MODEED by considering various levels of DSM participation. In addition, Figure 10 shows the set of different Pareto values with different levels of DSM participation. The finest adaptive quantity of fuel cost and distribution using the fuzzy membership ap-

Case 3: DSM Bundled with Economical and Emission Dispatch
The DEED model by considering the DSM is considered in this test case through the MOABS algorithm. The obtained numerical analysis is presented in Figure 9 based on MODEED by considering various levels of DSM participation. In addition, Figure 10 shows the set of different Pareto values with different levels of DSM participation. The finest adaptive quantity of fuel cost and distribution using the fuzzy membership approach is found.  It is seen from the results that the proposed model has the ability to present the best solution compared to the single objective model for the similar participation level in DSM. For instance, the production companies could have saved a partial cost of daily fuel of 60.91%, 180%, and 256% for participation levels of 5%, 10%, 15%, and 20% of inhabitant demand, respectively. At this time, the production companies can decrease their sets of distribution levels. The decreases of distribution levels for different levels of participation are 0.032, 0.056, 0.074, and 0.121, respectively. The values of production units for timing programs of the proposed model in the participation level of DSM are achieved at 20% with or without DSM and are presented in Figure 11a,b.  It is seen from the results that the proposed model has the ability to present the best solution compared to the single objective model for the similar participation level in DSM. For instance, the production companies could have saved a partial cost of daily fuel of 60.91%, 180%, and 256% for participation levels of 5%, 10%, 15%, and 20% of inhabitant demand, respectively. At this time, the production companies can decrease their sets of distribution levels. The decreases of distribution levels for different levels of participation are 0.032, 0.056, 0.074, and 0.121, respectively. The values of production units for timing programs of the proposed model in the participation level of DSM are achieved at 20% with or without DSM and are presented in Figure 11a,b. It is seen from the results that the proposed model has the ability to present the best solution compared to the single objective model for the similar participation level in DSM. For instance, the production companies could have saved a partial cost of daily fuel of 60.91%, 180%, and 256% for participation levels of 5%, 10%, 15%, and 20% of inhabitant demand, respectively. At this time, the production companies can decrease their sets of distribution levels. The decreases of distribution levels for different levels of participation are 0.032, 0.056, 0.074, and 0.121, respectively. The values of production units for timing

Test Function for Multi-Objective Algorithm
This section presents a new test system to evaluate the proposed multi-objective ABC algorithm and compare its results with other well-known multi-objective optimization algorithms. Therefore, the ZDT2 test function issue is selected [29]: If n = 30, the problem of non-dominated Pareto front will be 219 in the regional Pareto front. To compare the proposed algorithm, consider two powerful MOPSO and NSGA-I algorithms with the same initial conditions such as the number of populations at 100 and a maximum iteration of 1000. The best Pareto distribution for the proposed MOABC is arrived after 350, while this number for NSGA-I is 650, and for MOPSO, it is about 1000 iterations. The final Pareto solutions are shown in Figure 12. Figure 12 shows that the superiority of scattering solution points (Pareto front) is achieved through the proposed algorithm compared to the NSGA and MOPSO models.

Test Function for Multi-Objective Algorithm
This section presents a new test system to evaluate the proposed multi-objective ABC algorithm and compare its results with other well-known multi-objective optimization algorithms. Therefore, the ZDT2 test function issue is selected [29] If n = 30, the problem of non-dominated Pareto front will be 219 in the regional Pareto front. To compare the proposed algorithm, consider two powerful MOPSO and NSGA-I algorithms with the same initial conditions such as the number of populations at 100 and a maximum iteration of 1000. The best Pareto distribution for the proposed MOABC is arrived after 350, while this number for NSGA-I is 650, and for MOPSO, it is about 1000 iterations. The final Pareto solutions are shown in Figure 12. Figure 12 shows that the superiority of scattering solution points (Pareto front) is achieved through the proposed algorithm compared to the NSGA and MOPSO models.

Conclusions
This paper combined the MODEED with the DSM model in order to use the advantages of DSM in the public and production industry side. During the process of DSM, the load shift technique based on the next day for controlling residential loads and minimizing the public facilities energy bill energy is considered. Considering the regular manner of consumers, all essential household devices are included. The objective functions of emission and power dispatch were minimized for a system with six thermal units. Four various levels of inhabitant participation for investigating the impacts of DSM in the production side were considered. The multi-objective method is utilized for minimizing the MODEED based on a non-dominated arrangement. The simulation results demonstrated that the proposed model has the ability to present benefits for public companies and the production side. Furthermore, it is observed that the benefits of a 20% participation level are high.

Conclusions
This paper combined the MODEED with the DSM model in order to use the advantages of DSM in the public and production industry side. During the process of DSM, the load shift technique based on the next day for controlling residential loads and minimizing the public facilities energy bill energy is considered. Considering the regular manner of consumers, all essential household devices are included. The objective functions of emission and power dispatch were minimized for a system with six thermal units. Four various levels of inhabitant participation for investigating the impacts of DSM in the production side were considered. The multi-objective method is utilized for minimizing the MODEED based on a non-dominated arrangement. The simulation results demonstrated that the proposed model has the ability to present benefits for public companies and the production side. Furthermore, it is observed that the benefits of a 20% participation level are high.