Power System Economic Dispatch with Spatio-temporal Wind Forecasts

Transcription

1 Power System Economic Dispatch with Spatio-temporal Wind Forecasts Le Xie*, Yingzhong Gu*, Xinxin Zhu t, and Marc G. Genton t *Department of Electrical and Computer Engineering, Texas A&M University, College Station, TX, USA Lxie@ece.tamu.edu, gyzdmgqy@tamu.edu t Department of Statistics, Texas A&M University, College Station, TX, USA {xzhu,genton}@stat.tamu.edu Abstract-In this paper spatio-temporal wind forecast is incorporated in power system economic dispatch models. Compared to most existing power system dispatch models, the proposed formulation takes into account both spatial and temporal wind power correlations. This in turn leads to an overall more costeffective scheduling of system-wide wind generation portfolios. The potential economic benefits are manifested in the systemwide generation cost savings, as well as the ancillary service cost savings. We illustrate in a modified IEEE 24 bus system that the overall generation cost can be reduced by 12.7% by using spatiotemporal wind forecasts compared with only using a persistent forecast model. I. INTRODUCTION This paper is motivated by the fact that renewable wind energy is increasingly being integrated into electric power systems throughout the world [1] [2]. By 21 several regions in the US as well as around the globe have already achieved more than 1% wind installed capacity [3] [4]. Compared to conventional fossil fuel generation, wind generation exhibits limited predictability and high inter-temporal variations (i.e ramping) [5]. In order to fully extract the environmental and economic potential of renewable energy resources, novel power system operating paradigm is needed to reduce the balancing and ancillary service burden on by fast responsive fossil fuel units [6]. The basic objective of power system scheduling is to maintain the supply and demand balance at a minimum cost subject to transmission constraints and plausible contingencies. Before high penetration of renewable energy such as wind and solar, the uncertainty in power system scheduling primarily comes from the demand side [5]. With the high presence of intermittent wind power, the uncertainty now comes from both demand and supply (wind and solar) sides. Load forecasting has been an active area of research for more than four decades. State-of-the-art load forecasts could achieve high accuracy in the day-ahead stage. In contrast to load forecasting, variable wind generation is very hard to be accurately predicted at hours ahead, which is typically based on the physical weather model [7]. Only near term (lo-minute to 2 hours ahead) prediction is of acceptable confidence. The short-term forecast is usually based on statistical models, a good set of references can be found in [8]-[11]. The main subject of this paper is to exploit a novel economic dispatch model by leveraging short-term spatio- temporal wind forecast information for regions with high wind penetration. Wind generation is driven by wind patterns, which tend to follow certain geographical spatial correlations. For large regions of wind farms, the wind generation forecast in the downstream of wind could significantly benefit from the upstream wind power generation. Enabled by the technological advances in communications, spatially correlated wind data could be leveraged for higher system-wide short-term wind forecasts. This is potentially very applicable to large-scale wind farms, e.g offshore wind generation. Starting from our recent work on look-ahead dispatch, which utilizes inter-temporal wind forecasts in economic dispatch [15], we proposed in this paper a first attempt to incorporate both spatial and temporal wind forecasts in power system dispatch. In summary, the main contribution of this paper is twofold: 1) We propose to use a spatio-temporal correlated forecast model for short-term wind generation in a region. This forecast model takes into account both local and nearby wind farms' historical data. This method is tested with realistic wind data, and has shown improved forecast accuracy. 2) We formulate an economic dispatch model which incorporates the spatio-temporal wind forecast information. Numerical study in a modified IEEE 24-bus test system shows the improved economic benefits compared to both static economic dispatch and look-ahead dynamic dispatch. This paper is organized as follows. In Section II we propose an economic dispatch model which incorporates available short-term spatio-temporal wind power forecast data. The detailed wind forecast method is discussed in Section III, in which both spatial and temporal historical wind data is utilized for improved forecast accuracy. An illustrative power system economic dispatch example is presented in Section IV, which quantifies the potential savings in both generation cost and ancillary services in the proposed dispatch model. Concluding remarks and future work are discussed in Section V. II. ECONOMIC DISPATCH MODEL WITH SPATIO-TEMPORAL WIND FOREC AST In this section we formulate the power system economic dispatch problem considering both spatially and temporally /11/$ IEEE

2 G D W T n CGi CWi CRi G D p t P Q' Fk Fmax w.6.t prmp Pa.in, Pa ax Pw in, p w ax w Pw Pw TABLE I NOTATIONS Set of conventional power plants Set of inelastic loads Set of wind farms MPC window size (look-ahead horizon) number of wind farms installed in system Generation cost function of power plant i Generation cost function of wind farm i Reserve cost function of power plant i Scheduled generation of power plant i at time k Scheduled generation of wind farm i at time k Forecasted load level of bus i at time k Scheduled reserve capacity of power plant i at time k Scheduled regulation capacity of power plant i at time k Vector of branch flow at time k Vector of transmission line limits Energy Market scheduling interval Ramping constraints of power plant i Minimax generation limit for unit i Minimax generation limit for wind farm i Forecasted wind availability for wind farm i at time k The vector of forecasted wind availability The vector of historical wind data correlated wind forecast. The notations are summarized in Table I. With spatio-temporal wind forecast, the economic dispatch model is presented as follows: s.t. T min: L [ LCGi(P8J + L CWi(PJ kko iew + L CR8i(P8J + CR9i(P9J] iew ied LP9i Rg(aD,aw,PLS) (3) IF k l Fmax, k ko,..., T (4) 1P8i - pll pi R D..T, i E G U W, k ko,..., T (5) pe; i n P8i Pe;ax, k ko,..., T (6) o pti Pe;ax, k ko,..., T (7) o P9i Pe;ax, k ko,..., T (8) pe; i n P8i + pti + P9i Pe;ax, k ko,.. T (9) Pwi i n Pi PWiax, k ko,.. T (1) PWi k Ak PWi' k ko,.. T (11) Pw f(pw). (12) In the proposed formulation, the objective function is to minimize the power system operating costs which include costs of generation and costs of providing reserve services. Constraints of this problem are system and individual units operating constraints for security and reliability purposes. (1) is the energy balance equations, requiring the generated energy to be always equal to the energy consumed in steady state. (2) is the system reserve requirement. (3) is the system regulation requirement balancing equation. Given the system reliability requirement, the wind forecast accuracy and load forecast accuracy, the functions used to determine the reserve requirement and regulation requirement of the system are assumed to be linear. (4) is the transmission line capacity limitation, which contributes to network transmission congestion. (5) is the ramping constraints of units. (6) is the upper bound and lower bound of conventional generators' output. (7) is the available reserve capacity constraints. (8) is the available regulation capacity. (9) is the capacity constraints of each generator for providing energy reserve and regulation services. (1) is the mechanical upper and lower constraints of wind farms' output. In our study, wind resources are assumed to not participate into ancillary services market such as reserve and regulation markets. (11) is the wind forecast for each wind farm, which is regarded as the nature availability for each wind farm at specific times. (12) is the spatial and temporal correlation between historical wind data and wind forecast data, which will be explained in detail in section III. III. SPATIO-TEMPORAL WIND SPEED FORECASTING In this section, we focus on wind speed forecasting, which is a hot topic due to the large scale penetration of wind energy. In [11], some statistical models for short-term wind speed forecasting were reviewed, including conventional time series forecasting models and spatio-temporal forecasting models. In our case, we illustrate equation (12) in detail with a wind speed dataset from three locations in the northwest of the U.S.. Based on the dataset, the short-term future wind speed at each location is forecast with two models: persistent forecasting (PSS) and a trigonometric direction diurnal (TDD) model [9] which takes the spatio-temporal information into the forecasting model. And then wind forecasts are converted into wind power with a certain power curve. Specifically, 1 minutes average wind speed and wind direction data from 22 and 23 were collected from three meteorological towers along the Columbia River Gorge at the boundary of Washington (WA) to the south and Oregon (OR) to the north in Northwest U.S.: Vansycle (OR), Kennewick (WA) and Goodnoe Hills (WA). Vansycle is close to the Stateline wind energy center; Kennewick is 39 kin northwest of Vansycle; and Goodnoe Hills is located 146 km west of Vansycle. Missing data were imputed with the method mentioned in [17]. Wind speed in this area mainly blows from west or east due to its special geographic features. The Columbia River Gorge runs from east to west passaging through the Cascade Mountains, with high terrain in the north and south and this condition results in west wind or east wind. Therefore, wind

3 speed at these three locations are correlated with each other, and to predict near future wind speed at one of the three locations, spatial correlation should be considered in addition to correlation with historical wind speed observations of itself. This can be seen in Fig. I in which auto-correlations and cross-correlations or spatial correlations of the three locations up to time lag 4 min. with 23 data are displayed. Autocorrelations at each location are strong. Auto-correlation coefficients are larger than.7 up to 4 min. lag. And crosscorrelation coefficients are also significantly large. All the values are around.6 up to 4 min. lag. First, we aim at forecasting k min. ahead wind speed at Vansycie, Kennewick and Goodnoe Hills respectively, based on the 1 min. average wind speed and wind direction data, with k 1,...,12, or 1 min. to 2 hours, and a time horizon 24 hours is of interest. In the following, persistent forecasting and a most recent space-time forecasting model will be applied after being briefly introduced Lag Lag Lag,.,q K, V K, G Lag Lag Lag G, K Lag Lag Lag Fig. I. Auto-correlations and cross-correlations of 1 min. wind speed at Vansycle (V), Kennewick (K) and Goodnoe Hills (G). A. Persistent Forecasting PSS assumes the future wind speed is the same as the current one. For example, if vt is wind speed at time t at Vansycie, then k min. future wind speed is predicted as vt, or Vt+k vt. PSS works very well for short-term forecasting, such as 1 min. due to the nature of wind. PSS is usually taken as a reference and an advanced forecasting model is thought to be good if it outperforms PSS. B. Trigonometric Direction Diurnal Model The TDD model proposed by [9] is one of the most advanced current space-time statistical forecasting models which generalizes the Regime-Switching Space-Time model [17] by taking wind direction into the model. We describe the TDD model with 1 min. ahead wind speed forecasting at Vansycie as an example Let vt, Kt, and Ct be 1 min. average wind speed at Vansycie, Kennewick and Goodnoe Hills at time t, and forecast vt+l, or 1 step ahead (1 min. ahead) wind speed. The TDD model assumes vt+l follows a truncated normal distribution on the nonnegative real domain, that is vt+l rv N + (f-lt+l, at+l) (this can be detected by the density plots in Fig. 2), with center parameter f-lt+l and scale parameter at+l. The key lies in modeling these two parameters appropriately. Vansycle Kennewick Goodnoe Hills Z, " 3, " 3, " 3 Fig. 2. Density plot of 1 min. wind speed at Vansycle (V), Kennewick (K) and Goodnoe Hills (G). The center parameter f-lt+l is modeled as f-lt+l Dt+1 + f-lr+l' where Dt+l is made of trigonometric functions to fit the diurnal pattern of wind speed. Specifically, Ds do + d1 sin (2S) + d2 cos (2S) + d3 sin (:S) + d4 cos (:S), where s 1,2,...,24. The residual series after removing the diurnal pattern, f-lr+l' is modeled as a linear function of current and past (up to time lag h) wind speed and trigonometric functions of wind direction at Vansycie, as well as Kennewick and Goodnoe Hills as follows: f-lr+l f(vt,..., Vt;': h,k[,...,k[_h' K[,...,Cr-h, cos(ev,t), sin(ev,t), cos(ev,t_h), sin(ev,t_h) cos(ek,t), sin(ek,t),, cos(ek,t_h), sin(ek,t_h) cos(ec,t), sin(ec,t),, cos(ec,t_h), sin(ec,t_h) ) ' (13) where ev,t is wind direction at Vansycie at time t and Vt, K[, Cr, ev t' ex t' and ec t are residual series of wind speed and wi d dir ction at th three locations without diurnal component. The scale parameter at+l is modeled as where bo, b1 > and Vt is the volatility value: 1 v -( '" ((v:r _ v:r ) 2 + (Kr _ Kr ) 2 io t - 6 t-i t-i-l t-i t-i-l ) 1/ 2 + (Cr-i - Cr_i_l) 2 ). (14)

4 The coefficients in (13) along with bo, b1 in (14) are estimated by the continuous ranked probability score (CRPS) method, see [18] for more details. Predictors in (13) are selected with the Bayesian Information Criteria (BIC), see [9] for details. This linear model is trained by data in 22 with a 45 days sliding window, and for different k values, the models' structure are different, since different predictors would be selected into the linear equation. For example, for k 5, table II displays correlation coefficients between 5 min. ahead wind speed at Vansycle and other variables at current time and up to past 4 min.. Numbers in bold indicate variables that are selected into the linear equation (13) for k 5 at Vansycle. A. Simulation Platform Setup The power grid of the Columbia River Basin where Vansycle (OR), Kennewick (WA) and Goodnoe Hills (WA) are located, are operated by Bonneville Power Administration (BPA). BPA is a self-funded non-profit agency which markets wholesale electrical power in the northwest area [13]. Modified from the IEEE Reliability Test System (RTS-24), the simulation system of 4 MW total installed capacity is configured according to the practical characteristics of BPA. The simulation system network is presented in Fig. 3. TABLE II CORRELATION COEFFICIENTS Variable V cos(ov) sin(ov) K COS(OK) sin(ok) G cos(oa) sin(oa) t t-i t-2 t t So far, the predictive distribution N + (ILt+!, O't+l) is estimated, and here we take the median of the truncated normal distribution as wind speed forecast, which is defined as Z2 ILt+! + O't+1 <I>-1 [1/2. + (1/2)<I>(-ILt+1)/O't+1], where <I>(.) is the cumulative distribution function of a standard normal distribution. C. Results Wind speed forecasting results of the TDD model are transfered to wind power based on a particular manufacturer's power curve [12]. The results are compared with PSS based on the Mean Absolute Error (MAE) on 23 data. For example, each hourly average wind speed during one day at Vansycle is predicted with 6 step ahead forecasting based on the TDD model and PSS, denoted by ft, t 6,12,...,144 and lit-i, while the true values are lit. Then the MAE for the TDD model is (lit - ft) 2 and (lit - lit_l) 2 for PSS. Forecasting results for 1 and 2 hour ahead are compared in table III for July 7, 23. TABLE III MEAN ABSOLUTE ERROR OF WIND POWER FORECAS TS Locations Vansycle Kennewick Goodnoe Hills Models PSS TDD PSS TDD PSS TDD 1 hr ahead hr ahead IV. ILLUSTRATI VE EXAMPLE In this section, details of simulation platform setup, illustrative examples, simulation results and analysis are provided. Fig. 3. System Network Diagram The generators are assigned as different kinds of power sources including hydropower, coal power, wind power, natural gas power and nuclear power. The generator capacity portfolio (installed capacity percentage of different technologies) is configured according to the generation portfolio of the practical BPA system [13]. The load profile is scaled from the historical load profile of the BPA system [14]. Simulation duration is 24 hours. Both energy balancing market and ancillary market including reserve market and regulation market are simulated. The market operating interval is set to be 1 minutes for both energy balancing market and ancillary market. Wind profile of 24 hours are collected and scaled from BPA System [17]. As is shown in Fig. 4, wind generation potentials at the three locations are respectively calculated according to the manufacturer's power curve. Wind energy which is over the maximum generation capability of wind turbines has been curtailed for security purpose. One hour ahead forecast data is used as the inputs of the established simulation system. The size of the MPC optimization window T is 6 (I hour). The generator parameters are factored out and modified according to [16]. Ramping rates and marginal costs are applied as is shown in Table IV. Bus number, type, capacity (Cap), marginal cost (MC), and ramp rate (RP) of each generator are listed. The computation environment is in Matlab 29 under a PC Windows 7 system. B. Market Results Analysis According to the results of the simulation system, the TDD model can increase the actual wind resources utilization, reduce system-wide generation cost, system ancillary services (including regulation and reserve services) cost, and the total system operating cost.

5 12 i 1 :! c 8 J!! Q. c 6 :;:: I!! c 4 Cl." C 2 c I!! c Cl." C iii 4 c( 2 I--TD PS&-.:... AvailabilitYI Time Sample (1 Omins) Time Samle (1 mins) Fig. 4. The Wind Generation Potential Fig. 5. Actual Wind Generation at Goodnoe Hills Bus I TABLE IV GENERATORS CONFIGURATION Type Cap (MW) MC ($/MWh) Hydro 4 6 Coal 2 37 Coal Wind(GH) \ 3 Nuclear Hydro 7 5 Hydro Natural Gas 5 79 Hydro Wind(VS) 11 2 Wind(KW) 8 I RP (p.u.lmin) The results of the electrical power system operating and power market are listed in Table V. TABLE V SYS TEM OPERATING RES ULTS TOO PSS PSSITOO Total system cost $1,94,235 $1,233, % Generation cost $526,39 $535, % Reserve cost $278,4 12 $336, % Regulation cost $289,5 12 $36 1, % Ancillary Services $567,925 $698, % Total wind generation 3,98 MWh 2,764 MWh 89.2% It can be observed in Table V that the system-wide wind generation using the TDD model is 12.1 % higher than using the PSS model. Given the same wind pressure pattern and system load pattern, the TDD model enables a higher wind resources utilization and a higher wind generation ratio than other technologies. This is because the increased accuracy of the TDD model decreases the wind resources which are going to be wasted by underestimation of wind generation potential in forecast. And therefore, a lower generation cost (1.7%) can be achieved with more wind resources utilized than conventional resources such as coal and natural gas. In Fig. 5, the actual wind generation output at Goodnoe Hills is presented. The MAE value of actual wind generation output by using the TDD model is 22.4 while the MAE value of actual wind generation output by using the PSS model is A lower MAE value indicate a lower forecast error (or average wind forecast deviation). And hence, wind forecast produced by the TDD model is more accurate than the PSS model. It can also be observed that due to the better wind forecast quality of the TDD model, less underestimation of wind potential is introduced. The energy generated by a wind farm in Goodnoe Hills by using the TDD model is much higher (44%) than the wind energy generated by using the PSS model. Besides, the improvement of wind forecast quality in the TDD model will also reduce the uncertainty of wind generation scheduling. A lower mean absolute error (MAE), or average wind deviation between forecast value and actual value, applies on the simulation day for the TDD model (21.), which is compared with the PSS model (229.2). It can be verified in Table V that the lower uncertainty of wind generation scheduling costs less in power system ancillary services (decreased by 23.%) including reserve (decreased by 2.8%) and regulation services (decreased by 25.%). The system overall reserve requirement takes account of uncertainty in both wind generation uncertainty (or wind forecast errors) and load level uncertainty (or load forecast errors). The selected reserve capacity is used to compensate the energy imbalance within time frame of half an hour to 2 hours. In Fig. 6, the total system reserve requirement is compared between using the TDD model for wind forecast and the PSS model for wind forecast. It can be observed that by using the TDD model, due to the improved forecast accuracy, the overall reserve requirement can be reduced. Regulation services also help to compensate the energy imbalance of the system in order to keep the system frequency within a secure range. Different from reserve services, regulation capacity is used to smooth out the short-time (l min. to 1 mins) frequency fluctuation and energy imbalance. It can

6 45 FPsSI C :; <T ii i"!. C " E i!? -:; g- " : c 1; -:; CI " :... O r--r- -;-- Fig. 6. o Time Sample (1 mins) Total System Reserve Service Requirement Fig. 7. Time Sample (1 mins) Total System Regulation Service Requirement be seen from Fig. 7, the TDD model can decrease the system requirement for regulation services and therefore reduce the corresponding regulation cost. Because of the reduction in both energy balancing market (total generation cost) and ancillary services market by using the TDD model, the total system operating cost has been reduced by 12.7% compared with using the PSS model V. CONCLUSIONS In this paper, we incorporated a spatio-temporal wind forecast model (TDD model) in power system economic dispatch. Compared with conventional temporal-only wind forecast models such as the PSS model, the TDD model considers both the local and geographically correlated wind farms for wind forecast. By leveraging both temporal and spatial wind historical data, more accurate wind forecasts can be obtained. We illustrated the potential economic benefits of more accurate wind forecast in a modified IEEE RTS 24 bus system. It was observed that the TDD model can increase the wind resources utilization, and reduce the system costs in both energy balancing and ancillary services. In our future work, we plan to investigate the potential applicability of the proposed dispatch model to large-scale wind farms such as offshore wind farms. Given the more consistent wind pattern over larger geographical area, the potential benefits of the proposed method could be higher. Another important direction for future research is to analyze the tradeoff between communication/computation burden and the improved economic benefits by incorporating more spatially correlated wind data into power system dispatch models. ACKNOWLEDGEMENTS The work of the first two authors was supported in part by Texas Engineering Experiment Station and National Science Foundation ECCS Grant # The work of the last two authors was supported in part by NSF Grant DMS REFERENCES [I) U.S. Department of Energy, 2% Wind Energy by 23: Increasing Wind Energy's Contribution to US. Electricity Supply, May 28. [2) M. B. McElroy, X. Lu, C. P. Nielsen, and Y. Wang, "Potential for windgenerated electricity in China," Science, Vol. 325, No. 5946, pp , Sep 29. [3) Electric Reliability Council of Texas, ERCOT Quick Facts, July 21. [4) H. Weigt, "Germanys Wind Energy: The Potential for Fossil Capacity Replacement and Cost Saving," Applied Energy, Vol. 86, Issue 1, pp , October 29. [5) L. Xie, P.M.S. Carvalho, L.A.F.M. Ferreira, J. Liu, B.H. Krogh, N. Popli, and M.D. llie, "Wind energy integration in power systems: operational challenges and possible solutions," Special Issue of the Proceedings of IEEE on Network Systems Engineering for Meeting the Energy and Environment Dream, Vol. 99, No. I, pp , January 211. [6) Y. V. Makarov, C. Loutan, J. Ma and P. de Mello, "Operational impacts of wind generation on California power systems," IEEE Transactions on Power Systems, Vol. 24, Issue 2, pp , May 29. [7) IBM, Deep Thunder-P recision Forecasting for Weather-Sensitive Business Operations, available online at [8) C. Monteiro, R. Bessa, V. Miranda, A. Botterud, J. Wang, and G. Conzelmann, Wind Power Forecasting: State-of-the-Art 29, Argonne National Laboratory, US Department of Energy, 29. [9) A.S. Hering, and M.G. Genton, "Powering up with space-time wind forecasting," Journal of the American Statistical Association, vol. 15, pp , 2 1. [1) N. Abdel-Karim and M.D. Ilie, "Short Term Wind Speed Prediction by Finite and Infinite Impulse Response Filters: A State Space Model Representation Using Discrete Markov Process," Proceedings of IEEE PowerTech Conference, Romania June 29. [II) X. Zhu and M.G. Genton, "Short-term wind speed forecasting for power system operations," IAMCS technical report, [12) GE Energy, 1.5 MW Wind Turbine, available online at [13) Bonneville Power Administration, "29 BPA facts," Portland, Novermber 21. [14) Bonneville Power Administration, "BPA total transmission system 1ad(TTSL)," Portland, July 27. [15) L. Xie and M. llie, "Model Predictive EconomiclEnvironmental Dispatch of Power Systems with Intermittent Resources," IEEE Power and Energy Society General Meeting, Calgary, July 29. [16) Y. Gu, and L. Xie, "Look-ahead Coordination of Wind Energy and Electric Vehicles: A Market-based Approach," North American Power Symposium 21, The University of Texas, Arlington, Sept.21O. [17) T. Gneiting, K. Larson, K. Westrick, M.G. Genton, and E. Aldrich, "Calibrated probabilistic forecasting at the Stateline wind energy center: The regime-switching space-time method," Journal of the American Statistical Association, vol. 1 1, pp , 26. [18) T. Gneiting and A. E. Raftery, "Strictly proper scoring rules, prediction, and estimation," Journal of the American Statistical Association, vol. 21, pp , 27.