Colegio de Ingenieros Electricistas y Electrónicos

Feria de Ingeniería Eléctrica – 1era. Versión

El uso del Python en la Industria Eléctrica.
Gonzalo Oscar Eulate Choque - Luis Adolfo Pozo Barrera - Daniel Ramos Huancani - Cristhian Coyo Llanque - Limbert Omar Apaza Yahuasi

Abstract. utiliza para el análisis de los sistemas eléctricos de

potencia.
This paper presents a summary of the work
presented at the 1st Technology Fair of the En el software DIgSILENT POWER FACTORY, es
electrical engineering of the San Andrés Major posible la automatización de tareas mediante el uso
University of La Paz Bolivia, in December 2016. del rutinas de programación utilizando un lenguaje
propio del mismo, conocido como DPL (DIgSILENT
Keywords – Script, Digsilent, Sistemas Eléctricos Programming Language) o programación a través
de Potencia (SEP), Python. de software externo, Python para este caso de
estudio, el cual permite acceder a la modificación
de sus parámetros, donde por ejemplo se puede
1. Resumen. cambiar el estado de una línea de transmisión de un
estado fuera de servicio a un estado de servició o
El presente documento, es un resumen del trabajo
en operación, mediante código escrito en el Python
presentado en la 1era Feria de Tecnología de la
[11].
carrera de ingeniería eléctrica de la Universidad
Mayor de San Andrés, realizada en el mes de
De la misma forma se puede acceder a los objetos
diciembre 2016.
de los sistemas eléctricos de potencia mediante
código, para obtener una lista de los elementos
Se presenta el uso del software Python en el análisis
específicos como son: transformadores, líneas,
del crecimiento de la demanda, donde se considera
demandas en nodos, motores, etc., mediante el uso
el ejemplo que dispone PowerFactory, de la red
de comandos en Python.
IEEE de 39 nodos, para un horizonte de 5 años.
De manera general, la utilización de un lenguaje
Así también se presenta el uso del software Python
de programación (Python) coadyuva al especialista
en el ejemplo de la red IEEE 39 nodos, para evaluar
a generar códigos y rutinas personalizadas
las contingencias que afectan a este sistema
que permiten optimizar y estudiar diferentes
eléctrico de potencia (SEP) mediante el criterio
comportamientos en los proyectos eléctricos
determinístico (N-1).
analizados bajo esta lógica.
Asimismo se presenta como un ejemplo de
aplicación al Sistema Interconectado Nacional (SIN) 2. Antecedentes.
de Bolivia, considerando la base datos que dispone
el Comité Nacional de Despacho de Carga para Un sistema eléctrico de potencia es el conjunto
estudios eléctricos del periodo 2016 – 2020, para de elementos que constituyen la red eléctrica de
el análisis de contingencias mediante el criterio potencia siendo su función; generar, transmitir y
determinístico (N-1). distribuir, la energía eléctrica hasta los clientes
finales a costo mínimo, con una determinada
calidad de servicio.
1. Introducción.
En Bolivia la base de datos del Sistema
Los sistemas eléctricos de potencia y las redes
Interconectado Nacional (SIN), es actualizada en
de distribución, se encuentran en constante
forma semestral por el Comité Nacional de Despacho
crecimiento, requiriendo análisis exhaustivos de la
de Carga (CNDC) en el software PowerFactory.
operación en diferentes horarios, bajo determinados
niveles de demanda y bajo contingencias en los
El software DIgSILENT PowerFactory, es una
elementos de los sistemas eléctricos de potencia,
herramienta integrada para el análisis de sistemas
para el conocimiento y comportamiento de las
eléctricos de potencia, caracterizando técnicas
diferentes variables eléctricas.
confiables y flexibles de modelado y algoritmos. El
programa utiliza un ambiente de trabajo muy similar
En la actualidad existen herramientas como el
al que utiliza Windows.
software DIgSILENT POWER FACTORY que se

11
 Colegio de Ingenieros Electricistas y Electrónicos

Una alternativa de automatización de diversas Se presenta el problema del desarrollo

tareas en el software PowerFactory es el uso del
lenguaje de programación Python. Así mismo se
de programas en Python para la
hace notar que se trabajó rutinas con DPL para automatización de tareas repetitivas y
tareas de automatización similares a Python. la presentación de resultados en una
planilla en Excel de forma automática.
3. Planteamiento del Problema.
5. Justificación.
Para la operación óptima de los sistemas eléctricos
de potencia, se requiere contar con la mayor El presente trabajo reducirá el uso de recursos en
información acerca de las variables, como son: las tareas de automatización para el análisis de
voltaje en los nodos del sistema eléctrico de los sistemas eléctricos de potencia en diferentes
potencia, corrientes en las líneas de transmisión y niveles de demanda y contingencias, así como,
transformadores de potencia, pérdidas en líneas de reducir los errores en la actualización de datos, el
transmisión, pérdidas en transformadores, estado tiempo de actualización de datos y proporcionar
de carga en diferentes elementos del sistema mayor cantidad de escenarios para el análisis.
eléctrico de potencia para diferentes niveles de
demanda, etc.
6. Hipótesis.
Para conocer las variables mencionadas, es
necesario realizar tareas de manera repetitiva y Una alternativa para los análisis de los sistemas
automatizada, para la optimización de recursos, eléctricos de potencia en diferentes estados de
estas tareas en el caso del software PowerFactory, carga, condiciones, etc., es “la automatización de
pueden ser realizadas mediante el trabajo en un tareas en el software PowerFactory mediante el uso
script en Python [11]. de programación libre Python”.

Se presenta el problema del desarrollo de

programas en Python para la automatización de
7. Metodología.
tareas repetitivas y la presentación de resultados en
una planilla en Excel de forma automática. Debido a que el proyecto se investigación considera
un desarrollo de software a través de un lenguaje
Los problemas secundarios que se presentan, de programación, se adoptó una metodología
son la selección de la versión de Python y las clásica para su conformación, tomando como base
librerías asociadas para la interfaz con el Excel, las premisas de un proyecto del tipo “espiral” [12],
para la interfaz con cálculos científicos de apoyo en donde:
diferentes análisis, debido a que dependen en gran
medida de la plataforma donde se alojará Python. • A través del diagnóstico de la información
Así mismo, una problemática radica en contar disponible, se desarrollaron diferentes
con un conocimiento previo y manejo del software prototipos, que evolucionaron en su concepción
PowerFactory en el análisis de los sistemas a medida del nivel de programación, integración
eléctricos de potencia. y modificación de los parámetros en DIGSILENT
• Se registró la información inicial para el
ingreso de datos en DIGSILENT a través de la
4. Objetivos. programación en planillas Excel
• Se evaluó los progresos de la programación en
Se presenta a continuación el objetivo general y los
objetivos específicos. Python considerando los resultados de flujos de
potencia mediante una comparación directa con
Objetivo General. los modelos normales.
La automatización de tareas en el software • Se retroalimento los parámetros de entrada
PowerFactory, mediante el uso del lenguaje de con los resultados parciales del prototipo
programación Python. inicial, considerando la necesidad de
información complementaria para cubrir con los
Objetivos Específicos. requerimientos del modelo en DIGSILENT.
• Manejo del software PowerFactory en el análisis
de los sistemas eléctricos de potencia, en 8. Solución Propuesta.
particular flujos de carga.
• Realizar la programación en Python, para
Se trabaja en el entorno del lenguaje de
el intercambio de datos entre el Excel y el
programación Python, debido a que el mismo
PowerFactory.
12
 Colegio de Ingenieros Electricistas y Electrónicos

permite el acceso y modificación de elementos Tabla Nº 1: Datos de las demandas para un

de los sistemas eléctricos de potencia que son horizonte de 5 años.
modelados en el software PowerFactory.

El script en el lenguaje Python puede ser usado en

el PowerFactory para lo siguiente:

• Automatización de tareas.
• Creación de comandos para el uso de funciones.
• Integración del PowerFactory en otras
aplicaciones.

La funcionalidad del PowerFactory es proporcionada Fuente: Elaboración propia.

en Python mediante un módulo dinámico
“PowerFactory” que realiza la interfaz con el API del Los datos que se consideran para el análisis del
DIgSILENT. De esta forma con el uso de scripts de crecimiento de la demanda son los siguientes:
Python se puede acceder a un rango de datos en el
PowerFactory. • Estado de carga de las líneas.
• Estado de carga de los transformadores.
• Todos los objetos. • Voltaje en p.u. en todos los nodos.
• Todos los atributos (Datos de los elementos, • Estado de carga de los autotransformadores.
tipos de datos, resultados). • Potencias generadas, demandas, pérdidas de
• Todos los comandos (cálculo de flujos de carga, potencia por zonas.
etc.).
Se presenta a continuación el diagrama de flujo
• Funciones especiales creados en DPL.
utilizado para la programación en Python, a la
derecha se presenta parte del programa utilizado
Así también scripts creados en Python pueden ser
para los diferentes procesos del diagrama de flujos.
importados por otros scripts en Python, en forma
externa puede ser utilizado el DIgSILENT. Figura N° 1: Diagrama de flujo para el Análisis del
Crecimiento de la Demanda.
El problema de flujos de carga comprende el cálculo
de los flujos de potencia y los voltajes en un sistema
bajo condiciones normales de operación. La
mayor parte de los sistemas de transmisión están
balanceados y una representación monofásica de
la red puede ser usada en este caso.

La solución del flujo de cargas es esencial para las

distintas evaluaciones de los sistemas de potencia
durante los períodos de planificación y operación.

Las alternativas y escenarios son analizados

usando en forma repetitiva el cálculo de flujos de
carga en condiciones normales y de contingencia.

8.1 Análisis de Crecimiento de la

Demanda.
Para el crecimiento de la demanda, se considera
el ejemplo que dispone el PowerFactory, de la red
IEEE de 39 nodos. Se presenta a continuación el
comportamiento de la demanda para un horizonte
de 5 años.

Fuente: Elaboración propia.

13
 Colegio de Ingenieros Electricistas y Electrónicos

Como ejemplo se tiene el caso del proceso (bloque) Tabla N° 2: Casos que se considera la salida de servicio.
EXPORTACIÓN RESULTADOS A EXCEL, en Caso Nombre Nodo i Nodo j
esta parte desde del cálculo de flujos de carga,
1 Line 02 - 25 2 25
los resultados son exportados a Excel en el año
2 Line 03 - 04 3 4
respectivo, para el año 1, los resultados se exportan
a la hoja “1” del archivo denominado RED39. 3 Line 03 - 18 3 18
xlsx. A continuación se incluye el código fuente 4 Line 04 - 05 4 5
correspondiente a esta rutina: 5 Line 04 - 14 4 14
6 Line 08 - 09 8 9
for linea in lineas: #En esta parte se recorre todos 7 Line 09 - 39 9 39
los objetos (líneas) 8 Line 10 - 11 10 11
if linea.outserv == 0: #En esta parte se verifica que 9 Line 10 - 13 10 13
el objeto línea esta en servicio 10 Line 13 - 14 13 14
linea_c = round(linea.GetAttribute(‘c:loading’),2)
#En esta parte se obtiene el estado de carga de la Fuente: Elaboración propia.
línea analizada
ws.Cells(j+1,1).Value=str(linea.loc_name) Los datos que se consideran para el análisis n-1 de
#En esta parte se obtiene el nombre de la línea y contingencias son los siguientes:
se exporta a Excel
ws.Cells(j+1,2).Value=linea_c • Estado de carga de las líneas.
#En esta parte se exporta el estado de carga en • Estado de carga de los transformadores.
Excel • Voltaje en p.u. en todos los nodos.
j=j+1 • Estado de carga de los autotransformadores.
• Potencias generadas, demandas, pérdidas de
En el anterior código de ejemplo, se tiene el caso de potencia por zonas.
obtener los datos de las líneas de transmisión.
Se presenta a continuación el diagrama de flujo
utilizado para la programación en Python:
8.2 Análisis del Sistema Eléctrico
de Potencia frente a Contingencias. Figura N° 2: Diagrama de flujo para el análisis N-1 de contingencias.

En el ejemplo de la red IEEE 39 barras, se utiliza

el criterio N-1 para evaluar las contingencias que
afectan a este sistema eléctrico de potencia (SEP).

Criterio N-1. Un SEP cumple con el criterio del N-1

si al aplicar la contingencia simple más severa, el
sistema continua en condiciones aceptables de
funcionamiento considerando que los flujos en las
líneas se mantienen dentro de límites normales de
operación, los voltajes en las barras no superan
sus niveles mínimos, máximos de variación, no
existen actuaciones de protecciones y no existen
desconexiones forzadas de cargas o equipos.

La aplicación del criterio N-1 consiste en la

simulación de una pérdida de un componente del
SEP que puede ser; una línea de transmisión,
transformador de potencia o un generador.

Para el análisis de contingencias se considera

la salida de servicio simulada para 10 líneas de
transmisión en el SEP de 39 barras, que se presenta
a continuación:

Fuente: Elaboración propia.

14
 Colegio de Ingenieros Electricistas y Electrónicos

Como ejemplo se tiene el caso del proceso (bloque) carga de los transformadores de potencia, tomando
SELECCIÓN LINEA CAMBIAR A FUERA DE en cuenta el crecimiento de la demanda para un
SERVICIO, en esta parte se realilza la selección del horizonte de 5 años.
elemento a simular el estado fuera de servicio, se
realiza el cálculo de flujos de carga, los resultados Tabla N° 4: Estado de carga de los
son exportados a Excel en el caso respectivo, para transformadores de la red 39 nodos.
el año caso 1, los resultados se exportan a la hoja
“1” del archivo RED39N-1.xlsx.

8.3 Aplicación del análisis de

contingencias al SIN.
Como un ejemplo de aplicación al SIN de Bolivia,
se considera la base datos que dispone el Comité
Nacional de Despacho de Carga para estudios Fuente: Elaboración propia.
eléctricos del periodo 2016 – 2020. Se considera
la zona de Oruro como está clasificado en la base En la anterior tabla se observa que se presenta
de Digsilent, en la siguiente tabla se presenta el un problema de sobrecarga en el transformador
resumen de los casos a ser analizados: Trf 06-31, a partir del año 1 en el horizonte de 5
años considerado en el ejemplo. Para el año 5 la
Tabla N° 3: Nombres de las líneas de la zona de sobrecarga está por encima del 200%.
Oruro.
Caso Nombre Nodo i Nodo J VOLTAJE kV LONGITUD km En la siguiente tabla se presentan el estado de carga
1 lne_CSG_BLV069 CSG069 BLV069 69 27,90 de las líneas de transmisión del SEP de 39 nodos,
2 lne_SOC_TNO069 SOC069 TNO069 69 1,81
tomando en cuenta el crecimiento de la demanda
3 lne_CAT_CSG069 CAT069 CSG069 69 3,93
4 lne_TES_VIN069 TES069 VIN069 69 1,02 para un horizonte de 5 años.
5 lne_TET_TES069 TET069 TES069 69 2,64
6 lne_TES_VIN069(1) TES069 VIN069 69 0,001 Tabla N° 5: Estado de carga de las líneas de la red
7 lne_HUY_MIG069 HUY069 MIG069 69 14,46 39 nodos.
8 lne_NOR_TNO069 NOR069 TNO069 69 1,20 Nombre Año Base Año 1 Año 2 Año 3 Año 4 Año 5
9 lne_ANG_MIG069 ANG069 MIG069 69 3,90 Linea % % % % %
10 lne_SUD_TET069 SUD069 TET069 69 5,33 Line 03 - 18 8,75 7,09 6,71 10,6 16,34 23,23
Line 04 - 05 22,88 41,06 59,07 78,93 101,25 127,34
Fuente: CNDC – Elaboración propia. Line 04 - 14 45,75 44,11 43,55 43,29 43,58 44,98
Line 05 - 06 76,22 97,74 121,99 148,62 178,6 214,13
Line 05 - 08 54,18 57,39 63,21 69,75 77,37 86,88
En la anterior se presentan los 10 casos a ser Line 06 - 07 71,88 78,6 88,79 100,17 113,25 129,23
Line 06 - 11 60,67 48,98 37,82 27,83 24,46 36,21
analizados, donde se simula que la línea en cada Line 07 - 08 31,21 40,39 47,64 55,69 64,84 75,81
caso se encuentra fuera de servicio. Line 08 - 09 18,11 15,74 21,83 30,37 41,06 54,71
Line 09 - 39 11,84 12,75 16,88 26,06 37,03 50,77
Line 10 - 11 61,43 50,46 41,15 33,67 32,21 42,68
El Sistema Eléctrico de la zona Oruro comprende Line 10 - 13 47,45 58,01 69,22 81,54 95,41 111,77
66 nodos, 36 líneas y 32 transformadores de Line 13 - 14 45,69 56,99 69,24 82,66 97,7 115,34
Line 17 - 27 7,7 11,33 13,8 16,92 20,64 25,07
potencia, no se considera el comportamiento en Line 21 - 22 99,08 100,5 100,47 100,7 101,33 102,7
los autotransformadores y los generadores en esta
zona, como parte del análisis. Fuente: Elaboración propia.

En la anterior tabla se observa que se presenta

9. Resultados. problemas de sobrecarga en las líneas Line 04-
05, Line 05-06, etc. en el horizonte de 5 años
Se presenta a continuación los resultados del uso considerado en el ejemplo. Nos muestra la posible
del Python en las tareas de automatización para el sobrecarga en varias líneas si no se realizan
análisis del crecimiento de la demanda y el análisis proyectos en el SEP, para cubrir el crecimiento de la
N-1 de contingencias. demanda en el horizonte de 5 años en este ejemplo.

9.1 Resultados del análisis del 9.2 Resultados del análisis N-1 de
crecimiento de la demanda. contingencias.
En la siguiente tabla se presentan el estado de En la siguiente tabla se presentan los estados de

15
 Colegio de Ingenieros Electricistas y Electrónicos

carga de los transformadores del SEP de 39 nodos, Tabla N° 8: Voltajes en p.u. en los nodos de la red 39 nodos.
para los 10 casos analizados. Barra Año Base Caso 1 Caso 2 Caso 3 Caso 4 Caso 5 Caso 6 Caso 7 Caso 8 Caso 9 Caso 10
p.u. p.u. p.u. p.u. p.u. p.u. p.u. p.u. p.u. p.u.
En la siguiente tabla se observa que en ningún Bus 01 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05
caso analizado se presenta sobrecarga en los Bus 02 1,05 1,05 1,06 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05
Bus 24 1,04 1,03 1,04 1,04 1,04 1,04 1,03 1,04 1,04 1,04 1,04
transformadores, sin embargo se puede considerar
Bus 25 1,06 1,04 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06
los casos próximos a un estado de carga del 90%, Bus 26 1,05 1,04 1,06 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05
para la toma de acciones. Bus 35 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05 1,05
Bus 36 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06 1,06
Tabla N° 6: Estado de carga de los Bus 37 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03
Bus 38 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03
transformadores de la red 39 nodos.
Bus 39 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03 1,03
Nombre Xfo Año Base Caso 1 Caso 2 Caso 3 Caso 4 Caso 5 Caso 6 Caso 7 Caso 8 Caso 9 Caso 10
% % % % % % % % % % Fuente: Elaboración propia.
Trf 02 - 30 27,65 26,23 26,03 28,11 27,92 28,36 28,16 27,97 28,08 27,9 27,86
Trf 06 - 31 79,58 80,01 82,31 79,9 79,66 81,03 82,45 81,39 82,1 80,56 81,03 En la anterior tabla se observa, que en el caso 2 que
Trf 10 - 32 86,67 86,79 88,72 86,82 86,9 86,66 88,3 87,71 86,15 86,55 87,19 es el caso en que se asume que la línea Line 03-
Trf 11 - 12 14,18 14,26 14,53 14,23 14,97 13,88 13,7 13,86 40,89 35,68 16,05 04, como el estado de fuera de servicio, es el que
Trf 13 - 12 15,53 15,61 15,59 15,56 14,97 18,35 16,32 16,06 46,73 31,1 19,15 presenta que en 4 nodos la tensión se encuentra
Trf 19 - 20 12,43 12,46 12,42 12,43 12,45 12,42 12,47 12,45 12,45 12,46 12,46 por encima del valor máximo de 1,05 p.u.
Trf 19 - 33 80,41 80,55 80,4 80,44 80,49 80,4 80,57 80,51 80,5 80,54 80,52
Trf 20 - 34 87,98 88,13 87,96 88,01 88,06 87,97 88,15 88,09 88,07 88,12 88,1
Trf 22 - 35 81,46 81,73 81,43 81,52 81,61 81,44 81,76 81,66 81,63 81,72 81,67 9.3 Aplicación del análisis de
Trf 23 - 36
Trf 25 - 37
76,44 76,54
75,06 75,75
76,43
75,1
76,46
75,06
76,5
75,06
76,43
75,06
76,55
75,06
76,51
75,06
76,5
75,06
76,53
75,06
76,52
75,06
contingencias al SIN.
Trf 29 - 38 80,89 81,02 80,87 80,89 80,89 80,89 80,9 80,9 80,89 80,9 80,9
Se utilizó el programa desarrollado en Python que
Fuente: Elaboración propia. tiene nombre SIN_N-1_V1.py, donde se realiza
la lectura de los 10 casos a ser analizados que
En la siguiente tabla se presentan el estado de corresponde a la líneas de la zona de Oruro.
carga de las líneas de transmisión del SEP de 39
nodos, para los 10 casos analizados. Se hace notar que el caso que corresponde a la
simulación de fuera de servicio de la línea lne_CAT_
Tabla N° 7: Estado de carga de las líneas de la red CSG069, el cálculo de flujos de carga no converge
39 nodos. (corresponde al caso N° 3). No se considera los
resultados para este caso, pero da una alarma que
Nombre Año Base Caso 1 Caso 2 Caso 3 Caso 4 Caso 5 Caso 6 Caso 7 Caso 8 Caso 9 Caso 10 si esta línea se encuentra fuera de servicio debe
Linea % % % % % % % % % % ser analizada en detalle la operación del SIN.
Line 04 - 14 45,75 52,08 58,66 49,41 62,87 46,55 46,25 79,4 19,81 11,95
Line 05 - 06 76,22 77,82 82,18 77,43 57,97 103,93 77,55 76,93 41,75 104,12 116,31 En la siguiente tabla se presentan los estados de
Line 05 - 08 54,18 56,61 51,53 54,28 58,39 51,6 59,76 58,29 55,76 53,23 52,87 carga de los transformadores de potencia en la
Line 06 - 07 71,88 74,16 71,08 72,22 71,91 75,09 76,87 75,52 66,85 76,34 78,28 zona de Oruro.
Line 06 - 11 60,67 64,19 63,14 61,92 42,22 90,31 63,3 62,7 20,53 92,43 106,51
Line 07 - 08 31,21 33,53 29,9 31,52 31,24 34,26 34,8 33,82 25,87 35,67 37,64
En la siguiente tabla se observa que los casos con
Line 08 - 09 18,11 18,3 22,27 18,33 18,04 19,53 6,28 21,22 18,9 19,44
mayor estados de sobrecarga en transformadores
Line 09 - 39 11,84 11,89 16,45 12,02 11,61 13,24 20,84 15,4 12,47 13,05
Line 10 - 11 61,43 64,59 63,99 62,58 44,87 88,49 64,37 63,64 108,49 102,95
se presenta para el caso N° 7 analizado (lne_HUY_
Line 10 - 13 47,45 44,43 47,59 46,5 64,5 20,46 46,91 46,76 107,97 9,61 MIG069) y el caso N° 9 (lne_ANG_MIG069).
Line 13 - 14 45,69 42,24 45,21 44,56 64,39 17,8 45,07 44,89 86,52 15,62
Tabla N° 9: Estado de carga de los
Fuente: Elaboración propia. transformadores en la zona de Oruro.
Nombre Xfo Año Base Caso 1 Caso 2 Caso 3 Caso 4 Caso 5 Caso 6 Caso 7 Caso 8 Caso 9 Caso 10
En la anterior tabla se observa, que en el caso 10 % % % % % % % % % %
que es el caso en que se asume que la línea Line trf_ANG069 58,8 58,8 58,66 58,58 58,58 58,59 58,81 58,93 58,91 58,61
13-14, con el estado de fuera de servicio, es el que trf_AVI06901 0 0 4,46 4,48 4,48 4,47 4,47 4,41 4,46 4,47
trf_AVI06902 0 0 114,52 114,96 114,93 114,81 114,8 113,23 114,49 114,67
presenta que tres líneas presentan sobrecarga. trf_CRB069 55,64 55,64 55,68 55,77 55,77 55,76 134,26 55,67 113,29 55,73
trf_CRC069 16,16 16,16 16,22 16,44 16,43 16,41 16,59 16,06 16,44 16,38
En la siguiente tabla se presentan los valores en trf_CRQ069 31,41 31,41 32,02 32,51 32,48 32,32 31,9 31,24 31,76 32,24
p.u. de los voltajes en nodos del SEP de 39 nodos, trf_CSG069 99,48 99,48 102,26 102,62 102,6 102,5 102,49 101,2 102,24 102,38
trf_EST069 102,89 102,89 104,31 107,7 107,56 106,7 104,47 102,41 104,02 103,24
para los 10 casos analizados.
trf_SUD06901 105,07 105,07 107,22 108,92 108,82 108,27 106,78 104,49 106,3 107,99
trf_SUD06902 106 106 108,16 109,88 109,77 109,22 107,72 105,41 107,24 108,94
trf_TES069 64,83 64,83 65,48 68,03 65,26 67,32 65,75 64,57 65,49 65,24
trf_VLC069 71,29 71,29 71,53 71,75 71,74 71,73 74,6 71,11 73,26 71,68

Fuente CNDC – Elaboración propia.

16
 Colegio de Ingenieros Electricistas y Electrónicos

En la siguiente tabla se presentan los estados de de las tareas de la simulación de cambiar el estado
carga en las líneas para los casos analizados en la de las líneas de estado en servicio a estado fuera
zona de Oruro. de servicio en la red de 39 barras, mediante el script
ej39N-1.py en Python.
Tabla N° 10: Estado de carga de las líneas en la
zona de Oruro. Así también se concluye que se cumplieron los
Nombre Xfo Año Base Caso 1 Caso 2 Caso 3 Caso 4 Caso 5 Caso 6 Caso 7 Caso 8 Caso 9 Caso 10 objetivos específicos, tomando en cuenta que se
% % % % % % % % % % realizó una revisión en el manejo y modelamiento
lne_MCH_PAI069 4,85 4,85 8,06 4,22 4,5 5,04 9,41 15,82 10,1 6,54 de sistemas eléctricos de potencia en el DIgSILENT
lne_MCH_SUD069 9,27 9,27 12,51 8,83 9,11 9,39 13,55 20,3 14,26 11,12 PowerFactory, se puede realizar el cálculo de flujos
lne_NOR_TNO069 42,25 42,25 43,38 43,27 43,25 43,14 42,91 42,72 43,05 de carga mediante las opciones que presenta el
lne_SOC_TNO069 59,69 59,69 145,82 142,92 140,24 56,97 36,91 57,95 121,66 PowerFactory o mediante la programación de script
lne_SUD_CRQ069 3 3 3,03 3,05 3,05 3,04 3,02 2,99 3,02 3,04 en Python.
lne_SUD_TET069 35,39 35,39 75,02 15,74 13,87 15,7 40,28 30,22 39,81
De igual forma se cumplió con el objetivo de efectuar
lne_SUD_TNO069 18,3 18,3 43,37 102,61 99,74 97,63 14,22 36,94 15,36 78,68
la programación en Python, para el intercambio de
lne_TAB_TCO069 20,94 20,94 19,89 23,25 23,18 23,14 7,46 20,32 8,76 22,65
datos entre el Excel y el PowerFactory. Para el caso
lne_TAC_TCO069 11,41 11,41 10,27 11,26 11,24 11,53 17,53 11,67 14,03 11,13
del análisis de contingencias se presentan los casos
Fuente CNDC – Elaboración propia. a ser analizados en el archivo RED39N-1.XLSX
en la hoja CASOS, para luego mediante el script
En la tabla anterior se presenta en los casos 4, 5, 6 programado en Python se realiza el cálculo de flujos
y 10, donde se presentan estados de sobrecarga en de carga para el caso analizado y los resultados
el orden de 40%. son exportados a las hojas correspondientes en el
archivo RED39N-1.XLSX.
10. Conclusiones.
Finalmente se realizó el análisis de contingencias
Se concluye en el presente trabajo que se verifica tomando en cuenta la base de datos del SIN,
la automatización de tareas en el software que dispone el Comité Nacional de Despacho de
PowerFactory, mediante el uso del lenguaje de Carga en su página oficial. Se consideró el análisis
programación libre Python. para la zona de Oruro, en este ejemplo
de aplicación de un caso real del
Se realizó el análisis del crecimiento de la demanda
Sistema Eléctrico
para un horizonte de 5 años en el SEP de 39
Nacional de
barras del ejemplo que dispone el DIgSILENT
PowerFactory, para este análisis se utilizó el
script desarrollado en Python ej39GetV2.
py, que se explica en detalle
el funcionamiento
en el numeral
10 del presente
trabajo. La
actualización de
las demandas
y el cálculo de
flujos de carga
para cada año
se realizaron
mediante el
script en Python
de manera
automatizada.

A s i m i s m o
se realizó la
automatización

17
 Colegio de Ingenieros Electricistas y Electrónicos

Bolivia, se concluye que los scripts programados es/~pdo/Seminarios/Python2005.pdf.

en Python funcionaron correctamente, es decir no [11] DIgSILENT GmbH Heinrich-Hertz-StraBe 9
se tuvieron que realizar modificaciones a los scripts “DIgSILENT PowerFactory Advanced Tutorial
utilizados en la red de 39 barras. Python”, Alemania 2016, 42 páginas.

Bibliografía Consultada. Acerca de los autores.

[1] Harper, Enriquez, Análisis Moderno de Sistemas
Eléctricos de Potencia: Soluciones de flujos de Gonzalo Oscar
potencia, Operación económica de sistemas Eulate Choque
de potencia. Impreso en México: MacGraw-Hill/
Limusa, 1981, 574 páginas. Ingeniero Electricista,
[2] Allen J. Wood, Bruce F. Wollemberg, Power Magister Sc en Ingeniería
Generation Operation and Control: Economic del Software, Docente
Dispatch of Thermal Units and Methods of Universidad Mayor de San
solution. Impreso en Estados Unidos de América: Andrés, Trabaja en el Sistema Nuevo de la Empresa
John Wiley & Sons Inc, 1996, 592 páginas. de Electricidad La Paz S.A., DELAPAZ, áreas de
[3] Gonzalo Oscar Eulate Choque, “Análisis de interés: Regulación de Electricidad, Inteligencia
Flujos de Carga mediante el uso de Técnicas Artificial, Sistemas Eléctricos de Distribución,
Heurísticas”, Tesis de Maestría UMSA, La Paz Senior Member IEEE. Estudios de Postgrado
Bolivia, marzo 2013. en Administración de Empresas, Planificación,
[4] Gonzalo Oscar Eulate Choque, “Despacho Protección de Redes de Distribución, Métodos
económico de carga mediante el uso de Numéricos de la Simulación y Modelación, Análisis
Algoritmos Genéticos en el software Digsilent”, de los Regímenes de las Redes de Distribución,
Revista Electromundo No 74, ISSN 2220-5276. Economía Informática, Ingeniería del Software.
La Paz Bolivia, julio 2015.
[5] Dr. Adrian Will “Algoritmos Genéticos y Página web: http://geulate.okicode.com
Optimización Heurística”, apuntes de clases. Email: geulate@gmail.com
Grupo de Aplicaciones de Inteligencia Artificial
Universidad Nacional de Tucumán, 290 páginas.
[6] Digsilent GmbH, “Python Tutorial, Digsilent
Technical Documentaction”, DIgSILENT
PowerFactory, Aplication Guide Power Factory Luis Adolfo
Version 15.2. 2014, 42 páginas. En línea: http:// Pozo Barrera
documents.tips/documents/python-tutorial-
560f11ef2caba.html. Ingeniero Electricista,
[7] 2016 GitHub, Inc., “Python scripting for Universidad Mayor de San
PowerFactory”, Python for PowerFactory, Andrés, con Maestría en
2015, en línea: https://github.com/susantoj/ Gestión de Sistemas de
powerfactory_python. Información, Especialidad
[8] José María Herrera Fernández, Luis Miguel en Tecnologías de Información (TI), Diplomados
Sánchez Brea “Computación científica con en Gestión de Redes y Seguridad Informática,
Python para módulos de evaluación continua en Gerencia de Proyectos Empresariales, Preparación,
asignaturas de ciencias aplicadas”, Esta página gestión y Evaluación de Proyectos de Inversión y
web es un extracto del Proyecto de Innovación Educación Superior. Actualmente trabaja en el Área
y Mejora de la Calidad Docente, en línea: de Planificación, Control de Planes de Inversión y
http://pendientedemigracion.ucm.es/info/aocg/ Sistemas de Información Geográfico de la empresa
python/modulos_cientificos/matplotlib/index. DELAPAZ.
html.
[9] Eugenia Bahit “Curso: Python para Principiantes”, Email: luis.pozo@delapaz.bo
safeCreative 1207302042960. Buenos Aires
Argentina 2012, 136 páginas. En línea: www.
eugeniabahit.com.
[10] Francisco Barranco Espósito, Antonio Guerrero
Galindo, Manuel entrena Casas y Álvaro
González Nonay “Introducción a Python”,
España 2012, 56 páginas. En línea: http://lsi.ugr.

18
 Colegio de Ingenieros Electricistas y Electrónicos

Daniel Ramos
Cristhian Coyo Llanque
Huancani
Estudiante de 9no. semestre
Estudiante último semestre
de la carrera de Ingeniería
de la carrera de Ingeniería
Eléctrica de la Universidad
Eléctrica de la Universidad
Mayor de San Andrés.
Mayor de San Andrés.
Trabajó como Supervisor de
Campo y Supervisor de Digitalización y Codificación
Miembro del equipo ganador de la Primera feria
de Datos en el Instituto Nacional de Estadística, en
de Ingeniería eléctrica. Áreas de interés: Sistemas
varias gestiones. Sus áreas de interés contemplan;
Eléctricos de Distribución, Planificación, Protección
Los Sistemas Eléctricos de Potencia, Automatización
de Sistemas Eléctricos de Potencia y Redes de
de Procesos y Programación. Participó en Talleres
Distribución, programación en Python y Digsilent.
de formación en; Programación con Matlab y la
Estabilidad en los Sistemas Eléctricos de Potencia.
Email: dmjdan@gmail.com
Email: coyosky_3.16@hotmail.com

Limbert Omar Apaza Yahuasi

Estudiante de 10mo. semestre de la carrera de Ingeniería Eléctrica de la

Universidad Mayor de San Andrés. Miembro del equipo ganador de la Primera
feria de Ingeniería eléctrica. Participó en seminarios de PLC’s y Sistemas SCADA,
Sistemas ESCADA en Sistemas de Distribución, Operación y Mantenimiento de
Redes de Distribución, Confiabilidad en los Sistemas Eléctricos de Potencia.

Áreas de interés: Sistemas Eléctricos de Distribución, Planificación, Protección de Sistemas Eléctricos

de Potencia y Redes de Distribución, Informática y Economía.

Email: limbert_17pts@hotmail.com

19