CAPÍTULO 6

FLUJO DE CARGA LINEALIZADO

6.1 Introducción

El flujo de carga linealizado llamado también flujo de carga C.C. es una forma aproximada
de resolver ecuaciones de flujo de carga no lineal. En la deducción del flujo de carga C.C. es
llevada en cuenta la relación estrecha entre el flujo de potencia activa en la lı́nea de transmisión
y la abertura angular entre los extremos de la lı́nea.

Pkm ⇒ δ = θk − θm

6.1.1 Motivaciones para el uso del Flujo de Carga C.C.

• Los valores obtenidos con este método simplificado presentan resultados de flujos de
potencia activa en las lı́neas con valores muy próximos a los obtenidos usando el modelo
no lineal (≈ 5 − 10%)

• El esfuerzo computacional es muy pequeño y robusto.

• En determinados tipos de problemas, como en el planeamiento de sistemas de transmisión,

la estructura de los datos lleva con frecuencia a la no convergencia del flujo de carga no
lineal convencional. En estos casos es mejor una información de flujos de potencia activa
que permitirı́a encontrar puntos crı́ticos en el sistema (entregada por el flujo de carga
c.c.); en estos casos el problema no lineal puede no converger y la información entregada
por este no serı́a útil.

• La aproximación es válida solamente para sistemas de transmisión; esto es, para niveles
elevados de tensión. La aproximación es mejor cuanto mayor fuese el nivel de tensión.

• El modelo de flujo de carga C.C. no tiene en cuenta las magnitudes de las tensiones
nodales y el flujo de potencia reactiva.

124
6.2 Linealización de las Ecuaciones No Lineales del Flujo
de Carga

6.2.1 El flujo de Potencia Pkm de una Lı́nea de Transmisión es:

Pkm = Vk2 gkm − Vk Vm gkm Cosθkm − Vk Vm bkm Senθkm (6.1)

Pmk = Vm2 gkm − Vk Vm gkm Cosθkm + Vk Vm bkm Senθkm (6.2)

De (6.1) y (6.2) se determinan las pérdidas ası́:

Pe = Pkm + Pmk = gkm [Vk2 + Vm2 − 2Vk Vm Cosθkm ] (6.3)

Despreciando las pérdidas en la lı́nea; esto es, considerando la conductancia de la lı́nea

despreciable, entonces

gkm Vk [Vk − Vm Cosθkm ] ≈ 0

Pkm = −Pmk = −Vk Vm bkm Senθkm (6.4)

Introduciendo además las siguientes aproximaciones:

Vk ≈ Vm ≈ 1p.u. (6.5)

Senθkm ≈ θkm (en radianes) (6.6)

−xkm 1
bkm = 2 2
≈− (6.7)
rkm+ xkm xkm

Tomando en cuenta las ecuaciones (6.5), (6.6), y (6.7) en la ecuación (6.4), se llega a:

θkm θk − θm
Pkm = = (6.8)
xkm xkm

La ecuación (6.8) tiene la forma de la Ley de Ohm aplicada a una resistencia, tal como se
presenta en la figura 6.1

125
 i km P km
Vk θk

Vk - Vm θk - θm
i km = r km r km P km = x km x km

Vm θm

Figura 6.1: Comparación de modelos de corriente continua y D.C.

6.2.2 El flujo Pkm de un Transformador en fase es:

Pkm = (akm Vk )2 gkm − (akm Vk )Vm gkm Cosθkm − (akm Vk )Vm bkm Senθkm (6.9)

Pmk = Vm2 gkm − (akm Vk )Vm gkm Cosθkm + (akm Vk )Vm bkm Senθkm (6.10)

Con las simplificaciones introducidas en (6.5) - (6.7), se tiene:

akm θkm
Pkm = −Pmk =
xkm

y haciendo la aproximacón akm ≈ 1, se tiene

θkm θk − θm
Pkm = −Pmk = = (6.11)
xkm xkm

Que tiene la misma forma de la deducida para una lı́nea de transmisión.

6.2.3 El flujo Pkm de un desfasador puro es dado por la relación:

Pkm = Vk2 gkm − Vk Vm gkm Cos(θkm + ϕkm ) − Vk Vm bkm Sen(θkm + ϕkm ) (6.12)

Pmk = Vm2 gkm − Vk Vm gkm Cos(θkm + ϕkm ) + Vk Vm bkm Sen(θkm + ϕkm ) (6.13)

Introduciendo las simplificaciones (6.5) - (6.7), se tiene:

Vk2 gkm − Vk Vm gkm Cos(θkm + ϕkm ) = Vk gkm [Vk − Vm Cos(θkm + ϕkm )] ≈ 0

126
 Vk ≈ 1; Vm ≈ 1; Cos(θkm + ϕkm ) ≈ 1

Entonces (6.13) se reduce a:

Sen(θkm + ϕkm )
Pkm = (6.14)
xkm

La abertura angular (θkm + ϕkm ) entre los extremos de Ykm es relativamente pequeña;
entonces es posible la aproximación: Sen(θkm + ϕkm ) ≈ (θkm + ϕkm ).
Ası́, (6.14) se reduce a lo siguiente:

(θkm + ϕkm )
Pkm = −Pmk = (6.15)
xkm

Pkm tiene entonces dos componentes, una variable que depende del estado del sistema
θkm /xkm y otra fija e igual a ϕkm /xkm

En la siguiente figura se representa adecuadamente el modelo linealizado con un desfasador

puro:

θk θm
k m

θp = θk + ψ km

(θ km + ψ km )
P km = x km

Esquema equivalente

θk θm
k m

x km

ψ km ψ km
- x km x km

Figura 6.2: Modelo del Transformador Desfasador Puro

127
 Donde la parte del flujo invariante ϕkm /xkm aparece como una carga adicional en la barra
k, y una generación adicional en la barra m.

6.3 Formulación Matricial del Flujo de Carga C.C.


P =Bθ

Considerando que el sistema eléctrico tiene solamente lı́neas de transmisión (sin transfor-
madores en fase o desfasadores). Entonces, el flujo de potencia Pkm es dado por:

Pkm = x−1
km θkm (6.16)

La inyección de potencia activa en la barra k es igual a la suma de los flujos que salen de la
barra k; ası́:

1
k
P Gk P k1

P k2 2

P kn
n

P Dk

Figura 6.3: Ecuación de balance de potencia nodal

De la figura (6.3) se tiene que:

Pk = PGK − PDK = Pk1 + Pk2 + · · · + Pkn

Expresados los flujos de potencia activos en función de (xkm y θkm )


Pk = x−1
km θkm k = 1, 2, · · · , nb (6.17)
m∈Ωk

Descripción de la ecuación (6.17)

1 1 1
Pk = (θk − θ1 ) + (θk − θ2 ) + · · · + (θk − θn )
xk1 xk2 xkn

128
 Que puede ser separada en dos componentes:

 
Pk = x−1
km θk + − x−1
km θm (6.18)
m∈Ωk m∈Ωk

Pk = − xθk1
1
− θ2
xk2
− θ3
xk3
−··· − θn
xkn
+ θk
xk1
+ θk
xk2
+ θk
xk3
+ ···+ θk
xkn


Pk = − xθk1
1
− θ2
xk2
··· + 1
xk1
+ 1
xk2
+···+ 1
xkn
θk · · · − θn
xkn

Pk = − xθk1
1
− θ2
xk2
··· + θk
xkk
··· − θn
xkn

Expresarlo en forma vectorial es como sigue:

  ⎡ ⎤

n θ1
− x1k1 − 1
···+ 1
···− 1 ⎢ ⎥
xk2 xkl xkn ⎢ ⎥
⎢ ⎥
l=1 ⎢ θ2 ⎥
⎢ ⎥
⎢ .. ⎥
Pk = ⎢ . ⎥
⎢ ⎥
⎢
⎢ θk ⎥ ⎥
⎢ .. ⎥
⎢ ⎥
⎣ . ⎦
θn

(6.18) Puede ser representada en forma matricial; ası́: P = B  θ (6.19)

Sistema de 4 nodos
⎡ ⎤
⎡ ⎤ 
4
⎡ ⎤
P1 ⎢
⎢
1
x1n
− x112 − x113 − x114 ⎥
⎥ θ1
⎢ ⎥ ⎢ n=1 ⎥⎢ ⎥
⎢ ⎥ ⎢ 
4 ⎥⎢ ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ P2 ⎥ ⎢ − 1 1
− x123 −− 1 ⎥⎢ θ2 ⎥
⎢ ⎥ ⎢ x21 x2n x24 ⎥⎢ ⎥
⎢ ⎥ ⎢ n=1 ⎥⎢ ⎥
⎢ ⎥ = ⎢ ⎥⎢ ⎥
⎢ ⎥ ⎢ 
4 ⎥⎢ ⎥
⎢ P3 ⎥ ⎢ − 1 − x132 1
− x134 ⎥⎢ θ3 ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ ⎥ ⎢ x31 x3n ⎥⎢ ⎥
⎣ ⎦ ⎢ n=1 ⎥⎣ ⎦
⎢ 
4 ⎥
P4 ⎣ ⎦ θ4
− x41
1
− x142 − x143 1
x4n
n=1

θ ⇒ Vector de los ángulos de las tensiones nodales

Siendo que: P ⇒ Vector de las inyecciones lı́quidas de potencia activa

B ⇒ Matriz del tipo admitancia nodal y cuyos elementos son:

 1
Bkm = −x−1
km = − (6.20)
xkm

129
   1

Bkk = x−1
km = (6.21)
m∈Ωk m∈Ωk xkm


La matriz B en (6.19) es singular y se debe eliminar una de las ecuaciones y adoptar para
la barra correspondiente como referencia angular θk = 0
La matriz es singular, ya que las pérdidas son despreciadas y la suma de las componentes
de P es nula; esto es, la inyección de potencia en una barra cualquiera puede ser obtenida a
partir de la suma algebraica de las demás. Ası́, ahora se transforma en un sistema no singular
de dimensión nb − 1. La solución de este sistema algebraico entrega los ángulos de las nb − 1
barras restantes.
Si el sistema eléctrico tiene transformadores en fase o desfasadores, entonces el sistema de
ecuaciones (6.19) continúa válido con las siguientes observaciones:


• La formación de la matriz B cuando existen transformadores en fase o desfasadores, es
realizada exactamente igual que en el caso de las lı́neas de transmisión.

• En la formación del vector P, se debe considerar la representación equivalente de los

desfasadores; esto es, se debe adicionar el flujo fijo ϕkm /xkm en las barras k y m para
cada transformador desfasador puro.

Ejemplo: Dado el siguiente sistema de Tres (3) Barras, calcular los flujos de potencia activa
por las lı́neas.

1 2
P 12 X 12 = 1/3
P 1 = 1.5 P 2 = -0.5

P 13 P 23

X 13 = 1/2 X 23 = 1/2

P 3 = -1.0

Figura 6.4: Sistema de Tres (3) Barras

130
 ⎡ ⎤ ⎡ ⎤⎡ ⎤
1, 5 5 −3 −2 θ1
 ⎢ ⎥ ⎢ ⎥⎢ ⎥
P = Bθ ⇒ ⎣ −0, 5 ⎦ = ⎣ −3 5 −2 ⎦ ⎣ θ2 ⎦
−1, 0 −2 −2 4 θ3

Considerando la referencia angular en la barra 1, θ1 = 0 se tendrá el siguiente sistema

reducido:
    
−0, 5 5 −2 θ2
=
−1, 0 −2 4 θ3

   −1     
θ2 5 −2 −0, 5 1/4 1/8 −0, 5
= = =
θ3 −2 4 −1, 0 1/8 5/16 −1, 0

   
θ2 −1/4
=
θ3 −3/8

Cálculo de Flujos
P12 = θ12
x12
= 0−(−1/4)
1/3
= 3
4
⇒ P12 = 0, 75

P13 = θ13
x13
= 0−(−3/8)
1/2
= 3
4
⇒ P13 = 0, 75

−1/4−(−3/8)
P23 = θ23
x23
= 1/2
= 1
4
⇒ P23 = 0, 25

6.4 Representación de las pérdidas en el modelo de Flujo

de Carga C.C.

En el modelo C.C. se puede representar el efecto de las pérdidas en el sistema de transmisión.

En este las pérdidas son modeladas como cargas distribuidas en todo el sistema y son atendidas
por la barra de referencia.

6.4.1 Flujo de carga A.C.:

Potencia generada en la barra de referencia = pérdidas de transmisión + carga neta de todas

las demás barras del sistema.
El flujo de potencia activa inyectada en la barra k, es dada por:


Pk = Vm (Gkm Cosθkm + Bkm Senθkm ) (6.22)
m∈k

131
En que k es el conjunto de las barras vecinas de la barra k incluida la propia barra k.
Considerando Vk = Vm = 1p.u. de (6.22) tenemos:


Pk = Gkk + (Gkm Cosθkm + Bkm Senθkm ) (6.23)
m∈Ωk

Donde Ωk es el conjunto de barras vecinas de la barra k, excluida la barra k.

Para una lı́nea de transmisión o transformadores se tiene:

Gkm = −gkm (6.24)


Gkk = gkm (6.25)
m∈Ωk

1
Bkm ≈ (6.26)
xkm

Empleando (6.24), (6.26) en (6.23), se tiene:

   1
Pk = gkm − gkm Cosθkm + Senθkm
m∈Ωk m∈Ωk m∈Ωk xkm

 
Pk = (1 − Cosθkm )gkm + x−1
km Senθkm (6.27)
m∈Ωk m∈Ωk

Usando aún la aproximación:

2
θkm
Cosθkm ≈ 1 − 2

Senθkm ≈ θkm (en radianes)

Entonces (6.27) asume la siguiente forma:

1  
Pk = 2
gkm θkm + x−1
km θm (6.28)
2 m∈Ωk m∈Ωk

132
 Comparando las ecuaciones (6.28) y (6.17), se observa que la diferencia está en el término:

1  2
gkm θkm (6.29)
2 m∈Ωk

Que representa las pérdidas aproximadas de las lı́neas de transmisión que salen de la barra k.
Este hecho puede ser verificado analizando las pérdidas de transmisión en la lı́nea k-m.

Pe = Pkm + Pmk = gkm (Vk2 + Vm2 − 2Vk Vm Cosθkm ) (6.30)

Haciendo las siguientes aproximaciones en (6.30), se tiene:

2
θkm
Vk = Vm = 1p.u. y Cosθkm = 1 −
2

θ2
Pe = gkm 2 − 2 1 − km ⇒ Pe = gkmθkm
2
(6.31)
2

Ası́, (6.29) representa la mitad de las pérdidas activas de todas las lı́neas adyacentes a esta
barra. Ası́, el efecto de las pérdidas pueden ser representadas de manera aproximada como
cargas adicionales obtenidas al dividir las pérdidas de cada lı́nea del sistema entre sus barras
terminales (mitad para cada lado).
Por lo tanto, cuando en el modelo de flujo de carga C.C., se consideran las pérdidas en
forma aproximada; este asume la siguiente forma:


P + Ppérdidas = B θ (6.32)

Forma aproximada para resolver (6.32):


• Resolver el sistema P = B θ̃ sin tener en cuenta las pérdidas.

• Calcular las pérdidas aproximadas empleando (6.31) y la solución θ̃ obtenida anterior-

mente. Distribuir estas pérdidas en todas las barras del sistema como cargas adicionales.

• Resolver el sistema (6.32) con los valores de pérdida obtenidas en el paso enterior. Con
los θ conocidos, se calculan los flujos en las lı́neas.

Ejemplo:
En el ejemplo anterior, recalcular los flujos considerando las pérdidas.
Asumir que rkm = 15 xkm

133
Cálculo aproximado de las pérdidas:

rkm
gkm = 2 +x2
rkm km

1
15
g12 = 1 2
15
2 = 26
( ) ( )
15
+ 13

1
5
g13 = 1 2
10
2 = 13
( ) ( )
10
+ 12

1
5
g23 = 1 2
10
2 = 13
( ) ( )
10
+ 12
  2
1
Pe12 = 2
g12 θ12 = 15
26
0− − = 15
26 × 26
= 0, 036
4
  2
3 5×9
2
Pe13 = g13 θ13 = 5
13
0− − = 13 × 64
= 0, 054
8
  2
3
2
Pe23 = g23 θ23 = 5
13
− 14 − − = 5
13 × 64
= 0, 006
8

1 2
P 1 = 1.5 P 2 = -0.5

-0.045 -0.021

P 3 = -1.0 -0.030

Figura 6.5: Sistemas de 3 nodos incluido el efecto por pérdidas

Al plantear el problema se obtiene:

       
θ2 1/4 1/8 −0, 521 θ2 −0, 259
= ⇒
θ3 1/8 1/16 −1, 030 θ3 −0, 387

134
Ahora es posible calcular una mejor aproximación de flujos:

θ12 0 − (−0, 259)

P12 = = 1 ⇒ P12 = 0, 777 p.u.
x12 3

θ13 0 − (−0, 387)

P13 = = 1 ⇒ P13 = 0, 774 p.u.
x13 2

θ23 −0, 259 − (−0, 387)

P23 = = 1 ⇒ P23 = 0, 256 p.u.
x23 2

Nuevo valor de P1:

P1 = P12 + P13 = 1, 551 p.u.

Ejemplo:
En el sistema de 3 barras anterior, sustituir la lı́nea 1-2 por un desfasador con las siguientes
caracterı́sticas:

ϕ12 = 5◦ , x12 = 1/3

Recalcular los flujos en las lı́neas: ϕ12 = 5◦ = 0, 0873rad

ϕ12 0, 0873
P12 = = 1 = 0, 2618
x12 3

1 2
P 1 = 1.5 X 12 = 1/3 P 2 = -0.5

-0.2618 0.2618

X 13 = 1/2 X 23 = 1/2

P 3 = -1.0

Figura 6.6: Sistema de 3 nodos incluido el efecto del desfasador

135
 ⎡ ⎤ ⎡ ⎤⎡ ⎤
1, 5 − 0, 2618 5 −3 −2 θ1
⎢ ⎥ ⎢
5 −2 ⎦ ⎣ θ2 ⎥
⎥⎢

P =Bθ ⇒ ⎣ −0, 5 + 0, 2618 ⎦ = ⎣ −3 ⎦
−1, 0 −2 −2 4 θ3
Considerando la referencia angular en la barra 1, entonces θ1 = 0
    
−0, 2382 5 −2 θ2
=
−1, 0 −2 4 θ3

   
θ2 −0, 1846
=
θ3 −0, 3423

 θ12 0,1846
P12 = X12
= 1 = 0, 5538 p.u.
3

 ϕ12
P12 = P12 + X12
= 0, 5538 + 0, 2618 = 0, 8156 p.u.

θ13 0,3423
P13 = X13
= 1 = 0, 6846 p.u.
2

θ23 −0,1846−(−0,3423)
P23 = X23
= 1 = 0, 3154 p.u.
2

Ası́, el flujo entre 1-2 pasó de 0, 75 [Link] 0, 8156 p.u.

Antes del desfasador P12 = 0, 75 pu.

Colocando el desfasador P12 = 0, 8156 p.u.

1 2
0.8156

P 1 = 1.5 P 2 = -0.5

0.6844 0.3156

P 3 = -1.0

Figura 6.7: Sistema de 3 nodos mostrando los flujos incluido el transformador desfasador

136