XXV Jornadas de Automática ISBN: 84-688-7460-4

SISTEMA PARA LA GENERACIÓN AUTOMÁTICA DE CURVAS

I-V, P-V Y MONITORIZACIÓN DE MÓDULOS FOTOVOLTAICOS

Andújar Márquez, José M.

Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 017380, Fax: +34 959 017304, andujar@[Link]

Enrique Gómez, Juan M.

Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 017364, Fax: +34 959 017304, [Link]@[Link]

Durán Aranda, Eladio

Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 017388, Fax: +34 959 017304, aranda@[Link]

Martínez Bohórquez, Miguel A.

Dpto. de Ing. Electrónica, de Sist. Infor. y Aut.- Universidad de Huelva.
Tlf: +34 959 017388, Fax: +34 959 017304, Bohórquez@[Link]

Resumen Los convertidores DC/DC son ampliamente

utilizados en sistemas de generación fotovoltaicos
El control de instalaciones fotovoltaicas de modo como interface entre el panel fotovoltaico y la carga,
eficiente exige un conocimiento preciso de las permitiendo el seguimiento del punto de máxima
características I-V y P-V de los módulos transferencia de potencia (MPPT). Principalmente, su
fotovoltaicos. De hecho, estas curvas permiten situar misión es la de acondicionar la energía producida por
el punto de máxima transferencia de potencia (MPPT el conjunto de células siguiendo una estrategia de
) del panel fotovoltaico, cuyo seguimiento es un control determinada [1,2,3]. Existe una característica
problema de control aún no resuelto intrínseca asociada al proceso de conversión DC/DC,
satisfactoriamente. Hasta ahora, la obtención de las relacionada con la capacidad que presentan estos
curvas I-V y P-V se viene realizando mediante sistemas para reflejar o adaptar la resistencia vista
cargas electrónicas variables de elevado coste, lo desde su entrada y la conectada a su salida. En este
cual hace económicamente inviable la sentido, los convertidores son similares a los
monitorización en instalaciones de tipo doméstico. transformadores, cuando son empleados como
En este trabajo se propone una nueva metodología adaptadores de impedancia, excepto que en los
basada en el uso de convertidores DC/DC convertidores el parámetro de adaptación no es la
controlados mediante un sistema electrónico relación de espiras entre el secundario y el primario,
automático que genera ciclos de trabajo variables en sino el ciclo de servicio δ, que puede ser gobernado
todo el rango, lo cual permite la generación del electrónicamente [4,5,6,7]. Una nueva aplicación de
conjunto de valores (V,I) necesarios para la esta capacidad es propuesta en el presente trabajo
caracterización del módulo en uso. Respecto del para determinar la característica I-V de un panel
método tradicional de la carga activa, el propuesto fotovoltaico.
en este trabajo presenta las siguientes ventajas
principales: disminución de tamaño y coste, así como Los procedimientos utilizados experimentalmente
la posibilidad de monitorización y análisis de fallos para obtener la característica I-V de un panel
de la instalación. fotovoltaico, están basados en la conexión de una
impedancia variable en los terminales del panel y
seguidamente medir la tensión y la corriente
Palabras Clave: Control, monitorización, análisis de
proporcionadas por el módulo, mientras dicha
fallos, trazador I-V, panel fotovoltaico.
impedancia varía entre cortocircuito y circuito
abierto [8,9,10,11,12]. Esta misma estrategia será
implementada en este trabajo empleando un
convertidor DC/DC haciendo variar su ciclo de
1 INTRODUCCIÓN servicio, pero en este caso el barrido de los distintos
puntos se realiza, con reguladores no disipativos, lo
XXV Jornadas de Automática ISBN: 84-688-7460-4

cual cobra mayor importancia a medida que aumenta ideal, e IS es la corriente inversa de saturación, que
la potencia suministrada por el panel. presenta una dependencia con la temperatura del
panel. IL representa la corriente (foto corriente)
En la sección 2 se describe brevemente el principio generada por la radiación solar (S). Dicha corriente
de funcionamiento de una célula fotovoltaica así exhibe una relación sensiblemente lineal respecto de
como sus características I-V y P-V. la radiación y la temperatura.

En la sección 3 se realiza una descripción de la La expresión (1) (considerando la dependencia de sus

capacidad que presentan los distintos convertidores parámetros con T y S), proporciona las denominadas
DC/DC para realizar un barrido de las curvas I-V y curvas I-V teóricas de un panel fotovoltaico, y el
P-V de un panel fotovoltaico, y se demuestra que la producto de ambas magnitudes la potencia
estructura Reductora-Elevadora es la más adecuada suministrada. En la Figura 1 se muestran las curvas
para este fin. de un panel determinado, así como la dependencia de
estas variables con la temperatura y la radiación
En la sección 4, la metodología propuesta es solar.
validada mediante simulación. En la sección 5 se
aplica el diseño anterior a una instalación real en la Cada curva presenta un punto de máxima potencia
que se monitorizan la corriente, tensión, temperatura (Vp), el cual proporciona el punto óptimo de
y radiación, con objeto de efectuar un análisis de operación, para un uso eficiente del panel. Otros dos
fallos de la instalación. Por último se aportan algunas puntos importantes de esta curva, son la tensión a
conclusiones. circuito abierto (Voc) y la corriente de cortocircuito
(Isc). La tensión en circuito abierto representa la
máxima tensión proporcionada por el panel a
2 CARACTERÍSTICAS I-V Y P-V corriente cero (sin carga o en vacío), mientras que la
DE PANELES SOLARES. corriente en cortocircuito representa la máxima
corriente extraíble del panel (carga en cortocircuito).
El principio físico que rige el funcionamiento de una
célula fotovoltaica, puede ser explicado tomando
como base una unión p-n de material semiconductor
(normalmente silicio) sensible a la luz solar. Cuando
ésta incide sobre la unión (generalmente sobre el lado
n), los fotones que la constituyen suministran la
cantidad de energía necesaria a los electrones de
valencia del semiconductor para romper el enlace que
los mantenía unidos a sus átomos respectivos.

Una única célula es capaz de proporcionar una

tensión de apenas 0,5V y una potencia de entre 1 y
2W. Para elevar la potencia suministrada, las células
se conectan en serie y en paralelo, dando lugar a lo
Figura 1: Características I-V y P-V de un panel
que se denomina panel o módulo fotovoltaico.
fotovoltaico.
La expresión exponencial simplificada propuesta por
Gow y Manning [13], describe la relación entre la 3 CONVERTIDORES DC/DC
tensión (V) y la corriente (I) proporcionada por un COMO EMULADORES DE
módulo, ecuación (1). RESISTENCIA VARIABLE.

 
Los convertidores DC/DC son muy empleados en
 q  V + I RS   V IR
   nS nP   + S  aplicaciones donde se requiere obtener una tensión
nS nP
I = n P  I L − I s e − 1 − 
AKT
(1) media de salida que puede ser mayor o menor que la
   RP  aplicada a su entrada, gobernando los tiempos en que
   
  el interruptor principal del convertidor conduce o no-
conduce (técnica PWM), generalmente a frecuencia
donde nP y nS indican el número de células constante. A la relación entre el intervalo de tiempo
conectadas en paralelo y en serie; RP y RS las en el que el interruptor esta cerrado (TON) y el
resistencias intrínsecas paralelo y serie asociadas a la periodo de conmutación (T), se le denomina ciclo de
célula, K es la constante de Boltzman (1.38·10-23 J/K) servicio (δ).
y q es la carga del electrón. El factor A determina la
desviación de las características de una unión p-n
 XXV Jornadas de Automática ISBN: 84-688-7460-4

Tanto en el modo de funcionamiento con corriente no se alcanzan valores próximos a la corriente de

continuada (CCM, Continuous Conduction Mode), cortocircuito del panel.
como en el modo discontinuo (DCM, Discontinuous
Conduction Mode), las tres topologías básicas pueden
ser asemejadas a un transformador de corriente
continua, donde la relación de transformación puede 3.2 CONVERTIDOR ELEVADOR.
ser controlada electrónicamente variando el ciclo de
servicio δ del convertidor en el rango [0,1]. La figura De la tabla 1 se deduce:
2 muestra el diagrama de un panel solar conectado a  
un convertidor DC/DC, donde la resistencia reflejada  4⋅ R 
Ri (CCM ) ∈ [ 0, R ] ; Ri ( DCM ) ∈  , 2R 
( )
2
a la entrada del convertidor viene representada por Ri  1+ 1+ 4 
(R es la resistencia de carga).  K 

lo cual indica que un convertidor elevador es incapaz

de reflejar a su entrada impedancias mayores que la
de carga, y que por tanto no se alcanzan valores
próximos a la tensión de circuito abierto del panel.

3.3 CONVERTIDOR REDUCTOR-ELEVADOR

Figura 2: Conexión Panel-Convertidor. Las topologías reductora-elevadoras Zeta, Cuk y

SEPIC proporcionan la misma resistencia de entrada
Suponiendo convertidores sin pérdidas, la relación mostrada en la tabla 1; y por tanto que
entre la resistencia de entrada y la carga conectada a
los convertidores se muestra en la Tabla 1, tanto para
Ri (CCM ) ∈ [ 0, ∞ ) y Ri ( DCM ) ∈  2 L , ∞
 T )
CCM, como DCM. lo cual implica que estas estructuras son capaces de
caracterizar toda la curva I-V de un módulo en CCM,
Convertidor Kcrit. Ri (CCM) Ri (DCM)
desde la tensión a circuito abierto hasta la corriente
de cortocircuito.
R R  2
Reductor 1−δ ⋅ 1 + 1 + 4 K 2 
δ2 4  δ 
4⋅ R
δ ⋅ (1 − δ )2 R ⋅ (1 − δ )
2 2
Elevador  
1 + 1 + 4δ K 
2

 
R ⋅ (1 − δ )
2
K ⋅R
Reductor-
(1 − δ )2 a) Convertidor Reductor.
elevador.
δ 2
δ2 d)Resistencia de
2L ; entrada vs. δ en CCM.
Con K= DCM ocurre para K<Kcrit
RT

Tabla 1: Valores de Ri para los convertidores de la

Figura 3.
b) Convertidor Elevador. e) Resistencia de
En la figura 3 se muestran los tres convertidores
básicos que proporcionan las diferentes relaciones de entrada vs. δ en CCM.
conversión dadas en la Tabla 1, junto con una
representación gráfica de la resistencia de entrada
reflejada en función del ciclo de servicio δ para
CCM.

3.1 CONVERTIDOR REDUCTOR. c) Convertidor Red-Elev. f ) Resistencia de

entrada vs. δ en CCM.
De la tabla 1 se deduce: Figura 3: Convertidores DC/DC usados comúnmente
R
(
Ri (CCM ) ∈ [ R, ∞ ) ; Ri ( DCM ) ∈  ⋅ 1 + 1 + 4 K )  y su resistencia de entrada.
2
,∞
4 
lo cual indica que un convertidor reductor presenta 4. DISEÑO DEL SISTEMA Y
una incapacidad para reflejar a su entrada METODOLOGÍA.
impedancias menores que la de carga, y que por tanto
 XXV Jornadas de Automática ISBN: 84-688-7460-4

El análisis teórico realizado En la sección 3, ha sido triangular y un comparador, como se observa en la

verificado mediante PSpice® utilizando un modelo de figura 5. El oscilador superior proporciona una señal
panel implementado mediante bloques ABM (Analog triangular de unos 25Khz. El inferior genera otra
Behavioral Modeling) y basado en la ec. (1). Los señal de 25 Hz que podríamos denominar de
parámetros del modelo corresponden al módulo “BP referencia y de amplitud superior a la primera, para
Saturno” (n s= 60 y np = 1) [8]. garantizar los estados δ = 0 y δ = 1. La relación de
frecuencias indica de forma aproximada el nº de
PANEL FOTOVOLTAICO
PARAMETERS: puntos por barrido (aproximadamente 103 en este
Rv ar = 1k

R2
caso).
out

I1 R1
0.768

R4
En las figuras 6 y 7, se presentan los resultados
21
2.904 511.72 {Rv ar} obtenidos en la simulación de las topologías
Elevadora y SEPIC, para las mismas condiciones de
(2.26e-11)*(exp(V(%IN1,%IN3)/((55.82)*(8.625e-5)*(273.15+V(%IN2))))-1)
radiación, temperatura y carga (R=20Ω). En estas
figuras se observan las predicciones del análisis
0 teórico realizado en la sección 3, en cuanto a la
4.0A
imposibilidad de reproducir los puntos cercanos a
I
VOC por la topología Elevadora, y la capacidad de la
2.0A
estructura Reductora-Elevadora para representar
todos los puntos de la característica I-V de un panel,
0A incluido el VOC y el ISC. En la Tabla 2 se muestran los
100W valores de los puntos más representativos de dichas
P curvas. Nótese la precisión de la metodología
50W
empleada, la cual permite obtener resultados muy
próximos a los deseados.
0W
0V 5V 10V 15V 20V 25V 30V 35V 40V

V(OUT)
R4
Figura 4: Modelo y curvas I-V y P-V del modulo “BP 9.5k
C1 HI
LO

R5
Saturno”. 100k
5n
15Vdc
V1
-15Vdc
V2

8
HI
V+ 0 0
La figura 4 presenta las curvas I-V y P-V del módulo 3
+ R1
LO
4

en cuestión a 21ºC y 1000W/m2. Para obtener un TL082

OUT
V-
1
15k
2
-
V-

2 TL082 1
conjunto de valores (V, I) a la salida del módulo, es -
4 3
+
OUT
V+
R52

R6 50k LO
necesario implementar una Ri variable ya que ésta 100k
LO

100k
R2
8
Vtriang 4

R53 V-
fija el punto de operación del módulo. Según lo 0
0
HI
100k 2
-

explicado en la sección 3, esto se puede conseguir R41

0
TL082OUT
V+
1

utilizando un convertidor DC/DC con δ variando

3 PWM
C11 +
10k 8
R51
10u Vref HI
entre 0 y 1. 100k

HI 8
LO 4
3 V+
4.1 OBTENCIÓN DE UNA SEÑAL PWM +
TL082 1
R11
2
V-

CON δ VARIABLE.
OUT -
10k
2 V- TL082 1
- OUT
V+
3
R61 4 +
LO 8
R21
100k
La metodología de elaboración de curvas I-V y P-V 0
100k HI

0
propuesta, se implementa en base a los montajes de
las figuras 5, 6 y 7. La generación de la señal PWM .
se consigue básicamente comparando una señal 1.0V
triangular Vtriang con otra señal de referencia Vref. El -1.0V
comparador satura a nivel alto o bajo dependiendo de 0s 1.0ms 2.0ms
VTriang
los valores instantáneos de ambas señales. Time

Manteniendo Vref constante se obtiene una señal

PWM con δ constante. Variando lentamente Vref -2.0V
0V

obtendremos una señal PWM con δ variando de 0s

Vref
20ms 40ms
Time
forma suave. La frecuencia de la señal PWM viene
fijada por la señal Vtriang (señal más rápida), mientras 10V

que δ depende de la relación entre la amplitud de

0s 400us 800us
Vtriang y el valor de Vref en ese instante. VPWM
Time

En nuestro caso, el sistema de generación PWM-δ Figura: 5. Sistema de generación de la señal PWM
variable, se compone de dos generadores de señal con δ variable.
 XXV Jornadas de Automática ISBN: 84-688-7460-4

adquisición de datos (SAD) tal como el mostrado en

la figura 8. Además, añadiendo al SAD las señales
procedentes de las medidas de temperatura y
radiación luminosa del panel, es posible efectuar de
forma automática un análisis de fallos en el módulo o
detección de anomalías en su funcionamiento, debido
a depósitos de polvo en su superficie, por ejemplo.
Las señales adquiridas de tensión y corriente
(mediante un sensor de efecto Hall para no interferir
en el sistema) se acondicionan antes de ser
introducidas en una tarjeta de adquisición de datos
(TAD), con objeto de poder utilizar todo el margen
de entrada del canal analógico. La TAD puede ser
una de tipo estándar alojada en un slot de expansión
de un PC que actúa a modo de controlador de toda la
instalación. Las señales de tensión y corriente
adquiridas junto con el control del ciclo de trabajo
del convertidor DC/DC se emplean para, en base a la
Figura 6: Curvas I-V y P-V obtenidas con un metodología presentada en el apartado 4,
convertidor elevador. implementar las característica I-V y P-V reales del
panel.

PANEL FOTOVOLTAICO

R2 L2 C2 D2
out 1 2
2
100uH MBR3520
0.768 330uF R2
R1 M1 L1 C1
20
I1 21 R1
100uH 330uF
511.72 10
2.904 IRF150

(2.26e-11)*(exp(V(%IN1,%IN3)/((55.82)*(8.625e-5)*(273.15+V(%IN2))))-1)
GENERADOR 0
PWM
δ-VARIABLE
0

Figura 8:Esquema del circuito de monitorización y

análisis de fallos

Además de las señales de tensión y corriente, el SAD

adquiere los valores de temperatura y radiación
presente en el módulo, lo cual permite implementar
Figura 7: Curvas I-V y P-V obtenidas con un las características I-V y P-V teóricas del panel.
convertidor SEPIC. Una comparativa entre las característica I-V y P-V
reales y teóricas permite analizar el comportamiento
correcto del panel y detectar cualquier fallo o
anomalía en su funcionamiento, bien por problemas
VOC (v) ISC (a) VMPP (v) IMPP (a) PMPP (w) en las propias células o por la acción de agentes
Real 36.18 2.90 29.80 2.71 80.66 externos.
Elevad. - 2.90 29.80 2.71 80.66
SEPIC 36.19 2.90 29.81 2.71 80.66
El control de toda la instalación se efectúa mediante
Tabla 2: Comparativa de los valores más un PC y software diseñados para ello. Para la
significativos de las curvas I-V y P-V obtenidas con implementación del software de instrumentación y
distintos convertidores control se puede diseñar una aplicación específica en
un lenguaje de programación como C, por ejemplo, o
bien utilizar programas comerciales de los usuales en
5. MONITORIZACIÓN Y ANALISIS los sistemas de adquisición y distribución de señales.
DE FALLOS En nuestro caso hemos optado por lo segundo y el
software del instrumentación y control se ha
El sistema de generación automática de curvas I-V y implementado en Labview®.
P-V diseñado en el apartado anterior permite una
implementación práctica inmediata. Para ello se
puede utilizar un simple y económico sistema de
XXV Jornadas de Automática ISBN: 84-688-7460-4

6. CONCLUSIONES. Power Electronics, vol. 17, no. 6, November

2002.
En el presente trabajo se ha propuesto y demostrado
mediante simulación la posibilidad de caracterizar las [5] K.K. Tse, M.T. Ho, H.S.-H. Chung, [Link]:
curvas I-V y P-V de un panel fotovoltaico mediante "A novel maximum power point tracker for PV
el uso de convertidores CC/CC con δ variable. Esta panels using switching-frequency modulation
metodología abre una nueva alternativa en el diseño scheme". IEEE Trans. on Industrial Elec., vol.
de trazadores I-V fotovoltaicos, pudiendo tener 51, no. 2, April 2004.
consecuencias importantes en el proceso de
simplificación de éstos. Además, como se ha [6] S. Singer: "The application of loss-free resistors
demostrado, la topología Elevadora-Reductora es la in power processing circuits". Power Elect.,
única eficaz para este fin, ya que es capaz de barrer la IEEE Trans., vol. 6, no. 4, pp. 595-600, 1991.
curva completa.
Respecto de los métodos tradicionales, el propuesto [7] C. Jingquan, R. Erickson, D. Maksimovic:
en este trabajo proporciona las ventajas siguientes: "Averaged switch modeling of boundary
1) Mínima pérdida de potencia al contrario de los conduction mode DC-to-DC converters".
sistemas que operan en zona lineal (zona Industrial Electronics Society, IECON '01. The
activa). Esto implica ventajas en tamaño y 27th Annual Conference of the IEEE, vol. 2, pp.
coste. 844 – 849, 29 Nov.-2 Dec. 2001.
2) Adaptación automática del intervalo de
interpolación. [8] CIEMAT: "Fundamentos, dimensionado y
La metodología propuesta, una vez verificada aplicaciones de la energía solar fotovoltaica",
mediante simulación, se ha implementad en un Cap. 11, Editorial CIEMAT, Madrid 2000.
sistema real, el cual permite además, mediante la
adquisición de los valores de radiación luminosa y [9] D,L, King: "Photovoltaic module and array
temperatura en el panel, monitorizarlo y detectar performance characterization methods for all
cualquier fallo o anomalía en su funcionamiento. system operating conditions". Proceedings of
En la actualidad estamos trabajando en la NREL/SNL Photovoltaic Program Review
estandarización de la metodología y en el Meeting, pp. 347-368, November 1996.
dimensionado del sistema de instrumentación y
control, con objeto de adaptarlo a las necesidades de [10] C.M. Whitaker, T.U. Townsend, J.D.
los sistemas comerciales. Newmiller, D.L. King, W.E. Boyson, J.A.
Kratochvil, D.E. Collier, D.E. Osborn:
Referencias "Application and validation of a new PV
performance characterization method".
[1] C. Hua, C. Shen: "Study of maximum power Photovoltaic Specialists Conference,
tracking techniques and control of DC/DC Conference Record of the Twenty-Sixth IEEE,
converters for photovoltaic power system". 29th pp. 1253-1256, 29 Sept.-3 Oct. 1997.
Annual IEEE PESC, IEEE Computer Soc.
Press, pp. 86-93, New York 1998. [11] M. A. de Blas, J. L. Torres, E. Prieto, A. García:
"Selecting a suitable model for characterizing
[2] M.A.S. Masoum, H. Dehbonei, E.F. Fuchs: photovoltaic devices". Renewable Energy, vol.
"Theoretical and experimental analyses of 25, no. 3, pp. 371-380, March 2002.
photovoltaic systems with voltage-and-current-
based maximum power-point tracking". IEEE [12] B. G. Bahgat, N. H. Helwa, G. E. Ahamd, E. T.
Transactions on Energy Conversion, vol. 17, no. El Shenawy: "Estimation of the maximum
4, December 2002. power and normal operating power of a
photovoltaic module by neural networks".
[3] K.H. Hussein, I. Muta, T. Hoshino, M. Osaka: Renewable Energy, vol. 29, no. 3, pp. 443-457,
"Maximum photovoltaic power tracking: an March 2004.
algorithm for rapidly changing atmospheric
conditions". IEE Proc. Gener. Transnm. [13] J.A. Gow, C.D. Manning: "Development of a
Distrib., vol. 142, no.1, pp. 59-64, January photovoltaic array model for use in power-
1995. electronics simulation studies". IEE Proc. on
Electric Power Applications, vol. 146, no.2, pp.
[4] K.K. Tse, M.T. Ho, H.S.-H. Chung, [Link]: 193-200, March 1999.
"A comparative study of maximum-power-point
tracker for photovoltaic panels using switching
frequency modulation". IEEE Transactions on