RESUMEN

La celda y el módulo fotovoltaico se representan típicamente mediante un circuito equivalente. La correcta modelación del circuito equivalente es hoy en día una herramienta esencial para diseñadores y proyectistas, para la realización de estudios y optimizar la energía en los sistemas fotovoltaicos. Sin embargo, la determinación de los parámetros del circuito equivalente, es aún una tarea difícil y normalmente requiere de datos experimentales y procedimientos de cálculo, generalmente no disponibles por el usuario. Este artículo presenta una nueva metodología, que posibilita el cálculo exacto de los parámetros y magnitudes eléctricas, del modelo de simple diodo con cinco parámetros, obtenidos a partir de los datos aportados por los fabricantes en las hojas de datos de los módulos fotovoltaicos. Finalmente se comprueba la validez y precisión de la metodología propuesta, comparando datos y curvas experimentales obtenidas por el fabricante de un módulo comercial, con los obtenidos mediante el modelo matemático propuesto.

Palabras clave: fotovoltaico, simulación, arreglo, modelo, celdas solares.

ABSTRACT

The cell and photovoltaic module is typically represented by an equivalent circuit. The correct equivalent circuit modeling is nowadays an essential tool for designers and planners, for studies and optimizes energy in photovoltaic systems. However, the extraction of the equivalent circuit parameters is still a difficult task and typically requires experimental data and calculation procedures, generally not available for the user. This paper presents a new methodology, which enables the exact calculation of the parameters and electrical quantities, simple diode model with five parameters, obtained from the data provided by manufacturers in the data sheets of photovoltaic modules. Finally is showed the validity and accuracy of the proposed methodology, comparing experimental data and curves obtained by the manufacturer of a commercial module, with those obtained by the proposed mathematical model is checked.

Key words: photovoltaic, simulation, array, model, solar cells.

INTRODUCCIÓN

La energía solar constituye una fuente de energía limpia e inagotable. La conversión directa de esta energía en electricidad se realiza a través de las celdas fotovoltaicas. Las celdas se agrupan para formar los módulos fotovoltaicos, estos a su vez se agrupan en serie y/o paralelo para formar arreglos fotovoltaicos, conformándose así los denominados generadores fotovoltaicos. Los generadores fotovoltaicos pueden alimentar directamente cargas eléctricas de corriente directa. Aplicaciones más sofisticadas requieren de convertidores electrónicos para procesar la energía eléctrica proveniente del generador fotovoltaico. Estos convertidores pueden ser usados para regular la tensión y la corriente en las cargas, para controlar el flujo de potencia en los sistemas conectados a las redes eléctricas y principalmente para seguir el punto de máxima potencia (en inglés MPPT), del generador fotovoltaico para lograr la mejor eficiencia del sistema.

La rápida expansión de la energía solar fotovoltaica, como consecuencia de las ventajas que representa su uso, esencialmente en los sistemas interconectados a las redes eléctricas, por su carácter totalmente descentralizado, tanto en grandes parques fotovoltaicos o en cubiertas y fachadas de edificios, han hecho posible que se amplíen cada vez más las investigaciones en este campo [1]. Los dispositivos fotovoltaicos iluminados sobre la superficie de la tierra dan una respuesta que depende de la calidad e intensidad de la radiación recibida y de la temperatura de las celdas [2]. La respuesta obtenida es el resultado en su conjunto de los parámetros eléctricos de los módulos y de la estructura del generador fotovoltaico. Por tanto los diseñadores seleccionan diferentes opciones de módulos fotovoltaicos para obtener el mejor, antes de dimensionar el generador fotovoltaico y por otra parte estiman la dimensión del sistema requerido para suministrar la energía necesaria [3-4]. Por tanto, la predicción de la respuesta eléctrica del generador fotovoltaico en cualquier condición de operación es vital para predecir el comportamiento del sistema fotovoltaico desde su etapa de diseño y para lograr los mejores indices energéticos en explotación.

El comportamiento eléctrico de un dispositivo fotovoltaico, bajo determinadas condiciones de operación, se caracteriza por sus parámetros eléctricos y por la curva Volt-Ampere (I=f(V)), que describe en operación. Numerosas investigaciones se han desarrollado encaminadas a lograr una elevada precisión de los modelos. En todos los casos los parámetros de los modelos desarrollados, se obtienen a partir de los escasos datos suministrados por los fabricantes. Sin embargo, para el cálculo de los parámetros se realizan consideraciones como: se estima el factor de idealidad del diodo [5-8], se considera la fotocorriente igual a la corriente de cortocircuito [9-11], se desprecia la corriente de saturación inversa para el cálculo de la fotocorriente [12-14], se desprecia la resistencia paralelo para calcular la corriente de saturación inversa [15-17], además de las consideraciones anteriores en [18-19], los algoritmos propuestos para el cálculo de los parámetros requieren  de una elevada carga computacional.   

En el presente artículo se presenta una metodología que permite calcular de manera exacta los parámetros del circuito equivalente de simple diodo con cinco parámetros, a partir de los datos aportados por el fabricante, logrando desacoplar en las ecuaciones propuestas los parámetros del circuito equivalente y con un mínimo de esfuerzo computacional calcular cada uno de ellos.

MATERIALES Y MÉTODOS

Circuito equivalente del módulo fotovoltaico

Desde el siglo XX y lo que va del XXI, se han desarrollado diferentes circuitos equivalentes para celdas y módulos fotovoltaicos. Tipicamente la opción de uno u otro circuito equivalente depende de la precisión deseada de la curva I=f(V) modelada con respecto a la curva obtenida de forma experimental. La figura 1, muestra el circuito equivalente más usado para modelar el módulo fotovoltaico, por garantizar elevada precisión y desempeño en estudios de sistemas fotovoltaicos.

Aplicando las leyes de Kirchoff  y el modelo matemático del diodo Shockley para representar la corriente I_D del diodo [12], la característica Volt-Ampere del circuito I=f(V), se expresa a través de la ecuación (1).

(1)

Donde, I₀ es la corriente de saturación inversa del diodo, N_c es el número de celdas del módulo fotovoltaico, n es el factor de idealidad del diodo, desde la teoría de difusión de Shockley. El valor de n para celdas de siliciose encuentra en valores comprendidos en el intervalo 1≤n≤2 y su elección tiene una marcada incidencia en la forma que adopta la característica I_m=f(V_m) [12; 19]. La mayoría de los modelos adoptan para este factor un valor arbitrario dentro del intervalo, recomendándose por varios autores elegir un valor por el entorno de 1,3 [12; 14]. Como se verá posteriormente en la metodología propuesta este valor se calcula de manera automática garantizando el mejor ajuste de la característica I_m=f(V_m), a partir de los datos del fabricante. Finalmente V_t, es la tensión térmica del diodo, una constante que depende de la temperatura T (K), la constante de Boltzman (K=1,3806503 10^-23 J/K ) y la carga del electrón (q=1,.60217646  10^-19 C). La dependencia de V_t, de estos parámetros se expresa como V_t=K T/q.

La ecuación (1), resultado del circuito equivalente propuesto en la figura 1, contiene cinco parámetros desconocidos, I_phm, R_sm, R_pm, I₀ y n, cuyos valores se determinan a partir de los datos aportados por el fabricante en la hoja de datos de sus módulos fotovoltaicos, que básicamente representan tres puntos de la característica: 1) el punto de cortocircuito (0,I_ccm),  2) el de circuito abierto (V_ocm,0) y 3) el de máxima potencia (V_max,I_max) .Los datos de estos tres puntos corresponden a las condiciones estándar de 1000W/m², una temperatura de 25 ^ᵒC y una masa de aire (AM=1.5).

La figura 2.a, muestra la característica externa del módulo fotovoltaico I_m=f(V_m), en la que aparecen señalados los puntos antes mencionados, mientras en la figura 2.b, aparece la característica de potencia del módulo fotovoltaico en función de la tensión, P_m=f(V_m).

Es poco probable el trabajo del módulo fotovoltaico en las condiciones estándar de operación, por lo que dos parámetros de la ecuación (1), son dependientes de la variación de la radiación y temperatura. La fotocorriente posee una elevada dependencia de la radiación (G) y en menor medida de la temperatura y por su parte la corriente de saturación inversa (I₀) posee una elevada dependencia de la temperatura. A efecto de los cálculos para condiciones de operación no estándar, el valor de estas magnitudes eléctricas se obtiene a partir de las ecuaciones (2) y (3), [5; 8; 11; 19].

(2)

(3)

Donde I_phmn e I_0n  representan la fotocorriente y la corriente de saturación inversa en condiciones estándar respectivamente, ΔT=T_x-T_n (T_x y T_nson las temperaturas a la condición real y la estándar respectivamente en K, K_i (A/K), representa el factor de variación de la fotocorriente por grados Kelvin, G_xy G_nom representan las respectivas radiaciones en la condición real y estándar. Por su parte E_g es la energía de activación de la banda prohibida del semiconductor, para el silicio poli-cristalino E_g=1.12eV y para el silicio amorfo E_g=1.7eV [5].

La tensión de circuito abierto V_ocmes otro de los parámetros que sufre una marcada variación con la temperatura [5; 8;11;19], la ecuación (4), para su cálculo en condiciones no estándar es la siguiente:

(4)

Donde V_ocmn y V_ocmx en (V), representan las correspondientes tensiones de circuito abierto en condiciones nominales y en la condición no estándar, K_v (V/K), representa el factor de variación de la tensión de circuito abierto con el incremento de la temperatura, su valor es negativo, el valor de V_ocmx disminuye con el incremento de la temperatura.

Método de obtención de I₀ a partir de los datos del fabricante

A partir del circuito equivalente del módulo (fig.1), y tomando en cuenta los datos nominales del fabricante, se plantean los tres circuitos representado en la figura 3, para facilitar la comprensión de los cálculos y los posteriores análisis.

En la tabla 1, aparecen los valores que toman las magnitudes eléctricas para cada una de las condiciones planteadas en función de los datos aportados por el fabricante.

(5)

Aplicando el teorema a las combinaciones posibles de los circuitos de la figura 3, sustituyendo los valores de la tabla 1, según corresponda, desarrollando manipulaciones algebraicas y despejando en cada uno de los casos I₀, se obtienen tres expresiones que permiten calcular el valor de I₀, para las condiciones estándar en función de los datos aportados por el fabricante en la hoja de datos de los módulos. La ecuación (6), corresponde con la aplicación del teorema a la condición 1 y 3, la ecuación (7), con la aplicación del teorema a la condición 1 y 2, y la ecuación (8), con la aplicación del teorema a la condición 2 y 3.

(6)

(7)

(8)

Método de obtención de la fotocorriente (I_phm), a partir de los datos del fabricante

Los circuitos de las figuras 3.a, 3.b y 3.c se corresponden con las condiciones del circuito equivalente del módulo fotovoltaico bajo las condiciones estándar de los datos aportado por los fabricantes. Como consecuencia de ello, la fotocorriente (I_phm), en los tres circuitos es la misma. Las tres ecuaciones anteriores (6; 7; 8), se obtuvieron bajo esta condición, por lo que igualando dos de ellas y despejando I_phm, se obtiene la ecuación (9), donde la fotocorriente solo queda en función de dos parámetros de modelo, I_phm=f(n,R_sm).

(9)

Donde E_Vocm, E_Iccm y E_Pmax se determinan por las ecuaciones (10), (11) y (12):

(10)

(11)

(12)

Método de obtención de la resistencia paralelo del módulo (R_pm), a partir de los datos del fabricante

El valor de la resistencia paralelo (R_pm), se obtiene evaluando en la ecuación (1), obtenida del circuito equivalente del módulo, los valores correspondientes con los datos del fabricante en el punto de máxima potencia (V_max,I_max), el que constituye un punto contenido en la característica Volt-Ampere. Después de realizarse la sustitución de los valores del punto de máxima potencia en la ecuación y despejando R_pm, se obtiene la ecuación (13), para su cálculo, quedando expresada en función de los restantes cuatro parámetros del circuito equivalente, R_pm=f(I_phm,I₀,n,R_sm).

(13)

Consideraciones sobre la potencia del módulo fotovoltaico

La potencia de salida del módulo se calcula a través de la ecuación (14), (ver fig.1).

(14)

(15)

De la ecuación (15), se obtiene la conductancia del módulo fotovoltaico en el punto de máxima potencia mediante la ecuación (16).

(16)

Método de obtención de la pendiente de la característica Volt-Ampere del módulo fotovoltaico

La pendiente de la característica Volt-Ampere del módulo fotovoltaico, expresa la conductancia del módulo, se obtiene diferenciando la ecuación (1), respecto a la tensión del módulo, ecuación (17).

(17)

(18)

Para el caso particular en que V_m=V_max e I_m=I_max, cuando se cumpla que dI_m/dV_m= -(I_max/V_max), se garantiza la pendiente que debe poseer la característica Volt-Ampere del módulo fotovoltaico en el punto de máxima potencia. Para este caso particular dI_m/dV_m (V_max,I_max)=f(I₀,n,R_sm,R_pm), cuatro de los parámetros del circuito equivalente.

Procedimiento de cálculo de I_phm, I₀, n, R_pm y R_sm

El procedimiento de cálculo se realiza explorando el intervalo en que puede variar n. Para cada valor de n se asume inicialmente el valor de R_sm igual a cero (R_sm=0), y se calculan en el siguiente orden los restantes tres parámetros I_phm, I₀ y R_pm. Para cada valor de n se incrementa R_sm en un valor ΔR_sm, se va obteniendo el valor de la pendiente y se va calculando el error absoluto, mediante la ecuación (19), hasta lograr que el error sea menor a un error mínimo (E_min), definido por el usuario, cumpliendo esta condición quedan calculado para el valor n, los restantes cuatro parámetros del circuito equivalente.

(19)

Así se calculan de manera sucesiva los valores de los parámetros del circuito equivalente, para cada nuevo valor de n, explorando el intervalo con un incremento de Δn, cada vez que se realiza un cálculo se compara el resultado del error absoluto, con el obtenido anteriormente, con el objetivo de que en caso que sea menor desechar los resultados anteriores y admitir como más precisos los valores recién calculados. De manera que al final de la exploración realizada se obtiene el mejor ajuste de la curva. Como el valor de n influye marcadamente en la forma de la curva, cada vez que se realiza un cálculo se tiene una forma diferente de la característica Volt-Ampere, y finalmente se obtiene el resultado de los parámetros de la curva más precisa. La figura 4, muestra el algoritmo de cálculo implementado.

Expresión de cálculo de I_0x

(20)

Validación de la metodología

El algoritmo de cálculo de los cincos parámetros del modelo se implementó en Matlab. El valor de Δn  y ΔR_sm, se asumió igual a 0,01 y el valor de E_min se asumió igual a 10^-5. Los valores asumidos garantizan una elevada precisión de los cálculos.  La metodología se comprobó con más de cien módulos fotovoltaicos de diferentes fabricantes y en todos los casos se obtuvieron excelentes resultados. Para ilustrar la validez de la metodología se comparan los resultados obtenidos mediante el modelo con los de la hoja de datos del módulo KC200GT, del fabricante Kyocera [23].

La tabla 2, muestra los datos ofrecidos por el fabricante en su hoja de datos y los obtenidos mediante el modelo propuesto, donde se aprecia que los valores coinciden con los del fabricante.

Además, en la figura 5, se muestran, curvas de corriente en función de la tensión I_m=f(V_m), obtenidas a diferentes irradiancias y temperaturas con los parámetros calculados  mediante la metodología desarrollada, así como se adicionan puntos marcados en la  figura  mediante círculos, que representan puntos experimentales, extraídos de las curvas de la hoja de datos del fabricante [23]. Como se puede observar, para todos  los casos, se obtienen magníficos resultados, los puntos experimentales se encuentran sobre las curvas obtenidas mediante el modelo propuesto.

CONCLUSIONES

Se desarrolló y validó una metodología que permite obtener de forma exacta los parámetros del circuito equivalente del módulo fotovoltaico a partir de los datos aportados por el fabricante en la hoja de datos de los módulos fotovoltaicos.

La expresión propuesta para obtener la corriente de saturación inversa a temperatura no estándar (I_0x), se obtiene a partir del propio circuito equivalente del modelo propuesto y presenta la ventaja que permite obtener una mejor interrelación de las magnitudes eléctricas del circuito. Los datos requeridos para su cálculo, se obtienen a partir de los datos aportados por el fabricante en la hoja de datos de los módulos fotovoltaicos y el factor de idealidad calculado.

El artículo provee al lector de toda la información necesaria para obtener los cinco parámetros del modelo de simple diodo de manera exacta, así como de un conjunto de expresiones de cálculo de gran utilidad para proponer modelos de simulación dinámica del módulo y/o arreglo fotovoltaico.

REFERENCIAS

Recibido: 1/10/2016