  Electricity load forecasting is crucial for effectively managing and
optimizing power grids. Over the past few decades, various statistical and deep
learning approaches have been used to develop load forecasting models. This
paper presents an interpretable machine learning approach that identifies load
dynamics using data-driven methods within an operator-theoretic framework. We
represent the load data using the Koopman operator, which is inherent to the
underlying dynamics. By computing the corresponding eigenfunctions, we
decompose the load dynamics into coherent spatiotemporal patterns that are the
most robust features of the dynamics. Each pattern evolves independently
according to its single frequency, making its predictability based on linear
dynamics. We emphasize that the load dynamics are constructed based on coherent
spatiotemporal patterns that are intrinsic to the dynamics and are capable of
encoding rich dynamical features at multiple time scales. These features are
related to complex interactions over interconnected power grids and different
exogenous effects. To implement the Koopman operator approach more efficiently,
we cluster the load data using a modern kernel-based clustering approach and
identify power stations with similar load patterns, particularly those with
synchronized dynamics. We evaluate our approach using a large-scale dataset
from a renewable electric power system within the continental European
electricity system and show that the Koopman-based approach outperforms a deep
learning (LSTM) architecture in terms of accuracy and computational efficiency.
The code for this paper has been deposited in a GitHub repository, which can be
accessed at the following address github.com/Shakeri-Lab/Power-Grids.
