Estabilidad en El Sistema Multi Maquina Kundur

Estabilidad en un sistema
multimaquina
Presentado por:
Andrs Felipe Ros Tamayo.
Jorge Alexander Rodrguez Ros.
presentado a: Phd. Sandra Milena Prez Londoo
Universidad Tecnolgica de Pereira

Facultad de ingenieras
2016
Enunciado del problema:
Para el sistema de dos reas de Kundur:
1. Calcular los valores propios del sistema. Puede emplear cualquier paquete
como Neplan, PSAT, Digsilent, entre otros.
2. Calcular los valores propios del sistema usando Matlab. Comparar los
resultados con los obtenidos en el tem anterior.
3. Calcular los vectores propios, factores de participacin y relacin de

amortiguamiento del sistema.
4. Identificar los modos de oscilacin del sistema.
5. Analizar los resultados, graficar los valores propios y determinar cules
generadores tienen ms influencia sobre cada uno de los modos de

oscilacin.
calculo de los valores propios del sistema usando Neplan:
el sistema de dos areas del kundur lo trae Neplan por defecto:
las maquinas sncronas estn configuradas por defecto para trabajar con el
modelo subtransitorio, pero como vamos a trabajar el caso de estabilidad de
pequea seal, escogemos el modelo clsico de la maquina en la ventana de
configuraciones de cada generador:
los valores propios arrojados por Neplan son los siguientes:

calculo de los valores propios del sistema usando Matlab:
se obtiene la Ybus del sistema:
se corre un flujo de potencia para conocer el voltaje en todos los nodos y a

partir de estos calcular las demandas y condensadores vistos como
admitancias:
L7=(967-100*1i)/100;% se calcula la demanda total en el nodo 7 en P.u

Y_L7=conj(L7)/Voltajes(7,1)^2;
L9=(1767-250*1i)/100;% se calcula la demanda total en el nodo 9 en P.u

Y_L9=conj(L9)/Voltajes(9,1)^2;
las demandas y condensadores conectados a los nodos 7 y 9 son agregados a

la Ybus como admitancias.
ybus_inicial(7,7)=ybus_inicial(7,7)+Y_L7;
ybus_inicial(9,9)=ybus_inicial(9,9)+Y_L9;
se aplica la reduccion de kron:
Ybus=ybus_inicial(1:4,1:4)-ybus_inicial(1:4,5:end)...
*inv(ybus_inicial(5:end,5:end))*ybus_inicial(5:end,1:4);
el sistema resultante es el siguiente:
ya teniendo una Ybus reducida, se agregan las admitancias de los

generadores a la Ybus reducida:
Yg=1/(1i*Xdp);
ybus(1,1)=ybus(1,1)+Yg;
aux=[Yg,Yg,Yg,Yg];
El sistema resultante es el siguiente:

y se aplica de nuevo una reduccin de kron:
YbusN=[diag(aux),-diag(aux);-diag(aux),Ybus];
obtenindose una Ybus que adems de contener la informacin de las

admitancias de las lneas y transformadores, tambin contiene las
admitancias internas de los generadores y las cargas y condensadores vistos
como admitancia.
(las barras de este sistema estn al voltaje interno de cada generador)
se hallan los voltajes internos de los generadores (E1,E2,E3 y E5) a partir de

las corrientes inyectadas por cada generador y su respectivo voltaje en
terminales:
I1=conj((700+185*1i)/(100*V(1)));
E1=V(1)+I1*Xdp*1i;
E(1,1)=abs(E1);
E(1,2)=angle(E1);
I2=conj((700+185*1i)/(100*V(2)));
E2=V(2)+I2*Xdp*1i;
E(2,1)=abs(E2);
E(2,2)=angle(E2);
I3=conj((700+185*1i)/(100*V(3)));
E3=V(3)+I3*Xdp*1i;
E(3,1)=abs(E3);
E(3,2)=angle(E3);
I4=conj((700+185*1i)/(100*V(4)));
E4=V(4)+I4*Xdp*1i;
E(4,1)=abs(E4);
E(4,2)=angle(E4);
con el procedimiento anterior hemos pasado de la red que se muestra en la

figura (a) a la red que se muestra en la figura (b).
asumiendo una variacin pequea del ngulo, la ecuacin de oscilacin

puede ser escrita de la siguiente manera para cada generador:
la potencia activa en cada nodo es una funcin no lineal de los voltajes del
sistema, por lo que es necesario linealizar estas ecuaciones en la vecindad de
su punto de operacin, esta linealizacion se lleva a cabo aplicando la
expansin de taylor de primer grado, el cambio en la potencia activa
inyectada en todos los nodos puede ser expresada como se muestra
acontinuacin:
de la matriz anterior nicamente necesitamos los trminos H, los cuales se

obtienen de la siguiente manera:
la siguiente ecuacin matricial describe el comportamiento dinmico del

sistema de potencia :
hallamos los valores propios:
en Matlab este clculo se hace con el comando: [lambda]=eig(A) el cual

arroja los siguientes valores:
el siguiente diagrama de polos muestra grficamente la ubicacin de los

valores propios:
a partir de los valores propios, calculamos el damping ratio el cual determina
la tasa de decaimiento de la amplitud de la oscilacin.
Luego calculamos la frecuencia de oscilacin para cada valor propio con la

siguiente ecuacin:
Factores de participacin:
los factores de participacin son obtenidos a partir de un anlisis de
sensibilidad de cada valor propio:
diferenciando a ambos lados con respecto a a kj se obtiene el siguiente

termino:
pre multiplicando por a ambos lados y notando que y adems

que llegamos a la siguiente ecuacin:
todos los elementos de sern iguales a cero, excepto por los

elementos en la K_esima fila y la J_esima columna, por lo tanto:
Por lo tanto, la sensibilidad del valor propio con respecto al elemento

de la matriz de estados, es igual al producto del elemento del vector propio
izquierdo y el elemento del vector propio derecho .
Para calcular los factores de participacin, se hace uso de los vectores propios
por derecha (kundur) 0 [W] (Machowsky) y los vectores propios por
izquierda (kundur) 0 [U] (Machowsky)
Los factores de participacin se hallan de la siguiente manera:
Los coeficientes Pki=U ikW ki se refieren a los factores de participacin,

cada factor de participacin es el producto del K_esimo elemento del i_esimo
vector propio izquierdo y derecho, este cuantifica la sensibilidad del i_esimo
valor propio con respecto al K_esimo elemento de la diagonal de la matriz
[A].
Los resultados que arroja Neplan son los siguientes:

Los resultados calculados en Matlab son:
Conclusin:
*La Ybus juega un papel importante a la hora de hacer una anlisis de
pequea seal, ya que esta debe contener la informacin de todos los
elementos del sistema, incluyendo las demandas , los elementos
compensadores como los condensadores y las impedancias internas de los
generadores.
*el numero de nodos del sistema debe ser reducido al nmero de

generadores que lo componen.
*se debe tener cuidado con el cambio de base de los parmetros de los
generadores, ya que el cambio de base del valor H (almacenamiento de
energa cintica) suele pasarse por alto.
*para que la cantidad de valores propios que arroja Neplan sea igual a los
calculados en Matlab, se debe trabajar la maquina con su modelo clsico, y el
AVR debe ser eliminado de cada mquina.
*El anlisis de modos de oscilacin permite obtener las frecuencias de

oscilacin y su impacto en el sistema, as como los generadores a los que hay
que controlar para evitar que aparezcan oscilaciones en el sistema.

Estabilidad en El Sistema Multi Maquina Kundur

Cargado por

Copyright:

Formatos disponibles

Estabilidad en El Sistema Multi Maquina Kundur

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Estabilidad en El Sistema Multi Maquina Kundur

Cargado por

Copyright:

Formatos disponibles

Estabilidad en un sistema

Universidad Tecnolgica de Pereira

Para el sistema de dos reas de Kundur:

3. Calcular los vectores propios, factores de participacin y relacin de

4. Identificar los modos de oscilacin del sistema.

5. Analizar los resultados, graficar los valores propios y determinar cules

generadores tienen ms influencia sobre cada uno de los modos de

los valores propios arrojados por Neplan son los siguientes:

se corre un flujo de potencia para conocer el voltaje en todos los nodos y a

L7=(967-100*1i)/100;% se calcula la demanda total en el nodo 7 en P.u

L9=(1767-250*1i)/100;% se calcula la demanda total en el nodo 9 en P.u

las demandas y condensadores conectados a los nodos 7 y 9 son agregados a

se aplica la reduccion de kron:

el sistema resultante es el siguiente:

ya teniendo una Ybus reducida, se agregan las admitancias de los

El sistema resultante es el siguiente:

obtenindose una Ybus que adems de contener la informacin de las

(las barras de este sistema estn al voltaje interno de cada generador)

se hallan los voltajes internos de los generadores (E1,E2,E3 y E5) a partir de

con el procedimiento anterior hemos pasado de la red que se muestra en la

asumiendo una variacin pequea del ngulo, la ecuacin de oscilacin

de la matriz anterior nicamente necesitamos los trminos H, los cuales se

la siguiente ecuacin matricial describe el comportamiento dinmico del

en Matlab este clculo se hace con el comando: [lambda]=eig(A) el cual

el siguiente diagrama de polos muestra grficamente la ubicacin de los

Luego calculamos la frecuencia de oscilacin para cada valor propio con la

diferenciando a ambos lados con respecto a a kj se obtiene el siguiente

pre multiplicando por a ambos lados y notando que y adems

todos los elementos de sern iguales a cero, excepto por los

Por lo tanto, la sensibilidad del valor propio con respecto al elemento

Los coeficientes Pki=U ikW ki se refieren a los factores de participacin,

Los resultados que arroja Neplan son los siguientes:

*el numero de nodos del sistema debe ser reducido al nmero de

*El anlisis de modos de oscilacin permite obtener las frecuencias de

También podría gustarte