See discussions, stats, and author profiles for this publication at: [Link]
net/publication/266047110
SISTEMA DE CONTROL DE BAJO COSTO PARA EL SEGUIMIENTO DEL PUNTO
DE MÁXIMA POTENCIA EN SISTEMAS FOTOVOLTAICOS
Article
CITATION READS
1 562
3 authors:
Jose Manuel Andujar Marquez Juan Manuel Enrique
Universidad de Huelva Universidad de Huelva
406 PUBLICATIONS 3,524 CITATIONS 48 PUBLICATIONS 622 CITATIONS
SEE PROFILE SEE PROFILE
Miguel Angel Martínez Bohórquez
Universidad de Huelva
59 PUBLICATIONS 382 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Energy retrofitting and social housing instrumentation attending passive criteria. Winter performance View project
Proyecto Aurora View project
All content following this page was uploaded by Jose Manuel Andujar Marquez on 10 January 2015.
The user has requested enhancement of the downloaded file.
SISTEMA DE CONTROL DE BAJO COSTO PARA EL
SEGUIMIENTO DEL PUNTO DE MÁXIMA POTENCIA EN
SISTEMAS FOTOVOLTAICOS
Andújar Márquez, José M.
Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 217379, Fax: +34 959 217304, andujar@[Link]
Enrique Gómez, Juan M.
Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 217374, Fax: +34 959 217304, [Link]@[Link]
Martínez Bohórquez, Miguel A.
Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 217656, Fax: +34 959 217304, bohórquez@[Link]
Durán Aranda, Eladio
Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 217655, Fax: +34 959 217304, aranda@[Link]
Resumen Una única célula, es capaz de proporcionar una
tensión de apenas 0,5V y una potencia de entre 1 y
Aunque la fuente de energía sea limitada y gratuita, 2W. Para elevar la potencia suministrada, las células
el bajo rendimiento de un panel fotovoltaico se conectan en serie y en paralelo, dando lugar a lo
comercial (en torno al 15%) aconseja que éste esté que se denomina panel o módulo fotovoltaico.
siempre lo más próximo posible a su valor de
rendimiento máximo. Para ello, es fundamental que La expresión (1) describe la relación entre la tensión
la transferencia de potencia entre el grupo de (V) y la corriente (I) proporcionada por un módulo
paneles y la instalación sea la mejor posible. Un [8], donde nP y nS indican el número de células
factor que beneficia de forma acusada la eficiencia conectadas en paralelo y en serie; RP y RS las
de la instalación, es que el sistema de paneles pueda resistencias intrínsecas paralelo y serie asociadas a la
estar la mayor parte posible del tiempo funcionando célula, K es la constante de Boltzman (1.38·10-23 J/K)
en su punto de máxima potencia (Maximum Power y q es la carga del electrón. El factor A determina la
Point, MPP). En este trabajo, se presentan dos desviación de las características de una unión p-n
sistemas de control de bajo costo para el eficaz ideal, e IS es la corriente inversa de saturación, que
seguimiento del MPP ante cambios en la carga y en presenta una dependencia con la temperatura del
las condiciones ambientales. panel. IL representa la corriente (foto corriente)
generada por la radiación solar (G). Dicha corriente
exhibe una relación respecto de la radiación y la
Palabras Clave: algoritmo MPPT, convertidor DC- temperatura.
DC, Generador fotovoltaico, “P&O”, rendimiento.
⎡ ⎡ q⎛⎜ V + IRS ⎞⎟ ⎤ V + IRS ⎤
⎢ ⎢ ⎜⎝ nS nP ⎟⎠ AKT ⎥ nS nP ⎥
I = nP ⎢IL − Is ⎢e −1⎥ − ⎥ (1)
⎢ ⎢ ⎥ RP ⎥
1 INTRODUCCIÓN ⎢ ⎢⎣ ⎥⎦ ⎥
⎣ ⎦
El principio físico que rige el funcionamiento de una
célula fotovoltaica, puede ser explicado tomando P = I ·V (2)
como base una unión p-n de material semiconductor
(normalmente silicio) sensible a la luz solar. Cuando ⎡ ⎡ q ⎛⎜ V + P ·RS ⎞⎟ ⎤ V + P·RS ⎤
⎢ ⎢ ⎜⎝ nS V · nP ⎟⎠ ⎥ nS V ·nP ⎥
ésta incide sobre la unión (generalmente sobre el lado P = nP ·V ·⎢ I L − I s ⎢ e
AKT
− 1⎥ − ⎥ (3)
n), los fotones que la constituyen suministran la ⎢ ⎢ ⎥ RP ⎥
⎢ ⎢⎣ ⎥⎦ ⎥
cantidad de energía necesaria a los electrones de ⎣ ⎦
valencia del semiconductor, para romper el enlace ⎛ dP ⎞
que los mantenía unidos a sus átomos respectivos. ⎜ ⎟ =0 (4)
⎝ dV ⎠ MPP
La expresión (1) (considerando la dependencia de sus Como se ha indicado anteriormente, las curvas I-V y
parámetros con T y G), proporciona las denominadas P-V dependen de los valores de radiación incidente y
curvas I-V de un panel fotovoltaico, y el producto de temperatura en el panel. Esto hace que dichas curvas
ambas magnitudes la potencia suministrada. En la no sean estáticas y, consecuentemente, tampoco lo
figura 1 se muestran las curvas de un panel sea el punto de máxima potencia. Así, dicho punto
determinado, así como la dependencia de las recorre una trayectoria a lo largo del día, que
variables con la temperatura y la radiación solar. depende de los valores de radiación y temperatura en
Cada curva presenta un punto de máxima potencia cada instante. Por tanto, para obtener un rendimiento
(Vmpp), el cual proporciona el punto óptimo de óptimo del sistema fotovoltaico necesitamos disponer
operación, para un uso eficiente del panel. Otros dos de un sistema de seguimiento del punto de máxima
puntos importantes de esta curva, son la tensión a potencia (MPPT, Maximum Power Point Tracking).
circuito abierto (Voc) y la corriente de cortocircuito Básicamente, estos sistemas se componen de un
(Isc). La tensión en circuito abierto representa la convertidor DC-DC conectado entre el generador
máxima tensión proporcionada por el panel a fotovoltaico y la carga, y un sistema de control en el
corriente cero (sin carga o en vacío), mientras que la que se ejecuta un algoritmo de seguimiento del MPP.
corriente en cortocircuito representa la máxima Tal sistema de control genera una señal PWM con un
corriente extraíble del panel (carga en cortocircuito). ciclo de trabajo (δ) apropiado, que es insertada en el
convertidor DC-DC (ver Fig.3).
El convertidor presenta una impedancia de entrada
(Ri) que depende básicamente de la impedancia de la
carga (RL) y del ciclo de trabajo (δ)[2,3,4]. Por tanto,
el algoritmo MPPT buscará la δ óptima para que el
punto de trabajo del array fotovoltaico (intersección
de la característica I-V del array con la recta de carga
Ri) coincida con el punto de máxima potencia MPP.
En [4] se realiza un análisis detallado de las distintas
topologías de convertidores DC-DC para el
seguimiento del MPP, concluyendo que la topología
Buck-Boost (utilizada en este trabajo) es la más
apropiada para tal fin.
Figura 1: Características I-V y P-V de un panel
fotovoltaico.
El punto de máxima potencia (MPP, Maximum
Power Point) se obtiene resolviendo las ecuaciones
(3) y (4) (ver Fig.2). Este cálculo es lento y tedioso,
ya que ambas expresiones no poseen una solución
analítica, y por tanto han de ser resueltas mediante
métodos numéricos.
Figura 3: Esquema básico de un sistema MPPT.
Si bien la fuente de energía de un sistema
fotovoltaico es ilimitada y gratuita, un
funcionamiento ineficiente del sistema provoca un
sobredimensionamiento del mismo y un consiguiente
Figura 2. dP/dV frente a V. Cuando dP/dV=0 el encarecimiento de su precio. Hay que tener en cuenta
sistema está trabajando en el punto de máxima que el rendimiento de un panel fotovoltaico
potencia (MPP). comercial es muy bajo (en torno al 15%), con lo cual
el inadecuado seguimiento del MPP hace que este
rendimiento sea aún menor y, por consiguiente, la
superficie necesaria de paneles aumente.
Si bien en la literatura pueden encontrarse desarrollos
sobre algoritmos para el seguimiento del MPP (en
[5,6,7] se realiza un análisis detallado de los
algoritmos más usados), nuestra experiencia nos
muestra que si bien sus funcionamientos son
adecuados en horas centrales del día, sus eficiencias
disminuyen en horas tempranas y tardías, haciéndose
más patente en días de movimientos constantes de
nubes.
El principal objetivo de este trabajo es el desarrollo
de dos algoritmos para la mejora del seguimiento del
punto de máxima potencia. En la sección 2 se
detallan tales algoritmos. En la sección 3 se muestran
los resultados obtenidos bajo simulación del
funcionamiento de dichos algoritmos y, por último, Tabla 1: Casos posibles y ley de control para
se presentan algunas conclusiones. aproximar el sistema fotovoltaico al MPP.
bi = − sign ⎡⎣( ∆V panel ) ⋅ ( ∆Ppanel ) ⎤⎦
2 ALGORITMOS (5)
DESARROLLADOS PARA EL
SEGUIMIENTO DEL MPP. δ i = δ i −1 + bi ⋅ ∆δ (6)
En esta sección se desarrollan dos algoritmos para el
seguimiento del MPP. Éstos se derivan del conocido El valor de ∆δ se ha fijado a 0.004 (precisión de un
método “Perturbación y Observación (P&O)”. Los convertidor digital/analógico de 8 bits). En cada
algoritmos los hemos denominado “P&O Orientado”, iteración “i” se obtienen medidas de ∆Ppanel y ∆Vpanel.
por tener en cuenta si el sistema está trabajando ala Seguidamente, se ajusta el valor de δ para
izquierda o a la derecha del MPP, y “P&O aproximarlo al MPP. Obsérvese que no se necesitan
Modificado” que necesita un “conocimiento” de medidas de temperatura y radiación solar para hacer
cambios producidos en la radiación solar que llega al el seguimiento, lo cual abarata mucho el sistema de
panel. control. El diagrama de flujo de este algoritmo de
seguimiento se muestra en la Fig.5.
2.1 “P&O” ORIENTADO.
Considérese la Fig. 4 donde se muestra la
característica P-V de un panel fotovoltaico a una
temperatura y radiación determinadas. Supongamos
que debido a una modificación en el ciclo de trabajo
(δ) del convertidor, el sistema evoluciona desde VA
hasta VB (∆Vpanel>0 y ∆Ppanel>0). Como se observa en
la Fig. 4, Vmpp es mayor que VB, por tanto la tensión
de salida debe continuar aumentando. Supongamos
ahora que la perturbación ha movido al panel desde
VB hasta VA. En este caso, la tensión del panel debe
volver a aumentar para aproximarse a Vmpp . En la
Tab. 1 pueden observarse todas las combinaciones
posibles.
La única variable a la que el sistema de control tiene
acceso es el ciclo de trabajo (δ). Como se indica en Figura 4: Característica P-V de un panel
fotovoltaico. El algoritmo MPPT desarrollado
[2,3,4], cualquier incremento en δ implica una
es capaz de discernir si el sistema está
disminución en la resistencia de entrada Ri del
operando a derecha o izquierda del MPP y
convertidor DC-DC, y por tanto una disminución en
actuar en consecuencia.
la tensión de operación del generador fotovoltaico y
viceversa. El sistema parte de un valor inicial de δ=
0.5 que varía a incrementos constantes (∆δ) según las
expresiones (5) y (6).
radiación producirá un aumento en el valor de la
corriente del MPP (Impp) y viceversa (ver Fig. 6).
Cuando el algoritmo detecta una variación en la
radiación (G) por encima de un determinado umbral,
responde con un aumento inmediato deδ, de esta
forma, el convertidor DC-DC decrementa su
impedancia de entrada Ri, y por consiguiente obliga
al generador fotovoltaico a desplazarse a puntos de
más alta corriente cercanos al MPP. En la Fig.7 se
muestra el diagrama de flujo de este algoritmo de
seguimiento. Obsérvese que este algoritmo no
requiere medidas precisas de radiación (hecho que
encarecería bastante su implementación), sino que
necesita disponer del signo del incremento de
radiación en un intervalo de medidas. Un simple
fotodiodo puede servir para tal fin.
1200
14
radiación
1000 Impp
12
800 10
Radiación (W/m2)
Impp (A)
8
600
6
400
4
200
2
0 0
07:12 09:36 12:00 14:24 16:48
Hora local
Figura 6: Representación simultanea de la radiación
luminosa (G) e intensidad del punto de máxima
potencia (Impp). Obsérvese que ambas curvas varían
en la misma forma.
3 RESULTADOS DE SIMULACIÓN.
Figura 5: Diagrama de flujo del algoritmo MPPT
P&O Orientado. Para comprobar el funcionamiento de los dos
algoritmos propuestos se han realizado varias
2.2 ALGORITMO “P&O” MODIFICADO. simulaciones usando MATLAB con datos de
temperatura y radiación correspondientes a un día
Si se analiza con detenimiento el algoritmo anterior, completo (08-10-2000) (ver Fig. 8), adquiridos a
se observa que ante una variación brusca en la intervalos de 1 minuto en una instalación real con un
radiación incidente (ocasionada por ejemplo por el sistema de adquisición de datos “Hydra” (FLUKE)
paso de una nube) y por tanto en la potencia previo aislamiento galvánico.
suministrada por el generador fotovoltaico, el sistema
es incapaz de discernir la dirección apropiada del Como generador fotovoltaico para la simulación se
cambio en δ [5]. Para correguir este efecto, se ha elegido el modelo “SX60” (BP), cuyos parámetros
propone una ligera modificación al algoritmo se muestran en la tabla 2.
anterior. Como es sabido, las variaciones en la
radiación luminosa repercuten principalmente y de
manera directa en la corriente que suministra el
generador fotovoltaico [1]. Así, un incremento en la
‘SX60’ (BP):
% definición de los parámetros para una célula bp sx60
A = 1.2; % factor de idealidad de la unión PN.
Eg = 1.12; % energía de la banda de salto
np = 1; % número de módulos conectados en
paralelo
ns = 36; % número de células conectadas en series
Pmax = 60; % potencia máxima
Vmax = 16.8; % tensión de Pmax a 25ºC y 1sol
Imax = 3.56; % intensidad de Pmax a 25ºC y 1sol
Tr = 273 + 25; % temperatura de referencia de la
célula
Iscr = 3.80; % corriente de la célula en cortocircuito a
la temperatura y radiación de referencias
Voc = 21.06 / ns; % Tensión a circuito abierto por célula
a la temperatura de referencia
kv = -80e-3; % Coeficiente de temperatura de Voc
ki = 0.065/100; % Coeficiente de temperatura de Iscr
Tabla 2: Parámetros del módulo “SX60” (BP).
El sistema completo bajo simulación se compone de
un generador fotovoltaico conectado a un convertidor
DC-DC y a una carga de 20Ω, todo ello controlado
por los dos algoritmos de seguimiento detallados
anteriormente.
En la Fig. 9 se presentan simultáneamente la
trayectoria del punto real de máxima potencia y la
seguida por el sistema bajo el control del algoritmo
“P&O Orientado”, denominado en la gráfica P1.
Obsérvese que el algoritmo tiene un comportamiento
correcto cuando no existen variaciones bruscas en la
radiación (parte derecha de la Fig. 9), sin embargo,
este funcionamiento no es tan bueno cuando estas
variaciones son rápidas y de cierta magnitud (parte
izquierda de la Fig. 9).
En la Fig. 10 se vuelve a mostrar simultáneamente la
trayectoria del punto real de máxima potencia y la
seguida por el sistema bajo el control del algoritmo
“P&O Modificado”, denominado en la gráfica P2.
Obsérvese que este algoritmo tiene un
comportamiento mejor que el anterior, ya que
persigue de forma eficiente la trayectoria del MPP
incluso para variaciones bruscas en la radiación.
60 radiación 1200
Temperatura
Temperatura
1000
Radiación (W/m )
Figura 7: Diagrama de flujo del algoritmo MPPT 50
P&O Modificado. 800
Temperatura (ºC)
40 600
2
400
30
200
20
0
07:12 09:36 12:00 14:24 16:48
Hora local
Figura 8: Valores de temperatura y radiación.
Referencias
Pmpp
60
Pmpp
[1] Alonso G. M. C., 1998: El Generador
Fotovoltaico. CIEMAT, Serie de Ponencias,
50
Editorial CIEMAT, Madrid.
40
Potencia (W)
30
[2] Andujar J. M.,. Enrique J. E., Durán E.,
Bohórquez M. A. M., 2004. Sistema para la
20
generación automática de curvas I-V, P-V y
10 monitorización de módulos fotovoltaicos. Proc.
De las XXV Jornadas de Automática. ISBN:
0
84-688-7460-4. Ciudad Real (España).
07:12 09:36 12:00 14:24 16:48
Hora local
[3] Enrique J. M., Durán E., Sidrach M., Andújar J.
Figura 9: Trayectoria real del MPP y la obtenida con M., 2005. A new approach to obtain I-V and P-
el algoritmo “P&O Orientado”. V curves of PV panels by using DC-DC
converters. 31st IEEE Photovoltaic Specialists
Conference, January 3-7, 2005.
Pmpp
P2
[4] Enrique J.M., Sidrach-de-Cardona M., Durán
60
E., Bohórquez M.A, Carretero, J.E., Andújar
50 J.M.. An Optimun Configuration of DC/DC
Converters to Use in Photovoltaics Facilites.
40 20th European Photovoltaic Solar Energy
Potencia (W)
Conference and Exhibition, EUPVSEC 2005.
30
Junio 2005. Barcelona.
20
[5] Hohm D. P., Ropp M. E., 2002. Comparative
10 Study of Maximum Power Point Tracking
Algorithms. Progress in Photovoltaics:
0 Research and applications, November 2002.
07:12 09:36 12:00 14:24 16:48
Hora local
[6] Hua C., Lin J., Tzou H., 2003. MPP Control of
a Photovoltaic Energy System. ETEP Vol. 13,
No. 4, July/August 2003.
Figura 10: Trayectoria real del MPP y la obtenida
con el algoritmo “P&O Modificado”
[7] Hua C., Shen C., 1998: Study of maximum
power tracking techniques and control of
DC/DC converters for photovoltaic power
4 CONCLUSIONES. system. 29th Annual IEEE PESC, IEEE
Computer Soc. Press, pp. 86-93, New York.
En este trabajo, se han discutido los problemas de
merma del rendimiento que puede causar en una [8] J.A. Gow, C.D. Manning: "Development of a
instalación fotovoltaica, el inadecuado seguimiento photovoltaic array model for use in power-
del punto de máxima potencia del arreglo de paneles. electronics simulation studies". IEE Proc. on
Para conseguir la máxima eficiencia de la instalación Electric Power Applications, vol. 146, no.2, pp.
es necesario que durante el mayor tiempo posible, el 193-200, March 1999.
MPP sea el punto de trabajo. Para conseguir esto, se
han presentado dos algoritmos de muy bajo costo en
cuanto a su implementación práctica. El primero,
“P&O Orientado”, presenta un adecuado
funcionamiento en días con poco movimiento de
nubes, y además no necesita de medidas de
radicación y temperatura. El segundo, “P&O
Modificado”, aunque necesita de medidas “gruesas”
de radiación, presenta un excelente comportamiento
para cualquier condición ambiental.
View publication stats