$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Coleta e descrição do conjunto de dados
O conjunto de dados (Figura 1) utilizado nesta pesquisa compreende registros históricos cruciais para a previsão de energia solar. O conjunto de dados compreende dados operacionais diários de uma usina solar térmica concentrada de 50 MW operada pela Megha Engineering and Infrastructures Limited (MEIL), localizada próxima a Anantapur, Andhra Pradesh, Índia, utilizando tecnologia de Energia Solar Concentradora (CSP) por canal parabólico que captura Irradiância Normal Direta (DNI) e transfere calor por meio de um Fluido de Transferência de Calor (HTF) para gerar eletricidade. O conjunto de dados foi coletado de 01 de janeiro de 2015 a 3 de outubro de 2025 e contém sete atributos-chave que capturam informações temporais, medições de irradiância solar e produção de energia. Os atributos temporais incluem 'Data', que fornece a data do calendário em formato padrão, 'Ano' indicando o ano da coleta dos dados, 'Mês' representando o número do mês, 'Dia' denotando o dia do mês, e 'Dia Juliano', oferecendo um sistema sequencial de numeração de dias ao longo do ano para análise temporal contínua. A principal variável meteorológica de entrada é a 'SUM' do DNI, medida em kWh/m², que representa a Irradiância Normal Direta total (DNI), a energia solar acumulada recebida por metro quadrado da superfície do coletor, servindo como fator crítico que influencia a eficiência da conversão térmica da usina CSP. A variável alvo 'Geração Real', medida em kWh, captura a produção de energia elétrica produzida pela usina CSP, refletindo o resultado do processo de conversão de energia solar para térmica para elétrica. Esses atributos permitem coletivamente uma análise abrangente do desempenho da planta, incluindo determinação de eficiência térmica, modelagem de conversão DNI-para-energia, identificação de influências atmosféricas e de cobertura de nuvens por meio de análise espectral, implementação de protocolos de correção de erro zero para calibração de sensores e desenvolvimento de modelos avançados híbridos de previsão de aprendizado de máquina para otimizar o planejamento operacional em tempo real e melhorar a eficiência e confiabilidade geral da central CSP. Detalhes da planta disponíveis em: https://solarpaces.nrel.gov/project/megha-solar-plant

Figura 1: Cinco principais linhas do conjunto de dados. Dados de exemplo mostrando as entradas iniciais do conjunto de dados de geração de energia solar, exibindo características de entrada e variáveis-alvo usadas para treinamento e avaliação do modelo. Por favor, clique aqui para ver uma versão ampliada desta figura.
Preparação de dados
O estudo utiliza dados de séries temporais de geração solar que vão de 01 de janeiro de 2015 a 10 de março de 2025. Para considerar possíveis problemas de qualidade dos dados nos primeiros anos e focar em padrões mais recentes, os registros foram filtrados a partir de 01 de janeiro de 2017. Colunas temporais (Data, Ano, Dia) foram removidas com base em análises preliminares de correlação que mostraram valor preditivo negligenciável. Valores ausentes foram imputados usando uma técnica de média móvel para manter a continuidade temporal enquanto minimizava a distorção dos padrões subjacentes. Três recursos de lag foram criados usando a variável alvo (Geração Atual (kW/h)) para capturar dependências temporais.
Divisão de conjuntos de dados
Para estabelecer coortes equilibradas e representativas de treinamento, validação e teste, o conjunto de dados pré-processado foi segmentado usando um método de amostragem estratificada. Essa abordagem garantiu que 70% dos dados (2091 dados) fossem alocados para treinamento, enquanto tanto os conjuntos de validação quanto de teste representavam 15% (448 dados por conjunto).
Normalização de dados
As funcionalidades foram padronizadas usando o StandardScaler, enquanto os valores-alvo foram normalizados via MinMaxScaler para a faixa [0,1] para estabilidade da rede neural.
Treinamento em modelos
Cinco modelos híbridos (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU e CNN-Transformer) foram implementados para previsão de geração solar. Os dados de entrada foram reestruturados em um formato sequencial, reformulado como (amostras, passos de tempo, características) onde passos de tempo = 1 para a maioria dos modelos, exceto CNN-LSTM, que usou uma janela deslizante de 15 passos. Os conjuntos de treinamento, validação e testes foram escalados preservando a ordem temporal para evitar vazamento de dados. Todos os modelos são treinados com 32 lotes e 30 épocas.
Random Forest-BiLSTM (Figura 2): O modelo híbrido proposto combina uma rede Bidirecional de Memória de Curto Prazo (BiLSTM) com um regressor de Floresta Aleatória (RF) para melhorar a precisão da previsão. Primeiro, o modelo BiLSTM é treinado com os dados de séries temporais de entrada para capturar padrões temporais e gerar previsões iniciais. Após isso, os erros residuais (diferenças entre valores reais e previstos) do BiLSTM são calculados. Um modelo de Floresta Aleatória é então treinado nas características originais de entrada para aprender e prever esses resíduos. Para melhorar o desempenho do modelo RF, as seis características mais importantes são selecionadas com base nas pontuações de importância das características. Por fim, a previsão corrigida é obtida adicionando os resíduos previstos por RF às saídas do BiLSTM. Essa abordagem híbrida aproveita a capacidade de modelagem de sequências do BiLSTM e a força de aprendizado em conjunto da Random Forest para alcançar melhor generalização e desempenho preditivo.
Seja
a sequência de entrada no passo de tempo t.
Previsão do BiLSTM:
,
Cálculo Residual:
Aprendizado Residual usando Random Forest: Seja Z⊂X as k principais características selecionadas usando a importância da característica.
Previsão final:

Figura 2: Arquitetura do modelo de memória de curto prazo bidirecional de floresta aleatória. Diagrama esquemático ilustrando a arquitetura híbrida RF-BiLSTM, mostrando a integração do processamento de características em Random Forest com camadas LSTM bidirecionais para aprendizado de sequências temporais. Por favor, clique aqui para ver uma versão ampliada desta figura.
CNN-LSTM (Figura 3): O modelo híbrido CNN-LSTM começa processando a sequência de entrada usando uma camada convolucional 1D para extrair características espaciais locais, seguida por uma ativação LeakyReLU, normalização em lote e max pooling. As características extraídas são então passadas por uma pilha de três camadas de LSTM para aprender dependências temporais, com normalização de camada e abandono aplicados após os dois primeiros LSTMs para regularização. A saída final do LSTM é passada por camadas densas totalmente conectadas com ativação e dropout e, finalmente, mapeada para a saída usando um único neurônio.
Seja
a sequência de entrada, onde T é a janela de tempo e F é o número de características.
Operação da CNN:
Máximo pooling:
Célula LSTM:




Saída:

Figura 3: Arquitetura do modelo CNN-LSTM. Representação estrutural do modelo híbrido de Rede Neural Convolucional-Memória Longa de Curto Prazo, demonstrando extração convolucional de características seguida pelo processamento de sequência temporal unidirecional. Por favor, clique aqui para ver uma versão ampliada desta figura.
CNN-BiLSTM (Figura 4): O modelo híbrido CNN-BiLSTM primeiro extrai características espaciais usando uma camada convolucional 1D com 32 filtros, seguida pela normalização em lote e max pooling para reduzir a dimensionalidade. A saída é então passada por uma pilha de duas camadas bidirecionais de LSTM para capturar dependências temporais de longo prazo tanto para frente quanto para trás. A regularização é aplicada por meio de dropout e normalização em lote. Uma camada densa com 128 neurônios refina a representação aprendida antes que a camada final de saída a mapeie para um único valor previsto.
Operação da CNN:
Máximo pooling: 
LSTM bidirecional:
,
Saída:

Figura 4: Arquitetura do modelo CNN-BiLSTM. Diagrama de arquitetura do modelo Convolutional Neural Network-Bidirectional Long Short-Term Memory, destacando a combinação de camadas convolucionais com processamento bidirecional recorrente para captura aprimorada de dependência temporal. Por favor, clique aqui para ver uma versão ampliada desta figura.
CNN-GRU (Figura 5): O modelo híbrido CNN-GRU começa com uma camada Conv1D usando um kernel tamanho 1 para extrair características espaciais do único passo de tempo. O max pooling reduz as dimensões espaciais. Isso é seguido por uma pilha de camadas GRU: a primeira retorna sequências para capturar dependências temporais, e a segunda resume a sequência em uma representação compacta. Uma camada densa final gera o valor previsto. A regularização dropout é aplicada entre as camadas de GRU para evitar o sobreajuste.
Operação da CNN:
Máximo pooling:
Célula GRU:



Saída:

Figura 5: Arquitetura do modelo CNN-GRU. Esquema do modelo híbrido Convolutional Neural Network-Gated Recurrent Unit, mostrando pré-processamento convolucional integrado com camadas GRU para modelagem temporal eficiente. Por favor, clique aqui para ver uma versão ampliada desta figura.
CNN-transformador (Figura 6) O modelo CNN-transformador começa com uma camada Conv1D para extrair características locais da sequência de entrada, seguida por uma camada de max pooling. Essas características são passadas por um bloco codificador Transformer composto por um mecanismo de autoatenção multi-cabeça, normalização de camadas e uma rede densa feed-forward. O agrupamento de média global é então aplicado antes que uma camada densa final produza a previsão. Essa arquitetura foi projetada para capturar tanto padrões espaciais (via CNN) quanto dependências globais (via atenção Transformer).
Operação da CNN:
Autoatenção com Múltiplas Cabeças:

Onde: Q, K, V = XWQ, XWK,XW V e dk é a dimensão das chaves.
Rede Feed Forward:

Camadas de Adicionar & Normalizar:


Saída:

Figura 6: Arquitetura do modelo CNN-Transformer. Visão estrutural do modelo híbrido Convolutional Neural Network-Transformer, incorporando extração convolucional de características com mecanismos de atenção multi-cabeça para reconhecimento avançado de padrões temporais. Por favor, clique aqui para ver uma versão ampliada desta figura.
Desenvolvimento de modelos de conjunto
Para aumentar a precisão da previsão e a robustez do modelo, implementamos uma abordagem de conjunto com média ponderada usando previsões dos cinco modelos híbridos de aprendizado profundo: RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU e CNN-Transformer. O conjunto foi construído atribuindo pesos otimizados às previsões de cada modelo, com maiores ponderações atribuídas a modelos que demonstrassem desempenho individual superior, medido por seus escores R². Essa estratégia de ponderação garante que modelos mais precisos contribuam de forma mais significativa para a previsão final, ao mesmo tempo em que aproveitam as forças complementares de todas as arquiteturas. A saída do conjunto foi então avaliada usando métricas padrão de desempenho: R², erro absoluto médio (MAE) e erro quadrático médio (MSE) para avaliar sua precisão preditiva, consistência e capacidade de generalização. Esse conjunto de aprendizado profundo visa integrar extração de características temporais de múltiplas perspectivas, alcançando assim maior precisão e robustez do que qualquer modelo híbrido isolado.
Formulação matemática da técnica de conjunto:
Seja
representam o conjunto de modelos base correspondentes a CNN-RF-BiLSTM, CNN-LSTM, CNN-BILSTM, CNN-GRU e CNN-Transformer.
Cada modelo base Mi produz uma previsão:
A matriz meta-feature para empilhamento é formada como:
O meta-aprendiz da Regressão de Ridge estima a previsão final como:
onde:
-- wi são os pesos empilhados aprendidos
-- w₀ é o termo de viés
Para evitar o sobreajuste, a Regressão de Crista minimiza a seguinte função de perda regularizada:

onde:
-- yj = alvo verdadeiro para a jésima amostra
-- N = número total de amostras
-- α = parâmetro de regularização controlando a redução do peso
A previsão do conjunto é obtida como:
onde os pesos wi são automaticamente aprendidos minimizando a função de perda de crista.