Table of Contents
Fetching ...

Symmetry-Constrained Forecasting of Periodically Correlated Energy Processes

Cyril Voyant, Candice Banes, Luis Garcia-Gutierrez, Gilles Notton, Milan Despotovic, Zaher Mundher Yaseen

Abstract

Time series in energy systems, such as solar irradiance, wind speed, or electrical load, are characterized by strong diurnal and seasonal periodicities. Accurate forecasting requires accounting for time varying statistical properties that stationary or classical persistence models cannot capture. \cyr{A family of analytical forecasting operators for cyclostationary processes is introduced, extending persistence through a closed form coefficient $\tildeλ(t,τ)=\tfrac{1}{2}\bigl(1+ρ(t,τ)\bigr)$, where $ρ(t,τ)$ denotes the local correlation between the current observation and its phase aligned counterpart}. This formulation preserves periodic variance and covariance, achieving a symmetry induced reduction of effective degrees of freedom. The resulting operator defines a training free analytical limit of persistence under periodic non stationarity. Validation on synthetic cyclostationary signals and empirical renewable energy datasets demonstrates consistent accuracy gains over classical persistence, particularly at multi hour horizons. By embedding temporal symmetry into the prediction process, the framework provides a physically interpretable, reproducible, and computationally minimal baseline for forecasting periodic processes across energy and complex systems.

Symmetry-Constrained Forecasting of Periodically Correlated Energy Processes

Abstract

Time series in energy systems, such as solar irradiance, wind speed, or electrical load, are characterized by strong diurnal and seasonal periodicities. Accurate forecasting requires accounting for time varying statistical properties that stationary or classical persistence models cannot capture. \cyr{A family of analytical forecasting operators for cyclostationary processes is introduced, extending persistence through a closed form coefficient , where denotes the local correlation between the current observation and its phase aligned counterpart}. This formulation preserves periodic variance and covariance, achieving a symmetry induced reduction of effective degrees of freedom. The resulting operator defines a training free analytical limit of persistence under periodic non stationarity. Validation on synthetic cyclostationary signals and empirical renewable energy datasets demonstrates consistent accuracy gains over classical persistence, particularly at multi hour horizons. By embedding temporal symmetry into the prediction process, the framework provides a physically interpretable, reproducible, and computationally minimal baseline for forecasting periodic processes across energy and complex systems.
Paper Structure (15 sections, 11 theorems, 51 equations, 8 figures, 3 tables)

This paper contains 15 sections, 11 theorems, 51 equations, 8 figures, 3 tables.

Key Result

Theorem 1

Let $I(t)$ be a random process exhibiting statistical periodicity of period $T$. The process is periodically stationary if its first- and second-order moments satisfy where $C(I(t),I(t+n\Delta t))$ denotes the covariance between $I(t)$ and $I(t+n\Delta t)$. The variance and covariance are periodic functions of the phase $t \bmod T$, and the correlation coefficient is defined by

Figures (8)

  • Figure 1: Flowcharts illustrating the forecasting process for the BLEND Persistence operator in three cases: (a) stationary, (b) cyclostationary, and (c) simplified cyclostationary. The stationary case assumes constant statistical properties, while the cyclostationary case integrates periodic variations in mean, variance, and correlations. The simplified cyclostationary case uses an empirical approximation to further reduce complexity while retaining predictive efficiency.
  • Figure 2: Geographical distribution of the 68 $\mathtt{SIAR}$ weather stations and kriging surface of daytime averaged $I$ (2017-2020)
  • Figure 3: Cyclic parameters (mean $\mu(t)$, standard deviation $\sigma(t)$, and correlation coefficient $\rho(t)$) shown on the left and prediction model comparisons (right) are applied to a noisy periodic signal ($T=24$, $n=1$ and $\Delta t=1$). The left side illustrates the smoothed and original cyclic statistics, while the right side displays the performance of models.
  • Figure 4: Normalized Root Mean Square Error (nRMSE) plotted against variability metrics for synthetic time series across multiple prediction horizons. Each row represents a specific horizon (1, 3, 6, and 9 hours), and each column corresponds to a variability metric: Coefficient of Variation (CV), Mean Absolute Return (MAR), Root Mean Square Error (RMSE), and Autocorrelation at Lag 1. Trend lines using third order polynomials highlight the relationship between nRMSE and the variability characteristics of the data, illustrating the impact of variability on model performance.
  • Figure 5: Boxplots of the Normalized Root Mean Square Error (nRMSE) for different models across prediction horizons (1h, 3h, 6h, 9h). The x-axis represents the forecasting models, while the y-axis shows the distribution of nRMSE values across synthetic time series simulations.
  • ...and 3 more figures

Theorems & Definitions (49)

  • Definition 1: Simple Persistence Operator $[\emptyset ; \mathcal{O}(1)]$
  • Remark 1
  • Definition 2: Cyclic Persistence Operator $[\emptyset ; \mathcal{O}(1)]$
  • Remark 2
  • Remark 3
  • Definition 3: Smart Persistence Operator $[\mathcal{O}(\mathtt{ref}); \mathcal{O}(1)]$
  • Remark 4
  • Remark 5
  • Definition 4: $\mathtt{CLIPER}$ Operator under Stationary Assumptions $[\mathcal{O}(n); \mathcal{O}(1)]$
  • Remark 6
  • ...and 39 more