01 - Ejercicios de Estabilidad PDF
01 - Ejercicios de Estabilidad PDF
DIgSILENT GmbH
Heinrich-Hertz-Strasse 9
D-72810 Gomaringen
Tel.: +49 7072 9168 - 0
Fax: +49 7072 9168- 88
http://www.digsilent.de
e-mail: [email protected]
2 October 2010
Indice de Contenidos
1 Introducción ..................................................................................................................................... 4
3.1 AVR................................................................................................................................................... 9
3.1.1 Definición del elemento DSL ............................................................................................................ 9
3.1.2 Definición del modelo de planta (modelo compuesto) ........................................................................ 9
3.1.3 Respuesta del regulador al escalón ................................................................................................ 10
3.1.4 Comportamiento durante falla ....................................................................................................... 13
1 Introducción
El propósito de los siguientes ejercicios es introducir los métodos fundamentales para el estudio de los diferentes
fenómenos asociados a la estabilidad de sistemas eléctricos de potencia (SEPs). Los diferentes tipos de cálculos y
las herramientas disponibles en PowerFactory se ejercitan de manera de permitir una familiarización con los
métodos de análisis para estabilidad de SEPs, visualización e interpretación de los resultados.
Durante los ejercicios el instructor le brindará ayuda para la realización de las tareas y estará a su disposición
para contestar preguntas y atender dudas. Por favor, no dude en consultarlo.
• Este proyecto contiene el sistema de una máquina que se muestra en el diagrama unifilar de la Figura 1.
La potencia del generador sincrónico se transmite a través de un sistema de 2 líneas (de diferente
longitud) hacia una red de potencia infinita, representada en este caso por una fuente de tensión con
impedancia de salida =0. En los próximos ejercicios se analizarán diferentes aspectos referidos a la
estabilidad de la operación de esta máquina.
• Antes de comenzar, calcular un flujo de cargas y verificar mensajes de error y/o advertencias en la
ventana de salida. Corregir cualquier error antes de proseguir.
• En la página de opciones básicas seleccionar el método de simulación RMS, utilizar una representación
de la red balanceada, y seleccionar la opción para verificar las condiciones inciales.
• En la página “Tamaño de Paso” ajustar el paso de integración a 1 ms y el tiempo inicial a – 200 ms.
• Definir a continuación las variables a monitorear: picar con el botón derecho sobre el elemento, Definir -
-> Conjunto de variables (Sim).
− Generador: potencia activa y reactiva, corriente y ángulo del rotor (con referencia a la tensión de la
barra de referencia del sistema firot).
− Generador: variable “outofstep” (ver Simulación RMS\Señales). Esta variable lógica cambia de 0 a 1
cuando el generador pierde sincronismo con la red.
− Cortocircuito trifásico a la mitad (50%) de la línea ‘CCT2’ para t=0. Picar con el botón derecho la
línea, y en el menú contextual Definir -> Evento de Cortocircuito. Ajustar allí el tiempo de ejecución
y el tipo de falla (trifásica). La impedancia de falla en este caso será 0.
− Observar que la ubicación de la falla en la línea (ej. al 50%) se define directamente en la línea
sobre la página Simulación RMS. Al elegir la ubicación activar la opción Disponible.
− Apertura tripolar de los interruptores de la línea para t=100 ms a los efectos de despejar la falla.
Picar con el botón derecho la línea, y en el menú contextual Definir -> Evento de Switcheo.
• Para generar diagramas (instrumentos virtuales) adicionales picar y seleccionar alguno de los
siguientes tipos disponibles:
• Repetir la simulación modificando el tiempo de despeje de la falla hasta encontrar de manera iterativa el
tiempo crítico de despeje de la falla.
• Picar con el botón derecho en la línea y en el menú contextual elegir Ejecutar escritos DPL.
− minTclear: 0.05 s; tiempo mínimo (el escrito no iterará si el sistema se vuelve inestable para un
tiempo menor a minTclear)
− maxTclear: 1 s; tiempo máximo (el escrito no iterará si el sistema se mantiene estable para una
duración de la falla de maxTclear)
− tstep: paso de la iteración, se recomienda usar aquí el período de la frecuencia del sistema.
• El escrito reporta el resultado a la ventana de salida. Verificar este valor con el calculado anteriormente
de forma manual.
• Ejecutar el escrito nuevamente para diferentes ubicaciones de la falla (10% y 90%). Comparar los
resultados y sacar conclusiones.
− Remover la falla 100 ms después de la aperture de los interruptores de línea. Picar con el botón
derecho la línea, y en el menú contextual Definir -> Evento de Cortocircuito. Ajustar allí el tiempo
de ejecución y en tipo de falla seleccionar Despejar Cortocircuito.
− Cierre tripolar de los interruptores de la línea para t=1 s. Picar con el botón derecho la línea, y en el
menú contextual Definir -> Evento de Switcheo.
100. km
186. km
Type CCT
Type CCT
DIgSILENT
PowerFactory 14.0.512
1299.40 698.60
453.37 243.75
1.50 0.81
HV
500. kV 528.75
1.06
16.07
-1998.00
-697.12
2.31
Trf
Trf 500kV/24kV/2220MVA
1998.00
967.68
48.14
Seminario PowerFactory
Date:
G1
Annex:
Gen 2220MVA/24kV
Graphic: Grid
8/7/2009
Project: PF Seminar
DIgSILENT
-8-Estabilidad Angular Transitoria
-9-Sistema de Excitación
3 Sistema de Excitación
A continuación se incluirá un regulador de tensión (AVR) y un estabilizador (PSS) para el generador. De esta
manera se incluirá en el estudio los reguladores del generador y se analizará su influencia sobre la estabilidad del
sistema que está siendo analizado.
3.1 AVR
• El diagrama de bloques del regulador se muestra en la Figura 2. Para mostrarlo en PF picar con el botón
derecho sobre el elemento DSL del regulador y seleccionar “Mostrar Gráfica”. Identificar sobre el
diagrama de bloques las variables de entrada y de salida del regulador.
• Mediante el administrador de datos crear un nuevo modelo compuesto (Composite Model) dentro de la
red y nombrarlo “Planta G1”.
• Ejecutar un flujo de carga y calcular las condiciones iniciales para verificar la no existencia de errores.
• Ejecutar un flujo de carga y calcular las condiciones iniciales de la simulación dinámica para verificar la
no existencia de errores. Observar que pegado al símbolo del regulador deberá aparecer la tilde
indicando que el modelo se encuentra activo.
• Definir variables de resultado para el regulador de tensión del mismo modo que se realizó en el ejercicio
anterior para los demás elementos de la red. Seleccionar todas las señales de entrada y salida del
regulador.
• Recalcular las condiciones iniciales, editar nuevamente las variables de resultado , luego picar allí con
botón derecho sobre el símbolo del regulador y seleccionar Salida -> Variables Seleccionadas. PF
imprimirá en la ventana de salida los valores iniciales de las señales de entrada y salida del regulador.
• Definir un evento de parámetro sobre la variable de entrada usetp para simular un escalón en la tensión
de referencia. Para ello:
− Editar la lista de los eventos de simulación . De ser necesario, poner eventualmente fuera de
servicio los eventos de la simulación del caso anterior.
− Nuevo Valor: salto del 5% sobre el valor inicial de la señal de entrada ustep. (ej. Nuevo valor = 1.05
para el caso de ser 1 el valor inicial de esta señal).
• Opcionalmente se pueden cambiar los ajustes de los parámetros del regulador de tensión (AVR) y
observar su influencia en la respuesta al escalón.
usetp
0
Vs
Vrmax
O
Vrmin
upss
2
Vf
K o18 KeSe
Se(efd)+Ke
Kd
Ke,E1,Se1,E2,..
0
_Fex_
curex Kc
3 1
sK/(1+sT) uerrs..
• Analizar los resultados obtenidos y compararlos con aquellos previos a la conexión del regulador de
tensión.
• Determinar el tiempo crítico de despeje de falla. ¿Se ha modificado respecto del caso §2.2.1?
• Seleccionar el diagrama de bloques “pss_STAB2A” de entre los modelos IEEE standard disponibles en la
librería global de PowerFactory.
• Ejecutar un flujo de carga y calcular las condiciones iniciales de la simulación dinámica para verificar la
no existencia de errores. Observar que pegado al símbolo del regulador deberá aparecer la tilde
indicando que el modelo se encuentra activo.
• Recalcular las condiciones iniciales, editar nuevamente las variables de resultado , luego allí picar con
botón derecho sobre el símbolo del regulador y seleccionar Salida -> Variables Seleccionadas. PF
imprimirá a la ventana de salida el valor inicial de las señales de entrada y de salida. Analizar los valores
iniciales de las variables del estabilizador.
• Comparar la respuesta de la planta (generador) tras el agregado del estabilizador de potencia con
aquella obtenida en §3.1.4 y sacar conclusiones sobre la influencia del estabilizador en el sistema.
• Determinar el tiempo crítico de despeje de falla para l línea CCT2 (misma ubicación que en §2.2.2) y
comparar resultados respecto del caso anterior.
0
-15-Sistema de Excitación
paux
pgt KsTd/(1+sT1) KsTd/(1+sT1) KsTd/(1+sT1) KB
1
K2,T2,T2 K2,T2,T2 K2,T2,T2 inv
KB
K4
A los efectos de brindar una mejor comprensión del análisis y ejercitar el uso de la herramienta en PF, se
verificarán los resultados del análisis de autovalores con aquellos correspondientes a una simulación transitoria de
una pequeña perturbación en el sistema.
Antes de comenzar, crear un nuevo caso de estudio de manera de no perder los resultados de los ejercicios
anteriores.
• Editar los eventos de simulación y asegurarse de eliminar o poner fuera de servicio los eventos de
simulación de los ejercicios anteriores.
− Picar con el botón derecho sobre el generador sincrónico en el diagrama unifilar y en el menú
contextual seleccionar Definir -> Evento de Máquina Sincrónica.
− Definir un torque adicional de 0.01 p.u. (pequeña señal o perturbación). Este torque se adicionará al
torque inicial de la máquina, produciendo una leve aceleración del rotor.
• Mediante un nuevo evento para la máquina sincrónica volver el torque adicional nuevamente a cero para
t=0.2 s.
− Ubicar el cursor sobre un punto de la curva, botón derecho para acceder al menú contextual y
seleccionar Mostrar dx/dy para medir la distancia entre dos puntos.
− De ser necesario también se podrá definir normalizar la curva (factor de escala para el eje Y). Para
ello editar el gráfico, y en la tabla de las curvas, activar en la penúltima columna “Normalizar” y
ajustar el valor normalizado según se desee.
• Poner a continuación el estabilizador (PSS) fuera de servicio, ejecutar nuevamente la simulación, y medir
la nueva frecuencia de la oscilación y el factor de amortiguamiento.
El resultado de las simulaciones transitorias se verificará a continuación mediante la función a análisis modal
(cálculo de autovalores).
• En la barra principal del programa, cambiar la barra de herramientas por (análisis de autovalores).
− Opciones avanzadas: asegurase que esté seleccionada la opción de cálculo de los factores de
participación. La variable de estado deberá ser la velocidad (s:speed).
− Valores Propios -> Participaciones -> Detallado para reportar la información detallada incluidos los
factores de participación.
− Selección de valores propios: elegir la opción de Filtro para todos los valores propios con un período
menor que 30 segundos y un amortiguamiento máximo de 10 1/s.
• Identificar el período y el factor de amortiguamiento de la oscilación principal y comparar con los valores
obtenidos de la simulación transitoria en §4.1.
• Para representar los resultados del análisis de autovalores gráficamente en el plano complejo,
seleccionar y luego . En la ventana de diálogo seleccionar Mostrar Todo y escalar los ejes X e Y.
− Picar doble sobre cualquiera de los autovalores del diagrama para acceder a sus datos (partes real e
imaginaria, factor de participación etc.).
− Picar con el botón derecho en cualquiera de los autovalores y en el menú contextual seleccionar
Generar Diagrama Fasorial o General Diagrama de Barras para acceder a una representación
fasorial o de barras de los factores de participación de ese modo de oscilación.
• Repetir el cálculo de autovalores para el caso con estabilizador (PSS en servicio) y comparar los
resultados nuevamente.
5 Estabilidad de tensión
A continuación se analizará la estabilidad de tensión del sistema de una máquina. Esto puede ser llevado a cabo
mediante el uso de curvas tensión-potencia reactiva (V-Q) y curvas potencia activa-tensión (P-V). Para agilizar el
cambio de los parámetros del sistema (por ejemplo el P de las cargas) y generar las curvas P-V y V-Q se dispone
de escritos DPL que automatizan esta tarea.
• Sobre el diagrama unifilar, marcar el generador y la barra LV (presionar la tecla control para la selección
múltiple).
• Botón derecho sobre la selección, y en el menú desplegable seleccionar Ejecutar Escritos DPL.
• En la ventana de diálogo seleccionar el escrito “V-Q Curve” que se encuentra dentro de la carpeta
“Scripts” del proyecto.
• En la ventana de diálogo del escrito ajustar los parámetros de entrada como sigue:
El escrito DPL generará automáticamente tantas curvas V-Q como pasos para al potencia active del generador se
hayan elegido. Observar además que el escrito DPL cambia el despacho del generador de P-Q a P-V (tensión en
la barra “LV”) para la generación de las curvas V-Q.
Analizar los resultados y determinar los puntos estables/inestables del sistema. Notar que también es posible
mostrar todas las curvas sobre un mismo diagrama editando la tabla de curvas de cualquiera de las gráficas.
• Opcional: crear una variación de red a los efectos de guardar las modificaciones que se introducirán a
continuación.
• Poner fuera de servicio el generador, y conectar una carga en el Terminal “LV”. Definir esta carga como
carga PQ, vale decir con una potencia activa y un factor de potencia constante. Asignarle un valor de
100 MW (el valor inicial no es relevante ya que será modificado durante el cálculo) y un factor de
potencia de 0.8 inductivo.
• Sobre el diagrama unifilar, marcar la carga, y las barras LV, HV y la barra de potencia infinita (presionar
la tecla control para la selección múltiple).
• Botón derecho sobre la selección, y en el menú contextual seleccionar Ejecutar Escritos DPL.
• En la ventana de diálogo seleccionar el escrito “P-V Curve” que se encuentra dentro de la carpeta
“Scripts” del proyecto.
• Ejecutar para iniciar la simulación. El escrito generará automáticamente las curvas P-V para las 3 barras
seleccionadas.
• Para una tensión de operación de 1 p.u. +/- 5% determinar cual es la máxima potencia que puede ser
sumistrada a una carga en la barra LV.
• Ejecutar nuevamente el escrito DPL “P-V curves2 y comparar los resultados con el caso anterior. ¿Qué
influencia tiene el conmutador bajo carga sobre la estabilidad de tensión en este caso? Determinar los
nuevos límites para la máxima carga que se puede suministrar en la barra LV.
6 Estabilidad de Frecuencia
− Importa: el área 2 está importando aprox. 500 MW del resto del sistema
• Calcular un flujo de cargas y verificar la condición inicial del sistema (perfiles de tensiones, sobrecarga
de elementos, etc.). Imprimir para ello un reporte de verificación (ventana Salidas en el comando de
flujo de cargas).
• Verificar las potencias de intercambio entre las diferentes áreas del sistema. Utilizar para ello un reporte
de resumen por área. Alternativamente podrá definir un Boundary para acceder a los resultados en
forma tabular.
• Verificar que no existan modos inestables de operación (polos en el semiplano derecho, marcados en
rojo en el diagrama).
• Cambiar las propiedades del eje Y para mostrar la frecuencia. ¿En qué rangos de frecuencia se
encuentran distribuidos los modos de oscilación?
• Chequear el modo de oscilación alrededor de 0.8 Hz con el mínimo coeficiente de atenuación. Generar
un diagrama con los factores de participación. Identificar los generadores que participan de esta
oscilación.
• Identificar en cada área los generadores con el mayor despacho por escenarios. Estos generadores son
candidatos a producir la mayor desviación de frecuencia en el sistema de producirse su desconexión.
• Simular a continuación la desconexión (n-1) de dichos generadores. Considerar los siguientes casos: 1
generador por area desconectado, dos escenarios de operación por generador.
• Ajuste de la simulación:
− Monitorear la velocidad (o frecuencia) de todos los generadores en cada una de las áreas.
• De acuerdo al código de operación del sistema, la frecuencia no puede bajar de 49.7 Hz para una
contingencia (n-1) en el sistema. ¿Se excede este límite para la desconexión de alguno de los
generadores en el sistema?
Al igual que en el caso anterior, el análisis se llevaráa cabo para ambos escenarios de operación.
• Analizar el comportamiento de la frecuencia en ambas islas tras la separación. ¿Se mantienen ambas
islas estables? Determinar la desviación de frecuencia de estado estacionario en ambas islas del sistema.
− Como consecuencia de la mala operación de una protección, ambas líneas son desconectadas.
• ¿Resulta afectada por la falla la desviación transitoria de la frecuencia en este caso con relación a la
simulación anterior?
• ¿Cómo resulta la desviación de frecuencia de estado estacionario en ambas islas del sistema en este
caso con relación al caso anterior?
• Identificar aquellos casos previamente analizados (pérdida de una unidad generadora o del interconector
con el área 2) donde se produce una desviación transitoria de la frecuencia de hasta 49 Hz o mayor.
• Editar el relé de alguna de las cargas y observar su definición y los diferentes pasos de rechazo de
carga.
• Verificar los mensajes en la ventana de salida. ¿A qué frecuencia se inicia el rechazo de carga? Verificar
los ajustes de los diferentes pasos de los relés en función de los resultados de la simulación.
• Simular nuevamente el disparo del interconector del área 2 para t=0 sin existencia de falla previa,
mediante un evento de switcheo para ambas líneas.
• Determinar la máxima y mínima desviación transitoria de frecuencia en ambas islas del sistema.
• Cambiar el droop de los reguladores de velocidad de los generadores del área 2 del valor 0.05 a 0.025
p.u. Notar que:
− Como el droop del regulador es un parámetro del modelo de la red, se recomienda aquí la definición
de una variación de red a los efectos de no modificar los valores anteriores.
• Editar el tipo de carga general y en la hoja de la simulación RMS cambiar la dependencia de P con la
frecuencia de 0 a 0.1 y luego tambien a 1.
• Es correcta la hipótesis a veces asumida de despreciar este efecto en caso de que no se conozco esta
característica de la carga?