FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO

ESTACIONARIO USANDO LA METODOLOGÍA DE REFORMULACIÓN DE

POTENCIAS ESPECIFICADAS
 UNIVERSIDAD POLITÉCNICA SALESIANA
SEDE QUITO

CARRERA:
INGENIERÍA ELÉCTRICA

Trabajo de titulación previo a la obtención del título de

INGENIERO ELÉCTRICO

TEMA:
FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO
ESTACIONARIO USANDO LA METODOLOGÍA DE REFORMULACIÓN DE
POTENCIAS ESPECIFICADAS

AUTOR:
DANIELA FERNANDA AGUAS BELTRÁN

TUTOR:
CARLOS ANDRÉS BARRERA SINGAÑA

Quito, septiembre 2020

 . Datos de Catalogación Bibliográfica .
Daniela Fernanda Aguas Beltrán

FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO

ESTACIONARIO USANDO LA METODOLOGÍA DE REFORMULACIÓN DE
POTENCIAS ESPECIFICADAS

Universidad Politécnica Salesiana, Quito – Ecuador 2020

Ingeniería Eléctrica

Breve reseña histórica e información de contacto.

Daniela Fernanda Aguas Beltrán (Y’1995). Realizó sus estudios de

nivel secundario en el Colegio Particular “Francisca de las Llagas” de
la ciudad de Quito. Egresada de Ingeniería Eléctrica de la Universidad
Politécnica Salesiana. Su trabajo se basa en el estudio de flujo de
potencias con modelos de carga en estado estacionario usando la
metodología de reformulación de potencias especificadas.
daguasb@est.ups.edu.ec

Dirigido por:

Carlos Andrés Barrera Singaña (Y’1989). Realizó sus estudios

superiores en la Escuela Politécnica Nacional de Quito, donde se graduó
de Ingeniero Eléctrico en 2013. Además, cursó estudios de posgrado en
la Universidad de Zaragoza, España, donde obtuvo el título en Máster
en Energías Renovables y Eficiencia Energética con Especialidad en
Sistemas Eléctricos. Actualmente es profesor ocasional a tiempo
completo de la Universidad Politécnica Salesiana. cbarrera@ups.edu.ec

Todos los derechos reservados:

Queda prohibida, salvo excepción prevista en la ley, cualquier forma de reproducción,

distribución, comunicación pública y transformación de esta obra para fines comerciales,
sin contar con la autorización de los titulares de propiedad intelectual. La infracción de
los derechos mencionados puede ser constitutiva de delito contra la propiedad intelectual.
Se permite la libre difusión de este texto con fines académicos o investigativos por
cualquier medio, con la debida notificación a los autores.

DERECHOS RESERVADOS
©2020 Universidad Politécnica Salesiana
QUITO - ECUADOR

I
 DECLARATORIA DE COAUTORÍA DEL DOCENTE TUTOR

Yo, Carlos Andrés Barrera Singaña declaro que bajo mi dirección y asesoría fue
desarrollado el trabajo de titulación “FLUJOS DE POTENCIA CON MODELOS DE
CARGA EN ESTADO ESTACIONARIO USANDO LA METODOLOGÍA DE
REFORMULACIÓN DE POTENCIAS ESPECIFICADAS” realizado por Daniela
Fernanda Aguas Beltrán, obteniendo un producto que cumple con todos los requisitos
estipulados por la Universidad Politécnica Salesiana para ser considerados como trabajo
final de titulación.

Quito D.M., septiembre de 2020

………………………………….
Carlos Andrés Barrera Singaña
C.C.: 0503503336

II
 CESIÓN DE DERECHOS DE AUTOR

Yo, Daniela Fernanda Aguas Beltrán, con documento de identificación N° 1724346836,

manifiesto mi voluntad y cedo a la Universidad Politécnica Salesiana la titularidad sobre
los derechos patrimoniales en virtud de que soy autora del trabajo de grado/titulación
intitulado: “FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO
ESTACIONARIO USANDO LA METODOLOGÍA DE REFORMULACIÓN DE
POTENCIAS ESPECIFICADAS”, mismo que ha sido desarrollado para optar por el título
de: Ingeniero Eléctrico, en la Universidad Politécnica Salesiana, quedando la Universidad
facultada para ejercer plenamente los derechos cedidos anteriormente.
En aplicación a lo determinado en la Ley de Propiedad Intelectual, en mi condición de
autor me reservo los derechos morales de la obra antes citada. En concordancia, suscribo
este documento en el momento que hago entrega del trabajo final en formato digital a la
Biblioteca de la Universidad Politécnica Salesiana.

Quito, D.M., septiembre de 2020

………………………………….
Daniela Fernanda Aguas Beltrán
C.C.: 1724346836

III
1. ÍNDICE GENERAL

1 Introducción ........................................................................................................... 2
2 Marco teórico ......................................................................................................... 3
2.1 Modelos de carga ................................................................................................... 3
2.2 Componente de carga ............................................................................................. 3
2.3 Clase de carga......................................................................................................... 3
2.4 Composición de la carga ........................................................................................ 3
2.5 Mezcla de clase de carga. ....................................................................................... 3
2.6 Característica de carga ........................................................................................... 3
2.7 Modelos en estado estacionario.............................................................................. 5
2.8 Modelo de impedancia constante ........................................................................... 5
2.9 Modelo de corriente constante ............................................................................... 5
2.10 Modelo de potencia constante ................................................................................ 5
2.11 Modelo ZIP ............................................................................................................ 5
2.12 Modelo lineal.......................................................................................................... 6
2.13 Modelo exponencial ............................................................................................... 6
2.14 Modelo polinomial ................................................................................................. 7
2.15 Motor de inducción (MI) ........................................................................................ 7
3 Flujo de potencia (FP) ............................................................................................ 7
3.1 Generadores sincrónicos ........................................................................................ 7
3.2 Líneas de transmisión ............................................................................................. 8
3.3 Líneas cortas ........................................................................................................... 8
3.4 Líneas medias ......................................................................................................... 8
3.5 Líneas largas ........................................................................................................... 9
3.6 Transformadores..................................................................................................... 9
3.7 Compensación de potencia reactiva en paralelo .................................................... 9
3.8 Matriz Ybarra ....................................................................................................... 10
3.9 Solución de flujos de potencia ............................................................................. 10
4 Metodología de reformulación potencias especificadas....................................... 11
5 Implementación en Matlab ................................................................................... 12
5.1 Sistema de prueba 1: 4 barras Grainger y Stevenson ........................................... 13
5.2 Sistema de prueba 2: IEEE 9 Nodos .................................................................... 13
5.3 Sistema de prueba 3: IEEE 14 nodos ................................................................... 14
6 Análisis de resultados ........................................................................................... 14
6.1 Sistema 1 .............................................................................................................. 14
6.2 Sistema 2 .............................................................................................................. 15
6.3 Sistema 3 .............................................................................................................. 15
7 Conclusiones ........................................................................................................ 17
7.1 Trabajos futuros.................................................................................................... 17
8 Referencias ........................................................................................................... 17
8.1 Matriz de Estado del Arte..................................................................................... 21
8.2 Resumen de Indicadores....................................................................................... 26
9 Anexos.................................................................................................................. 28

IV
2. ÍNDICE DE FIGURAS

Figura 1. Clases de carga y componentes [6]. .................................................................... 4

Figura 2. Aportación porcentual de los componentes de carga [12]. ................................. 4
Figura 3. Categorización de modelados de carga [9]. ........................................................ 4
Figura 4. Diagrama parejo al motor de inducción, modelo de carga estacionario [7]. ...... 7
Figura 5. Dirección de un FP de 4 barras en Power World. .............................................. 7
Figura 6. FP del modelamiento de un generador síncrono P [MW]. ................................. 8
Figura 7. FP del modelamiento de un generador síncrono Q [Mvar]. ............................... 8
Figura 8. Modelo de línea corta [19]. ................................................................................. 8
Figura 9. Modelo de línea media π [19]. ............................................................................ 8
Figura 10. Modelamiento de implementación en líneas muy largas y largas [17]. ........... 9
Figura 11. Transformador[7]. ............................................................................................. 9
Figura 12. Taps del transformador que modifican su voltaje de salida [33]. ..................... 9
Figura 13. Taps del transformador con modelo π [11]. ..................................................... 9
Figura 14. Esquema del flujo, restructuración de potencias especificadas. ..................... 12
Figura 15. Sistema de 4 barras Grainger [33]. ................................................................. 13
Figura 16. Modelo IEEE con 9 nodos [49]. ..................................................................... 13
Figura 17. Modelado IEEE 14 nodos [50] ....................................................................... 14
Figura 18. Porcentaje de error respecto a voltajes en el sistema de 4 barras. .................. 15
Figura 19. Porcentaje de error respecto ángulos en el sistema de 4 barras. ..................... 15
Figura 20. Porcentaje de error respecto a voltajes en el sistema de 9 barras. .................. 16
Figura 21. Porcentaje de error respecto ángulos en el sistema de 9 barras. ..................... 16
Figura 22. Porcentaje de error respecto a voltajes en el sistema de 14 barras. ................ 16
Figura 23. Porcentaje de error respecto ángulos en el sistema de 14 barras. ................... 17
Figura 24. Resumen e indicador de temática - matriz estado del arte.............................. 26
Figura 25. Indicador de formulación del problema-matriz estado del arte ...................... 26
Figura 26. Indicador de solución -matriz estado del arte ................................................. 27

V
3. ÍNDICE DE TABLAS

Tabla 1. Exponentes de carga [6]. ...................................................................................... 6

Tabla 2. Demanda de cargas modelo de 4 nodos. ............................................................ 13
Tabla 3. Despacho de generadores del sistema de 4 barras. ............................................ 13
Tabla 4. Valores p.u en líneas del modelo de prueba 1.................................................... 13
Tabla 5. Demanda de energía del modelo de 9 nodos...................................................... 13
Tabla 6. Despacho de generadores del sistema de 9 barras. ............................................ 13
Tabla 7. Datos de las líneas del modelado de 9 nodos. .................................................... 13
Tabla 8. Demanda de energía del modelo 14 nodos. ....................................................... 14
Tabla 9. Despacho, generadores modelo de 14 nodos. .................................................... 14
Tabla 10. Valores en p.u en líneas del modelamiento con 14 nodos. .............................. 14
Tabla 11. Porcentaje de error en voltajes y ángulos del sistema 4 barras entre
PF_CARGAS y PowerFactory. ........................................................................................ 15
Tabla 12. Porcentaje de error en voltajes y ángulos del sistema 9 barras entre
PF_CARGAS y PowerFactory. ........................................................................................ 16
Tabla 13. Porcentaje de error en voltajes y ángulos del sistema 14 barras entre
PF_CARGAS y PowerFactory. ........................................................................................ 17
Tabla 14. Datos de implementación de modelos de carga en sistema 4 barras Grainger y
Stevenson. ........................................................................................................................ 28
Tabla 15. Datos de implementación de modelos de carga en sistema IEEE 9 barras. ..... 28
Tabla 16. Datos de implementación de modelos de carga en sistema IEEE 14 barras. ... 29

VI
Resumen Abstract
El estudio de flujos de potencia es The power flow study is undoubtedly a
indudablemente un punto muy very important factor considered in
importante en los SEP, el cual permite SEP which allow to determine
determinar voltajes, pérdidas, niveles de voltages, losses, load levels, and other
carga entre otras variables para variables for the deployment of both
aplicaciones de operación y operation and planning. Hence, the
planificación. Por este motivo la accuracy in this type of electrical
precisión de este tipo de estudio eléctrico studies ought to be the highest possible
debe ser la mayor posible en cuanto a la regarding the elements that are
representación de los elementos que represented in the electrical grid. An
conforman la red eléctrica. Una mejora improvement to electrical system
al modelo del sistema eléctrico es la modeling is the consideration of steady
consideración de modelos de carga en state load models whose power flow
estado estacionario, y cuyo resultado de results are closer to reality. Having said
flujos de potencia tiene resultados más that, the following study will focus on
cercanos a la realidad. Por lo antes incorporating steady state load models
descrito, el presente trabajo se enfoca en into the Newton-Raphson method for
incorporar modelos de carga power flows that use specified power
estacionarios en el método Newton- restructuring methods. This is
Raphson para flujos de potencia usando developed on Matlab, and the results
la metodología de reestructuración de are compared to DIgSILENT
potencias especificadas. Esto se PowerFactory concerning three
desarrolla en Matlab y se contrastan los different test power systems.
resultados con el software DIgSILENT
PowerFactory sobre tres sistemas de
potencia de prueba.

Palabras Clave: Análisis de errores,

Distribución de potencia, Flujo de carga, Keywords: Error analysis, Power
Método de Newton, Modelado de carga, distribution, Load flow, Newton
Potencia reactiva, Sistema de potencias, method, Load modeling, Reactive
Simulación del sistema de potencia. power, Power systems, Power system
simulation.

1
1 Introducción modelos ya establecidos, pero su
optimización estima dificultades [4], [5]:
Los ingenieros eléctricos generalmente 1. Las cargas involucradas son
toman decisiones en la operación y numerosas.
planificación, apoyados por resultados de 2. El consumo está relacionado con
estudios actuales de energía y de los usuarios, el ingreso es
estabilidad [1]. Dependiendo de la enrevesado debido al dominio y
precisión que se desee dar a los resultados localización.
de un estudio, estos deben tener la mayor 3. El proceder de la carga que presenta
cantidad de modelos de los distintos es voluble frente tiempo, debido a
componentes, los cuales incorporen su necesidad o consumo.
efecto en los resultados. Por ejemplo, las 4. Carece de datos sobre la
cargas son introducidas en el problema de constitución de cargas en el sistema
flujos de potencia de modo tradicional eléctrico.
como elementos que consumen potencia de 5. Ante los diversos cambios en
manera constante; lo cual en forma práctica componentes del sistema, las
no es verdadero, ya que la potencia que cualidades y comportamientos no
consume una carga puede estar sujeto entre se comprimen lo suficiente
otros aspectos al voltaje de alimentación de especialmente cuando ocurren
dicha carga. cambios de alto voltaje y
La expresión carga permite valerse de periodicidad.
diferente modo para cumplir el propósito Las cargas se pueden catalogar de la
del contexto establecido, según la NEC siguiente forma [4]:
(Norma Ecuatoriana de la Construcción), 1. Cargas que contienen
es la potencia instalada o demandada en un características eléctricas y
circuito eléctrico [2]. También se mecánicas “dinámicas rápidas”.
especifica como el intercambio de potencia 2. Cargas cuya respuesta a las
dada o recibida en un punto sobre un excursiones de voltaje exhiben
intervalo de tiempo. discontinuidades significativas.
El firme progreso de la tecnológica 3. Cargas cuya respuesta a las
sobresalta claramente la evolución y variaciones de voltaje no exhibe
composición de las cargas. A pesar de los discontinuidades significativas o
avances notorios, se utilizan modelos retraso de tiempo.
diseñados de hace tiempo atrás. Analizar 4. Carga con características de
un SEP es el requisito fundamental la “dinámica lenta”.
representación precisa de las En trabajos similares [4], [6] se enfocan
características dinámicas y de estado en modelar y monitorear parámetros de
estable de las cargas dentro del sistema [3]. energía agregando modelos de carga para
La potencia del elemento es una de las un rendimiento dinámico de los sistemas a
partes más importantes del sistema de prueba.
energía, y el modelado apropiado genera En [7] se insertan los modelos de carga
objetivos más adecuados. Diferentes para flujos de potencia (FP) usando los
indagaciones científicas evidencian que un métodos alternativos entre ellos Newton
algoritmo de cargas influye Raphson (NR) estado estable, y realizar el
significantemente en el perfil de cálculo del flujo denominado: Nodos de
estabilidad. De este modo perfeccionar inyección, reorganización de potencia
calculada, flujo continuo.

2
 El presente trabajo se enfoca en formar mayor número de tipos de cargas existentes
un método de flujos de potencia por en el sistema para su análisis [4], [8].
Newton-Raphson que pueda resolver flujos En los siguientes estudios es
de potencia usando el método de sumamente importante incluir modelos de
reformulación de las potencias carga, ya que sus resultados pueden
especificadas, para cualquier modelo de asemejarse más a la realidad:
carga en estado estacionario. Este método a. Estabilidad: Transitoria, Pequeña
de flujos de potencia se aplica en diversos señal, Voltaje, Frecuencia
sistemas con distinto número de barras y (desequilibrio carga-generación).
sometido a distintos modelos de carga. b. Flujos de potencia (FP).
El artículo se estructura de la manera
siguiente: Sección 1: introducción y Para entender de mejor manera la
generalidad de cargas. Sección 2: marco connotación de un modelo de carga es
teórico, tipos de cargas, modelos de cargas factible representar la estructura de la carga
y elementos de un flujo de potencia. [9], [10], para lo cual se indican las
Sección 3: Metodología de reformulación siguientes definiciones:
potencias especificadas. La sección 4:
publica la realización del diseño 2.2 Componente de carga
matemático en los distintos sistemas de Un componente de carga es el agregado
potencia. La sección 5: detalla todos los equivalente de todos los dispositivos de un
resultados adquiridos del ensayo, la tipo específico o similar.
sección 6: presenta las conclusiones
alcanzadas por la elaboración del artículo y 2.3 Clase de carga
trabajos a desarrollarse en futuras Las clases de carga para fines de
investigaciones, la sección 7: despliega las simulación es posible categorizarlas como:
referencias adquiridas y finalmente la residencial, comercial o industrial.
sección 8: que son los anexos que es un
soporte fundamental para la elaboración 2.4 Composición de la carga
del presente. La composición fraccional de una carga,
puede aplicarse a la carga de una barra o a
una clase de carga específica.
2 Marco teórico
2.1 Modelos de carga 2.5 Mezcla de clase de carga.
Un modelo se describe como la La composición fraccional de la carga en
representación matemática de un una barra por clases de carga.
componente o fenómeno físico que puede
ser tan sencilla como un número complejo, 2.6 Característica de carga
o tan difícil como una ecuación diferencial. Un conjunto de parámetros, como factor de
Las modelaciones de carga logran potencia, variación de P con V, etc., que
representar matemáticamente una identifican la conducta de una carga
correlación entre la P y Q de la carga unida especificada. Este término puede aplicarse
a un bus, su diferencia de potencial de a un dispositivo de carga, un componente
sistema y frecuencia. de carga, una clase de carga o la carga total
Un modelo de carga debe figurar la de la barra [11].
realidad en todo instante de tiempo y La siguiente terminología se usa
además deben ser idóneos para incluir el comúnmente para representar diferentes
tipos de modelos de carga:

3
 1. Composición de carga típica de
Barra de transmisión categoría de carga.
2. Combinación categórica de carga de
P+jQ cada barra.
Los modelos de carga se dividen en dos
Barra de distribución grandes conjuntos, los cuales son
Clases de carga
dinámicos y estacionarios como se indican
Residencial Comericial Industrial Agricultura en la Figura 3.

Componentes de carga Exponencial Dinámica exponencial

Sistema de Calefacción Refrigeración Aparatos electrónico Iluminación Polinomial

Dinámica IM

Lineal

Figura 1. Clases de carga y componentes [6]. Compuesto

Modelos de Carga
ZIP

Componentes de carga
Motor de inducción Modelos Estacionario Modelos dinámicos Distribución

Iluminación
Potencia constante

Características de carga Transporte de energía

Características de carga
Corriente constante

100%
95% 5% 0% 0% 0%
Impedancia constante Modelo Dess
LM SM Z P I
LM SM Z P I
Figura 3. Categorización de modelados de carga [9].

Refrigeración
Los modelados de carga estacionarios y
Aparatos electrónicos
activos deben facilitar información sobre
Características de carga Características de carga las características de carga relevantes,
0% 0%
generalmente las demandas de potencia
activa y reactiva como también las
LM SM Z P I LM SM Z P I
respuestas ante variaciones de voltaje y
Figura 2. Aportación porcentual de los componentes
frecuencia dentro del sistema.
de carga [12]. La principal diferencia al utilizar los
modelos de carga dinámica es que
En la Figura 2 se observa los
proporcionan esta información en función
componentes de carga y su contribución
del tiempo, mientras que los modelos de
porcentual. Dentro de las clases de carga
carga estática son invariables en el tiempo.
coexisten componentes de carga típicas
Los modelos de carga estacionarios se
que figuran la mayoría de la potencia
emplean para figurar las cargas que
consumida por los beneficiarios finales.
responden instantáneamente a un cambio
Para conseguir esta información se
en el voltaje o frecuencia. Los modelos de
debe tomar en cuenta que difiere para todo
carga dinámica se utilizan para representar
sistema eléctrico, suele ser un
las cargas con una respuesta dependiente
procedimiento confuso que conlleva
del tiempo a un cambio de voltaje o
mucho tiempo y esencialmente se
frecuencia, que se basa en los estados
fundamenta en encuestas a clientes. Para
anteriores de ambos sistemas y carga;
superar estos problemas, se idearon un
además influido por las interacciones de
enfoque más estructurado, simple y fácil de
carga del sistema durante y después de la
realizar, la información necesaria es [13]:

4
transición del estado anterior al siguiente el voltaje también puede surgir como
[5], [11], [14], [12]. admitancia constante. Esto quiere decir
Además, se incorpora varias que, ante un incremento del voltaje,
modelaciones de carga al cálculo de FP; obligatoriamente debe crecer la corriente
siendo este estudio eléctrico en estado para conservar la impedancia [15], [16].
estacionario. A continuación, se realizan
una revisión de los principales modelos 2.9 Modelo de corriente constante
asociados a esta característica no temporal. El modelado es reconocido por el
comportamiento en el que la potencia se
2.7 Modelos en estado estacionario afecta solamente con el voltaje. Esto quiere
Los modelos estacionarios son idóneos indicar que, ante una dilatación del voltaje,
para la representación de cambios sencillos la corriente no se modifica mientras la
y relativamente breves a lo largo del impedancia y energía se alteran [15], [16].
tiempo en cargas del sistema. La
representación habitual del modelo es que 2.10 Modelo de potencia constante
depende de P. reactiva (Q) y real (P) con la Se define como un modelo estático en el
frecuencia y el voltaje para un mismo cual la potencia no se modifica con el
instante de tiempo [15]. voltaje, debido que, para el caso de un
incremento de voltaje, la corriente se
𝑃 = 𝑃(𝑃𝑛 , 𝑄𝑛 , 𝑓𝑛 , 𝑉, 𝑓) (1) reduce con el fin de conservar la misma, lo
que se interpreta como incremento de
𝑄 = 𝑄(𝑃𝑛 , 𝑄𝑛 , 𝑓𝑛 , 𝑉, 𝑓) (2) impedancia. También es llamado como
modelo MVA constante [15], [16].
Donde:
2.11 Modelo ZIP
P Potencia activa al voltaje de El modelo de carga ZIP corresponde al
operación. acoplamiento de distintos modelados de
Pn Potencia activa a voltaje nominal. carga previamente descritos como son:
Q Potencia reactiva al voltaje de potencia, impedancia y corriente constante,
operación. los que se configuran como polinomios de
Qn Potencia reactiva a voltaje nominal. segundo grado como se muestra a
V Voltaje en la barra. continuación [17]:
Vn Voltaje de operación.
2
f Frecuencia del voltaje en la barra. 𝑉 𝑉
fn Frecuencia nominal. 𝑃 = 𝑃𝑛 [𝑝1 ∗ ( ) + 𝑝2 ∗ ( ) + 𝑝3 ] (3)
𝑉𝑛 𝑉𝑛

2
En los FP se considera que la 𝑉 𝑉 (4)
𝑄 = 𝑄𝑛 [𝑞1 ∗ ( ) + 𝑞2 ∗ ( ) + 𝑞3 ]
frecuencia es única e invariante, razón por 𝑉𝑛 𝑉𝑛
la cual, los siguientes modelos desacoplan
a la frecuencia como elemento que actúe
como variador de P y Q que consume una Donde:
carga. p1 y q1 Fracción de la carga de tipo
impedancia constante.
2.8 Modelo de impedancia p2 y q2 Fracción de la carga de tipo
constante impedancia corriente.
Se define como un modelo estático debido p3 y q3 Fracción de la carga de tipo
a que la potencia cambia directamente con impedancia potencia

5
 El modelamiento estacionario es a2 y b2 Elementos añadidos a carga
diseñado de elementos: potencia cte (p3 y con impedancia fija.
q3), componentes de impedancia cte (p1 y
q1), componentes de corriente cte, (p2 y q2). 2.13 Modelo exponencial
Las variables de modelamiento ZIP son la Es una modelación con carga estacionario
fracción de cada modelo de carga dentro de que agrega una dependencia de potencia
la representación ZIP, pueden tomar con su voltaje como un exponencial y la
valores positivos y deben desempeñar las ecuación, habitualmente se compone de la
siguientes limitaciones [18], [6]: siguiente manera:

𝑝1 + 𝑝2 + 𝑝3 = 1 (5) 𝑃 = 𝑃𝑛 ∗ (𝑉)𝑛𝑝 (9)

𝑞1 + 𝑞2 + 𝑞3 = 1 (6) 𝑄 = 𝑄𝑛 ∗ (𝑉)𝑛𝑞 (10)

La contribución de carga de los Donde:

elementos (Z, I, P) en la carga integral varía 𝑉
𝑉= (11)
desde 0 y 1 p.u. tal que su suma total es 1 𝑉𝑛
p.u. La contribución de los elementos La P y Q frente a diferenciaciones de
puede ser mayor a 1 p.u. o menor a 0 p.u, voltaje podría ser apartada e personificada
pero su suma debe ser igual a 1 p.u. a este por un modelamiento exponencial. Los
modelo se lo conoce como “Modelo ZIP parámetros de este modelo son: Pn, Qn y Vn,
restringido”. estos representan los datos iniciales del
modelo de carga y los exponentes, np y nq
2.12 Modelo lineal dependen de la carga y de los datos que se
Cuando el voltaje se modifica en rangos han encontrado mediante ensayos [6].
angostos cercano al nominal, entonces se La conducta veraz en sus cargas no está
usa el modelamiento de carga lineal, sea en asociada puntualmente con cifras
el caso de existir pequeñas perturbaciones particulares como 0, 1 o 2 por motivo que
debería analizar estabilidad. su estructura es compleja y depende de
Al variar elevadamente el voltaje, no se otros componentes que se describe en
recomienda utilizar el modelamiento Tabla 1.
descrito, caso contrario al aplicarlo
provocara resultados erróneos. El Tabla 1. Exponentes de carga [6].
modelamiento de carga lineal se Componente np Nq
representa, tanto como para potencia activa
Aire Acondicionado 0.5 2.5
y potencia reactiva [18] . Freidora, calentador
2.0 0
de agua
𝑉 Lavadora 0.08 1.6
𝑃 = 𝑃𝑛 ∗ [𝑎𝑜 + 𝑎2 ∗ ( )] (7)
𝑉𝑛 Secadora 2.0 3.2
𝑉 Incandescentes 1.55 0
𝑄 = 𝑄𝑛 ∗ [𝑏𝑜 + 𝑏2 ∗ ( )] (8) Fluorescentes 0.96 7.4
𝑉𝑛 Motor industrial 0.07 0.5
Motor del ventilador 0.08 1.6
Donde: Bomba de agua 1.4 1.4
Transformador en
a0 y b0 Intervención relativa de 3.4 11.5
vacío
carga con energía firme.

6
2.14 Modelo polinomial 𝑅2 𝑉2
𝑃 = ( 𝑅1 + )( )
Este modelamiento pertenece al 𝑆
( 𝑅1 +
𝑅2 2
) + (𝑋 1 + 𝑋 2 )2 (14)
𝑆
acoplamiento de modelados de carga,
corriente, impedancia, potencia constante, 2 2
𝑉 𝑉
los que se figuran con polinomios distintos 𝑄 = (𝑋 1 + 𝑋 2 ) ( 𝑅2 2 )+ (15)
(𝑋 1 + 𝑋 2 ) 𝑋
según sea la necesidad del sistema. (𝑅 1
+
𝑆
) +
2

𝑉 𝛼 𝑉 𝛽 𝑉 𝛾
𝑃 = 𝑃𝑛 [𝑝1 ( ) + 𝑝2 ( ) + 𝑝3 ( ) ]
𝑉𝑛 𝑉𝑛 𝑉𝑛
(12) 3 Flujo de potencia (FP)
𝑉 𝛼 𝑉 𝛽 𝑉 𝛾 Los FP son significativos para un diseño o
𝑄 = 𝑄𝑛 [𝑞1 ( ) + 𝑞2 ( ) + 𝑞3 ( ) ]
𝑉𝑛 𝑉𝑛 𝑉𝑛
(13) expansión futura del SEP, así como para
establecer una mejor ejecución del sistema
Donde: existente.
Cuando un sistema de potencia está
𝛼 , 𝛽, 𝛾 Son los índices de trabajando bajo circunstancias específicas
potencias. de generación, carga y topología de red, la
solución del sistema se llama flujo de
2.15 Motor de inducción (MI) potencia [21]. El propósito del resultado
Los modelos de carga estacionario de del FP es formar los valores de ángulo,
motor de inducción o conocidos como MI, magnitud del voltaje, potencia real y
son considerablemente usados en estudios reactiva en la red eléctrica [22]. Así mismo,
de sistemas eléctricos de países con una se obtienen la carga de los componentes,
demanda total de estas máquinas rotativas las pérdidas en potencia siempre para las
que oscila entre el 60% y 70%. peores condiciones que pueda enfrentar el
La fórmula del modelamiento tiene sistema eléctrico.
base en el circuito semejante a un motor de Es importante conocer los modelos de
inducción, expuesto en el transcurso los elementos en un SEP que lo conforman,
reconociendo que los elementos están como lo son: generadores, líneas de
referidos al lado del estator [19], [20]. transmisión, transformadores y
X1+X2
compensación reactiva.
R1

R2/S
X

Figura 4. Diagrama parejo al motor de inducción,

modelo de carga estacionario [7].

Donde:
R1 Resistencia asociada al estator.
R2 Resistencia asociada al rotor. Figura 5. Dirección de un FP de 4 barras en Power
X Reactancia shunt. World.
X1 Reactancia asociada al estator. 3.1 Generadores sincrónicos
X2 Reactancia asociada al rotor. En el SEP el generador síncrono es una
s Deslizamiento de MI. parte importante, transforman la potencia

7
mecánica en potencia eléctrica [23], [24]. c. Capacitancia, C, expresada en F.
Para su indagación, se lo representa al Las líneas de transmisión se clasifican
mismo de acuerdo con la Figura 6, Figura según su longitud, lo cual permite
7, siendo las representaciones SL (Slack), formular modelos de FP.
PV (voltaje controlado) y PQ (carga),
aquellas que modelan un generador. 3.3 Líneas cortas
Se deduce como línea corta la que consta
1
con una longitud hasta de 80 km o menos,
Pi1
i
Logra propagar incluso 1.5 veces su
Pgi
Pik
potencia establecida. La capacitancia en
G k paralelo es despreciable por ser muy
Pin pequeña y solo se valora en serie la
resistencia y la inductancia, que al
multiplicarlo con su distancia resulta la
n
impedancia global, su modelamiento de
Pd1
Figura 6. FP del modelamiento de un generador síncrono línea corta es constituido en la Figura 8
P [MW]. [28].
1
Is R XL Ir
Qi1
i
Qgi
Qik
G k
Vs Vr
Qin

n
Qd1 Figura 8. Modelo de línea corta [19].

Figura 7. FP del modelamiento de un generador síncrono Is R XL Ir

Q [Mvar].

3.2 Líneas de transmisión

La línea de transmisión (LT) es un Vs Y/2 Y/2 Vr
componente más usual que conforma el
SEP. Es el medio donde fluye energía
eléctrica a partir de los ejes de generación
hasta el usuario final. Los sistemas Figura 9. Modelo de línea media π [19].
eléctricos usualmente se encuentran 3.4 Líneas medias
sometidos a cambios constantes Se deduce como una línea media la que
induciendo el incremento de la demanda, constan con una longitud de 80 hasta 240
esto provoca un cambio en el FP que fluye kilómetros, misma que puede transportar
por la LT, el cual conduce un cambio de su potencia nominal entre 1.2 a 1.3 veces.
voltaje que excede su cargabilidad y Este modelo de línea es factible debido a
estimular más allá del rango decretado que representa con exactitud con R y L de
[25], [26]. la línea, si la admitancia se encuentra
Los principales componentes de las LT centrada en un punto medio se denomina
son [27]: circuito T , si la admitancia se fracciona en
a. Resistencia, R, expresada en Ω. partes similares al final de la línea, se
b. Inductancia, L, expresada H.

8
denomina circuito π [26], [28]. Este carga, mientras que la ramificación de
modelo se indica en la Figura 9. magnetización del modelo semejante se
puede ignorar durante el funcionamiento
3.5 Líneas largas
en estado estable. [31], [32], resultando el
Se define como líneas largas a las que
circuito de la Figura 12.
constan con una longitud correspondiente
mayor o igual a 250 km. El modelo de la V2 θ2
V1 θ1
línea larga y muy larga está constituido por
parámetros distribuidos [29], y se muestra
en Figura 10. r + jx
Z

n1:1

Figura 12. Taps del transformador que modifican su

Y /2 Y /2 voltaje de salida [33].
Vs Vr
Donde:
n1 Tap en el lado primario.
Figura 10. Modelamiento de implementación en líneas r+jx Resistencia y reactancia del
muy largas y largas [17]. primario y secundario
referidos al lado primario.
3.6 Transformadores
Los transformadores son elementos Entonces, el transformador se puede
eléctricos con dos devanados plegados obtener por medio el modelo de admitancia
sobre un núcleo magnético como presenta π, como se indica la Figura 13.
la Figura 11 [30].
Una función primordial del
transformador en SEP es rebajar o 1
aumentar su nivel de voltaje, así establecer 𝑦= 2 (𝑟
(16)
condiciones de carga o transmisión. De 𝑛1 + 𝑗𝑥)
manera común en los sistemas eléctricos,
los transformadores tienen cambio de
tomas (taps), esto permite ser posible que V1 θ1 V2 θ2

la Q (Mvar) sea redistribuida para el

mejoramiento en la curva de voltaje de la
red. n1y

PRIMARIO SECUNDARIO (1-n1) y n1 (n1-1) y

Ip FLUJO MAGNETICO Is

Vp
Np Ns
Vs Figura 13. Taps del transformador con modelo π [11].
a.c a.c
3.7 Compensación de potencia
NUCLEO reactiva en paralelo
Para la compensación se ha figurado
Figura 11. Transformador[7].
componentes acoplados al bus de la barra
de forma paralela, tiene como función
En el modelo del transformador doble utilizar inductores o condensadores para
devanado su corriente de excitación es compensar la misma. La representación
inferior que el flujo de electrones en la

9
para esta compensación reactiva es la y todas las admitancias mutuas signo
impedancia acoplada a tierra, y su positivo [35], [33].
admitancia se detalla:
3.9 Solución de flujos de potencia
𝑌𝑠ℎ = (𝑔𝑠ℎ + 𝑗𝑏𝑠ℎ ) (17) Las de flujos de energía resultan en un
sistema de ecuaciones de característica no
Donde: lineal, por lo que es preciso emplear
Ysh Admitancia equivalente. métodos iterativos.
gsh Conductancia equivalente. El método de Newton Raphson (NR) es
bsh Susceptancia equivalente. un algoritmo abierto e iterativo, de
convergencia rápida en forma general
Cabe recalcar la admitancia se puede sobre todo en sus adaptaciones
emplear en el cálculo de la admitancia desacopladas, el método permite encontrar
Ybarra respetando su posición del raíces de funciones mediante un valor
componente de compensación de Q. En numérico próximo a la raíz. Debido a que
algunos casos, la compensación tanto el método Gauss-Seidel posee un método
capacitiva como inductiva se produce en el de convergencia menor comparada con la
mismo bus, por lo que las dos deben ser convergencia cuadrática, que es muy
equivalentes. eficiente. La cifra de iteraciones para que
el método NR con su resultado no tiene
3.8 Matriz Ybarra nada que ver con la dimensión y
La matriz Ybarra es también llamada propiedades del sistema, al contrario para
matriz de admitancia o Ynodo es la cantidad de iteraciones de método
considerada una matriz cuadrada con Gauss-Seidel es directamente proporcional
tamaño (n×n) siendo “n” número de nodos a la dimensión del sistema eléctrico [36],
del sistema eléctrico. Se puede definir [37], [38], [39].
como una matriz simétrica es decir los Este método y sus ramificaciones son
elemento Yij serán iguales a los elementos usados para solucionar flujos de potencia
Yji, donde las “i” son las filas y “j” las en el SEP a gran escala y es parte del
columnas, los elementos de la diagonal son iterativos de métodos numéricos, básico
positivos y se calcula sumando las para calcular aplicaciones de cálculo de
admitancias de los elementos conectados a flujo de potencia. Aunque este método
los nodos, la formación de la matriz tiene varias particularidades de
consiste en [34]: convergencia, hay muchas metodologías
Primordialmente tener el diagrama para mejorarlo, entre ellas las más
unifilar del sistema con las reactancias reconocidas [29]:
referidas al valor base asignado. a. Limitar los aumentos de voltaje y
Las admitancias propias son elementos ángulo en cada iteración.
que son parte de la diagonal principal, se b. Escoger buenos valores de arranque.
deducen sumando todas las admitancias de Generalmente se emplea el arranque
los elementos unidos a la barra mientras las plano 1+j0 para los nodos PQ y V+j0
admitancias mutuas serán las admitancias para los PV.
similares de los elementos existentes entre c. Los controles de reactivos en los nodos
el par de barras. PV son agregados al algoritmo de flujo
Es importante tener en cuenta que las de carga después de alcanzar una
admitancias propias tienen signo negativo convergencia inicial.

10
 Cuando se resuelve por el método NR, 𝐻 𝑁
el estudio del FP las ecuaciones no lineales 𝐽=[ ] (27)
𝑀 𝐿
se transforman en ecuaciones lineales. La
ecuación no lineal utilizada representa la P Para la formulación fuera de la
y Q según su voltaje [40], [41]. diagonal de la matriz del Jacobiano “J”.
Se representa la formula 𝐻𝑖𝑗 = 𝑉𝑖 𝑉𝑗 (𝐺𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 − 𝐵𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 ) (28)
matemáticamente del método NR,
empezando con las potencias calculas [41],
[42]: 𝑁𝑖𝑗 = −𝑉𝑖 𝑉𝑗 (𝐺𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 − 𝐵𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 ) (29)
𝑛
𝑐𝑎𝑙
𝑃𝑖 = 𝑉𝑖 ∑ 𝑉𝑗 (𝐺𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 + 𝐵𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 ) (18) 𝑀𝑖𝑗 = 𝑉𝑖 𝑉𝑗 (𝐺𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 − 𝐵𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 ) (30)
𝑗=1

𝑛
(19) 𝐿𝑖𝑗 = −𝑉𝑖 𝑉𝑗 (𝐺𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 − 𝐵𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 ) (31)
𝑄𝑖 𝑐𝑎𝑙 = 𝑉𝑖 ∑ 𝑉𝑗 (𝐺𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 − 𝐵𝑐𝑜𝑠𝜃𝑖𝑗 ) (39)
𝑗=1
Para la formulación de la matriz del
Jacobiano “J” dentro de la diagonal:
Para el cálculo de los errores de P y Q,
se resta la potencia especificada en la barra
Con los resultados adquiridos de ángulo
menos la potencia calculada:
y amplitud del voltaje, son sujetos a
∆𝑃𝑖 = 𝑃𝑖 𝑒𝑠𝑝 − 𝑃𝑖 𝑐𝑎𝑙 (20) modificaciones por medio de las siguientes
formulaciones:

∆𝑄𝑖 = 𝑄𝑖 𝑒𝑠𝑝 − 𝑄𝑖 𝑐𝑎𝑙 (21) 𝐻𝑖𝑖 = 𝑉𝑖 2 𝐵𝑖𝑖 + 𝑄𝑖 (32)

∆𝑃𝑖 = 𝑃𝑖,𝑒𝑠𝑝 − 𝑉𝑖 ∑ 𝑉𝑗 (𝐺𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 + 𝐵𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 ) (22) 𝑁𝑖𝑖 = −𝑉𝑖 2 𝐺𝑖𝑖 − 𝑃𝑖 (33)
𝑗=1

∆𝑄𝑖 = 𝑄𝑖,𝑒𝑠𝑝 − 𝑉𝑖 ∑ 𝑉𝑗 (𝐺𝑖𝑗 𝑠𝑖𝑛𝜃𝑖𝑗 − 𝐵𝑖𝑗 𝑐𝑜𝑠𝜃𝑖𝑗 ) (23) 𝑀𝑖𝑖 = 𝑉𝑖 2 𝐺𝑖𝑖 − 𝑃𝑖 (34)
𝑗=1

La potencia especificada corresponde

a: 𝐿𝑖𝑖 = 𝑉𝑖 2 𝐵𝑖𝑖 − 𝑄𝑖 (35)

𝑃𝑖 𝑒𝑠𝑝 = 𝑃𝑖 𝑔𝑒𝑛 − 𝑃𝑖 𝑐𝑎𝑟𝑔𝑎 (24)

𝜃 𝐾+1 = 𝜃 𝐾 + ∆𝜃 𝐾 (36)
𝑒𝑠𝑝 𝑔𝑒𝑛 𝑐𝑎𝑟𝑔𝑎
𝑄𝑖 = 𝑄𝑖 − 𝑄𝑖 (25)
∆𝑉 𝐾
Para encontrar la actualización voltajes 𝑉 𝐾+1 = 𝑉 𝐾 (1 + ) (37)
𝑉
y ángulos.
4 Metodología de
[
∆𝜃 ∆𝑃
] = [𝐽]−1 [ ] (26)
reformulación potencias
∆𝑉/𝑉 ∆𝑄 especificadas
La matriz Jacobiana “J” se muestra a Esta metodología se enfoca en incorporar
continuación. en las potencias especificadas tanto activa
como reactiva del método Newton-

11
Raphson los modelos de cargas Inicio

estacionarios (exponencial, polinomial,

lineal, ZIP, motor de inducción, potencia Análisis
de datos
cte , impedancia cte ,corriente cte) en el
cálculo para FP [43], [44]. Por ejemplo, si Cambiar a p.u todos los datos
se incluye el modelo lineal, las ecuaciones Realizar de matriz de admitancias
de las potencias especificadas: Ybarra

𝑔𝑒𝑛 𝑉
𝑃𝑖 𝑒𝑠𝑝 = 𝑃𝑖 − 𝑃𝑛 ∗ [𝑎0 + 𝑎2 ∗ ( )] (38) Calculo de potencias
𝑉𝑛 Especificadas, tomando en cuenta
el efecto de los modelos de carga
en estado estacionario
𝑔𝑒𝑛 𝑉
𝑄𝑖 𝑒𝑠𝑝 = 𝑄𝑖 − 𝑄𝑛 ∗ [𝑏0 + 𝑏2 ∗ ( )] (39)
𝑉𝑛

En la ecuación (38) se analiza las Calculo de errores de potencia

potencias activas especificadas para un

si
modelo de carga lineal y en la ecuación Calculo de
ΔP<Toleracia potencias en Resultados
(39) las potencias reactivas especificadas. ΔQ<Toleracia barras
Las ecuaciones generales para no
establecer potencias especificas se fin
Calculo de Jacobiano
muestran en la fórmula general, que agrega
la influencia de la carga modela. Determinación de las variación de
las variables de estado
𝑃𝑖 𝑒𝑠𝑝 = 𝑃𝑖 𝑔𝑒𝑛 − 𝑃𝑖 𝑚𝑜𝑑𝑒𝑙𝑜 𝑑𝑒 𝑐𝑎𝑟𝑔𝑎 (40)
Actualización de variables de
estado
(41)
𝑄𝑖 𝑒𝑠𝑝 = 𝑄𝑖 𝑔𝑒𝑛 − 𝑄𝑖 𝑚𝑜𝑑𝑒𝑙𝑜 𝑑𝑒 𝑐𝑎𝑟𝑔𝑎
Figura 14. Esquema del flujo, restructuración de
potencias especificadas.
Las ecuaciones (40) y (41) respectivamente 5 Implementación en Matlab
muestran formulaciones para calcular P y
Q especificadas, debido a estos cambios en La implementación de logaritmos
las expresiones se identifica la presentados en el artículo se desarrolló
modificación directa para desarrollar una mediante la plataforma de Matlab, y a cuyo
matriz en base a sus errores de energía. código resultante se le ha llamado
PF_CARGAS. El método se aplica en
∆𝑃𝑖 𝑃𝑖 𝑒𝑠𝑝 − 𝑃𝑖 𝑐𝑎𝑙 distintos sistemas de potencia, estos son
[ ] = [ 𝑒𝑠𝑝 ] tomados de las distintas literaturas. El
∆𝑄𝑖 𝑄𝑖 − 𝑄𝑖 𝑐𝑎𝑙 (42)
primero es tomado de [33], el cual consiste
en un SEP de cuatro barras, dos
En esta metodología, los elementos generadores y cuatro cargas. El segundo
del Jacobiano no sufren cambios [45]. Y modelo IEEE con 9 nodos [47], el cual
deben agregarse completamente las consta de 9 nodos, 3 generadores y 3
modelaciones de energía del estado cargas, finalmente el tercer modelo es
estacionario descritos en la sección 2.7. IEEE con 14 nodos [48], que consta de 5
Para entender de mejor manera el proceso generadores y 11 cargas.
de dicha metodología, se lo detalla en el Para cada SEP se usan varios
esquema de flujo de la Figura 14 [7], [46]. modelos de carga, como se expone en el
apartado de anexos de este trabajo. A
continuación, se encuentran los sistemas de

12
potencia antes descritos con sus 5.2 Sistema de prueba 2: IEEE 9
respectivos datos técnicos cuyos Nodos
parámetros son descritos de forma El modelado de la Figura 16, se indica el
universal de 100 MVA. método tradicional IEEE de 9 nodos, el
cual es empleado de manera común para el
5.1 Sistema de prueba 1: 4 barras desarrollo de estudios eléctricos de distinta
Grainger y Stevenson índole [47], [49].
Nótese la Figura 15, un esquema unifilar
del modelo de prueba 1 [33]. 2 7 Carga 2 9 3

8
G2 G3

carga1 carga 2
5 6
1
2

Carga 1 Carga 3

4 3 1

G1
carga 3
carga 4
Figura 16. Modelo IEEE con 9 nodos [49].
Figura 15. Sistema de 4 barras Grainger [33].
Tabla 5. Demanda de energía del modelo de 9 nodos.
Tabla 2. Demanda de cargas modelo de 4 nodos.
Barra Carga P (MW) Q (Mvar)
Barra Carga P (MW) Q (Mvar)
5 Carga 1 125 50
1 Carga 1 50 30.99 6 Carga 2 90 30
2 Carga 2 170 105.35 8 Carga 3 100 35
3 Carga 3 200 123.94
Tabla 6. Despacho de generadores del sistema de 9
4 Carga 4 80 49.58
barras.
Tabla 3. Despacho de generadores del sistema de 4
barras. Barra Generador P (MW) Q (Mvar)
1 Generador 1 N.A. N.A.
Barra Generador P (MW) Q (Mvar) 2 Generador 2 163.0 1.02
3 Generador 3 85 1.02
1 Generador 1 N.A. N.A.
2 Generador 2 0 0 Tabla 7. Datos de las líneas del modelado de 9 nodos.
3 Generador 3 0 0
4 Generador 4 318 N.A. Barras R (p. u) X (p.u) B (p.u.)

Tabla 4. Valores p.u en líneas del modelo de prueba 1. 4-5 0.01 0.08 0.08
4-6 0.01 0.09 0.07
Barras R (p. u) X (p.u) B (p.u.) 5-7 0.03 0.16 0.15
6-9 0.039 0.17 0.17
1-2 0.010 0.050 0.051 7-8 0.008 0.07 0.07
1-3 0.007 0.037 0.038 8-9 0.01 0.10 0.10
2-4 0.007 0.037 0.038
3-4 0.012 0.063 0.063

13
5.3 Sistema de prueba 3: IEEE 14 Tabla 10. Valores en p.u en líneas del modelamiento con
14 nodos.
nodos
La Figura 17, detalla el esquema unifilar
Barras R (p. u) X (p.u) B (p.u.)
del modelo IEEE 14 nodos. Así como el
modelo IEEE con 9 nodos, este es 1-2 0.01 0.05 0.05
empleado a nivel académico para la 1-5 0.05 0.22 0.04
realización de estudios eléctricos. 2-3 0.04 0.19 0.04
2-4 0.05 0.17 0.03
13
2-5 0.05 0.17 0.03
14
3-4 0.06 0.17 0.03
12
4-5 0.01 0.04 0.01
11
10
6-11 0.09 0.19 0.00
6-12 0.12 0.25 0.00
9
6 6-13 0.06 0.13 0.00
8 G5
9-10 0.03 0.08 0.00
7
4 9-14 0.02 0.27 0.00
G1 10-11 0.08 0.19 0.00
1 12-13 0.22 0.19 0.00
5
13-14 0.17 0.34 0.00

2
6 Análisis de resultados
3
G2 Dado que el modelo propuesto se ha
G3
implementado en diferentes sistemas de
Figura 17. Modelado IEEE 14 nodos [50] prueba para obtener el flujo de energía
incluyendo modelamientos de carga, se
Tabla 8. Demanda de energía del modelo 14 nodos. procede a realizar la comparación de los
valores numéricos con aquellos
Barra Carga P (MW) Q (Mvar) encontrados con el programa DIgSILENT
PowerFactory, para así determinar el
2 Carga 2 21.7 12.7 porcentaje de error tanto en ángulo del
3 Carga 3 94.2 19.0
4 Carga 4 47.8 -3.9
voltaje y magnitud. Los parámetros de
5 Carga 5 7.6 1.6 cada modelo empleado son presentados en
6 Carga 6 11.2 7.5 la sección de anexos.
9 Carga 9 29.5 16.6 Los valores encontrados del modelo de
10 Carga 10 9.0 5.8 prueba permiten considerar y profundizar
11 Carga 11 3.5 1.8 la influencia de distintos métodos de carga
12 Carga 12 6.1 1.6
en el sistema eléctrico de forma teórica y
13 Carga 13 13.5 5.8
14 Carga 14 14.9 5.0 práctica.

Tabla 9. Despacho, generadores modelo de 14 nodos. 6.1 Sistema 1

Tomando en cuenta el esquema unifilar en
Barra Generador P (MW) Q (Mvar) la Figura 15, aquí ostentan los valores
alcanzados del modelo 4 nodos, para lo
1 Generador 1 N.A N.A cual se usan los siguientes modelos de
2 Generador 2 40.0 N.A carga:
3 Generador 3 0.0 N.A a. Carga barra 1: ZIP.
6 Generador 6 0.0 N.A
b. Carga barra 2: lineal.
8 Generador 8 0.0 N.A
c. Carga barra 3: exponencial.

14
 d. Carga barra 4: ZIP.
Los datos de voltaje y ángulo
conseguidos mediante el programa
desarrollado en Matlab y comparado con el
software DIgSILENT PowerFactory se
exponen en la Tabla 11. Como se observa
el modelo propuesto PF_CARGAS tiene
un porcentaje de error tanto en voltaje
como en ángulo del 0%.

Tabla 11. Porcentaje de error en voltajes y ángulos del

sistema 4 barras entre PF_CARGAS y PowerFactory. Figura 18. Porcentaje de error respecto a voltajes en el
sistema de 4 barras.

Pf_cargas PowerFactory Error (%)

Ángulo (°)

Ángulo (°)
V (p. u)

Ángulo
Voltaje
V (p.u)
Barra

1 1 0 1 0 0.00 0.00
2 0.98 -0.95 0.98 -0.95 0.00 0.00
3 0.97 -1.87 0.97 -1.87 0.00 0.00
4 1.02 1.51 1.02 1.51 0.00 0.00

Por medio de los valores resultantes de Figura 19. Porcentaje de error respecto ángulos en el
la Tabla 11, se puede concluir que la sistema de 4 barras.
integración de modelos de carga en estado Como se observa a continuación el
estacionado mediante reformulación de modelado PF_CARGAS tiene un
potencias especificadas no refleja errores porcentaje de error tanto en voltaje como
para el sistema de 4 barras tanto en voltajes en ángulo del 0%.
y ángulos como se observa en la Figura 18, Basado en los valores resultantes de la
Figura 19 al ser contrastados con los Tabla 12, se puede concluir que el
resultados del software comercial. programa desarrollado funciona de manera
adecuada y eficiente para sistemas
6.2 Sistema 2 pequeños y medianos; su porcentaje de
En el SEP de la Figura 16, se muestran los error es nulo de manera ilustrativa se
resultados alcanzados en el sistema de 9 grafica en la Figura 20, Figura 21 el voltaje
barras. Para la simulación se distinguen los en p.u de cada una de las barras y el ángulos
diferentes modelos de carga: de igual forma.
1. Carga A: ZIP.
2. Carga B: lineal. 6.3 Sistema 3
3. Carga C: exponencial. A continuación, se describe los valores
En la Tabla 12 se exponen los valores alcanzados en el modelo de 14 nodos IEEE,
ángulo y voltaje conseguidos mediante el mediante la simulación de distintos
programa desarrollado en Matlab. modelos de carga:

15
 a. Cargas 2, 5, 10, 12 y 14: ZIP. Tabla 12. Porcentaje de error en voltajes y ángulos del
sistema 9 barras entre PF_CARGAS y PowerFactory.
b. Cargas 3, 6, 11: lineal.
c. Cargas 4, 9 y 13: exponencial.
Pf_cargas PowerFactory Error (%)
En la Tabla 13 describe datos de ángulo
y voltajes conseguidos mediante el

Ángulo (°)

Ángulo (°)
programa desarrollado en Matlab y

V (p. u)

Ángulo
Voltaje
V (p.u)
Barra
PowerFactory. En este caso el error en
voltaje oscila entre el 0% hasta el 1.50% y
en ángulos desde 0% hasta 9.96%. 1 1.04 0 1.04 0 0.00 0.00
La aplicación de PF_CARGAS para 2 1.025 9.25 1.02 9.25 0.00 0.00
sistemas grandes y complejos, arroja un 3 1.025 4.63 1.02 4.63 0.00 0.00
4 0.00 0.00
error moderado, como muestra la Tabla 13 1.025 -2.22 1.02 -2.22
0.00 0.00
5 0.99 -3.98 0.99 -3.98
y se ilustra en la Figura 21 y Figura 22. En 6 1.01 -3.71 1.01 -3.71 0.00 0.00
la barra 3 se presenta un porcentaje de error 7 1.02 3.69 1.02 3.69 0.00 0.00
0.00 0.00
notorio en el ángulo, pero nulo en voltaje 8 1.01 0.69 1.01 0.69
9 1.03 1.93 1.03 1.93 0.00 0.00
debido a la implementación del modelo de
carga lineal, mientras que en la barra 4 el
modelo de carga exponencial arroja un
error tanto en voltaje como en ángulo al
igual que la barra 5 con un modelo ZIP.
La implementación de modelamientos
de carga resulta una solución a FP sin
errores en sistemas pequeños, mientras que
en sistemas de mayor cantidad de barras y
cargas aparece un ligero error, el cual se
debe a la característica de convergencia del
método Newton-Raphson, que con la
reestructuración de las potencias
Figura 21. Porcentaje de error respecto ángulos en el
especificadas se modifica el vector de sistema de 9 barras.
errores haciendo que la convergencia sea
en valores levemente distintos a los
entregados por el software comercial.

Figura 22. Porcentaje de error respecto a voltajes en el

sistema de 14 barras.

Figura 20. Porcentaje de error respecto a voltajes en

el sistema de 9 barras.

16
Tabla 13. Porcentaje de error en voltajes y ángulos del La plataforma de Matlab permite la
sistema 14 barras entre PF_CARGAS y PowerFactory.
fácil implementación de la metodología y
análisis, por tal razón, los valores
Pf_cargas PowerFactory Error (%)
alcanzados por el flujo, bajo el modelo de
carga y el modelo sin carga encuentran una
combinación de modelamientos de carga

Ángulo (°)

Ángulo (°)
V (p. u)

Ángulo
Voltaje
V (p.u)
Barra

con el objetivo de aporta una mayor

precisión para los resultados del FP.
Las investigaciones de flujos de
1 1.06 0.00 1.06 0.00 0.00 0.00 potencia no solamente para fines
2 1.04 -5.28 1.04 -5.29 0.00 0.08
académicos sino comerciales o industriales
3 1.01 -14.58 1.01 -13.26 0.00 9.96
4 1.02 -11.12 1.01 -10.95 1.14 1.51
se debe en lo posible emplear modelos de
5 1.03 -9.56 1.01 -9.34 1.50 2.40 carga ya que esto le da una mejor
6 1.07 -15.76 1.07 -15.27 0.00 3.23 representación de la realidad a las
7 1.05 -14.58 1.06 -14.35 0.53 1.56
8
simulaciones.
1.09 -14.58 1.09 -14.35 0.00 1.56
9 1.04 -16.36 1.05 -16.11 0.59 1.53
10 1.04 -16.55 1.04 -16.26 0.48 1.77 7.1 Trabajos futuros
11 1.05 -16.29 1.05 -15.90 0.24 2.42 Con los resultados de este trabajo es
12 1.05 -16.68 1.05 -16.20 0.04 2.96
13 1.048 -16.78 1.04 -16.33 0.09 2.79
posible derivar en investigaciones
14 1.02 -17.55 1.03 -17.23 0.38 1.86 relacionadas a:
1. Estabilidad de voltaje de largo
plazo tomando en cuenta la
dependencia del voltaje
conjuntamente con la carga.
2. Estabilidad en frecuencia
considerando modelamientos de
carga dependientes del voltaje y la
frecuencia.
3. Planificación de la expansión de la
generación y transmisión
considerando modelos de carga el
crecimiento de la demanda.
Figura 23. Porcentaje de error respecto ángulos en el
sistema de 14 barras. 8 Referencias
7 Conclusiones [1] I. N. G. Juan and C. Morataya, “ADIEL
IZAHU MENDEZ ACEVEDO 2002-
La modelación matemática propuesta ha 19926 ANALISIS DE SISTEMA DE
permitido lograr la consolidación de POTENCIA 1,” 2002.
modelamientos de carga para el estado [2] I. Adrián, D. Sandoya, A. Gustavo, and
estacionario permitiendo que el método R. Ordóñez, “Arq . Cristóbal Leonel
NR logre resolver problemas del FP. La Chica Martínez Arq . Jenny Lorena
metodología propuesta modifica solamente Arias Zambrano.”
una parte de las ecuaciones del método [3] J. V. Milanović, K. Yamashita, S.
Newton-Raphson, por lo cual su cálculo no Martínez Villanueva, S. Ž. Djokić, and
L. M. Korunović, “International
involucra una mayor demora en la
industry practice on power system load
convergencia. modeling,” IEEE Trans. Power Syst.,

17
 vol. 28, no. 3, pp. 3038–3046, 2013, Power Networks (566), no. February.
doi: 10.1109/TPWRS.2012.2231969. 2014.
[4] D. L. Lubkeman, S. Member, and A. A. [14] K. Yamashita et al., Modelling and
Girgis, “12, No. 2,” Power, vol. 12, no. aggregation of loads in flexible power
2, pp. 472–482, 1997. networks - Scope and status of the work
[5] D. Q. J. Dqj et al., “Research on of cigre WG C4.605, vol. 8, no. PART
adaptive synthesis dynamic load model 1. IFAC, 2012.
based on multiple model ideology,” vol. [15] G. Olguín Parada, “Modelos De Carga
3, pp. 5–9, 2015. Y Sus Efectos En Estudios,” p. 70,
[6] L. Representation, F. O. R. Dynamic, 2013, [Online]. Available:
and P. Analysis, “Load representation http://repositorio.uchile.cl/bitstream/ha
for dynamic performance analysis ieee ndle/2250/114820/cf-
task force on load representation for quezada_ml.pdf?sequence=1.
dynamic performance,” IEEE Trans. [16] M. Est and M. Din, “Modelo estático y
Power Syst., vol. 8, no. 2, pp. 472–482, dinámico de carga.”
1993, doi: 10.1109/59.260837. [17] Y. Ge, A. J. Flueck, D. K. Kim, J. B.
[7] C. Javier Jiménez Salazar, Ahn, J. D. Lee, and D. Y. Kwon, “An
“Incorporación de Modelos de Carga en Event-Oriented Method for Online
Estado Estable Para el Cálculo de Load Modeling Based on
Flujos de Potencia,” p. 170, 2018, Synchrophasor Data,” IEEE Trans.
[Online]. Available: Smart Grid, vol. 6, no. 4, pp. 2060–
http://bibdigital.epn.edu.ec/bitstream/1 2068, 2015, doi:
5000/19639/1/CD-9042.pdf. 10.1109/TSG.2015.2405920.
[8] Anon, “Load representation for [18] “Walter Brokering, Rodrigo Palma,
dynamic performance analysis,” IEEE Luis Vargas-Ñom Lüfke (El rayo
Trans. Power Syst., vol. 8, no. 2, pp. domado) o los sistemas eléctricos de
472–482, 1993, doi: potencia .pdf.” .
10.1109/59.260837. [19] V. F. Syvokobylenko and S. N.
[9] M. Jin, H. Renmu, and D. J. Hill, Tkachenko, “The mathematical model
“Composite load modeling via of an induction machine in terms of the
measurement approach,” 2006 IEEE skin effect in the rotor and the
Power Eng. Soc. Gen. Meet. PES, p. saturation of magnetic circuits,” 2018
4244, 2006, doi: 10th Int. Conf. Electr. Power Drive
10.1109/pes.2006.1708962. Syst. ICEPDS 2018 - Conf. Proc., no. 2,
[10] A. A. Sallam and O. P. Malik, Power pp. 0–4, 2018, doi:
system stability: Modelling, analysis 10.1109/ICEPDS.2018.8571608.
and control. 2015. [20] M. Rivera and J. De Dios, “Tarea # 2 de
[11] M. A. Merkle and A. M. Miri, conversión de energía,” pp. 1–7.
“Modelling of industrial loads for [21] N. Muñoz, F. Villada, and D. R.
voltage stability studies in power Cadavid, “Análisis de los flujos de
systems,” Can. Conf. Electr. Comput. potencia de sistemas eléctricos y su
Eng., vol. 2, pp. 881–886, 2001, doi: interpretación física,” Inf. Tecnol., vol.
10.1109/ccece.2001.933557. 23, no. 5, pp. 57–68, 2012, doi:
[12] S. Palm and P. Schegner, “Static and 10.4067/S0718-07642012000500007.
transient load models taking account [22] J. Caicedo, “Análisis del flujo de
voltage and frequency dependence,” potencia en líneas de transmisión
19th Power Syst. Comput. Conf. PSCC largas, usando el banco de pruebas para
2016, no. 1, 2016, doi: líneas de transmisión de la UPS-Sede
10.1109/PSCC.2016.7540887. Guayaquil,” 2017.
[13] J. V Milanović et al., Modelling and [23] F. Martínez and P. Andrada,
Aggregation of Loads in Flexible “Generadores síncronos de alta

18
 velocidad. Una evaluación crítica,” Inf. [37] “Flujo de Carga Newton-Raphson
Tecnol., vol. 28, no. 2, pp. 47–58, 2017, Acoplado Rápido con Técnicas para
doi: 10.4067/S0718- Orientar el Análisis en Caso de,” 1998.
07642017000200006. [38] J. A. Terán Tarapués and C. M. Rúa
[24] B. Gomez and M. Javier, “Control Alvarez, “El Método de Newton para
digital de la excitación de un generador raíces complejas. Fractales en el
síncrono,” 2013. problema de Cayley.,” Rev. EIA, vol.
[25] J. Horacio and T. Hernández, 15, no. 29, p. 97, 2018, doi:
“Modelado y Análisis de Sistemas 10.24050/reia.v15i29.1131.
Eléctricos de Potencia en Estado [39] E. L. Metodo, D. E. N. L. A.
Estacionario,” pp. 2–5, 2006. Alternativa, and D. E. L. Ingeniero, “El
[26] S. Estimation, lJ Nagrath. . Metodo De Newton-Raphson - La
[27] L. D. E. Transmision, “Capítulo 9 Alternativa Del Ingeniero Para
Líneas de Transmisión,” pp. 279–320. Resolver Sistemas De Ecuaciones No
[28] I. Samila Mat Zali, Member, IEEE, and Lineales,” Sci. Tech., vol. XI, no. 27,
Jovica V. Milanović, Fellow, “Generic pp. 221–224, 2005.
Model of Active Distribution Network [40] S. Plaza, J. Manuel, and G. Jiménez,
for Large Power System Stability Dinámica del método de Newton. .
Studies,” Samila Mat Zali, Member, [41] M. D. N. Raphson and M. D. N.
IEEE, Jovica V. Milanović, Fellow, Raphson, “Método de Newton
IEEE, 1397. Raphson,” no. October, pp. 0–3, 2018.
[29] U. C. Bandile and K. O. Woodford, [42] W. F. Tinney, S. Member, and C. E.
“ESTUDIO COMPARATIVO DE Hart, “Power Flow Newton ’ s,” Proc.
MÉTODOS PARA EL CÁLCULO IEEE, no. 11, 1967, doi:
DELPUNTO DE CARGA MÁXIMA 10.1177/1053825913518895.
EN UN SISTEMA ELÉCTRICO [43] emadwiandr, Electric Energy Systems
DEPOTENCIA CONSIDERANDO 2nd edition. Exposito Cañizares
COORDENADAS POLARES Conejo, vol. 53, no. 9. 2013.
YRECTANGULARES,” no. [44] E. Dlala and A. Arkkio, “General
December, pp. 2–4, 2015. formulation for the Newton-Raphson
[30] R. M. Mujal Rosas, Cálculo de líneas y method and the fixed-point method in
redes eléctricas. 2002. finite-element programs,” 19th Int.
[31] JOSE ALCON, “ESTUDIO DE LA Conf. Electr. Mach. ICEM 2010, no. 1,
CAPACIDAD DE CARGA SEGURA 2010, doi:
EN TRANSFORMADORES DE 10.1109/ICELMACH.2010.5607928.
POTENCIA,” vol. 1, no., pp. 1–476, [45] S. Kamel, M. Abdel-Akher, and M. K.
2010, doi: El-Nemr, “A new technique to improve
10.1017/CBO9781107415324.004. voltage controlled nodes (PV nodes) in
[32] P. F. Analysis and F. Analysis, “Power the current injection Newton-Raphson
System Analysis,” no. September, power-flow analysis,” Proc. Univ.
2012. Power Eng. Conf., pp. 7–10, 2010.
[33] J. Grainger and W. J. Stevenson, [46] M. Abokrisha, A. Diaa, A. Selim, and
“Analisis de Sistemas de Potencia,” S. Kamel, “Development of Newton-
Cenace, Uanl. p. 743, 1996. raphson power-flow method based on
[34] “vnsnoo GEOVANNY PARDO second order multiplier,” 2017 19th Int.
SALAZAR,” 1998. Middle-East Power Syst. Conf.
[35] C. De, M. Admitancia, and D. Barra, MEPCON 2017 - Proc., vol. 2018-
“Construcción de la Matriz Admitancia Febru, no. December, pp. 976–980,
de Barra del Sistema,” pp. 12–22, 2013. 2018, doi:
[36] B. Curso, D. S. El, and A. E. Zuluaga, 10.1109/MEPCON.2017.8301299.
“Matriz admitancia Y BU S,” 2011. [47] D. Powerfactory, “Nine-bus System,”

19
 pp. 1–10.
[48] D. Powerfactory, “14 Bus System,”
DIgSILENT PowerFactory, pp. 1–9,
2015.
[49] “PSCAD IEEE 09 Bus System,” pp. 0–
4, 2014.
[50] IEEE, “Data Sheets for IEEE 14, 30, 62
Bus System,” 2003.

20
8.1 Matriz de Estado del Arte
FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO ESTACIONARIO
USANDO LA METODOLOGÍA DE REFORMULACIÓN DE POTENCIAS
ESPECIFICADAS
PROPUESTAS
FORMULACIÓN PARA
RESTRICCIONES SOLUCIÓN
DATOS TEMÁTICA DEL RESOLVER
DEL PROBLEMA PROPUESTA
PROBLEMA EL
PROBLEMA

Metodología de Reestructuración
Potencia Instalada o demandada
Perfil de voltaje de acuerdo a los
Simulación de modelos de carga

Implementación de modelos de
Modelos en estado estacionario

Modelos en estado estacionario

Sistemas eléctricos de potencia
NOMBRE DEL ARTICULO

Método de Newton Raphson

Modelación matemática NR
Compensación de potencia

Compensación de potencia

de Potencias Especificadas
Modelos y tipos de carga

limites admisibles
Flujo de Potencia

Modelos de carga

Modelos de carga
reactiva

reactiva
CITAS
ITEM

(SEP)

carga
AÑO

Arq. Cristóbal Leonel Chica Martínez Arq. Jenny Lorena Arias

1 2018
Zambrano.
54 x x x x x
“The mathematical model of an induction machine in terms of the skin
2 2018
effect in the rotor and the saturation of magnetic circuits,”
410 x x x
“Incorporación de Modelos de Carga en Estado Estable Para el
3 2018
Cálculo de Flujos de Potencia,”
31 x x x x x x x x
“El Método de Newton para raíces complejas. Fractales en el
4 2018
problema de Cayley.,”
7 x x x x x
S. Plaza, J. Manuel, and G. Jiménez, Dinámica del método de
5 2018
Newton. .
15 x x x x
M. D. N. Raphson and M. D. N. Raphson, “Método de Newton
6 2018
Raphson,”
111 x x x x x
“Development of Newton-Raphson power-flow method based on
9 2018
second order multiplier,”
32 x x x x x x
M. Rivera and J. De Dios, “Tarea # 2 de conversión de energía,” pp.
7 2017
1–7.
69 x x x x x x x x
J. Caicedo, “Análisis del flujo de potencia en líneas de transmisión
8 2017 largas, usando el banco de pruebas para líneas de transmisión de la 148 x x x x x x x x x x
UPS-Sede Guayaquil,” 2017.
“Generadores síncronos de alta velocidad. Una evaluación crítica,”
10 2017
Inf. Tecnol., vol. 28, no. 2, pp. 47–58, 2017
65 x x x x x x x x

21
 FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO ESTACIONARIO
USANDO LA METODOLOGÍA DE REFORMULACIÓN DE POTENCIAS
ESPECIFICADAS
PROPUESTAS
FORMULACIÓN PARA
RESTRICCIONES SOLUCIÓN
DATOS TEMÁTICA DEL RESOLVER
DEL PROBLEMA PROPUESTA
PROBLEMA EL
PROBLEMA

Metodología de Reestructuración
Potencia Instalada o demandada
Perfil de voltaje de acuerdo a los
Simulación de modelos de carga

Implementación de modelos de
Modelos en estado estacionario

Modelos en estado estacionario

Sistemas eléctricos de potencia
NOMBRE DEL ARTICULO

Método de Newton Raphson

Modelación matemática NR
Compensación de potencia

Compensación de potencia

de Potencias Especificadas
Modelos y tipos de carga

limites admisibles
Flujo de Potencia

Modelos de carga

Modelos de carga
reactiva

reactiva
CITAS
ITEM

(SEP)

carga
AÑO

“Static and transient load models taking account voltage and

11 2016
frequency dependence,”
3 x x x x x x
D. Q. J. Dqj et al., “Research on adaptive synthesis dynamic load
12 2015
model based on multiple model ideology,” vol. 3, pp. 5–9,
17 x x x x x x x x x
A. A. Sallam and O. P. Malik, Power system stability: Modelling,
13 2015
analysis and control.
8 x x x x x x
“An Event-Oriented Method for Online Load Modeling Based on
14 2015
Synchrophasor Data,”
18 x x x x x x x x
U. C. Bandile and K. O. Woodford, “ESTUDIO COMPARATIVO
15 2015 DE MÉTODOS PARA EL CÁLCULO DELPUNTO DE CARGA 23 x x x x x x x x
MÁXIMA…
D. Power factory, “14 Bus System,” Dig SILENT Power Factory,
16 2015
pp. 1–9, 2015.
7 x x x x
J. V Milanović et al., Modelling and Aggregation of Loads in
17 2014
Flexible Power Networks (566), no. February.
150 x x x x x x x
18 2014 “PSCAD IEEE 09 Bus System,” pp. 0–4, 2014. 12 x x x x x x
“International industry practice on power system load modeling,”
19 2013
IEEE Trans. Power Syst
56 x x x x x x x
G. Olguín Parada, “Modelos De Carga Y Sus Efectos En Estudios,” p.
20 2013
70,
16 x x x x x
21 2013 M. Est and M. Din, “Modelo estático y dinámico de carga.” 3 x x x x
“Generic Model of Active Distribution Network for Large Power
22 2013
System Stability Studies,”
2 x x x

22
 FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO ESTACIONARIO
USANDO LA METODOLOGÍA DE REFORMULACIÓN DE POTENCIAS
ESPECIFICADAS
PROPUESTAS
FORMULACIÓN PARA
RESTRICCIONES SOLUCIÓN
DATOS TEMÁTICA DEL RESOLVER
DEL PROBLEMA PROPUESTA
PROBLEMA EL
PROBLEMA

Metodología de Reestructuración
Potencia Instalada o demandada
Perfil de voltaje de acuerdo a los
Simulación de modelos de carga

Implementación de modelos de
Modelos en estado estacionario

Modelos en estado estacionario

Sistemas eléctricos de potencia
NOMBRE DEL ARTICULO

Método de Newton Raphson

Modelación matemática NR
Compensación de potencia

Compensación de potencia

de Potencias Especificadas
Modelos y tipos de carga

limites admisibles
Flujo de Potencia

Modelos de carga

Modelos de carga
reactiva

reactiva
CITAS
ITEM

(SEP)

carga
AÑO

23 2013 emadwiandr, Electric Energy Systems 2nd edition. 1 x x x x

B. Gomez and M. Javier, “Control digital de la excitación de un
24 2013
generador síncrono,” 2013.
8 x x x x
C. De, M. Admitancia, and D. Barra, “Construcción de la Matriz
41 2013
Admitancia de Barra del Sistema,” pp. 12–22, 2013.
10 x x x x x
K. Yamashita et al., Modelling and aggregation of loads in
25 2012
flexible power networks -
59 x x x x x x x x
P. F. Analysis and F. Analysis, “Power System Analysis,” no.
26 2012
September, 2012.
33 x x x x
Analisis de los flujos de potencia de sistemas eléctricos y su
39 2012
interpretación fisica"
140 x x x x x x
B. Curso, D. S. El, and A. E. Zuluaga, “Matriz admitancia Y BU S,”
28 2011
2011.
32 x x x
ESTUDIO DE LA CAPACIDAD DE CARGA SEGURA EN
27 2010
TRANSFORMADORES DE POTENCIA
19 x x x x x x x x
“General formulation for the Newton-Raphson method and the fixed-
29 2010
point method in finite-element programs,”
27 x x x x x
“A new technique to improve voltage-controlled nodes (PV nodes) in
30 2010
the current injection Newton-Raphson power-flow analysis,”
23 x x x x x x x
“Walter Brokering, Rodrigo Palma, Luis Vargas-Ñom Lüfke (El rayo
31 2008
domado) o los sistemas eléctricos de potencia .pdf.” .
327 x x x
32 2006 Composite Load Modeling via Measurement Approach 49 x x x x x x x x

23
 FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO ESTACIONARIO
USANDO LA METODOLOGÍA DE REFORMULACIÓN DE POTENCIAS
ESPECIFICADAS
PROPUESTAS
FORMULACIÓN PARA
RESTRICCIONES SOLUCIÓN
DATOS TEMÁTICA DEL RESOLVER
DEL PROBLEMA PROPUESTA
PROBLEMA EL
PROBLEMA

Metodología de Reestructuración
Potencia Instalada o demandada
Perfil de voltaje de acuerdo a los
Simulación de modelos de carga

Implementación de modelos de
Modelos en estado estacionario

Modelos en estado estacionario

Sistemas eléctricos de potencia
NOMBRE DEL ARTICULO

Método de Newton Raphson

Modelación matemática NR
Compensación de potencia

Compensación de potencia

de Potencias Especificadas
Modelos y tipos de carga

limites admisibles
Flujo de Potencia

Modelos de carga

Modelos de carga
reactiva

reactiva
CITAS
ITEM

(SEP)

carga
AÑO

J. Horacio and T. Hernández, “Modelado y Análisis de Sistemas

33 2006
Eléctricos de Potencia en Estado Estacionario,” pp. 2–5, 2006.
21 x x x x x x
“El Metodo De Newton-Raphson - La Alternativa Del Ingeniero Para
34 2005
Resolver Sistemas De Ecuaciones No Lineales,”
32 x x x x x
35 2003 IEEE, “Data Sheets for IEEE 14, 30, 62 Bus System,” 2003. 39 x x x x
ADIEL IZAHU MENDEZ ACEVEDO 2002-19926 ANALISIS DE
36 2002
SISTEMA DE POTENCIA 1
20 x x x x
37 2002 R. M. Mujal Rosas, Cálculo de líneas y redes eléctricas. 2002. 6 x
M. A. Merkle and A. M. Miri, “Modelling of industrial loads for
38 2001
voltage stability studies in power systems,”
201 x x x x x x x
40 1998 “vnsnoo GEOVANNY PARDO SALAZAR,” 1998. 1 x x x x x
“Flujo de Carga Newton-Raphson Acoplado Rápido con Técnicas
42 1998
para Orientar el Análisis en Caso de,”
17 x x x x
43 1997 D. L. Lubkeman, S. Member, and A. A. Girgis, “12, No. 2,” 3 x x x x x
44 1996 J. Grainger and W. J. Stevenson, “Análisis de Sistemas de Potencia,” 879 x x x x x x x x x x x
45 1995 B. On, “For power,” vol. 10, no. 1, 29 x x x x
“Load representation for dynamic performance analysis ieee task force
46 1993
on load representation for dynamic performance,”
6 x x x x x x x x
47 1967 W. F. Tinney, S. Member, and C. E. Hart, “Power Flow Newton’ s,” 35 x x x x
48 S. Estimation, lJ Nagrath. . 2 x x x x x

24
 50
49
ITEM

AÑO

320.
NOMBRE DEL ARTICULO
DATOS

D. Power factory, “Nine-bus System,” pp. 1–10.

L. D. E. Transmisión, “Capítulo 9 Líneas de Transmisión,” pp. 279–

8
12

CITAS
CANTIDAD
x
x

22

Flujo de Potencia
16

Modelos de carga
16

Modelos en estado estacionario

TEMÁTICA

17

Compensación de potencia
ESPECIFICADAS

reactiva
17

Modelos y tipos de carga

x
x

24

Sistemas eléctricos de potencia

(SEP)
DEL

11

Compensación de potencia
reactiva
PROBLEMA

x
21

Simulación de modelos de carga

FORMULACIÓN

x
19

Modelos de carga

Perfil de voltaje de acuerdo a los

x
x

20

limites admisibles
13

Potencia Instalada o demandada

DEL PROBLEMA
RESTRICCIONES

15

Modelación matemática NR
EL
PARA

Metodología de Reestructuración
19
RESOLVER

PROBLEMA
USANDO LA METODOLOGÍA DE REFORMULACIÓN DE POTENCIAS

de Potencias Especificadas
PROPUESTAS

18

Método de Newton Raphson

x

Implementación de modelos de
17
FLUJOS DE POTENCIA CON MODELOS DE CARGA EN ESTADO ESTACIONARIO

carga
SOLUCIÓN

17

Modelos en estado estacionario

PROPUESTA

25
8.2 Resumen de Indicadores

TEMÁTICA

Flujo de Potencia
25
20
15
10
5
Compensacion de
0 Modelos de carga
potencia reactiva

Modelos en estado
estacionario
Figura 24. Resumen e indicador de temática - matriz estado del arte

FORMULACIÓN DEL PROBLEMA

30

25

20

15

10

0
Modelos y tipos Sistemas Compensación de Simulación de
de carga eléctricos de potencia reactiva modelos de carga
potencia (SEP)
Figura 25. Indicador de formulación del problema-matriz estado del arte

26
 SOLUCIÓN PROPUESTA

18

15

12

0
Método de Newton Implementacion de Modelos en estado
Raphson modelos de carga estacionario
Figura 26. Indicador de solución -matriz estado del arte

27
9 Anexos

Tabla 14. Datos de implementación de modelos de carga en sistema 4 barras Grainger y Stevenson.

INDICES DE MODELOS DE CARGA EN INDICES DE MODELOS DE CARGA EN

DATOS SISTEMA GRAINGER 4 BARRAS
POTENCIA ACTIVA POTENCIA REACTIVA
Angulo
Voltaje

Ycomp
Barra

Tipo

ZQ

PQ
Qg

Qc

ZP

PP
Pg

IQ
Pc

bp

bq
ap

cp

aq

cq
IP
1 1 0 0 0 0.5 0.30 SLACK 0 0.25 0.5 0.25 2 1 0 0.1 0.65 0.25 2 1 0
2 1 0 0 0 1.7 1.05 PQ 0 0.3 0.7 0 0 1 0 0.25 0.75 0 0 1 0
3 1 0 0 0 2 1.23 PQ 0 1 0 0 0.08 1 1 1 0 0 1.05 1 1
4 1.02 0 3.18 0 0.8 0.49 PV 0 0.3 0.35 0.35 2 1 0 0.5 0.15 0.35 2 1 0

Tabla 15. Datos de implementación de modelos de carga en sistema IEEE 9 barras.

INDICES DE MODELOS DE CARGA EN INDICES DE MODELOS DE CARGA EN
DATOS SISTEMA IEEE 9 BARRAS
POTENCIA ACTIVA POTENCIA REACTIVA
Angulo

Ycomp
Voltaje
Barra

Tipo

ZQ

PQ
Qg

Qc

ZP

PP
Pg

IQ
bp

bq
Pc

ap

cp

aq

cq
1 1.04 0 0 0 0 0 SLACK 0 0
IP
0 0 0 0 0 0 0 0 0 0 0
2 1.02 0 1.63 0 0 0 PV 0 0 0 0 0 0 0 0 0 0 0 0 0
3 1.02 0 0.85 0 0 0 PV 0 0 0 0 0 0 0 0 0 0 0 0 0
4 1 0 0 0 0 0 PQ 0 0 0 0 0 0 0 0 0 0 0 0 0
5 1 0 0 0 1.25 0.5 PQ 0 0.4 0.05 0.55 2 1 0 0.7 0.12 0.18 2 1 0
6 1 0 0 0 0.9 0.3 PQ 0 0.7 0.3 0 0 1 0 0.3 0.7 0 0 1 0
7 1 0 0 0 0 0 PQ 0 0 0 0 0 0 0 0 0 0 0 0 0
8 1 0 0 0 1 0.35 PQ 0 1 0 0 0.25 1 1 1 0 0 2.09 1 1
9 1 0 0 0 0 0 PQ 0 0 0 0 0 0 0 0 0 0 0 0 0

28
 Tabla 16. Datos de implementación de modelos de carga en sistema IEEE 14 barras.
INDICES DE MODELOS DE CARGA EN INDICES DE MODELOS DE CARGA EN
DATOS SISTEMA IEEE 14 BARRAS
POTENCIA ACTIVA POTENCIA REACTIVA

Angulo

Ycomp
Voltaje
Barra

Tipo

ZQ

PQ
Qg

Qc

ZP

PP
Pg

IQ
bp

bq
Pc

ap

cp

aq

cq
IP
1 1.06 0 0 0 0 0 SLACK 0 1 0 0 0 0 0 1 0 0 0 0 0
2 1.04 0 0.4 0 0.21 0.12 PV 0 0.11 0.5 0.39 2 1 0 0.5 0.1 0.4 2 1 0
3 1.01 0 0 0 0.94 0.19 PV 0 0.31 0.69 0 0 1 0 0.5 0.5 0 0 1 0
4 1 0 0 0 0.47 -0.03 PQ 0 1 0 0 3.01 1 1 1 0 0 0.09 1 1
5 1 0 0 0 0.07 0.01 PQ 0 0.11 0.5 0.39 2 1 0 0.5 0.1 0.4 2 1 0
6 1.07 0 0 0 0.11 0.07 PV 0 0.31 0.69 0 0 1 0 0.5 0.5 0 0 1 0
7 1 0 0 0 0 0 PQ 0 1 0 0 0 0 0 1 0 0 0 0 0
8 1.09 0 0 0 0 0 PV 0 1 0 0 0 0 0 1 0 0 0 0 0
9 1 0 0 0 0.29 0.16 PQ 0.19 1 0 0 3.01 1 1 1 0 0 0.09 1 1
10 1 0 0 0 0.09 0.05 PQ 0 0.11 0.5 0.39 2 1 0 0.5 0.1 0.4 2 1 0
11 1 0 0 0 0.03 0.01 PQ 0 0.31 0.69 0 0 1 0 0.5 0.5 0 0 1 0
12 1 0 0 0 0.06 0.01 PQ 0 0.11 0.5 0.39 2 1 0 0.5 0.1 0.4 2 1 0
13 1 0 0 0 0.13 0.05 PQ 0 1 0 0 3.01 1 1 1 0 0 0.09 1 1
14 1 0 0 0 0.14 0.05 PQ 0 0.11 0.5 0.39 2 1 0 0.5 0.1 0.4 2 1 0

29