Low-carbon generation expansion planning considering uncertainty of renewable energy at multi-time scales

0 Introduction

For carbon emission reduction,the installed capacities of wind turbines (WTs)and photovoltaics (PVs)have grown rapidly in China in recent years [1-2].As of the end of 2019,the cumulative installed capacity of renewable energy exceeded 400 million kW,accounting for 34% of the global installed capacity of renewable energy [3-4].Although wind and PV power have low carbon emissions,they pose inherent uncertainty and fluctuation challenges for the lowcarbon planning and operation of a power system [5-6].Therefore,it is necessary to fully consider the uncertainty of wind and PV power in a low-carbon generation expansion planning model.

Considerable research has been conducted on the correlation and fluctuation of renewable energies to capture their characteristics.In [7],a synergy coefficient based on the Pearson linear correlation coefficient was proposed to analyze the correlation between wind and PV power.Based on the probability distributions of wind and PV power,the Copula function was used to describe the spatial correlation in [8] and [9]; however,the temporal correlation was not considered.In [10],a wind power sequence model was presented based on wind speed,considering the probability,autocorrelation,and season characteristics.In [11],a hybrid forecasting method based on a neural network was proposed by analyzing the forecasting methods of wind power and wind speed.In [12],a probability density model of wind power fluctuation was presented considering the statistical laws of wind power fluctuation at multi-time scales.In [13],a stochastic differential equation was used to capture the multi-time dependence structure of wind speed to derive predictive distributions and time-path trajectories.Although the correlation and fluctuation of renewable energies have been studied at certain scales,their uncertainty at multi-time scales still needs further study.

An integrated low-carbon generation expansion planning model was proposed in [14],which integrates and formulates the impacts of various low-carbon factors.A comprehensive generation planning model with a suitably modified objective function and additional constraints was proposed in [15] considering carbon emission trade and carbon tax.In [16] and [17],a generation planning model with carbon capture equipment was proposed,and the model results demonstrated that 90% carbon emission reduction can be achieved at the cost of a small amount of energy consumption.In [18],a generation planning model based on the operating characteristics of carbon capture equipment was presented,which considered carbon emission and primary coal supply constraints.A lowcarbon generation expansion planning model was proposed by introducing the green certificate trade mechanism and the carbon trade mechanism [19] with consideration of the maximum income and renewable energy quota,carbon trade,and number of green certificates.A model considering economic and environmental benefits [20] was proposed based on investment decision optimization and multi-time scale optimization operation,and was found to be suitable for large-scale renewable energy integration.Although the method of low-carbon generation expansion planning was studied in the aforementioned research,the integration of low-carbon elements and the characteristics of renewable energy at multi-time scales should be further considered.For large-scale renewable energy integration in a power system,effect of the uncertainty at multi-time scales on the generation expansion planning cannot be considered simply.

The contributions of this paper are as follows.Aiming to solve the aforementioned problem,a two-layer lowcarbon generation expansion planning model considering the uncertainty of renewable energy at multi-time scales is proposed based on the concept of decomposition coordination.The uncertainty at multi-time scales which is match investment decisions and multi-time scale operation assessment is considered in detail.In the upper layer of the model,the investment decision is optimized considering the uncertainty of renewable energy at a monthly scale.In the lower layer,the optimal operation is simulated to verify the peak regulation ability and flexibilities of the investment decision optimized in the upper layer.The peak regulation ability is verified by daily operation simulation at an hourly scale,whereas the flexibility is verified at a 15-minute scale.The model promotes carbon emission reduction while maintaining flexibility.

The remainder of this paper is organized as follows: an uncertainty model of renewable energy,which considers the spatio-temporal correlation and fluctuation at multi-time scales,is described in section 2.The proposed low-carbon generation expansion planning model comprising investment decision optimization and the operation simulation models is presented in section 3.A case study and simulation results are discussed in section 4,including the results and model comparisons.Finally,section 5 concludes the paper.

1 Uncertainty model of renewable energy at multi-time scales

1.1 Spatio-temporal correlation model of renewable energy

1.1.1 Probability density function of renewable energy

The kernel density estimation model is used in this study to describe wind power.The kernel density estimation is expressed as

where h is the smoothing coefficient (h＞ 0),K(·)is the kernel function,and xi is the ith sample of the variable,and the cumulative distribution function (CDF)can be obtained by integrating f(x).

PV power presents an approximately linear correlation[21] with irradiance.The beta distribution can be used to describe the irradiance in certain periods.Subsequently,the probability density function of PV power can be obtained as follows:

where α and β are the shape parameters of the beta distribution and ymax is the maximum PV power output.

1.1.2 Spatial correlation model of renewable energy

Based on the probability distribution functions of wind and PV power,the Copula function [22] is used to describe the spatial correlation.The modeling steps are as follows:

1) The parameters of the probability distributions expressed in (1)and (2)are obtained by the maximum likelihood estimation.

2)The probability distribution functions obtained in 1)are regarded as the marginal distribution functions of the Copula function.Subsequently,the Kendall correlation coefficient and the Euclidean distance are employed to fit the Copula function indices [23-24].

3)The optimal Copula function obtained in 2)is used to generate the sequences considering the spatial correlation of wind and PV power using the Monte Carlo method.These sequences considering spatial correlation can be expressed asandrespectively,where N is the sequence length.

1.1.3 Spatio-temporal correlation model based on spatial correlation model

The sequences considering spatio-temporal correlation can be generated based on the model of temporal correlation and the sequences considering spatial correlation,which are described in section 2.1.2.

The temporal correlation model is as follows:

Let ht-1 and ht be the outputs of wind or PV power at times t-1 and t, respectively.Thus,the probability distributions of wind and PV power at the corresponding times can be expressed as Ft-1(ht-1)and Ft (ht),respectively.Based on the Copula theory,the joint probability density function at these adjacent times can be expressed as follows:

where c(·)is the probability density function of the optimal copula function,F(ht-1,ht)is the joint probability distribution,and f t-1(ht-1)and f t (ht)are the probability density functions of renewable energy at times t-1 and t,respectively.

Based on the conditional probability distribution theory,the probability density of renewable energy at time t can be obtained using ht-1.

whereis the conditional probability density function of renewable energy at time t.

Based on the temporal correlation model,the wind power sequence considering the spatio-temporal correlation can be generated by the following steps:

1)u1w is substituted for h1 in (4).

2)The wind power sequence at time t + 1 is updated based on the joint probability density in (4).

3)t = t + 1 is set,and step 2)is repeated until t = N.

Finally,the wind power sequence considering the spatiotemporal correlation is expressed aswhere hiw is the wind power output at time i.Similarly,the PV power sequence,is also updated.

1.2 Fluctuation model of renewable energy

The spatio-temporal correlation is considered in the aforementioned generated sequences,whereas the fluctuation is not.The sequences still need to be modified by considering the fluctuation.The fluctuation distributions of renewable energy at multi-time scales are different.The normal distribution [8] has been commonly used to describe the fluctuation of renewable energy.t-location-scale and logistic distributions are introduced in this study to describe the fluctuation of renewable energy,which are more flexible than the normal distribution.

The probability density function of t-location-scale is expressed as

where v,μ,and σ are the shape,position,and scale parameter,respectively.

The probability density function of the logistic can be expressed as

where u1 is the position parameter and σ1 is the scale parameter.

The statistic,S,in the K-S test [25] and fitting error I are introduced as indices to test the performance of fitting.

where I is the fitting error,f is the fitting probability density function, is the height of the ith histogram,yi is the value of the fitting probability density function at position ai,and ai is the center of the ith histogram.

1.3 Sequence considering uncertainty at multitime scales

Based on the data demand of the low-carbon generation expansion planning model proposed in this paper,fluctuation models at multi-time scales are combined with the models considering the spatio-temporal correlation of wind and PV power to generate sequences at multi-time scales.The procedure of generating sequences at multi-time scales is shown in Fig 1.The method of generating wind power sequence is similar to that of PV power sequence.Taking wind power as an example in this section,the method is as follows:

Step 1: Based on the historical data of wind power,the hourly sequence considering the spatio-temporal correlation is generated by the model presented in section 2.1 and denoted as

Step 2: Based on the hourly sequence in step 1 and the fluctuation model at a 15-minute level,the sequence of the 15-minute level is generated.The wind power at period t is htw in an ideal condition when the fluctuation is not considered.

Based on the fluctuation model described in section 2.2,the value of the fluctuation at the 15-minute scale during period t is generated,which is denoted as

where and are the output of the wind power and the value of the wind power fluctuation at the ath 15-minute interval at period t,where a=1,2,3,4.The wind power sequence generated here is used for a 15-minute operation simulation in the generation planning model.

Step 3: Similarly,the hourly sequence of wind power is generated using the hourly fluctuation model to update the hourly sequence considering the spatio-temporal correlation.

where htwr and ΔPtw are the output of the wind power and the value of the wind power fluctuation at time t.The generated sequence here is used for hourly operation simulation on typical days in the generation planning model.

Step 4: Based on the hourly sequence in step 1,the capacity credit [19] is used to describe the monthly sequence.The capacity credit is calculated using the dichotomy method.The monthly sequence is used for investment decision optimization in the generation planning model.

2 Low-carbon generation expansion planning model

2.1 Framework of low-carbon generation expansion planning model

The uncertainty of renewable energy at multi-time scales and low-carbon elements are introduced into the traditional generation expansion planning model.Subsequently,the two-layer low-carbon generation expansion planning model considering the uncertainty of renewable energy at multi-time scales is proposed.Based on the concept of decomposition coordination,the scale of the problem can be simplified and the uncertainty at multi-time scales can be considered more accurately using the two-layer model than that by the traditional model.The upper layer model,i.e.,the investment decision optimization,mainly considers the planning scheme,load growth,and total cost of the generation planning after the introduction of a carbon trading model and carbon capture power plants.The lower layer one,i.e.,the operation simulation,comprises the typical daily operations considering peak regulation verification and a 15-minute operation considering flexibility verification.In addition,typical days are selected by clustering the generated sequences.By coordinated iteration and rolling correction between investment decision optimization and multi-time scale operation simulation,the optimized planning scheme with comprehensive benefits of economy,low-carbon emission,and flexibility is obtained.The overall framework is shown in Fig.2.

Fig.1 Flowchart of generating sequence at multi-time scales

Fig.2 Framework of generation planning model

2.2 Model of investment decision optimization

The investment decision is made monthly to consider the fluctuation and uncertainty of renewable energy accurately.

2.2.1 Objective function of investment decision optimization

where Fiv is the investment cost,

where N is the number of planning months; ζi is the set of candidate power plants; dy is the discount coefficient in month y, and d y=(1+ r/12)1-y,where r is the annual discount interest rate; py is the discount coefficient in month y; Ui is the capital cost; and Pimax is the capacity of power plant i.

Fop is the operating cost of a power plant,comprising its fixed operating cost and the fuel cost.

where ζe is the set of existing power plants,Gi,y is the fixed cost of power plant i in month y,Ei,y is the power generation in month y,and fi is the power generation cost of power plant i.

During the planning stage,the carbon trading cost,Fco2,is expressed as

where Pco2 is the carbon trading cost in month y,

where Ed is the carbon emission quota,which is denoted asis the total power generation of all power plants at time t,η is the coefficient of the carbon emission quota,WCO2 is the carbon trading price,Eh is the upper limit purchase of carbon emissions,and Wpc is the penalty fee for the excess of carbon emission quota.

2.2.2 Investment decision optimization constraints

1)Power balance constraint

wher e Pyi is the available capacity of power plant i in month y,which consists of the installed capacity of coal-fired power plants,the expected output of hydropower power plants,and monthly capacity credit for renewable energy considering multi-time scales.Lmyax is the maximum load in month y and RyD is the coefficient of the reserve.

2) Electricity balance constraint

where Eym is the demand of electricity in month y and Rye is the coefficient of electricity reserve in month y.

3)Primary energy supply constraint

The introduction of carbon capture technology may exacerbate the primary energy shortage because of the additional energy consumption.Therefore,the upper limit of the primary energy supply is set as

where qi is the fuel consumption characteristics of fuel m and Qi,y is the total supply of fuel m in month y.Primary coal resources are considered in this study because of inclusion of a coal-based system.

4)Carbon trade constraint

where is the carbon emission coefficient (t/MWh)of power plant i,Erye is the actual carbon emission, is the carbon emissions quota,Ehy denotes the carbon emission purchases,and is the carbon emissions sale in month y.

5)Penetration of renewable energy constraint

where is the increased installed capacity of renewable energy power plant z in month y and λy is the penetration coefficient of renewable energy installed capacity in month y.

6)Negative reserve constraint

where θi is the negative reserve coefficient of power plant I and Rydown is the coefficient of the negative reserve demand in month y.

2.3 Model of operation simulation

The operation simulation is used to verify the peak regulation and ramping abilities at hourly and 15-minute scales on typical days.

2.3.1 Hourly operation simulation on typical days

The objective function comprises the operating costs,penalty fees for curtailment,penalty fees for the loss of load,and carbon trading costs.

where T is the total number of operating hours on typical days,Ωg is the set of all operable power plants,Pi,t is the power output of power plant i at time t,ρres is the penalty fee of renewable energy curtailment,Eres,t is the renewable energy curtailment at time t,ρloss is the penalty fee of the loss of load,Eloss,t is the loss of load of the system,and Pco,2t is the carbon trading cost of the system.

The constraints on typical daily hourly operation simulation are as follows:

1) Coal-fired power plants constraints

where Pi,min is the minimum output of coal-fired power plant i,Pi,max is the maximum output,Ii,t is the binary state status(1 for running status and 0 otherwise)at time t,ΩThe is the set of coal-fired power plants,Ti,on is the minimum uptime,Ti,down is the minimum downtime; Ui,t is the start-up binary variable,and Di,t is the binary shutdown variable.The gasfired and nuclear power plants are similar to the coal-fired power plants.The aforementioned constraints comprise the upper and lower output limits of the coal-fired power plants,minimum up-down time constraints,and logic constraints of the statuses.

2) Renewable energy constraint

where Pi,t is the available output of a renewable energy power plant and Ωres is the set of renewable energy power plants.

3)Positive and negative reserve constraints

where Rtu and Rtd are the posi tive and negative reserves of the system at time t respectively,and Lt is the load at time t.

4) Power balance constraint

2.3.2 Fifteen-minute operation simulation on typical days

The schedule of each power plant is obtained from the hourly operation simulation on typical days,as described in section 2.3.1.

The objective function of the 15-minute operation simulation comprises the operating costs,penalty fees for renewable energy curtailment,and penalty fees for insufficient ramping capacity.

where Tm is the total number of operating periods,ρramp is the penalty fee for insufficient ramping capacity,and and are the relaxation of the upward and downward ramping capacities at time t,respectively.

Note that the constraints are similar to those of the hourly operation simulation presented in section 1.3.1.The added 15-minute ramping constraints of the system are as follows:

where and are the upward and downward ramping rates,respectively.

2.4 Solving procedure

The proposed generation expansion planning model proposed in this paper is a mixed-integer linear programming model.Solver Cplex or Gurobi can be used to solve the model.The solving procedure is shown in Fig.3.

Fig.3 Flowchart of solving procedure

3 Case study

The generation expansion models are implemented in MATLAB using solver CPLEX 12.5 on a PC with an Intel Core(TM)i7-8750H CPU and 8 GB RAM,and the model is established using YALMIP.

3.1 Test system and data

The proposed low-carbon generation expansion planning model is applied to a provincial power grid covering 2020 to 2026.The installed capacity in 2020,in this case,is shown in Fig.4.

Fig.4 Installed capacity in 2020

The full life cycle carbon emissions of wind and PV power are considered in this study.The carbon emissions of wind and PV power are 0.252 and 0.246 kg/kWh [26-28],respectively.The overall carbon emission of coal-fired power is 0.89 kg/kWh.

The main operating parameters of the different types of power plants are listed in Table 1.

Table 1 Operating parameters of power plants

Carbon emission(t/MWh)Types Generation cost(yuan/MWh)Coal cost(t/MWh)Coal-fired 350 0.89 0.34 Carbon Capture 410 0.09 0.40 Gas-fired 610 0.4 0 Nuclear 84 0 0 Pump 105 0 0 Hydro 0 0 0 WT 0 0.252 0 PV 0 0.246 0

The data of the candidate plants are summarized in Table 2.

Table 2 Data of candidate plants

Total capacity(MW)Coal-fired 375 100 6400 Types Investment(million yuan/MW)Fixed cost(million yuan/year)Carbon Capture 393 115 12600 Gas-fired 315 113 7650 Nuclear 1150 450 8150 Pump 440 205 1800 Hydro 1136 295 480 WT 710 220 11220 PV 782 240 12100

The carbon trading price is set as 150 yuan/ton.The carbon emission quota is set as 0.75 ton/MWh.RyD is set as 10%.Eh is set as 0.4 times the carbon emission quota.The discount rate is set as 10%.Wpc is set as 600 yuan/ton.λ is set as 30%.The threshold of peak regulation is set as 0.01 times the value of load in the corresponding period.The maintenance of the plants is arranged based on the maintenance schedule of the province.

3.2 Scenarios generation at different scales

The hourly fluctuation fitting indices of wind power and photovoltaic are shown in Table 3.

Table 3 Fitting indexes of hourly fluctuation

Types index normal t-location-scale logistic WT S 0.0832 0.0197 0.0461 I 0.0051 0.0016 0.0031 PV S 0.1733 0.1467 0.1550 I 0.1060 0.079 0.1024

Based on Table 3 the t-location-scale distribution is more accurate than the normal and logistic distributions to fit the hourly fluctuation.

Furthermore,the renewable energy fluctuations at other time scales are analyzed.The results demonstrate that the t-location-scale distribution is more suitable to fit the hourly fluctuation,whereas the logistic distribution is more accurate to fit the fluctuation at the 15-minute scale.

Using the historical data of 2019,the fitting results of different Copula functions are provided in Table 4.

Table 4 Parameters of candidate Copula functions

Euclidean distance Gaussian 0.9796 0.8713 0.0057 t 0.9868 0.8964 0.0042 Gumbel 9.0747 0.8988 0.0041 Clayton 11.5878 0.8528 0.0099 Frank 35.7888 0.9134 0.0026 Types Model coefficient Kendall coefficient

From Table 4 the optimal Copula function is Frank-Copula because of the maximum Kendall coefficient and the minimum Euclidean distance.

Using the proposed simulation method considering the uncertainty of renewable energy at multi-time scales and the provincial historical data,wind and PV power sequences are generated at the 15-minute and hourly scales.

The generated sequences are validated by the verification of the spatial correlation,autocorrelation,and fluctuation.Taking the generated sequences in a week at an hourly time scale as an example,the spatial correlation coefficients for the original and generated sequences are -0.3962 and-0.3621,respectively.Using the coefficients,the proposed method can capture the spatial correlation accurately.The verification of autocorrelation and fluctuation are shown in Fig.5 and 6,respectively.The figures show that the autocorrelation function (ACF)and CDF for the generated sequence are similar to those for the original sequence.The proposed method can also capture the autocorrelation and fluctuation accurately.

Fig.5 ACF comparison for renewable energy sequence at hour scale

Fig.6 CDF for renewable energy fluctuation at hour scale

3.3 Simulation results and analysis

The investment decision results of the proposed model and their costs are shown in Fig.7 and listed in Table 5 respectively.

Fig.7 Results of investment decision optimization

Table 5 Costs of the expansion plan using the proposed model

Investment(billion yuan)Coal demand(billion tons)152.79 854.2 -70.75 3.641 1.37 Operating cost (billion yuan)Carbon trad cost(billion yuan)Carbon emission(billion tons)

The added installed capacities of the coal-fired,carbon capture,gas-fired,hydro,pump,nuclear,wind power,and PV power are 3670 MW,8400 MW,5390 MW,480 MW,1800 MW,8150 MW,and 15900 MW respectively.The power plants with low generation cost and low carbon emissions,i.e.,pump,hydro,renewable energy,and nuclear,are placed into production in the early stages.Owing to the limitation of the primary coal supply,carbon capture power plants with low carbon emissions are placed into production in large quantities until 2024 to deal with the increase in the load and system carbon emissions.At the end of the planning period,the installed capacity of the gas-fired plant gradually increases to alleviate the burden of the primary coal supply and the carbon emissions.

In addition,the effects of the introduction of carbon capture power plants and large-scale renewable energy integration on peak regulation are considered by hourly operation simulation of typical days.Taking the day of maximum load in 2026 as an example,the results of the simulation are shown in Fig.8.

Fig.8 Hourly simulation of typical daily

Based on Fig.8 the nuclear power plants,coal-fired power plants with low generation costs,and carbon capture power plants with better peak regulation performance have priority for power generation.The proposed generation planning scheme has sufficient capacity for peak regulation and non-curtailment,which fully accommodates wind and PV power and promotes carbon emission reduction.

3.4 Comparison with traditional model

The proposed model in this paper is compared with the traditional model without considering the uncertainty of renewable energy at multi-time scales and the carbon capture power plants.The comparison results of the installed capacity in 2026 are shown in Fig.9.

Fig.9 Comparison of installed capacity in 2026

The key indices of the two low-carbon generation expansion planning schemes are listed in Table 6.

Table 6 Comparisons of key indices

Total cost(billion yuan)Proposed 152.79 854.24 -70.75 3.641 936.28 Traditional 140.13 862.68 -55.21 3.746 947.66 Model Investment cost(billion yuan)Operation cost (billion yuan)Carbon trading cost(billion yuan)Carbon emissions(billion tons)

The proposed model performs better than the traditional model.The installed capacity of the coal-fired power plants is significantly reduced and the carbon capture power plants with low carbon emissions are placed into production on priority.Based on Table 6 compared to the traditional model,the proposed scheme reduces the total cost by 13.24 billion yuan and the carbon emissions by 82 million ton,which ensures the economic benefits of the system while achieving carbon emission reduction.

3.5 Sensitivity analysis of proposed model

1)Effects of carbon emission quota on planning results

Based on the proposed scheme,the carbon emissions quota is increased by 10% and decreased by 10% as comparison cases.The low-carbon generation expansion planning schemes in 2026 under different carbon emission quotas are shown in Fig.10.

Fig.10 Planning results with different carbon emission quotas

From Fig 10,the carbon emission quota has a significant impact on the installed capacities of the coal- and gasfired plants in the planning period.The carbon-free power plants,such as nuclear,pump,and hydro plants,are placed into production on priority to alleviate the burden of carbon emissions under the condition of a low-carbon emission quota.Large-scale coal-fired power plants with low generation costs are placed into production to ensure economic benefits when the carbon emission quota is high.

2)Effects of carbon trading price on planning results

Based on the proposed scheme,the carbon trading price is increased by 20% and decreased by 20% as comparison cases.The low-carbon generation expansion planning schemes under the different carbon emission quotas are shown in Fig.11.

Fig.11 Planning results with different carbon trading prices

Carbon trading price has a significant effect on coalfired and carbon capture plants.The coal-fired power plants are placed into production on priority with lowcarbon trading prices,whereas the carbon capture and gasfired power plants with higher generation costs have lower installed capacities.With the increase in the carbon trading price,the carbon capture,gas-fired,and nuclear power plants are placed into production on priority to alleviate the carbon emission burden.

3)Effects of penetration of renewable energy on planning results

Based on the proposed scheme,the penetration of renewable energy is increased by 10% and decreased by 10% as comparison cases.the low-carbon generation expansion planning schemes under the different proportions of renewable energy are shown in Fig.12.

Fig.12 Planning results with different proportions

With the increase in the renewable energy penetration,system tends toward low carbonization.The carbon capture,gas-fired,and nuclear power plants with low carbon emissions are placed into production,whereas the proportion of coal-fired power plants with high carbon emissions is increasingly reduced.

4 Conclusion

In this paper,a two-layer low-carbon generation expansion planning approach considering the uncertainty of renewable energy at multi-time scales is proposed.Based on the comparison and sensitivity analysis,the following conclusions are obtained:

1)A low-carbonization power structure is achieved by considering the life cycle carbon emissions of renewable energy and introducing low-carbon technology and market.

2)Compared with the traditional model,the uncertainty of renewable energy at multi-time scales and carbon capture power plants are introduced into the model proposed in this paper.The renewable energy outputs are simulated more accurately,and the operation costs and carbon emissions are reduced correspondingly.

3)The installed capacity of the coal-fired power plants with high carbon emissions is decreased when the carbon quota is low and the carbon trading price is high.in comparison,the power plants with low-carbon emissions,such as carbon capture,nuclear,and renewable energy plant,are placed into production on priority.

Acknowledgements

This study was supported partly by the National Key R＆D Program of China (2018YFA0702200)and the Science and Technology Project of State Grid Shandong Electric Power Company (520604190002).

Declaration of Competing Interest

We declare that we have no conflict of interest.