Two-Stage Optimization Model for Two-Side Daily Reserve Capacity of a Power System Considering Demand Response and Wind Power Consumption

: Today, wind power is becoming an important energy source for the future development of electric energy due to its clean and environmentally friendly characteristics. However, due to the uncertainty of incoming wind, the utilization e ﬃ ciency of wind energy is extremely low, which means the problem of wind curtailment becomes more and more serious. To solve the issue of wind power large-scale consumption, a two-stage stochastic optimization model is established in this paper. Di ﬀ erent from other research frameworks, a novel two-side reserve capacity mechanism, which simultaneously takes into account supply side and demand side, is designed to ensure the stable consumption of wind power in the real-time market stage. Speciﬁcally, the reserve capacity of thermal power units is considered on the supply side, and the demand response is introduced as the reserve capacity on the demand side. At the same time, the compensation mechanism of reserve capacity is introduced to encourage generation companies (GENCOs) to actively participate in the power balance process of the real-time market. In terms of solution method, compared with the traditional k-means clustering method, this paper uses the K-means classiﬁcation based on numerical weather prediction (K-means-NWP) scenario clustering method to better describe the ﬂuctuation of wind power output. Finally, an example simulation is conducted to analyze the inﬂuence of reserve capacity compensation mechanism and system parameters on wind power consumption results. The results demonstrate that with the introduction of reserve capacity compensation mechanism, the wind curtailment quantity of the power system has a signiﬁcant reduction. Besides, the income of GENCOs is gradually increasing, which motivates their enthusiasm to provide reserve capacity. Furthermore, the reserve capacity mechanism designed in this paper promotes the consumption of wind power and the sustainable development of renewable energy.


Background and Motivation
With the excessive consumption of fossil energy, energy transformation has been paid attention by all countries worldwide. Renewable energy is favored because of its cleanness, environmental protection, and renewability [1]. Therefore, renewable energy generation will become the main method of electric energy production in the future. However, these clean energy sources also have some disadvantages. For instance, some kinds of renewable energy supply such as PV (photovoltaic) and wind power, due to the randomness and fluctuation during the generation process, will increase the Sustainability 2019, 11, 7171 2 of 22 volatility of system power generation output. What is more, with their large-scale penetration, the uncertainty of market balance will also increase, resulting in the power system operation facing serious challenges. For this reason, countries around the world have not been very effective in large-scale use of clean energy. In China, the utilization efficiency of wind energy is very low, and wind curtailment is a serious problem [2]. In 2017, the wind curtailment rate in Xinjiang and Gansu province reached 33% and 29% respectively [3]. Consequently, how to use renewable energy such as wind energy efficiently and reduce spillage of wind power is the main problem to be solved in this paper. Specifically, we further discuss the root cause of wind spillage, and we found that it is mainly the uncertainty of wind in nature that leads to the failure of accurate prediction of wind power output, which leads to the deviation of real-time wind power output, and finally affects the scheduling and production arrangement of the real-time spot market. In order to realize large-scale consumption of wind power, reserve capacity in the real-time stage is necessary, but a free reserve makes power generation companies gradually lose interest [4]. Therefore, the establishment of a compensation mechanism for reserve capacity will encourage more and more power generation manufacturers to provide more reserve capacity in the real-time stage, offer sufficient reserve space for wind power, and promote wind power consumption.

Literature Review
With large-scale wind power penetrating into the power system, the significant increase of the frequency modulation (FM) reserve unit is an inevitable trend [5]. However, considering the uncertainty of wind power, how to reasonably arrange the production of these reserve units is the main issue we need to focus on. In view of this problem, many studies have been done in previous research. The following is an overview of the methodology used by predecessors' studies to solve this problem.
For the uncertainty processing of wind power, most of the processing methods adopted in the previous research can be divided into two categories, namely robust optimization and stochastic optimization. References [6][7][8][9][10] adopt the idea of robust optimization to cope with the volatility of wind power. Reference [6] studies the wind power uncertainty of gas-electric integrated system coordination optimization and puts forward a two-stage distributionally robust optimization (DRO) model with the goal of minimum system total operating cost to analyze the impact of residents of gas load, abandon the wind punishment cost coefficient, and climb rate parameters on the final optimization results. Based on robust optimization, the literature [7] establishes an optimization model for joint scheduling of energy and auxiliary services markets that considered wind power uncertainty in the real-time power market and proposes redundant constraint reduction strategy to improve the computational efficiency of robust joint scheduling. Reference [8] utilizes distributionally robust chance constrained and interval optimization to handle the impact of wind power fluctuation. In reference [9], a robust unit commitment model with multi objectives (minimizing the operating cost and maximizing the peak shaving capacity of the power system) is proposed to resist the disturbance of wind power stochastic fluctuation. In addition, reference [10] utilizes robust optimization to deal with the uncertainty of wind power and photovoltaic power output, and the Gaussian process method is used to predict the confidence interval of uncertain parameters, which support a day-ahead market scheduling strategy formulation. However, both robust optimization and distributionally robust optimization have drawbacks, which are that scheduling decisions are too conservative and at the expense of economic interests.
Therefore, stochastic optimization, another method to solve uncertain problems, is favored by researchers. Reference [11] formulates a stochastic optimization model to accomplish the day-ahead economic dispatch fully considered wind power volatility. In reference [12], a scheduling model for the coordinated operation of hydro and wind power is established by utilizing the complementary characteristics of hydro and wind energy, and a stochastic programming model is adopted to cope with the uncertainty of wind power output. In addition, reference [13], in considering the influence of solar radiation intensity and wind power uncertainty on the whole scheduling process, establishes a stochastic optimal scheduling model for centralized wind-solar power stations. Moreover, a day-ahead stochastic scheduling model of thermal-hydro-wind-photovoltaic power system is presented in reference [14], which adopts Monte Carlo simulation to generate the output scenarios of wind power and photovoltaic power to describe their uncertainty. Yet all of the above studies did not consider the impact of wind power uncertainty on the power system concretely. Specifically, wind power in the electricity spot market is divided into two stages, i.e., the first stage is day-ahead market dispatch phase, and the second stage is a real-time market, whereas the actual deviation of the wind power output often occurs in the second stage. Therefore, the established stochastic programming model divided into two stages should be taken into account, which can restore the real-time situation of actual scheduling.
Generally, the consumption of wind power not only considers the reserve capacity market but also other methods to realize the safe and reliable connection of wind power to the grid. For example, demand response management (DRM) [15], joint operation of energy storage units, and other cooperative operation mechanism can balance the deviation of wind power real-time output. However, many researchers tend to consider these in terms of reserve capacity and ignore other approaches that can also promote wind power consumption. In reference [16], a dynamic economic scheduling model of integrated power system contained wind power is proposed to optimize the generation and reserve capacity of all units, and the capacity and ramp rate of the generator have an important impact on the optimal scheduling of reserve capacity, but the demand response of the user side is not mentioned. What is more, reference [17] introduces the concept of interruptible load and presents a cost-benefit tradeoff-based algorithm for optimal reserve capacity in the power market, which emphasized that both generator sets and interruptible load could participate in the reserve market. In addition, some of the literature also adopts the joint operation of energy storage device and wind turbine to reduce the interference of wind power instability on the system [18][19][20]. For instance, reference [18] proposes a wind power, pumped storage, and thermal plant joint operation mechanism with the objective of maximizing the operator's profit gained from the day-ahead market and reserve market, and a heuristic optimization algorithm has been used to solve this problem. In reference [19], a combination of wind turbines and air compression energy storage devices is introduced. Similar to reference [18], reference [20] also presents a decision-making mechanism for wind-energy storage joint system to participate in day-ahead market and reserve market. In summary, it can be found that the research on the scheduling mechanism of DRM as a flexible resource to balance wind power is relatively weak. Consequently, it is necessary to fill the research gap in this field.
According to the summary of the above literature review, it can be obviously concluded that previous research still has three insufficiencies. First of all, in terms of research methods, most of them adopt the idea of robust optimization, which leads to the over-conservative results. Second, the traditional stochastic optimization process cannot fully reflect the actual scheduling process when dealing with the uncertainty of wind power. Finally, on the demand side, the utilization of flexible demand response resources is insufficient.

Contributions and Organization
To fill the research gap, this paper proposes a novel two-stage stochastic optimization model for a wind power system with the participation of spot market, which considers both the flexibility of supply side and demand side, that is, the flexibility capacity includes interruptible load of demand response and reserve capacity of thermal power unit. The promotion effect of adding flexible capacity market in the real-time stage on the consumption of clean energy is fully studied, and the consumption situation of wind power is further analyzed. Additionally, the main contributions of this paper can be divided into three aspects. First, the model is constructed to establish a compensation mechanism when standby power is invoked in the real-time stage. Second, a two-stage stochastic optimization method is used to set up the scheduling optimization model, and the clustering method K-means-NWP is adopted to describe the uncertainty of wind power output. Finally, the flexibility of adjusting resources based on demand response is considered, which can increase the reserve power of the system and improve the consumption capacity of wind power. Besides, the organization of the remaining sections are as follows: Section 2 elaborates the main problems studied in this paper; in the Section 3, the mathematical model is built and the solution method is given; Section 4 uses a numerical example to verify the rationality of the established model; and the final conclusions are given in Section 5.

Problem Description
In this section, the main problem that needs to be addressed in this paper is described in detail, which contains the main problem introduction and the proposed optimization method process.

Wind Power Consumption
With the rapid development of renewable energy sources (wind power, photovoltaic, etc.) in China, wind power installed capacity increases year by year. In order to implement the large-scale consumption of wind power and ensure the safe and stable operation of the power system, the power system needs to arrange a certain rotation reserve based on the original operation mode to cope with the uncertainty of wind power output [21]. The influence of wind power on system dispatch and operation is mainly caused by the uncertainty of wind power output, which changes rapidly and violently and is unpredictable. The previous method of reserve configuration cannot meet the new requirements brought by large-scale wind power grid connection [22]. With the increasing proportion of wind power being grid-connected, the demand for a spinning reserve is increasing. It is not enough to rely on conventional power supply to adjust and provide a spinning reserve. Therefore, how to fully exploit the other regulatory resources in the system and rationally optimize the rotating standby according to the characteristics of the spinning reserve requirements of a wind power grid-connected system is an urgent problem that needs to be solved [23]. At present, flexible capacity resources (like the reserve capacity of thermal power units, demand response, energy storage, etc.) are the key to solving the problem of wind power consumption. This paper mainly considers the impact of thermal power unit flexibility, demand response, and wind penalties on wind power consumption.
In order to ensure the safety and stability of the system, accept more clean energy, and maintain a good system operation economy, it is urgent to make a reasonable and effective decision on the optimization of spinning reserve capacity of the system after large-scale wind power integration.

Proposed Capacity Optimization Configuration Process
In this paper, a two-stage stochastic optimization model for dual-side reserve capacity of power system considering wind power consumption is proposed. The overall goal of two-side reserve capacity focuses on minimizing the total cost of system operation and improving wind power consumption. Specifically, the stochastic optimization model can be divided into two stages. In the first stage, in the day-ahead market, the dispatching operator sorts the operating costs and spare capacity costs reported by each generator unit and determines the spare capacity that can be provided in the market on the next day. Then, according to the forecast value of wind power output, the dispatching operator calculates the reserve capacity that needs to be reserved for each generator unit in the second stage. In the second stage, the method of randomly generating the scene is used to simulate the operation of the market on the next day. When the wind power has error fluctuation, each generator unit provides backup power for its system and at the same time considers the interruptible load on the demand side to provide the unit with backup power. Finally, the impact on wind power consumption is analyzed by introducing abandonment penalty cost, demand response, and thermal power unit flexibility (climbing rate, minimum technical output) in the objective function.

Mathematical Modeling and Solution
In this section, we present a two-stage stochastic optimization model that mainly focuses on optimizing the dual-side reserve capacity of the power system to improve wind power consumption. In the first stage, the wind power output is predicted to determine the system backup capacity that needs to be generated in the second stage. In the second stage, the wind power output is clustered according to historical meteorological data by the K-means-NWP method in order to simulate the wind power output on the next day, then according to the error between the wind power output simulation value and the actual value. The remaining reserve capacity in the day-ahead phase can be utilized to balance the error.

Mathematical Modeling
This paper presents a two-stage stochastic optimization model. In the first stage, the independent system operator (ISO) ranks the operation cost and reserve capacity cost reported by each unit in the day-ahead market. With the aim of minimizing the total operating cost of the system, the reserve capacity available on the intra-day market the next day can be determined by optimization [24]. Then according to the typical output scenario of wind power output, the reserve capacity provided by other units is calculated. In the second stage, the K-means method based on meteorological data is used to simulate the operation of the intra-day-market in the second day. In case of fluctuations in wind power, thermal units provide reserve capacity for the system, and meanwhile, interruptible load on the demand side is considered to provide reserve electricity (interruptible load refers to the load that can be interrupted by the user side when the power grid needs it, which is usually realized by signing a contract with the user). Considering the compensation price of reserve electricity supply from thermal unit and interruptible load in real-time market, the scheduling and production arrangement in the second stage is completed with the objective function of minimum system equilibrium in the real-time stage [25]. Finally, the effect on wind power consumption is analyzed by introducing penalty coefficient of wind curtailment, demand response, and flexibility of thermal units (ramp rate and minimum technical output) into the objective function.
where a i , b i , c i respectively represent the coefficients of the generating cost function of thermal units; µ on g,t , µ o f f g,t denote the start and stop status of thermal units; ρ u g,r , ρ d g,r refer to the compensation price of reserve capacity of thermal units; ρ ω is the probability of wind power combination scenario ω; Ω represents the total number of wind power combination scenarios, which equals to the product of the number of output scenes of all wind turbines; λ t is the user electricity price; ρ u c , ρ d c represent the compensation price of standby power provided by thermal units in the real-time stage; c s denotes the penalty coefficient of wind curtailment; W spill ω represents the wind abandon quantity of wind power combination scenario ω.

The Constraints of the Day-ahead Stage
(1) Day-ahead system power balance where N G , N W respectively represent the number of thermal power and wind power units. P Da g,t is the power of thermal power unit g at time t in the day-ahead stage. P Da w,t denotes the power of the wind turbine g at time stage t in the day-ahead stage. L t refers to the total load of the system at time t. (2) Day-ahead unit output constraint where µ g,t is the start and stop state of thermal unit g at time stage t. P max g , P min g separately denote the maximum and minimum power of thermal unit g.
(3) Day-ahead wind power output constraints where P w,t is the predictive power of wind turbine w at time stage t. (4) Day-ahead wind power ramping constraints where r u w , r d w refer to the ramp rate of wind turbine w. (5) Day-ahead thermal units startup and shut off constraints [26] µ (6) Day-ahead thermal units ramping constraints where r u g , r d g refer to the ramp rate of thermal unit w. (7) Day-ahead spinning reserve capacity constraints Upper reserve capacity provided by thermal units at time stage t: Lower reserve capacity provided by thermal units at time stage t: Lower reserve capacity required by wind turbine: Upper reserve capacity needed by wind turbine: Sustainability 2019, 11, 7171 where R up,max g,t and R down,max g,t respectively refer to the maximum of the upper reserve capacity and lower reserve capacity provided by thermal unit g at time stage t. R up w,t , R down w,t separately refer to the maximum of the upper reserve capacity and lower reserve capacity required by wind turbine w at time stage t.
where P Rt w,t is the power of wind turbine w at time stage t in the real-time stage; W spill w refers to the wind abandon quantity of wind turbine w. (9) Real-time units output constraints where P Rt g,ω,t denotes the power of thermal unit g at time stage t in the wind power combination scene w at the real-time stage. R up g,ω,t and R down g,ω,t separately denote the upper reserve capacity and lower reserve capacity provided by thermal unit g at time stage t in the wind power combination scene w at the real-time stage. (10) Real-time reserve capacity constraints where u up g,t and u down g,t respectively refer to the state provided by thermal unit g at time stage t in the real-time stage.
where R IL ω,t denotes the interruptible load provided by users at time stage t in the wind power combination scenario ω in the real-time stage. θ represents the proportion of the maximum of interruptible load that can be provided by the power side to the total load [27]. (11) Real-time reserve service (only one in use at any given time) Sustainability 2019, 11, 7171 8 of 22 (12) Real-time wind power output constraints (13) Real-time wind power ramping constraints

K-Means Classification Based on Numerical Weather Prediction (K-Means-NWP)
The K-means clustering method is the most classic one in the dynamic clustering algorithm. Its basic idea is to divide each sample into the category closest to the mean, which is clustered by the distance and near distance.
The K-means clustering algorithm generally includes the following processing steps: (a) Dividing all data into K initial classes, selecting K sample points as the initial cluster center, and record as z 1 (l), z 2 (l) . . . . . . z k (l); Where, the initial value l = 1; (b) All samples are assigned to the K ω j (K) class represented by each cluster center according to the nearest neighbor rule, The number of samples included in each category is N j (l); (c) Calculating the mean vector of each type and use the vector as the new cluster center: where j = 1, 2 . . . . . . k, i = 1, 2 . . . . . . N j (l); (d) When z j (l + 1) z j (l), it means the clustering result is not optimal, then returns to step (b) and continues the iterative calculation; (e) When z j (l + 1) = z j (l), it means that the clustering result is optimal at this time and the iteration process is over.
According to the literature [28], there is a close relationship between wind power generation and meteorological data, especially with wind speed. However, the existing statistical forecasting methods based on numerical weather forecast do not classify meteorological information according to different characteristics; that is, they do not consider the different characteristics of wind power output corresponding to different meteorological data. Therefore, it is necessary to solve the following problems: 1) different characteristics of meteorological information correspond to large differences in wind power output, and 2) various factors in meteorological data are related to wind power output; that is, the factors that play a leading role in wind power output. Based on this, this paper comprehensively considers various meteorological factors and proposes a method of wind power output prediction based on fuzzy k-means of numerical weather forecast (K-means-NWP). The specific steps are as follows: Step 1 The meteorological factors affecting wind power output are analyzed, including daily atmospheric pressure average (Pav), daily wind speed minimum (Vim), daily wind speed maximum (Vmax), daily temperature minimum (Tmin), daily temperature maximum (Tmax), and daily wind direction, the sinusoidal mean (Dsin), and the daily wind cosine mean (Dcos). Step 2 Select the historical meteorological data of the previous month, x i = [Pav Vmin Vmax Tmin Tmax Dsin Dcos], then X = [x 1 ; x 2 . . . . . . ; x n ]; n = 30. Since the dimensions of the components in the NWP vector are different, normalization is required. The air pressure, wind speed, and temperature are divided by their respective historical maximum values. The wind direction sine and cosine values are normalized values and need no further processing.
Step 3 Set the initial value k of clustering number according to the number of target clusters.
Step 4 According to Equation (27), calculate the mean vector of each type and use the vector as the new cluster center.
Step 5 According to the fuzzy k-means criterion, judge the final result. If it is not the optimal result, return to step 4 to continue the iteration. Otherwise, the iteration ends, and the clustering result is output.
Step 6 Observe the meteorological data of the target day, calculate the distance between them and the clustering results and normalize them as the membership degree of the corresponding clustering results.
The flow chart of K-means-NWP steps is as shown in Figure 1.
power output prediction based on fuzzy k-means of numerical weather forecast (K-means-NWP). The specific steps are as follows: Step 1 The meteorological factors affecting wind power output are analyzed, including daily atmospheric pressure average (Pav), daily wind speed minimum (Vim), daily wind speed maximum (Vmax), daily temperature minimum (Tmin), daily temperature maximum (Tmax), and daily wind direction, the sinusoidal mean (Dsin), and the daily wind cosine mean (Dcos Step 3 Set the initial value k of clustering number according to the number of target clusters. Step 4 According to Equation (27), calculate the mean vector of each type and use the vector as the new cluster center.
Step 5 According to the fuzzy k-means criterion, judge the final result. If it is not the optimal result, return to step 4 to continue the iteration. Otherwise, the iteration ends, and the clustering result is output.
Step 6 Observe the meteorological data of the target day, calculate the distance between them and the clustering results and normalize them as the membership degree of the corresponding clustering results.
The flow chart of K-means-NWP steps is as shown in Figure 1. Step 1 Step 2 Step 3 Step 4 Step 5 Step 6 In this paper, in order to simplify the calculation, we take n = 3 to get three scenarios of high, medium, and low wind power output, which is expressed as [K h ; K m ; K l ], and then calculate the membership degree of the target day meteorological data corresponding to the three scenarios of high, medium and low wind power output, which is expressed as [ρ h , ρ m , ρ l ]. Thus, the predictive value of target daily wind power output in each scenario can be obtained. [ min max min max sin cos]

Solution Methodology
In this paper, in order to simplify the calculation, we take n = 3 to get three scenarios of high, medium, and low wind power output, which is expressed as  

Numerical Example
In the numerical example, this paper adopts the basic data of related physical components in the previous literature, and on this basis, the model algorithm proposed in Section 3 is used to realize the simulation, which studies the effect of reserve capacity compensation mechanism on GENCOs' revenue and quantity of wind spillage. Meanwhile, the impact of changes in related parameters on spare capacity allocation, reserve costs, and wind power consumption is investigated, and the calculated results also meet the expectations, verifying the effectiveness and applicability of the proposed method.

Parameters Setting
In this section, the power system contained thermal units, wind power generators and interruptible load is set up, and necessary input information and data is presented [29,30]. We build a power system that aggregates four thermal units and two wind turbines. The install capacity of the two wind turbines is 150 MW and 100 MW. Table 1 demonstrates the key parameters of the GENCOs, and Table 2 shows the parameters of electricity market. In addition, we divide spot trading into 96 time periods to simulate scheduling. And the above parameters are substituted into the model, and the model can be solved by Matlab R2016b and CPLEX Optimization Studio v 12.8.

The Influence of Compensation Mechanism
A two-stage stochastic optimization model proposed in Section 3.1 introduces the reserve capacity compensation mechanism, which provides compensation for the demand response provided by thermal units and the user side. In this section, two different scenarios are set up. By analyzing the operation results of the model, the effect of the introduction of reserve capacity compensation mechanism on wind power consumption and generator revenue is discussed. In order to observe more clearly the influence of compensation mechanism on wind power consumption, the penalty coefficient of wind curtailment is 0 in the two scenarios discussed in this section.

Case 1
In the traditional power system without reserve capacity compensation mechanism, the total system cost is composed of the operating cost and startup and shut off cost of thermal units. From the perspective of thermal power plants, their revenue is generated by day-ahead dispatching minus operation cost and startup and shut off cost.
Case 2 After the reserve capacity compensation mechanism is added into the conventional system, the power plant can obtain the compensation by reserving capacity in the day-ahead market, and meanwhile, the reserve power compensation can be obtained according to the net adjustment amount in the real-time market under the compensation mechanism. In addition to thermal units, the real-time reserve power provided by interruptible load will also obtain compensation. The total system cost increases the reserve capacity cost of thermal units, the compensation cost of calling thermal power plant and the reserve power of interruptible load on the original basis. Similarly, thermal power plants increase the income of reserve capacity and reserve power compensation.

Result and Discussion (High-Medium-Low)
In the day-ahead market, the reserved reserve capacity of thermal units 1-4 and the predicted output of wind turbines 1 and 2 can be obtained according to the calculation results in the first stage of the model, among which the wind power output is the predicted value. According to the clustering method in Section 3.2, the actual output scenarios of wind power are classified as high, medium, and low. The probabilities of each scenario are shown in Table 3. Different dispatching results can be obtained under different wind power output combinations. This paper selects five typical combinations-High-High (HH), High-Low (HL), Medium-High (MH), Medium-Medium (MM), and Low-Low (LL). The influence of compensation mechanism on wind power consumption and power generation benefit is analyzed under different wind power output combination. Figure 3 shows the comparison of wind curtailment results in scenario 1 and scenario 2 under different wind power output combinations. It can be seen that the wind curtailment quantity in scenario 1 is significantly less than that in scenario 2, which intuitively illustrates the influence of the reserve capacity compensation mechanism on promoting wind power consumption. In the real-time stage, the thermal unit provides upper reserve capacity R t,up g,ω and lower reserve capacity R t,down g,ω , while the interruptible load IL provides upper reserve capacity R t,IL g,ω . Then, to promote system balance and wind power consumption, the net reserve power generated in the real-time stage can be expressed as  Table 4 and Figure 4, when there is a reserve capacity compensation mechanism, thermal units are more active in providing reserve capacity to maintain system balance compared with the net reserve capacity in scenario 1 and scenario 2.
stage, the thermal unit provides upper reserve capacity  Table 4 and Figure 4, when there is a reserve capacity compensation mechanism, thermal units are more active in providing reserve capacity to maintain system balance compared with the net reserve capacity in scenario 1 and scenario 2.    Table 5 demonstrates the composition of the total system cost and the benefits of the power plant under the two scenarios. According to the wind curtailment result in Figure 3, the wind curtailment penalty cost in scenario 2 is also greatly reduced compared with scenario 1 due to the small quantity of wind curtailment. In addition, since scenario 2 absorbs more wind power, and the power generation cost of the wind turbines is 0, its power generation economy is obviously better than that of thermal power units, so the total system cost of scenario 2 is far less than scenario 1. In terms of the income of the GENCOs, the GENCOs in scenario 1 can only obtain the electricity revenue of pre-dispatching in the day-ahead stage. In scenario 2, due to the compensation mechanism of reserve capacity, GENCOs can obtain the compensation of reserve capacity reserved in the day-ahead stage and also gain revenue when the upper reserve power is provided in the real-time stage. Consequently, in scenario 2, the profit of the GENCOs is greatly increased, which can more effectively promote their enthusiasm to provide reserve capacity. Sustainability 2020, 12, x FOR PEER REVIEW 14 of 23   Table 5 demonstrates the composition of the total system cost and the benefits of the power plant under the two scenarios. According to the wind curtailment result in Figure 3, the wind curtailment penalty cost in scenario 2 is also greatly reduced compared with scenario 1 due to the small quantity of wind curtailment. In addition, since scenario 2 absorbs more wind power, and the power generation cost of the wind turbines is 0, its power generation economy is obviously better than that of thermal power units, so the total system cost of scenario 2 is far less than scenario 1. In terms of the income of the GENCOs, the GENCOs in scenario 1 can only obtain the electricity revenue of predispatching in the day-ahead stage. In scenario 2, due to the compensation mechanism of reserve capacity, GENCOs can obtain the compensation of reserve capacity reserved in the day-ahead stage and also gain revenue when the upper reserve power is provided in the real-time stage. Consequently, in scenario 2, the profit of the GENCOs is greatly increased, which can more effectively promote their enthusiasm to provide reserve capacity.

Impact Factors Analysis on the Wind Power Consumption
From the analysis of Section 4.2, it can be seen that the existence of compensation mechanism in the system is conducive to maintaining the balance of the system and promoting the consumption of the system's wind power. From the perspective of the GENCOs, the benefits are also increased. This section mainly discusses the influence of the change of relevant parameters in the system on wind power consumption, and the analysis is based on the market with backup capacity compensation mechanism.
In the base scenario, the tariff discount rate is 0.2, the maximum proportion of interruptible load IL in the total load is 0.005, and the penalty coefficient of wind curtailment is 800. Figure 5 draws the interruptible load provided by the power side at various time nodes in different wind power scenarios under different tariff discount rates in the real-time stage, As can be seen from Figure 5, in the five different wind power combination scenarios, when ξ = 0.7, the power side provides the most interruptible load to balance the market. Table 6 illustrates the total amount of interruptible load provided by the power side in one day. Obviously, the interruptible load provided by the power side in the real-time stage will increase with the increase of the tariff discount rate. Compared with the scenario of a low tariff discount rate, the incentive of a high tariff discount rate can promote the power side to make up for the negative error of wind power output to a greater extent. The wind curtailment in various wind power scenarios with different discount rates are illustrated in Table 7. From the perspective of the overall operation of the system, the increase of the tariff discount rate is conducive to wind power consumption and reduce wind curtailment. Figure 6 denotes the pre-output of wind turbines at various time nodes in the day-ahead stage of the system under different tariff discount rates, and the total output of day-ahead wind power under different discount rates is shown in Table 8. The increase of tariff discount rate affects the output of units in the equilibrium stage and also changes the pre-output of wind turbines in the day-ahead stage. It can be seen that, as opposed the situation in Section 4.1, with the increase of the tariff discount rate, the pre-output of wind turbines in the day-ahead stage is gradually increased under the condition that the actual output of wind power in the real-time stage is the same, which is equivalent to reducing the pressure of balancing wind power fluctuations in the real-time stage. In this process, the results in Table 9 are produced because more wind power has completed pre-dispatching in the day-ahead stage. In addition, with the increase of the tariff discount rate, the net reserve power provided by other units in the real-time phase gradually decreases.       Another parameter that affects the demand response is the interruptible load limit that can be provided by the power side, represented by the proportion of interruptible load to the total load. The system operation results under different θ values are shown in Tables 10-12. Similar to the tariff discount rate, with the increase of the θ value, the wind curtailment quantity in each scenario decreases, and as the output value of pre-dispatch of wind power in the day-ahead stage increases with the increase of the θ value, the total reserve capacity provided by other units in the real-time stage also decreases.

Penalty Coefficient for Wind Curtailment
The change of penalty coefficient of wind curtailment also has a great influence on the result of the system operation. Different system operation results are obtained when the penalty coefficient of wind curtailment is the various values shown in Table 13. Figure 7 denotes the comparison of wind curtailment curves of the system in one day under different penalty coefficient of wind curtailment. Meanwhile, Table 14 demonstrates the sum of wind curtailment quantity in each wind power scene under different penalty coefficient of wind curtailment. It is clear that, in each wind power scenario, the higher the penalty coefficient of wind curtailment, the less the wind curtailment quantity of the system will be. It is obvious that, in the case of a high penalty coefficient of wind curtailment, the system will try to use the reserve capacity of thermal units to absorb wind power in order to reduce the total cost of the system. The change of penalty coefficient of wind curtailment also has a great influence on the result of the system operation. Different system operation results are obtained when the penalty coefficient of wind curtailment is the various values shown in Table 13. Figure 7 denotes the comparison of wind curtailment curves of the system in one day under different penalty coefficient of wind curtailment. Meanwhile, Table 14 demonstrates the sum of wind curtailment quantity in each wind power scene under different penalty coefficient of wind curtailment. It is clear that, in each wind power scenario, the higher the penalty coefficient of wind curtailment, the less the wind curtailment quantity of the system will be. It is obvious that, in the case of a high penalty coefficient of wind curtailment, the system will try to use the reserve capacity of thermal units to absorb wind power in order to reduce the total cost of the system.   Table 15 shows the sum of the net reserve capacity provided by the thermal unit and the load side in one day under different penalty coefficient of wind curtailment in the real-time stage. While the penalty coefficient of wind curtailment increases, the sum of net reserve capacity in the real-time stage gradually decreases. This is because, in the case that the actual output value of wind power is the same, the higher the penalty coefficient of wind curtailment, the more pre-dispatch the system will make to wind turbines in the day-ahead stage, as shown in Table 16. should be noted that, compared with the data in Table 9, Table 11, and Table 16, the change of penalty coefficient of wind curtailment has a more significant impact on wind power consumption than the change of wind turbine pre-output caused by the change of relevant parameters of demand response. The change of penalty coefficient of wind curtailment can affect wind power consumption to a greater extent.

Conclusions
In this paper, we introduce a reserve capacity compensation mechanism and propose a two-stage stochastic optimization model based on demand response and wind power consumption on the basis of considering the total cost of the day-ahead stage and the real-time stage. Since the actual output of a wind turbine is uncertain, this paper uses K-means-NWP scenario clustering method to process. The wind power output is divided into three scenarios-high, medium, and low-and the actual unit data and market data are used to analyze the system operation result and wind power consumption situation with or without a reserve capacity compensation mechanism and then to compare the influence of different parameters in the system on the wind power consumption result. The conclusions are as follows: (1) Operating results in a variety of wind power output combination scenarios all prove that after the introduction of reserve capacity compensation mechanism, the system's wind abandoning capacity is significantly reduced, and thermal units and interruptible loads are also significantly more active in providing reserve capacity, which is conducive to maintaining system balance and promoting wind power consumption.
(2) From the perspective of GENCOs, the reserve capacity compensation mechanism greatly increases the income of GENCOs and thus their enthusiasm to provide reserve capacity. This also plays an active role in promoting wind power consumption, stabilizing wind power output instability, and maintaining system balance.
(3) As the tariff discount rate increases, the power side provides more interruptible load in the balance stage. Compared with the scenario of a low tariff discount rate, the incentive of a high tariff discount rate can promote the power side to stabilize the negative error of wind power output to a greater extent. The increase of the tariff discount rate is helpful for wind power consumption and reduce wind abandoning. Similarly, with the increase of the proportion of the total load of interruptible load stations, the analysis of the system operation results also reached the same conclusion.
(4) The increase of the penalty coefficient of wind curtailment has a great influence on the system operation result. With the increase of the penalty coefficient of wind curtailment, the total wind abandoning quantity of the system is greatly reduced. The higher the penalty coefficient of wind curtailment is, the more pre-dispatching the system will make to the wind turbine in the day-ahead stage. Compared with the change of parameters ξ and θ, the change of C S has a more significant impact on the wind abandoning quantity and wind power consumption of the whole system.
In conclusion, this paper puts forward targeted suggestions for further promotion of wind power consumption, hoping to provide reference for subsequent research. However, there are other ways to promote wind power consumption such as energy storage. The energy storage system has a response speed of at least minute, and it can store the power in the valley load period and release it in the peak load period. Considering the contribution of an energy storage system in promoting wind power consumption is a worthy subject for further research, which will be further discussed in a future paper. In addition, the construction of the core model is still relatively ideal, which does not consider the blocking effect of the physical power grid in the actual operation process. Future research can be extended in this direction.

Nomenclature
Indexes ω wind power combination scenario index, running from 1 to Ω w wind turbine index, running from 1 to N W g thermal turbine index, running from 1 to N G t Time step index, running from 1 to T Constants a i , b i , c i the coefficients of the generating cost function of thermal units λ t the user electricity price P max g the maximum power of thermal unit g P min g the minimum power of thermal unit g r u w , r d w the up and down ramp rate of wind turbine w R up,max g,t the maximum of the upper reserve capacity provided by thermal unit g at time stage t R down,max g,t the maximum of the lower reserve capacity provided by thermal unit g at time stage t R up w,t the maximum of the upper reserve capacity required by wind turbine w at time stage t R down w,t the maximum of the lower reserve capacity required by wind turbine w at time stage t R up g,ω,t the upper reserve capacity provided by thermal unit g at time stage t in the wind power combination scene w at the real-time stage R down g,ω,t the lower reserve capacity provided by thermal unit g at time stage t in the wind power combination scene w at the real-time stage u up g,t , u down g,t the state provided by thermal unit g at time stage t in the real-time stage ρ u g,r , ρ d