Abstract

To alleviate the issues of global warming and energy crisis, countries are vigorously developing renewable energy technology. The integration of large-scale renewable energy, including wind energy, hydropower, and photovoltaic (PV), has a great impact on system operation scheduling and economic dispatch. This paper presents an economic dispatching method of wind-PV-CSP-hydro-battery system with wind and photovoltaic power generation as the main energy sources. The long short-term memory (LSTM) neural network is applied to predict wind and PV power, besides, the Latin Hypercube Sampling (LHS) method and the synchronous reduction algorithm are used to obtain 10 typical wind and PV power scenarios. A day-ahead economic dispatch model of wind-PV-CSP-hydro-battery mathematical model is established, and relevant constraints are considered. Concentrated solar power (CSP), hydropower stations, batteries, and transferable loads are used as flexible resources to increase the penetration rate of wind and photovoltaic power generation. Finally, three cases are tested to demonstrate the feasibility of the proposed model. The results show that: (1) LSTM neural network can well predict the output power of wind and photovoltaic power generation with a small root mean square error (RMSE). (2) The introduction of transferable loads and CSP power station into the renewable energy power system can effectively reduce the fluctuation and curtailment rates of wind power and PV power generation.

1. Introduction

Climate warming has become an indisputable fact in the world. The Chinese government has promised to reach carbon peak and carbon neutralization in 2030 and 2060, respectively [13]. For the power industry, developing large-scale renewable energy and reducing the use of thermal power units are effective ways to reduce carbon emissions. The renewable energy power generation reached 22148 billion kWh in 2020, including 466.5 billion kWh for wind power and 260.5 billion kWh for photovoltaic power [4]. For large-scale renewable energy, the uncertainty and randomness of its output power limit its energy utilization, especially for the economic dispatching of renewable energy power systems [5].

To reduce the uncertainty of renewable energy in economic dispatch, scholars have conducted a lot of research on the power side and load side. On the power side, the intermittence of wind and photovoltaic power generation is considered in the renewable energy obligation model, which improves the penetration of wind and photovoltaic power generation [6]. The complementarity of wind and PV power generation is used to reduce its randomness [7]. In [8], the direct cost, under estimation penalty cost and over estimation penalty cost are applied to quantify the usage of wind power and photovoltaic power generation. To improve the accuracy of power generation calculation, illumination and wind speed are predicted in [9, 10]. In [1113], wind power and PV power plants are jointly dispatched, and economic dispatch is carried out considering the objectives of the maximum wind-solar joint benefit. In summary, the utilization rate and output characteristics of wind power and photovoltaic power generation can directly describe the uncertainty of their behavior. On the load side, users can follow the wind and photovoltaic power generation adaptively according to the power consumption curve. The flexible loads are employed in the scheduling model to suppress the uncertainty brought by large-scale wind power grid connections [14]. To decrease the uncertainty in wind power generation and improve the penetration rate of wind power generation, transferable loads on the demand side are investigated [15]. The interruptible load is regarded as a virtual unit to tackle the challenges of having large-scale wind power integration into the power system [16]. The demand-side response resources function is fully tapped by analyzing the difference of different loads’ responses to electricity price for solving the random fluctuation of wind power in [17]. In general, the rational use of flexible loads can reduce the uncertainty of wind power and photovoltaic power generation.

A short-term PV power generation output prediction method based on particle swarm optimization-deep belief network (PSO-DBN) is proposed in [18], although the PSO algorithm can optimal the DBN network, the optimization results are inevitably affected by the setting of parameters. References [19, 20] use the back propagation neural networks to predict wind power. However, this method may cause over-learning and aggravate prediction errors due to the increase in the number of nodes in the hidden layer. Therefore, it is necessary to find a suitable method to predict the output of wind and PV power based on the wind and PV power characteristics. The long short-term memory (LSTM) neural network is a kind of recurrent neural network (RNN) model that can deal with time series forecasting problems [21]. It can predict the information for the next time by learning the data of the previous time. Meanwhile, the transferable load has the best flexibility and the better users’ satisfaction with electricity consumption than reducible load [22]. It should be noted that the high penetration of wind power and PV in power system cannot be ignored. Concentrated solar energy (CSP) has an excellent performance in improving the penetration of renewable energy due to its thermal energy storage (TES) [2325]. References [2628] incorporate hydropower stations and batteries into wind power generation and PV power generation systems, which can maximize the system's ability to accept wind and PV power generation, as well as the economy and reliability of power system operation. In the economic dispatch of power system with renewable energy, there are still two problems: (1) the prediction method of renewable energy cannot effectively reflect the output characteristics of renewable energy; (2) to reduce the fluctuation of renewable energy, how to reasonably combine flexible resources.

To address these issues, this paper investigates a day-ahead economic dispatch model of the pure renewable energy power system. It takes wind power generation and photovoltaic power generation as the main power supply, and CSP power plant, hydropower station, battery, and transferable load as flexible resources. Meanwhile, LSTM neural network is used to predict the wind and PV power. Besides, the Latin Hypercube Sampling (LHS) method and the synchronous back reduction algorithm are applied to obtain 10 typical wind and PV power scenarios. Ultimately, the CPLEX is employed to get the optimal solutions for the proposed model. Moreover, three cases are used to prove the effectiveness of the model. Compared with the literature, this paper possesses the following contributions:(1)A day-ahead economic dispatch model of wind-PV-CSP-hydro-battery is established. Compared with references [710] including thermal power generation, this paper is more meaningful to the economic scheduling of purely renewable energy systems in the future.(2)The LSTM neural network is used to learn the output characteristics of wind energy and photovoltaic power generation, which can achieve a better prediction effect than the traditional neural network [19, 20].(3)The simulation results show that compared with hydro-battery [9] and hydro-CSP-battery [7], hydro-CSP-battery-translatable load as flexible resources can better reduce the fluctuation of wind power and photovoltaic.

The rest of this paper is arranged as follows: In Section 2, a day-ahead economic dispatch model of wind-PV-CSP-hydro-battery is established. Section 3 shows the forecasting method of wind and PV output and introduces the solving steps of the day-ahead economic dispatch model of wind-PV-CSP-hydro-battery. The relevant cases are set up for proving the validity of the model, and the corresponding results are obtained by case simulation in Section 4. Finally, some concluding remarks are given in Section 5.

2. Model Formulation

To improve the flexibility of power system and reduce the fluctuation of renewable energy, this paper establishes a day-ahead economic dispatch model with the minimum sum of the economic costs of wind power, PV power generation, CSP power station, hydropower station, storage battery, and transferable load. In the meantime, the relevant constraints for each power generation unit also are considered, as shown in Figure 1.

2.1. Objective Function

In our work, the cost objective function includes wind power generation cost, PV power generation cost, CSP power station operation cost, CSP power station spinning reservation cost, hydropower generation cost, battery cost, and transferable load cost. The mathematical expression of economic cost objective function is given as follows:where is the total cost, is wind power generation cost, is PV power generation cost, is CSP power station cost, is hydropower generation cost, is the battery cost, and is the compensation cost of the transferable load.

The cost of wind power generation and PV power generation is shown as follows:where is wind farms index, is the total number of wind farm, is the cost coefficient of the j-th wind farm, is the output power of the j-th wind farm at time . is PV power station index, is the total number of PV power station, is the cost coefficient of the k-th PV power station, and is the output power of the k-th PV power station at time .

The cost of CSP power station has two parts, the first part is the operation cost of CSP power station, and the second part is the spinning reverse cost, which is stated as follows:where is CSP power station index, is the total number of CSP power stations, is the operation cost of the l-th CSP power station, is the output power of the l-th CSP power station at time , is the spinning reverse cost of the l-th CSP power station, and is the spinning reverse power of the l-th CSP power station at time .

The cost of hydropower station and the cost of battery are as follows:where is the hydropower station index, is the total number of hydropower stations, is the cost coefficient of the m-th hydropower station, is the output power of the m-th hydropower station at time , is the battery, is the total number of batteries, is the cost coefficient of the n-th battery, and is the output power of the n-th battery at time .

The compensation cost of the transferable load is given as follows:where is the original operation time set of transferable load, is the transferable period of load, , are the start time and the end time of transferable period, is the compensation cost of transferable load at time , is the power transferred from time to , is the status of the transferable load, means the load transferred to time and otherwise.

2.2. Constraint Condition
2.2.1. Power Balance Constraints

To maintain the power balance of the system, the power balance constraint is established as Eq. (8). It is defined as the sum of the output power of j-th wind farm, the output power of k-th PV power station, the output power of l-th CSP power station, the output power of the m-th hydropower station, and the output power of the n-th battery equals the load of power system considering the transferable load. The specific formulas are as follows:where is the load of power system considering the transferable load, is the load of power system not considering the transferable load, and is the load transferred to time .

2.2.2. Operation Constraints of Wind Farm

Wind farm operation constraints include the upper and lower limits constraint of the output power of j-th wind farm and forecast power constraint. The specific formulas are as follows:where is the maximum power of the j-th wind farm, , are the forecast power and forecast error power of the j-th wind farm at time .

2.2.3. Operation Constraints of PV Power Station

PV power station operation constraints include the upper and lower limits constraint of the output power of k-th PV power station and forecast power constraint. The specific formulas are as follows:where is the maximum power of the k-th PV power station, and , are the forecast power and forecast error power of the k-th PV power station at time .

2.2.4. Operation Constraints of CSP Power Station

The CSP power station can provide a spinning reserve because it has thermal energy storage (TES) units. The constraints of a CSP power station include upper and lower limits of the output power, upper and lower limits of TES discharge and charging power, positive and negative spinning reserve constraints provided by TES of CSP power station, and capacity constraints of TES of CSP power station. The specific formulas are as follows [9]:where is the maximum power of the l-th CSP power station, , are the TES discharge power and the TES charge power of the l-th CSP power station at time , , are the positive spinning reserve and the negative spinning reserve provided by TES of the l-th CSP power station at time , , are the TES maximum discharge power and the TES maximum charge power of the l-th CSP power station at time , is state variable, means TES is discharging, while means TES is charging.

2.2.5. Operation Constraints of Hydropower Station

The constraints of hydropower station include upper and lower limits of the output power of hydropower station, storage capacity, and discharge capacity constraints of hydropower station, neglecting the spillage of hydropower station, which is shown as follows:where are the hydropower generation coefficients, , are the maximum output power and the minimum output power of the m-th hydropower station, , are the storage capacity of the reservoir of the m-th hydropower station at time and time , , are the maximum and the minimum storage capacity of the reservoir of the m-th hydropower station, is the outflow of the m-th hydropower station at time , , are the maximum and the minimum outflow of the m-th hydropower station, is the inflow of the m-th hydropower station at time , is an upstream hydropower station of the m-th hydropower station, is the time delay from the o-th hydropower station to the m-th hydropower station, and is the outflow of the o-th hydropower station by time .

2.2.6. Constraints of Battery

The constraints of battery include capacity constraints, charge and discharge power constraints, and charge and discharge state constraints. The specific formulas are as follows:where , are the capacities of the n-th battery at time and time , , are the maximum and the minimum capacities of the n-th battery, , are the charge and discharge power of the n-th battery at time , , are the maximum charge power and the maximum discharge power of the n-th battery, , are the states of charge and discharge of the n-th battery at time , , means not charge and charge of the n-th battery at time , , means not discharge and discharge of the n-th battery at time , , are the charge and the discharge efficiencies, respectively, and is the unit time interval.

2.2.7. Spinning Reserve Constraints of Power System

Spinning reserve constraints of power system include up spinning reserve and down spinning reserve. The specific formulas are as follows:where , are the up spinning reserve and down spinning reserve of the power system.

2.2.8. Constraints of Transferable Load

Based on not affecting the user's satisfaction with power consumption, the transferable load shall be limited to a certain time area, that is, only the transferable load within the allowable time shall be considered. The specific formula is as follows:(a)Transferable load within allowable time:where is transfer power of the transferable load.(b)Transferable load in disallowed time:

3. Solution Steps

3.1. Forecast Wind Power and PV Power

Wind power and PV power are uncertain and random. If the deterministic wind power and PV power data are used to solve the economic dispatch, it will influence the veracity of the economic dispatch results. Therefore, it is necessary to use historical data to predict wind power and PV power.

The LSTM neural network has a special gate structure, which solves the gradient vanishing problem in the RNN neural network [29]. It can predict the information of the next time by learning data of the previous time. Therefore, LSTM neural network is used to predict wind power and PV power. The specific structure of LSTM neural network is shown in Figure 2 [30].

In Figure 2, is the input variable of the LSTM neural network at time , , are the output variables of the LSTM neural network at time and , , are the memory units status at time and , is forget gate, is input gate, is output gate, is memory unit status, and is the sigmoid function. The mathematical formula of each gate is as follows [31]:where , are the recursive weight matrixes of , , , . , , , are the input weight matrixes of , , , , respectively, and , , , are the offset matrixes of , , , , respectively [32].

The explanation of each gate in Figure 2 is as follows: the forget gate determines the information to be retained and discarded through the sigmoid function, the input gate determines the information to be updated through the sigmoid function. is generated by the tanh function. The output gate firstly obtains the initial output through the sigmoid function and then multiplies the initial output with which through the tanh function layer at the time to obtain the final output .

3.2. Prediction Error Fitting of Wind Power and PV Power

Most literature regards the probability distribution of prediction error of wind power and PV power generation as having a normal distribution. However, according to references [33, 34], the t-Location Scale function may better reflect the probability distribution of prediction error of wind power and PV power generation.where is the mean, is the variance, () is the free degree.

3.3. Scenario Generation and Reduction of Wind Power and PV Power

LHS method has the characteristics of stratified sampling and disordered sorting [35], which avoids the disadvantages of clustering and repeated sampling. Therefore, the LHS method is adopted to generate 1000 wind power and PV power scenarios.

To simulate economic dispatch for 1000 scenarios could be highly time consuming. It is therefore necessary to reduce the number of scenarios. In this paper, the synchronous back reduction algorithm [36] is used to solve the problem. The solution flow chart is shown in Figure 3. The specific steps are given inTable 1.

4. Case Simulation

In this paper, the data of wind power and PV power generation are from Elia company in Belgium. The cost coefficients of wind power and PV power generation are adopted from reference [10], the relevant parameters of hydropower station are from reference [17], and the coefficient of CSP station comes from [35]. The relevant parameters of the storage battery are shown in Table 2. The system load data are shown in Table 3, and 5% of the system reserve is considered. The case simulation adopts a PV power station, a wind farm, two hydropower stations, two batteries and a CSP station, and three different cases as shown below are set to verify the mathematical model. Matlab R2019b is used to predict wind power and PV power generation, and the mathematic model is solved by CPLEX.

Case 1. Economic dispatch of wind-PV-hydro-battery.

Case 2. Economic dispatch of wind-PV-hydro-CSP-battery.

Case 3. Economic dispatch of wind-PV-hydro-CSP-battery considering transferable load.

4.1. Forecast of Wind Power and PV Power

The data of wind power and PV are from the measured data from July 1, 2020 to September 1, 2020 provided by Elia company in Belgium. The sampling time of wind power and PV is 1 hour. A total of 24 data sets are sampled a day, and the historical data of wind power and PV are 2208 points, respectively.

The LSTM neural network is used to predict wind power and PV. The hidden layer of LSTM neural network is taken 200 layers. The first 90% of the historical data sequence is used for training and the last 10% is used for testing. The solver is set to ‘Adam’. According to the different output characteristics of PV and wind power, PV shall carry out 250 rounds of training. We specify the initial learning rate of 0.005 and reduce the learning rate by multiplying the factor of 0.2 after 125 rounds of training. We conducted 300 rounds of wind power training and specified an initial learning rate of 0.005. After 150 rounds of training, multiplying by a factor of 0.2 reduces the learning rate.

To enhance the diagnostic accuracy of wind power output, K-means clustering method is adopted for processing the historical data of wind power, and the historical data of wind power with similar output characteristics are used for prediction three days before the prediction date. The forecasted wind power output data is shown in Figure 4. The output characteristics of PV power generation are obvious, and the historical data are not processed. The forecasted PV power generation is shown in Figure 5. The RMSE is adopted to measure the accuracy of scenery prediction. The expression of RMSE is as follows:where is the observed value, is the true value, and is the number of observations.

In Figures 4 and 5, ‘Observed’ is the real observation of the data to be predicted, ‘Forecast’ is the predicted result. The RMSE of wind power forecast is 190.6511 MW, and the RMSE of PV power generation forecast is 313.0239 MW. By observing the ‘Observed’ and ‘Forecast’ curves in Figures 4 and 5, it is found that the ‘Forecast’ curve can well reflect the power size and trend of the ‘Observed’ curve. And the RMSE is also small, so the LSTM neural network has good prediction results.

4.2. Scenario Generation and Reduction of Wind Power and PV

The instability and randomness of clean energy have negative effects on the safe and stable operation of power systems with access to clean energy, thus it is necessary to discuss the uncertainty of clean energy. In this paper, the LHS method is used to generate 1000 scenery scenarios to describe their uncertainty. To reduce the calculation pressure, the synchronous back reduction algorithm is used to reduce the scenery scenarios, and 10 representative scenarios are generated for economic scheduling research. 1000 wind power and 1000 PV power generation scenarios generated by LHS are shown in Figure 6(a) and Figure 6(b), respectively. The 10 typical wind power and PV power generation scenarios generated by synchronous back generation reduction algorithm are shown in Figure 6(c) and Figure 6(d), respectively, in which s1, s2 … s10 represent 10 scenarios.

The LHS data of wind power is the first 24 points in Figure 4, and the PV LHS data is the 25th to 48th points in Figure 6. 1000 samples are taken for measuring the prediction errors of wind power and PV, respectively, and 1000 groups of sampling data are superimposed on the prediction curve to obtain 1000 wind power output scenarios as shown in Figure 6(a) and Figure 6(b). To reduce the pressure of calculation, the 10 typical wind and PV output scenarios generated by the synchronous back reduction method are shown in Figure 6(c) and Figure 6(d). The 10 typical scenarios generated can well represent 1000 scenarios of wind power and PV power generation from the output range and output direction.

4.3. Influence of CSP Power Station on Wind Power and PV

The output results of each generation unit in Case 1 are shown in Figure 7(a), and the charge and discharge of the battery are shown in Figure 7(b). The results of Case 2 are shown in Figure 8(a); Figure 8(b) is the output of the storage battery and CSP power station. In the simulation picture, ‘WIND’ represents wind power, ‘PV’ is photovoltaic, ‘H1’ and ‘H2’ mean two hydropower stations, ‘BD’ and ‘BC’ represent battery discharge and charging, respectively, ‘CSPD’ and ‘CSPC’ are CSP power station discharge and charging, respectively, and ‘load’ represents the load of the system. If ‘BD’ and ‘CSPD’ are positive, it means that the battery and CSP power station are discharging, and ‘BC’ and ‘CSPC’ are negative, it means that the battery and CSP power station are charging. The abscissa is time and the ordinate is output power.

To reflect the influence of CSP power stations on the output of wind power and PV, the fluctuation rate of wind power and PV and the curtailment rate of wind power and PV are defined. The specific expression is as follows:(1)The fluctuation rate of wind power and PV:where , is the output power of wind power and PV in and , is the fluctuation rate of wind power and PV.(2)The curtailment rate of wind power and PV:where PV0 and WIND0 are the output power of PV station and wind power station, PV1 and WIND1 are the actually used PV and wind power in an economic dispatch, is the curtailment rate of wind power and PV.

To verify the effect of CSP power stations on system regulation performance, the wind power and PV output curves in Figure 7(a) and Figure 8(a) are compared. It can be found that the wind power and PV output curves in Figure 8(a) are relatively stable compared with those in Figure 7(a). Most obviously the load curve fluctuates greatly from 14 : 00 to 15 : 00, specifically, the PV output changes by 1066 MW from 14 : 00 to 15 : 00 in Case1, and by 703 MW in Case 2, which is due to the addition of CSP power station in Case 2. The output of the hydropower station is increased in Case 2 from 14 : 00 to 15 : 00, and stored the excess PV output in CSP power station and storage battery, making the PV output curve relatively stable. At the same time, it can be found from Table 4 that the fluctuation rate of wind power and PV is 5.4926 in Case 2, which is smaller than that of Case 1, which is 7.4328.

It can be seen from Figure 7(b) and Figure 8(b) that after adding CSP power station, excess wind power and PV are stored in CSP power station and battery at 1 : 00–6:00 and 12 : 00–17 : 00 with sufficient wind and PV output. At 20 : 00–23 : 00 with smaller wind and solar output, CSP power station and battery provide more output to make up for the lack of wind and solar output. At the same time, it can be found in Table 5 that adding CSP power station reduces the curtailment rate of wind power and PV.

4.4. Influence of Transferable Load on Wind Power and PV

To verify the influence of transferable load on system operation, transferable load is considered in Case 3. The fluctuation rate of wind power and PV, and economic cost in three cases are compared. The translatable intervals of translatable load are 9 : 00–10 : 00, 10 : 00–11 : 00, 11 : 00–12 : 00, and the acceptable intervals of translatable load are 1 : 00–2:00, 15 : 00–16 : 00. Figure 9 shows the load curve after translation, in which ‘load0’ is the load curve before translation and ‘load’ is the load curve after translation.

The load curve is the load curve after taking the translational load in Case 3. The output curve of each generation unit is shown in Figure 10(a), and the output of CSP power station and the battery is shown in Figure 10(b).

We can see from Figure 10 that after considering the transferable load, the peak-to-valley difference of the system load is reduced, which reduces the flexibility requirements for the output of each unit in the system. At the same time, it can be found from Table 5 that the fluctuation rates of the wind power and PV decrease after the introduction of the translatable load, which is 5.3891, and which is the lowest among the three cases, but its curtailment rates of wind power and PV rates are not much different from those in Case 2. Also, we can see from Table 6 that the economic cost of 4513 $ in Case 3 is higher than that of 4444 $ in Case 2, which is due to the introduction of transferable load compensation cost, but the economic cost is lower than that of Case 1.

In conclusion, the curtailment rate of wind power and PV rate and the fluctuation rate of wind and PV in Case 2 considering CSP power station are better than those of Case 1 without CSP power station. The fluctuation rate of wind power and PV of Case 3 considering CSP power station and transferable load is better than that of Case 2 considering only CSP power station and Case1 not considering CSP power station and translational load, its curtailment rate of wind power and PV is also better than that of Case 1.

5. Conclusion

In this paper, taking the carbon peak and carbon neutrality as the background, we investigated the economic dispatch of power systems with wind and PV generation as main power supplies. To improve the permeability of wind power and PV power generation and reduce the output uncertainty, this paper established an economic dispatching model including wind-PV-hydro-CSP-battery and transferable load. The model took the economic cost of the system as the objective function and considered the relevant constraints of wind, PV, hydropower station, CSP power plants, batteries, and transferable load. The following conclusions are obtained:(1)The wind and PV power generation predicted by LSTM neural network can well reflect their power output, and the RMSE is small in which wind power forecast is 190.6511 MW, PV power generation forecast is 313.0239 MW.(2)Compared with the traditional renewable energy system with wind power and photovoltaic as energy source, the proposed wind-PV-hydro-CSP-battery system considering transferable load reduces the fluctuation rate of wind and PV power by 27.5%, the reduction of wind and photovoltaic power generation curtailment rate by 2.6% and the economic cost by 1.1%, which prove the effectiveness of the combination of hydropower, CSP power station, battery and transferable load resources in reducing the fluctuation of wind and PV power generation.

In this paper, we find that only considering the flexible resource of transferable load has little impact on reducing the fluctuation of wind and PV power generation. Therefore, we will consider more demand-side response policies in future research.

Nomenclature

C:Total cost
C:Cost of the thermal power unit
C:Wind power generation cost
Cpv:PV power generation cost
Ccsp:CSP power station cost
Ch:Hydropower generation cost
Cb:Battery cost
Ctl:Compensation cost of the transferable load
j:Wind farms index
N:Total number of wind farm
Kw,j:Cost coefficient of the j-th wind farm5
Pj,t:Output power of the j-th wind farm at time
k:PV power station index
Npv:Total number of PV power station
l:CSP power station index
Ncsp:Total number of CSP power stations
KCSP,l:Operation cost of the l-th CSP power station
Pl,t:Output power of the l-th CSP power station at time t
KCSP,R,l:Spinning reverse cost of the l-th CSP power station
Pl,t,r:Spinning reverse power of the l-th CSP power station at time t
m:Hydropower station index
Nh:Total number of hydropower stations
Nb:Total number of batteries
KH,m:Cost coefficient of the m-th hydropower station
KB,n:Cost coefficient of the n-th battery
Pm,t:Output power of the m-th hydropower station at time
tts tl:Start time of transferable period
Ttl:Original operation time set of transferable load
Ctl,t:Compensation cost of transferable load at time t
PL,t:Load of power system considering the transferable load
PL0,t:Load of power system not considering the transferable load
Pj,t,f:Forecast power of the j-th wind farm at time t
cn,t:States of charge of the n-th battery at time t
Pn,t:Output power of the n-th battery at time t
Ptl:Transfer power of the transferable load
Pt,tl:Power transferred from time t to tl
Utl:Status of the transferable load
Pmax j:Maximum power of the j-th wind farm
Δtn:Unit time interval
Pj,t,e:Forecast error power of the j-th wind farm at time t
Pmax k:Maximum power of the k-th PV power station
Pk,t,f:Forecast power of the k-th PV power station at time t
Pk,t,e:Forecast error power of the k-th PV power station at time t
Pmax l:Maximum power of the l-th CSP power station
PTES l,t,d:TES discharge power of the l-th CSP power station at time t
PTES l,t,c:TES charge power of the l-th CSP power station at time t
PTES l,t,r+:Positive spinning reserve provided by TES of the l-th CSP power station at time t
PTES l,t,r-:Negative spinning reserve provided by TES of the l-th CSP power station at time t
PTES,max l,t,d:TES maximum discharge power of the l-th CSP power station at time t
PTES,max l,t,c:TES maximum charge power of the l-th CSP power station at time t
RUpv:Up spinning reserve of the power system
RDs:Down spinning reserve of the power system
ηc n:Charge efficiencies
ηd n:Discharge efficiencies
dn,t:States of discharge of the n-th battery at time t
te tl:End time of transferable period
Wh,f:Recursive weight matrixes of ft
Wh,i:Recursive weight matrixes of it
Wh,o:Recursive weight matrixes of ot
Wx,f:Input weight matrixes of ft
Wx,i:Input weight matrixes of it
Wx,o:Input weight matrixes of ot
bf:Offset matrixes of ft
bi:Offset matrixes of it
bo:Offset matrixes of ot
μt:The mean
σt:The variance
t:Free degree
Kpv,k:Cost coefficient of the k-th PV power station
Pk,t:Output power of the k-th PV power station at time t
Yl,t:State variable
Pmax m:Maximum output power of the m-th hydropower station
Pmin m:Minimum output power of the m-th hydropower station
Vm,t:Storage capacity of the reservoir of the m-th hydropower station at time t
Vm,t+1:Storage capacity of the reservoir of the m-th hydropower station at time t+1
Vmax m:Maximum storage capacity of the reservoir of the m-th hydropower station
Vmin m:Minimum storage capacity of the reservoir of the m-th hydropower station
Qm,t:Outflow of the m-th hydropower station at time t
Qmax m:Maximum outflow of the m-th hydropower station
Qmin m:Minimum outflow of the m-th hydropower station
Im,t:Inflow of the m-th hydropower station at time t
o:Upstream hydropower station of the m-th hydropower station
τom :Time delay from the o-th hydropower station to the m-th hydropower station
Qo,(t-τom):Outflow of the o-th hydropower station by time (t-τom)
En,t:Capacities of the n-th battery at time t
En,t-1:Capacities of the n-th battery at time t-1
Emax n:Maximum capacities of the n-th battery
Emin n:Minimum capacities of the n-th battery
Pc n,t:Charge power of the n-th battery at time t
Pd n,t:Discharge power of the n-th battery at time t
Pc,max n:Maximum charge power of the n-th battery
Pd,max n:Maximum discharge power of the n-th battery.
Abbreviations
PV:Photovoltaic
CSP:Concentrating solar power
LSTM:Long short-term memory
CFCEP:Chaotic fast convergence evolutionary programming
LHS:Latin hypercube sampling
GA:Genetic algorithm
SQP:Sequential quadratic programming
AGC:Automatic generation control
TES:Thermal energy storage
RNN:Recurrent neural network
RMSE:Root means square error.

Data Availability

Detailed data can be found within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant no U2003110 and in part by the Key Laboratory Project of Shaanxi Provincial Department of Education (No. : 20JS110) and Shaanxi Youth Science and Technology New Star Project (No. : 2020KJXX-094) and Xi’an Science and Technology Planning Project (No. : 2020KJRC0084).