Consumption–Production Profile Categorization in Energy Communities
Abstract
:1. Introduction
- Cornwall LEM was a four-year trial (2016–2020), aligned with the 2050 net-zero economy objective. The project was based on a flexibility solution. The partners involved in the project were Centrica N-side Western Distribution National Grid ESO, University of Exeter, and Imperial College London.
- The Brooklyn Microgrid case is a P2P Microgrid where quarter neighbors sell their surpluses to other neighbors. It started in 2012 because of hurricane Sandy, during which of networks lost electricity. Whole-community residential panels produce 250–400 MWh each unit of small-scale capacity.
- Great Manchester Local Energy Market (GM LEM). This was a two-year project whose mission was to reach a net-zero emission target by 2038. The project was funded by 11 partners (Bruntwood, OVOEnergy, Carbon Co-Op, Hitachi, Daikin, Regent, et al.).
- Gothenburg Communities. The project was financed by the Swedish Energy Agency and f3 (includes universities and research institutes. Chalmers Industrileknik (CIT) is the host of the f3 organization. The Swedish Knowledge Centre for Renewable Transportation Fuels. Gothenburg Communities comprise 16 municipalities, members of the Klimatkommunerna (Climate Municipalities).
- Empower Hvaler, Norway. Empower was the European Union’s Horizon 2020 Research and Innovation program based on the idea that the local market maximizes social welfare. The consortium included Smart Innovation Østfold AS, (NO) Schneider Electric Norge AS (NO), eSmart Systems AS (NO), Fredrikstad Energi Nett AS (NO), University of St. Gallen (CH), Universitat Politècnica de Catalunya (ES), Malta Intelligent Energy Management Agency (MT), and NewEn Projects GmbH (DE) to develop a neighborhood energy using rooftop PV (3–5 kW) on about 100 houses, micro wind turbines, other de-central energy production, and EVs.
- Los Molinos del Rio Aguas (LMRA). Ecological Community in Almería, Southeast of Spain. Off-grid community-owned microgrid interconnecting solar home systems. This LEM is based on Blockchain to exchange energy and contribute to the Community’s social welfare, an example of prosumption. They have achieved a Leveraged Cost of Energy (LCOE) reduction of around 30% (around 0.08 €/kWh).
2. Materials and Methods
- 1.
- Scalability (from small to significant Energy Communities).
- 2.
- LEM Type (namely Transactive Energy, Energy Community, or P2P).
- 3.
- Price Scheme (i.e., how the pricing is computed).
- 4.
- Bidding Strategies chosen (i.e., the market mechanism used to form the price, typically Long Term Auction, Double auction Auction, Intraday Market, Power Purchase Agreement, and Continuous Double Auction/Two-sided bidding market).
- 5.
- The Forecasting Models used.
- 6.
- Storage optimization.
- 7.
- Smart charging (to execute upward and downward flexibility).
- 8.
- Automation and Blockchain technologies implemented.
2.1. Time Series Forecasting
2.2. ARIMA, SARIMA, and SARIMAX Methods
- 1.
- AR(p) Autoregression (AR) is a regression model that uses past values as input to predict future values.
- 2.
- I(d) Integration is the process of differencing the time series to make it stationary, a necessary condition for using ARIMA models. They suppose the error is white noise. Differencing involves subtracting the current values of a series from its previous values d number of times. A Dickey–Fuller test may be run to check the null hypothesis of the series to be stationary.
- 3.
- MA(q) Moving Average the moving average component depicts the model’s error as a combination of previous error terms. The order q represents the number of terms included in the model.
- 1.
- Select the dependent variable: Consumption kWh, Production kWh, or BESS Headroom %. The Consumption kWh measures the energy consumption required to run appliances and heat a particular site at a given time. Production kWh measures the PV generation. Headroom% measures the storage available as a market opportunity.
- 2.
- Prepare the dependent variable (data cleansing/pre-processing activities): handle null values with linear interpolation, detect and replace all kinds of outliers with cut-off values, and reject variables with a quality score less than 25%.
- 3.
- Select exogenous variables and prepare for modeling (run all data cleansing/pre-processing activities). The power flows, storage headroom%, and meteorological measures have been selected as exogenous variables for all models. The list includes the following: Discharge kWh, Charge kWh, Grid Export kWh, PV Charge kWh, PV Consumption kWh, PV Export kWh, Grid Discharge kWh, Grid Charge kWh, Grid Consumption kWh, Consumption Discharge kWh, precipitation, precipitation probability, sunshine hours, solar irradiation, wind speed, and wind direction. These variables are described in detail in the Experimentation and Results Section 3.
- 4.
- Check dependent variable stationarity using residual analysis and Augmented Dickey–Fuller test. The Augmented Dickey–Fuller test (ADF Test) checks the null hypothesis of a series as stationaryADF Test checks if . If so, then is a random walk. It is a stationary process if −1 < 1 + < 1. This test helps us to determine integration (d and D) parameters. Figure 3 shows the Augmented Dickey–Fuller Test applied to the Hourly Production kWh. In this case, the natural series is stationary. The same experiment applied to Daily Production kWh data showed non-stationarity and the need to integrate with one difference to make it stationary.
- 5.
- Run SARIMAX models.
- 6.
- Plot ACF and PACF correlograms for SARIMAX models to determine optimal seasonal and non-seasonal autoregressive and moving average (p, q, P, and Q) parameters. Analyze the overall fitness, cause and effects, and impact diagram in TCM models.
- 7.
- Diagnose and assess random residuals model using the Ljiung–Box Q-test. The Ljiung–Box Q-test contrasts the null hypothesis that the autocorrelations of a time series are different from zero. The null hypothesis tests’ residual errors are not random, which implies that there is a structure in the observed series that the model does not explain. The more random the errors, the more likely it is to be a good model. Instead of testing for randomness on each lag, the total randomness is tested on all study lags. Not passing this test implies not properly decomposing the time series in its trend and seasonality components from the stochastic component.
- 8.
- Analyze all correlogram lags inside the 95% confidence interval to determine more optimal autoregressive and moving average parameters order.
- 9.
- Compute validation metrics to determine model fitness. The validation metrics are RMSE, MSE, MAE, and R2.
- 10.
- Graph residuals and residuals Q-Q Plot to verify residuals’ stationarity and normality.
- 11.
- Compute incumbent forecasting model validation metrics (RMSE, MSE, MAE, and R2) to determine model fitness.
- 12.
- Benchmark SARIMAX performance model improvement with incumbent forecasting validation metrics and plot measured dependent variable vs. incumbent forecasting model estimates vs. SARIMAX model estimates.
Exponential Smoothing Average Method (ESA)
2.3. Temporal–Causal Models
- 1.
- Select and prepare the dependent variable. The dependent variables are the Local Energy Market Consumption kWh, Production kWh, and BESS Headroom %.
- 2.
- Select exogenous variables and prepare for modeling. The variables selected are: Discharge kWh, Charge kWh, Grid Export kWh, PV Charge kWh, PV Consumption kWh, PV Export kWh, Grid Discharge kWh, Grid Charge kWh, Grid Consumption kWh, Consumption Discharge kWh, precipitation, precipitation probability, sunshine hours, solar irradiation, wind speed, and wind direction.
- 3.
- Check dependent variable stationarity using residual analysis and Dickey–Fuller test and determine d and D parameters.
- 4.
- Analyze the overall quality model by computing validation metrics (RMSE, MSE, MAE, and R2). Analyze the test of model effects (how exogenous variables affect the dependent variable and vice versa). Analyze inputs’ impact on the dependent variable. Analyze input variables lags’ significance.
- 5.
- Graph residuals and residuals Q-Q Plot to verify residuals normality.
- 6.
- Compute incumbent forecasting model validation metrics to determine model fitness. Validation metrics have been RMSE, MSE, MAE, and R.
- 7.
- Compare TCM performance model improvement against incumbent forecasting validation metrics and plot measured dependent variable vs. incumbent forecasting model estimates vs. TCM model estimates.
- 1.
- Data preparation and selectionThe Cornwall LEM source dataset crosses the Power Flows and Storage State of Charge of a certain site measured at minute granularity level with the site metadata information, weather forecasts data, and incumbent forecasting model production and consumption estimates. All analyses were created for the annual period of 1 April 2019 to 31 March 2020. Information was then aggregated daily and hourly.
- 2.
- Stationarity Analysis of Dependent VariableThe Augmented Dickey–Fuller Test has been applied to check stationarity and identify the level of integration needed. Figure 3 shows the Augmented Dickey–Fuller Stationarity Test for hourly Production kWh.
- 3.
- Modeling and Assessment
- (a)
- Model IdentificationProduce ESA-SARIMAX correlograms. Figure 4 shows a time series correlogram with 24 lags for an hourly time series. The correlograms present the autocorrelation and partial autocorrelation functions with a 95% confidence interval.
- (b)
- ESA-SARIMAX random residual diagnosis with Ljung–Box Q-Test. Check the absence of serial autocorrelation.
- (c)
- Model AssessmentValidation metrics include Root Mean Square Error and Mean Square Error (which show how close the forecasts are to the actual values), Mean Absolute Error or an average of the absolute values of the errors across all records (it indicates the average magnitude of error, independent of the direction), and the coefficient of determination R, which determines the proportion of variance in the dependent variable that the exogenous variables can explain.
- 4.
- Residual Normality VisualizationThe Q-Q Plot is utilized to test Normality in Residuals of selected forecasting models. The Q-Q plot for Consumption kWh residual is presented in Figure 5. The Quantile–Quantile Plot helps us to compare two distributions. In the Figure, the points represent the residual quantile distribution, while the straight line/45º reference line represents the standard normal distribution. If the quantile distribution points follow a straight line, the time series is normally distributed. In our case, that means the model residuals are normally distributed.
- 5.
- Performance BenchmarkComparison Table with validation metrics from ESA-SARIMAX-TCM models and incumbent forecasting model for Consumption kWh and Production kWh. Data come from the incumbent forecasting models at postcode and quarter-hourly granularity levels.
- 6.
- Error Reduction ReportComparison Table of MSE and MAE Reduction from selected ESA-SARIMAX-TCM forecasting models and the incumbent forecasting models
3. Experimentation and Results
- 1.
- Verify the stationarity of all time series with the Augmented Dickey-Fuller Test (it checks if the dependent variable series is stationary).
- 2.
- Check if the residual autocorrelation is random; that is, if the model has captured the data dynamic process component, e.g., trend, seasonality, and stochastic/irregular, using the Ljiung–Box Q-Test and analyzing the residual normality with graphical analysis and residual Q-Q Plots.
- 3.
- In SARIMAX models, ADF and PACF functions of the model that pass the residual autocorrelation randomness determine the optimal p and q parameters. Significance tests of these parameters measure their suitability. In the Time Causal Model, the Cause and Effects test measures the parameters’ aptitude using significance tests.
- 4.
- Validate the best-selected model with RMSE, MSE, MAE, R metrics.
- 5.
- Benchmark with incumbent validation metrics (using RMSE, MSE, MAE, R metrics) and present the error reduction in terms of MSE and MAE.
3.1. Models for the Cornwall LEM Average Site
- a.
- Average site daily granularity models
- b.
- Average site hourly granularity models
3.2. Models for Meaningful Cornwall LEM Site Clusters
- a.
- Selected cluster daily granularity models
- Most meaningful inputs: Heating Electric/Non-Electric, TV, Freezer, Fridge, Ground source heating pump, Air source heating pump, Solar PV, Floor Area, Dwelling Type, Consumption kWh
- Dwelling Type: bungalow, house
- Floor Area: mostly big (350–400 m2) with a small portion of small dwellings (0–100 m2)
- DER: they have mostly PV panels, air source heat pump, and ground source heat pump
- Appliances: fridge, freezer, TV
- Heating: all of them have Electric Heating
- Power Flows and BESS Headroom% (measured at site and date level):
- –
- Consumption kWh: It ranges between 0 and 157 kWh. The highest frequency is at 110 kWh.
- –
- Grid Consumption kWh: it ranges between 0 and 153 kWh. The highest frequency at 90 kWh
- –
- Grid Import kWh: it ranges between 0 and 153 kWh. The highest frequency at 90 kWh
- –
- PV Consumption kWh: it ranges between 0 and 28.65 kWh. The highest frequency at 7 kWh
- –
- Production kWh: it ranges between 0 and 84.7 kWh. The highest frequency at 40 kWh
- –
- Headroom: it ranges between 0.047 and 100%. The highest frequency is between 40 and 90
- –
- Grid Charge kWh: it ranges between 0 and 27.96 kWh. The highest frequency at 27 kWh
- –
- Discharge kWh: it ranges between 0 and 23.61 kWh. The highest frequency at 24 kWh
- –
- Consumption Discharge kWh: it ranges between 0 and 23.6 kWh. The highest frequency is at 23 kWh.
- –
- Grid Export kWh: it ranges between 0 and 56.51 kWh. The highest frequency at 12 kWh
- –
- PV Export kWh: it ranges between −0.57 and 56.51 kWh. The highest frequency at 12 kWh
- –
- PV Charge kWh: it ranges between 0 and 17.07 kWh. The highest frequency at 11 kWh
- –
- Charge kWh: it ranges between 0 and 28.11 kWh. The highest frequency at 27 kWh
- –
- Grid Discharge kWh: it ranges between 0 and 27.96 kWh. The highest frequency at 5 kWh
- b.
- Selected cluster hourly granularity models
- Most meaningful inputs: Heating Electric/Non-Electric, TV, Freezer, Fridge, Ground source heating pump, Air source heating pump, Solar PV, Floor Area, Dwelling Type, Headroom.
- Dwelling Type: bungalow, house.
- Floor Area: the highest frequency with 300–350 m2.
- DER: they have mostly PV panels and ground source heat pump.
- Appliances: fridge, freezer, TV.
- Heating: they have mostly Electric Heating.
- Power Flows and BESS Headroom% (measured at site and hour level).
- 1.
- Consumption kWh: It ranges between 0 and 16 kWh. The highest frequency at 16 kWh.
- 2.
- Grid Consumption kWh: It ranges between 0 and 16 kWh. The highest frequency is at 16 kWh
- 3.
- Grid Import kWh: It ranges between 0 and 16 kWh. The highest frequency at 16 kWh
- 4.
- PV Consumption kWh: It ranges between 0 and 5.3 kWh. The highest frequency is at 4.2 kWh
- 5.
- Production kWh: It ranges between 0 and 11.2 kWh. The highest frequency is at 7 kWh
- 6.
- Headroom: It ranges between 0.047 and 100%. The highest frequency is between 20 and 90%
- 7.
- Grid Charge kWh: It ranges between 0 and 3.47 kWh. The highest frequency is at 1.5 kWh
- 8.
- Discharge kWh: It ranges between 0 and 3.34 kWh. The highest frequency is at 3.2 kWh
- 9.
- Consumption Discharge kWh: It ranges between 0 and 3.33 kWh. The highest frequency is at 3.2 kWh
- 10.
- Grid Export kWh: It ranges between 0 and 9.95 kWh. The highest frequency is at 6.2 kWh
- 11.
- PV Export kWh: It ranges between -0.038 and 9.95 kWh. The highest frequency is at 6.2 kWh
- 12.
- PV Charge kWh: It ranges between 0 and 3.03 kWh. The highest frequency is at 3.2 kWh
- 13.
- Charge kWh: It ranges between 0 and 3.5 kWh. The highest frequency is at 1.5 kWh
- 14.
- Grid Discharge kWh: It ranges between 0 and 3.21 kWh. The highest frequency is at 2.2 kWh
4. Discussion
- 1.
- Dataset variability: different data attributes collected at every Energy Community may contribute to the model comparison.
- 2.
- Weather forecast information: Consumption and production depend heavily on meteorology. Weather forecast quality and granularity are critical issues in forecasting power flows.
- 3.
- Price signals information: It is necessary to count on a fourth forecasting model: the price signals forecasting model. That way, the LEM will understand when to charge the BESS inexpensively and when to discharge with the highest profit for the Energy Community.
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Braunholtz-Speight, T.; Sharmina, M.; Manderson, E.; Mclachlan, C.; Hannon, M.; Hardy, J.; Mander, S. Business models and financial characteristics of community energy in the UK. Nat. Energy 2020, 5, 169–177. [Google Scholar] [CrossRef]
- Okwuibe, G.C.; Gazafroudi, A.; Hambridge, S.; Dietrich, C.; Trbovich, A.; Shafie-khah, M.; Tzscheutschler, P.; Hamacher, T. Evaluation of Hierarchical, Multi-Agent, Community-Based, Local Energy Markets Based on Key Performance Indicators. Energies 2022, 15, 3575. [Google Scholar] [CrossRef]
- Faia, R.; Pinto, T.; Vale, Z.; Corchado, J.M. A Local Electricity Market Model for DSO Flexibility Trading. In Proceedings of the International Conference on the European Energy Market (EEM), Ljubljana, Slovenia, 18–20 September 2019; pp. 1–6. [Google Scholar] [CrossRef]
- Cheng, L.; Yu, T. Game-Theoretic Approaches Applied to Transactions in the Open and Ever-Growing Electricity Markets from the Perspective of Power Demand Response: An Overview. IEEE Access 2019, 7, 25727–25762. [Google Scholar] [CrossRef]
- Hu, J.; Wu, J.; Ai, X.; Liu, N. Coordinated Energy Management of Prosumers in a Distribution System Considering Network Congestion. IEEE Trans. Smart Grid 2021, 12, 468–478. [Google Scholar] [CrossRef]
- Capper, T.; Gorbatcheva, A.; Mustafa, M.A.; Bahloul, M.; Schwidtal, J.M.; Chitchyan, R.; Andoni, M.; Robu, V.; Montakhabi, M.; Scott, I.J.; et al. Peer-to-peer, community self-consumption, and transactive energy: A systematic literature review of local energy market models. Renew. Sustain. Energy Rev. 2022, 162, 112403. [Google Scholar] [CrossRef]
- Khatib, H. Electricity trading. Economic Evaluation of Projects in the Electricity Supply Industry; Institution of Engineering and Technology: London, UK, 2014; pp. 197–211. [Google Scholar] [CrossRef]
- Kane, D.; Peacock, A.; McCallum, P. LEM Residential BESS Utilisation Summary Report; Technical Report; Trilemma Consulting Limited: Windsor, UK, 2020; Available online: https://www.centrica.com/media/4637/lem-residential-bess-utilisation-summary-report.pdf (accessed on 4 October 2023).
- Kane, D.; Peacock, A. Cornwall Local Energy Market Residential Project a Whistle Stop Tour; Technical Report November; Trilemma Consulting Limited: Windsor, UK, 2020. [Google Scholar]
- Nicholls, D.; Kane, D. Cornwall LEM Residential Electricity Dataset with Solar Production and Battery Storage, 2018–2020; UK Data Service: Colchester, UK, 2021. [Google Scholar] [CrossRef]
- Kane, D.; Peacock, A.; McCallum, P. LEM Residential Data Dictionary PUBLIC; Technical Report; Trilemma Consulting Limited: Windsor, UK, 2020. [Google Scholar]
- Kane, D.; Peacock, A.; McCallum, P. LEM Residential MetaData Summary Report; Technical Report; Trilemma Consulting Limited: Windsor, UK, 2020. [Google Scholar]
- Kane, D.; Peacock, A.; McCallum, P. LEM Residential Fleet Self-Consumption Summary Report; Technical Report; Trilemma Consulting Limited: Windsor, UK, 2020. [Google Scholar]
- Zendehboudi, A.; Baseer, M.A.; Saidur, R. Application of support vector machine models for forecasting solar and wind energy resources: A review. J. Clean. Prod. 2018, 199, 272–285. [Google Scholar] [CrossRef]
- Das, U.K.; Tey, K.S.; Seyedmahmoudian, M.; Mekhilef, S.; Idris, M.Y.I.; Van Deventer, W.; Horan, B.; Stojcevski, A. Forecasting of photovoltaic power generation and model optimization: A review. Renew. Sustain. Energy Rev. 2018, 81, 912–928. [Google Scholar] [CrossRef]
- Wang, H.; Lei, Z.; Zhang, X.; Zhou, B.; Peng, J. A review of deep learning for renewable energy forecasting. Energy Convers. Manag. 2019, 198. [Google Scholar] [CrossRef]
- Lai, J.P.; Chang, Y.M.; Chen, C.H.; Pai, P.F. A survey of machine learning models in renewable energy predictions. Appl. Sci. 2020, 10, 5975. [Google Scholar] [CrossRef]
- Voyant, C.; Notton, G.; Kalogirou, S.; Nivet, M.L.; Paoli, C.; Motte, F.; Fouilloy, A. Machine learning methods for solar radiation forecasting: A review. Renew. Energy 2017, 105, 569–582. [Google Scholar] [CrossRef]
- Yin, S.; Ai, Q.; Li, J.; Li, Z.; Fan, S. Energy Pricing and Sharing Strategy Based on Hybrid Stochastic Robust Game Approach for a Virtual Energy Station With Energy Cells. IEEE Trans. Sustain. Energy 2020, 12, 772–784. [Google Scholar] [CrossRef]
- Box, G.W.; Jenkins, G. Time Series Analysis: Forecasting and Control; Holden-Day: San Francisco, CA, USA, 1970. [Google Scholar]
- Holt, C. Forecasting seasonals and trends by exponentially weighted moving averages. Int. J. Forecast. 2004, 20, 5–10. [Google Scholar] [CrossRef]
- Winters, P.R. Forecasting Sales by Exponentially Weighted Moving Averages. Manag. Sci. 1969, 6, 324–342. Available online: http://www.jstor.org/stable/2627346 (accessed on 4 October 2023). [CrossRef]
- Song, X.; Taamouti, A. A Better Understanding of Granger Causality Analysis: A Big Data Environment. Oxf. Bull. Econ. Stat. 2019, 81, 911–936. [Google Scholar] [CrossRef]
- Granger, C. Investigating Causal Relations by Econometric Models and Cross-Spectral Methods. Econom. Econom. Soc. 1969, 37, 424–438. [Google Scholar] [CrossRef]
- Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 1990. [Google Scholar] [CrossRef]
- Mellit, A.; Pavan, A.M.; Ogliari, E.; Leva, S.; Lughi, V. Advanced methods for photovoltaic output power forecasting: A review. Appl. Sci. 2020, 10, 487. [Google Scholar] [CrossRef]
- Mosavi, A.; Salimi, M.; Ardabili, S.F.; Rabczuk, T.; Shamshirband, S.; Varkonyi-Koczy, A.R. State of the art of machine learning models in energy systems, a systematic review. Energies 2019, 12, 1301. [Google Scholar] [CrossRef]
- Pérez-Ortiz, M.; Jiménez-Fernández, S.; Gutiérrez, P.A.; Alexandre, E.; Hervás-Martínez, C.; Salcedo-Sanz, S. A review of classification problems and algorithms in renewable energy applications. Energies 2016, 9, 607. [Google Scholar] [CrossRef]
- Bermejo, J.F.; Fernández, J.F.; Polo, F.O.; Márquez, A.C. A review of the use of artificial neural network models for energy and reliability prediction. A study of the solar PV, hydraulic and wind energy sources. Appl. Sci. 2019, 9, 1844. [Google Scholar] [CrossRef]
- Khare, V.; Nema, S.; Baredar, P. Solar-wind hybrid renewable energy system: A review. Renew. Sustain. Energy Rev. 2016, 58, 23–33. [Google Scholar] [CrossRef]
- Ahmed, A.; Khalid, M. A review on the selected applications of forecasting models in renewable power systems. Renew. Sustain. Energy Rev. 2019, 100, 9–21. [Google Scholar] [CrossRef]
- Alkabbani, H.; Ahmadian, A.; Zhu, Q.; Elkamel, A. Machine Learning and Metaheuristic Methods for Renewable Power Forecasting: A Recent Review. Front. Chem. Eng. 2021, 3, 1–21. [Google Scholar] [CrossRef]
- Bracale, A.; Caramia, P.; Carpinelli, G.; Di Fazio, A.R.; Ferruzzi, G. A Bayesian method for Short-Term probabilistic forecasting of photovoltaic generation in smart grid operation and control. Energies 2013, 6, 733–747. [Google Scholar] [CrossRef]
- Ritter, H.; Karaletsos, T. TyXe: Pyro-based Bayesian neural nets for Pytorch. arXiv 2021, arXiv:2110.00276. [Google Scholar] [CrossRef]
- Oreshkin, B.N.; Carpov, D.; Chapados, N.; Bengio, Y. N-Beats: Neural Basis Expansion Analysis for Interpretable Time Series Forecasting. In Proceedings of the 8th International Conference on Learning Representations, ICLR 2020, Addis Ababa, Ethiopia, 26–30 April 2020; pp. 1–21. [Google Scholar]
- Guanoluisa-Pineda, R.; Arcos-Aviles, D.; Flores, M.; Ibarra, A.; Motoasca, E.; Martinez, W.; Guinjoan, F. Photovoltaic power forecast using deep learning techniques with hyperparameters based on Bayesian optimization. Sustainability 2023, 15, 12151. [Google Scholar] [CrossRef]
KPI | Description |
---|---|
Discharge kWh | Energy discharged from the battery |
Charge kWh | Energy used to charge the battery |
Consumption kWh | Energy consumption on the site |
Grid Export kWh | Energy exported to the grid |
PV Charge kWh | Energy from Solar PV system that is diverted instantaneously to BESS Charge |
PV Consumption kWh | Energy from solar PV system that is used instantaneously for Consumption |
PV Export kWh | Energy from Solar PV system that is spilled instantaneously to Grid Export |
Grid Discharge kWh | Energy from BESS Discharge that is spilled instantaneously to Grid Export |
Grid Charge kWh | Energy for BESS Discharge that is supplied instantaneously by Grid Import |
Grid Consumption kWh | Energy from Grid Import that is used instantaneously for Consumption |
Grid Consumption Discharge | Energy from BESS Discharge that is used instantaneously for Consumption |
BESS State of Charge % (SOC) | Percentage that indicates relative storage capacity |
Precipitation probability % | Probability of precipitation |
Precipitation mm | Precipitation measured in mm |
Wind direction | Wind direction in degrees: 0 is north, 90 is east, 270 is west) |
Wind speed | Wind speed measured in knots (one knot is 1.852 km/h) |
Solar radiation | energy density measured in J/cm2 (the energy density is the total energy delivered per unit area in Joules per square centimeter, J/cm2), and sunshine duration measured in minutes (calculated using the geographical coordinates to compute the sun gradient and, hence, the energy density). |
Model | Method | RMSE(1) | MSE(1) | MAE(1) | R2(1) | RMSE(2) | MSE(2) | MAE(2) | R2(2) | ∇ MSE | ∇ MAE |
---|---|---|---|---|---|---|---|---|---|---|---|
Daily Consumption kWh | ESA Winters Multiplicative | 0.909 | 0.826 | 0.672 | 0.944 | 3.173 | 10.071 | 2.173 | 0.301 | 91.80% | 69.08% |
SARIMAX(1,1,1) × (1,0,1) | 0.909 | 0.827 | 0.662 | 0.946 | 3.173 | 10.071 | 2.173 | 0.301 | 91.79% | 69.54% | |
SARIMAX(0,1,1) × (0,1,1) | 0.562 | 0.316 | 0.427 | 0.979 | 3.173 | 10.071 | 2.173 | 0.301 | 96.86% | 80.35% | |
Time Causal Model | 0.890 | 0.792 | 0.653 | 0.95 | 3.173 | 10.071 | 2.173 | 0.301 | 92.14% | 69.95% | |
Daily Production kWh | ESA Winters Multiplicative | 3.775 | 14.250 | 3.041 | 0.59 | 4.632 | 21.451 | 2.983 | 0.374 | 33.57% | −1.94% |
SARIMAX(0,1,1) × (0,0,0) | 0.119 | 0.014 | 0.087 | 1.00 | 4.632 | 21.451 | 2.983 | 0.374 | 99.93% | 97.08% | |
Time Causal Model | 3.790 | 14.364 | 2.911 | 0.630 | 4.632 | 21.451 | 2.983 | 0.374 | 33.04% | 2.41% | |
Daily Headroom % | ESA Winters Multiplicative | 9.929 | 98.592 | 7.734 | 0.724 | N/A | N/A | N/A | N/A | N/A | N/A |
SARIMAX(1,0,4) × (0,1,1) | 2.076 | 4.311 | 1.582 | 0.989 | N/A | N/A | N/A | N/A | N/A | N/A | |
SARIMAX(1,0,4) × (1,1,1) | 9.112 | 83.020 | 6.739 | 0.781 | N/A | N/A | N/A | N/A | N/A | N/A | |
Time Causal Model | 9.310 | 86.676 | 7.031 | 0.780 | N/A | N/A | N/A | N/A | N/A | N/A |
Model | Method | RMSE(1) | MSE(1) | MAE(1) | R2(1) | RMSE(2) | MSE(2) | MAE(2) | R2(2) | ∇ MSE | ∇ MAE |
---|---|---|---|---|---|---|---|---|---|---|---|
Hourly Consumption kWh | ESA Winters Multiplicative | 0.082 | 0.007 | 0.059 | 0.874 | 0.190 | 0.036 | 0.143 | 0.301 | 80.56% | 58.74% |
SARIMAX(0,0,8) × (1,1,1) | 0.000 | 0.000 | 0.000 | 1.000 | 0.190 | 0.036 | 0.143 | 0.301 | 100.00% | 100.00% | |
SARIMAX(9,0,9) × (2,1,2) | 0.065 | 0.004 | 0.049 | 0.918 | 0.190 | 0.036 | 0.143 | 0.301 | 88.89% | 65.73% | |
Time Causal Model | 0.085 | 0.007 | 0.064 | 0.864 | 0.190 | 0.036 | 0.143 | 0.301 | 80.56% | 55.24% | |
Hourly Production kWh | ESA Brown’s Linear Trend | 0.16 | 0.026 | 0.087 | 0.928 | 0.344 | 0.118 | 0.171 | 0.668 | 77.97% | 49.12% |
SARIMAX(3,0,2) × (2,0,0) | 0.032 | 0.001 | 0.017 | 0.997 | 0.344 | 0.118 | 0.171 | 0.668 | 99.15% | 90.06% | |
SARIMAX(3,0,2) × (1,0,0) | 0.036 | 0.001 | 0.019 | 0.996 | 0.344 | 0.118 | 0.171 | 0.668 | 99.15% | 88.89% | |
SARIMAX(24,0,24) × (1,0,0) | 0.033 | 0.001 | 0.017 | 0.997 | 0.344 | 0.118 | 0.171 | 0.668 | 99.15% | 90.06% | |
SARIMAX(3,0,1) × (2,0,0) | 0.033 | 0.001 | 0.017 | 0.997 | 0.344 | 0.118 | 0.171 | 0.668 | 99.15% | 90.06% | |
SARIMAX(3,0,1) × (2,0,2) | 0.033 | 0.001 | 0.017 | 0.997 | 0.344 | 0.118 | 0.171 | 0.668 | 99.15% | 90.06% | |
Time Causal Model | 0.109 | 0.012 | 0.065 | 0.966 | 0.344 | 0.118 | 0.171 | 0.668 | 89.83% | 61.99% | |
Hourly Headroom % | ESA Brown’s Linear Trend | 2.454 | 6.023 | 1.518 | 0.991 | N/A | N/A | N/A | N/A | N/A | N/A |
SARIMAX(0,1,6) × (2,0,1) | 0.358 | 0.128 | 0.247 | 1.000 | N/A | N/A | N/A | N/A | N/A | N/A | |
SARIMAX(5,1,6) × (2,0,1) | 0.757 | 0.573 | 0.520 | 0.999 | N/A | N/A | N/A | N/A | N/A | N/A | |
SARIMAX(0,1,6) × (3,0,2) | 2.150 | 4.623 | 1.449 | 0.993 | N/A | N/A | N/A | N/A | N/A | N/A | |
Time Causal Model | 0.860 | 0.740 | 0.526 | 1.000 | N/A | N/A | N/A | N/A | N/A | N/A |
Model | Method | RMSE(1) | MSE(1) | MAE(1) | R2(1) | RMSE(2) | MSE(2) | MAE(2) | R2(2) | ∇ MSE | ∇ MAE |
---|---|---|---|---|---|---|---|---|---|---|---|
Daily Cluster 4 Consumption kWh | ESA Winters Multiplicative | 0.958 | 0.917 | 0.730 | 0.771 | 5.332 | 28.434 | 4.634 | −6.21 | 96.77% | 84.25% |
SARIMAX (1,0,6) × (0,1,1) | 0.498 | 0.248 | 0.394 | 0.942 | 5.332 | 28.434 | 4.634 | −6.21 | 99.13% | 91.50% | |
SARIMAX (2,1,6) × (1,0,1) | 0.871 | 0.758 | 0.663 | 0.822 | 5.332 | 28.434 | 4.634 | −6.21 | 97.33% | 85.69% | |
Time Causal Model | 0.980 | 0.960 | 0.718 | 0.790 | 5.332 | 28.434 | 4.634 | −6.210 | 96.62% | 84.51% | |
Daily Cluster 4 Production kWh | ESA Winters Multiplicative | 3.307 | 10.937 | 2.680 | 0.567 | 3.880 | 15.056 | 2.657 | 0.396 | 27.36% | −0.87% |
SARIMAX(0,1,2) × (0,0,0) | 0.185 | 0.034 | 0.123 | 0.999 | 3.880 | 15.056 | 2.657 | 0.396 | 99.77% | 95.37% | |
SARIMAX(0,1,2) × (1,0,0) | 2.247 | 5.047 | 1.759 | 0.809 | 3.880 | 15.056 | 2.657 | 0.396 | 66.48% | 33.80% | |
Time Causal Model | 3.230 | 10.433 | 2.445 | 0.633 | 3.880 | 15.056 | 2.657 | 0.396 | 30.71% | 7.98% | |
Daily Cluster 4 Headroom % | ESA Winters Multiplicative | 11.053 | 122.179 | 8.677 | 0.701 | N/A | N/A | N/A | N/A | N/A | N/A |
SARIMAX(0,1,1) × (0,0,0) | 4.554 | 20.739 | 3.555 | 0.950 | N/A | N/A | N/A | N/A | N/A | N/A | |
SARIMAX(0,1,1) × (1,0,1) | 7.685 | 59.06 | 6.606 | 0.860 | N/A | N/A | N/A | N/A | N/A | N/A | |
Time Causal Model | 9.971 | 99.416 | 7.882 | 0.759 | N/A | N/A | N/A | N/A | N/A | N/A |
Model | Method | RMSE(1) | MSE(1) | MAE(1) | R2(1) | RMSE(2) | MSE(2) | MAE(2) | R2(2) | ∇ MSE | ∇ MAE |
---|---|---|---|---|---|---|---|---|---|---|---|
Hourly Cluster 1 Consumption kWh | ESA Winters Multiplicative | 0.095 | 0.009 | 0.070 | 0.870 | 0.213 | 0.045 | 0.163 | 0.339 | 80.00% | 57.06% |
SARIMAX(0,0,2) × (1,1,1) | 0.000 | 0.000 | 0.000 | 1.000 | 0.213 | 0.045 | 0.163 | 0.339 | 100.00% | 100.00% | |
SARIMAX(9,0,9) × (2,1,2) | 0.065 | 0.004 | 0.049 | 0.918 | 0.213 | 0.045 | 0.163 | 0.339 | 91.11% | 69.94% | |
Time Causal Model | 0.098 | 0.010 | 0.075 | 0.864 | 0.213 | 0.045 | 0.163 | 0.339 | 77.78% | 53.99% | |
Hourly Cluster 1 Production kWh | ESA Brown’s Linear Trend | 0.163 | 0.027 | 0.090 | 0.927 | 0.350 | 0.123 | 0.174 | 0.667 | 78.05% | 48.28% |
SARIMAX(3,0,4) × (2,0,1) | 0.027 | 0.001 | 0.015 | 0.998 | 0.350 | 0.123 | 0.174 | 0.667 | 99.19% | 91.38% | |
SARIMAX(22,0,22) × (2,0,1) | 0.027 | 0.001 | 0.015 | 0.998 | 0.350 | 0.123 | 0.174 | 0.667 | 99.19% | 91.38% | |
SARIMAX(8,0,8) × (2,0,1) | 0.032 | 0.001 | 0.017 | 0.997 | 0.350 | 0.123 | 0.174 | 0.667 | 99.19% | 90.23% | |
Time Causal Model | 0.112 | 0.013 | 0.067 | 0.966 | 0.350 | 0.123 | 0.174 | 0.667 | 89.43% | 61.49% | |
Hourly Cluster 1 Headroom % | ESA Brown’s Linear Trend | 2.529 | 6.397 | 1.559 | 0.991 | N/A | N/A | N/A | N/A | N/A | N/A |
SARIMAX(2,1,2) × (2,0,1) | 0.426 | 0.182 | 0.291 | 0.994 | N/A | N/A | N/A | N/A | N/A | N/A | |
SARIMAX(3,1,3) × (2,0,2) | 0.835 | 0.698 | 0.580 | 0.999 | N/A | N/A | N/A | N/A | N/A | N/A | |
Time Causal Model | 0.932 | 0.868 | 0.579 | 0.999 | N/A | N/A | N/A | N/A | N/A | N/A |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Rozas, W.; Pastor-Vargas, R.; García-Vico, A.M.; Carpio, J. Consumption–Production Profile Categorization in Energy Communities. Energies 2023, 16, 6996. https://doi.org/10.3390/en16196996
Rozas W, Pastor-Vargas R, García-Vico AM, Carpio J. Consumption–Production Profile Categorization in Energy Communities. Energies. 2023; 16(19):6996. https://doi.org/10.3390/en16196996
Chicago/Turabian StyleRozas, Wolfram, Rafael Pastor-Vargas, Angel Miguel García-Vico, and José Carpio. 2023. "Consumption–Production Profile Categorization in Energy Communities" Energies 16, no. 19: 6996. https://doi.org/10.3390/en16196996
APA StyleRozas, W., Pastor-Vargas, R., García-Vico, A. M., & Carpio, J. (2023). Consumption–Production Profile Categorization in Energy Communities. Energies, 16(19), 6996. https://doi.org/10.3390/en16196996