$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Datensatzsammlung und -beschreibung
Der in dieser Forschung verwendete Datensatz (Abbildung 1) umfasst historische Aufzeichnungen, die für die Vorhersage von Solarenergie entscheidend sind. Der Datensatz umfasst tägliche Betriebsdaten eines 50-MW-konzentrierten Solarthermiekraftwerks, das von Megha Engineering and Infrastructures Limited (MEIL) betrieben wird und sich in der Nähe von Anantapur, Andhra Pradesh, Indien, befindet, das Parabolrinn-Konzentrierende Solarenergietechnologie (CSP) nutzt, die Direct Normal Irradiance (DNI) erfasst und Wärme über eine Wärmeübertragungsflüssigkeit (HTF) zur Stromerzeugung überträgt. Der Datensatz wurde vom 1. Januar 2015 bis zum 3. Oktober 2025 gesammelt und enthält sieben Schlüsselattribute, die zeitliche Informationen, Messungen der Sonneneinstrahlung und die Stromerzeugungsleistung erfassen. Zu den zeitlichen Attributen gehören 'Datum', das das Kalenderdatum im Standardformat angibt, 'Jahr', das das Jahr der Datenerhebung angibt, 'Monat' für die Monatsnummer, 'Tag', das den Tag des Monats angibt, und 'Julianischer Tag', der ein fortlaufendes Tagesnummerierungssystem für kontinuierliche zeitliche Analyse über das ganze Jahr hinweg bietet. Die primäre meteorologische Eingabevariable ist die 'DNI SUM', gemessen in kWh/m², die die gesamte Direct Normal Irradiance (DNI) darstellt, also die kumulierte Sonnenenergie pro Quadratmeter der Kollektorfläche, die als kritischer Faktor für die Wärmeumwandlungseffizienz des CSP-Kraftwerks dient. Die Zielvariable 'Tatsächliche Erzeugung', gemessen in kWh, erfasst die vom CSP-Kraftwerk erzeugte elektrische Leistung und spiegelt das Ergebnis des Umwandlungsprozesses von Solar-Wärme-Strom zu Strom wider. Diese Eigenschaften ermöglichen zusammen eine umfassende Analyse der Anlagenleistung, einschließlich der Bestimmung der thermischen Effizienz, der Modellierung der DNI-zu-Energie-Umwandlung, der Identifizierung atmosphärischer und wolkenbezogener Einflüsse durch Spektralanalyse, der Implementierung von Fehlerfrei-Korrekturprotokollen für die Sensorkalibrierung sowie der Entwicklung fortschrittlicher hybrider Machine-Learning-Prognosemodelle zur Optimierung der Echtzeitbetriebsplanung und zur Verbesserung der Gesamteffizienz und Zuverlässigkeit der CSP-Anlagen. Pflanzendetails sind verfügbar unter: https://solarpaces.nrel.gov/project/megha-solar-plant

Abbildung 1: Die obersten fünf Zeilen des Datensatzes. Beispieldaten, die die ersten Einträge des Solarstromerzeugungsdatensatzes zeigen, Eingabemerkmale und Zielvariablen darstellen, die für Modelltraining und -auswertung verwendet werden. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Datenvorbereitung
Die Studie verwendet Zeitreihendaten zur Solarerzeugung, die vom 01. Januar 2015 bis zum 10. März 2025 reichen. Um mögliche Qualitätsprobleme in den Anfangsjahren zu berücksichtigen und sich auf neuere Muster zu konzentrieren, wurden die Datensätze ab dem 1. Januar 2017 gefiltert. Zeitliche Spalten (Datum, Jahr, Tag) wurden auf Grundlage einer vorläufigen Korrelationsanalyse entfernt, die einen vernachlässigbaren Prädiktionswert zeigte. Fehlende Werte wurden mittels einer gleitenden Durchschnittstechnik imputiert, um zeitliche Kontinuität aufrechtzuerhalten und gleichzeitig Verzerrungen der zugrunde liegenden Muster zu minimieren. Drei Lag-Features wurden mit der Zielvariablen (Actual Generation (kW/h)) erstellt, um zeitliche Abhängigkeiten zu erfassen.
Datensatzaufteilung
Um ausgewogene und repräsentative Trainings-, Validierungs- und Testkohorten zu etablieren, wurde der vorverarbeitete Datensatz mittels einer stratifizierten Stichprobenmethode segmentiert. Dieser Ansatz stellte sicher, dass 70 % der Daten (Daten von 2091) für das Training zugewiesen wurden, während sowohl Validierungs- als auch Testsätze jeweils 15 % ausmachten (448 Daten pro Satz).
Datennormalisierung
Funktionen wurden mit StandardScaler standardisiert, während Zielwerte über MinMaxScaler auf den [0,1]-Bereich für die Stabilität neuronaler Netzwerke normalisiert wurden.
Modelltraining
Fünf Hybridmodelle (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU und CNN-Transformer) wurden für die Prognose der Solarerzeugung implementiert. Die Eingabedaten wurden in ein sequentielles Format umstrukturiert, das als (Beispiele, Zeitschritte, Merkmale) umgestaltet wurde, wobei Zeitschritte = 1 für die meisten Modelle sind, außer CNN-LSTM, das ein gleitendes Fenster von 15 Schritten verwendete. Trainings-, Validierungs- und Testsets wurden skaliert, wobei die zeitliche Ordnung erhalten blieb, um Datenlecks zu vermeiden. Alle Modelle werden mit 32 Chargen und 30 Epochen trainiert.
Random Forest-BiLSTM (Abbildung 2): Das vorgeschlagene Hybridmodell kombiniert ein bidirektionales Long Short-Term Memory-Netzwerk (BiLSTM) mit einem Random Forest (RF)-Regressor, um die Vorhersagegenauigkeit zu verbessern. Erstens wird das BiLSTM-Modell mit den Eingabezeitreihendaten trainiert, um zeitliche Muster zu erfassen und erste Vorhersagen zu generieren. Anschließend werden die Restfehler (Unterschiede zwischen tatsächlichen und vorhergesagten Werten) des BiLSTM berechnet. Ein Random Forest-Modell wird dann auf den ursprünglichen Eingabemerkmalen trainiert, um diese Residuen zu lernen und vorherzusagen. Um die Leistung des RF-Modells zu verbessern, werden die sechs wichtigsten Merkmale anhand von Merkmalsbedeutungswerten ausgewählt. Schließlich wird die korrigierte Vorhersage erhalten, indem die RF-vorhergesagten Residuen zu den BiLSTM-Ausgaben hinzugefügt werden. Dieser hybride Ansatz nutzt die Sequenzmodellierungsfähigkeit von BiLSTM und die Ensemble-Lernstärke von Random Forest, um eine bessere Verallgemeinerung und prädiktive Leistung zu erzielen.
Sei
die Eingabesequenz im Zeitschritt t.
BiLSTM-Vorhersage:
,
Residualberechnung:
Residual Learning mit Random Forest: Lass Z⊂X die Top-K Features sein, die anhand der Merkmals-Importance ausgewählt werden.
Endgültige Prognose:

Abbildung 2: Architektur des zufälligen wald-bidirektionalen Langzeitgedächtnismodells. Schematisches Diagramm, das die RF-BiLSTM-Hybridarchitektur illustriert und die Integration der Random Forest-Feature-Verarbeitung mit bidirektionalen LSTM-Schichten für das temporale Sequenzlernen zeigt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
CNN-LSTM (Abbildung 3): Das CNN-LSTM-Hybridmodell beginnt mit der Verarbeitung der Eingabesequenz mit einer 1D-Faltungsschicht, um lokale räumliche Merkmale zu extrahieren, gefolgt von einer LeakyReLU-Aktivierung, Batch-Normalisierung und maximalem Pooling. Die extrahierten Features werden dann durch einen Stapel von drei LSTM-Schichten geleitet, um zeitliche Abhängigkeiten zu lernen, wobei nach den ersten beiden LSTMs Schichtnormalisierung und Dropout zur Regularisierung angewendet werden. Der endgültige LSTM-Ausgang wird durch vollständig verbundene, dichte Schichten mit Aktivierung und Dropout geleitet und schließlich mit einem einzelnen Neuron auf den Ausgang abgebildet.
Sei
die Eingabesequenz, wobei T das Zeitfenster und F die Anzahl der Merkmale ist.
CNN-Betrieb:
Max-Pooling:
LSTM-Zelle:




Ausgabe:

Abbildung 3: Architektur des CNN-LSTM-Modells. Strukturelle Darstellung des Hybridmodells des konvolutionellen neuronalen Netzwerks-Langzeit-Kurzzeitgedächtnisses, das eine konvolutionelle Merkmalsextraktion gefolgt von unidirektionaler temporaler Sequenzverarbeitung demonstriert. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
CNN-BiLSTM (Abbildung 4): Das CNN-BiLSTM-Hybridmodell extrahiert zunächst räumliche Merkmale mittels einer 1D-Faltungsschicht mit 32 Filtern, gefolgt von Batch-Normalisierung und maximalem Pooling zur Reduzierung der Dimensionalität. Die Ausgabe wird dann durch einen Stapel von zwei bidirektionalen LSTM-Schichten geleitet, um langfristige zeitliche Abhängigkeiten sowohl in Vorwärts- als auch in Rückwärtsrichtung zu erfassen. Die Regularisierung erfolgt über Dropout und Batch-Normalisierung. Eine dichte Schicht mit 128 Neuronen verfeinert die gelernte Darstellung, bevor die endgültige Ausgabeschicht sie auf einen einzigen vorhergesagten Wert abbildet.
CNN-Betrieb:
Max-Pooling: 
Bidirektionales LSTM:
,
Ausgabe:

Abbildung 4: Architektur des CNN-BiLSTM-Modells. Architekturdiagramm des Convolutional Neural Network-Bidirectional Long Short-Term Memory-Modells, das die Kombination von Faltungsschichten mit bidirektionaler rekurrierender Verarbeitung zur verbesserten zeitlichen Abhängigkeitserfassung hervorhebt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
CNN-GRU (Abbildung 5): Das CNN-GRU-Hybridmodell beginnt mit einer Conv1D-Schicht, die eine Kerngröße von 1 verwendet, um räumliche Merkmale aus dem einzelnen Zeitschritt zu extrahieren. Max-Pooling reduziert die räumlichen Dimensionen. Darauf folgt ein Stapel von GRU-Schichten: die erste gibt Sequenzen zurück, um zeitliche Abhängigkeiten zu erfassen, und die zweite fasst die Sequenz zu einer kompakten Darstellung zusammen. Eine endgültige dichte Schicht liefert den vorhergesagten Wert. Zwischen den GRU-Schichten wird eine Dropout-Regularisierung angewendet, um Überanpassungen zu verhindern.
CNN-Betrieb:
Max-Pooling:
GRU-Zelle:



Ausgabe:

Abbildung 5: Architektur des CNN-GRU-Modells. Schaltplan des Convolutional Neural Network-Gated Recurrent Unit Hybrid-Modells, das eine konvolutionelle Vorverarbeitung mit GRU-Schichten für effiziente temporale Modellierung zeigt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
CNN-Transformator (Abbildung 6) Das CNN-Transformator-Modell beginnt mit einer Conv1D-Schicht, um lokale Merkmale aus der Eingabesequenz zu extrahieren, gefolgt von einer maximalen Pooling-Schicht. Diese Funktionen werden durch einen Transformer-Encoder-Block geleitet, der aus einem Mehrkopf-Selbstbeobachtungsmechanismus, Schichtnormalisierung und einem vorwärtsdichten Netzwerk besteht. Anschließend wird globales Durchschnittspooling angewendet, bevor eine letzte dichte Schicht die Vorhersage liefert. Diese Architektur ist darauf ausgelegt, sowohl räumliche Muster (über CNN) als auch globale Abhängigkeiten (über Transformer Attention) zu erfassen.
CNN-Betrieb:
Mehrköpfige Selbstaufmerksamkeit:

Wobei: Q, K, V = XWQ, XWK, XWV und dk die Dimension der Schlüssel ist.
Feedforward-Netzwerk:

Hinzufügen und Norm-Schichten:


Ausgabe:

Abbildung 6: Architektur des CNN-Transformer-Modells. Struktureller Überblick über das Convolutional Neural Network-Transformer-Hybridmodell, das konvolutionale Merkmalextraktion mit Multi-Head-Attention-Mechanismen für fortgeschrittene zeitliche Mustererkennung integriert. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Entwicklung des Ensemblemodells
Um die Prognosegenauigkeit und Modellrobustheit zu verbessern, implementierten wir einen gewichteten Durchschnittsensemble-Ansatz unter Verwendung von Vorhersagen aus den fünf hybriden Deep-Learning-Modellen: RF-BiLSTM, CNN-LSTM, CNN-BILSTM, CNN-GRU und CNN-Transformer. Das Ensemble wurde konstruiert, indem jedem Modell optimierte Gewichtungen zugewiesen wurden, wobei höhere Gewichtungen für Modelle mit überlegener individueller Leistung vergeben wurden, gemessen an ihren R²-Werten. Diese Gewichtungsstrategie stellt sicher, dass genauere Modelle wesentlich zur endgültigen Prognose beitragen und gleichzeitig die komplementären Stärken aller Architekturen nutzen. Das Ensemble-Ergebnis wurde anschließend mit Standardleistungskennzahlen bewertet: R², mittlerer absoluter Fehler (MAE) und mittlerer quadratischer Fehler (), um seine prädiktive Genauigkeit, Konsistenz und Verallgemeinerungsfähigkeit zu bewerten. Dieses Deep-Learning-Ensemble zielt darauf ab, zeitliche Merkmalextraktion aus mehreren Perspektiven zu integrieren und so eine größere Genauigkeit und Robustheit zu erreichen als jedes einzelne hybride Modell isoliert.
Mathematische Formulierung der Ensemble-Technik:
Lass
repräsentieren die Menge der Basismodelle, die CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU und CNN-Transformer entsprechen.
Jedes Basismodell Mi liefert eine Vorhersage:
Die Meta-Feature-Matrix für das Stapeln wird gebildet wie folgt:
Der Ridge-Regression-Meta-Learner schätzt die endgültige Vorhersage wie folgt:
wobei:
-- wi sind die erlernten Stapelgewichte
-- w₀ ist der Verzerrungsterm
Um Überanpassungen zu vermeiden, minimiert die Ridge Regression die folgende regularisierte Verlustfunktion:

wobei:
-- yj = wahres Ziel für diej-te Probe
-- N = Gesamtzahl der Stichproben
-- α = Regularisierungsparameter, der die Gewichtsschrumpfung steuert
Die Ensemble-Vorhersage wird wie folgt erhalten:
wobei die Gewichte wi automatisch durch Minimierung der Ridge-Verlustfunktion ermittelt werden.