Open access peer-reviewed chapter

Time Series and Renewable Energy Forecasting

By Mahmoud Ghofrani and Anthony Suherli

Submitted: March 29th 2017Reviewed: October 5th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71501

Downloaded: 295

Abstract

Reliability is a key important criterion in every single system in the world, and it is not different in engineering. Reliability in power systems or electric grids can be generally defined as the availability time (capable of fully supplying the demand) of the system compared to the amount of time it is unavailable (incapable of supplying the demand). For systems with high uncertainties, such as renewable energy based power systems, achieving a high level of reliability is a formidable challenge due to the increased penetrations of the intermittent renewable sources such as wind and solar. A careful and accurate planning is at the utmost importance to achieve high reliability in renewable energy based systems. This chapter will assess wind-based power system’s reliability issues, and provide a case study that proposes a solution to enhance the reliability of the system.

Keywords

  • availability
  • energy storage
  • renewable energy
  • reliability
  • wind

1. Introduction

The world is moving forward in technology as power systems lean toward renewable energy more and more each year. While the idea of using renewable energy has long been the focus of numerous researches from all over the world, the implementation itself is more complicated than said. Dealing with renewable energy proposes new challenges that must be carefully addressed and solved. The uncertainty of renewable energy sources, such as wind speed (for wind turbines) or solar radiation (for solar photovoltaic (PV) panels), and the fact that it is unreliable from time to time due to said uncertainty, are two of the major issues that rough up the transition from fossil based energy sources to renewable energy sources. The main objective of operational and planning strategies is to enable power systems to constantly and continuously meet the consumers’ demand or the system load. The volatility of renewable energy sources jeopardizes the power system’s ability to reliably meet this objective.

Climate change concerns, and new state initiatives are some of the factors that contributed in pushing and escalating the number of wind power based technology deployment during recent years. The stochastic nature of wind power resources makes it difficult to perform a reliable operation. While this issue has been frequently studied and numerous methods have been developed, a flawless solution for every case has not yet been designed. One of the proposed solutions is to use fast-responding units like gas generators as the operating reserves to keep up with the demand [1, 2], although doing that reduces the system’s efficiency and increases its operating costs [3]. Another possible solution is to install energy storage systems, which store wind power during low-demand periods and release power during periods when the system cannot provide sufficient power to meet the load [4]. This increases the flexibility of the power system as the energy storage system counterbalances the unexpected wind power fluctuations to more efficiently utilize the smoothened wind power for supplying the system demand.

In the upcoming future, energy storage systems are expected to be an essential part of electric grids. However, its deployment depends heavily on its economic advantages when compared to the more conventional operational practices. To come up with the most economically beneficial plan, a cost–benefit analysis must be done for each possible technology, especially in regulated utilities where the limited market opportunities diminish the potential economic benefits of storage technologies over gas-fired generators [3]. To assure the effectiveness of an energy storage system, we must approach the problem with an appropriate strategy [5, 6]. An optimal storage sizing strategy furnishes the system with the capability to stabilize against forecast uncertainty and integrate wind power more reliably [7, 8, 9], added with optimal scheduling, it also improves the system’s transmission capacity utilization.

For high wind penetrations, fast-response thermal units are used as reserve capacities to provide the fast ramping capability required to deal with wind power fluctuations. Recent developments in storage technologies have advanced its energy efficiency and enhanced its capability in dealing with fast ramping. Additionally, storage systems bring forth several benefits when compared to fast-response thermal units, such as efficiency enhancement of renewable integration, reduced emission, and improved utilization of grid assets.

Several applications have been proposed for energy storage systems, which include but are not limited to renewable capacity firming and reliability enhancement of renewable integration. Each application requires a case-by-case optimal allocation strategy, which might result in different solutions. Particularly, the matter of optimally sizing, siting, scheduling and operating storage systems to address the reliability issues of intermittent renewable integration is of great importance. A solid, probabilistic optimization framework is needed to supplement grid operability and reliability while at the same time reduces overall costs for systems with high wind penetrations. The framework developed by the author in Ref. [10] is adopted for this chapter and is provided in the next section. A case study is presented in this chapter to analyze the reliability of renewable energy based systems and compare storage technologies and conventional gas-fired alternatives for reliably integrating different wind penetrations. An economic analysis is also provided to calculate costs and benefits associated with each technology to determine the most economical solution.

2. Methodology

The following methodology is one possible solution example to model an intermittent renewable energy-based power system.

2.1. Wind, load, and equipment availability modeling

We use probability distribution functions (PDFs) to model the stochastic nature of load and wind generation, which parameters are calculated using 10 years of historical hourly data for load and wind speed [11]. The produced model will then be used to generate hourly samples for the planning period. We use Fuzzy C-Means (FCM) clustering to capture a statistical model that takes into account seasonal variations [12]. We grade each of the sample points with a value within the range of [0, 1], then we minimize the weighted distance between any sample point and a cluster center by using an iterative algorithm. The elbow method determines the total number of clusters [13]. By combining the FCM clustering and the elbow method, we categorize our planning days into 40 clusters of 24-hour wind speed and load samples. We utilize the maximum likelihood method to find the parameters of the PDFs for the samples. Two sets of 24 individual PDFs will represent each of the clusters for a 24-hour period.

2.1.1. Wind power modeling

The total power generated by a wind turbine can be calculated by the product of a simple kinetic energy equation through a cross sectional area A as follows [14].

P=12ρv3A=12ρv3πd24E1

where vis the wind speed in meters per second (m/s), ρrepresents air density in kg/m3 and dis the rotor diameter in meters.

The power output of a wind turbine depends heavily on the speed of the wind, and the wind speed itself can be best characterized by using the Weibull PDF, which formula is the following [11, 15]:

fvck=kcvck1evckE2

where kis a shape vector, cis the scale vector, and vis a vector of the measured wind speed. The average width of the wind speed distribution is determined by the shape vector k, while the scale vector indicates where the majority of the distribution lies and how wide the distribution-stretch is.

The wind power output can be calculated by using the power-speed curve [16]:

GW=0vvi,vvovvivrviGWrvivvrGWrvrvvoE3

where GWis the output wind power, and vi,v0,vr,vrepresents cut-in speed, cut-out speed, rated speed, and wind speed respectively.

2.1.2. Load modeling

The variation of the load is described by the Gaussian distribution [11]:

flμσ=12πσ2expLμ22σ2E4

where σ and μ represent the standard distribution, and mean of the Gaussian distribution respectively, and L represents the load demand.

2.1.3. Equipment availability modeling

Forced outage rate (FOR) of an equipment is the unavailability of the equipment estimated for a long-time period [17]. FOR models the availability of the equipment stochastically by the binomial PDF as follows:

fqnp=nqpq1pnqE5

where nis the number of units for each power plant and q = 0, 1, 2, …., n. The availability of each unit pis:

p=1FORE6

where FOR is basically the probability of the system’s unavailability. For systems with long operating cycles, FOR can adequately estimate the unavailability probability of units that operates under similar conditions. On the other hand, it is not an adequate estimator for systems with short demand cycles. The most important period in the operation of a unit is the start-up period, and a peaking unit (example of a system with short demand cycles) will have less operating hours with more start-up and shut-down periods [17].

2.2. Energy storage modeling

The model of a storage system must be able to handle the energy balance between the sum of the stored and generated energy and the load, where it stores excess energy gained from wind generation and releases the energy to supply the peak demand. We can use compressed air energy storage (CAES) to enhance the wind integration performance in a transmission network due to its beneficial features such as large power capacity, long lifetime, and low operation costs [18]. The charging and discharging equations of the storage system are as follows:

St=1dsSt1+ηcsLsttTE7
St=1dsSt1ηdsGsttTE8

where Strepresents the energy stored in the storage system at hour t, ηcsand ηdsrepresents the charging and discharging efficiencies for the CAES, Lstrepresents the storage loading capacity at hour t, and Gstrepresents the storage generating capacity at hour t, and dsrepresents the self-discharge rate for CAES.

The state of charge of the storage system at any time t is within the minimum and maximum storage capacity requirements:

SminStSmaxtTE9

where Sminand Smaxare the minimum and maximum storage capacities.

The stored power must not exceed the maximum power rating at any given time as follows:

PtPmaxtTE10

where Ptand Pmaxare the storage power at time t and the maximum storage power respectively.

The following is the ramping constraints for the storage:

GStGSt1RUstTE11
GSt1GStRDstTE12

where RUsand RDsare the ramp up and ramp down of the turbine for the storage system respectively.

CAES has an expected lifetime of 30 years [19].

2.3. Economic modeling

The storage cost is the sum of the energy and power costs associated with each energy storage technology. The storage cost is described by the following equation [20, 21]:

ICS=CS.Smax+CP.PmaxE13

where ICSis the cost of investment for the storage system, CSis the energy cost for the storage system, and CPis the power cost for the storage system.

The energy cost for CAES is 53 $/kWh, which includes the combined reservoir and the balance of plant costs. The power cost of CAES is around 425 $/kW [20], which includes turbine, compressor, and other power related costs.

The operation expenses are the sum of operation and maintenance (O&M) and fuel costs, which can be described by the following equation [22]:

OCSt=HR.GSt.CNGt+COM.PmaxtTE14

where OCStis the operation cost of the storage system, HRis the turbine heat rate for the storage system, CNGtis the cost of natural gas of the storage system, and COMrepresents the cost of operation and maintenance for the storage.

CNGt, COM, and HRare 4300 Btu/kWh, 5 $/MBtu and 2.5 $/kW-year [20, 22].

For a gas-fired conventional generator, the investment cost and heat rate are 695 $/kW and 8000 Btu/kWh respectively.

The total annual cost can be calculated by uniformly distributing the investment costs over the lifetime as follows:

A=d1+dN1+dN+11ICE15

where Ais the annual equivalent cost for the investment, dis the discount rate, Nis the life cycle of the investment, and ICis the investment cost.

We assume a discount rate of 10% and a lifetime of 30 years for the investment.

2.4. DC optimal power flow

We use optimal power flow (OPF) to find the steady state condition that at the same time minimizes the total operation and reliability costs. The objective function of the deterministic OPF is as follows:

Obj.Function=Mini=1ngaiPgi,t2+biPgi,t+ci+i=1nbIEARi×ILi,t=MinOCt+ILCttTE16

where IEARiis the interrupted energy assessment rate at each bus. ai, biand ciare the coefficients of the cost function for the ith generator, Pgi,trepresents the power output of the i-th generator at hour t, and ILC is the interrupted load cost.

The objective function above is subject to each of the following constraints:

Power balance equation:

i=1nbPgi,t=i=1nbPdi,ttTE17

where Pdi,tis the supplied load at bus i at hour t, and nbis the bus number.

Power generation and load limitations:

Pgi,tminPgi,tPgi,tmaxtTE18
Pgi,tminPgi,tPgi,tmaxtTE19

where Pgi,tminand Pgi,tmaxare the lower and upper generation limits for the i-th generator at hour t respectively.

Interrupted load:

ILi,t=PDi,tPdi,ttTE20

where PDi,tis the load demand at bus i at hour t.

Generation ramp up and ramp down:

Pgi,tPgi,t1RUitTE21
Pgi,t1Pgi,tRDitTE22

And the transmission line limitation:

i=1nbHri×Pgi,tPdi,tfr¯rΩ&tTE23

where Hriis the generalized distribution factor of line rwith respect to bus i, fr¯is the maximum transmission capacity for line r, and Ωis the set of transmission lines.

2.5. Probabilistic optimal power flow

A probabilistic OPF is a more appropriate approach when dealing with uncertainties of loads and wind power fluctuations, which process includes running the deterministic power flow continuously to account for the majority of possible system states. This chapter utilizes an approximate method called Hong’s point estimate method (2 m + 1 scheme) to characterize uncertainties, which uses the first few front most statistical moments of stochastic variables to approximate the probability functions [23].

K is the number of concentration points that we use to represent the statistical information of the random input variable in our K×m scheme. A location xi,kand a weight wi,k(xi,k, wi,k) represent the kth concentration of the random variable xi. In order to relate the input and output variables to each other, we apply the non-linear function Fx1x2xixm. The location of the kth value of variable xiis determined by the following equation:

xi,k=μxi+ξi,kσxiE24

where μxiis the mean for the input variable xi, ξi,krepresents the standard location for the input variable xi, and σxiis the standard deviation for the input xi. We assign a weighing factor wi,kto the current random output variable of the kth concentration. To determine ξi,kand wi,k, for the kth concentration of xi, we use the following equations [23]:

k=1Kwi,k=1mE25
k=1Kwi,kξi,kj=λi,jj=1,,2K1E26

λi,jin the equation above represents the jth standard central moment for the random variable xi, and its probability density function fxican be described as:

λi,j=MjxiσxijE27

The jth central moment of the random variable xiis given by:

Mjxi=xiμxijfxidxiE28

Once we obtain every concentration (xi,k, wi,k), we use the nonlinear function F to calculate the vector of random output variables Zikfor each point μx1μx2xi,kμxmas follows:

Zik=Fμx1μx2xi,kμxmE29

By using the values from Zik, and the weighing factors, the jth moments of the random output variables can be approximated by:

EZji=1mk=1Kwi,kZikjE30

We can extract the desired statistical information of our random output variable using a 2 m + 1 scheme by solving (25) for K = 3 and ξi,3=0. The standard locations and weight produced by the equation are:

ξi,k=λi,32+13kλi,434λi,32k=1,2ξi,3=0E31
wi,k=13kξi,kξi,1ξi,2k=1,2E32
wi,3=1m1λi,4λi,32E33

λi,3and λi,4are the skewness and kurtosis of xi.

The scheme above sets up ξi,3=0, which results in xi,k=μxiin (25), and yields m of the 3 m locations at the same point. By that done, it only requires one additional function evaluation for this particular location to complete 1 iteration of our probabilistic OPF. We update the corresponding weight to w0as follows:

w0=i=1mwi,3=1i=1m1λi,4λi,32E34

The deterministic DC-OPF is executed 2 m + 1 times in order to take all the random variables into account.

2.6. Reliability analysis

Reliability analysis provides an index to measure the degree of supply availability to meet the system demand. In the times when generated and stored energy is insufficient to supply the load, load is interrupted to maintain the power balance in the system. Load and generation variations as well as equipment failures are among the system uncertainties that could contribute to the load interruption in a power system. For wind turbines, the reliability model is a combination of a two-state model and power output model defined by (3). This combination is illustrated in Figure 1 to provide the reliability model for wind generators.

Figure 1.

Reliability model for wind generator.

The interrupted load in the system is equivalent to the amount of energy that is not supplied for each hour of the scheduling period, which can be described as:

ENSt=i=1nbILi,ttTE35

where ENStis the energy not supplied at hour t, and ILi,tis the interrupted load at bus iat hour t.

The interrupted load is defined as a random output variable whose first moment is calculated by (30), with j = 1. Expected energy not supplied (EENS) is then calculated for a one-year planning duration to provide a probabilistic index for our reliability analysis.

EENS=c=1Ct=1Ti=1nbEILi,tc.ncE36

where ncis the number of days within cluster c, C is the total number of clusters, c is the cluster number, and E is the average function.

We use the energy index of reliability (EIR) to estimate the reliability of the system, which can be calculated as follows:

EIR=1EENSEEE37

EE represents the expected energy demand of the system during the planning interval and is defined as:

EE=c=1Ct=1Ti=1nbPDi,tc.ncE38

2.7. Genetic algorithm optimization

We use a Genetic Algorithm (GA)-based optimization to install the energy storage with its optimal location and size. The GA begins by initially taking a set of randomly selected solutions, and then ranking the solutions based on their fitness values. We then perform recombination, crossover, selection, and mutation, to evolve the solution population. Once the satisfaction criterion is satisfied, we put the process into a halt. We assign a large penalty factor to the violated constraint to ensure satisfying constraints.

2.8. Proposed method

We model the storage system into our POPF as a load that stores excess, unconsumed energy generated by the system during off-peak periods. The storage system is modeled as a generator to release the stored energy to meet the peak load when sufficient transmission capacity is available. The location and scheduling of the storage systems are then optimized using GA. The optimized solution is the most cost efficient as it minimizes the total operation and interrupted-load costs for the span of the planning period. In order to optimally enhance the grid operability for wind integration, the storage technologies must possess an adequate capacity per the system’s need. The fitness function that we use for the proposed method is the total weighted sum of the system’s cost for each cluster over the planning period, as follows:

Fit.Function=Minc=1Ct=1TOCt+ILCtcncE39

where OCis the operation cost of the system, ILCis the interrupted-load cost of the system, and nc represents the number of days within the cluster.

The proposed GA-based POPF can be described in the following steps:

1. Input wind speed, loads, and FOR data

2. Initialize the first population.

A. For t = 1, until t = T:

3. Initialize the first input variable by setting i = 1 and EZ=0& EZ2=0

B. For i = 1, until i = m:

4. Select input random variable xi

5. Calculate ξi,k, wi,k, λi,j

6. Initialize k=1.

C. For k= 1, until k= 3:

7. Calculate xi,k

8a. If GWt>Lt, model the storage as a variable load with the following constraints:

0LStminSmax1dsSt1Pmax

8b. If GWt<Lt, model the storage as a generator with the following constraints:

0GStmin1dsSt1SminPmax

9. Run Deterministic OPF using Zik=Fμx1μx2xi,kμxm

10. Calculate S for charging-discharging by using Eqs. (7 and 8)

11. Calculate OCt, ILCt, and EENSt

12. Update raw moments using the following equations:

EZ=EZ+wi,kZik;EZ2=EZ2+wi,kZik2

13. If k= 3, go to step 14, if not, go to step C with k=k+1.

14. If i= m, go to step 15, if not, go to step B withi=i+1.

15. If t = T, go to step 16, if not, go to step A with t=t+1.

16. Evaluate the fitness function and constraints

17. Generate children by using crossover and mutation

18a. If termination criteria are not met, produce next generation by selection and combination and go to step 2.

18b. If termination criteria are met, calculate the statistical output information.

3. Case study

We evaluate our proposed method by applying it on the IEEE 24-bus system with the goal of solving for the optimal size and location for the storage units [10, 24]. In our simulations, we test our method for different situations by inputting different wind penetrations. To take into account possible geological restrictions for CAES deployment in a real-world situation, we excluded busses 2, 7, 8, 11, and 17 from candidate locations on purpose. A wind farm is also pre-determinedly installed at bus 14 in each of the case studies. A diagram of the IEEE 24-bus system is shown in Figure 2. Tables 13 provide more information regarding the flow limitations, generators’ cost functions, and IEAR values of the IEEE 24-bus system in use for the case study.

Figure 2.

IEEE 24-bus system.

From busTo busFlow limit (MW)From busTo busFlow limit (MW)
121751113500
131751114500
151751213500
241751223500
261751323500
391751416500
3244001516500
491751521500
5101751524500
6101751617500
781751619500
891751718500
8101751722500
9114001821500
9124001920500
10114002023500
10124002122500

Table 1.

Transmission flow limitations.

GeneratorCost function coefficients
ai$/MW2hbi$/MWhci$
G10.10371.051313.6
G20.10871.041168.1
G30.09066.191078.8
G40.09167.26969.8
G50.07871.60958.2
G60.07871.60958.2
G70.10073.90471.6
G80.09073.90471.6
G90.09869.70445.4
G100.10166.51702.7

Table 2.

Coefficients of the cost function for the generators.

BusNo.12345678
IEAR6.204.895.305.626.115.505.415.40
BusNo.910111213141516
IEAR2.304.145.393.413.013.54
BusNo.1718192021222324
IEAR3.752.293.64

Table 3.

IEAR ($/kWh) values at each bus in the IEEE bus-system.

20% WP & 2500 MW PL*30% WP & 3000 MW PL40% WP & 3500 MW PL50% WP & 4000 MW PL
CS**CA***CSCACSCACSCA
Optimal Placement (Bus #)14142323
Smax(MWh)77.81778.715241858
Pmax(MW)59.6318.00258.6138.1652.3415.5947.5807.6
EIR%99.8099.8099.2899.2899.0399.0399.1999.19
Wind utilization (%)94.4986.4288.6978.5880.2670.9972.0662.22
OC(109$)0.5880.6030.6200.6610.6560.7110.6960.774
GC(109$)0.0020.0120.0370.077
ILC(109$)0.0410.0410.1940.1940.3240.3240.2950.295
SC(109$)0.0050.0240.0530.072

Table 4.

Simulation results for different wind penetrations (WP).

Peak load.


Centralized storage.


Conventional alternative.


We define wind penetration (WP) as the ratio between the wind capacity installation and the system maximum load. We use real-world historical data obtained from the BPA for the system load [25] and from Mesonet (Ames Station) for the wind speed [26], to create a more realistic simulation environment. To calculate the cost of electric service reliability in the IEEE 24-bus system that we are going to run our demo on, we use the values of IEAR for our load busses [27]. The cost of the storage system, storage cost (SC), is equal to the sum of the cost of investment (A), and its cost of operation (OCs) for the planning period. We can then calculate the cost of conventional generation (OC) by excluding the storage operation cost from the total operation cost.

The objective is to achieve the maximum possible reliability level. Our scenario’s goal is to solve for the optimal placement and sizing for the storage system to meet the reliability objective. Our control strategy is to use the available wind energy to supply the load first, followed by utilizing the existing conventional generation capacity, and last, if necessary, to discharge power from the storage system to satisfy the load. The result of the simulations including the comparison with other conventional alternatives is shown in Table 4. Same reliability level is considered for both technologies to make a fair economic comparison. Our storage system, which enhances the reliability of wind integration, can be economically evaluated by comparing the sum of its associated costs with the total cost for the conventional alternative. The cost–benefit analysis from our simulation results shows the economic merits of the CAES, which can be found in Table 4.

4. Conclusions

Our case study concludes that energy storage technologies are more economic and technically sound options than fossil-fuelled generators to reliably and efficiently integrate intermittent renewable energy such as wind. The merits of energy storage application for reliability enhancement of renewable integration become even more highlighted when the emission costs associated with fossil-fuelled generators are included in the evaluation. This provides the subject of future studies.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Mahmoud Ghofrani and Anthony Suherli (December 20th 2017). Time Series and Renewable Energy Forecasting, Dependability Engineering, Fausto Pedro García Márquez and Mayorkinos Papaelias, IntechOpen, DOI: 10.5772/intechopen.71501. Available from:

chapter statistics

295total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Introductory Chapter: Introduction to Dependability Engineering

By Fausto Pedro García Márquez and Mayorkinos Papaelias

Related Book

First chapter

Digital Filters for Maintenance Management

By Fausto Pedro García Márquez and Diego José Pedregal Tercero

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us