#### Abstract

Annual electricity consumption forecasting is one of the important foundations of power system planning. Considering that the long-term electricity consumption curves of developing countries usually present approximately exponential growth trends and linear and accelerated growth rate trends may also appear in certain periods, this paper first proposes a small-sample adaptive hybrid model (AHM) to extrapolate the above curves. The iterative trend extrapolation equation of the proposed model can simulate the linear, exponential, and steep trends adaptively at the same time. To estimate the equation parameters using small samples, the partial least squares (PLS) and iteration starting point optimization algorithms are suggested. To evaluate forecasting performance, the artificial neural network (ANN), grey model (GM), and AHM are used to forecast electricity consumption in China from 1991 to 2014, and then the results of these models are compared. Analysis of the forecasting results shows that the AHM can overcome stochastic changes and respond quickly to changes in the main electricity consumption trend because of its specialized equation structure. Overall error analysis indicators also show that AHM often obtains more precise forecasting results than the other two models.

#### 1. Introduction

Given that electricity cannot be stored at a large scale, forecasting electricity consumption is important for power enterprises to balance supply and demand. Considering the differences among forecasting indicators, researches on electricity consumption forecasting are usually divided into two categories, namely, short-term and mid/long-term. Annual electricity consumption forecasting, which is usually used for power system planning, is categorized as mid/long-term.

For most developed countries that have achieved economic development at the postindustrialization stage and have fulfilled most of the electricity needs of people, their electricity consumption changes slowly and can be easily forecasted. By contrast, with the rapid development of the social economy, electricity consumption of developing countries usually changes rapidly and is difficult to forecast [1]. Since 2001, the average yearly growth rate of the top 15 developed electricity consumers of the world is only 1.1%, whereas that of 15 highest electricity consumers in developing nations or countries is as high as 7.4% [2]. This study mainly focuses on the annual electricity consumption forecasting for developing countries.

At present, two ideas for annual electricity consumption forecasting have been proposed. Some researchers had focused on the relationship between annual electricity consumption and its influencing factors and designed multivariate equations to simulate this relationship, obtaining the forecasting results by inputting the values of the explanatory variable in the forecasting period into the equation [3–5]. However, as the exact values of the explanatory variable in the forecasting period are unknown, most models are usually only used for scenario forecasting. Many other researchers concentrated only on the change regularity of the annual electricity consumption trend. Equations are often used to simulate the annual electricity consumption curve and obtain the forecasting results by trend extrapolation [6–8]. These methods can obtain the “natural” change results of future electricity consumption. The adaptive hybrid model (AHM) proposed in this paper is a trend extrapolation method.

The long-term annual electricity consumption curves for most developing countries usually present increasing trends with different curvatures [2]. That is, these annual electricity consumption curves are often approximately exponential growth curves. The artificial neural network (ANN) and other similar models have been extensively used in trend extrapolation of short-term or monthly electricity consumption time series [9–11]. However, compared with these time series, annual electricity consumption has fewer data and, more importantly, is easily affected by socioeconomic environment changes. As such, many minor unexpected waves usually exist in its curve. In this manner, if ANN is used to extrapolate the annual electricity consumption, then the generalization ability of the trained network will easily be affected significantly. Furthermore, time series algorithms [12, 13] cannot easily handle situations wherein the sample size is small. The grey models (GM), which are represented by GM(1, 1), have been successfully used in the trend extrapolation of the annual time series [14–16]. The merits of GM(1, 1) are in the following two aspects. First advantage is that GM(1, 1) can use small samples to estimate the equation parameters. Second, by adjusting the parameters, GM(1, 1) can simulate the approximately exponential growth curve using any slope and curvature. However, when used to forecast annual electricity consumption, GM(1, 1) still has the following limitations.

The forecasting equation of GM(1, 1) is essentially a homogeneous exponential equation, which can only perfectly simulate the time series with a steady growth rate. This makes it unable to perfectly adapt the characteristics of the real annual electricity consumption data. For one thing, the feature of steady growth rate does not necessarily exist in annual electricity consumption data. For another, with the changes of the socioeconomic environment, the annual electricity consumption trend usually cannot remain steady for a long time. Except for the approximately exponential growth trend, the linear trend also usually appears in specific periods. Evidently, GM(1, 1) cannot simulate this linear trend. Furthermore, as mentioned previously, the forecasting equation of GM(1, 1) is essentially a homogeneous exponential equation. The structure of this equation has the upper limit of growth speed. As the social economy develops, the electricity consumption share in end-use energies generally increases gradually. As a result, the annual electricity consumption trend is usually steeper than the total energy consumption and CO_{2} emissions trends [17]. In the period of economic boom, the electricity consumption trend will be very steep and the growth rate of it may increase gradually. Evidently, GM(1, 1) cannot handle this situation.

In this paper, the AHM is first proposed to forecast annual electricity consumption by trend extrapolation. The forecasting equation of this model contains linear, nonhomogeneous exponential, and steep items. By adjusting the equation parameters using small samples, this equation can simulate the above three trends adaptively at the same time.

#### 2. Hybrid Equation

This section discusses the hybrid equation which can simulate linear, nonhomogeneous exponential, and steep trends at the same time.

The fitting and forecasting equation of GM(1, 1) to [14–16] is then written aswhere and are the equation parameters, respectively.

Notably, (1) is a homogeneous exponential equation. Based on its analytical expression, the following relationship between and should exist:

As a result, (1) can be written as

Furthermore, as previously discussed, the steep trend that has an increasing growth rate may also appear in the annual electricity consumption time series. Considering that a linear growth rate can ensure a highly steep trend, the simulation equation of the steep trend can be written as follows:where and are the equation parameters.

That is,

Based on (3) and (5) and with the assumption that the linear trend may also appear in the annual electricity consumption time series, the hybrid equation can be written as

The first part of (6) shows the iterative relationship between two adjacent values, while the second part ascertains the starting point of iteration. The iterative equation obviously contains linear, exponential, and steep items. This equation structure has the ability to simulate the three trends at the same time. Therefore, (6) is more applicable than (1) in forecasting the annual electricity consumption time series.

#### 3. Parameter Estimation

##### 3.1. Parameter Estimation of the Iterative Equation

Influenced by the changes in the socioeconomic environment, the annual electricity consumption time series cannot easily maintain a steady trend for a long time. That is, collecting a large number of samples to estimate the parameters of the iterative equation is nearly impossible. Moreover, considering the unique structure of the iterative equation, multicollinearity between the independent variables is serious, especially when the sample size is small. Multicollinearity will render the estimations obtained by the traditional ordinary least squares (OLS) as highly unstable. As such, the forecasting results will be with significant risks. Similar to OLS, the partial least squares (PLS) can also estimate the parameters of a linear equation. However, unlike OLS, PLS has the ability to determine the regression parameters using few observations with multicollinearity [18]. Hence, this study adopts the idea of PLS to design the parameter estimation algorithm of the previously presented iterative equation.

The key idea of PLS is to extract components to overcome the influence of multicollinearity. Before component extraction, preprocessing of the observations is essential.

Given the structure of (6), letthe extracted component is obtained using the following equation [19, 20]: where and .

Using as the independent variable to explain and , the residuals are written as and , respectively. That is, where and .

Replacing and in (8) with and , respectively, the second component is obtained. Using a similar method, and are also obtained.

More extracted components do not necessarily result in better forecasting precision. In fact, if the current components already contain useful information, more components will introduce stochastic information into the final equation and accordingly reduce the generalization capability of the equation, although simulation precision may be improved. From the second component , the importance of each component should be tested first before being introduced into the final equation. The indicator to evaluate the importance of a component is defined as follows:where is the forecasting result of by the forecasting equation, which is established using the first components and samples (the th sample is excluded), and is the forecasting result of by the forecasting equation, which is established using the first components and all samples.

Based on statistical experience, when , the th component should be introduced into the final equation [21]. Given that the importance of the extracted components decreases one after the other, if the th component is identified as unnecessary to be added into the equation, then the subsequent components ( and ) will unnecessarily be tested and should be abandoned directly.

Given that components remained, the explanatory equation to is written as

As the extracted components have only slight multicollinearity, the regression parameters of (11) can be obtained by OLS. Using the inverse operation of preprocessing and component extraction, (11) can be transformed into the iterative equation of (6).

##### 3.2. Parameter Estimation of the Starting Point

Equation (6) shows that the starting point of the iterative equation will affect the long-term simulation and forecasting results. To obtain the optimum estimation of which can decide the starting point, the transformation to (6) is necessary.

Let (6) is further written as

Let (13) can be written as

Based on the general optimization idea, the optimum estimation of should ascertain the residual sum of squares (RSS) to reach its minimum value.

Let be the RSS. The first derivative of with respect to should be equal to 0 to obtain the minimum value. That is,

By solving (18), the following can be obtained:

As such, the estimations of , , , , and are all obtained. Let in (6), the simulation and forecasting results of the annual electricity consumption time series will be obtained.

#### 4. Model Test

##### 4.1. Data Selection and Preliminary Analysis

China is the largest electricity consumer among the developing countries. At present, the share of global electricity consumption in China is as high as 23.18% [2]. More importantly, electricity consumption in China is still in the phase of rapid growth. The precise forecasting of electricity consumption in China is important for the production planning of its power system. In this study, the electricity consumption data in China are used to test the performance of AHM.

Figure 1 shows the electricity consumption curve from the time of the founding of the People’s Republic of China.

As shown in Figure 1, the long-term electricity consumption trend in China is an approximately exponential growth curve, except for some minor irregular waves. The growth speed of the curve is unstable, with linear or accelerated growth rate trends in certain periods. For example, from 1991 to 1997, the annual increase amounts are 56.18, 76.70, 82.24, 91.45, 79.07, 72.41, and 54.85 TWh. These relatively stable increase amounts imply that the electricity consumption curve during this period nears a linear trend. However, the electricity consumption curve suddenly becomes very steep in the next several years. From 1998 to 2003, the annual growth rates are 2.07%, 6.52%, 10.98%, 8.43%, 10.06%, and 16.53%. That is, the electricity consumption curve roughly presents a steep trend with accelerated growth rates in this period. These above characteristics preliminary prove the necessity of the equation structure of AHM.

##### 4.2. Forecasting Results

To evaluate forecasting performance, ANN, GM(1, 1), and the proposed AHM were all used to one-step-ahead forecast China’s annual electricity consumption from 1991 to 2014, and then their results were compared.

Among the ANNs, the back-propagation (BP) and radial basis function (RBF) are the most represented models. While the latter is usually better than the former in terms of learning efficiency and stability [10], the RBF ANN was selected for this study. Before forecasting, 20 consecutive previous values of the forecasting point were selected to train the network. Among the 20 selected values, the 11th to 20th values were selected as the output values, respectively, while ten previous values of each output value were selected as its input vector. After training, the forecasting result was obtained by inputting the 11th to 20th values as the input vector into the trained network.

In forecasting using GM(1, 1) and AHM, 10 consecutive previous values of each forecasting point were selected to build the forecasting equation. For example, in forecasting electricity consumption in China in 1991, the electricity consumption values from 1981 to 1990 were selected to estimate the equation parameters of GM(1, 1) [ and in (1)] and AHM [ to in (6)]. Let in (1) and (6), the forecasting results of GM(1, 1) and AHM on the electricity consumption in China in 1991 were obtained.

The real data (RD, unit: TWh) of the annual electricity consumption in China from 1991 to 2014, the forecasting results of the three models (unit: TWh), and the relative error (RE, unit: %) of each forecasting result are all listed in Table 1.

Table 1 shows that, compared with GM(1, 1) and AHM, the large errors of ANN appear more irregular. This finding implies that ANN is more easily affected by the stochastic changes of the annual electricity consumption curve than the other two models. GM(1, 1) and AHM usually present large errors at the same time. Furthermore, compared with GM(1, 1), AHM can always more quickly reduce the forecasting errors. This finding proves the adaptive ability of AHM. The curves of the RD and the forecasting results of the three models are drawn in Figure 2 to clearly demonstrate the forecasting performance.

As shown in Figure 2, the rate of increase of the electricity consumption curve decreased unexpectedly in 1998 because of the Southeast Asian financial crisis and returned to its main trend after 2000. This can be identified as a stochastic change. Figure 2 shows that this stochastic change significantly affected the forecasting precision of ANN and resulted in continuous large errors. By contrast, GM(1, 1) and AHM were less affected by this unexpected change, especially AHM, which showed an immediate adjustment of the forecasting result and smaller errors were obtained since 1999.

In 2008, the rate of increase of the electricity consumption decreased again. Unlike in 1998, the RD curve in Figure 2 shows that the change caused a downshift to the main trend. The forecasting precision of ANN was again significantly affected. The forecasting results of AHM, by contrast, exhibited an immediate response, with the values becoming close to the RD since 2010. Compared with AHM, the response of GM(1, 1) is evidently slower, and this makes it have larger forecasting errors.

The performance of AHM can be explained by its specialized forecasting equation. The basic items comprising the forecasting equation of AHM are all constant trend equations, thus allowing AHM to overcome the stochastic effects to the electricity consumption curve. Moreover, at the change in the main electricity consumption trend, AHM can quickly adapt through the adjustment of the parameters of the different items.

##### 4.3. Error Analysis

To comprehensively evaluate the performance of the three forecasting models, the mean absolute percentage error (MAPE), median absolute percentage error (MdAPE), and maximal absolute percentage error (MaxAPE) were used as the indicators to measure the forecasting precision. These indicators are defined as follows:

These indicators have similar functions in evaluating the performance of each forecasting model, but they still have fine distinctions. MAPE is the most widely used indicator, which can reflect the general closeness of the forecasting results to the RD. In addition to having functions similar to those of MAPE in measuring the general closeness of the forecasting results to the RD, MdAPE can overcome the influence of several outliers. MaxAPE records the worst forecasting result and reflects the maximal forecasting risk.

Using the data in Table 1 and (20) the results for the error analysis based on the three indicators were obtained and are shown in Table 2.

As shown in Table 2, AHM is better than ANN and GM(1, 1) based on each indicator. This finding proves that compared with ANN and GM(1, 1), the forecasting results of AHM are not only generally closer to the real data, but also have less forecasting risk.

#### 5. Conclusions

This study first presented a small-sample AHM for annual electricity consumption forecasting. The proposed iterative equation, which contained linear, exponential, and steep items, was used for trend extrapolation of the possible changes of the electricity consumption curve. The PLS algorithm was employed to estimate the parameters of the iterative equation to improve forecasting stability using small samples. Furthermore, the parameter of the iteration starting point was obtained by minimizing the RSS. To evaluate forecasting performance, ANN, GM(1, 1), and the proposed AHM are used to forecast China’s electricity consumption from 1991 to 2014, and then their results are compared. With the forecasting equation of the proposed AHM composed of three constant trend terms, the forecasting results of AHM are less affected by stochastic change than those of ANN. Moreover, at the change in the main electricity consumption trend, the forecasting results of the proposed AHM showed a more immediate adjustment than GM(1, 1) because of the specialized equation of AHM that can simulate three trends at the same time. To comprehensively evaluate the forecasting precision, MAPE, MdAPE, and MaxAPE of each model were calculated. All indicators showed that AHM is better than ANN and GM(1, 1).

#### Notations

: | The th values of the annual electricity consumption time series |

: | The fitting/forecasting result to |

: | The standardized vector of observations of independent variable |

: | The standardized matrix of observations of dependent variables |

: | The extracted component which is the weighted combination of the columns of |

: | The residual vector obtained by using to explain |

: | The residual matrix obtained by using to explain |

: | The indicator to evaluate the importance of a extracted component. |

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The authors are grateful for the financial support provided by the National Natural Science Foundation of China (NSFC) (71471061 and 71201057).