Combined Economic Emission Dispatch with Environment-Based Demand Response Using WU-ABC Algorithm

Owing to the growing interest in environmental problems worldwide, it is essential to schedule power generation considering the effects of pollutants. To address this, we propose an optimal approach that solves the combined economic emission dispatch (CEED) with maximum emission constraints by considering demand response (DR) program. The CEED consists of the sum of operation costs for each generator and the pollutant emissions. An environment-based demand response (EBDR) program is used to implement pollutant emission reduction and facilitate economic improvement. Through the weighting update artificial bee colony (WU-ABC) algorithm, the penalty factor that determines the weighting of the two objective functions is adjusted, and an optimal operation solution for a microgrid (MG) is then determined to resolve the CEED problem. The effectiveness and applicability of the proposed approach are demonstrated via comparative analyses at a modified grid-connected MG test system. The results confirm that the proposed approach not only satisfies emission constraints but also ensures an economically superior performance compared to other approaches. These results present a useful solution for microgrid operators considered environment issues.


Introduction
Economic dispatch (ED) is one of the most important issues pertaining to the operation and control of power systems [1,2]. It involves finding an optimal scheduling for all committed generators to minimize fuel costs, while satisfying various constraints such as load demand balance and generation capacity constraints. However, owing to increasing global concerns regarding the environmental issues caused by the combustion of fossil fuels, ED has become a significant concern from an economical perspective and for dealing with pollutant emissions from fossil fuel combustion [3]. Some countries have set maximum emission limits and impose fines when these limits are exceeded [4]. Thus, emission constraints should also be considered during operational scheduling. Economic emission dispatch (EED) has been proposed for scheduling to account for the emission of harmful gases such as carbon dioxide (CO 2 ), sulfur oxides (SO x ), and nitrogen oxides (NO x ), which pollute the air and exacerbate global warming [5]. The objective function of the EED usually adds emission criteria to the fuel cost of the thermal unit. To integrate the emission component with ED, Dhillon et al. [6] and Kulkarni et al. [7] presented a method whereby a penalty factor is multiplied to the emission term in the objective function. This technique allows both components to become commensurately involved in the optimization.
Demand response (DR) is expected to help solve environmental problems, as it is environmentally friendly and offers several benefits to the entire system [8,9]. This is because DR resources are one of Here, the penalty factor is used to convert the emission into an equivalent cost value. Various constraints are considered in solving the problem. In particular, the maximum emission constraint is additionally considered. To solve this problem optimally, we propose an environment-based demand response (EBDR) program that determines volume by optimizing economics and reducing the degree of pollutant emissions from the system. A weighting update artificial bee colony algorithm (WU-ABC) is applied to solve the CEED problem under a maximum emission constraint. Here, the penalty factor is updated using the weighting update method until all the constraints are satisfied. By using WU-ABC algorithm, the proposed approach becomes capable of finding the optimal CEED solution of the MG.
The primary contributions of this paper can be summarized as follows: • The CEED problem, with a maximum emission constraint on the MG, involves the minimization of operation costs and pollutant emissions. A penalty factor is used by the MGO to control the importance of both objective functions with respect to pollutant emissions. This approach is sufficiently effective to satisfy the maximum emission constraints on MG operation. • EBDR is applied and both economic and environmental factors are considered. Through the concept of elasticity, the MGO can command demand shifts and reduction, and participants in the DR program are provided appropriate incentives. Through this EBDR program, economics can be improved and environmental hazards can be reduced.

•
To solve the CEED problem, we propose the use of the WU-ABC algorithm. This algorithm can efficiently update the penalty factor according to the maximum emission constraint by using the weighting update method in combination with the ABC algorithm, which has been widely used for the optimal operation of power systems of late.
The proposed approach provides an optimal power scheduling using EBDR for MG systems in CEED problems with maximum emission constraints. In this regard, MGOs can be provided with more reasonable and flexible solutions for optimal operation with environmental constraints.
The remainder of this paper is organized as follows. Section 2 introduces the grid-connected MG system model. Section 3 presents the formulation of the problem for the CEED with a maximum emission constraint. Section 4 outlines the proposed EBDR strategy. Section 5 presents the solution of the WU-ABC algorithm and summarizes the overall process. Section 6 presents the simulation results, and Section 7 presents the conclusions of this study. Figure 1 shows the structure of the proposed grid-connected MG that includes DG, photovoltaic (PV), wind turbine (WT) units and controllable demands; it is connected to the main grid. MGOs can sell surplus power or buy insufficient power through interactions with the main grid, and they can request demand reduction from the customers using the DR program [28,29]. emissions. Here, the penalty factor is used to convert the emission into an equivalent cost value. Various constraints are considered in solving the problem. In particular, the maximum emission constraint is additionally considered. To solve this problem optimally, we propose an environment-based demand response (EBDR) program that determines volume by optimizing economics and reducing the degree of pollutant emissions from the system. A weighting update artificial bee colony algorithm (WU-ABC) is applied to solve the CEED problem under a maximum emission constraint. Here, the penalty factor is updated using the weighting update method until all the constraints are satisfied. By using WU-ABC algorithm, the proposed approach becomes capable of finding the optimal CEED solution of the MG. The primary contributions of this paper can be summarized as follows:

Grid-Connected MG System
• The CEED problem, with a maximum emission constraint on the MG, involves the minimization of operation costs and pollutant emissions. A penalty factor is used by the MGO to control the importance of both objective functions with respect to pollutant emissions. This approach is sufficiently effective to satisfy the maximum emission constraints on MG operation. • EBDR is applied and both economic and environmental factors are considered. Through the concept of elasticity, the MGO can command demand shifts and reduction, and participants in the DR program are provided appropriate incentives. Through this EBDR program, economics can be improved and environmental hazards can be reduced.

•
To solve the CEED problem, we propose the use of the WU-ABC algorithm. This algorithm can efficiently update the penalty factor according to the maximum emission constraint by using the weighting update method in combination with the ABC algorithm, which has been widely used for the optimal operation of power systems of late.
The proposed approach provides an optimal power scheduling using EBDR for MG systems in CEED problems with maximum emission constraints. In this regard, MGOs can be provided with more reasonable and flexible solutions for optimal operation with environmental constraints.
The remainder of this paper is organized as follows. Section 2 introduces the grid-connected MG system model. Section 3 presents the formulation of the problem for the CEED with a maximum emission constraint. Section 4 outlines the proposed EBDR strategy. Section 5 presents the solution of the WU-ABC algorithm and summarizes the overall process. Section 6 presents the simulation results, and Section 7 presents the conclusions of this study. Figure 1 shows the structure of the proposed grid-connected MG that includes DG, photovoltaic (PV), wind turbine (WT) units and controllable demands; it is connected to the main grid. MGOs can sell surplus power or buy insufficient power through interactions with the main grid, and they can request demand reduction from the customers using the DR program [28,29].

DG Model
The power generation costs of a diesel generator are generally expressed as the sum of a quadratic and sinusoidal function.
Here, the sinusoidal function represents a valve point effect to practically account for the generation cost function. Figure 2 illustrates the valve point effect for a conventional generator. It represents a sharp increase in losses due to the wire drawing effect caused by the opening of each steam admission valve [30].

DG Model
The power generation costs of a diesel generator are generally expressed as the sum of a quadratic and sinusoidal function.
Here, the sinusoidal function represents a valve point effect to practically account for the generation cost function. Figure 2 illustrates the valve point effect for a conventional generator. It represents a sharp increase in losses due to the wire drawing effect caused by the opening of each steam admission valve [30].

PV Model
The output power of a PV is dependent on solar irradiation. Forecasted solar irradiation is commonly determined using the beta probability distribution function (PDF) expressed as follows [31]. ( 1) ( 1) The shape parameters a and b are determined according to the mean (μs) and standard deviation (σs) of solar irradiation data.
Mathematically, the generated power from the PV array is represented as The generated solar power is represented by the product of PV panel efficiency, size of the PV array, and solar irradiation. β denotes temperature coefficients of the maximum power of the PV array.

WT Model
The Weibull PDF has been regularly used to model wind speed and can be expressed as

PV Model
The output power of a PV is dependent on solar irradiation. Forecasted solar irradiation is commonly determined using the beta probability distribution function (PDF) expressed as follows [31].
The shape parameters a and b are determined according to the mean (µ s ) and standard deviation (σ s ) of solar irradiation data.
Mathematically, the generated power from the PV array is represented as The generated solar power is represented by the product of PV panel efficiency, size of the PV array, and solar irradiation. β denotes temperature coefficients of the maximum power of the PV array.

WT Model
The Weibull PDF has been regularly used to model wind speed and can be expressed as The hourly output of a WT unit is highly dependent on wind speed. The mathematical expression for converting the hourly wind speed to power is as follows:

Utility Model
For MG operation, MGOs purchase deficit energy from the main grid depending on price and internal generation levels. Conversely, they sell excess energy to the main grid. The amount of power and exchange cost traded between the utility and the MG in each hour is expressed as: Here, the power generation state of the MG b u (t) takes the value 1 when there is a deficit in power, and 0 when there is a surplus of power.

Problem Formulation
In the optimization for the CEED problem, two competing objective functions are mathematically formulated by nonlinear functions, which minimize both generation costs and emission function by fulfilling the equality and inequality constraints. In other words, the problem is formulated for the multiobjective function as the minimization of the operation costs and pollutant emissions. The emission function is converted to cost by a multiplying with a weight called the penalty factor.

Generation Cost Function
The generation cost function is the sum of costs for each generator, including renewable energy source (RES). Conventional generators such as thermal generating units are generally expressed as the sum of a quadratic and sinusoidal function. Costs for PV and WT units, such as the available RES, are determined by considering investment and maintenance costs. The operation costs can be formulated as: P = [P G,1 , . . . , P G,I , P s,1 , . . . , P s,J , P w,1 , . . . , P w,K ] T (11) Energies 2020, 13, 6450 6 of 20 Equation (11) denotes the vector of the real power output of each generator. Equations (12) and (13) represent the costs of power generation with renewable energy sources, which include installation and maintenance costs. AC is an annuitization coefficient, and it can be calculated using the following equation [32]: 3.

Emission Function
The quantity of atmospheric pollutants, such as CO 2 , SO x , and NO x emitted by a conventional generator can be expressed as the sum of a quadratic and exponential function: In Equation (15), the emission dispatch function is considered as a convex polynomial, similar to the operation cost function. In addition, it is assumed that RES is not considered in the emission function because it does not release air pollutants [33].

Combined Economic Emission Dispatch
As discussed above, the generation costs and emission dispatch function are two different objectives. Hence, a compromise solution is required to solve the CEED problem that minimizes generation costs and the emitted quantities of pollutants. To solve a multiobjective function, a penalty factor is used to reform the emission criteria into the generation cost. Mathematically, the penalty factor is a multiplied value that transforms two different functions into a single objective function. In our work, a DR program is considered for microgrid operation. Therefore, the costs of participating in the DR program are included in the CEED problem, as follows.
Here, the penalty factor h i indicates the ratio ($/kg) of the generation costs and emission quantity of each generator i. This factor is initially determined by considering the minimum and maximum of each objective function, and then it is updated to account for environmental constraints. Equation (13) indicates the costs of participating in the DR program; it includes incentive and penalty costs. In this paper, it is determined based on the proposed EBDR program.

Generation Capacity Constraint
The power output of each generation unit is restricted by lower and upper limits for stable operation: Energies 2020, 13, 6450 7 of 20 Equations (18)-(21) indicate the minimum and maximum power limits of DG, utility, PV, and WT, respectively.

Power Balance Constraint
The total power generation must cover the total power demand in the presence of transmission line loss: Equation (23) represents the power demand that is finally determined through the DR.

DR Capacity Constraints
The DR capacity is determined by the contract between the MGO and the customers, and it cannot exceed the maximum capacity:

Transmission Line Constraints
For stable operation, MGO should consider the transmission line constraint. The upper limit of the transmission line is as follows:

Emission Constraints
The proposed CEED problem takes into account maximum emission constraints to prevent air pollution. If the MGOs violate emission constraints, they will have to pay a fee for the volume of excess.

Environment-Based Demand Response Program
This paper proposes a new type of incentive DR program called EBDR with emission control measures by including emission constraints. Generally, aiming to reduce peak demand, the MGO uses DR programs to adjust consumer behavior and manage power demand. The DR program includes the modification of electricity consumption patterns and incentives to promote this change. Demand and pollutant emission are directly related during MG operation. Thus, when solving the CEED problem, the EBDR is used to control pollutant emissions as well as reduce peak demand for economic benefits. By utilizing EBDR, MGOs can operate the MGs in an ecofriendly manner.

Price Elasticity of Demand
The EBDR includes the concept of elasticity based on market price. Elasticity is a measure used in economics to assess the percentage of change in demand caused by price fluctuations. In terms of power consumption, this percentage changes as the power demand varies with the changes in market prices, which is defined as an increase in price over time and not an absolute value. Elasticity is expected to be negative because higher power prices can cause demand reductions. The elasticity of the demand for electricity is calculated as follows: Energies 2020, 13, 6450 8 of 20 Two types of elasticity of demand exist: self-price elasticity and cross-elasticity [34,35]. Self-price elasticity disregards the shift in demand from one period to another and considers changes in consumption over a given period. In this case, increments in market price lead to demand reduction, and self-elasticity always has a negative value. Meanwhile, cross-elasticity is the transferred demand from the peak period to the off-peak demand period within a day. When t 1 is not equal to t 2 , the price drop at t 2 causes a demand reduction at t 1 . Thus, cross-elasticity always has a positive value.

Modeling of EBDR
The MGO enters into contracts with the customers through the EBDR program and mentions the values of incentive and penalty. The customers change their demand from D 0 (t) to D(t) according to the contract. The values of incentive and penalty for DR are given by Equation (29) represents the shift in demand of a customer due to the contract. Equations (30) and (31) denote the total values of incentive and penalty, respectively.
The EBDR program shifts the peak demand to reduce operation costs using cross-elasticity based on market prices. The relationship between market prices and demand is represented as follows: Equation (32) represents the power demand at t 1 according to the market price at t 2 , considering a 24-h interval. When incentive and penalty are included in the price, the formula is modified as follows: Equation (33) indicates only the change in demand due to the market price, without considering emission constraints. The change in demand when emission constraints are included is calculated by utilizing the concept of self-elasticity as follows: Here, emission constraints are imposed when emissions are greater than the set emission limit. The above equation yields the 24-h interval consumption of customers participating in the EBDR program, in order to minimize MG operating costs while considering emission constraints.

WU-ABC Optimization
In this study, we propose a solution through the WU-ABC algorithm, a hybrid of the weighting update method and the ABC algorithm, to solve the CEED optimization problem. The ABC algorithm is a swarm-based optimization method that simulates the foraging behavior of a bee swarm; it has recently been used in many studies to solve optimization problems [36,37]. However, it is difficult to solve the proposed CEED problem under emission constraints using the ABC algorithm alone. Therefore, the weighting update method is also employed to solve the emission problem by adjusting the penalty factor and EBDR capacity. The WU-ABC algorithm involves five main steps, as described below: • Initialization step: In this step, the initial population of food source is placed randomly in a D-dimensional problem space.
Equation (35) represents the m th random food source of dimension n in the CEED problem. • Employed bees step: Each bee repeatedly explores a food source to find the optimal solution, and then it chooses a new optimal position (v mn ) close to the reference position.
When an employed bee finds a food source that is better than the reference one (x mn ), the reference food source is superseded with the new candidate. • Onlooker bees step: Onlooker bees look for new positions that are close to the old position through a greedy search method. The bees consider a fitness value, the amount of nectar, and probabilistically determine the space for their next exploration, as follows: In Equation (37), fit m is a value that is proportional to the amount of nectar. • Scout bee step: If the reference solution cannot be improved through the number of predetermined trials, this food source is abandoned, and the scout bee tries to find a new food source to replace the abandoned food source. The new food source in the solution space is determined using Equation (35). This process is repeated for maximum cycles. • Weighting update step: This step is executed when an optimal solution violates the emission constraint. The penalty factor is updated as much as the actual pollutant emission exceeds the maximum emission, through a weight adjustment method. This process is repeated until the pollutant emission constraints are satisfied for optimal MG operation. Multiplication by a linear function is generally used as the weight update rule. In our work, to improve the convergence speed and accuracy, we use multiplication with an exponential function for the weighting update, as follows: After updating the penalty factor using Equation (38), optimization is performed again through the ABC algorithm. Step 1: Construct the MG model comprising DG, PV, and WT units and controllable demand, as shown in Figure 1, and initialize the input data for Equations (1)-(9).
Step 4: Set up parameters of the ABC algorithm, including the number of populations, maximum cycle, and maximum exploration.
Step 5: Initialize the power generation capacity of each unit randomly based on the initialization step in Equation (35).
Step 6: Search the neighborhood of the refence position through the employed bees, using Equation (36).
Step 7: Subject to Equation (37), apply a greedy search and choose one position to subsequently search for new neighboring positions.
Step 8: Replace the abandoned solution with the new solution found by the scout bees.
Step 9: Check that the result does not violate the emission constraint.
Step 10: If the emission constraint is violated, adjust the penalty factor in Equation (38).
Step 11: Repeat this process until the MG has satisfied all constraints.

System Data
To verify the effectiveness and validity of the proposed approach, a simulation is conducted on the modified grid-connected MG system consisting of three DG units, one 40-MW PV system, and one 30-MW wind farm [25]. Table 1 lists each DG characteristic, including the min/max power limit,

System Data
To verify the effectiveness and validity of the proposed approach, a simulation is conducted on the modified grid-connected MG system consisting of three DG units, one 40-MW PV system, and one 30-MW wind farm [25]. Table 1 lists each DG characteristic, including the min/max power limit, generation costs, and emission coefficients [32]. The operation costs for the wind farm and PV system include maintenance and investment costs of 30.8$/MWh and 23.4$/MWh, respectively [32].  Table 2 depicts the hourly demand and generation power of RES, which are calculated for wind speeds and solar irradiation in the east coast of the USA [32]. Table 3 shows the market price and DR incentive over the 24-h period [38]. The pollutant emission violation fee is set to 6.34$/kg [39]. The maximum DR capacity is assumed to be 10% of the total demand for each hour, while the maximum power flow capacity between MG and the main grid is 30 MW. The proposed WU-ABC algorithm has been applied to solve the associated CEED problem. To that end, we used MATLAB R2020a installed on a personal laptop with a 2.90 GHz core i5-9400F processor and 16 GB RAM. A simulation is run with 30 food sources and 300 iterations for 100 repeated trials, and the weighting update rate is set to 1.

Comparison Analysis
The U.S. plans to reduce CO 2 emissions from fossil fuel power plants by 32% before 2030 [39]. Thus, in this study, we assumed that emissions should ideally reduce by 32% when compared with an optimal operation situation of DG without consideration of emissions. Table 4 shows the optimal operation results when operating the MG with only DG units. The total operation costs and pollutant emissions are $101,065.60 and 5064.76 kg, respectively. Thus, maximum emission is set at 3444.04 kg, which is 32% of the pollutant emissions from DG units. To reduce operation costs and to avoid the violation of the maximum emission constraint, the proposed EBDR and WU-ABC algorithms are used for solving the CEED problem. The penalty factor should be set taking generation costs and the emissions of each generator into account. Various types of penalty factors are calculated and listed in Table 5. These factors vary according to the characteristics of each generator. In multiobjective problems, the penalty factor affects the relative importance of different objective functions, because it affects the optimal solution. In other words, a penalty factor that is excessively large may overestimate the effect of one objective function in a multiobjective problem. Thus, a min-max relationship with a small average penalty factor is chosen to prevent predominant effects of only one objective function, and the penalty factor is iteratively updated from the initial value to find the best solution that satisfies all constraints.  In this study, three cases are considered to validate the proposed approach. Case 1 is a conventional CEED proposition without a DR program and penalty factor update. Case 2 involves a conventional economic DR program but does not include penalty factor updates. In Case 3, the proposed EBDR is used and the penalty factor is updated through the WU-ABC algorithm. All cases are simulated under the same constraints, except for the DR programs and penalty factor update. Figure 4 shows the hourly scheduling of optimal power generation for all cases. To reduce operation costs, the MGO sells power for utility when market prices are high, in all cases. Although the generation cost of G1 is higher than that of the two other generators, it generates the most power due to low pollutant emissions per MWh. As shown in Figure 4a, Case 1 does not shift the peak demand; operation mainly occurs through the power generated by DG units, because no DR program is considered. Compared to Case 1, Case 2 utilizes a conventional economic DR program and demonstrates demand shifts in consideration of elasticity according to market prices. As shown in Figure 4b, the MGO uses the DR program to shift demand from periods of high market price to periods when power is relatively inexpensive. Thus, the power generation of DG units in the MG is reduced, and the operation costs are reduced despite considering incentives for the DR participants. However, in Case 2, no maximum emission constraints are imposed; to remedy this problem, penalty factor update and EBDR are considered in Case 3. The results for Case 3 are shown in Figure 4c. Compared to Case 1 and Case 2, it can be noted that the period taken by the MGO to sell power is reduced and the number of participants in the DR program is increased. Moreover, the power generation capacity of G2 and G3 (which have relatively large pollutant emissions per unit power generated) is also decreased. Case 3 could satisfy the maximum emission constraint. Table 6 shows the optimal CEED results for each case. Compared to the results in Table 4, which are essentially optimal operation propositions without environmental considerations, it can be seen that Case 1 has marginally higher operation costs (by 1.1%) but shows significant reduction in pollutant emissions (by 26.1%). These results violate the maximum emission constraint and there is no significant reduction in operation costs. For Case 2, the MGO reduced the operation costs of DG by using an economic DR to reduce peak demand. Thus, additional DR incentive costs are incurred, but the total operation costs are reduced by selling more power during high market price periods. The operation costs are reduced to $98,187.99, and pollutant emissions are decreased to 3571.14 kg. These results show that scheduling under a DR program for MG operation can have a positive impact, leading to reductions in peak demand, operating costs, and pollutant emissions. Compared to Case 1, Case 2 reduces both operation costs and pollutant emissions, but it does not satisfy the maximum emission constraints. Therefore, the MGO pays $805.81 as a violation fee, which is proportional to the excess pollutant emissions. In Case 3, the proposed EBDR and penalty factor update method are used to satisfy the maximum emissions constraint. Compared to Case 2, the operation costs of each DG unit are reduced due to the reduction of the power generation capacity of the units. The quantity of power that is sold to the utility decreased. The cost of participating in the DR program increased with the use of EBDR. Consequently, operating costs increased by approximately 0.05% compared to Case 2, but the maximum emission constraint is satisfied, curtailing emissions within 3444.03 kg. As can be seen from these results, the proposed approach is effective in solving the CEED problem with a maximum emission constraint.

Max-Min
Energies 2020, 13, 6450 14 of 20 penalty factor update and EBDR are considered in Case 3. The results for Case 3 are shown in Figure 4c. Compared to Case 1 and Case 2, it can be noted that the period taken by the MGO to sell power is reduced and the number of participants in the DR program is increased. Moreover, the power generation capacity of G2 and G3 (which have relatively large pollutant emissions per unit power generated) is also decreased. Case 3 could satisfy the maximum emission constraint.   Table 6 shows the optimal CEED results for each case. Compared to the results in Table 4, which are essentially optimal operation propositions without environmental considerations, it can be seen that Case 1 has marginally higher operation costs (by 1.1%) but shows significant reduction in pollutant emissions (by 26.1%). These results violate the maximum emission constraint and there is no significant reduction in operation costs. For Case 2, the MGO reduced the operation costs of   Figure 5 indicates the penalty factor and pollutant emissions with respect to the number of updates. The penalty factor is updated proportionally to the excess pollutant emissions; consequently, the penalty factor update is performed 11 times, and the penalty factor is finally determined to be 23.55$/kg. Table 7 summarizes the optimal CEED solution for all cases. Case 1 without a DR program presents a greater total operation cost and higher emissions than the other cases. Meanwhile, for Case 2 with the conventional economic DR program, the MGO is able to reduce operation costs and pollutant emissions relative to Case 1 but pays a violation fee for exceeding the maximum emission constraint. In Case 3, the maximum emission constraint is satisfied through the use of penalty factor update and the EBDR program; although it slightly increases the operating cost compared to Case 2, the constraint violation cost of $805.81 is not incurred. Therefore, it can be concluded that the proposed approach exhibits the most promising results for the CEED problem, when considering a violation fee.  To demonstrate the superiority of WU-ABC, performance tests have been conducted through comparison with various algorithms [40,41]. For a fair comparison, the maximum emission constraints and violation fee are equally considered in grid-connected MG operation. A total of five algorithms are considered, such as ABC, improved multilayer artificial bee colony (IML-ABC) [41], weighting update genetic algorithm (WU-GA), weighting update particle swarm optimization (WU-PSO), and proposed WU-ABC. WU-GA and WU-PSO are the addition of the weighting update method to the conventional GA and PSO. The simulation process for all algorithms is repeated 10 times, and then the results are calculated as average value. Table 8 shows the CEED results of each algorithm. It can be observed that WU-ABC reaches a best solution compared with other algorithms. The algorithms without the weighting update method do not satisfy the emission  To demonstrate the superiority of WU-ABC, performance tests have been conducted through comparison with various algorithms [40,41]. For a fair comparison, the maximum emission constraints and violation fee are equally considered in grid-connected MG operation. A total of five algorithms are considered, such as ABC, improved multilayer artificial bee colony (IML-ABC) [41], weighting update genetic algorithm (WU-GA), weighting update particle swarm optimization (WU-PSO), and proposed WU-ABC. WU-GA and WU-PSO are the addition of the weighting update method to the conventional GA and PSO. The simulation process for all algorithms is repeated 10 times, and then the results are calculated as average value. Table 8 shows the CEED results of each algorithm. It can be observed that WU-ABC reaches a best solution compared with other algorithms. The algorithms without the weighting update method do not satisfy the emission constraint, and as a result, the total cost is high due to imposing violation fees. In addition, WU-ABC indicates the lowest cost compared to other algorithms considered weighting update methods and takes the least CPU time. Therefore, these results reveal that the proposed WU-ABC algorithm is appropriate for solving CEED problems considering maximum emission constraints.

Conclusions
This paper proposed a CEED approach with the EBDR program by using the WU-ABC algorithm. The objective function was constructed in consideration of generation costs and pollutant emission. The emission function took the emitted quantities of CO 2 , SOx, and NOx, into account and converts their sum into operation costs through unit conversion, using a penalty factor. To meet the maximum emission constraint considered in addition to the CEED in our work, the EBDR program was proposed to balance economics and environmental issues. The EBDR program is based on the concept of elasticity, and the participation volume changes depending on the violation of emission constraints. To solve the multiobjective optimization problem, we proposed the use of the WU-ABC algorithm, a combination of the weighting update method and the conventional ABC algorithm. To demonstrate the effectiveness of the proposed CEED approach, simulations were conducted on the modified grid-connected MG system. Three cases were considered depending on the maximum emission constraint and the application of a DR program. The proposed approach was able to reduce operation costs compared to the conventional solution without DR program. Furthermore, it was easy to satisfy the maximum emission constraints by updating the penalty factor. In summary, the use of the EBDR program reduced operation costs and pollutant emissions, and the maximum emission constraints were satisfied through the WU-ABC algorithm. Although operation costs were slightly increased compared to that of a conventional economic DR, an emission violation fee was not imposed. In addition, the performance test results indicated that the WU-ABC algorithm outperformed other algorithms in terms of cost and CPU time. Therefore, the proposed CEED approach can assist MGOs in optimizing MG operations under environmental constraints. Our future work will be under way to focus on solving the sensitivity problem in MG when the consumers do not comply with the command of MGO in the CEED problem.   Installation cost of j PV panels G w,k Installation cost of k WT units Penalty factor Variables C G,i (t) Generation cost function of i DG units at t P G,i (t) Generation power of i DG units at t P min G,i (t) Minimum generation power of i DG units at t P max G,i (t) Maximum generation power of i DG units at t PDF B Beta PDF Γ Gamma function T t Temperature at t P s, j (t) Generation power of j PV panels at t P min s, j (t) Minimum generation power of j PV panels at t P max s, j (t) Maximum generation power of j PV panels at t PDF w Weibull PDF v t Wind speed at t P w,k (t) Generation power of k WT units at t P min w,k (t) Minimum generation power of k WT units at t P max w,k (t) Maximum generation power of k WT units at t P u (t) Power exchanged with utility at t P min u (t) Minimum power exchanged with utility at t P max u (t) Maximum power exchanged with utility at t P buy (t) Power purchased from utility at t P sell (t) Power sold to utilities at t b u (t) Unit function C u (t) Power exchanged cost with utility at t λ sell (t) Power purchased cost from utility at t λ buy (t) Power sold costs to utilities at t F(P) Generation costs function C s, j Solar power cost C w,k Wind power cost E(P) Emission function em(t) Pollutant emissions at t em max (t) Maximum pollutant emissions at t C(P) CEED function C DR (t) Total DR cost at t I(∆D(t)) Total Incentive cost P(∆D(t)) Total Penalty cost P D (t) Total power demand at t P loss (t) Total line loss at t D 0 (t) Initial power demand at t D(t) Power demand at t ∆D(t) Change in power demand S l Power flow at line l S max l Maximum power flow at line l E(t 1 , t 2 ) Elasticity for t 1 and t 2 MP 0 (t) Initial market price at t MP(t) Market price at t inc(t) DR incentive cost at t pen(t) DR penalty cost at t CP(t) Contracted power capacity at t D e (t) Power demand considering emission E e (t 1 , t 1 ) Emission elasticity