A Joint Scheduling Optimization Model for Wind Power and Energy Storage Systems considering Carbon Emissions Trading and Demand Response

To reduce the influence of wind power random on system operation, energy storage systems (ESSs) and demand response (DR) are introduced to the traditional scheduling model of wind power and thermal power with carbon emission trading (CET). Firstly, a joint optimization schedulingmodel for wind power, thermal power, and ESSs is constructed. Secondly, DR and CET are integrated into the joint scheduling model. Finally, 10 thermal power units, a wind farm with 2800MW of installed capacity, and 3 × 80MW ESSs are taken as the simulation system for verifying the proposed models. The results show backup service for integrating wind power into the grid is provided by ESSs based on their charge-discharge characteristics. However, system profit reduces due to ESSs’ high cost. Demand responses smooth the load curve, increase profit from power generation, and expand the wind power integration space. After introducing CET, the generation cost of thermal power units and the generation of wind power are both increased; however, the positive effect of DR on the system profit is also weakened.The simulation results reach the optimumwhen both DR and CET are introduced.


Introduction
The implementation of China's energy-saving and pollutant emission reduction strategies has prompted large-scale wind power development.In 2014, the installed capacity of wind power reached 115 million kW, ranking first in the world.However, influenced by the intermittency that is characteristic of wind power, the growth rate of wind power grid integration is smaller than the growth rate of the installed capacity.This phenomenon leads to a high rate of curtailed wind power in China.The average rate of curtailed wind power is approximately 12.8%.Especially in the "three north areas," the rate of curtailed wind power has already reached 15.4%.In order to solve the problem of curtailed wind power, suitable backup service should be provided on the generation site for wind power connected to the grid.Energy storage systems (ESSs) could flexibly provide backup service by charging and discharging.This property gives ESSs the most potential as a means to provide backup service for wind power integration.Additionally, demand response (DR) could optimize customers' power consumption behavior and incentivize customers to participate in system scheduling for wind power consumption.
Currently, China is implementing carbon emission trading (CET) pilot projects and plans to establish a CET market to prompt energy-saving and emissions reduction in the "thirteenth five-year plan" period.CET could affect the generation cost of thermal power units and make clean energy generation more advantageous.Hetzer et al. [1] regard carbon emission as a virtual network flow to subsequently build a theoretical framework for carbon emissions from the power system based on analyses of carbon emission trading and developing trends in the power industry.References [2][3][4] study carbon emission right definition problems in crossregional power trading and build a fair allocation principle based on the construction of a mathematical model that tracks carbon flow.The above research shows that CET can affect the cost of thermal power generation and make grid 2 Mathematical Problems in Engineering integration of clean energy more advantageous.Therefore, a study of how to integrate ESSs, DR, and CET into the traditional generation scheduling of a power system with wind power has important practical significance.
Selecting a better backup method is the most effective way to overcome the random nature of wind power.Currently, backup service for wind power grid integration consists of three main parts: thermal power, pumped storage power plants, and ESSs.Jiang et al. [5] develop a wind-thermal economic emission scheduling model considering coordination of the power allocation of thermal units and wind turbines.Yuan et al. [6] establish a multiobjective economic scheduling model for the hydrothermal-wind problem, considering the uncertainty cost of wind power.Wang et al. [7] propose a novel stochastic constraint model to solve for the uncertainty cost of wind power and present an improved particle swarm optimization (PSO) algorithm to solve the proposed model.Thermal power units can provide backup service for wind power by adjusting their start-stop condition, but their fuel consumption and pollution emission are not environmentally friendly.Therefore, pumped storage power stations are a better way to provide backup service.Papaefthymiou et al. [8] achieve high penetration levels of renewable energy in power systems by combining wind power and pumped storage power plants.Papaefthymiou and Papathanassiou [9] propose a novel unit commitment problem model and binary PSO algorithm to find the optimal schedule scheme.Ming et al. [10] calculate the effect of pumped storage power stations on wind power regulation and develop an economic evaluation model for combined wind power and pumped storage systems.Pumped storage power plants have the advantages of saving energy and reducing pollutant emission, but they are not suitable for large-scale application because they are restricted by geographic location and conditions.
The operational theory of ESSs is similar to that of pumped storage power plants, but ESSs have more flexible installation requirements with better prospects for largescale application.Wu et al. [11] take constraints on power generation units and energy storage units into consideration and build a static model for joint operation of wind power and energy storage.Hu et al. [12] combine opportunity and constraint theory and build a joint scheduling model for wind power and energy storage systems considering the uncertainty in wind power output, which is applied with good results.Ding et al. [13] build a wind-storage joint scheduling model considering risk constraints and use a Monte Carlo simulation method to simulate wind power output.García-González et al. [14] improve the power ramp mathematic expression for wind power and build a joint scheduling optimization model for the control of curtailed wind power and energy storage systems.
The researchers cited above have achieved good results in actual applications; however, the high initial investment cost limits the scale of application of energy storage systems.Therefore, other routes are needed to optimize the application of ESSs.DR can optimize customers' power consumption behavior, smooth the demand load curve, and increase the consumption of wind power.Greening [15] puts forward the basic concept of demand response.Niu et al. [16] classify DR into price-based demand response (PBDR) and incentivebased demand response (IBDR).López et al. [17] propose an optimization model for performing load shifting in the context of a smart grid.Nwulu and Xia [18] integrate game theory into dynamic economic emission scheduling considering DR.A framework for optimizing the bidding strategy of a smart distribution company that contains wind farms and responsive loads in the day-ahead energy market is proposed by Ghasemi et al. [19].Wang et al. [20] construct a modeling framework for an integrated electricity system in which loads become an additional resource.
CET can highlight the environmental friendliness of wind power and increase the advantages of wind power generation [21].Zhu et al. [22] developed a full-infinite fuzzy stochastic programming method for planning municipal electric power systems associated with greenhouse gas control under uncertainty.Khalid and Savkin [23] developed and tested a methodology for controlling the emissions from a group of microgenerators aggregated in a virtual power plant with wind power.The above papers discussed the impact of CET on system scheduling; the effects of collaborative optimization of CET, ESSs, and DR should be analyzed.
The rest of this paper is organized as follows.Section 2 puts forward a demand response mathematical model.Section 3 presents a model for charging and discharging ESSs.Section 4 establishes the joint scheduling optimization model for wind power and ESSs with and without CET.Section 5 takes 10 thermal units, a wind farm with 2800 MW of installed capacity, and 3 × 80 MW ESSs as the simulation system and comparatively analyzes the influence of ESSs, DR, and CET on system operation.Section 6 presents the primary conclusions.

Demand Response Model
Power demand response refers to a situation in which customers dynamically adjust their power consumption behavior according to price, which should guarantee a balance between power supply and demand.From an economic point of view, when electricity price increases, the demand for power should decrease.Part of the power demand during the peak load period will be transferred to other periods, while the other part will be reduced.Thus, for the peak load period, the load reduction consists of three parts: one is load transfer due to an increased electricity price, the second is load reduction, and the third is load transfer due to a reduced electricity price during the valley load period.For the float load period, the load change consists of two parts: one part comes from the peak load period and the other part goes to the valley load period.For the valley load period, the load increase consists of three parts: one is load transfer due to a decreased electricity price, the second part comes from the peak load period, and the third part is new power demand due to the price reduction.
This study defines the power demand during the peak load period, float load period, and valley load period before implementing time-of-use (TOU) price as  peak ,  flat , and  valley .Power demand at period  is   ; therefore, The proportion of power demand reduction during the peak load period is given by , the proportion of the load that transfers to other periods is  1 and the load loss proportion is 1 −  1 , and the proportion of the load that transfers to the float load period is  2 and the proportion that transfers to the valley load period is 1 −  2 .The proportion of power demand increase during the valley load period is , the proportion of the load transferred to the valley load period is  2 , and the new demand proportion is 1 −  2 .The proportion of the load that comes from the valley load period is  2 , and the proportion from the float load period is 1 −  2 .Thus, the demand during the peak, float, and valley load periods is calculated as follows: If the proportion of load change is the same at each point in time in the same period, the load at each time point is The system load period will change from   to    with the introduction of DR; load demand reduces during the peak load period and increases during the valley load period.The demand load curve becomes smoother after peak load shifting.

Charging and Discharging Model for ESSs
ESSs can be regarded as both power resources and load demand.When wind power output is high at night, ESSs are regarded as load demand.In the daytime, they are regarded as power resources to meet load demand.Charge and discharge of ESSs are limited by the system capacity.Assuming the storage energy of ESSs at time  is  , , the charging and discharging power balance should obey where  + , is the charging power at time ,  − , is the discharging power at time , and   is the loss coefficient for charging and discharging power.
Charging and discharging power from ESSs is limited as follows: where   is the upper limit for charging and discharging power.
In addition, the energy storage capacity of ESSs is limited: where  max  is the maximum storage capacity of the ESSs.

Mathematical Model without CET
4.1.1.Objective Function.Wind farm operators hope for higher consumption of wind power to gain more profits; however, this will cause more frequent adjustment of thermal power units for peak regulation, improving wind power grid integration but also increasing system coal consumption.To achieve the optimum energy efficiency, a joint optimization scheduling model for wind power and thermal power is built.
The maximum total profit of ESSs, wind power, and thermal power is taken as the optimization objective: where   ,   , and   are the profits of wind farms, thermal power units, and ESSs, respectively.The wind farm profit is calculated as follows: The profit of thermal power units is calculated as follows: where   is the benchmark price of thermal power,  , is the real-time power generation of unit  at time ,  , is the power consumption rate of unit ,  fuel is the fuel cost for power generation, OM , is the operation and maintenance cost of unit , and  , is the depreciation cost of unit .
The fuel cost of a thermal power unit is calculated by where  coal is the standard coal purchase price,  ,   ( , ) is the standard coal consumption, and  , is a binary variable.When a thermal power unit is shut down, coal consumption is zero, when the thermal power unit is operating, coal consumption is determined by the consumption characteristics function   (⋅) and the real-time power output  , : where   ,   , and   are coal consumption parameters of unit ,  , (1 −  ,−1 )SU  is the start-up cost of unit  at time , SU  is the start-up cost of unit ,  ,−1 (1 −  , )SD  is the shutdown cost of unit  at time , and SD  is the shutdown cost of unit : where  ,char is the electricity price when charging ESSs,  ,disc is the electricity price when discharging, and   is the fixed cost of the ESSs.

Constraint Conditions.
In the joint optimization model, the constraints of demand and supply balance, thermal power unit operation, wind power operation, and ESS operation should be comprehensively considered.
(1) System Demand and Supply Balance Constraint.Before DR, the system demand and supply balance constraint is described by After DR, the constraint is described by (2) Thermal Power Unit Power Generation Constraint.Realtime generation output is limited by the installed capacity and the minimum generation output: (3) Unit Ramp Rate Constraint.Depending on the technology level, unit generation output change is constrained by the adjacent period.The real-time output increment and decrement should obey (4) Unit Start-Stop Time Constraints.Frequent start-stop affects the performance of a thermal unit.The continuous unit start-stop constraint is shown as follows: Equation ( 17) is the shortest time constraint on unit . on ,−1 is the continuous running time of unit  at time  − 1.  on  is the shortest running time of unit .Equation ( 18) is the shortest shutdown time constraint on unit . off ,−1 is the continuous shutdown time of unit  at time  − 1.  off  is the unit shortest shutdown time.
(5) Wind Power Output Constraint.The real-time wind power output is constrained by the wind farm capacity: where   is the equivalent utilization efficiency and   is the total installed capacity of wind farm .
(6) Charging and Discharging Power Constraint on ESSs.For ESSs, the cumulative charging and discharging power should obey Therefore, if ESS operators hope to profit, the charging and discharging prices should obey (7) System Generation Reserve Constraints.When the power system is in operation, fluctuations may happen on both the generation side and the demand side.To ensure real-time balance, the supply of power should be adjusted to fall within a certain margin by increasing or reducing the power output: Equations ( 22)-(24) are the system upward spinning reserve constraints. max , is the maximum possible output of unit  at time . usr  is the upward spinning reserve demand, depending on thermal and wind generation power in the corresponding period.  is the maximum possible energy generation by unit  at time , adjusted for the installed capacity.Δ +  is the upward ramp rate, namely, the maximum power generation in the adjacent period.  is the thermal power unit reserve coefficient.  is the power reserve coefficient for wind turbines: Equations ( 25)-(27) are the system downward spinning reserve constraints. min , is the minimum possible output of unit  at time , restricted by two factors, namely, the minimum possible generation capacity under operation and the unit downward ramp rate constraint. dsr  is the system downward spinning reserve demand, depending on wind power in the corresponding period.  is the minimum possible generation capacity of unit , adjusted for the realtime minimum power output.Δ −  is the unit downward ramp rate, that is, the maximum power reduction generation unit in the adjacent period.

Mathematical Model with CET.
Currently, China is performing pilot construction and planning to establish a CET market in the "thirteenth five-year plan" period.CO 2 emission from the thermal power industry accounts for approximately 40% of the total.Generation rights displacement and the CET mechanism are both market mechanisms to optimize the thermal industry structure and reduce energy consumption and emission, which are consistent in purpose and results.
The marginal generation cost of thermal power changes under a CET mechanism, and carbon emission parameters are different due to different unit technologies, so the generation scheduling plan also changes.To maximize system profit under a carbon trading mechanism, this study builds an optimization model with the objective of maximizing thermal and wind power profit: Thermal power profit should meet the following conditions: where  CO 2 is the cost of carbon emission: where  CO 2 is the actual carbon emission of thermal units during the operation period,  0 is the total initial carbon emission right, and  CO 2 is the carbon trading price, which is related to the carbon trading demand.To simplify the model, this study assumes that the price does not change with the carbon trading demand.
The actual carbon emission of thermal units is related to the power load rate.Generally speaking, the actual carbon emission of units can be expressed as a quadratic function, similar to (22): where  CO 2 , ,  CO 2 , , and  CO 2 , are parameters of the carbon emission function.Then, total system emissions are as follows: Scheduling and operation constraint conditions for wind power and ESSs should be considered comprehensively in carbon emissions trading.The system demand and supply constraints, wind power unit operation constraints, and ESSs operation constraints are shown in ( 13) to ( 27).

Case Descriptions.
To analyze the impact of ESSs, DR, and CET on system operation, four cases are set up as follows.
Case 1 (baseline case).Self-scheduling of the system without DR and CET: DR and CET are not considered, so the impact of ESSs on wind power grid integration is analyzed alone in this case.The ESS capacity is 3 × 80 MW, the charging and discharging power of a single ESS unit is 20 MW, and the charging and discharging loss coefficient is 15%.
Case 2 (self-scheduling of the system with DR).Demand response is introduced into the joint scheduling.The demand load curve is divided into peak, valley, and float load periods according to the literature [24], which are listed in Table 1.The values of  1 and  2 are 0.95 and 0.7, respectively,  1 and  2 are 0.90 and 0.40, respectively, and both  and  are 5%.
Case 3 (self-scheduling of the system with CET).90% of CO 2 emissions from Case 1 are taken as the initial baseline for carbon emission trading, with a carbon emissions trading price of 80 ¥/t.Case 4 (self-scheduling of the system with DR and CET).Collaborative optimization of CET and DR is analyzed with joint scheduling of wind power and ESSs.

Basic Data.
This study uses 10 thermal power units and a wind farm with an installed capacity of 2800 MW as the simulation system.Coal consumption and carbon emission parameters of the thermal power units are listed in Table 2, and the operation coefficients are listed in Table 3.The grid purchase price of wind power is assumed to be 540 ¥/MW⋅h, the total operation and depreciation cost is 600 million/year, the grid purchase price of thermal power is 380 ¥/MW⋅h, and the price of standard coal is 800 ¥/t.The wind power unit equivalent utilization rate and system load distribution are set according to the literature [22] and listed in Table 4.

Simulation Results.
The simulation was implemented in GAMS optimization software using the CPLEX 11.0 linear solver from ILOG solver.The CPU time required for solving the problem for different case studies with an idea pad450 series laptop computer powered by a core T6500 processor and 4 GB of RAM was less than 10 s.

Case 1: Self-Scheduling of the System without DR or CET.
This case mainly analyzes the impact of ESSs on wind power grid integration.The scheduling result is shown in Figure 1, in which the maximum load is 2860 MW, the minimum load is 1230 MW, and the peak-valley ratio is 2.33.The maximum load for the thermal power equivalent output curve is 2458 MW, the minimum load is 204 MW, and the peak-valley ratio is 12.07.Figure 1 shows the wind and thermal power outputs.
After introducing ESSs, the peak-valley ratio is 2, the coal consumption rate is reduced from 326 kg/MWh to 322.5 kg/MWh, the system profit is enhanced by 160000 yuan, the electricity from wind power delivered to the grid is increased from 16840.5 MWh to 18620.6 MWh, and the curtailed wind power rate is reduced by 8.1%.ESSs can smooth the demand load curve, provide backup service for wind power, and reduce the start-stop cost of thermal power units.The scheduling optimization results for the power system with and without ESSs are listed in Table 5.
If ESS operators hope to maximize their economic benefit in the optimization period, they should discharge all stored energy to gain more economic benefit.However, to reduce the impact of wind power output fluctuation, ESSs make charging and discharging decisions based on wind power output to reduce thermal power peak regulation, as shown in Figure 2. The ESSs charge when wind power output increases and discharge when wind power output decreases.
Total profit declines with addition of ESSs, due to their high investment cost and the lack of large-scale commercial production.China has gradually become concerned with large-scale ESS development; in the long term, it has great potential with the establishment of price mechanisms and mature technology.

Case 2:
Self-Scheduling of the System with DR.The impact of DR on joint scheduling is analyzed in this case.DR can optimize customers' power consumption behavior and smooth the demand load curve.With  and  assumed to equal 3%, the system load curves for the two cases are shown in Figure 3.
According to Figure 3, the peak-valley differences for the two cases are 1400 MW and 1252 MW, respectively.The peakvalley ratios are 2.33 and 2. DR can smooth the demand load curve and reduce the curtailed wind power rate, which is 15.9% in Case 1 and 13.8% in Case 2. The utilization rate of wind power increases.The system scheduling optimization results for Case 2 are outlined in Table 6.

Case 3:
Self-Scheduling of the System with CET.The impact of CET on the joint scheduling problem is analyzed in this case.The CET price is 100 yuan/t.98% of total emissions in Case 1 are taken as the initial carbon emission permit, namely, 26171.73.The output of the thermal power units with different carbon trading prices is shown in Figure 4.
With a carbon trading price of 100 ¥/t, the output of wind power increased by 470.21 MWh, and the curtailed wind power rate decreased to 13.78%.CET would increase the cost of thermal power generation and change the system scheduling plan.The output of units with high carbon emission coefficient decreased, for example, unit 2 and unit 3. The output of units with low carbon emission coefficient rose, for example, unit 5.The optimization scheduling result in Case 3 is given in Table 6.
Compared with Cases 1 and 3, CET can increase the cost of thermal power generation and enhance the competitive advantage of wind power.System coal consumption decreased by 3.2072 million yuan, which is lower than that for Case 2. Profit increased to 2.9445 million yuan with increased wind power integration, which is lower than that for Case 2.

Case 4:
Self-Scheduling of System with DR and CET.The synergistic effects of DR and CET on the joint scheduling problem are analyzed in Case 4. The maximum wind power output was 19354.72 MWh.The curtailed wind power rate declined from 19.5% to 12.58%.After introducing DR and CET, system coal consumption decreased to 314.82 kg/MWh.The maximum system profit was 3.1272 million yuan.The system optimization scheduling results under different cases are shown in Table 6.
In summary, DR leads to proper power consumption behavior, while system operation obtains an optimal result when introducing DR and CET at the same time.The comprehensive operation result of Case 4 is better than Cases 2 and 3, which shows that a collaborative optimization effect exists for DR and CET.The curtailed wind power rate and coal consumption reach the minimum when DR and CET are both introduced.

Conclusions
The randomness of wind power output has become the primary impediment to wind power grid integration.To improve the system's capability to consume wind power, proper backup service should be provided on the generation side.Demand response should also be introduced to incentivize customers to respond to system scheduling.This paper constructs a joint optimization scheduling model for wind power and energy storage systems with CET and DR.The simulation results show the following: (1) ESSs can provide backup service for wind power and increase the system's ability to consume wind power.However, total system generation profits would decrease due to the high fixed cost of ESSs.
In the long term, China's large-scale ESS market has great potential with the establishment of a reasonable price mechanism and the development of mature energy storage technology.
(2) Demand response can optimize customers' power consumption behavior, smooth the demand load curve, and improve wind power grid integration.After introducing DR, wind power grid integration is enhanced, the system power structure becomes reasonable, thermal power costs are reduced, and the total system generation profit increases.
(3) Carbon emissions trading can increase the generation cost of thermal power units, enhance the advantages of integrating clean energy into the power grid, reduce the curtailed wind power rate, and optimize the system power structure.However, the total system generation profit is reduced due to the increased cost of thermal power generation.
(4) With the introduction of DR and CET, system operation and scheduling obtained an optimal result.The total system profit increased and the high cost of ESSs was reduced.In addition, the demand load curve was smoothed and the curtailed wind power rate and system coal consumption were reduced to the minimum.System profit increased to the maximum, which shows that a collaborative effect exists for DR and CET.
(5) Though the proposed model examined the optimization effect of ESSs, CET, and DR on wind power grid integration, some constraints related to practical applications are ignored for the sake of convenient analysis and reaching a conclusion.

Nomenclature
: Time index : Thermal power unit index : Wind farm index : Binary variable:

Figure 2 :
Figure 2: Charge and discharge optimization result for ESSs.

Table 1 :
Period division of TOU.

Table 2 :
Coal consumption and carbon emission parameters of thermal power units.

Table 3 :
Operation coefficients of coal-fired power units.

Table 4 :
Equivalent utilization of wind power units.

Table 5 :
Scheduling optimization result of power system in different cases.
Generation fuel cost OM , : Operation and maintenance cost of unit   , : Depreciationcostofunit  coal : Standardcoalpurchaseprice   ( , ): Standard coal consumption of unit  at time    ,   ,   : Coal consumption parameters of unit  SU  : Start-upcostofunit SD  : Shutdown cost of unit   ,char : Electricity price when charging ESSs  ,disc : Electricity price when discharging ESSs   : Fixed cost of ESSs  on ,−1 : Running time of unit  at time  − 1  off ,−1 : Shutdown timeof unit  at time  − 1  on  : Shortest running time of unit   off  : Unit shortest shutdown time of unit    : Equivalent utilization efficiency at time    : Total installed capacity of wind farm   max , : Maximum possible output of unit  at time   usr  : Upward spinning reserve demand at time    : Maximum possible energy generated by unit  at time  Δ +  : Upward ramp rate of unit    : Thermal power unit reserve coefficient   : Power reserve coefficient for wind turbines  min , : Minimum possible output of unit  at time   dsr  : Downward spinning reserve demand at time    : Minimum possible generation capacity of unit  Δ −  : Downward ramp rate of unit   CO 2 : Carbon emission cost  CO 2 : Actual carbon emission of thermal unit during operation period  0 : Total initial carbon emission right  CO 2 : Carbon trading price.