$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Raccolta e descrizione del dataset
Il dataset (Figura 1) utilizzato in questa ricerca comprende i documenti storici fondamentali per la previsione dell'energia solare. Il dataset comprende dati operativi giornalieri provenienti da un impianto solare termico concentrato da 50 MW gestito da Megha Engineering and Infrastructures Limited (MEIL), situato vicino ad Anantapur, Andhra Pradesh, India, che utilizza la tecnologia a soffitto parabolico Concentrated Solar Power (CSP) che cattura l'Irradianza Normale Diretta (DNI) e trasferisce calore tramite un Fluido di Trasferimento di Calore (HTF) per generare elettricità. Il dataset è stato raccolto dal 1° gennaio 2015 al 3 ottobre 2025 e contiene sette attributi chiave che catturano informazioni temporali, misurazioni di irradianza solare e output di generazione di energia. Gli attributi temporali includono 'Data', che fornisce la data del calendario in formato standard, 'Anno' che indica l'anno di raccolta dei dati, 'Mese' che rappresenta il numero del mese, 'Giorno' indica il giorno del mese, e 'Giorno Giuliano', che offre un sistema di numerazione sequenziale dei giorni durante tutto l'anno per un'analisi temporale continua. La principale variabile meteorologica di input è la 'DNI SUM' misurata in kWh/m², che rappresenta l'Irradianza Normale Diretta totale (DNI), l'energia solare cumulativa ricevuta per metro quadrato della superficie del collettore, che funge da fattore critico che influenza l'efficienza di conversione termica dell'impianto CSP. La variabile target 'Generazione Effettiva', misurata in kWh, cattura la potenza elettrica prodotta dall'impianto CSP, riflettendo il risultato del processo di conversione dell'energia solare a termica a elettrica. Questi attributi complessivamente permettono un'analisi completa delle prestazioni dell'impianto, inclusa la determinazione dell'efficienza termica, la modellazione di conversione DNI-potenza, l'identificazione delle influenze atmosferiche e della copertura nuvolosa tramite analisi spettrale, l'implementazione di protocolli di correzione a zero errori per la calibrazione dei sensori e lo sviluppo di modelli avanzati ibridi di previsione del machine learning per ottimizzare la pianificazione operativa in tempo reale e migliorare l'efficienza e l'affidabilità complessiva dell'impianto CSP. Dettagli sull'impianto disponibili su: https://solarpaces.nrel.gov/project/megha-solar-plant

Figura 1: Le prime cinque righe del dataset. Dati di esempio che mostrano le voci iniziali del dataset di generazione di energia solare, mostrando le caratteristiche di input e le variabili target utilizzate per l'addestramento e la valutazione del modello. Clicca qui per visualizzare una versione più grande di questa figura.
Preparazione dei dati
Lo studio utilizza dati di serie temporali sulla generazione solare che coprono dal 1° gennaio 2015 al 10 marzo 2025. Per tenere conto di potenziali problemi di qualità dei dati nei primi anni e concentrarsi su modelli più recenti, i record sono stati filtrati dal 1° gennaio 2017 in poi. Le colonne temporali (Data, Anno, Giorno) sono state rimosse sulla base di un'analisi preliminare di correlazione che ha mostrato un valore predittivo trascurabile. I valori mancanti sono stati imputati utilizzando una tecnica di media mobile per mantenere la continuità temporale minimizzando la distorsione dei modelli sottostanti. Tre caratteristiche di lag sono state create utilizzando la variabile target (Actual Generation (kW/h)) per catturare dipendenze temporali.
Divisione dei dataset
Per stabilire coorti di addestramento, validazione e test bilanciate e rappresentative, il dataset pre-elaborato è stato segmentato utilizzando un metodo di campionamento stratificato. Questo approccio ha garantito che il 70% dei dati (2091 dati) fosse allocato per l'addestramento, mentre sia i set di validazione che quelli di test rappresentavano ciascuno il 15% (448 dati per insieme).
Normalizzazione dei dati
Le funzionalità sono state standardizzate utilizzando StandardScaler, mentre i valori target sono stati normalizzati tramite MinMaxScaler nell'intervallo [0,1] per la stabilità della rete neurale.
Addestramento al modello
Cinque modelli ibridi (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU e CNN-Transformer) sono stati implementati per la previsione della generazione solare. I dati di input sono stati ristrutturati in un formato sequenziale riformulato come (campioni, passi temporali, caratteristiche) dove passi temporali = 1 per la maggior parte dei modelli, eccetto CNN-LSTM, che utilizzava una finestra scorrevole di 15 passi. Addestramento, validazione e set di test sono stati scalati mantenendo l'ordine temporale per evitare perdite di dati. Tutti i modelli sono addestrati con 32 unità di lotto e 30 epoche.
Foresta Casuale-BiLSTM (Figura 2): Il modello ibrido proposto combina una rete Bidirezionale Long Short-Term Memory (BiLSTM) con un regressore Random Forest (RF) per migliorare l'accuratezza della previsione. Innanzitutto, il modello BiLSTM viene addestrato sui dati delle serie temporali in ingresso per catturare schemi temporali e generare previsioni iniziali. Successivamente, vengono calcolati gli errori residui (differenze tra valori effettivi e previsti) derivanti dalla BiLSTM. Un modello di Foresta Casuale viene quindi addestrato sulle caratteristiche di input originali per apprendere e prevedere questi residui. Per migliorare le prestazioni del modello RF, le sei caratteristiche più importanti sono selezionate in base ai punteggi di importanza delle caratteristiche. Infine, la predizione corretta si ottiene aggiungendo i residui previsti per RF agli output del BiLSTM. Questo approccio ibrido sfrutta la capacità di modellazione delle sequenze di BiLSTM e la forza di apprendimento in ensemble di Random Forest per ottenere una migliore generalizzazione e performance predittiva.
Sia
la sequenza di input al passo temporale t.
Previsione BiLSTM:
,
Calcolo residuo:
Apprendimento residuo usando Random Forest: Sia Z⊂X le prime k caratteristiche selezionate usando l'importanza delle caratteristiche.
Previsione finale:

Figura 2: Architettura del modello di memoria a breve termine bidirezionale a lungo termine in foresta casuale. Diagramma schematico che illustra l'architettura ibrida RF-BiLSTM, mostrando l'integrazione dell'elaborazione delle caratteristiche Random Forest con strati LSTM bidirezionali per l'apprendimento temporale delle sequenze. Clicca qui per visualizzare una versione più grande di questa figura.
CNN-LSTM (Figura 3): Il modello ibrido CNN-LSTM inizia elaborando la sequenza di input utilizzando uno strato convoluzionale 1D per estrarre le caratteristiche spaziali locali, seguito da un'attivazione LeakyReLU, normalizzazione batch e max pooling. Le caratteristiche estratte vengono poi passate attraverso uno stack di tre strati LSTM per apprendere le dipendenze temporali, con normalizzazione dei livelli e dropout applicati dopo i primi due LSTM per la regolarizzazione. L'output finale dell'LSTM viene passato attraverso strati densi completamente connessi con attivazione e dropout e infine mappato all'uscita usando un singolo neurone.
Sia
la sequenza di input, dove T è la finestra temporale e F è il numero di caratteristiche.
Operazione CNN:
Pooling massimo:
Cellula LSTM:




Output:

Figura 3: Architettura del modello CNN-LSTM. Rappresentazione strutturale del modello ibrido Convolutional Neural Network-Long Short-Term Memory, che dimostra l'estrazione convoluzionale delle caratteristiche seguita dall'elaborazione della sequenza temporale unidirezionale. Clicca qui per visualizzare una versione più grande di questa figura.
CNN-BiLSTM (Figura 4): Il modello ibrido CNN-BiLSTM estrae inizialmente le caratteristiche spaziali utilizzando uno strato convoluzionale 1D con 32 filtri, seguito da normalizzazione batch e max pooling per ridurre la dimensionalità. L'output viene poi passato attraverso uno stack di due strati LSTM bidirezionali per catturare dipendenze temporali a lungo termine sia in direzione avanti che indietro. La regolarizzazione viene applicata tramite dropout e normalizzazione batch. Uno strato denso con 128 neuroni affina la rappresentazione appresa prima che lo strato finale di output la mappi a un singolo valore previsto.
Operazione CNN:
Pooling massimo: 
LSTM bidirezionale:
,
Output:

Figura 4: Architettura del modello CNN-BiLSTM. Diagramma architettonico del modello Convolutional Neural Network-Bidirectional Long Short-Term Memory, che evidenzia la combinazione di strati convoluzionali con l'elaborazione ricorrente bidirezionale per una migliore cattura delle dipendenze temporali. Clicca qui per visualizzare una versione più grande di questa figura.
CNN-GRU (Figura 5): Il modello ibrido CNN-GRU parte da uno strato Conv1D che utilizza un kernel dimensione 1 per estrarre le caratteristiche spaziali dal singolo passo temporale. Il max pooling riduce le dimensioni spaziali. Segue una pila di strati GRU: il primo restituisce sequenze per catturare dipendenze temporali, e il secondo riassume la sequenza in una rappresentazione compatta. Uno strato denso finale produce il valore previsto. La regolarizzazione dropout viene applicata tra gli strati GRU per prevenire il sovrafitting.
Operazione CNN:
Pooling massimo:
Cellula GRU:



Output:

Figura 5: Architettura del modello CNN-GRU. Schema del modello ibrido Convolutional Neural Network-Gated Recurrent Unit, che mostra preprocessing convoluzionale integrato con strati GRU per una modellazione temporale efficiente. Clicca qui per visualizzare una versione più grande di questa figura.
CNN-transformer (Figura 6) Il modello CNN-transformer parte da uno strato Conv1D per estrarre le caratteristiche locali dalla sequenza di input, seguito da uno strato di max pooling. Queste caratteristiche vengono passate attraverso un blocco encoder Transformer costituito da un meccanismo di auto-attenzione multi-testa, normalizzazione degli strati e una rete densa feed-forward. Viene quindi applicato il pooling della media globale prima che uno strato denso finale produca la previsione. Questa architettura è progettata per catturare sia i pattern spaziali (tramite CNN) sia le dipendenze globali (tramite l'attenzione dei Transformer).
Operazione CNN:
Autoattenzione con più teste:

Dove: Q, K, V = XWQ, XWK,XW V e dk è la dimensione delle chiavi.
Rete Feed Forward:

Aggiungi & Norma i livelli:


Output:

Figura 6: Architettura del modello CNN-Transformer. Panoramica strutturale del modello ibrido Convolutional Neural Network-Transformer, che incorpora l'estrazione convoluzionale delle caratteristiche con meccanismi di attenzione multi-head per il riconoscimento avanzato dei pattern temporali. Clicca qui per visualizzare una versione più grande di questa figura.
Sviluppo del modello di ensemble
Per migliorare l'accuratezza delle previsioni e la robustezza del modello, abbiamo implementato un approccio a ensemble a media ponderata utilizzando previsioni dai cinque modelli ibridi di deep learning: RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU e CNN-Transformer. L'ensemble è stato costruito assegnando pesi ottimizzati alle previsioni di ciascun modello, con pesi più alti assegnati ai modelli che dimostrano prestazioni individuali superiori, misurate dai loro punteggi R². Questa strategia di ponderazione garantisce che modelli più accurati contribuiscano in modo più significativo alla previsione finale, pur sfruttando i punti di forza complementari di tutte le architetture. L'output dell'ensemble è stato poi valutato utilizzando metriche standard di prestazione: R², errore assoluto medio (MAE) ed errore quadratico medio (MSE) per valutarne l'accuratezza predittiva, la coerenza e la capacità di generalizzazione. Questo insieme di deep learning mira a integrare l'estrazione temporale delle caratteristiche da molteplici prospettive, raggiungendo così maggiore accuratezza e robustezza rispetto a qualsiasi singolo modello ibrido isolato.
Formulazione matematica della tecnica dell'ensemble:
Sia
rappresentano l'insieme dei modelli base corrispondenti a CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU e CNN-Transformer.
Ogni modello base Mi produce una previsione:
La matrice meta-feature per l'impilamento si forma come:
Il meta-apprendista della Ridge Regression stima la previsione finale come:
dove:
-- wi sono i i appresi a sovrapporre i pesi
-- w₀ è il termine di pregiudizio
Per evitare l'overfitting, la Regressione Ridge minimizza la seguente funzione di perdita regolarizzata:

dove:
-- yj = vero target per il jesimo campione
-- N = numero totale di campioni
-- α = parametro di regolarizzazione che controlla la riduzione del peso
La previsione dell'insieme si ottiene come:
dove i pesi wi vengono automaticamente appresi minimizzando la funzione di perdita di cresta.