Crite Rio Igual Dad Are As

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 15

CONTROL Y ESTABILIDAD DE

SISTEMAS ELÉCTRICOS DE
POTENCIA CON MATLAB
CRITERIO DE IGUALDAD DE ÁREAS

𝑉̇ − 𝑉̇
𝑆 = 𝑉̇ 𝐼 ∗̇ = 𝑉̇
POTENCIA ENTREGADA 𝑗𝑥

𝑉 𝑐𝑜𝑠𝛿 + 𝑗𝑉 𝑠𝑒𝑛𝑜𝛿 − 𝑉
Generador – Barra Infinita 𝑆 = 𝑉̇
𝑗𝑥
𝑉 𝑠𝑒𝑛𝑜𝛿 + 𝑗 𝑉 𝑐𝑜𝑠𝛿 − 𝑉
𝑆 = 𝑉̇
𝑥
𝑉 𝑠𝑒𝑛𝑜𝛿 + 𝑗 𝑉 𝑐𝑜𝑠𝛿 − 𝑉
𝑉 𝑐𝑜𝑠𝛿 + 𝑗𝑉 𝑠𝑒𝑛𝑜𝛿
𝑥
𝑉 𝑐𝑜𝑠𝛿𝑠𝑒𝑛𝑜𝛿 − 𝑗𝑉 𝑐𝑜𝑠𝛿𝑠𝑒𝑛𝑜𝛿 + 𝑉 𝑉 𝑠𝑒𝑛𝑜𝛿
𝑉𝑉 𝑥
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
𝑥
𝑉 𝑐𝑜𝑠 𝛿 − 𝑉 𝑉 𝑐𝑜𝑠𝛿 + 𝑉 𝑠𝑒𝑛𝑜 𝛿
+
𝑥
CURVA POTENCIA -
ÁNGULO Potencia - Angulo de Decalaje
12

𝑃 =𝑃 𝑠𝑒𝑛𝑜 𝛿 Pmax
10

𝐸𝑛 𝑒𝑞𝑢𝑖𝑙𝑖𝑏𝑟𝑖𝑜: Pm

P [MW]
6
𝑃 = 𝑃 = 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
4

𝛿=𝛿
2

0
0 20 40 60 80 90⁰100 120 140 160 180
δo Δ
ECUACIÓN DE OSCILACIÓN
Tm Te
𝑑𝜔 GAS/AGUA
TURBINA
𝐽
𝑑𝑡
= 𝑇 −𝑇 G

• J, momento de inercia Generador- GOBERNADOR PL


Turbina

• 𝜔 , velocidad angular del rotor


2𝐻 𝑑𝜔
= 𝑃 − 𝑃 (1)
𝜔 𝑑𝑡
La constante de Inercia esta definida.
1 𝐽𝜔
2𝐻 𝑑 𝛿
𝐻= = 𝑃 − 𝑃 (2)
2 𝑆 𝜔 𝑑𝑡
2𝐻𝑆 𝑑𝜔 • 𝑃 , potencia mecánica en p.u.
=𝑇 −𝑇
𝜔 𝑑𝑡
2𝐻 𝑑𝜔 𝑇 −𝑇 𝜔 • 𝑃 , potencia eléctrica en p.u.
=
𝜔 𝑑𝑡 𝑆
• Recordando:
DIAGRAMAS DE BLOQUE 𝑑
DE MATLAB ℒ 𝑓(𝑡) = 𝑠𝐹 𝑠 − 𝑓(0)
𝑑𝑡
De la ecuación (1):
𝑑 1
∆𝜔 = 𝑠 ∆𝜔 (𝑠) = 𝑃 −𝑃
𝑑𝑡 2𝐻
2𝐻 𝑑𝜔
=𝑃 −𝑃 1
𝜔 𝑑𝑡 ∆𝜔 (𝑠) = 𝑃 −𝑃
𝑑 𝜔 1
2𝐻 𝑠
= 𝑃 −𝑃
𝑑𝑡 𝜔 2𝐻
𝑑 1
∆𝜔 = 𝑃 −𝑃 𝑓( 𝑃 sin 𝛿)
𝑑𝑡 2𝐻
∆𝜔 , variación de velocidad del rotor 𝑃 𝑑∆𝜔 ∆𝜔 𝑑𝛿
𝑃 - 𝛿
𝑑𝛿 1 𝑑𝑡 1 𝜔 𝑑𝑡 1
= 𝜔 = ∆𝜔 𝜔 + 2H S
𝑑𝑡 S
CAMBIO EN POTENCIA
MECÁNICA
Potencia - Angulo de Decalaje
1.2

Si en un momento dado se produjese un


incremento de la potencia mecánica 1

entrante al generador desde Pm0 = 0,5 𝑃 0.8

hasta Pm1 = 0,8. La relación en la curva

PE [PU]
0.6
Potencia – Ángulo permite interpretar la 𝑃
variación de potencia eléctrica y 0.4

ángulo. 0.2

𝛿𝑜 𝛿 𝛿
0
0 20 40 60 80 100 120 140 160 180
Δ

𝑃𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝐸𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑎 𝑒𝑛 𝑝. 𝑢.
CAMBIO EN POTENCIA
MECÁNICA

La modelación del sistema en formato


SIMULINK, reproduce las oscilaciones
de potencia, ángulo y variación de
velocidad.

𝑃𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝐸𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑎 𝑒𝑛 𝑝. 𝑢.
• Integrando entre dos puntos

CRITERIO DE IGUALDAD 𝑑𝛿 𝑑𝛿 𝜔 𝑃 −𝑃
DE ÁREAS | − | = 𝑑𝛿
𝑑𝑡 𝑑𝑡 𝐻
Considerando la ecuación (2): Considerando puntos de equilibrio
2𝐻 𝑑 𝛿
= 𝑃 −𝑃 𝜔 𝑃 −𝑃
𝜔 𝑑𝑡 𝑑𝛿 = 0
𝐻
𝑑 𝛿 𝜔 𝑑𝛿
= 𝑃 −𝑃 |2
𝑑𝑡 2𝐻 𝑑𝑡 Integrando por partes.
𝑑𝛿 𝑑 𝛿 𝜔 𝑑𝛿
2 = 𝑃 −𝑃 2
𝑑𝑡 𝑑𝑡 2𝐻 𝑑𝑡 𝜔 𝑃 −𝑃 𝜔 𝑃 −𝑃
𝑑𝛿 + 𝑑𝛿 = 0
𝑑 𝑑𝛿 𝜔 𝑃 − 𝑃 𝑑𝛿 𝐻 𝐻
=
𝑑𝑡 𝑑𝑡 𝐻 𝑑𝑡
𝜔 𝑃 −𝑃 𝜔 𝑃 −𝑃
𝑑𝛿 = 𝑑𝛿
𝐻 𝐻
CRITERIO DE IGUALDAD
DE ÁREAS
Potencia - Angulo de Decalaje
1.2

𝜔 𝑃 −𝑃
𝐴 = 𝑑𝛿 1

𝐻 𝑨𝟐
0.8
𝑃
𝜔 𝑃 −𝑃 𝑨𝟏
𝐴 = 𝑑𝛿

PE [PU]
𝐻 0.6

𝑃
Considerando la igualdad de áreas la 0.4

energía ganada es igual a la energía 0.2

perdida.
0
0 20 40 60 80 100 120 140 160 180
• El sistema es estable si el área 2 es al 𝛿 𝛿 𝛿 Δ

menos igual al área 1.


RESPUESTA A UNA FALLA
DE CORTOCIRCUITO

Un sistema conformado por cuatro


unidades de 13.9 MVA, 10 kV y 50 Hz se
conecta a
una barra infinita a través de un sistema
de transmisión.

• La reactancia de la red mostradas en


la figura, están en por unidad sobre
una base de 55.6 MVA,
10 kV (referido al lado de baja
tensión).
RESPUESTA A UNA FALLA
DE CORTOCIRCUITO

• Condiciones iníciales:
𝑃 = 0.9 ; 𝑄 = 0.436 ; 𝑉 = 1/28.34°

𝐸 = 0.90081/0° 𝑥′ = 0.3 𝐻 = 3.5

• La f.e.m. del generador se determina


como:

𝐸 = 𝑉 + 𝑗𝑥′ 𝐼 = 1.1626/41.77°
RESPUESTA A UNA FALLA
DE CORTOCIRCUITO

La potencia entregada por el


Generador es:
𝑉𝑉
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
𝑥
La impedancia equivalente es:
0.5 ∗ 0.93
𝑥 = 0.3 + 0.15 +
0.5 + 0.93
𝑥 = 0.7752

La potencia es
1.1626 0.90081
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
0.7752
𝑃 = 1.351 𝑠𝑒𝑛𝑜 𝛿
RESPUESTA A UNA FALLA
DE CORTOCIRCUITO

La potencia entregada por el


Generador es:
𝑉𝑉 V=0
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
𝑥
1.1626 0
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
0.45

𝑃 =0
RESPUESTA A UNA FALLA
DE CORTOCIRCUITO

La potencia entregada por el


Generador es:
𝑉𝑉
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
𝑥

1.1626 0.90081
𝑃 = 𝑠𝑒𝑛𝑜 𝛿
0.95

𝑃 = 1.1024 𝑠𝑒𝑛𝑜 𝛿
SIMULACIÓN EN MATLAB

Modelando en Simulink de Matlab

También podría gustarte