Spinning Reserve Requirements Optimization Based on an Improved Multiscenario Risk Analysis Method

This paper proposes a novel security-constrained unit commitment model to calculate the optimal spinning reserve (SR) amount. The model combines cost-benefit analysis with an improved multiscenario risk analysis method capable of considering various uncertainties, including load and wind power forecast errors as well as forced outages of generators. In this model, cost-benefit analysis is utilized to simultaneously minimize the operation cost of conventional generators, the expected cost of load shedding, the penalty cost of wind power spillage, and the carbon emission cost. It remedies the defects of the deterministic and probabilistic methods of SR calculation. In cases where load and wind power generation are negatively correlated, this model based onmultistep modeling of net demand can consider the wind power curtailment tomaximize the overall economic efficiency of system operation so that the optimal economic values of wind power and SR are achieved. In addition, the impact of the nonnormal probability distributions of wind power forecast error on SR optimization can be taken into account. Using mixed integer linear programming method, simulation studies on a modified IEEE 26-generator reliability test system connected to a wind farm are performed to confirm the effectiveness and advantage of the proposed model.


Introduction
Wind power generation has been used in many countries as the most promising sustainable energy source to reduce the consumption of fossil fuels and curb emissions of carbon dioxide.However, because of the stochastic volatility and limited predictability of wind speed, power system operations are significantly challenged by the largescale wind power integration.The influences include power quality, operating reserve, transient stability, and frequency and voltage control [1,2].In particular, the traditional unit commitment (UC) with deterministic spinning reserve (SR) requirements cannot adequately cope with uncertain wind power generation.Additional SR needs to be provided to ensure operational reliability [3,4].Therefore, various new methods taking into account wind power uncertainty in [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20] have been presented to determine SR requirements.These methods can mainly be divided into three types: (i) The deterministic methods [5][6][7][8][9]: in [5][6][7][8][9], the minimum amount of SR is set to be the largest capacity of online generators, or some proportion of the standard deviation of the net demand or wind power forecast error, or their combination.Although these deterministic methods are easy to implement, they are inadequate in dealing with the stochastic nature of various uncertainties in power systems and the economic efficiency of system operation.
(ii) The probabilistic methods [10][11][12][13]: in the reliabilityconstrained UC [10,11,13] and method [12], probabilistic reliability criteria that are defined by the loss of load probability (LOLP) or the expected energy not served (EENS) implicitly determine the SR requirements.In each time period, the probabilistic method can ensure that the system will maintain a specified reliability level under which the optimal economic efficiency of system operation can be achieved.However, this method neglects the setting of the reliability metrics and the rationality of such metrics.
(iii) The methods based on cost-benefit analysis [14][15][16][17][18][19][20]: by introducing the notion of value of lost load (VOLL), the optimization of the SR requirements based on cost-benefit analysis conducts an economic value analysis for the cost and benefit of providing SR so as to achieve a tradeoff between economics and reliability.Thus, the experiential settings of the required SR amount [5][6][7][8][9] and reliability metrics [10][11][12][13] are avoided.
It is well known that wind power generation needs to be controlled by "spilling wind" owing to the transmission congestion.Due to the fact that there is no fuel consumption in wind power generation, wind power spillage, which is also called "wind power curtailment," increases the load supplied by conventional generators and thereby increases the fuel cost and emission cost of these generators.On the other hand, in cases where load and wind power generation are negatively correlated, wind power spillage will diminish the peak-valley difference of net demand so that the start-up and shutdown costs of peaking generators can be reduced [5].Moreover, wind power spillage reduces the underestimation of wind power output, which may decrease the system EENS and thus reduce additional SR amount.Therefore, besides security issues, these two conflicting and opposite effects of wind power spillage should be considered while optimizing SR requirements from the overall economic and low-carbon efficiency of system operation.
The UC in [5,6,[14][15][16][17][18] comprehensively consider the effects of wind power spillage on the total cost of system operation and SR quantification.Nevertheless, these studies have other various limitations in system EENS calculation.The following are specifically mentioned.
(i) Modeling: a stochastic planning method based on a scenario tree is adopted in [14][15][16] to simulate the probability distribution of forecast errors of net demand or wind power so that the EENS of every scenario can be determined.However, the forced outage rates (FOR) of generators are not taken into account, which may lead to an underestimation of SR.Similarly, Lee [17] used a capacity outage probability table to calculate the system EENS based on the simplified assumptions of wind power and load forecast errors.
(ii) Optimization process: Ortega-Vazquez [18] considered the load and wind power forecast errors as well as FOR of generators in determining the optimal SR amount for each time period.In reserve-constrained UC, these SR requirements are then set as constraints.
However, based on a time-decoupled UC problem, the preprocessing may result in a suboptimal solution.
The aforementioned conservative methods [14][15][16][17][18] can be improved by using a multiscenario risk analysis method [20] that is capable of taking into account the probabilistic distribution of load and wind power forecast errors, as well as FOR of generators.However, wind power spillage and the nonnormal probability distributions of load and wind power forecast errors cannot be taken into account in the system EENS calculation in [20].This may affect the optimal determination of SR requirements.
Therefore, based on [20] and by combining costbenefit analysis with an improved multiscenario risk analysis method, this paper proposes a novel security-constrained unit commitment (SCUC) model to determine the optimal SR amount.Specifically, in order to consider wind power spillage and the nonnormal probability distribution of wind power forecast error in system EENS calculation, the single discretization of net demand forecast error in [20] is extended into the respective discretization of wind power and load forecast errors.After the EENS under each scenario is calculated, the system EENS can be determined by probability-weighted aggregation.The proposed UC model can be formulated as mixed integer linear programming (MILP) and solved with a commercial solver like CPLEX to improve the computational efficiency and robustness.
The contributions and added value of this paper compared with the respective cited methods are threefold: (1) an improved multiscenario risk analysis method capable of considering various uncertainties including load and wind power forecast errors as well as forced outages of generators is proposed to optimize the spinning reserve requirements; (2) compared with the respective cited methods including the deterministic methods [5][6][7][8][9] and probabilistic methods [10][11][12][13] as well as the cost-benefit analysis [14][15][16][17][18][19][20], the proposed method can consider the impacts of various uncertainties including the nonnormal probability distributions and curtailments of wind power on the overall economic efficiency of system operation via combining cost-benefit analysis with the proposed improved multiscenario risk analysis method so that the optimal economic values of wind power and SR are achieved; (3) the proposed UC model is formulated as a MILP which can be solved effectively via commercial solvers.However, it should be noted that the practical limitation of the proposed method is its computational efficiency.Because of the huge number of decision variables caused by various uncertainties including forecast errors of load and wind power as well as FOR of generators, the simulation times of [20] and the proposed method are larger than the ones of the deterministic methods and probabilistic methods.Hence, the method in [20] and the proposed method are not suitable for online applications of a large-scale power system.
The remainder of this paper is laid out as follows.First of all, Section 2 presents the procedure of improved EENS calculation method with consideration of wind power spillage and the nonnormal probability distribution of wind power forecast error.Then, the SCUC model is proposed in Section 3.After that, in Section 4, the case studies and simulation results analysis are presented.At last, in Section 5, main conclusions are summarized.

Formulation of EENS
In this section, based on multistep modeling of net demand, an improved multiscenario risk analysis method capable of taking into account different probability distributions of the load and wind power forecast errors, as well as FOR of generators, is proposed to determine the system EENS with consideration of wind power spillage.The formulation of EENS is later transformed into an equivalent MILP form.

Improved EENS Calculation Method with Consideration
of Wind Power Spillage.The system EENS is caused by the available SR not being able to meet the excess of load and wind power forecast errors plus shortfalls in the power output of online generators [20].In this paper, the unit of EENS is MWh.
In [20], the generator outages are a series of binary variables, whereas the net demand forecast error is a continuous variable.In general, the EENS formulation can be divided into three steps.
Step 2. The normal distribution of net demand forecast error is discretized into NL (NL = 7) intervals and then combined with the contingency events in each scenario.
Step 3. By summing the expectations of NL intervals resulting in some loss of load, the EENS in a scenario is determined.The total system EENS is the summation of all EENS of each scenario weighted by probabilities of corresponding scenarios.
Unlike [20], wind power spillage and the nonnormal distribution of load and wind power forecast errors are considered in the system EENS calculation of the proposed method.The method in [20] is unable to combine wind power spillage with net demand forecast error and deal with different probability distributions.Consequently, the single discretization of net demand forecast error in Step 2 of [20] needs to be extended into the respective discretization of wind power and load forecast errors.Therefore, the three steps in [20] are extended into four steps in the proposed method as shown in Figure 1.
In Step 1, only single-outage random events are considered, because multiple-outage random events have relatively small probabilities in the scheduling horizon while consuming far more computational resources.In addition, the probability distribution of wind power and load forecast errors in Steps 2 and 3 can be divided into smaller intervals.Though Step 1 Contingency events Step 2 Discretization of net demand forecast error Step 1 Contingency events Step 2 Discretization of wind power forecast error considering wind power spillage Step 3

Discretization of load forecast error
Method of [20] e proposed method Step 3 Total system EENS calculation Step 4 Total system EENS calculation Figure 1: Comparison of EENS calculation process between the method of [20] and the proposed method.
higher accuracy of the result can be achieved with small intervals, more computational resources will be required.

Procedure of Improved EENS Calculation Method with Consideration of Wind Power Spillage and the Nonnormal Probability Distribution.
The improved EENS calculation process with consideration of wind power spillage and the nonnormal probability distribution can be described as follows.
Step 1.Note that the random outage events among NI generators are independent of each other.The probability    of all scheduled generators available except generator  is NI + 1 fault scenarios can be constructed according to no contingency event and single-order contingency events.During period t, taking scenario  as an example, the deficient or redundant SR amount    of other generators in accordance with the outage generator  under scenario  can be calculated as Note that the scenario of no contingency event is the base scenario, when  = 0.
Step 2. Wind power forecast error is considered to follow a nonnormal distribution in [21,22], similar to load forecast error.But, for the convenience of description and without loss of generality, the normal distribution is still chosen as an example in this paper.The normal distribution of wind Seven-interval discretization of normal distribution of wind power forecast error.
power forecast error is approximately discretized into NL (NL = 7) intervals.Figure 2 shows typical seven-interval discretization of wind power forecast error modeled as a continuous normally distributed random variable with zero mean    and a standard deviation    .The probability of interval  1 is   1 , and the width of each interval is    .For each interval  1 ( 1 = 1, 2, . . ., 7), the midvalue represents the value of corresponding whole interval.Figure 2 shows typical seven-interval discretization of continuous normal distribution of wind power forecast error.
Combining the realization of uncertainties of generators in fault scenarios with NL intervals of wind power forecast error, (NI + 1) × NL new scenarios can be constructed.Then, after adding   1 to    , the system forecast error   , 1 caused by the uncertainties of generators and wind power forecast error can be formulated as In (3), a binary variable    1 is introduced to judge whether    has effects on the value of wind power forecast error ( 1 −(NL+ 1)/2)   ( 1 ≤ 4).When  1 ≥ 5, wind power spillage has no effect on the overestimation of wind power output.The binary variable    1 ( 1 ≤ 4) satisfies The nonlinear conditional expression in (4) can be equivalent to the following linear inequalities: Step 3. Load forecast error is similarly modeled as a normally distributed random variable with zero mean    and a standard deviation    .The probability distribution of load forecast error can be similarly discretized into NL intervals, and the probability of interval  2 is denoted as   2 .In each interval  2 ,   , 1 is subtracted from ( 2 − (NL + 1)/2)   to calculate the total system forecast error.The intervals only need to satisfy In order to differentiate the intervals satisfying ( 6) from others, it is necessary to introduce another new binary variable   , 1 , 2 which satisfies The nonlinear formulation above can be linearized to Step 4.During period t, the EENS in each scenario constructed in Step 2 can be determined by summing the expectations of all intervals causing some loss of load as expressed in Then, during period , the total system EENS is the summation of EENS  , 1 weighted by probabilities of corresponding scenarios It is worth noting that when the assumptive normal distribution of wind power forecast error is not correct, the aforementioned EENS calculation method will still be valid.By discretizing the new probability distribution of wind power forecast error into several intervals in Step 2,   1 in (3), (4), and (5) and   1 in (10) can be replaced by the values of wind power forecast error and probability in each interval, respectively.In the same way, the load forecast error which follows the nonnormal distribution can be similarly treated.

EENS Linearization.
In (10), EENS  are formulated as the sum of products of a continuous variable and three binary variables.According to the linearization method in [23], EENS  can be transformed into a series of linear inequalities.
The linearization process can be described as follows. Step Step 2.
, 1 , 2 , and   , 1 , 2 can be formulated as where   , is a nonlinear formulation consisting of the product of a bounded continuous variable and a binary variable, which can be equivalent to the following linear constraints:

Problem Formulation of SCUC Based on Cost-Benefit Analysis
In this section, the mathematical formulation of the SCUC problem is presented.The expected cost of load shedding (ECLS) is expressed as the approximation of EENS given in Section 2 multiplied by VOLL and then added to the objective function of the proposed SCUC model so that the optimal wind power spillage and the optimal quantification of SR are achieved.

Objective Function.
The objective function of conventional UC usually only considers the generation cost (GC) of conventional generators over all scheduling periods.This generation cost contains the fuel cost, start-up cost, and SR cost of generators.However, in the proposed formulation based on cost-benefit analysis, the SR requirements are determined by minimizing the sum of generation cost of generators and ECLS.In addition, though wind power spillage results in additional fuel cost and carbon emission cost (CC) of conventional generators, the penalty cost (PC) of wind power spillage is still considered in the objective function from the wind power's priority in scheduling policy and other environmental benefits of it.Thus, the objective function of the proposed SCUC is defined as where VOLL is the generator outage loss value for users, that is, the loss for outage of 1 MWh electricity, usually estimated via consumer survey [24].

Carbon Emission Cost of Conventional Generators.
In this paper, it is supposed that Clean Development Mechanism (CDM) in Kyoto Protocol [25] is the carbon trading mechanism employed.The carbon emission cost during period  CC  can be calculated by The calculation of carbon emission cost ( 18) can be divided into two cases: (i) When    >    , the carbon emission of power systems exceeds the credit and    −    is the excessive part.Supplementary carbon emission allowance obtained by CDM is then required, thereby increasing the cost for system operation.
(ii) When    ≤    , the carbon emission of power systems is lower than the credit; that is, there is surplus which can be sold through emission transaction and the system operation cost is then decreased equivalently.

The total carbon emission 𝑀 𝑡
during period  is commonly calculated by Because ( 19) is too complex for calculation, it can be replaced by the following equation given in [25]: where the unit of  , is tCO 2 /tce where tce refers to the ton of standard coal equivalent.The carbon emission credit    during period  and the total carbon emission credit   of power systems are determined by the load of each period; that is, where    satisfies 0 ≤    ≤    .

SR Constraints of Generators.
Here, one has where  10    is the maximum 10-minute SR supplied by generator  [17].the early complete enumeration, priority list method, and dynamic programming to the current branch-and-bound method and intelligent optimization algorithm in [8,17].In recent years, the commercial software CPLEX has utilized branch-cutting on the basis of branch-and-bound method and combined various techniques, such as the heuristic method and cutting-plane algorithm, to solve the MILP problem.So far, CPLEX has been widely used to solve UC problems in [5,7,[14][15][16][18][19][20] by approximately linearizing the UC problem to a MILP model.Therefore, this paper also adopts this solver for the solution of the proposed SCUC problem where the operation cost in objective function given by ( 15) and constraints ( 23)-( 29) are linearized in accordance with the method proposed in [18] which is presented in the Appendix.

Upper and the Lower Limits for the Power Output of Generators. One has
The deterministic methods and probabilistic methods are incomplete without the inclusion of a discussion on the impact of FOR.From a computational point of view, the addition of generator contingencies would render the solution process of the UC problem of [20] and the proposed UC problem much more challenging.The addition of a single generator contingency, also considering its possible times of failure, would multiply the number of scenarios by the number of periods of the scheduling horizon.It can be recalled that, with each extra scenario, there are corresponding extra variables and constraints.As a result, realistically sized problems such as the European interconnected power system may be very hard to handle with current computing tools in a reasonable amount of time required for day-ahead UC.However, some modeling simplifications could be considered.One possibility is to make use of the scenario reduction techniques [26].Secondly, decomposition techniques [27] are promising because they exploit the intrinsic decomposable structure of the problem whereby each scenario is optimized individually under the command of a master coordinating problem.These aspects need to be further investigated in subsequent researches.

Case Studies and Simulation Results Analysis
The proposed method for determining SR in this section is tested and verified on a modified IEEE 26-generator reliability test system.Simulation studies are carried out to draw comparisons of UC results between the proposed method and method of [20] to verify the effectiveness and advantage of the proposed SCUC model.

Modified IEEE 26-Generator Reliability Test System.
A modified IEEE 26-generator reliability test system is adopted in this paper to verify the effectiveness of the proposed model.In this system, there are 26 thermal generators with a total capacity of 3105 MW.The hydro generating units have been omitted.The transmission limits, failure rates, ramp rates, cost coefficients, and minimum up and down time are obtained from [18,28].The power output of the generators committed at  = 0 is obtained by the economic dispatch of the committed generators at a load level of  1 lists the forecast load and power output of wind farm used for testing.It can be observed from the table that the wind power output and load are negatively correlated.The standard deviations of load and wind power forecast error can be approximately calculated according to [18].VOLL is set at 4000 $/MWh [20].  is set at 10 $/MWh.The CDM transaction price  CDM is 50 $/tCO 2 ,  , is 3 tCO 2 /tce, and   is 100 $/tce [25].
Considering the effects of maximizing carbon emission cost on SR optimization, the carbon emission credit    is set at 0.
The proposed model is solved on a 2.66 GHz Windowsbased workstation with 8 GB of RAM using a MILP solver CPLEX 12.5 under MATLAB environment.When the duality gap tolerance for CPLEX solver is set at 0.5%, the computation time is about 11.2 s.

Comparison of the Proposed Method and
the Method of [20]

Necessity Analysis with Consideration of Wind Power
Spillage.Wind power has been given priority in scheduling.However, the development of Chinese power grid in regions with abundant wind energy resources (Xinjiang, Ningxia,  Gansu, etc.) is far behind the increase of installed wind power capacity at present, and thereby the large-scale wind power accommodation capacity of power systems is insufficient.Wind power spillage must be adopted to ensure the security and stability of systems operation.
When the capacity limit 500 MW of the transmission lines between buses 14 and 11 and 14 and 16 is modified to 250 MW in the test system, it can be seen from Figure 3 that the transmission flow capacity is insufficient and thereby the wind power cannot be completely accommodated in the proposed method.At this time, the network security constraints in [20] which has no solution are not satisfied.This verifies the necessity of utilizing multistep modeling of net demand to improve the traditional multiscenario risk analysis method of SR optimization.

Economic Analysis with Consideration of Wind Power
Spillage.Compared with the method presented in [20], wind power spillage is taken into account.Table 2 compares different costs (fuel cost, start-up cost, SR cost, ECLS, penalty cost, carbon emission cost, and total cost) between the proposed method and the method of [20] over the whole time periods.It can be observed from this table that although the fuel cost of conventional generators, the penalty cost, and the carbon emission cost increase when wind power spillage is considered, reductions in ECLS, SR cost, start-up cost, and total cost are achieved.Therefore, the optimal determination

Wind power
The proposed method Method of [20]   of SR requirements and the optimal economic value of wind power are realized.This verifies the advantage of the proposed model.The lower total cost of the proposed method can be explained in detail as in Table 2.
Firstly, Figure 3 shows the variation of wind power spillage during all periods in the proposed method.Figures 4 and 5 compare the SR and EENS, respectively, between the proposed method and the method of [20].
(i) When the load level is low and wind power output is large during periods 1-5, 16, and 17, the amount of wind power spillage is not equal to zero.This is because the cost of providing additional SR is higher than the sum of increased fuel cost of conventional generators, the penalty cost, and the carbon emission cost caused by wind power spillage during these periods.Consequently, by increasing a small amount of wind power spillage, the amount of SR is reduced and the system EENS is reduced simultaneously.
Take periods 3 and 16 as an example.On the one hand, this method curtails wind power values 11.87 MW and 24.89 MW, respectively, which leads to the penalty cost increase by 118.70$ and 248.90$, respectively, and the carbon emission cost increase by 256.40$ and 575.70$, respectively, as well as the fuel cost of conventional generators increase by 152.62$ and 353.37$, respectively.On the other hand, the SR cost decreases by 539.23$ and 1109.52$,respectively, and ECLS decreases by 30.06$ and 55.37$, respectively.In summary, during period 3, the total scheduling cost can reduce by 41.57$.During period 16, the amount of wind power spillage is larger, and the sum of these costs increases by 13.08$ caused by wind power spillage.However, Figure 6 compares the generation schedules between the two methods.It can be seen from this figure that when the system load level is low during periods 16 and 17, the outages of generators 13 and 15 are avoided in the proposed method due to additional wind power curtailments.Consequently, the start-up cost of generators during the next period can be reduced (see the comparison of start-up cost shown in Table 2) and the overall economic efficiency of system operation will be improved.
(ii) When the load level is high and wind power output is small during other periods, more expensive generators are committed.As a result, the sum of increased fuel cost and the penalty cost caused by wind power spillage is higher than the cost of providing additional SR.Hence, the wind power will be completely accommodated during high load levels.

Simulation Analysis considering Different Probability Distributions of Wind Power Forecast
Error.The study in [21] shows that the tail of the actual day-ahead wind power forecast error data curve is situated between the Laplace distribution and the normal distribution curve.Kurtosis is chosen as the statistical parameter to evaluate the tail of the studied probability density function, and the kurtosis of the actual wind power forecast error ranges from 3 to 6. Therefore, in this paper, the kurtosis value of the new wind power forecast error is set at 4.8, which has the same mean and standard deviation as the normal distribution in [18].This new probability density function is  similarly divided into seven probability intervals, and the median, range, and probability of each interval are shown in Table 3.
1 in (3), (4), and (5) and   1 in (10) are taken from the median and probability values of each wind power forecast error interval in Table 3, respectively.Figure 7 compares the SR capacity between the normal distribution and the probability distribution with kurtosis 4.8.It can be seen from this figure that the tail region of forecast error will affect the determination of SR and using the normal distribution to model the wind power forecast error will lead to an increase of the probability of SR insufficiency.

Conclusions
Combining cost-benefit analysis with an improved multiscenario risk analysis method capable of taking into account the load and wind power forecast errors as well as FOR of generators, a novel SCUC model with consideration of wind power spillage is proposed in this paper to determine the optimal amount of SR.
In order to consider wind power spillage and the nonnormal distribution of wind power forecast error in system EENS calculation, based on multistep modeling of net demand, the single discretization of net demand forecast error in [20] is extended into the respective discretization of wind power and load forecast errors.From the results of simulation studies, in cases where load and significant wind power

Figure 3 :
Figure 3: Wind power curtailments with different capacity limits of transmission lines.

Figure 4 :
Figure4: Comparison of wind power curtailments and SR between the proposed method and the method in[20].

Figure 5 :
Figure5: Comparison of EENS between the proposed method and the method of[20].

Figure 6 :
Figure6: Comparison of the generation schedules between the proposed method and the method of[20].

Figure 7 :
Figure 7: Comparison of SR between different probability distributions.
1.A binary variable   , 1 , 2 and a continuous variable   , 1 are introduced.Let   , 1 , 2 be the product of the binary variable   , 1 , 2 and the binary variable    , and let   , 1 be the product of the bounded continuous variable    and the binary variable    1 .Then, EENS  will be transformed into the summation of products of a binary variable and a continuous variable.The above new variables   , 1 , 2 and   , 1 can be equivalent to the following linear constraints: (   ,    ) is the fuel cost of generator  during period  which is described by the consumption characteristic curve of a quadratic function;  su,     (1 −  −1  ) is the start-up cost of generator  during period ;  ,     is the reserve cost of generator  to supply SR amount    during period t.

Table 1 :
Forecast data of load and power output of wind farm.

Table 2 :
[20]arison of different costs between the proposed method and the method presented in[20].

Table 3 :
Parameters of seven-interval approximation for probability distribution with kurtosis 4.8.
na n ds t a n d a r dd e v i a t i o no fl o a d forecast error   1 ,   1 : V a l u ea n dp r o b a b i l i t yo fw i n dp o w e r forecast error in interval  1   2 : Probability of load forecast error in interval  2    : Probability of all scheduled generators available except generator     ,    : Ramp-up and ramp-down rates of generator  in MW/h  on ,min ,  off ,min : Minimum up and down time of generator   on ,−1 ,  off ,−1 : Up and down time which have been accumulated up to period  − 1   , : Power output of conventional generators connected to bus  during period    , ,   sw, : P o w e ro u t p u to fw i n df a r mc o n n e c t e d to bus  and the corresponding amount of wind power spillage during period    , : L o a d a t b u s  during period    :   ,   ,   and   : Carbon emission coefficients of thermal generator   , :