Potencia 1 Curso (Yuri)
Potencia 1 Curso (Yuri)
Potencia 1 Curso (Yuri)
ESTUDIANTES:
DOCENTE:
AREQUIPA
2009
INDICE
2
2.12 Compensación de una Línea
2.12.1 Reactores Shunt
2.12.2 Compensación capacitiva Shunt
2.12.3 Compensación capacitiva en serie
3
CONCEPTOS
FUNDAMENTALES
4
Capitulo I
Conceptos Fundamentales
1.1 Introducción
5
En general un sistema eléctrico presenta cierta estructura
que se utiliza no solo para describirlo sino para
modelarlo, simularlo, controlarlo y gestionarlo.
a) Generación:
Generador trifásico AC Generador Síncrono. Alternador.
Transformadores: transfieren la energía eléctrica a una
alta eficiencia.
b) Transmisión:
Transfiere energía desde los centros de producción hasta
los de carga.
c) Subtransmisión:
Niveles de tensión 33kV, 10kV, 22.9kV, 60kV, 138kV, 220kV.
d) Consumo, carga:
Comerciales, Residenciales, Industriales
6
Comerciales, Residenciales: Independientes de la frecuencia
y consumen poca energía reactiva.
Industriales: Cargas compuestas, motores de inducción.
Carga Promedio
Factor de Carga
Carga Pico
Maxima Demanda
F.U.
Capacidad Instalada
7
v t Vmcos ω t θV
i t Imcos ω t θi
Operando:
1 1
p t Vm.Im. cos θV θi cos 2ω t θV θi
2 2
1
p t .Vm.Im.cos θV θi cos 2ω t θV θi
2
1
p t .Vm.Im.cos θV θi cos 2. ω t θV θV θi
2
1
p t .Vm.I m.cos θ V θi cos 2. ω t θ V .cosθ V θi sen 2. ω t θ V .senθ V θi
2
Vm V. 2
Im I. 2
Además θ θV θi
p t V.I.cosθ1 cos 2ω t θV V.I.senθsen 2ω t θV (1.1)
V.I.cosθ1 cos 2ω t θV
Flujo de Energía en el circuito
Pr(t), en merced a la resistencia que ofrece el circuito.
8
px t V.I.senθ.sen 2ω t θV
Q V.I.senθ (1.4)
S P2 Q2 V.I
V V θv , I I θi
V.I* V e jθV .I e jθi V .I.ej θV θi V .I θv θi
S V.I* V.Iθ
Donde θ θv θi
S V .I cos θ j.V .I senθ (1.5)
S P j.Q
Triángulos de Potencias
9
Q es negativa cuando la carga es capacitiva y la corriente
adelanta al voltaje.
S S1 S2 S3 60 j80
S 10053.13º kVA
S * 100 53.13º
I 71.43 53.13º
V * 1400 0º
FP cos φ cos 53.13º 0.6 Retraso
Se conecta un capacitor de resistencia despreciable en
paralelo a las cargas para mejorar el FP. a 0.8 en atrás.
Calcular la potencia del capacitor y sus capacitancias.
P=60KW
V
36.87º
53.13º
QN
Q=80K
S=100KVA
I
Qc
QC Q QN 80 45
I
10
QC 35kVAR
2
V 14002
XC
QC 35 103
XC 56
Luego:
1
56 , f 60Hz , C 47.37μF
2πfC
S * 60 j45
I'
V * 1.4 j0
I' 53.57 36.87º
Ejemplo 1:
Solución:
760
Demanda máxima la carga= 950kVA
0.8
950
Demanda máxima del suministrador= 1000kVA
0.95
Perdidas en la distribución = 5%
Costos de los dispositivos para el caso=60x1000= $60000
Cargo anual de depreciación=0.12x60000= $7200
11
Y el costo anual por los kilo watts consumidos será:
Ejemplo 2:
Calcular:
Solucion:
Demanda Promedio
Factor de carga
Demmanda Maxima
Luego:
12
1.7 Flujo de Potencia Compleja
V1 V1 δ1
V2 V2 δ2
V1 V2 V δ1 V2 δ2
I12 1
Z Zγ
V V
I12 1 δ1 γ 2 δ2 γ
Z Z
S12 V1.I12 *
V V
S12 V1 δ1 1 γ δ1 2 γ δ2 (1.6)
Z Z
2
V V .V
S12 1 γ 1 2 γ δ1 δ2 Expresion Base
Z Z
2
V V .V
P12 1 cos γ 1 2 cosγ δ1 δ2 (1.7)
Z Z
2
V1 V .V
Q12 senγ 1 2 senγ δ1 δ2 (1.8)
Z Z
Z X90º
Luego:
V1.V2
P12 senδ1 δ2 (1.9)
X
13
1.- La expresion (1.9) muestra que P12 depende de V1, V2,
X, S1 y S2 pero X es un parametro del sistema (fijo) que no
se puede modificar, el angulo si se puede modificar ya que
la tension y la corriente en todo momento varian, este
angulo influye drasticamente en la potencia P12 ya que el
voltaje se trata de mantener constante.
V1.V2
P12
X
Ejemplo 3:
14
Solucion:
En la linea:
Luego:
15
alcanzan su valor pico en el orden secuencial ABC se dice
que el generador posee una secuencia positiva.
VAN 2Vfcos ω t θV
2π
VBN 2Vfcos ω t - θV 1
3
4π
VCN 2Vfcos ω t - θV
3
16
Si la carga es balanceada
ia 2Ifcos ω t θi
ib 2Ifcos ω t -
2π
θi 2
3
ic 2Ifcos ω t - θi
4π
3
Luego:
17
Y si mantenemos la simetria, la potencia trifasica compleja
sera:
S3φ P3φ jQ3φ
S3φ 3VfIf *
En una conexión estrella En una conexión triangulo
If I L Vf VL
VL 3Vf IL 3If
Vca
120 Van
120
120
Vbn
18
Vba
Van V 0º
Vbn V 120º
Vcn V 240º
Vab Van Vbn V 10º 1 120º 3 V 30º
Vbc 3V 90º
Vca 3 V 150º
Ejemplo:
Solucion:
480
30º
3
IA 25.83 73.78º
185º1040º
I A I1 I3 480
150º
IB I2 I1 3
IB 25.83166.22º
IC I3 I2 185º1040º
480
90º
3
IC 25.8346.22º
185º1040º
19
Luego:
4800º 185º. 25.83 73.78º I1 3040º 185º. 25.83166.22º 0
I1 14.91 43.78
I2 IB I1 25.83166.22º 14.91 43.78 14.91 163.79
I3 I1 IA 14.91 43.78 25.83 73.78º 14.9176.23
Va'b' 3040º. 14.91 43.78 447.3 3.78
Vb'c' 3040º. 14.91 163.79 447.3 123.79
Vc'a' 3040º. 14.91 76.23 447.3 116.23
20
PARAMETROS Y
MODELOS DE LINEAS DE
TRANSMISION
21
Capitulo II
2.1 Introduccion
22
AAAC: Conductor de aleacion de aluminio.
L
Rdc ρ
A
23
A 60 hz la resistencia Rac es alrededor del 2% mayor que la
resistencia Rdc, como un conductor de hebras es espiralado,
entonces, cada hebra es mas larga que el mismo conductor,
po lo tanto Rdc es mayor. La resistencia del condcutor
aumenta con el incremento de la temperatura y puede
calcularse mediante la siguiente expresion:
T t2
R2 R1.
T t1
Donde:
R1 = resistencia del conductor a un temperatura t1
R2 = resistencia del conductor a un temperatura t2
T = constante de temperatura que depende del material
dφ d
L L (2.1)
dt dt
24
De la ley de Ampere:
H.dl Ix (2.2)
Luego:
2
Ix .I
X
R
En (2.2)
2
X
H.2π.X .I
X
H .I
R 2πR 2
B μH
X
B μ. .I (2.3)
2πR 2
1
W .L.i2 (2.4)
2
25
2
1 μ. X .I .2πxdx
dW
2μrμo
2πR
2
Integrando:
μI 2
W (2.5)
16π
1 μI 2 μ
.L.I 2 Lint (2.6)
2 16π 8π
Para μr 1 entonces o 4 10 7
1
Lint 10 7
2
26
(2.7)
(2.8)
1
Donde: r' r .e 4
1 1
(2.9)
27
(2.10)
1 L11I1 L12I2
2 L22I2 L21I1
Agrupando terminos
1
L11 2 10 7. ln
r'1
1
L22 2 10 7. ln
r'2
1
L12 L21 2 10 7. ln
D
Espaciamiento Simetrico
28
Asumiendo corrientes trifasicas balanceadas:
Ia I b Ic 0
Luego:
1 1 1
a 2 10 7 Ia ln I b ln Ic ln
r'a Dab Dac
Dab Dbc Dca D ; ra rb rc
1 1 1
a 2 10 7 Ia ln I b ln Ic ln (**)
r'a D D
Ia I b Ic
1 1
a 2 10 7 I a ln I a ln
r' D
D
a 2 10 7.I a ln
r'
L 2 10 7 ln
D
r'
D
L 2 10 7 ln
RMG (2.11)
μ D
L ln
2π r'
Por lo tanto la inductancia por fase por Kilometro es :
29
D
L 0.2 ln
Ds (2.11’)
Ds RMG
Espaciamiento Asimetrico
1 1 1
a 2 10 7 Ia ln I b ln Ic ln
r' Dab Dac
1 1 1
b 2 10 7 I b ln Ic ln Ia ln (α)
r' Dbc Dba
1 1 1
c 2 10 7 Ic ln I b ln Ia ln
r' Dcb Dca
L.I
1 1 1
ln ln ln
r' Dab Dac
1 1 1
L 2 10 7 ln ln ln
Dba r' Dbc
1 1 1
ln ln ln
Dac Dbc r '
I b Ia240º a2Ia
Ic Ia120º aIa
Sustituyendo en (α):
30
1 1 1
La 2 10 7 ln a2 ln a ln
r' Dab Dac
1 1 1
Lb 2 10 7 a ln ln a2 ln (2.12)
r' Dbc Dba
1 1 1
Lc 2 10 7 a2 ln a ln ln
r' Dcb Dca
L La Lb Lc
Reemplazando
1 1 1 1
L 2 10 7 3 ln ln ln ln
r' Dab Dac Dcb
1 1
L 2 10 7 ln ln
r' 1
Dab.Dac.Dcb 3
D .D .D 13
L 2 10 ln ab ac cb
7
(2.13)
r'
31
DMG
L 0.2 ln
RMG
1
DMG Dab.Dac.Dcb 3
radio rx radio ry
I I
corriente corriente
n m
(2.14)
32
Como todos los subconductores del arreglo “x” son
electricamente paralelas la inductancia de x sera:
(2.15)
Ejemplo:
6
14 14 6
RMGx
e r.2r.2 3r.4r.2 3r.2r.2r
e r. 2r
49
6
14
1
RMG x 49 e .2.2 3.4.2 3.2.2 r 42.e 4r.64.r 6
RMGx 2.17r
33
Ejemplo:
L 2 10 7 ln
0.608 7
8.025 10
0.011
X L 2π60.8.025 10 7 3.025 10 4
X L 0.3025 km
3 RMGagrup 16
RMG .d.d. 2
c / subconductor
4
(2.20)
34
2.4.9 Inductancia de lineas trifasicas de doble circuito
RMG A 4 RMG .D A 1A 2 2
agrup
2 RMGagrup.D A 1A 2
DMG mH
L 0.2 ln
RMG L km
35
Ejemplo:
IA 30 j50
IB 25 j55 IC IA IB
IC 55 j105
a) Utilizando (α)
1 1 1
N 2 10 7 Ia ln I b ln Ic ln
Dan Dbn Dcn
N 2 10 7 1.504Ia 1.098I b 0.405Ic
N 2 10 7 1.099Ia 0.693I b
N 2.1157 10 5 118.39
Luego:
VN jωN 15km
VN 87.712 j47.402
VN 99.701 151.61
b) Aplicamos (α)
1 1 1
a 2 10 7 Ia ln I b ln Ic ln
r' Dab Dac
a 2 10 7 Ia ln I b ln Ic ln
1 1 1
r' D 2 D
Ic Ia Ib ra rb rc r
36
Va jωa 15km
Va j 314.15 2 10 7 aln
2D
I b ln 2 15km
ra'
Va 404.51 149.67º
q
C
V
q q
D
A 2πxL
D
E
εo
37
q
I
2πεox
D2 C
q
V12 Edx
D1 D1
2πεox
dx
q D
V12 ln 2 (2.23)
2πεo D1
q1 D
Vab q1 ln
2πεo r
q2 D
Vba q1 ln
2πεo r
38
Y como V12=-V21
(2.24)
2πεo
C
D (2.25)
ln
r
0.0556 F
C
D km (2.25)’
ln
r
39
2.5.3 Diferencia de Potencial de 1 conductor empaquetado
q 1 q 2 q 3 ... q n 0
1 n
D
Vij
2πεo
q k ln kj
k 1 Dki
q a1 q b q c 0
40
Asumimos que son sistemas trifasicos, entonces decimos que
despreciando en efecto de la tierra y cubierta de los
cables y asumiendo que la linea de transmision es
transpuesta, obtenemos de (*) que la Diferencia de
Potencial entre la fase a y la fase b para el primer tramo
aplicando (*).
(2.26)
El valor promedio:
(2.27)
2πεo
C
DMG
ln
ragrup
El efecto del agrupamiento es introducir una “r”
equivalente “ra”; el radio equivalente es similar al RMG
41
calculado para la inductancia con la diferencia que se usa
el radio “r” de cada subconductor en lugar del RMG “Ds”,
entonces si “d” es la distancia de separacion de 2
subconductores, el radio del agrupamiento va a ser:
2 subconductores: r A r.d
3 subconductores: r A 3 r.d 2
4 subconductores: r A 1.09 4 r.d 3
2πεo 0.0556
C C
DMG DMG
ln ln
RMGc RMG
RMGc 3 rA.rB.rC
Ejemplo:
DMG 3
35.35.70
0.5328"
RMGL 0.0444'
12
42
1.345"
RMGC 0.056'
2 12
L 0.2 ln
44.097' mH
1.38
0.0444' km
0.0556 μF
C 0.00834
44.097 km
ln
0.056
Ejemplo:
DMG 3
0.9 0.9 1.8 1.134m
43
El efecto de la tierra es incrementar la capacitancia pero
la altura del conductor es mucho mas grande comparada con
la distancia entre conductores y entonces el efecto de la
tierra resultada despreciable. Por tal motivo para todos
los modelos de lineas para analasis de estado estacionario,
el efecto de la tierra en la capacitancia puede ser
despreciable. Sin embargo para analisis del caso
desbalanceado (caso fallas asimetricas) el efecto de la
tierra asi como el de cables recubiertos debe ser
considerado.
X L 2πfL
mH
X L 213.37
km
XL 0.213
km
1.134 100 H
L 2 10 7 ln 1 7.436 10 7
m
2 2.5e 4
mH
L 0.7436
km
X L 2πfL
X L 2π 60 0.7436
X L 0.44
km
44
parametros circuitales apropiados bajo una representacion
por fase.
Z r jωL
SR3φ *
IR
3VR
VE VR IRZ (2.28)
IE IR
A 1
B Z
VE AVR BIR C 0
IE CVR DIR D 1
VE
VR vacio (2.29)
A
45
Carga Capacitiva
SR3φ 3VEIE *
Perdidas totales en la linea
SL3φ SE 3φ SR3φ
Eficacia de la Regulacion
PR3φ
η
PE 3φ
46
Ejemplo:
mH 40km
ZL 0.15 .40km 2πf.1.3263 . 6 j20
km km 1000mH
2200º
VR 1270º
3
SR3φ 380 cos1 0.8 304 j228
SR3φ 304 j228
IL 997.375 36.87º
3VR3φ 3 1270º
VE VR ZIR
VE 1270º6 j20997.375 36.87º 144.2874.917º kV
VE L L 3 144.2874.917º 249.912kV
Y G jωC l
47
Aplicando la primera ley de Kirchoff
Y
IE IR VR. (1)
2
VE Z.IL VR 2 (2)
Reemplazando:
Y Y
IE IR VE. VR.
2 2
IE IR Z.IR Z 1.VR .VE
Y
2
Agrupando terminos
Y
A 1Z
2
B Z
C Y 1 Z
Y Y
D 1Z
4 2
VR A B VE
I C D I
R E
48
Resolviendo la matriz de transmision se puede obtener las
magnitudes reciprocas.
V x x V x ZxI x
V x x V x
ZI x (2.30)
x
V x x V x
x 0 lim
x 0 x
dV x
ZI x (α)
dx
I x x I x V x x .Yx
I x x I x
V x x .Y (2.31)
x
dI x
Si x 0 YV x (β)
dx
d 2V x dI x
2
Z (α’)
dx dx
49
d 2I x dV x
2
Y (β’)
dx dx
d 2V x
ZYV x
dx 2
d 2V x dI x
2
Z
dx dx
(2.32)
d V x
2
ZYV x
dx 2
ZY γ 2
d 2V x
2
γ 2V x 0 (2.32’)
dx
γ = Constante de propagacion
α = constante de atenuacion
β = constante de fase
De la expresion (2.30)
1 dV x
I x
Z dx
1
Z
γA1eγx γA2e γx
I x
γ
Z
A1eγx A2e γx
I x
1
ZC
A1eγx A2e γx
Z
ZC Zc = Impedancia caracteristica
Y
50
VR ZcIR
A1
2
V ZcIR
A1 R
2
Si x=l
Ve A B Vr
Ie C D Ir
Donde:
1
A cosh γl C senhγl
Zc
B Zc.senγl D cosh γl
A partir de (2.34*) podemos encontrar un modelo equivalente
חpara una linea larga, como el modelo de la figura.
Asimismo podemos reemplazar los coeficientes del cuadripolo
parecidas al del modelo de linea media.
Ve 1 Z'Y'.Vr Z'.Ir
(2.35)
Ie Y' 1
Z'Y' Z'Y'
.Vr 1 .Ir
4 2
51
Comparando 2.35 con A y B de la línea media para encontrar
Z’ y Y’ a partir de (2.34)*
cosh γl 1
tanh
γl
2 senhγl
Zsenhγl
Z' Zc.senhγl
γl
γl
tanh
tanh
Y' 1 γl Y 2
2 Zc 2 2 γl
2
1
Velocidad de propagación v
LC
1
Longitud de onda λ , γ jβ
f LC
52
Ve cos βl.Vr j.Zc.senβl.Ir
1
Ie j. senβl.Vr cos βl.Ir
Zc
En regimen de vacio
Ve
Vrvacio
cos βl
Vr
Ir
Zc
Como:
Vno min al
Vr
3
2
kVno min al
SIL
Zc
Para una linea de modelo largo que esta cargado con un SIL
el voltaje y la corriente seran.
53
componentes reactivos, es decir, no tienen potencia
reactiva Qe=Qr=0, esto indica que para una SIL las perdidas
reactivas en la inductancia se compensan con la potencia
reactiva generada por la capacitancia shunt de la linea.
Ejemplo:
a)
rad
β ω LC 2πf LC 2π 60 0.97 0.0115 109 1.259 10 3
km
L 0.97 103
Zc 290
C 0.0115 10 6
1 km
v 299409.25
LC s
b)
5000º
Vr 288.6750º kV
3
800
Sr cos 1 0.8 100036.87 800 j600
0.8
Sr * 1000 36.87
Ir 1154.70 36.87A
3Vr * 3 288.675 0º kV
Ve cos βl.Vr jZc.senβl.Ir
Ve cos 21.77 288.6750ºj290sen 21.77 1154.7
Ve 356.8616.19kV
54
Ve L L 3 356.86 618kV
1
Ie j senβl.Vr cos βl.Ir
Zc
1
Ie j sen 21.77 288.675 cos 21.77 1154.7
290
Ie 900.82 17.66
Luego:
Se 3Ve.Ie* 3 356.8616.19900.8217.66
Se 800 j537MVA
Regulacion de voltaje:
356.86
288.6
cos 21.77
%RV 100%
288.6
%RV 33.15%
Ve A.Vr B.Ir
55
Ve A.Vr
Ir
B
A A α
B B β
Ve Ve δ
Vr Vr 0º
Ve δ A α.Vr 0º
Ir
B β B β
Ve A Vr
Ir δ β α β
B B
2
Vr Ve A Vr
Pr cosδ β cos α β
B B
2
Vr Ve A Vr
Qr senδ β sen α β
B B
PL Pe Pr
QL Qe Qr
56
2
A Vr Ve Vr
Sr (A)
B B
2
D Vr Ve Vr
Se β α β δ (B)
B B
Ve Ve δ Vr Vr 0º
A
Vr β α
2
(2.38)
B
Equivalentemente
A
Vr cos β α
2
Coordenada horizontal del Centro
B
A
Vr sen β α
2
Coordenada vertical del Centro
B
57
El centro del extremo se localiza al trazar el segmento Ocr
en un angulo de (β-α) en direccion positiva respecto al eje
negativo del eje de los MW. Del centro Cr se traza la
circunferencia del extremo receptor con radio:
Ve Vr
B
Ve Vr
B
Para una linea sin perdidas resulta que α=0, B=jX, β=90º
A A0º
58
La potencia activa trifasica transmitida por la linea
resulta:
Ve L L Vr L L
Pr3φ senδ
X'
(2.39)
L L L L L L 2
Ve Vr Vr
Qr3φ cos δ cos βl
X' X'
X' Zc.senβl
senδ
P3φ Ve pu.Vrpu.SIL. (2.40)
senβl
59
Ejemplo:
a)
2π 2π
βl l 350 0.44rad 25.2º
λ 5000
sen36.85
900 1.0 0.9SIL
sen25.2
SIL 709.96MW
2
kV
SIL
Zc
b)
60
Una liea cargada con su impedancia caracteristica no
presenta flujo neto de potencia reactiva y los perfiles de
tension de toda la linea son aproximadamente los mismos. En
lineas largas cargadas con una potencia considerando
potencias menores que la SIL dan como resultados un aumento
del voltaje en la recepcion; y cargas pesadas
considerablementes mayores que la SIL dan como resultado
fuertes caidas de tension en la recepcion.
IR
Vr
+ + Ir
VE VR SR jX LSHUNT
- -
senβl
X LSHUNT Zc
Ve
cos βl
Vr
senβl
X LSHUNT Zc (2.41)
1 cos βl
Y la corriente a mitad de la linea resulta igual a 0.
61
Vr
Vm
βl Im 0
cos
2
Ejemplo:
a)
500
Ve 0º
3
Ve 288.70º
0.95 103
Zc 289.9
0.0113 10 6
βl 24.77º
Ve 288.70º
Vr
cos βl cos 24.77º
Vrnom 550.7kW
b)
sen24.77º
X LSHUNT 289.9 1320
1 cos 24.77º
kV 2 5002
Q3φ 189.39MVAR
X LSHUNT 1320
62
Se usa la compensacion capacitiva shunt para retrasar el
factor de potencia de lineas cargadas pesadamente. El
efecto es proveer la capacitancia necesaria para mantener
Vr en nivel adecuado. Los capacitores se conectan
directamente a una barra o al enrollado terciario del
transformador y tambien puede estar dispuestos a lo largo
de la linea para minimizar perdidas y caidas de voltaje.
Ve.Vd.senδ
P Potencia transferida en la linea
X'XcSERIE
Z' j X'Xcserie
Xcserie
Porcentaje de compensacion (25-75%)
X'
63
a)
L
Zc 289.9
C
X' Zc.senβl
X' 289.9 sen 24.76º
X' 121.41
P 800MW
500.500
800 senδ
121.41
δ 22.86º
500.500 5002
Qr cos 22.86º cos 24.76º
121.41 121.41
Qr 22.57MVAR
2
V 5002
Xc 436.75
Scap j572.4
1
Xc
2πfC
1
C 6.07μF
2π60 436.75
b)
Z' j X'Xserie
Z' j X'X serie 40% j 121.41 48.58 j72.8
tan j tan
Z βl 2 24.76
Y' j
Zc 2 289.9 2
Y' j0.0015
64
Z'Y' j72.80.0015
A 1 1 j
2 2
A 0.945
5000º
Vr 288.67kV
3
Sr * 1000 36.87
Ir 1.153 36.87kA
3Vr * 3 288.670º
Ve A.Vr B.Ir
Ve 0.945 288.67 j 72.8 1.155 36.87º
Ve 330.17 11 .76kV
L L
Ve 571.8711 .76kV
Regulacion
571.87
500
%RV 0. 945 100% 21.0%
500
65
ANALISIS DE FLUJO DE
POTENCIA
66
Capitulo III
67
3.1 Ecuaciones de un sistema de dos nodos
2
S1 PG 1 PD1 j QG 1 PD1 ZL ZL θ
V1
XL
V2 V2 δ2
Si SGi SDi
Si ViIi *
La corriente Ii tiene 2 componentes, la que fluye por la
resistencia o impedancia de linea serie, y la admitancia en
paralelo.
S1 *
I1 V1Ysh V1 V2 Ys
V1 *
(3.1)
S2 *
I2 V2Ysh V2 V1 Ys
V2 *
Ecuaciones generales de flujo de carga o de flujo de
potencia
I1 Y11V1 Y12V2
68
(3.2)
I2 Y21V1 Y22V2
Y11 YSHUNT YS
Y22 YSHUNT YS
1
i) YSH j Yshunt es capacitivo(*)
2XC
R
ii) Se define un factor de perdidas α
XL
π R
j
iii) Z
serie R jXL X Le
2 XL
1 j α 2
π
iv) La admitancia en serie YS e (**)
XL
De la expresion (3.1)
1 1 j α 2
π
S1* jV1 V1 V2 e V1 *
2X C X L
1 1 j α 2
π
S2* jV2 V2 V1 e V2 *
2X C X L
Simplificando
69
π
j α
2 2 π δ2 δ1 2
V1 V j α V1 V2 e e
S1* j 1 e 2
2XC XL XL
π
j α
V1 V2 e δ1 δ2 e
2 2
V2 V2 j α π2 2
S2* j e
2XC XL XL
2 2
V V π π
P1 jQ1 1 90º 1 cos α 2 jsenα 2
2XC XL
V1 V2
cosδ1 δ2 jsenδ2 δ1 . cosα π .senα π
XL 2 2
2
V V V
P1 PG 1 PD1 1 senα 1 2 sen α δ1 δ2
XL XL
(3.4)
2
V2 V V
P2 PG 2 PD 2 senα 1 2 sen α δ2 δ1
XL XL
2 2
V1 V V V
Q1 QG 1 QD1 cos α 1 2 cos α δ1 δ2 1
XL XL 2XC
(3.5)
2 2
V2 V V V
Q1 QG2 QD 2 cos α 1 2 cos α δ2 δ1 2
XL XL 2XC
QG 1 QG 2 QD1 QD 2
cos α
ZL
2 2
V1 V2 2V1 V2 cosδ1 δ2
(3.7)
V1 2
V1 .YSH
2
70
De (3.6) y (3.7)
71
P1 PD1
P P
2 D2
P Pn PDn
Pn 1 PDn 1 QD1
P P Q
2n D 2n Dn
ii) Variable de control o independientes
U1 PG 1
U P
2 G2
U U n PGn
U n 1 QG 1
U Q
2n Gn
X1 δ1
X δ
2 2
X X n δn
X n 1 V1
X V
2n n
i Pi,Qi
En este tiepo de barra se conoce previamente
la potencia ctiva y reactiva neta de inyeccion
en la barra. Las incognitas son siempre el
72
modulo y el argumento de la tension. Son barras que no
poseen generadores por lo tanto, no generan potencia activa
aunque si pueden llegar a generar potencia reactiva
mediante banco de condensadores, maquinas sincronas o
dispositivos electronicos.
Vimin Vi Vimax
73
La diferencia entre angulos de tensiones nodales conectados
por una linea no debe sobrepasar cierto limite, esto para
evitar sobrecargas de lineas y transformadores, asi como
favorecer las estabilidad del sistema.
δi δj δi δj max
ii) Variables de control [U]
74
Aplicando la primera ley de Kirchoff
I4 Y43 V4 V3
Admitancias
Uniformizando
0 Y34V3 Y34V4
Luego:
75
Y44 Y34 Y34 Y43 y 34
Generalizando
IBARRA YBARRA.VBARRA
Donde:
Entonces:
1
VBARRA YBARRA.I BARRA
76
1
ZBARRA YBARRA
f Xº , Uº , Pº 0
77
determinada configuracion nominal de cargas Pº se
especifican exactamente “2n” componentes y de los dos
vectores Xº y Uº, luego se resuelven los restantes “2n”
componentes no especificados usando las ecuaciones SLFE.
78
a) Se desprecian las resistencias de lineas
79
Ejemplo
SG1
SG2 XL 0.959
km
1 2 V 220kV
PG 2 100MW
SL 400MVA
l 55km
SD1 100 j50MVA
SD 2 100 j100MVA
SB 100MVA
3
VB 220kV
SG3 SD3
VB2 2202
ZB 484
SB 100
XL 0.959 55
ZLpu l j j0.109
ZB 484
YLpu j9.17
Valores Iniciales
Para la Barra 1
80
P1 jQ1 Y11V1V1 * Y12V2V1 * Y13V3V1 *
PG1 PD1 j QG 1 QD1 j9.1762 j9.1762 1.00º 1.0 0º j9.1762 1.0δ2 1.0
j9.1762 1.0δ3 1.0 0º
PG1 1.0 9.1762senδ2 9.1762senδ3
QG1 0.5 18.3525 9.1762 cos δ2 9.1762 cos δ3
Para la Barra 2
Para la Barra 3
Luego:
δ2 2.0813º
δ3 2.0813º
QG1 0.5122
QG 2 0.0304
QG3 1.0304
81
V1 pu V2pu
P12pu sen δ2 0.3332
X pu
V2pu V3 pu
P23pu senδ2 δ3 0.6659
X pu
V1 pu V3 pu
P13pu senδ1 δ3 0.3332
X pu
82
3.6. Aspectos computacionles de flujos de carga
Pi jQi
Yi1V 1 Y12 ... Y1nVn
Vi *
3.7 Metodos Numerico Iterativos para la solucion de
ecuaciones algebraicas no lineales
F(x)=0 (3.9)
F1(x1,x2,x3,…,xn)=0
F2(x1,x2,x3,…,xn)=0
F3(x1,x2,x3,…,xn)=0
83
El 2do resultado se usa para encontrar para encontrar el
3er estimado.
f(x) x 2 5x 4 0 y x 2 5x 4
x1 4
x2 1
1 2 4
x x
5 5
xu 1 g(x v )
x 3
Iteracion 1
x 1 g x
84
1 2 4 13
x 1 g 3 3
5 5 5
Iteracion 2
x 2 g 2.15
13
5
Y asi sucesivamente
xu 1 x v ζ
Ejemplo
2x1 x1x2 1 0
2x2 x1x2 1 0
x1x2
x1 0.5
2 x1 1
x1x2 x2 1
x2 0.5
2
Por Gauss
Iteracion “0”
x(10) 0
x(20) 0
Iteracion “1”
Iteracion “2”
0.5 0.5
x(12) 0.5 0.625
2
0.5 0.5
x(22) 0.5 0.625
2
85
Iteracion “3”
0.625 0.625
x(13) 0.5 0.695
2
0.625 0.625
x(23) 0.5 0.695
2
Iteracion “0”
x(10) 0
x(20) 0
Iteracion “1”
Iteracion “2”
0.5 0.5
x(12) 0.5 0.625
2
0.625 0.5
x(22) 0.5 0.656
2
Iteracion “3”
0.625 0.656
x(13) 0.5 0.705
2
x(23) 0.5
0.705 0.656 0.731
2
f1 x1, x2 0
f2 x1, x2 0
86
El metodo puede ser extendido a sistemas de n ecuaciones,
nosotros asumimos lo valores inciales.
x1(0),x2(0)
x1(0), x2(0)
x1(0) x1(0)
Son las soluciones buscadas
x2(0) x2(0)
Por lo que estas estas soluciones buscadas deben satisfacer
el sistema de ecuaciones.
f1 x1(0) x1(0), x2(0) x2(0) 0
f x
2 1
(0)
0
x1(0), x2(0) x2(0)
x2 (0)
(0)
x1 (0)
2
x2 2 x1(0) ... 0
x2 (0)
(0)
x1 (0)
2
x2 2 x1(0) ... 0
2f
Los terminos orden superior son despreciables.
x 2
En forma matricial
F(0) J(0)X(0) 0
87
f1 f1 f1
f1 x(0) X1(0) x1 (0) x1 (0) x1 (0)
(0) (0) f1 f1
f x X f1
2 2 J x1
F(0) X(0) x1 (0) x1 (0)
(0)
(0) (0)
fn x
X n f f1 f1
1
x1 (0) x1 (0) x1 (0)
X(0) J 1f(0)
88
Capitulo IV
Solución del flujo de potencia
4.1 Introducción:
89
Pi jQi ViIi*
Pi jQi (*)
Ii
Vi*
Pi jQi n n
Vi yin yijVj) (4.2)
Vi* i0 j1
j 0
n
(k)
Pi ReVi (Vi i0 yin yijVj )
(k 1) (k)* (k) n
j 1
(4.4)
n
(k)
Qi ImVi (Vi j1 yin yijVj )
(k 1) (k)* (k) n
j 1
90
En la expresión (4.3) merced a que los elementos fuera de
la diagonal de la matriz admitancia son (Yij=-yij) el
negativo d la admitancia de ambas barras,(4.3) se pude
escribir como:
n
Vi (Vi Yii
(k)* (k)
YijVj )
j1
(k)
Vi(k 1)
Yii
Por lo tanto (4.4) se pude escribir como:
n
Pi(k 1) ReVi(k)*(Vi(k).Yii YijVj(k))
j1
(4.4)’
j i
n
Qi(k 1) ImVi(k)*(Vi(k).Yii YijVj(k))
j1
ji
| ViR(k 1) ViR(k) | ξ
| ViI(k 1) ViI(k) | ξ
91
Ejemplo:
P2prog jQ2prog
(k)*
y12V1(k) y 23V3(k)
V2
V 2(k 1)
y12 y 23
P3 prog jQ3 prog
(k)*
y13V1(k) y 23V2(k)
V3
V 3(k 1)
y13 y 23
N V2 E V3 E
0 1 ∟0º -- 1 ∟0º --
92
3 0.9591∟-3.7850 0.0223 0.9783∟-3.3231 0.0115
Flujo de potencia:
S12 V1.I12*
S21 V2.I21*
S13 V1.I13*
S31 V3.I31*
S23 V2.I23*
S32 V3.I32*
93
S12 1.8640 j0.8903
S21 1.7833 j0.6892
S13 2.2604 j1.6153
S31 2.1881 j1.3609
S23 0.3725 j0.1541
S32 0.3760 j0.1498
Ejemplo 2:
94
Sbase 100MVA
S2 2.5 1.09
P3 1.2
Valores iniciales:
Hallando V2
P2prog jQ2prog
(k)*
y12V1(k) y 23V3(k)
V2
V 2(k 1)
y12 y 23
P3 prog jQ3 prog
(k)* y13V1(k) y 23V2(k)
V3
V 3(k 1)
y13 y 23
Hallando Q3:
Q3 Im V (V .(y y ) y V1 y V2 )
(1) (0)* (0) (0) (1)
3 3 13 23 13 23
Hallando V3:
95
P3prog jQ3prog
(0)*
y13V1(0) y23V2(1)
V3
V 3(1)
y13 y 23
N V2 Q3 V3
0 1 +j0º -- 1.02 +j0º
1 0.9831-j0.0327 0.2481 1.02-j0.0054
2 0.9816-j0.0354 0.3771 1.02-j0.0079
3 0.9815-j0.0369 0.4118 1.02-j0.0090
4 0.9815-j0.0376 0.4277 1.02-j0.0095
5 0.9815-j0.0379 0.4366 1.02-j0.0097
6 0.9815-j0.0380 0.4377 1.02-j0.0098
7 0.9815-j0.0380 0.4390 1.02-j0.0098
V2pu=0.9815-j0.0380;Q3pu=0.4390;V3pu=1.02-j0.0098
Luego :
96
Las Corrientes :
Flujo de potencia :
S12 V1.I12*
S21 V2.I21*
S13 V1.I13*
S31 V3.I31*
S23 V2.I23*
S32 V3.I32*
97
El flujo de potencia reactiva en una línea está determinada
por la diferencia angular de los voltajes en sus
terminales, mientras que el flujo de potencia reactiva
es originado principalmente por la diferencia de los
módulos de voltajes en sus terminales.
1
Vx Vj
a
Ii a*Ij
98
Aplicando la primera le de Kirchhoff
1
Ii (Vi Vx)Y YVi YVj
a
Ii
Ij
a*
1 1 1
Ij * [Y(Vi Vx)] * [Y(Vi Vj)]
a a a
1 1
Ij * YVi YVj
a | a |2
1
Y Y Vi
Ii a
Ij 1 1 Vj
* Y Y
a | a |2
99
SOLUCION DEL FLUJO
DE POTENCIA NEWTON
RAPHSON
100
Capitulo V:
n n
Ii Vi yij yijVj (5.1)
j0 j1
j1
n
Ii YijVj
j1
n
(5.2)
Ii | Yij || Vj | (θij δj)
j1
Pi jQi Vi*Ii
n
(5.3)
Pi jQi | Vi | (δi) | Yij || Vj |(θij δj)
j1
101
n
Pi | Yij || Vj | | Vi | cos(θij δi δj)
j 1
n
(5.3)’
Qi | Yij || Vj | | Vi | sen(θij δi δj)
j1
n
Pi | Vi |2 Gii | Vi | | Vj || Yij | cos(θij
j1
δi δj)
ji
n
Pi [| Vi |2 Bii | Vi | | Vj || Yij | cos(θij δi
j1
δj)]
ji
------(5.3’’)
102
P2 P P2 P2
P2(0) δ2(0) 2 δ4(0) | V2 |(0) | V4 |(0) 0
δ2 (0) δ4 (0) | V2 |(0) | V4 |(0)
P3 P P3 P3
P3(0) δ3(0) 3 δ4(0) | V3 |(0) | V4 |(0) 0
δ3 (0) δ4 (0) | V3 |(0) | V4 |(0)
P4(0)...
Q2(0)...
Q3(0)...
103
Ejemplo:
S2=-2.5-j1.09
P3=1.2
V1=1.03∟0º
V2=1.0+j0.0
V3=1.02+j0.0
Solución:
104
Hallando la matriz Jacobiana:
P 2
V 2 V 1 Y 12 cos(θ12 δ2 δ1) V 2 V 3 Y 23sen(θ23 δ2 δ3)
δ2
P 2
V 2 V 3 Y 23sen(θ23 δ2 δ3)
δ3
P 2
V 1 Y 12 cos(θ12 δ2 δ1) 2V 2 Y 22 cos(θ22) V 3 Y 23 cos(θ23 δ2 δ3)
V2
P 2
47.7586
δ2
P 2
30.0000
δ3
P 2
23.9838
V 2
P3
V 3 V 2 Y 23sen(θ23 δ3 δ2)
δ2
P3
V 3 V 1 Y 13sen(θ13 δ3 δ1) V 3 V 2 Y 23sen(θ23 δ3 δ2)
δ3
P3
V 3 Y 23 cos(θ23 δ3 δ2)
V 2
P 3
30.0000
δ2
P 3
57.7517
δ3
P 3
18.0000
V2
Q2
V 2 V 1 Y 12 cos(θ12 δ2 δ1) V 2 V 3 Y 23 cos(θ23 δ2 δ3)
δ2
Q2
V 2 V 3 Y 23 cos(θ23 δ2 δ3)
δ3
Q2
V 1 Y 12sen(θ12 δ2 δ1) 2V 2 Y 22sen(θ22) V 3 Y 23sen(θ23 δ2 δ3)
V2
Q2
25.1034
δ2
Q2
18.0000
δ3
Q2
45.5477
V 2
105
47.7586 30.0000 23.9838
J 30.0000 57.7517 18.0000
25.1034 18.0000 45.5477
0.0273 0.0166 0.0078
J 1
0.0168 0.0256 0.0013
0.0084 0.0010 0.0171
P 2 δ2
P3 J δ3
Q2 V 2
δ2 P 2
1
δ3 J P3
V 2 Q2
Valores iniciales:
106
δ2(1) 0.0378rad
δ3(1) 0.0090rad
(1)
V 2 0.09831rad
P 2(2) 2.4544
P3(2) 1.1929
Q2(2) 1.0664
P 2 P 2 P 2(2)
P3 P3 P3(2)
Q2 Q2 Q2(2)
P 2 0.0456
P 3 0.0071
Q2 0.0246
P 2 δ2
P3 J δ3
Q2 V 2
δ2 P 2
1
δ3 J P3
V 2 Q2
Luego:
V2=0.9815-j0.0381;Q3=0.44;V3=1.02-j0.0099
107
I12 y12(V1 V2)
I21 I12
I13 y13(V1 V3)
I31 I13
I23 y 23(V2 V3)
I32 I23
S12 V1.I12*
S21 V2.I21*
S13 V1.I13*
S31 V3.I31*
S23 V2.I23*
S32 V3.I32*
108
J 1 0
J
0 J 2
Q 0 J2V
O la última ecuación matricial se pude desacoplar en 2
ecuaciones.
P
P J 1δ δ
δ
Q
(5.6)
Q J 4 V V
V
Pi n
Vi Vj Yijsen(θij δj δi)
2
Vi Yijsen(θii)
δi j1
Pi 2
Qi Vi Yii sen(θii)
δi
Pi 2
Qi Vi Bii
δi
Pi
Vi Bii (5.8)
δi
109
En condiciones de operación normales δi- δj es bastante
pequeño por lo que la ecuación (5.3)’’ θii- δi+δj
resulta θii, es decir
Pi
Vi VjBij (5.9)
δi
Pi
Vi Bij (5.9’)
δi
Qi n
Vi
ViYiisen(θii) Vi Vj Yijsen(θij δj δi)
j1
Qi
ViYiisen(θii) Qi
Vi
Luego Bii>>Qii
Qi
Vi Bii
Vi
Qi
Vi Bij Para los electos fuera de la diagonal (5.11)
Vi
P
B' δ (5.12)
Vi
Q
B'' Vi (5.13)
Vi
110
DESPACHO ÓPTIMO DE
GENERACION
111
Capitulo VI
112
derivada de F respecto a las variables de equivalencia a
cero.
f
0
xi
i 1,2,3,...n
f f f f
f , , ,...,
x1 x2 x3 xn
2 f
H
xi xi j
Ejemplo:
2 2 2
f ( x1, x 2, x3,..., xn) x1 2 x2 3 x3 x1 x2 x2 x3 8 x1 16 x2 32 x3 110
f
2x1 x2 8 0
x1
f 2 1 0 x1 8
4x2 x1 x3 16 0 1 4 1 x2 16
x2
f 0 1 6 x3 32
6x3 x2 32 0
x3
f(x1, x2, x3) f(3,2,5) 2
113
El 2 representa el mínimo de la función F para este caso el
mínimo local es igual al mínimo global.
2 f 2 f 2 f
2
x1 x1 x2 x1 x3 2 1 0
2 f 2 f 2 f
H 1 4 1
x22 x1 x2
2
x2 x3
0 1 6
f 2 f 2 f
x X x3 X 2 x3
2
3 1
La matriz hessiana es positiva por lo tanto el vector
a) Restricciones de Igualdad:
k
L f igi
i 1
L f k
gi
i 0
xi xi i 1 xi
L
gi 0
i
114
Ejemplo:
( x 8) 2 ( y 6) 6 25
f ( x, y ) x 2 y 2
L x 2 y 2 [( x 8) 2 ( y 6) 6 25]
L
2 x [2( x 8)] 2 x 2x 16 0
x
2 x( 1) 16 (1)
L
2 y [ 2( y 6)] 2 y 2y 12 0
y
2 y ( 1) 12 ( 2)
L
( x 8) 2 ( y 6) 6 25 (3)
x x 8 0
y y 6 0
x 4, y 3
x 12, y 9
(4,3) 1
(12,9) 3
115
6.13 Optimización de parámetros restringidos: Restricciones
de desigualdad
Ejemplo:
Solución:
L x 2 y 2 λ[(x 8)2 (y 6)2 25] μ(2x y 12)
L
2x 2λx 16λ 2μ 0 (A)
x
L
2y 2λy 12λ μ 0 (B)
y
L
(x 8)2 (y 6)2 25 0 (C)
λ
L
2x y 12 0 (D)
μ
116
Resolviendo el sistema:
2 x 4 y 2x 4y 8 0 ( E )
y 12 2 x ( F )
4 4 .8
x (*)
1
4 2 .4
y (**)
1
(*)y(**) en (C)
2 2
4 4.8 4 2.4
8 6 25 0
1 1
2 0.36 0
2
1 1.8
2 0.2
117
6.2 Costo Operativo de una central térmica:
118
Usualmente en la práctica el costo del combustible del
generado ”i” s representa como una función cuadrática de
generación de potencia real.
Ct C1 ( P1 ) C2 ( P2 ) ... Cn ( Pn )[$ / h]
Pt=P1+P2+…+Pt
119
pequeños para 2 a 10min, por lo tanto, el problema de
despacho económico se pude plantear de la siguiente manera:
Nuestro problema es determinar los valores de las salidas
de las unidades P1, P2, P3, P4,…Pn, que minimice el costo
total dado por (6.1) sujeto a las restricciones de igualdad
dadas por (6.2)
dC1 dC dC
dP1 2 dP2 .... n dPn 0 (6.7)
dP1 dP2 dPn
120
6.4 Solución Del Despacho Económico Si Considerar Los
Límites Del Generador Ni La Pérdidas de línea.
Ejercicio:
2
C1 8 P1 10.10 3 P1
2
C2 7 P2 8.10 3 P2
Solución:
dC1
8 20.10 3 P1
dP1
dC2
7 16.10 3 P2
dP2
dC1 dC2
8 20.103 P1 7 16.103 P2 (1)
dP1 dP2
121
Pt P1 P2 (2)
1 20.10 3 P1
P2
16.10 3
1 20.10 3 P1
Pt P1
16.10 3
1 5
Pt 3
P1 1
16.10 4
4
P1 Pt 27.28MW
9
dC1 dC 2 4
8 20.10 3 Pt 27.28
dP1 dP2 9
dCi
7.44 0.0089 Pt ($ / MW h)
dPi
Costo de operación total:
Pt P1 P2 dCi/dPi Ct
500 194.4 305.6 11.89 4041
600 238.9 361.1 12.78 5197
700 283.3 416.7 13.67 7335
800 327.8 472.2 14.56 7775
900 372.2 527.8 15.45 9197
1000 416.6 583.4 16.34 11875
1100 461.1 638.0 17.23 12308
1200 505.6 694.4 18.12 13998
1300 550.0 750.0 19.01 15775
1400 594.4 805.6 19.90 17641
1500 638.9 861.1 20.79 19597
Pi min Pi Pi max
i 1,2,3,..., n
122
Cuando se incluyen restricciones de desigualdad, la
solución del flujo de potencia se modifica de la siguiente
manera.
100≤P1≤600 (MW)
400≤P2≤1000 (MW)
dC1
13.4 8 20.10 3 P1
dP1
P1 270 MW
270 400 670 MW
dC1
8 20.10 3 (600) 20$ / Mw h
dP1
dC2
20 7 16.10 3 P2
dP2
P2 812.5Mw
Pt 812.5 600 1412.5MW
123
670 Pt 1412.5
Pt P1 P2 dCi/dPi
500 100.0 400.0 10
600 200.0 400.0 12
670 270.0 400.0 13.4
700 283.3 416.7 13.67
800 327.8 472.2 14.56
900 372.2 527.7 15.45
1000 412.5 583.4 16.34
1100 461.1 638.9 17.23
1200 505.6 694.4 18.12
1300 549.9 750.1 19.01
1400 594.4 805.6 19.90
1412.5 600.0 812.5 20.00
1500 600.0 900.0 21.40
1600 600.0 1000 23.00
Pt=P1+P2+…+Pn-Pl (6.9)
P P P
( dP1 dP2 ... dPn ) L dP1 L dP2 ... L dPn 0
P1 P2 Pn
Multiplicando por
124
P P P
dP1 dP2 ... dPn L dP1 L dP2 ... L dPn 0
P1 P2 Pn
(6.10)
C1 P C P C P
L dP1 2 L dP2 ... n L dPn 0
P1 P1 P2 P2 Pn Pn
(6.11)
dCi P
L 0
dPi P1
dCi 1
dPi 1 PL
P1
P 1,2,..., n
dCi
Li (6.12)
dPi
n
Ct Ci
i 1
n
Ct i iPi iPi 2 (A)
i 1
125
Esta ecuación está sometida a la siguiente restricción de
igualdad.
Pi PD ----------(B)
i 1
De la primera condición:
k n
L Ct PD Pi 0 (C)
i 1 i 1
Entonces:
dL
0 1º condición
dPi
dL
0 2º condición
d
dL Ct Ct
(0 1) 0
dPi Pi Pi
Ct C1 C2 .... Ct
Ct Ci
Pi Pi
De la segunda condición:
Ci
Pi
i 1,2,3,...n (D)
i 2iPi
Pi PD (B)
i 1
i
Pi (E)
2i
dCi
para Pi min Pi Pi max
dPi
126
dCi
para Pi Pi max
dPi
dCi
para Pi Pi min
dPi
n
i
i 1 2i
PD (F)
n
i
PD
2i
n
i 1
(G)
1
i 1 2i
f ( ) PD
(k )
df ( )
f ( ) ( k ) ( k ) PD
d
PD f ( ) ( k ) (I)
( k ) (k )
PD ( K )
df ( )
d
P ( k )
( k ) (k )
n
dPi (J)
i 1 d
127
P ( k )
(k )
n
1 (K)
i 1 2i
( k 1) ( k ) ( k )
n
P ( k ) PD PD ( k ) PD Pi ( k )
i 1
dCi
para Pi min Pi Pi max
dPi
dCi
para Pi Pi max
dPi
dCi
para Pi Pi min
dPi
Ejemplo:
2
C1 500 5.3P1 0.004 P1
2
C2 400 5.5 P2 0.006 P2
2
C3 200 5.8 P3 0.009 P3
200 P1 450
150 P 2 350
100 P3 225
PD 975MW
128
Solución:
Asumiendo
( 0 ) 6
Usando la ecuación (E)(ecuación de coordinación)
6 5 .3
P1 87.5MW
2(0.004)
6 5 .5
P2 41.6667 MW
2(0.006)
6 5.8
P3 11 .1111 MW
2(0.009)
P (1) PD ( P1 P 2 P3 )
P (1) 834.7222MW
834.7222
(1) 3.1632
1 1 1
2(0.004) 2(0.006) 2(0.009)
( 2 ) (1) (1)
( 2 ) 6 3.1632
( 2 ) 9.1632
9.1632 5.3
P1 482.8947MW
2(0.004)
9.1632 5.5
P2 305.2632MW
2(0.006)
9.1632 5.8
P3 186.8421MW
2(0.009)
P (2) 975 (482.8947 305.2632 186.8421)
P (1) 0 MW
129
P (2) 975 (450 305.2632 186.8421)
P (1) 32.8947
32.8947
( 2 ) 0.2368
1 1
2(0.006) 2(0.009)
( 3) ( 2 ) ( 2 )
( 3) 9.4
9.4 5.5
P2 325MW
2(0.006)
9.4 5.8
P3 200 MW
2(0.009)
P1 450 MW
P 2 325MW
P3 200 MW
Luego el costo total de operación
Ct=8236.25$/h
Problema:
Para el ejercicio 6.5.1 las pérdidas de transmisión totales
del sistema en cuestión están dadas por:
2
PL 1.5.10 4 P1 2.10 5 P1 P2 3.105 P22 MW
dCi 1
dPi 1 PL
P1
130
8 20.10 3 P1
16 (1)
1 (3.10 4 P1 2.10 3 P2 )
7 16.10 3 P2
16 (1)
1 ( 2.10 5 P1 6.10 5 P2 )
P1 315.8103MW
P2 524.7017 MW
Pérdidas de transmisión
PD P1 P2 PL 813.972 MW
CT=9399.25$~/MW-h
dCi 1
dPi 1 PL
P1
Cuando se incluyen las pérdidas
131