Guzmandolores Matlab
Guzmandolores Matlab
Guzmandolores Matlab
Modalidad
Tesis
Presenta:
Guzmn Dolores Misael.
Asesores:
Dr. Villafuerte Daz Rubn
M.I. ngeles Manual Mario Cruz.
AGRADECIMIENTOS.
"A ti, Dios de mis padres, doy yo gracias y alabo, porque me has
dado sabidura y poder". (Daniel 2:23).
II
FLUJOS DE POTENCIA CON MATLAB
DEDICATORIA.
A mi esposa Isabel:
III
FLUJOS DE POTENCIA CON MATLAB
RESUMEN.
Se incluye adems, las 5 versiones ms recientes del software libre PSAT con
indicaciones para su instalacin y el manual de uso general del software.
IV
FLUJOS DE POTENCIA CON MATLAB
OBJETIVO GENERAL.
Aplicacin del programa MATLAB para determinar los voltajes en una red
elctrica y el flujo de potencia en lneas de transmisin ante cualquier
condicin de demanda de energa.
OBJETIVOS ESPECFICOS.
V
FLUJOS DE POTENCIA CON MATLAB
CONTENIDO.
AGRADECIMIENTOS. .................................................................II
DEDICATORIA. ...........................................................................III
RESUMEN.................................................................................. IV
OBJETIVO GENERAL. ............................................................... V
OBJETIVOS ESPECFICOS. ...................................................... V
CONTENIDO. ............................................................................. VI
INTRODUCCIN. ..................................................................... VIII
VI
FLUJOS DE POTENCIA CON MATLAB
CONCLUSIONES ..................................................................................................79
EJEMPLO 3.4.1 .....................................................................................................79
EJEMPLO 3.4.2 [5] ...............................................................................................80
EJEMPLO 3.4.3 .....................................................................................................80
CONCLUSIONES GENERALES ...........................................................................81
Notas .....................................................................................................................82
Bibliografa: ..........................................................................................................83
VII
FLUJOS DE POTENCIA CON MATLAB
INTRODUCCIN.
VIII
FLUJOS DE POTENCIA CON MATLAB
CAPITULO 1
EL ESTUDIO DE FLUJOS DE POTENCIA.
1.1 INTRODUCCIN.
1.2 ANTECEDENTES.
1
FLUJOS DE POTENCIA CON MATLAB
Ahora con el desarrollo de las computadoras se puede lograr hallar las soluciones del
estudio de flujos de potencia de sistemas complejos, se pueden manejar sistemas de
ms de 200 barras, 300 lneas y 500 transformadores. Los resultados completos son
obtenidos de manera rpida y econmica, simplemente con la impresin de algunas
hojas. Los listados impresos de los resultados que provee la computadora consisten
en una serie de tablas; por lo general, la informacin ms importante a ser
considerada es la tabla, que incluye todos los resultados de la red.
2
FLUJOS DE POTENCIA CON MATLAB
3
FLUJOS DE POTENCIA CON MATLAB
[1]
por John Grainger se estudia la solucin al problema de flujo de potencia por medio
de estos mtodos iterativos.
(1.1)
4
FLUJOS DE POTENCIA CON MATLAB
(1.2)
1
v
i p
0.5
-0.5
-1
(1.3)
5
FLUJOS DE POTENCIA CON MATLAB
v
1
0.5
p
-0.5
-1
(1.4)
6
FLUJOS DE POTENCIA CON MATLAB
v
1
0.5
p
-0.5
-1
7
FLUJOS DE POTENCIA CON MATLAB
( )
(1.5)
[ ( ) ]
(1.6)
| || | ( ) | || |
(1.7)
8
FLUJOS DE POTENCIA CON MATLAB
nunca se hace negativa, mientras que el segundo trmino tiene un valor promedio
cero. Definiendo entonces las siguientes cantidades:
| || | Potencia activa.
| || | Potencia reactiva.
(1.8)
( )
(1.9)
1.6
1.5
1.4
p
p
1.2
1 v
1
i P
0.5 0.8
0.6
0 0.4
0.2
Q
-0.5
0
-0.2
-1
-0.4
En la Figura 1.5 (a) y (b) se observa que la potencia instantnea toma valores
negativos durante ciertos periodos de tiempo, indicando con esto que la energa fluye
en esos momentos de la carga al generador.
9
FLUJOS DE POTENCIA CON MATLAB
| | | |( )
| | | |( )
(1.10)
10
FLUJOS DE POTENCIA CON MATLAB
| | | | | || | ( )
(1.11)
| || | | || |
(1.12)
La relacin que existe entre potencia aparente, reactiva y activa puede ser visto en
forma grfica utilizando lo que se conoce como tringulo de potencia, el cual se
muestra en la siguiente Figura:
11
FLUJOS DE POTENCIA CON MATLAB
| |
(1.13)
En donde representa una medida de la cantidad de potencia til que esta siendo
consumida por el elemento, por lo que al se le conoce como factor de potencia,
el cual al multiplicarlo por la potencia aparente, resulta en la potencia activa que el
elemento consume.
Una vez que se dispone de los modelos de los elementos que componen el SEP,
este debe representarse interconectado de alguna manera los modelos
correspondientes.
12
FLUJOS DE POTENCIA CON MATLAB
una de tales magnitudes. El valor por unidad de una magnitud cualquiera se define
como la razn de su valor al valor base:
(1.14)
El valor base siempre tiene las mismas unidades que el valor real, forzando al valor
unitario a ser adimensional. El valor en por ciento es igual a cien veces el valor por
unidad. Los mtodos de clculo que utilizan los valores por unidad o por ciento son
mucho ms sencillos que aquellos que emplean los valores reales en Volts, Ohms,
KVA, etc.
13
FLUJOS DE POTENCIA CON MATLAB
1
Vase Smbolos Grficos para Diagramas Electrnicos y Elctricos, IEEE Estndar 315-1975
14
FLUJOS DE POTENCIA CON MATLAB
15
FLUJOS DE POTENCIA CON MATLAB
CAPITULO 2
FORMACIN DE LA MATRIZ DE ADMITANCIAS
2.1 INTRODUCCIN.
Para el anlisis de sistemas elctricos es necesario disponer de todos los datos para
llevar a cabo una gran cantidad de estudios que permiten determinar sus condiciones
de operacin tanto en estado estacionario como en estado transitorio. Para ello es
importante conocer la matriz de admitancias de la red, debido a que as es posible,
mediante estudios de flujos de potencia calcular los voltajes de cada nodo de la red,
as como la potencia real y reactiva que circula a travs de los sistemas de
transmisin. Se presentan diferentes formas de calcular la matriz de admitancias de
una red elctrica.
16
FLUJOS DE POTENCIA CON MATLAB
La ley de corrientes establece que: la suma algebraica de las corrientes que entran
a un nodo en un nodo es igual a cero y se expresa matemticamente por medio de
la siguiente ecuacin:
(2.1)
Que tambin puede expresarse como: la suma de las corrientes que entran a un
nodo es igual a la suma de las corrientes que salen del mismo nodo. La ecuacin
es:
(2.2)
Dnde:
nce Corrientes que entran al polo.
ncs Corrientes que salen del polo.
EJEMPLO.
17
FLUJOS DE POTENCIA CON MATLAB
SOLUCION:
Figura 2.2 Nodo uno mostrando las corrientes que entran y salen a l.
En la Figura 2.2, se han dibujado las corrientes que entran y salen del nodo,
haciendo uso de la notacin con doble subndice para indicar que el primero tiene un
potencial mayor que el segundo y que la corriente se supone que fluye en la
direccin mostrada en la Figura del nodo uno con los dems nodos con los cuales
tiene conexin. As, la aplicacin de la ley de corrientes de Kirchhoff aplicada al nodo
uno permite establecer la ecuacin siguiente:
18
FLUJOS DE POTENCIA CON MATLAB
(2.3)
(2.4)
[ ]
(2.5)
Figura 2.3 Nodo dos mostrando las corrientes que entran y salen de l.
19
FLUJOS DE POTENCIA CON MATLAB
Que:
(2.6)
Igual a:
(2.7)
[ ]
(2.8)
Figura 2.4 Nodo tres mostrando las corrientes que entran y salen de l.
20
FLUJOS DE POTENCIA CON MATLAB
(2.9)
(2.10)
[ ]
(2.11)
En las ecuaciones (2.5), (2.8) y (2.11) las variables a determinar son voltajes en los
nodos 1, 2 y 3, por lo que se pueden relacionar matricialmente de la forma siguiente:
[ ] [ ]
[ ]
(2.12)
(2.13)
21
FLUJOS DE POTENCIA CON MATLAB
Dnde:
Matriz de admitancias.
As, para el sistema de tres nodos de la figura 2.1 la ecuacin (2.13) se puede
expresar en forma general de la forma siguiente:
[ ][ ] [ ]
(2.14)
[ ][ ] [ ]
(2.15)
>> YBUS=[-5.5i,2.5i,2i;2.5i,-8i,5i;2i,5i,-8.25i]
YBUS =
22
FLUJOS DE POTENCIA CON MATLAB
>> IBUS=[0.6-1.039230485i;0;0.2778370843-1.575692405i]
IBUS =
0.6000 - 1.0392i
0
0.2778 - 1.5757i
>> VBUS=inv(YBUS)*IBUS
VBUS =
0.9519 + 0.3537i
0.9032 + 0.2981i
0.9691 + 0.3001i
>> radian=57.2955951
23
FLUJOS DE POTENCIA CON MATLAB
24
FLUJOS DE POTENCIA CON MATLAB
(2.16)
Dnde:
Nmero de elementos conectados al nodo .
Impedancia conectada al nodo .
Admitancia propia del nodo .
(2.17)
Dnde:
ndice del nodo de inicio.
Nodo final.
admitancia mutua entre el nodo y el nodo .
Nota: los nodos y deben de ser diferentes al nodo de referencia.
EJEMPLO.
Analizar la tabla 2.1 sin considerar el nodo de referencia (0) y usando el mtodo de
inspeccin de la red, encontrar la matriz de admitancias.
25
FLUJOS DE POTENCIA CON MATLAB
SOLUCION:
El nodo uno est formado por los elementos 1, 2 y 3 de tal manera que:
(2.18)
De igual manera para el nodo dos, la admitancia esta formada por los elementos
2, 4 y 5, y es igual a:
(2.19)
Finalmente para el nodo tres, su admitancia est formada por los elementos 3, 5 y 6,
que es igual a:
(2.20)
(2.21)
Para los elementos tres y siete se tiene P=1 y Q=3, por lo tanto:
26
FLUJOS DE POTENCIA CON MATLAB
(2.22)
Para los elementos seis y ocho se tiene P=2 y Q=3, por lo tanto:
(2.23)
[ ] [ ][ ][ ]
27
FLUJOS DE POTENCIA CON MATLAB
(2.24)
Dnde:
[ ] Matriz de admitancias nodales.
[ ] Matriz de incidencia elemento-bus.
[ ] Transpuesta de la matriz de incidencia elemento-bus.
[ ] Matriz de admitancias primitiva (el inverso de la impedancia de
cada elemento).
EJEMPLO.
SOLUCIN:
Primero, se genera la grfica orientada a partir de la red original (figura 2.6 a). La
orientacin de la red generalmente se hace tomando como punto de partida el nodo
de referencia.
28
FLUJOS DE POTENCIA CON MATLAB
29
FLUJOS DE POTENCIA CON MATLAB
[ ]
(2.25)
30
FLUJOS DE POTENCIA CON MATLAB
[ ]
(2.26)
[ ]
(2.27)
31
FLUJOS DE POTENCIA CON MATLAB
[ ]
[ ] [ ]
(2.28)
[ ]
[ ]
[ ]
(2.29)
32
FLUJOS DE POTENCIA CON MATLAB
CAPITULO 3
SIMULACIONES CON MATLAB.
3.1 INTRODUCCIN.
33
FLUJOS DE POTENCIA CON MATLAB
decir, no existe la posibilidad de modificar sus rutinas y/o agregar nuevos modelos de
dispositivos elctricos, de modo que cuando se investiga dichos tipos de software no
cuentan con la flexibilidad necesaria para simular y evaluar el impacto de nuevas
tecnologas, en desarrollo, y aun no implementadas en los sistemas de potencia
reales; reduciendo as su aplicabilidad en centros de enseanza. Para propsitos
acadmicos y de investigacin es ms importante la flexibilidad que la eficiencia, all
radica la importancia de contar con herramientas computacionales de cdigo abierto
dentro de los cuales se encuentra el software libre PSAT.
MatEMTP X X X X
Matpower X X
PAT X X X X
PSAT X X X X X X X
PST X X X X
SPS X X X X X X
VST X X X X X
34
FLUJOS DE POTENCIA CON MATLAB
3.3 PSAT
PSAT: Sus siglas en ingls corresponden a Power System Analysis Toolbox, que en
espaol significa Caja de Herramientas para el Anlisis de Sistemas de
Potencia.
35
FLUJOS DE POTENCIA CON MATLAB
PSAT constituye una herramienta tan verstil que permite al estudiante responder
preguntas como: Qu pasa si?. Esto es posible porque su interfaz grfica de
usuario permite realizar fcilmente cambios en la topologa de la red o en los
componentes y parmetros elctricos de estos; de este modo se puede aprender a
aplicar los conceptos tericos de forma dinmica mediante prcticas que se acercan
un poco ms a la realidad.
36
FLUJOS DE POTENCIA CON MATLAB
Esta GUI proporciona fcil acceso a todas las herramientas del PSAT. Tambin
proporciona la posibilidad de asignar las configuraciones principales, tales como:
nmero de iteraciones del mtodo NR (Newton Raphson), valor base del sistema,
etc. Adems, PSAT no depende totalmente de esta GUI y usa variables globales
para almacenar las configuraciones de los parmetros de las rutinas y datos para los
modelos; logrando correr igualmente en la versin de lneas de comando.
37
FLUJOS DE POTENCIA CON MATLAB
38
FLUJOS DE POTENCIA CON MATLAB
Figura 3.3 Librera Simulink del PSAT. Se muestran las diferentes categoras para el
diseo de redes.
RUTINAS
Flujo de Potencia
39
FLUJOS DE POTENCIA CON MATLAB
El PSAT incluye una funcin de CPF que es novedad entre los paquetes basados en
MATLAB de anlisis de sistemas de potencia. El algoritmo CPF consiste en un
indicador que calcula un vector tangente normalizado y un corrector, cualquiera
puede ser obtenido por medio de una parametrizacin local o una interseccin
perpendicular.
PSAT es solo un programa Basado en Matlab que proporciona un algoritmo IPM para
resolver OPF basados en mecanismos de equilibrio de mercado. Albergando varias
funciones objetivo:
- La maximizacin del beneficio social
- La maximizacin de la distancia para la mxima condicin de carga
- Un enfoque multi-objetivo
De aqu solo ser examinado el OPF estndar la cual usa la maximizacin del
beneficio social como funcin objetivo.
Sin embargo, cabe sealar que el VSC (Voltage Stability Constrained) OPF que
maximiza la distancia para la mxima condicin de carga y el enfoque multiobjetivo
tambin estn disponibles para realizar anlisis en el PSAT.
PSAT es capaz de realizar anlisis SSS, este tiene la opcin de realizar anlisis
dinmicos y anlisis de sensibilidad QV.
40
FLUJOS DE POTENCIA CON MATLAB
PSAT calcula y grafica los eigenvalores y los factores de participacin del sistema
una vez que el flujo de potencia se ha llevado a cabo. Los eigenvalores son
calculados en el anlisis dinmico de la matriz de estado del sistema dinmico, para
el anlisis de sensibilidad QV son calculados de la matriz jacobiana. Una
caracterstica notable de estas opciones es que los eigenvalores son calculados
usando matrices jacobianas analticas, asegurando resultados con gran precisin.
PSAT provee la opcin de realizar simulaciones en el dominio del tiempo. Para esto
usa dos mtodos de integracin diferentes (regla trapezoidal y Euler retrasado).
PSAT es capaz de introducir perturbaciones comunes por medio de funciones
incrustadas. Estas funciones incrustadas son tiles para simular perturbaciones
comunes para anlisis transitorio tales como fallas y operaciones de interrupcin. Las
perturbaciones pueden ser obtenidas al cambiar los parmetros o valores de las
variables despus de completar el flujo de potencia.
PSAT posee los siguientes modelos estticos y dinmicos para hacer un completo
anlisis de sistemas de potencia:
Otros modelos
Instalacin.
Ejecutando PSAT.
42
FLUJOS DE POTENCIA CON MATLAB
Como se puede ver, aparecer una nueva ventana del entorno de trabajo del PSAT,
dentro de ella podemos ver los datos por defecto para el clculo del flujo de potencia,
la potencia base, la frecuencia, tiempos de anlisis para los anlisis en el dominio del
tiempo, etc. los cuales pueden ser modificados segn nuestro requerimiento de
simulacin solicitados
43
FLUJOS DE POTENCIA CON MATLAB
Cargando datos.
Casi todas las operaciones requieren que un archivo sea cargado. El nombre de este
archivo siempre se muestra en el espacio de edicin de texto del Data File de la
ventana principal, para cargar un archivo simplemente doble clic en este espacio, o
siguiendo la ruta en la barra de herramientas File / Open / Data File o por el atajo
<ctrl+D>. Si la fuente es de un formato diferente soportado por PSAT se puede
utilizar la conversin de formatos que provee este. Tambin es posible cargar
resultados previamente guardados tras utilizar el segundo botn de la izquierda de la
barra de herramientas, en el men File / Open/ Saved System o por el atajo <ctrl+y>.
Corriendo el programa.
Una vez cargado un archivo se debe asignarle una rutina, por ejemplo flujos de
potencia, ya que solo por el mero hecho de cargarlo no implica haber realizado
clculos sobre este. Las ltimas actualizaciones de los archivos son ledas cada vez
que el flujo de potencia es ejecutado. En caso de haber modificado los modelos
Simulink cargados originalmente los datos sern actualizados.
Visualizacin de resultados.
Los resultados pueden ser generalmente visualizados varias formas, mediante una
interfaz grfica de usuario en Matlab o como un archivo de texto en cdigo ASCII, o
como un documento de Excel. Por ejemplo, los resultados del flujo de potencia,
cualquiera que sea la solucin de del sistema actual, puede ser inspeccionado por
medio de una interfaz grfica de usuario que se encuentra en la direccin View
44
FLUJOS DE POTENCIA CON MATLAB
/Static Report o por el atajo <ctrl+v>. Entonces la interfaz permite guardar los
resultados en un archivo de texto.
Guardando resultados.
En cualquier momento el men File / Save / Current System, o por el atajo <ctrl+a>,
puede ser utilizado para guardar el sistema actual en un archivo con extensin .mat.
Todas las estructuras globales usadas por PSAT son almacenadas en este archivo el
cual est ubicado en la carpeta del actual directorio que tiene la extensin .out.
Configuraciones.
El ambiente del Simulink de Matlab y sus caractersticas grficas son utilizadas por
PSAT para crear una herramienta CAD o de dibujo asistido por computador capaz de
disear redes de sistemas de potencia; el Simulink permite visualizar y modificar la
topologa de la red as como almacenar los cambios realizados a esta, sin necesidad
de guardar los cambios directamente desde la lista de archivos.
Por lo tanto las rutinas en el dominio del tiempo no utilizan las capacidades para
construir diagramas de bloque de control que posee el Simulink.
Para acceder a la librera del Simulink desde la ventana principal se hace clic en el
cono de Simulink ubicado en la direccin Edit / Network / Edit Netwok / Simulink
Library o por el atajo <Ctrl+s>.
Herramientas
PSAT viene con varias herramientas, como por ejemplo la conversin de formato de
datos y la posibilidad de que el usuario defina nuevos modelos y/o rutinas. La
conversin de formato de datos permite importar archivos desde otros programas
utilizados para el anlisis de sistemas de potencia. Debe tenerse en cuenta que en
algunos casos la conversin no puede ser realizada desde archivos elaborados con
programas que poseen ms capacidades que PSAT. Los archivos de PSAT pueden
ser convertidos a formatos IEEE.
Interfaces
46
FLUJOS DE POTENCIA CON MATLAB
PSAT provee interfaces para GAMS y UWPFLOW lo cual extiende las capacidades
flujos ptimos de potencia y flujos de potencia continuos. El software General
Algebraic Modeling System (GAMS) cuenta con herramientas tiles para resolver
flujos de potencia ptimos de redes complejas y de grandes dimensiones.
47
FLUJOS DE POTENCIA CON MATLAB
Figura 3.7 Diagrama unifilar en el que se muestran los nombres y nmeros de las barras.
48
FLUJOS DE POTENCIA CON MATLAB
SOLUCIN:
El sistema que muestra la figura 3.7 debe ser inicialmente creado a travs de las
libreras simulink, para ello seguimos la ruta en la ventana principal del PSAT: Edit /
Simulink Library, o dando clic en el icono , o bien por el atajo <Ctrl+S>, despus
de visualizarse la ventana que se muestra en la figura 3.3 vamos a File / New / Model
o por el atajo <Ctrl+N> se nos abrir una nueva ventana para el diseo de las redes,
en esta ultima se debe crear el modelo del ejemplo anterior. Primero se dibujan las
barras del sistema, para ello se apertura el primer smbolo de la figura 3.3 llamado
librera de barras y conectores (ver figura 3.8) y luego se copian o se arrastran los
bloques necesarios.
Como se puede ver esta librera contiene los bloques de las barras y otros
conectores, por ejemplo la librera de flujo de potencia contiene los bloques
necesarios para la construccin de diseos de sistemas de potencia.
49
FLUJOS DE POTENCIA CON MATLAB
Para modificar el nombre de la barra se da clic en el nombre por default, ya que esta
el cursor hay, se borra el nombre actual y se le pone el que se requiera. Para
ingresar los datos de las barras como la tensin, el nmero de estradas y salidas
basta con dar doble clic sobre el smbolo de la barra respectiva.
La figura 3.10 muestra los valores de la barra 1 Abedul que se han introducido de los
datos de la tabla 3.4.
Tabla 3.4 Valores preparados para la entrada fcil de los buses al PSAT.
A B C D E F
1 Abedul 2 2 230 [1.00 0.00] 1 1
2 Olmo 1 2 230 [1.00 0.00] 1 1
3 Pino 2 1 230 [1.00 0.00] 1 1
4 Arce 2 2 230 [1.02 0.00] 1 1
50
FLUJOS DE POTENCIA CON MATLAB
Lo siguiente es agregar las lneas, para ello se abre la librera de flujo de potencia y,
desde all, se arrastra el bloque Transmission Line hasta la ventana CAD, ver figura
3.11. En caso que sea necesario alguna modificacin visual en el diseo, en la
ventana CAD en la pestaa formato de la barra de herramientas se encuentran
opciones como: Flipe Block (Voltear Bloque), Rotate Block (Rotar Bloque), etc.
51
FLUJOS DE POTENCIA CON MATLAB
Para unir los elementos (Lneas con Barras) cada bloque muestra un borne en sus
extremos que sirven para unirse por medio de lneas con otros (figura 3.12). Para
modificar sus parmetros de la misma forma como se hizo con las barras. Todos los
datos mostrados en las tablas 3.2 y 3.3 deben ser ingresados para el modelo.
La figura 3.13 muestra los valores de lnea de 1 a 2 que se han introducido de los
datos de la tabla 3.5.
Tabla 3.5 Valores preparados para la entrada fcil de las lneas al PSAT.
Lnea de
barra a
barra G H I J K L
1-2 [100 230 60] 0 0.01008 0.05040 2*0.05125 [0.0 0.0 0.0]
1-3 [100 230 60] 0 0.00744 0.03720 2*0.03875 [0.0 0.0 0.0]
2-4 [100 230 60] 0 0.00744 0.03720 2*0.03875 [0.0 0.0 0.0]
3-4 [100 230 60] 0 0.01272 0.06360 2*0.06375 [0.0 0.0 0.0]
52
FLUJOS DE POTENCIA CON MATLAB
Lo que prosigue es arrastrar las cargas a nuestro modelo y unirlas a las barras
(figura 3.14). Las potencias activa y reactiva se deben ingresar en por unidad, la
ecuacin 1.14 muestra cmo se obtiene el valor en p.u. dividiendo el valor real entre
el valor base.
La figura 3.15 muestra los valores de carga de la barra 1 que se han introducido de
los datos de la tabla 3.6.
53
FLUJOS DE POTENCIA CON MATLAB
Tabla 3.6 Valores preparados para la entrada fcil de las cargas al PSAT.
Carga de la barra
No. M N O
1 Abedul [100 230] [0.50 0.3099] [1.2 0.8]
2 Olmo [100 230] [1.70 1.0535] [1.2 0.8]
3 Pino [100 230] [2.00 1.2394] [1.2 0.8]
4 Arce [100 230] [0.80 0.4958] [1.2 0.8]
Figura 3.15 Ventana donde se modifican los parmetros de la carga de la barra 1 Abedul.
La figura 3.16 muestra los valores de los generadores de las barra 1 y 2 que se han
introducido de los datos de la tabla 3.7.
54
FLUJOS DE POTENCIA CON MATLAB
Tabla 3.7 Valores preparados para la entrada fcil de los generadores al PSAT.
Generacin de
la barra No. P Q R S T U V
1 Abedul. [100 230] 1.0 0.00 [5 -2] [1.1 0.9] 0 0.85
4 Arce. [100 230] 1.02 ---- [5 -2] [1.1 0.9] 3.18 0.85
Figura 3.16 ventana donde se modifican los parmetros de los elementos Slack y PV de
las barras 1 y 4 respectivamente.
55
FLUJOS DE POTENCIA CON MATLAB
56
FLUJOS DE POTENCIA CON MATLAB
57
FLUJOS DE POTENCIA CON MATLAB
Despus, dar doble clic en la carpeta tests ubicada en Folders in current path y
definimos la extensin del archivo que queramos cargar (para este caso PSAT
Simulink (.mdl)) en el men Filters (ver figura 3.20)
a cargar y modificarlo si es necesario por algn error, para culminar la carga del
archivo hacemos clic en el botn Load (figura 3.21) y si el sistema no tiene ningn
error en sus datos, se cargar en la ventana principal del PSAT tal como se muestra
en la figura 3.22.
59
FLUJOS DE POTENCIA CON MATLAB
60
FLUJOS DE POTENCIA CON MATLAB
Figura 3.24 Reporte esttico donde se muestran datos de tensin, ngulo y potencias
programadas
Tambin se puede tener una grfica comparativa de cada parmetro, haciendo clic
en los iconos de color azul que se encuentran a un lado de las unidades (ver figura
3.25).
61
FLUJOS DE POTENCIA CON MATLAB
Haciendo uso de Excel para una mejor visualizacin de los resultados en tablas, se
da clic en el botn Report o alguna de las formas antes mencionadas y se abrir un
archivo de Excel mostrando en la parte superior los datos de PSAT, los datos del
autor, la localizacin del archivo y la fecha de creacin. Los resultados de flujo de
potencia se organizan en tablas como se muestra a continuacin:
62
FLUJOS DE POTENCIA CON MATLAB
Tambin tenemos como dato extra para que se incluya en nuestro reporte: los
controles de violacin de lmite, basta con seleccionar primer recuadro ubicado bajo
el botn Report para que lo incluya en el reporte, en este caso debido a que el
sistema opera en ptimas condiciones al incluirlo en el reporte nos indica que
63
FLUJOS DE POTENCIA CON MATLAB
Hay que tomar en cuenta que todos los clculos realizado y obtenidos por MATLAB
los podemos ver en cualquier momento, como por ejemplo la matriz o las
matrices de los datos del sistema, en la ventana principal del MATLAB, dentro de la
ventana del WORKSPACE, quedaron almacenados todas las matrices de clculo o
data del sistema, por ejemplo ubicamos la variable Line en ella se encuentra la matriz
hay dos maneras para localizarla: primera, escribiendo en el prompt de MATLAB
el nombre de la variable respetando maysculas y minsculas y se desplegara un
listado de las sub-variables que se encuentran en ella, la matriz compleja esta
denotada por la letra Y, en el caso simulado la matriz que se almacen est
compuesta por 4 filas y 4 columnas al ser un sistema de 4 barras (Y: [4x4]). La
segunda forma es haciendo doble clic sobre la variable Line la cual abrir el editor
de variables, y hay tambin la podemos localizar. Para poder llamar los datos
almacenados en Y correspondientes a la matriz , en el promtp de MATLAB se
digita el nombre de la variable (Line) y su extensin (Y) unidas a travs del punto,
recordando que debemos llamarlo tal como est escrito respetando sus maysculas
y minsculas, es as que se puede obtener los datos o resultados del sistema
simulado.
>> Line.Y
Finalmente se muestra los datos almacenados en esta matriz:
64
FLUJOS DE POTENCIA CON MATLAB
65
FLUJOS DE POTENCIA CON MATLAB
66
FLUJOS DE POTENCIA CON MATLAB
67
FLUJOS DE POTENCIA CON MATLAB
Figura 3.30 Vista previa del sistema de 11 barras para el ejemplo 3.4.2.
Como se puede ver en la figura 3.30, el ejemplo prediseado del sistema de potencia
compuesto de 11 barras corresponde a la figura 3.31 del libro [5]. Cada elemento del
sistema puede ser modificado segn sea la necesidad, como por ejemplo fijar nuevos
datos de tensin y potencia activa en las barras tipo PV, en la barra de referencia fijar
la tensin y su ngulo, modificar la potencia activa y reactiva en las barras de carga,
etc.
68
FLUJOS DE POTENCIA CON MATLAB
Ahora para cargar el archivo y analizarlo, se retorna a la ventana de Load Data File y
dar clic en Load, como se puede ver, si el sistema no tiene ningn error en sus datos,
se cargar en la ventana principal del PSAT. Para el caso de flujo de potencia dar
clic en el botn Power Flow, y si el sistema converge hasta obtener su error se
mostrar una grfica de la convergencia en la parte inferior izquierda, en el caso de
que el sistema reporte algn error de datos, se abrir el grfico prediseado
coloreando los elementos que tengan algn error de datos o conexin.
69
FLUJOS DE POTENCIA CON MATLAB
Para tener una grfica comparativa de cada parmetro, se hace clic en los iconos de
color azul que se encuentran al lado de las unidades.
70
FLUJOS DE POTENCIA CON MATLAB
71
FLUJOS DE POTENCIA CON MATLAB
Para poder ver los clculos realizados y obtenidos como por ejemplo la matriz
se sigue el proceso descrito anteriormente. Para el caso simulado la matriz que se
almacen est compuesta por 11 filas y 11 columnas al ser un sistema de 11 barras
(Y: [11x11]). Ahora se hace el llamado de la matriz escribiendo en el promt de
MATLAB:
>> Line.Y
Finalmente se muestra los datos almacenados en esta matriz:
72
FLUJOS DE POTENCIA CON MATLAB
El flujo de potencia continuo puede ser calculado usando la rutina CPF del PSAT. En
este ejemplo se calcular el CPF para el caso IEEE de 14 barras. El modelo Simulink
para este caso es presentado en la figura 5.27. Este archivo es proporcionado por la
actual distribucin del PSAT y tiene el nombre de d_014.mdl, se encuentra dentro de
la carpeta Tests en el directorio de PSAT.
73
FLUJOS DE POTENCIA CON MATLAB
Para visualizar las curvas de estabilidad de tensin, o de nariz, en cualquier barra del
sistema se cuenta con la opcin Plot, que se encuentra en el GUI principal, dando
esta opcin se despliega una interfaz para grficos donde se podr trazar la curva de
nariz para una barra en cuestin.
A continuacin se muestran las curvas de estabilidad de tensin para las barras 4,9 y
14 del ejemplo en desarrollo.
74
FLUJOS DE POTENCIA CON MATLAB
Despus de haber cargado el archivo dar clic en el botn OPF (ver figura 3.2),
automticamente se correr la rutina. Para conocer resultados tales como potencias
demandas y ofertadas o prdidas u otros resultados despus de la optimizacin, se
procede a lo habitual, se ingresa a la static report y luego se apertura el reporte
detallado de resultados, aquellos resultados se pueden observar en las tablas 3.15 a
3.23.
75
FLUJOS DE POTENCIA CON MATLAB
76
FLUJOS DE POTENCIA CON MATLAB
77
FLUJOS DE POTENCIA CON MATLAB
78
FLUJOS DE POTENCIA CON MATLAB
CONCLUSIONES
EJEMPLO 3.4.1
Observando y comparando las tablas 3.12 y 3.13, que son elaboradas por PSAT a
partir de la tabla 3.10, se hace notorio que los valores de generacin y carga, para
las potencias activa y reactiva son prcticamente los mismos, debido a que los
sistemas de potencia no pueden almacenar energa y por ende toda la potencia que
generan a cada momento debe ser consumida.
Al comparar la tabla 3.14, la cual es construida por PSAT basndose en la tabla 3.11,
con sus dos predecesoras, se puede calcular que el porcentaje de prdidas de
potencias activa y reactiva son inferiores al uno y cinco por ciento, respectivamente.
Hecho que deja ver un sistema con buena eficiencia.
79
FLUJOS DE POTENCIA CON MATLAB
Un detalle importante que se puede ver en las figuras 3.33 y 3.35, es que las barras
5, 6, 8, 10 y 11 no generan ni consumen potencias activas y reactivas; esto se debe
a que ninguna de ellas poseen generadores o cargas asociadas, constituyendo
barras de maniobra. Como es de esperarse la barra Slack (Barra 3) es la que ms
genera potencia activa, y la 2 corresponde a la que ms compensa mediante
reactivos.
Es sencillo inferir que las cargas conectadas a la barra 9 son ms grandes que las
conectadas a la 7 debido a que la primera es la que consume mayores cantidades de
potencias activa y reactiva.
EJEMPLO 3.4.3
CONCLUSIONES GENERALES
81
FLUJOS DE POTENCIA CON MATLAB
Notas
82
FLUJOS DE POTENCIA CON MATLAB
Bibliografa:
Libros:
[1] John J. Grainger, William D. Stevenson Jr., Anlisis de sistemas de potencia
McGraw-Hill.
[2] Enrquez Harper, Introduccin al anlisis de los sistemas elctricos de
potencia, LIMUSA
[3] Charles A. Gross, Anlisis de sistemas de potencia, interamericana.
[4] Dr. Villafuerte Daz Rubn Anlisis de redes
[5] P. Kundur Power System Stability and Control Mc-Graw-Hill, Inc.
[6] Stagg and El-Abian.computer Methods in Power System Analysis McGraw-
Hill Book Company.
[7] Ogata, Katsuhiko. Problemas de Ingeniera de Control Utilizando Matlab
Madrid: Prentice Hall, 1999.
Paginas de internet:
[8] Dr. Federico Milano, Webpage PSAT
http://www.uclm.es/area/gsee/web/Federico/psat.htm
[9] Dr. Federico Milano, Power System Analysis Toolbox documentation for
PSAT version 1.3.4, July 14, 2005.
http://seit.unsw.adfa.edu.au/staff/sites/hrp/research/PSAT/psat-1.3.4.pdf
[10] Dr. Federico Milano, Power System Analysis Toolbox Quick Reference
Manual for PSAT version 2.1.2, June 26, 2008
http://www.uclm.es/area/gsee/web/Federico/archive/psat-2.1.2-ref.pdf
[11] Luigi Vanfretti A Brief Introduction to the Power System Analysis
Toolbox (PSAT) for Power System Analysis Undergraduate Courses
83
FLUJOS DE POTENCIA CON MATLAB
November 9, 2006
http://idisk.me.com/vanfretti/Public/pdfs/2010_LV_PEF_Compilation.pdf
[12] Luigi Vanfretti Introduction to PSAT: A toolbox for Education and
Research in Electric Power Engineering
http://kavir3.persiangig.com/IntroToPSAT.pdf
[13] H. Lotfizad PSAT Power System Analysis Toolbox
kavir3.persiangig.com/PSAT.ppt
[14] Daniel Luna Rivera Manual Bsico de uso del Toolbox Power Systems
Analysis Toolbox Psat
http://es.scribd.com/doc/17037847/Manual-de-Uso-Basico-del-Power-System-
Analysis-Toolbox
[15] Pedro Luis Torres Chacon, Pedro A. Torres C. Power System Analysis
Toolbox uso bsico
http://es.scribd.com/doc/64680278/MATLAB-EN-EL-ANALISIS-Y-
SIMULACION-DE-SISTEMAS-DE-POTENCIA
84