Next Article in Journal
Utilization of Electric Vehicles for Vehicle-to-Grid Services: Progress and Perspectives
Next Article in Special Issue
On the Use of Quantum Reinforcement Learning in Energy-Efficiency Scenarios
Previous Article in Journal
Switching Frequency Effects on the Efficiency and Harmonic Distortion in a Three-Phase Five-Level CHBMI Prototype with Multicarrier PWM Schemes: Experimental Analysis
Previous Article in Special Issue
Use of Artificial Neural Networks to Predict Fuel Consumption on the Basis of Technical Parameters of Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector

by
Felipe Leite Coelho da Silva
1,
Kleyton da Costa
2,
Paulo Canas Rodrigues
3,
Rodrigo Salas
4 and
Javier Linkolk López-Gonzales
5,6,*
1
Department of Mathematics, Federal Rural University of Rio de Janeiro, Seropédica 23890-000, Brazil
2
Department of Economics, Federal Rural University of Rio de Janeiro, Seropédica 23890-000, Brazil
3
Departament of Statistics, Federal University of Bahia, Salvador 40170-110, Brazil
4
Escuela de Ingeniería C. Biomédica, Universidad de Valparaíso, Valparaíso 2362905, Chile
5
Facultad de Ingeniería y Arquitectura, Universidad Peruana Unión, Lima 15, Peru
6
Instituto de Estadística, Universidad de Valparaíso, Valparaíso 2360102, Chile
*
Author to whom correspondence should be addressed.
Energies 2022, 15(2), 588; https://doi.org/10.3390/en15020588
Submission received: 10 December 2021 / Revised: 30 December 2021 / Accepted: 10 January 2022 / Published: 14 January 2022

Abstract

:
Forecasting the industry’s electricity consumption is essential for energy planning in a given country or region. Thus, this study aims to apply time-series forecasting models (statistical approach and artificial neural network approach) to the industrial electricity consumption in the Brazilian system. For the statistical approach, the Holt–Winters, SARIMA, Dynamic Linear Model, and TBATS (Trigonometric Box–Cox transform, ARMA errors, Trend, and Seasonal components) models were considered. For the approach of artificial neural networks, the NNAR (neural network autoregression) and MLP (multilayer perceptron) models were considered. The results indicate that the MLP model was the one that obtained the best forecasting performance for the electricity consumption of the Brazilian industry under analysis.

1. Introduction

In recent years, projections of electricity consumption for the Brazilian industrial sector have been studied, both for short and long term [1]. This interest is related to the development of the sector, energy planning, and energy efficiency [2]. Furthermore, electricity has economic and social importance for a country or region. Thus, considering the industrial sector as one of the largest electricity consumers, studies must be carried out to ensure a minimum of predictability for legislators and consumers’ decision-making processes [3].
In this context, several models have been used to obtain electricity predictions, such as the regression models using only weather variables for predicting load demand in England and Wales [4]; linear regression models for electricity consumption projections in Italy [5]; the Box and Jenkins models as well as the exponential smothing models for electricity demand in European countries [6]; the neural network models for power load forecast in Brazil [7]; the Bayesian dynamic linear model for short-term forecasting of Brazilian industry electricity consumption [8]; additive semi-parametric models for energy load forecasting in Australia [9]; bottom–up model for electricity consumption in Taiwan’s cement industry [10]; bottom–up approach for electricity consumption forecasting of the pulp and paper sector of the Brazilian industry [1]; and bottom–up stochastic approach for electricity consumption forecasting of a sector of the Brazilian industry [11]. Martínez-Álvarez et al. [12] used data-mining techniques to electricity demand forecasting; a comparative study of different time-series models for energy consumption forecasting of smart buildings in a university campus in the south of Spain [13]; there was also a study about energy consumption forecasting to an industrial building using an artificial neural network (ANN) algorithm [14], and artificial intelligence techniques were used for energy demand planning in smart homes [15]. Sulandari et al. [16] used singular spectrum analysis, fuzzy systems, and neural networks for electricity load time-series forecasting in Indonesia. Sulandari et al. [17] presented an SSA-based hybrid forecasting method for a complex seasonal time series of daily electricity load in Indonesia. Sulandari et al. [18] presented a study comparative with the methods of Singular Spectrum Analysis, fuzzy systems, and neural networks for Indonesian electricity load demand forecasting.
The time-series forecasting models are mathematical and computational modeling strategies used in academic research and to help in public policies, which are oriented by evidence. The literature concerning time-series forecasting has to a heterogeneous and dynamic degree taken into account an extensive amount of scientific production and competitions for evaluating models applied to observed data in different fields of knowledge [19,20,21]. It is possible to find in the literature prediction comparisons in which univariate time-series models are superior to (or as good as) multivariate models (as in the discussion proposed in [22,23,24,25]). A possible interpretation for this result is a sparse representation, in large-scale models, of the dynamic interactions in a system of variables [26].
The purpose of this work is to perform a comparative study between two classes of models for time-series forecasting (statistical and artificial neural networks) applied to the electricity consumption in the Brazilian industry.
To achieve the goal of this study, the Holt–Winters method, the seasonal autoregressive integrated moving average model (SARIMA), the dynamic linear model, and TBATS (Trignometric Box–Cox transform, ARMA errors, Trend, and Seasonal components) were considered as part of the statistical approaches. For the artificial neural networks approach, we consider the neural network autoregression (NNAR) and the multilayer perceptron (MLP). It is noteworthy that the use of these classes of models within the same comparative study investigating the consumption of electricity in the Brazilian industry was not found in the literature, although, as noted, there are studies that apply prediction models to the data under analysis. Therefore, the development of this work contributes to the literature available in this area, opening space for further discussions and applications.
This study is structured as follows. Section 2 describes the methodology, Section 3 presents the main results and discussion. Finally, Section 4 gives the conclusion and introduces problems for future research.

2. Methodology

The empirical strategy adopted in this work takes into account statistical and artificial neural networks models to forecast the time series of monthly industrial electricity consumption in the Brazilian energy system. The data were extracted from the Time-Series Management System of the Central Bank of Brazil [27]. We split the dataset into the training set (January 1979 to December 2018) for model fit and the test set (January 2019 to December 2020) to assess the predictive ability of the models under consideration. It is important to note that our analysis seeks to measure the predictive capacity of models for short-term predictions. Thus, the test set with 24 (5%) observations is reasonable for investigation.
The most suitable model was selected through the precision metric mean absolute percentage error (MAPE). All statistical analysis and graphical representations were made using R programming language [28].

2.1. Statistical Models

Statistical models have applications in several areas of knowledge and are considered as established models in the forecasting literature. Here, the models that are considered as part of statistical approach are the Holt–Winters method, the SARIMA model, the dynamic linear model, and the TBATS algorithm.

2.1.1. Holt–Winters Method

The Holt–Winters method was proposed with the contributions of [29,30] using exponentially weighted moving averages to update those needed for seasonal adjustment of the mean (trend) and seasonality. This method can be built in an additive or multiplicative way. The additive method is an extension of Holt’s exponential smoothing that captures seasonality and produces exponentially smoothed values for the level of the forecast, the trend of the forecast, and the seasonal adjustment to the forecast, adding the seasonality factor to the trended forecast. The multiplicative method multiplies the trended forecast by the seasonality, producing the Holt–Winters’ multiplicative forecast. The seasonal adjustment for the additive method subtracts a seasonality component from the level equation. For the multiplicative method, there is a division of the series by its seasonality component. Table 1 describes the three smoothing equations (level, trend, and seasonality) and the forecast equation to the Holt–Winters method.

2.1.2. SARIMA

Following the contributions of [31], we consider as the second model used in this work the SARIMA multiplicative model (seasonal autoregressive integrated moving average model), considering level, trend, and seasonality components from simple and seasonal operators (see Equation (1)).
ϕ ( B ) Φ B s d s D Y t = θ ( B ) Θ B s α t ,
where ϕ ( B ) is the simple autoregressive operator, Φ ( B s ) is the seasonal autoregressive operator, d is the simple differenciated operator, s D is the seasonal differenciated operator, θ ( B ) is the simple moving average operator, Θ ( B s ) is the seasonal moving average operator, and α t is a random noise.

2.1.3. Dynamic Linear Model

The dynamic linear model was the third model considered in this study. From a state-space structure, the general form of the dynamic linear model can be written using two main equations: the so-called observational equation (Equation (2)) and the system evolution equation (Equation (3)) [32].
Y t = F t θ t + v t , v t N m 0 , V t
θ t = G t θ t 1 + w t , w t N p 0 , W t
By definition, θ t is the state vector in time t, F t is a known regression matrix of the observational equation, and G t is a known matrix of the evolution of the system evolution equation. It is assumed that V t is a vector of variance associated with the errors of the observational equation and W t is a matrix of covariance associated with the errors of the equation of evolution of the system.
As mentioned, the matrices F t and G t are known matrices, being used to build the structure of the model according to the components found in the series to be analyzed—level, trend, and seasonality. However, in this study, we performed the estimation of the observational variance matrices V t and the covariance of the system W t using the Monte Carlo Markov Chain (MCMC) method, considered the Gibbs sampler algorithm.

2.1.4. Trignometric Box–Cox Transform, ARMA Errors, Trend, and Seasonal Components (TBATS)

Finally, the TBATS model was the fourth model considered as part of the statistical approach. This method uses a combination of Fourier terms with an exponential smoothing state-space model and Box–Cox transformation, being useful for adjusting the seasonality change over time. Proposed by [33], the model can be written in its reduced form as
φ p ( L ) η ( L ) y t ( ω ) = θ q ( L ) δ ( L ) ε f ,
where L is a lag operator, η ( L ) = det I F * L , δ ( L ) =   w * adj I F * L g * L + det I F * L , φ p ( L ) , and θ q ( L ) are polynomials with degrees p e q, w * = ( 1 , ϕ , a ) , g * = ( α , β , γ ) , with the matrix F * defined by
F * = 1 ϕ 0 0 ϕ 0 0 0 A .

2.2. Artificial Neural Networks Approach

An artificial neural networks seeks to model the relationship between a set of input signals and an output signal, which was inspired by the working mechanism of a biological brain. In this study, we use autoregressive neural networks (NNAR) and multilayer perceptron (MLP).

2.2.1. Autoregressive Neural Networks (NNAR)

The N N A R ( p , P , k ) m model takes into account a feedforward network with a single hidden layer, p inputs, k nodes in the hidden layer, P seasonal lags, and m periods [34,35]. We consider the algorithm proposed by [36] that defines the number of nodes in the hidden layer (k) as an average of the number of inputs and the number of outputs, that is, ( p + P + 1 ) / 2 . Thus, the model is capable of capturing the time-series components (level, trend, and seasonality). In this NNAR model, the inputs into each hidden layer neuron are combined linearly (Equation (6)) to give weight and produce output from artificial neural networks and the activation function as the binary sigmoid, which is a nonlinear function (Equation (7)).
z j = α j + i = 1 N w i , j x i ,
z j represents the jth hidden layer neuron, N represents the number of input layer neurons, α j represents the intercept of the j-th hidden neuron, w i , j denotes the weights assigned to the connection between the input and the hidden layer, x i are the observations (covariates or neurons) of the input layer, and the activation function is given by
g ( z ) = 1 1 + e z .

2.2.2. Multilayer Perceptron (MLP)

Multilayer perceptron or feedforward deep network is a mathematical function mapping a sort of input values to output values. According to [37], the main objective of a feedforward neural network is to approximate any function f * , defining a mapping y = f ( x ; θ ) , and learn the value of the parameter θ = ( w j , v 1 , , v N ) that makes the best function approximation.
According to [38], the MLP is given by
y ^ = v 0 + j = 1 N v j g ( w j T x )
where x is the input vector x, augmented with 1, i.e., x = ( 1 , x T ) T , w j is the weight vector for the jth hidden node, v 0 is the output layer intercept, v 1 , , v N are the weights for the output node, y ^ is the network output, and g is the activation function and is used to allow a possible nonlinearity at the hidden layer. We used the logistic activation function (Equation (7)).

2.3. Mean Absolute Percentage Error (MAPE)

In this study, MAPE was used to assess the ability of models fitted and forecast the time series. The MAPE was defined by
M A P E = 100 × 1 N i = 1 N y i y ^ i y i ,
where N is the number of fitted points, y i is the observed value, and y ^ i is the forecast value both for instant i.

3. Results and Discussion

Figure 1a shows the behavior of electricity monthly consumption of the Brazilian industry, and the within-year variability can be seen in the yearly box-plots of Figure 1b. A greater variability can be observed in years associated with crisis, such as 2001, 2009, and 2020. In 2001, there was a crisis in the Brazilian electrical system, and in 2009, there were the consequences of the world economic crisis that started in 2008 and occurred in Brazil in 2009. The year 2020 had the greatest variability due to the COVID-19 pandemic.
For its part, Table A1 shows some descriptive statistics by year of the industrial electricity consumption. It can be seen in the table that there are several levels of variability in electricity consumption over the years.
Using the mean of squared errors between additive and multiplicative methods, the Holt–Winters multiplicative method obtained a better result than the additive method. Thus, it was decided to consider the multiplicative method for comparison with the other proposed models, and there is a division of the series by its seasonality component.
Then, we applied a selection algorithm to the SARIMA model that considers the candidate models through the principle of parsimony. Thus, the structure with the best result for the Akaike Information Criterion metric was the SARIMA ( 1 , 1 , 1 ) × ( 1 , 1 , 1 ) 12 . We considered the dynamic linear model (DLM) that captures trend and seasonality components. In this work, the dynamic regression matrix F t and the evolution matrix G t of the model are
F t 1 × 13 = 1 0 1 0 0
G t 13 × 13 = 1 1 0 0 0 0 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 .
For DLM, it was assumed that the observational variance V t = σ 2 , and the covariance matrix of the system W t is a diagonal matrix introduced by W t = d i a g ( σ μ 2 , σ β 2 , σ γ 2 , 0 , 0 ) . These unknown variances were also estimated using Bayesian inference. Thus, to complete the specification of the model, we assumed independent inverse gamma priors distributions with means a , a θ 1 , a θ 2 , a θ 3 and variances b , b θ 1 , b θ 2 , b θ 3 , respectively, fixed in known values. Therefore, by using the unobservable states as latent variables, a Gibbs sampler can be run on the basis of the following full conditional densities:
   σ 2 I G a 2 b + n 2 , a b + 1 2 S S y , σ μ 2 I G a θ 1 2 b θ 1 + n 2 , a θ 1 b θ 1 + 1 2 S S θ 1 , σ β 2 I G a θ 2 2 b θ 2 + n 2 , a θ 2 b θ 2 + 1 2 S S θ 2 , σ γ 2 I G a θ , 3 2 b θ 3 + n 2 , a θ 3 b θ 3 + 1 2 S S θ 3 ,
with S S y = t = 1 n y t F t θ t 2 and S S θ i = t = 1 T θ t , i ( G t θ t 1 ) i 2 , for i = 1 , 2 , 3 . The full conditional density of the states is a normal distribution, and it is covered in the used dlm package [39].
Regarding the approach of artificial neural networks, the applied NNAR model has the results as an average of 20 networks, each of which is a NNAR(2,2,1) with nine weights. We also considered an MLP model with five nodes in the input layer, four nodes in the first hidden layer, two nodes in the second hidden layer, and one node in the output layer (Figure 2).
In that sense, Figure 3 presents the fitted results for all models considered. The models are able to fit the observed data even in periods of economic crisis. It can be seen in this figure that the models are able to capture the behavior of the time series of industrial electricity consumption.
Indeed, Figure 4 shows the forecasting results for the models under consideration. It can be seen in the figure that the models were unable to give an accurate forecast for the months of April, May, and June of 2020. The electricity consumption in these months was impacted by the COVID-19 pandemic. This period marked the beginning of the pandemic in Brazil.
From this perspective, Table 2 shows the mean absolute percentage errors (MAPE) for all models considering the training and testing data of the industrial electricity consumption time series. The MLP model provided the best MAPE results for the training data (or model fit) and test data (or forecasting). The fitted models to electricity consumption data provided a percent mean square error of less than 2.5 % .
Finally, Table 3 gives the MAPE for each forecast horizon yielded by all models under consideration. In this step, the fitted models were used to obtain the forecast up to h steps ahead ( h = 1 , 2 , 3 , , 24 ) , and then, the MAPE was calculated for each forecast horizon. The bold entries identify the model that performs best for the corresponding level and forecast horizon, based on the smallest value of MAPE. The last row presents the average MAPE considering all forecast horizons. It can be seen in the table that the Holt–Winters model presented better results for MAPE up to the horizon of eight steps ahead, with the exception of the horizon of three steps ahead, in which the TBATS model resulted in a better result. However, from the horizon of nine steps ahead, the MLP model presented the best results for the MAPE. The MLP model presented the best average MAPE among all the models considered.
To complement Table 3, Figure 5 shows the behavior of MAPE for the forecast horizon of all six models considered in this study. It can be seen in this figure that all models showed a significant increase in the MAPE values after the forecast horizon of 16 steps ahead. This increase can be associated to the beginning of the COVID-19 pandemic in Brazil. Another relevant feature of this figure was the decrease of the MAPE after July 2020 for all models.

4. Conclusions

The exercise of forecasting the industry’s electricity consumption is essential for energy planning in a given country or region. In this way, this study applied time-series forecasting models (statistical approaches and artificial neural network approaches) to the industrial electricity consumption in the Brazilian system.
The results of the study indicate the following: (i) the models considered have a satisfactory ability to adjust to the data; (ii) the models managed to capture the complex structure of the data involving the crises (peaks in the series) in the years 2001 and 2009; and (iii) the MLP model presented the best predictive capacity for the group of proposed models, among which the Holt–Winters method was the overall best for short-term forecasting. In addition, the results found are useful as instruments to support decision making by economic agents and legislators of the Brazilian energy system.
For future research, we can apply other univariate models, such as Recurrent Neural Network (RNN), Long Short-Term Memory (LSTM) [40,41], and bootstrap-based models, and develop a hierarchical time-series forecast strategy to compare with the classical structure applied in this study. Finally, another avenue for future work would be to study the multivariate models (such as vector autoregressive, Bayesian vector autoregressive, artificial neural networks models with multiple features in the input layer) and compare with the univariate approach.

Author Contributions

Conceptualization, F.L.C.d.S., K.d.C., P.C.R. and J.L.L.-G.; methodology, F.L.C.d.S.; software, F.L.C.d.S. and K.d.C.; validation, F.L.C.d.S., P.C.R., J.L.L.-G. and R.S.; formal analysis, F.L.C.d.S. and K.d.C.; investigation, F.L.C.d.S., K.d.C. and J.L.L.-G.; resources, R.S., J.L.L.-G. and P.C.R.; data curation, J.L.L.-G. and F.L.C.d.S.; writing—original draft preparation, F.L.C.d.S., J.L.L.-G., R.S. and P.C.R.; writing—review and editing, J.L.L.-G., R.S., F.L.C.d.S. and P.C.R.; supervision, J.L.L.-G. and F.L.C.d.S.; project administration, J.L.L.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Kleyton da Costa acknowledges financial support from the Brazilian National Council for Scientific and Technological Development (CNPq) grant Programa Institucional de Bolsas de Iniciação Científica (PIBIC). Javier Linkolk López-Gonzales acknowledges financial support from the ANID scholarship.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Summary of statistical description.
Table A1. Summary of statistical description.
YearMeanVarianceSt. Dev.AmplitudeMin.Max.
19794616.8370,024.88264.62725.004215.004940.00
19805123.8365,639.24256.20699.004806.005505.00
19815095.8311,991.61109.51329.004906.005235.00
19825324.0874,528.27273.00809.004836.005645.00
19835669.92142,872.81377.991076.004994.006070.00
19846704.58254,915.17504.891578.005834.007412.00
19857570.00100,401.82316.86890.007084.007974.00
19868094.83286,936.52535.661529.007190.008719.00
19878116.9267,643.36260.08797.007749.008546.00
19888377.2545,443.66213.18695.008073.008768.00
19898583.08246,853.36496.841705.007595.009300.00
19908322.58262,576.99512.421949.007145.009094.00
19918550.08391,295.36625.541814.007466.009280.00
19928610.5893,964.45306.541139.007953.009092.00
19938915.08145,204.81381.061083.008290.009373.00
19948921.92136,503.90369.461082.008368.009450.00
19959305.5044,580.09211.14577.009070.009647.00
19969709.67217,869.70466.771637.008753.0010,390.00
199710,143.08185,066.81430.191272.009455.0010,727.00
199810,164.83152,053.79389.941178.009545.0010,723.00
199910,324.33300,464.79548.151607.009257.0010,864.00
200010,940.00163,874.18404.811398.0010,024.0011,422.00
200110,211.50609,427.18780.662160.009178.0011,338.00
200210,635.67247,841.33497.841609.009431.0011,040.00
200310,852.67114,157.70337.871186.0010,345.0011,531.00
200412,846.83291,560.70539.961585.0011,829.0013,414.00
200513,217.33133,968.42366.021105.0012,496.0013,601.00
200613,598.42149,381.17386.501313.0012,851.0014,164.00
200714,530.67222,010.42471.181433.0013,592.0015,025.00
200814,652.83404,200.70635.771995.0013,417.0015,412.00
200913,483.17710,177.42842.722628.0011,924.0014,552.00
201014,956.58380,298.81616.682031.0013,425.0015,456.00
201115,298.00187,278.18432.761386.0014,467.0015,853.00
201215,285.42112,891.36335.991061.0014,567.0015,628.00
201315,390.25197,482.39444.391516.0014,370.0015,886.00
201414,925.4268,056.63260.88723.0014,537.0015,260.00
201514,071.50127,615.73357.231238.0013,327.0014,565.00
201613,687.75174,712.57417.991598.0012,538.0014,136.00
201713,903.92138,837.36372.611211.0013,105.0014,316.00
201814,121.92119,368.63345.501014.0013,525.0014,539.00
201913,858.1771,473.42267.35864.0013,442.0014,306.00
202013,802.421,019,275.901009.592936.0012,173.0015,109.00

References

  1. Silva, F.L.C.; Souza, R.C.; Oliveira, F.L.C.; Lourenco, P.M.; Calili, R.F. A bottom-up methodology for long term electricity consumption forecasting of an industrial sector—Application to pulp and paper sector in Brazil. Energy 2018, 144, 1107–1118. [Google Scholar] [CrossRef]
  2. López-Gonzales, J.L.; Castro Souza, R.; Leite Coelho da Silva, F.; Carbo-Bustinza, N.; Ibacache-Pulgar, G.; Calili, R.F. Simulation of the Energy Efficiency Auction Prices via the Markov Chain Monte Carlo Method. Energies 2020, 13, 4544. [Google Scholar] [CrossRef]
  3. Ardakani, F.J.; Ardehali, M.M. Novel effects of demand side management data on accuracy of electrical energy consumption modeling and long-term forecasting. Energy Convers. Manag. 2014, 78, 745–752. [Google Scholar] [CrossRef]
  4. Taylor, J.W.; Buizza, R. Using weather ensemble predictions in electricity demand forecasting. Int. J. Forecast. 2003, 19, 57–70. [Google Scholar] [CrossRef]
  5. Bianco, V.; Manca, O.; Nardini, S. Electricity consumption forecasting in italy using linear regression models. Energy 2009, 34, 1413–1421. [Google Scholar] [CrossRef]
  6. Taylor, J.; Mcsharry, P. Short-term load forecasting methods: An evaluation based on european data. Power Syst. IEEE Trans. 2007, 22, 2213–2219. [Google Scholar] [CrossRef] [Green Version]
  7. Santana, A.L.; Conde, G.B.; Rego, L.P.; Rocha, C.A.; Cardoso, D.L.; Costa, J.C.; Bezerra, U.H.; Francs, C.R. Predict decision support system for load forecasting and inference: A new undertaking for brazilian power suppliers. Electr. Power Energy Syst. 2012, 38, 33–45. [Google Scholar] [CrossRef]
  8. Sadownik, R.; Barbosa, E.P. Short term forecasting of industrial electricity consumption in brazil. J. Forecast. 1999, 18, 215–224. [Google Scholar] [CrossRef]
  9. Fan, S.; Hyndman, R.J. Short-term load forecasting based on a semi-parametric additive model. IEEE Trans. Power Syst. 2012, 27, 134–141. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, Y.-H.; Chang, Y.-L.; Fleiter, T. A critical analysis of energy efficiency improvement potentials in taiwan’s cement industry. Energy Policy 2016, 96, 14–26. Available online: https://www.journals.elsevier.com/energy-policy (accessed on 25 November 2021). [CrossRef]
  11. Silva, F.L.; Oliveira, F.L.; Souza, R.C. A bottom-up bayesian extension for long term electricity consumption forecasting. Energy 2019, 167, 198–210. [Google Scholar] [CrossRef]
  12. Martínez-Álvarez, F.; Troncoso, A.; Asencio-Cortés, G.; Riquelme, J.C. A Survey on Data Mining Techniques Applied to Electricity-Related Time Series Forecasting. Energies 2015, 8, 13162–13193. [Google Scholar] [CrossRef] [Green Version]
  13. Divina, F.; García Torres, M.; Goméz Vela, F.A.; Vázquez Noguera, J.L. A Comparative Study of Time Series Forecasting Methods for Short Term Electric Energy Consumption Prediction in Smart Buildings. Energies 2019, 12, 1934. [Google Scholar] [CrossRef] [Green Version]
  14. Ramos, D.; Faria, P.; Vale, Z.; Mourinho, J.; Correia, R. Industrial Facility Electricity Consumption Forecast Using Artificial Neural Networks and Incremental Learning. Energies 2020, 13, 4774. [Google Scholar] [CrossRef]
  15. Rocha, H.R.; Honorato, I.H.; Fiorotti, R.; Celeste, W.C.; Silvestre, L.J.; Silva, J.A. An Artificial Intelligence based scheduling algorithm for demand-side energy management in Smart Homes. Appl. Energy 2021, 282, 116–145. [Google Scholar] [CrossRef]
  16. Sulandari, W.; Subanar, S.; Lee, M.H.; Rodrigues, P.C. Time series forecasting using singular spectrum analysis, fuzzy systems and neural networks. MethodsX 2020, 7, 101015. [Google Scholar] [CrossRef] [PubMed]
  17. Sulandari, W.; Subanar, S.; Suhartono, S.; Utami, H.; Lee, M.H.; Rodrigues, P.C. SSA-based hybrid forecasting models and applications. Bull. Electr. Eng. Inform. 2020, 9, 2178–2188. [Google Scholar] [CrossRef]
  18. Sulandari, W.; Lee, M.H.; Rodrigues, P.C. Indonesian electricity load forecasting using singular spectrum analysis, fuzzy systems and neural networks. Energy 2020, 190, 116408. [Google Scholar] [CrossRef]
  19. Makridakis, S.; Hibon, M. The m3-competition: Results, conclusions and implications. Int. J. Forecast. 2000, 16, 451–476. [Google Scholar] [CrossRef]
  20. Makridakis, S.; Spiliotis, E.; Assimakopoulos, V.; Chen, Z.; Gaba, A.; Tsetlin, I.; Winkler, R. The m5 Uncertainty Competition: Results, Findings and Conclusions. 2020. Available online: https://www.researchgate.net/publication/346493740_The_M5_Uncertainty_competition_Results_findings_and_conclusions (accessed on 15 November 2021).
  21. Makridakis, S.; Fry, C.; Petropoulos, F.; Spiliotis, E. The Future of Forecasting Competitions: Design Attributes and Principles. 2021. Available online: https://arxiv.org/abs/2102.04879 (accessed on 1 November 2021).
  22. Du Preez, J.; Witt, S. Univariate versus multivariate time series forecasting: An application to international tourism demand. Int. J. Forecast. 2003, 19, 435–451. Available online: https://www.sciencedirect.com/science/article/pii/S0169207002000572 (accessed on 1 December 2021). [CrossRef]
  23. Todorov, H.; Searle-White, E.; Gerber, S. Applying univariate vs. multivariate statistics to investigate therapeutic efficacy in (pre)clinical trials: A Monte Carlo simulation study on the example of a controlled preclinical neurotrauma trial. PLoS ONE 2020, 15, 798. [Google Scholar] [CrossRef]
  24. Rana, M.; Koprinska, I.; Agelidis, V. Univariate and multivariate methods for very short-term solar photovoltaic power forecasting. Energy Convers. Manag. 2021, 121, 380–390. Available online: https://www.sciencedirect.com/science/article/pii/S0196890416303934 (accessed on 2 December 2021). [CrossRef]
  25. Ivanova, M.; Herron, T.; Dronkers, N.; Baldo, J. An empirical comparison of univariate versus multivariate methods for the analysis of brain–behavior mapping. Hum. Brain Mapp. 2021, 42, 1070–1101. Available online: https://onlinelibrary.wiley.com/doi/abs/10.1002/hbm.25278 (accessed on 4 December 2021). [CrossRef] [PubMed]
  26. Lütkepohl, H. Introduction to Multiple Time Series Analysis; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  27. Central Bank of Brazil. Time Series Management System—v2.1. 2021. Available online: https://www3.bcb.gov.br/sgspub/localizarseries/localizarSeries.do?method=prepararTelaLocalizarSeries (accessed on 10 November 2021).
  28. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 9 November 2021).
  29. Holt, C.C. Forecasting seasonals and trends by exponentially weighted moving averages. Int. J. Forecast. 1957, 20, 5–10. [Google Scholar] [CrossRef]
  30. Winters, P.R. Forecasting sales by exponentially weighted moving averages. Manag. Sci. 1960, 6, 324–342. [Google Scholar] [CrossRef]
  31. Box, G.E.P.; Jenkins, G.M.; Reinsel, G.C.; Ljung, G.M. Time Series Analysis: Forecasting and Control; Wiley: Hoboken, NJ, USA, 2015. [Google Scholar]
  32. West, M.; Harrison, J. Bayesian Forecasting and Dynamic Models; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  33. Livera, A.M.D.; Hyndman, R.J.; Snyder, R.D. Forecasting time series with complex seasonal patterns using exponential smoothing. J. Am. Stat. Assoc. 2011, 106, 1513–1527. [Google Scholar] [CrossRef] [Green Version]
  34. Hyndman, R.J.; Athanasopoulos, G. Forecasting: Principles and Practice; OTexts: Melbourne, Australia, 2018. [Google Scholar]
  35. Veloz, A.; Salas, R.; Allende-Cid, H.; Allende, H.; Moraga, C. Identification of lags in nonlinear autoregressive time series using a flexible fuzzy model. Neural Process. Lett. 2016, 43, 641–666. [Google Scholar] [CrossRef]
  36. Hyndman, R.; Athanasopoulos, G.; Bergmeir, C.; Caceres, G.; Chhay, L.; O’Hara-Wild, M.; Petropoulos, F.; Razbash, S.; Wang, E.; Yasmeen, F. Forecast: Forecasting Functions for Time Series and Linear Models. R Package Version 8.15. 2021. Available online: https://pkg.robjhyndman.com/forecast/ (accessed on 10 November 2021).
  37. Goodfellow, I.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
  38. Ahmed, N.K.; Atiya, A.F.; Gayar, N.E.; El-Shishiny, H.E. An Empirical Comparison of Machine Learning Models for Time Series Forecasting. Econom. Rev. 2010, 29, 594–621. [Google Scholar] [CrossRef]
  39. Petris, G. An R Package for Dynamic Linear Models. J. Stat. Softw. Artic. 2020, 36, 1–16. Available online: https://www.jstatsoft.org/v036/i12 (accessed on 3 November 2021).
  40. Vivas, E.; Allende-Cid, H.; Salas, R. A systematic review of statistical and machine learning methods for electrical power forecasting with reported mape score. Entropy 2020, 22, 1412. [Google Scholar] [CrossRef] [PubMed]
  41. Cordova, C.H.; Portocarrero, M.N.L.; Salas, R.; Torres, R.; Rodrigues, P.C.; López-Gonzales, J.L. Air quality assessment and pollution forecasting using artificial neural networks in Metropolitan Lima-Peru. Sci. Rep. 2021, 11, 24232. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Electricity consumption (a) in GWh and (b) box-plots for the Brazilian industry. Source: Central Bank of Brazil.
Figure 1. Electricity consumption (a) in GWh and (b) box-plots for the Brazilian industry. Source: Central Bank of Brazil.
Energies 15 00588 g001
Figure 2. MLP network architecture.
Figure 2. MLP network architecture.
Energies 15 00588 g002
Figure 3. Model fit for the six considered models applied to the Brazilian industrial electricity consumption.
Figure 3. Model fit for the six considered models applied to the Brazilian industrial electricity consumption.
Energies 15 00588 g003
Figure 4. Model forecasting for the six considered models applied to the Brazilian industrialelectricity consumption.
Figure 4. Model forecasting for the six considered models applied to the Brazilian industrialelectricity consumption.
Energies 15 00588 g004
Figure 5. Mean absolute percentage error, considering the six models, for h = 1 , , 24 steps-ahead out-of-sample forecasts applied to the electricity consumption of the Brazilian industry.
Figure 5. Mean absolute percentage error, considering the six models, for h = 1 , , 24 steps-ahead out-of-sample forecasts applied to the electricity consumption of the Brazilian industry.
Energies 15 00588 g005
Table 1. The smoothing and forecast equations of the Holt–Winters method.
Table 1. The smoothing and forecast equations of the Holt–Winters method.
EquationsAdditive MethodMultiplicative Method
Level ( t ) α y t s t m + ( 1 α ) t 1 + b t 1 α y t s t m + ( 1 α ) t 1 + b t 1
Trend ( b t ) β t t 1 + ( 1 β ) b t 1 β t t 1 + ( 1 β ) b t 1
Seasonal ( s t ) γ y t t 1 b t 1 + ( 1 γ ) s t m γ y t t 1 + b t 1 + ( 1 γ ) s t m
Forecast ( y ^ t + h t ) t + h b t + s t + h m ( k + 1 ) t + h b t s t + h m ( k + 1 )
Table 2. Mean absolute percentage error for the six models under consideration for model fit and model forecasting considering the training and testing data, respectively.
Table 2. Mean absolute percentage error for the six models under consideration for model fit and model forecasting considering the training and testing data, respectively.
ModelFittedForecast
Holt–Winters2.514.09
SARIMA1.886.17
TBATS1.993.77
DLM1.874.09
NNAR2.404.77
MLP1.483.41
Table 3. Mean absolute percentage error, considering the six models, for h = 1 , , 24 steps-ahead out-of-sample forecasts applied to the electricity consumption of the Brazilian industry.
Table 3. Mean absolute percentage error, considering the six models, for h = 1 , , 24 steps-ahead out-of-sample forecasts applied to the electricity consumption of the Brazilian industry.
StepHWSARIMATBATSDLMNNARMLP
10.270.860.502.283.290.96
20.341.560.952.673.461.96
31.001.970.952.622.961.82
40.843.081.713.063.252.49
51.112.771.462.532.812.17
61.432.861.762.652.912.03
71.532.961.792.452.791.82
81.563.221.922.292.651.75
92.033.652.292.362.851.81
101.933.622.132.182.661.73
112.063.802.202.132.691.63
122.444.122.372.333.101.62
132.544.042.292.243.421.53
142.434.062.152.243.541.45
152.374.102.042.183.481.41
163.135.292.973.154.442.44
173.866.213.733.935.253.24
184.526.864.284.485.833.55
194.576.914.214.385.763.38
204.406.834.004.225.473.31
214.286.693.854.155.253.26
224.286.433.894.215.183.39
234.156.343.774.124.983.39
244.006.173.784.094.773.42
Average2.544.352.543.043.872.32
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Leite Coelho da Silva, F.; da Costa, K.; Canas Rodrigues, P.; Salas, R.; López-Gonzales, J.L. Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector. Energies 2022, 15, 588. https://doi.org/10.3390/en15020588

AMA Style

Leite Coelho da Silva F, da Costa K, Canas Rodrigues P, Salas R, López-Gonzales JL. Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector. Energies. 2022; 15(2):588. https://doi.org/10.3390/en15020588

Chicago/Turabian Style

Leite Coelho da Silva, Felipe, Kleyton da Costa, Paulo Canas Rodrigues, Rodrigo Salas, and Javier Linkolk López-Gonzales. 2022. "Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector" Energies 15, no. 2: 588. https://doi.org/10.3390/en15020588

APA Style

Leite Coelho da Silva, F., da Costa, K., Canas Rodrigues, P., Salas, R., & López-Gonzales, J. L. (2022). Statistical and Artificial Neural Networks Models for Electricity Consumption Forecasting in the Brazilian Industrial Sector. Energies, 15(2), 588. https://doi.org/10.3390/en15020588

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