Next Article in Journal
Photovoltaic Array Fault Detection by Automatic Reconfiguration
Previous Article in Journal
A High-Power-Density Single-Phase Rectifier Based on Three-Level Neutral-Point Clamped Circuits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Hybrid Model Based on Extreme Learning Machine, k-Nearest Neighbor Regression and Wavelet Denoising Applied to Short-Term Electric Load Forecasting

School of Mathematics and Statistics, Lanzhou University, Lanzhou 730000, Gansu, China
*
Author to whom correspondence should be addressed.
Energies 2017, 10(5), 694; https://doi.org/10.3390/en10050694
Submission received: 24 March 2017 / Revised: 23 April 2017 / Accepted: 9 May 2017 / Published: 16 May 2017
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
Electric load forecasting plays an important role in electricity markets and power systems. Because electric load time series are complicated and nonlinear, it is very difficult to achieve a satisfactory forecasting accuracy. In this paper, a hybrid model, Wavelet Denoising-Extreme Learning Machine optimized by k-Nearest Neighbor Regression (EWKM), which combines k-Nearest Neighbor (KNN) and Extreme Learning Machine (ELM) based on a wavelet denoising technique is proposed for short-term load forecasting. The proposed hybrid model decomposes the time series into a low frequency-associated main signal and some detailed signals associated with high frequencies at first, then uses KNN to determine the independent and dependent variables from the low-frequency signal. Finally, the ELM is used to get the non-linear relationship between these variables to get the final prediction result for the electric load. Compared with three other models, Extreme Learning Machine optimized by k-Nearest Neighbor Regression (EKM), Wavelet Denoising-Extreme Learning Machine (WKM) and Wavelet Denoising-Back Propagation Neural Network optimized by k-Nearest Neighbor Regression (WNNM), the model proposed in this paper can improve the accuracy efficiently. New South Wales is the economic powerhouse of Australia, so we use the proposed model to predict electric demand for that region. The accurate prediction has a significant meaning.

Highlights:

  • A novel hybrid model named ELM-WA-KNN is proposed for electric load forecasting in New South Wales, Australia.
  • The wavelet analysis (WA) is introduced to eliminate the noise of the electric load time series.
  • k-Nearest Neighbor regression (KNN) is used to get the input-output relationship in the hybrid model.
  • The kernel function of KNN is established by extreme learning machine (ELM).
  • The proposed ELM-WA-KNN model has the best performance among all the considered models.

1. Introduction

Electric load prediction is one of the major tasks in power management departments which carry out electric power dispatching, usage and planning. Improving the accuracy of load forecasting is helpful for planning power management and is advantageous to the arrangement of power grid operations. Meanwhile, it is also advantageous for the reduction of power generation costs and is beneficial to establish the power for construction plans. Therefore, electric load forecasting has become an important part of power system management modernization.
In recent years, a large number of techniques for time series forecasting have been used in electric load forecasting. In traditional predictive methods, the linear regression model and grey model have been widely used. Goia et al. proposed a linear regression model using heating demand data to forecast the short-term peak load in a district-heating system [1]. Zhou et al. presented a trigonometric grey prediction approach by combing the traditional GM (1, 1) with the trigonometric residual modification technique for forecasting electricity demand [2]. Akay et al. proposed GPRM to predict Turkey’s total and industrial electricity consumption [3]. Meanwhile, artificial neural networks and support vector machines have a variety of applications in electric load forecasting. Chang et al. proposed the so-called EEuNN framework by adopting a weighted factor to calculate the importance of each factor among the different rules to predict monthly electricity demand in Taiwan [4]. Kavaklioglu used SVR to model and predict Turkey’s electricity consumption [5]. Wang et al. presented a combined ε-SVR model considering seasonal proportions based on development trends from history data to forecast the short-term electricity demand [6]. Other predictive techniques have also been proposed to deal with the electric load forecasting problem. Kucukali et al. attempted to forecast Turkey’s short-term gross annual electric demand by applying fuzzy logic methodology based on the economical, political and electricity market conditions of the country [7]. Dash et al. presented the development of a hybrid neural network to model a fuzzy expert system for time series forecasting of electric loads [8]. Taylor showed that for predictions up to a day-ahead the triple seasonal methods outperform the double seasonal methods in predicting electricity demand [9].
However, no matter which techniques are utilized to predict electric load, with only a single model it is difficult to achieve high precision because of the various shortcomings in the models [10]. Therefore, hybrid models are built with different combinations of data mining technology to improve the prediction accuracy in load forecasting research. Wu et al. proposed a Particle Swarm Optimization-Supporter Vector Machine (PSO-SVM) model based on cluster analysis techniques and data accumulation pretreatment in short-term electric load forecasting [11]. Chen et al. proposed a new electric load forecasting model by hybridizing FTS and GHSA with LSSVM [12]. Huang presented a SVR-based load forecasting model which hybridized the chaotic mapping function and QPSO with SVM [13]. Zhang et al. successfully established a novel model for electric load forecasting by combination of SSA, SVM and CS [14]. Hence, it is apparent that lots of novel hybrid models combined with different data mining techniques could improve the prediction precision in the research of load forecasting.
In this paper we also propose an outstanding hybrid model to predict electric loads. In non-parametric model research, KNN regression has been widely used in time series prediction. The KNN regression is one of the historical approximation methods in machine learning. The main idea of the algorithm is that we get the output by calculating the degree of similarity between the independent variables [15]. Poloczek proposed the KNN regression as a geo-imputation preprocessing step for pattern-label-based short-term wind prediction of spatio-temporal wind data sets [16]. Ban proposed a new multivariate regression approach for financial time series based on knowledge shared from referential nearest neighbors [17]. Hu proposed a conjunction model named EMD-KNN which was based on an EMD and KNN regressive model for forecasting annual average rainfall [18]. Li established a novel hybrid BAMO with KNN to predict apoptosis protein sequences using statistical factors and dipeptide composition [19]. Wang selected k-nearest neighbors to correct the commonly used precipitation data on the Qinghai-Tibetan plateau by establishing the relationship between daily precipitation and environmental as well as other meteorological factors [20]. Meanwhile in the neural network research field, Huang established an innovative algorithm called extreme learning machine (ELM) which was based on a traditional single-hidden-layer feed-forward neural network [21]. Generally speaking, ELM not only reduces the training time of neural networks but also has a great predictive performance. Meanwhile the speed of ELM is faster than traditional machine learning algorithms. Ma built an adaptive prediction model based on ELM to predict traffic flows [22]. Masri targeted predicting the functional properties of soil samples by establishing a hybrid model based on a Savitzky-Golay filter for preprocessing and ELM for obtaining non-linear relationship [23]. Liao used ELM and economic theory to study stock price forecasting and the results showed that ELM had the highest prediction accuracy [24]. Zhang researched the short-term prediction of wind by a proposed model based on wavelet decomposition and ELM [25]. Shamshirband predicted the horizontal global solar radiation by using and extreme learning machine method and the comparative results clearly specified that ELM can provide reliable predictions with better precision compared to the traditional techniques [26]. In the field of signal processing, compared to traditional methods such as EMD and EEMD, the wavelet transform (WT) has a great ability to obtain the characteristics of data from the time and frequency domains [27]. Hence we can eliminate the noise in signals and grab the main information using WT. Wang targeted forecasting future stock prices by utilizing wavelet analysis to denoise the time series applying a neural network to obtain the nonlinear relationships [28]. Wang proposed a novel approach for short-term load forecasting by applying wavelet denoising in a combined model that is a hybrid of SARIMA and a neural network [29]. Qin predicted chaotic time series based on wavelet denoising, phase space reconstruction and LSSVM, and the proposed model had better performance than the traditional models [30]. Abbaszadeh proposed a new hybrid model based on a wavelet denoising technique to denoise hydrological time series and applied ANN to acquire the best prediction of hydrological data [31].
In this paper, because the Australian region of New South Wales (abbreviated as NSW) is an active economic region in the Asia-Pacific region and it has a large population, its economic development and peoples’ lives have become inextricably tied to electricity, so the electricity department is concerned with ensuring sufficient power production. Accurate electric load forecasting could be particularly important for this, as surplus electricity production will lead to environmental pollution and resource wastage. Thus, effective electric load forecasting could be a useful indicator for decision-makers, and effective early warning of an increase in electricity demand is important to ensure supply-demand balance. Based on these facts, we propose a new hybrid model based on the wavelet denoising technique, k-nearest neighbor regression and extreme learning machine to forecast the short-term electricity load in New South Wales.
This paper mainly consists of three parts: (1) the first part introduces the wavelet denoising technique, the k-nearest neighbor regression (KNN), extreme learning machine (ELM) and the establishment of the proposed hybrid model; (2) the second part presents the data set of the experiment, the evaluation criteria of models, the predictive values of the electric load and the analysis of comparative models; (3) the third part is a summary of the proposed model to illustrate the great predictive performance we could achieve for electric loads through the proposed ELM-WA-KNN hybrid model (EWKM).

2. Materials and Methods

The proposed hybrid model presented in this paper for short-term electric load forecasting mainly consists of three basic data mining technologies: wavelet denoising technique, k-nearest neighbor regression (KNN) and extreme learning machine (ELM). Detailed descriptions are given below.

2.1. Wavelet Denoising Technique

The wavelet transform proposed by Mallat has become one of the strongest mathematical tools for providing a time-frequency representation of an analyzed signal in the time series domain. Detailed coefficients are produced by high-pass filters and approximation series are produced by low-pass filters [32]. With the development of the wavelet analysis theory, the dyadic wavelet transform (DWT) for discrete time series y t is as shown in Equation (1):
W m , n = 2 m 2 t = 0 N 1 ψ ( t 2 m n 2 m ) y t
where s = 2 m and τ = 2 m n are the scale and location of the discrete wavelet, respectively. N stands for the integer power of 2, and ψ ( ) for the wavelet function. By the transform method, DWT could eliminate the white nose of the time series and acquire the useful information of the time series on a different scale.

2.2. Extreme Learning Machine (ELM)

Extreme learning machine (ELM) is a type of the single-hidden-layer feed-forward neural network (SLFN) which cannot adjust the parameters of the neural network, but ELM has the threshold of the hidden layer and the weights between each layers. The model of the ELM is as shown in Figure 1.
Given a data set ν = { ( x i , y i ) | i = 1 , 2 , , N ;   x i R n ;   y i R m } , the output of the ELM is:
f ( x i ) = j = 1 L β j g ( α j x i + b j )
where i = 1 , 2 , , N ; b j stands for the learning parameters of the hidden layers and α j = [ α 1 j , α 2 j , , α r j ] T represents the weights between the input layer and hidden layer. Meanwhile the β j = [ β j 1 , β j 2 , , β j m ] T ( j = 1 , 2 , , L ) is the weights between the hidden layer and output layer and g ( x ) is the activation function of the hidden layer. Specifically, the radial basis function is selected as the active function in the hidden node in the experiment. So we could obtained the formula of the ELM as follows:
H β = Y H = [ g ( α 1 x 1 + b 1 ) g ( α L x 1 + b L ) g ( α 1 x N + b 1 ) g ( α L x N + b L ) ] N × L
where the weight matrix between the hidden layer and output layer is β = [ β 1 , β 2 , .. , β L ] T and Y = [ y 1 , y 2 , .. , y N ] T stands for the ELM output. It is apparent that H represents the result of the hidden layer. Here, two important ELM theorems by Liang must to be mentioned [33]:
Theorem 1.
If SLFN with L additive nodes and with activation function g ( x ) which is differentiable in any interval of R is given, then the output matrix of hidden layer H is invertible and H β T = 0 .
Theorem 2.
If small positive value ε > 0 and activation function g ( x ) : R R which is differentiable in any interval is given, then there exists L N such that N arbitrary distinct input vectors randomly produced based upon any continuous probability distribution with probability one.
Based on Theorems 1 and 2, we could solve the equation by the least squares method and the result is:
β = H 1 Y
where H 1 is the Moore-Penrose pseudo inverse of the hidden layer.

2.3. k-Nearest Neighbor Regression (KNN)

KNN is a non-parametric technology which derives from pattern recognition studies [34]. With the development of the study of nonlinear dynamics, many researchers have utilized KNN to frequently predict time series because the algorithm has a great ability to get the nonlinear properties of a time series. The main idea of KNN is that the similarity (neighborhood) between the independent variable of the predictors is used and the independent variable in the historical observations is calculated to acquire the best estimators for the predictor [15].
KNN applies a metric on the predictors to seek the set of k past nearest neighbors in the historical data set for the current condition. To deal with the regressive problem, Lall and Sharma proposed the kernel function and we could get the result of the KNN regression as follows [35]:
Y i = f ( Y i ( 1 ) , Y i ( 2 ) , , Y i ( k ) )
where Y i represents the value of the prediction; Y i ( j ) stands for the magnitude of nearest neighbor j in the above formula. It must be noticed that j is the order of the nearest neighbors based on the distance from the current condition i (j = 1 to k ). The similarity between predictor and historical label is depended on the distance as following formula:
r i j = t = 1 q ( d j t d i t ) 2
where d j t is the t th independent variable of Y j ; d i t stands for the t th independent variable of Y i and q represents the number of independent variables in the formula.
Establishing the kernel function of KNN is therefore the main concern to predict the time series. Many researchers have tried different methods to build fitter kernel functions in dealing with the regression. However almost all the kernel functions are established based on the linear relationship and these methods have failed to get the nonlinear property, so in this paper the kernel function is obtained by ELM, which is a popular data mining technique, to search for the nonlinear properties.

2.4. The Proposed Hybrid Model

In this section, the hybrid model of ELM-WA-KNN is proposed as follows:
Suppose the historical data set is O H = { oh 1 , oh 2 , , oh p } , the training data set is O T R = { ox 1 , ox 2 , , ox p } and the testing data set is O T E = { ox n + 1 , ox n + 2 , , ox n + m } .
Firstly, the electric load time series including the historical data set, the training data set and the testing data set is decomposed by DWT. The one low-frequency signal and one high-frequency signal are regarded as the available time series and white noise, respectively. It could be expressed as follows:
O H = { oh 1 , oh 2 , , oh p } W A H = { h 1 , h 2 , , h p }
O T R = { ox 1 , ox 2 , , ox p } W A T R = { x 1 , x 2 , , x p } O T E = { ox n + 1 , ox n + 2 , , ox n + m } W A T E = { x n + 1 , x n + 2 , , x n + m }
Secondly, we select the last six electric load data as the input variable and the following one as the output variable of the hybrid model. The formula of the system is as follows:
h i s t o r i c a l d a t a { h 1 h 2 h 3 h 4 h 5 h 6 h 7 h 2 h 3 h 4 h 5 h 6 h 7 h 8 h p 6 h p 5 h p 4 h p 3 h p 2 h p 1 h p t r a i n d a t a { x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x n 6 x n 5 x n 4 x n 3 x n 2 x n 1 x n t e s t d a t a { x n 5 x n 4 x n 3 x n 2 x n 1 x n x n + 1 x n 4 x n 3 x n 2 x n 1 x n x n + 1 x n + 2 x n + m 6 x n + m 5 x n + m 4 x n + m 3 x n + m 2 x n + m 1 x n + m
Thirdly, the distances between the training (or testing) target and one of historical targets is represented by the following formula:
Suppose x i stands for any one of the training (or testing) targets and { x i 6 , x i 5 , x i 4 , x i 3 , x i 2 x i 1 } are the corresponding characteristic indicators of the training (or testing) data. In addition, h j is the historical data target and { h i 6 , h i 5 , h i 4 , h i 3 , h i 2 , h i 1 } are the corresponding characteristic indicators of the historical data. In this paper, the distance is calculated by Euclidean distance:
d i , j = ( x i 6 h j 6 ) 2 + ( x i 5 h j 5 ) 2 + + ( x i 1 h j 1 ) 2
Furthermore, the historical target of x i and the distance corresponding to the historical target can be expressed as follows:
{ ( d i , 7 , h 7 ) , ( d i , 8 , h 8 ) , ( d i , 9 , h 9 ) , , ( d i , p , h p ) }
Next, the distances can be listed in ascending order and the first k historical targets can be obtained as h i ( 1 ) , h i ( 2 ) , , h i ( k ) .
Then, the kernel function can be obtained by ELM:
x i = f ( h i ( 1 ) , h i ( 2 ) , h i ( 3 ) , , h i ( k ) )
where the kernel function f ( ) can be trained by ELM. It cannot be ignored that a simple linear relationship is employed as the traditional method to build the kernel function, but the method cannot deal with the nonlinear relationship. Hence, ELM is selected to establish the kernel function because of its great accuracy and high speed. Finally, the prediction value of the innovative EWKM hybrid model can be obtained, and the basic structure of the proposed model is shown in Figure 2.

3. Empirical Study

3.1. Study Area Description

New South Wales is a state on the east coast of Australia (Figure 3). The estimated population of NSW at the end of March 2016 was 7.7 million, making it Australia’s most populous state. NSW is Australia’s economic powerhouse, and also one of the most active economic regions in the Asian-Pacific region. Among its industrial sectors, the most outstanding are the iron and steel industries. According to the data of the Australian immigration information network, the state’s GDP accounted for one third of Australia’s gross domestic product and more than 35% of the country’s products and services are manufactured in NSW. Its steel production accounted for about 85% of Australia’s total output, centering on the port of Newcastle and Ken Blah. Coal and related products are the state’s biggest exports. With an A$5 billion value, they account for about 19% of all exports. Electricity is inseparable from both industry and people’s life, and the electricity department must therefore ensure adequate power production.

3.2. Data Description

To verify the effectiveness of the proposed model, the data sets of electric load (Mw) from NSW (Australia) are used as the experimental data. They were obtained from the Australian Energy Market Operator (http://www.aemo.com.au/) [36]. Since the proposed hybrid model EWKM contains KNN and ELM, the sample data are divided into three groups: the first group is the historic subset of KNN, which contains 17,568 data points (from 1 January 2016 00:30 to 1 January 2017 00:00); the second is the training subset of ELM, which contains 2832 data points (from 1 January 2017 00:30 to 1 March 2017 00:00); the remaining 2158 data points (from 1 March 2017 00:30 to 14 April 2017 23:30) as the testing subset of ELM can be seen in Table 1 and Figure 4.

3.3. ELM-KNN as a Simulation Tool

The purpose of this section is to examine the fitting capacity of the combination of ELM and KNN for solving complex simulation problems. One thing that must be mentioned is that complex linear systems always appear in the application of projects, which can make it difficult to establish a suitable model, so there is a need for more research on how to find the reasonable and efficient methods to establish the nonlinear relationship between the input and the output. This section proposes a hybrid method that combines ELM and KNN. Three functions with different characteristics are employed as the benchmarks to prove its high fitting capacity. The structure of the three benchmarks considered in the experiment are as follows:
  • F1: Sphere function (d = 2):
    y = x 1 2 + x 2 2
  • F2: Rosenbrock function (d = 2):
    y = 20 + x 1 2 10 cos ( 2 π x 1 ) + x 2 2 10 cos ( 2 π x 2 )
  • F3: Ackley function (d = 2):
    y = 20 e 0.2 1 30 ( x 1 2 + x 2 2 ) e 1 30 [ cos ( 2 π x 1 ) + cos ( 2 π x 1 ) ]
The range of variables x1 and x2 are all from −2 to 2, with the step size of 0.04. That means the variables x1 and x2 have 101 possible values, respectively, and the algorithm can produce 10,201 groups of experimental data. After random ordering of these experimental data, all data are divided into three groups: the first group is the historical subset of KNN, which contains 8201 data points (from No. 1 to No. 8201); the second group is the training subset of ELM, which contains 1500 data points (from No. 8202 to No. 9701); the remaining 500 data points (from No. 9702 to No. 10,201) as the testing subset of ELM.
In order to prove the high fitting capacity of the proposed method ELM-KNN, three other method: KNN, BPNN-KNN and SVM-KNN are employed for comparison. After implementing these algorithms using Matlab 2014(b), we have carried out extensive simulations and each algorithm has been run 100 times so as to perform an average statistical analysis. Table 2 shows the results with the average value of sum error and time are taken as the criteria.
Multiple studies show that ELM-KNN can outperform KNN and other hybrid algorithm for solving non-linear simulation problems. The traditional KNN algorithm is suitable for solving linear problems, so although the run time is short, the sum error of this single algorithm is large. Therefore, this section take other algorithms which are suitable for non-linear fitting to combine with KNN. These hybrid algorithms are ELM-KNN, BPNN-KNN and SVM-KNN. It can be observed from Table 2 that while the accuracy of these hybrid algorithm is similar, the accuracy of ELM-KNN is slightly better than that of BPNN-KNN and SVM-KNN. More importantly, a marked improvement can be seen in the run speed with ELM-KNN. Each function evaluation is virtually instantaneous on a modern PC. For example, the computation time with ELM-KNN on a 2.13 GHz desktop is between 4 and 5 s, which was much superior to BPNN-KNN and SVM-KNN. These results show the high accuracy and efficiency of ELM-KNN that make it a very powerful tool for fitting.

3.4. Evaluation Criteria

The mean absolute error (MAE), the mean relative error (MRE) and the correlation coefficient (R) are used to evaluate the reliability of EWKM model. MAE and MRE measure residual errors, which give a global idea of the difference between the observed and forecasted values. MAE is used to evaluate the absolute error range of the predicted value, while MRE is used to reflect the specification of the predicted value on average. The lower the values of MAE and RMSE, the better the model is. The proportion of the total variance in the observed data can be described by the correlation coefficient (R). R is better when it is close to one.
MAE, MRE and R are calculated as follows:
M A E = i = 1 n | y i y ^ i | / n
M R E = 1 n i = 1 n | y ^ i y i y i |
R = i = 1 n ( y i y ¯ ) ( y ^ i y ^ ¯ ) i = 1 n ( y i y ¯ ) i = 1 n ( y ^ i y ^ ¯ )
where y i is the observed value, y ^ i is the predictive value to y i , y ¯ is the average of the observed value and n is the number of the observations of the validation set.

3.5. The Process of the Proposed Hybrid Model

3.5.1. Wavelet De-noising

The electrical demand is affected by a variety of factors, so electric load time series are usually accompanied by high noise. The direct forecasting of electric load with noisy data usually results in large errors. In this section, the DWT is executed for efficiently removing the noise from the observed data. In general, a normative procedure to select the decomposition level does not exist, so the selection is based on multiple experiments. In this study, the Daubechies wavelet of order 3 (db3) has a better performance, so the db3 is adopted in the wavelet denoising process [29,37]. Considering the characteristics of the experimental data, after testing the effect of different levels with db3, level 1 works best. The decomposition figure includes the approximation coefficients at level 1 (a1) and the detail coefficients at level 1 (d1). The decomposition process of the experimental data is shown in Figure 5. The a1 is a smoothed version of the original series and it represents the low frequency signal. Thus, a1 is selected to build forecasting model.

3.5.2. The Process of ELM-KNN

In Section 2.4, the EWKM hybrid model is established. The parameters of the hybrid model k = 12 are selected by lots of experiment tests. According to the choice of the electric load time series, there are 2832 training data from 1 January 2017 00:30 to 1 March 2017 00:00 (n = 2832) and 2158 testing data from 1 March 2017 00:30 to 14 April 2017 23:30 (m = 2158). Meanwhile, the half-hour electric load of the year 2016 (p = 17568) is selected as the historical data set in the experiment. Finally, the prediction of electric load could be calculated

3.6. Results and Analysis

3.6.1. Results of the Proposed Model

The results obtained from the modified hybrid model EWKM agree well with the original values. As is shown in Figure 6, the forecasting curve of EWKM closely approaches the original one. This figure confirms that the EWKM model has great performance in predicting the electric load series.
As Table 3 shows the MRE, R and MAE of the proposed model EWKM are 0.0262, 0.9660 and 196.7408, respectively.
The numerical value of the mean relative error (MRE) is close to zero and the correlation coefficient (R) is close to one. It clearly confirms that the EWKM model can capture the non-stationary and highly noisy features of the electric load series, and the new model can effectively improve the forecasting accuracy.

3.6.2. Model Comparisons

This section provides a comparison between the proposed EWKM model and three other benchmark models: EKM, WKM, and WNNM. Note that random selection of the parameters (ELM and BPNN) may cause different results. To prevent this uncertain phenomenon, fifty runs for each method are applied and the results of every run are recorded in Table 4.
The estimation performance of EWKM, EKM and WNNM are assessed by statistical indicators of MRE, MAE and R. The values are presented in Table 5. What needs special explanation is the fact that WKM is a non-parameter method and the result of this model is a fixed value, so it makes no sense to do lots of experiments with WKM. The mark “/” in Table 5 stands for the non-existent statistical values. To verify the EWKM model, this experiment compares it with EKM, WKM and WNNM, using the same electric load data. It must be noticed this analysis is done according to the mean value of three indicators. As listed in Table 5, EWKM has the lowest MRE, MAE and the highest R among the four models. Comparing EKM with EWKM, after introducing wavelet denoising into the model, the MRE and MAE have obviously decreased. Unlike the EKM model which is directly constructed from the original data, we decompose the original data into a low frequency and high frequencies parts, and the results indicate that the EWKM can capture the highly noisy and non-stationary features of electric load data after the wavelet denoising, so it can be concluded that the denoising process is crucial to predict electric loads. Then, comparing WKM with EWKM, the three average indicators MRE, R, MAE of the WKM are 0.0290, 0.9607 and 219.6566, respectively, which are all inferior to those of EWKM. This indicates that optimizing the kernel function by using the ELM has a wonderful effect, so the ELM is necessary to predict the electric load series. As for the WNNM, the BPNN is an artificial neural network which is widely used to predict time series. When using the BPNN to forecast electric loads, the two statistical errors MAE and MRE are clearly larger than the errors of the EWKM model. The comparison of the results between the four models shows that the proposed model EWKM can effectively improve forecasting accuracy. Therefore, the conclusion can be drawn that every part of the proposed model is suitable and reasonable. The proposed model EWKM that includes denoising processing, KNN and ELM constitutes a significant improvement in electric load forecasting.

4. Conclusions and Future Work

In electricity demand forecasting, noise signals, caused by various unstable factors, often corrupt electric load series. The contribution of this paper is a method that uses a hybrid model based on wavelet denoising processing. In previous studies, models had been usually established with the original data, however, this paper takes the low-frequency signal to modeling, so that it can reduce errors caused by noise signals. Moreover, in the construction of the kernel function, the traditional way is to construct a linear relationship, but this paper introduces the ELM into the establishment of the kernel function, so that it can optimize the KNN algorithm. Through the analysis of experimental results, a conclusion can be drawn that the every part included in the new hybrid model is necessary to predict future electric loads.
However, this paper only takes electric load data as the research subject, without taking other related variables into consideration. To resolve such limitations, future research should aim to include other factors which may influence the electric demand, such as the population and GDP, and there’s a lot of room for improvement.

Acknowledgments

The authors are grateful the editor and anonymous reviewers for their suggestions in proving the quality of the paper. This research is supported by the National Natural Sciences Foundation of China (No. 41571016).

Author Contributions

Weide Li designed research; Demeng Kong wrote the paper; Jinran Wu performed research and analyzed data.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANNartificial neural network
BAMObinary animal migration optimization
BPNNback propagation neural network
CScuckoo search algorithms
DWTdiscrete wavelet transform
EEMDensemble empirical mode decomposition
EEuNN frameworkevolving fuzzy neural network framework
ELMextreme learning machine
EMDempirical mode decomposition
FTSfuzzy time series
GHSAglobal harmony search algorithm
GMgrey model
GPRMgrey prediction with rolling mechanism
KNN regressionk-nearest neighbor regression
LSSVMleast squares support vector machines
MAEmean absolute error
MREmean relative error
PSOparticle swarm optimization
QPSOquantum particle swarm optimization algorithm
Rcorrelation coefficient
SARIMAseasonal auto-regressive integrated moving average
SLFNsingle hidden-layer feed-forward network
SSAsingular spectrum analysis
SVMsupporter vector machine
SVRsupport vector regression
WTwavelet transform

References

  1. Goia, A.; May, C.; Fusai, G. Functional clustering and linear regression for peak load forecasting. Int. J. Forecast. 2010, 26, 700–711. [Google Scholar] [CrossRef]
  2. Zhou, P.; Ang, B.W.; Poh, K.L. A trigonometric grey prediction approach to forecasting electricity demand. Energy 2006, 31, 2839–2847. [Google Scholar] [CrossRef]
  3. Akay, D.; Atak, M. Grey prediction with rolling mechanism for electricity demand forecasting of Turkey. Energy 2007, 32, 1670–1675. [Google Scholar] [CrossRef]
  4. Chang, P.C.; Fan, C.Y.; Lin, J.J. Monthly electricity demand forecasting based on a weighted evolving fuzzy neural network approach. Int. J. Electr. Power Energy Syst. 2011, 33, 17–27. [Google Scholar] [CrossRef]
  5. Kavaklioglu, K. Modeling and prediction of Turkey’s electricity consumption using Support Vector Regression. Appl. Energy 2011, 88, 368–375. [Google Scholar] [CrossRef]
  6. Wang, J.; Zhu, W.; Zhang, W.; Sun, D. A trend fixed on firstly and seasonal adjustment model combined with the ε-SVR for short-term forecasting of electricity demand. Energy Policy 2009, 37, 4901–4909. [Google Scholar] [CrossRef]
  7. Kucukali, S.; Baris, K. Turkey’s short-term gross annual electricity demand forecast by fuzzy logic approach. Energy Policy 2010, 38, 2438–2445. [Google Scholar] [CrossRef]
  8. Dash, P.K.; Liew, A.C.; Rahman, S.; Ramakrishna, G. Building a fuzzy expert system for electric load forecasting using a hybrid neural network. Expert Syst. Appl. 1995, 9, 407–421. [Google Scholar] [CrossRef]
  9. Taylor, J.W. Triple seasonal methods for short-term electricity demand forecasting. Eur. J. Oper. Res. 2010, 204, 139–152. [Google Scholar] [CrossRef]
  10. Moghram, I.; Rahman, S. Analysis and evaluation of five short-term load forecasting techniques. IEEE Trans. Power Syst. 1989, 4, 1484–1491. [Google Scholar] [CrossRef]
  11. Wu, H.; Niu, D.; Song, Z. Short-term electric load forecasting based on data mining. Int. J. Eng. Technol. 2017, 8, 250–253. [Google Scholar] [CrossRef]
  12. Chen, H.Y.; Hong, W.-C.; Shen, W.; Huang, N.N. Electric load forecasting based on a least squares support vector machine with fuzzy time series and global harmony search algorithm. Energies 2016, 9, 70. [Google Scholar] [CrossRef]
  13. Huang, M.L. Hybridization of chaotic quantum particle swarm optimization with SVR in electric demand forecasting. Energies 2016, 9, 426. [Google Scholar] [CrossRef]
  14. Zhang, X.; Wang, J.; Zhang, K. Short-term electric load forecasting based on singular spectrum analysis and support vector machine optimized by Cuckoo search algorithm. Electr. Power Syst. Res. 2017, 146, 270–285. [Google Scholar] [CrossRef]
  15. Karlsson, M.S. Nearest Neighbor Regression Estimators in Rainfall-Runoff Forecasting. Ph.D. Thesis, University Of Arizona, Tucson, AZ, USA, 1985. [Google Scholar]
  16. Poloczek, J.; Andre Treiber, N.; Kramer, O. KNN regression as geo-impuation method for spatio-temporal wind data. In Proceedings of the International Joint Conference SOCO’14-CISIS’14-ICEUTE’14, Bilbao, Spain, 25–27 June 2014; Volume 299. [Google Scholar] [CrossRef]
  17. Ban, T.; Zhang, R.; Pang, S.; Sarrafzadeh, A.; Inoue, D. Referential KNN regression for financial time series forecasting. In Proceedings of the International Conference on Neural Information Processing, Daegu, Korea, 3–7 November 2013; Volume 8226, pp. 601–608. [Google Scholar]
  18. Hu, J.; Liu, J.; Liu, Y.; Gao, C. Emd-KNN model for annual average rainfall forecasting. J. Hydrol. Eng. 2013, 18, 1450–1457. [Google Scholar] [CrossRef]
  19. Li, X.; Ma, S.; Wang, Y. BamoKNN: A novel computational method for predicting the apoptpsis protein locations. In Proceedings of the IEEE International Conference on Bioinformatic & Biomedicine, Shenzhen, China, 15–18 December 2016; pp. 743–746. [Google Scholar]
  20. Wang, Y.; Nan, Z.; Chen, H.; Wu, X. Correction of daily precipitation data of ITPCAS dataset over the Qinghai-Tibetan plateau with KNN model. In Proceedings of the IGARSS IEEE International Geoscience & Remote Sensing Symposium, Beijing, China, 10–15 July 2016; pp. 593–596. [Google Scholar]
  21. Huang, G.B.; Rong, S.; You, K. Extreme learning machine:a new learning scheme of feedforward neural networks. In Proceedings of the International Joint Conference on Neural Networks, Budapest, Hungary, 25–29 July 2004; Volume 2, pp. 985–990. [Google Scholar]
  22. Ma, Z.; Luo, G.; Huang, D. Short term traffic flow prediction based on on-line sequential extreme learning machine. In Proceedings of the Eighth International Conference on Advanced Computational Intelligence, Chiang Mai, Thailand, 14–16 February 2016; pp. 143–149. [Google Scholar]
  23. Masri, D.; Woon, W.L.; Aung, Z. Soil property prediction: An extreme learning machine approach. In Proceedings of the International Conference on Neural Information Processing: Neural Information Processing, Istanbul, Turkey, 9–12 November 2015; pp. 18–27. [Google Scholar]
  24. Liao, H.Y.; Wang, H.; Hu, Z.J.; Wang, K.; Li, Y.; Huang, D.S.; Ning, W.H.; Zhang, C.X. Stock price forecasting based on extreme learning machine. Comput. Mod. 2014, 12. (In Chinese) [Google Scholar]
  25. Zhang, Y.H.; Wang, H.; Hu, Z.J.; Wang, K.; Li, Y. A hybrid short-term wind speed forecasting model based on wavelet decomposition and extreme learning machine. Adv. Mater. Res. 2013, 860–863, 361–367. [Google Scholar]
  26. Shamshirband, S.; Mohammadi, K.; Yee, D.L.; Petković, D.; Mostafaeipour, A. A comparative evaluation for identifying the suitability of extreme learning machine to predict horizontal global solar radiation. Renew. Sustain. Energy Rev. 2015, 52, 1031–1042. [Google Scholar] [CrossRef]
  27. Dabuechies, I. The wavelet transform, time-frequency loaclization and signal analysis. IEEE Trans. Inf. Theory 1990, 36, 6–7. [Google Scholar]
  28. Wang, L.P.; Shekhar, G. Neural networks and wavelet denoising for stock trading and prediction. Time Ser. Anal. Model. Appl. 2013, 47, 229–247. [Google Scholar]
  29. Wang, J.; Wang, J.; Li, Y.; Zhu, S.; Zhao, J. Techniques of applying wavelet denoising into a combined model for short-term load forecasting. Int. J. Electr. Power Energy Syst. 2014, 62, 816–824. [Google Scholar] [CrossRef]
  30. Qin, Y.; Huang, S.; Zhao, Q. Prediction model of chaotic time series based on wavelet denoising and LS-SVM and ITS application. J. Geod. Geodyn. 2008, 28, 96–100. [Google Scholar]
  31. Abbaszadeh, P. Improving hydrological process modeling using optimized threshold-based wavelet denoising technique. Water Resour. Manag. 2016, 30, 1701–1721. [Google Scholar] [CrossRef]
  32. Kalteh, A.M. Improving forecasting accuracy of streamflow time series using least squares support vector machine coupled with data-preprocessing techniques. Water Resour. Manag. 2016, 30, 747–766. [Google Scholar] [CrossRef]
  33. Liang, N.Y.; Huang, G.B.; Saratchandran, P.; Sundararajan, N. A fast and accurate on-line sequential learning algorithm for feedforward networks. IEEE Trans. Neural Netw. 2006, 17, 1411–1423. [Google Scholar] [CrossRef] [PubMed]
  34. Cover, T.M.; Hart, P.E. Nearest neighbor pattern classification. IEEE Trans. Inf. Theory 1967, 13, 21–27. [Google Scholar] [CrossRef]
  35. Lall, U.; Sharma, A. A nearest neighbor bootstrap for resampling hydrologic time series. Water Resour. Res. 1996, 32, 679–694. [Google Scholar] [CrossRef]
  36. Australian Energy Market Operator (AEMO). Available online: http://www.aemo.com.au (accessed on 12 May 2017).
  37. Ali, Y.M.K. Forecasting Power and Wind Speed Using Artificial and Wavelet Neural Networks for Prince Edward Island (PEI); Electrical and Computer Engineering at Dalhouse University: Halifax, NS, Canada, 2010. [Google Scholar]
Figure 1. The structure of ELM.
Figure 1. The structure of ELM.
Energies 10 00694 g001
Figure 2. The basic structure of the proposed model.
Figure 2. The basic structure of the proposed model.
Energies 10 00694 g002
Figure 3. The study area location.
Figure 3. The study area location.
Energies 10 00694 g003
Figure 4. The electric load series of New Souths Wales, (a) historical data set; (b) train data set; (c) test data set.
Figure 4. The electric load series of New Souths Wales, (a) historical data set; (b) train data set; (c) test data set.
Energies 10 00694 g004
Figure 5. The process of wavelet denoising, (a) the decomposition of historical data set; (b) the decomposition of train and test data sets; note: a1, approximation signal; d1, detail signal.
Figure 5. The process of wavelet denoising, (a) the decomposition of historical data set; (b) the decomposition of train and test data sets; note: a1, approximation signal; d1, detail signal.
Energies 10 00694 g005
Figure 6. Plot of the original data and predictive values.
Figure 6. Plot of the original data and predictive values.
Energies 10 00694 g006
Table 1. Statistical parameters in each data set.
Table 1. Statistical parameters in each data set.
Data SetNMinMaxMeanStdSK
StatisticsStdStatisticsStd
Historical data set17,56817,568535013,459797812570.4730.0180.185
Train data of ELM28322832570413,986861217480.7850.0460.234
Test data of ELM21582158548911,000776810440.0000.053−0.557
Min is the minimum; Max is the maximum; Std is the standard deviation; S is the skewness; K is the kurtosis.
Table 2. The Fitting performance of the four methods.
Table 2. The Fitting performance of the four methods.
FunctionCriteriaKNNELM-KNNBPNN-KNNSVM-KNN
F1sum_error26.64617.94118.32917.806
time1.18754.78378.07386.8382
F2sum_error429.717378.105380.040378.502
time1.26154.67388.12566.9026
F3sum_error19.5114.6854.9364.824
time1.17014.81508.13736.8376
Table 3. Indicators of the EWKM.
Table 3. Indicators of the EWKM.
ModelMRERMAE
EWKM0.02620.9660196.7408
Table 4. The performance of the EWKM model compared to EKM and WNNM.
Table 4. The performance of the EWKM model compared to EKM and WNNM.
Run No.EWKMEKMWNNM
MREMAERMREMAERMREMAER
10.0272203.93540.96430.0310234.71190.95590.0340258.22240.9484
20.0254191.52160.96770.0330247.79880.95130.0361271.34610.9431
30.0262197.22020.96580.0306232.15430.95610.0336255.23620.9487
40.0256193.30320.96680.0331248.13380.94950.0361271.02280.9412
50.0269200.09460.96520.0298226.24270.95780.0328249.49730.9509
60.0255190.67210.96700.0324244.74480.95190.0354268.03750.9439
70.0267199.55740.96540.0308234.32780.95680.0342260.60110.9491
80.0278209.63880.96220.0306232.24620.95750.0337255.60750.9505
90.0256191.86920.96690.0305231.06650.95630.0336254.90350.9487
100.0258193.79430.96620.0312237.25510.95400.0344261.80850.9459
110.0264198.31340.96450.0328245.29840.95220.0355266.47260.9442
120.0265197.85500.96540.0308233.81560.95570.0338256.79170.9483
130.0261195.97840.96630.0311236.24160.95450.0344261.12750.9463
140.0265198.51580.96560.0302229.52130.95750.0334253.48240.9501
150.0253190.11620.96760.0336255.49140.95000.0362275.45770.9433
160.0254190.73570.96740.0323243.81350.95230.0349263.80250.9451
170.0274205.23060.96330.0309232.53880.95570.0342257.98390.9472
180.0271200.27890.96570.0315237.71040.95430.0346261.75180.9460
190.0250187.60620.96840.0330246.91600.95120.0356267.74340.9431
200.0257192.49590.96710.0310234.32470.95560.0342258.76660.9475
210.0261196.60630.96660.0334251.68970.94930.0365275.87470.9407
220.0256193.24920.96680.0337251.44580.94910.0367274.84530.9409
230.0267200.59480.96460.0315237.08840.95530.0345260.62530.9478
240.0258195.81410.96640.0298226.94380.95810.0329250.60530.9508
250.0266199.72600.96540.0302230.20090.95720.0334254.97890.9497
260.0253191.29830.96760.0304230.22000.95700.0336255.19380.9493
270.0271204.07360.96410.0306232.71970.95640.0336256.56390.9490
280.0265199.21130.96540.0336253.11850.94940.0365276.14220.9412
290.0278207.29880.96390.0308232.18030.95590.0337254.50380.9483
300.0265196.99440.96660.0312236.75810.95450.0344261.53990.9468
310.0250186.95440.96810.0319241.28990.95240.0351266.22530.9437
320.0280208.44720.96300.0311234.45060.95480.0343258.73300.9469
330.0254191.09990.96740.0330249.98620.95090.0361273.80560.9431
340.0255191.70660.96700.0303229.54050.95700.0334252.92840.9495
350.0259194.02020.96650.0305231.51810.95650.0332252.73510.9498
360.0262196.79850.96610.0341256.52690.94900.0370279.29290.9411
370.0258193.27830.96630.0317241.17370.95300.0348265.55010.9448
380.0265199.25080.96460.0324244.22650.95150.0352266.19420.9442
390.0270201.38690.96510.0306232.20340.95630.0336255.27850.9493
400.0265200.27360.96520.0321243.90470.95260.0352267.38540.9448
410.0278207.91200.96310.0321241.45780.95260.0350264.28430.9442
420.0251188.31510.96830.0321242.78730.95360.0349264.51000.9462
430.0260194.47050.96620.0313238.45540.95460.0345263.16480.9468
440.0257192.40900.96780.0304231.63170.95750.0336256.31010.9503
450.0249187.36780.96860.0312237.91460.95380.0345262.73470.9457
460.0273205.13110.96380.0320243.18070.95300.0352267.25740.9451
470.0250187.81900.96820.0306231.53290.95680.0338256.10220.9485
480.0290214.22470.96220.0328247.64850.94960.0359272.00240.9411
490.0252189.68220.96740.0315238.62180.95270.0345261.94170.9450
500.0256192.89320.96660.0313237.39660.95470.0347264.26950.9465
Table 5. Comparing the three criteria of four models.
Table 5. Comparing the three criteria of four models.
IndicatorsModelsMeanMinimumMaximumStandard DeviationMedianUpper QuantileLower Quantile
MREEWKM0.02620.02490.02900.00090.02610.02550.0267
EKM0.03160.02980.03410.00110.03130.03060.0324
WNNM0.03460.03280.03700.00110.03450.03370.0352
WKM0.0290//////
MAEEWKM196.7408186.9544214.22476.4440196.2924191.7473200.2289
EKM238.8433226.2427256.52697.7976237.3259232.2141244.1461
WNNM262.4248249.4973279.29297.4146261.7802256.1542267.0612
WKM219.6566//////
REWKM0.96600.96220.96860.00160.96630.96510.9671
EKM0.95400.94900.95810.00270.95450.95220.9563
WNNM0.94630.94070.95090.00300.94640.94420.9487
WKM0.9607//////

Share and Cite

MDPI and ACS Style

Li, W.; Kong, D.; Wu, J. A Novel Hybrid Model Based on Extreme Learning Machine, k-Nearest Neighbor Regression and Wavelet Denoising Applied to Short-Term Electric Load Forecasting. Energies 2017, 10, 694. https://doi.org/10.3390/en10050694

AMA Style

Li W, Kong D, Wu J. A Novel Hybrid Model Based on Extreme Learning Machine, k-Nearest Neighbor Regression and Wavelet Denoising Applied to Short-Term Electric Load Forecasting. Energies. 2017; 10(5):694. https://doi.org/10.3390/en10050694

Chicago/Turabian Style

Li, Weide, Demeng Kong, and Jinran Wu. 2017. "A Novel Hybrid Model Based on Extreme Learning Machine, k-Nearest Neighbor Regression and Wavelet Denoising Applied to Short-Term Electric Load Forecasting" Energies 10, no. 5: 694. https://doi.org/10.3390/en10050694

APA Style

Li, W., Kong, D., & Wu, J. (2017). A Novel Hybrid Model Based on Extreme Learning Machine, k-Nearest Neighbor Regression and Wavelet Denoising Applied to Short-Term Electric Load Forecasting. Energies, 10(5), 694. https://doi.org/10.3390/en10050694

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