Tecnológico Nacional de

México
Centro Nacional de Investigación
y Desarrollo Tecnológico

Avance de tercer semestre

Control de un convertidor bidireccional para la
óptima carga/descarga de baterías de un
sistema fotovoltaico.

presentada por
Ing. Ximena Olivares Pérez

como requisito para la obtención del grado de

Maestra en Ciencias en Ingeniería
Electrónica

Director de tesis
Dr. Jesús Darío Mina Antonio

Revisores
Dr. Abraham Claudio Sánchez
M. en C. Alfredo González Ortega
Cuernavaca, Morelos, México. Diciembre de 2024.

I
Índice

1. Introducción................................................................................................................................- 1 -
1.1. Objetivos..............................................................................................................................- 2 -
1.1.1. General.............................................................................................................................- 2 -
1.1.2. Específicos.......................................................................................................................- 2 -
2. Dimensionamiento del sistema...........................................................................................- 2 -
2.1. Dimensionamiento del conjunto de paneles...........................................................- 4 -
2.2. Dimensionamiento de baterías....................................................................................- 6 -
2.3. Dimensionamiento convertidor Boost.......................................................................- 7 -
2.4. Algoritmo MPPT..................................................................................................................- 9 -
2.5. Dimensionamiento convertidor bidireccional DC-DC Buck+Boost................- 11 -
3. Cronograma de actividades................................................................................................- 15 -
Referencias.......................................................................................................................................- 15 -
Anexo A..............................................................................................................................................- 16 -

II
 LISTA DE ACRÓNIMOS
1. PV panel fotovoltaico
2. SVF sistema fotovoltaico
3. BESS sistema de almacenamiento de energía en baterías (por sus siglas en
inglés Battery Energy Storage System)

III
 Índice de Figuras
FIGURA 1. ESQUEMA GENERAL DE UN SISTEMA FOTOVOLTAICO CON BESS CONECTADO A LA RED. -1-
FIGURA 2. MUESTREO DEL CONSUMO DE ENERGÍA -3-
FIGURA 3. COMPARATIVA DEL CONSUMO INDIVIDUAL Y EN CONJUNTO DE LA CARGA SELECCIONADA. -3-
FIGURA 4. BLOQUE DE ARREGLO DE PANELES SIMULINK -5-
FIGURA 5. A) CURVAS DE CORRIENTE-VOLTAJE, B) CURVAS DE POTENCIA-VOLTAJE DEL ARREGLO DE PANELES. - 5 -
FIGURA 6. CONSUMO DE CARGAS Y POTENCIA ENTREGADA POR EL ARREGLO DE PANELES. -6-
FIGURA 7. CIRCUITO ELÉCTRICO DEL CONVERTIDOR BOOST -7-
FIGURA 8. SIMULACIÓN DEL CONVERTIDOR BOOST. -8-
FIGURA 9. SEÑALES OBTENIDAS DEL CONVERTIDOR BOOST. TENSIÓN DE SALIDA(ROJO), TENSIÓN DE
ENTRADA(AMARILLO), CORRIENTE DE INDUCTOR (AZUL), CORRIENTE DE SALIDA (VERDE). -9-
FIGURA 10. ALGORITMO INC. - 10 -
FIGURA 11. CIRCUITO CONVERTIDOR BOOST CON MPPT INC. - 10 -
FIGURA 12. CIRCUITO CONVERTIDOR BUCK+BOOST. - 11 -
FIGURA 13. CONVERTIDOR OPERANDO EN MODO BOOST. - 11 -
FIGURA 14. CONVERTIDOR OPERANDO EN MODO BUCK - 11 -
FIGURA 15. SIMULACIÓN DEL CONVERTIDOR EN MODO BOOST - 13 -
FIGURA 16. SEÑALES OBTENIDAS DEL CONVERTIDOR MODO BOOST. TENSIÓN DE SALIDA(ROJO), TENSIÓN DE
ENTRADA(AMARILLO), CORRIENTE DE INDUCTOR (AZUL), CORRIENTE DE SALIDA (VERDE). - 13 -
FIGURA 17. SIMULACIÓN DEL CIRCUITO EN MODO BUCK. - 14 -
FIGURA 18. SEÑALES OBTENIDAS DEL CONVERTIDOR MODO BUCK. TENSIÓN DE ENTRADA(AMARILLO), TENSIÓN DE
SALIDA(VERDE), CORRIENTE DE INDUCTOR (AZUL), CORRIENTE DE SALIDA (ROJO). - 14 -

IV
1. Introducción
Debido a la creciente demanda de energía eléctrica que se tiene en la actualidad se
presenta un impacto negativo en nuestro ambiente debido a que al generar energía
eléctrica de una manera convencional se hace el uso de hidrocarburos. Es por ello
que se buscó nuevas alternativas para generar energía que sean más amigables
con en ambiente, dando así paso al uso de las energías renovables, actualmente las
energías renovables más utilizadas son la eólica y la solar, siendo la energía solar la
tecnología que mayor margen de progreso tiene debido a su mayor eficiencia y bajo
mantenimiento [1], [2].
Sin embargo, la irradiancia del sol y la fuerza del viento se mantienen en constante
cambio. En consecuencia, tanto la energía eólica como la solar poseen
intermitencias en su generación de energía debido a sus características
naturales [1] , [2] .

Teniendo en cuenta esto, es necesario encontrar alternativas para asegurar que

siempre haya energía disponible para la demanda energética que se tenga. En este
contexto una solución prometedora a estos problemas de intermitencia es la
implementación de sistemas de almacenamiento de energía en baterías (BESS por
sus siglas en inglés) en los sistemas fotovoltaicos como el que se observa en la
figura 1, los cuales no solo ayudan a los paneles a entregar la energía que se es
requerida, sino que también permiten tener una mayor independencia de la red
eléctrica.

Figura 1. Esquema general de un sistema fotovoltaico con BESS conectado a la red.

1
 1.1. Objetivos

[Link]
Implementar un controlador para el proceso de carga y descarga del banco
de baterías de un sistema fotovoltaico interconectado a red, ante las condiciones
cambiantes de la irradiancia y demanda de energía.

[Link]íficos
 Revisión del estado del arte para definir el caso de estudio y las condiciones
de operación.
 Diseñar los diferentes subsistemas del caso de estudio.
 Diseñar un sistema de gestión de energía.
 Diseñar el control de carga/descarga de baterías.
 Evaluar el desempeño del sistema.

2. Dimensionamiento del sistema

En el siguiente apartado se presenta el dimensionamiento que se lleva realizado
de los diversos elementos que integran el sistema propuesto.
Para definir los valores que serán utilizados en el consumo de cargas en el sistema
propuesto se basó en los sistemas realizados en (). En los cuales implementaban los
sistemas fotovoltaicos integrados en edificios (BIPV por sus siglas en inglés) y
sistemas de energía doméstica, en donde el consumo de cargas se promediaba
entre los valores de 12kWh/día y 11kWh/día.
Por lo que, para las cargas utilizadas en el sistema propuesto se basaron en datos
obtenidos de un muestreo por dos semanas en el mes de septiembre de manera
real de una casa, la cual presenta una demanda de energía promedio de 3.82 kW/h
al día, en la figura 2 se presenta el muestreo del consumo de energía de dicha casa.

2
 Consumo de energía
6
5.5
5
4.5
4
3.5
3
kWh
2.5
2
1.5
1
24 24 24 24 24 24 24 24 24 24 24 24 24 24 24
/ 20 /20 /20 /20 /20 /20 /20 /20 /20 /20 /20 /20 /20 /20 /20
2 3 4 5 6 7 8 9 0 1 2 3 4 5 6
9/ 9/ 9/ 9/ 9/ 9/ 9/ 9/ 9/1 9/1 9/1 9/1 9/1 9/1 9/1

Día

Figura 2. Muestreo del consumo de energía

Teniendo en consideración estos valores y las cargas propuestas en () se seleccionó

un día típico y se asumió un conjunto de tres casas que conforman un edificio
resultando en un consumo total de 11.8 kWh por día, en la figura 3 se muestra el
consumo de cargas del día seleccionado tanto de la casa individualmente, así como
el conjunto de ellas.

Gráfica del consumo de cargas

1
0.9
0.8
0.7
0.6
0.5
kW

0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

Hora

Consumo individual Consumo del conjunto de casas

Figura 3. Comparativa del consumo individual y en conjunto de la carga seleccionada.

3
 2.1. Dimensionamiento del conjunto de paneles
Una vez definidos los valores de las cargas se continuo con el dimensionamiento de
los demás elementos del sistema, empezando con el dimensionamiento del
conjunto de paneles.
Para ello se utiliza la siguiente formula

Et
N p=
HSP (0.9∗PPv ) max

N p : Número de panelesET :Cantidad de cargas con una estimación del75 %

HSP : Número de horas solares pico P Pv : Potenciamáxima
max

Para calcular las HSP de la región seleccionada, es decir, el promedio de horas

solares pico en condiciones estándar de prueba de los módulos fotovoltaicos
2
5.21 kW /m /día
HSP= 2
=5.21
1 k W /m /día
La potencia máxima ( P Pv ¿ que puede entregar el panel, varía según la marca y el
max

modelo que se escoja, en este caso se seleccionó la marca JA Solar y el modelo

JAM54S30-420/MR. En la tabla 1 se muestran los parámetros eléctricos que este
W
pude proporcionar en condiciones donde se presenta una irradiancia de 1000 2 y
m
una temperatura de 25°C.
Tabla 1. Características eléctricas del modelo de panel JAM54S30-420/MR.

Modulo JAM54S30-420/MR
Potencia máxima 420 W
Voltaje de circuito abierto (Voc) 37.58 V
Corriente de corto circuito (Isc) 14.10 A
Eficiencia del módulo 21.5%

Ya definidos los parámetros se puede conocer el número de paneles necesarios

para satisfacer las cargas, en este análisis se obtiene el siguiente número de
paneles.

Et 15.85 k Wh
N p= = =8.04 → 8 paneles
HSP (0.9∗PPv ) 5.21(0.9∗420 W )
max

Para la simulación del conjunto de paneles para el sistema se utilizó el software

Simulink de MATLAB® en el cual se comprobó que los valores de tensión y potencia
fuesen los esperados del conjunto de 8 paneles solares. Para ello se utilizó el bloque
de arreglo de paneles de Simulink, figura 4, en el cual se observaron las
características de las curvas de corriente-voltaje (I-V) así como de potencia-voltaje

4
 W
(P-V) cuando se asume una irradiancia de 1000 2 y una temperatura de 25C°, en la
m
figura 3 se observan dichas curvas del conjunto de paneles y en la tabla 2 se
muestran los parámetros del arreglo fotovoltaico simulado.

Figura 4. Bloque de arreglo de paneles Simulink

a)

b)

Figura 5. a) curvas de corriente-voltaje, b) curvas de potencia-voltaje del arreglo de paneles.

Tabla 2. Valores obtenidos de la simulación del arreglo de paneles.

Valores obtenidos
Tensión máxima 254 V
Corriente máxima 14.1 A
W 3360 W
Potencia a 1000 2 25°C
m
W 2687 W
Potencia a 800 2 25°C
m
W 1669 W
Potencia a 500 2 25°C
m
W 1329 W
Potencia a 400 2 25°C
m

5
 W 989 W
Potencia a 300 2 25°C
m
W 317 W
Potencia a 100 2 25°C
m

2.2. Dimensionamiento de baterías

Para el análisis del dimensionamiento de la batería para el sistema de
almacenamiento (BESS) se tomó en cuenta la diferencia existente entre la demanda
de potencia y la potencia que puede entregar el arreglo de paneles cuando este es
capaz de suministrar la energía requerida de manera independiente.
En la figura 6 se muestra la cantidad de potencia entregada por el arreglo de
paneles y el consumo de cargas por las tres viviendas, el área sombreada
representa la potencia que se utilizará para almacenar la energía sobrante y
suministrarla cuando el SFV no sea capaz de hacerlo.

Gráfica con el consumo de las tres viviendas

2.3
2.1
1.9
1.7
1.5
1.3
kW

1.1
0.9
0.7
0.5
0.3
0.0999999999999999
-0.1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

Horas

Potencia Panel Consumo de cargas

Figura 6. Consumo de cargas y potencia entregada por el arreglo de paneles.

De acuerdo a (“ENERGÍA SOLAR FOTOVOLTAICA”” RID CONNECTED PV SYSTEMS

WITH BATTERY ENERGY STORAGE SYSTEMS DESIGN GUIDELINES”) para determinar
la tensión para el almacenamiento de energía se recomienda que para cargas
mayores a los 5 kW se utilice una batería igual o mayor a 48V.

6
Para la capacidad que debe tener la batería se determinó con la siguiente formula

Ee 7.73 kW
C= = =57 . 51 Ah
V Bat∗DoD 192V ∗0.7

C :Capacidad Ee : Energía excedente DoD : Profundidad de descarga ( DoD por sus siglas en inglés )
V Bat :Voltaje de la batería
Teniendo estos valores definidos se buscó una batería que presentara las
características especificadas, siendo así que se seleccionó una batería con los
siguientes parámetros eléctricos.

Parámetros batería
Voltaje nominal 48 V
Capacidad 60 Ah
Para el análisis de este sistema se plantea el uso de 8 baterías en serie para
obtener un voltaje total de 192V, esto para mantener una adecuada proporción en
el ciclo de trabajo (D) y así el convertidor trabaje de manera eficiente, pues bien, la
ganancia tiene una relación con el ciclo de trabajo, si se tiene una ganancia alta el
ciclo de trabajo deberá serlo también.

2.3. Dimensionamiento convertidor Boost

Debido que es necesario elevar los niveles de tensión y potencia del arreglo de
paneles se hace uso de un convertidor elevador Boost (figura 4) para obtener los
valores de tensión del Bus de CD el cual está a la salida de este. Para realizar los
cálculos de los elementos pasivos se asumieron las características del arreglo de
paneles cuando se presenta el punto de máxima potencia, es decir, cuando se tiene
W
una irradiancia de 1000 2 y una temperatura de 25°C.
m
En la tabla 3 se muestran los parámetros iniciales del diseño.

Figura 7. Circuito eléctrico del convertidor Boost

Tabla 3. Parámetros iniciales del diseño del convertidor Boost

Parámetros Boost
Frecuencia de
10 kHz
conmutación ( f s)

7
 Ciclo de trabajo (D) 0.435
V_Pv 254 V
I_Pv 13.77 A
P_Pv 3.5 kW
V_BusDC 450 V
I_BusDC 7.7 A
P_BusDC 3.5 kW

El valor del inductor se calcula como

V Pv∗D 254 V ∗0.435

LB = = =16.059 mH
∆ i LB∗f s 0.7 A∗10 k Hz
El valor del capacitor se calcula como

I BusDC∗D 7.46 A∗0.435

C B= = =14.434 µF
∆V C B∗f s 22.5 V ∗10 k Hz

Para los valores de ∆ i LB y ∆ V C B se asumió un rizo de corriente y un rizo de voltaje

del 5%.
Se realizó una simulación para comprobar que los valores calculados de los
elementos pasivos del convertidor DC-DC sean los adecuados para elevar los
valores de tensión, corriente y potencia del arreglo de paneles, en la figura 8 se
observa la simulación realizada en Simulink de MATLAB®, en la figura 9 se
presentan las gráficas obtenidas.

Figura 8. Simulación del convertidor Boost.

8
Figura 9. Señales obtenidas del convertidor Boost. Tensión de salida(rojo), tensión de entrada(amarillo), corriente de inductor
(azul), corriente de salida (verde).

2.4. Algoritmo MPPT

Debido al constante cambio que presenta la irradiancia del sol, existen
variaciones en la cantidad de energía que puede proporcionar el arreglo de paneles,
sin embargo, existe un punto en donde se entrega la máxima potencia que el panel
puede entregar (MPP por sus siglas en inglés), el cual varía según la irradiancia que
se presente en el panel solar.
Es por ello que es necesario el seguimiento de punto de máxima potencia (MPPT por
sus siglas en inglés) este monitorea continuamente los valores de tensión y
corriente en los paneles para mantenerlos operando en ese punto óptimo, para así
extraer la máxima cantidad de energía posible. Existen diversas maneras en las
cuales se puede obtener el MPPT, los métodos más empleados son el método de
perturbar y observar (P&O) y el método de conductancia incremental (INC por sus
siglas en inglés).
Para este análisis fue seleccionado el método INC el cual se basa en el diagrama de
flujo presentado en la figura 10. En la figura 11 se observa el circuito realizado en
Simulink de MATLAB® que incluye el algoritmo MPPT.

9
 Figura 10. Algoritmo INC.

Figura 11. Circuito convertidor Boost con MPPT INC.

Para la obtención del MPPT en la simulación se utilizó el bloque de función MATLAB

en el cual se introdujo un código con el algoritmo INC (Anexo A).

2.5. Dimensionamiento convertidor bidireccional DC-DC Buck+Boost

Para el BESS se utiliza un convertidor bidireccional que se encarga de realizar
la carga/ descarga de la batería cuando sea requerido. En este análisis se propone

10
el uso del convertidor bidireccional Buck+Boost (figura 12), el cual posee dos
modos de operación dependiendo el interruptor que este activo.

Figura 12. Circuito convertidor Buck+Boost.

Para operar el convertidor en modo Boost (figura 13) se activa Q BB 1 , el cual

funciona como switch y Q BB 2 funciona como diodo, en este modo se realiza la
descarga de la batería. Por otro lado, cuando Q BB 2 es activado se opera en modo
Buck (figura 14) y se realiza la carga de la batería.

Figura 13. Convertidor operando en modo Boost.

Figura 14. Convertidor operando en modo Buck

Se realizó el análisis en los dos modos de operación para así obtener las ecuaciones
necesarias para realizar el diseño del convertidor. Para el modo de operación Boost
se tienen los siguientes parámetros iniciales del diseño.

11
Tabla 4. Parámetros de diseño inicial modo Boost.

Parámetros Modo Boost

Frecuencia de
10 kHz
conmutación ( f s)
Ciclo de trabajo ( D B t) 0.573
V_Bat 192 V
I_Bat 18.22 A
P_Bat 3.5 kW
V_BusDC 450 V
I_BusDC 7.7 A
P_BusDC 3.5 kW

El valor del inductor se calcula como

V B at ∗DB t 192 V∗0.57

LB B= = =12.07 mH
∆ i LB B∗f s 0.91 A∗10 kHz
El valor del capacitor se calcula como

I BusDC∗DB t 7.46 A∗0. 57

C B B 1= = =1 9. 81 µF
∆ V C B B 1∗f s 22.5V ∗10 kHz
Para el modo de operación Buck se tienen los siguientes parámetros iniciales del
diseño.
Tabla 5. Parámetros de diseño inicial modo Buck.

Parámetros Modo Buck

Frecuencia de
10 kHz
conmutación ( f s)
Ciclo de trabajo ( D B k ) 0.426
V_Bat 192 V
I_Bat 18.22 A
P_Bat 3.5 kW
V_BusDC 450 V
I_BusDC 7.7 A
P_BusDC 3.5 kW

El valor del inductor se calcula como

(V ¿ ¿ BusDC −V Bat )∗D Bk (450 V −192V )∗0. 426

LBB = = =12.07 mH ¿
∆ i LBB∗f s 0.91 A∗10 kHz
El valor del capacitor se calcula como

∆ i LBB∗D Bk 0.91 A∗0.4 26

C BB 2= = =1.1 8 µF
∆ V C BB 2∗f s 22.5V ∗10 kHz

12
Se considero para ambos modos un rizo de corriente y un rizo de voltaje del 5%.
Para determinar que los valores obtenidos son adecuados para el correcto
funcionamiento del convertidor bidireccional se hicieron dos pruebas en Simulink de
MATLAB®, para comprobar el correcto funcionamiento del modo Boost y para el
correcto funcionamiento del modo Buck

Figura 15. Simulación del convertidor en modo Boost

Figura 16. Señales obtenidas del convertidor modo Boost. Tensión de salida(rojo), tensión de entrada(amarillo), corriente de
inductor (azul), corriente de salida (verde).

13
 Figura 17. Simulación del circuito en modo Buck.

Figura 18. Señales obtenidas del convertidor modo Buck. Tensión de entrada(amarillo), tensión de salida(verde), corriente de
inductor (azul), corriente de salida (rojo).

3. Cronograma de actividades
A continuación se presenta el cronograma de actividades a realizar, en donde las actividades en las
cueles se esta trabajando son marcadas en azul, las actividades realizadas en verde y las actividades
por realizar en naranja.

14
 Tercer semestre Cuarto semestre
Actividades
Agosto Septiembre Octubre Noviembre Diciembre Enero Febrero Marzo Abril Mayo Junio Julio

Revisión del estado

del arte

Redacción de tesis

Definir las
condiciones de
operación para el
caso de estudio

Dimensionamiento
del arreglo de
paneles

Implementación de
MPPT

Dimensionamiento
de batería

Dimensionamiento
de elementos
pasivos convertidor
bidireccional CD-CD
Modelado del
convertidor
Bidireccional

Selección del sistema

de control para la
carga/descarga de
baterías

Revisión del
desempeño del
sistema
Presentación de
avance

Referencias
[1] M. Yekini Suberu, M. Wazir Mustafa, and N. Bashir, “Energy storage systems for renewable energy
power sector integration and mitigation of intermittency,” 2014, Elsevier Ltd. doi:
10.1016/[Link].2014.04.009.

[2] Larbi Tebessi University, Institute of Electrical and Electronics Engineers. Algeria Section, and Institute of
Electrical and Electronics Engineers, Proceedings, 2019 1st International Conference on Sustainable
Renewable Energy Systems and Applications (ICSRESA) : Tebessa, Algeria, 04 to 05 December 2019.

15
[12] M. Argyrou, C. Spanias, C. Marouchos, S. Kalogirou, and P. Christodoulide, Energy management and
modeling of a grid-connected BIPV system with battery energy storage. IEEE, 2019.

[13] M. M. Hasan, S. Jaman, T. Geury, and O. Hegazy, “Design and Simulation of a Grid-Connected Two-Stage
Bidirectional Converter for a Combined PV-Stationary Energy Storage System,” in 2022 2nd International
Conference on Sustainable Mobility Applications, Renewables and Technology, SMART 2022, Institute of
Electrical and Electronics Engineers Inc., 2022. doi: 10.1109/SMART55236.2022.9990149.

[14] P. Athira Anil and M. Rajesh, “PV-ESS connected to a grid tied Inverter with a High gain DC-DC
converter,” in ICCISc 2021 - 2021 International Conference on Communication, Control and Information
Sciences, Proceedings, Institute of Electrical and Electronics Engineers Inc., Jun. 2021. doi:
10.1109/ICCISc52257.2021.9484974.

[15] Y. Shao and X. Wu, “Modeling and Control of Bi-directional Converter based on Battery Management in
Photovoltaic System,” in 2022 IEEE International Conference on Advances in Electrical Engineering and
Computer Applications, AEECA 2022, Institute of Electrical and Electronics Engineers Inc., 2022, pp. 335–
341. doi: 10.1109/AEECA55500.2022.9919027.

[16] D. Xu, G. Wang, W. Yan, and X. Yan, “A novel adaptive command-filtered backstepping sliding mode
control for PV grid-connected system with energy storage,” Solar Energy, vol. 178, pp. 222–230, Jan.
2019, doi: 10.1016/[Link].2018.12.033.

[17] G. Sharma and R. Gupta, “Solar PV Based Grid Scale Battery Energy Storage System with IBC and DAB,” in
2023 7th International Conference on Computer Applications in Electrical Engineering-Recent Advances:
Sustainable Transportation Systems, CERA 2023, Institute of Electrical and Electronics Engineers Inc.,
2023. doi: 10.1109/CERA59325.2023.10455215.

[18] S. Pamu, R. Thogaru, and A. Mitra, Resiliency Analysis of Grid connected Solar PV System with Battery
Support. IEEE, 2020.

[19] S. Marhraoui, A. Abbou, Z. Cabrane, and S. Krit, “Sliding Mode Control for PV Grid-Connected System
With Energy Storage,” 2020, pp. 168–199. doi: 10.4018/978-1-7998-4381-8.ch008.

[20] Metwally S, Zhang X, Ali S, and Krishnamoorthy H, Solar PV and BESS based Home Energy System.
Institute of Electrical and Electronics Engineers, 2019.

Anexo A

16
17
18