A Bilateral Tradeo ff Decision Model for Wind Power Utilization with Extensive Load Scheduling

In this paper, we present the extensive load scheduling problem with intermittent and uncertain wind power availability. A chance-constrained bilateral tradeoff decision model is established to solve the problem. Our model aims at maximizing the wind power utilization and minimizing the system operation cost simultaneously by means of responsive loads, which are precisely divided into shiftable loads and high-energy loads. The chance constraint is applied to restrict the system imbalance with a small probability. Then, a revised sample average approximation (SAA) algorithm is developed to transform the chance constraint into sample average reformulations. Furthermore, the multi-objective differential evolution (MODE) method combined with SAA is proposed to solve the problem. Experiments enabling an effectiveness analysis of the two kinds of responsive loads are performed on the power system in Yancheng. The research of parameters of MODE, the sensitivity of different risk levels and the influence of iteration numbers are discussed. Finally, computational results prove that the combination of shiftable loads and high-energy loads have a better effect than adopting shiftable loads and high-energy loads separately, and the proposed method is convergent and valid in solving the problem.


Introduction
Wind power, which is clean and renewable, has been integrated with the grid on a large scale to replace fossil fuels [1,2].A lot of countries have set goals for high-penetration wind power levels in the grid to increase the wind power utilization.China aims to achieve an average wind power utilization rate of 95% in 2020.However, the randomness and intermittence of wind power have increased the tremendous challenges in wind energy utilization as well as the power system operators and planners [3,4].The traditional dispatching operating mode of a power system, designed to address limited uncertainty in the system, encounters challenges in accommodating the high penetration of large-scale wind power.Thus, new methods should be developed to cope with the situation [5].
The chance-constrained optimization technique has been adopted to manage the problems involved stochastic variables in some previous studies [6,7].In [7], the chance constraints in the optimization formulations guarantees that the failure probability of the EV charging plan fulfilling the driving requirement is below the predetermined confidence parameter.In [6], a chance constrained programing approach was presented to ensure desired confidence levels of meeting future stochastic power and natural gas demands while minimizing the investment cost.In the present study, the chance constraint is applied to ensure the probability of system imbalance caused by too little system reserve capacity is less than a certain probability.These types of constraints could not be inverted to obtain equivalent deterministic equalities.Thus, the method developed in [7] could not be directly applied here to solve our problem.In addition, the joint chance constrained problem in [6] are formulated through mixed-integer quadratic programming (MIQP), which is not efficient method to solve our problem related to not joint chance constraints.Therefore, we propose the sample average approximation (SAA) algorithm [8] to solve the problem.The approach can provide a solution that converges to the optimal one as the number of samples increases.
During the Spring Festival, the load is expected to drop sharply which will decrease the usage of power.In 2018, to ensure the safe and stable operation of the power grid, Shandong Electric Power Company took measures to limit renewable energy output, which means that a large amount of wind power would be abandoned.However, Jiangsu Electric Power Company actively encourages enterprises to produce, which not only increases the load of 9000 MW for trough load periods but also promotes the consumption of renewable energy.In the first three days of the Spring Festival, a total of 72000 MWh of renewable energy consumption was added.In the smart grid, the load, as an important responsive resource [9,10], is an excellent alternative for wind power consumption [11,12].
Extensive research has been done to increase wind power utilization through the participation of responsive loads.In [13], a model for solving the combined effect of PEVs and wind power integration with a DR (demand response) program on static transmission network expansion planning was proposed, but the uncertainty of wind power has not been considered.A two-stage stochastic model which handling the forecast errors of renewable generations in a microgrid with responsive loads was implemented in [14], however, the objective of the model is only operation cost not including the renewable energy utilization.Hajibandeh Neda [15] investigated the potential of DR as an emerging alternative in systems with significant amounts of wind power.A comprehensive set of DR programs including tariff-based, incentive-based and combinational DR programs are considered in a stochastic network-constrained market clearing framework.Kalavani Farshad [16] provided a stochastic method to conduct the optimal scheduling of the combination of wind power and energy storage with considering the DR program in the electricity market.The DR program was adopted to increase the total expected profit and decrease the total operational cost.
Most models in this context have optimized wind power utilization by responsive loads.However, the scheduling characteristics of loads were not considered, and responsive loads were not divided finely.Different types of responsive loads have different characteristics [17].Therefore, different loads can participate in the optimized operation in different ways [18], which may lead to better optimized results.
In this study, according to their regulating characteristics, loads are grouped into two types, including the high-peneration energy load (HL) and the shiftable load (SL).The HL possesses both adjustable and interruptible characteristics, which can achieve power control in the range of 0%-100%.The regulate action can be completed in an instant.The HL is suitable for increasing wind power consumption because we can raise its power when the wind power is high and reduce its power in the opposite case.The use of an SL implies shifting the load from one period to another to maintain the electricity consumption.Loads can be cut out or put in at any periods of the whole cycle.Generally, an SL is put in during the peak wind power periods and cut out in trough wind power periods to increase the wind power utilization.HLs and SLs are shown in Figure 1.The two types of responsive load participate in the optimal scheduling with conventional power, developing a new scheduling mode based on source-load coordinate optimal operation.However, the dispatch of responsive loads will damage the users' benefit.Economic compensation needs to be made for them, which incurs additional costs.Thus, besides the amount of consumption of wind power, the cost of system operation should be considered as well.The model contains two objectives: to obtain the maximum capacity of wind power accommodation and to minimize the system operation costs.A multi-objective optimal model is established in this study.
There are many algorithms that can be used to solve the multi-objective optimal problem, for example: the Non-dominated sorting genetic algorithm II (NSGA-II) [19], the Strength pareto evolutionary algorithm (SPEA) [20], the Pareto archived evolution strategy (PAES) [21] and the Multiobjective differential evolution (MODE) [22].In [19], NSGA-II, SPEA2 and PAES were designed for multi-objective problems, and NSGA-II achieved the best results in terms of Pareto front spread.In [22], NSGA-II, MODE and multi-objective particle swarm optimization (MOPSO) algorithms are compared in accuracy and computational time.The results show that the obtained Pareto Front of MODE is more accurate and faster.It can be concluded that MODE performs better among these algorithms.So we adopted MODE algorithm to solve the multi-objective model in this paper.
The remaining part of this paper is organized as follows.Section 2 describes the problem with chance constraints.Section 3 presents the solution methods which including SAA and MODE algorithms.Section 4 reports the computational experiments for the power system of Yancheng City of Jiangsu province in China.Section 5 concludes the study.

Problem Description
In this section, we develop a chance-constrained multi-objective formulation to address uncertain wind power availability.The two objectives are to maximize the utilization of wind power W E and to minimize the system operation cost GL C , which including the costs of thermal generator and the operation costs of SLs and HLs.HLs and SLs are adopted in the optimization.The formulation is to determine the regulate or shift operations of the responsive loads.The detailed formulation is shown as follows: . .s t  The two types of responsive load participate in the optimal scheduling with conventional power, developing a new scheduling mode based on source-load coordinate optimal operation.However, the dispatch of responsive loads will damage the users' benefit.Economic compensation needs to be made for them, which incurs additional costs.Thus, besides the amount of consumption of wind power, the cost of system operation should be considered as well.The model contains two objectives: to obtain the maximum capacity of wind power accommodation and to minimize the system operation costs.A multi-objective optimal model is established in this study.
There are many algorithms that can be used to solve the multi-objective optimal problem, for example: the Non-dominated sorting genetic algorithm II (NSGA-II) [19], the Strength pareto evolutionary algorithm (SPEA) [20], the Pareto archived evolution strategy (PAES) [21] and the Multi-objective differential evolution (MODE) [22].In [19], NSGA-II, SPEA2 and PAES were designed for multi-objective problems, and NSGA-II achieved the best results in terms of Pareto front spread.In [22], NSGA-II, MODE and multi-objective particle swarm optimization (MOPSO) algorithms are compared in accuracy and computational time.The results show that the obtained Pareto Front of MODE is more accurate and faster.It can be concluded that MODE performs better among these algorithms.So we adopted MODE algorithm to solve the multi-objective model in this paper.
The remaining part of this paper is organized as follows.Section 2 describes the problem with chance constraints.Section 3 presents the solution methods which including SAA and MODE algorithms.Section 4 reports the computational experiments for the power system of Yancheng City of Jiangsu province in China.Section 5 concludes the study.

Problem Description
In this section, we develop a chance-constrained multi-objective formulation to address uncertain wind power availability.The two objectives are to maximize the utilization of wind power E W and to minimize the system operation cost C GL , which including the costs of thermal generator and the operation costs of SLs and HLs.HLs and SLs are adopted in the optimization.The formulation is to determine the regulate or shift operations of the responsive loads.The detailed formulation is shown as follows: where In the above formulation, the functions (1) represents the two objectives.One is to maximize utilization of wind power; the other is to minimize the system operation cost, which is composed of the generation cost of thermal generators C gen , the operation cost of shiftable loads C SL and the operation cost of high-energy loads C HL .Constraint (2) represents the power balance of the system.P t L is the total amount of conventional system load over time period t which doesn't include SL and HL.Constraint (3) describes that the reserve needed for loads forecast error and wind power forecast error should be satisfied.Constraint (4) limits the upper and lower bounds of power generation of each unit.Constraints ( 5) and ( 6) mean the status (on or off) of each unit must last for a minimum time once it is started up or shut down.Constraints (7) and ( 8) are the unit ramping up constraints and ramping down constraints respectively.Constraint (9) represents that the utilization of wind power cannot exceed the maximum value.Chance-constraint (10) means that the probability of system imbalance caused by the too little system reserve capacity is less than risk level ε.Constraints ( 11) and ( 13) indicate the SLs and HLs must be less than or equal to the schedulable capacity.Constraint (12) show that the total power consumption of each shiftable load can't change.Constraints ( 14) and ( 15) describe the use of SLs and HLs should exceed a certain minimum time.Constraints ( 16) and ( 17) describe the limitation of the dispatched numbers of each SL and HL.Constraint (18) represents the transmission capacity constraints of each line.Constraint (19) is the formula of DC power flow calculation.
Equalities ( 20), ( 21) and ( 22) describe the costs of thermal generator and the operation costs of SLs and HLs.
To sum up, the bilateral tradeoff model can be described as (23).F = (−E W , C GL ) is the objective function for this study.x is a decision variable.f j (x) is the equality constraint function, and g k (x) is the inequality constraint function.

Revised Sample Average Approximation
SAA is an effective method to solve the chance-constrained problem.The basic idea of SAA is to approximate the true distribution of random variables with an empirical distribution through Monte Carlo sampling technology.Multiple theoretical research studies and computational studies of SAA have been developed for chance-constrained stochastic problems [23].In this section, we develop a combined SAA algorithm to solve the chance-constrained program.The combined SAA framework contains three parts: scenario generation, convergence analysis and solution validation.

Scenario Generation
The wind power availability is closely related to the wind speed.Recently, a variety of theoretical distributions were proposed to fit the probability distribution of wind speed [24].In the present study, we adopt two-parameter Weibull distribution to fit the probability distribution of actual wind speed.The Monte Carlo simulation is run to generate N scenarios of the wind power random variable ζ t,n wi , and each scenario has the same probability of 1/N. 1 (0,∞) (•) is an indicator function, i.e., The chance constraints (10) can be replaced by First, if δ = ε, as the sample size N goes to infinity, we can prove that the solution of the SAA problem converges to the solution of the true problem.Since N cannot be infinity, it is possible that the solution of the SAA problem is infeasible with respect to the true problem.To increase the confidence level for the solution of SAA problem to be feasible, we let δ < ε.Generally, and N = 2/ε, δ = ε/2.

Solution Validation
Assume that x is the solution of the SAA problem and that v is the corresponding value of objective function.To test the feasibility of x, we construct the upper bound of v.
First, S iterations were taken.For each iteration, we run the SAA problem M times with N scenarios.As described in [23], we let L be the largest integer to meet the constraint (27).
With B denoting the binomial distribution.Then, selecting the Lth smallest optimal value denoted as vLs , the confidence level for vLs to be the upper bound of the chance-constrained problem is more than 1 − ω.The average of vL1 , vL2 , . . ., vLS is the upper bound; the smallest one of vL1 , vL2 , . . ., vLS is the optimal value v o .

Summary of the SAA Algorithm
The proposed SAA algorithm is summarized in the following steps. (1 scenarios of the random vector ζ; ii.Solve the proposed model by MODE to obtain the solution as x and the objective value as v. 2 Sort the M optimal values in an ascending order, and denote the Lth value vL as vLs . (2) The upper bound is defined as v = 1 S S s=1 vLs , and the minimum value from vL1 to vLs is taken as the optimal objective value v o .The corresponding optimal solution x can also be obtained.(3) Calculate the optimality gap by

Multi-Objective Differential Evolution
The problem is a multi-objective mixed-integer linear program (MILP) [25], which cannot be solved by CPLEX.CPLEX can provide flexible high-performance optimization programs to solve the single-objective problems, which including linear programming problem, quadratic programming problem and mixed-integer programming problem.However, it is not applicable in multi-objective problems.Hence, MODE algorithm was proposed in this paper.Because of the dual objectives, a set of Pareto optimal solutions would be obtained.Therefore, a novel fuzzy membership function [26] is introduced to select the optimal compromise solution from Pareto solution set.

Multi-Objective Differential Evolution
MODE is a type of random parallel direct global search algorithm with the advantages of simple and easy-to-use global optimization.Its details are as follows: (1) Population initialization: The initial population X 0 is chosen arbitrarily in the solution space, and the population scale number is N p .(2) Mutation operation: The variation operation is based on the difference of random sampling individuals as follows: where G is the evolution times, v i,G+1 is a variation individual and x r1,G is the father individual.
(3) Crossover operation: Once the perturbed individual v i,G+1 is generated, this individual is subjected to crossover operation with a target individual X i,G+1 that ultimately generates the trial solution,u i,G+1 , as follows: (4) Selection operation: We make a temporary population X * and then let the test individual z i,G+1 compete with the father individual x i,G .The better one would be selected to X * [27].

Selection of the Optimal Compromise Solution
To obtain the balance of wind power cost and economical operation of the system, we need to measure the membership of each Pareto solution and use this measurement to find the optimal compromise solution.In this study, we introduce a novel fuzzy membership function as follows: where f imax and f imin are the upper and lower bound of the ith objective function, respectively.The standard membership degree is estimated as u = 2 i=1 u i .The largest membership is corresponding to the optimal compromise solution.The flowchart of the SAA algorithm combined with MODE is shown in Figure 2.

Results and Discussion
In this section, we conduct a series of experiments on the power system of Yancheng in Jiangsu, China to illustrate the proposed approach.The advantages of dual-objective model and the impact of SLs and HLs are manifested.We also discuss the sensitivity of risk level and the convergence of the proposed SAA algorithm combined with MODE.The time horizon was 24 hours at 0.5 hours step.The relevant information can be obtained from the state grid corporation.

Introduction of the Yancheng Power System
Yancheng, with six wind farms, has abundant wind power resources.The topological structure of the Yancheng power system is shown in Figure 3.The specific information of HLs and SLs is shown in the Table 1 and Table 2.

Results and Discussion
In this section, we conduct a series of experiments on the power system of Yancheng in Jiangsu, China to illustrate the proposed approach.The advantages of dual-objective model and the impact of SLs and HLs are manifested.We also discuss the sensitivity of risk level and the convergence of the proposed SAA algorithm combined with MODE.The time horizon was 24 h at 0.5 h step.The relevant information can be obtained from the state grid corporation.

Introduction of the Yancheng Power System
Yancheng, with six wind farms, has abundant wind power resources.The topological structure of the Yancheng power system is shown in Figure 3.The specific information of HLs and SLs is shown in the Tables 1 and 2. M Hb and M Sa are the maximum schedule times for HLs b and SLs a, respectively.The reserve capacity of the thermal generators needed for the load forecast error is 10% of the maximum system loads.The reserve of thermal units needed for the fluctuation of wind power is 15% of the wind power generation.The maximum ramping speed of thermal units is 1% of the rated capacity every minute.Other information for the thermal generators is shown in Table 3.The iteration numbers of SAA are kept as S × M = 5 × 5 and risk level is ε = 0.2.The population size of the MODE is 100, the maximum iteration number is 500, the crossover probability is 0.3, the mutation probability is 0.5.In the study, we adopt two-parameter Weibull distribution to fit the probability distribution of actual wind speed.All the experiments are conducted on a computer workstation with 4

Research for Parameters of MODE
There are some important parameters of MODE like population size, iteration number, crossover probability, and mutation probability.In order to study the impacts of population size and iteration number on the execution time, they were set as different values.The results are shown in Figure 4.

Research for Parameters of MODE
There are some important parameters of MODE like population size, iteration number, crossover probability, and mutation probability.In order to study the impacts of population size and iteration number on the execution time, they were set as different values.The results are shown in Figure 4. From Figure 4, it is apparent that population size and iteration number have great influences on the execution time.With the increase of population size, the execution time increased lightly at first stage, then there was a significant decline, afterwards it increased again.When the population size becomes larger, on the one hand, the calculation speed slows down; on the other hand, the global search ability of the algorithms becomes stronger, the optimal solution can be found in less generations.When the population size is 100, the execution time is least.As for the iteration number, the execution time became larger when it increases.Nevertheless, the smaller the iteration number, the less accurate the results.In order to take both execution time and accuracy into account, the iteration number was set as 500.
To verify the above analysis of population size and iteration number, we performed 10 runs for each value of population size and iteration number to count the run time.Besides, the hypervolume of each Pareto solution set was also calculated.The reference point is (41,739, 3281,867.2).The results are shown in Table 4 and Table 5.From Figure 4, it is apparent that population size and iteration number have great influences on the execution time.With the increase of population size, the execution time increased lightly at first stage, then there was a significant decline, afterwards it increased again.When the population size becomes larger, on the one hand, the calculation speed slows down; on the other hand, the global search ability of the algorithms becomes stronger, the optimal solution can be found in less generations.When the population size is 100, the execution time is least.As for the iteration number, the execution time became larger when it increases.Nevertheless, the smaller the iteration number, the less accurate the results.In order to take both execution time and accuracy into account, the iteration number was set as 500.
To verify the above analysis of population size and iteration number, we performed 10 runs for each value of population size and iteration number to count the run time.Besides, the hypervolume of each Pareto solution set was also calculated.The reference point is (41,739, 3281,867.2).The results are shown in Tables 4 and 5. Hypervolume can be used to assess the overall performance of Pareto solution set.And the bigger the hypervolume, the better the solution.With the increase of population size or iteration number, the corresponding hypervolume became larger as well.Because it is more likely to find the global optimal solution with larger population sizes or more iterations numbers.From Tables 4 and 5, it can be known that there are some deviations in the execution time and solutions obtained by the same parameters.It is because the initial population of MODE is chosen randomly.
Different crossover probabilities or mutation probabilities correspond to different optimization paths, and different optimization paths correspond to different calculation time or optimal solution.In order to explore the influences of crossover probability and mutation probability, we also performed 10 runs for each value of them.The results are shown in Table 6.It can be seen that the influences of crossover probability and mutation probability on the execution time and solutions are lighter than the influences of population size and iteration number.Although different crossover probability or mutation probability would lead to the different optimization path, the optimization path is also related to the initial population, which is chosen arbitrarily in the solution space.Thus, based on empirical values and Table 6, the crossover probability is set as 0.3 and the mutation probability is set as 0.5.

Effectiveness of Responsive Loads
To research the effect of SLs and HLs, three types of optimization model were established.Model 1 was to optimize the problem with SLs only.Model 2 was to optimize the problem with HLs only, and model 3 employed both SLs and HLs.A conventional model without responsive loads was also built to make a contrast.For the multi-objective optimization problems, the optimal compromise solutions are in Table 7, E W represents the wind energy utilization and C GL is the cost of system operation.Obviously, compared with conventional type, all the three model types not only contribute to the increasement of wind energy utilization but also facilitate the decrease of system operation cost.This finding is reasonable.Despite the dispatch of responsive loads incurs expenses, the increased consumption of wind power will lead to the decrease of thermal units' output.And the cost saved by thermal units are more than the expenses of responsive loads.Thus, it can be concluded that the adoption of responsive loads can bring better effect.According to Table 4, it is not difficult to find that model 3 is the best one for it can lead to most increase of wind power utilization and most decrease of system operation costs.In other words, the effect of combining two types of responsive load is better than using the two types of load separately.

Influence on the Load Curve
The scheduling of extensive loads, which is related with the trend of forecasted wind power availability, will influence the load curve of system.The results are shown in Figure 5.

Influence on the Load Curve
The scheduling of extensive loads, which is related with the trend of forecasted wind power availability, will influence the load curve of system.The results are shown in Figure 5. From Figure 5, it is obvious that the daily peak load is cut, and the daily valley load is increased after optimization.It promotes the daily load curve leveling and is conducive to the system security and stability.This kind of impact is reasonable because wind power owns anti-peak characteristics, which means the trough periods of forecasted wind power are the peak periods of daily load and the peak periods of forecast wind power are the trough periods of daily load.And to maximize the wind power utilization, the scheduling of responsive loads should make the valley load increased and the peck load decreased, leading the daily load curve consistent with the forecasted wind power curve.Contrasting the 3 wind power curves, it is easy to find that the optimal wind power utilization is larger than the conventional wind power utilization during 2:00 to 6:30 and 14:00 to 15:00, which verifies the effect of responsive loads.
To sum up, the combination of HLs and SLs in the optimization model is of great significance, for it can not only lead better effect in maximizing wind power utilization and minimizing system operation cost, but also contribute to leveling the load curve and maintaining the system security.

Influence of Risk Level ε
The availability of wind power is uncertain.Therefore, there is some risk that the system reserve capacity cannot surmount the prediction errors of wind power.The probability of system imbalance caused by too little system reserve capacity is the risk level.With different risk levels, this optimization model produces different results.To test the influence of the risk level, several levels of risk are tried in the model with both SLs and HLs involved in the optimization.The iteration number is 5 5 S M × = × .The results are shown in Table 8.From Figure 5, it is obvious that the daily peak load is cut, and the daily valley load is increased after optimization.It promotes the daily load curve leveling and is conducive to the system security and stability.This kind of impact is reasonable because wind power owns anti-peak characteristics, which means the trough periods of forecasted wind power are the peak periods of daily load and the peak periods of forecast wind power are the trough periods of daily load.And to maximize the wind power utilization, the scheduling of responsive loads should make the valley load increased and the peck load decreased, leading the daily load curve consistent with the forecasted wind power curve.Contrasting the 3 wind power curves, it is easy to find that the optimal wind power utilization is larger than the conventional wind power utilization during 2:00 to 6:30 and 14:00 to 15:00, which verifies the effect of responsive loads.
To sum up, the combination of HLs and SLs in the optimization model is of great significance, for it can not only lead better effect in maximizing wind power utilization and minimizing system operation cost, but also contribute to leveling the load curve and maintaining the system security.

Influence of Risk Level ε
The availability of wind power is uncertain.Therefore, there is some risk that the system reserve capacity cannot surmount the prediction errors of wind power.The probability of system imbalance caused by too little system reserve capacity is the risk level.With different risk levels, this optimization model produces different results.To test the influence of the risk level, several levels of risk are tried in the model with both SLs and HLs involved in the optimization.The iteration number is S × M = 5 × 5.The results are shown in Table 8.As the risk level increases from 0.05 to 0.5, the total wind power utilization increases, and the cost of system operation decreases.This is reasonable because a lower risk level corresponds to a more restrictive chance constraint.Correspondingly, the generation of thermal units is decreased, which results in a decrease in total cost.An extreme case is ε = 1, in which the chance constraint can be neglected.In such a case, the optimal wind power utilization is larger than that at any other risk level; furthermore, the cost of system operation is the smallest.

The Influence of Iteration Numbers
The SAA algorithm is a type of approximate method.The result is not the solution of the true problem.Changes in iteration numbers result in changes in the accuracy of the results.Thus, experiments are done to study the impact of the iteration numbers.In the experiments, SLs and HLs are involved in the optimization, and the risk level is kept as 0.2.The results are shown below.
Tables 9 and 10 show that with S = 5 and M increasing from 5 to 10, the gap between the upper bound and optimal value of wind power utilization decreases from 0.61% to 0.42%, and the gap of the system operation cost decreases from 0.34% to 0.21%.With M = 5 and S increasing from 5 to 10, the gap between the upper bound and optimal value of wind power utilization decreases from 0.61% to 0.31%, and the gap of the system operation cost decreases from 0.34% to 0.17%.The decrease in the gaps indicates that the solution became more precise; with accumulating iterations, the proposed SAA algorithm converges to the true values.What's more, an increase in S can lead to a greater decrease of the gap than an increase in M. In this study, the compensation for the operation of interactive loads is given by the demand response special fund of the Jiangsu power company, which is the main motivation of the schedule of extensive load.The focus of this study is to discuss how to increase the wind power utilization through the refined controllable loads.As for the power market mechanism incorporating renewable energy and interactive loads, we take it as our next researcus.

Conclusions
In this study, a chance-constrained multi-objective program considering the uncertain wind power availability has been studied.SLs and HLs are innovatively adopted to optimize the problem, and a revised SAA algorithm combined with MODE is developed toch fo solve it.The experiments performed on the power system of Yancheng show that the combination of SLs and HLs leads to better results in maximizing the wind power utilization and minimizing the system operation cost; and responsive loads are also helpful in leveling the load curve.The impacts of parameters of MODE on the execution time and hypervolume are taken into account.Additionally, the sensitivity of different risk level and the influence of iteration scales of SAA are discussed to verify the accuracy and validity of the proposed method.Start-up and shut-down cost for thermal unit j n j , m j , l j Operating cost factors for thermal generator j. λ Sa Cost factors for shiftable load b to regulate a unit load λ Hb Cost factors for high-energy load b to regulate a unit load P t L Total amount of system load over time period t P max

Gj
Maximum output of thermal generator j.P min

Gj
Minimum output of thermal generator j.R t L Amount of spinning reserve needed for the system over time period t. γ Reserve capacity factor of wind power R up Gj Ramp-up output limit of generator j.

Figure 1 .
Figure 1.Introduction of an HL and an SL.(a) The scheduling characteristic of an HL.(b) The scheduling characteristic of an SL.

Figure 1 .
Figure 1.Introduction of an HL and an SL.(a) The scheduling characteristic of an HL.(b) The scheduling characteristic of an SL.

Figure 2 .
Figure 2. Proposed SAA algorithm combined with MODE.

M
are the maximum schedule times for HLs b and SLs a,

Figure 2 .
Figure 2. Proposed SAA algorithm combined with MODE.

Figure 4 .
Figure 4. Execution time of different population size and iteration number.

Figure 4 .
Figure 4. Execution time of different population size and iteration number.

Figure 5 .
Figure 5. System load and wind power utilization.

Figure 5 .
Figure 5. System load and wind power utilization.

Table 1 .
Intel Cores and 8 GB RAM.Parameters of high-energy load users.

Table 2 .
Parameters of Shiftable load users.

Table 1 .
Parameters of high-energy load users.

Table 2 .
Parameters of Shiftable load users.

Table 3 .
Parameters of conventional units.

Table 4 .
Run time and hypervolume under different population sizes.

Table 5 .
Run time and hypervolume under different iteration numbers.

Table 4 .
Run time and hypervolume under different population sizes.

Table 5 .
Run time and hypervolume under different iteration numbers.

Table 6 .
Run time and hypervolume under different crossover probability and mutation probability.

Table 7 .
Results of different optimization models.

Table 8 .
Results of different risk levels ε.

Table 9 .
Wind power utilization of different iteration numbers.

Table 10 .
Cost of system operation of different iteration numbers.