2022 Flujo de Carga Part1

Flujos de carga en Sistemas de Potencia
Prof. Juan Bautista R.
Asignatura: Análisis de Sistemas de Potencia I

Facultad de Ingeniería Eléctrica y Electrónica
Bibliografía
Prof. Juan Bautista R. (FIEE - UNI)

Introducción: ECUACIONES NODALES
Sea el sistema:

Su representación será:
Generador Transformador Transformador Generador
Transformador
Motor

Generador Transformador Transformador Generador
Transformador
Motor

Fuente de corriente en paralelo

con la admitancia



Transformación de las fuentes y de impedancias a admitancias:
0 Voltios

0 Voltios

0 Voltios

0 Voltios

0 Voltios

0 Voltios

Corrientes entrantes y corrientes salientes:
0 Voltios

Corrientes entrantes y corrientes salientes:
0 Voltios
Ecuación Nodal:
Corriente total saliente

Ecuaciones Nodales:
Esta última es linealmente dependiente de las otras tres ecuaciones

Ecuaciones Nodales:
V0=0 referencia V0=0
V0=0
V0=0

Agrupando:
En forma matricial:
Vector de Tensiones
Vector de Inyección de
respecto a la Barra de
Corrientes
Prof. Juan Bautista R. (FIEE - UNI) referencia
Agrupando convenientemente:
Agrupando convenientemente:
Número de barras
Tendremos entonces:
Tendremos entonces:
Y además:
Vector de Tensiones
Corrientes
referencia
Ejemplo: Plantear la Matriz de admitancias de barras, en el sistema:
Y=
34 Y=
43 j10.0
Matriz de admitancias de barras del sistema:
Vector de Tensiones
Corrientes
referencia
I = V
Vector de Tensiones
Corrientes
referencia
 I1  V1 
   
I2  = V2 
 I3  V3 
   
I4  V4 
Determinar las Tensiones V1 y V2
Determinar las Tensiones V1 y V2
V1
V2
Determinar la Matriz de Impedancias de Barra del sistema:
Z=
Z= Y =
En el problema (anterior), representar la red con fuentes de tensión,
sabiendo que las corrientes Inyectadas son:
Respuesta:
Hasta aquí hemos aprendido a plantear las
ecuaciones de corrientes inyectadas a cada
una de las barras.:
Y además:
Vector de Tensiones
Corrientes
referencia
Formulación del Problema de Flujo de carga

Se analizarán los sistemas de potencia
en régimen estacionario, equilibrados,
sinusoidales, sin anomalías.

Flujos de carga: importancia
• Permite determinar los flujos de potencia activa y

reactiva en una red eléctrica.
• Permite determinar las tensiones en las barras de una
red eléctrica (Sistema eléctrico de potencia - SdP).
• Permite calcular las pérdidas en el SdP.
•Permite estudiar las alternativas para la
planificación de nuevos sistemas o ampliación de
los ya existentes.
• Permite evaluar los efectos de pérdidas
temporales de generación o de circuitos de
transmisión.
• Permite evaluar los efectos de reconfigurar los
circuitos de un SdP (por ejemplo ante la pérdida
de una línea de transmisión).
• Permite evaluar las mejoras que se producen
ante el cambio en la sección de los conductores
de un SEP.
Sea la barra k
Generación
Carga

Inyección líquida de potencia en la barra k
Inyección líquida de
Potencia en la barra “k”
Generación
Carga
Formulación del Problema de Flujo de carga…
Inyección de Corriente en la barra “k”:

Inyección de
Potencia líquida en
la barra k

Inyección de
la barra k

Inyección de
la barra k

Inyección de
la barra k

De las ecuaciones nodales:


Entonces:

Entonces:

Ecuaciones de Flujo de Potencia
Las partes real e imaginaria de la potencia en la barra k son:
Para la potencia activa:

Sean:

Sean:

Sean:
Barra k Barra m

Sean:
Como:

Sean:
Entonces

Sean:
Entonces

Sean:
Entonces

Sean:
Entonces

Sean:
Entonces

Sea: θk − θm =
θ km
Entonces:

Para la Potencia Reactiva:
Y como:

Ejemplo: escribir las ecuaciones para el Sistema:

Como:

Como:
En forma compleja:

En forma polar:

En forma polar:

En forma polar:

En forma polar:

En forma polar:

En forma polar:

En forma polar:

Barra Slack (Barra oscilante, swing bus o Barra de referencia)
En los sistemas de transmisión las pérdidas no son conocidas antes del

flujo de potencia.
Es necesario mantener una barra en donde P no es especificada.
Si se especifica en esta barra: V como δ
δ es mantenido constante durante la solución del flujo de potencia.

δ es el ángulo de referencia del sistema.

Es barra de referencia, debido a que la potencia activa P y la potencia reactiva Q, no son
especificadas y “cubre” las pérdidas del sistema de transmisión.
La situación de la barra oscilante o flotante puede influir en la complejidad de los cálculos;

deberá siempre utilizarse la barra que se aproxima más a una barra de potencia infinita.
Normalmente a la barra oscilante se le asigna el número 1 como referencia, esta barra tiene
conectada la generación y es necesaria para llegar a la solución del flujo de carga.

Tipos de Barras
Barra de Tensión controlada,
Son el 95% Son el 5% del total de las
del total de barras
las barras
Slack
Datos
Por calcular

Los sistemas de potencia generan matrices Y altamente esparcidas.
Si una red tiene NB = 1000 barras y NR = 2000 ramas (lineas)
Número total de elementos de la matriz Y: NB 2 1000

= = 2
1'000, 000
Número total de elementos no nulos: NE =NB + 2 NR =1000 + 2 ( 2000 ) =5000
Número de elementos nulos: NN = NT − NE = 1'000, 000 − 5000 = 995, 000
Grado de esparcidad:  NN   995, 000 

=GE 100
=   100  =  99.5%
 NT   1'000, 000 

Sistemas de Ecuaciones del Flujo de Potencia

Una red eléctrica cualquiera tiene:
Una barra slack (datos V y θ) incógnitas P y Q
NPQ barras tipo PQ (P y Q datos Incógnitas V y θ)

NPV barras tipo PV (P y V datos incógnitas Q y θ)
Número total de barras de la red: NPQ + NPV + 1
Por tanto:
2(NPQ+NPV+1) son datos
2(NPQ+NPV+1) son incógnitas

Son dos tipos de incógnitas:
y Incógnitas asociadas al estado de la red (Variables de estado)
Que pueden ser obtenidas una vez conocidas las variables de estado
y
(utilizando las ecuaciones de potenciales de nodos).

En función de la existencia de dos tipos de incógnitas, el problema del flujo de carga (FdeC)
puede ser descompuesto en dos subsistemas de ecuaciones algebraicas:
SUBSISTEMA 1: Dimensión 2.NPQ + NPV

Determinación de las variables de estado desconocidas:




Resulta entonces un total de 2.NPQ+NPV incógnitas variables de estado.

En términos de potencia, son datos (especificados):
y datos
datos

y datos
datos
Para cada potencia dada, se puede escribir una ecuación de flujo de carga:

y datos
datos
Para barras PQ y PV

y datos
datos
Para barras PQ y PV
Para barras PQ

y datos
datos
Para barras PQ y PV
Para barras PQ
Resultando un sistema de 2.NPQ+NPV ecuaciones y el mismo

número de incógnitas = SISTEMA DETERMINADO

SUBSISTEMA 2: Dimensión NPV + 2
Determinación de las Potencias Nodales desconocidas.

Resuelto el sistema de ecuaciones del Subsistema 1, todas las tensiones de la red son
conocidas, y las incógnitas restantes son:
Para la barra SLACK 01 incógnita
Para las barras PV y la barra Slack: NPV+1 incógnitas
Por tanto se determinarán: NPV + 2 incógnitas

Ecuaciones del Subsistema 2:
Como el estado de la red es conocida, basta aplicar directamente las ecuaciones

de potencias nodales para las respectivas barras:
Para la Barra SLACK:

Ecuaciones del Subsistema 2:
Como el estado de la red es conocida, basta aplicar directamente las ecuaciones

de potencias nodales para las respectivas barras:
Para la Barra SLACK y Barras PV:

Las incógnitas del Subsistema 1 pueden ser escritas como:
θ Es el vector de ángulos de tensiones de las barras PQ y PV

V Es el vector de las magnitudes de tensiones de las Barras PQ

Las ecuaciones del FdeC para el Subsistema 1 pueden ser reescritas como:
Para las barras PQ y PV

Las ecuaciones del FdeC para el Subsistema 1 pueden ser reescritas como:
Para las barras PQ y PV
Para las barras PQ

Son los valores de las inyecciones de potencia activa y
reactiva especificados para las barras (considerados
constantes en principio) y dados por:

Son calculados a través de las ecuaciones de
potencias nodales.
Son llamados residuos u errores de

Potencia activa y reactiva

Las ecuaciones de FdeC en forma vectorial son:
P Vector de inyecciones de potencia activa en las barras PQ y PV
Q Vector de inyecciones de potencia reactiva en las barras PQ

Considere la función vectorial:

Como una solución del Subsistema 1 es obtenida cuando los errores son iguales
a Cero, las ecuaciones del Subsistema 1 pueden ser escritas en la forma:

Por ejemplo en el sistema:
Las variables de estado

desconocidas son:


desconocidas son:
 
  NPQ+NPV
  2+1
x= 
 
  NPQ
  2


desconocidas son:
θ 2 
θ  NPQ+NPV
 3 2+1
x = θ 4 
 
V2  NPQ
V4  2

El Subsistema 1 consiste en la resolución del siguiente sistema de ecuaciones:
 ∆P2 
 ∆P 
 3
g ( x) =
= ∆P4  0
  En el que:
∆
 2Q
 ∆Q4 

El Subsistema 2 consiste en la obtención de:


En el ejemplo:
Barra Slack
Barra Slack
Barra PV

En el ejemplo:
Barra Slack
Barra Slack
Barra PV
Como la red tiene 4 Barras, podemos escribir 8 ecuaciones de potenciales de nodos.
Es decir:
El sistema de ecuaciones del subsistema 1 = 5 de ellas.
para el Subsistema 2 = 3 restantes.
Total = 8 ecuacio. de potenciales de nodos.

Para el sistema adjunto, escribamos las ecuaciones de FdeC del subsistema 1

PQ
Slack

PQ
Slack
θ 2  NPQ + NPV = 1 + 0 = 1
x= 
Variables de estado que
corresponden al Subsistema
V2 
1
NPQ = 1

subsistema 1
PQ
Slack
θ 2  NPQ + NPV = 1 + 0 = 1
Variables de estado que corresponden al
x= 
V2 
Subsistema 1
NPQ = 1
Entonces las ecuaciones del Subsistema 1 son:

subsistema 1
PQ
Slack
θ 2  NPQ + NPV = 1 + 0 = 1
x= 
V2  NPQ = 1
Datos:

Slack PQ
Matriz de admitancias:

Slack PQ

Slack PQ

Slack PQ
Sustituyendo valores se obtiene el sistema completamente determinado:

Slack PQ
Sustituyendo valores se obtiene el sistema completamente determinado:
Resolviendo obtenemos los

valores de las variables de
estado incógnitas:

Ejemplo: Escribir las ecuaciones de FdeC para el sistema:

PQ
Slack V
Y14=-j2.0
PV
PQ

Datos de las 4 barras:
Datos de las ramas del sistema:

Considerando los datos:

Separando:
Y= G + jB
Y= + j

Variables de Estado Incógnitas:
 
 
  NPQ + NPV = 2 + 1 = 3
x= 
 
  NPQ = 2
 

Variables de Estado Incógnita:
θ 2 
θ 
 3 NPQ + NPV = 2 + 1 = 3
x = θ 4 
 
V2  NPQ = 2
V4 

Dimensión del subsistema 1: 2.NPQ + NPV = 2 x 2 + 1 = 5
( θ , θ3 , θ, V2, V4)
2 4
Ecuaciones:
Para barras PQ y PV

( θ , θ3 , θ, V2, V4)
2 4
Ecuaciones:
Para barras PQ y PV
Para barras PQ

( θ , θ3 , θ, V2, V4)
2 4
Ecuaciones:
Para barras PQ y PV
Para barras PQ
Barra PQ
Barra PQ
Barra PV
Barra PQ
Barra PQ

Las ecuaciones restantes constituyen el Subsistema 2, de dimensión NPV + 2 = 1
+2=3

Las ecuaciones restantes constituyen el Subsistema 2, de dimensión NPV + 2 = 1
+2=3
En el ejemplo:
01 para la Slack
NPV+1 = 1+1 = 2 = 01 para la

Slack y el resto para las barras
PV
Ya no son incógnitas (del subsistema 1): 2 , θ, θV ,θV

3 4 , 2 4

Reescribiendo las ecuaciones del subsistema 1,

Reescribiendo las ecuaciones del subsistema 1,

Colocando en forma de diferencias de potencias:
y considerando además que:
Entonces:

Sustituyendo los valores conocidos de la matriz de admitancias, sobre las ecuaciones del
subsistema 1, tenemos:
Las variables (incógnitas) del subsistema 1

corresponden al vector x:

Agrupando las ecuaciones del subsistema 1 y utilizando Pk(x) y Qk(x) para representar las
inyecciones de potencia activa y reactiva calculadas en función de variables x, entonces:
ecuaciones
ecuaciones
En forma compacta, las expresiones de las ecuaciones del subsistema 2 son:
barra Slack
   
barras Slack y PV  =  1, 3
   

Algunos conceptos deben quedar claros
Diferencia angular de las tensiones de barras
Sistema de Sistema de
potencia A potencia B
Sincronizando dos áreas
Sistema de
Sistema de
Potencia 1
Potencia 2
Sincronizando dos áreas
Area n° 1 Area n° 2
Estableciendo una segunda línea de enlace
Area n° 1 Area n° 2
TIPICO ACSR CONDUCTOR
Typical Transfer Limits
Capacitancia natural de la línea de transmisión
Modelo de la línea de transmisión
Proceso de apertura del
interruptor
Interruptor
Inductive and Capacitive Phase Angle Examples
Line-to-Line & Line-to-Ground Voltages
Power System Current Flow Path
Falla línea - línea
Voltage Phase Angle Difference
Model of Basic Governor Control System
Measurement of the Power Angle (u)
Illustration of Phase & Power Angles
Transmission Line PI Model
Use of the Active Power Transfer Equation
Use of the Active Power Transfer Equation
The Power-Angle Curve
The Power-Angle Curve
Evaluating Power Transfer Limits
Evaluating Power Transfer Limits
Total and Available Transfer Capability (TTC)
Simple Frequency Control System
Direction of Reactive Power Flow
A generator supplies reactive power to an inductive load. A generator
absorbs reactive power from a capacitive load
Direction of Reactive Power Flow
A generator supplies reactive power to an inductive load. A generator
absorbs reactive power from a capacitive load
Radial Power System
MW Required Versus Mvar Transferred
Illustration of Surge Impedance Loading
Note that this figure assumes the VS and
VR buses are similar strength.
Illustration of Surge Impedance Loading
Note that this figure assumes the VS and
VR buses are similar strength.
P-V Curve
P-V Curve
BENEFICIOS DE UNA ADECUADA ASIGNACION DE COMPENSACIÓN
REACTIVA
•Reducción de pérdidas de transmisión.

•Mejora del perfil de tensiones.
•Mejorar la utilización de la capacidad de transmisión.
•Postergar inversiones por necesidades de refuerzos.
•Aumento de los márgenes de reserva de potencia reactiva frente a eventuales
perturbaciones.
•Minimizar costos de operación del sistema; evitar restricciones al despacho
económico; prevenir colapsos de tensión y mantener márgenes de seguridad.
•Mejora del margen de estabilidad estacionaria y transitoria del sistema.
INFLUYEN EN EL BALANCE DE POTENCIA REACTIVA
•La variación de la tensión de generación.

•La variación de los taps en los transformadores.
•Conexión y desconexión de bancos de capacitores y
reactores.
•Configuración de los compensadores estáticos (SVC)
•La configuración de la red y el despacho de generación.
INTRODUCCION A LAS RESTRICCIONES EN LOS SISTEMAS ELÉCTRICOS DE POTENCIA
Demanda del sistema eléctrico
restricciones de igualdad se encuentran en las ecuaciones de flujo de carga.

Restricciones de Capabilidad del generador (polos lisos)
Estator: +Q Límite de corriente

de campo
S = 3 Vnom I nom
S = 3 Vnom I nom Límite de corriente
Rotor: en el estator
VE VE Límite mínimo de la
P= sen δ fuente de energía
Xs Xs r mecánica
V
Q= [E cos δ − V ] m MVAnom
Xs fp = 0.9 +
Recordando: φ P
0 1 pu
P2 + Q2 = S 2
2
 V2 VE 2 Límite máximo de la
2
P + Q +  = ( ) δ m′
j fuente de energía
 X s X s
V2 mecánica
( x − a ) 2 + ( y − b) 2 = r 2 (0,− )
Xs
Entonces: a=0 − 0.6 Q Límite de
i calentamiento de
2
V e cabezales o de
b=− h subexcitación
Xs
−Q Límite práctico de
estabilidad
Restricciones de capabilidad del generador (polos salientes)
Estator: Ángulo δ :
S = 3 Vnom I nom I X q cos(φ )

δ = tan −1
V + I X q sen(φ )
Rotor:
Voltaje interno:
VE V2 1 1
P= sen δ + ( − ) sen(2δ )
Xd 2 Xq Xd
E = V cos(δ ) + I X d sen(δ + φ )
VE V Xd − Xq
2
V 2
1 1
Q= cos δ + ( ) cos(2δ ) − ( + )
Xd 2 Xd Xq 2 Xd Xq
Ángulo máximo δ m:
Recordando:  C   C  2 1  
δ m = cos  −
−1
+   + 
S = P + jQ 8B   8B  2  
  
V2 1 1 V 2 Xd − Xq VE
S =−j ( + )+ ( )( sen(2δ ) + j cos(2δ )) + ( senδ + j cos δ )
2 Xd Xq 2 Xd Xq Xd
A B C
S = A+ B +C
Restricciones de capabilidad del generador (polos salientes)
+Q
S = A+ B +C S = 3 Vnom I nom Límite de corriente
de campo
V2 1 1 Límite de corriente
A=−j ( + ) en el estator
2 Xd Xq
MVAnom
m
V 2 Xd − Xq
B= ( )(sen(2δ ) + j cos(2δ ))
2 Xd Xq 1 pu P
0
VE
C= (sen δ + j cos δ ) C
Xd
Límite practico de estabilidad

(margen de 10%) δ
10% Límite práctico de

2δ estabilidad
Circulo de
reluctancia B permanente
A
δm
2δ m Límite teórico de
estabilidad
−Q
permanente
Por tanto la potencia de un generador no
puede estar fuera del rango de potencias
mínima y máxima
RESTRICCIONES DE TENSION
En una red de transmisión es muy importante conservar las magnitudes de los voltajes en los diferentes nodos,
los valores de voltaje y ángulo de fase pueden variar, pero solo dentro de límites establecidos en donde no afecte
al sistema ni a los equipos conectados.
Por otro lado sería incosteable colocar reguladores de voltaje en cada lugar donde existiera un cambio en las
magnitudes, por lo tanto existe un rango de voltaje de operación donde operan eficientemente los instrumento
conectados, el cual es fácil de obtener esta restricción es utilizada en flujos óptimos.
RESTRICCIONES DE RESERVA RODANTE.
Estas restricciones son requeridas para enfrentar a contingencias como la salida de alguna de las unidades,
o variaciones de carga inesperada en el sistema. Es muy común que la reserva rodante de los generadores
en un sistema sea suficiente para suplir la salida de la unidad más grande, con esto se garantiza el
suministro de energía ante esta contingencia. La total generación debe ser mayor que la demanda más las
pérdidas, más un cierto valor preestablecido de reserva rodante.
RESTRICCIONES DE LÍNEAS DE TRANSMISIÓN.
El flujo de potencia en una línea de transmisión está restringido principalmente por cuestiones de
seguridad de estabilidad y de intercambio preestablecido, en este tipo de restricciones intervienen las
características de la red y su capacidad:
Controles en un sistema de
potencia
Estructura de un sistema de
potencia
Estructura de un sistema de
potencia
Método de Newton Raphson

para la solución de las ecuaciones
de flujo de carga
Bibliografía :

Características del método:
1. Velocidad de convergencia ‘cuadrática’ (el
número de cifras significativas se duplica luego de
cada iteración)
2. Confiable, no sensible a la elección de la barra
flotante.
3. Solución precisa obtenida luego de 4-6
iteraciones.
4. J debe ser re-calculada e invertida luego de
cada iteración. (J es una matriz esparsa, tiene
estructura simétrica, pero los valores no son
simétricos)
Las ecuaciones básicas del Sistema S1 que serán resueltas son:
Errores (power mismatches):

Sistema a ser resuelto por el método Newton Raphson
Siendo n el número de barras tenemos:

 ∆P1   ∆θ1 
 ∆P   ∆θ 
 2   2 
 ...   ...  NPQ + NPV

   
 ...   ... 
Ecuación
Matricial a
 ∆Pn −1   ∆θ n −1 
resolver:   = −J  
 ∆Q1   ∆V1 
 ∆Q   ∆V 
 2
  2

 ...   ...  NPQ
 ...   ... 
    Que tiene la
forma:
 ∆Ql   ∆Vl 
[b] = [ A][ x ]
(i ) (i )
 ∆P1   ∆θ1 
 ∆P   ∆θ 
 2   2 
 ...   ... 
   
 ...   ... 
Para el proceso
iterativo:  ∆Pn −1  (i )
 ∆θ n −1 
  = −J  
 ∆Q1   ∆V1 
 ∆Q   ∆V 
 2
  2

 ...   ... 
 ...   ... 
   
 ∆Ql   ∆Vl 
El Jacobiano (Matriz J) General
Sea:  f1 
f 
F =  2
:
 
Y las variables:
 fn 
El Jacobiano es:

Jacobiano aplicado a la solución del Flujo de Potencia
No. Barras PQ

Resumen del Método de Newton Raphson
Ecuaciones del Subsistema 1 (S1)
Sistema que será resuelto:

Sistema matricial
Que tiene la forma:

Entonces las variables son:

J=

Orden de la matriz Jacobiana

Algoritmo:
1.- Escribir
2.- Establecer las condiciones

iniciales de las variables de
estado:
Y hacer: i=0
3.- Calcular:
Y verificar
convergencia, si:
Y si:
Parar.
4.- Hacer:
Escribir el
Jacobiano:
5.- Resolver el sistema

linealizado:

6.- Actualizar la solución del
problema
7.- Volver al paso 3.-

Por ejemplo, en el sistema:
PV
Slack
PQ

Por ejemplo, en el sistema:
PV
Slack
θ 2 

x = θ3  
PQ
V3 
Ecuación del subsistema 1
1+1
θ 2 

x = θ3  
V3 

Jacobiano de dimensión n - 1 + 1 = 3 -1 +1
1+1
1+1 1

Elementos de las submatrices del Jacobiano: H, N, M, y L
kk se refiere a un término de la diagonal

km se refiere a un término fuera de la
diagonal. Línea k columna n
kk se refiere a un término de la diagonal
km se refiere a un término fuera de la
diagonal. Línea k columna n
Estructura del Jacobiano
1.- Los elementos fuera de la diagonal principal correspondientes a las barras no directamente
conectadas son nulos. Es decir, el Jacobiano es una matriz altamente esparsa.
Si las barras k y m no están directamente conectadas, entonces:
y por tanto:

2.- Las matrices H, N, M y L tienen una estructura semejante a la Matriz YBARRA, excepto por las lineas
y columnas no representadas.
Si todas las Barras fueran tipo PQ , la estructura del jacobiano será semejante a la
estrucutura de YBARRA y las submatrices H, N, M y L serían cuadradas.
Las matrices H, N, M y L son simétricas en estructura. Si existe H12 entonces existe H21

3.- El Jacobiano es simétrico en valores, asi como las matrices H, N, M y L
Además que:

Por ejemplo, escribamos el Jacobiano del Sistema:

2+1
2+1 2

Escribir el JACOBIANO del SdeP:

Slack PQ
PQ
PV

Slack PQ
PQ
PV
− −
− − NPV+NPQ=
    1+2=3
− = − J −
   
 − − NPQ = 2
 −   − 
Slack PQ
PQ
PV
−  ∆δ 2 
−  ∆δ  NPV+NPQ=
   3 1+2=3
 −  = − J  ∆δ 4 
   
 − ∆
 2V
NPQ = 2
 −   ∆V4 
Slack PQ
PQ
PV
 ∆P2   ∆δ 2 
 ∆P   ∆δ  NPV+NPQ=
 3   3 1+2=3
 ∆P4  = − J  ∆δ 4 
   
∆
 2Q ∆
 2V
NPQ = 2
 ∆Q4   ∆V4 
JACOBIANO = J
NPV+NPQ=
1+2=3
NPQ = 2
NPV+NPQ= 1+2=3 NPQ = 2

JACOBIANO = J
NPV+NPQ=
1+2=3
NPQ = 2
NPV+NPQ= 1+2=3 NPQ = 2

1- 3 y 2-4 no conectados

Escribir el Jacobiano del Sistema:

PQ PV
Slack
PQ PV

PQ PV
Slack
NPV = 2
NPQ = 2
PQ PV

PQ PV
Slack
NPV = 2
NPQ = 2
PQ PV
−
−
  NPV+NPQ variables d
−
x= 
−
−
  NPQ variables V
−
PQ PV
Slack
NPV = 2
NPQ = 2
PQ PV
δ 2 
δ 
 3 NPV+NPQ variables d
δ 4 
x= 
δ 5 
V2 
NPQ variables V
 
V4 
PQ PV
Slack
NPV = 2
NPQ = 2
PQ PV
δ 2 
δ 
δ 4 
El JACOBIANO es de 6x6 x= 
δ 5 
V2 
NPQ variables V
 
V4 
PQ PV
Slack
NPV = 2
NPQ = 2
PQ PV
δ 2 
δ 
δ 4 
El JACOBIANO es de 6x6 x= 
Barras no conectadas: δ 5 
1-2, 1-3, 1-5, 2-5, 3-4
V2 
NPQ variables V
 
V4 
Ecuación teórica a  ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 
resolver:  ∂δ ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
 2 
 ∂P3 ∂P3 ∂P3 ∂P3 ∂P3 ∂P3 
 ∆P2   ∂δ ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4   ∆δ 2 
 ∆P   2  
 3   ∂P4 ∂P4 ∂P4 ∂P4 ∂P4 
∂P4  3 ∆δ
 
 ∆P4   ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4   ∆δ 4 
  = −   
 ∆P5  ∂P5 ∂P5 ∂P5 ∂P5 ∂P5 ∂P5  ∆δ 5 
 
 ∆Q2   ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4  ∆V2 
   ∂Q2   
 ∆Q4  ∂Q2 ∂Q2 ∂Q2 ∂Q2 ∂Q2  ∆V4 
 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
 ∂Q ∂Q4 ∂Q4 ∂Q4 ∂Q4 ∂Q4 
 4 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 

Ecuación teórica a  ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 
resolver, los  ∂δ ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
marcados en rojo
deberán ser CEROS
 2 
 ∂P3 ∂P3 ∂P3 ∂P3 ∂P3 ∂P3 
 ∆P2   ∂δ ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4   ∆δ 2 
 ∆P   2  
 3   ∂P4 ∂P4 ∂P4 ∂P4 ∂P4 
∂P4  3 ∆δ
 
 ∆P4   ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4   ∆δ 4 
  = −   
 ∆P5  ∂P5 ∂P5 ∂P5 ∂P5 ∂P5 ∂P5  ∆δ 5 
 
 ∆Q2   ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4  ∆V2 
   ∂Q2   
 ∆Q4  ∂Q2 ∂Q2 ∂Q2 ∂Q2 ∂Q2  ∆V4 
 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
Barras no conectadas:  ∂Q ∂Q4 ∂Q4 ∂Q4 ∂Q4 ∂Q4 
1-2, 1-3, 1-5, 2-5, 3-4  4 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂δ 5 ∂V2 ∂V4 

 ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 
 ∂δ 0
∂δ 3 ∂δ 4 ∂V2 ∂V4 
 2 
 ∂P3 ∂P3 ∂P3 ∂P3 
 ∂δ 0 0 
 ∆P2  ∂δ 3 ∂δ 5 ∂V2 ∆δ
 ∆P   2   2 
 3   ∂P4 ∂P4 ∂P4 ∂P4 ∂P4   ∆δ 3 
 0 
 ∆P4   ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4   ∆δ 4 
  = −   
 ∆P5  ∂P5 ∂P5 ∂P5 ∂P5  ∆δ 5 
 0 0 
 ∆Q2   ∂δ 3 ∂δ 4 ∂δ 5 ∂V4  ∆V2 
   ∂Q2   
 ∆Q4  ∂Q2 ∂Q2 ∂Q2 ∂Q2  ∆V4 
 0 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂V2 ∂V4 
Barras no conectadas:  ∂Q ∂Q4 ∂Q4 ∂Q4 ∂Q4 
1-2, 1-3, 1-5, 2-5, 3-4  4 0 
 ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4 

El sistema analizado fue: Eliminemos ahora la línea: 2-3
Barras no conectadas: Barras no conectadas:

1-2, 1-3, 1-5, 2-5, 3-4 1-2, 1-3, 1-5, 2-5, 3-4 y 2-3

Barras no conectadas:
1-2, 1-3, 1-5, 2-5, 3-4 y 2-3
 ∂P2 ∂P2 ∂P2 ∂P2 ∂P2 
 ∂δ 0
Los marcados de
∂δ 3 ∂δ 4 ∂V2 ∂V4 
MAGENTA deberán ser  2 
CEROS también  ∂P3 ∂P3 ∂P3 ∂P3 
 0 0 
∆
 2P ∂δ 2 ∂δ 3 ∂δ 5 ∂V2  ∆δ 2 
 ∆P    
 3  ∂P4 ∂P4 ∂P4 ∂P4 ∂P4   ∆δ 3 
 0 
 ∆P4   ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4  ∆δ 4 
  = −   
 ∆P5  ∂P5 ∂P5 ∂P5 ∂P5  ∆δ 5 
 0 0 
 ∆Q2   ∂δ 3 ∂δ 4 ∂δ 5 ∂V4  ∆V2 
   ∂Q2   
 ∆Q4  ∂Q2 ∂Q2 ∂Q2 ∂Q2  ∆V4 
 0 
 ∂δ 2 ∂δ 3 ∂δ 4 ∂V2 ∂V4 
 ∂Q ∂Q4 ∂Q4 ∂Q4 ∂Q4 
 4 0 
 ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
Barras no conectadas:
1-2, 1-3, 1-5, 2-5, 3-4 y 2-3
 ∂P2 ∂P2 ∂P2 ∂P2 
 ∂δ 0 0
∂δ 4 ∂V2 ∂V4 
 2 
 ∂P3 ∂P3 
 0 0 0 0 
 ∆P2  ∂δ 3 ∂δ 5  ∆δ 2 
 ∆P    
 3  ∂P4 ∂P4 ∂P4 ∂P4 ∂P4   ∆δ 3 
 0 
 ∆P4   ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4  ∆δ 4 
  = −   
 ∆P5  ∂P5 ∂P5 ∂P5 ∂P5  ∆δ 5 
 0 0 
 ∆Q2   ∂δ 3 ∂δ 4 ∂δ 5 ∂V4  ∆V2 
   ∂Q2   
 ∆Q4  ∂Q2 ∂Q2 ∂Q2  ∆V4 
 0 0 
 ∂δ 2 ∂δ 4 ∂V2 ∂V4 
 ∂Q ∂Q4 ∂Q4 ∂Q4 ∂Q4 
 4 0 
 ∂δ 2 ∂δ 4 ∂δ 5 ∂V2 ∂V4 
Por ejemplo:
Datos de Barras:
Datos de Líneas:

1.- Determinamos YBARRA
Donde:

2.- Prueba de Convergencia en relación a las condiciones iniciales i = 0
No converge, el proceso comienza.

3.- Proceso iterativo de Newton Raphson,
Primera iteración: i = i + 1 = 1

No converge. El proceso continúa

Segunda iteración: i = i + 1 = 2

Converge

Solución encontrada para todas las variables de estado o sea V y : θ
4.- Solución del Subsistema 2, sustitución de variables
Expandir:

Sustituyendo los valores fijos:
Sustituyendo los valores encontrados en el proceso iterativo:

Expandiendo:

Expandiendo:

Solución del Flujo de carga:

Ejemplo Problema expuesto en el libro de SERGIO HAFFNER (Brasil)
Teniendo en cuenta las tolerancias en potencias activa y reactiva respectivamente: P

θ θ
=0.001 y Q =0.001 , realizar los cálculos de FdeC del sistema:

Variables conocidas para el subsistema 1

Matriz de admitancias

Ecuaciones del subsistema 1
Observar que son ecuaciones no lineales.

Resolviendo el sistema:
Siendo:

Jacobiano
Como:

Como la solución inicial es:
Entonces los cálculos son:
Iteración

Para la tolerancia
Tendremos:
Las ecuaciones del subsistema 2 son:
Sustituyendo los valores ya

conocidos:

Ejemplo:
En el sistema eléctrico de potencia que se muestra:
a) Determinar la matriz de admitancias Y del sistema.

b) El sistema de ecuaciones del Subsistema 1 que al resolver, obtendremos las incógnitas
variables de estado.
c) Determinar por el Método de NEWTON RAPHSON la Tensión E2 (Módulo y ángulo) para dos
iteraciones (llenar la tabla adjunta sustentando con los cálculos de cómo se obtienen los
resultados).

Matriz de admitancias
Ecuaciones a ser utilizadas
Entonces las ecuaciones de flujo de carga son:

Solución:

Analizar el sistema DEL PROBLEMA 6.9 del SARMA:

Datos de Barras
Datos de Líneas

Datos de Transformadores
Solución: Datos de entradas e incógnitas

De los cálculos en pu se obtiene:
0.0015 + j 0.02 =Z 0.00225 + j 0.025 0.00075 + j 0.01

Y = j 0.44
Carga
0.0045 + j 0.05 0.009 + j 0.10

j 0.88 j1.72
Carga

0.0015 + j 0.02 =Z 0.00225 + j 0.025 0.00075 + j 0.01

Y = j 0.44
Carga
PQ PQ
Slack 0.009 + j 0.10 PV

0.0045 + j 0.05
j 0.88 j1.72
PQ
Carga

0.0015 + j 0.02 =Z 0.00225 + j 0.025 0.00075 + j 0.01

Y = j 0.44
Carga
PQ PQ
Slack 0.009 + j 0.10 PV

0.0045 + j 0.05
j 0.88 j1.72
θ 2 
θ 
PQ  3
θ 4 
Variables de  
Estado x = θ5 
Carga
V2 
 
V4 
V 
 5
0.0015 + j 0.02 =Z 0.00225 + j 0.025 0.00075 + j 0.01

Y = j 0.44
Carga
PQ PQ
Slack 0.009 + j 0.10 PV

0.0045 + j 0.05
j 0.88 j1.72
 ∆P2   ∆θ 2 
 ∆P   ∆θ 
 3   3
PQ Ecuación a  ∆P4   ∆θ 4 
resolver:    
∆
 5P =− J ∆
 5θ
 ∆Q2   ∆V2 
Carga    
∆Q
 4 ∆V
 4
 ∆Q   ∆V 
 5  5
Formación de la Matriz de Admitancias
pu
pu
pu
pu
pu

La Matriz de admitancias será:
 1 1 
 0.0015 + j 0.02 0 0 0 − 
0015 + j 0.02
 
 1 1 1 1 
 0 + j 0.88 + + j1.72 0 − − j1.72 − − j 0.88 
0.0045 + j 0.05 0.009 + j 0.10 0.009 + j 0.10 0.0045 + j 0.05
 
 1 1 
Y = 0 0 − 0 
0.00075 + j 0.01 0.00075 + j 0.001
 
 1 1 1 1 1 
 0 − − j1.72 − + + j 0.44 − − j 0.44 
0.009 + j 0.10 0.00075 + j 0.001 0.00075 + j 0.02 0.00225 + j 0.025 0.00225 + j 0.025
 
 1 1 1 1 1 1 
 − 0015 + j 0.02 − − j 0.88 0 − − j 0.44 + + j 0.44 + + j 0.88
 0.0045 + j 0.05 0.00225 + j 0.025 0.0015 + j 0.02 0.00225 + j 0.025 0.0045 + j 0.05 
Haciendo operaciones obtenemos:
 3.73 − j 49.72 0 0 0 −3.73 + j 49.72 

 0 2.68 − j 28.46 0 −0.89 + j9.92 −1.79 + j19.84 

Y = 0 0 7.46 − j99.44 −7.46 + j99.44 0 
 
 0 −0.89 + j9.92 −7.46 + j99.44 11.92 − j147.96 −3.57 + j39.68
 −3.73 + j 49.72 −1.79 + j19.84 0 −3.57 + j39.68 9.09 − j108.58 

 3.73 − j 49.72 0 0 0 −3.73 + j 49.72 
 0 2.68 − j 28.46 0 −0.89 + j9.92 −1.79 + j19.84 

Y = 0 0 7.46 − j99.44 −7.46 + j99.44 0 
 
 0 −0.89 + j9.92 −7.46 + j99.44 11.92 − j147.96 −3.57 + j39.68
 −3.73 + j 49.72 −1.79 + j19.84 0 −3.57 + j39.68 9.09 − j108.58 
El software muestra:
El software muestra la MATRIZ DE ADMITANCIAS:

Ojo: El resultado es solo una “fotografía” de un instante
El resultado del FdeC será: de operación del sistema.
Problema analizado en el
libro

Ejemplo
Encontrar las magnitudes de ángulos de fase y los flujos en las líneas para el SdeP que se
muestra.
Parámetros en pu

Slack PV
NPQ =1
PQ
NPV =1
P2=0.5
Slack Z=0.02+j0.3 PV
V=1.02
B=j0.15 V2=1.0
Z=0.01+j0.1
B=j0.1 Z=0.01+j0.1
B=j0.1
NPQ =1
PQ
NPV =1
1.2+j0.5
P2=0.5
V=1.02
B=j0.15 V2=1.0
Z=0.01+j0.1
B=j0.1 Z=0.01+j0.1
B=j0.1
NPQ =1
PQ
NPV =1
1.2+j0.5
−
Variables de   NPQ + NPV = 2 incógnitas δ
estado: x = −
 −  NPQ = 1 incógnita V
P2=0.5
V=1.02
B=j0.15 V2=1.0
Z=0.01+j0.1
B=j0.1 Z=0.01+j0.1
B=j0.1
NPQ =1
PQ
NPV =1
1.2+j0.5
δ 2 
Variables de   NPQ + NPV = 2 incógnitas δ
estado: x = δ 3 
V3  NPQ = 1 incógnita V
No olvidar las variables de estado: d2, d3 y V3

Sustituyendo las cantidades conocidas: V1=1.02, V2=1.0 y d1=0

Para aplicar las iteraciones NEWTON RAPHSON tendremos:

Donde:

En general tener en cuenta que:

En general tener en cuenta que:
Para iniciar los cálculos pueden tomarse: Vi = 1.0 y di=0

Iteración 1 Evaluando el Jacobiano y la Matriz de errores: Ojo:
Vi = 1.0 y di=0

Iteración 1 Evaluando el Jacobiano y la Matriz de errores:



Error = El mayor valor

absoluto.

Resolviendo la ecuación



Resolviendo la ecuación Matriz incógnita

Resolviendo la ecuación Matriz incógnita
Resolviendo:

Entonces:
Radianes
Error

Iteración 2
Evaluando el Jacobiano con: y

Iteración 2
En general los elementos de la diagonal son mayores en

magnitud (son sumas de términos).

Iteración 2
Error = El mayor valor absoluto.

Resolviendo
Y además
Con el error

Iteración 3
Evaluando el Jacobiano con los valores: y
Error = El mayor valor absoluto.

Resolviendo obtenemos
Y entonces
Suficientemente pequeño:
Con el nuevo error
HAY CONVERGENCIA en tres
iteraciones.

La última tarea en el cálculo del FdeC es determinar las potencias reactivas generadas, la
potencia activa en la barra SLACK y los flujos en las líneas. Las potencias generadas pueden ser
calculadas directamente de las ecuaciones de flujo de carga:
Por tanto:

El total de la potencia de pérdidas en el sistema es:
Las pérdidas en la línea i - j son calculadas en el envío y recepción de las líneas.

Las potencia aparente en la línea de i a j es:

Las potencia aparente en la línea de j a i es:
En forma similar hay que calcular las potencia Pji y Qji

Ejemplo:
Y=
PQ
PV
NPQ =1
NPV =1
−
x =  − 
Variables de NPQ + NPV = 2 incógnitas
estado:
 −  NPQ = 1 incógnita
PQ
PV
NPQ =1
NPV =1
δ 2 
Variables de
  NPQ + NPV = 2 incógnitas d
estado: x = δ 3 
V3  NPQ = 1 incógnita V
Resumen del Método de Newton Raphson
Ecuaciones del Subsistema 1 (S1)
Sistema que será resuelto:

Y=
Y=
Y=
Sistema que será resuelto
Y las variables
Valores iniciales
Valores fijos
Sustituyendo ahora susceptancias (B) y conductancias (G)
Y=
Es decir
Elementos de la Matriz JACOBIANA
Elementos de la Matriz JACOBIANA
Simplificando
Continúa
Y=
Para aplicar las iteraciones NEWTON RAPHSON tendremos:
Ejemplo
Resolver el problema de flujo de carga usando el método de NR:
V=1, θ=0 j0.1

1.5+j0.8
2
1
j0.2 3 j0.25
P=1, V=1
Método de Newton-Raphson
Ejemplo
Barra 1: Flotante (V1 y θ1 dados) V=1, θ=0 j0.1

1.5+j0.8
1 2
Barra 2: Barra PQ j0.2 3 j0.25
(V2 y θ2 desconocidos) P=1, V=1
2 ecuaciones - balance de
potencia activa y reactiva.
Barra 3: Barra PV - θ3 desconocido

(V3 especificado)
1 ecuación: balance de
potencia activa.
Ejemplo
− j15 j10 j5 
Y = G + jB =  j10 − j14 j 4 
 j 5 j 4 − j 9
θ 2 θ 3 V2
P2  H 22 H 23 N 22 
J=
P3 H
 32 H 33 N 32 
Q2  M 22 M 23 L22 
Ejemplo
k =n
P2 = ∑ V2 Vk b2 k sen θ 2 k = V2 [10V1 sen θ 2 + 4V3 sen(θ 2 − θ 3 )]
k =1
k =n
P3 = ∑ V3 Vk b3k sen θ 3k = V3 [5V1 sen θ 3 + 4V2 sen(θ 3 − θ 2 )]
k =1
k =n
Q2 = −∑ V2 Vk b2 k cosθ 2 k = 14V22 − V2 [10V1 cosθ 2 + 4V3 cos(θ 2 − θ 3 )]
k =1
Ejemplo
V10 = 1, V20 = 1, V30 = 1,θ10 = 0,θ 20 = 0,θ 30 = 0
P2 = V2 [10V1 sen θ 2 + 4V3 sen(θ 2 − θ 3 )] = 1[10 ⋅1sen 0 + 4 ⋅1sen 0] = 0
P3 = V3 [5V1 sen θ 3 + 4V2 sen(θ 3 − θ 2 )] = 1[5 ⋅1sen 0 + 4 ⋅1sen 0] = 0
Q2 = 14V22 − V2 [10V1 cosθ 2 + 4V3 cos(θ 2 − θ 3 )] =

14 ⋅1 − 1[10 ⋅1 cos 0 + 4 ⋅1 cos 0] = 0
Ejemplo
k =n
∆Pi = Pi − ∑ Vi Vk ( g ik cosθ ik + bik sen θ ik )
sp
k =1
k =n
∆Qi = Q − ∑ Vi Vk ( g ik sen θ ik − bik cosθ ik )
i
sp
k =1
 ∆P2   − 1.5 − 0   − 1.5 

 ∆P  =  1.0 − 0  =  1.0 
 3    
∆Q2  − 0.8 − 0 − 0.8
Ejemplo
θ2 θ3 V2
P2  − Q2 + 14V22 − 4V2V3 cos(θ 2 − θ 3 ) P2 
J=  2 
P3  − 4V V
2 3 cos(θ 2 − θ 3 ) − Q 3 + 9V 3 − 4V V
2 3 sen(θ 2 − θ )
3 
Q2  P − 4V V sen( θ − θ ) Q + 14V 2 
 2 2 3 2 3 2 2 
θ 2 ........ θ 3 ........ V2
P2  14.000 − 4.000 0.000 
J=  
P3  − 4 .000 9 .000 0 .000 
Q2  0.000 0.000 14.000
Ejemplo
0.0818 0.0364 0.0000
J −1 = 0.0364 0.1273 0.0000
0.0000 0.0000 0.0714
 ∆θ 2  0.0818 0.0364 0.0000  − 1.5 

 ∆θ  = 0.0364 0.1273 0.0000 ⋅  0 
 3     
∆V2 / V2  0.0000 0.0000 0.0714 − 0.8
 ∆θ 2  − 0.0864
 ∆θ  =  0.0727 
 3   
∆V2 / V2   − 0.0571
Ejemplo
θ = θ + ∆θ 2 = 0 − 0.0864 = −0.0864
1
2
0
2
θ 31 = θ 30 + ∆θ 3 = 0 + 0.0727 = 0.0727
1 ∆V2
0 0
V = V +V
2 2 = 1 − 1 ⋅ 0.0571 = 0.9429
2
V2
Esto completa la primer iteración.
Ahora re-calculamos las potencias de la barra con los nuevos
valores de las variables de estado:
Ejemplo
V11 = 1, V21 = 0.9429, V31 = 1,θ11 = 0,θ 21 = −0.0864,θ 31 = 0.0727
P2 = V2 [10V1 sen θ 2 + 4V3 sen(θ 2 − θ 3 )] = −1.4107

P3 = V3 [5V1 sen θ 3 + 4V2 sen(θ 3 − θ 2 )] = 0.9608
Q2 = 14V22 − V2 [10V1 cosθ 2 + 4V3 cos(θ 2 − θ 3 )] = −0.6715
 ∆P2   − 1.5 + 1.4107   − 0.0893

 ∆P  =  1.0 − 0.9608  =  0.0392 
 3    
∆Q2  − 0.8 − 0.6715 − 0.1285
Ejemplo
θ2 θ3 V2
P2  − Q2 + 14V22 − 4V2V3 cos(θ 2 − θ 3 ) P2 
J=  2 
P3  − 4V V
2 3 cos(θ 2 − θ 3 ) − Q 3 + 9V 3 − 4V V
2 3 sen(θ 2 − θ )
3 
Q2  P − 4V V sen( θ − θ ) Q + 14V 2 
 2 2 3 2 3 2 2 
θ2 θ3 V2
P2  13.1172 − 3.7238 − 1.4107 
J=  
P3  − 3. 7238 8. 7106 0. 5975 
Q2  − 1.4107 0.5975 11.7742 
Ejemplo
0.0876 0.0369 0.0086 
J −1 = 0.0369 0.13707 − 0.0022
0.0086 − 0.0022 0.0861 
 ∆θ 2  0.0876 0.0369 0.0086   − 0.0893

 ∆θ  = 0.0369 0.13707 − 0.0022 ⋅  0.0392 
 3     
∆V2 / V2  0.0086 − 0.0022 0.0861  − 0.1285
 ∆θ 2   − 0.075 
 ∆θ  =  0.021 
 3   
∆V2 / V2  − 0.0119
Ejemplo
θ = θ + ∆θ 2 = −0.0864 − 0.0075 = −0.09385

2
2
1
2
θ 32 = θ 31 + ∆θ 3 = 0.0727 + 0.0021 = 0.07485

2 1 ∆V2 1
V = V +V
2 2 = 0.9429 − 0.0119 ⋅ 0.9429 = 0.9316
2
V2
Esto completa la segunda iteración.
Ahora re-calculamos las potencias de la barra con los nuevos
valores de las variables de estado:
Ejemplo
V12 = 1, V22 = 0.9316, V32 = 1,θ12 = 0,θ 22 = −0.09385,θ 32 = 0.07485
P2 = V2 [10V1 sen θ 2 + 4V3 sen(θ 2 − θ 3 )] = −1.4987

P3 = V3 [5V1 sen θ 3 + 4V2 sen(θ 3 − θ 2 )] = 0.9995
Q2 = 14V22 − V2 [10V1 cosθ 2 + 4V3 cos(θ 2 − θ 3 )] = −0.7979
 ∆P2   − 1.5 + 1.4987  − 0.0013

 ∆P  =  1.0 − 0.9995  =  0.0005 
 3    
∆Q2  − 0.8 − 0.7979  − 0.0021
Ejemplo
θ2 θ3 V2
P2  − Q2 + 14V22 − 4V2V3 cos(θ 2 − θ 3 ) P2 
J=  2 
P3  − 4V V
2 3 cos(θ 2 − θ 3 ) − Q 3 + 9V 3 − 4V V
2 3 sen(θ 2 − θ )
3 
Q2  P − 4V V sen( θ − θ ) Q + 14V 2 
 2 2 3 2 3 2 2 
θ2 θ3 V2
P2  12.9488 − 3.7736 − 1.4987 
J=  
P3  − 3 .6736 8 .6596 0 .6257 
Q2  − 1.4987 0.6257 11.3529 
Ejemplo
0.0888 0.0370 0.0097 

J −1 = 0.0370 0.1313 − 0.0024
0.0097 − 0.0024 0.0895 
 ∆θ 2  0.0888 0.0370 0.0097   − 0.0893

 ∆θ  = 0.0370 0.1313 − 0.0024 ⋅  0.0392 
 3     
∆V2 / V2  0.0097 − 0.0024 0.0895  − 0.1285
 ∆θ 2  − 0.00012
 ∆θ  =  0.00002 
 3   
∆V2 / V2  − 0.00020
Ejemplo
θ = θ + ∆θ 2 = −0.09385 − 0.00012 = −0.09397

3
2
2
2
θ 33 = θ 32 + ∆θ 3 = 0.07485 + 0.00002 = 0.7486

3 2∆V2 2
V = V +V
2 2 = 0.9316 − 0.0002 ⋅ 0.9316 = 0.9314
2
V2
Esto completa la tercera iteración.
El método ha convergido ya que el vector de apartamiento es
casi cero.
En el sistema:

Jacobiano

Jacobiano

Iniciamos el proceso con:
Resultados
Corriente linea 1-2:

Resultados
Corriente línea 1-2:

Resultados

Resultados
Corriente línea 1-2: Potencia en la barra Slack:

Resultados
Pérdidas en el Sistema = Pot Slack –

Corriente línea 1-3: Pot Carga

Resultados
Corriente línea 1-3: Pérdidas en el Sistema = Pot Slack – Pot

Carga

Gracias

MODELADO DEL GENERADOR
SÍNCRONO Y CURVA DE CAPABILIDAD
- CIRCUITO EQUIVALENTE.
- CURVA DE CAPABILIDAD.
Gómez Palacio, Durango.

23 de Julio de 2009
Instructor: Raúl García Kasusky
CIRCUITO EQUIVALENTE
Introducción
El generador síncrono es uno de los elementos mas importantes

de un sistema de potencia, ya que éste se encarga de generar la
energía eléctrica que será transmitida a grandes distancias para
ser posteriormente utilizada por los usuarios.
El modelado del generador depende del tipo de análisis que se

pretenda realizar, el enfoque dado en esta presentación será el
correspondiente a un modelo simplificado para el análisis de
estado estable.
MÁQUINA SÍNCRONA TRIFÁSICA
Eje de campo
(CD) Eje de la fase A
c’ b Rotor
Entrehierro
Estator
a a’
If
Eje de la fase C Eje de la fase B
b’ c
Maquina síncrona de dos polos

Generador síncrono (Tipos de rotor)
GENERADOR DE POLOS GENERADOR DE POLOS
LISOS (ROTOR CILÍNDRICO) SALIENTES
S S
N N N N
S S
4 polos 4 polos
Generador síncrono de polos lisos (Rotor Cilíndrico)
Circuito equivalente Diagrama fasorial


Ra jX s E
 δ  
I  V jX s I
 V
E φ 
 Ra I
I
La referencia es el voltaje de terminales:
Ecuación de voltaje del circuito: 
   V = V∠0
E = V + ( Ra + jX s ) I 
E = E∠δ
Donde:
 La impedancia de la maquina:
E = Voltaje interno. 
Ra = Resistencia de armadura. Z = Ra + jX s = Z∠α
X s = Reactancia síncrona. La magnitud del voltaje interno es proporcional a la corriente
 de campo:
V = Voltaje en terminales. ωMf If
I = Corriente E=
2
Generador síncrono de polos lisos (Rotor cilíndrico)
Potencia de generación:
  ∗ Si Ra se desprecia:
S =V I 
Z = jX s y α = 90
Donde:
  Donde:
 E − V E∠δ − V∠0
I =  = 
Z Z VE
P= sen δ
 ∗ Ee − jδ − V Xs
I = V
Ze − jα Q= [E cos δ − V ]
Xs
Entonces:
 VE V2 VE V2
S= cos(α − δ ) − cos α + j sen(α − δ ) − j sen α
Z Z Z Z
Donde:
VE V2
P= cos(α − δ ) − cos α
Z Z
VE V2
Q= sen(α − δ ) − sen α
Z Z
Generador síncrono de polos salientes
S Diagrama fasorial

Iq E Eje q
φ
jI q X q
N
Id 
I  jI d X d
V
Eje d
VE V2 1 1
P= sen δ + ( − ) sen(2δ )
Xd 2 Xq Xd
Potencia de generación:
VE V 2 Xd − Xq V2 1 1
Q= cos δ + ( ) cos(2δ ) − ( + )
Xd 2 Xd Xq 2 Xd Xq
Condensador Síncrono
Suponiendo: Ra = 0

E

Generador sobreexcitado: jX d I I X d cos φ
δ 
φ V I X d sen φ

I
 Lugares
E geométricos
Generador con excitación de potencia
normal: jI X d constante
para :
δ E cos δ 
 V E, I
I
I X d sen φ

Generador subexcitado: E 
I  I X d cos φ
jX d I
φ δ
Condensador Síncrono
Suponiendo: Ra = 0

E

Generador sobreexcitado: jX d I I X d cos φ
δ 
φ V I X d sen φ

I
 Lugares
E geométricos de
Generador con excitación potencia
normal: jI X d constante para :
δ E, I
E cos δ 
 V
I
I X d sen φ

Generador subexcitado: E 
I  I X d cos φ
jX d I
φ δ

2022 Flujo de Carga Part1

Cargado por

Copyright:

Formatos disponibles

2022 Flujo de Carga Part1

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

2022 Flujo de Carga Part1

Cargado por

Copyright:

Formatos disponibles

Flujos de carga en Sistemas de Potencia

Prof. Juan Bautista R.

Asignatura: Análisis de Sistemas de Potencia I

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Generador Transformador Transformador Generador

Prof. Juan Bautista R. (FIEE - UNI)

Generador Transformador Transformador Generador

Prof. Juan Bautista R. (FIEE - UNI)

Fuente de corriente en paralelo

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Corriente total saliente

Esta última es linealmente dependiente de las otras tres ecuaciones

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

• Permite determinar los flujos de potencia activa y

Prof. Juan Bautista R. (FIEE - UNI)

Inyección de Corriente en la barra “k”:

Prof. Juan Bautista R. (FIEE - UNI)

Inyección de Corriente en la barra “k”:

Prof. Juan Bautista R. (FIEE - UNI)

Inyección de Corriente en la barra “k”:

Prof. Juan Bautista R. (FIEE - UNI)

Inyección de Corriente en la barra “k”:

Prof. Juan Bautista R. (FIEE - UNI)

Inyección de Corriente en la barra “k”:

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Las partes real e imaginaria de la potencia en la barra k son:

Para la potencia activa:

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)

Prof. Juan Bautista R. (FIEE - UNI)