Utilización de las curvas PV y QV como herramienta de análisis de

planificación.

Ing. Claudio Saldivar

Facultad Politécnica
Universidad Nacional del Este
Ciudad del Este-Paraguay
saldivar.claudio@gmail.com

Resumen
En el presente trabajo se aborda el problema de seguridad del voltaje en los sistemas de potencia. En este
sentido, se analiza los márgenes de carga activo y reactivo de un determinado sistema con la curvas PV y QV
respectivamente. Para la utilización de estas herramientas dentro de un determinado escenario de planificación,
con algunas técnicas de estabilidad de tensión conocidas, se ha procedido a seleccionar las barras candidatas para
la instalación de generadores y equipos de compensación de potencia reactiva. Para este propósito, a partir de un
caso base, se calculan el margen de carga activo y reactivo del sistema. Además, con el margen de potencia
reactiva de cada barra, la cual es obtenida con el uso de la curva de QV, es posible recabar informaciones
relacionadas al soporte de reactivos, con lo cual se visualiza futuros problemas de inestabilidad de tensión. En
particular, la aplicación de estas técnicas para el sistema paraguayo es muy interesante, debido a que este sistema
contiene a su principal centro de cargas, sistema metropolitano, conectado de forma radial desde los centros de
generación y con distancias bien significativas, con lo cual podrían surgir algunos problemas asociados a la
estabilidad de la tensión.
Abstract

In this paper, the security issue is addressed in voltage power systems. In this sense, the margins of active
and reactive load of a given system with PV and QV curves were analyzed respectively. For the use of
these tools within a given scenario planning techniques with some voltage stability, we proceeded to select
candidates buses for the installation some equipment such as generators and reactive compensation shunt.
For this purpose, from a base case, the margin of active and reactive load of the system is calculated.
Moreover, reactive power margin of each bus, which is obtained using QV curve, it is possible to obtain
information related reagent carrier, whereby future voltage instability problems is displayed. In
particular, the application of these techniques to the Paraguayan system is very interesting, because this
system contains a main load, metropolitan system, connected radially from the generating centers with
well significant distances, which could arise some problems associated with the voltage stability.

lo tanto, identificar al sistema de barra crítica (o

1. INTRODUCCIÓN barras), esto es, la barra donde se origina el problema,
es muy importante. Estas barras críticas constituyen
La planificación de la expansión de un sistema de el conjunto de barras candidatas para el refuerzo
potencia es una tarea compleja que involucra contra el colapso de tensión. Además, estas barras
diferentes temas como la estabilidad transitoria, son las más probables, para la reducción de pérdidas
estabilidad de tensión, la fiabilidad y los esquemas de cuando se considera la compensación de potencia
protección. Las herramientas computacionales reactiva, otro conjunto de barras de interés.
desarrolladas pueden incluir aspectos estáticos y Otra forma de evaluar la seguridad de la tensión del
dinámicos del sistema, que abarca una amplia gama sistema es por medio de la reducción de pérdidas.
de posibles estudios. En cuanto a los estudios de Esta idea proviene del hecho de que el colapso de
estabilidad de tensión, en general, un punto de tensión puede estar asociado con altos valores de
colapso de tensión se asocia con un Jacobiano pérdida. En este sentido, la reducción de las pérdidas
singular. En este punto de singularidad, conocida puede mantener un sistema alejado de la caída de
como bifurcación silla-nodo, la matriz jacobiana flujo tensión. Como se ha discutido en la literatura, esto no
de potencia tiene un valor propio cero. No se puede es correcto, ya que la reducción de pérdidas y colapso
obtener una solución más allá de este punto crítico, de tensión pueden ser diferentes problemas. Sin
que es el punto de capacidad de carga máxima del embargo, la reducción de la pérdida en el área crítica
sistema. También puede ocurrir bifurcación inducida, del sistema puede producir buenos resultados, y la
debido al agotamiento de las fuentes repentinas de reducción de pérdidas se debe considerar en los
potencia reactiva. El fenómeno Colapso de Tensión estudios de seguridad de tensión.
es local o, como mucho, un fenómeno regional. Por
Combinando los estudios de estabilidad de tensión y
reducción de pérdidas nos proporcionan resultados
interesantes que pueden evaluarse mediante el cálculo
del margen de carga, por lo que la robustez del
sistema es fundamental. Otra forma de analizar la
robustez del sistema es mediante el uso de la curva
QV [1]. Esta curva proporciona el margen de donde [J] representa la matriz Jacobiana del flujo de
potencia reactiva de cada barra y se obtiene mediante carga linealizado, y V las variables de estado
la transformación de una barra PQ en una PV y (ángulo de fase y magnitud de la tensión,
variando su nivel de tensión. Para cada nivel de respectivamente), y P0 y Q0 son las potencias netas
tensión considerado la potencia reactiva generada se activa y reactiva conectadas a cada barra. Es
almacena, generando una curva representativa de la importante notar que λestá relacionada con la
robustez de las barras con respecto al soporte de variación del parámetro del sistema, mientras que las
potencia reactiva. constantes Kp y Kq indican una dirección particular
En este trabajo se propone una estrategia para el de crecimiento de la carga asociada con cada barra. Si
escenario de planificación basado en las herramientas el crecimiento de las cargas es idéntico en todas las
teóricas descritas anteriormente. Por este motivo, se barras con factor de potencia constante en relación al
emplea el sistema paraguayo real, con todos sus caso base, entonces los valores de K asumen el valor
límites tomados en consideración. Los resultados se unitario. La longitud del paso predictor está dada por:
obtienen en las posibilidades reales de operación y se
discute la efectividad de las acciones propuestas.

2. TÉCNICAS DE MARGEN DE CARGA

en el que || . || representa la norma de vector tangente.
En esta sección se describe cómo funciona el método De esta expresión se deduce que, cuanto más
de continuación, así como la importancia del vector pronunciada sea la curva, más pequeño es el paso
tangente calculada durante el proceso. Tras la predictor, y viceversa. El método avanza con grandes
presentación del método de continuación, una pasos mientras el sistema se encuentre lejos del punto
discusión del área crítica tiene lugar. de bifurcación, y con pasos menores en la medida que
Se describe el método de la continuación utilizado se acerca a dicho punto. El punto de operación del
para el trazado de las Curvas PV, y la importancia del sistema es obtenido posteriormente con el paso
vector tangente calculado durante el proceso. corrector.
Posteriormente, en las siguientes secciones, se realiza
una discusión acerca de las pérdidas del sistema y la 2.3 Paso Corrector
determinación del área crítica. El Paso Corrector se obtiene mediante la inclusión de
2.1 Método de la Continuación una ecuación adicional. Tal ecuación surge del hecho
de que el paso predictor y el paso corrector son
El Método de la Continuación puede ser utilizado vectores perpendiculares entre sí. Una alternativa a
para trazar el camino andado por un sistema de este paso es la obtención del nuevo punto operativo
potencia desde un punto de equilibrio estable hasta el mediante el paso predictor como la estimación inicial
punto de bifurcación [4], [5], obteniéndose para un programa de flujo de carga. En general, con
finalmente el Margen de Carga del Sistema para la esta metodología, el sistema converge en muy pocas
dirección de crecimiento de carga considerada. Esta iteraciones.
metodología se basa en el siguiente modelo del Fig.1: Curva PV tipica.
sistema:

(1)
en el que x representa las variables de estado, siendo
λ un parámetro del sistema utilizado para llevar al
sistema desde un punto de equilibrio a otro. Este tipo
de modelo se ha empleado en numerosos estudios de
colapso de tensión, siendo λ normalmente el factor de
crecimiento carga / generación del sistema. Dos pasos 3. METODO DE LA CURVA QV
mueven al sistema a lo largo del camino de la
bifurcación: el paso predictor y el paso corrector. El método de la Curva QV [7], puede ser utilizado
como una herramienta de planificación de sistemas
2.2 Paso Predictor de potencia, la cual podría ser complementada con
otros estudios estáticos y dinámicos [8], [9]. El
El Paso Predictor es utilizado para definir una
análisis de las curvas QV debe ser realizado en
dirección de movimiento. El Vector Tangente puede
conjunto con las curvas PV, lo cual permite
ser utilizado para este propósito, y viene dado por:
identificar las barras críticas del sistema así como los

2
requerimientos de potencia reactiva en las barras para
garantizar la seguridad de tensión.
La posibilidad de reducir el costo computacional
asociado con el cálculo de las curvas QV es abordada
en esta sección. Al igual que para las curvas PV, no
toda la curva resulta de interés. Más bien lo que se
busca, es aquél punto de mínima inyección de
potencia reactiva, proporcionando como resultado el
margen de potencia reactiva para la barra analizada.
Para este propósito, un método relativamente
novedoso llamado método de la continuación QV, se
emplea tal como se describe a continuación.
3.1 Método de la Continuación QV
La idea es trazar la curva QV mediante pasos
controlados hasta llegar al punto mínimo de la Curva. Fig.2: Curva QV tipica.
El tamaño del paso está regulado por el vector
tangente según: 4. METODOLOGÍA

Debido a que el sistema paraguayo presenta algunos

puntos débiles con respecto a la estabilidad de
donde k es un valor escalar utilizado para acelerar o tension, la idea es evaluar cómo responde el sistema
desacelerar el cómputo; || . || denota la norma de un cuando se instala una unidad de pequeña generación.
vector tangente modificado VTQ, el cual se calcula Se evalúan los efectos en el margen de carga del
según: sistema (MC) y el margen de potencia reactiva
(MCR). En este sentido, las siguientes opciones de
barras se consideran para la instalación:
- Barras críticas identificados por el vector tangente,
En la ecuación (7), [J] es el jacobiano del flujo de
- Los barras con el menor y el margen de potencia
potencia; Q1 es un vector con todos sus componentes
reactiva mayor, calculado por la curva QV,
iguales a cero, excepto por aquella asociada a la
Los resultados obtenidos se comparan, por lo que se
potencia reactiva de la barra en estudio. Es
obtiene una buena estrategia de planificación.
importante mencionar que los límites de potencia
reactiva para todas las demás barras de tensión 5. RESULTADO DEL TEST
controlada (barras PV) son tenidos en cuenta. Con la
longitud del paso calculada con la ecuación (6), las 5.1 Test Sistemas IEEE
tensiones en las barras PQ se actualizan según: Las metodologías propuestas fueron puestas a prueba
mediante el uso de los sistemas de prueba IEEE.
Debido a que los resultados cualitativos obtenidos
son similares para todos los sistemas de prueba, se
discute sólo el sistema de 30 barras. La idea es
La actualización de las variables de estado a través de considerar pequeñas unidades generadoras y obras de
la ecuación (8) (paso corrector) mantiene al problema transmisión en las barras indicadas por los siguientes
bien condicionado durante el proceso de índices:
convergencia del problema de flujo de potencia,
durante el cual, la barra analizada es considerada - Vector Tangente, caso base (VT): La barra con la
como una barra del tipo PV. Sin embargo, el cálculo entrada más grande en el vector tangente calculada en
de la longitud del paso y la corrección del nivel de el caso base.
tensión se ejecutan considerando a esta barra como - Margen reactiva de potencia (MCR): Barra con el
una del tipo PQ. Esto ocurre independientemente del margen de potencia reactiva por lo menos en el caso
tipo de barra original en estudio. El costo base.
computacional para la obtención de VTQ, y - Vector tangente en el punto de colapso de tensión
consecuentemente de la longitud del paso y del (VTC): La barra más crítica en el punto de colapso
término de corrección, es sin embargo muy reducido, de voltaje indicado por el vector tangente.
ya que su cómputo sólo se ejecuta luego de la - Mezcla las técnicas (Mix): Se obtiene un conjunto
convergencia del método de Newton-Raphson, por lo de barras mediante la mezcla de las técnicas
que no se requiere de un nuevo cálculo de la matriz anteriores. Esto produce una reducción de las barras a
jacobiana. ser considerado. Se elige la barra que produce la
menor MCR.
La Figura I muestra los resultados obtenidos: A seguir en la tabla II se presenta el Listado de obras
de expansión tomadas en consideración para
Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del
sistema IEEE 30 Barras en 1˚ vecindad (Flow Tradicional) comparación de los márgenes de carga obtenidos en
1.2000
los diferentes métodos utilizados.
Aumento del Margen Carga (pu)

Tabla II: Listado de Obras de Expansión del sistema Este de

1.1000
ANDE.
Caso Obra realizada - Descripción Sistema
Pybz1 - base Sin obras de ampliación (caso base) Este
Pybz2 Margen Derecha - Acaray (nueva LT 220 kV) Este
1.0000
Pybz3 Acaray - Presidente Franco (nueva LT 220 kV Este
30b8z - GB30 30b12z - GB30 30bz9 - GB29 30bz8 - GB29 30bz9 - GB27

Casos analizados - Sistema IEEE 30 Barras

30bz10 - GB27 30bz11 - GB27 30bz12 - GB27
Pybz4 TR 500/220 kV en SE Margen Derecha (Nuevo TR) Este
Casos Anarede - Flow - Expansión Nuevas LTs Casos Anarede - Flow - Expansión Generación Local 5 MW Pybz5 Paranambú - Naranjal (nueva LT 66 kV) Este
Casos Anarede - Flow - Expansión Generación Local 10 MW Casos Anarede - Flow - Expansión Generación Local 15 MW
Pybz6 Transformador 66/23 kV en Naranjal (Nuevo TR) Este
Fig. 1: Comparación de los resultados obtenidos para el IEEE 30 – Pybz7 Transformador220/66 kV en Paranambú (Nuevo TR) Este
Método Flujo Tradicional (Anarede). Pybz8 Generación de 5 MW en la barra de 23 kV de Naranjal Este
Pybz9 Generación de 10 MW en la barra de 23 kV de Naranjal Este
Pybz10 Generación de 15 MW en la barra de 23 kV de Naranjal Este
En la figura I se observa los valores de MC obtenidos Pybz11 Generación de 20 MW en la barra de 23 kV de Naranjal Este
Pybz12 Margen Derecha - Itakyry (nueva LT 220 kV) Este
utilizando un de programa de flujo de potencia Pybz13 Transformador 220/66 kV en Itakyry (Nuevo TR) Este
(Anarede), y en la Figura 2 utilizando el método de la Pybz14 Itakyry - Catueté (nueva LT 66 kV) Este
continuación. Puede observarse también, que por el Pybz15 Transformador 66/23 kV en Catueté (Nuevo TR) Este
Pybz16 Generación de 5 MW en la barra de 23 kV de Catueté Este
método de la continuación los MC obtenidos son Pybz17 Generación de 10 MW en la barra de 23 kV de Catueté Este
superiores a los de método tradicional. Pybz18 Generación de 15 MW en la barra de 23 kV de Catueté Este
Pybz19 Generación de 20 MW en la barra de 23 kV de Catueté Este

Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del

sistema IEEE 30 Barras en 1˚ vecindad (Método de la continuación)

2.0000
Los valores de márgenes de potencia reactiva cerca
1.9000 de la zona de Naranjal son indicados a continuación:
1.8000
Aumento del Margen Carga (pu)

1.7000

1.6000
Margen de Carga Reactiva (MCR) - Expansión cerca de la barra de mayor y menor MCR (Zona de Naranjal)
1.5000

1.4000 100
1.3000 0
Itakyry Barra 1 Campo Dos Presidente Franco Kilometro 30 Itakyry Barra 2 Margen Derecha Carlos A. López Paranambú Acaray
1.2000 -100
1.1000
-200
1.0000
30b8z - GB30 30b12z - GB30 30bz9 - GB29 30bz8 - GB29 30bz9 - GB27 30bz10 - GB27 30bz11 - GB27 30bz12 - GB27 -300
MVAr

Casos analizados - Sistema IEEE 14 Barras

-400
Casos Matlab - Flow - Expansión Nuevas LTs Casos Matlab - Flow - Expansión Generación Local 5 MW
Casos Matlab - Flow - Expansión Generación Local 10 MW Casos Matlab - Flow - Expansión Generación Local 15 MW
-500

Fig. 2: Resultados obtenidos para el sistema IEEE 30 – Método de -600

la continuación. -700
-800
-900
Barras del SIN - Subsistema ESTE
Esta particularidad es algo esperado debido a que con VT Pybz1 - Base VT pybz5 Pybz11 Pybz1-Base Pybz2 Pybz3 Pybz4 Pybz5 Pybz6 Pybz7 Pybz8 Pybz9 Pybz10 Pybz11

el método de la continuación el punto de bifurcación

Fig. 3: Márgenes de carga reactiva – Zona Naranjal.
es alcanzado, mientras que con el método tradicional
la misma no es atingida debido a problemas de
convergencia que surgen cuando el aumento de la Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del sistema real
estudiado 1˚ y 2˚ vecindad (Flow Tradicional) - Zona Naranjal

carga tiende hacia el punto de la bifurcación y es 1.10

utilizado el método de Newton Raphson.

Aumento del Margen Carga (pu)

5.2 Resultados para el sistema de ANDE

En esta sección se presentan los resultados obtenidos

cuando se emplea el Subsistema 1 (SS1) de Paraguay 1.00
pybz (Base) MD-ACY ACY - PFO Otro TR MD PBU-NAR Otro TR NAR Otro TR PBU

y específicamente dentro del mismo el Sistema Este. Casos analizados - Sistema IEEE 30 Barras
Obra de Transmisión Gen 5 MW Barra 250 Gen 10 MW Barra 250 Gen 15 MW Barra 250 Gen 20 MW Barra 250

Las Subestaciones pertenecientes a este sistema así

Fig. 4: Margen de Carga (MC) para expansiones cerca de la barra
como sus respectivos niveles de tensión son crítica en el sistema real estudiado con flujo de potencia Tradicional
indicados en la tabla I a continuación: – Zona Naranjal.

Tabla I: Subestaciones que componen el subsistema Este de la

ANDE y sus respectivos niveles de tensión. Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del sistema real
estudiado 1˚ y 2˚ vecindad (Método de la continuación) - Zona Naranjal
1.10
Subestaciones Niveles de tensión (kV)
Margen Derecha 500 - 220 - 66
Aumento del Margen Carga (pu)

Acaray 220 - 66 - 23
Kilometro 30 220 - 66 - 23
Presidente Franco 220 - 66 - 23
Itakyry 220 - 66 - 23
Campo Dos 220 - 23
Alto Paraná 66 - 23
Hernandárias 66 - 23
1.00
Kilometro 15 66 - 23 pybz (Base) MD-ACY ACY - PFO Otro TR MD PBU-NAR Otro TR NAR Otro TR PBU

Curuguaty 66 - 23 Obra de Transmisión Gen 5 MW Barra 250

Casos analizados - Sistema IEEE 30 Barras
Gen 10 MW Barra 250 Gen 15 MW Barra 250 Gen 20 MW Barra 250

Catueté 66 - 23
Minga Porá 66 - 23 Fig. 5: Margen de Carga (MC) para expansiones cerca de la barra
Paranambú 220 - 66 - 23 crítica - Método de la continuación – Zona Naranjal.
Naranjal 66 - 23
Carlos Antonio López 220 - 23

4
De igual manera los valores de márgenes de potencia actuaciones de control propuestas. Comparación con
reactiva cerca de la zona de Catueté son indicados a los criterios utilizados en la actualidad muestra una
continuación: metodología que resulto ser muy eficaz para el
análisis.
Margen de Carga Reactiva (MCR) - Expansión cerca de la barra de mayor y menor MCR (Zona de Catueté)

0
Itakyry Barra 1 Campo Dos Presidente Franco Kilometro 30 Itakyry Barra 2 Margen Derecha Carlos A. López Paranambú Acaray

-200 7. REFERENCIAS
MVAr

-400

-600
[1] Fritz W. Mohn, A. C. Zambroni de Souza,
-800
Tracing PV and QV Curves with the Help of
-1000
Barras del SIN - Subsistema ESTE
a CRIC Continuation Method. IEEE
VT Pybz1 - Base
Pybz15
VT pybz19-Gen 20 MW
Pybz16-Gen 5 MW
Pybz1-Base
Pybz17-Gen 10 MW
Pybz12
Pybz18-Gen 15 MW
Pybz13
Pybz19-Gen 20 MW
Pybz14
Transactions on Power Systems, v.21,
Fig. 6: Márgenes de carga reactiva – Catueté. p.1104 - 1114, 2006.
[2] D. A. Alves, L. C. P. da Silva, C. A. Castro
Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del
sistema real estudiado en 1˚y 2˚ vecindad (Flow Tradicional) - Zona Catueté
e V. F. da Costa., Continuation fast
1.1000
decoupled power flow with secant predictor,
IEEE Transactions on Power Systems, Vol.
Aumento del Margen Carga (pu)

18, No. 3, Aug. 2003, pp. 1078 - 1085.

[3] P. Kundur, Power System Stability and

Control. Palo Alto: McGraw-Hill, 1994.
1.0000
pybz (Base) Otro TR MD MD - IRY Otro TR IRY LT IRY-CAT Otro TR CAT

Obras Transmisión
Casos analizados - Sistema IEEE 30 Barras
Gen 5 MW Barra 291 Gen 10 MW Barra 291 Gen 15 MW Barra 291 Gen 20 MW Barra 291
[4] T. V. Cutsem, C. Vournas, Voltage Stability
Fig. 7: Margen de Carga (MC) para expansiones cerca de la barra
of Electric Power Systems, Kluwer
crítica en el sistema real estudiado con flujo de potencia Academic Publishers, 1998.
Tradicional – Zona Catueté.
[5] B. H. Chowdhury, C. W. Taylor, “Voltage
Comparación de Margen de Carga (pu) - Expansión cerca de la barra crítica del
stability analysis: V-Q power flow
sistema real estudiado en 1˚y 2˚ vecindad (Método de la continuación) - Zona Catueté simulation versus dynamic simulation”,
1.1000
IEEE Transactions on Power Systems, vol.
15, No. 4, pp. 1354-1359, Nov. 2000.
Aumento del Margen Carga (pu)

[6] Show-Kang Chang, Vladimir Brandwajn.,

“Adjusted solutions in fast decoupled load
1.0000
flow”, IEEE Transactions on Power
pybz (Base) Otro TR MD MD - IRY Otro TR IRY LT IRY-CAT Otro TR CAT

Casos analizados - Sistema IEEE 30 Barras

Systems, Vol.3, N. 2, pp. 726-733, May
Obras Transmisión Gen 5 MW Barra 291 Gen 10 MW Barra 291 Gen 15 MW Barra 291 Gen 20 MW Barra 291

1988.
Fig. 8: Margen de Carga (MC) para expansiones cerca de la barra
crítica - Método de la continuación – Zona Catueté.
[7] A . C. Zambroni de Souza, Fritz W. Mohn,
Isabella F. Borges, Tito R. Ocariz, “Using
6. CONCLUSION PV and QV Curves with the Meaning of
Static Contingency Screening and
En este trabajo se propone una herramienta de Planning”, accepted for publication by the
planificación en base a unos criterios de estabilidad Electric Power Systems Research.
de tension. A tal efecto, una combinación del margen
de potencia reactiva y aumento del margen de carga [8] V. Venkatasubramanian, “Singularity
se centran. Pequeñas unidades de generación se Induced Bifurcation and the van der Pol
consideran como acciones de control en algunos Oscillator”, IEEE Transactions on Circuits
grupos de barras, como así también la inserción de and Systems I: Fundamental Theory and
obras de transmisión. En primer lugar, los sistemas de Applications, Vol. 41, No. 11, pp. 765-769,
IEEE académicas han sido probados y los resultados November 1994.
alentadores estimulado para emplear la técnica
propuesta en el sistema de Paraguay. Este sistema ha [9] M. K. Pal, “Voltage Instability Considering
sido elegido debido a su fragilidad en seguridad de Load Characteristics”, IEEE Transactions on
tensión, por lo que los resultados podrían ser mejor Power Systems, Vol. 7, No. 1, pp. 243-249,
evaluados. Una vez más, se han obtenido buenos February 1992.
resultados y algunas características no observadas en
los sistemas académicos podrían ser exploradas en el [10] R. A. Schlueter, “A Voltage Stability
sistema eléctrico de potencia de Paraguay ANDE. Security Assessment Method”, IEEE
La metodología que aquí se propone puede ser Transactions on Power Systems, Vol. 13,
utilizada para otros sistemas sin restricción, ya que No. 4, pp. 1423-1438, November 1998.
los diferentes grupos de barras se identifican y se
ensayaron como candidatos para acoger las