1 Flujo de carga en Sistemas de Potencia.Modelado de Sistemas de Potencia Flujo de carga en Sistemas de Potencia. Dr. Ing. Mario Vignolo

2 CONTENIDO: GeneralidadesModelado del sistema y planteo del problema del flujo de carga Solución del flujo de carga Método de Newton Raphson para la resolución del flujo de carga Método Desacoplado rápido

3 PROPÓSITO DEL FLUJO DE CARGA:Determinación de voltajes, intensidades y potencias activas y reactivas en distintos puntos de una red eléctrica.

4 HIPÓTESIS DE TRABAJO: Sistemas en régimen, equilibrados, sinusoidales, sin anomalías.

5 Importancia de los flujos de cargaPermite determinar los flujos de potencia activa y reactiva en una red eléctrica. Permite determinar los voltajes en las barras de una red eléctrica. Permite calcular las pérdidas en una red eléctrica. Permite estudiar las alternativas para la planificación de nuevos sistemas o ampliación de los ya existentes. Permite evaluar los efectos de pérdidas temporales de generación o de circuitos de transmisión.

6 Importancia de los flujos de cargaPermite evaluar los efectos de reconfigurar los circuitos de un SEP (por ejemplo ante la pérdida de una línea de transmisión). Permite evaluar las mejoras que se producen ante el cambio en la sección de los conductores de un SEP.

7 Carga, generación y modelado de la red en análisis de flujo de carga.

8 Modelado de los componentes del sistema.Líneas de transmisión - circuito Pi Transformadores - impedancia Generadores - Potencia activa constante con capacidad de control (limitado) de voltaje del primario (P = cte, V= cte). Cargas - Potencia compleja constante (P = cte, Q= cte).

9 Línea de transmisión. i k i k

10 Generadores y Cargas. GeneradoresPotencia Activa - inyección constante Potencia reactiva - regulación de voltaje Demanda de carga Inyección constante de potencia activa y reactiva

11 Flujo de carga & Balance de potencia1 i k Carga n

12 Análisis Voltaje - Corriente Análisis voltaje - potencia.versus Análisis voltaje - potencia. 1 i k Carga n

13 Análisis Voltaje - Corrientey la Matriz Ybus Carga i 1 k n Sistema de ecuaciones lineales Vtierra=0

14 Análisis Voltaje - PotenciaInyección en la red 1 i G k n Sistema de ecuaciones no lineales

15 Forma de las ecuaciones de flujo de carga.Voltaje en forma polar Voltaje en forma rectangular Admitancia en forma polar Admitancia en forma rectangular

16 Forma polar de las ecuaciones de flujo de cargaEl voltaje está expresado en coordenadas polares, mientras que la admitancia está expresada en coordenadas rectangulares.

17 Balance de potencia activa y reactiva.1 k n G i 1 k n G

18 Ecuaciones de flujo de cargai=1,2,3...n funciones de voltajes complejos desconocidos especificado balance de pot. activa y reactiva

19 Ecuaciones de flujo de cargaSi la potencia activa o reactiva para la barra i no es especificada, la ecuación de balance de energía no puede ser definida. (si la barra i no tiene generación o carga, la potencia especificada es igual a cero.) Potenciales variables desconocidas:

20 Tipos de barras Barras de carga (PQ): No hay generaciónPotencia activa y reactiva especificada Barras de generación (PV): Voltaje constante y especificado Potencia activa especificada

21 Número de incógnitas y número de ecuacionesHipótesis: Sistema de n barras Ng - cantidad de barras de generación y voltaje controlado Nd - cantidad de barras de carga n = Ng + Nd

22 Número de incógnitas y número de ecuacionesPara cada barra de generación tengo: una ecuación de balance de potencia activa el voltaje de la barra especificado Para cada barra de carga tengo: una ecuación de balance de potencia reactiva

23 Número de incógnitas y número de ecuacionesCuatro variables por cada barra: Las potencias reactivas Qi de las barras de generación pueden ser calculadas una vez determinados los voltajes de las barras (módulos y fases)

24 Barra flotante ¿Es posible especificar la potencia activa inyectada por todos los generadores y la potencia activa consumida por las cargas en forma independiente? Las pérdidas RI2 no son conocidas inicialmente

25 Barra flotante ¿Es este criterio razonable?Una barra del sistema puede realizar el balance de potencia activa demandada y potencia activa consumida (BARRA FLOTANTE) ¿Es este criterio razonable? La potencia activa se transmite “bien” a través del sistema

26 Barra flotante ¿Cómo se realiza el balance de potencia reactiva en el sistema? ¿Es posible utilizar una única barra para realizar el balance de reactiva en el sistema? La potencia reactiva no se transmite “bien” a través del sistema (produce caídas de tensión importantes) Cada barra PV realiza el balance de reactiva en forma local

27 Modelado de sistemas de potencia.Resolviendo el problema de flujo de carga.

28 Ejercicio: Ecuaciones de flujo de carga.Formar Matriz Ybus del sistema. Determinar tipos de barras. Listar variables conocidas y desconocidas. Escribir las ecuaciones de flujo de carga. 1 2 3 P=0.5 V=1 P=1, V=1 j0.1 j0.2 j0.25 1.5+j0.8

29 Ybus.

30 Tipos de barras. Barra 1: Flotante (V1 y 1 dados)Barra 2: Barra PQ (V2 y 2 desconocidos) 2 ecuaciones - balance de potencia activa y reactiva. Barra 3: Barra PV - 3 desconocido (V3 especificado) 1 ecuación: balance de potencia activa. 1 2 3 P=0.5 V=1 P=1, V=1 j0.1 j0.2 j0.25 1.5+j0.8

31 Ecuaciones.

32 Métodos para resolver las ecuaciones de flujo de carga.Sistema de ecuaciones algebraicas no lineales. Métodos: Método de Gauss-Seidel. Método de Newton-Raphson. Algoritmo de desacoplado rápido de flujo de carga.

33 Método de Newton Raphson.Idea básica. 1 4 6

34 Método de Newton - Raphson.Ejemplo ¿Qué tan buena es esta aproximación?

35 Método de Newton Raphson.Ejemplo

36 Método de Newton Raphson.Ejemplo

37 Método de Newton-Raphson.Ejemplo

38 Método de Newton-Raphson.Resumen El caso de una dimensión:

39 Sistemas de ecuaciones no lineales.Sistema general de ecuaciones algebraicas no lineales simultáneas. f1,...fn, son funciones dadas, x1,...xn, son incógnitas.

40 Método de Newton-RaphsonAproximación lineal por Taylor:

41 Método de Newton-RaphsonSupongamos que tomamos una estimación inicial de la solución x=xr

42 Método de Newton-RaphsonEstimación del error x:

43 Método de Newton-RaphsonMatriz Jacobiana Vector de apartamiento estimador lineal del error

44 Método de Newton-Raphsonestimador lineal del error

45 Método de Newton-RaphsonEstimador mejorado del valor supuesto inicialmente

46 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia Elegir las variables de estado (x): (a) Para barras PQ, elegir la magnitud del voltaje de barra y su ángulo de fase asociado. (b) Para barras PV, elegir el ángulo de fase (la magnitud del voltaje es fija) Para barra flotante (referencia), tanto magnitud de voltaje como ángulo de fase son cantidades especificadas. PQ&PV PQ

47 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia especificado funciones de x desconocidas

48 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia PQ&PV PQ PQ&PV PQ

49 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia

50 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia

51 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia

52 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia PQ&PV PQ

53 Método de Newton Raphson. Aplicación al flujo de carga del sistema de potencia Características del método: 1. Velocidad de convergencia ‘cuadrática’ (el número de cifras significativas se duplica luego de cada iteración) 2. Confiable, no sensible a la elección de la barra flotante. 3. Solución precisa obtenida luego de 4-6 iteraciones. 4. J debe ser re-calculada e invertida luego de cada iteración. (J es una matriz esparsa, tiene estructura simétrica, pero los valores no son simétricos)

54 Método de Newton RaphsonEjemplo Resolver el problema de flujo de carga usando el método de NR: V=1, =0 j0.1 1.5+j0.8 2 1 j0.25 j0.2 3 P=1, V=1

55 Método de Newton-RaphsonEjemplo V=1, =0 Barra 1: Flotante (V1 y 1 dados) Barra 2: Barra PQ (V2 y 2 desconocidos) 2 ecuaciones - balance de potencia activa y reactiva. Barra 3: Barra PV - 3 desconocido (V3 especificado) 1 ecuación: balance de potencia activa. j0.1 1.5+j0.8 2 1 j0.2 3 j0.25 P=1, V=1

56 Método de Newton-RaphsonEjemplo

57 Método de Newton-RaphsonEjemplo

58 Método de Newton-RaphsonEjemplo

59 Método de Newton-RaphsonEjemplo

60 Método de Newton-RaphsonEjemplo

61 Método de Newton-RaphsonEjemplo

62 Método de Newton-RaphsonEjemplo Esto completa la primer iteración. Ahora re-calculamos las potencias de la barra con los nuevos valores de las variables de estado:

63 Método de Newton-RaphsonEjemplo

64 Método de Newton-RaphsonEjemplo

65 Método de Newton-RaphsonEjemplo

66 Método de Newton-RaphsonEjemplo Esto completa la segunda iteración. Ahora re-calculamos las potencias de la barra con los nuevos valores de las variables de estado:

67 Método de Newton-RaphsonEjemplo

68 Método de Newton-RaphsonEjemplo

69 Método de Newton-RaphsonEjemplo

70 Método de Newton-RaphsonEjemplo Esto completa la tercera iteración. El método ha convergido ya que el vector de apartamiento es casi cero.

71 Método de Newton-RaphsonEjemplo

72 Desacoplado rápido del flujo de carga (FD)Desacoplando las ecuaciones PQ&PV PQ

73 Desacoplado rápido del flujo de carga (FD)Desacoplando las ecuaciones PQ&PV PQ Las ecuaciones están desacopladas pero los coeficientes de las matrices H y L son interdependientes: H depende del módulo del voltaje, L depende del ángulo de fase. Este esquema requiere evaluación de las matrices en cada iteración.

74 Simplificaciones de Stott & Alsac1. Las diferencias entre los ángulos de fase de barras típicas del sistema son usualmente pequeñas: 2. Las susceptancias de línea Bikson mucho mayores que las conductancias de línea Gik: 3. La potencia reactiva inyectada en cualquier barra es mucho menor que la potencia reactiva que circularía si todas las líneas que parten de esa barra se corticircuitaran al neutro del sistema:

75 Elementos Jacobianos Potencia activa

76 Elementos Jacobianos Potencia reactiva

77 Modificaciones posterioresPQ&PV PQ PQ&PV PQ

78 Modificaciones posterioresPQ&PV PQ Desacoplado rapido de las ecuaciones. PQ&PV PQ

79 Método de desacoplado rápidoCaracterísticas PQ&PV PQ 1. B’ y B’’ son matrices esparsas reales. 2. B’ y B’’ son aproximaciones del Jacobiano con gradiente constante. (El resultado final es el correcto!) 3. Aunque FD requiere más iteraciones, la solución se puede obtener mucho más rápido. 4. FD es más robusto que NR (puede encontrar soluciones donde NR falla) 5. Problemas potenciales en redes con R>X.

80 Método de desacoplado rápidoEjemplo Resolver el problema de flujo de carga usando el método FD: V=1, =0 j0.1 1.5+j0.8 2 1 j0.25 j0.2 3 P=1, V=1

81 Método de desacoplado rápidoEjemplo Barra 1: Flotante (V1 y 1 dados) Barra 2: Barra PQ (V2 y 2 desconocidos) 2 ecuaciones - balance de potencia activa y reactiva. Barra 3: Barra PV - 3 desconocido (V3 especificado) 1 ecuación: balance de potencia activa. j0.1 1.5+j0.8 2 1 j0.2 3 j0.25 P=1, V=1

82 Método de desacoplado rápidoEjemplo

83 Método de desacoplado rápidoEjemplo

84 Método de desacoplamiento rápidoEjemplo Apartamiento de potencia activa

85 Método de desacoplado rápidoEjemplo

86 Método de desacoplado rápidoEjemplo Apartamiento de potencia reactiva

87 Método de desacoplado rápidoEjemplo