$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Collecte et description des ensembles de données
L’ensemble de données (Figure 1) utilisé dans cette recherche comprend des archives historiques cruciales pour la prévision de l’énergie solaire. L’ensemble de données comprend des données opérationnelles quotidiennes provenant d’une centrale solaire thermique concentrée de 50 MW exploitée par Megha Engineering and Infrastructures Limited (MEIL), située près d’Anantapur, Andhra Pradesh, Inde, utilisant la technologie de concentration solaire (CSP) à auge parabolique qui capture l’irradiance normale directe (DNI) et transfère la chaleur via un fluide de transfert de chaleur (HTF) pour produire de l’électricité. Le jeu de données a été collecté du 1er janvier 2015 au 3 octobre 2025 et contient sept attributs clés qui capturent l’information temporelle, les mesures d’irradiance solaire et la production d’énergie. Les attributs temporels incluent « Date », qui indique la date du calendrier au format standard, « Year » indiquant l’année de collecte des données, « Month » représentant le numéro du mois, « Day » désignant le jour du mois, et « Julian Day » offrant un système de numérotation séquentielle tout au long de l’année pour une analyse temporelle continue. La principale variable météorologique d’entrée est la « SUME DNI » mesurée en kWh/m², qui représente l’irradiance normale directe totale (DNI), l’énergie solaire cumulative reçue par mètre carré de surface du collecteur, servant de facteur critique influençant l’efficacité de conversion thermique de la centrale CSP. La variable cible « Production réelle », mesurée en kWh, capture la puissance électrique produite par la centrale CSP, reflétant le résultat du processus de conversion d’énergie solaire-thermique-électrique. Ces attributs permettent collectivement une analyse complète des performances de l’usine, incluant la détermination de l’efficacité thermique, la modélisation DNI-puissance de conversion, l’identification des influences atmosphériques et nuageuses via l’analyse spectrale, la mise en œuvre de protocoles de correction zéro erreur pour l’étalonnage des capteurs, ainsi que le développement de modèles avancés hybrides de prévision par apprentissage automatique pour optimiser la planification opérationnelle en temps réel et améliorer l’efficacité et la fiabilité globales de l’usine CSP. Détails de l’usine disponibles à : https://solarpaces.nrel.gov/project/megha-solar-plant

Figure 1 : Les cinq premières lignes du jeu de données. Données d’exemple montrant les entrées initiales de l’ensemble de données de production d’énergie solaire, affichant les caractéristiques d’entrée et les variables cibles utilisées pour l’entraînement et l’évaluation du modèle. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Préparation des données
L’étude utilise des séries temporelles de production solaire allant du 1er janvier 2015 au 10 mars 2025. Pour tenir compte des problèmes potentiels de qualité des données dans les premières années et se concentrer sur des tendances plus récentes, les enregistrements ont été filtrés à partir du 1er janvier 2017. Les colonnes temporelles (date, année, jour) ont été supprimées sur la base d’analyses préliminaires de corrélation montrant une valeur prédictive négligeable. Les valeurs manquantes ont été imputées à l’aide d’une technique de moyenne mobile afin de maintenir la continuité temporelle tout en minimisant la distorsion des motifs sous-jacents. Trois caractéristiques de latence ont été créées en utilisant la variable cible (Génération réelle (kW/h)) pour capturer les dépendances temporelles.
Répartition des jeux de données
Pour établir des cohortes d’entraînement, de validation et de test équilibrées et représentatives, l’ensemble de données pré-traité a été segmenté à l’aide d’une méthode d’échantillonnage stratifiée. Cette approche garantissait que 70 % des données (2091 données) étaient allouées à l’entraînement, tandis que les ensembles de validation et de test représentaient chacun 15 % (448 données par ensemble).
Normalisation des données
Les fonctionnalités ont été standardisées via StandardScaler, tandis que les valeurs cibles ont été normalisées via MinMaxScaler à la plage [0,1] pour la stabilité du réseau de neurones.
Entraînement des modèles
Cinq modèles hybrides (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU et CNN-Transformer) ont été mis en œuvre pour la prévision de la production solaire. Les données d’entrée ont été restructurées dans un format séquentiel remodelé en (échantillons, pas de temps, caractéristiques) où les pas de temps = 1 pour la plupart des modèles, sauf CNN-LSTM, qui utilisait une fenêtre glissante de 15 étapes. Les ensembles d’entraînement, de validation et de test ont été mis à l’échelle tout en préservant l’ordre temporel pour éviter la fuite de données. Tous les modèles sont entraînés avec 32 lots et 30 époques.
Forêt aléatoire-BiLSTM (Figure 2) : Le modèle hybride proposé combine un réseau bidirectionnel de mémoire à court terme long terme (BiLSTM) avec un régresseur de forêt aléatoire (RF) pour améliorer la précision des prédictions. Tout d’abord, le modèle BiLSTM est entraîné sur les données de série temporelle d’entrée pour capturer les schémas temporels et générer des prédictions initiales. Ensuite, les erreurs résiduelles (différences entre les valeurs réelles et prédites) du BiLSTM sont calculées. Un modèle de forêt aléatoire est ensuite entraîné sur les caractéristiques d’entrée originales pour apprendre et prédire ces résidus. Pour améliorer les performances du modèle RF, les six caractéristiques les plus importantes sont sélectionnées en fonction des scores d’importance des caractéristiques. Enfin, la prédiction corrigée est obtenue en ajoutant les résidus prédits par RF aux sorties BiLSTM. Cette approche hybride exploite la capacité de modélisation de séquences de BiLSTM et la force d’apprentissage d’ensemble de Random Forest pour obtenir une meilleure généralisation et une meilleure performance prédictive.
Soit
la séquence d’entrée au pas de temps t.
Prédiction BiLSTM :
,
Calcul résiduel :
Apprentissage résiduel en utilisant la forêt aléatoire : Soit Z⊂X les k principales caractéristiques sélectionnées selon l’importance des caractéristiques.
Prédiction finale :

Figure 2 : Architecture du modèle de mémoire à court terme bidirectionnel aléatoire foreste. Schéma illustrant l’architecture hybride RF-BiLSTM, montrant l’intégration du traitement des caractéristiques en forêt aléatoire avec des couches LSTM bidirectionnelles pour l’apprentissage temporel des séquences. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
CNN-LSTM (Figure 3) : Le modèle hybride CNN-LSTM commence par traiter la séquence d’entrée à l’aide d’une couche convolutionnelle 1D pour extraire les caractéristiques spatiales locales, suivie d’une activation LeakyReLU, d’une normalisation batch et d’un max pooling. Les caractéristiques extraites sont ensuite passées à travers une pile de trois couches LSTM pour apprendre les dépendances temporelles, avec une normalisation des couches et un abandon appliqués après les deux premiers LSTM pour régularisation. La sortie finale du LSTM est passée à travers des couches denses entièrement connectées avec activation et démission, puis finalement mappée à la sortie à l’aide d’un seul neurone.
Soit
la séquence d’entrée, où T est la fenêtre temporelle et F le nombre de caractéristiques.
Opération CNN :
Pooling maximal :
Cellule LSTM :




Sortie :

Figure 3 : Architecture du modèle CNN-LSTM. Représentation structurelle du modèle hybride réseau neuronal convolutionnel à mémoire courte durée long, démontrant l’extraction convolutionnelle de caractéristiques suivie d’un traitement unidirectionnel de séquences temporelles. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
CNN-BiLSTM (Figure 4) : Le modèle hybride CNN-BiLSTM extrait d’abord les caractéristiques spatiales à l’aide d’une couche convolutionnelle 1D avec 32 filtres, suivie d’une normalisation par lots et d’un max pooling pour réduire la dimensionnalité. La sortie est ensuite transmise à travers une pile de deux couches LSTM bidirectionnelles pour capturer les dépendances temporelles à long terme dans les directions avant et arrière. La régularisation est appliquée via le dropout et la normalisation par lots. Une couche dense avec 128 neurones affine la représentation apprise avant que la couche finale ne l’associe à une seule valeur prédite.
Opération CNN :
Pooling maximal : 
LSTM bidirectionnel :
,
Sortie :

Figure 4 : Architecture du modèle CNN-BiLSTM. Diagramme architectural du modèle Convolutional Neural Network-Bidirectional Long Short-Term Memory, mettant en lumière la combinaison des couches convolutionnelles avec un traitement récurrent bidirectionnel pour une capture accrue des dépendances temporelles. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
CNN-GRU (Figure 5) : Le modèle hybride CNN-GRU commence avec une couche Conv1D utilisant une taille de noyau de 1 pour extraire les caractéristiques spatiales de la seule étape de temps. Le pooling max réduit les dimensions spatiales. Cela est suivi d’une pile de couches GRU : la première retourne des séquences pour capturer les dépendances temporelles, et la seconde résume la séquence en une représentation compacte. Une couche dense finale produit la valeur prédite. La régularisation du dropout est appliquée entre les couches GRU pour éviter le surajustement.
Opération CNN :
Pooling maximal :
Cellule GRU :



Sortie :

Figure 5 : Architecture du modèle CNN-GRU. Schéma du modèle hybride Convolutional Neural Ret-Gated Recurrent Unit, montrant le prétraitement convolutionnel intégré aux couches GRU pour une modélisation temporelle efficace. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
CNN-transformateur (Figure 6) Le modèle CNN-transformateur commence par une couche Conv1D pour extraire les caractéristiques locales de la séquence d’entrée, suivie d’une couche de max pooling. Ces caractéristiques sont transmises à travers un bloc encodeur Transformer composé d’un mécanisme d’auto-attention multi-têtes, de la normalisation des couches et d’un réseau dense en avance. Un pooling global de la moyenne est alors appliqué avant qu’une couche dense finale ne produise la prédiction. Cette architecture est conçue pour capturer à la fois des motifs spatiaux (via CNN) et des dépendances globales (via l’attention Transformer).
Opération CNN :
Attention personnelle à plusieurs têtes :

Où : Q, K, V = XWQ, XWK, XWV et dk est la dimension des clés.
Réseau Feed Forward :

Couches d’ajout et de normalisation :


Sortie :

Figure 6 : Architecture du modèle CNN-Transformer. Aperçu structurel du modèle hybride réseau neuronal convolutionnel-transformer, incorporant l’extraction convolutionnelle des caractéristiques avec des mécanismes d’attention multi-têtes pour une reconnaissance avancée des motifs temporels. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Développement de modèles d’ensemble
Pour améliorer la précision des prévisions et la robustesse des modèles, nous avons mis en place une approche d’ensemble moyenne pondérée utilisant des prédictions issues des cinq modèles hybrides d’apprentissage profond : RF-BiLSTM, CNN-LSTM, CNN-BILSTM, CNN-GRU et CNN-Transformer. L’ensemble a été construit en attribuant des poids optimisés aux prédictions de chaque modèle, avec des poids plus élevés attribués aux modèles démontrant une performance individuelle supérieure, mesurée par leurs scores R². Cette stratégie de pondération garantit que des modèles plus précis contribuent de manière significative à la prévision finale tout en tirant parti des forces complémentaires de toutes les architectures. La sortie de l’ensemble a ensuite été évaluée à l’aide de métriques de performance standard : R², erreur absolue moyenne (MAE) et erreur quadratique moyenne (MSE) afin d’évaluer sa précision prédictive, sa cohérence et sa capacité de généralisation. Cet ensemble d’apprentissage profond vise à intégrer l’extraction de caractéristiques temporelles sous plusieurs perspectives, obtenant ainsi une plus grande précision et robustesse que tout modèle hybride isolé.
Formulation mathématique de la technique d’ensemble :
Soit
représentent l’ensemble des modèles de base correspondant à CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU et CNN-Transformer.
Chaque modèle de base Mi produit une prédiction :
La matrice méta-caractéristique pour l’empilement est formée comme suit :
Le méta-apprenant de la régression de crête estime la prédiction finale ainsi :
où :
-- wi suis les poids d’empilage appris
-- w₀ est le terme de biais
Pour éviter le surapprentissage, la régression de crête minimise la fonction de perte régularisée suivante :

où :
-- yj = cible vraie pour le jième échantillon
-- N = nombre total d’échantillons
-- α = paramètre de régularisation contrôlant le retrait du poids
La prédiction de l’ensemble s’obtient ainsi :
où les poids wi sont automatiquement appris en minimisant la fonction de perte de crête.