$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Colección y descripción de conjuntos de datos
El conjunto de datos (Figura 1) utilizado en esta investigación comprende registros históricos cruciales para la predicción de la energía solar. El conjunto de datos comprende datos operativos diarios de una planta solar térmica concentrada de 50 MW operada por Megha Engineering and Infrastructures Limited (MEIL), situada cerca de Anantapur, Andhra Pradesh, India, utilizando tecnología de energía solar concentrada (CSP) por canal parabólico que captura la irradiancia normal directa (DNI) y transfiere calor mediante un fluido de transferencia de calor (HTF) para generar electricidad. El conjunto de datos se recopiló desde el 1 de enero de 2015 hasta el 3 de octubre de 2025 y contiene siete atributos clave que capturan información temporal, mediciones de irradiancia solar y producción de energía. Los atributos temporales incluyen 'Fecha', que proporciona la fecha del calendario en formato estándar, 'Año' que indica el año de recogida de datos, 'Mes' que representa el número del mes, 'Día' que denota el día del mes, y 'Día Juliano', que ofrece un sistema secuencial de numeración de días a lo largo del año para un análisis temporal continuo. La variable meteorológica principal de entrada es la 'SUMA DNI' medida en kWh/m², que representa la Irradiancia Normal Directa (DNI) total, la energía solar acumulada recibida por metro cuadrado de la superficie del colector, que sirve como factor crítico que influye en la eficiencia de conversión térmica de la planta de CSP. La variable objetivo 'Generación Real', medida en kWh, captura la potencia eléctrica producida por la planta CSP, reflejando el resultado del proceso de conversión de energía solar a térmica a eléctrica. Estos atributos permiten en conjunto un análisis exhaustivo del rendimiento de la planta, incluyendo la determinación de la eficiencia térmica, modelado de conversión de DNI a energía, identificación de influencias atmosféricas y de la cubierta de nubes mediante análisis espectral, implementación de protocolos de corrección de cero errores para la calibración de sensores y desarrollo de modelos avanzados híbridos de predicción por aprendizaje automático para optimizar la planificación operativa en tiempo real y mejorar la eficiencia y fiabilidad general de la planta CSP. Detalles de la planta disponibles en: https://solarpaces.nrel.gov/project/megha-solar-plant

Figura 1: Cinco filas principales del conjunto de datos. Datos de muestra que muestran las entradas iniciales del conjunto de datos de generación de energía solar, mostrando características de entrada y variables objetivo utilizadas para el entrenamiento y evaluación del modelo. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
Preparación de datos
El estudio utiliza datos de series temporales de generación solar que abarcan desde el 1 de enero de 2015 hasta el 10 de marzo de 2025. Para tener en cuenta posibles problemas de calidad de los datos en los primeros años y centrarse en patrones más recientes, los registros fueron filtrados a partir del 1 de enero de 2017. Las columnas temporales (Fecha, Año, Día) se eliminaron basándose en un análisis preliminar de correlación que mostró un valor predictivo despreciable. Se imputaron valores faltantes utilizando una técnica de media móvil para mantener la continuidad temporal minimizando la distorsión de los patrones subyacentes. Se crearon tres características de lag utilizando la variable objetivo (Generación Actual (kW/h)) para capturar dependencias temporales.
División de conjuntos de datos
Para establecer cohortes de entrenamiento, validación y pruebas equilibradas y representativas, el conjunto de datos preprocesado se segmentó mediante un método de muestreo estratificado. Este enfoque aseguraba que el 70% de los datos (2091 datos) se asignaran para entrenamiento, mientras que tanto los conjuntos de validación como de prueba representaban el 15% (448 datos por conjunto).
Normalización de datos
Las características se estandarizaron usando StandardScaler, mientras que los valores objetivo se normalizaron mediante MinMaxScaler en el rango [0,1] para la estabilidad de la red neuronal.
Formación en modelos
Se implementaron cinco modelos híbridos (Random Forest-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU y CNN-Transformer) para la previsión de generación solar. Los datos de entrada se reestructuraron en un formato secuencial remodelado como (muestras, pasos de tiempo, características) donde pasos de tiempo = 1 para la mayoría de los modelos, excepto CNN-LSTM, que usaba una ventana deslizante de 15 pasos. Los conjuntos de entrenamiento, validación y pruebas se escalaron preservando el orden temporal para evitar fugas de datos. Todos los modelos se entrenan con 32 unidades de lote y 30 épocas.
Random Forest-BiLSTM (Figura 2): El modelo híbrido propuesto combina una red Bidireccional Long Short-Term Memory (BiLSTM) con un regresor de Random Forest (RF) para mejorar la precisión de la predicción. Primero, el modelo BiLSTM se entrena con los datos de series temporales de entrada para capturar patrones temporales y generar predicciones iniciales. Después de esto, se calculan los errores residuales (diferencias entre valores reales y predichos) de la BiLSTM. A continuación, se entrena un modelo de Bosque Aleatorio sobre las características de entrada originales para aprender y predecir estos residuos. Para mejorar el rendimiento del modelo RF, se seleccionan las seis características más importantes en función de las puntuaciones de importancia de las características. Finalmente, la predicción corregida se obtiene sumando los residuos predichos por RF a las salidas de BiLSTM. Este enfoque híbrido aprovecha la capacidad de modelado de secuencias de BiLSTM y la fuerza del aprendizaje conjunto de Random Forest para lograr una mejor generalización y rendimiento predictivo.
Sea
la secuencia de entrada en el paso de tiempo t.
Predicción BiLSTM:
,
Cálculo residual:
Aprendizaje residual usando Random Forest: Sea Z⊂X las k características principales seleccionadas usando la importancia de la característica.
Predicción final:

Figura 2: Arquitectura del modelo de memoria a corto plazo bidireccional de bosque aleatorio. Diagrama esquemático que ilustra la arquitectura híbrida RF-BiLSTM, mostrando la integración del procesamiento de características en Bosque Aleatorio con capas LSTM bidireccionales para el aprendizaje temporal de secuencias. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
CNN-LSTM (Figura 3): El modelo híbrido CNN-LSTM comienza procesando la secuencia de entrada utilizando una capa convolucional 1D para extraer características espaciales locales, seguida de una activación LeakyReLU, normalización por lotes y max pooling. Las características extraídas se pasan entonces por una pila de tres capas LSTM para aprender dependencias temporales, aplicándose la normalización de capa y el dropout tras los dos primeros LSTM para regularización. La salida final del LSTM se pasa a través de capas densas totalmente conectadas con activación y dropout, y finalmente se mapea a la salida usando una sola neurona.
Sea
la secuencia de entrada, donde T es la ventana temporal y F es el número de características.
Operación de CNN:
Máximo pooling:
Celda LSTM:




Salida:

Figura 3: Arquitectura del modelo CNN-LSTM. Representación estructural del modelo híbrido Red Neuronal Convolucional-Memoria a Corto Plazo Largo, que demuestra la extracción de características convolucionales seguida de procesamiento de secuencias temporales unidireccionales. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
CNN-BiLSTM (Figura 4): El modelo híbrido CNN-BiLSTM extrae primero características espaciales utilizando una capa convolucional 1D con 32 filtros, seguida de la normalización por lotes y el max pooling para reducir la dimensionalidad. La salida se pasa entonces a través de una pila de dos capas LSTM bidireccionales para capturar dependencias temporales a largo plazo tanto hacia adelante como hacia atrás. La regularización se aplica mediante dropout y normalización por lotes. Una capa densa con 128 neuronas refina la representación aprendida antes de que la capa final de salida la mapea a un único valor predicho.
Operación de CNN:
Máximo pooling: 
LSTM bidireccional:
,
Salida:

Figura 4: Arquitectura del modelo CNN-BiLSTM. Diagrama arquitectónico del modelo Convolutional Neural Network-Bidirectional Long Short-Term Memory, que destaca la combinación de capas convolucionales con procesamiento recurrente bidireccional para una captura de dependencias temporales mejorada. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
CNN-GRU (Figura 5): El modelo híbrido CNN-GRU comienza con una capa Conv1D que utiliza un tamaño de núcleo de 1 para extraer características espaciales de ese único paso de tiempo. El max pooling reduce las dimensiones espaciales. Esto es seguido por una pila de capas GRU: la primera devuelve secuencias para capturar dependencias temporales, y la segunda resume la secuencia en una representación compacta. Una capa densa final produce el valor predicho. Se aplica regularización de dropout entre capas GRU para evitar el sobreajuste.
Operación de CNN:
Máximo pooling:
Celda GRU:



Salida:

Figura 5: Arquitectura del modelo CNN-GRU. Esquema del modelo híbrido Convolutional Neural Network-Gated Recurrent Unit, mostrando preprocesamiento convolucional integrado con capas GRU para un modelado temporal eficiente. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
CNN-transformador (Figura 6) El modelo CNN-transformador comienza con una capa Conv1D para extraer características locales de la secuencia de entrada, seguida de una capa de max pooling. Estas características pasan a través de un bloque codificador Transformer que consiste en un mecanismo de auto-atención multicabezal, normalización de capas y una red densa de avance de avance. A continuación, se aplica un agrupamiento de media global antes de que una capa densa final produzca la predicción. Esta arquitectura está diseñada para capturar tanto patrones espaciales (vía CNN) como dependencias globales (mediante atención Transformer).
Operación de CNN:
Autoatención con múltiples cabezas:

Donde: Q, K, V = XWQ, XWK,XW V y dk es la dimensión de las tonalidades.
Red de Feed Forward:

Capas de Añadir y Normalizar:


Salida:

Figura 6: Arquitectura del modelo CNN-Transformador. Visión estructural del modelo híbrido Red Neuronal Convolucional-Transformer, que incorpora la extracción de características convolucionales con mecanismos de atención multicabeza para el reconocimiento avanzado de patrones temporales. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
Desarrollo de modelos de conjunto
Para mejorar la precisión de las predicciones y la robustez del modelo, implementamos un enfoque de conjunto de media ponderada utilizando predicciones de los cinco modelos híbridos de aprendizaje profundo: RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU y CNN-Transformer. El conjunto se construyó asignando pesos optimizados a las predicciones de cada modelo, con pesos más altos a los modelos que demostraban un rendimiento individual superior, medido por sus puntuaciones R². Esta estrategia de ponderación garantiza que los modelos más precisos contribuyan de forma más significativa a la previsión final, aprovechando al mismo tiempo las fortalezas complementarias de todas las arquitecturas. La salida del conjunto se evaluó entonces utilizando métricas estándar de rendimiento: R², error absoluto medio (MAE) y error cuadrático medio (MSE) para evaluar su precisión predictiva, consistencia y capacidad de generalización. Este conjunto de aprendizaje profundo pretende integrar la extracción temporal de características desde múltiples perspectivas, logrando así mayor precisión y robustez que cualquier modelo híbrido individual en aislamiento.
Formulación matemática de la técnica de conjunto:
Sea
representan el conjunto de modelos base correspondientes a CNN-RF-BiLSTM, CNN-LSTM, CNN-BiLSTM, CNN-GRU y CNN-Transformer.
Cada modelo base Mi produce una predicción:
La matriz meta-característica para apilamiento se forma como:
El meta-aprendiz de Regresión de Ridge estima la predicción final como:
Dónde:
-- wi soy el aprendido a apilar pesos
-- w₀ es el término de sesgo
Para evitar el sobreajuste, la regresión de crestas minimiza la siguiente función de pérdida regularizada:

Dónde:
-- yj = objetivo verdadero para la jésima muestra
-- N = número total de muestras
-- α = parámetro de regularización que controla la contracción del peso
La predicción del conjunto se obtiene como:
donde los pesos wi se aprenden automáticamente minimizando la función de pérdida de cresta.