Título: Balanceo de los circuitos de distribución primaria 32 y 131 de Santa Clara Autor: Domingos Isabel Filipe Tutor: Dr. C. Ignacio Pérez Abril , Junio del 2018 , Junio del 2018 , Junio del 2019 Departamento de Electroenergética Electroenergetic Department , Junio del 2019 Author: Domingos Isabel Filipe Title: Balance of primary distribution circuits 32 and 131 of Santa Clara Thesis Director: Dr. C. Ignacio Pérez Abril i Este documento es Propiedad Patrimonial de la Universidad Central “Marta Abreu” de Las Villas, y se encuentra depositado en los fondos de la Biblioteca Universitaria “Chiqui Gómez Lubian” subordinada a la Dirección de Información Científico Técnica de la mencionada casa de altos estudios. Se autoriza su utilización bajo la licencia siguiente: Atribución- No Comercial- Compartir Igual Para cualquier información contacte con: Dirección de Información Científico Técnica. Universidad Central “Marta Abreu” de Las Villas. Carretera a Camajuaní. Km 5½. Santa Clara. Villa Clara. Cuba. CP. 54 830 Teléfonos.: +53 01 42281503-1419 ii PENSAMIENTO Para tener éxito, es necesario encontrar algo a lo que aferrarse, algo que te motive, algo que te inspire. (Tony Dorsett) iii DEDICATORIA Esta tesis está dedicada a mis padres Mario Filipe y Isabel Augusta quienes con su amor, paciencia y esfuerzo me han permitido llegar a cumplir hoy un sueño más, gracias por inculcar en mi el ejemplo de esfuerzo y valentía, de no temer las adversidades porque Dios está conmigo siempre. A mis hermanos por su cariño y apoyo incondicional, durante todo este proceso, por estar conmigo en todo momento gracias. A toda mi familia porque con sus oraciones, consejos y palabras de aliento hicieron de mí una mejor persona y de una u otra forma me acompañaron en todos mis sueños y metas. A mi novia, en el camino encuentras personas que iluminan tu vida, que con su apoyo alcanzas de mejor manera tus metas, a través de sus consejos, de su amor, y paciencia me ayudó a concluir esta meta. iv AGRADECIMIENTOS Son muchas las personas que han contribuido al proceso y conclusión de este trabajo. En primer lugar, quiero agradecer a la Empresa Eléctrica de Santa Clara que me pudieron ayudar en los datos necesarios de este proyecto. Agradezco al profesor Reinier Herrera Casanova quien me brindó un apoyo incondicional de manera personal y me alentó para que concluyera esta investigación. v RESUMEN Dentro de las redes trifásicas más desbalanceadas se encuentran los circuitos de distribución primaria, esto se debe fundamentalmente a la distribución de cargas monofásicas no equilibradas a lo largo de la red y al uso de ramales de una o dos fases y neutro. Si el circuito opera en condiciones de desbalance, se incrementan las pérdidas de energía en los alimentadores primarios y se produce la circulación de elevadas corrientes por el conductor neutro en condiciones normales de operación, lo que puede provocar la operación incorrecta de las protecciones de sobrecorriente de tierra. El objetivo del presente trabajo, consiste en reducir la circulación de corriente por el conductor neutro y minimizar las pérdidas de energía a lo largo del circuito de distribución primaria, aplicando el programa para el balance de fases basado en el algoritmo genético por ordenamiento no dominado NSGA-II, en dos circuitos de la ciudad de Santa Clara, realizando las reconexiones mínimas necesarias de los ramales de una o dos fases y neutro y de los transformadores de distribución monofásicos. vi ÍNDICE PENSAMIENTO ................................................................................................................. ii AGRADECIMIENTOS ....................................................................................................... iv RESUMEN ......................................................................................................................... v INTRODUCCIÓN ............................................................................................................... 1 CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCIÓN PRIMARIA............... 4 1.1 Principales definiciones de desbalance ................................................................ 4 1.1.1 Definición de la IEC ....................................................................................... 5 1.1.2 Definición de la NEMA .................................................................................. 6 1.1.3 Definición de la IEEE .................................................................................... 6 1.2 Normas referentes al grado de desbalance .......................................................... 7 1.3 Principales causas del desbalance ...................................................................... 7 1.4 Principales efectos del desbalance ...................................................................... 8 1.5 Características fundamentales de las redes de distribución primaria en el mundo y en Cuba ....................................................................................................................... 8 1.6 Situación actual del desbalance en los circuitos analizados ................................. 9 CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA ....................................................................................................................... 12 2.1 Introducción ....................................................................................................... 12 2.2 Formulación del problema .................................................................................. 14 2.2.1 Variables independientes ............................................................................ 14 2.2.2 Funciones objetivo ...................................................................................... 17 2.3 Implementación del NSGA-II .............................................................................. 20 2.3.1 Representación de las variables (cromosoma) ............................................ 20 2.3.2 Cálculo de las funciones objetivo ................................................................ 20 2.4 Procedimiento de estimación de cargas ............................................................. 21 2.4.1 Ajuste de las cargas .................................................................................... 22 2.4.1.1 Transformadores trifásicos ............................................................................ 23 2.4.1.2 Bancos de tres transformadores Y-∆ ............................................................ 23 vii 2.4.1.3 Bancos de dos transformadores ................................................................... 24 2.4.1.4 Transformadores monofásicos ...................................................................... 25 2.4.2 Algoritmo de ajuste de las cargas ............................................................... 26 2.4.3 Algoritmo de estimación .............................................................................. 27 CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS ......................................................................................................................................... 29 3.1 Introducción ....................................................................................................... 29 3.2 Análisis del Circuito 32 ....................................................................................... 30 3.2.1 Estimación de cargas y balanceo del circuito 32 ......................................... 31 3.3 Análisis del Circuito 131 ..................................................................................... 35 3.3.1 Estimación de cargas y balanceo del circuito 131 ....................................... 37 CONCLUSIONES ............................................................................................................ 40 RECOMENDACIONES .................................................................................................... 41 REFERENCIAS BIBLIOGRÁFICAS ................................................................................ 42 ANEXOS .......................................................................................................................... 46 Anexo I: Planos atualizados de los circuitos ................................................................. 46 Anexo II: Base de datos para el circuito 32................................................................... 47 Anexo III: Base de datos para el circuito 131 ................................................................ 48 INTRODUCCIÓN 1 INTRODUCCIÓN El Sistema Electroenergético Nacional (SEN) es el encargado de brindar el suministro de energía eléctrica a los consumidores de forma continua dentro de los parámetros de calidad establecidos. El SEN está compuesto de cuatro subsistemas o partes fundamentales: generación, transmisión, subtransmisión y distribución, esta última a su vez se divide en distribución primaria y distribución secundaria y es la parte del sistema más afectada por los desbalances de carga. El desbalance tanto de voltaje como de corriente trae consigo un incremento de las pérdidas de energía y la circulación de elevadas corrientes por el conductor neutro en condiciones normales de operación, provocando disparos de los dispositivos de protección. Lograr el balance en los circuitos de distribución primaria contribuye a brindar un mejor servicio a los usuarios en cuanto a calidad y continuidad del servicio. Un balanceo adecuado de las cargas por fase en el circuito de distribución durante las 24 horas del día y en cada uno de los nodos de éste, constituye la condición de máxima eficiencia en un circuito trifásico, pero desde el punto de vista operativo, es prácticamente imposible de lograr. El uso de transformadores monofásicos, si bien tiene ciertas ventajas desde el punto de vista de la explotación de las redes en cuanto a inversiones y flexibilidad, al poder conformar bancos de diferentes características, contribuye de forma importante a acentuar el fenómeno del desbalance de carga entre las fases del circuito. Por otra parte, la práctica de explotación de estos circuitos y el continuo incremento de la carga, provocan la adición de nuevos bancos de transformadores o la modificación de los ya existentes, lo que corrientemente se hace sin atender a consideraciones de balance en el circuito [1]. El balance de las redes de distribución primaria se puede considerar como un problema muy actual y se relaciona fundamentalmente con la necesidad del ahorro de energía y la mejora de la calidad del servicio eléctrico. En varios de los circuitos de Santa Clara, se aprecian grandes porcientos de desbalance, esto provoca un incremento de las pérdidas de energía en el circuito de distribución primaria y contribuye al desbalance de los voltajes en los distintos nodos. Además, aparece la circulación de una elevada corriente por el conductor neutro en condiciones normales de operación que complica la detección de las fallas a tierra a partir de protecciones de sobrecorriente de neutro lo que conlleva a la desconexión de este tipo de protección. Con la ayuda de las mediciones tomadas por los interruptores principales de los circuitos 26 y 170 del municipio Santa Clara, se pudo conocer la situación real en que se encuentra el desbalance en los mismos. Para dar INTRODUCCIÓN 2 solución a este problema se deben seleccionar de manera óptima las fases del circuito primario a que se deben conectar los distintos bancos de transformadores de distribución, así, como los ramales bifásicos y monofásicos, de manera que se minimicen las pérdidas de energía y se logre un mejor balance de las fases en toda la longitud del circuito. Por lo que se plantea como problema de investigación: ¿Cómo lograr el balanceo efectivo de los circuitos de distribución primaria en Santa Clara? De acuerdo al problema de investigación, se establece como objetivo general:  Proponer variantes de balanceo de los circuitos de distribución primaria 32 y 131 de la ciudad de Santa Clara. Para dar cumplimiento al mismo se trazan como objetivos específicos:  Describir los métodos de balanceo de los circuitos de distribución primaria.  Obtener los datos necesarios a partir del celaje de los circuitos y actualizar sus bases de datos.  Realizar la estimación y el ajuste de las cargas.  Aplicar el programa de balanceo a varios circuitos y seleccionar las variantes más adecuadas. Las tareas científicas que contribuyen al cumplimiento de los objetivos específicos son:  Obtención de los datos necesarios mediante el celaje de los circuitos de distribución primaria.  Recopilación de las mediciones de los interruptores principales de cada circuito (NULEC).  Actualización de las bases de datos de los circuitos.  Estimación de las cargas a partir de las mediciones disponibles.  Modelación de los circuitos en el Radial.  Determinación de las variantes de balanceo más adecuadas. El informe de la investigación se estructura de la siguiente manera: En la introducción se dejará definida la necesidad, importancia y actualidad del tema que se aborda y se dejarán explícitos los elementos del diseño teórico. En el primer capítulo se realiza un acercamiento al problema del desbalance y sus efectos negativos, las características de los circuitos de distribución primaria y se describe la situación real en que se encuentra el desbalance en los circuitos analizados. En el segundo capítulo se describen las principales características del programa de balanceo mediante NSGAII y se aborda lo referente al programa de estimación de cargas y en el tercer capítulo se describe brevemente cada uno de los circuitos, se le aplica el balanceo a cada uno de ellos y INTRODUCCIÓN 3 posteriormente se muestran los resultados obtenidos. Finalmente, se presentan las conclusiones de la investigación desarrollada y las recomendaciones para trabajos futuros. CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 4 CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCIÓN PRIMARIA En el presente capítulo se realiza un acercamiento al problema del desbalance de las redes de distribución primaria y a la situación actual en que se encuentra este tema. También se analizan algunas normas existentes sobre esta temática. 1.1 Principales definiciones de desbalance El desbalance de una variable trifásica es una condición en la que las tres fases presentan diferencias en el módulo y/o desplazamientos angulares entre fases distintas de 120 grados eléctricos. La variable trifásica indicada de manera genérica puede ser indistintamente, la tensión fase-neutro, la tensión de línea o la corriente de un sistema trifásico [2]. Los circuitos de distribución son los encargados de recibir la energía eléctrica de las líneas de transmisión o subtransmisión y entregarlas a los usuarios a la tensión adecuada, con la conveniente continuidad y calidad de suministro para los distintos usos. Es en esta parte donde se producen los porcentajes más grandes de pérdidas de energía en todas sus manifestaciones debido al gran volumen de elementos que los conforman (alimentadores, transformadores, líneas y cables, capacitores y equipos de protección), y a los bajos niveles de tensión que se manejan [2, 3]. Aunque no todos los circuitos de distribución tienen las mismas configuraciones y longitudes, sí comparten características comunes. En la Figura 1.1 se muestra un típico circuito de distribución [4]. CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 5 Figura 1.1. Típica subestación de distribución con uno de los varios alimentadores mostrados. (Algunos de los taps laterales no se muestran). Los sistemas de distribución son típicamente radiales, donde el flujo de potencia nace solo de un nodo. El nodo principal se reconoce como la subestación que alimenta al resto de la red, en la misma se reduce la tensión del nivel de alta tensión al de media tensión [5]. La distribución se hace en el nivel de media tensión o en baja tensión. Los clientes residenciales o comerciales se alimentan en baja tensión y los clientes industriales, de acuerdo a sus requerimientos particulares, se alimentan en media tensión o baja tensión. Se consideran baja tensión los valores menores que 1 kV, media tensión a los comprendidos entre 1kV y 35 kV, y alta tensión a valores superiores a 35 kV [5]. 1.1.1 Definición de la IEC La IEC (International Electrotechnical Commission) de acuerdo a lo expresado en su norma 61000-4-30 [6], define el desbalance de tensión utilizando el método de las componentes simétricas, como la magnitud de la relación entre la componente de secuencia negativa con respecto a la componente de secuencia positiva expresado en por ciento. 𝑉𝑈𝐹 = 𝑉2 𝑉1 ∗ 100 (1.1) Donde: CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 6 𝑉2- voltaje de secuencia negativa 𝑉1- voltaje de secuencia positiva El voltaje 𝑉2 se determina mediante la siguiente expresión: 𝑉2 = 100 ∗ √ 1−√3−6𝛽 1+√3−6𝛽 (1.2) β se calcula como: β = 𝑈𝑎𝑏𝑓𝑢𝑛𝑑4+𝑈𝑏𝑐𝑓𝑢𝑛𝑑4+𝑈𝑐𝑎𝑓𝑢𝑛𝑑4 (𝑈𝑎𝑏𝑓𝑢𝑛𝑑2+𝑈𝑏𝑐𝑓𝑢𝑛𝑑2+𝑈𝑐𝑎𝑓𝑢𝑛𝑑2)2 (1.3) Donde las tensiones 𝑈𝑖𝑗𝑓𝑢𝑛𝑑 corresponden a la componente fundamental de la tensión entre la fase 𝑖 y la fase 𝑗. 1.1.2 Definición de la NEMA La definición de la NEMA (National Electric Manufacturing Association) [7] sobre desbalance de voltaje, conocida también como índice de desbalance de voltaje de línea (LVUR), es dada por: %LVUR = 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 𝑚á𝑥𝑖𝑚𝑎 𝑑𝑒𝑙 𝑣𝑜𝑙𝑡𝑎𝑗𝑒 𝑑𝑒 𝑙í𝑛𝑒𝑎 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑣𝑜𝑙𝑡𝑎𝑗𝑒 𝑑𝑒 𝑙í𝑛𝑒𝑎 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 ∗ 100 (1.4) Donde la desviación máxima del voltaje de línea promedio se determina mediante la siguiente expresión: 𝑑𝑒𝑠𝑣 𝑚á𝑥 𝑉𝑙í𝑛 𝑝𝑟𝑜𝑚=Máx[|Vab − Vprom|, |Vbc − Vprom|, |Vca − Vprom|] (1.5) El voltaje de línea promedio se determina como: 𝑉𝑙í𝑛𝑒𝑎 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 = Vab+Vbc+Vca 3 (1.6) Esta definición aportada por la NEMA solo trabaja con magnitudes, los ángulos de fase no se toman en cuenta. 1.1.3 Definición de la IEEE La definición de la IEEE (Institute of Electrical and Electronics Engineers) [8] sobre desbalance de voltaje, conocida también como índice de desbalance de voltaje de fase (PVUR), es dada por: CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 7 %𝑃𝑉𝑈𝑅 = 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 𝑚á𝑥𝑖𝑚𝑎 𝑑𝑒𝑙 𝑣𝑜𝑙𝑡𝑎𝑗𝑒 𝑑𝑒 𝑓𝑎𝑠𝑒 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑣𝑜𝑙𝑡𝑎𝑗𝑒 𝑑𝑒 𝑓𝑎𝑠𝑒 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 ∗ 100 (1.7) La IEEE usa la misma definición que la NEMA, diferenciándose en que la IEEE usa voltajes de fase en lugar de voltajes de línea a línea. Los ángulos de fase son despreciados y se consideran solamente las magnitudes. 1.2 Normas referentes al grado de desbalance La norma europea EN 50160 [9] establece que, tanto para baja tensión como para media tensión, en condiciones normales de explotación, para cada período de una semana, el 95% de los valores eficaces promediados en 10 minutos de la componente inversa de la tensión de alimentación debe situarse entre el 0% y el 2% de la componente directa. Se establece la salvedad, en baja tensión que, en algunas regiones equipadas con líneas parcialmente monofásicas o bifásicas, los desequilibrios pueden alcanzar el 3% en los puntos de suministros trifásicos. La norma IEC 61000-2-2 [10] al igual que la norma IEEE 1159, recomienda que el índice de desbalance en un suministro eléctrico no debe superar el 2%. En la norma IEC 61000-4-30 [11] se establece que el período de medida debe ser de una semana con valores cada 10 minutos y/o cada 2 horas. Evidentemente estos límites que establece la normativa internacional son adoptados generalmente por las reglamentaciones de calidad de servicio vigentes en cada país. La norma cubana NC 365: 2011 toma todos los elementos aplicables de la norma IEC 60038: 2002, dicha norma establece que los sistemas de suministro eléctrico deberán ser diseñados y operados para limitar el desbalance máximo de tensión (fase con mayor desviación con respecto a la tensión promedio / tensión promedio de las fases) al 2% medido en los terminales de entrega del suministro en régimen de vacío [12]. 1.3 Principales causas del desbalance Lograr el balance perfecto de tensiones y corrientes en las redes de distribución primaria es técnicamente inalcanzable. El motivo principal para la existencia de redes trifásicas desbalanceadas lo constituye la distribución de cargas monofásicas no equilibradas a lo largo de la red, dichas cargas varían en el tiempo dependiendo de la hora del día [13]. Existen otras causas que provocan desbalances en estas redes y deben ser consideradas, entre estas se encuentran: la transposición incompleta de líneas, transformadores CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 8 conectados en delta-abierta o estrella-abierta, impedancias asimétricas en las redes de alimentación, fallas monofásicas, desperfectos en empalmes, uniones o contactos y asimetrías propias de las fuentes de suministro eléctrico [14]. Pueden existir también otros elementos, como son: hornos monofásicos, hornos de inducción y hornos trifásicos de arco eléctrico, máquinas de soldadura eléctrica, aparatos de rayos x, anomalías en el sistema eléctrico tales como: apertura de un conductor, falla en el aislamiento de los equipos, corrientes de magnetización de transformadores trifásicos debido a las características magnéticas propias de su construcción y banco trifásico de capacitores con una fase fuera de servicio [15]. 1.4 Principales efectos del desbalance Los voltajes desbalanceados pueden dar lugar a efectos nocivos sobre los equipos y sobre el sistema de energía, que es intensificado por el hecho de que un desbalance pequeño en los voltajes de fase puede causar un desbalance desproporcionado más grande en las corrientes de fase [16]. El sistema de energía, bajo condiciones de desbalance, incurrirá en mayores pérdidas y efectos térmicos y será menos estable, porque cuando las fases son equilibradas, el sistema está en una posición mejor para responder a las transferencias de la carga de emergencia. Además, bajo estas condiciones aparece una corriente de desbalance que circula por el conductor neutro en condiciones de operación normal, esto complica la detección de las fallas a tierra a partir de protecciones de sobrecorriente de tierra [17, 18]. El efecto del desbalance de voltaje puede también ser severo en equipos tales como motores de inducción, convertidores electrónicos de potencia y accionamientos de velocidad variable (ASD). 1.5 Características fundamentales de las redes de distribución primaria en el mundo y en Cuba Las redes de distribución primaria ocupan un lugar importante en el sistema electroenergético, siendo su función tomar la energía eléctrica de la fuente y distribuirla o entregarla a los consumidores. La efectividad con que las redes de distribución realizan esta función se mide en términos de regulación de voltaje, continuidad del servicio, flexibilidad, eficiencia y costo. El costo de las redes de distribución representa aproximadamente el 50% del costo del sistema eléctrico en su conjunto. Estas redes deben ser capaces de brindar servicio con un mínimo de variaciones de voltaje y un mínimo de http://www.ecured.cu/index.php?title=Redes_de_distribuci%C3%B3n&action=edit&redlink=1 CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 9 interrupciones, y en caso de que ocurran interrupciones del servicio, estas deben ser de corta duración y afectar al menor número posible de consumidores. Las redes de distribución presentan características muy particulares que los diferencian de las redes de transmisión. Entre estas se distinguen:  topologías radiales,  múltiples conexiones,  cargas de distintas naturalezas,  alta razón R/X (líneas de resistencias comparables a la reactancia), líneas sin transposiciones [19]. 1.6 Situación actual del desbalance en los circuitos analizados Los registros históricos obtenidos a partir de los interruptores principales de dos circuitos de distribución primaria del municipio de Santa Clara (circuitos 32 y 131), muestran la existencia de elevados niveles de desbalance en los mismos. Esto, como se explicó anteriormente, dificulta la operación adecuada de los circuitos, por lo tanto, se hace necesario reducir estos niveles de desbalance [1]. Con los datos obtenidos a partir de los recerradores NULEC, se puede determinar el desbalance definido por la NEMA (Desbalance NEMA o %NEMA) y el porciento de corriente por el conductor neutro (%In) que presentan los circuitos analizados: Desbalance NEMA = Máx[|Ia−Iprom|,|Ib−Iprom|,|Ic−Iprom|] Iprom ∗ 100 (1.8) Iprom = Ia+Ib+Ic 3 (1.9) % In = In Iprom ∗ 100 (1.10) Donde: Máx[|Ia − Iprom|, |Ib − Iprom|, |Ic − Iprom|] es el valor máximo de la diferencia entre las corrientes de fase y la corriente promedio. In- corriente por el conductor neutro. Con el uso de las expresiones anteriores y los datos obtenidos a partir de los interruptores principales, se puede determinar el desbalance definido por la NEMA y el porciento de corriente por el conductor neutro, para diferentes estados de operación de los circuitos. Los resultados obtenidos se muestran en la tabla 1.1. Tabla 1.1 Valores para el desbalance definido por la NEMA y el %In. CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 10 Circuito Estado Ia (A) Ib (A) Ic (A) In (A) %NEMA %In 32 Medio 44,30 39,00 35,90 8,50 11,50 7,20 Máximo 70,00 58,00 55,00 18,00 14,80 29,50 131 Medio 74,40 83,10 81,60 22,20 8,65 27,85 Máximo 86 112 106 41 15,13 40,46 En las figuras 1.1 y 1.2 se puede apreciar el comportamiento que presentan en un día promedio el desbalance definido por la NEMA (%NEMA) y el porciento de corriente por el conductor neutro (%In) en los dos circuitos de distribución primaria analizados. Figura 1.1: Comportamiento del %NEMA y %In en el circuito 32. CAPÍTULO I. DESBALANCE EN LAS REDES DE DISTRIBUCION PRIMARIA 11 Figura 1.2: Comportamiento del %NEMA y %In en el circuito 131. En las figuras anteriores se puede observar que en los circuitos de distribución primaria analizados los niveles de desbalance son bastante elevados y se encuentran muy por encima de los niveles recomendados por las diferentes instituciones internacionales. 0 5 10 15 20 25 30 35 40 45 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Desbalances en el circuito 131 %Nema %In CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 12 CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA En el presente capítulo, se realiza un acercamiento a las principales características del programa de balanceo de circuitos de distribución primaria mediante el algoritmo genético NSGA-II, y debido a su gran importancia, se describe el procedimiento utilizado para realizar la estimación de las cargas eléctricas suministradas por los distintos transformadores del circuito. 2.1 Introducción En los circuitos de distribución primaria pueden aplicarse dos estrategias relacionadas, con el objetivo de reducir las pérdidas de energía y la circulación de corriente por el conductor neutro, estas son: la reconfiguración y el balance de fases. La reconfiguración consiste en el cambio de la topología del circuito mediante la apertura o cierre de desconectivos, mientras que el balance de fases se refiere a la redistribución de las cargas entre las fases del circuito para mejorar el equilibrio entre estas [20]. Los circuitos de distribución primaria se encuentran entre las redes trifásicas más desbalanceadas, el desbalance entre las corrientes de fase provoca dos problemas fundamentales: el incremento de las pérdidas en los conductores primarios y la presencia de una alta corriente de neutro en operación normal que dificulta la detección de las fallas a tierra por los dispositivos de protección. Debe tenerse en cuenta que la corriente de desbalance no se manifiesta de igual forma en toda la extensión del circuito. Un circuito puede tener un balance perfecto en su extremo de la fuente y estar muy desbalanceado en otras secciones o ramales, de esta forma, las pérdidas debidas al desbalance pueden ser altas en un circuito aparentemente balanceado [20]. Los circuitos de distribución siempre están protegidos en la subestación, no obstante, otras protecciones pueden ubicarse aguas abajo en algunas secciones o ramales del circuito. Adicionalmente, un circuito, o algunas secciones del mismo, pueden ser suministrados desde otro circuito adyacente en condiciones de emergencia. En todos los puntos donde exista una protección de falla a tierra, se debe reducir al mínimo la corriente de desbalance considerando todas las condiciones de operación posibles que pueden cambiar la topología CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 13 del circuito. Normalmente, la reconexión de un reducido número de ramales o de transformadores de distribución puede mejorar grandemente el balance de fases en los puntos de interés y reducir las pérdidas de energía en los conductores primarios [20]. Varias técnicas de optimización como: redes neuronales artificiales [21, 22], lógica difusa [23], evolución diferencial [24], colonia de hormigas [25], etc. se han utilizado para la reconfiguración de la distribución mediante la operación de los desconectivos. También, el balance de fases se ha resuelto mediante diferentes métodos como: programación entera mixta [26], simulación del recocido [27, 28], algoritmo genético (AG) [29], algoritmos heurísticos [30, 31], algoritmo inmune [32, 33], etc. El AG de la referencia [29] busca el balance de fases y la reducción de las pérdidas mediante la optimización de las conexiones de los transformadores de distribución. Sin embargo, el número de reconexiones necesarias resulta muy alto para ser práctico. La referencia [30] presenta un algoritmo heurístico para la reconexión de los ramales de dos fases y neutro o de una fase y neutro con el objetivo de reducir las pérdidas y mejorar el balance en varios puntos del circuito. Además, otra heurística de búsqueda hacia atrás se propone en [31] para ajustar las conexiones de los ramales, pero solo se considera la minimización de un índice de desbalance en cada segmento de línea. En las referencias [32, 33] se propone un algoritmo inmune para reconectar los ramales y transformadores de distribución. La función multiobjetivo considera el desbalance de las corrientes, el costo de las interrupciones del servicio y el costo de la mano de obra para implementar las soluciones. Un sistema experto para obtener la estrategia de balance de fases se propone en [34]. Esta herramienta trata de reducir la corriente de neutro y consecuentemente evitar los disparos indeseables de la protección de falla a tierra. Las referencias [35, 36] proponen una herramienta de optimización combinatoria basada en lógica difusa, redes neuronales y un método heurístico rápido para balancear las corrientes. La aproximación propuesta en [37] efectúa el balance de fases mediante la introducción de nuevas conexiones en los devanados de los transformadores de distribución. El balance de fases se realiza en [38, 39] por: optimización por enjambre de partículas, forraje de bacterias y técnicas difusas. La función multiobjetivo considera: la corriente de neutro, el costo de reconexión, las caídas de tensión y las pérdidas. En la referencia[40] se investiga el uso de simulación del recocido caótica, mientras que en [41] y [42] se emplea una técnica auto-adaptativa de evolución diferencial y un algoritmo heurístico híbrido CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 14 respectivamente. Seis algoritmos para el balance de fases se estudian en [43], donde se sugiere a la programación dinámica como la técnica más apropiada para este problema. En [44] se propone el uso de un transformador balanceador que reduce el desbalance y mitiga las corrientes de armónicos en el neutro de los sistemas de distribución secundaria. Como se ha explicado, varios métodos han sido aplicados para resolver el balance de fases. Sin embargo, no todos los trabajos referidos consideran todas las condiciones pertinentes, ni emplean un método eficiente de optimización multiobjetivo que permita obtener la frontera de Pareto del problema [45]. A partir de un trabajo previo sobre este problema [46] se desarrolló una aplicación computacional basada en el algoritmo genético por ordenamiento no-dominado (NSGA-II) que ahora se extiende para minimizar el desbalance en múltiples puntos y para todas las posibles condiciones de operación de la red eléctrica. Este trabajo formula el problema de balance de fases como la minimización de: 1) la corriente de desbalance en los puntos deseados del circuito; 2) las pérdidas de energía en los conductores primarios; y 3) el número de elementos que tienen que reconectarse para lograr estos objetivos [47]. 2.2 Formulación del problema Para formular el problema de optimización del balance de fases, deben definirse las variables independientes y las funciones objetivo del problema. 2.2.1 Variables independientes Las variables independientes de este problema son las conexiones a las fases del circuito primario de: los ramales de dos fases y neutro, los ramales de una fase y neutro, los bancos trifásicos de dos o tres transformadores monofásicos y los transformadores monofásicos de distribución. La determinación de las corrientes en los terminales primarios de los distintos bancos de transformadores o transformadores de distribución se desarrolla seguidamente. Los bancos de tres transformadores (3T) generalmente utilizan la conexión (Y-Δ) como se muestra en la figura 2.1. CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 15 Figura 2.1: Banco de transformadores con conexión Y-Δ. La carga monofásica de 120/240V (S1φ) se alimenta por el transformador de alumbrado mientras que dos transformadores de fuerza iguales completan el banco para suministrar la carga trifásica de 240 V (S3φ). Para simplificar el análisis, la carga monofásica se considera distribuida entre los dos secundarios del transformador de alumbrado. De la figura 2.1, las tensiones de línea en el secundario están en fase con las tensiones de fase en el primario (V1n, V2n y V3n). Entonces, aplicando el principio de superposición, las corrientes primarias (I1, I2, I3) se obtienen mediante:            *333 1* 113 1 3 * 233 1* 113 1 2 * 133 1* 113 2 1 nn nn nn VSVSI VSVSI VSVSI       (2.1) Los bancos de dos transformadores (2T) utilizan la conexión (Y abierta - Δ abierta) como se muestra en la figura 2.2. Figura 2.2: Banco de transformadores con conexión Y abierta - Δ abierta. La carga monofásica de 120/240 V se alimenta del transformador de alumbrado mientras que un solo transformador de fuerza completa el banco para suministrar la carga trifásica CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 16 de 240 V. Las corrientes del primario (I1, I2, I3) se obtienen de:       03 * 33233 1 2 * 11 * 33133 1 1    I VSVSI VSVSVSI nn nnn   (2.2) Los transformadores monofásicos de distribución (1T) suministran solo carga monofásica de 120/240 V. En este caso solo existe el terminal 1 y:   0 0 3 2 * 111    I I VSI n (2.3) En la tabla 2.1 se muestra un resumen para el cálculo de las corrientes en los terminales primarios de los distintos bancos de transformadores. Tabla 2.1 Corrientes en el primario de los transformadores. Corrientes Banco de tres transformadores Banco de dos transformadores Transformador monofásico 1I * 1 3 3 1 * 1 1 3 2              nn V S V S  * 1 1 * 3 3 1 3 3 1              nnn V S V S V S  * 1 1       nV S  2I * 2 3 3 1 * 1 1 3 1              nn V S V S  * 3 3 2 3 3 1        nn V S V S  0 3I * 3 3 3 1 * 1 1 3 1              nn V S V S  0 0 En forma general, el circuito de distribución contiene (Nlat) ramales de dos fases y neutro y de una fase y neutro, así como (Ntrf) bancos trifásicos de transformadores monofásicos y transformadores monofásicos de distribución cuyos terminales 1, 2, 3 pueden conectarse según resume la tabla 2.2. La posición de cada fase se corresponde con los terminales 1, 2 y 3 en cada caso. Se debe tener cuidado al reconectar los bancos de dos transformadores para no modificar la secuencia de fases en el secundario. Además, si un transformador monofásico está conectado a una sección de fase y neutro no hay que considerar su posible reconexión. Tabla 2.2 Posibles conexiones de ramales y transformadores. Fases disponibles en el circuito CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 17 Elementos ABC AB BC CA Ramales de dos fases y neutro AB, BC, CA Ramales de una fase y neutro A, B, C A, B B, C C, A Bancos de tres transformadores ABC, BCA, CAB Bancos de dos transformadores AB, BA, BC, CB, CA, AC AB, BA BC, CB CA, AC Transformadores monofásicos A, B, C A, B B, C C, A Para representar las conexiones de los ramales y los transformadores de distribución, se define el arreglo (xcon) de (Nlat+Ntrf) elementos que está formado por los subvectores (xlat) y (xtrf). ][ trflatcon xxx  (2.4) Los elementos de (xcon) son enteros restringidos a los valores mostrados en la tabla 2.2. Desde el punto de vista práctico, se requiere que el número de reconexiones sea limitado a un porciento reducido de los ramales y transformadores. Sin embargo, todos los (Nlat+Ntrf) elementos pueden modificarse en el proceso de optimización, lo cual debe evitarse de alguna manera. Para reducir el número de reconexiones, se utiliza el arreglo adicional (xvar) formado por (Nvar) enteros acotados entre 1 y la dimensión del arreglo (xcon). La función del arreglo (xvar) es determinar el subconjunto de los elementos del arreglo (xcon) que pueden variar su valor original. Es decir, solo los elementos de (xcon) seleccionados en el arreglo (xvar) pueden variar su conexión original. El número máximo de reconexiones (Nvar) se escoge por el usuario. Añadiendo los elementos de (xvar), las variables independientes del problema se representan por el arreglo (x) de (Nlat+Ntrf+Nvar) elementos formado por los subvectores (xcon) y (xvar). ][ var xxx con (2.5) 2.2.2 Funciones objetivo Diversos objetivos pueden lograrse mediante la reconexión de los ramales y transformadores de distribución. En este trabajo se van a considerar tres objetivos fundamentales para la optimización: CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 18 1) Mínimas corrientes de neutro en los puntos deseados. 2) Mínimas pérdidas de energía. 3) Mínimo número de reconexiones. Para calcular las funciones objetivo de este problema, se considera que la red eléctrica radial (que puede estar formada por varios circuitos interconectados), tiene un estado de operación normal (t=0) y varios estados posibles de contingencia (t=1...S). Las corrientes de neutro deben reducirse al mínimo en todos los estados de operación de la red, para posibilitar la correcta operación de las protecciones de fallas a tierra en todas las condiciones. Sin embargo, las pérdidas de energía se van a determinar solo en el estado de operación normal, ya que las contingencias son de corta duración y ocurren esporádicamente. Para un cierto vector solución (x) (conexión de los elementos) y para el estado de operación (t), las corrientes primarias de todos los ramales, bancos de transformadores y transformadores de distribución se calculan por las expresiones de la tabla 2.1. Además, las corrientes primarias de otros elementos como transformadores trifásicos o bancos de capacitores trifásicos se calculan y se añaden a los resultados de las corrientes por las fases determinados anteriormente. A la hora (h), la corriente (Ji,k,h) en la fase (i) de la sección (k) del circuito se calcula como la suma de las corrientes primarias (Ii,n,h) del conjunto (Ck) de todas las cargas suministradas desde esta sección del circuito.    kCn hnihki txItxJ ),(),( ,,,, (2.6) La corriente de neutro (JNk,h) es la suma de las corrientes de las (Mk) fases de la sección (k).    kM i hkihk txJtxJN 1 ,,, ),(),( (2.7) Mientras que el máximo de la corriente de neutro en dicha sección (k) se define como:   241 ),(),( ,max  hwheretxJtxJ hkh N k Nmax (2.8) Las pérdidas de potencia a la hora (h) en la sección (k) se calculan mediante: CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 19 2 , 1 2 ,,, ),(),(),( txJRtxJRtxP hkkN M i hkikhk N k    (2.9) Donde (Rk) y (RNk) son las resistencias de los conductores de fase y de neutro respectivamente. Como se considera la minimización de la corriente de neutro en varios puntos y para varios estados de operación de la red, se ha diseñado una arquitectura abierta para el problema formulado que permite la declaración de varias funciones objetivo (1…W) para minimizar las corrientes de neutro. Cada función determina el máximo de un conjunto de corrientes de neutro en diferentes puntos y diferentes condiciones de operación. Cada función objetivo se define como:   TitKikwheretxJxJ ktki NmaxNfun  , ),()( ,max (2.10) Donde los conjuntos: (Ki) y (Ti) representan las localizaciones y estados de operación que se consideran para evaluar la función. El segundo objetivo a minimizar son las pérdidas de energía en el estado normal (t = 0), o sea, la suma de las pérdidas en las N secciones del circuito para las 24 horas del ciclo de carga diario:     N k h hk xPxE 1 24 1 , )0,()( (2.11) Finalmente, el tercer objetivo a minimizar es el número (NR(x)) de elementos a reconectar. De esta forma, el problema de optimización se expresa por (W+2) funciones objetivo como sigue:                )()( )()( )()( )()( min 2 1 11 xNxf xExf xJxf xJxf RW W WW Nfun Nfun  (2.12) Los elementos del vector (x), son enteros acotados por las posibles conexiones de los elementos (xcon) (como se establece en la tabla 2.2), así como acotados por los límites para el arreglo auxiliar (xvar). CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 20 2.3 Implementación del NSGA-II El problema presentado es un problema de optimización combinatoria que requiere un algoritmo especial de optimización multiobjetivo como el NSGA-II. Este es un algoritmo muy efectivo que puede optimizar simultáneamente varias funciones de tipo conflictivo para obtener el conjunto de soluciones óptimas del problema llamado frontera de Pareto. Una descripción simplificada de este método se reproduce seguidamente. 1) Se genera una población inicial aleatoria de tamaño N. 2) La población se ordena mediante un procedimiento de ordenamiento no-dominado. 3) Se repiten los siguientes pasos hasta obtener un número de generaciones deseado. a. Se seleccionan por torneo los padres de la nueva generación. b. Se obtienen los descendientes por cruzamiento y mutación de los padres. c. Se forma una generación intermedia de tamaño 2N que incluye los padres y los descendientes. d. La nueva generación de tamaño N se obtiene por ordenamiento y selección de la generación intermedia. Como todo algoritmo genético, la solución de un problema por medio del NSGA-II implica algunas adaptaciones y la programación de ciertas partes del algoritmo. En este caso, se ha adaptado una versión del NSGA-II en Matlab codificado en números reales [48] para resolver el problema presentado. Además, se han modificado los operadores genéticos para incluir el cruzamiento y mutación de variables de tipo entero. 2.3.1 Representación de las variables (cromosoma) Debido a que la versión del NSGA-II utilizada emplea un cromosoma codificado en números reales; no hay necesidad de codificar las variables. De esta forma, el cromosoma definido para este problema contiene un arreglo de enteros capaz de representar el vector x de las variables que se define en la ecuación (2.5). 2.3.2 Cálculo de las funciones objetivo Una parte especial del NSGA-II que debe programarse para resolver el problema presentado es el procedimiento que calcula todas las funciones objetivo. Este procedimiento se invoca por el NSGA-II cada vez que se genera un nuevo individuo (solución) por medio de cruzamiento o mutación. En esta implementación, el procedimiento referido considera tensión nominal (módulo) en todos los nodos. Las corrientes primarias de todos los bancos de transformadores y CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 21 transformadores de distribución se calculan mediante las expresiones de la tabla 2.1. Con estas corrientes se determinan las corrientes de fase y de neutro en todas las secciones del circuito y se evalúan las pérdidas. Este procedimiento reduce el esfuerzo de cálculo para evaluar las soluciones y si se realiza una buena caracterización de las cargas, se obtienen resultados suficientemente precisos. 2.4 Procedimiento de estimación de cargas La base para cualquier tipo de estudio a realizar en un circuito de distribución primaria, es el conocimiento con la mayor precisión posible de las cargas eléctricas suministradas por los distintos bancos de transformadores o transformadores monofásicos del circuito. La carga se caracteriza normalmente por su magnitud pico (kW y kVAr) y por su comportamiento horario (gráfico de carga). Además, debido a las características de los bancos de transformadores que se emplean en estos circuitos, para determinar su corriente de carga en cada fase del circuito primario, es preciso conocer en todo momento la magnitud de la carga monofásica y trifásica del banco [45]. La única forma de conocer la magnitud de las cargas es medirlas, e incluso las mediciones solo serán válidas para el momento en que estas se realizan. Sin embargo, los circuitos de distribución en Cuba carecen de instrumentación para monitorear la carga en los bancos de transformadores de distribución y la experiencia dice que normalmente no se dispone de tomas de carga de los mismos [45]. De esta manera, se puede concluir que comúnmente no existen mediciones actualizadas de las cargas y se hace casi imposible realizarlas antes de hacer cualquier tipo estudio que se requiera. Esto lleva muchas veces a los especialistas a estimar las cargas en los diferentes bancos de transformadores (potencia activa y reactiva máximas de la carga monofásica y trifásica, así como el comportamiento horario de dicha carga). Esta estimación se realiza mediante diferentes criterios que pueden ser más o menos exactos y que muchas veces introducen errores apreciables. Además, en el caso de estudios en que es importante la caracterización de la carga por fase, como es el caso de los análisis para el balance de fases, se complica más aún esta estimación [45]. Normalmente, todos los circuitos de distribución primaria utilizan recerradores NULEC o dispositivos ION que monitorean constantemente el comportamiento del circuito, por lo cual se disponen de mediciones de los parámetros totales del circuito: potencia activa, reactiva y aparente, factor de potencia, tensiones, corrientes de fase y corriente de neutro entre otros. Por lo tanto, aunque no se dispone de mediciones en los bancos de transformadores, CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 22 sí se conocen los valores de los parámetros totales del circuito y su comportamiento horario. Estas mediciones totales pueden emplearse como base de comparación para determinar si una estimación dada de las cargas puede estar cerca o lejos de la realidad [45]. A partir de esta consideración, se desarrolla una herramienta computacional que dadas las mediciones totales del circuito en un día característico y contando con una base de datos actualizada de la estructura del mismo y las potencias de los distintos bancos de transformadores, pueda determinar una estimación de las cargas que produzca resultados compatibles con las mediciones disponibles. Evidentemente, esta estimación solo puede considerarse como una aproximación posible a la realidad y puede estar más o menos distante de los datos reales del circuito [45]. 2.4.1 Ajuste de las cargas A partir de las expresiones de la tabla 2.1, las cuales se utilizan para el cálculo de las corrientes en los terminales primarios de los distintos bancos de transformadores, pueden determinarse las potencias aparentes por fase multiplicando estas corrientes por las tensiones correspondientes. En la tabla 2.3 se muestran las expresiones para el cálculo de dichas potencias. Tabla 2.3 Potencias en el primario de los transformadores. Potencias Banco de tres transformadores Banco de dos transformadores Transformador monofásico * 111 IVS n   33 1 13 2 SS   1 3 1 33 1 1 S V V S n n        1S * 222 IVS n   33 1 1 2 13 1 S V V S n n               n n V V S 3 2 33 1 1 0 * 333 IVS n   33 1 1 3 13 1 S V V S n n        0 0 En un circuito que suministra un conjunto de bancos de transformadores o transformadores monofásicos de distribución, la suma de todas las potencias por fase de dichas cargas debe igualar a la potencia total en cada fase del circuito (Stot) si se desprecian las pérdidas. De esta manera, para un tiempo (t), en cada fase (k), debe cumplirse que:    N i tkitk SStot 1 ,,, (2.13) Donde N es el número de cargas del circuito. CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 23 Si no se tiene ninguna indicación sobre la potencia que puede estar llevando cada transformador, es razonable suponer que esta sea proporcional a dichas potencias, con lo cual se puede obtener una primera aproximación de la potencia de carga de cada transformador, de acuerdo a la fase primaria a que está conectado, es decir:          N i kikitktki StrfStrfStotS 1 ,,,,, (2.14) Esta expresión permite determinar un conjunto de potencias que sumadas igualan a la potencia medida en cada fase en el tiempo (t). Sin embargo, cada tipo de banco de transformadores tiene cierta distribución de la carga entre las fases, que es propia de las conexiones de dicho banco de transformadores según se determina en la tabla 2.3. Por lo tanto, es necesario desarrollar un procedimiento para que dadas las potencias estimadas en los distintos transformadores del banco (S1, S2 y S3), pueda determinar unos valores de la carga monofásica (S1φ) y trifásica (S3φ) del banco, tales que las potencias por fase en el primario se acerquen a los valores estimados, pero siempre cumpliendo las relaciones propias de cada tipo de carga. En el desarrollo de las expresiones de estimación se va a considerar igual factor de potencia para la carga trifásica y la monofásica. Además, para evitar la sobrecarga de las expresiones, en el análisis siguiente se van a obviar los subíndices (i,t) que representan la carga y el tiempo para los que se hace el ajuste. 2.4.1.1 Transformadores trifásicos En este caso están tanto los transformadores trifásicos como los bancos de transformadores con conexión (Y-Y) que se emplean para suministrar carga balanceada. Para el caso de estos transformadores, solo existe carga trifásica que se calcula como la suma de las cargas de cada fase:     3 1 3 1 0 k kSS S   (2.15) Obtenida la carga trifásica, la nueva estimación de carga en cada fase es: 33 1 321 SSSS  (2.16) 2.4.1.2 Bancos de tres transformadores Y-∆ Este banco tiene tanto carga monofásica como trifásica y si no hay una indicación de la CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 24 carga de cada transformador, se puede considerar que la carga se reparta proporcional a la potencia de los transformadores, es decir, que operen a igual coeficiente de carga. A partir de la tabla 2.3 y suponiendo que los tres transformadores tienen una carga igual a su potencia nominal, las potencias por fase tienen que cumplir:    33 1 1 3 13 1 3 33 1 1 2 13 1 2 33 1 13 2 1 S V V SStrf S V V SStrf SSStrf n n n n                (2.17) Además, la suma de las potencias de los tres transformadores tiene que igualar a la suma de las cargas trifásicas y monofásicas, por lo que puede llegarse a:   31 33 1 13 2 2 SSSS SSS FA A   (2.18) Resolviendo este sistema de ecuaciones se obtiene que: AF FA SSS SSS   4 )(2 3 1   (2.19) La relación entre la carga monofásica y la trifásica para que el banco esté plenamente cargado o que el coeficiente de carga sea igual para los tres transformadores es: )4()(231 AFFA SSSSSS   (2.20) Obtenida esta relación, puede efectuarse la estimación de la siguiente manera:     13 11 )/12(3 SS SS   (2.21) Esta forma de estimación da iguales coeficientes de carga, pero solo respeta la carga en el transformador de alumbrado. Otra forma que pudiera usarse para estimar la carga que no conduce a igual coeficiente de carga de los tres transformadores, pero respeta ambas cargas en alumbrado y fuerza sería tomar como potencia en los transformadores de fuerza al promedio de la carga en las fases 2 y 3 y calcular: ))(2())(2( 13232131 SSSSSSSS   (2.22) 2.4.1.3 Bancos de dos transformadores A partir de la tabla 2.3 y suponiendo que los dos transformadores tienen una carga igual a su potencia nominal, las potencias por fase tienen que cumplir: CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 25               n n n n V V SStrf S V V SStrf 3 2 33 1 2 1 3 1 33 1 1 1 1   (2.23) Como se quiere igualar el coeficiente de carga de ambos transformadores se tiene que cumplir la igualdad: F n n A n n S V V S S S V V S                3 2 33 1 1 3 1 33 1 11  (2.24) Si se sustituye (S1φ = α S3φ) se obtiene: F n n A n n S V V S V V                3 2 3 1 3 1 3 1 11  (2.25) Cuya solución para alfa es:                             n n n n n n F A V V real V V imag V V S S 3 1 2 3 1 2 3 2 2 3 1 111 (2.26) Una vez obtenida la relación entre la carga monofásica y trifásica correspondiente a las potencias nominales de ambos transformadores, se determinan las cargas como:       13 3111 1)3/()/1( SS VVSS nn   (2.27) Esta forma de estimación da iguales coeficientes de carga, pero solo respeta la carga en el transformador de alumbrado. La otra forma posible de hacer esta estimación es calcular las cargas a partir de las cargas estimadas para ambos transformadores:               n n n n V V SSS V V SS 3 1 33 1 11 3 2 23 1 13   (2.28) 2.4.1.4 Transformadores monofásicos En este caso, la carga monofásica es igual a la carga del transformador directamente: 03 11     S SS (2.29) CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 26 Es necesario destacar, que el factor de potencia de todas las cargas monofásicas y trifásicas se considera igual al de la carga total medida para el circuito, por lo cual este factor de potencia se impone a todas las cargas calculadas. 2.4.2 Algoritmo de ajuste de las cargas Para realizar el ajuste de las cargas a las mediciones totales del circuito para un tiempo (t), que corrientemente puede ser el horario pico, es preciso considerar que no todas las cargas tienen igual comportamiento horario, por lo que pueden existir bancos de transformadores que tienen una carga reducida en dicho tiempo. Para resolver este problema, se parte de considerar que el gráfico de carga que se estima represente cada carga del circuito. Si (Pgt) es la potencia en por unidad del gráfico en el tiempo (t), las potencias nominales de los transformadores del banco se afectarán por dicho coeficiente: tktk PgStrfStrf , (2.30) De esta forma se obtienen los valores de potencia nominal referidas al tiempo (t), que son las que se emplearán para distribuir la carga en dicho tiempo. Así, si un transformador no tiene carga en este tiempo o es muy baja, recibirá cero o un valor reducido de la carga medida que se repartirá principalmente entre los bancos cuya carga tiene el pico en un tiempo cercano o igual a (t). Utilizando estos valores como una aproximación inicial de las potencias aparentes por fase en cada carga se realiza el siguiente algoritmo hasta lograr la convergencia. end S SS SSStotS nn errorSStotwhile StrfS n n jki N i n tki n tkitk n tki N i n tkitk tki n tki 2.2 Tabla lapor nuevasobtienen se 2.29)-(2.15 sexpresione las utilizando ,obtienen se 1 max 0 ,, 31 1 1 ,, 1 ,,,,, 1 ,,, ,,,,                  (2.31) Una vez terminado el algoritmo de ajuste, se determinan los valores máximos de las cargas monofásicas y trifásicas por el proceso inverso al utilizado para adaptar las potencias nominales al tiempo (t), es decir: CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 27 tt tt PgSS PgSS       3max3 1max1 (2.32) 2.4.3 Algoritmo de estimación Para realizar la estimación de las cargas en un circuito cuya carga es homogénea (todas las cargas tienen igual gráfico horario), normalmente basta con ajustar la carga en el horario pico para obtener buenos resultados. Sin embargo, la mayoría de los circuitos están compuestos por cargas heterogéneas y si se ajusta para el horario pico de un tipo de carga por ejemplo residencial, el ajuste para las cargas no residenciales no es bueno. Para intentar resolver este problema, el algoritmo de estimación realiza el ajuste de la carga a dos horas diferentes (t1) y (t2) que se corresponden con el horario pico nocturno y el horario pico del circuito y se determina (t1) como el horario pico entre las 6:00 pm y las 12:00 pm y (t2) el pico en el resto del día. Escogidos los gráficos de carga que se van a emplear en los diferentes bancos de transformadores, se realiza un ajuste de la carga para las mediciones de (t1), de donde se obtienen las cargas monofásicas y trifásicas para este horario S1φ(t1), S3φ(t1), así como otro ajuste en el horario (t2) para obtener S1φ(t2), S3φ(t2). Se supone que el mejor ajuste integral en todo el horario del día se obtiene por unos valores de carga monofásica y trifásica que están dados por una combinación lineal de los valores de carga obtenidos por el ajuste en ambos horarios, es decir: )()1()( )()1()( 23133 21111 tStSS tStSS       (2.33) Donde (λ) es una constante entre 0 y 1 que puede determinarse sencillamente por búsqueda directa. Hasta este momento se ha realizado la determinación de las cargas con dos objetivos fundamentales: 1) Mejor ajuste a la corriente de neutro. 2) Mejor ajuste a las corrientes de fase. Esta metodología permite automatizar la estimación de las cargas en los circuitos de distribución asegurando un buen ajuste a las mediciones registradas. Sin embargo, no hay que olvidar que esta estimación puede estar lejos de las cargas reales si no se cumplen los presupuestos considerados. Si se cuenta con tomas de carga en los bancos de transformadores, mediciones de energía en los transformadores monofásicos o los registros actualizados de facturación a los consumidores de cada banco o transformador, pueden CAPÍTULO II: PROGRAMA PARA EL BALANCEO DE CIRCUITOS DE DISTRIBUCIÓN PRIMARIA 28 reducirse los posibles errores. Por otra parte, la selección correcta de los gráficos de carga para los diferentes bancos de transformadores o transformadores monofásicos es fundamental [45]. CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 29 CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS En el presente capítulo, se realiza un análisis de las principales características de cada uno de los circuitos seleccionados y se explica brevemente lo referente a la estimación de cargas en cada uno de ellos. Posteriormente se muestran los resultados emitidos por el programa de balanceo y la descripción de la variante seleccionada. 3.1 Introducción Para lograr un balanceo efectivo de los circuitos de distribución primaria, se deben considerar algunas cuestiones fundamentales. En primer lugar, se debe realizar un celaje del circuito que permita conocer y actualizar todos los datos del mismo, luego se introduce toda esta información en el programa Radial y se confecciona una base de datos, posteriormente se procede a realizar un proceso de estimación y ajuste de las cargas del circuito. Después de obtenidos todos estos datos se corre el programa de balanceo y se selecciona la variante más adecuada para cada circuito, por último, deben aplicarse estos resultados en la práctica. Dentro de las cuestiones analizadas anteriormente, el proceso de celaje de los circuitos posee una gran importancia, para su correcta realización deben tenerse en cuenta los pasos siguientes: 1) Identificar la disposición de las fases a la salida de la subestación con la ayuda de los datos que se visualizan en el interruptor principal de cada circuito. Se debe tener presente esta disposición durante todo el recorrido e identificar a que fase(s) del tronco del circuito se conectan cada uno de los ramales. 2) Identificar el calibre de los conductores de fase y del conductor neutro, así como la distancia aproximada de nodo a nodo. 3) Verificar la potencia de cada transformador e identificar el tipo de conexión, la fase a la que se conecta, identificar los transformadores de fuerza y de alumbrado (en el caso de bancos de dos o tres transformadores), el tipo de carga que alimenta cada banco y por ende el gráfico de carga característico que le corresponde. CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 30 4) Caracterizar los bancos de capacitores, de los cuales se debe verificar su potencia y el estado en que se encuentran. 3.2 Análisis del Circuito 32 El circuito 32 tiene su origen en la subestación (RF) ubicada en la carretera a Manicaragua, se alimenta a través de un transformador de 1600 kVA, con conexión (Δ-Y aterrada) y relación de voltajes de 34,5 - 13,8 kV. Las líneas que conforman el tronco del circuito ocupan una extensión de 2,87 km aproximadamente. El circuito cuenta con 25 transformadores monofásicos con conexión entre fase y neutro, 8 bancos de dos transformadores con conexión estrella abierta-delta abierta y 3 bancos de tres transformadores con conexión estrella-delta, para un total de 36 bancos de transformadores conectados en el circuito, además posee 4 ramales bifásicos y 9 ramales monofásicos. La carga predominante en este circuito es la residencial, pero también existen cargas estatales de importancia, entre estas se encuentran: la torre de televisión, los talleres del Estado Mayor de las FAR, una imprenta, la base de taxis, la base de ómnibus Transmetro, la Empresa de Comunales, la Empresa Cimex y la Escuela de Preparación para la Defensa. En la tabla 3.2 se muestra la cantidad de transformadores que se encuentran conectados en cada fase, las potencias instaladas por fase y la potencia total instalada en el circuito. Tabla 3.2 Potencias instaladas en el circuito 32. Potencias instaladas Fase A Fase B Fase C Cantidad de transformadores por fase 16 18 16 Potencia instalada por fase (kVA) 1019,5 1072 977 Potencia total instalada en el circuito (kVA) 3068,5 La operación del circuito 32 es controlada en su totalidad por el interruptor NULEC (V-224), los registros históricos obtenidos a partir de dicho interruptor permiten conocer el comportamiento de los parámetros fundamentales del circuito en un período de tiempo determinado. Las figuras 3.1 y 3.2 muestran el comportamiento de las corrientes y de las potencias de este circuito en el período de tiempo comprendido entre el 11 y el 17 de marzo de 2019. En estas figuras se observa que el circuito analizado presenta un gráfico mixto- residencial, hay evidencias de que la fase más cargada es la (a) y la menos cargada es la (c). Se debe señalar que, aunque los valores de corriente de neutro parecen pequeños, en realidad para este circuito se pueden considerar como significativos pues las corrientes de fase también son pequeñas, por lo tanto, es obvio que existe un desbalance considerable. CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 31 Figura 3.1: Comportamiento de las corrientes en el circuito 32. Figura 3.2: Comportamiento de las potencias en el circuito 32. 3.2.1 Estimación de cargas y balanceo del circuito 32 Para realizar cualquier tipo de estudio en un circuito de distribución es necesario conocer con la mayor precisión posible las cargas eléctricas suministradas por los distintos bancos de transformadores o transformadores monofásicos del circuito. Al no existir mediciones actualizadas de las cargas y la imposibilidad de realizarlas, se hace necesario estimar las cargas en los diferentes bancos de transformadores (potencia activa y reactiva máximas de la carga monofásica y trifásica, así como el comportamiento horario de dicha carga). Con los datos obtenidos en el celaje realizado y los registros históricos descargados del interruptor principal del circuito, se procede a realizar la estimación de las cargas eléctricas. Para esto se utiliza un programa de estimación elaborado en el software MATLAB, este programa trabaja con las mediciones aportadas por el interruptor principal para un día 11-03-2019 12-03-2019 13-03-2019 14-03-2019 15-03-2019 16-03-2019 17-03-2019 11-03-2019 12-03-2019 13-03-2019 14-03-2019 15-03-2019 16-03-2019 17-03-2019 CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 32 promedio (se analiza el comportamiento del circuito durante el mes de marzo de 2019 y se escoge como día promedio el 15 de marzo de 2019), con una base de datos confeccionada a partir del celaje realizado y con una base de datos que contiene los posibles gráficos de carga para cada uno de los bancos de transformadores en función de la carga servida. Disponiendo de todos estos elementos, se corre el programa de estimación y se selecciona el tipo de ajuste deseado, es decir: mejor ajuste a las corrientes de fase o mejor ajuste a la corriente de neutro. Una vez seleccionado el tipo de ajuste, el programa calcula las potencias monofásicas y trifásicas (P1, Q1, P3 y Q3) de los distintos bancos de transformadores y estima un gráfico para cada corriente de fase y para la corriente que circula por el conductor neutro. Evidentemente, esta estimación solo puede considerarse como una aproximación posible a la realidad y puede estar más o menos distante de los datos reales. Las potencias monofásicas y trifásicas (P1, Q1, P3 y Q3) que se obtienen en este proceso se muestran en el Anexo II. La estimación realizada para este circuito, utilizando un ajuste por corriente de neutro, se muestra en la figura 3.3, puede verse que el ajuste de las cargas es bueno respecto a los gráficos medidos en un día característico por el NULEC. Figura 3.3: Estimación realizada para el circuito 32. Al concluir con el proceso de estimación y ajuste de las cargas, se conforma una base de datos con estos elementos y con los datos recopilados en el celaje. Posteriormente se carga esta base de datos desde el programa de balanceo, se ajustan los datos referentes al voltaje nominal del circuito, los máximos cambios permitidos y se seleccionan los elementos que se desean variar (ramales bifásicos, ramales monofásicos, bancos de tres transformadores, bancos de dos transformadores y transformadores monofásicos) , se ajustan los puntos CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 33 donde se desea minimizar la corriente de neutro y por último se selecciona el tamaño de la población y la cantidad de generaciones deseadas. Luego se procede a realizar varias corridas del programa de balanceo mediante el algoritmo genético de ordenamiento no dominado NSGA-II. En el caso de este circuito se decide variar solamente los ramales de una y dos fases y los transformadores monofásicos, se utiliza una población de 100 individuos, se realizan 100 generaciones y se permiten 8 cambios como máximo. Los resultados iniciales emitidos por el programa de balanceo para este circuito, en sus condiciones originales, muestran que las pérdidas de energía en las líneas primarias alcanzan los 13,68 kWh/día y que la corriente de neutro representa un 26,26 % de las corrientes de fase. Las variantes de balanceo que se obtienen se encuentran ordenadas por corriente de neutro, pérdidas de energía y número de reconexiones respectivamente. Se debe escoger una variante, que, con un pequeño número de reconexiones, pueda lograr una buena reducción de la corriente de neutro y de las pérdidas de energía, atendiendo a esto se seleccionó la variante # 45. En la figura 3.4 se muestran las variantes de balanceo obtenidas para este circuito y los resultados de la variante seleccionada. Figura 3.4: Variantes de balanceo para el circuito 32. Al analizar los resultados obtenidos para la variante seleccionada, puede verse que la corriente de desbalance debe experimentar una disminución de un 19,92 %, las pérdidas de energía deben reducirse en 2,20 kWh/día y solo deben efectuarse tres cambios. Se debe señalar que, aunque los ahorros de energía que se consiguen en este caso son pequeños, CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 34 debe notarse que las pérdidas de energía del circuito también lo son. En la figura 3.5 se muestran los puntos del circuito donde deben efectuarse dichos cambios. Figura 3.5: Cambios propuestos para el circuito 32. En la figura 3.6 se observan los gráficos estimados por el programa de balanceo para las corrientes de fase y de neutro, antes y después de que se ejecuten las reconexiones. CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 35 Figura 3.6: Corrientes estimadas antes y después de efectuar las reconexiones. Se observa en estos resultados teóricos la disminución que debe experimentar la corriente de neutro y como tienden a valores muy cercanos las corrientes por las tres fases. Por tanto, se reducen las pérdidas de potencia y energía en toda la longitud del circuito primario y se pueden ajustar las protecciones de sobrecorriente de tierra, lo que permite elevar la fiabilidad del circuito y mejorar la calidad del servicio eléctrico que se presta a los clientes. 3.3 Análisis del Circuito 131 El circuito 131 parte de la subestación de la Universidad. Se alimenta a través de un transformador de 1600 kVA con conexión (Δ-Y aterrada) y relación de voltajes de 34,5 - 4,16 kV. Este circuito tiene una longitud aproximada de 7 km y cuenta con 5 ramales bifásicos, 3 ramales monofásicos, 4 bancos de tres transformadores, 2 transformador triásico, 7 bancos de dos transformadores y 6 transformadores monofásicos. La mayor parte de las cargas alimentadas por este circuito son monofásicas, aunque existen otros tipos de cargas, por lo que su gráfico de carga es mixto-residencial. Entre las cargas estatales más importantes de este circuito se encuentran: - Correo de la Universidad - Facultad de Construcciones y Facultad de Ciencias Agropecuarias - Almacén de la Universidad - CEDER - Pollera - Cámara 1,2 y 3 CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 36 En la tabla 3.2 se muestra la cantidad de transformadores que se encuentran conectados en cada fase, las potencias instaladas por fase y la potencia total instalada en el circuito. Tabla 3.2 Potencias instaladas en el circuito 131. Potencias instaladas Fase A Fase B Fase C Cantidad de transformadores por fase 12 11 13 Potencia instalada por fase (kVA) 400 462.5 447.5 Potencia total instalada en el circuito (kVA) 1310 La operación del circuito 131 es controlada por el interruptor NULEC (V-1321), los registros históricos obtenidos a partir de dicho interruptor permiten conocer el comportamiento de los parámetros fundamentales del circuito en un período de tiempo determinado. En las figuras 3.7 y 3.8 se muestra el comportamiento de las corrientes y las potencias de este circuito en el período de tiempo comprendido entre el 1 y el 7 de abril de 2019. Figura 3.7: Comportamiento de las corrientes en el circuito 131. Figura 3.8: Comportamiento de las potencias en el circuito 131. 1-04-2019 2-04-2019 3-04-2019 4-04-2019 5-04-2019 6-04-2019 7-04-2019 1-04-2019 2-04-2019 3-04-2019 4-04-2019 5-04-2019 6-04-2019 7-04-2019 CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 37 3.3.1 Estimación de cargas y balanceo del circuito 131 Con los datos obtenidos en el celaje realizado y los registros históricos descargados del interruptor principal se procede a realizar la estimación de las cargas del circuito. Como se explicó en el caso del circuito anterior, para realizar esta estimación se debe disponer de las lecturas del interruptor principal en un período de tiempo determinado para entonces seleccionar un día promedio, para este circuito se analiza su comportamiento en el período de tiempo comprendido entre el 17 de marzo del 2019 y el 17 de abril del 2019 y se seleccionó como día promedio el miércoles 3 de abril del 2019. Las potencias monofásicas y trifásicas (P1, Q1, P3 y Q3) que se obtienen en este proceso se muestran en el Anexo III. La estimación realizada para este circuito, con un ajuste por corriente de neutro, se muestra en la figura 3.9, puede verse que el ajuste de las cargas es bueno respecto a los gráficos medidos en un día característico por el NULEC. Figura 3.9: Estimación realizada para el circuito 131. Cuando se concluye con el proceso de estimación y ajuste de las cargas, se conforma una base de datos con estos elementos y con los datos recopilados en el celaje. Posteriormente se carga esta base de datos desde el programa de balanceo, se ajustan todos los elementos necesarios y se procede a realizar varias corridas del programa. En el caso de este circuito se decide variar solamente los ramales de una y dos fases y los transformadores monofásicos, se utiliza una población de 100 individuos, se realizan 100 generaciones y se permiten 8 cambios como máximo. Además, se decide realizar el balanceo en dos puntos diferentes, el primero de ellos corresponde a la salida de la CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 38 subestación y el segundo es un punto del circuito que se encuentran un gran número de ramificaciones monofásicas y bifásicas. Los resultados iniciales emitidos por el programa de balanceo para este circuito, en sus condiciones originales, muestran que las pérdidas de energía en las líneas primarias alcanzan los 72,34 kWh/día y que la corriente de neutro en el primer nodo del circuito representa un 43,06 % de la de fase. Las variantes de balanceo que se obtienen se encuentran ordenadas por corriente de neutro, pérdidas de energía y número de reconexiones respectivamente. Se debe escoger una variante, que, con un pequeño número de reconexiones, pueda lograr una buena reducción de la corriente de neutro y de las pérdidas de energía, atendiendo a esto se seleccionó la variante # 54. En la figura 3.10 se muestran las variantes de balanceo obtenidas para este circuito y los resultados de la variante seleccionada. Figura 3.10: Variantes de balanceo para el circuito 131. Al analizar los resultados obtenidos para la variante seleccionada, puede verse que la corriente de desbalance debe experimentar una disminución de un 39,43 %, las pérdidas de energía deben reducirse en 7,53 kWh/día y solo deben efectuarse cuatro cambios. Aunque los ahorros de energía que se consiguen en este caso son pequeños, debe notarse que las pérdidas de energía del circuito también lo son. Además, se puede lograr una buena reducción de la corriente de neutro con un reducido número de reconexiones en el circuito. CAPÍTULO III: RESULTADOS DEL BALANCEO DE LOS CIRCUITOS SELECCIONADOS 39 En la figura 3.11 se muestran los puntos del circuito donde deben efectuarse dichos cambios. Figura 3.11: Cambios propuestos para el circuito 131. En la figura 3.12 se observan los gráficos estimados por el programa de balanceo para las corrientes de fase y de neutro, antes y después de que se ejecuten las reconexiones. Figura 3.12: Corrientes estimadas antes y después de efectuar las reconexiones. En la figura anterior se puede ver que, al efectuar los cambios propuestos por el programa para este circuito, debe mejorar el equilibrio entre las corrientes de fase y debe disminuir significativamente la corriente de neutro. Por tanto, puede decirse que al aplicar los cambios propuestos por el programa este circuito debe operar de forma más eficiente y con menores pérdidas de energía. CONCLUSIONES 40 CONCLUSIONES A partir de los resultados obtenidos en el presente trabajo, puede llegarse a las siguientes conclusiones:  Mediante el celaje realizado a los circuitos en estudio se ha detectado que los datos de los circuitos se encontraban desactualizados.  Las bases de datos del Radial para todos los circuitos se han ajustado de forma bastante precisa a los gráficos de carga de las corrientes medidas por los recerradores NULEC.  En los circuitos analizados solo se tuvo acceso a mediciones de los interruptores principales, por lo tanto, fue necesario realizar una estimación de las cargas utilizando un programa de estimación.  El programa de balanceo basado en el algoritmo NSGA-II ha tenido un comportamiento favorable en todas las corridas en cuanto a tiempo de ejecución y convergencia. Brindando múltiples soluciones para balancear la red, lo que permite escoger la variante más adecuada consultando siempre con la dirección técnica de la Empresa Eléctrica.  Las soluciones obtenidas en los diferentes circuitos permiten estimar un buen comportamiento en cuanto a balance y pérdidas de los mismos cuando estas se logren aplicar. RECOMENDACIONES 41 RECOMENDACIONES  Aplicar los cambios propuestos por el programa de balanceo para los circuitos 32 y 131.  Continuar perfeccionando los métodos computacionales de estimación y ajuste de cargas en los circuitos de distribución primaria.  Realizar tomas de cargas en los transformadores para apoyar la estimación correcta de estas cargas.  Extender el estudio a otros circuitos de la provincia en aras de poder balancear la mayor parte de los circuitos primarios, y de esta forma lograr un servicio más eficiente.  Utilizar al máximo las potencialidades del programa de balanceo y realizar el balance de los circuitos en diferentes puntos cuando sea pertinente. REFERENCIAS BIBLIOGRÁFICAS 42 REFERENCIAS BIBLIOGRÁFICAS [1] I. P. Abril, "Genetic Algorithm for the Load Balance on Primary Distribution Cir cuits," IEEE Latin América Transactions, vol. 8, September 2010. [2] S. R. Castaño. Redes de distribución de Energía. Universidad Nacional de Colombia [online]. [3] S. R. Castaño. (2004, Redes de distribución de energía. Universidad Nacional de Colombia [online]. [4] E. P. D. Handbook. (2004). [5] N. C. N. O. N. d. Normalización. (2004. [6] I. 61000-4-30. (2003, Electromagnetic compatibility (EMC) - Part 4-30: Testing and measurement techniques - Power quality measurement methods. IEC 61000-4-30. [7] A. S. MG1. (1993, Motors and Generators. [8] T. P. f. P. I. M. a. G. I. Standard, IEEE Standard 112. (1991. [9] E. CENELEC, ed. [10] I.-.-. I. E. Commission, " in Compatibility levels for low-frequency conducted disturbances and signalling in public low-voltaje power supply systems, ed.>>, . (1990. [11] I.-.-. I. E. Commission, " in Part 2-2 Compatibility levels for low-frequency conducted disturbances and signalling in public low-voltaje power supply systems, ed.>>,, 2003. [12] O. N. d. N. O. N. de N. (NC), "Norma Cubana NC 365: 2011," in Tensiones Normalizadas, ed:>>, 2011. [13] D.-E. d. a. p. s. e. G. Casaravilla y V. Echinope. [14] M. p. e. B. d. C. d. D. P. M. A. Cuellar, Universiidad Central <> de Las Villas, Santa Clara. (2014. [15] D. d. T. e. S. T. ECAMEC Tecnología. (2009. [16] E. o. u. v. o. t. o. p. o. a. t.-p. i. m. C. Y. Lee. (1999, vol. 14, pp. 202-208. [17] I. P. Abril. (2012, Optimización Multiobjetivo del Balanceo de Circuitos de Distribución Primaria en in 16 Convención Científica de Ingeniería y Arquitura. [18] I. J. V. P. J. C. G. T. M. I. M. Piumetto. (2012, Análisis de la Revolución de los Factores de Desbalance en Sistema de Distribución de MT por la Inserción de la Generación Distribuida. en II Congreso Latinoamericano de Distribución Eléctrica. REFERENCIAS BIBLIOGRÁFICAS 43 [19] R. C. V. Sanabria. (1993, Análise Interactivo Gráfico de Sistemas Eléctricos de Distribución Primaria. Para completar las exigencias del grado de Magister en Ciencias de la Ingeniería. [20] A. I. P, "Fundamentación Técnica del Programa de Balance de Fases," Universidad Central "Marta Abreu" de Las Villas, Santa Clara, Mayo de 2016 2016. [21] L. D. A. Auguglio, and R. Romero. (2003, Minimum Losses reconfiguración of MV distribution networks through local control of tie-switches. IEEE Transactions on Power delivery 18, 762-771. [22] R. G. H. Salazar, and R. Romero. (2006, Artificial neural networks and clustering technique applied in the reconfiguration of distribution systems. IEEE Transactions on Power delivery vol. 21, pp. 1735-1742. [23] D. Das. (2006, Reconfiguration of distribution system using fuzzy multi- objective approach. Electric Power and Energy Systems vol. 28, pp. 331-338. [24] C.-F. C. J. -P. Chiou, and C. -T. Su. (2005, Variable scaling hybrid differential evolution for solving network reconfiguration of distribution systems IEEE Latin América Transactions on Power Systems vol. 20, pp. 668-674. [25] C.-F. C. C. -T. Su, and J. -P Chiou. (2005, Distribution network reconfiguration for loss reduction by ant colony search algorithm. EPRI vol. 75, pp. 190-199. [26] M. C. Z. Jinxiang, and Z. Chang. (1999, Phase Balancing using Mixed Integer Linear Programming. IEEE Transactions on Power Systems 13, 1487-1492. [27] B. G. Z. Jinxiang, and MY. Chow. (1999, "Phase Balancing using Simulated Annealing". IEEE Transactions on Power Systems 14, 1508-1513. [28] J. C. G. A. Ruiz, and R. Gallego. (2006, Solution to the phase balancing problem and the primary feeders reconfiguration by three phase modeling using simulated annealing. Scientia et Technica vol. 12, pp. 1-6. [29] T.-H. C. a. J.-T. Cherng. (2000, “Optimal phase arrangement of distribution transformers connected to a primary feeder for system unbalance improvement and loss reduction using a genetic algorithm”. IEEE Transactions on Power Systems 15. [30] R. P. B. M. Dilek, J. C. Thompson, and R. Sequin. (2001, “Simultaneous phase balancing at substations and switches with time-varying load patterns,” IEEE Transactions on Power Systems 16, 922-928. [31] C.-S. C. C.-H. Lin, H.-J. Chuang, and C.-Y. Ho. (2005, “Heuristic Rule-Based Phase Balancing of Distribution Systems by Considering Customer Load Patterns,”. IEEE Transactions on Power Systems 20, 709-716. [32] C.-S. C. C.-H. Lin, M.-Y. Huang, H.-J. Chuang, M. S. Kang, C.-Y. Ho, and C.W. Huang. (2007, “Optimal Phase Arrangement of Distribution Feeders Using Immune Algorithm,”. 14th International Conference on Intelligent System Applications to Power Systems, ISAP 2007, Kaohsiung, Taiwan, November 4 - 8. REFERENCIAS BIBLIOGRÁFICAS 44 [33] C.-S. C. M.-Y. Huang, C.-H. Lin, M. S. Kang, H.-J. Chuang, and C.-W. Huang. (2008, “Three-Phase Balancing of Distribution Feeders using Immune Algorithm,”. IET Generation, Transmission and Distribution 2, 383-392. [34] C.-S. C. C.-H. Lin, H. J. Chuang, M.-Y. Huang, and C.-W. Huang. (2008, “An Expert System for Three-Phase Balancing of Distribution Feeders,”. IEEE Transactions on Power Systems 23, 1488-1496. [35] M. W. S. A. Ukil, and J. Jordaan. (2006, "Feeder Load Balancing using Neural Network,” Lecture Notes in Computer Science, Springer 39, 1311-1316. [36] A. U. a. M. W. Siti. (2008, “Feeder Load Balancing using Fuzzy Logic and Combinatorial Optimization-based Implementation,”. Electric Power Systems Research 78(11), 1922-1932. [37] A. S. N. Gupta, and K. R. Niazi. (2011, A novel strategy for phase balancing in three-phase four-wire distribution systems. IEEE Power and Energy Society General Meeting, San Diego (CA), 1-7. [38] R. A. H. a. S. H. Soltani. (2012, Fuzzy Optimal Phase Balancing of Radial and Meshed Distribution Networks Using BF-PSO Algorithm. IEEE Transactions on Power Systems 27, 47-57. [39] a. S. H. S. R. A. Hooshmand. (2012, Simultaneous optimization of phase balancing and reconfiguration in distribution networks using BF–NM algorithm. International Journal of Electrical Power & Energy Systems 41, 76- -86. [40] F. C. Guo. (2012, Solving Load Phase Balancing Problem in LV Distribution Networks by Chaotic Simulated Annealing. Advanced Materials Research 463, 689-693. [41] L. L. M. Sathiskumar, and S. Thiruvenkadam. (2012, A self adaptive hybrid differential evolution algorithm for phase balancing of unbalanced distribution system. International Journal of Electrical Power & Energy Systems 42, 91- 97. [42] M. S. G. Mahendran, S. Thiruvenkadam, and L. Lakshminarasimman. (2013, Multi-objective Unbalanced Distribution Network Reconfiguration through Hybrid Heuristic Algorithm. Journal of Electrical Engineering & Technology 8, 215-222. [43] S. S. K. Wang, and T. G. Robertazzi. (2013, Phase balancing algorithms. Electric Power Systems Research 96, 218-224. [44] A. S. N. Gupta, and K. R. Niazi. (2014, A novel method for simultaneous phase balancing and mitigation of neutral current harmonics in secondary distribution systems. International Journal of Electrical Power & Energy Systems 55, 645-656. [45] I. P. Abril. (2016, Fundamentación Técnica del Programa de Estimación de Cargas. Universidad Central “Marta Abreu” de Las Villas, Santa Clara. REFERENCIAS BIBLIOGRÁFICAS 45 [46] I. P. Abril, "NSGA-II phase balancing of primary distribution circuits by the reconnection of their circuit laterals and distribution transformers," Electric Power Systems Research, vol. 109, pp. 1-7, 2014. [47] I. P. Abril, "Fundamentación Técnica del Programa de Balance de Fases," Universidad Central “Marta Abreu” de Las Villas, Santa Clara, mayo de 2016. [48] A. Seshadri. NSGA-II source code. available in: «http://www.mathworks.com /matlabcentral/fileexchange/10429-nsga-ii-a-multi-objective-optimization- algorithm/content/NSGA-II/». http://www.mathworks.com/ ANEXOS 46 ANEXOS Anexo I: Planos atualizados de los circuitos ANEXOS 47 Anexo II: Base de datos para el circuito 32 ANEXOS 48 Anexo III: Base de datos para el circuito 131