373185c TesisEfraínPaucar
373185c TesisEfraínPaucar
373185c TesisEfraínPaucar
TESIS
ELABORADO POR
ASESOR
LIMA - PERÚ
2014
DEDICATORIA
Dedico esta tesis a: Mi madre María por su verdadero amor, a la memoria de mi padre
Gonzalo, a mi esposa Gladys y a mis hijos Aldrin y Greiz por su comprensión y constante
apoyo.
ii
AGRADECIMIENTOS
iii
ÍNDICE DE CONTENIDOS
INTRODUCCIÓN ............................................................................................................... 1
CAPITULO I
PLANTEAMIENTO DE INGENIERÍA DEL PROBLEMA ..................................................... 2
1.1 Objetivos ...................................................................................................................... 2
1.1.1 Objetivo general ....................................................................................................... 2
1.1.2 Objetivos específicos ................................................................................................ 2
1.2 Delimitación del estudio ............................................................................................... 3
1.2.1 Alcances ................................................................................................................... 3
1.2.2 Limitaciones.............................................................................................................. 3
1.3 Descripción del problema............................................................................................. 3
1.4 Antecedentes de los trabajos realizados con anterioridad............................................ 6
1.5 Hipótesis básica ........................................................................................................... 7
1.6 Metodología utilizada para la solución de la problemática planteada 8
CAPITULO II
MARCO TEÓRICO ............................................................................................................ 9
2.1 Pérdidas de energía en los sistemas eléctricos de distribución .................................... 9
2.1.1 Cálculo de pérdidas de energía en conductores eléctricos de distribución.............. 11
2.1.2 Cálculo de pérdidas de energía en transformadores de distribución ....................... 13
2.1.3 Cálculo de pérdidas de energía en contadores de energía ..................................... 14
2.1.4 Cálculo de pérdidas de energía en aisladores ........................................................ 16
2.2 Flujo de potencia en sistemas de distribución radiales ............................................... 18
2.2.1 Método escalera ..................................................................................................... 18
2.2.2 Método de suma de potencias ................................................................................ 20
2.2.3 Método de suma de corrientes................................................................................ 23
2.2.4 Otros métodos de flujo de carga para distribución .................................................. 25
2.3 Métodos de cálculo de pérdidas técnicas y no técnicas de energía en sistemas de
distribución eléctrica ........................................................................................................ 30
2.3.1 Flujo de carga ......................................................................................................... 30
iv
2.3.2 Procedimiento de cálculo de pérdidas de potencia y energía en distribución de
ANEEL – Brasil ................................................................................................................ 33
2.3.3 Otros métodos aproximados utilizados para la estimación de pérdidas de potencia y
energía en Baja Tensión .................................................................................................. 42
2.4 Datos del sistema eléctrico de distribución................................................................. 44
2.4.1 Datos de las redes eléctricas de distribución .......................................................... 45
2.4.2 Datos de las cargas en el sistema de distribución eléctrica .................................... 47
2.4.3 Estudios de caracterización de carga ..................................................................... 47
2.4.4 Estudios de cargabilidad de conductores y transformadores MT/BT desde la óptica
de pérdidas mínimas........................................................................................................ 49
CAPITULO III
METODOLOGÍA PARA LA SOLUCIÓN DEL PROBLEMA .............................................. 56
3.1 Diseño de base de datos para estudio de pérdidas técnicas de energía .................... 58
3.1.1 Datos de redes eléctricas de media tensión............................................................ 58
3.1.2 Datos de transformadores MT/BT ........................................................................... 60
3.1.3 Datos de redes eléctricas de baja tensión .............................................................. 60
3.1.4 Datos de acometidas .............................................................................................. 61
3.2 Procedimiento para la estimación de pérdidas técnicas y no técnicas de energía
mediante el método propuesto......................................................................................... 61
3.2.1 Importación de información de GIS y consumos de clientes ................................... 63
3.2.2 Cálculo de demanda máxima de cargas BT a partir de los consumos de energía .. 64
3.2.3 Obtención de curva de carga de una subestación MT/BT para análisis de flujo de
carga en baja tensión....................................................................................................... 65
3.2.4 Cálculo de pérdidas técnicas en MT y BT ............................................................... 65
3.2.5 Cálculo y localización de pérdidas no técnicas en MT y BT .................................... 67
3.2.6 Análisis en mapas temáticos de las pérdidas de energía ........................................ 67
3.3 Viabilidad del desarrollo aplicativo ............................................................................. 72
CAPITULO IV
ANÁLISIS Y PRESENTACIÓN DE RESULTADOS ......................................................... 73
4.1 Areas de influencia y características técnicas de los elementos y cargas que
conforman la topología de la red ...................................................................................... 73
4.2 Resultados obtenidos ................................................................................................ 76
v
4.2.1 Pérdidas técnicas de energía obtenidas en media tensión ..................................... 76
4.2.2 Pérdidas técnicas de energía obtenidas en baja tensión ........................................ 77
4.2.3 Pérdidas técnicas de energía total en distribución MT y BT .................................... 78
4.2.4 Pérdidas no técnicas de energía en distribución ..................................................... 79
4.3 Análisis e interpretación de los resultados ................................................................. 81
CONCLUSIONES ............................................................................................................ 82
RECOMENDACIONES .................................................................................................... 84
BIBLIOGRAFÍA ................................................................................................................ 85
ANEXOS.......................................................................................................................... 88
vi
ÍNDICE DE TABLAS
Tabla 4.1: Magnitud de líneas de distribución de media tensión por alimentadores ......... 75
Tabla 4.2: Magnitud de líneas de distribución de baja tensión por alimentadores ............ 75
Tabla 4.3: Número de usuarios por alimentadores ........................................................... 76
Tabla 4.4: Porcentaje de pérdidas técnicas de energía en media tensión ........................ 76
Tabla 4.5: Porcentaje de pérdidas técnicas de energía en baja tensión ........................... 77
Tabla 4.6: Porcentaje de pérdidas técnicas totales en distribución .................................. 79
Tabla 4.7: Porcentaje de pérdidas totales, técnicas y no técnicas de energía en BT ....... 80
vii
ÍNDICE DE FIGURAS
ix
RESUMEN
El presente trabajo tiene por objetivo resolver el problema que tienen las empresas
eléctricas de distribución para definir sus políticas y estrategias de reducción y control de
pérdidas de energía, debido a la dificultad en conocer cuáles son los niveles de pérdidas
de energía no técnicas y técnicas en cada parte de los sistemas de distribución.
Para ello, se parte de la premisa que se conocen las pérdidas totales de energía a
través de las mediciones de los alimentadores de media y baja tensión y en los puntos de
consumo (clientes). Luego mediante métodos simplificados sin perder precisión se
obtienen las pérdidas técnicas. Por último, las pérdidas no técnicas se calculan de
manera indirecta, restando las pérdidas técnicas de las pérdidas totales.
Por ello, en esta tesis se plantea una metodología de bajo costo para el cálculo de
pérdidas técnicas en sistemas eléctricos de distribución. La misma se basa en el Flujo de
Potencia en Distribución, usa datos fáciles de obtener y toma en cuenta resultados de
estudios anteriores de caracterización de carga y pérdidas en conductores y
transformadores de distribución.
x
ABSTRACT
This paper aims to solve the problem with the electricity distribution companies to
define their policies and strategies to reduce and control energy losses due to the difficulty
in knowing what the levels of energy losses in technical and non-technical each part of the
distribution systems.
To do this, we start from the premise that we know the total energy losses through
measurements feeders and low voltage and consumption points (customers). Then using
simplified methods without losing accuracy obtained technical losses. Finally, non-
technical losses are calculated indirectly, by subtracting the technical losses of total
losses.
This task is complicated by the large amount of information to handle the uncertainty
of the data, the high costs associated with measurements, the high level of imbalance, the
complexity of the distribution network and the lack of standardized methods.
xi
INTRODUCCIÓN
Con los métodos tradicionales aplicados hasta ahora para determinar las pérdidas de
energía en las redes de distribución del sistema eléctrico peruano resulta bastante
oneroso y difícil para las empresas de distribución eléctrica identificar y localizar los
excesos de pérdidas de energía que superan las pérdidas técnicas. Esto conlleva a que
en los últimos años los reportes de pérdidas de energía de las empresas eléctricas de
distribución del país, muestren pérdidas económicas permanentes debido al exceso de
pérdidas de energía respecto de las pérdidas estándares.
La delimitación de este estudio viene dada por los alcances y las limitaciones que se
encuentran en la elaboración de la presente tesis y que se mencionan a continuación.
1.2.1 Al cances
- El presente trabajo se enfoca en la aplicación de la metodología de determinación de
pérdidas de energía para empresas de distribución eléctrica ubicadas en el territorio
nacional.
- La metodología será aplicada a un sistema de distribución eléctrica de
aproximadamente 70,000 usuarios.
1.2.2 Limitaciones
Debido a limitaciones de tiempo y recursos económicos, la implementación de la
metodología propuesta se hace a nivel académico, por lo que el desarrollo de un software
comercial no está contemplado en esta etapa de la investigación.
De otra parte, en los últimos años los reportes de pérdidas de energía de las
empresas eléctricas de distribución del país, muestran pérdidas económicas
permanentes debido al exceso de pérdidas de energía respecto de las pérdidas
estándares. Esta tendencia se inicia desde que en el año 2005 el Organismo Regulador
del Mercado Eléctrico de Distribución (OSINERGMIN – GART), retira los incentivos
económicos vía tarifa que tenían las empresas de distribución. Dichos incentivos
consistían en la incorporación en el VAD, desde el 9.0% de excesos de pérdidas de
energía respecto a las pérdidas estándares en el año 1993 hasta llegar al 0% para el año
2005.
3
En nuestro país en estos últimos diez años se ha incrementado la demanda con una
tasa de crecimiento jamás prevista y la gestión eficiente de la distribución eléctrica es
muy importante para disponer de recursos para atender las inversiones necesarias para
atender la demanda actual y futura.
Esta importante reducción en las pérdidas se produjo porque las mismas tenían un
alto componente de pérdidas no técnicas (conexiones clandestinas, alteraciones a las
mediciones del medidor de energía, y otras formas de hurto y fraude) fáciles de localizar
e identificar; asimismo, se detectaron las pérdidas técnicas de energía más evidentes
como transformadores y conductores sobrecargados, cargas desbalanceadas, etc.
TOTAL PERÚ
(1993-2011)
25%
Excesos de Pérdidas incluido en VAD
21.9%
Pérdidas Estándares
20.6%
19.7% Pérdidas Reales
20%
17.1%
14.6%
15%
12.4%
11.5%
9.0%
8.7%
8.5%
8.4%
10.3%
7.7%
9.7%
9.1% 9.0% 8.8%
10%
4.8%
8.6% 8.6%
4.3%
8.2% 8.0%
3.6%
3.3%
7.8% 7.8%
2.5%
7.7%
2.0%
1.5%
0.8%
0%
0%
0%
0%
0%
0%
5%
7.91%
7.59%
7.45%
7.39%
7.33%
7.28%
7.25%
7.12%
7.10%
7.05%
7.04%
6.95%
6.88%
6.85%
6.83%
6.84%
6.82%
6.81%
6.74%
0%
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
4
Este comportamiento se explica, porque al encontrarse en estos niveles de
eficiencia, es difícil y costoso con los métodos tradicionales aplicados hasta ahora,
identificar y localizar los excesos de pérdidas de energía por encima de las pérdidas
técnicas.
Con relación a las pérdidas estándares, en los años 2003, 2005 y 2011 estas fueron
de 7.10%, 6.95% y 6.74%, respectivamente. El reto para las empresas eléctricas es
alcanzar niveles de eficiencia con pérdidas de energía menor al estándar de 6.74%.
Por otro lado, las herramientas computacionales utilizadas por las empresas para
realizar estudios de esta naturaleza se basan en programas de flujos de potencia
trifásicos, los que no son adecuados para modelar el gran volumen de datos necesario
para el análisis de pérdidas de energía en redes de distribución. Esto, debido a que estas
aplicaciones solo proporcionan las pérdidas totales de potencia en las redes, no genera
reportes referidos a las pérdidas de potencia por cargas, no modelan redes
desbalanceadas ni tienen información de la ubicación física de la carga. Asimismo,
debido a que el modelamiento de las cargas se hace de forma tradicional (muestreo por
un grupo de alimentadores de media tensión y circuitos de baja tensión con los que se
obtiene indicadores promedio de la empresa o sistema analizado), no identifican donde
se encuentran los excesos de pérdidas de energía por lo que no ayudan a tomar
decisiones.
5
En resumen, se puede decir que con respecto a las pérdidas de energía, existe la
necesidad por parte de las empresas de distribución eléctrica de estimar en forma
detallada, en cada componente o parte los sistemas eléctricos de distribución, sus
pérdidas técnicas y no técnicas.
6
distribución. Debido a esto el proceso de determinación de pérdidas era un proceso que
demandaba varios meses de trabajo.
De esta forma lo que en esta tesis se probará será la siguiente hipótesis básica: Es
posible determinar con precisión la ubicación de las pérdidas no técnicas de energía en
las redes de distribución del sistema eléctrico peruano, a pesar de lo voluminosa que
resulta la información y por otra parte, a veces escasez de algunos datos de los sistemas
de distribución.
7
1.6 Metodología utilizada para la sol ución de la problemática
planteada
Figura 1.2 Flujograma de la metodología propuesta en esta tésis para determinar las
pérdidas de energía en sistemas de distribución
(Fuente: Elaboración propia)
8
C API TULO II
M ARCO TEÓRICO
Clasificación de
Pérdidas de Energía
-Líneas
-Transformadores
Administrativas Accidentales Fraudulentas
-Aisladores
-Medidores
Las pérdidas no técnicas se definen como la diferencia entre las pérdidas totales de
un sistema eléctrico y las pérdidas técnicas estimadas para el mismo [1]. Se pueden
clasificar en tres tipos:
a) Accidentales, las cuales tienen su origen en el mal uso u operación de los elementos y
equipos de los circuitos eléctricos, tal es el caso de errores en el conexionado.
b) Administrativas, energía que por algún motivo no se contabiliza: usuarios sin
medidores (toma directa), ferias, etcétera.
10
c) Fraudulentas, referidas a la energía que toman algunos consumidores evitando
mediante algún mecanismo pasar por los contadores de energía de la compañía de
electricidad.
Para fines del presente estudio las pérdidas técnicas se clasifican en relación a la
variabilidad de la carga (o corriente eléctrica) en dos tipos [1]:
Para calcular las pérdidas en los conductores eléctricos de distribución una línea de
distribución se puede representar de manera simplificada como se muestra en la Figura
2.2. El exceso de pérdidas de energía al utilizar cables especificados sin considerar la
minimización de los costos totales del sistema (costos de inversión y de operación a lo
11
largo de la vida útil de la instalación) se traduce en mayores costos para las empresas y
los usuarios.
i j
R + jX
Vi Vj
Donde:
Vi : Tensión en el punto i (V).
Vj : Tensión en el punto j (V).
R : Resistencia de la línea (Ω).
X : Reactancia de la línea (Ω).
Donde:
Plínea max : Pérdidas de potencia en el conductor en condiciones de máxima demanda
T : Tiempo en horas
Llínea : Pérdidas de energía en el conductor en el periodo T
12
2.1.2 Cálculo de pérdi das de energía en transformadores de
distribución
ii ij
i j
Rs + jXs
Vi ip Rm Xm im Vj
Donde:
13
im : Corriente magnetizante (A).
Vi 2
Pfe , ( 2.3)
Rm
Pcu Rs I 2 , ( 2.4)
Donde:
Donde:
T : Tiempo en horas
Las pérdidas en el circuito de tensión se conocen como pérdidas en vacío o fijas (son
proporcionales al cuadrado de la tensión) y las pérdidas en el circuito de corriente como
14
las pérdidas en carga (proporcionales al cuadrado de la corriente que transporta el
contador de energía).
Los valores de los consumos internos en vacío y con carga de los contadores de
energía ya sean electromecánicos o electrónicos se obtienen de forma práctica al realizar
las mediciones en vacío y con carga. Asimismo, son datos proporcionados por los
fabricantes.
ii ij
i j
Rs + jXs
Vi ip Rm Xm im Vj
Donde:
15
Vi 2
Pvm , (2.6)
Rm
Pcm Rs I 2 , (2.7)
Donde:
Donde:
T : Tiempo en horas
16
Aunque no pueden darse valores estándares de pérdidas, es frecuente adoptar para
las mismas los siguientes valores de corriente permanente en ambiente seco y otro valor
para ambiente con humedad>85%:
- Corriente de fuga de (5 A) a (50 A) por aislador de silicona.
- Corriente de fuga de (400 A) a 1.20 mA por aislador de porcelana.
- Por encima de 1.20 mA se superponen a la onda de corriente de fuga, impulsos de
muy corta duración modificando la forma de onda con predominancia en un valor
medio de 500 mA [18].
La pérdida de potencia estimada oscila entre los 0.5 W y los 10 W por aislador, las
mayores pérdidas ocurren donde hay mayor contaminación y humedad.
Es importante tener en cuenta la influencia de las estaciones (primavera, verano,
otoño, invierno) en la humedad del ambiente, ya que con una humedad relativa entre 85 y
100% se pueden triplicar los promedios de corrientes de fuga respecto a las de estación
seca.
Las normas aplicables para los aisladores son: IEC 61109, IEC 91952, ANSI C29.1,
C29.2, C29.5, C29.6, C29.7, C29.11 y C29.13.
La fórmula para determinar el efecto aislador por km y fase será la dada para la
conductancia en general.
Las pérdidas de potencia totales simplemente serán las pérdidas producidas por un
aislador, multiplicadas por el número de aisladores totales:
PTaisl Paisl Naisl , (2.9)
Donde:
Este cálculo se realiza para los aisladores de media tensión, en baja tensión las
corrientes de fuga en los aisladores son prácticamente iguales a cero.
17
2.2 Flujo de potencia en sistemas de distri bución radiales
Aquí se describen algunos métodos de flujo de carga utilizados para sistemas
eléctricos de distribución. Los cuales son típicamente radiales, desbalanceados, con una
alta relación R/X y un gran número de líneas y cargas balanceadas y desbalanceadas,
por lo que no se puede utilizar el análisis de flujo de potencia con modelos matriciales
aplicados a los sistemas eléctricos de potencia.
Entre los métodos de flujo de potencia más conocidos para sistemas de distribución
se tienen:
- Método Escalera.
- Método de suma de potencias.
- Método de suma de corrientes (seleccionado en la presente tesis para desarrollar el
método práctico para el análisis de pérdidas de energía en los sistemas eléctricos de
distribución).
- Otros métodos de flujo de carga para distribución.
Para las ilustraciones de los métodos más conocidos, en esta tesis se utiliza el
siguiente diagrama simplificado de un sistema de distribución.
Vs I1 I2
Pc1+jQc1 Pc2+jQc2
Figura 2.5 Diagrama unifilar simplificado de un sistema de distribución
(Fuente: Tesis de R. Vilcahuamán, PUCC-Chile 1993[3])
Este método resuelve la red aguas arriba, asigna valores iniciales de perfil de tensión
y aplica directamente las leyes de Kircchoff, hasta llegar al nodo fuente. La diferencia de
tensiones (error) que se obtenga entre la tensión de este nodo y el calculado en el
proceso, se suma al perfil de tensión previamente supuesto a fin de obtener un nuevo
perfil de tensiones para la próxima iteración. El proceso termina cuando la tensión
calculada en el nodo fuente es el especificado o la diferencia entre ellos es pequeña. Este
método fue presentado inicialmente por Kersting y otros [3].
Se compara las tensiones de la iteración actual de los nodos contra los valores de la
iteración anterior:
Donde:
: Es la tolerancia de error aceptable para que converjan las ecuaciones del
sistema (ejemplo 0.001).
VS : Módulo de la tensión conocida en la fuente (valor especificado).
VS ( k 1) : Módulo de la tensión calculada en la iteración “k+1” en el nodo fuente.
I 1( k ) , I 2( k ) : Corrientes en los tramos 1 y 2 correspondiente a la iteración “k”.
Z 1 , Z 2 : Impedancias de los tramos 1 y 2.
19
2.2.2 Método de suma de potencias
Dado un perfil de tensiones iniciales asignados, se determinan las potencias activas
y reactivas en los consumos (vista desde los terminales de los transformadores de
distribución o punto de conexión de la carga), además de las potencias reactivas
inyectadas por los condensadores. En un proceso simultáneo aguas arriba de la red, se
calculan las potencias equivalentes en las barras y las pérdidas de potencia en las líneas
[4].
En cada uno de los nodos, inicialmente se asumen las tensiones equivalentes o muy
cercanas a 1.0 en p.u. Se calcula las potencias en el nodo 1 y en el nodo Fuente:
2
S c2
S 1( k ) ( Pc1 jQc1 ) Z 2 ( Pc 2 jQc 2 ) (2.14)
V 2( k )
2
S 1( k )
S s(k ) S 1( k ) Z 1
V 1( k )
(2.15)
Se calcula las tensiones de nodos aguas abajo:
*
S
V 1( k 1)
V s Z1
1( k )
V (2.16)
1( k 1)
* *
V1( k 1) 2
V 1( k 1) Vs Z 1 S 1( k ) (2.17)
20
Al desarrollar la ecuación anterior se tiene:
V1( k 1) 2 V1( k 1)Cos jV1( k 1) Sen Vs ( R1 P1( k ) X1Q1( k ) ) j ( X1 P1( k ) R1Q1( k ) ) (2.18)
V1( k 1) 4 AV
1 1( k 1) B1 0
2
(2.21)
Donde:
La ecuación con estos coeficientes no depende del ángulo de desfase (β) entre las
2
A A
V1( k 1) 1 1 B1
2
(2.24)
2
Si se considera que en condiciones normales de operación, las caídas de tensión y las
pérdidas del sistema deben ser pequeñas y que los módulos de las tensiones son
cercanos a la unidad, se tiene que A1 será siempre negativo y B1 positivo y de valor
2
A1
reducido. Por lo tanto, el término B1 siempre será pequeño y cercano a
2
2
A1
. De esta forma la única solución que tiene sentido físico es de la forma:
2
21
2
A A
V1( k 1) 1 1 B1 (2.25)
2 2
Una vez que se tienen las magnitudes de las tensiones, se pueden calcular los
ángulos de éstos a partir de la ecuación:
( X 1 P1( k ) R1Q1( k ) )
Cos1( k 1) tg 1 2
(R P X Q ) V (2.26)
1 1( k ) 1 1( k ) 1( k 1)
Tensiones en el nodo 2:
*
S
V 1( k 1) Z 2
c2
V 2( k 1)
V (2.27)
2( k 1)
Para obtener las tensiones en el Nodo 2, se procede de la misma forma que para el
nodo 1.
Vi ( k 1) Vi ( k )
(2.28)
Ps ( k 1) Ps ( k ) (2.29)
Qs ( k 1) Qs ( k ) (2.30)
Donde:
22
VS ( k 1) : Módulo de la tensión calculada en la iteración “k+1” en el nodo fuente.
V1( k 1) ,V2( k 1) : Módulo de las tensiones calculadas en los nodos 1 y 2, iteración “k+1”.
I 1( k ) , I 2( k ) : Corrientes en los tramos 1 y 2 correspondiente a la iteración “k”.
Las variables conocidas son la tensión de la fuente (Vs = 1.0) y las cargas (Pi+jQi), las
variables a calcular son las corrientes (Ii) y las tensiones (Vi) en los nodos de carga o
ramales.
En cada uno de los nodos, inicialmente se asumen las tensiones equivalentes o muy
cercanas a 1.0 en p.u. Seguidamente, se calculan las corrientes en los tramos 1 y 2:
*
S2
I 2( k )
V (2.31)
2( k )
*
S1
I 1( k ) I 2( k )
V (2.32)
1( k )
Se calcula las tensiones de nodos aguas abajo:
23
V 1( k 1) V s Z 1 I 1( k ) (2.33)
V 2( k 1) V 1( k 1) Z 2 I 2( k ) (2.34)
Se comparan las tensiones de los nodos de la iteración actual contra los valores de
Vi ( k 1) Vi ( k )
la iteración anterior. Si se cumple que se culmina el cálculo de
flujo de carga y se reportan los valores obtenidos en la última iteración.
V1( k 1) ,V2( k 1) : Módulo de las tensiones calculadas en los nodos 1 y 2, iteración “k+1”.
I 1( k ) , I 2( k ) : Corrientes en los tramos 1 y 2 correspondiente a la iteración “k”.
Pc1 , Qc1 , Pc 2 , Qc 2 : Potencia activa y reactiva de las cargas en los nodos 1 y 2.
Z 1 , Z 2 : Impedancias de los tramos 1 y 2.
Se debe tener presente que para el flujo de carga trifásico para sistemas
desbalanceados con este método u otro similar se parte con tensiones especificadas por
fase en la fuente cuyos valores en p.u. son: 1.00 1.0 120 1.0120 y las
24
tensiones iniciales de los nodos aguas abajo también presentan estos mismos valores o
muy cercanos.
V a ( n 1) V an Z aa Z ab Z ac Z ad I a ( n 1)
V b ( n 1) V bn Z ba Z bb Z bc Z bd I b ( n 1)
(2.35)
V c ( n 1) V cn Z ca Z cb Z cc Z cd I c ( n 1)
V d ( n 1) V dn Z da Z db Z dc Z dd I d ( n 1)
Aquí se representa la influencia de las impedancias mutuas entre las distintas fases y
la corriente de desbalance I d ( n1) I a ( n1) I b ( n1) I c ( n1) que hay en el conductor
neutro o retorno por tierra.
Las cargas en los nodos vienen a ser la suma de las cargas de las laterales
(ramales) incluida las cargas que se conectan en el nodo.
25
i
V t 1 2 n
s n-1
…
Z Z Z
1 2 n
i i i i
1 2 n
(n-1)
Siendo la tensión en el nodo “n” respecto al neutro igual a Vn y en el nodo “n-1” igual
a Vn-1, se tiene la siguiente ecuación que representa la tensión en el nodo “n”:
(2.36)
n n 1 n
Vn Vs it Z j i j ( Zi ) (2.40)
j 1 j 1 i j 1
26
Para el caso de redes de distribución con cargas desbalanceadas y disposición
asimétrica de los conductores el modelo de red a utilizar se muestra en la Figura 2.7.
n n 1 n
itb Z abk ibk ( Z abi )
k 1 k 1 i k 1
n n 1 n
itc Z ack ick ( Z aci )
k 1 k 1 i k 1
n n 1 n
Vbn Vbs itb Zbbk ibk ( Zbbi ) (2.42)
k 1 k 1 i k 1
n n 1 n
ita Z abk iak ( Z abi )
k 1 k 1 i k 1
n n 1 n
itc Zbck ick ( Zbci )
k 1 k 1 i k 1
n n 1 n
Vcn Vcs itc Z cck ick ( Z cci ) (2.43)
k 1 k 1 i k 1
n n 1 n
ita Z ack iak ( Z aci )
k 1 k 1 i k 1
n n 1 n
itb Zbck ibk ( Zbci )
k 1 k 1 i k 1
Vi ( k 1) Vi ( k )
El proceso de cálculo culmina cuando se cumple que
27
Figura 2.7 Red equivalente de sistemas trifásicos desbalanceados
(Fuente: Masoud Aliakbar GOLKAR y otros, Tehrán-IRAN, 2007[16])
28
Tomando como referencia la red de distribución mostrada en la Figura 2.6 y
reemplazando el modelo de cada tramo de red por el equivalente mostrado en la Figura
2.7; vemos que la tensión de la fase “c” en el nodo “n” es igual a la tensión Vcs menos la
caída de tensión que se produce entre el nodo “s” y “n”. Aplicando ley de Kirchhoff
calculamos la tensión en “n”:
Vcn Vcs itc Z cc1 ita Z ac1 itb Z bc1 (itc ic1 )Z cc 2 (ita ia1 )Z ac 2 (itb ib1 )Z bc 2
(itc ic1 ic 2 ) Z cc 3 (ita ia1 ia 2 )Z ac 3 (itb ib1 ib 2 )Z bc 3
(itc ic1 ic 2 ... icn 1 ) Z ccn (ita ia1 ia 2 ... ian1 ) Z acn
(2.44)
(itb ib1 ib 2 ... ibn1 ) Zbcn
Vcn Vcs itc Z cc1 ita Z ac1 itb Z bc1 itc Z cc 2 ic1Z cc 2 ita Z ac 2 ia1Z ac 2
itb Zbc 2 ib1Zbc 2 itc Z cc 3 (ic1 ic 2 )Z cc 3 ita Z ac3 (ia1 ia 2 )Z ac3
itb Zbc 3 (ib1 ib 2 ) Z bc 3 itc Z ccn (ic1 ic 2 ... icn 1 ) Z ccn
ita Z acn (ia1 ia 2 ... ian 1 ) Z acn itb Z bcn (ib1 ib 2 ... ibn 1 ) Z bcn (2.45)
Reordenando los términos de la ecuación anterior de forma que se agrupen las caídas
de tensión debido a impedancias propias y aparte las que se deben a las impedancias
mutua se tiene:
Vcn Vcs itc Zcc1 itc Zcc 2 itc Zcc3 itc Zccn ic1Zcc 2 (ic1 ic 2 )Z cc3 (ic1 ic 2 ... icn1 )Z ccn
ita Z ac1 ita Z ac 2 ita Z ac3 ita Z acn ia1Z ac 2 (ia1 ia 2 )Z ac3 (ia1 ia 2 ... ian1 )Z acn
itb Zbc1 itb Zbc 2 itb Zbc3 itb Zbcn ib1Zbc 2 (ib1 ib 2 )Zbc3 (ib1 ib 2 ... ibn1 )Zbcn
(2.46)
29
n n 1 n
Vcn Vcs itc Z cck ick ( Z cci ) (2.47)
k 1 k 1 i k 1
n n 1 n
ita Z ack iak ( Z aci )
k 1 k 1 i k 1
n n 1 n
itb Zbck ibk ( Zbci )
k 1 k 1 i k 1
Donde:
Vcn : Tensión de fase “c” con referencia al neutro en el nodo “n”
Vcs : Tensión de fase “c” con referencia al neutro en el nodo “s” (fuente)
: Corriente total en cada fase que circula en el tramo de línea inicial
: Corriente de las fases “a”, “b” y “c” del ramal que parte del nodo “1”
: Corriente de las fases “a”, “b” y “c” del ramal que parte del nodo “n”
: Impedancia propia de cada fase de un tramo de línea
: Impedancia mutua de un tramo de línea
: Impedancia propia por fase del tramo de línea que llega al nodo “n”
: Impedancia mutua del tramo de línea que llega al nodo “n”
30
El éxito de este método depende principalmente de la cantidad y calidad de
información de la infraestructura eléctrica de distribución y de las cargas (cuya
información es escasa), así como del modelo de flujo de carga utilizado [1].
El valor del factor de pérdidas depende de cada sistema en particular. Una primera
aproximación del valor se puede estimar a partir del factor de carga. Se demuestra que:
f c FL f c
2
(2.49)
Donde:
fc
: Factor de carga del sistema.
31
Figura 2.8 Diagrama de carga horario de un periodo de un día
(Fuente: Elaboración propia)
Una relación empírica entre el factor de pérdidas y el factor de carga desarrollado por
Bullery Woodrow es:
FL X fc (1 X ) fc 2 (2.50)
32
coincidentes con la máxima demanda) hasta cumplir con el requisito de balance de
potencias.
Datos de
Redes
Eléctrica y Proceso de cálculo Estimación de
Inicio ¿Ps = Pc+L? SI
Demandas de de flujo de carga Pérdidas de Energía
Cargas y
Fuente
NO
Modificar Pc
Fin
(Pc*Fs)
33
Figura 2.10 Diagrama de Flujo energía requerida y transferida
(Fuente: ANEEL-Brasil, 2007[10])
34
Transformadores A3/A3a e A2/A3a;
Redes A3;
Transformadores A2/A3;
Redes A2;
Transformadores A1/A2;
Redes A1.
fp %MerccUC
nUCs
fp c
UC
(2.51)
c1
Donde:
fp: factor de potencia típico
nUCs: Número de clientes o puntos de consumos medidos
%MerccUC: % de participación del mercado por este tipo de consumos
fpcUC: factor de potencia del consumo medido UC
El factor de carga típico será calculado con la expresión siguiente:
fc %MerccUC
nUCs
fc c
UC
(2.52)
c1
Donde:
fc: factor de carga típico
nUCs: Número de clientes o puntos de consumos medidos
%MerccUC: % de participación del mercado por este tipo de consumos
fccUC: factor de carga del consumo medido UC
El factor de pérdidas típico será calculado con la expresión siguiente:
fpe %MerccUC
nUCs
fpe c
UC
(2.53)
c1
Donde:
fpe: factor de pérdidas típico
nUCs: Número de clientes o puntos de consumos medidos
%MerccUC: % de participación del mercado por este tipo de consumos
36
fpecUC: factor de pérdidas del consumo medido UC
k: Constante.
Nt
d Snom
i i
dceq i 1
Nt (2.55)
Snom
i 1
i
Donde:
dceq: Distancia de carga equivalente
Nt: Número total de transformadores conectados a la red
di: Distancia de la fuente (origen de red) hasta el transformador “i”
Snomi: Potencia nominal del transformador “i”
P max
2 2
ltot Vb Cos b
2
p (2.56)
mp V Cos
Donde:
p: Pérdidas de Potencia del alimentador (kW)
Pmax: Potencia máxima del alimentador
di: Distancia de la fuente (origen de red) hasta el transformador “i”
ltot: Longitud total del circuito
mp: Momento de pérdidas del alimentador (MW 2km/kW)
Vb: Tensión de referencia del alimentador utilizado para mp
V: Tensión de operación del alimentador
b: Angulo de referencia del alimentador correspondiente al factor de carga
utilizado para mp
: Angulo de factor de potencia del alimentador
37
mp (nd ) (n p ) (2.57)
Donde:
nd: 360/
: Angulo del sector circular del alimentador MT
np: Número de transformadores propios y particulares conectadas
a(rt rr )b (2.58)
r
c d ln( t ) (2.59)
rr
: e
rt: Resistencia del conductor de la troncal del alimentador MT
rr: Resistencia del conductor de la ramal del alimentador MT
a,b,c,d,e : Constantes definidas de acuerdo al valor de rt, rr y distancia de carga
equivalente dceq.
nt
Pfe p fet (2.60)
t 1
Las pérdidas con carga del transformador se calculan utilizando las expresiones
siguientes:
pcu t ( fut )2 pcu t (2.61)
nt
Pcu pcu t (2.62)
t 1
Donde:
38
Cálculo de pérdidas de potencia en redes eléctricas de baja tensión:
Las ecuaciones a utilizar para el cálculo de pérdidas de potencia son diferentes para
cada tipología. En el caso de la Tipología 1 se utiliza la siguiente ecuación:
i 2l 2
pe f (r , l , i, I j ) rl iI j l I j 2 (2.63)
3
nf
p c 3 2pe(r0 , l0 , i0 , I j0 ) (2.64)
10
Donde:
39
l0: Longitud del tramo elemental (lcir/2)
10Snom fu
i0 :Corriente de la línea(A/m)
3Vnomlcir
Ij0: 0 (A)
Las energías que ingresan al segmento provienen de Generación aguas abajo (Ereqs),
de fuentes de transformación (Etransfs); las salidas son las energías de clientes libres
(EformSRAs) en el punto de inicio del segmento y las salidas al final del segmento son
Energía de clientes formales (Eforms) e informales (EformIs) y la energía transferida al
siguiente segmento (Etransfs+1) ó (Einjs+1). Las pérdidas de energía del segmento se
representan como (Es).
2
Ereginj
Ereg Efio
Efio (2.65)
inj
40
Figura 2.13 Balance de energía en un segmento de red
(Fuente: ANEEL-Brasil, 2007[10])
Donde:
Ereginj : Energía regular (consumos registrados) que pasa por el segmento
Este método de cálculo de pérdidas técnicas tiene como propósito obtener los
valores de las pérdidas de potencia y energía de manera global en las redes de media
tensión, transformadores MT/BT, redes de baja tensión, conductores de acometida y
medidores de energía.
Los indicadores que se obtienen con este método son importantes pero no son
suficientes para definir los planes y programas con el objetivo de mejorar el nivel de
eficiencia de la distribución eléctrica.
Los estudios de pérdidas de energía con esta metodología pueden tomar entre seis a
ocho meses y sus resultados son goblales y no proporcionan la localización e
41
identificación de los componentes del sistema que contribuyen al incremento de las
pérdidas de energía.
Figura 2.14 Diagrama unifilar de una red de baja tensión con un tramo único
(Fuente: Elaboración propia)
La idea es que solamente a través de las mediciones de tensión en el nodo final (V1),
tensión en el origen (Vs) y potencia activa en el origen del sistema de distribución (Ps),
sea posible estimar las pérdidas de potencia en condiciones de máxima demanda.
42
De la ecuación anterior la resistencia R en función de las mediciones de Potencia y
Tensión será:
Vs Vs V1
R (2.67)
Ps
P V V V P V V
2
Plínea max
s s s 1 s s 1 (2.68)
Vs Ps Vs
% Plínea max
Vs V1 100%
(2.69)
Vs
Una vez que se tienen los valores de las Pérdidas de Potencia en horas punta, se
complementa esta información con los factores de carga y factores de pérdidas de la red
de distribución de baja tensión. Por último, se calculan las pérdidas de energía y los
valores porcentuales con respecto a la energía distribuida en baja tensión.
Si bien este método permite conocer el nivel de las pérdidas de potencia y energía
en un sistema de distribución de baja tensión, no permite identificar las partes del sistema
(tramos de red) que contribuyen al incremento o excesos de pérdidas técnicas. Asimismo,
para obtener una mejor precisión es necesario realizar mediciones de tensiones en gran
parte de los nodos de carga y ello no es práctico por que incrementa el tiempo requerido
para el estudio y un consecuente aumento del costo.
Esta es una de las metodologías más utilizadas por las empresas para la estimación
de pérdidas de técnicas de energía y generalmente les toma entre ocho a doce meses
concluir dichos estudios dependiendo del tamaño de los sistemas eléctricos de
distribución.
43
2.4 Datos del sistema eléctrico de distribución
En esta parte se hace una revisión general de los datos que son de utilidad para
determinar las pérdidas técnicas y no técnicas de las empresas de distribución en el
sector eléctrico peruano, así como las fuentes de donde pueden ser obtenidos estos
datos.
Las empresas eléctricas del Perú, en general, cuentan con sistemas de información
georeferenciado (GIS) y sistemas de información comercial. Entre los sistemas de
información georeferenciado que se utilizan se cuentan Small World, ArcMap, ArcView,
MapInfo y también Autocad. Asimismo, sus bases de datos de información comercial los
tienen en SQL Server, Oracle, FoxPro, etc. Lo más importante, es que el Osinergmin
tiene los datos georeferenciados de las redes de las empresas de distribución en MapInfo
y SQL Server. Y más importante aún, para quien necesite usar dicha información, es que
el Osinergmin los mantiene actualizados como parte de sus funciones de regulación de
tarifas para tener control del cambio en el Valor Nuevo de Reemplazo (VNR) y la
facturación de consumos de energía a todos los usuarios de electricidad.
Para conocer los datos disponibles de los sistemas de distribución y cómo es que
están organizados se puede optar, entre otras opciones, por lo siguiente:
- Solicitar dicha información directamente a las empresas eléctricas de distribución.
- Buscar información pública de las empresas en OSINERGMIN que como ente
regulador de tarifas eléctricas tiene información necesaria para cumplir sus funciones.
Gran parte de dicha información es de libre acceso en el portal web de Osinergmin.
44
AVNR COyM
VAD (2.70)
D max
Donde:
La Ley de Concesiones Eléctricas aprobada con Decreto Ley N° 25844 del año 1992,
en el título “Precios Máximos de Distribución” y sus respectivos artículos del N° 63 al 81,
establecen las obligaciones, los considerandos y la forma de definir los precios máximos
de distribución.
Es así que, el Osinergmin para definir el VAD y cumplir sus funciones de fijación de
tarifas cuenta con información de los sistemas eléctricos de distribución de las empresas
eléctricas mínimas requeridas para fijar precios, por lo que es evidente que las empresas
cuentan como mínimo con esta información además de tener otros datos de sus sistemas
eléctricos.
La Figura 2.15 muestra las instalaciones que componen los grupos de media y baja
tensión del sistema eléctrico de distribución.
46
- Salidas de Media tensión, tramos de media tensión, equipos de protección y
seccionamiento, nodos de enlace, estructuras y coordenadas de ubicación geográfica.
- Subestaciones de distribución en media y baja tensión (SED MT/BT), subestaciones
elevadora/reductora y coordenadas de ubicación geográfica.
- Salidas de baja tensión, tramos de baja tensión, acometidas, puntos de conexión de los
suministros y coordenadas de ubicación geográfica.
A diciembre del año 2011 el mercado eléctrico de distribución estuvo compuesto por
más de 5.5 millones de usuarios finales de los cuales más del 99% corresponden al
sector de baja tensión y más del 92.3% son usuarios con tarifa BT5B con contadores que
registran únicamente la energía total consumida [6].
47
En la misma ocasión de cada cuatro años, la GART realiza estudios de
caracterización de carga en los sectores típicos en las que están divididos las diferentes
regiones en función de la densidad de carga y otros factores establecidos por el
Osinergmin [8]. El último estudio realizado fue el año 2008-2009 para la fijación de tarifas
2009-2013 y como muestra se presenta en el anexo III la información resumida del sector
típico 1 correspondiente a la zona de Lima (en la concesión de EDELNOR) [9].
Para usuarios con otras tarifas se utilizan los registros de potencia horaria
disponibles por la Distribuidora, a fin de minimizar los errores de muestreo para el
desarrollo del estudio de caracterización de carga.
En la Tabla A3.2 del Anexo III se detalla el número de mediciones realizadas en cada
una de las tarifas y estratos estudiados para fijar las tarifas 2009-2013.
De esta forma se obtienen los factores típicos (factor de carga, factor de potencia,
factor de pérdidas, etc.) de cada usuario de acuerdo a su opción tarifaria, que permiten
caracterizar la carga (demanda) de cualquier sistema eléctrico de distribución para
cualquier estudio que requiere de datos de las demandas.
48
Los resultados del estudio del año 2009 se muestran en el anexo III, figuras A3.2 al
A3.5.
Con el objeto de garantizar una adecuada representatividad de los factores calculados se
compara la curva obtenida a partir de los datos de las mediciones realizadas versus la
curva de compra de la Empresa Distribuidora. En la Figura A3.4 del Anexo III se presenta
un detalle de las curvas comparadas en el cual se aprecia que los datos inferidos a partir
de las mediciones reproduce de manera aceptable la curva de compra anual de la
empresa.
Siendo los factores típicos de interés para los estudios de pérdidas de energía las
correspondiente a FC (factor de carga), Fperdidas (Factor de pérdidas), FCPP (factor de
coincidencia en horas de punta), NHUBTPP (número de horas de utilización en punta),
NHUBTFP (número de horas de utilización en fuera de punta), etc.
49
2.4.4.1 Conductor Óptimo
Existen diversos estudios referentes a cuál es el calibre óptimo de conductor para
una determinada corriente de carga que pasa a través del conductor, la más cercana y
que se pueden encontrar en el portal de Osinergmin, son los resultados de estudios
realizados con motivo de Estudios de Valor Agregado de Distribución (VAD) que se hace
en el Perú cada 4 años para fijas las tarifas eléctricas en Distribución [2]. Estos estudios
indican que el conductor óptimo es aquel que minimiza la siguiente ecuación de costos:
Donde:
Costo_Inv: Costos de Inversión (materiales, mano de obra y gastos indirectos)
Costo_OyM: Valor presente de costos de operación y mantenimiento de la red para
cada año, durante la vida útil del conductor.
Costo_Pérdidas: Valor presente de los costos de pérdidas de potencia y energía que se
determinan a partir de:
Donde:
I: Corriente Máxima por el conductor
R: Resistencia del Conductor
Tep: Tiempo equivalente de Pérdidas (equivalente a T*Fpérdidas)
T: Tiempo en horas.
CE: Costo promedio de la energía
CP: Costo promedio de la potencia
K1: Constante representativa de la corriente máxima para una densidad
uniforme de carga. (Igual a 1.73 para una carga distribuida
uniformemente).
50
Parte de los resultados del estudio VAD del año 2009 [2] se muestran en el Gráfico
A3.5 del Anexo III de la tesis.
% Carga
80%
60%
40%
20%
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 25 26 27 28 29 30
Año
Aluminio Cobre
Figura 2.16 Evolución de carga del conductor requerido a lo largo de su vida útil para
que no genere excesos de pérdidas de energía
(Fuente: Elaboración propia)
Este análisis nos muestra que para no caer en excesos de pérdidas de energía la
carga inicial deberá ser superior al 20% de la capacidad nominal del conductor en su
puesta en operación y en los últimos años de su vida útil terminar con una carga menor al
70% de su capacidad nominal.
51
2.4.4.2 Optimización de Transformador de Distribución (MT/BT)
Existen diversos estudios para encontrar la solución al tamaño o capacidad de
potencia de los transformadores de distribución (MT/BT) para la atención de una
determinada carga, y se pueden encontrar fácilmente en el portal de Osinergmin, que son
los resultados de estudios realizados con motivo de Estudios de Valor Agregado de
Distribución (VAD) que se hace en el Perú cada 4 años para fijas las tarifas eléctricas en
Distribución [2].
Donde:
CCT_MTB: Costo total mínimo de la subestación MT/BT (todos los costos
expresados en valor presente).
CINI: Costos iniciales de la subestación de distribución MT/BT.
CFUTURO: Costos de instalaciones futuras en la subestación MT/BT.
COyM: Costos de operación y mantenimiento de la subestación MT/BT.
CPERDIDAS: Costos por pérdidas de energía y potencia en la subestación MT/BT.
Donde:
CCENTRO: Costo de instalación o montaje de la subestación.
CTRAFOS: Costo del transformador MT/BT.
Donde:
CFTRAFOi: Valor presente de los costos futuros de transformadores a instalar.
CFCENTROi: Valor presente de los costos futuros de subestaciones a construir.
CFCAMBIOi: Valor presente de los costos futuros de cambios de transformador a
realizar.
52
Nota: El incremento de la demanda se cubre mediante el reemplazo de
transformadores por otros de mayor potencia y/o mediante la instalación de nuevas SED.
El valor presente de los costos futuros por cambios de transformadores, está dado
por el valor presente del costo inicial de cambio de un transformador y de los costos de
desmontaje y montaje asociados al mismo.
Donde:
PCU: Pérdidas de Potencia en el devanado del transformador.
PFE: Pérdidas de Potencia en el núcleo del transformador.
Tep: Tiempo equivalente de Perdidas (equivalente a T*Fpérdidas)
T: Tiempo en horas.
CE: Costo promedio de la energía
CP: Costo promedio de la potencia
Estos costos se determinan para cada año a lo largo de la vida útil del trasformador y
se llevan al valor presente con una tasa de retorno del 12%. El factor de utilización
cambiará año tras año con el incremento de la demanda.
El Tep depende del factor de pérdidas y este último depende del factor de carga.
Este proceso de optimización se realiza para cada sistema de distribución en particular.
53
% Pérdidas de Energía enTransformadores
(Según % de Carga)
% Pérdidas
de Energía
2.5%
2.0%
1.5%
1.0%
0.5%
0.0%
65%
67%
69%
71%
73%
98%
75%
78%
80%
82%
85%
87%
90%
93%
95%
101%
104%
107%
111%
114%
117%
121%
125%
128%
132%
% Carga
Pérdidas Fijas (%) Pérdidas Variables (%) Pérdidas Totales (%)
140%
120%
100%
80%
60%
40%
20%
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Figura 2.18 Porcentaje de carga del transformador, con factor de carga de 0.5, tasa
de crecimiento de 2.5%, que optimiza las pérdidas de energía en el equipo.
Este análisis nos muestra que para no caer en excesos de pérdidas de energía la
carga inicial deberá ser superior al 60% de la capacidad nominal del transformador de
distribución en su puesta en operación y en los últimos años de su vida útil terminar con
una carga menor al 130% de su capacidad nominal.
55
C API TULO III
METODOLOGÍ A P AR A L A SOLUCIÓN DEL PROBLEM A
Para este efecto, se parte de la premisa que la data común que se puede conseguir
en todas las empresas eléctricas es el correspondiente a la base de datos de VNR (Valor
Nuevo de Reemplazo) y a los reportes mensuales de consumos facturados que son
reportados periódicamente al Osinergmin con motivo de brindar la información para el
FOSE (Fondo de Compensación Social de Electricidad).
Como solución para realizar este análisis, se plantea utilizar el método que se
muestra de forma general en la Figura 6.1 y el flujograma completo en la Figura 6.4. Para
ello se toman como datos entrada para el modelo de flujo de carga en sistemas
desbalanceados, los parámetros y características técnicas de los conductores,
transformadores, contadores de energía y cargas; los resultados de estudios de
caracterización de carga, los consumos de energía de los clientes, datos de protocolos de
prueba de transformadores de distribución, contadores de energía, aisladores, etc.; y la
representación de toda la infraestructura eléctrica con sus coordenadas geográficas. Los
resultados de este análisis se incorporan al balance de energía utilizados para determinar
las pérdidas no técnicas de energía.
Infraestructura
Parámetros de Consumos
Conductores, Facturados
Transformadores,
Medidores, etc. Caracterización
de Carga
57
Los sistemas eléctricos de distribución tienen contadores de energía en los puntos de
transferencia (entradas y salidas) de energía al sistema y en los puntos de consumo, por
lo que para determinar las pérdidas totales de energía se realizan balances de energía,
según se esquematiza en la Figura 3.2.
A partir del conocimiento de los datos que son posibles de obtener de las redes
eléctricas, se definen los datos mínimos requeridos para el análisis de pérdidas técnicas.
Un aspecto importante es que se modelan las redes eléctricas en forma detallada cada
vano por vano para facilitar la identificación de los puntos críticos o los elementos que
más inciden en el exceso de pérdidas de energía y su visualización gráfica
georeferenciado sobre el catastro urbano a diferencia de los tradicionales diagramas
unifilares.
58
- Código del centro de transformación o de la central de generación aislada.
- Código del sistema eléctrico al que pertenece.
- Tensión nominal lado AT, MT, BT.
- Potencia Instalada.
- Demanda máxima.
- Coordenadas geográficas de ubicación.
Estructura (M_POSTEMT):
- Código de la estructura.
- Código del tramo de media tensión y/o tramo de baja tensión al cual soporta.
- Etiqueta de la estructura asignada por la empresa.
- Tipo de Aisladores y cantidad.
- Coordenadas de ubicación.
59
3.1.2 Datos de transf ormadores MT/ BT
60
- Código del neutro que identifica el material, sección y fase del conductor o cable.
- Longitud del tramo.
- Fechas de puesta en servicio.
- Código del tramo anterior o tramo padre.
- Coordenadas geográficas de ambos extremos del tramo.
Acometida (M_ACOMETID):
- Código de la acometida.
- Tramo de baja tensión al que pertenece.
- Coordenadas de ambos extremos del tramo de acometida.
- Longitud del cable de acometida.
- Código del suministro.
- Marca, Modelo de Medidor.
- Código de Tipo de Medidor (electrónico, electromecánico, etc.).
- Tarifa que corresponde al suministro.
- Consumo de energía mensual (kWh).
- Máxima demanda kW si dispone.
- Sector Típico al que pertenece el suministro.
El flujograma de la Figura 3.3 nos muestra el método práctico que de manera general
consiste en:
61
Figura 3.3 Diagrama de flujo de datos que utiliza el método
62
Figura 3.4 Diagrama de procedimiento que utiliza el método práctico
(Fuente: Elaboración propia)
Seguidamente se realiza el cálculo de flujo de carga para las redes de media tensión
de forma similar que para la baja tensión y se guardan los resultados y se emiten los
reportes en formatos tabulares y mapas temáticos.
Las pérdidas de energía calculadas con el método de flujo de carga nos
proporcionan las pérdidas técnicas de energía en el sistema eléctrico de distribución.
Finalmente tomando los datos de balance de energía obtenidos de los registros de
los contadores de energía que nos proporcionan las pérdidas de energía totales se
procede a realizar el cálculo de las pérdidas no técnicas como la diferencia entre las
pérdidas totales y las pérdidas técnicas.
63
transferencia de datos hacia la herramienta de cálculo implementado para la utilización
del método práctico desarrollado. Este proceso se realiza con conexión directa a la Base
de Datos GIS y la Base de Datos de Clientes que está típicamente en SQL Server u otro
motor de base de datos, para que la importación de información sea rápida y
transparente. Igualmente se puede exportar los datos del GIS y los datos de los clientes a
formatos Excel o texto para luego organizar la información de la forma como lo requieren
las herramientas de análisis de flujo de carga.
(3.1)
Donde:
(3.2)
Donde:
Las corridas de flujo de carga se validan con datos de medición que se disponen a la
entrada del sistema de distribución. En términos resumidos se tiene que cumplir:
(3.3)
(3.4)
Donde:
Ps: Potencia a la entrada al sistema de distribución.
Pc: Potencia correspondiente a la sumatoria de las cargas.
L: Pérdidas de Potencia en el sistema de distribución.
Zred: Impedancia equivalente del sistema de distribución.
f(Pc,Zred): L es función de una variable conocida Zred y otra desconocida Pc.
65
A los datos que se tienen de cada tramo de red, transformadores, contadores de
energía, cables de acometida, etc.; se le adicionan los flujos de corriente obtenidos por el
flujo de carga, seguidamente se obtienen las pérdidas de potencia y energía en cada
parte del sistema de distribución aplicando los modelos matemáticos presentados en el
capítulo 2.
Cuando se corre el flujo de carga en media tensión se considera como carga todas
las subestaciones MT/BT, donde dicha carga en lado MT es igual a:
Donde:
Pcu max : Pérdidas de Potencia en el devanado de transformador en condiciones de
máxima demanda
Pfe : Pérdidas de Potencia en el Núcleo del transformador
¿Se disponen
Tomar estos
de datos de este
Inicio SI parámetros para
elemento en
este elemento
particular?
NO
Datos a tomar en
¿Se disponen datos de la Tomar parámetros
cuenta para el Salida hacia
marca al que pertenece el SI de la marca para
cálculo de pérdidas proceso de cálculo
elemento? este elemento
de energía
NO
Tomar parámetros
típicos para este tipo
de elementos
(3.6)
Siendo:
(3.7)
Los resultados se organizan en tablas que guardan además de los resultados las
coordenadas geográficas de cada elemento de red para los reportes de mapas temáticos.
Los mapas temáticos son utilizados por profesionales de todas las disciplinas
sociales y ambientales porque presentan distintos tipos de datos como información visual,
que resulta de más rápida lectura y más sencilla interpretación.
67
Luego que se tienen los valores de pérdidas técnicas y no técnicas de energía en las
redes de distribución, se procede al análisis en entorno GIS, consistente en la
visualización de los resultados en mapas temáticos típicamente pintando de color rojo los
tramos de líneas de media y baja tensión y subestaciones MT/BT, etc. con excesos de
pérdidas de energía.
68
A partir de este análisis se determinan las estrategias de reducción y control de las
pérdidas de energía tanto técnicas como no técnicas.
69
Figura 3.8 Visualización de Líneas MT en KML (Keyhole Markup Language)
(Fuente: Elaboración propia)
70
Figura 3.9 Visualización de Líneas BT en KML
(Fuente: Elaboración propia)
71
3.3 Viabilidad del desarrollo aplicati vo
Dado que el análisis de pérdidas de energía requiere del procesamiento de datos
numerosos y en diversos formatos, es necesaria la programación de un aplicativo que
procese la información de los sistemas eléctricos de distribución conforme a la
metodología planteada.
72
C API TULO IV
AN ÁLISI S Y PRESENTACIÓN DE RESULTADO S
74
En las tablas 4.1 y 4.2 se muestra la longitud de líneas de distribución de media y
baja tensión, respectivamente. Los datos que se consignan por alimentador consisten en
el nombre del alimentador, la subestación a la que pertenecen, su longitud total y su
composición por tramos en conductores de cobre, aluminio o AUT (autosoportados).
La Tabla 4.3 muestra el número de usuarios por alimentadores segmentados por el tipo
de tarifa que poseen. Se aprecia que la mayoría de ellos pertenecen a la tarifa BT5B.
75
Subestación Alimentador BT5B (N° MT2 (N° MT3 (N° MT4 (N° TOTAL (N°
AT/MT MT Clientes) Clientes) Clientes) Clientes) Clientes)
SET002 ALI003 9698 1 9699
SET003 ALI004 6663 1 2 1 6667
ALI005 6311 6311
SET004 ALI006 1350 1350
SET001 ALI001 2728 1 1 2730
ALI002 7098 7098
SET005 ALI007 6409 2 1 6412
ALI008 7560 2 1 7563
ALI009 2017 2017
SET006 ALI010 6728 1 1 6730
ALI011 7729 7729
ALI012 9470 2 9472
Total 73761 5 7 5 73778
% de Pérdida de % Pérdida de
Sub Estación Alimentador % Pérdida de
Energía en la Energía
AT/MT MT Energía Total MT
Línea MT Aisladores
SET001 ALI001 0.93% 0.06% 0.98%
SET001 ALI002 1.65% 0.09% 1.74%
SET002 ALI003 1.24% 0.68% 1.92%
SET003 ALI004 0.95% 0.92% 1.88%
SET003 ALI005 1.11% 0.15% 1.27%
SET004 ALI006 0.61% 0.25% 0.86%
SET005 ALI007 0.63% 0.18% 0.81%
SET005 ALI008 0.88% 0.22% 1.09%
SET005 ALI009 0.31% 0.23% 0.54%
SET006 ALI010 2.93% 0.22% 3.16%
SET006 ALI011 0.45% 0.09% 0.54%
SET006 ALI012 1.88% 0.08% 1.96%
TOTAL 1.28% 0.20% 1.48%
76
En promedio las pérdidas técnicas de energía en líneas de media tensión se
encuentran en el orden de 1.48%, valor coherente con los niveles que son reconocidos
en las tarifas eléctricas para sistemas de distribución ubicadas en las zonas semi-urbanas
y rurales.
0.20%
% Pérdida de
% de Pérdida de % Pérdida de Energía % de Pérdida Técnicas
Sub Estación AT/MT Alimentador MT Energía en
Energía en la Línea BT en Acometidas de Energía en BT
Medidores
SET001 ALI001 1.61% 0.08% 0.51% 2.20%
SET001 ALI002 3.42% 0.06% 0.75% 4.23%
SET002 ALI003 0.82% 0.03% 2.41% 3.25%
SET003 ALI004 0.55% 0.01% 2.31% 2.87%
SET003 ALI005 3.35% 0.03% 1.16% 4.55%
SET004 ALI006 0.96% 0.02% 1.82% 2.80%
SET005 ALI007 1.24% 0.05% 0.96% 2.26%
SET005 ALI008 1.49% 0.05% 0.93% 2.47%
SET005 ALI009 1.25% 0.03% 1.65% 2.92%
SET006 ALI010 1.77% 0.04% 0.97% 2.77%
SET006 ALI011 2.85% 0.03% 0.82% 3.71%
SET006 ALI012 2.20% 0.05% 0.67% 2.92%
TOTAL 2.08% 0.05% 0.99% 3.11%
77
En baja tensión las pérdidas técnicas de energía en líneas están alrededor del
2.08%, ella es debido a que las redes se encuentran con poca carga, y este hecho
contribuye también a que las pérdidas de energía en los medidores de energía sean
significativas. En las zonas rurales los consumos de energía típicamente ocurren en las
horas de la noche y en la mañana, pero en el día el consumo de la electricidad se reduce
drásticamente, por lo que en dicho periodo lo que existe en las redes es básicamente el
consumo de los medidores de energía.
2.08%
0.99%
0.05%
Los valores mostrados en forma porcentual vienen a ser el resultado de dividir las
pérdidas técnicas de energía en media tensión (“conductores y aisladores”),
transformadores de distribución (“núcleo y devanado”) y baja tensión (“conductores de
red de baja tensión, conductores de acometidas y contadores de energía”) entre la
energía total que ingresa al sistema de distribución en estudio.
78
% de Pérdida de % Pérdida de % Pérdida de % de Pérdida
Sub Estación Alimentador
Energía en la Energía en Energía en Redes Técnicas de
AT/MT MT
Línea MT Transformadores BT Energía Total
SET001 ALI001 0.98% 1.57% 1.37% 3.93%
SET001 ALI002 1.74% 1.68% 2.77% 6.20%
SET002 ALI003 1.92% 5.12% 2.25% 9.29%
SET003 ALI004 1.88% 5.73% 1.75% 9.35%
SET003 ALI005 1.27% 2.15% 3.04% 6.46%
SET004 ALI006 0.86% 2.27% 1.00% 4.13%
SET005 ALI007 0.81% 2.05% 1.38% 4.24%
SET005 ALI008 1.09% 2.64% 1.76% 5.49%
SET005 ALI009 0.54% 1.29% 0.82% 2.65%
SET006 ALI010 3.16% 2.40% 1.84% 7.40%
SET006 ALI011 0.54% 1.78% 2.67% 4.99%
SET006 ALI012 1.96% 1.76% 2.14% 5.86%
TOTAL 1.48% 2.27% 2.04% 5.78%
Nota: % de Pérdidas en BT calculados con referencia a la energía inyectada en la salida MT.
2.27% 2.04%
1.48%
79
práctico visto en el capítulo 3 y las pérdidas no técnicas obtenidas se muestran en la
Tabla 4.7 y la Figura 4.5.
6.00%
4.00% 3.11%
2.00%
0.00%
% Pérdidas de Energía % Pérdidas de Energía % Pérdidas de Energía
Total BT Técnica BT No Técnica BT
80
4.3 Análisis e interpret ación de los resultados
En el sistema estudiado se observa que el nivel de Pérdidas No Técnicas está entre
el 5% y 9%, muy por encima de las Pérdidas Técnicas, por lo que se observa que se
tienen exceso de Pérdidas No Técnicas de Energía, especialmente en los alimentadores
ALI004, ALI005, ALI007, ALI008, ALI009, ALI010 y ALI012.
Analizando los resultados obtenidos en la tabla 4.4 vemos que las pérdidas técnicas
en media tensión se deben principalmente por las pérdidas producidas en los
conductores. En general los valores obtenidos de pérdidas técnicas en media tensión se
encuentran alrededor de los niveles reconocidos por las tarifas eléctricas.
Analizando los resultados obtenidos en baja tensión en la tabla 4.5, vemos que tiene
mayor participación en las pérdidas técnicas las pérdidas producidas en los conductores
y contadores de energía. Especialmente en las zonas rurales se observa que se presenta
mayores niveles de pérdidas técnicas en los contadores de energía.
En la tabla 4.6 se tienen los resultados totales que nos muestran que los excesos de
pérdidas técnicas se encuentran en los transformadores de distribución principalmente en
los alimentadores con cargas dispersas.
Finalmente al analizar la tabla 4.7 vemos que los excesos de pérdidas de energía
para este sistema se deben a las pérdidas no técnicas.
Con los resultados mostrados, que son los mismos que se presentaron a la Alta
Dirección (Gerencia de Distribución) de la empresa donde se realizó el trabajo, la
organización tomó decisiones de acciones a ejecutar para la reducción de pérdidas de
energía. Antes del estudio, los funcionarios de esta empresa tenían la percepción que las
pérdidas de energía totales elevadas se debían a las pérdidas técnicas por lo que tenían
la intención de realizar inversiones en cambio de conductores y montaje de nuevos
circuitos de distribución para la reducción de pérdidas técnicas de energía. Luego de ver
los resultados obtenidos de la estimación de pérdidas técnicas se encontró que el exceso
de pérdidas de energía se debe a las pérdidas no técnicas, por lo que la dirección de la
empresa tomó la decisión de aplicar las acciones de reducción de pérdidas de energía no
técnicas tales como auditorias a los procesos de toma de lectura de contadores de
energía (partes del proceso que son realizados por empresas externas), contrastación de
contadores de energía, etc.
81
CONCLUSIONES
82
que el análisis en entorno GIS es clave para entender este problema y también
para plantear sus soluciones.
83
RECOMEND ACIO NES
3. A Osinergmin-GART: consignar de forma explícita las fases (R, S, T, RS, ST, RT,
SRT, etc.) a las que están conectadas las líneas MT, BT, transformadores MT/BT
y cargas (consumidores finales) en la información de VNR de las instalaciones de
distribución de media y baja tensión. Con ello se espera que la información del
VNR tenga mayor relevancia para estudios desbalanceados de flujo de carga,
pérdidas de energía, calidad de tensión, etc.
84
BIBLIOGR AFÍ A
[1] OLADE (1990). Manual Latinoamericano y del Caribe para el Control de Pérdidas
Eléctricas. 1ra. Edición. Bogotá, Colombia. PP 256.
[2] EDELNOR (2001). Informe Final del Estudio de Costos del VAD del Sector Típico 1 –
Volumen 1 Capítulo 5, Resultados. Lima, Perú. PP 70.
[4] Muchoz, C.M. (1989). Flujo de Potencia Trifásico para Sistemas de Distribución.
Tesis de grado de ingeniero. Pontificia Universidad Católica de Chile, Chile. PP 155.
[5] Mok, H. M., Elangovan, S., Salama, M.M.A., Longjian C. (2000). Power Flow
Analysis for Balanced and Unbalanced Radial Distribution Systems [en línea]. [Fecha
de acceso 20 de marzo de 2013]. URL disponible en:
http://www.researchgate.net/publication/252701407_Power_Flow_Analysis_for_Bala
nced_and_Unbalanced_Radial_Distribution_Systems
85
[8] Osinergmin-GART (2005). Estudio de Caracterización de la Carga de los Sistemas
de Distribución Eléctrica – Sectores Típicos 1, 2, 3, 4, 5 y E. [Fecha de acceso 25 de
febrero de 2013]. Resolución de Fijación del VAD. URL disponible en:
http://www2.osinerg.gob.pe/ProcReg/VAD/VAD2005/Descarga/Publicacion/InformeV
AD%20Anexo10-Parte%201%20de%206.pdf
[11] YUSOFF*, M., BUSRAH*, A., MOHAMAD*, M., TENG, M. (2009). A Simplified
Approach in Estimating Technical Losses in TNB Distribution Network Based on Load
Profile and Feeder Characteristics. [Fecha de acceso 08 de enero de 2013].
http://www.wseas.us/e-library/conferences/2010/Penang/MMF/MMF-15.pdf
[14] Troy (1999). Reference Guide MapX 4.0. MapInfo Corporation, New York, USA. PP
341.
[16] GOLKAR, M. A. (2007). A Novel Method for Load Flow Analysis of Unbalanced
Three–Phase Radial Distribution Networks. [Fecha de acceso 23 de octubre de
86
2012]. URL disponible en: http://journals.tubitak.gov.tr/elektrik/issues/elk-07-15-3/elk-
15-3-2-0607-4.pdf
[18] Zhicheng, G., Yingke, M., Liming, W., Ruihai L., Hua, W., Yi, M (2009). Leakage
Current and Discharge Phenomenon of Outdoor Insulators. [Fecha de acceso 10 de
mayo de 2013]. URL disponible en: http://www.ijeei.org/docs-
17922250284b44498900d74.pdf
87
ANEXO A
Sistemas de Información Geo gráfico ( GIS) y formatos más utilizados
Una herramienta GIS que está al alcance de todos viene a ser la aplicación de
Google Earth, que nos permite situarnos en cualquier lugar de la Tierra para ver
88
imágenes de satélite, mapas, imágenes de relieve y edificios 3D, o para explorar desde
vista aéreas, etc. El formato de archivo GIS basado en XLM que utiliza esta aplicación es
el formato KML, que facilita el intercambio de datos GIS. En la tesis se utiliza este formato
para visualizar los resultados de las pérdidas de energía en mapas temáticos en KML
(del acrónimo en inglés Keyhole Markup Language), éste fue desarrollado para ser
manejado con Keyhole LT, precursor de Google Earth (Google adquirió Keyhole LT en
octubre de 2004 tras lanzar su versión LT 2).
Otras aplicaciones GIS como es el caso de MapInfo, utilizan archivos con extensión
TAB, en la que se guardan los datos de los objetos geográficos o georeferenciados y sus
respectivos datos tabulares.
89
ANEXO B
Breve descripción de los p rogramas principales en Visual Studio 6.0
Permite entrar a la base de datos del sistema que ya cuenta con los datos obtenidos
con el procedimiento precedente para elegir la red eléctrica de media o baja tensión que
será utilizado para el Análisis de Pérdidas de Energía. De esta forma se visualiza en la
pantalla solamente la data de la red eléctrica objeto de estudio para fines de revisión,
actualización y análisis de flujo de carga, etc. Para la visualización gráfica en entorno
georeferenciado este procedimiento utiliza otras funciones adicionales denominadas
“DibujaRedMT” y “DibujaRedBT”.
Este procedimiento se encarga del cálculo iterativo de flujo de carga para la red de
media o baja tensión seleccionada.
91
Procedimiento: Iteración (Código en Visual FoxPro 6.0)
92
ANEXO C
Información de los Sistemas Eléctricos de Distribución del Perú
A continuación se listan los datos que son posibles de conseguir con facilidad de los
sistemas eléctricos de distribución en el Perú.
93
- Coordenadas de los vértices de la polilínea.
Nodo de Enlace:
- Código del nodo.
- Código del tramo al cual está conectado.
- Código de la subestación a la que pertenece.
- Tensión nominal.
- Tipo de propietario.
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas del vértice de la polilínea del tramo que llega o sale de la subestación.
Estructura:
- Código de la estructura.
- Código del tramo de media tensión y/o tramo de baja tensión al cual soporta.
- Etiqueta de la estructura asignada por la empresa.
- Código VNR que identifica el tipo de estructura y altura.
- Tipo de propietario.
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas de ubicación.
94
- Tensiones nominales.
- Demandas máximas del servicio particular y alumbrado público.
- Potencia instalada.
- Etiqueta de la subestación asignada por la empresa.
- Código del subcomponente asociado.
- Código VNR que identifica el tipo de subestación, potencia y fases.
- Tipo de propietario.
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas de ubicación.
Subestaciones Elevadora/Reductora:
- Código de la subestación.
- Código de la salida de MT a la que pertenece.
- Tensiones nominales.
- Potencia instalada.
- Etiqueta de la subestación asignada por la empresa.
- Código del subcomponente asociado.
- Código VNR que identifica el tipo de subestación, potencia y fases.
- Tipo de propietario.
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas de ubicación.
- Código de la salida.
- Código de la subestación de distribución MT/BT a la que pertenece o de la central de
generación aislada de ser el caso.
- Tensión nominal.
- Nombre de la salida asignado por la empresa.
- Tipo de servicio (servicio particular o alumbrado público).
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas de ubicación.
95
Tramo de Baja Tensión:
Acometida:
- Código de la acometida.
- Tramo de baja tensión al que pertenece.
- Estado del elemento respecto a la fijación anterior.
- Fechas de puesta en servicio y retiro.
- Coordenadas de los vértices de la polilínea.
96
Tablas y Gráficos del Estudio VAD Estudios de Caracterización de la Carga de
los Sectores Típicos 1, 2, 3, 4, 5, Especial y SER [9] y de otros estudios:
97
TARIFA ESTRATO MEDICIONES
BT5B 1 113
BT5B 2 146
BT5B 3 137
BT5B 4 95
BT2 1 178
BT3 1 1254
BT4 1 1413
BT5-A 1 173
MT1 1 77
MT2 1 116
MT3 1 619
MT4 1 303
98
Tabla A3.3 Evolución número de usuarios de electricidad y ventas de energía (Perú)
(Fuente: Osinergmin-GART, Tarifas y mercado eléctrico)
99
Figura A3.1 Diagrama de carga por tarifas en días hábiles
(Fuente: Osinergmin, estudio de caracterización de carga sector típico 1, 2009)
100
Figura A3.2 Diagrama de carga por tarifas en días sábados
(Fuente: Osinergmin, estudio de caracterización de carga sector típico 1, 2009)
101
Figura A3.3 Diagrama de carga por tarifas en días domingos
(Fuente: Osinergmin, estudio de caracterización de carga sector típico 1, 2009)
102
En resumen los factores típicos necesarios para caracterizar la carga obtenidos
fueron los siguientes:
103
Tabla A3.4 Factores típicos obtenidos para cada tipo de Tarifa
(Fuente: Osinergmin, estudio de caracterización de carga sector típico 1, 2009)
Donde:
104
Figura A3.5 Costos totales de conductores, estudio VAD 2009
(Fuente: Osinergmin, informe final del estudio VAD sector típico 1, 2009)
105