Journal of Electrical and Computer Engineering

Volume 2016 (2016), Article ID 2165324, 10 pages

http://dx.doi.org/10.1155/2016/2165324

## A Novel Hybrid Method for Short-Term Power Load Forecasting

Department of Economic Management, North China Electric Power University, Baoding, China

Received 20 March 2016; Revised 21 June 2016; Accepted 11 July 2016

Academic Editor: Jit S. Mandeep

Copyright © 2016 Huang Yuansheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Influenced by many uncertain and random factors, nonstationary, nonlinearity, and time-variety appear in power load series, which is difficult to forecast accurately. Aiming at locating these issues of power load forecasting, an innovative hybrid method is proposed to forecast power load in this paper. Firstly, ensemble empirical mode decomposition (EEMD) is used to decompose the power load series into a series of independent intrinsic mode functions (IMFs) and a residual term. Secondly, genetic algorithm (GA) is then applied to determine the best weights of each IMF and the residual term named ensemble empirical mode decomposition based on weight (WEEMD). Thirdly, least square support vector machine (LSSVM) and nonparametric generalized autoregressive conditional heteroscedasticity (NPGARCH) are employed to forecast the subseries, respectively, based on the characteristics of power load series. Finally, the forecasted power load of each component is summed as the final forecasted result of power load. Compared with other methods, the forecasting results of this proposed model applied to the electricity market of Pennsylvania-New Jersey-Maryland (PJM) indicate that the proposed model outperforms other models.

#### 1. Introduction

In the operation of power system, power load forecasting is not only an important part of power system planning but also one of the most influential factors for the improvement of social economy developing, which has a significant impact on generation, transmission, and distribution. Therefore, it is crucial to design an efficient method to forecasting power load accurately by improving both economic efficiency and power supply quality and enhancing reliability of power system operation [1, 2].

However, the power load forecasting results are influenced by many uncertain random factors, including changes of nature environment and electricity price and factors influencing power load [3, 4], causing the complex volatility characteristics of power load, which lead the time-series of power load series to inherent nonstationary and discrete in practice, even in the short-term dynamics. Notably, due to the factors, influencing the short-term power load, previous studies [5–7] mainly focus on the establishment of scientific forecasting method to enhance the adaptability for these multifactors. Therefore, the time-varying components of power load are hardly separated to extract the information in them.

Under this circumstance, a novel hybrid method for power load forecasting has been proposed in this paper. Based on artificial algorithm, the advantages of econometric models and signal processing are emphasized by this research. Specifically, this hybrid method for load power forecasting includes ensemble empirical mode decomposition based on variable weights (WEEMD), genetic algorithm (GA), least square support vector machine (LSSVM), and nonparametric generalized autoregressive conditional heteroscedasticity (NPGARCH). Notably, in order to test the validity and feasibility of the model, amount of historical data of power load in American Electric Power (AEP) has been adopted to apply this new hybrid method and compare it with previously well-known methods for the power load forecasting accuracy. Results indicate that this hybrid method outperforms the compared methods with the forecasting accuracy.

The rest of the paper is organized as follows: Section 2 summaries previous literature over methodology for power load forecasting; Section 3 briefly describes EEMD, GA, nonparametric GARCH, and LSSVM and proposes the improved WEEMD, as well as presenting the procedure of the newly hybrid method; Section 4 evaluates the forecasting accuracy and tests the robustness of the proposed method, by comparing with other forecasting methods. Section 5 draws conclusions.

#### 2. Related Literature Review

Current power load forecasting methods can be broadly divided into two categories: statistical analysis method [8–11] and artificial algorithm [12–15]. The traditional statistical analysis methods are always under the framework of time-series analysis [16]. By using ARMA and ARMA with weather inputs, Bolzern and Fronza [8] propose two short-term predictors of the winter power load. Respectively, the power load forecast results are satisfied with both predictors and similar in these two cases. In order to model certain behavior of energy consumption, Wang [9] employs conventional fuzzy systems based on the integrated algorithm, and results show that the monthly electricity consumption of Iran is accurately forecasted with this proposed method. In [10], considering the nonlinearity and volatility of power load series, analyzing the threshold characteristics, Wang et al. adopt a novel double-threshold generalized autoregressive conditional heteroscedasticity (DT-GARCH) model for short-term power load forecasting. To forecast the power load efficiently and scientifically, Liu [11] synthesizes the time-series modeling with the regression modeling. According to the residual error from regression forecasting, the accuracy of power load forecasting can be improved by correcting the error of GARCH model.

As the artificial algorithm applied in the power load forecasting is prevalent in recent decades, there are plenty of researches related to that. Taking the characteristics of randomness, tendency, and periodicity of short-term power load into account, Sun and Ye [12] propose a model based on LSSVM and fruit fly algorithm (FOA) for short-term load forecasting. In order to reduce the nonlinear power load sequence and improve the accuracy of forecasting result, Hu and Chang [13] decompose the time-series into several components by local wave method and use optimal parameters to establish the DEMD-LSSVM model for component forecasting. Lauret et al. [14] utilize Bayesian techniques to design an optimal neural network (NN) based model for power load forecasting, where Bayesian technique offers great advantages on traditional neural network learning methods. Based on an artificial neural network (ANN), Shao et al. [15] proposed an approach combined with a fuzzy system for short-term power load forecasting. The load error is obtained from the historical information and past forecasted load errors are caused by fuzzy systems, and the final forecasted load can be obtained by adding the load error to preliminary load forecasted by ANN. In [16], to improve the forecasting performance by searching a suitable parameters combination, the paper presents an SVR-based electric load forecasting model by applying a novel algorithm named chaotic ant swarm optimization (CAS). Combined with the chaotic behavior of single ant and self-organization behavior of ant colony in the foraging process, the CAS is proposed to overcome premature local optimum. In [17], a combination of the wavelet transform (WT) and gray model is proposed for short-term power load forecasting, which is improved by particle swarm optimization (PSO). With this proposed method taking mean temperature, mean relative humidity, mean wind speed, and previous days load data into consideration and eliminating the high frequency of historical data built by the WT, the accuracy can be largely improved.

With the improvement of the forecasting accuracy, Ghelardoni et al. [18] decompose time power load series with empirical mode decomposition (EMD) into two sets of components, respectively, describing the trend of energy consumption values. LSSVM is built to forecast these two components. The hybrid method proposed in [19] enhances the capability for forecasting. This proposed method is modeled based on supporting vector regression (SVR), EMD, and regression (AR), which can simultaneously provide forecasting with high accuracy and interpretability. In order to solve the problem of mode mixing and high frequency random components, Liu et al. [19] proposed an optimized method based on ensemble empirical mode decomposition (EEMD) and subsection particle swarm optimization (SS-PSO). By extracting and reconstructing intrinsic mode function (IMF), the power load series movement is well forecasted. In [20], the proposed method based on complementary ensemble empirical mode decomposition (CEEMD) fuzzy entropy and echo state network (ESN) with leaky integrator neurons (LiESN) enhances the forecasting accuracy of power load.

In summary, although the characteristics of power load series are complex, such as being rich multidimensional, nonlinear, time-varying, and nonstationary, the previous literature has accumulated a great deal of experience about power load forecasting. Particularly, the forecasting methods have been constantly proposed and the forecasting performance has been continuously improved, all of which provide a significant foundation for the present study. However, we can still find that the method of power load forecasting can be problematical based on the analysis and research of previous literature, while there are also some uncertain factors of power load forecasting. For example, both EEMD and EMD decompose the original power load series into several components, and other forecasting methods are directly built based on the decomposed components in most of the literature, without analysis of the weights of each component in original power load series. Besides, as the power load is affected by great uncertainties, the power load series movement is difficult to capture, and there will be nonlinear and time-varying feature in power load series caused by no matter the external factors or the internal factors. But the existing methods for power load forecasting are usually not effective to separate nonlinear and time-varying components of power load and unable to extract their inherent moving mechanisms, which consequently affects the forecasting accuracy.

Based on the potential faultiness mentioned above, a hybrid method for power load forecasting is developed and applied to PJM by this study. Specially, the ensemble empirical mode decomposition (EEMD) model decomposes power load into series of intrinsic mode functions (IMFs) and residual. And genetic algorithm (GA) is used to determine the weights of each component, which is based on deviation between fitted values and actual values. For weighted components with the heteroscedasticity character, GARCH is more favorable for forecasting. However, the LSSVM method is used in the forecasting of the subseries with the characteristics of nonstationary and nonlinearity subseries. On the one hand, power load forecasting can be more convenient and accurate, after the decomposition of the power load components because of the separation of other influential factors. On the other hand, with other neural network models, which cannot fundamentally solve problems of the local minimum, the difficulty in determination of hidden layer, and the slow training rate, the LSSVM can not only get over these disadvantages, but also can improve the accuracy of forecasting. Therefore, based on the components weighted by GA method, LSSVM method can be more in line with the power load study in this research.

#### 3. Power Load Forecasting Methodology

##### 3.1. Genetic Algorithms (GA)

Genetic algorithms code the candidate solutions of an optimization algorithm as a string of characters which are usually binary digits [21]. In accordance with the terminology that is borrowed from the field of genetics, this bit string is usually named as chromosome. The solution represented by its chromosome is considered as an individual. The algorithm starts with the initial generation of the population. The fitness of the individuals within the population is assessed, and new individuals are generated for the further generation. A number of genetic operators, containing selector operator, crossover, and mutation, are available for this purpose. By a number of fixed generations which is the termination condition, the best individual will be obtained with the max fitness, which is the global optimal solution of the issues.

##### 3.2. The WEEMD Method

EMD, an effective method for signal processing, is gradually replaced by EEMD which overcomes the mixing model problem [22–25]. The essence of EEMD is to decompose a time-series into a set of independent intrinsic mode functions (IMFs) and the residue obtained by adding a random Gaussian white noise sequence, which is different from EMD. While the time-series is decomposed by EEMD, the IMF is a function, satisfying the following two conditions:(1)In the whole data set, the number of extreme and the number of zero crossings must either differ or differ at most by one.(2)At any point, mean values of the envelope, defined by both the local maxima and minima, are zero.

For an arbitrary time-series, , procedures of EEMD method can be described as follows:(1)Add the white noise to power load series , with , and set the number of ensemble (): where denotes the th added white noises series and represents the noise-added power load of the th trial.(2)Decompose the noise-added series into IMFs by using EMD, where is the th IMF of the th trial and is the number of IMFs.(3)Repeat Steps and until .(4)Calculate the ensemble mean of trails for each IMF, then , where is the th IMF component by using EEMD.

In this paper, weights assignment, for each IMF, is proposed based on EEMD, and the rationality of this weights assignment proved by this study is as follows: for an arbitrary time-series, , setting any value of each IMF satisfies the function , and by the conditions mentioned above, the number of extreme and zero crossings of each IMF is , which are obtained by the following equations:Assigning weight to each IMF, thenwhere , a constant, is the weights of th IMF.

Obviously, the weights have no effect on the values, which satisfies the first conditions. Specially, mean values of the envelope will not be changed by weight assignment, and the deduction can be expressed as follows:Assigning weight to it,where is the mean value of the envelope of the local maxima and is the mean values of the envelope of the local minima.

GA method is built to determine the weights of each IMF, and individuals are first randomly generated as initial population . By using the genetic operators, selector operator, crossover, and mutation, the best individual can be obtained. Specially, in the selector operator, the fitness function is defined as follows:where is error and is the size of initial population.

As the GA operator is designed to maximize the fitness function, the above minimization problem can be solved by using the following transformation:where is the fitness of th individual.

The election probability of each individual is :

As described above, the process of EEMD is indeed like sifting, which has an effect to eliminate riding waves. The IMFs are extracted from the power load series and contain “information” about the time-series. This paper uses weights as the contribution of IMF to ; the greater the weight is, the larger the amount of “information” IMF contains.

##### 3.3. The Nonparametric GARCH (1, 1) Method

Nonparametric GARCH (1, 1) model for the error fluctuation [26], which requires less assumptions, is defined aswhere is estimation of random processes; is specific transformation sequence of zero mean and unit variance. is the conditional variance. is the volatility.

Nonparametric model can be used to estimate the conditional variance, and (9) can be rewritten as follows:where is the martingale difference sequence.

According to the equations above, the function is regressed between lagged variables and can be estimated by using nonparametric smoothing method and the autoregression function. The calculation of nonparametric GARCH [27–30] is as follows:(1)Using parameters GARCH (1, 1) model fits the volatility , with estimation using maximum likelihood estimate parameter being employed as parameter estimation, where .(2) is as weight, and is estimated by using and . Smooth nonparametric estimation method is applied to obtain the estimation of autoregression function .(3)Standard deviation is obtained by using(4)Increment and return to step if , where is a prespecified maximum number of iterations.

The nonparametric estimation is an improvement over the parametric GARCH estimation of volatility. By means of step continuous iteration, there is little to pick and choose between volatility estimates for various values of . However, the algorithm can often be improved by averaging over the final estimates to obtain

##### 3.4. The LSSVM Method

LSSVM is proposed as an improved algorithm based on support vector machine (SVM) [31–35], with the given training data set , with the input and the output . The following regression model is constructed by using nonlinear mapping function :

With the given training data set , the optimization problem of LSSVM is defined as follows:

According to the Kuhn-Tucker conditions, the LSSVM regression model can be expressed as

is the kernel function, which can map variables to the feature space and avoid high dimensional complex difficulties. This paper applies RBF as the kernel function, which is defined aswhere is the kernel function parameter. The LSSVM method can be used by establishing the parameters and .

##### 3.5. The Hybrid Method for Power Load Forecasting

Considering the complex volatility characteristics of power load series, much more scientific forecasting models are required to address the nonlinearity and time variations. Under this circumstance, the WEEMD, based on EEMD, is proposed to extract different components of power load series and assign weights to each component (IMF) according to its standard deviation, where the LSSVM is presented to forecast the subseries with the characteristics of nonstationary and nonlinearity, and the nonparametric GARCH (1, 1) is used to forecast the subseries with heteroscedasticity. With this hybrid model, the power load movement can be well forecasted. The procedures of the improved model can be described as Figure 1 and the concrete steps are given as follows:(1)The power load series is first decomposed by EEMD into intrinsic mode functions (IMFs) and one residual series, and then where is the original power load series and and are decomposed from the series.(2)Each IMF series and the residual series are assigned to be weighted by GA, which can be represented as follows:(3) To verify the existence of ARCH effect and heteroscedasticity, ARCH-LM is used to test the subseries which is related to the stochastic error, and then where is the subseries with heteroscedasticity and represents the subseries without heteroscedasticity.(4)The LSSVM is built to forecast the future values of ; meanwhile, the nonparametric GARCH model is presented to forecast the future values of , and their forecasted results are and , respectively, which can be represented as follows: where is the forecasted value of power load series.