Abstract

Wind speed high-accuracy forecasting, an important part of the electrical system monitoring and control, is of the essence to protect the safety of wind power utilization. However, the wind speed signals are always intermittent and intrinsic complexity; therefore, it is difficult to forecast them accurately. Many traditional wind speed forecasting studies have focused on single models, which leads to poor prediction accuracy. In this paper, a new hybrid model is proposed to overcome the shortcoming of single models by combining singular spectrum analysis, modified intelligent optimization, and the rolling Elman neural network. In this model, except for the multiple seasonal patterns used to reduce interferences from the original data, the rolling model is utilized to forecast the multistep wind speed. To verify the forecasting ability of the proposed hybrid model, 10 min and 60 min wind speed data from the province of Shandong, China, were proposed in this paper as the case study. Compared to the other models, the proposed hybrid model forecasts the wind speed with higher accuracy.

1. Introduction

In the past few decades, with environmental degradation and resource depletion, renewable energy [1] has received more attention. Wind energy, one of the cleanest forms of renewable energy, is developing rapidly throughout the world. With the rapid increase in the utilization of wind energy, the primary concern is the security and stability of feeding electricity into the grid [2]. Wind speed high-accuracy forecasting is an important part of the electrical system monitoring and control. However, due to the instability of wind energy and the inherent complexity, transferring electricity into the power grid is limited and costly [3, 4]. To improve the efficiency of wind power and reduce the comprehensive cost of wind energy, accurate prediction of wind speed is necessary.

Many methods have been proposed to improve the forecasting accuracy of wind speed in recent decades. Based on the computational mechanism, these forecasting models can be grouped into four main categories: (i) physical models, (ii) statistical models, (iii) intelligence models, and (iv) hybrid forecasting models [5].

Physical methods [6, 7], which are based on the lower atmosphere or numerical weather prediction (NWP), can accurately forecast the wind speed. However, the physical methods require long running times and are not applicable for short-term forecasting. Statistical models [811], which are known as time-series-based models, do not apply historical data. These models are trained with measurement data and the differences between forecasted and actual wind speed are used to adjust the model parameters. The ARMA and the ARIMA models are the most popular models used to forecast future wind speed. Many forecasting results based on statistical models show that these models are useful in the wind speed forecasting fields [1215]. These models have numerous advantages, and the approaches need only historical wind speed data and are easy to implement. However, if the nonlinear characteristics of the wind speed series are prominent, the prediction accuracy of these methods decreases rapidly. The intelligent methods adopt artificial intelligence (AI) theories or evolutionary algorithms to forecast wind speed. Many intelligent methods are used to forecast wind speed, such as ANN (Artificial Neural Network) [1618], FLM (Fuzzy Logic Method) [19, 20], and SVM (Support Vector Machine) [18, 2126]. Unlike the single methods, the hybrid methods proposed by experts and scholars always combine several models to improve the accuracy of wind speed forecasting. Most of the recently proposed forecasting methods are hybrids, and the decomposition algorithms are often used to enhance their precision [27, 28].

Chaotic theory has been used to handle time series in many fields [2931]. Considering chaotic characteristics of the wind speed series, a hybrid prediction model was introduced in [32] using largest Lyapunov exponent prediction method to predict. Due to the inherent complexity of wind speed, describing the moving trend of wind speed and accurate prediction is difficult. Therefore, many studies use other methods to enhance the forecasting capacity of the original series. The specific methods include hybrid models that employ different approaches or combine different forecasting models to extract the inner traits of the original series in different aspects to perform wind speed forecasting. For the former type of hybrid models, the most common methods, such as wavelet transform (WT), singular spectrum analysis (SSA), and empirical mode decomposition (EMD), are used to preprocess the original series and forecast the wind speeds [3235]. These data processing methods are used to eliminate the influence of outliers on the forecasting accuracy, thereby improving the forecasting accuracy.

In this paper, a novel algorithm is proposed that hybridizes SSA (singular spectrum analysis), FAPSO (Firefly Algorithm and Particle Swarm Optimization), and RENN (rolling Elman neural network), to forecast wind speed. To verify the performance of the model, several hybrid models and single models are also used to forecast wind speed. In this model, besides the multiple seasonal patterns used to reduce interferences from the original data, the rolling model is utilized to forecast the multistep wind speed. To verify the forecasting ability of the proposed hybrid model, 10 min and 60 min wind speed data from the province of Shandong, China, were used as the case study.

The details of the algorithm are described below, and the flow diagram is shown in Figure 1.

Step 1. The SSA is used to decompose the original wind speed datasets into several subseries. Then, the new series is reconstructed. The wind speed data used in this paper is typically a chaotic time series, and the use of SSA can eliminate the influence of outliers and improve the prediction accuracy of the wind speed forecast model.

Step 2. The hybrid optimization algorithm (FAPSO) that combines the FA with the PSO is utilized to optimize the weights and thresholds of the ENN model. The optimization algorithm can provide better initial weights and thresholds to the ENN and improve the search ability. Compared with the single optimization model, the hybrid optimization model has better optimization effects.

Step 3. Construct the ENN model for the reconstructed series. Then, use the established model to forecast the one-step wind speed. The optimized ENN model can avoid getting trapped into local optimum and the global searching ability of the algorithm is enhanced.

Step 4. The rolling ENN model is used to forecast the multistep results. Multistep wind speed forecasting with high-precision is helpful for electricity production to produce various benefits, such as avoiding a power-grid collapse, reducing production costs, and reducing the spinning reserve capacity of thermal power units.

Step 5. The Diebold-Mariano test is used to validate the accuracy and stability of the proposed model.

2. Methodology

In this paper, numerous methods are involved. In this section, the relative algorithms including singular spectrum analysis, the firefly algorithm, particle swarm optimization, and the hybrid model are described in detail.

2.1. Singular Spectrum Analysis

The singular spectrum analysis [3638] is a signal processing technique capable of capturing the intrinsic oscillation modes of a signal. The SSA has two main stages: decomposition and reconstruction.

2.1.1. Wind Speed Time Series Decomposition

The decomposition has two main steps: embedding and SVD (singular value decomposition). During the embedding stage, given a time series , N is the effective length. Define L   called “window length” and .

To perform the embedding, the original time series is mapped into a sequence of lagged vectors of size L by forming lagged vectors

Then, the trajectory matrix is derived:

From the trajectory matrix, both the rows and columns of are subseries of the original series, and is a Hankel matrix, which means it has equal elements on antidiagonals.

In the singular value decomposition step, the singular value decomposition of the matrix can be computed through eigenvalues and eigenvectors of the matrix . Suppose () as the eigenvalues of , and suppose as the corresponding eigenvector. The singular value decomposition of the trajectory matrix is shown below:

The matrices are elementary matrices. The collection is called the th eigentriple of the SVD. Each eigentriple consists of an eigenvector, a factor vector, and a singular value. and are th left singular vectors and right singular vectors of , respectively.

2.1.2. Wind Speed Series Reconstruction

The decomposition has two main steps: grouping and diagonal averaging. During the grouping stage, is split into disjoint groups . Let ; then, the resultant matrix corresponding to the group is defined as

These matrices are computed for and the expansion of (3) leads to the decomposition as follows:

In the diagonal averaging step, the grouped decomposition in (6) will be transformed into a new series of length . Let be a matrix with elements   (). Diagonal averaging transfers the matrix to the series by the following formula:where ; ; ; if , ; and if not .

The original time series can be decomposed into series by SSAwhere is the result of diagonal averaging of the matrix .

2.2. Elman Neural Network (ENN)

The Elman recurrent neural network, proposed by Elman, is a partial recurrent network model [39]. Compared with classic feed-forward perception and pure recurrent network, ENN has a context layer that feeds back the hidden layer outputs in the previous time-steps. The context layer can enhance the ability of processing dynamic information and improve the forecasting accuracy.

The neurons contained in each layer are used to disseminate information from one layer to another. The nonlinear state space expression of Elman networks is as follows:where is the -dimension output node vector; is -dimension hidden layer node vector; is -dimension input layer vector; is the feedback state vector; , , and denote their corresponding weights; is the transfer function of the output neuron; and is the transfer function of the hidden neuron.

Then, adjust the weights of the network to minimize the squared error between the actual values and forecasting results:where is a target output vector.

Although the ENN has strong predictive power, the limitations are obvious. The initial weights and threshold values of ENN are randomly generated, the training speed is slow, and ENN is susceptible to falling into the local optimal value. The intelligent optimization algorithm can effectively overcome these shortcomings.

2.3. Firefly Algorithm and Particle Swarm Optimization (FAPSO)

The optimization algorithm is composed of the firefly algorithm and particle swarm optimization. Compared with a single optimization algorithm, the proposed optimization algorithm avoids many shortcomings and determines a better solution.

2.3.1. Firefly Algorithm (FA)

The firefly algorithm, proposed by Yang, is a multimodal nature inspired metaheuristic algorithm based on the flashing behavior of fireflies [40, 41]. The algorithm has proved effective in solving linear design problems and multimodal optimization problems.

The firefly algorithm has two stages, which are described as follows.

Step 1. The brightness is dependent on the intensity of light emitted by the firefly. Suppose there are a group of fireflies and the position for an th firefly is , where indicates the fitness value of the firefly. The brightness of a firefly is chosen to reflect the fitness value of its current position

Step 2. All fireflies have a unique attractiveness , which indicates the ability to attract other fireflies. The attractiveness is related to the distance factor at locations and , and between the two corresponding fireflies, and are given by

The attractiveness function of the firefly is computed aswhere is the largest attraction and is the coefficient of light absorption.

The movement of the less bright firefly toward the most bright firefly is computed aswhere is the randomization parameter and rand is a randomly selected number between the interval []. A rudimentary FA algorithm is outlined in Algorithm 1.

Parameters:
—the fitness function of firefly  .
—the maximum number of iterations.
—the number of fireflies.
—light intensity of firefly  .
—coefficient of light absorption.
(1) Set the parameters of FA.
(2) Generate initial population of fireflies randomly.
(3) FOR EACH  :   DO
(4) Evaluate the corresponding fitness function  .
(5) END FOR
(6) Determine light intensity.
(7) FOR EACH  :    DO
(8) Determine light intensity    using  
(10) WHILE     DO
(11) FOR EACH    DO
(12) FOR EACH    DO
(13) Move firefly    towards    in all    dimensions.
(14) IF  THEN
(15)
(16)
(17) END IF
(18) Evaluate the new solution and update the new light intensity  .
(19) END FOR  
(20) END FOR  
(21) END WHILE
2.3.2. Particle Swarm Optimization (PSO)

Particle Swarm Optimization (developed by Kennedy and Eberhart) [42, 43] is based on the behavior of birds. The principle of PSO is assuming a location has no mass or volume, flying like a bird in multidimensional space, and not only adjusts its position but also exchanges information about its current position in search space according to its own earlier experience and that of its neighbors. In this mechanism, members of a swarm communicate their information and modify their positions and velocities using the group information according to the best position appearing in the current movement of the swarm. The particles of the swarm find the optimal point by cooperation. The particle velocity and position updating formulas are shown below:where is the inertia weight; is the velocity of the th particle; and are acceleration coefficients (nonnegative constant); rand() is a random value between 0 and 1; is the position of the th particle; and denote local best particle of the th particle and global best particle among local bests at time .

A rudimentary PSO algorithm is outlined in Algorithm 2.

Parameters:
—the number of particles.
—the maximum number of iterations.
—the parameters of PSO.
(1) Set the parameters of PSO.
(2) Generate initial population of particles randomly.
(3) FOR EACH  :   DO
(4) ;
(5) ;
(6) END FOR
(7) Find the best value of particles.
(8) WHILE    DO
(9) Find the best fitness value for each candidates.
(10) FOR EACH  :   DO
(11) Calculate each particle fitness function.
(12) IF    THEN
(13)
(14) END IF
(15) END FOR
(16) Choose the candidate with the best fitness value of all the candidates
(17) FOR EACH :   DO
(18) IF    THEN
(19)
(20)
(21) END IF
(22) END FOR
(23) Update the values of all the candidates by using PSO’s evolution equations.
(24) FOR EACH :    DO
(25)
(26)
(27) END FOR
(28) END WHILE
2.3.3. Hybrid Optimization Algorithm

In this section, a modified optimization model hybridizing firefly algorithm and particle swarm optimization is proposed to improve the accuracy of wind speed forecasting. The specifics of the FAPSO are described below.

Step 1. The firefly algorithm was used to optimize fireflies.(1.1)Generate initial population of fireflies randomly.(1.2)Compute the brightness of each firefly by using objective function.(1.3)Move firefly and evaluate new fireflies.(1.4)Rank the fireflies and find the current best as the firefly researched.(1.5)Optimize fireflies in total.

Step 2. The particle swarm algorithm is used to search the best particle.(2.1)The fireflies, searched in the first step, are the initial population.(2.2)Calculate the fitness values of particles.(2.3)Find the global best position.(2.4)Update the positions and velocities of the agents.(2.5)Determine whether the termination conditions have been satisfied; if so, proceed to step (2.1) and continue to search; otherwise, continue to step (2.6).(2.6)Find the best global position.

2.3.4. Test of FAPSO

To verify the optimization performance and convergence speed of the modified algorithm, four benchmark functions are selected in this paper. These benchmark functions have different characteristics, which are used to fully investigate the optimization ability of the algorithm. The four common test functions are shown in Table 1. The experimental parameters of PSO and FAPSO are shown in Table 2.

Thirty experiments, searching for the minimum value point by 2000 iterations, were carried out independently. The results, including the maximum, minimum, average, and standard deviations, are displayed in Table 3. The results show that the optimization effect of FAPSO is better than PSO.

2.4. The Proposed Hybrid Model

In this paper, a novel algorithm, SSA-FAPSO-RENN, is proposed to forecast wind speed. SSA is used to acquire the moving tendency of wind speed and enhance the forecasting abilities. The hybrid optimization algorithm (FAPSO) that combines the FA and the PSO is utilized to optimize the parameters of the ENN model. To forecast the multistep wind speed, the rolling Elman neural network (RENN) model is used.

For the convenience of narrative, the proposed hybrid model is named the SSA-FAPSO-RENN model.

3. Experimental Simulation

In this section, the details of experimental simulation will be introduced. Wind speed series of 10 min and 60 min are used to verify the effect of the model.

3.1. Performance Metric

The primary concern is to determine whether the prediction model is superior to other models. The performance of the model is usually evaluated using statistical criterions.

To estimate the forecasting performance, the Diebold-Mariano (DM) test and three error criterions are adopted, including MAE, MAPE, and MSE. DM test [44] is a comparison test that focuses on the predictive accuracy and can be used to evaluate the forecasting performance of the proposed hybrid model and other comparing models. The details of DM test are given as follows: where is the loss function. and are the forecast errors from two models. is an estimator of the variance of . The hypothesis test is defined asThe null hypothesis is that the two forecasts have the same accuracy. Under the null hypothesis, the test statistics DM are asymptotically distributed. If , the null hypothesis will be rejected.

The detailed equations of these three error criterions are given as follows.

MSE (Mean Squared Error)

MAE (Mean Absolute Error)

MAPE (Mean Absolute Percentage Error)where and denote the real and predicted values at time , respectively.

3.2. Wind Speed Datasets

To verify the forecasting ability of the proposed hybrid model, 10 min and 60 min wind speed data (January 1, 2011, to November 9, 2011) from the province of Shandong, China, are proposed as the case study in this paper. In the two tests, multiple seasonal patterns are used to reduce interferences from the original data, March 1 to May 31 (spring), June 1 to August 31 (summer), September 1 to November 9 (fall), and January 1 to February 28 (winter), and the wind speed datasets are randomly selected.

To further assess the forecasting accuracy, every wind speed series is divided into a training set and a validation set. In addition, an entire day of data will be used as a test set to test the forecasting ability of the models.

The first case study is 10 min wind speed forecasting. The total number of available samples is 1152. The training set also includes 806 wind speed datasets and the validation set includes 140 wind speed datasets. The remaining data are used to calculate the predictive ability of these models. Figure 2 shows four wind speed datasets from three wind observation sites corresponding to the four seasons.

The second case study is 60 min wind speed forecasting. The total number of available samples is 1032. The training set includes 806 wind speed datasets, and the validation set includes 202 wind speed datasets, and the remaining data are the test set. Figure 3 shows four wind speed datasets (60 min wind speed) from three wind observation sites corresponding to the four seasons.

From Figures 2 and 3, several features can be summarized:(a)The data for four seasons are quite different.(b)There are three wind observation sites. The wind speed data from the same site is similar.(c)The intensity of the wind in winter is large but small for the wind in summer.(d)The experimental datasets reveal the chaotic nature and intrinsic complexity of wind speed.

3.3. Set Parameters of the Related Models

In this paper, a novel algorithm that hybridizes SSA, FAPSO, and RENN is proposed to forecast wind speed. Parameters have a large influence on prediction accuracy.

3.3.1. SSA Decomposition and Reconstruction

SSA performs well with complex irregular time series. The parameter setting of SSA is very important for the forecasting effect. The window length is the only parameter in the SSA decomposition process. The window length was chosen as an integer fulfilling the conditions: and , where is the data length. In this paper, the data length of 10 min and 60 min wind speed is 1152 and 1032, respectively; therefore, the window length chosen is 400.

Choosing allows the trend to be extracted simultaneously. Since the trend of wind speed series is complex, many eigentriples are required to reconstruct it. In this paper, the trend is reconstructed by eigentriples 1–150 in the 10 min wind speed forecasting experiments, and the 60 min wind speed trend is reconstructed by eigentriples 1–100. The graph in Figure 4 depicts the initial series and trend of the 10 min wind speed series from the wind observation site A.

3.3.2. Parameters of the Hybrid Model

Setting the parameters is very important for the prediction of wind speed. To compare the prediction effect of the model and attain a scientific conclusion, the initial parameters of these models need to be unified. The details are shown in Table 4.

3.4. Case Study One: 10 Min Wind Speed Forecasting

In this section, 10 min wind speed series, which are from four datasets of three wind observation sites, are used to test the forecasting capacity of the proposed hybrid model.

The forecasting results of the SSA-FAPSO-RENN model are compared with the forecasting results of the BPNN, ARIMA, RENN, SSA-RENN, and SSA-PSO-RENN. The BPNN, ARIMA, and RENN models are the single models, and the others are combination models. The parameters of BPNN are the same as those of ENN. The MAE, MAPE, and MSE values are the evaluation standard.

Figure 5 shows the multistep predicted results, the original wind speed data from spring datasets of wind observation site A, using the different involved models. The forecasting results are given in Tables 57. The values in bold indicate the smallest values of MSE, MAPE, and MAE. The smallest values of MSE, MAPE, and MAE are not all the forecasting results of the proposed model, but the minimum and the predictive values of the proposed model are very close. The results may differ with different error criterions. The proposed model is shown to have better prediction accuracy for most of the sample wind speed series.

To reflect the forecasting results more directly, the results of Tables 57 can be averaged. The average results were calculated by the results of three wind observation sites and four seasons. They are shown in Table 8.

Table 8 indicates the following:(1)The forecasting results of 1-step are better than 2-step, 3-step, and 5-step. For example, the MAPE values of the SSA-FAPSO-RENN model change from 5.06% to 6.43%, 7.69%, and 9.45% at 2-step, 3-step, and 5-step. This conclusion can be reached through other models.(2)Among all involved single models, the RENN model has the best performance except for the 1-step forecasting result, and the ARIMA model has the worst performance in every forecasting step.(3)Compared with combined models, the single model forecasting effect is relatively poor.(4)SSA-FAPSO-RENN offered the most accurate forecast value, with MAPE values of 5.06%, 6.43%, 7.69%, and 9.45% at 1-step, 2-step, 3-step, and 5-step, respectively.(5)The above conclusion can also be reached with MSE and MAE.(6)The SSA can improve the forecasting performance of the RENN model.

Wind speed high-accuracy forecasting, an important part of electrical system monitoring and control, is crucial to protect the safety of wind power utilization but is always a difficult and arduous task. Compared with the other forecasting models involved in this paper, the proposed hybrid model has better forecasting ability in the 10 min wind speed forecasting study.

3.5. Case Study Two: 60 Min Wind Speed Forecasting

In this case, one-hour wind speed series were used to test the forecasting capacity of the proposed hybrid model. Figure 6 depicts the initial series and trend of the 60 min wind speed series from the wind observation site A.

The forecasting results of proposed hybrid model, SSA-FAPSO-RENN, are compared with the forecasting results of BPNN, ARIMA, RENN, SSA-RENN, and SSA-PSO-RENN.

Figure 7 shows the multistep predicted results with the one-hour wind speed data from spring datasets of wind observation site A, using the different involved models. The estimated results of these predictions are given in Tables 911.

From Tables 911, the values in bold indicate the smallest values of MSE, MAPE, and MAE, and the minimum values and the predictive value of the proposed model are approximately the same; compared to 10 min forecasting results, one-hour wind speed forecasting has more error; with the increase of prediction steps, precision will rise; some wind speed series are more adaptable to other models, but overall, the proposed model has better prediction ability.

To reflect the forecasting results more directly, the results of Tables 911 can be averaged. The average results are shown in Table 12.

Table 12 indicates the following:(7)The forecasting results of 1-step are better than 2-step, 3-step, and 5-step. The MAPE values of the proposed model change from 10.23% to 11.78%, 13.39%, and 15.50% at 2-step, 3-step, and 5-step. Other models can also reach this conclusion.(8)Among all involved single models, the RENN model has the best performance in every forecasting step, and the ARIMA model has the worst performance. The MAPE values of the ARIMA model are 15.65%, 17.61%, 19.03%, and 21.02% at 1-step, 2-step, 3-step, and 5-step, respectively. However, the values of RENN reduce to 12.92%, 14.40%, 15.36%, and 16.67%.(9)SSA-FAPSO-RENN offered the most accurate forecast value. In 1-step forecasting experiment, the MAPE values of the combined model were reduced by 20.8%, 3.8%, and 2.1% compared with the Elman model, SSA Elman, and SSAFA Elman, respectively. Similar conclusions can be obtained from multistep experiments.(10)The above conclusion can also be achieved with MSE and MAE.

The forecasting results are generally good. The proposed hybrid model can be used to forecast 60 min wind speed. Compared with traditional single models and other models involved in this paper, the proposed model has the best forecasting ability. The forecasting results show that the model has a better performance in the 10 min wind speed forecasting study than the 60 min study.

The average results of study one and study two are shown in Table 13. The proposed model has the best performance as evaluated by MAPE, MAE, and MSE. The Diebold-Mariano values of the SSA-PSO-RENN and SSA-RENN models are larger than the upper limits at the 10% significance level, the DM values of the RENN and the BP models are larger than the upper limits at the 5% significance level, and the DM value of the ARIMA models is larger than the upper limits at the 1% significance level. The result indicates that the proposed hybrid model is significantly better than other models.

4. Conclusions

Wind power systems need to further develop accurate and reliable technology for short-term wind speed forecasting. Due to the influence of various meteorological factors, wind speed series are intermittent and randomly characterized, making it difficult to forecast wind speed using a single model. The focus of recent research has been the development of new methods and combinations of methods. However, individual models do not always achieve desirable performance. Hybrid models can decrease negative influences that are intrinsic in each of the individual models. These models can use the advantages of each individual model and are less sensitive, in certain cases, to the factors that cause the individual models to perform undesirably. Therefore, the hybrid model is more effective than individual models for wind speed forecasting.

To forecast the 10 min and one-hour wind speed more accurately, a new hybrid model, SSA-FAPSO-RENN, is proposed, which can overcome many limitations of single models, such as poor prediction accuracy and artificial parameters. The forecasting results show that the proposed model can improve the accuracy of 10 min and 60 min wind speed forecasting. Compared with other models involved in this paper, the prediction precision of the proposed model is the largest.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Social Science Foundation of China (Grant no. 13&ZD171).