$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Zbieranie i opis zbiorów danych
Zbiór danych (Rysunek 1) użyty w tych badaniach obejmuje historyczne zapisy kluczowe dla prognozowania energii słonecznej. Zbiór danych obejmuje codzienne dane operacyjne z elektrowni słonecznej o mocy 50 MW, obsługiwanej przez Megha Engineering and Infrastructures Limited (MEIL), zlokalizowanej w pobliżu Anantapur, Andhra Pradesh, Indie, wykorzystującej technologię skoncentrowanej energii słonecznej (CSP) z parabolicznego trumfu (CSP), która wychwytuje bezpośrednie normalne promieniowanie (DNI) i przekazuje ciepło za pomocą płynu transferowego ciepła (HTF) do wytwarzania energii elektrycznej. Zbiór danych został zebrany od 1 stycznia 2015 do 3 października 2025 i zawiera siedem kluczowych atrybutów obejmujących informacje czasowe, pomiary promieniowania słonecznego oraz produkcję energii. Atrybuty temporalne obejmują 'Date', podająca datę kalendarzową w standardowym formacie, 'Year' oznaczającą rok zbierania danych, 'Month' oznaczającą numer miesiąca, 'Day' oznaczającą dzień miesiąca oraz 'Julian Day' oferującą sekwencyjny system numeracji dni w ciągu roku do ciągłej analizy czasowej. Główną zmienną wejściową meteorologiczną jest 'DNI SUM' mierzona w kWh/m², która reprezentuje całkowite bezpośrednie normalne napromieniowanie (DNI), czyli skumulowaną energię słoneczną otrzymywaną na metr kwadratowy powierzchni kolektora, będącą kluczowym czynnikiem wpływającym na efektywność konwersji cieplnej elektrowni CSP. Zmienna docelowa "Rzeczywista generacja", mierzona w kWh, obejmuje moc elektryczną wyprodukowaną przez elektrownię CSP, odzwierciedlając wynik procesu konwersji energii słonecznej na cieplno-elektryczną. Te cechy łącznie umożliwiają kompleksową analizę wydajności elektrowni, w tym określanie efektywności cieplnej, modelowanie konwersji DNI na zasilanie, identyfikację wpływów atmosfery i chmur poprzez analizę spektralną, wdrożenie protokołów korekcji zerowych błędów do kalibracji czujników oraz rozwój zaawansowanych hybrydowych modeli prognozowania uczenia maszynowego dla optymalizacji planowania operacyjnego w czasie rzeczywistym oraz zwiększenia ogólnej efektywności i niezawodności elektrowni CSP. Szczegóły dotyczące rośliny dostępne na: https://solarpaces.nrel.gov/project/megha-solar-plant

Rysunek 1: Pięć górnych wierszy zbioru danych. Przykładowe dane pokazujące początkowe wpisy zbioru danych o wytwarzaniu energii słonecznej, pokazujące cechy wejściowe i zmienne docelowe używane do trenowania i oceny modelu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
Przygotowanie danych
Badanie wykorzystuje dane z szeregów czasowych z produkcji energii słonecznej obejmujące okres od 1 stycznia 2015 do 10 marca 2025 roku. Aby uwzględnić potencjalne problemy z jakością danych w pierwszych latach i skupić się na nowszych wzorcach, dane zostały przefiltrowane od 1 stycznia 2017 roku. Kolumny czasowe (Data, Rok, Dzień) zostały usunięte na podstawie wstępnej analizy korelacji, która wykazała znikomą wartość predykcyjną. Brakujące wartości zostały sprowadzone za pomocą techniki średniej kroczącej w celu utrzymania ciągłości czasowej przy jednoczesnym minimalizowaniu zniekształceń ukrytych wzorców. Utworzono trzy funkcje opóźniające przy użyciu zmiennej docelowej (Rzeczywista Generacja (kW/h)) do rejestrowania zależności czasowych.
Dzielenie zbiorów danych
Aby ustalić zrównoważone i reprezentatywne kohorty treningowe, walidacyjne i testowe, zbiór danych przygotowany wcześniej został podzielony na segmenty metodą stratyfikowanego próbkowania. Takie podejście zapewniło, że 70% danych (2091 danych) było przeznaczonych do treningu, podczas gdy zarówno zestawy walidacyjne, jak i testowe stanowiły po 15% (448 danych na zestaw).
Normalizacja danych
Cechy zostały ustandaryzowane za pomocą StandardScaler, natomiast wartości docelowe zostały znormalizowane za pomocą MinMaxScaler do zakresu [0,1] dla stabilności sieci neuronowej.
Trening modelu
Do prognozowania energii słonecznej zaimplementowano pięć modeli hybrydowych (Random Forest-BiLSTM, CNN-LSTM, CNN-LSTM, CNN-GRU oraz CNN-Transformer). Dane wejściowe zostały przekształcone w format sekwencyjny przekształcony w (próbki, kroki czasowe, cechy), gdzie kroki czasowe = 1 dla większości modeli, z wyjątkiem CNN-LSTM, które korzystało z przesuwanego okna 15 kroków. Zestawy treningowe, walidacyjne i testowe były skalowane, zachowując porządek czasowy, aby uniknąć wycieku danych. Wszystkie modele są trenowane z 32 rozmiarami partii i 30 epochami.
Losowy Las – BiLSTM (Rysunek 2): Proponowany model hybrydowy łączy dwukierunkową sieć pamięci długotrwałej (BiLSTM) z regressorem lasu losowego (RF), aby poprawić dokładność predykcji. Po pierwsze, model BiLSTM jest trenowany na danych szeregów czasowych wejściowych, aby wychwycić wzorce czasowe i generować wstępne prognozy. Następnie obliczane są błędy resztkowe (różnice między rzeczywistymi a przewidywanymi wartościami) z BiLSTM. Model lasu losowego jest następnie trenowany na oryginalnych cechach wejściowych, aby nauczyć się i przewidzieć te reszty. Aby zwiększyć wydajność modelu RF, sześć najważniejszych cech wybiera się na podstawie wskaźników ważności cech. Na koniec skorygowane przewidywanie uzyskuje się przez dodanie rezydualów prognozowanych RF do wyjść BiLSTM. To hybrydowe podejście wykorzystuje zdolność modelowania sekwencji BiLSTM oraz siłę uczenia się zespołowego Random Forest, aby osiągnąć lepszą generalizację i wydajność predykcyjną.
Niech
będzie sekwencją wejściową w kroku czasowym t.
Prognoza BiLSTM:

,
Obliczenia resztkowe:
Uczenie rezydualne z wykorzystaniem lasu losowego: Niech Z⊂X będzie cechami wybranych z najwyższych k, biorąc pod uwagę ważność cech.
Ostateczna prognoza: 

Rysunek 2: Architektura modelu pamięci długoterminowej krótkotrwałej lasu i dwukierunkowej pamięci. Schemat ilustrujący hybrydową architekturę RF-BiLSTM, pokazujący integrację przetwarzania cech Random Forest z dwukierunkowymi warstwami LSTM do uczenia sekwencji czasowej. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
CNN-LSTM (Rysunek 3): Model hybrydowy CNN-LSTM rozpoczyna się od przetwarzania sekwencji wejściowej za pomocą warstwy 1D splotowej do wyodrębnienia lokalnych cech przestrzennych, po czym następuje aktywacja LeakyReLU, normalizacja wsadowa i maksymalna pula. Wyodrębnione cechy są następnie przekazywane przez stos trzech warstw LSTM w celu poznania zależności czasowych, a normalizacja warstw i dropout są stosowane po pierwszych dwóch LSTM do regularizacji. Ostateczne wyjście LSTM jest przekazywane przez w pełni połączone, gęste warstwy z aktywacją i przerwą, a ostatecznie mapowane na wyjście za pomocą pojedynczego neuronu.
Niech
będzie ciągiem wejściowym, gdzie T to okno czasowe, a F to liczba cech.
Operacja CNN:
Maksymalna pula:
Komórka LSTM:




Efekt:

Rysunek 3: Architektura modelu CNN-LSTM. Strukturalna reprezentacja hybrydowego modelu hybrydowego sieci neuronowej i długiej pamięci krótkotrwałej, demonstrująca ekstrakcję cech splotowych, po której następuje jednokierunkowe przetwarzanie sekwencji czasowej. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
CNN-BiLSTM (Rysunek 4): Model hybrydowy CNN-BiLSTM najpierw wyodrębnia cechy przestrzenne za pomocą 1D warstwy splotowej z 32 filtrami, następnie następuje normalizacja wsadowa i maksymalna pula w celu zmniejszenia wymiarowości. Następnie wyjście jest przekazywane przez stos dwóch dwukierunkowych warstw LSTM, aby uchwycić długoterminowe zależności czasowe zarówno w kierunku do przodu, jak i do tyłu. Regularizacja jest stosowana poprzez dropout i normalizację wsadową. Gęsta warstwa z 128 neuronami udoskonala wyuczoną reprezentację, zanim ostateczna warstwa wyjściowa odwzoruje ją na pojedynczą przewidywaną wartość.
Operacja CNN:
Maksymalna pula: 
Dwukierunkowy LSTM: 
,
Efekt: 

Rysunek 4: Architektura modelu CNN-BiLSTM. Diagram architektury modelu Convolutional Neural Network-Bidirectional Long Short-Term Memory, podkreślający połączenie warstw splotowych z dwukierunkowym przetwarzaniem rekurencyjnym dla lepszego przechwytywania zależności czasowej. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
CNN-GRU (Rysunek 5): Model hybrydowy CNN-GRU zaczyna się od warstwy Conv1D, używającej jądra o rozmiarze 1 do wyodrębniania cech przestrzennych z pojedynczego kroku czasowego. Maksymalna pula zmniejsza wymiary przestrzenne. Następnie następuje stos warstw GRU: pierwsza zwraca sekwencje do przechwytywania zależności czasowych, a druga podsumowuje sekwencję do zwartej reprezentacji. Końcowa warstwa gęsta generuje wartość przewidywaną. Regularizacja dropoutu jest stosowana pomiędzy warstwami GRU, aby zapobiec nadmiernemu dopasowaniu.
Operacja CNN:
Maksymalna pula:
Komórka GRU:



Efekt:

Rysunek 5: Architektura modelu CNN-GRU. Schemat hybrydowego modelu hybrydowego Convolutional Neural-Network-Gated Recurrent Unit, pokazujący konwolucyjne wstępne przetwarzanie zintegrowane z warstwami GRU dla efektywnego modelowania czasowego. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
Transformator CNN (Rysunek 6) Model transformatora CNN zaczyna się od warstwy Conv1D do wyodrębnienia lokalnych cech z sekwencji wejściowej, po której następuje warstwa maksymalnej puli. Funkcje te są przekazywane przez blok enkodera transformatora, składający się z wielogłowowego mechanizmu samouwagi, normalizacji warstw oraz gęstej sieci feed-forward. Następnie stosuje się globalne ułożenie średniej przed ostateczną warstwą gęstą wygenerowaną prognozę. Ta architektura została zaprojektowana tak, aby przechwytywać zarówno wzorce przestrzenne (za pośrednictwem CNN), jak i globalne zależności (za pomocą uwagi transformera).
Operacja CNN:
Wielogłowa samouwaga:

Gdzie: Q, K, V = XWQ, XWK, XWV i dk to wymiar kluczy.
Sieć Feed Forward:

Warstwy Dodaj i Normuj:


Efekt:

Rysunek 6: Architektura modelu CNN-Transformer. Przegląd strukturalny hybrydowego modelu splotowej sieci neuronowej-transformatora, łączący ekstrakcję cech splotowych z mechanizmami uwagi wielogłowej dla zaawansowanego rozpoznawania wzorców czasowych. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
Rozwój modelu zespołowego
Aby zwiększyć dokładność prognozowania i odporność modelu, wdrożyliśmy podejście zespołu ważonego średniego, wykorzystując prognozy z pięciu hybrydowych modeli głębokiego uczenia: RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU oraz CNN-Transformer. Zespół został skonstruowany poprzez przypisanie zoptymalizowanych wag do przewidywań każdego modelu, przy czym wyższe wagi przyznawano modelom wykazującym lepsze indywidualne wyniki, mierzone ich wynikami R². Ta strategia ważenia zapewnia, że dokładniejsze modele znacząco przyczyniają się do ostatecznej prognozy, jednocześnie wykorzystując komplementarne zalety wszystkich architektur. Następnie wynik zespołu został oceniony za pomocą standardowych wskaźników wydajności: R², średniego błędu bezwzględnego (MAE) oraz średniego błędu kwadratowego (MSE), aby ocenić jego dokładność predykcyjną, spójność i zdolność uogólnienia. Ten zespół głębokiego uczenia ma na celu integrację ekstrakcji cech czasowych z różnych perspektyw, osiągając tym samym większą dokładność i odporność niż jakikolwiek pojedynczy model hybrydowy w izolacji.
Matematyczne sformułowanie techniki zespołowej:
Niech
reprezentują zestaw modeli bazowych odpowiadających CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU oraz CNN-Transformer.
Każdy model bazowy Mi daje prognozę:
Macierz meta-cech dla układania składów powstaje jako:
Meta-uczący się regresji grzbietu szacuje ostateczną prognozę jako:
gdzie:
-- wi są wyuczonymi wagami stosowania
-- w₀ to wyraz stronniczości
Aby uniknąć nadmiernego dopasowania, regresja grzbietowa minimalizuje następującą uregulowaną funkcję strat:

gdzie:
-- yj = prawdziwy cel dlaj-tej próbki
-- N = całkowita liczba próbek
-- α = parametr regularizacji kontrolujący kurczenie się masy
Predykcja zespołowa uzyskuje się jako:
gdzie wagi wi są automatycznie uczone przez minimalizację funkcji straty grzbietu.