En esta investigación, se desarrolla un sistema autónomo de carga para vehículos eléctricos (VE) alimentado por un sistema fotovoltaico bifacial y un sistema de almacenamiento de energía por baterías (SAEB), como se muestra en la Figura 1, utilizando un método de seguidor del punto de máxima potencia (MPPT) controlado por una red neuronal artificial (ANN). Se utiliza un circuito equivalente de un solo diodo para modelar el arreglo fotovoltaico bifacial, teniendo en cuenta la temperatura, la resistencia en serie, la resistencia en derivación, la irradiancia solar y la irradiancia en el lado posterior. La salida del sistema fotovoltaico se suministra a la estación de carga para VE mediante un convertidor reductor-elevador de corriente continua (CC–CC), mientras que el SAEB contribuye a la regulación del voltaje del bus de corriente continua. El algoritmo de optimización basado en el comportamiento del ave secretaria (SBOA) se utiliza para optimizar las ganancias del controlador PI y los parámetros del filtro. La modelización de los componentes del sistema y la metodología propuesta se describen en las siguientes secciones.
Modelado de Componentes
1. Sistema fotovoltaico bifacial
Los módulos fotovoltaicos ayudan a alcanzar la corriente y el voltaje requeridos. En el presente estudio, se utilizaron dos módulos fotovoltaicos monofaciales para emular las características de generación de energía de matrices bifaciales bajo diferentes condiciones de albedo superficial9. La ganancia bifacial de la irradiación BGg se determina mediante la Ecuación (1)
(1)
Aquí, la irradiación frontal es Gf, y la irradiancia posterior es Gr. La corriente máxima fotovoltaica bifacial es Imp, y el voltaje es Vmp. Finalmente, la salida producida por el sistema bifacial con el número elegido de módulos en serie Ns y en paralelo Np se da mediante la Ecuación (2). El control del sistema fotovoltaico se muestra en Figura 2
(2)
2. Sistema de almacenamiento de batería (BSS)
Las baterías son un medio crucial para almacenar energía eléctrica18. Este trabajo utiliza baterías de iones de litio como sistema de almacenamiento de energía. Entre las opciones disponibles de baterías recargables, este estudio se centra en las baterías de iones de litio, que se consideran óptimas para vehículos eléctricos debido a su eficiencia y densidad de potencia. Las ecuaciones (3) y (4) describen la batería de iones de litio con restricciones.
(3)
Las especificaciones del PV, la batería de almacenamiento y el convertidor se resumen en la Tabla 2, mientras que la Tabla 3 presenta la distribución del flujo de potencia entre el arreglo fotovoltaico, la unidad de almacenamiento de energía y las cargas conectadas.
(4)
3. Convertidor de CC a CA
El inversor con punto neutro钳ado (NPC, por sus siglas en inglés) es ampliamente utilizado en sistemas modernos de conversión de energía para energías renovables. En comparación con la estructura convencional de dos niveles, la configuración NPC ofrece una menor distorsión armónica total (THD), una menor interferencia electromagnética y un esfuerzo de voltaje reducido en los dispositivos semiconductores. El enlace de corriente continua (CC) de esta topología se divide en tres niveles de potencial: +Vdc/2, 0 y -Vdc/2, mediante capacitores emparejados y diodos钳adores que forman la unión neutra. Al controlar adecuadamente la conmutación, el inversor genera tres estados de voltaje discretos en los terminales de salida, mejorando así la calidad de la forma de onda y la eficiencia de la conversión de potencia. El punto neutro del bus de CC en un inversor NPC desempeña un papel clave en el equilibrio de los voltajes de los capacitores y en garantizar el funcionamiento adecuado del inversor. El esquema de control del inversor se muestra en Figura 3.
4. Convertidor elevador de corriente continua a corriente continua
Los convertidores CC-CC elevan el voltaje de entrada en corriente continua a un voltaje de salida más alto controlando el ciclo de trabajo de un dispositivo de conmutación de alta frecuencia, manteniendo así la eficiencia. Tiene dos modos de operación. Durante el estado ON, el diodo está polarizado en inversa, y la fuente de entrada suministra energía al inductor, que la almacena como un campo magnético, mientras que el condensador de salida suministra energía a la carga. Durante el estado OFF, el diodo queda polarizado en directa. Aquí, la energía almacenada en el inductor se libera a través del diodo hacia el condensador conectado en la salida y hacia la carga, combinándose con la fuente de entrada para producir un voltaje de salida mayor que el voltaje de entrada. Al variar el ciclo de trabajo de la modulación por ancho de pulso (PWM), el voltaje de salida del convertidor varía según la relación
, donde D es el ciclo de trabajo. Se utiliza en sistemas solares, de celdas de combustible, de vehículos eléctricos y de almacenamiento.
Esquema de control de ANN para MPPT
En este trabajo, se selecciona un MPPT basado en una red neuronal artificial (ANN) para recopilar la potencia máxima generada por el sistema solar17. Los conjuntos de datos de irradiación y temperatura utilizados para el entrenamiento de la ANN se generaron en MATLAB/Simulink bajo condiciones ambientales variables de operación, correspondientes a niveles de irradiancia solar de 800–1000 W/m2 y rangos de temperatura de 20–25 °C, con un ciclo de trabajo (D). La ANN se entrena para minimizar el error cuadrático medio (MSE) entre las salidas obtenidas y requeridas (Op,
) con el fin de extraer la máxima potencia del sistema fotovoltaico bifacial. Aquí, la irradiación solar y la temperatura se utilizan como entradas para la ANN, y el ciclo de trabajo es la salida para controlar el convertidor elevador. La estructura del modelo ANN desarrollado para el MPPT, tal como se muestra en Figura 4, se considera en este trabajo. Tabla Suplementaria 1 presenta las ventajas y desventajas del método propuesto en comparación con otros métodos estándar.
Una red neuronal artificial consta de tres capas principales: una capa de salida (OPL), una capa de entrada (IPL) y una capa oculta (HIL). Esto ayuda a transferir los datos entre la IPL y la HIL. Posteriormente, se genera mediante los pesos en los enlaces que conectan la IPL y la HIL. En este contexto, los cálculos se realizan con un sesgo determinado aplicado a la variable HL, y los resultados obtenidos se acumulan en la variable OPL. Aquí, se selecciona una red neuronal artificial de tipo LMBP [17]. Los pesos de conexión se ajustan durante el entrenamiento mediante la medición del error para alcanzar la salida deseada. En este caso, se utiliza LMBP para entrenar redes neuronales artificiales con el MSE como función de rendimiento. El algoritmo LMBP utiliza las derivadas obtenidas para actualizar los pesos, lo que muestra propiedades de aprendizaje eficaz y convergencia acelerada.
Cada neurona en una red de perceptrón multicapa tiene una función de sumación y una función de activación. Sin embargo, existen ciertos pesos numéricos (wpk) que conectan estas neuronas a través de los niveles. Cuando las entradas se multiplican por los pesos según se especifica en la Ecuación (5). Se considera la función sigmoide no lineal con minimización del MSE dada en la Ecuación (6).
(5)
(6)
SBOA optimizó los parámetros del filtro y de control del sistema propuesto
Hoy en día, los algoritmos metaheurísticos desempeñan un papel clave en la resolución de problemas de ingeniería. Figura 5 muestra una clasificación de los algoritmos. El algoritmo SBOA23 imita el comportamiento del SB para sobrevivir en su hábitat natural. Las fases de exploración y explotación del comportamiento de caza del Ave Secretaria se modelan en el SBOA. La fase de exploración del algoritmo imita el comportamiento del SB al atrapar serpientes, mientras que su fase de explotación imita su comportamiento para esquivar depredadores como águilas. El proceso de solución comienza generando valores aleatorios dentro de sus límites respectivos para cada SB en la población, y calculando el valor de la función objetivo mediante la ejecución del modelo Simulink, tratando el valor de cada SB como un parámetro de diseño. En función del valor objetivo, cada SB cambia su posición para representar una solución mejor. A continuación se describe la modelización matemática del SBOA:
Fase inicial de preparación
Primero, es necesario identificar las soluciones iniciales que se utilizaron para comenzar la búsqueda de una minimización típica de la función objetivo F(Y). En este caso, la población aleatoria inicial Y = [Y1, Y2, ......YN] de SB se forma mediante N soluciones iniciales. La ecuación (7) inicializa la población y Yi representa una solución.
(7)
Donde los límites inferior y superior de las variables de elección se denotan mediante lb y ub. En [0, 1], r es un número aleatorio. N es la dimensión del problema. Además, el valor de aptitud de la solución Fi = F(Yi) se utiliza para cuantificar la calidad de Yi.
Estrategia de caza de los SB
Encontrar, comer y atacar a la presa son las tres fases principales de la caza del SB. Todo el proceso de caza se ha dividido en tres intervalos de tiempo iguales: t < 1/3T, 1/3T < t < 2/3T y 2/3T < t < T. Aquí, t es la iteración actual y T es el número máximo de iteraciones. Estos intervalos corresponden a las tres fases de la depredación del SB: búsqueda, consumo y ataque a la presa. Estas divisiones se basan en las estadísticas biológicas de las fases de caza del ave y en la duración temporal de cada fase. Por consiguiente, cada fase del SBOA se modela de la siguiente manera:
En busca de presas
El SB debe cazar presas desde una distancia segura durante esta fase. Para recopilar datos suficientes en toda la región de búsqueda, el primer paso en los algoritmos de optimización requiere una exploración más intensa. Utilizando las posiciones de los otros dos SB como guía, el SB puede investigar regiones posibles adicionales. Para aumentar aún más la diversidad del algoritmo, se introducen así procesos de mutación diferencial. Cuando t es menor que 1/3 de T, se utiliza la Ecuación (8) para actualizar la posición de cada individuo. Buscará, consumirá y atacará a las presas. Por consiguiente, cada fase del SBOA se modela de la siguiente manera:
(8)
Donde, Yi(t) es la posición actual del iº SB, Ynewi(t) proporciona la solución actualizada, dos individuos Yi(t), Yr2(t) se seleccionaron aleatoriamente de la población existente. Un vector aleatorio, para la estrategia de selección de parámetros, R1 tiene 1 × N seleccionado aleatoriamente del intervalo [0, 1] para lograr un equilibrio efectivo entre las fases de explotación y exploración. La solución actualizada se considera únicamente si produce un valor objetivo mejorado en comparación con la solución obtenida previamente.
Consumo de presas
SB utiliza una agilidad ágil y se mueve para mantenerse cerca de la serpiente tras detectar una posible presa. La paciencia de la presa se aprovechará para disminuir sus defensas, observándola y atrayéndola mientras giran a su alrededor. Otros SB ajustan su postura cerca de la presa. El éxito en la caza aumentará significativamente con este enfoque. La ecuación (9) ilustra cómo se utiliza el movimiento browniano (B) para simular el movimiento aleatorio de SB cuando 1/3 T < t < 2/3 T.
(9)
Donde B se implementa en MATLAB utilizando la distribución normal estándar mediante el vector generado aleatoriamente r(1,M). Luego, los SB utilizan la Ecuación (10) para actualizar sus posiciones.
(10)
Donde, Ymejor(t) es la mejor solución y el factor de escala exponencial exp((t/T)4) se elige para aumentar la capacidad de explotación en cada iteración.
Atacar a la presa
Después de consumir continuamente, la víctima se quedará sin energía. Las aves secretarias deben lanzar el ataque ahora. Aquí, se simulan diferentes métodos de ataque, como pasos incessantes o saltos largos esporádicos en un breve período de tiempo, utilizando la técnica de vuelo de Lévy. Las características se describen mediante las ecuaciones (11) y (12). En este momento, las opciones candidatas están cerca de la respuesta óptima porque las SB se acercan rápidamente a la presa. Este plan se implementará cuando t > 2/3T.
(11)
(12)
Donde, LF representa la estrategia de vuelo de Levy y
, s= 0.01, y φ = 1.5 son valores fijos. μ y υ son valores arbitrarios que se encuentran en [0, 1]. Donde, τ denota gamma y η = 0.5, el factor de ponderación 0.5 fue adoptado, según se propuso en la formulación original del SBOA, para regular el tamaño del paso y mantener un equilibrio entre grandes saltos exploratorios y refinamientos locales. Tabla 4 presenta los valores utilizados para la ANN y el algoritmo SB en el sistema desarrollado.
Estrategia de escape para SB
Al cazar otras presas en la naturaleza, los SB corren el riesgo de ser cazados. Águilas, halcones, zorros y chacales son los principales adversarios con los que deben enfrentarse. Deben emplear una variedad de técnicas de evasión para protegerse a sí mismos o a sus presas cuando perciben peligro. Este algoritmo simula los métodos de escape mediante la modelación de modos de carrera (D2) y camuflaje (D1). Al enfrentar adversarios, los pájaros secretario inicialmente se mezclan con su entorno para mantenerse a salvo. Aquí, los pájaros secretario procuran evitar óptimos locales en los algoritmos al actualizar sus posiciones alrededor de la presa (el mejor individuo). Utilizarán tácticas de vuelo o de carrera rápida para mantenerse a salvo si no logran evadir al enemigo. Como referencia, se selecciona un individuo aleatorio Yrand como líder para evitar quedar restringido a un óptimo local. Las ecuaciones (13) y (14) son utilizadas por los pájaros secretario para actualizar sus posiciones. El diagrama de flujo del SBOA se muestra en Figura 6.
(13)
(14)
Representación de las variables de diseño
Las variables del problema en este estudio incluyen los parámetros de diseño, como las ganancias del controlador PI de los convertidores CC-CC, así como las resistencias y las inductancias de los filtros y convertidores. La ecuación (15) muestra la representación de las variables de control. Los límites de los parámetros de control se enumeran en la Tabla Suplementaria 2.
(15)
Función de ajuste (FF)
Se seleccionó la minimización de THD como la función objetivo (Obj) en este estudio. El valor máximo de FF corresponde al valor mínimo de Obj, dado en las ecuaciones (16) y (17).
(16)
(17)