$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Datasetverzameling en beschrijving
De dataset (Figuur 1) die in dit onderzoek wordt gebruikt, bestaat uit historische gegevens die cruciaal zijn voor het voorspellen van zonne-energie. De dataset bevat dagelijkse operationele gegevens van een 50 MW geconcentreerde zonne-thermische installatie, geëxploiteerd door Megha Engineering and Infrastructures Limited (MEIL), gelegen nabij Anantapur, Andhra Pradesh, India, die gebruikmaakt van parabolische trog Concentrating Solar Power (CSP) technologie die Direct Normal Irradiance (DNI) opvangt en warmte overdraagt via een Heat Transfer Fluid (HTF) om elektriciteit op te wekken. De dataset werd verzameld van 1 januari 2015 tot 3 oktober 2025 en bevat zeven belangrijke attributen die temporele informatie, metingen van zonnestraling en energieproductie vastleggen. De temporele attributen omvatten 'Datum', dat de kalenderdatum in standaardformaat aangeeft, 'Jaar' dat het jaar van gegevensverzameling aangeeft, 'Maand' dat het maandnummer aangeeft, 'Dag' dat de dag van de maand aangeeft, en 'Julian Day' dat een sequentieel dagnummeringssysteem biedt gedurende het hele jaar voor continue temporele analyse. De primaire meteorologische invoervariabele is 'DNI SUM', gemeten in kWh/m², die de totale Direct Normal Irradiance (DNI) vertegenwoordigt, de cumulatieve zonne-energie die per vierkante meter van het collectoroppervlak wordt ontvangen, en dient als de kritieke factor die de thermische conversieefficiëntie van CSP-installatie beïnvloedt. De doelvariabele 'Werkelijke Opwekking', gemeten in kWh, omvat de elektrische stroomproductie van de CSP-centrale en weerspiegelt het resultaat van het energieconversieproces van zon-naar-thermisch naar elektriciteit. Deze eigenschappen maken samen een uitgebreide analyse van de prestaties van de installatie mogelijk, waaronder het bepalen van thermische efficiëntie, DNI-naar-vermogen-conversiemodellering, identificatie van atmosferische en wolkenbedekking door middel van spectrale analyse, implementatie van nul-foutcorrectieprotocollen voor sensorkalibratie, en ontwikkeling van geavanceerde hybride machine learning-voorspellingsmodellen voor het optimaliseren van realtime operationele planning en het verbeteren van de algehele efficiëntie en betrouwbaarheid van CSP-installaties. Plantendetails beschikbaar op: https://solarpaces.nrel.gov/project/megha-solar-plant

Figuur 1: De bovenste vijf rijen van de dataset. Voorbeeldgegevens die de eerste vermeldingen van de zonne-energieopwekkingsdataset tonen, met invoerkenmerken en doelvariabelen die worden gebruikt voor modeltraining en evaluatie. Klik hier om een grotere versie van deze figuur te bekijken.
Datavoorbereiding
De studie maakt gebruik van tijdreeksgegevens uit zonne-energie die zich uitstrekken van 1 januari 2015 tot 10 maart 2025. Om rekening te houden met mogelijke kwaliteitsproblemen in de vroege jaren en zich te richten op recentere patronen, werden de records vanaf 1 januari 2017 gefilterd. Temporele kolommen (Datum, Jaar, Dag) werden verwijderd op basis van voorlopige correlatieanalyse die een verwaarloosbare voorspellende waarde aantoonde. Ontbrekende waarden werden geïmplicteerd met behulp van een voortschrijdend gemiddelde om de temporele continuïteit te behouden en tegelijkertijd vervorming van onderliggende patronen te minimaliseren. Er werden drie vertragingsfuncties gemaakt met behulp van de doelvariabele (Werkelijke Generatie (kW/h)) om temporele afhankelijkheden vast te leggen.
Datasetsplitsing
Om evenwichtige en representatieve trainings-, validatie- en testcohorten op te zetten, werd de vooraf verwerkte dataset gesegmenteerd met behulp van een gestratificeerde steekproefmethode. Deze aanpak zorgde ervoor dat 70% van de data (data uit 2091) werd toegewezen aan training, terwijl zowel validatie- als testsets elk 15% omvatten (448 data per set).
Datanormalisatie
Functies werden gestandaardiseerd met StandardScaler, terwijl doelwaarden via MinMaxScaler werden genormaliseerd naar het [0,1]-bereik voor stabiliteit van neurale netwerken.
Modeltraining
Vijf hybride modellen (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU en CNN-Transformer) werden geïmplementeerd voor de voorspelling van zonne-energie. De invoergegevens werden herstructureerd in een sequentieel formaat dat werd hervormd als (samples, timesteps, features) waarbij timesteps = 1 voor de meeste modellen, behalve CNN-LSTM, dat een schuifvenster van 15 stappen gebruikte. Training, validatie en testsets werden opgeschaald terwijl temporele orde behouden bleef om datalekken te voorkomen. Alle modellen worden getraind met 32 batch-grootte en 30 epochs.
Random Forest-BiLSTM (Figuur 2): Het voorgestelde hybride model combineert een Bidirectional Long Short-Term Memory (BiLSTM) netwerk met een Random Forest (RF) regressor om de voorspellingsnauwkeurigheid te verbeteren. Ten eerste wordt het BiLSTM-model getraind op de invoertijdreeksgegevens om temporele patronen vast te leggen en initiële voorspellingen te genereren. Daarna worden de restfouten (verschillen tussen werkelijke en voorspelde waarden) van de BiLSTM berekend. Een Random Forest-model wordt vervolgens getraind op de oorspronkelijke inputkenmerken om deze residuen te leren en te voorspellen. Om de prestaties van het RF-model te verbeteren, worden de zes belangrijkste kenmerken geselecteerd op basis van feature-importance scores. Ten slotte wordt de gecorrigeerde voorspelling verkregen door de RF-voorspelde residuen toe te voegen aan de BiLSTM-uitgangen. Deze hybride benadering maakt gebruik van het sequence-modelleringsvermogen van BiLSTM en de ensemble-leerkracht van Random Forest om betere generalisatie en voorspellende prestaties te bereiken.
Laat
de invoerreeks zijn op tijdstap t.
BiLSTM-voorspelling:
,
Residuele berekening:
Residueel leren met Random Forest: Laat Z⊂X de top-k features zijn die geselecteerd zijn met feature-belangrijkheid.
Eindvoorspelling:

Figuur 2: Architectuur van het willekeurig bos-bidirectioneel langetermijngeheugenmodel. Schematisch diagram dat de RF-BiLSTM hybride architectuur illustreert, met de integratie van Random Forest-featureverwerking met bidirectionele LSTM-lagen voor temporele sequentieleren. Klik hier om een grotere versie van deze figuur te bekijken.
CNN-LSTM (Figuur 3): Het CNN-LSTM hybride model begint met het verwerken van de invoersequentie met behulp van een 1D Convolutionele laag om lokale ruimtelijke kenmerken te extraheren, gevolgd door een LeakyReLU-activatie, batchnormalisatie en maximale pooling. De geëxtraheerde features worden vervolgens door een stapel van drie LSTM-lagen gestuurd om temporele afhankelijkheden te leren, waarbij laagnormalisatie en dropout worden toegepast na de eerste twee LSTM's voor regularisatie. De uiteindelijke LSTM-uitgang wordt door volledig verbonden, dichte lagen geleid met activatie en uitval en uiteindelijk met behulp van één neuron op de uitgang toegewezen.
Laat
de invoerreeks zijn, waarbij T het tijdsvenster is en F het aantal kenmerken.
CNN-operatie:
Maximale pooling:
LSTM-cel:




Output:

Figuur 3: Architectuur van het CNN-LSTM-model. Structurele representatie van het Convolutioneel Neuraal Netwerk-Lang Kortetermijngeheugen hybride model, waarbij convolutionele feature-extractie gevolgd wordt door unidirectionele temporele sequentieverwerking. Klik hier om een grotere versie van deze figuur te bekijken.
CNN-BiLSTM (Figuur 4): Het CNN-BiLSTM hybride model extraheert eerst ruimtelijke kenmerken met behulp van een 1D convolutionele laag met 32 filters, gevolgd door batchnormalisatie en maximale pooling om de dimensionaliteit te verminderen. De output wordt vervolgens door een stapel van twee bidirectionele LSTM-lagen gestuurd om langdurige temporele afhankelijkheden zowel in voorwaartse als achterwaartse richting vast te leggen. Regularisatie wordt toegepast via dropout en batchnormalisatie. Een dichte laag met 128 neuronen verfijnt de aangeleerde representatie voordat de uiteindelijke uitvoerlaag deze aan één voorspelde waarde koppelt.
CNN-operatie:
Maximale pooling: 
Bidirectionele LSTM:
,
Output:

Figuur 4: Architectuur van het CNN-BiLSTM-model. Architectuurdiagram van het Convolutional Neural Network-Bidirectioneel Long Short-Term Memory-model, dat de combinatie van convolutionele lagen met bidirectionele recidiverende verwerking voor verbeterde temporele afhankelijkheidsvastlegging benadrukt. Klik hier om een grotere versie van deze figuur te bekijken.
CNN-GRU (Figuur 5): Het CNN-GRU hybride model begint met een Conv1D-laag waarbij een kernelgrootte van 1 wordt gebruikt om ruimtelijke kenmerken uit de enkele tijdstap te halen. Max pooling vermindert de ruimtelijke dimensies. Hierop volgt een stapel GRU-lagen: de eerste geeft sequenties terug om temporele afhankelijkheden vast te leggen, en de tweede vat de sequence samen tot een compacte representatie. Een laatste dichte laag geeft de voorspelde waarde uit. Dropout-regularisatie wordt toegepast tussen GRU-lagen om overfitting te voorkomen.
CNN-operatie:
Maximale pooling:
GRU-cel:



Output:

Figuur 5: Architectuur van het CNN-GRU-model. Schema van het Convolutional Neural Network-Gated Recurrent Unit hybride model, met convolutionele preprocessing geïntegreerd met GRU-lagen voor efficiënte temporele modellering. Klik hier om een grotere versie van deze figuur te bekijken.
CNN-transformator (Figuur 6) Het CNN-transformatormodel begint met een Conv1D-laag om lokale kenmerken uit de invoerreeks te halen, gevolgd door een maximale poolinglaag. Deze functies worden doorgestuurd door een Transformer-encoderblok dat bestaat uit een multi-head zelf-aandacht mechanisme, laagnormalisatie en een feed-forward dicht netwerk. Global average pooling wordt vervolgens toegepast voordat een laatste dichte laag de voorspelling geeft. Deze architectuur is ontworpen om zowel ruimtelijke patronen (via CNN) als globale afhankelijkheden (via Transformer attention) vast te leggen.
CNN-operatie:
Zelf-aandacht met meerdere hoofden:

Waarbij: Q, K, V = XWQ, XWK, XWV en dk de dimensie van sleutels is.
Feed Forward-netwerk:

Voeg toe en norm lagen:


Output:

Figuur 6: Architectuur van het CNN-Transformer model. Structureel overzicht van het Convolutioneel Neural Network-Transformer hybride model, waarbij convolutionele feature-extractie wordt gecombineerd met multi-head attention mechanismen voor geavanceerde temporele patroonherkenning. Klik hier om een grotere versie van deze figuur te bekijken.
Ontwikkeling van ensemblemodellen
Om de voorspellingsnauwkeurigheid en modelrobuustheid te verbeteren, hebben we een gewogen gemiddelde ensemblebenadering geïmplementeerd met voorspellingen van de vijf hybride deep learning-modellen: RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU en CNN-Transformer. Het ensemble werd samengesteld door geoptimaliseerde gewichten toe te kennen aan de voorspellingen van elk model, waarbij hogere gewichten werden gegeven aan modellen die superieure individuele prestaties aantonen, gemeten aan hun R²-scores. Deze wegingsstrategie zorgt ervoor dat nauwkeurigere modellen aanzienlijk bijdragen aan de uiteindelijke prognose, terwijl ze toch de complementaire sterke punten van alle architecturen benutten. De ensemble-output werd vervolgens geëvalueerd met behulp van standaard prestatie-metrics: R², gemiddelde absolute fout (MAE) en gemiddelde kwadratische fout (MSE) om de voorspellende nauwkeurigheid, consistentie en generalisatiecapaciteit te beoordelen. Dit deep learning-ensemble heeft als doel temporele feature-extractie vanuit meerdere perspectieven te integreren, waardoor een grotere nauwkeurigheid en robuustheid wordt bereikt dan bij elk enkel hybride model op zichzelf.
Wiskundige formulering van de ensembletechniek:
Laat
Vertegenwoordigen de verzameling basismodellen die overeenkomen met CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU en CNN-Transformer.
Elk basismodel Mi levert een voorspelling op:
De meta-feature-matrix voor stapeling wordt gevormd als:
De Ridge Regression meta-learner schat de uiteindelijke voorspelling als:
waarbij:
-- wi zijn de geleerde stapelgewichten
-- w₀ is de biasterm
Om overfitting te voorkomen, minimaliseert Ridge Regressie de volgende geregulariseerde verliesfunctie:

waarbij:
-- yj = ware target voor dej-de steekproef
-- N = totaal aantal monsters
-- α = regularisatieparameter die de gewichtskrimp regelt
De ensemblevoorspelling wordt verkregen als:
waarbij de gewichten wi automatisch worden geleerd door de Ridge-verliesfunctie te minimaliseren.