UNIVERSIDAD NACIONAL EXPERIMENTAL POLITÉCNICA
“ANTONIO JOSE DE SUCRE”
VICERRECTORADO PUERTO ORDAZ
DEPARTAMENTO DE INGENIERIA ELECTRICA
SISTEMAS DE POTENCIA I
SEMESTRE 2020-1
Profesor:
Ing. Aurilse Ortiz
Alumno:
Germán Manrique
CI 26549026
Octubre, 2020
Tabla de Contenido
Pág.
Introducción…………………………………………………………………………………….…….3
Contenido……………………………………………………………………………………………..4
1. Parte I……………………………………………………………………………………….....4
1.1 Análisis de Flujo de Cargas…………………………………………………………....4
1.2 Análisis de Cortocircuito………………………………………………………………13
2. Parte II (Descripción de la Actividad)…………………………………………………….16
2.1 Análisis de Flujo de Cargas Mediante ETAP 12.6.0……………………………….23
2.2 Análisis de Cortocircuito Mediante ETAP 12.6.0…………………………………...26
Conclusiones……………………………………………………………………………………….30
Referencias…………………………………………………………………………………………31
2
Introducción
Cuando se tiene un sistema eléctrico de potencia, es necesario conocer su estado y
comportamiento, y así poder tomar las decisiones correspondientes para garantizar el
suministro de la energía eléctrica con calidad, seguridad y confiabilidad. Por tal motivo se
realizan los análisis de flujo de carga.
El principal objetivo de realizar un análisis de flujo de Carga es que partiendo de una
condición de demanda y de generación conocida, se determina el comportamiento de la red
en términos de tensiones y flujos de potencia.
Los métodos utilizados comúnmente para realzar dicho estudio son los de Gauss-
Seidel, Newton Raphson y Newton Raphson Desacoplado Rápido; el primero y el último se
utilizarán para esta experiencia.
Para la experiencia a realizar se tiene un Sistema Eléctrico de Potencia, el cual se
estudiara con el fin de conocer los niveles de tensión en las barras y el nivel de corriente de
cortocircuito (simétrica y asimétrica) cuando ocurren fallos en la barra de carga.
La herramienta a utilizar es el ETAP 12.6.0 el cual es uno de los software más
eficiente que se utilizan hoy en día para realizar dichos procedimientos. El análisis se
realizará analíticamente y después se realizara una comparación con los datos extraídos del
programa.
3
Contenido
Parte I. (Valor 5ptos)
1.1 Aplicando los conocimientos de análisis de flujo de carga. Determine de
forma analítica, considerando 2 iteraciones y comparando con una Tolerancia de 10-3:
a) Vector de Tensiones mediante el método de GAUSS-SEIDEL
b) Vector de Tensiones mediante el método de NEWTON RAPHSON DESACOPLADO
RAPIDO
Pasando los datos a por unidad (pu):
( ⁄ ⁄ ) ⁄
( ⁄ ⁄ ) ⁄
( ) ( )
( ) ( )
Buscando la Matriz de Admitancia (YBUS)
4
La matriz YBUS de secuencia positiva, negativa y cero será:
La matriz de impedancia es la inversa de la matriz YBUS (ZBUS=Y-1BUS):
5
Gauss-Seidel
Barra Tipo (Barra) Voltaje Inicial
1 Slack
6 PQ
3 PQ
4 PQ
5 PV
Vector de Tensiones Iniciales:
Potencia Inyectada en las Barras:
La potencia inyectada en las barras esta denotada por la potencia generada menos la
potencia demandada, esto es:
Iteración 1:
*( ) ( )+
[( ) [ ]]
*( ) ( )+
[( ) [ ]]
6
*( ) ( )+
[ [ ]]
*( ) ( )+
[ ]
Iteración 2:
*( ) ( )+
[( ) [ ]]
*( ) ( )+
[( ) [ ]]
*( ) ( )+
[ [ ]]
*( ) ( )+
7
[ ]
Resultados Gauss-Seidel
[ ] [ ] [ ] [ ]
[ ] [ ] [ ] [ ]
Newton Raphson Desacoplado Rápido (NRDR)
Se utilizaran los mismos valores iniciales utilizados en el método de Gauss-Seidel, el
sistema a resolver será:
[ ] [ ][ ] [ ] [ ][ ]
Dónde:
[ ]
Iteración 1:
La potencia activa calculada será:
∑| || || |
∑| || || |
8
∑| || || |
∑| || || |
[ ] [ ] [ ]
Los errores serán:
Ahora la a potencia reactiva calculada:
∑| || || |
∑| || || |
∑| || || |
[ ] [ ] [ ]
Los errores serán:
9
Ahora se debe resolver el sistema:
[ ][ ] [ ][ ]
[ ] [ ]
[ ] [ ][ ]
[ ] [ ] [ ] [ ]
El voltaje será:
[ ][ ] [ ][ ]
[ ] [ ]
[ ] [ ][ ]
[ ] [ ]
( )
( )
10
( )
( )
[ ] [ ] [ ] [ ]
Iteración 2:
La potencia activa calculada será:
∑| || || |
∑| || || |
∑| || || |
∑| || || |
[ ] [ ] [ ]
Los errores serán:
11
Ahora la a potencia reactiva calculada:
∑| || || |
∑| || || |
∑| || || |
[ ] [ ] [ ]
Ahora:
[ ] [ ][ ]
[ ] [ ] [ ] [ ]
El voltaje será:
[ ] [ ][ ]
[ ] [ ]
( )
12
[ ] [ ] [ ] [ ]
1.2 Determine de forma analítica, considerando análisis de cortocircuito:
a) Corriente de Falla en Secuencia y Fase considerando una Falla Trifásica en la
Barra de Carga
Se debe saber que se tomará en cuenta la impedancia de thevenin de cortocircuito
vista desde la barra 3 al momento del fallo. Esta impedancia no es más que la impedancia
de posición (3,3) de la matriz de impedancia ZBUS (ZThe33). Para el voltaje Prefalla o de
thevenin (VThe33) se considerará el voltaje del vector inicial, para este caso (1 ángulo
cero).
13
La corriente será:
( )
Para calcular las corrientes de secuencia se debe multiplicar por la matriz de Fortesscue:
[ ] [ ][ ] [ ][ ]
[ ] [ ]
Para pasar a Real se multiplica por la corriente base:
√ √
{ }
{ }
b) Corriente de Falla en Secuencia y Fase considerando una Falla Monofásica a
tierra en la Barra de Carga
El circuito correspondiente de la falla, quedará:
𝑍𝑇 𝑍𝑇 𝑍𝑇 𝑗
𝐼𝐹
𝑍𝑇 𝑗
𝐼𝐹
𝑝𝑢 𝐼𝐹
𝐼𝐹 pu
𝐼𝐹
𝐼 𝐼 𝐼 𝑝𝑢 𝑆𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎
Las corrientes de fase serán:
14
[ ] [ ][ ] [ ] [ ][ ]
[ ] [ ]
[ ] [ ] [ ] [ ]
c) Corriente de Falla en Secuencia y Fase considerando una Falla Bifásica en la
Barra de Carga
El circuito equivalente de la falla será:
𝑍𝑇 𝑍𝑇 𝑗
𝐼
𝑍𝑇 𝑗
𝐼 𝑝𝑢 𝑆𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎
𝐼 𝐼 𝑝𝑢 𝑆𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎
𝐼 𝑝𝑢 𝑆𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎
Las corrientes de fase serán:
[ ] [ ][ ] [ ] [ ][ ]
[ ] [ ]
[ ] [ ] [ ] [ ]
d) Corriente de Falla en Secuencia y Fase considerando una Falla Bifásica a tierra
en la Barra de Carga
El circuito equivalente de la falla será:
15
Las corrientes de fase serán:
[ ] [ ][ ] [ ] [ ][ ]
[ ] [ ]
16
[ ] [ ] [ ] [ ]
II Parte. (Valor 10ptos)
Para realizar la simulación del sistema eléctrico de potencia (SEP) asignado,
usaremos el software ETAP 12.6.0 Este es un programa que es utilizado de manera
frecuente por estudiantes de ingeniería eléctrica e inclusive de otras carreras, consta de un
alto repertorio de funciones que ayudan a realizar estudios de alta gama.
A continuación se describen una serie de pasos que se realizaron para simular el
sistema eléctrico en estudio:
1.) Lo primero es abrir el simulador ETAP 12.6.0 y crear un nuevo proyecto (file).
2.) Una vez creado el nuevo proyecto, se empieza a diseñar el diagrama unifilar
seleccionando los equipos que componen el sistema de potencia (transformadores, líneas de
transmisión, cargas, generadores y buses) (ver figura).
17
3.) Luego se cargan los datos de los equipos que aparecen en el sistema eléctrico, como se
muestra a continuación:
El Generador 1 está El Generador 2 está
conectado a la barra conectado a la barra de
slack voltaje controlado PV
Los transformadores tendrán las siguientes características (no poseerán cambiadores de
tap):
18
Las lineas de transmisión tendrán distinta distancia,pero igual impedancia:
19
Los buses operaran bajo ciertas condiciones (se considerará el vector inicial de
tensiones), esto es:
20
4.) Terminado la carga de datos, se procede a realizar el estudio de flujo de carga,
para esto se debe seleccionar en la parte superior el modulo flujo de carga. Además, se debe
seleccionar el método a utilizar (para este caso se utilizará Gauss-Seidel y Newton Raphson
Desacoplado Rápido); como se muestra a continuación:
21
5.) Los resultados obtenidos del flujo de carga mediante el método de Gauss-Seidel
fueron los siguientes:
Las barras están operando bajo los niveles adecuados, una muestra de ello se puede
apreciar a continuación:
22
Resumiendo se tiene:
[ ] [ ]
23
6.) Los resultados obtenidos del flujo de carga mediante el método de Newton
Raphson Desacoplado Rápido (NRDR), fueron los siguientes:
[ ] [ ]
24
7.) Ahora se debe realizar el estudio de cortocircuito, se considerarán cortocircuitos
trifásicos, monofásicos a tierra, bifásicos entre fases y bifásicos a tierra. Para esto se debe
seleccionar el bus en donde se va a realizar la falla (para este caso la barra de carga “3”)
como se muestra en la figura.
8.) Además de seleccionar el bus de falla, se debe seleccionar el tipo de cortocircuito a
realizar (ver figura) y la corriente que quiere mostrar después de la falla (secuencia o de
fase).
25
9.) Para el cortocircuito trifásico se obtuvieron los siguientes resultados:
[ ] [ ] [ ] [ ]
26
10.) Cuando se realizó el cortocircuito monofásico a tierra, se obtuvieron los siguientes
valores:
Corrientes de Secuencia:
Corrientes de Fase:
[ ] [ ] [ ] [ ]
27
11.) Al realizar el cortocircuito bifásico, se obtuvieron los siguientes valores:
Corrientes de Secuencia:
Corrientes de Fase:
[ ] [ ] [ ] [ ]
28
12.) Al realizar el cortocircuito bifásico tierra, se obtuvieron los siguientes valores:
Corrientes de Secuencia:
Corrientes de Fase:
[ ] [ ] [ ] [ ]
29
Conclusiones
El problema del flujo de consiste en obtener las tensiones en todos las barras y los
flujos de potencia por todos los elementos de un sistema eléctrico de potencia, es decir, en
determinar las condiciones de operación en régimen permanente de un sistema de potencia.
Para ello se toma como datos de partida el consumo, la inyección de potencia en todas las
barras de la red y el modelo eléctrico de la misma.
El estudio del problema de flujo de carga se realiza a través de métodos iterativos
como Gauss-Seidel, Newton Raphson Desacoplado Rápido (NRDR), entre otros. Tiene la
característica de que si se tiene un sistema eléctrico muy grande lo reduce a uno muy
pequeño a través de la matriz de admitancia nodal (YBUS).
Para la experiencia realizada se manejó el método de Gauss-Seidel y Newton
Raphson Desacoplado Rápido (NRDR). Además del programa computacional ETAP 12.6.0,
el cual es muy utilizado actualmente para desarrollar actividades de índole eléctrico.
De la experiencia realizada se puede decir que los valores obtenidos analíticamente
mediante los métodos de Gauss-Seidel y Newton Raphson Desacoplado Rápido, no
concuerdan con los de la simulación porque apenas se realizó 2 iteraciones mientras que el
programa utilizo 99 iteraciones para su convergencia. Otro factor importante de no
concordancia fue que para el cálculo analítico no se tomó en cuenta las impedancias de los
generadores, mientras que para el programa desarrollado en ETAP si se tomó en cuenta.
Además, se pudo constatar que en el sistema eléctrico estudiando las barras están
operando fuera de los rangos establecidos (la norma establece +/- 5%), lo cual se pudiese
corregir agregando compensación reactiva (por medio del generador o por medio de banco
de capacitores); para que así el sistema opere de manera óptima.
En cuanto al estudio de cortocircuito se da el mismo escenario, para el cálculo
analítico no se tomaron en cuanto las impedancias transitorias y subtransitorias de los
generadores. Aunque en este estudio las corrientes de secuencia y de falla obtenidas
mediante el simulador y el cálculo analítico parecen semejantes; están en el orden de
diferencia de dos o tres por ciento (2% o 3%).
30
Referencias
1. Zurita, Sady. Análisis de Flujo de Carga. Universidad Nacional Experimental Politécnica
“Antonio José de Sucre” Vice-Rectorado Puerto Ordaz.
2. Díaz, Alexis. Análisis de Flujo de Carga. Universidad Nacional Experimental Politécnica
“Antonio José de Sucre” Vice-Rectorado Puerto Ordaz.
3. Grainger/Stevenson.2001.AnálisisdeSistemasdePotencia.México.McGrawHill.
31
