Abstract
To mitigate the increase of anxiety resulting from the depletion of fossil fuels and destruction of the ecosystem, wind power, as the most common renewable energy, is a flourishing industry. Thus, accurate wind speed forecasting is critical for the efficient function of wind farms. However, affected by complicated influence factors in meteorology and volatile physical property, wind speed forecasting is difficult and challenging. Based on previous research efforts, an intelligent hybrid model was proposed in this paper in an attempt to tackle this difficult task. First, wavelet transform was utilized to extract the main components of the original wind speed data while eliminating noise. To make better use of the backpropagation artificial neural network, the initial parameters of the network are substituted with optimized ones, which are achieved by using the artificial fish swarm algorithm (AFSA), and the final combination model is employed to conduct wind speed forecasting. A series of data are collected from four different observation sites to test the validity of the proposed model. Through comprehensive comparison with the traditional models, the experiment results clearly indicate that the proposed hybrid model outperforms the traditional single models.
1. Introduction
With fossil fuels gradually drying up and the aftermath of the polluted environment becoming increasingly obvious, the world is sparing no effort to explore and exploit renewable energy for future power generation [1]. However, intermittence is the nature of the most renewable energy sources, which presents a great challenge to maintain and ensure the stability and reliability of power network [2]. Thus, effective energy storage technologies play a paramount role in mitigating the volatility [3]. As key ingredients, electric vehicles (EVs) would be widely regarded as an integral part of such a renewable energy system if they use electricity generated by a renewable energy source [4]. Besides, oil and coal supply uncertainty, growing mobility demand, and increasingly stringent regulations on pollutants and carbon footprints are expediting a paradigm shift towards sustainable power system and transportation [5, 6]. Until now, even though numerous clean energy technologies, such as plugin hybrid electric vehicles (PHEVs) [7], battery management system (BMS) [8], and battery electric vehicles (BEVs) [9], are being widely deployed and developed to significantly reduce the fuel consumption and the carbon emission throughout the world, they also increase the demand of power in the meantime. Thus, expanding investment on exploration and exploitation of the supplement of traditional energy is badly needed.
Wind power, the most recognized renewable energy source (RES), definitely possesses the prominent potential to become the resource of electricity, which plays a fundamental and indispensable role in social modernization [10]. Among different types of renewable energy, such as solar, wind, geothermal, and tidal energy, wind power is the fastestdeveloping and the most widely used one and also meets the demand for large scale commercial exploitation, owing to the maturity of the technology [11].
Currently, the technology of wind power has attracted worldwide attention and occupied a considerable share in the proportion of energy construction. Many countries take the development of wind power as an effective action to improve their domestic energy structures, preserve the ecological environment, alleviate environmental pollution, and even integrate it into the strategic development planning for the entire power system.
China, with such a great population, will soon replace America as the leading energy consuming country [12]. Depending on a series of positive policies, the wind power industry has stepped into rapid development and the present scale is exciting. The capacity of exploitable wind power is on rise year by year as technology progresses. According to a study started by the National Energy Administration (NEA), in the national wind power longterm plan, it is explicitly proposed that, to 2020, the national wind power installed capacity goal is 20,000,000 kilowatts [13]. In order to achieve such a capacity of power, nearly 708 million tons coal are needed according to the internationally accepted conversion standard, and the accompanying carbon emission even reaches to 480 million tons. Therefore, developing wind power efficiently does not only have economic significance, but also high ecological value.
With such a large capacity of wind power, the integration of wind power becomes an issue. Furthermore, high randomness of wind brings a big problem for a wind farm to provide electric grids to meet its necessary need of power, and this characteristic of wind power is a big challenge to traditional power system which has a relatively mature operation mode. Hence, accurate wind speed forecasting, as a necessary step in wind power integration, is particularly important, and these are the contributions of this paper.
There are three acknowledged methods for forecasting wind speed, including longterm prediction, mediumterm prediction, and shortterm prediction [14]. Briefly, longterm forecasting is used in the process of site selection for wind farms and in the maintenance and overhaul of the power network. When researching the dynamic economic dispatch scheme in a gridconnected wind power system, the mediumterm forecasting of wind that usually spans a few days or weeks comes in handy. As for the latter method, it mostly works to meet the control requirements of wind turbines. Overall, apparently, forecast error will become greater as the forecast horizon increases. Thus, the accuracy of shortterm forecasting is particularly important, for it lays the foundation of the renewable energy program, security evaluation, optimal dispatching model, and other critical decisions in power network dispatching.
Until now, many scientists have devoted research efforts to shortterm wind speed forecasting, and various achievements have been successfully applied. Numerical simulation, as early results, is usually used to construct mathematical models to simulate the dynamical physical behaviour of wind. These models include many linear, nonlinear, difference, and differential equations related with many conditions of wind, and those numerical solutions are usually achieved by calculation and iterative methods. However, because of the astatic and chaotic characteristics of the factors involved in the numerical simulation, as well as many uncertainties in the forecasting process, the results of these methods usually cannot meet the requirements of precision [15].
To solve these difficulties, many new models have been applied, which can roughly be categorized into physical, such as the Numerical Weather Prediction systems (NWPs), statistical, containing linear methods such as the Auto Regressive model (AR), Moving Average (MA), Auto Regressive Integrated Moving Average (ARIMA) model, and intelligent methods, which come from artificial intelligence and machine learning fields, such as the artificial neural networks (ANN) [16], and even by hybrid approach methods, which are combination of statistical and physical methods that use the analysis of a time series [17].
ARIMA, as a traditional statistical approach, has been proposed by some researchers to forecast wind speed, though it should be improved when dealing with the volatility and chaotic characteristics of wind speed. Ernst et al. supplied a method using past power measurements and meteorological predictions of the data of wind speed and direction interpolated at a wind farm [18]. The Grey prediction model GM, as a novel method for wind speed forecasting, as well as a wind power predictor, that reflects a Grey model that has a differential equation with a single variable, was provided by ElFouly. The structure of the artificial neural network, which was applied for forecasting the mean monthly wind speed in regions of Cyprus, is discussed by Baran, though the data were not sufficient [19]. In addition, a hybrid model that combined the Kalman filter (KF) with the artificial neural network was proposed based on an ARIMA model to further enhance the forecasting accuracy of wind speed. Sideratos and Hatziargyriou [20] employed a hybrid method to forecast wind speed in which three independent algorithms had been adopted. They were the Seasonal Adjustment Algorithm, the Exponential Smoothing Method (ESM), and the Radial Basis Function (RBF) neural network, and the final results indicated that the suggested hybrid method had satisfactory performance.
Backpropagation, as one of the most common learning algorithms among the artificial neural networks, has obvious advantages in nonlinear model fitting and predication because of its unique characteristic in terms of high tolerance with data errors. From the operational principle of the learning algorithm of the BPANN model, it can be observed that the initial weights and thresholds of the backpropagation neural network have a serious impact on the accuracy of the model [21].
In addition to the model, the characteristic of the initial data can also affect the final forecasting of the wind speed. In fact, the original data from the wind farm usually cannot be used in the model directly because there is much chaotic and scrambled information [22]; then, the WT is employed in this work to preprocess the data of tenminute wind speed to eliminate the high frequency disturbance. For the sake of improving the accuracy of wind speed forecasting further, a new optimization algorithm called the artificial fish swarm algorithm (AFSA) is proposed in this paper. The AFSA adapted here is used to find the proper parameters of the initial backpropagation artificial neural network, whose input is the processed tenminute wind speed data. The final results of the data showed that the hybrid optimal AFSA and BPANN method can perform better in shortterm wind speed forecasting compared with previous models.
The rest of this paper is organized as follows. The analysed data are described in Section 2. A general description of relative methods is provided in Section 3. In Section 4, the hybrid optimal model is discussed in detail; specifically, the process of selection of the optimal parameters of the propagation and the overall structure of the WAFSABPANN model will be presented and introduced explicitly. In Section 5, the simulation of the model and the evaluation of its final performance of wind speed forecasting are discussed in detail. The conclusions of this paper are presented in Section 6.
2. The Sources and Characteristics of the Data Set
To meet the demands of the model proposed in this paper, a variety of data were collected from turbines of wind farms from the Shandong Peninsula of China. Owing to its abundant sources of wind from the sea, which are richly endowed by nature, the Shandong Peninsula plays a more and more important role within the entire large sector of the Chinese wind power structure. However, it is difficult to forecast and evaluate wind speed because there are too many factors that can affect the changing of the wind speed. Our effective and reliable method is to find out the valuable and useful information, which the proposed model can then utilize to display the forecasting through the historical data.
The data used in this paper are collected from three different turbines of a wind farm. To ensure the fairness and consider the potential periodicity of seasons, the involved data series of wind speed cover four seasons in 2011. To be specific, data from four representative months in different seasons are chosen to implement the simulation and forecasting, namely, January, May, July, and September. For instance, tenminute wind speed data series from January 1 to January 31 in 2011 are collected. Each day, the data series span from 00:00 to 23:50, and the wind speed data are recorded every ten minutes. Figure 1 shows the whole profile of the wind speed data series spanning four seasons in one of the three chosen observation sites.
From Figure 1, the obvious feature of fluctuation and undulation of the original wind speed data series can be observed, and there are no rules that may be utilized to conduct forecasting through simple estimation. Besides, the results of parameter estimation and testing indicate that wind speed data series fit the Weibull distribution. Though many models can tolerate the instability of the original data, the accuracy of the forecasting models can still be impaired by the noise and disturbance suffused in them. Thus, there are active imperious demands to preprocess the initial data series to achieve a relatively smooth time series. Herein, wavelet transform was proposed to dispose the wind speed data series firstly, and preprocessed data were applied to execute the hybrid model, which will be discussed in the next section.
3. Description of Relative Models
Several methods are referred to in this paper, such as wavelet transform, backpropagation artificial neural network, and artificial fish swarm algorithm; each of them will be explained briefly below.
3.1. Wavelet Transform (WT)
Wavelet transform is a type of denoising method that can decompose signals with different frequencies. The essence of WT is to transform a onedimensional signal in the time domain into a twodimensional timefrequency signal. WT can carry out the multiresolution analysis, which indicates that WT has a high time resolution and low frequency resolution in the high frequency part; in the low frequency part, WT has a high frequency resolution and low time resolution [23]. The basic idea of WT is to approach the primitive function by using two of the above functions. Assume that the signal is the square integrable function; then, the WT of the signal is the inner product with :where and are the scaling and shift factor, respectively, and is the adjoin function of .
The basic steps of WT are as follows.
Step 1 (Decompose the original signal). Choose the wavelet basis and decomposition hierarchy and calculate coefficients of wavelet decomposition in each layer.
Step 2 (Address the threshold value of high frequency coefficients). Choose a threshold value for each decomposition hierarchy and process the high frequency coefficients to eliminate the noise centred on the high frequency part.
Step 3 (Carry out wavelet reconstruction of signal). Carry out the wavelet reconstruction for low frequency coefficients and high frequency coefficients such that the threshold value is processed quantitatively [24].
Through the wavelet transform, the high frequency wind speed data series with significant noise is eliminated, and the residual of the subseries is combined to substitute the original data. From Figure 2, it is obvious that the denoising data series is more smooth and stable than the original data. The simulation and forecasting results upon the preprocessed data employed in this paper also indicate the necessity and superiority of the wavelet transform.
It should be noted here that the entire preprocessed data are only used to train the hybrid model involved hereinafter. To be specific, wind speed data during the whole month are employed to use the wavelet transform, and the final forecasting of the next week is achieved from the original real wind speed data.
3.2. BackPropagation Artificial Neural Network (BPANN)
The schematic of the BP artificial neural network is plotted in Figure 3(a). Applying the classic backpropagation optimal algorithm, the BPANN is widely used in multidomains. This particular BP neural network has three layers: the input layer, the hidden layer, and the output layer [25]. The number of the nodes of the input layer is determined by our concrete data of wind speed. Herein, our data series is the tenminute wind speed time series, and three nodes were selected as the input layer [26, 27]. This choice is based on extensive experiments: the prediction results obtained with three nodes in the input layer are much better than other cases. According to the data series of wind speed, these three input nodes represent the history data at continuous time epochs. Similarly, only one node of the output layer is also used to adapt our data series. As to the hidden layer, it is not difficult to find that it has a great effect on the robustness and stability of the neural network. According to HechtNielsen [28], the proper number of the nodes of the hidden layer, that is, , was selected.
Before implementing the network, it is necessary to normalize the data sets, with both the input and the output data included. Thereafter, the outputs of all hidden layer nodes are calculated aswhere is the weight connected from input nodes to hidden nodes , while the bias of neural is , is the sigmoid function, and represents the value of each input node. Then, calculate the output of neural network:where is the connection weight from hidden layer to output layer, is the corresponding bias of neural cell, and the activation function is . At last, minimize the error via the training process. In this paper, the network was trained firstly; then, the network was validated, and finally we performed model testing. To see more details of the BPANN, [29] can be referred to.
3.3. Artificial Fish Swarm Algorithm (AFSA)
Artificial fish swarm algorithm is an effective optimal algorithm for increasing the exploitation capability during the process of searching. As we all know, fish can easily find the area with more nutrition through individual searching or following after the ambient fishes. Li et al. [30] introduced this efficient algorithm in 2002, which models the behaviours of fishes to get the global optimal solution.
The state of the artificial fish can be described as , and represents the initial variable. As the objective function, is the consistency of the current position where the fish is in. represents the distance of individual fish. Visual is the cognitive distance, while Step and represents the maximum step size and factor of congestion degree, respectively. There are four types of behaviours in the AFSA:(1)Foraging behaviour: the current status of an artificial fish is , and the fish selects a state in its sensing range randomly; if the problem is approximately maximum and (or if the problem is approximately minimum; because problems that are approximately maximum or minimum can be converted to each other; in the following, we only discuss the problem that is approximately maximum), then make a step forward in this direction; otherwise, randomly select state again and determine whether it meets the condition to move forward; after attempting for try_number times, if the condition still cannot be met, then choose a direction randomly to move a step:(2)Swarming behaviour: represents the current position of an artificial fish; the task is to explore the numbers of partners and the central location in the neighbourhood of the current position. If , which indicates that there is more nutrition at the centre of partners and less competitiveness, then make a step forward towards the direction of the central of partner; otherwise, execute foraging behaviour:(3)Following behaviour: the current status of the artificial fish is ; explore the number of partners and the biggest partner in the neighbourhood of the current position. If , which indicates that the surroundings of have more nutrition and are not so congested, then make a step forward towards the direction of ; otherwise, execute foraging behaviour:(4)Random behaviour: implementation of random behaviour is relatively simple; the action is to move a step in a direction that is randomly picked in the field. In fact, this is a default behaviour during foraging behaviour.
There is a bulletin board in the algorithm, which is also defined as an artificial fish, to record the information of the optimal individual [31]. Each of the artificial fish would compare its current state with the state recorded in the bulletin board after every movement, and then substitute it if the current state of the individual is superior to the one recorded in bulletin board before.
According to the nature of the problem to be solved, evaluate the current environment of the fish and then choose a corresponding behaviour. A common method is to choose a behaviour from the alternative behaviours that can make the most progress toward the optimal direction; if there is no choice that can make the next state better than the current state, then the fish has to take a random action.
4. The Structure of Hybrid WAFSABPANN Model
Just as in the aforementioned references, the initial data of wind speed usually contain much interference, which can affect the accuracy of forecasting; it is therefore a critical procedure for optimal models to select valuable information while filtering the disturbance. Therefore, a relatively mature WT model was introduced to preprocess the original data series, which is nonstationary wind speed data. The sublayers wind speed is a much more stable and smooth time series. Meanwhile, because the initial weights and thresholds of backpropagation artificial neural networks are initialized randomly, it is difficult to reach the best state in most cases. Therefore, finding the best initial parameters of the initial network is greatly needed [32]. Herein, the artificial fish swarm algorithm was applied to finish this work. Based on the WAFSABPANN model, the steps of training are depicted as follows.
Step 1. Given the structure of the backpropagation artificial neural network, there are three layers in the network, and the number of nodes is in the input layer, the hidden layer, and the output layer, respectively [33].
Step 2. Set the dimensions of the artificial fish, including the weights and thresholds of the artificial neural network:
Step 3. Initialize the basic parameters of the artificial fish swarm algorithm, such as the size of the population pop_size, the cognitive distance Visual, the maximum step size Step, and the factor of congestion degree , while Max_Gen represents the maximum number of iterations and the objective value.
Step 4. Initialize the iteration step . Produce artificial fish individuals with a population of pop_size randomly in the feasible region as the initial fish; meanwhile, the component of each individual fish should strictly be random, with digits between .
Step 5. Calculate the food concentration of the current position of each individual fish, choose the best one, and mark it on the bulletin board.
Step 6. Simulate each individual fish to execute the following behaviour and swarming behaviour, and then calculate ; after comparing and determining the bigger , execute the corresponding behaviour (the default action is to execute the foraging behaviour).
Step 7. Update the bulletin board after each iteration, and always substitute the exiting with more optimal in fish individuals.
Step 8. Judge if the iteration times Gen have already reached Max_Gen or if the error has met the requirement, and output the digit on the bulletin board once the termination condition is met; otherwise, make and go back to Step 6.
The pseudocode of the artificial fish swarm algorithm is listed in Algorithm 1.

Specifics parameters of the AFSA algorithm are presented in Table 1.
The hybrid method has better performance than the traditional backpropagation artificial neural network and other single models, and the structure of the WAFSABPANN is briefly shown in Figure 3, from which a visualized understanding of the hybrid optimum model can be got.
5. Simulation and Analysis
By means of the aforementioned hybrid model and the time series wind speed data, a series of simulations and forecasting have been operated using Matlab R2012b, and the results are discussed as follows.
5.1. Related Index
To ensure fairness and real forecasting results, the data volume, training data set, and forecasting data set are kept unchanged. Herein, three common indexes of error, such as the mean square error (MSE), mean absolute percentage error (MAPE), and the mean absolute error (MAE), are adopted to evaluate the performance of related models, and their specific formulations are shown below:where represents the concrete value of forecasting series and is the length of the series, while the corresponding value of original data series is .
5.2. Discussion of the Results from Different Seasons
Affected by many meteorological factors which are usually erratic, wind speed fluctuates all of the time, making it difficult to forecast. Based on a previous study, a hybrid new model is proposed in this paper. To achieve intuitive results of the hybrid model, the specific prediction data with a certain time interval on the first weeks of February, June, August, and October are collected, while the training data sets span the whole month of January, May, July, and September, respectively. Thus, the length of the prediction series is 1008 because of the time intervals, which makes it improper to list all these data under the limited space. In this case, forecasting data from several comparative models are recorded every 6 hours. The results from our hybrid model are compared with those from other traditional hybrid models. In this experiment, preprocessed data from four months in different seasons are utilized to train the neural network, and upon training the neural network, the final forecasting is made. Herein, data collected from January 1 to January 31 are utilized as the training set to train the network, upon which the prediction of wind speed from February 1 to February 7 is achieved. Next, preprocessed data acquired from May 1 to May 31 are used to train the artificial neural network, and the data from June 1 to June 7 are the relevant testing set. According to the same principle, the remaining training data series are formed. For example, the testing data from August 1 to August 7 are obtained from the neural network, which are trained by wind speed data spanning July 1 to July 31.
Tables 2, 3, 4, and 5 display the prediction results of four seasons, respectively. To be specific, forecasting results in the first week of February obtained from different models are presented in Table 2. Similarly, Tables 3, 4, and 5 represent the forecasting results in other three seasons.
As is depicted in Tables 2, 3, 4, and 5 and Figure 4, the specific forecasting results are enumerated vividly based on four different models: the initial backpropagation artificial neural network, the ARIMA model, the ARIMA with data from wavelet transforming (WTARIMA), and the hybrid BP network with optimal parameters as well as the preprocessed data. Through the comparison of data, it can be found that the MAPE of WAFSABPANN is relatively more smooth and smaller than the other three models. Here, we can know that the wind speed in February has a remarkable fluctuation, which results in a relatively big error in all four models. Through further analysis, a rough speculation can be formed: compared with other times, the falling of temperature in the winter may lead to the fluctuation of the wind speed.
5.3. Discussion of the Results from Different Observation Sites
To test the stability of the experiments as well as the effectiveness of the proposed hybrid model, the scope of the training data and testing data was enlarged, including four whole months in different seasons. At the same time, horizontal comparison is also employed to strengthen the validity of the proposed hybrid models. Furthermore, to eliminate random error in wind speed data, data fields from three different observation sites are applied in the process of simulation and forecasting. All these efforts aim to reduce the randomness of the experiments as well as testify the effectiveness of the proposed models, and the analysis results of data series gathered from three different observation sites are presented in Tables 6–9.
Obvious improvements in the performance of the hybrid models can be noticed, which can be read from all three error criteria. And the superiority of our WAFSABPANN model exists in observation site 1, site 2, and site 3. As is depicted in Table 6, it is not hard to find that there are 1.28%, 1.6%, 1.31%, 1.18%, and 1.13% reductions in MAPE, 0.1419, 0.1586, 0.1447, 0.1113, and 0.1199 reductions in MSE, and 0.0871, 0.1048, 0.0888, 0.073, and 0.0779 reductions in MAE, respectively, of the hybrid model proposed in this paper when compared with BPANN, ARIMA, AFSABPANN, WTBPANN, and WTARIMA model.
On other aspects, through the horizontal comparison of the data shown in Tables 6, 7, 8, and 9, WT’s superior performance on the three error indexes, compared with the performances of the other traditional models, can be seen. The average index of errors in all tables covering four seasons also indicates this rule: the hybrid model indeed has much more superiority when dealing with the original data set, which contains much noise and fluctuation. In other words, the process of denoising of WT has a critical effect on the final forecasting when analysing the comprehensive comparison, from which we can find an obvious improvement for the model when the wavelet transform is added. Synthesizing these diagrams, a rough rule can be observed, showing that the wavelet transform has an impact on simulation and forecasting: the tendency of the error index is more flat compared with the original error. The hybrid model (WAFSABPANN) can at least be confirmed to finish forecasting the wind speed with a tolerated error.
Moreover, it is necessary to consider the running time when implementing these models. On average, the BPANN model is the fastest which usually costs less than two seconds, and this can be seen in Tables 6, 7, 8, and 9. It is still very quick even the process of wavelet transform is added. The ARIMA and WTARIMA model are relatively fast with the running time about forty seconds. However, the speed of model always slows down when it involves the AFSA algorithm, which usually lasts 150 seconds. Through further analysis, it is not hard to explain the longer period of hybrid model proposed in this paper since there are so many loops and iterations in the AFSA algorithm.
Though much more time is needed when applying the hybrid model, it is still a reasonable choice to make wind speed forecasting for wind farms, for the improvement of prediction accuracy can compensate the costs in time and bring more profits both in technology and economy.
5.4. Comprehensive Comparison between Different Models
Table 10 shows the final comprehensive comparison between the traditional and hybrid models. Based on the metrics tabularized in Table 10, the overall perspective of each model can be seen distinctly. It should be noted here that all of those digits are the means of the whole experiment. For instance, the MAE of the BPANN model is 0.4284, which is the mean of the aforementioned forecasting results from these three observations sites covering four seasons. The MSE and MAPE are formed in the similar way.
It can be found that, when compared with BPANN, ARIMA, AFSABPANN, WTBPANN, and WTARIMA, there are 24.7%, 25.5%, 24.78%, 24.64%, and 25.53% improvements in MAE, 45.25%, 46.2%, 45.01%, 45.58%, and 46.93% improvements in MSE, and 21.28%, 21.17%, 21.87%, 20.1%, and 20.52% improvements in MAPE, respectively, of the hybrid model proposed in this paper. Similarly, reductions in MAPE, MSE, and MAE occurred at all three observation sites in four seasons which are presented in aforementioned tables.
Through the aforementioned table and analysis, a general rule could be summarised now: the proposed WAFSABPANN can always perform better than other referred models, which is obviously embodied in MAPE, MAE, and MSE. Above all, a visualized conclusion can be reached: the proposed hybrid model performs better than the others. In summary, the proposed optimal hybrid model can, to some extent, make relatively better forecasts than the traditional prediction models. In other words, the optimal algorithm introduced in this paper is suitable to be applied to optimizing the traditional backpropagation artificial neural network based on the wind speed data.
6. Conclusions
With the ongoing need for wind power, accurate wind speed forecasting has a strong impact on wind farm management, even on the entire wind power dispatch hybrid system. Therefore, a new hybrid model is proposed in this paper to solve this tough task. Firstly, based on WT, extra noise in the original wind speed data series which are nonstationary is eliminated while the remaining data contain most of the effective information, which is critical for forecasting wind speed. Next, to achieve the most proper parameters of the backpropagation artificial neural network, AFSA, as an effective method, is applied in this paper, and with the optimal parameters, the initial network is formed. Then, the preprocessed wind speed data are utilized to train the BPANN artificial neural network; in the final process, the forecasting of wind speed is obtained from the trained network. A number of sets of wind speed data from four different observation sites are gathered to execute the simulation and forecasting. Meanwhile, three error metrics are built up to evaluate the performance of the models. Through comprehensive comparison, the conclusion elicited is that the proposed intelligent hybrid model (WAFSABPANN) outperforms the other traditional models, such as BPANN, ARIMA, AFSABPANN, WTBPANN, and WTARIMA. Hence, WAFSABPANN is an effective and reliable model for shortterm wind speed forecasting in wind farms of China.
Though relatively superior precision of wind speed prediction can be obtained applying the hybrid model proposed in this paper, there is still room for improvements. Investigating these forecasting approaches, such as the combination of numerical simulation and intelligence algorithms, could be a future research topic.
Abbreviations
AFSA:  Artificial fish swarm algorithm 
ANN:  Artificial neural network 
AR:  Auto Regressive Model 
ARIMA:  Auto Regressive Integrated Moving Average 
BEVs:  Battery electric vehicles 
BMS:  Battery management system 
BPANN:  Backpropagation artificial neural network 
ESM:  Exponential Smoothing Method 
EVs:  Electric vehicles 
GM:  Grey prediction model 
KF:  Kalman filter 
MA:  Moving Average Model 
MAE:  Mean absolute error 
MAPE:  Mean absolute percentage error 
MSE:  Mean square error 
NWPs:  Numerical Weather Prediction Systems 
NEA:  National energy administration 
PHEVs:  Plugin hybrid electric vehicles 
RBF:  Radial Basis Function 
RES:  Renewable energy source 
WT:  Wavelet transform. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was funded by the Support Plan for Leaders of the Firstclass Discipline with Characteristics in Colleges and Universities of Liaoning Province of China under Grant (XKRC201403) and the project of DUFE research base (DUFE2014J29), the modeling and analysis of the inflow and outflow of fullcaliber regional water sources of China.