Economic Emission Dispatch for Wind Power Integrated System with Carbon Trading Mechanism

Nowadays, the power system is faced with some new changes from low-carbon approaches, though these approaches have proved to be effective in developing low-carbon electricity. Specifically, wind power integration and carbon trading influence the traditional economic emission dispatch (EED) mode, allowing for the disturbance of wind power uncertainties and the fluctuation of carbon trading price. Aiming at the above problems, this study firstly builds a stochastic EED model in the form of chance-constrained programming associated with wind power reliability. Next, wind power features are deduced from the statistic characteristics of wind speed, and thus the established model is converted to a deterministic form. After that, an auxiliary decision-making method based on the technique for order preference by similarity to an ideal solution (TOPSIS) is designed to draw the optimal solution based upon the specific requirements of carbon emission control. The simulation results eventually indicate that the minimization of fuel costs and carbon emissions comes at the expense of wind power reliability. Meanwhile, carbon emission reduction can be effectively realized by carbon trading rather than a substantial increase in fuel costs, and carbon trading may help to improve power generation efficiency. Furthermore, carbon trading prices could be determined by the demands of carbon emission reduction and power generation efficiency improvement.


Introduction
With the rapid development of the social economy and the continuous progress of science and technology, people's demand for electric power resources is increasingly urgent [1]. However, electrical production not only needs to consume a large amount of traditional energy, but also releases carbon dioxide (CO 2 ), which causes the greenhouse effect [2]. Against this background, low-carbon electricity development naturally becomes one of the important measures to accelerate the sustainable development of the electric power industry [3]. In recent years, many low-carbon approaches have been carried out in power systems [4].
In the process of power generation, wind power as a kind of renewable energy has the advantages of being renewable, environmentally friendly and cheap, so a large amount of wind power integration into the power system has become an effective low-carbon approach [5,6]. China's installed wind power capacity had already reached 281 million kilowatts by the end of 2020, of which offshore wind power generation occupied about 9 million kilowatts [7]. While in the process of power dispatching, based on guaranteeing the safe operation of the power grid, economic emission dispatch (EED), synthesizing fuel costs and CO 2 emissions simultaneously, has become a commonly used power dispatching technology [8]. It not only needs to meet the safety and economic requirements of traditional power dispatching, but also to reduce CO 2 emissions as much as possible [9]. Obviously, EED, with its many advantages, is identified as a more practical approach for carbon electricity [23]. Specifically, wind power integration and carbon trading influence traditional EED mode, allowing for the disturbance of wind power uncertainties and carbon trading price fluctuations. As one of the essential processes in low-carbon electricity development, exploring the effects of wind power integration and carbon trading on EED becomes especially significant. To achieve this, it is quite meaningful to investigate the EED problem in a wind power integrated system with a carbon trading mechanism. Based on the above considerations, this paper tries to establish how wind power reliability affects the optimal dispatching solution, how carbon emission reduction and power generation efficiency improvements are influenced by carbon trading price and how to set the rational price of carbon trading according to the specific requirements of carbon emission reduction and power generation efficiency improvement. Concentrating on wind power integration, EED and carbon trading simultaneously, together with the relationships between them, this paper may help the decision-makers to choose low-carbon power dispatching strategies more scientifically and rationally.
This paper is organized as follows: Section 2 builds a stochastic EED model considering wind power reliability; Section 3 studies the statistical characteristics of wind energy, and the effect of decision-making methods upon carbon trading price; Section 4 presents a case study to investigate the relationships between wind power integration, EED and carbon trading; and Section 5 summarizes the full text and provides some conclusions.

Problem Formulation
In consideration of the effect of wind power uncertainties on the optimal dispatch strategy for wind power integrated systems, the traditional EED model is required to be modified to describe this effect. To be specific, the probabilities of real power balance constraints related to wind power uncertainties are defined as confidence levels to reflect wind power reliability. In other words, the present formulation treats this EED problem as a chance-constrained programming model for minimizing both fuel costs and carbon emissions under uncertainty, and to satisfy both load demand and system operation constraints.

Cost Function
Assuming there is no cost once the wind power generators are up and running, the major consideration of generation cost is the thermal power generator's cost in the wind power integrated system, and the cost function is constructed as follows [24]: where F 1 is the total costs of the power system; C i is the thermal power generator's cost function; P i is the scheduling output for the thermal power generator; and N is the number of thermal power generators. Generally, the cost function is a quadratic function for a single thermal power unit. In the actual operation of the system, the phenomenon of wire drawing superimposes a pulsating effect on the characteristic curve of unit consumption; that is, the valve point effect. In order to improve the accuracy of the EED model, a sine function is usually added to the cost function which is further developed as shown below [25]: where a i , b i , c i , e i and f i are the cost function coefficients for i, the thermal power generator; and P min i is the thermal power generator's lower limit of power generation.

Emission Function
The wind power generators do not produce carbon emissions directly, as wind power is clean energy, while the thermal power units do produce some carbon emissions, thus the emission function is shown below [26]: where F 2 is the total carbon emissions of the power system; and E i is the emission function for the thermal power generator. The total carbon emissions arising from thermal power generators could be depicted as a quadratic function in the following form [27]: where α i , γ i and λ i are the thermal power generator's emission coefficients.

Real Power Output Constraint
In order to ensure the normal operation of the thermal power generator, the requirements of the upper and lower limits of power generation ought to be satisfied [28]: where P max i is the power generation upper limit for the thermal power generator.

Real Power Balance Constraint
Due to the wind power uncertainties, wind power output ought to be considered as a random variable, and the constraints of the power equilibrium which the wind power integrated system should satisfy under certain confidence levels can be formulated as the probability forms [28]: where P w is the wind power output; P load is the total load demand; P loss is the transmission line losses; and α is the confidence level, thus the reliability of the wind power integrated system can be defined as 1 − α.
In particular, P loss in Equation (6) based on Kron's loss equation can be described as follows [29]: where B ij , B i0 and B 00 are B-coefficients of transmission network power loss.

Uncertainty Processing
With respect to the proposed EED model in the form of chance-constrained programming, a real power balance constraint appears as a random event on account of the wind power uncertainties, and this model should be switched to a deterministic form for the ease of the optimization process [30]. Specifically, the statistical characteristics of wind power output could be deduced from the probability distribution function (PDF) of wind speed.
There are mainly five kinds of wind speed PDF: Weibull distribution, gamma distribution, Rayleigh distribution, log-normal distribution and Burr distribution, which are widely applied in wind speed forecasting.
The Weibull PDF can be defined as shown below [31]: where v is the wind speed; c denotes the scale parameter (the same dimension of wind speed); and k represents the shape parameter (dimensionless). The gamma PDF can be written as indicated below [31]: where ζ is the shape parameter; β is the scale parameter; and Γ is the gamma function. The Rayleigh PDF can be formulated as follows [32]: where τ is the scale parameter. The log-normal PDF can be defined as follows [32]: where µ is the location parameter; and σ is the scale parameter. The Burr PDF can be written as follows [32]: where ε is the shape parameter; and ω and ϕ are the scale parameters. The Weibull distribution in a relatively simple form has the widest application range, which is explicitly for wind speed, and convenient to calculate. However, its fitting effect is only good for ordinary wind speed, for it cannot fit some wind speed characteristics at extreme values. The gamma distribution is the first distribution model used to depict wind speed. The Rayleigh distribution sets a high standard for data, which requires not only mass data but also for the dispersion degree of the data to be as small as possible, and the fitting error is usually around 10%. The fitting effect of the log-normal distribution is generally good, but its effect is poor when the wind speed frequency is too low or high. The Burr distribution also has excellent results, but this model is a little complex and its solution seems a bit complicated. In summary, every distribution has its own advantages and disadvantages, so it is necessary to further explore which one is more suitable for the statistical characteristics of wind speed in the actual application.
Based on the wind speed PDF, the statistical characteristics of wind power output may further be determined, and the corresponding transformation could be formulated as a linear relationship [33]: where W is wind power output; V is wind speed; v r represents the rated wind speed; v i denotes the cut-in wind speed; and T signifies a transformation, and where w rated is rated wind power. Relying on Equation (13), the PDF of wind power output together with its cumulative distribution function (CDF) can be induced as below: After discussing the statistical characteristics of wind power output, the constraint in Equation (6) can be converted into a deterministic form: Taken together, uncertain constraints in the proposed EED model are quantitatively deduced, and finally, a deterministic model is established to obtain the optimal dispatch strategy for a wind power integrated system.

Decision-Making
With regards to the above deterministic EED model, such a multi-objective optimization problem has multiple objectives and each objective function is contradictory. That is, the increase in proximity of the cost function to the optimal value will cause the emission function to move away from its optimal value. The solutions of the proposed EED problem are a set of non-inferior solutions, since no solution is better than any of the others with respect to multi-objectives [26]. For obtaining the optimal dispatch strategy, a set of "best" solutions could be selected from these non-inferior solutions via specific evaluation methods, such as inferior solution distance method (TOPSIS), analytic hierarchy process, reordering method, superior and hierarchical sequence method, etc. [34,35]. By virtue of TOPSIS, this section manages to draw the optimal solution as the decision-making scheme, and therein the definitions of the weight of each objective function become a key problem. The critical decision steps which can provide convenience for the selection of an optimal dispatch strategy, and specifically the relationship between carbon trading prices and the optimal objectives' weights are devised as follows [36]: Step 1: Supposing that the obtained Pareto-efficient solution set is , m is the number of the objective function corresponding to the solution. Calculating the normalized decision matrix ( f im ) K×2 for the purpose of simply calculation, where f im is the normalization of F m (x i ): Step 2: Assuming that the weight of the cost function F 1 is ω 1 , and the weight of the emission function F 2 is ω 2 , these two weights may fluctuate with the carbon trading price P CDM : the higher the P CDM , the more benefits (thermal power generators with lower emissions) or fees (thermal power generators with higher emissions) the power system has to afford. That is, P CDM is positively correlated with ω 2 , which is inclined to carbon emission control. ω 2 cannot be defined as P CDM directly, for the latter has to be normalized: , where ρ is a compensation coefficient relying on the current average P CDM , and F Unit 1 , F Unit 2 stand for per unit cost and per unit emission respectively.
Step 3: Defining the ideal solution A + and the non-ideal solution A − : Step 4: Generating the separation measures d + i and d − i for each solution based on A + and A − , together with ω 1 and ω 2 : Step 5: Figuring out each solution's relative distance R i , which reflects the distance from each solution to the non-ideal solution, and ultimately extracting the Pareto-efficient solution with the largest relative distance R j 0 :

Case Study
From the discussions in Sections 2 and 3, the modeling and implementation of power dispatching problems have been adequately presented, and the proposed model with general forms has been identified as suitable for current low-carbon development by setting various parameters of wind power output and carbon trading price. By virtue of the multi-objective optimization algorithm and multi-attribute decision-making method, this study tries to further study how wind power reliability affects the optimal dispatching solution, how carbon emission reduction and power generation efficiency improvement are influenced by carbon trading price, and how to set the rational price of carbon trading according to the specific requirements of carbon emission reduction and power generation efficiency improvement.
In this study, together with one wind power generator, a classic six-unit system with an emission level and non-smooth fuel cost was applied to test the performance of the proposed EED model [37]. The test system's operating coefficients are offered in Table A1 of the Appendix A, and the parameters of the wind power unit are supplied in Table A2. Furthermore, the system's demand was set as 900 MW, and the transmission loss matrix is also supplied in the Appendix A.

Description of Wind Power Output
In order to simulate the distribution function of wind power output relying on wind speed, one month's wind speed data from Rudong, East China's Jiangsu Province, was collected to describe the statistical characteristics of wind speed [38]. More specifically, the maximum likelihood method was used to estimate the parameters of different wind speed PDFs as discussed in Section 3.1, and the estimated values of each distribution's parameters are shown in Table 1: Based on Table 1, fitted curves for these five kinds of wind speed PDFs are plotted in Figure 1, where the histograms of actual wind speed are regarded as a reference picture. Specifically, the step size was set as 0.05 m/s and a series of equilateral points were selected from the wind speed interval (0 and 10 m/s). To compare the simulation effects from different wind speed PDFs, the Euclidean distance between the estimated probability density of each curve and the estimated probability density of the histogram was calculated, as shown in Table 2. Specifically, the probability density function values of various curves were calculated respectively at the same equal-point abscissa value derived from the estimated histogram probability density. It can be seen from Table 2 that the Euclidean distance between the log-normal distribution and histogram was the smallest, therefore the fitting effect of the log-normal distribution is considered to be the best one for the current wind speed data. At the same time, since the classical distribution of the wind speed forecast is the Weibull distribution, the linear weighting method combining the log-normal distribution and the Weibull distribution is formulated as the final form of wind speed distribution: To compare the simulation effects from different wind speed PDFs, the Euclidean distance between the estimated probability density of each curve and the estimated probability density of the histogram was calculated, as shown in Table 2. Specifically, the probability density function values of various curves were calculated respectively at the same equal-point abscissa value derived from the estimated histogram probability density. It can be seen from Table 2 that the Euclidean distance between the log-normal distribution and histogram was the smallest, therefore the fitting effect of the log-normal distribution is considered to be the best one for the current wind speed data. At the same time, since the classical distribution of the wind speed forecast is the Weibull distribu-tion, the linear weighting method combining the log-normal distribution and the Weibull distribution is formulated as the final form of wind speed distribution: where the weight is θ, which can be defined as follows: where d 1 is the distance between the Weibull distribution and histogram, and d 2 is the distance between the log-normal distribution and histogram, thus θ depending on d 1 and d 2 can be quantitatively fixed as 0.4952, based on Table 2. Furthermore, the CDF of the wind power output in Equation (17), which can provide convenience when dispatching optimization, can be finally deduced in the following form: With regards to the simplified restriction in Equation (17), it should be noted especially that the standard normal distribution function in Equation (28) could be approximated by an analytic expression of the scheduling output P i [39], and then the upper and lower bounds related to the restriction can be achieved the traversal way.

Impact Analysis of Wind Power Reliability
After simulating the wind power output, the proposed EED model was converted into a deterministic form. Next, with reduction in the total fuel costs and carbon emissions of the power system as optimization goals, the proposed EED model was minimized by NSGA-II through MATLAB software. To guarantee the diversity and convergence of the evolution of the population in the process of optimization, the optimal front-end individual coefficient was set as 0.6, and the population size was set to 100. Moreover, both the maximum evolutionary algebra and maximum number of iterations were set as 200, and the deviation of the fitness function value was fixed at 0.01. In order to explore how wind power reliability affected the optimal dispatching solutions, according to various reliability levels indicated by different confidence levels, the distributions of 60 non-dominated solutions representing three different levels of reliability were calculated and are depicted in Figure 2. gies 2021, 14, x FOR PEER REVIEW 10 o Furthermore, the CDF of the wind power output in Equation (17), which can provi convenience when dispatching optimization, can be finally deduced in the followi form: With regards to the simplified restriction in Equation (17), it should be noted esp cially that the standard normal distribution function in Equation (28) could be appro mated by an analytic expression of the scheduling output i P [39], and then the upp and lower bounds related to the restriction can be achieved the traversal way.

Impact Analysis of Wind Power Reliability
After simulating the wind power output, the proposed EED model was convert into a deterministic form. Next, with reduction in the total fuel costs and carbon emissio of the power system as optimization goals, the proposed EED model was minimized NSGA-II through MATLAB software. To guarantee the diversity and convergence of t evolution of the population in the process of optimization, the optimal front-end indiv ual coefficient was set as 0.6, and the population size was set to 100. Moreover, both t maximum evolutionary algebra and maximum number of iterations were set as 200, a the deviation of the fitness function value was fixed at 0.01. In order to explore how wi power reliability affected the optimal dispatching solutions, according to various reliab ity levels indicated by different confidence levels, the distributions of 60 non-dominat solutions representing three different levels of reliability were calculated and are depict in Figure 2.  From the above three figures, the distributions of the Pareto frontier which tracks both fuel costs and carbon emission fluctuations are always extensive and uniform. With the enhancement of the confidence level for chance constraint, the reliability of the wind power integrated system defined in Equation (6) becomes lower and lower, and directly impacts the optimal dispatching solutions. To illustrate the differences between the three figures in Figure 2 as far as two types of single-objective, Table 3 indicates the extreme values of fuel costs and carbon emissions in terms of different reliability levels. As indicated in Table 3, with an increase in confidence level, which leads to a decrease in the reliability of the wind power integrated system, the minimum fuel cost decreases from 57,337.63 to 56,927.77 dollars accordingly, and the minimum carbon emission simultaneously decreases from 5686.42 to 4529.59 lb. This is obviously the result of the improvement in the wind power generators' power output, and the reduction in the thermal power generators' power output. It suggests that the improvement in wind power proportion is conducive to the control of fuel costs and carbon emissions; however, it can reduce wind power reliability depending on its uncertainties. In other words, the minimization of fuel costs and carbon emissions comes at the cost of wind power reliability. To verify the accuracy of optimization, Table 4 lists all scheduled power generations for the final solution by taking α = 0.05 as an example. Table 4. Scheduled generation (MW), respective (total) cost ($) and carbon emissions (lb) (α = 0.05).

Impact Analysis of Carbon Trading Price
Having obtained the Pareto-optimal solutions under different reliability levels, it is necessary to select an appropriate method of decision-making via the specific evaluation method designed in Section 3.2, since no solution is superior to all the others with respect to multi-objectives. Furthermore, this section focuses on the impact analysis of the carbon trading price; hence, the reliability level as another parameter is assumed to be 0.95 (α = 0.05) in accordance with Table 4. With increasing carbon trading price P CDM , the satisfying solutions with the largest relative distance under different weights ω 2 , combined with the corresponding fuel costs and carbon emissions are summarized in Table 5. To clarify, ω 1 and ω 2 , as defined in Equation (19), are the normalizations of P CDM , and specifically, ω 2 is positively correlated with P CDM , thus the impact analysis of P CDM is roughly equivalent to the impact analysis of ω 2 .
As ω 2 increases, the fuel costs and carbon emissions of satisfying solutions vary monotonically in Table 5. For the purposes of normalized comparison, Figure 3 depicts the correlations between ω 2 and the rate of change of the fuel costs and carbon emissions.  increases, the fuel costs and carbon emissions of satisfying solutions vary monotonically in Table 5. For the purposes of normalized comparison, Figure 3 depicts the correlations between 2  and the rate of change of the fuel costs and carbon emissions. From the line charts in Figure 3, it can be found that the growth rate of fuel costs and the decrement rate of carbon emissions simultaneously increase as 2  increases, namely due to the improvement in carbon trading price, whereas these rates of change become almost stable as soon as 2  reaches 0.7. In particular, the decrement rate of carbon emissions fully exceeds the growth rate of fuel costs for any value of 2  between 0 and 1. For example, the fuel costs increased by 4.29 and 5.74%, while the carbon emissions separately reduced by 28.3 and 31.27%, respectively, when 2  varied from 0 to 0.3 and 0.7. In other words, carbon emission reduction can be effectively realized by carbon trading without significantly increasing fuel costs, considering that the decrement rate of carbon emissions is superior to the growth rate of fuel costs when the carbon trading price increases.
In addition, to find out how power generation efficiency improvement is influenced by carbon trading price, Figure 4 shows the proportions of power generation shared by six thermal power generators corresponding to different values of 2  . From the line charts in Figure 3, it can be found that the growth rate of fuel costs and the decrement rate of carbon emissions simultaneously increase as ω 2 increases, namely due to the improvement in carbon trading price, whereas these rates of change become almost stable as soon as ω 2 reaches 0.7. In particular, the decrement rate of carbon emissions fully exceeds the growth rate of fuel costs for any value of ω 2 between 0 and 1. For example, the fuel costs increased by 4.29 and 5.74%, while the carbon emissions separately reduced by 28.3 and 31.27%, respectively, when ω 2 varied from 0 to 0.3 and 0.7. In other words, carbon emission reduction can be effectively realized by carbon trading without significantly increasing fuel costs, considering that the decrement rate of carbon emissions is superior to the growth rate of fuel costs when the carbon trading price increases.
In addition, to find out how power generation efficiency improvement is influenced by carbon trading price, Figure 4 shows the proportions of power generation shared by six thermal power generators corresponding to different values of ω 2 .
With an increase in ω 2 , it is shown in Figure 4 that the power generation proportions shared by units 1 and 6 show a significant rising tendency, while the proportion shared by unit 2 appears to be an obviously decreasing trend. The reasons for this phenomenon are that the carbon emissions per-unit generating capacity from units 1 and 6 that are beneficial to carbon trading are smaller. As a consequence, ω 2 has an obvious impact on the power generation proportions of units 1 and 6: the higher the carbon trading price, the more output power they can afford. On the contrary, there is a negative correlation between ω 2 and the output power for unit 2, with a larger carbon emissions per-unit generating capacity. Specifically, the carbon emissions per-unit generating capacity among the six different generators is shown in Figure 5. Based on the above discussion and analysis, it shows that carbon trading may contribute to the improvement of power generation efficiency, as units with a lower carbon emissions per-unit generating capacity benefit far more than others from carbon trading. With an increase in 2  , it is shown in Figure 4 that the power generation proportions shared by units 1 and 6 show a significant rising tendency, while the proportion shared by unit 2 appears to be an obviously decreasing trend. The reasons for this phenomenon are that the carbon emissions per-unit generating capacity from units 1 and 6 that are beneficial to carbon trading are smaller. As a consequence, 2  has an obvious impact on the power generation proportions of units 1 and 6: the higher the carbon trading price, the more output power they can afford. On the contrary, there is a negative correlation between 2  and the output power for unit 2 ,with a larger carbon emissions per-unit generating capacity. Specifically, the carbon emissions per-unit generating capacity among the six different generators is shown in Figure 5. Based on the above discussion and analysis, it shows that carbon trading may contribute to the improvement of power generation efficiency, as units with a lower carbon emissions per-unit generating capacity benefit far more than others from carbon trading.   From what has been discussed above, it is found that carbon trading is beneficial for carbon emission reduction and power generation efficiency improvement. On the contrary, since carbon trading is implemented by carbon trading price setting, it is necessary to set a rational P CDM according to the specific requirements of carbon emission reduction and power generation efficiency improvement. Table 6 mainly illustrates the comparison between ω 2 and P CDM based on Equations (19) and (20), and Table A1, besides the relationships with P CDM , shows the decrement rate of carbon emissions and the proportions shared by units 1 and 6. The P CDM for ω 2 = 0.5, with a moderate view on environmental protection, is assumed to be 41.39 (¥/ton), which is an average price from the Beijing carbon trading market in 2019 [13], so that different P CDM corresponding to various ω 2 can be obtained after deriving the compensation coefficient ρ in Equation (19). In theory, P CDM depending on specific requirements could be set once the decrement rate of carbon emissions and the proportions shared by more environmental units are fixed in advance.

Conclusions and Recommendation
In order to solve the EED problem in wind power integrated systems with carbon trading mechanisms, this paper firstly builds a stochastic EED model in the form of chance-constrained programming. Next, the established EED model is converted to a deterministic form for the ease of the optimization process. After that, facing a set of non-inferior solutions for model optimization, an auxiliary decision-making method based on TOPSIS is designed to draw the optimal solution upon the specific requirements of carbon emissions control.
To find out how wind power reliability affects the optimal dispatching solution, how carbon emission reduction and power generation efficiency improvement are influenced by carbon trading price, and how to set the rational price of carbon trading according to the specific requirements of carbon emission reduction and power generation efficiency improvement, a classic six-unit system together with one wind power generator is introduced to settle this optimization problem in the case study. To begin with, the mixture of the log-normal and Weibull distribution functions with good forecasting performance is used to simulate wind power output, relying on actual wind speed data. Then, empirical analysis indicates that the improvement of wind power proportion is conducive to the control of fuel costs and carbon emissions; however, it will reduce wind power reliability depending on its uncertainties. In other words, the minimization of fuel costs and carbon emissions comes at the cost of wind power reliability. Moreover, carbon emission reduction can be effectively realized by carbon trading rather than a substantial increase in fuel costs, considering that the decrement rate of carbon emissions is superior to the growth rate of fuel costs when the carbon trading price increases. Additionally, carbon trading may contribute to the improvement of power generation efficiency, as units with a lower carbon emissions per-unit generating capacity will benefit far more than others from carbon trading. Furthermore, once the decrement rate of carbon emissions and the proportions shared by more environmental units are fixed in advance, carbon trading prices according to the specific requirements of carbon emission reduction and power generation efficiency improvement could be determined. In brief, within the framework of EED, the impacts of carbon trading and wind power integration on power dispatching are proved to be significant. With a lot of attention drawn to wind power integration, EED and carbon trading, the proposed model and the obtained conclusions can help decision-makers to choose low-carbon power dispatching strategies more scientifically and rationally. This paper focuses on the EED problem in wind power integrated systems with carbon trading mechanisms, and more attention is paid to the influences of wind power integration and carbon trading on EED. However, there might be some other low-carbon approaches to power systems which allow us to investigate its corresponding impacts. Additionally, it may further explore the present EED problem from various perspectives. Another potential direction is to pursue the other effective optimization algorithms to solve the EED problem in this paper.
Author Contributions: J.J.: conceptualization, methodology, investigation, validation and writingoriginal draft. Q.W.: investigation, data analysis, software and writing-review and editing. X.Z.: resources, data analysis, software and formal analysis. S.C.: writing-review and editing. X.G.: writing-review and editing. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
Variables: P i the scheduling output for the thermal power generator P w the wind power output P load the total load demand for the power system