Simulación de Una Planta Fotovoltaica en Atp-Emtp PDF
Simulación de Una Planta Fotovoltaica en Atp-Emtp PDF
Simulación de Una Planta Fotovoltaica en Atp-Emtp PDF
Captulo 1
Introduccin
El futuro est en la energa. Este sector est llamado a ser el nuevo motor de la economa
mundial, igualmente de la economa espaola. La Organizacin Internacional del Trabajo (OIT)
apunta a que en el horizonte de 2030 habr en todo el mundo 20 millones de puestos de trabajo
relacionados con las energas renovables.
Se calcula que para hacer frente al cambio climtico y minimizar sus consecuencias,
debemos conseguir una reduccin del 80 % de las emisiones de gases de efecto invernadero para el
2050 sobre niveles de 1990. Para conseguirlo la comunidad internacional ha fijado como Ley
Internacional de cumplimiento obligatorio, desde febrero de 2005, del Protocolo de Kyoto. Cabe
destacar, sin embargo, todava hoy la existencia de puntos de controversia en dicha ley entre los
diferentes pases, como la exclusin de la energa nuclear como una fuente de energa til en la
reduccin de emisiones de gases contaminantes.
Las consecuencias del cambio climtico que ya son visibles hoy en da incluyen (fuente:
Greenpeace):
Los siete aos ms calurosos de los ltimos 130 se han registrado en los ltimos 11 aos.
El nivel del mar se ha elevado 10-20 centmetros en los ltimos 100 aos, 9-12 cm de los
cuales en los ltimos cincuenta.
La energa nuclear es costosa, sus residuos son altamente contaminantes a largo plazo y
no cuenta con el visto bueno de la opinin pblica, aunque parece que cada da tiene menos
oposicin. Sin embargo, no emite gases contaminantes y, una vez construida la central que la
produce, el coste de produccin no resulta tan caro como otras fuentes de energa.
Sus detractores recuerdan los problemas de seguridad, la futura escasez de uranio, la difcil
gestin de los residuos radiactivos y aseguran que son un freno para las energas renovables.
Adems, sealan el alto coste para construir una central nuclear y para almacenar los residuos.
Francia ha apostado por la energa nuclear y se ha afanado en firmar acuerdos con otros
pases, sobre todo del mundo rabe. El sector atmico francs vende as centrales, uranio de
combustible y trata los residuos. Francia cuenta con 58 reactores nucleares y est construyendo
otro. Con ellos genera el 78 % de su electricidad, lo que le convierte en el pas del mundo con
mayor porcentaje de este tipo de energa.
Por otra parte, reducir el consumo de energa, a travs del ahorro y la eficiencia, es tan
necesario como sustituir las fuentes de energa contaminantes por limpias y renovables.
En 2008, el precio del petrleo alcanz los 140 dlares el barril antes de hundirse con el
resto de la economa, lo que indica los riesgos de supeditar el futuro a alto tan poco predecible
como el petrleo.
Figura. Camino ejemplar, energa primaria global. Fuente: German Advisory Council on Global
Change, www.wbgu.de.
Los proyectos puestos en servicio en Espaa hasta el ao 2010 supondrn una inversin de
23.598 millones de euros, y tan slo el 2.9 % de los fondos invertidos correspondern a fondos
pblicos (680.9 millones de euros).
Existen sin embargo ciertas barreras que dificultan la instalacin de energas renovables en
Espaa. Entre ellas encontramos barreras financieras: debido a que resulta imprescindible la
financiacin de los proyectos a travs de entidades financieras, es necesario demostrar una
rentabilidad atractiva. Para ello, es fundamental un marco normativo estable, adems de establecer
mecanismos basados en primas a la produccin.
Cuando se anunci en el verano de 2007, la prima en Espaa para la energa solar era la
ms generosa del mundo (58 cntimos por kilovatio-hora1). Con el tiempo, se ha visto que era
demasiado alta (se construyeron plantas de baja calidad y mal diseadas). Ahora se ajusta
trimestralmente: es de unos 39 cntimos por kilovatio-hora de electricidad de plantas de energa
solar independientes, y ligeramente superior para los paneles en tejados.
Otro de los factores que suponen una barrera a la promocin de las energas renovables en
Espaa, es el tratamiento de acceso de conexin a la red.
El presupuesto del gobierno de Barack Obama para 2010 prev duplicar en tres aos la
capacidad nacional de produccin de energas renovables. Las turbinas elicas y los
biocombustibles sern una aportacin importante. Pero ninguna forma de energa es ms
abundante que la del sol.
1
Un kilovatio-hora har brillar una bombilla de 100 vatios durante 10 horas.
Cada hora, el sol inunda la Tierra con un torrente de energa trmica equivalente a 21.000
millones de toneladas de carbn.
La energa total que necesitamos en todo el mundo es de unos 16 teravatios (un teravatio
equivale a un billn de vatios). Para 2020 necesitaremos 20 teravatios, y el sol derrama 120.000
teravatios sobre las tierras emergidas del planeta. O, dicho de otro modo, una energa 6.000 veces
superior al consumo mundial de electricidad se derrama constantemente sobre la Tierra. Desde
este punto de vista, podemos decir que la energa solar es virtualmente ilimitada.
Bastara cubrir con paneles fotovoltaicos tres milsimas de la superficie de Estados Unidos
(un cuadrado de 160 kilmetros de lado) para suministrar electricidad a todo el pas.
Los expertos predicen que, posiblemente el ao que viene, Italia ser el primer lugar donde
la electricidad de origen solar no necesitar subvenciones para competir con la electricidad
producida con combustibles fsiles.
Captulo 2
El Sol es una fuente de energa de 3,81020 MW, lo que equivale a una densidad de energa
por cada metro cuadrado de superficie solar de 62,5 MW. Sin embargo, como consecuencia de la
distancia que nos separa, solo una pequea parte, aproximadamente 1 kW, llega a la superficie de
la Tierra.
Llegados a este punto hay que decir que la potencia radiante que alcanza a nuestro planeta
es atenuada por diferentes fenmenos en la atmsfera. Estos fenmenos son:
Absorcin. Debida a las molculas que componen la misma: O3, H2O, O2, CO2,
etc.
Adems se ha de tener en cuenta que la energa emitida por el Sol no llega a la Tierra de
manera uniforme. Vara en funcin de la hora del da, la inclinacin estacional del globo terrqueo
respecto al Sol, la zona de la superficie terrestre, etc.
La intensidad de la radiacin solar sobre la superficie de la tierra est influenciada por tanto
por la forma de la Tierra, los movimientos de la misma, fenmenos atmosfricos o la propia
atmsfera, el ciclo da/noche e, incluso, la actividad humana.
Si se suma toda la radiacin global que incide sobre un lugar determinado en un periodo de
tiempo definido (hora, da, mes o ao), se obtiene la energa en kWh/m2. Se obtiene as un valor
diferente para cada regin.
En funcin de cmo reciben la radiacin solar los objetos situados en la superficie terrestre
se pueden distinguir los siguientes tipos de radiacin:
Directa. Es aquella que llega directamente del Sol sin haber sufrido cambio alguno
en su direccin. Se caracteriza por proyectar una sombra definida de los objetos
opacos que la interceptan.
Difusa. Parte de la radiacin que atraviesa la atmsfera es reflejada por las nubes o
absorbida por stas. Esta radiacin se orienta en todas las direcciones, como
consecuencia de las absorciones y reflexiones sucesivas, no solo de las nubes, sino
de las partculas de polvo atmosfrico, montaas, rboles, edificios, el propio suelo,
etc. Se caracteriza por no producir sombra alguna respecto a los objetos opacos
interpuestos.
Reflejada. La radiacin reflejada es, como su nombre indica, aquella reflejada por
la superficie terrestre. La cantidad de radiacin depende el coeficiente de reflexin
de la superficie, tambin llamado albedo. Las superficies horizontales no reciben
ninguna radiacin reflejada, porque no ven ninguna superficie terrestre y las
superficies verticales son las que ms radiacin reflejada reciben.
A la hora de estudiar la inclinacin ms adecuada con la que se deben orientar los mdulos
fotovoltaicos, es necesario precisar la posicin del Sol en cada instante, optimizando as el
rendimiento de la captacin. Para definir cada una de estas posiciones se recurre al sistema de
coordenadas polares.
La situacin del sol en un lugar cualquiera viene determinada por la altura y el azimut del
Sol.
Equipos:
Cada sistema tiene sus ventajas. En este momento, la generacin de vapor, tambin
conocida como energa solar termoelctrica, es ms eficiente que la fotovoltaica, ya que transforma
en electricidad un mayor porcentaje de la luz solar recibida. Sin embargo, requiere muchas
hectreas de terreno y largos tendidos de cables de alta tensin para poner la energa en el
mercado. Los paneles fotovoltaicos, en cambio, se pueden instalar en los tejados, en el lugar donde
se necesita la energa. Los dos procedimientos comparten un inconveniente obvio: su produccin
disminuye cuando el cielo est nublado y desaparece por la noche. Pero los ingenieros ya estn
trabajando para desarrollar sistemas que permitan almacenar la energa y utilizarla en las horas de
oscuridad.
Figura. Mayor planta fotovoltaica construida en Asia. Corea del Sur, 2008.
Las centrales elctricas termosolares pueden convertirse en los parques elicos marinos del
desierto. No hay barreras tcnicas, econmicas o de recursos para suministrar el 5 % de las
necesidades elctricas mundiales con slo electricidad solar trmica en el 2040, incluso asumiendo
que la demanda elctrica global se podra doblar para entonces. La industria termosolar podra as
convertirse en un negocio dinmico e innovador de 15.000 millones de euros anuales en 20 aos,
suministrando electricidad a 100 millones de personas que viven en las partes ms soleadas del
mundo.
En muchas regiones del mundo, un kilmetro cuadrado de tierra basta para generar unos
100-200 Gigavatios hora (GWh) de electricidad al ao usando la tecnologa solar termoelctrica.
Esto equivale a la produccin anual de una central trmica convencional de carbn o de gas de 50
MW. A nivel mundial, la explotacin de menos del 1 % del potencial solar trmico total sera
suficiente para estabilizar el clima mundial mediante reducciones masivas de CO2.
En el desarrollo actual se han abierto dos grandes vas para el suministro a gran escala de
electricidad generada con energa solar trmica. Una es el tipo ISCC, operacin hbrida de
coleccin solar combinada con una planta de ciclo combinado alimentada por gas. La otra es la
operacin solar nica, con una turbina de vapor convencional, y cada vez ms con el uso de un
medio como la sal fundida para el almacenamiento. Se elimina el aceite trmico, se genera vapor
directamente (GDV) en el tubo absorbedor, mejorando la eficiencia energtica y reduccin de
campo solar. Existe una planta piloto de este tipo en Almera, Proyecto DISS, que entr en
funcionamiento en 1998.
Figura. Planta DGV (Generacin Directa de Vapor) con agua-vapor. Se elimina el aceite trmico.
Figura. Diagrama de AndaSol: configuracin del campo solar, sistema de almacenamiento y ciclo
de vapor en el proyecto AndaSol-1, en el sur de Espaa.
1
Rgimen especial. Produccin de energa elctrica realizada en instalaciones cuya potencia instalada no supera
los 50 MW, a partir de cogeneracin u otras formas de produccin de electricidad asociadas a actividades no
elctricas, siempre que supongan un alto rendimiento energtico, o en grupos donde se utilicen como fuente de
energa primaria alguna de las energas renovables no consumibles, biomasa o cualquier tipo de biocarburante, o
residuos no renovables o procedentes de los sectores agrcola, ganadero y de servicios, con una potencia instalada
igual o inferior a 25 MW, cuando supongan un alto rendimiento energtico. La produccin en rgimen especial
est acogida a un rgimen econmico singular.
En el difcil equilibrio permanente que hay que tener entre generacin y demanda, el tener
unas energas con una proporcin muy alta de energas menos gestionables que las convencionales
dada las potencias que se han instalado y que se van a instalar en el futuro-, obliga a un gran reto
al operador del sistema. Sobre todo, si tenemos en cuenta nuestra escasa capacidad de
interconexin con Francia. El pas vecino es nuestra unin con un sistema diez veces mayor que el
nuestro y, por tanto, nos da la suficiente estabilidad y garanta de que an, con esta apuesta por las
energas renovables, vamos a seguir teniendo un sistema seguro.
Cabe destacar la creacin, en 1974, del Centro de Estudios de la Energa con el fin de
actuar ordenadamente en el campo de las energas alternativas. Desde entonces este organismo ha
promovido la utilizacin de energa solar, posibilitando as el surgimiento de una industria nacional
en este sector. Espaa es hoy el primer pas europeo productor de clulas y paneles fotovoltaicos,
con el 10 % de la produccin mundial.
Ley 54/1977.
RD 661/2007.
RD 1578/2008.
Rgimen jurdico y econmico para las instalaciones posteriores a la fecha lmite del
mantenimiento de la retribucin del RD 661/2007. Dispone la realizacin de
inspecciones aleatorias anuales con el objeto de verificar, entre otros aspectos, el
cumplimiento de los requisitos tcnicos.
RD 1663/2000.
Conexin a la red elctrica de baja tensin.
Retribucin.
Primera convocatoria:
Tipo I: 340 /MWh para el subtipo I.1 y 320 /MWh para el subtipo I.2
Durante el ao 2012, a la vista de la evolucin tecnolgica del sector y del mercado, y del
funcionamiento del rgimen retributivo, se podr modificar la retribucin de la actividad de
produccin de electricidad mediante tecnologa solar fotovoltaica.
En cuanto a la potencia total acumulada, Espaa (3.300 MW) se alza al segundo puesto
mundial, slo por detrs de Alemania (5.308 MW) y por delante de Japn (2.149 MW) y EEUU
(1.137 MW). Ya muy por detrs se sitan Corea del Sur (352 MW), Italia (350 MW), Francia (146
MW), China (145 MW), Australia (100 MW) e India (90 MW).
En 2008, gracias a las facilidades crediticias concedidas por las entidades financieras, en
Espaa se invirtieron en fotovoltaica ms de 16.000 millones de euros. La inmensa mayora de la
potencia que se conect corresponde a grandes plantas en suelo. El sistema elctrico y las redes de
distribucin pudieron absorber el enorme volumen de instalacin. El comportamiento del
mercado fotovoltaico fue totalmente inesperado.
El sobrecoste de la energa fotovoltaica en nuestro pas rond los 2.000 millones de euros
en 2009. El coste total de la energa del sistema elctrico espaol ascendi a 7.755 millones de
euros en dicho ao.
Por otro lado, el 37 % de las instalaciones tiene seguidores, de los cuales el 24 % son con
seguimiento a dos ejes y el 14 % a un solo eje. Estos sistemas se encuentran incluso en
instalaciones de unos pocos kW de potencia.
Las principales compaas fotovoltaicas espaolas son Isofoton, BP Solar, Atersa, Gamesa,
Solar, Siliken, Grupo Solar, Guascor, Solaria, Instalaciones Pevafersa. Mientras que las principales
compaas proveedoras de inversores son Atersa, Ingeteam y Enertron.
Se estima que con un 1 % de la superficie desrtica del planeta sera suficiente para cubrir
la demanda elctrica mundial con sistemas CST. En Espaa se mantendra dicha proporcin. El
potencial tcnico-econmico de produccin solar termoelctrica en Espaa, asumiendo
emplazamientos apropiados con irradiancias superiores a los 2.000 kW/m2/ao, es de 1.278
TWh/ao, siendo la demanda actual inferior a los 300 TWh/ao. Todo el sur peninsular situado
por debajo del paralelo 40, as como amplias zonas de las cuencas del Ebro y el Duero, presentan
excelentes condiciones para el aprovechamiento de la radiacin solar. Los costes de generacin
solar termoelctrica en Espaa estn todava por encima de los 0.18 /kWh y los costes de
inversin se sitan en torno a los 4.000 /kW. El objetivo espaol para la nueva generacin de
plantas termosolares estara en bajar de los 2.500 /kW y 0.12 /kWh hacia el ao 2020.
En marzo de 2007, fue inaugurada en Europa la primera planta de energa solar cerca de
Sevilla. La planta de 11 MW conocida como PS10, produce electricidad mediante 634 helistatos.
Cada uno de estos helistatos tiene una superficie de 121 metros cuadrados que concentra los
rayos del sol en lo alto de una torre de 114 metros de altura donde est instalado un receptor solar
y un caldern, desde ah, el vapor es enviado a una turbina de vapor situada en un edificio anexo a
la torre.
2
Por sus explicaciones sobre el efecto fotoelctrico y sus numerosas contribuciones a la fsica terica, en 1921
obtuvo el Premio Nobel de Fsica.
Las clulas solares han sido utilizadas frecuentemente para producir electricidad en lugares
donde no llega la red de distribucin elctrica, tanto en reas remotas de la Tierra como del
Espacio, haciendo posible el funcionamiento de todo tipo de dispositivos elctricos como satlites
de comunicaciones, radiotelfonos o bombas de agua. Ensambladas en paneles o mdulos y
dispuestas sobre los tejados de las casas, por medio de un inversor, pueden inyectar la electricidad
generada en la red de distribucin para el consumo, favoreciendo la produccin global de energa
primaria de un pas, de manera limpia y sostenible. Igualmente estos ltimos pueden configurarse
en amplias zonas no urbanas, dando lugar a plantas de generacin elctrica fotovoltaica.
Bifaciales: estas clulas estn constituidas por una doble unin N+-P-P+,
permiendo as aprovechar tanto la radiacin frontal como la reflejada del suelo (la
radiacin del albedo). El rendimiento de estas clulas puede llegar al 30 % siempre
que se tenga especial cuidado con la calidad de la superficie reflejada y las
condiciones mecnicas de colocacin del panel.
Puede situarse en 1954 el inicio del uso de clulas y paneles fotovoltaicos para producir
electricidad, al estar incorporados en los primeros satlites artificiales.
Hoy da, los investigadores del Laboratorio Nacional de Energas Renovables de Golden,
Colorado, Estados Unidos, intentan aprovechar la propiedad de diferentes semiconductores de
captar diferentes colores de los rayos de luz. Superponiendo capas de unos compuestos llamados
fosfuro de indio-galio y arseniuro de galio-indio, y utilizando una lente3 para concentrar la luz
solar, el ao pasado lograron fabricar una clula fotovoltaica con una eficiencia del 40.8 %. Pero la
produccin masiva de estas clulas an es un objetivo lejano: ahora mismo pueden fabricarse a un
precio de unos 7.000 euros por centmetro cuadrado. La tecnologa es de una complejidad
increble.
Otra posibilidad es renunciar a una mayor eficiencia a cambio de un coste menos elevado.
Los semiconductores de capa fina generan menos energa, pero requieren menos materia prima, lo
que los convierte en una alternativa econmica para las grandes instalaciones fotovoltaicas. Existen
ya clulas solares de capa fina a un coste de unos 80 cntimos de euro por vatio.
La cantidad de silicio necesario para producir un Wp se reduce a un ritmo del orden del 5
% anual.
Con la mirada puesta en un futuro ms lejano, los ingenieros del NREL estn trabajando
con lquidos fotovoltaicos. El objetivo es que cuesten lo que un litro de pintura. La eficiencia no
ser del 50 o el 40 %, sino del 10 %. Pero si es barato, la gente puede pintar las paredes de su casa,
conectarse y listo.
3
Lentes de Fresnel. Permiten concentrar hasta 500 veces los rayos del sol.
{
I = I ph I o e ko (V + Rs I ) 1 } V
RSH
Potencia mxima o pico, Pmx: punto de operacin de la clula (IMPP, VMPP; MPP
maximum power point) para el cual la potencia generada es mxima, esto es:
Pmax=IMPPVMPP
I MPP VMPP
FF = Pmx = FF I SC VOC
I SC VOC
V
I = I sc 1 C1 e 2 oc 1
C V
Donde
1
I C2 Voc
VMPP
V I MPP
C1 = 1 MPP e C2 = MPP 1 ln 1
I SC Voc I SC
Irradiancia.
Temperatura.
Los paneles o mdulos se fabrican en una amplia gama de tamaos, con el fin de servir a
diferentes aplicaciones.
Entre los paneles fotovoltaicos ha de existir una cierta distancia, para evitar el efecto de
sombra.
Para estructuras fijas la distancia entre paneles se calcula mediante la siguiente expresin:
h
d=
tg (61 latitud )
Se ha de hacer notar que los mdulos fijos pueden disponerse a distancias menores que los
colocados en estructuras de seguimiento solar, puesto que estos ltimos producen mayores
superficies de sombra.
Instalaciones pequeas de 3 kWp como planta tipo (con rango hasta 5 kW).
Son aplicaciones rurales aisladas por ser una solucin limpia y muchas veces econmica,
o aplicaciones conectadas a red sobre tejados, azoteas de casas, hechas por particulares
en zonas de su propiedad o influencia; la motivacin es generalmente medioambiental.
Con la generacin de 3 kWp se cubrira el consumo propio de una casa tipo medio en el
que vivan 2-3 personas, excluyendo el consumo de calefaccin y aire acondicionado.
Instalaciones grandes de 300 kWp como planta tipo (entre 100 kW y 1MWp).
Adems de los paneles fotovoltaicos que alimentan los satlites artificiales, ya a nivel
terrestre el sector de las telecomunicaciones es tambin uno de los mayores usuarios de este tipo
de energa, ya que sus equipos usualmente se encuentran en lo alto de las montaas, en lugares de
difcil acceso o en zonas deshabitadas.
Un tipo de aplicaciones distintas a las anteriores son las que no requieren la utilizacin de
acumuladores y que funcionan siempre que haya sol, consumiendo la energa que generan los
paneles. Se trata de sistemas relativamente simples, siendo los sistemas de bombeo de agua el caso
ms utilizado.
El generador fotovoltaico se estima que tiene una vida til superior a 30 aos, siendo la
parte ms fiable de la instalacin. La experiencia indica que los paneles nunca dejan de producir
electricidad, aunque su rendimiento pueda disminuir ligeramente con el tiempo.
Los paneles que forma el generador apenas requieren mantenimiento, basta limpiarlos con
algn producto no abrasivo cuando se detecte suciedad solidificada.
Por otro lado, las bateras con un correcto mantenimiento tienen una vida
aproximada de diez aos. Se debe controlar que el nivel de agua del electrolito est dentro de unos
lmites aceptables. Para reponerlo se utiliza agua desmineralizada o destilada. Se debe revisar su
nivel mensualmente en cada uno de los elementos y mantener los bornes de conexin libres de
sulfato. Actualmente existen bateras sin mantenimiento o de electrolito gelificado que no
necesitan reposicin de agua.
Entre las principales aplicaciones de los sistemas conectados a la red elctrica son los
tejados de vivienda. Son sistemas modulares de fcil instalacin donde se aprovecha la superficie
de tejado existente para sobreponer los mdulos fotovoltaicos. El peso de los paneles sobre el
tejado no supone una sobrecarga para la mayora de los tejados existentes. Una instalacin de unos
3 kWp que ocupa cerca de 30 metros de tejado inyectara a la red la energa necesaria para una
pequea vivienda. El mantenimiento se reduce a la limpieza de los paneles, cuando se detecte
suciedad solidificada.
Plomo-cido.
Nquel-cadmio.
Nquel-hierro.
Plata-cadmio.
Zinc-xido de plata.
Figura. Bateras de gel. Tambin son estacionarias. El electrodo est gelificado, de modo que no
precisan mantenimiento y pueden funcionar colocadas en cualquier posicin.
Para una planta de 1 MW que utilice estructuras fijas, la superficie de ocupacin ser
aproximadamente de 1,4 hectreas, es decir x = y =120 metros.
La superficie ocupada por una planta fotovoltaica de 1 MW con seguidores de un solo eje
ser aproximadamente de 4,6 hectreas (x = y = 215 metros), considerando condiciones de
radiacin en torno a 2.120 kWh/m2ao. Igualmente, con seguidores de dos ejes la superficie
ocupada sera de 5 hectreas (x = y = 225 metros).
La regulacin del seguimiento del sol por dos ejes se realiza bsicamente mediante los
siguientes modos:
Los mdulos tienen una vida esperada de ms de 40 aos. Realmente no se tienen datos
para saber con exactitud la vida real de un generador conectado a red porque no se tiene suficiente
perspectiva.
La vida til de los restantes elementos que componen la planta FV, inversores y medidores,
as como los elementos auxiliares, cableado, canalizaciones, cajas de conexin, etc., es la vida til
tpica de todo tipo equipo electrnico y material elctrico, la cual es compatible con la larga vida
til del generador FV, con el adecuado mantenimiento.
2.8 El inversor.
Se denominan inversores a los convertidores de corriente continua (Direct Current, DC) a
corriente alterna (Alternating Current, AC). La funcin de un inversor es cambiar una tensin de
entrada de DC a una tensin simtrica de salida de AC, con la magnitud y frecuencia deseadas.
Las formas de onda de la tensin de salida de los inversores prcticos no son senoidales y
contienen ciertas armnicas. En aplicaciones de potencia baja e intermedia se pueden aceptar
tensiones de onda cuadrada o de onda casi cuadrada, y para aplicaciones con alta potencia se
requieren formas de onda senoidal con poca distorsin. Con la disponibilidad de los dispositivos
semiconductores de potencia de alta velocidad, se pueden minimizar los contenidos de armnicos
en la tensin de salida, o al menos reducirlos de forma importante, mediante tcnicas de
conmutacin.
La entrada del inversor puede ser una batera, una celda de combustible, celda solar u otra
fuente de DC.
Un inversor se llama inversor alimentado por tensin (Voltage Source Inverter, VSI) si la
tensin de entrada permanece constante; inversor alimentado por corriente (Current Source Inverter,
CSI) si la corriente de entrada permanece constante, y convertidor enlazado con DC variable si la
tensin de entrada es controlable.
Para obtener tensiones negativas, es necesario contar con dos clulas de conmutacin. Tal
es el caso del siguiente puente inversor monofsico:
Estos inversores usan en general seales de control por modulacin por ancho de pulso
(PWM) para producir una tensin de salida de AC.
Elevador y puente en H
Inversores de dos niveles requieren conmutacin en alta frecuencia, junto con diversas
estrategias de modulacin por ancho de pulso (PWM), para obtener una forma de onda de tensin
o corriente de alta calidad. Sin embargo, en aplicaciones con alta potencia y alta tensin, tienen
algunas limitaciones para operar en alta frecuencia, principalmente a causa de las prdidas por
conmutacin y limitaciones de las especificaciones nominales de los dispositivos.
En trminos generales, puede ser ms fcil producir un inversor de alta potencia y alta
tensin con la estructura multinivel, al aumentar la cantidad de niveles de tensin sin necesidad de
mayores especificaciones nominales de los dispositivos semiconductores. Cabe destacar, adems,
que a medida que aumenta la cantidad de niveles de tensin se reduce en forma importante el
contenido de armnicos en la forma de onda de la tensin de salida, permitiendo el empleo de
filtros de menor tamao. Se debe sealar, sin embargo, que un elevado nmero de dispositivos
semiconductores aumenta el coste y complejidad de control del convertidor. Supone, adems, un
reto para este ltimo el equilibrado de tensiones debidas al elevado nmero de condensadores.
Conviene decir, por otra parte, que los sistemas fotovoltaicos (FV) conectados a red,
en particular de baja potencia e integrados en las arquitecturas o tejados de los edificios, son
reconocidos cada vez en mayor medida en todo el mundo por su contribucin a la generacin
limpia de energa elctrica. Este tipo de sistemas FV generalmente son propiedad de un particular
o particulares, llegando a alcanzar una potencia de hasta 10 kW. Lgicamente, el principal objetivo
de estos operadores privados es maximizar la produccin de energa. Cuestiones como la vida til
del sistema (20 aos o ms), elevada eficiencia y buenas condiciones medioambientales
(disponibilidad de radiacin solar) son de importancia para estos operadores. Otros requerimientos
importantes de estos sistemas FV son el cumplimiento de las normas de calidad de la energa
elctrica, compatibilidad electromagntica, limitacin de ruidos, o requerimientos de seguridad y
proteccin. A pesar de ello, su implantacin o aceptacin depender en buena medida de su
fiabilidad y bajo coste.
La tendencia actual en este campo se basa en unidades dispuestas en cadena (ver apartado
siguiente) de potencias de 1 kW sin transformador (en particular en inversores con tensiones de
entrada elevadas), debido a su elevada eficiencia. Evitar el transformador trae consigo beneficios
como la reduccin del coste, tamao, peso y complejidad del inversor. Sin embargo, prescindir del
transformador y por consiguiente del aislamiento galvnico que proporciona debe ser considerado
cuidadosamente.
Figura. (a) Mxima (~nF) y (b) mnima (~pF) capacidad a tierra de un mdulo FV.
Figura. Sistemas FV conectados a red con (a) inversor de medio puente con fijador
de nivel por diodos de tres niveles y (b) inversor de medio puente con fijador
de nivel por diodos de cinco niveles
La conexin a tierra del sistema FV en su punto medio supone una clara ventaja, ya que
elimina las corrientes capacitivas a tierra y su influencia negativa en la compatibilidad
electromagntica del circuito.
Figura. Sistemas FV conectados a red con (a) inversor de puente completo con fijador de nivel por
interruptores y (b) inversor de puente completo con fijador de nivel con diodos
Un sistema FV sin transformador con similares caractersticas puede ser realizado con un
inversor de puente completo con fijador de nivel con diodos (ver figura anterior). En ambas
topologas los subconjuntos FV estn simtricamente cargados.
2.8.1.1.3 Cascada.
Para la conversin de energa DC a AC se muestra a continuacin un sistema FV
conectado a red sin transformador con un inversor en cascada (cascaded inverter). La topologa se
compone de dos puentes completos con sus salidas AC conectadas en serie. Cada puente puede
crear tres niveles diferentes de tensin AC, permitiendo generar una salida de tensin global de
cinco niveles. La ventaja de esta topologa es su carcter modular. En la bibliografa especializada
este concepto se sugiere en sistemas con una tensin de 40 V en cada bus DC, adems de en
aplicaciones de alta potencia.
2.8.1.1.4 Paso.
Un convertidor de paso (step converter) conmuta subconjuntos FV de diferentes tensiones
para obtener una tensin AC de salida. J. Mller, en su artculo Assesment and Optimisation of Step-
Inverters for Photovoltaic Systems, expone una topologa compuesta por cinco mdulos con tensiones
nominales de 11 V, 22 V, 44 V, 88 V y 176 V para un sistema FV conectado a red como se
muestra en la siguiente figura:
La energa generada por cada subconjunto FV crece con la tensin del bus DC.
V0 n
HFn = n >1
V01
La distorsin armnica total, que es una medida de la coincidencia de formas entre una
onda y su componente fundamental, se define como
1 2
THD = V02n
V01 n =02, 3, ...
1 V0 n 2
DF = 2
V01 n = 2, 3, ... n
V0 n
DF =
V01 n 2
Cuando uno analiza los costes de los componentes de la energa fotovoltaica (ver
figura a continuacin), se muestra claramente que pagar un precio mayor por el
convertidor probablemente pueda ser justificado por una potencial superior reduccin
de los costes de la instalacin.
Entre todos los elementos activos y pasivos, los condensadores electrolticos tienen
una menor vida til, y con la desventaja de que incrementan el valor de la resistencia
serie equivalente y consecuentemente tambin las prdidas. Esta topologa no
requerir un alto valor de este condensador a fin de permitir el uso de la tecnologa de
pelcula, especialmente del tipo polipropileno, la cual tiene una larga vida til, buena
estabilidad trmica y elctrica y finalmente una baja resistencia serie.
Los grupos de mdulos no alcanzan una potencia superior a los 2 kW. Su potencia
total puede alcanzar varios MW, pero su implantacin no es an efectiva.
Es apropiado para pequeas instalaciones fotovoltaicas, como una nica vivienda familiar.
a0
f (t ) = + [a m cos( mwt ) + bm sen( mwt ) ]
2 m =1
donde
1
am =
f (t ) cos(mwt ) dwt
m = 0,1, ...
1
bm =
f (t ) sen(mwt ) dwt
m = 1, 2, ...
Para una mejor comprensin, se muestra la siguiente figura, donde una seal peridica es
suma de dos funciones seno:
Igualmente, esta teora puede desarrollarse para una forma de onda f(x,y) dependiente de
dos variables, esto es:
A00
f ( x, y ) = + [A0 n cos(ny ) + B0 n sen(ny )] + [Am 0 cos(mx) + Bm 0 sen( mx)]
2 n =1 m =1
+ [A
m =1 n =
mn cos(mx + ny ) + Bmn sen(mx + ny )]
( n0)
donde
1
Amn =
2 2 f ( x, y) cos(mx + ny) dxdy
1
Bmn =
2 2 f ( x, y)sen(mx + ny) dxdy
A00
f (t ) = + [ A0 n cos(n[ w0 t + 0 ]) + B0 n sen( n[ w0 t + 0 ])]+
2 n=1
DC Offset Armni cos de la componente fundamental de la onda senoidal
+ [Am 0 cos(m[ wc t + c ]) + Bm 0 sen( m[ wc t + c ])]+
m =1
El primer sumatorio, ... donde m=0, define la forma de onda de salida fundamental y
n=1
los armnicos asociados a ella (si existen). Estos ltimos deben ser minimizados o preferiblemente
eliminados con el proceso de modulacin.
El segundo sumatorio, ... donde n=0, corresponde a los armnicos de la onda
m =1
portadora, los cuales son relativamente componentes de alta frecuencia, ya que el trmino de ms
baja frecuencia es la propia frecuencia de la seal portadora.
El doble sumatorio final, ... donde m, n 0, se refiere a los armnicos
m =1 (nn=
0)
Para ilustrar esta idea, se considera la estructura general de un inversor alimentado por
tensin monofsico (single-phase VSI). El apagado y encendido de cada dispositivo semiconductor
se realiza mediante la comparacin de una seal moduladora triangular (portadora) y una seal
senoidal de referencia.
Ar
M=
Ac
Igualmente, para los interruptores T2 y T4 se obtiene la tensin Vbz desfasada 180 de Vaz:
Vdc
1
+ 4 mJ m M sen [ m + n] cos( m wc t + n w0 t )
m =1 n =
n
2 2
Vdc
1
+ 4 mJ m M sen [ m + n] cos( m wc t + n [ w0 t ])
m =1 n =
n
2 2
+ 8
Vdc
2m J
1
2 n 1 (mM ) sen ([m + n + 1] ) cos(2m wc t + [2n 1]w0 t )
m =1 n =
Esta expresin muestra que los armnicos impares y de banda laterales presentes en la
tensin de fase debidos a la seal triangular, son completamente cancelados en le tren de pulsos de
la tensin de lnea de salida. En ella quedan sin embargo los trminos impares (2n-1) debidos a la
seal de referencia y los trminos pares (2m) de la seal triangular.
La siguiente figura ilustra el espectro frecuencial de las tensiones de fase y de lnea de salida
expuestas anteriormente:
Figura. Espectro armnico terico de un inversor monofsico modulado por PWM: (a) Fase a y
(b) Tensin de lnea de salida, M=0.9, fc/fo=21.
Las ms utilizadas en aplicaciones industriales son la tercera y, sobre todo, la cuarta, debido
a su mejor factor de potencia y baja distorsin armnica, por lo que nos centraremos en ellas.
Empecemos por la modulacin por ancho de pulso senoidal, ya presentada anteriormente:
En lugar de mantener igual el ancho de todos los pulsos, como en el caso de la modulacin
de varios pulsos, se hace variar el ancho de cada pulso en proporcin con la amplitud de una onda
senoidal. El DF y la LOH se reducen en forma apreciable. Dichas seales de control (g1, g2, etc.)
de los dispositivos semiconductores, como se ve en la figura anterior, se generan comparando una
seal senoidal de referencia con una onda portadora triangular de frecuencia fc. La tensin de salida
es vo = Vs (g1-g4). La cantidad de pulsos por medio ciclo, p, depende de la frecuencia de la portadora.
fc m f
p= =
2 f0 2
Esta clase de modulacin elimina todas las armnicas menores o iguales a 2p-1. Para p=5,
la LOH es la novena.
Se pueden generar las mismas seales de disparo con una onda portadora triangular
unidireccional como se ve en el caso (d) de la figura anterior. Es ms fcil de implementar este
mtodo, y es preferible.
La representacin (c) de la figura anterior indica que los anchos de los pulsos ms cercanos
al pico de la onda senoidal no cambian mucho al variar el ndice de modulacin. De esta manera, la
tcnica anterior se puede modificar para que se aplique la onda portadora durante los primeros y
ltimos intervalos de 60 por medio ciclo. Esta modulacin por ancho de pulso senoidal
modificada (modified sinusoidal pulse-width modulation, MSPWM) se ve en la siguiente figura.
Aumenta la componente fundamental, y mejora sus caractersticas armnicas. Reduce la cantidad
de conmutacin de los dispositivos de potencia, y tambin reduce las prdidas por conmutacin.
Figura. Perfil de armnicos con modulacin por ancho de pulso senoidal modificado.
fc
= 6q + 3
fo
1. Modulacin trapezoidal.
Las seales de compuerta se generan comparando una onda portadora triangular con una
onda moduladora trapezoidal como se ve en la siguiente figura.
Esta clase de modulacin aumenta la tensin pico fundamental de salida hasta a 1.05 Vs,
pero la salida contiene armnicos de orden menor (LOH).
5. Modulacin delta.
En la modulacin delta se deja oscilar una onda triangular dentro de una ventana definida
V, arriba y abajo de la onda senoide de referencia vr. La funcin de conmutacin del inversor,
que es idntica a la tensin de salida v0, se genera desde los vrtices de la onda triangular vc, como
se ve en la figura. Tambin se llama modulacin por histresis o bang-bang. Si cambia la frecuencia de la
onda moduladora, manteniendo constante la pendiente de la onda triangular, cambia la cantidad de
pulsos y de anchos de pulso de la onda modulada.
La tensin fundamental de salida puede ser hasta 1Vs, y depende de la amplitud pico Ar y
de la frecuencia fr de la tensin de referencia. La modulacin delta puede controlar la relacin de
tensin a frecuencia, lo que es una caracterstica deseable, en especial para controlar motores.
Figura. Rizado de corriente de salida como una funcin del ndice de modulacin.
Entre todos ellos, los mtodos de perturbacin y observacin son los ms populares,
debido a su fcil implementacin.
A continuacin se muestra el diagrama de control bsico del convertidor elevador con lazo
cerrado de tensin de salida, lazo cerrado de corriente y control MPPT.
I p = I ph I o exp ( KV c / N s )
Pp = Vc I p = Vc (I ph I o exp ( KVc / N s ) )
El cambio de potencia de salida del sistema fotovoltaico bajo una perturbacin de tensin
puede ser dado por:
Wc C1 Vc V
Pc = =
T T
donde Vc es la tensin en C1 y T es el periodo MPPT. La precisin del MPPT puede ser mejorada
disminuyendo el paso de tensin o aumentado el periodo MPPT. Estas medidas , sin embargo,
disminuyen la velocidad de respuesta del MPPT.
En la siguiente figura se muestra una curva P-I de un sistema fotovoltaico compuesto por
8 mdulos de 130 Wp en serie bajo condiciones de radiacin solar 750 W/m2, 1.5 masa de aire y
25 C de temperatura. Adems, se representa PC, Pp y la curva de cambio de potencia en el bus
DC Pi(Pp- PC ), con paso de tensin 1.56 V y periodo MPPT de 160 ms.
PA
MPPT = x 100 %
Pm
La precisin del MPPT puede ser mejorada con alimentacin hacia delante de la potencia
del condensador del filtro. El algoritmo MPPT puede ser dado como sigue,
Vc = Vc ( k 1) + V (k 1)sign ( Pi + Pc )
Las teoras que tratan con potencias instantneas pueden ser clasificadas
fundamentalmente en dos grupos. El primero de ellos se basa en la transformacin de las fases abc
a tres ejes ortogonales, y el segundo opera directamente sobre las fases abc. La primera, que ser
llamada teora pq, se basa en la transformacin abc a o. En el mbito que nos ocupa, la teora pq
evidencia un tratamiento ms eficiente con corrientes no senoidales que la anterior,
particularmente cuando estn presentes componentes de secuencia cero. Adems, desde el punto
de vista de la compensacin, la teora pq permite la compensacin separada de las potencias activa
p, y reactiva q, adems de potencia activa de secuencia cero, po, siendo esta la razn por la que
esta teora puede compensar incluso bajo tensiones distorsionadas y/o desequilibradas.
1 1 1
2 2 2 v
v o a
2 1 1
v = 3 1
2
v b
2
v v
3 3 c
0
2 2
1
1 0
2 v o
v a
v = 2 1
1 3
v
b 3 2 2 2
vc v
1 1 3
2 2 2
Igualmente, las intensidades instantneas de lnea pueden ser transformadas a los ejes o,
esto es:
1 1 1
2 2 2
io ia
2 1 1
i = 3 1
2
2 ib
i ic
3 3
0
2 2
y su transformada inversa:
1
1 0
2 io
ia
i = 2 1
1 3
i
b 3 2 2 2
ic i
1 1 3
2 2 2
1 1
1 v a
v 2 2 2
v = vb
3 3 3
0 v
2 2 c
1 0
v a
v = 2 1 3 v
b
3 2 2 v
vc
1
3
2 2
Tres potencias instantneas (la potencia instantnea de secuencia cero po, la potencia activa
instantnea p, y la potencia reactiva instantnea q) se desprenden de los ejes o a partir de las
corrientes de lnea y tensiones de fase:
p o v o 0 0
p = 0 v v
q 0 v v
i 1 v v p
i = 2
v + v
2
v v q
Figura. Esquema de variables elctricas que forman parte del control pq.
Fotovoltaica.
Elica.
Clulas de combustible.
Microturbinas.
Cogeneracin.
Estas tecnologas juegan ya un papel importante en los picos de generacin, durante las
horas de mayor demanda, cuando los costes elctricos son altos. Junto a otros beneficios, y a pesar
de ello, hay lmites tcnicos asociados con la proliferacin de unidades de generacin distribuida,
como puede ser el funcionamiento en isla.
Qu es el funcionamiento en isla?
Se trata de la generacin de energa elctrica por parte de estas unidades cuando la red
elctrica se encuentra, aguas arriba, desconectada, debido a una operacin de falta, un evento
planeado de mantenimiento, o al fallo de algn equipo o dispositivo de la aparamenta del sistema
elctrico.
Como para todos los equipos de proteccin, el rendimiento de la proteccin en isla puede
ser cuantificada mediante los siguientes parmetros:
Tiempo de operacin.
UL 1741-2001
IEEE 1547.1-2005
5. Mtodos hbridos.
Si la potencia elctrica suministrada por el inversor, y por tanto por la planta fotovoltaica,
es similar a la potencia total consumida por las cargas aguas abajo, la variacin de tensin y
frecuencia en el PCC son menos evidentes cuando opera en isla. En tales casos, estos mtodos
pierden su efectividad, al incurrir en una zona de no deteccin (Non-Dectection Zone, NDZ).
IQ = 0
PLOAD = P + P
QLOAD = Q + Q
P
I=
V
V2 V2
Z LOAD = =
P + P + j (Q + Q ) P + P + j Q
'
Z LOAD Z LOAD
V
V ' = I Z LOAD
'
I Z LOAD =
P Q
2 2
1 + +
P P
P Q V
2 2 2
1 + + =
P P V '
La zona de no deteccin del mtodo OVP/UVP para DG controlada por corriente es:
2 2
P Q V
2 2
V
1 + +
Vmax P P Vmin
Sabemos de la norma IEEE Std. 1547 que la tolerancia respecto a la amplitud normal de
tensin es 88 % ~ 110 %:
Vmin = 0.88 V
Vmax = 1.10 V
P Q
2 2
0.91 1 +
2
+ 1.14
2
P P
Figura. NDZ del mtodo OVP/UVP para sistemas DG controlados por corriente.
V2
PLOAD = = P + P
R
1
QLOAD = Q = V 2 2fC
2fL
1
f'=
2 LC
Q C 1
=R 2f LC
PLOAD L 2f LC
Y sustituyendo:
1
Q Q P f' f
= 1 + =Qf
PLOAD P P f f '
C
Qf = R
L
1
f f Q P f f
Q f min 1 + Q f max
f f min P P f f max
donde fmax y fmin son los valores predeterminados del umbral de frecuencias antes mencionado. De
la norma IEEE Std. 1547, conocemos que la tolerancia de frecuencia de la tensin normal es 98.83
%~100.83 %, considerando una red de 50 Hz, tenemos
f min = 49.4 Hz
f max = 50.4 Hz
1
Q P
0.024 * Q f 1 + 0.016 * Q f
P P
Figura. NDZ del mtodo OFP/UFP para sistemas DG controlados por corriente.
Tiempo de disparo
P 50 % 10 ciclos
Fdp carga aislada <0.95 10 ciclos
P 50 % + Fdp carga aislada <0.95 10 ciclos
P < 50 % + Fdp carga aislada > 0.95 2 seg
Los mtodos activos tratan de superar las dificultades de deteccin en isla de los mtodos
pasivos mediante la inyeccin de una perturbacin en el sistema elctrico, causando una condicin
anormal de las magnitudes elctricas bajo condiciones de funcionamiento en isla. La
representacin en el espacio P-Q de la NDZ no es adecuada para los mtodos activos, ya que
para un determinado intercambio de potencia reactiva, Q, diferentes combinaciones de L y C son
posibles. Algunas de estas combinaciones resultan en funcionamiento en isla y otras no.
Todos los mtodos de deteccin en isla activos deben cumplir las normas internacionales a
tal efecto. As, la distorsin armnica total (THD) de la corriente suministrada por los inversores
distribuidos debe ser menor de un 5 %. Ello conduce a un incremento del tiempo de deteccin y
afecta a la fiabilidad de la operacin.
Segn el IEEE Std. 929-2000, el factor de calidad Qf est definido como 2 veces la
relacin entre la mxima energa almacenada y la energa disipada por ciclo a una frecuencia dada:
1
2 C R 2 I 2
Qf = 2 = w R C = R = R C
R I / w0
2 0
w0 L L
donde
1
w0 = 2f 0 =
LC
frecuencia resonante en rad/s.
La magnitud y fase de la impedancia de una carga paralela RLC a una frecuencia arbitraria f,
en funcin del factor de calidad Qf y la frecuencia resonante f0 de la carga, son dadas por:
1 R
Z LOAD = =
2 2
1 1 f f
+ wL 1 + Q f 0
R 2 wL f f0
1 w 2 LC f f
LOAD = tan R
1
= tan 1 Q f 0
wL f f 0
Cabe destacar, en primer lugar, que la magnitud de la impedancia de carga cambia muy
poco para cargas con Qf < 2.5. Por otra parte, debido a que el margen de frecuencias de operacin
es relativamente pequeo, 59.4 Hz a 60.5 Hz, impuesto, entre otras, por la norma IEEE Std. 929-
2000, se asume que la magnitud de tensin no vara con la frecuencia del sistema.
Figura. Variacin del ngulo de fase de una carga RLC con la frecuencia para diferentes valores de
frecuencia de resonancia (fo) y factor de calidad (Qf)
Las curvas muestran que a medida que aumenta el factor de calidad Qf, hay una mayor
variacin del ngulo de fase de la carga LOAD en torno a la frecuencia del sistema. Adems, la
interseccin de las curvas, ngulo de fase cero, tiene lugar en la frecuencia de resonancia fo. Por
otra parte, es interesante hacer notar que una carga cuya frecuencia de resonancia sea menor que la
frecuencia de la red (fo<fg) se muestra como una red capacitiva, ya que la corriente adelanta a la
tensin (LOAD>0). Igualmente, una carga con fo>fg es una red inductiva a la frecuencia de la red.
Este mtodo es resistente a falsos disparos ocasionados por saltos de fase aleatorios en la
tensin de la red. Por otro lado, se ha de tener en cuenta la posibilidad de interferencias entre
mltiples unidades y la habilidad de detectar fiablemente cambios de impedancia con valores
prcticos de la carga insertada.
Figura. Forma de onda de salida del inversor usando el mtodo SMS: (a) fase de adelanto de la
intensidad, (b) fase de atraso de la intensidad.
Este mtodo es relativamente fcil de implementar, ya que requiere solo una ligera
modificacin de un componente que ya es requerido, el PLL.
Figura. Marcos.
Figura. Esquema SRF-PLL (Synchronous Reference Frame PLL.). Este esquema es simple
y usado en casi todas las tcnicas PLL para sistemas trifsicos.
f fg
SMS = M sen
2 fm f g
donde fm es la frecuencia a la cual sucede el ngulo de fase mximo, M.
Despus de la desconexin de la red, la frecuencia del sistema aislado se desviar del valor
nominal fg si:
d LOAD d SMS
<
df f = fg
df f = fg
LOAD 12 Q f
fm f g 2
donde Qf es el valor mximo para el cual la condicin de isla ha de ser detectada, y fm-fg es tomada
normalmente como 3 Hz.
La NDZ del SMS IDM en el espacio Qf-fo se obtiene usando el criterio de fase (phase
criteria), esto es:
f o f is f fg
tan 1 Q f = M sen is
f f 2 f f
is o m g
Por tanto,
f is tan ( SMS ( f is ) )
f o2 f o f is2 = 0
Qf
Puede observarse que la NDZ es nula para Qf<2.5 cuando M=10 y fm-fg=3 Hz como
especificaciones de diseo. Adems, mientras M disminuye, el factor de calidad para el cual el
funcionamiento en isla no es detectado decrece. El rendimiento de este mtodo es similar al
mtodo pasivo UFP/OFP (M=0) para cargas con Qf elevado.
Figura. IPV para AFD (a) hacia arriba (drift up) y (b) hacia abajo (drift down).
La relacin entre el tiempo cero, tz, y la mitad del periodo de la onda de tensin, TVutil/2, se
denomina fraccin de corte (chopping fraction, cf):
2tz
cf =
TVutil
En rgimen permanente la magnitud Vis y frecuencia fis en un sistema aislado con una carga
local RLC pueden ser calculadas desde las siguientes ecuaciones:
Vis2
PLOAD ,is =
R
Vis2 f 2
QLOAD ,is = o 1 = Qinv
X c,g f is
AFD puede ser implementado forzando la frecuencia de la corriente a ser f (<1.5 Hz) por
encima de la frecuencia de la tensin en el ciclo previo (fik=fvk-1+ f), y manteniendo la corriente
del inversor igual a cero desde el final de su semiciclo negativo hasta el paso por cero de la tensin
positiva. Cuando la red elctrica est conectada, sta mantiene la frecuencia de la tensin. Cuando
no es as, la frecuencia de la tensin en el PCC tiende a aumentar (drift upward), alcanzando valores
por encima de fo. En ambos casos, la componente fundamental de la corriente adelanta a la
tensin un pequeo ngulo AFD.
ik = 2 I sen (2 ( f vk 1 + f ) t )
t z / 2 AFC
=
TV 2
donde
1 1 1 1 f
tz = = =
fV f i f f + f f ( f + f )
El ngulo de fase del inversor para el mtodo de deteccin en isla AFD es:
f
AFC = f t z =
f + f
f o f is f
tan 1 Q f =
f is f o f is + f
Por tanto,
f is tan ( AFD ( f is ) )
f o2 f o f is2 = 0
Qf
La siguiente figura muestra la NDZ de un AFD IDM (method detection islanding) para
diferentes valores de f:
Puede apreciarse, en la figura anterior, que la NDZ se desplaza hacia menores valores de
frecuencia de resonancia cuando f crece. Para Qf=2.5 y f=0.5 Hz, el mtodo falla en la
deteccin de funcionamiento en isla para cargas con 58.99 Hz fo 60.19 Hz. Para cargas con
factores de calidad por encima de 2.5 este mtodo activo no presenta ninguna mejora con respecto
a la proteccin pasiva UFP/OFP, ya que la capacidad de desplazamiento de la frecuencia es
significativamente reducida. Adems, para cargas capacitivas (fo<60 Hz) puede incurrir en fallos de
deteccin en isla, como se observa en la figura. Sin embargo, ha de tenerse en cuenta que este
mtodo es bastante efectivo para cargas inductivas.
Conviene decir, por otra parte, que bajo estas condiciones de anlisis (P=0 y Q=0,
Vis=Vg) el modo de control del inversor, por corriente constante o por potencia constante, no
tiene impacto sobre la NDZ.
cf = cf 0 + k ( f f g )
Donde cfo es la fraccin de corte cuando no hay un error de frecuencia y k es una ganancia
de aceleracin. Si la fraccin de corte aumenta, la frecuencia del inversor fotovoltaico se
incrementa, provocando la actuacin del rel de salida del inversor.
Para cfo=0.05, la distorsin armnica total (Total Harmonic Distortion, THD) de la corriente
del inversor es igual a un 5 % de la frecuencia nominal.
Este mtodo posee una de las ms pequeas zonas de no deteccin de todos los mtodos
de deteccin en isla activos. Adems, no pierde efectividad si es empleado por mltiples inversores
en isla, y puede ser extremadamente efectivo si es implementado en combinacin con el mtodo
Cambio Sandia de la Tensin (Sandia Voltage Shift, SVS).
Sin embargo, con este mtodo la calidad de la energa elctrica de salida del inversor DG se
reduce ligeramente. Tambin, la forma de onda de referencia puede ser estimulada por ruido o
armnicos.
wt z cf ( f )
SFS ( f ) = =
2 2
La NDZ de SFS IDM en el espacio Qf-fo se desprende de la expresin del criterio de fase
(LOAD =SFS):
f o f is
tan 1 Q f = [cf o + k ( f is f o )]
f is f o 2
Por tanto,
f is tan ( SFS ( f is ) )
f o2 f o f is2 = 0
Qf
A continuacin se muestra grficamente el resultado del clculo de la NDZ del SFS IDM.
Se observa que la NDZ es nula para cargas con Qf<4.8 cuando cfo=0.05 y k=0.1. Por otra
parte, a medida que k disminuye, el factor de calidad Qf para el cual no se detecta el
funcionamiento en isla tambin decrece. Para k=0 este mtodo llega a ser el AFD IDM, donde
habr cargas, mayormente capacitivas, para el cual no ser posible la deteccin del funcionamiento
en isla. Finalmente, cabe decir que el rendimiento de este mtodo es similar al mtodo pasivo
UFP/OFP para cargas con elevados factores de calidad.
Conviene decir, adems, que la NDZ en el espacio Qf-fo derivada para un inversor puede
ser directamente extendida para el caso de mltiples inversores, si emplean el mismo mtodo
activo y se encuentran razonablemente cerca. Sin embargo, uno debera usar un factor de calidad
modificado para la carga local (Qf_A&P), que es dado por:
Qf
Q f _ A& P =
1 PIDM
Este mtodo complementa el mtodo SFS anteriormente expuesto y puede ser usado en
caso de mltiples inversores.
Este mtodo presenta cierta dificultad de implementacin debida a la eleccin del valor
umbral. Un ajuste altamente sensible puede provocar molestos disparos de corte de suministro del
inversor DG, sin embargo la fiabilidad de deteccin en isla depende de la correcta eleccin de este
valor umbral.
Como punto dbil puede argirse una pequea reduccin de la calidad de energa elctrica
suministrada. La penetracin en el sistema elctrico de inversores DG que implementan este
mtodo debe mantenerse, por ello, baja.
Para detectar una isla accidental, se prueba la continuidad de la lnea de potencia enviando
una seal por el transmisor (T). Cuando la seal PLCC se pierde el receptor (R) puede ordenar el
cese da la operacin del inversor o inversores, o puede abrir su propio interruptor para aislar la
planta fotovoltaica de la carga del PCC.
Puede decirse que este mtodo posee mltiples ventajas: no tiene una NDZ, la calidad de
la energa elctrica suministrada por el inversor DG no est degradada, el nmero de inversores en
el sistema no afecta al rendimiento, y podra ser efectivo en cualquier nivel de penetracin de
inversores fotovoltaicos en la red. Sin embargo, se ha de considerar el coste del receptor y del
transmisor. Adems, existe la posibilidad de que una carga replique la seal PLCC siendo sta la
que detecte el dispositivo de deteccin, no permitiendo detectar la operacin anormal en isla.
Las compaas elctricas poseen aos de experiencia con este esquema para varias
aplicaciones de proteccin. Esta nueva aplicacin puede ser por tanto fcilmente aceptada. Este
mtodo, por otra parte, podra permitir la coordinacin de las unidades DG y los recursos de la
red elctrica. El mismo sistema podra ser usado tambin para la reconexin de los generadores
distribuidos una vez clareada la falta.
Su principal desventaja son los costes y potencial complejidad, ya que se hace necesario el
empleo de transmisores de seal para todos los posibles puntos de desconexin del sistema. Si hay
muchos interruptores automticos y la topologa de la red vara, un esquema de transferencia de
disparo puede llegar a ser bastante complicado.
Este mtodo presenta las siguientes ventajas: los bancos de condensadores requeridos
estn fcilmente disponibles en el mercado, las compaas elctricas tienen una gran experiencia en
ellos, pueden ser usados en caso de mltiples inversores, y no presenta zona de no deteccin. Por
otra parte, no est claro qu parte podra ser responsable del gasto del banco de condensadores.
Un nuevo mtodo propone la operacin del inversor como una inductancia virtual, ya que
la frecuencia de la tensin de salida del inversor es ligeramente superior a la frecuencia de la red.
Cuando se produce un corte de suministro de energa elctrica en la red, la amplitud o la
frecuencia de la tensin en la carga cambia. El funcionamiento en isla de la unidad DG se asocia a
este cambio, y la condicin de operacin en isla puede ser detectada.
La seal S2 es empleada para controlar el inversor conectado a red como una inductancia
virtual. La frecuencia de operacin de la inductancia virtual es ligeramente superior a la frecuencia
fundamental de la red.
El filtro pasa-banda (II) posee una estrecha banda de paso, y est sintonizado a una
frecuencia ligeramente superior a la frecuencia fundamental. Ya que la corriente de la inductancia
es proporcional a la integracin de la tensin, la salida de este filtro es enviada al integrador. El
filtro paso alto es empleado para bloquear la componente de continua generada por el inversor.
Bajo condiciones normales de operacin, S2 es muy pequea (casi cero). En este caso, el
sistema de generacin elctrica distribuida suministra solo potencia activa. Sin embargo, la
operacin de inductancia virtual proporcionar una retroalimentacin positiva. As, en caso de
funcionamiento en isla, provocar un cambio en la potencia reactiva general del sistema,
incrementndose con el tiempo. Se consigue as un cambio significativo en la amplitud de la
tensin en la carga. Ambos fenmenos, una variacin en la amplitud y la frecuencia de la tensin
en la carga, pueden ser detectados como condicin de funcionamiento en isla.
El concepto de gen escondido (hidden gene) se ha introducido para combinar los mritos de
los mtodos pasivos y activos. Un filtro de media movible es introducido en el controlador de
modulacin de frecuencia como un gen escondido, cuyo efecto no ser perceptible en una
situacin de red elctrica conectada. Sin embargo, la condicin de operacin en isla provocar que
el gen escondido se exprese distintivamente, proporcionado un marcado indicador de isla.
Esto implica que a igualdad de irradiacin solar incidente, un mismo sistema fotovoltaico
producir menos energa en un lugar clido que en un clima fro.
Un problema adicional puede surgir cuando hay sombras sobre el generador FV. En este
caso puede haber escalones en la curva I-V y el inversor operar en un punto que no es el de
mxima potencia.
Captulo 3
FILTRO.
CARGAS.
o CARGA 1.
o CARGA 2.
o CARGA 3.
3.2 Clculos.
3.2.1 Filtro.
P = 500 kW
P = 3 VL I L cos
P 500000
IL = = = 760 A
3 VL cos 3 420 0.95
di VDC i V
L = V primario (t ) L = DC
dt 2 t 2
VDC t
L= {VDC = 640 V } L 0.4 mH
2 i
P = 500 kW
P = 3 VL I L cos
P 500000
IL = = = 760 A
3 VL cos 3 420 0.95
760
8 LNEAS : = 95 A
8
X
X 0.25 R L = = 7.686 10 8 H / m
2f
U 4
100 = 100 = 0.95 %
U 420
20000 N1
m= = 82.48 ; m= N 1 = 82.48 50 = 4124
420 N2
3
S = 3 V L I L
S2 1000 10 3
I 2L = = = 1374.64 A
3 V2 L 3 420
S1 1000 10 3
I 1L = = = 28.86 A
3 V1L 3 20000
Ensayo de vaco :
1700
P0 = V2 N I 0 cos 0 cos 0 = = 0.3922
420
17.87
3
420
V2 V2 3
RFe = = = = 34.6
I Fe I 0 cos 0 17.87 0.3922
l = 0.6 m ; r = 0.07 m ; S = 0.0154 m 2
l 0.6
= r = 500 ( Ncleo de hierro); = = 62008.42 H 1
S 500 4 10 0.0154
7
7 H
0 = 4 10
m
N 2 I 0 50 17.87
N i = 0 = = = 0.0144 Wb
62008.42
Ensayo de cortocircuito :
V2 ccL 6 420
cc = 100 = 6 V2 ccL = = 25.2 V
V2 N 100
25.2
V1ccF = m V2 ccF = 82.48 =1200 V
3
Pcc 10500
Pcc = 3 V1ccL I 1NL cos cc cos cc = = = 0.175
3V1ccL I 1NL 3 1200 28.86
sen cc = 0.9845
X cc
X cc = X 1 + X 2' ; X 1 = X 2' = = 35.45 L1 = 0.1128 H
2
X 2' 35.45
X2 = 2 = = 5.21 m L2 =16.58 uH
m 82.48 2
Clculos:
VL = 20 kV LONG. = 2 km
P = 1 MW cos = 0.9
110 6
P = 3 VL I L cos I L = = 32 A
3 20 10 3 0.9
Re sistencia R
1.03 1.03
Rca = = = 0.245
K Ts S 28 150 km
Re ac tan cia X
0.123 H
Catlago General Cable : X L = 0.123 L = = 3.915 10 4
km 100 km
Tresbolillo
D
L ( H / km) = 210 4 ln 1.28
R
13.9 mm
R= ; D = 30.9 mm; L = 3.477 10 4 H / km; X L = 2 f L = 0.109 / km
2
uF
Cat log o General Cable : C = 0.179
km
LONG = 2 km
U Lnea (V ) = 3 32 2 (0.245 0.8 + 0.123 0.6) = 29.9 V < 100 V (0.5 %)
U 29.9
100 = 100 = 0.15 %
U 20000
VL = 20 kV LONG. = 6 km
P = 3 MW cos = 0.9
3 10 3
P = 3 I L V L cos IL = = 96.22 A
3 2010 3 0.9
LA110
Im ax. Adm = 313.4 A
S = 116.2 mm 2
R = 0.3066 / km
D
L = 0.05 + 0.46 ln 10 4 = 2.725 10 4
H
{D = 1904 mm, r = 5.67 mm }
r km
X L = 0.0856
km
LONG. = 6 km
e 312.52
x 100 = x 100 =1.56 %
V 20.000
66000 N1
m= = 5.715 ; m= N 1 = 5.715 1088 = 6218
20000 N2
3
S = 3 V L I L
S1 3000 10 3
I 1L = = = 26.24 A
3 V1L 3 66000
S2 3000 10 3
I 2L = = = 86.6 A
3 V2 L 3 20000
Ensayo de vaco :
5100
P0 = V2 N I 0 cos 0 cos 0 = = 0.1132
20000
1.3
3
20000
V2 V2 3
RFe = = = = 78465.65
I Fe I 0 cos 0 1.3 0.1132
l = 0.8 m ; r = 0.09 m ; S = 0.0254 m 2
l 0.8
= r = 500 ( Ncleo de hierro); = = 50127.54 H 1
S 500 4 10 0.0254
7
H
0 = 4 10 7
m
N2 I0 50 1.3
N i = 0 = = = 0.0129 Wb
50127.54
Ensayo de cortocircuito :
V 2 ccL 8 20000
cc = 100 = 8 V 2 ccL = = 1600 V
V2 N 100
Pcc 31500
Pcc = 3 V1ccL I 1NL cos cc cos cc = = = 0.1313
3V1ccL I 1 NL 3 5299.3 26.24
sen cc = 0.9913
X cc
X cc = X 1 + X 2' ; X 1 = X 2' = =173.385 L1 = 0.552 H
2
X 2' 173.385
X2 = = = 5.3085 L2 = 0.0169 H
m 2 5.715 2
3.2.7 Cargas.
3.2.7.1 Carga 1.
P = 0.8 MW
VL = 20 kV
cos = 0.87
Conexin en estrella
0.8 10 6
P = 3 VL I L cos I L = = 26.545 A
320 10 3 0.87
P 0.8 10 6
P = 3 R I 2 R = = 378
3 I 2 3 26.545 2
Z = R + jX ; tg = 0.567;
Q
tg =
P
Q 453600
Q = 3 X L I F2 X L = 2
= 214.44
3 I F 3 26.545 2
XL 214.44
L= = 0.6826 H
2 f 2 50
3.2.7.2 Carga 2.
P = 0.2 MW
VL = 20 kV
cos = 0.87
Conexin en estrella
0.2 10 6
P = 3 VL I L cos I L = = 6.64 A
320 10 3 0.87
P 0.2 10 6
P = 3 R I 2 R = = 1514
3 I 2 3 6.64 2
Z = R + jX ; tg = 0.567;
Q
tg =
P
Q 113400
Q = 3 X L I F2 X L = 2
= 857.35
3 I F 3 6.64 2
XL 857.35
L= = 2.73 H
2 f 2 50
3.2.7.3 Carga 3.
P = 1.5 MW
VL = 20 kV
cos = 0.87
Conexin en estrella
1.5 10 6
P = 3 VL I L cos I L = = 49.77 A
320 10 3 0.87
P 1.5 10 6
P = 3 R I 2 R = = 202
3 I 2 3 49.77 2
Z = R + jX ; tg = 0.567;
Q
tg =
P
Q 850500
Q = 3 X L I F2 X L = 2
= 114.45
3 I F 3 49.77 2
XL 114.45
L= = 0.3643 H
2 f 2 50
3.2.7.4 Carga 4.
P = 0.8 MW
VL = 20 kV
cos = 0.98
Conexin en estrella
0.8 10 6
P = 3 VL I L cos I L = = 23.56 A
320 10 3 0.98
P 0.8 10 6
P = 3 R I R =
2
= 480
3 I 2 3 23.56 2
Z = R + jX ; tg = 0.2030;
Q
tg =
P
Q 162400
Q = 3 X L I F2 X L = 2
= 97.52
3 I F 3 23.56 2
XL 97.52
L= = 0.3104 H
2 f 2 50
INPUT
Gate
OUTPUT
Rigbt
VAR
Ron[1..31]
Roff[1..19]
Rigbt
tup
toff
n
m
up
down
dt
INIT
Ron[1]:=1000
Ron[2]:=530
Ron[3]:=260
Ron[4]:=170
Ron[5]:=100
Ron[6]:=65
Ron[7]:=39
Ron[8]:=26
Ron[9]:=18
Ron[10]:=14
Ron[11]:=11
Ron[12]:=9
Ron[13]:=7
Ron[14]:=5.6
Ron[15]:=4.6
Ron[16]:=3.9
Ron[17]:=3.2
Ron[18]:=2.3
Ron[19]:=1.4
Ron[20]:=1.18
Ron[21]:=0.86
Ron[22]:=0.61
Ron[23]:=0.43
Ron[24]:=0.36
Ron[25]:=0.27
Ron[26]:=0.18
Ron[27]:=0.09
Ron[28]:=0.06
Ron[29]:=0.03
Ron[30]:=0.02
Ron[31]:=0.01
Roff[1]:=0.24
Roff[2]:=1.01
Roff[3]:=2.18
Roff[4]:=3.87
Roff[5]:=6.43
Roff[6]:=10.08
Roff[7]:=18.34
Roff[8]:=27.62
Roff[9]:=41.71
Roff[10]:=56.23
Roff[11]:=83.67
Roff[12]:=101.09
Roff[13]:=117.34
Roff[14]:=193.11
Roff[15]:=312.26
Roff[16]:=487.38
Roff[17]:=743.86
Roff[18]:=1822.76
Roff[19]:=100000000
Rigbt:=100000000
tup:=0
toff:=0
up:=0
down:=0
ENDINIT
EXEC
dt:=timestep
tup:=t
up:=1
down:=0
ENDIF
n:=(t-tup)/dt+1
IF n<=31 THEN
Rigbt:=Ron[n]
ELSE
Rigbt:=Ron[31]
ENDIF
ENDIF
toff:=t
down:=1
up:=0
ENDIF
m:=(t-toff)/dt+1
IF m<=19 THEN
Rigbt:=Roff[m]
ELSE
Rigbt:=Roff[19]
ENDIF
ENDIF
ENDEXEC
ENDMODEL
Captulo 4
Simulaciones
Inversor.
Planta fotovoltaica 1.
1 Inversor.
Funcionamiento en isla.
Carga 1.
Planta fotovoltaica 2.
1 Inversor.
Conectada a red.
Carga 1.
Planta fotovoltaica 3.
1 Inversor.
Conectada a red.
Carga 3.
Planta fotovoltaica 4.
1 Inversor.
Mtodo pasivo de deteccin en isla.
Carga 1/Carga 2/Carga 3.
Planta fotovoltaica 5.
1 Inversor.
Mtodo activo de deteccin en isla.
Carga 4.
Planta fotovoltaica 6.
1 Inversor.
Mtodo pasivo de deteccin en isla. Desconexin aguas arriba de la subestacin.
Carga 1.
Planta fotovoltaica 7.
2 Inversores.
Funcionamiento en isla.
Carga 1.
Planta fotovoltaica 8
2 Inversores.
Conectada a red.
Carga 1.
4.1 Inversor.
fpwm=1 kHz
Timestep=10 ns
Representaciones:
ton=2E-7 s 5 MHz
TRANSFORMADA DE FOURIER
ETC
Despus de la falta:
[V]
100
50
-50
-100
-150
-200
-250
0.1378 0.1378 0.1378 0.1378 0.1379 0.1379 [s] 0.1379
500
[V]
375
250
125
-125
-250
-375
-500
0.1472 0.1475 0.1477 0.1479 0.1481 0.1484 [s] 0.1486
f 49.87 Hz
Trabajo Fin de Mster 265 Juan Antonio Rodrguez Gonzlez
Energa Solar Fotovoltaica
[V]
375
250
125
-125
-250
-375
-500
0.1469 0.1469 0.1469 0.1469 0.1470 0.1470 [s] 0.1470
f 49.77 Hz
[V]
1875
1250
625
-625
-1250
-1875
-2500
0.1460 0.1465 0.1470 0.1475 0.1480 0.1485 0.1490 [s] 0.1495
1200
[V]
800
400
-400
-800
-1200
0.1572 0.1573 0.1574 0.1575 0.1576 0.1577 [s] 0.1578
f 50.15 Hz
Despus de la falta:
Captulo 5
Anlisis de resultados y
conclusiones
Se ha puesto de manifiesto la presencia de tres frecuencias en la forma de onda de la
tensin de salida del inversor fotovoltaico:
Cabe destacar que la presencia del filtro inductivo a la salida del inversor reduce
casi totalmente la amplitud de las ondas de tensin a estas estas frecuencias. Sin
embargo, puede que no lo suficiente para dejar de considerar la utilizacin de un
filtro EMI. En este sentido ha de tenerse en cuenta la siguiente consideracin:
= T =
f
100.000 10 3 m / s
= = 6.67 m
15 10 6 1 / s
donde
longitud de onda
velocidad de propagacin de la onda
T periodo
f frecuencia
Para la simulacin del comportamiento de los cables frente a estas elevadas frecuencias
pueden considerarse dos opciones: modelos de parmetros concentrados (circuitos en
cascada) o modelos de parmetros distribuidos.
La representacin de las intensidades en cada lnea nos muestra que la planta fotovoltaica
suministra la potencia activa requerida por las cargas aguas abajo. La red elctrica de
distribucin suministra la potencia reactiva.
Vmin=10.16 kV
Vmax=12.70 kV
Sin embargo, nos hayamos en zona de no deteccin para la carga 1 (0.8 MW,
cos=0.87), ya que la potencia activa suministrada por la planta fotovoltaica (1
MW) es muy parecida a la consumida por la carga. En tal caso la tensin
obtenida en el PCC fue de 11.34 kV.
Para el estudio del mtodo activo de deteccin en isla mediante la inyeccin de una
corriente armnica a la salida del inversor de la planta fotovoltaica, se ha considerado,
como sabemos, una onda de corriente de 1 A y 500 Hz (h=10). Adems, se ha tenido en
cuenta la carga 4 (0.8 MW, cos=0.98), ya que en tal caso la planta fotovoltaica opera en
zona de no deteccin para los mtodos pasivos OVP/UVP.
Captulo 6
Tendencias futuras
Una fuerte industria Europea apoyada por programas nacionales, con exportaciones
significativas y un fuerte mercado interno que tendr un gran retorno social de la
inversin en trminos de crecimiento econmico, creacin de empleo, reduccin de
emisiones de CO2 y,sobre todo, menor dependencia energtica.
Los escenarios de futuro, aunque auguran un muy importante descenso del ritmo de
crecimiento experimentado durante 2009, estn marcados por el alcance de la Paridad
de la Red. Por Paridad de la Red se entiende el punto de competitividad del kWh solar,
sin ningn tipo de ayudas, en relacin con el coste de consumo del kWh. Se puede
alcanzar en los prximos aos en amplias regiones del mundo, naciendo con ello un
mercado gigantesco y una nueva dimensin de la fotovoltaica.
Las tecnologas de capa fina o pelcula delgada irn ganando cuota de mercado de una
forma progresiva durante los prximos aos.
Generar electricidad durante la noche sigue siendo el gran desafo para la casa solar. Un
nuevo catalizador, barato y autorrenovable, descurbierto por el MIT, podra hacer
posible el uso del agua como medio de almacenamiento, lo que permitira tener las luces
encendidas por la noche e incluso recargar un coche elctrico o de hidrgeno. La casa se
convierte as en una central elctrica. Se convierte en una gasolinera.
Bibliografa
Electronic Controlled Device for the Analisis and Design of Photovoltaics Systems.
Pablo Sanchs, Member, IEEE, Alfredo Ursa, Member, IEEE, and Luis Marroyo, Member,
IEEE.
Analysis of Photovoltaic Charging System Based on MPPT. Lixin Pang, Hui Wang, Yuxia
Li, Jian Wang, Ziyu Wang.
A novel Utility Interactive Photovoltaic Inverter with Generation Control Circuit. Hisao
Watanabe, Toshihisa Shimizu, Gunji Kimura.
Toward Giga-watt production of silicon photovoltaic cells, modules and systems.
Takashi Tomita.
A New Topology for Photovoltaic 4 Series DC/DC Converter with High Efficiency
Under Wide Load Range. Jong Pil Lee, Byung Duk Min, Jeong Joong Kim, Kang Yui Ryul,
Tae Jin Kim, Dong Wook Yoo, Eui Ho Song and Ji Yoon Yoo.
Research on Testing Technology of High-power Photovoltaic Arrays. Yinglei He, Qing
Fu, Xiangfeng Li, Xiaogao Chen.
Research of an Evaluation Method of Photovoltaic Output and Generated Energy
Connected with a Power Grid. Atsushi Iga, Tomoyuki Keneko, and Yoshiyuki Ishihara,
Member, IEEE.
Small Photovoltaic Module with Rectangular Cells for Reducing Output Degradation
Caused by Spot Dirt. Shin-ichi Kobayashi, Tomonori Iino, Hironori Kobayashi, Kaxumasa
Yamada, Toshiaki Yachi.
Novel Grid-connected Photovoltaic Generation System. Li-Shiang Lai, Wen-Chieh Hou,
Ya-Tsung Feng, Yu-An Chen.
Electronic Converter for the Analysis of Photovoltaic Arrays and Inverters. P. Sanchs, I.
Echeverra, A. Ursa., O. Alonso, E. Guba and L. Marroyo.
Cloud effects on distributed photovoltaic generation: show transients at the gardner,
Massachussets photovoltaic experiment. Edward C. Kern Jr., Edward M. Gulachenski,
Gregory A. Kern.
Design Considerations and Topology Selection for DC-Module-Based Building
Integrated Photovoltaic System. Bangyin Liu, Chaohui Liang and Shanxu Duan.
Resonant DC-DC Converters for Photovoltaic Energy Generation Systems. A. Bellini, S.
Bifaretti and V. Iacovone.
A simple photovoltaic simulator for testing of power electronics. Ole-Morten Midtgard.
Simulation method for PV Module Power Generation with Dirt Spots and Reduction of
Output Degradation. Shn-ichi Kobayashi, Yuichi Aoyama, Masatoshi Kano, Toshiaki Yachi.
An Optimal Design of a Grid Connected Hybrid Wind/Photovoltaic/Fuel Cell System
for Distributed Energy Production. Debosmita Das, Longya Xu, Dave Nichols.
Apndice A
ATP-EMTP
A.1 Qu es ATP?.
Protegido por las leyes de Estados Unidos, ATP puede ser adquirido bajo licencia,
expedida y garantizada por los propietarios del material.
MS-Windows 9x/NT/2000/XP:
Linux:
Lneas 6000
Elementos lineales 10000
Interruptores 1200
Fuentes 900
Elementos no lineales 2250
Mquinas sncronas 90
Se muestra de esta manera la siguiente ventana de dilogo, compuesta por dos pestaas. La
primera de las cuales hace referencia a las dos versiones del programa ATP que posteriormente se
podrn ejecutar desde ATPCC:
Save on: salvar los resultados obtenidos en un fichero de extensin .lis (disk) o
presentarlos tambin en la pantalla (both).
Output Parameters:
*:se marca esta casilla cuando se desea que el nombre del fichero de
salida coincida con el nombre del fichero .atp utilizado como entrada. En caso
contrario se especifica el nombre del fichero de salida (incluyendo su extensin .lis
o .out).
-R:se selecciona esta opcin si se quiere que, cuando ya exista un fichero .lis
correspondiente al fichero .atp utilizado como entrada, el resultado de una nueva
simulacin se sobreescriba en ese mismo fichero.
Figura. Ventana de dilogo para la configuracin del resto de los programas principales.
Cut extension of input file: elimina la extensin de los ficheros de entrada, antes
de ejecutarlos con el programa adicional.
Default PLOT: se utiliza cuando se pretende que el programa adicional aadido
(por ejemplo: GTPPLOT o PlotXY) acte como programa grfico por defecto. En
este caso, la extensin de los ficheros de entrada debe estar definida como .pl4.
ATP Control Center Working Directory: permite definir el directorio de trabajo del
programa ATPCC.
Save last window position. Permite guardar la posicin y el tamao de la ventana del
programa para la prxima ocasin en que se utilice.
Conviene sealar que la estructura de los ficheros .atp es muy rgida, ya que toda la
informacin del sistema elctrico a simular debe ocupar una posicin especfica en cada fila del
fichero. Para facilitar la escritura de los datos en el fichero se incluye una serie de plantillas
(template). Adems, resulta indispensable tener a mano el manual de usuario ATP-EMTP Rule
Book.
Tal y como se ha comentado, puede adjuntarse una plantilla al editor de textos PFE32.
Ello puede hacerse mediante la opcin Template Attach File. Una vez adjuntada, se puede
utilizar sin ms insertndola en el fichero sobre el que se est trabajando, con la opcin Template
Edit.
Precediendo a las lneas MISCELLANEUS pueden utilizarse una serie de rdenes que
permiten la realizacin de clculos especiales, como, por ejemplo: flujos de carga, clculo de
armnicos, etc. stas pueden ser:
- CABLE CONSTANT
- DATA BASE MODULE
- FREQUENCY SCAN
- HARMONIC FREQUENCY SCAN
- HYSTERESIS
- LINE CONSTANT
- SATURATION
- XFORMER
- ZNO FITTER
- FIX SOURCE (esta opcin permite calcular un flujo de cargas)
Por otra parte, en cualquier lnea del programa pueden especificarse las siguientes tarjetas
especiales:
A continuacin se muestra un ejemplo, el cual ha sido generado por el editor grfico ATPDraw del
circuito elctrico que se muestra a continuacin:
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
C < n 1>< n 2><ref1><ref2>< R >< A >< B ><Leng><><>0
VA XX0002 33. 1. 0
XX0002 33. 1. 0
XX0021VA 33. 1. 0
XX0021 33. 1. 0
XX0002POS .01 1
POS XX0021 1.E3 3
XX0021POS 20. 3
VS XX0025 1. 0
VS XX0025 300. 0
/SWITCH
C < n 1>< n 2>< Tclose ><Top/Tde >< Ie ><Vf/CLOP >< type >
11VA XX0002 0
11 XX0002 0
11XX0021VA 0
11XX0021 0
XX0025VA MEASURING 1
/SOURCE
C < n 1><>< Ampl. >< Freq. ><Phase/T0>< A1 >< T1 >< TSTART >< TSTOP >
14VS 0 167.7 60. -90. -1. 1.
BLANK BRANCH
BLANK SWITCH
BLANK SOURCE
2POS 7.500000E+0001
2XX0021-7.500000E+0001
3POS XX0021 1.500000E+0002
VS
BLANK OUTPUT
BLANK PLOT
BEGIN NEW DATA CASE
BLANK
Fichero 1R.DAT:
El comando $ERASE clarifica el buffer donde se guardar el fichero .pch obtenido. DUM denota
el nombre de los nodos internos del componente (en este caso no hay).
KARD 3 3 3 3
KARG 1 2 3 4
KBEG 3 9 27 80
KEND 8 14 32 80
KTEX 1 1 0 0
$ERASE
/BRANCH
C<Bus M<Bus k<Bus 3<Bus 4<----R<---- L<---- C<----------------------->0
FRNODETONODE VALUER H
C . .^ . . ^ . . ^ . . ^ . . ^ . . ^xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx^
El fichero punch contiene los punteros numricos de todos los argumentos usados en la plantilla.
Igualmente puede obtenerse un fichero de extensin .lib si se escribe en el fichero .dat $PUNCH
name.lib.
El elemento creado, en este caso una resistencia, puede ser utilizado en la simulacin de un
circuito elctrico (.dat) haciendo uso de la orden $INCLUDE. En este caso sera:
$INCLUDE, 1R,FRNODE,TONODE,VALUER,H
H determina el tipo de salida: 0 para no hay salida, 1 para la corriente que circula por el elemento,
2 para la diferencia de potencial entre los terminales del elemento, 3 para disponer de ambas
magnitudes (intensidad y tensin) y 4 para el consumo de potencia y energa.
A.7.3.1 TACS.
El mdulo TACS (Transient Analysis of Control Systems) se puede usar para simular el control
de convertidores HVDC, sistemas de excitacin de mquinas sncronas, funciones de limitacin de
intensidad de pararrayos, cebado de arcos en interruptores y, en general, aquellos dispositivos o
fenmenos que no se pueden modelizar con los componentes elctricos existentes en ATP.
A.7.3.2 MODELS.
Es un lenguaje de programacin que se usa en ATP-EMTP para simular variables
dependientes del tiempo con caractersticas especiales. Los ficheros MODELS se crean de forma
independiente al programa principal y se pueden llamar desde cualquier programa cuantas veces se
desee. Su estructura es muy similar a la de otros lenguajes de programacin, como por ejemplo
FORTRAN, con sentencias del tipo FOR, IF, WHILE, etctera. Han de especificarse las
siguientes declaraciones:
MODEL DEFAULT
INIT -- Initialization
ENDINIT
EXEC -- Execution
ENDEXEC
ENDMODEL
Cada uso de un fichero MODEL se introduce con el comando USE, especificando las
entradas y salidas correspondientes a esa llamada. Los nombres de los elementos definidos en un
fichero MODEL no son visibles desde el exterior y se pueden elegir sin tener en cuenta las
posibles coincidencias con nombres usados fuera de este fichero MODEL.
Una vez construido el circuito, el paso final es dar nombre a los nodos del mismo.
ATPDraw puede automticamente dar nombre a todos los nodos, sin embargo el usuario puede
asimismo hacerlo, ATP Make Names, o dar nombre solo a los nodos de especial inters. Esto
ltimo se realiza simplemente haciendo click con el botn derecho del ratn en el nodo a
nombrar.
Antes de crear el fichero de entrada ATP o ejecutar la simulacin, se han de especificar los
parmetros de la simulacin (parameters miscellaneus), que se encuentran en el men ATP
Settings
Donde:
Print freq: un valor de 500 significa que cada 500 pasos de simulacin, se escribir
en el fichero .lis.
La opcin ATP Make File del men principal genera un fichero ATP del circuito. Podr
guardarse en el directorio que deseemos. Finalmente se procede a su ejecucin mediante la opcin
run ATP del men ATP.
El primer paso para la creacin de un nuevo modelo es la generacin del fichero de texto
donde se describe el nuevo elemento que se quiere modelizar. Este fichero de texto se puede
desarrollar mediante una de las siguientes opciones:
Escribiendo directamente desde el editor de texto, con la ayuda del ATP-EMTP Rule
Book para definir correctamente la estructura del fichero.
Utilizando la opcin Data Base Module (DBM), se quiere que tenga parmetros
variables.
El fichero de texto creado se tiene que compilar mediante alguna de las versiones de ATP,
de modo que se obtenga otro fichero con la extensin .lib. Este ltimo fichero es el que se
almacena en el directorio Usp de ATPDraw.
Por otro lado, tambin hay que definir un icono y una ventana que representen el nuevo
modelo de componente creado mediante el fichero .lib. Para ello, con la opcin Objects User
Specified New sup-file se crea un nuevo componente, cuya extensin es .sup.
En la ventana que aparece a continuacin hay que indicar los datos que definen dicho
componente y el nmero de nudos disponibles para su conexin externa.
El nmero de nodos y datos debe coincidir con las declaraciones ARG y NUM del fichero
escrito mediante Data Base Module. Permite la introduccin de hasta 36 datos, siendo el nmero
mximo de nodos 12.
Kind:
Nodo MODELS:
0: nodo de salida.
1: nodo de entrada de intensidad.
2: nodo de entrada de tensin.
3: nodo de entrada de estado de un interruptor.
4: nodo de entrada de una variable de una mquina.
5: variable TACS.
6: nodo de tensin de la parte imaginaria de estado-fijo (imssv).
7: nodo de intensidad de la parte imaginaria de estado-fijo (imssi).
8: salida de otro model.
Nodo TACS:
0: nodo de salida.
1: nodo de entrada de suma positiva.
2: nodo de entrada de suma negativa.
3: nodo de entrada desconectado.
En la ventana aparecen dos botones. Pulsando el botn donde aparece una lmpara
elctrica, se abre un editor de iconos que permite disear el icono con el que representar el nuevo
componente creado. Mediante el botn en el que aparece un interrogante, se accede a un editor de
ayuda donde se puede escribir el texto que se desea mostrar, cuando se solicite ayuda desde la
ventana de dilogo correspondiente a este componente.
Tras definir todos los datos, se pulsa la opcin save as y se guarda el fichero con extensin
.sup en el subdirectorio Usp, con lo que queda definido el nuevo componente. Para poder utilizar
este nuevo modelo como un componente ms del circuito elctrico a disear, se hace uso de la
opcin User specified Files y se selecciona el fichero .sup correspondiente al componente.
MODEL FLASH_1
comment ****************************************************************
* *
* Function: set or cancel the gap firing control signal *
* Inputs : voltage and current across ZnO resistor *
* Output : the firing signal to the electrical ZnO component *
* *
************************************************************* endcomment
OUTPUT trip
INIT trip:=0
tfire:=0
ENDINIT
EXEC
------------------------------------------------------------------
vcap:=V1-V2
power:=vcap*iczn
energy:=INTEGRAL(power)
------------------------------------------------------------------
IF trip>0 -- is already firing
AND t-tfire>fdur*1.e-3 -- has exceeded firing duration
THEN
trip:=0 -- cancel the firing signal
tfire:=0 -- null the firing time
ENDIF
------------------------------------------------------------------
IF trip=0 -- is not signaling to fire
AND tfire=0 -- firing condition not yet detected
AND ( power >= Pset * 1.e9 -- power setting exceeded
OR energy >= Eset * 1.e6 ) -- energy setting exceeded
THEN
tfire:=t -- set the firing detection time
ENDIF
------------------------------------------------------------------
IF trip=0 -- is not signaling to fire
AND tfire>0 -- firing condition has been detected
AND t-tfire>=fdel*1.e-3 -- firing delay exceeded
THEN
trip:=1 -- set the firing signal
ENDIF
------------------------------------------------------------------
ENDEXEC
ENDMODEL
Igualmente ha de definirse un fichero .sup, asociado a dicho modelo, tal y como se mostr
anteriormente. El nombre de los nodos ha de coincidir con las variables especificadas en las
secciones INPUT y OUTPUT del fichero .mod. Ambos ficheros habrn de almacenarse en la
carpeta /MOD.
Probe Volt
Probe Curr
Probe Tacs
Splitter
T
abla. Line Distributed.
A.8.5.6 SWITCHES.
Tabla. Switches.
A.8.5.7 SOURCES.
Tabla. Sources.
A.8.5.8 MACHINES.
Tabla. Machines.
A.8.5.9 TRANSFORMERS.
Tabla. Transformers.
A.8.5.10 TACS.
Coupling to circuit
TACS sources
Fortran statement
TACS devices
Al abrir un nuevo fichero aparece la ventana de dilogo de la figura, donde se eligen las
variables a representar y el tipo de representacin deseado.
Las variables escogidas se pueden dibujar en funcin del tiempo o en funcin de una de
ellas, segn se seleccione la opcin correspondiente. En el segundo caso, la primera de las variables
seleccionada constituye el eje x y aparece sealizada con una X a su izquierda.
di
Tensin en una inductancia: v = L
dt
dv
Corriente en un condensador: i =C
dt
dy
x=
dt
t
y (t ) = xdt + y(t t )
t t
t
y (t ) = y (t t ) + ( x(t ) + x(t t ))
2
Se resuelve por consiguiente esta ecuacin para x(t), que puede representar la intensidad a travs
de un condensador o la diferencia de potencial entre los terminales de una inductancia.
2
x(t ) = x(t t ) + ( y (t ) y (t t ))
t
2
x(t ) =
t
2 2 2
x(t + t ) = + (1.0 1.0) =
t t t
2
x(t + 2t ) =
t
Se aprecia en este caso la existencia de una oscilacin, como aparece en la siguiente figura. El valor
medio de x(t) corresponde a la forma de onda deseada.
Una transicin de y(t) de dos pasos de tiempo en lugar de uno resolvera esta situacin, como se
muestra a travs del siguiente conjunto de ecuaciones. Se asume en este caso y(t-t)=0, y(t)=0.5 e
y(t+t)=1.0.
2 1
x(t ) = 0.5 =
t t
x(t + 2t ) = 0.0 +
2
(1.0 1.0) = 0.0
t
x(t + 3t ) = 0.0
...
En la siguiente figura se ilustra los valores de x(t) e y(t) en consideracin de estas ecuaciones:
De todo ello, se deduce que reducir el paso de tiempo permite en algunos casos paliar o
eliminar la aparicin de oscilaciones numricas. Se debe sealar sin embargo que la aparicin de
oscilaciones debidas a operaciones de conmutacin no se vern reducidas de esta manera, ya que el
interruptor ideal siempre cambiar en un paso de tiempo.
A.11.3 Ejemplos.
Interruptores electrnicos de potencias con grandes reas de seguridad (Safe Operating Area,
SOA) tales como IGBTs y MOSFETs no necesitan snubbers. En estos casos el usuario del
programa puede aadir snubbers numricos. El valor del condensador puede estar comprendido
entre 1 y 2 nF, si no se considera la resistencia. Por el contrario, si se considera una rama RC
snubber, como se muestra en la figura siguiente, la constante de tiempo =RC habr de ser 2
3 veces el tiempo de paso elegido para la simulacin. Ntese, sin embargo, que los valores de R y
C podrn variar segn el sistema elctrico a simular.
Llegados a este punto se considera la simulacin del paso de corriente a travs de una
inductancia presentado anteriormente teniendo en cuenta una rama snubber. En la siguiente
figura se muestra los resultados de tensin e intensidad obtenidos, donde se aprecia la ausencia de
oscilaciones.
A.11.4.4 Amortiguamiento.
Alternativamente, puede considerarse la introduccin de resistencias de amortiguamiento
en el circuito elctrico. En el caso de una inductancia, la resistencia habra de disponerse en
paralelo con la misma, como se muestra en la siguiente figura, proporcionando un camino
alternativo a la intensidad cuando el interruptor se abre. La intensidad en el circuito RL decaer as
con el tiempo.
L
=
Rd
L
Rd =
t
Igualmente puede considerarse para el caso del condensador. En este caso la resistencia se
coloca en serie con el condensador, como se muestra en la siguiente figura. Para la eleccin de la
constante de tiempo puede considerarse la siguiente expresin:
Ks t
Rd =
2C
A.11.4.5 Interpolacin.
El esquema de interpolacin detecta el cambio de paso de la intensidad y la tensin cuando
se presentan oscilaciones. Entonces considera el tiempo anterior a la misma y realiza una
interpolacin lineal a la intensidad cero o en la operacin del interruptor. Las ecuaciones
diferenciales se resuelven en este punto y se comienza de nuevo en el siguiente tiempo de paso.