FLUJO DE POTENCIA

EN LOS SISTEMAS DISTRIBUCIN

DE ENERGA ELCTRICA (SDEE)

Caractersticas de los Sistemas de

Distribucin de Energa Elctrica
.

Mtodos de Clculo de Flujo de

Potencia
.

Modelacin de los Componentes para

estudios de flujo de potencia

Caractersticas de los Sistemas

de Distribucin Energa Elctrica
Estructura radial o dbilmente
enmallada
Sin transposicin
Elevados valores de R/X
Naturaleza multifsica
Carga distribuida
Existencia de una sola barra fuente
Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodos de clculo de flujo de potencia

en los sistemas de distribucin
Mtodos Bsicos y Trifsicos:
.

1. Mtodo escalonado (Ladder method)

2. Suma de corrientes (Current summation)
3. Mtodo suma de potencias
4. Newton Rapson desacoplado rpido
modificado
5. Mtodo de Gauss-Seidel implcito
Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D.

Mtodo Escalonado Monofsico (Ladder Method)

Se suponen conocidas las tensiones de nudos: 1 [pu]

S c 2
S c 3 V2 (k ) V3 (k ) Z 2 I3 (k )
I3 (k )

I2 (k )
I3 (k )
V2 (k )
V3 (k )
*

V1 (k ) V2 (k ) Z1 I2 (k )
E rror

Se corrige:

E rror V1 V1 (k )

V3 (k 1) V3 (k ) E rror
Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D.

Mtodo Escalonado: Sistema radial con laterales

Considrese la barra 3 como
extremo del alimentador :
a) Se estima un valor para
V3 y se calculan :

S c*3
I3 *
V3 V2 V3 Z 2 I3
b) Considerando inicialmente el lateral desconectado, se estima V4 :

V debe ser igual a V (dentro de cierto error)

2
2
I4 (Sc 4 / V4 ) * V2' V4 Z 3 I4 antes de reconectar el lateral

E 2 V2 V2' V4 (nueva) V4 (antigua ) E 2

I2 (Sc 2 / V2 )* I3 I4

donde

V1 V2 Z1 I2

E1 V1 (especifica da) V1 (calculada )

El proceso se repite hasta que

Cuando |E2|< se conecta

el lateral y se calcula I1 e I2
Se calculan:

V3 (nueva) V3 (antigua ) E1

E1
Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D: Mtodo Escalonado

Sistema radial con mltiples laterales

Sistema radial con laterales compuestos

Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D.

Mtodo Suma de Corrientes (Current Summation Method)

Proceso iterativo,
se inicia estimando
los valores para
todas las tensiones
de nudo, excepto
el extremo
transmisor,
en una iteracin k,
se tiene:

a) Partiendo desde el extremo receptor:

I3 (k ) Sc*3 / V3* (k )

I2 (k ) Sc*2 / V2* (k ) I3

b) Partiendo desde el extremo transmisor:

V2 (k ) V1 Z1 I2 (k )

V3 (k ) V2 (k ) Z 2 I3 (k )

En cada barra se comprueba si:

V2 (k ) V2 (k 1)
V3 (k ) V3 (k 1)
Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D.

Mtodo Suma de Potencias Monofsico

Durante una iteracin, se recorre la red en dos sentidos:(inverso y directo)

Partiendo de un perfil plano de tensin para todas las barras, las siguientes
ecuaciones describen el proceso iterativo en la k-sima iteracin:

a) Sentido inverso (aguas arriba): A partir de ltima barra, se calculan las

potencias perdidas y en cada tramo se suman los consumos
.

P2 (k ) R2 | I 2 | R2
2

| S C 2 |2

| V 2 (k ) | 2

A2 (k ) PC 2 PC 3 P2 (k )

X2
Q2 (k ) X 2 | I 2 |
P2 (k )
R2
2

B2 (k ) QC 2 QC 3 Q2 (k )
Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodo Suma de Potencias Monofsico (continuacin)

b) Sentido directo (aguas abajo): Se calculan los mdulos de las

tensiones usando una ecuacin bicuadrtica

V24 (k ) 2 A2 (k ) R1 B2 (k ) X 1 V12 V22 (k ) A22 (k ) B22 (k ) R12 X 12 0

X 1 A2 (k ) R1 B2 (k )
tg 2 (k )
R1 A2 (k ) X 1 B2 (k ) V22 (k )

V34 (k ) 2R2 PC 3 X 2 QC 3 V22 (k ) V32 (k ) PC23 QC2 3 R22 X 22 0

X 2 PC 3 R2 QC 3
'
tg 3 (k )
R2 PC 3 X 2 QC 3 V32 (k )
Criterio de convergencia:

3 (k ) 3' (k ) 2 (k )
(k )

( k 1)

|V p | |V p

Para los sistemas con laterales,el mtodo es aplicable en forma directa

Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D..

Mtodo de Newton-Raphson Desacoplado Rpido

El MNR-DR en su forma normal, est definido por:

P / V [ B ]

(1) donde:

Q / V [ B ]V

(2)

'

''

'
B pq
1

X pq

G pq

B pq
B pq

'
''
''
B pp
B pq B pp
B pp B pq
B pq

Durante los los ltimos aos se han desarrollados diversos procedimientos

para aplicar el MNR-DR a sistemas de distribucin:
Aumentar artificialmente la razn X/R mediante la rotacin de la impedancia
Modificar sustancialmente las matrices jacobianas B y B (ms reciente)
En redes de transmisin la razn R/X es suficientemente pequea y las
hiptesis del desacoplamiento son verdaderas
L Wang 1990, encontr en forma emprica, que el MNR-DR converge mejor
si los Gpq y Bpq son similares.
Esto indica que la hiptesis de desacoplamiento Gpq<<Bpq no afecta en forma
crtica la convergencia, a menos que esta no se cumpla en un amplio margen.
Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodo de Newton-Raphson Desacoplado Rpido modificado

Segn L.Wang.la convergencia puede ser mejorada incluyendo los efectos de la
conductancia en las iteraciones. Para esto conviene sumar las ecuaciones de
potencia neta, escritas en la siguiente forma:

Pp V pVq G pq cos( p q ) B pq sen( p q )

n

q 1

Q p V pVq G pq sen( p q ) B pq cos( p q )

n

q 1

de donde:

Pp Q p V pVq (G pq B pq ) cos( p q ) (G pq B pq sen( p q )

n

q 1

La ecuacin (2) es sustituida por:

P Q [ B ]V
'''

donde:

'''
'''
B pp
G pp B pp B pq
G pq B pq

Para la ecuacin (1), conviene modificar sus trminos por:

P / V [ B ]
'

donde:

'
'
B pp
B pq

donde: K = 0,5 cuando 0 R/X 1; K = 1,0 cuando R/X 1

'
B pq
B pq KG pq

B pq G pq
B pq G pq

L. Wang y otros, Novel decoupled power flow, IEEE Proc. Vol. 137, Pt. C. N1, January 1990, p.1-7.
Distribucin de la Energa Elctrica: Luis Ortiz N.

Flujo de potencia en S.D.

Mtodo de Newton-Raphson Desacoplado Rpido

modificado mediante rotacin en el plano de la impedancia
El MNR-DR es aplicable para elevados
valores de X/R
Una solucin para aplicarlo a los SED
consiste en aumentar artificialmente la
relacin X/R mediante una rotacin en
el plano de la impedancia

Z Z

Z ' Z

Impedancia de una lnea

cualquiera
Impedancia transformada

Se cumplen las siguientes relaciones:

R ' R cos X sen

X R sen X cos
'

Re
.

donde, X/R > X/R

Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodo de Newton-Raphson Desacoplado Rpido modificado

Una rotacin en el plano de la impedancia afecta
las potencias complejas en barras y al flujo de
potencia de la red, pero mantiene inalterable las
tensiones.
La potencia neta en una barra
p cualquiera considerando la
rotacin de impedancias, ser:

S Vp*Vq Ypq'
'
p

p 1

Re

'
'
Ypq
Y pq pq

elemento de la matriz
admitancia transformada

Ypq Y pq pq

elemento de la matriz admitancia original

'
pq
pq

donde:

'
Pp' V pVqY pq cos( p q pq
) o bien, Pp' Pp cos Q p sen

'
Q p' V pVqY pq sen( p q pq
)

Q p' Pp sen Q p cos

Para calcular el ngulo de rotacin se han propuestos diferentes criterios. El mtodo

presenta problemas cuando existen ramas cuya razn X/R vara en un rango muy elevado
Algunas lneas pueden generar impedancias transformadas de parte real negativa.
Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodos trifsicos de flujo de potencia en S.D.

Modelacin trifsica con desacoplamiento

El acoplamiento mutuo entre las fases de las lneas de cada seccin
de un alimentador, se modela empleando impedancias propias y
fuentes de tensin dependientes de la corrientes.

Como consecuencia del empleo de las fuentes de tensin dependientes: Vab, Vac, etc,
se hace necesario calcular tensiones en nudos intermedios en cada fase, como x e y
Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodos trifsicos de flujo de potencia en S.D.

Mtodo de Gauss-Seidel Implcito

a) Se forma la matriz impedancia de barras:

VB VF Z B I B

VB

Vector de tensiones de barras (n-1)x1

VF

Vector cuyos elementos son todos iguales a la tensin de la barra flotante (n-1)x1

Z B
IB

Matriz impedancia de barras, formada considerando la barra flotante como referencia y

excluyendo todos los elementos shunt conectados a las barras (n-1)(n-1)

Vector de corrientes inyectadas, calculadas considerando las potencias netas y los elementos
shunt como fuentes de corrientes (n-1)x1

b) Se adoptan valores iniciales para todas las tensiones barras, excepto la flotante
c) Se calculan las corrientes inyectadas a las barras:

Iq( k )

Pq jQ q
( k 1)

y
V
q q
V ( k )
q

d) Se calculan las nuevas tensiones:

( k 1)
p

p 1

V f Z pq Iq( k 1) Z pq Iq( k )
q 1

Los pasos b) y c) se repiten hasta que se cumpla el criterio de convergencia

Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodos trifsicos de flujo de potencia en S.D.

Mtodo de Gauss-Seidel Implcito

Se modelan todos los componentes de la red trifsica, se aplica el
MGS-ZB
La matriz impedancia de barras se puede obtener de la matriz admitancia
de barras bifactorizada
En los alimentadores de distribucin el caso tpico contiene slo una barra
de tensin controlada que corresponde a la SE de alimentacin

Distribucin de la Energa Elctrica: Luis Ortiz N.

Mtodos trifsicos de flujo de potencia en S.D.

Mtodo de Gauss-Seidel Implcito

El mtodo de Gauss-Seidel MGS-YB bifactorizada se puede interpretar, aplicando el
principio de superposicin: que la tensin de barra es el resultado de las componentes
de la tensin en la SE y las fuentes de corrientes inyectadas en las barras

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia

De acuerdo a las conexiones usuales de los transformadores, corresponde analizar:
M.T: Lneas trifsicas sin retorno por tierra; B.T: Lneas trifsicas con neutro conectado a tierra

a) Lnea trifsica sin neutro ni retorno por tierra

Va Va Va ' Z aa

Vb Vb Vb ' Z ba
Vc Vc Vc ' Z ca

Z ab
Z bb
Z cb

o bien:
V Z I

Z ac I a
donde: Zii=Rii+jXii

Z bc I b
Zij=iXij
i, j = a, b, c
Z cc Ic X ij 0 f n(1 / Dij ) (/Km)
Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

b) Lnea trifsica con neutro y retorno por tierra

Va Va Va '
Ia
V V V abcn I
b b b' Z b
Vc Vc Vc ' Ic


Vn Vn Vn '
In

o bien:

abc

V
Z

V
Z nf
n

abc
Z fn I

I
Z nn n

abc

Considerando Vn=0:

abc

Z n I abc
abc

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

Lnea trifsica con neutro y retorno por tierra

Z n Z Z fn Z nn1 Z nf
abc

abc

Z aa g

Z ba g
Z ca g

Z ab g
Z bb g
Z cb g

Z ac g

Z bc g
Z cc g

abc

Z n I abc
abc

De las ecuaciones anteriores el tramo de lnea con retorno por tierra, se representa por:

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

Lnea trifsica con neutro y retorno por tierra

Z
abcn

Z n
abc

Matriz de 4x4, con trminos Zij

Matriz de 3x3, con trminos Zij-g

}Z

ij g

Z ij

Z in Z nj

i,,j = a,b,c

Z nn

Las expresiones para las impedancias propias y mutuas, considerando el efecto de retorno por
tierra, en forma simple, con un grado de aproximacin aceptable de la serie de Carson, son:

Z ii g ( Rii 9,8696 10

Z ij g 9,8696 10

f ) j 2,8935 10

f j 2,8935 10

659 / f
f og10
(/Km)
RMG

659 / f
f og10
(/Km)
Dij

f = frecuencia, en Hz
= resistividad del terreno, en -m
RMG = radio medio geomtrico del conductor, en m
Dij = distancia entre los conductores i, j, en m
Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

c) Lnea bifsica con neutro a tierra
En una lnea bifsica que utilice por ejemplo las fases a y c, las ecuaciones bsicas conducen
a una matriz de 3x3 (a,c,n), la cual se puede reducirse a 2x2, empleando la condicin Vn=0.

Para la modelacin trifsica se debe expandir a una matriz de 3x3, introduciendo ceros en la
lnea y columna correspondiente a la fase no existente (b en este caso):

Va Vag Va ' g ' Z aan 0 Z acn Ia

0

0
0
0
b


Vc Vcg Vc ' g ' Z can 0 Z ccn Ic

Z 2
abc

Z aan 0 Z acn
0
0
0
Z can 0 Z ccn
Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

d) Lnea monofsica con neutro a tierra

In Ia

Va Vaa' Van Va 'n '

Va (Z aa Ia Z an In ) (Z nn In Z an Ia )

La modelacin trifsica queda:

Va (Z aa Z nn 2Z an ) Ia

Z1 (Z aa Z nn 2Z an )

Va Van Va 'n ' Z1

0
b
0
Vc
0
0

0 0 Ia

0 0 0
0 0 0

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de componentes para estudio de flujo de potencia en S.D.

Modelo trifsico de lnea empleando la matriz admitancia de barra
Para ello es necesario determinar la matriz nodal de cada lnea (de 6x6), a partir de la matriz
admitancia primitiva correspondiente

Z Y (Z )
abc

abc

abc

Matrices de impedancias propias y mutuas de la lnea trifsica (3x3)

y matriz admitancia primitiva de la lnea trifsica (3x3)

La ecuacin nodal de (6x6), es:

abc
abc
abc abc
Y V I
Yabc
a 'b 'c '
a 'b 'c '
abc
Y Y V I

Matriz nodal de la lnea

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelos de carga para estudio de flujo de potencia en SED

En los SED, existe una gran variedad de cargas que requieren: P y Q.
Ambos tipos de potencia dependen fundamentalmente de V y f
Las cargas pueden representarse por diversos modelos:

Modelo bsico: potencia constante, impedancia (o admitancia) constante y

corriente constante
Modelo incremental: considera variaciones de P y Q con la tensin y frecuencia.
Se realiza con coeficientes diferenciales semiempricos
Modelo exponencial: considera variaciones de tensin y frecuencia.
Los coeficientes para P y Q son caractersticos de cada tipo de carga
Modelo polinomial (2 o 3 orden), puede incluir variaciones de tensin y
frecuencia, por ejemplo refrigeracin:

P 0,798 0,606V 1,146V 2 0,418f 2,69Vf

Q 0,624 1,54V 3,37V 2 0,889f 7,37Vf
Los valores de V y f estn expresadas en pu.
Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelacin de cargas para estudio de flujo de potencia en S.D.

Concepto de Ventana de Carga

Distribucin de la Energa Elctrica: Luis Ortiz N.

Ventana de Carga

Distribucin de la Energa Elctrica: Luis Ortiz N.

Concepto de Ventana de Carga aplicado a un alimentador

Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelo de transformadores
Modelo por fase de un transformador:

Otra expresin, obtenida experimentalmente, para determinar las prdidas en

el ncleo de transformadores, es:

( KVA) nom
2
C V 2
Pn ( pu )
A V Be
SB
( KVA) nom
2
F V 2
Qn ( pu )
D V Ee
SB

Valores tpicos para las constantes:

A=0,0067;
B=0,73410-9
C=13,5
D=0,00167
E=0,268310-3
F=22,74
V = mdulo de la tensin
Distribucin de la Energa Elctrica: Luis Ortiz N.

Modelos de condensadores paralelo

Ventana de Carga

Condensadores paralelo

Vi
Qi Q 0
Vi
0
i

Modelo de admitancia constante

Qi0 Vi 0

: Valores iniciales; i = a,b,c

Distribucin de la Energa Elctrica: Luis Ortiz N.