Next Article in Journal
Effects of Processing Conditions on the Simultaneous Extraction and Distribution of Oil and Protein from Almond Flour
Next Article in Special Issue
A DEA Approach for Assessing the Energy, Environmental and Economic Performance of Top 20 Industrial Countries
Previous Article in Journal
New Aspects on the Modeling of Dithiolactone-Mediated Radical Polymerization of Vinyl Monomers
Previous Article in Special Issue
Investment Decisions of Fired Power Plants on Carbon Utilization under the Imperfect Carbon Emission Trading Schemes in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wind Power Short-Term Forecasting Hybrid Model Based on CEEMD-SE Method

1
School of Economics and Management, North China Electric Power University, Beijing 102206, China
2
Beijing Key Laboratory of New Energy and Low-Carbon Development, Beijing 102206, China
3
North China Power Dispatching and Control Centre, Beijing 100053, China
*
Author to whom correspondence should be addressed.
Processes 2019, 7(11), 843; https://doi.org/10.3390/pr7110843
Submission received: 20 October 2019 / Revised: 4 November 2019 / Accepted: 6 November 2019 / Published: 11 November 2019
(This article belongs to the Special Issue Energy, Economic and Environment for Industrial Production Processes)

Abstract

:
Accurately predicting wind power is crucial for the large-scale grid-connected of wind power and the increase of wind power absorption proportion. To improve the forecasting accuracy of wind power, a hybrid forecasting model using data preprocessing strategy and improved extreme learning machine with kernel (KELM) is proposed, which mainly includes the following stages. Firstly, the Pearson correlation coefficient is calculated to determine the correlation degree between multiple factors of wind power to reduce data redundancy. Then, the complementary ensemble empirical mode decomposition (CEEMD) method is adopted to decompose the wind power time series to decrease the non-stationarity, the sample entropy (SE) theory is used to classify and reconstruct the subsequences to reduce the complexity of computation. Finally, the KELM optimized by harmony search (HS) algorithm is utilized to forecast each subsequence, and after integration processing, the forecasting results are obtained. The CEEMD-SE-HS-KELM forecasting model constructed in this paper is used in the short-term wind power forecasting of a Chinese wind farm, and the RMSE and MAE are as 2.16 and 0.39 respectively, which is better than EMD-SE-HS-KELM, HS-KELM, KELM and extreme learning machine (ELM) model. According to the experimental results, the hybrid method has higher forecasting accuracy for short-term wind power forecasting.

1. Introduction

With the increasingly serious problem of energy shortage and environmental pollution, energy saving and emission reduction strategy is the key measure to promote green development [1,2], and renewable energy is increasingly favored by people [3]. As one of the most promising renewable energy sources, wind power has been paid worldwide attention owing to its clean, low-carbon and economic characteristics [4,5]. In recent years, China has vigorously developed renewable energy, and both the installed capacity and the grid-connected capacity of wind power are becoming larger and larger. According to the relevant data released by the National Energy Administration, in 2018, China added 21.27 million kW of wind power grid-connected installed capacity [6], and the cumulative grid-connected installed capacity reached 184 million kW, accounting for 9.7% of the total wind power installed capacity. Although wind power generation technology is gradually matured, due to the intermittence, randomness, fluctuation, and uncertainty of wind energy, the large-scale access of wind farms brings tough challenges to the safe and stable operation of power systems [7,8]. Therefore, accurate and effective wind power forecasting technology is necessary for the security and stability of the power grid. Firstly, accurately forecasting wind power can provide decision support for grid dispatching plan and power market transaction [9], strengthen the stability and safety of the power system [10]; secondly, it can be helpful to the reduction of systems rotating reserve capacity and power generation cost, and the enhancement of the wind farms’ economic benefits and the utilization rate of wind power by users [11]; in addition, with the advancement of power market reform, accurate short-term wind power forecasting can also provide relevant basis for the grid-connected wind power sales under power market conditions [3], promote inter-provincial priority to absorb renewable energy, and reduce the risk brought by wind power uncertainty for power market participants [12]. Thus, it follows that wind power forecasting is of great significance to the safe operation of power system and the sound development of the power industry.
At present, scholars all over the world have conducted extensive research on the models and methods of wind power forecasting, mainly including three categories: physical models, statistical models, and intelligent forecasting models [13,14,15]. The physical model is more complex and relies on meteorological and wind speed data provided by numerical weather prediction (NWP) system, which requires high accuracy and completeness of data, and is more geared for long-term prediction of wind power [16,17]. Although the physical model does not require the support of historical data, NWP data needs a lot of calculations in order to forecast wind power, resulting in fewer engineering applications in China. The statistical methods mainly include regression analysis, exponential smoothing, time series, and Kalman filtering, which are superior to physical methods when forecasting short-term wind power [18]. Statistical models are more mature in the application of wind power forecasting than physical methods by studying the information related to time and space in the data. However, due to the fact that the linear models are commonly used in statistical methods, the randomness of meteorological data cannot be accurately expressed, and the forecasting results of wind power series are poor because of the non-linear and non-stationary characteristics of the series [19,20]. In recent years, artificial intelligence has developed rapidly, and a variety of intelligent forecasting models and machine learning algorithms have been successfully used to wind power forecasting, and good results have been obtained. Common intelligent forecasting models include artificial neural network (ANN) [21], support vector machine (SVM) [22,23], least squares support vector machine (LSSVM) [24], etc. The single intelligent algorithm itself has certain defects [25]. It is difficult to ensure high forecasting accuracy, for it cannot accurately grasp the law of wind power variation, resulting in large errors at the forecasting points where wind power fluctuates sharply. Therefore, some hybrid models of intelligent algorithms have been proposed. And it can be seen from the empirical results that the hybrid model has higher forecasting accuracy and better fitting effect. In addition, because the time series of wind power data are non-stationary sequences with strong randomness and volatility, data pre-processing method can effectively improve the wind power forecasting accuracy [26]. Recently, various data preprocessing methods have been used in the forecasting the short-term wind power to achieve the purpose of feature dimension reduction and data noise reduction. In general, the results are satisfactory. There are some common data preprocessing methods used by researchers, including principal component analysis (PCA) [27], correlation analysis [28], cluster analysis [29], wavelet transform (WT) [30], empirical mode decomposition (EMD) [31] and other corresponding improved methods. For example, the empirical wavelet transform (EWT) is an innovative adaptive wavelet decomposition method. The core idea of EWT is to extract modal signal components with Fourier compact support by constructing orthogonal wavelet filter banks adaptively. Compared with EMD method, it can avoid the phenomenon of mode aliasing and false modes, and get fewer modal signal components.
In wind power forecasting, the common hybrid models usually combine data preprocessing methods with intelligent forecasting models. Through the data preprocessing method, data screening and dimension reduction can be performed to reduce data redundancy. Then, through the signal decomposition algorithm, multiple stable components’ sets of wind speed or wind power can develop. After obtaining the forecasting results of each component with intelligent forecasting model, the final forecasting outcome of wind speed or wind power can be achieved by reconstructing the results. Du et al. [32] used complementary ensemble empirical mode decomposition (CEEMD) to decompose wind power series, by doing so, the noise data can be eliminated, and the main features of original data can be extracted. Then, the improved wavelet neural network (WAN) was used to forecast wind power. Jiang et al. [33] considered that the existing wind speed forecasting models neglect the fluctuation influence between adjacent wind turbines, resulting in poor forecasting accuracy. For this purpose, they constructed a combined wind speed forecasting model, which uses gray correlation analysis to screen the fluctuation information and inputs the obtained fluctuation information into the SVM for wind speed forecasting. Yan et al. [34] proposed a new hybrid wind speed prediction method, which uses correlation-assisted discrete wavelet transform (DWT) for decomposition the wind speed sequence and utilizes generalized autoregressive conditional heteroscedastic (GARCH) to reflect the fluctuation of the subsequence and adjust the parameters of the prediction model based on LSSVM in real time, so as to reflect the change of wind speed better. This method has achieved good results in accuracy and stability. Fei [35] proposed an integrated model for wind speed forecasting using EMD and multi-kernel relevance vector regression (MkRVR) algorithm. The wind speed was decomposed into intrinsic mode functions (IMFs) with different frequency ranges by EMD method, and the forecasting model was established by using MkRVR for the forecasting of each decomposed signal. The empirical results proved that the EMD-MkRVR model is accurate and efficient in predicting wind speed. Khosravi et al. [36] studied the prediction effect of three machine learning algorithm models for wind power output, including multilayer feed-forward neural network (MLFFNN), support vector regression with a radial basis function (SVR-RBF), and adaptive neuro-fuzzy inference system (ANFIS). The input variables of the three models were temperature, pressure, relative humidity and local time. The comparison between the empirical results showed that SVR-RBF was superior to MLFFNN and ANFIS-PSO. Wu and Lin [37] firstly used variational mode decomposition (VMD) to decompose the original wind speed sequence into subsequences of various frequencies. Then bat algorithm (BA) was used to the parameters optimization of least squares support vector machine (LSSVM), and each subsequence was forecast by the improved LSSVM model. Finally, the results of wind speed forecasting were achieved by superimposing the forecasting results of each subsequence. Jiang et al. [38] developed a hybrid wind speed prediction approach based on EMD, VMD and sample entropy (SE). The performance of the hybrid approach was verified by the real measured wind speed data of two cases, which proved that it had better prediction effect. Tian et al. [39] constructed a hybrid wind speed prediction model consisting of data preprocessing, optimization model and prediction module. This model balanced the contradiction between prediction accuracy and stability. In the empirical study, the mean absolute percentage error (MAPE) of the model was less than 4%, which confirmed that it has high prediction accuracy. Yin et al. [40] constructed a hybrid model for short-term wind power forecasting. The extreme learning machine (ELM) optimized by the crisscross optimization algorithm (CSO) was used to forecast, which solved the problem of premature convergence of ELM and had high forecasting accuracy. Through the analysis of the existing literature, it is found that the existing literature pays insufficient attention to data preprocessing. The effective data preprocessing method can improve the input data quality of the model and improve the prediction accuracy. However, VMD, EEMD, CEEMD and other methods will get more subsequences when decomposing the original power series, which increases the complexity of forecasting. As an improved neural network, KELM has the advantage of fast training speed and strong generalization ability. The KELM method optimized by HS has better prediction performance and stronger global search capability.
On account of the above analysis, for the purpose of obtaining more accurate wind power forecasting results, this paper proposes a multi-step hybrid wind power forecasting model combining complete ensemble empirical mode decomposition (CEEMD), sample entropy and improved extreme learning machine with kernel (KELM). The hybrid forecasting models is based on the combination of mixed signal decomposition algorithm and optimized machine learning algorithm, which includes data preprocessing stage, optimizing stage and forecasting stage. Firstly, Pearson correlation analysis is adopted to screen out the input data which have high correlation with wind power, and extract the main features of the original data, so as to reduce data redundancy. Then, CEEMD method is utilized to decompose the wind power time series into a series of modal components, so as to eliminate noise and improve the quality of input data of the forecasting model. The sample entropy is used to classify and reconstruct subsequences, which can effectively reduce the amount of calculation. Secondly, the harmony search (HS) algorithm is used to optimize the parameters and improve the search ability of KELM model [41]. Finally, the improved KELM is applied to forecast each component, and the forecasting results of respective component are reconstructed to obtain the wind power forecasting value. The KELM model is faster in training and overcomes the shortcomings of easily falling into local optimal solution. In this paper, different Chinese wind farms are used as practical cases, and the multi-step hybrid wind power forecasting model is used to forecast wind power in ultra-short term and short term. Under the same conditions, the EMD-SE-HS-KELM, HS-KELM, KELM and ELM methods are compared to validate the effectiveness of the wind power forecasting model constructed in this paper. The innovations of this paper are as follows:
(1)
Ensure the quality of input data of the forecasting model. Most wind power forecasting models only adopt a single data preprocessing method. Thus, the accuracy is limited due to the fact that the complexity of original data is not properly dealt with. This paper first analyses the correlation of indicators, eliminates the indicators with low correlation degree to reduce data redundancy, and then uses CEEMD to decompose wind power to improve the quality of input data of the forecasting model, and use sample entropy (SE), which was proposed by Richman [42] to classify and reconstruct subsequences to reduce the computational complexity.
(2)
Realize multi-step forecasting of wind power. The hybrid wind power forecasting model constructed in this paper comprises three parts: data preprocessing, optimization and forecasting. Reasonable multi-level forecasting model can be more supportive for decision-making.
(3)
Balance model forecasting accuracy and stationarity. This paper uses a variety of methods to enhance the forecasting accuracy of wind power and improve the stability of the model. The KELM model with faster training speed is used, and the parameters of the KELM model are optimized by the HS algorithm to improve the search performance.
(4)
Verify the performance of the forecasting model comprehensively. According to the real measured data of wind farms in China, the ultra-short-term and short-term forecasting of wind power are carried out by adopting the model respectively, and the comprehensive performance of the forecasting model is investigated by calculating four forecasting accuracy evaluation indicators, which confirms the feasibility of multi-scenario application of the model.
The rest of the paper is organized as follows: Section 2 introduces the method used in the wind power forecasting model constructed in this paper. Section 3 establishes a hybrid wind power forecasting model. Firstly, the Pearson correlation coefficient between the original data set and the wind power is calculated, and the data with low correlation degree is eliminated. CEEMD- SE is used to decompose and reconstruct the wind power time series, which can eliminate noise and reduce the amount of computation; and the KELM model optimized by HS wind power is utilized to forecast wind power. Section 4 validates the effectiveness of the model proposed in this paper by comparing it with EMD-SE-HS-KELM, HS-KELM, KELM, and ELM methods. The conclusion is shown in Section 5.

2. Materials and Methods

2.1. Pearson Correlation Coefficient

The Pearson correlation coefficient is calculated as follows:
P x , y = cov ( X , Y ) σ X σ Y = E ( ( X μ X ) ( Y μ Y ) ) σ X σ Y = E ( X Y ) E ( X ) E ( Y ) E ( X 2 ) E 2 ( X ) E ( Y 2 ) E 2 ( Y )
where ( X , Y ) are two different variables. The value of P x , y whose range is [ - 1 , 1 ] reflects the degree of linear correlation between Y and X . The following conclusions can be drawn:
If | P x , y | 1 , it means that the correlation degree between Y and X is stronger; In contrast, if | P x , y | 0 , the correlation degree is weaker, or they are non-linear correlation, or even irrelevant.

2.2. Complementary Ensemble Empirical Mode decomposition (EMD)

EMD is widely applied into signal analysis, which was initially proposed by Hibert-Huang as an adaptive data mining methodology [42]. The complex signal is disintegrated into a finite number of intrinsic mode functions (IMFs) components and a residual component through EMD method. The EMD method can be theoretically employed to the decomposition of signals of any time type, as well as signals according to the temporal feature scale of the data itself. It is superior to other methods in handling non-stationary and nonlinear data [42]. However, it also is liable to produce serious defects of mode mixing problem [43], which will affect the accuracy of intrinsic mode functions. CEEMD is an improved algorithm of EMD and it can effectively reduce the occurrence of this phenomenon by using noise characteristics. CEEMD decomposition is based on EMD, and its implementation process is as follows:
(1) Add N pairs of white noise to original sequence to obtain a set containing 2 N signals. The auxiliary noise is Gaussian white noise n i ( t ) ( i = 1 , 2 , , N ) with a mean of 0 and a constant amplitude coefficient k . When N ranges from 100 to 300, the value of k is 0.001~0.5 times of the signal standard deviation.
[ x i 1 ( t ) x i 2 ( t ) ] = [ 1 1 1 1 ] [ x ( t ) n i ( t ) ]
where x ( t ) is original signal; n i ( t ) is auxiliary white noise; x i 1 ( t ) and x i 2 ( t ) are the signal pairs after the noise is added.
(2) EMD decomposition is performed on the obtained 2 N signals, and a set of I M F components is obtained for each signal. The j -th I M F component of the i -th signal is I M F i j , and the residual component is recorded as the last I M F component.
(3) The corresponding I M F components are averaged to obtain the components of each phase after the original sequence x ( t ) is decomposed by CEEMD:
I M F j = 1 2 N i = 1 2 N I M F i j
where I M F j denotes the j -th I M F obtained through the decomposing process of the original signal by CEEMD.

2.3. Sample Entropy Theory

After decomposing the data sequence by CEEMD, n subsequences can be obtained. If the forecasting model is used directly to model and predict each subsequence, the calculation scale will be large. Therefore, the sample entropy theory can be used to analyze the complexity of each subsequence, and the subsequences can be classified and reconstructed, which can lead to the significant reduction of the computational complexity. The SE is initially proposed by Richman [44]. SE is used to measure the possibility of generating new patterns in signals to assess the complexity of time series. A high sequence self-similarity indicates a low sample entropy value; while a complex sample sequence indicates a large sample entropy value.
For the time series { x ( n ) } = x ( 1 ) , x ( 2 ) , , x ( N ) , the calculation consists following steps:
(1) Construct a sequence of vectors X m ( 1 ) , , X m ( N m + 1 ) with dimension m , where:
X m ( i ) = { x ( i ) , x ( i + 1 ) , , x ( i + m 1 ) } , 1 i N m + 1
X m ( i ) is a vector, composed of m consecutive x values starting from the i th point.
(2) The distance d [ X m ( i ) , X m ( j ) ] between vectors X m ( i ) and X m ( j ) is defined as the maximum distance between the corresponding elements:
d [ X m ( i ) , X m ( j ) ] = max k = 0 , , m 1 ( | x ( i + k ) x ( j + k ) | )
(3) For a specific X m ( i ) , denote the number of distances between X m ( i ) and X m ( j ) ( 1 j N m , j i ) that are less than or equal to r as B i . For 1 j N m , define:
B i m ( r ) = 1 N m 1 B i
(4) Define B ( m ) ( r ) :
B ( m ) ( r ) = 1 N m i = 1 N m B i m ( r )
B ( m ) ( r ) is the probability that two sequences match m points under similar tolerance r.
(5) Expand the dimension to m + 1 . Similar to B i , repeat step (3) and denote the number of conditional distances as A i . Define A i m ( r ) :
A i m ( r ) = 1 N m 1 A i
(6) Define A m ( r ) :
A m ( r ) = 1 N m i = 1 N m A i m ( r )
A m ( r ) is the probability that two sequences match m + 1 points.
The Sample Entropy is defined as:
S a m E n ( m , r ) = lim N { ln [ A m ( r ) B m ( r ) ] }
When N is a finite value, it can be estimated as follows:
S a m E n ( m , r , N ) = ln [ A m ( r ) B m ( r ) ]

2.4. Harmony Search (HS) Algorithm

HS algorithm is a new optimization algorithm enlightened by musical performance process in bands [41]. It is featured by strong global optimization ability, simple structure and few parameters. The HS algorithm is also popular in the optimization of neural networks. All in all, the procedure of the HS algorithm include initialization, a new harmony generation, update of harmony memory (HM), and judgment of algorithm termination conditions. The calculation of HS algorithm is described as follows:
(1) Initialization
Setting parameters. Parameters comprise harmony memory size (HMS), harmony memory considering rate (HMCR), pitch adjusting rate (PAR), and bandwidth (BW). The meanings of parameters are explained as follows:
HMS: the music played by each instrument has a certain range. A solution space is constructed by the music playing range of each instrument, and then a harmony memory is randomly generated by this solution space.
HMCR: it is necessary to take a group of harmony from this harmony memory through a certain probability, and fine tune this group of harmony to get a group of new harmony, and then judge whether this group of new harmony is better than the worst harmony in the harmony memory, and a memory matching rate needs to be randomly generated.
PAR: select a group of harmonies in the harmony memory with a certain probability to fine tune.
BW: a group of harmonies taken from the memory are tuned with a certain probability, which is specified here.
X i = ( x i 1 , x i 2 , , x i d , x i D ) represents a harmony vector with D-dimension. The i th vector of HM is denoted as X i , and each dimension is generated by the following formula:
X i , d = x min , d + ( x max , d x min , d ) * r a n d ( )
where: d [ 1 , D ] and i [ 1 , H M S ] . r a n d ( ) represents a random number, which is uniformly distributed in [ 0 , 1 ] . x min , d and x max , d represent the maximum value and minimum value of the search range for each dimension variable respectively.
(2) New harmony Generation
Generate a new harmony vector X n e w , and each dimension of it is created as follows:
X n e w , d = { X i , d i f   r 1 < H C M R X n e w , d o t h e r w i s e
X n e w , d = { X i , d ± r a n d ( ) * b w d i f   r 2 < P A R X min , d + ( X max , d X min , d ) * r a n d ( ) o t h e r w i s e
where: X i , d is randomly chosen from HM; b w d is the d-dimensional of b w ; r 1 and r 2 are random numbers uniformly distributed in the range [ 0 , 1 ] .
(3) Update of HM
The objective function is used to evaluate the new harmony. Compare the new harmony X n e w with the worst harmony in the HM, if the former is better, then replace the worst harmony in the HM with the new harmony. For the KELM model, the kernel parameter δ and the penalty coefficient C are evaluated.
(4) Judgment of algorithm termination conditions
To determine whether the algorithm satisfies the termination condition, the algorithm outputs the best harmony vector, otherwise it returns to step 2.

2.5. Extreme Learning Machine with Kernel

Traditional learning algorithms, such as back propagation neural network (BPNN) algorithm, have some shortcomings, like slow training speed and easy to fall into local minimum points. Based on the single hidden layer feed-forward neural network (SLFN), the extreme learning machine (ELM) is a new algorithm, which has been widely used in many fields for its good learning ability. The connection weight between the input layer and the hidden layer is randomly generated as well as the threshold of the hidden layer neurons. The algorithm simply set the hidden layer neuron number in training, and a unique global optimal solution can be attained with quick learning speed and excellent generalization capability. The network structure of the extreme learning machine is shown in Figure 1.
In Figure 1, ω i j is the connection weight involving x i ( i = 1 , 2 , , n ) and v j ( j = 1 , 2 , , l ) , while v l is randomly generated as the threshold of the hidden layer; β i j represents the connection weight involving v i ( i = 1 , 2 , l ) and y j ( j = 1 , 2 , , m ) . For a single hidden layer neural network with a ( n l m ) structure, given Q samples ( x i , t i ) , the input data is x i = [ x i 1 , x i 2 , , x i m ] T , the expected output is t i = [ t i 1 , t i 2 , , t i m ] T . The output of the network ELM is [45]:
y j = i = 1 l β i σ i ( ω i x j + v j ) , j = 1 , 2 , , Q
where: σ i is the activation function.
H ω , v , x β = T , β = [ β 1 T β l T ] , T = [ t 1 T t m T ]
where: T is the expected output vector; H ω , v , x is denoted as the output matrix of the hidden layer, which can be expressed as:
H ω , v , x = ( σ ( ω 1 x 1 + v 1 ) σ ( ω l x 1 + v l ) σ ( ω 1 x Q + v 1 ) σ ( ω l x Q + v l ) )
Calculate the solution by means of the Moore-Penrose generalized inverse [46]:
β * = H + T
On the basis of ELM, Huang combines nuclear learning with ELM, replaces random mapping in ELM with nuclear mapping, and proposes KELM algorithm [47]. By introducing the kernel function to obtain better application potential, the definition of kernel function is as follows:
X is the input space, H is the feature space, if there is a mapping φ ( x ) : X H from X to H , so that for all x , y X , function K ( x , y ) = φ ( x ) · φ ( y ) , then K ( x , y ) is the kernel function, φ ( x ) is the mapping function, and φ ( x ) · φ ( y ) is the inner product of x , y mapped to the feature space [48]. This method achieves better application by introducing a kernel function. Defines the kernel matrix as:
{ Q E L M = H H T Q E L M i j = h ( x i ) h ( x j ) = K ( x i , x j )
where: h ( x ) denotes the output function of the hidden layer node; K ( x i , x j ) is a kernel function, namely:
K ( x i , x j ) = exp { x i x j / 2 δ 2 }
According to Equations (19) and (20) he output and output weights of KELM are as follows:
f ( x ) = [ K ( x , x 1 ) M ( x , x Q ) ] ( I C + Ω E L M ) 1 y
β = ( I C + Ω E L M ) 1 y
This paper utilize the HS algorithm to optimize the kernel parameter δ and the penalty coefficient C , which is crucial for the search ability of the KELM algorithm.
Extreme Learning Machine with Kernel (KELM) has strong generalization ability with a high learning speed. Combined with the SLFN of the kernel learning map, this method overcomes the shortcomings of the traditional neural network, which easily falls into the local optimal solution.

3. Wind Power Forecasting Model

3.1. Model Design

This section introduces the development of wind power forecasting model, including model design and model comprehensive evaluation. The flow chart of the CEEMD-SE-HS-KELM model is presented in Figure 2. The hybrid model is composed by three parts: data preprocessing, optimization and prediction phase.
The time series of wind power has strong non-linearity and non-stationarity. CEEMD is used to decompose wind power, reduce the fluctuation of original data, and overcome the modal mixture phenomenon of EMD to obtain multiple IMF components. Then, the complexity of each subsequence is analyzed by SE, and the reconstruction of the subsequence is carried out to reduce the amount of calculation, and HS-KELM model is utilized to forecast the IMF components. Eventually, the final wind power forecasting results are obtained through integrated processing. This paper establishes a CEEMD-SE-HS-KELM model to forecast wind power, as shown below:
(1)
Screen the original data set through correlation analysis to obtain data indicators with high correlation, which can used as the input data of CEEMD-HS-KELM forecasting model;
(2)
Decompose the original wind power sequence X through CEEMD method to obtain n sub sequences components from high frequency to low frequency, which are n 1 intrinsic mode function I M F i ( t ) and one approximately monotonous residual R ( t ) ;
(3)
Utilize SE theory to calculate the complexity of each subsequence and reconstruct the subsequence decomposed by CEEMD;
(4)
HS-KELM model is constructed for each subsequence and the forecast values of each subsequence are obtained;
(5)
Superimpose the wind power forecasting results of each subsequence to obtain the final forecasting results.

3.2. Evaluation Criteria

For the performance test of the proposed model, three evaluation criteria are of great importance, including Root Mean Square Error (RMSE), Mean Absolute Error (MAE) and Determination Factor R2. The calculation methods are as follows:
R M S E = 1 n i = 1 n ( p i p i ' ) 2
M A E = 1 n i = 1 N | p i ' p i |
R 2 = i = 1 n ( p i ' p i ¯ ) 2 i = 1 n ( p i p i ¯ ) 2
where: p i represents the true power value, kW; p i ' represents the anti-normalized power value after CEEMD-SE-HS-KELM output, kW; n represents the data number.

4. Short-Term Wind Power Forecasting

To comprehensively verify the performance of the CEEMD-SE-HS-KELM wind power hybrid forecasting model proposed in this paper, the following case is designed. Taking the measured data of a wind farm in China as an example, short-term wind power is made by the hybrid forecasting model, and compared with many intelligent forecasting models, multiple error indicators are calculated, and the fitting effect and forecasting accuracy of the model are analyzed.
The wind power data of Beijing wind farm in May 2017 is selected as experimental sample. The rated installed capacity of the wind farm is 49 MW and the sampling frequency is one point per 5 min. The data of 17 consecutive days from May 15th to 31st are taken, totaling 4608 sampling points, of which the first 4322 sampling points are training sets of the forecasting model. The last 288 sampling points are the test set, the original data is in the “Supplementary Materials”.

4.1. Data Set Screening

The original data sets are divided into two categories, one is historical meteorological data, including historical wind power data and historical weather data such as wind speed and wind direction. According to the Pearson correlation coefficients analysis of the original data sets, the observed correlations degrees between the meteorological data and wind power are shown in Table 1.
According to the Pearson correlation coefficients, the coupling degree between wind speed, wind direction and wind power are relatively high, so the wind speeds of different heights, the wind direction of 30 m, the wind direction of 50 m are selected as the input variables of the model. The original data is converted into a number between [0, 1], the purpose of which is to eliminate the magnitude difference between the various dimensions of data and avoid the influence on the forecasting result caused by the large magnitude difference between input and output data. The normalization method is as follows:
x i = x i x min x max x min
where: x min and x max are respectively the minimum and the maximum value in the sequence, x i is the initial input data, and x i is the normalized data.

4.2. CEEMD Decomposition of Wind Power Sequence and Subsequence Reconstruction

The original wind power data is decomposed by CEEMD decomposition method, and 12 IMF components and 1 residual component are obtained. The result of decomposition is shown in Figure 3. After CEEMD decomposition processing, the characteristic changes of the signal are extracted from high frequency to low frequency, and the components are relatively stable.
Using the HS-KELM model to directly forecast each subsequence will increase the amount of computation. SE method is used to calculate the complexity of each subsequence, and the results are shown in Figure 4.
The criterion of SE for IMF sequence classification is about 0.2 times of the original sequence standard deviation [49]. By calculating the standard deviation of subsequence complexity, the standard deviation is 0.25, so the similarity difference of subsequence reconstruction is obtained. The value is 0.12 and is rounded down to 0.1. As shown in Figure 5, the sample entropy of IMF2 and IMF3 subsequences has small difference, and the two subsequences can be classified into one class, and integrated and reconstructed as a new subsequence input into the HS-KELM for training and forecasting. All subsequences are categorized and the results are presented in Table 2 and the reconstructed results are presented in Figure 5.

4.3. Wind Power Forecasting by CEEMD-SE-HS-KELM Model

Following decomposition and reconstruction for the sequences, 8 HS-KELM forecasting models are constructed, and the above 8 subsequences are respectively trained and forecast. For the HS-KELM model, the common values of 3 parameters H M C R , P A R and B W of HS algorithm is [ 0.63 , 0.99 ] , [ 0.01 , 0.73 ] and [ 0.0004 , 0.3 ] respectively [41,50], the maximum iteration limit is detailed in reference [51], and the objective function is set as RMSE. Therefore, H M S = 100 , H M C R = 0.9 , P A R = 0.35 , B W = 0.25 , the number of new harmony vectors generated each time is 10.
For the HS-KELM model, the number of hidden layer neurons is set to 20, the data of a certain day is randomly selected and the wind power is forecast. RMSE is selected as the objective function, and the number of terminations is set as 100. As shown in Figure 6, the target function value tends to be gentle after 21 iterations, and ends at the 100th generation. The convergence speed is faster, and it is more consistent with the actual power curve. Therefore, the forecasting model structure of HS-KELM model is set as “6-20-1”. Train and forecast the eight subsequences respectively, and the final forecasting results of each subsequence are integrated to get the wind power forecasting output value, as shown in Figure 7.

4.4. Comparative Analysis of Forecasting Models

To confirm the validity of this model, ELM, KELM, HS-KELM, EMD-SE-HS-KELM, CEEMD-SE-HS-KELM models are constructed for comparative analysis, and they are respectively named configuration #1–#5. The KELM structure is in the form of “6-20-1”, and the maximum number of training time is 500. To avoid the influence of randomness on the forecasting results, the average value of each model is taken after 50 independent runs and compared with the actual wind power. The forecasting results are presented in Figure 8.
To compare and evaluate forecasting accuracy of different models, determination coefficient R2 and model error indexes such as Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) are calculated respectively, as shown in Table 3 and Figure 9.
It can be seen from Table 3 and Figure 8, Figure 9:
(1)
Comparing EMD-SE-HS-KELM and CEEMD-SE-HS-KELM, the RMSE and MAE of the latter are improved by 53.63% and 24.36% respectively compared with EMD-SE-HS-KELM, which indicates that the hybrid model of data preprocessing based on CEEMD-HS has better processing effect.
(2)
Comparing HS-KELM and EMD-SE-HS-KELM, the RMSE and MAE of the latter are improved by 52.10% and 86.58%, respectively, compared with HS-KELM, which indicates that for non-stationary wind power series, pre-processing can effectively eliminate noise, ensure data quality and improve forecasting accuracy.
(3)
Comparing KELM and HS-KELM, the RMSE and MAE of the latter are improved by 37.67% and 45.23%, respectively, compared with KELM, which indicates that the parameters of KELM algorithm are optimized by HS algorithm, which effectively improves the search ability and the forecasting accuracy.
(4)
Compared with ELM, the RMSE and MAE of the KELM are improved by 25.02% and 24.06% respectively compared with ELM, which indicates that the forecasting accuracy of KELM algorithm is better than ELM model, and KELM has a stronger generalization ability.
It can be seen from Figure 8, that a single KELM and ELM forecasting model can only roughly reflect the wind power trend, but the forecasting value at each time point is quite different from the actual. The CEEMD-SE-HS-KELM hybrid forecasting model constructed in this paper has a good fitting effect between the forecasting result and the actual value at each time point, and its forecasting accuracy is higher.

5. Conclusions

Aiming at short-term forecasting of non-linear and unsteady wind power time series, a hybrid forecasting model consisting of CEEMD-SE and HS optimized KELM is proposed in this paper. Firstly, the Pearson correlation coefficient is used to screen the input data to reduce data redundancy. Secondly, the combined data preprocessing strategy of CEEMD-SE is utilized to process the data. The wind power time series is decomposed and reconstructed to eliminate data noise and reduce the computational load. After that, the KELM model optimized by HS is used to forecast each subsequence after reconstructing, and the final wind power forecasting value is obtained after integration processing. Finally, a specific wind farm in China is taken as an example, the CEEMD-SE-HS-KELM, EMD-SE-HS-KELM, HS-KELM, KELM and ELM models are established to forecast wind power respectively. The case study shows that:
(1)
KELM model has higher forecasting accuracy than ELM model, and has broad application prospects in wind power forecasting.
(2)
Compared with the single KELM model, HS can optimize the kernel parameters and penalty function of KELM to obtain higher forecasting accuracy, which indicates that HS-KELM model has stronger global search ability and more stable forecasting performance.
(3)
Compared with EMD-SE, the data preprocessing strategy based on CEMD-SE has better performance and effectively improves the forecasting accuracy. The hybrid model proposed in this paper can be well applied to short-term wind power forecasting.

Supplementary Materials

The following are available online at https://www.mdpi.com/2227-9717/7/11/843/s1.

Author Contributions

All of the authors have contributed to this research. K.W., L.S. and H.Z. collected the data and wrote this paper; D.N. provided professional guidance; J.L. collected the data, G.D. and X.X. revised this manuscript. All authors have approved the submitted manuscript.

Funding

This work was supported by the 2018 Key Projects of Philosophy and Social Sciences Research, Ministry of Education, China (grant number 18JZD032); 111 Project, (grant number B18021); Natural Science Foundation of China (grant number 71804045).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

NWPNumerical weather prediction
ANNArtificial neural network
SVMSupport vector machine
LSSVMLeast squares support vector machine
PCAPrincipal component analysis
WTWavelet transform
EMDEmpirical mode decomposition
EWTEmpirical wavelet transform
CEEMDComplementary ensemble empirical mode decomposition
WASWavelet neural network
DWTDiscrete wavelet transform
GARCHGeneralized autoregressive conditional heteroscedastic
MkRVRMulti-kernel relevance vector regression
IMFIntrinsic mode functions
MLFFNNMultilayer feed-forward neural network
SVRSupport vector regression
RBFRadial basis function
ANFISAdaptive neuro-fuzzy inference system
PSOParticle swarm optimization
VMDVariational mode decomposition
BABat algorithm
ELMExtreme learning machine
CSOCrisscross optimization algorithm
KELMExtreme learning machine with kernel
HSHarmony search
SESample entropy
HMHarmony memory
HMSHarmony memory size
HMCRHarmony memory considering rate
PCRPitch adjusting rate
BWBandwidth
BPNNBack propagation neural network
RMSERoot mean square error
MAEMean absolute error
R2Determining factor

References

  1. Cai, W.; Lai, K.-H.; Liu, C.; Wei, F.; Ma, M.; Jia, S.; Jiang, Z.; Lv, L. Promoting sustainability of manufacturing industry through the lean energy-saving and emission-reduction strategy. Sci. Total Environ. 2019, 665, 23–32. [Google Scholar] [CrossRef] [PubMed]
  2. Cai, W.; Liu, C.; Lai, K.-H.; Li, L.; Cunha, J.; Hu, L. Energy performance certification in mechanical manufacturing industry: A review and analysis. Energy Convers. Manag. 2019, 186, 415–432. [Google Scholar] [CrossRef]
  3. Hong, C.; Lin, W.M.; Wen, B.Q. Review of wind speed and wind power prediction methods for wind farms. Power Syst. Clean Energy. 2011, 27, 60–66. [Google Scholar]
  4. Ren, Y.; Suganthan, P.N.; Srikanth, N. A novel empirical mode decomposition with support vector regression for wind speed forecasting. IEEE Trans. Neural Netw. Learn. Syst. 2016, 27, 1793–1798. [Google Scholar] [CrossRef] [PubMed]
  5. Zhao, X.; Wang, C.; Su, J.; Wang, J. Research and application based on the swarm intelligence algorithm and artificial intelligence for wind farm decision system. Renew. Energy 2019, 134, 681–697. [Google Scholar] [CrossRef]
  6. National Energy Administration. China Power Industry Annual Development Report: Renewable Energy Added Capacity Accounts for Over 50%. Available online: http://www.nea.gov.cn/2019-06/19/c_138155207.htm (accessed on 19 June 2019).
  7. Zhang, C.; Zhou, J.; Li, C.; Fu, W.; Peng, T. A compound structure of ELM based on feature selection and parameter optimization using hybrid backtracking search algorithm for wind speed forecasting. Energy Convers. Manag. 2017, 143, 360–376. [Google Scholar] [CrossRef]
  8. Tasnim, S.; Rahman, A.; Oo, A.M.T.; Haque, M.E. Wind power prediction in new stations based on knowledge of existing stations: A cluster based multi source domain adaptation approach. Knowl. Based Syst. 2018, 145, 15–24. [Google Scholar] [CrossRef]
  9. Wang, J.; Du, P.; Niu, T.; Yang, W. A novel hybrid system based on a new proposed algorithm-Multi-Objective Whale Optimization Algorithm for wind speed forecasting. Appl. Energy 2017, 208, 344–360. [Google Scholar] [CrossRef]
  10. Li, H.; Wang, J.; Lu, H.; Guo, Z. Research and application of a combined model based on variable weight for short term wind speed forecasting. Renew. Energy 2018, 116, 669–684. [Google Scholar] [CrossRef]
  11. Feng, S.L.; Wang, W.S.; Liu, C.; Dai, H.Z. Study on the physical approach to wind power prediction. Proc. CSEE 2010, 30, 1–6. [Google Scholar]
  12. Yang, W.; Wang, J.; Lu, H.; Niu, T.; Du, P. Hybrid wind energy forecasting and analysis system based on divide and conquer scheme: A case study in China. J. Clean. Prod. 2019, 222, 942–959. [Google Scholar] [CrossRef]
  13. Zhang, W.; Qu, Z.; Zhang, K.; Mao, W.; Ma, Y.; Fan, X. A combined model based on CEEMDAN and modified flower pollination algorithm for wind speed forecasting. Energy Convers. Manag. 2017, 136, 439–451. [Google Scholar] [CrossRef]
  14. Jiang, P.; Yang, H.; Heng, J. A hybrid forecasting system based on fuzzy time series and multi-objective optimization for wind speed forecasting. Appl. Energy 2019, 235, 786–801. [Google Scholar] [CrossRef]
  15. Hao, Y.; Tian, C. The study and application of a novel hybrid system for air quality early-warning. Appl. Soft Comput. 2019, 74, 729–746. [Google Scholar] [CrossRef]
  16. Wang, Y.; Wang, J.; Wei, X. A hybrid wind speed forecasting model based on phase space reconstruction theory and Markov model: A case study of wind farms in northwest China. Energy 2015, 91, 556–572. [Google Scholar] [CrossRef]
  17. Sun, W.; Liu, M. Wind speed forecasting using FEEMD echo state networks with RELM in Hebei, China. Energy Convers. Manag. 2016, 114, 197–208. [Google Scholar] [CrossRef]
  18. Lahouar, A.; Slama, J.B.H. Hour-ahead wind power forecast based on random forests. Renew. Energy 2017, 109, 529–541. [Google Scholar] [CrossRef]
  19. Li, C.; Zhu, Z. Research and application of a novel hybrid air quality early-warning system: A case study in China. Sci. Total Environ. 2018, 626, 1421–1438. [Google Scholar] [CrossRef] [PubMed]
  20. Xiao, L.; Shao, W.; Yu, M.; Ma, J.; Jin, C. Research and application of a hybrid wavelet neural network model with the improved cuckoo search algorithm for electrical power system forecasting. Appl. Energy 2017, 198, 203–222. [Google Scholar] [CrossRef]
  21. Li, G.; Shi, J.; Zhou, J. Bayesian adaptive combination of short-term wind speed forecasts from neural network models. Renew. Energy 2011, 36, 352–359. [Google Scholar] [CrossRef]
  22. Liu, D.; Niu, D.; Wang, H.; Fan, L. Short-term wind speed forecasting using wavelet transform and support vector machines optimized by genetic algorithm. Renew. Energy 2014, 62, 592–597. [Google Scholar] [CrossRef]
  23. Hu, J.; Wang, J.; Ma, K. A hybrid technique for short-term wind speed prediction. Energy 2015, 81, 563–574. [Google Scholar] [CrossRef]
  24. Zhou, J.; Shi, J.; Li, G. Fine tuning support vector machines for short-term wind speed forecasting. Energy Convers. Manag. 2011, 52, 1990–1998. [Google Scholar] [CrossRef]
  25. Li, H.; Wang, J.; Li, R.; Lu, H. Novel analysis-forecast system based on multi-objective optimization for air quality index. J. Clean. Prod. 2019, 208, 1365–1383. [Google Scholar] [CrossRef]
  26. Huang, G.-B.; Zhou, H.; Ding, X.; Zhang, R. Extreme learning machine for regression and multiclass classification. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2011, 42, 513–529. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, Y.; Zhang, C.; Zhao, Y.; Gao, S. Wind speed prediction with RBF neural network based on PCA and ICA. J. Electr. Eng. 2018, 69, 148–155. [Google Scholar] [CrossRef]
  28. Huang, Y.; Shen, L.; Liu, H. Grey relational analysis, principal component analysis and forecasting of carbon emissions based on long short-term memory in China. J. Clean. Prod. 2019, 209, 415–423. [Google Scholar] [CrossRef]
  29. Qi, J.; Xu, C.Z.; Liu, Y.Q.; Han, S.; Li, L. Short-term wind power prediction method based on wind speed cloud model in similar day. Autom. Electr. Power Syst. 2018, 42, 53–59. [Google Scholar]
  30. Gilles, J. Empirical wavelet transform. IEEE Trans. Signal Process. 2013, 61, 3999–4010. [Google Scholar] [CrossRef]
  31. Naik, J.; Satapathy, P.; Dash, P. Short-term wind speed and wind power prediction using hybrid empirical mode decomposition and kernel ridge regression. Appl. Soft Comput. 2018, 70, 1167–1188. [Google Scholar] [CrossRef]
  32. Du, P.; Wang, J.Z.; Yang, W.D.; Niu, T. A novel hybrid model for short-term wind power forecasting. Appl. Soft Comput. 2019, 80, 93–106. [Google Scholar] [CrossRef]
  33. Jiang, P.; Wang, Y.; Wang, J. Short-term wind speed forecasting using a hybrid model. Energy 2017, 119, 561–577. [Google Scholar] [CrossRef]
  34. Jiang, Y.; Huang, G.; Peng, X.; Li, Y.; Yang, Q. A novel wind speed prediction method: Hybrid of correlation-aided DWT, LSSVM and GARCH. J. Wind Eng. Ind. Aerodyn. 2018, 174, 28–38. [Google Scholar] [CrossRef]
  35. Fei, S.W. A hybrid model of EMD and multiple-kernel RVR algorithm for wind speed prediction. Int. J. Electr. Power Energy Syst. 2016, 78, 910–915. [Google Scholar] [CrossRef]
  36. Khosravi, A.; Koury, R.; Machado, L.; Pabon, J. Prediction of wind speed and wind direction using artificial neural network, support vector regression and adaptive neuro-fuzzy inference system. Sustain. Energy Technol. Assess. 2018, 25, 146–160. [Google Scholar] [CrossRef]
  37. Wu, Q.; Lin, H.X. Short-term wind speed forecasting based on hybrid variational mode decomposition and least squares support vector machine optimized by bat algorithm model. Sustainability 2019, 11, 652. [Google Scholar] [CrossRef] [Green Version]
  38. Jiang, Y.; Huang, G.; Yang, Q.; Yan, Z.; Zhang, C. A novel probabilistic wind speed prediction approach using real time refined. Energy Convers. Manag. 2019, 185, 758–773. [Google Scholar] [CrossRef]
  39. Tian, C.; Hao, Y.; Hu, J. A novel wind speed forecasting system based on hybrid data preprocessing and multi-objective optimization. Appl. Energy 2018, 231, 301–319. [Google Scholar] [CrossRef]
  40. Yin, H.; Dong, Z.; Chen, Y.; Ge, J.; Lai, L.L.; Vaccaro, A.; Meng, A. An effective secondary decomposition approach for wind power forecasting using extreme learning machine trained by crisscross optimization. Energy Convers. Manag. 2017, 150, 108–121. [Google Scholar] [CrossRef]
  41. Geem, Z.W.; Kim, J.H.; Loganathan, G.V. A new heuristic optimization algorithm: Harmony search. Simulation 2001, 76, 60–68. [Google Scholar] [CrossRef]
  42. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. A: Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  44. Yeh, J.R.; Shieh, J.S.; Huang, N.E. Complementary ensemble empirical mode decomposition: A novel noise enhanced data analysis method. Adv. Adapt. Data Anal. 2010, 2, 135–156. [Google Scholar] [CrossRef]
  45. Lin, T.; Cai, R.Q.; Zhang, L.; Yang, X.; Liu, G.; Liao, W.J. Prediction intervals forecasts of wind power based on IBA-KELM. Renew. Energy Resour. 2018, 7, 23. [Google Scholar]
  46. Bapat, K.M.P.B. The Generalized Moore-Penrose Inverses. Linear Algebra Appl. 1992, 165, 59–69. [Google Scholar]
  47. Huang, G.-B. An insight into extreme learning machines: Random neurons, random features and kernels. Cogn. Comput. 2014, 6, 376–390. [Google Scholar] [CrossRef]
  48. Li, D. Kernel Function. In Encyclopedia of Microfluidics and Nanofluidics; Springer: Boston, MA, USA, 2008. [Google Scholar] [CrossRef]
  49. Hornero, R. Optimal parameters study for sample entropy-based atrial fibrillation organization analysis. Comput. Methods Programs Biomed. 2010, 99, 124–132. [Google Scholar]
  50. Kim, J.H.; Geem, Z.W.; Kim, E.S. Parameter estimation of the nonlinear muskingum model using harmony search. J. Am. Water Resour. Assoc. 2001, 37, 1131–1138. [Google Scholar] [CrossRef]
  51. Omran Mahamed, G.H.; Mahdavi, M. Global-best harmony search. Appl. Math. Comput. 2008, 198, 643–656. [Google Scholar]
Figure 1. Extreme learning machine (ELM) training model with ( n l m ) structure.
Figure 1. Extreme learning machine (ELM) training model with ( n l m ) structure.
Processes 07 00843 g001
Figure 2. The flow chart of CEEMD-SE-HS-KELM model.
Figure 2. The flow chart of CEEMD-SE-HS-KELM model.
Processes 07 00843 g002
Figure 3. Decomposition results of CEEMD.
Figure 3. Decomposition results of CEEMD.
Processes 07 00843 g003
Figure 4. Sample entropy of each subsequence.
Figure 4. Sample entropy of each subsequence.
Processes 07 00843 g004
Figure 5. Wind power sub-sequences processed by CEEMD-SE.
Figure 5. Wind power sub-sequences processed by CEEMD-SE.
Processes 07 00843 g005
Figure 6. Iteration curve.
Figure 6. Iteration curve.
Processes 07 00843 g006
Figure 7. Wind power forecasting results.
Figure 7. Wind power forecasting results.
Processes 07 00843 g007
Figure 8. Wind power forecasting results (5 min). (a) configuration #1: ELM model; (b) configuration #2: KELM model; (c) configuration #3: HS-KELM model; (d) configuration #4: EMD-SE-HS-KELM model; (e) configuration #5: CEEMD-SE-HS-KELM model; (f) comparison of forecasting values of different models.
Figure 8. Wind power forecasting results (5 min). (a) configuration #1: ELM model; (b) configuration #2: KELM model; (c) configuration #3: HS-KELM model; (d) configuration #4: EMD-SE-HS-KELM model; (e) configuration #5: CEEMD-SE-HS-KELM model; (f) comparison of forecasting values of different models.
Processes 07 00843 g008
Figure 9. Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) and R2 of different forecasting models.
Figure 9. Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) and R2 of different forecasting models.
Processes 07 00843 g009
Table 1. Pearson correlation coefficient analysis results.
Table 1. Pearson correlation coefficient analysis results.
IndicatorsCorrelationIndicatorsCorrelationIndicatorsCorrelation
wind speed10 m0.788wind direction10 m−0.391temperature0.227
30 m0.79630 m−0.308humidity−0.51
50 m0.77750 m−0.025rainfall0.29
70 m0.76470 m0.289pressure−0.37
hub height0.764hub height0.289-
Table 2. Results of the new subsequence with merged intrinsic mode functions (IMF) components.
Table 2. Results of the new subsequence with merged intrinsic mode functions (IMF) components.
New Subsequence NumberInitial Subsequence NumberNew Subsequence NumberInitial Subsequence Number
1157
22,368,9
34710,11,12
45,68R
Table 3. The calculation results of model evaluation index.
Table 3. The calculation results of model evaluation index.
Evaluation IndexELMKELMHS-KELMEMD-SE-HS-KELMCEEMD-SE-HS-KELM
Configuration#1#2#3#4#5
RMSE20.8415.639.744.672.16
MAE9.297.053.860.520.39
R21.121.050.961.031.01

Share and Cite

MDPI and ACS Style

Wang, K.; Niu, D.; Sun, L.; Zhen, H.; Liu, J.; De, G.; Xu, X. Wind Power Short-Term Forecasting Hybrid Model Based on CEEMD-SE Method. Processes 2019, 7, 843. https://doi.org/10.3390/pr7110843

AMA Style

Wang K, Niu D, Sun L, Zhen H, Liu J, De G, Xu X. Wind Power Short-Term Forecasting Hybrid Model Based on CEEMD-SE Method. Processes. 2019; 7(11):843. https://doi.org/10.3390/pr7110843

Chicago/Turabian Style

Wang, Keke, Dongxiao Niu, Lijie Sun, Hao Zhen, Jian Liu, Gejirifu De, and Xiaomin Xu. 2019. "Wind Power Short-Term Forecasting Hybrid Model Based on CEEMD-SE Method" Processes 7, no. 11: 843. https://doi.org/10.3390/pr7110843

APA Style

Wang, K., Niu, D., Sun, L., Zhen, H., Liu, J., De, G., & Xu, X. (2019). Wind Power Short-Term Forecasting Hybrid Model Based on CEEMD-SE Method. Processes, 7(11), 843. https://doi.org/10.3390/pr7110843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop