SOLUCION A UN SISTEMA DE POTENCIA DE 10 NODOS
PROYECTO FINAL DE LA MATERIA SISTEMAS DE POTENCIA
Maestra: Alejandra Jiménez Vega
Carrera: Ingeniería en Energías Renovables
Jorge Daniel Jasso Rodriguez
3 de junio del 2023, Tijuana, México
pág. 1
INDICE
Introducción……………………………………………………………. 3
Planteamiento del problema de flujos de potencia………………… 3
Solución Analítica………………………………………………………5
Código en Matlab………………………………………………………16
Simulación……………………………………………………………… 20
Graficas de comparación……………………………………………. 25
Análisis de resultados de forma analítica y simulada………………27
Conclusión………………………………………………………………28
Referencias………………………………………………………………29
pág. 2
INTRODUCCION
En la actualidad, es indispensable contar con energía eléctrica, ya es un servicio básico,
una necesidad, pues es nuestra acompañante todos los días de nuestras vidas. Como
es de tal importancia, siempre se ha estudiado el tema de como generarla mas
eficientemente, sin contaminar, sin dejar impacto ambiental, como transmitirla sin tener
perdidas en el camino, como distribuirla sin tener perdidas en el cambio de tensión, etc.
Se ha intentado resolver todo esto, pero es un tema muy complicado por el hecho de que
siempre habrá perdida de diferentes tipos. Hay muchas personas que no saben cómo
llegar la energía eléctrica en su casa, trabajo, o algún sitio en específico, y este tema y
esta materia nos muestra la tal importancia de esto, lo complejo que es, lo que hemos
avanzado como ser humano para poder “controlar ” la electricidad, es sin dudas algo
increíble. Pero, ¿qué es un sistema de potencia?
Un sistema eléctrico de potencia es un conjunto de dispositivos que convierte energía
de una forma primaria a energía eléctrica, la transporta y la distribuye a los
consumidores finales. Por razones históricas y de eficiencia se adoptó en forma
generalizada el uso de la corriente alterna. Un sistema eléctrico de potencia está
compuesto por tres componentes principales: generación, transmisión y distribución [1]
.
PLANTEAMIENTO DEL PROBLEMA DE FLUJOS DE
POTENCIA
El estudio de flujos de potencia es importante para mantener siempre la estabilidad,
disponibilidad y sobre todo seguridad de los sistemas eléctricos, aunque estos pueden
ser bastante complejos, debido a que conlleva el estudio de a generación de electricidad,
hasta donde se encuentra el usuario consumiendo esa electricidad generada en la
central. Con el fin de poder plantear el problema de un análisis de flujos de potencia, se
demuestra el sistema lo mas simple posible sin perder su esencia, dado que el sistema
puede tener varios buses, distintas fuentes de generación y carga. Esto se debe a que
la demanda al crecimiento de la demanda que continuamente se registra en las redes de
distribución, y los problemas operativos, las compañías suministradoras de energía
deben contar con técnicas apropiadas de análisis para resolver de manera eficiente la
problemática.
Un sistema de potencia consta de N nodos, que como se mencionó anteriormente,
cuenta con nodos que son de unidades generadoras que inyectan potencia al sistema,
pág. 3
esta potencia inyectada es potencia activa (𝑃𝐺 ) y dependiendo de las condiciones de
operación del sistema, se puede inyectar o absorber potencia reactiva al sistema (𝑄𝐺 ), a
otros nodos o buses del sistema en los cuales se conectan las cargas las cuales
absorberán la potencia activa y reactiva.
Figura 1. Diagrama unifilar de un sistema de potencia de 3 nodos.
El planteamiento del problema inicia según el tipo de nodo que sea el que este en el
sistema, que pueden ser los siguientes:
• Bus de carga: Estos buses poseen generación de potencia activa y reactiva, las
cuales son desconocidas y las potencias activas y reactivas de carga son
conocidas (son cargas del sistema). Este tipo de barra se le denomina como P-
Q.
• Bus de voltaje controlado: En este bus el voltaje se mantiene constante, por
eso su nombre. Se le conoce como barra P-V. En cada bus donde existe
conectado un generador, se puede regular la generación de MW por medio del
ajuste de la fuente de energía mecánica y la magnitud de tensión por medio de
la excitación del generador.
• Bus slack: Suele ser un bus al cual no se le asigna potencia activa y reactiva de
generación, entonces esto permite absorber las pérdidas del sistema.
Método para la resolución de este problema
pág. 4
La solución a este tipo de problemas suele ser con un método iterativo, que tenga
ciclos y que, por lo tanto, se pueda programar, ya que, al tener sistemas tan grandes,
con tantos nodos y siendo un problema bastante complejo, nos apoyamos de los
métodos números para poder realizar este proceso. Para ser especifico, se utiliza el
método Newton-Raphson. Este es un método basado en la expresión de funciones de
la serie de Taylor, para funciones de dos o mas variables, que ayudan a resolver las
ecuaciones de flujos de potencia relativamente rápido. Las ecuaciones para flujos de
potencia activa y reactiva, son las siguientes:
Ec. (1)
Ec. (2)
Para el calculo del error, se utilizan las siguientes ecuaciones:
Ec. (3)
Ec. (4)
Estas serian las ecuaciones principales, aunque hay que saber que, para realizar estas
iteraciones, se logra hacer un Jacobiano dependiendo de los nodos que haya en el
sistema, así como la solución de este Jacobiano con derivadas parciales. Resolver este
tipo de problemas a mano, teniendo gran cantidad de buses, se vuelve una tarea
sumamente laboriosa, ahí la importancia de la programación y del desarrollo de
tecnología para resolver el problema de flujos de potencia.
SOLUCION ANALITICA
Cuando haremos la resolución a un problema de flujos de potencia, sabemos que nos
enfrentamos a un gran reto y mas en un tipo de sistema de potencia con es el siguiente:
pág. 5
Figura 2. Sistema de Potencia de 10 nodos
Sus datos son los siguientes:
pág. 6
Tabla 1. Datos de carga generación del Sistema de Potencia
Tabla 2. Datos de red
El ejercicio nos pide:
pág. 7
1. La matriz de admitancias nodales utilizando la matriz de incidencia de la red.
2. La matriz de impedancias nodales por el método directo.
3. Formular las ecuaciones de potencias calculadas para el método de Newton-Raphson.
4. Encontrar los voltajes complejos nodales con tolerancia de 0.00001 p.u. en las
desviaciones de potencia.
5. Determinar los flujos de potencia a través de las líneas de transmisión.
6. Determinar las potencias generadas y las pérdidas en el sistema.
7. Analizar sus resultados.
8. Con la finalidad de mejorar el perfil de voltaje en el sistema, se desea instalar un banco
de capacitores de 30 MVAR a 230 kV nominal. Proponga la ubicación del banco y
determine los nuevos voltajes en los nodos.
Entonces iniciamos con el proceso de la solución analítica. Esto se hará programando
un código en Matlab, y a partir de ahí obtener los resultados. De los datos que tenemos,
iniciamos declarando las impedancias para tener como primer objetivo, obtener Ybus.
Las impedancias del sistema son:
𝑍𝑏𝑢𝑠:
0.0000 + 0.2000𝑖
0.0400 + 0.1500𝑖
0.0300 + 0.1000𝑖
0.0500 + 0.2200𝑖
0.0400 + 0.2000𝑖
0.0900 + 0.2500𝑖
0.0600 + 0.2100𝑖
0.0000 + 0.2500𝑖
0.0400 + 0.1500𝑖
0.0600 + 0.3000𝑖
0.0000 + 0.1500𝑖
0.0300 + 0.1000𝑖
Obtenemos las admitancias del sistema, aplicando la siguiente ecuación a cada
impedancia del sistema:
1
𝑌= 𝐸𝑐. (5)
𝑍
𝑌 =
pág. 8
0.0000 − 5.0000𝑖
1.6598 − 6.2241𝑖
2.7523 − 9.1743𝑖
0.9823 − 4.3222𝑖
0.9615 − 4.8077𝑖
1.2748 − 3.5411𝑖
1.2579 − 4.4025𝑖
0.0000 − 4.0000𝑖
1.6598 − 6.2241𝑖
0.6410 − 3.2051𝑖
0.0000 − 6.6667𝑖
2.7523 − 9.1743𝑖
Una vez obtenido las admitancias del sistema, analizamos la Yshunt de cada nodo, para
al final colocarla en una matriz de este tipo:
𝑌𝑠ℎ𝑢𝑛𝑡
0.0000 + 0.0000𝑖
0.0000 + 0.3700𝑖
0.0000 + 0.0900𝑖
0.0000 + 0.2800𝑖
0.0000 + 0.2600𝑖
0.0000 + 2.0000𝑖
0.0000 + 2.0000𝑖
0.0000 + 0.0000𝑖
0.0000 + 0.0000𝑖
0.0000 + 0.0000𝑖
Una vez tenemos nuestra Yshunt con las admitancias del sistema, sumamos las
admitancias (Y) y las Yshunt de su diagonal. Esto lo declaramos así en Matlab:
𝑌𝑝 = 𝑑𝑖𝑎𝑔([𝑌 ; 𝑌𝑠ℎ𝑢𝑛𝑡])
Realizando esta operación en Matlab, solo tenemos que obtener nuestra matriz de
incidencia para obtener Ybus. La matriz de incidencia, es:
pág. 9
Figura 3. Matriz de incidencia
Para obtener nuestra Ybus, aplicamos la siguiente formula:
𝑌𝑏𝑢𝑠 = 𝐴 ∗ 𝑌𝑝 ∗ 𝐴′ 𝐸𝑐. (6)
Dándonos como resultado:
Columna 1 al 10 Fila del 1 al 5
𝑌𝑏𝑢𝑠
1.6598 − 11.2241𝑖 5𝑖 0 0 0
5𝑖 4.6962 − 22.9342𝑖 −2.7523 + 9.1743𝑖 −0.9823 + 4.3222𝑖 −0.9615 + 4.8077𝑖
0 −2.7523 + 9.1437𝑖 4.0271 − 12.6254𝑖 0 0
0 −0.9823 + 4.3222𝑖 0 2.2402 − 12.4447𝑖 −1.2579 + 4.4025𝑖
0 −0.9615 + 4.8077𝑖 0 −1.2579 + 4.4025𝑖 3.8792 − 15.1743𝑖
0 0 −1.2748 + 3.5411𝑖 0 0
0 0 0 0 −1.6598 + 6.2241𝑖0
0 0 0 0 0
−1.6598 + 6.2241𝑖 0 0 0 0
0 0 0 4𝑖 0
Columna del 1 a 10 Fila del 6 al 10
pág. 10
0 0 0 −1.6598 + 6.2241𝑖 0
0 0 0 0 0
−1.2748 + 3.5411𝑖 0 0 0 0
0 0 0 0 4𝑖
0 −1.6598 + 6.2241𝑖 0 0 0
1.9158 − 4.7462𝑖 −0.6410 + 3.2051𝑖 0 0 0
−0.6410 + 3.2051𝑖 2.3008 − 14.0959𝑖 6.6667𝑖 0 0
0 6.6667𝑖 −6.6667𝑖 0 0
0 0 0 4.4120 − 15.3984𝑖 −2.7523 + 9.1743𝑖
0 0 0 −2.7523 + 9.1743𝑖 2.7523 + 13.1743𝑖
En la tabla de datos de la red, se tienen valores que corresponden a los
transformadores regulantes en ese nodo, aunque no todos tiene este valor en la tabla.
En esta parte, con el transformador regular presente, seguiría un procedimiento para
seguir con el problema, pero en este caso, como todos los valores que se tienen es de
1.0, son unitarios, y los transformadores de este tipo no intervienen con los datos de la
red, por lo que no se le hace ninguna modificación a los datos de la red, según la sección
9.6 llamada Transformadores regulantes de Análisis de Sistemas de Potencia del
Grainger, Stevenson
La matriz de impedancias nodales por el método de método directo. En esta ocasión,
se utiliza:
𝑍𝑏𝑢𝑠 = 𝑌𝑏𝑢𝑠 −1
Figura 4. Matriz de incidencia
pág. 11
Las potencias demandadas son la siguientes (Sd) :
0.4000 + 0.1500𝑖
0 + 0𝑖
0.5400 + 0.2000𝑖
0.4000 + 0.3500𝑖
0.6500 + 0.4500𝑖
0.3000 + 0.1500𝑖
0.6000 + 0.2500𝑖
0.1000 + 0.0800𝑖
0.2500 + 0.1500𝑖
0.3000 + 0.2000𝑖
Las potencias generadas son las siguientes (Sg) :
0 + 0𝑖
0 + 0𝑖
0 + 0𝑖
0 + 0𝑖
0 + 0𝑖
0 + 0𝑖
0 + 0𝑖
1.100
0 + 0𝑖
0 + 0𝑖
Entonces como en ambos casos ya tenemos las potencias generadas y demandas, se
calculan las potencias especificadas, las cuales son:
𝑆𝑠 = 𝑆𝑔 − 𝑆𝑑
−0.4000 − 0.1500𝑖
0 + 0𝑖
−0.5400 − 0.2000𝑖
−0.4000 − 0.3500𝑖
−0.6500 − 0.4500𝑖
−0.3000 − 0.1500𝑖
−0.6000 − 0.2500𝑖
1.000 − 0.0800𝑖
−0.2500 − 0.1500𝑖
−0.3000 − 0.2000𝑖
Se declaran los voltajes iniciales para poder empezar a calcular el jacobiano y las
potencias calculas. Los voltajes son los siguientes:
𝑉 =
pág. 12
1.0200
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0100
1.0000
1.0000
Ya que tenemos las potencias especificadas, podemos obtener las potencias
calculadas, las cuales se obtienen a partir del voltaje:
𝑆𝑐 = 𝑉.∗ 𝑐𝑜𝑛𝑗(𝑌𝑏𝑢𝑠 ∗ 𝑉)
𝑆𝑐 =
0.0339 + 0.2290𝑖
0.0000 − 0.4700𝑖
0.0000 − 0.0900𝑖
0.0000 − 0.2800𝑖
0.0000 − 0.2600𝑖
−0.0000 − 2.0000𝑖
0.0000 − 2.0667𝑖
0.0000 + 0.0673𝑖
−0.0332 − 0.1245𝑖
0.0000 + 0.0000
En el código no viene todo el jacobiano como uno se imagina, una matriz enorme con
muchísimas derivadas parciales, si no que viene programado, en la forma en la que
pág. 13
uno ingresa los datos y ya nos da el Jacobiano del problema, hasta que cumpla el error
mínimo o hasta un numero de ciertas iteraciones que el usuario quiera.
El Jacobiano quedo:
Figura 5. Jacobiano
La potencia en el nodo 1 es de:
𝑆1 = 2.8016 − 1.4729𝑖
Después seguimos con el programa y obtenemos los flujos de potencia del sistema, los
cuales son:
−1.2173 − 0.4269𝑖
−0.9110 − 0.4659𝑖
−0.6658 − 0.1203𝑖
−0.2037 − 0.1491𝑖
−0.3482 − 0.1491𝑖
𝐹𝑝 = −0.10 − 0.050𝑖
0.3062 − 0.1154𝑖
0.2514 + 0.179𝑖
0.1904 + 0.198𝑖
−1.000 − 0.4748𝑖
−0.6207 − 0.1650𝑖
Para las pérdidas del sistema, queda de esta manera:
pág. 14
Figura 6. Perdidas en la línea
Después del jacobiano, se calculo el voltaje, el cual se detuvo en la iteración #15, al
haber cumplido con el error mínimo que se pedía. Los voltajes fueron estos:
Figura 7. Voltajes
pág. 15
CODIGO EN MATLAB
Código en el que realizaron estos cálculos:
Z = [2i; 0.04+0.15i; 0.03+0.1i; 0.05+0.22i; 0.04+0.2i; 0.09+0.25i; 0.06+0.21i; 0.25i;
0.04+0.15i; 0.06+0.3i; 0.15i; 0.03+0.1i];
Y=[(0.0000 + 0.2000i)^-1; (0.0400 + 0.1500i)^-1; (0.0300 + 0.1000i)^-1; (0.0500 +
0.2200i)^-1; (0.0400 + 0.2000i)^-1; (0.0900 + 0.2500i)^-1; (0.0600 + 0.2100i)^-1;
(0.0000 + 0.2500i)^-1; (0.0400 + 0.1500i)^-1; (0.0600 + 0.3000i)^-1; (0.0000 +
0.1500i)^-1; (0.0300 + 0.1000i)^-1];
Yshunt = [0+0i; 0+0.37i; 0+0.09i; 0+0.28i; 0+0.26i; 0+2i; 0+2i; 0+0i; 0+0i; 0+0i];
Yp = diag([Y ; Yshunt]);
A = [1 -1 0 0 0 0 0 0 0 0; 1 0 0 0 0 0 0 0 -1 0; 0 1 -1 0 0 0 0 0 0 0; 0 1 0 -1 0 0 0
0 0 0; 0 1 0 0 -1 0 0 0 0 0; 0 0 1 0 0 -1 0 0 0 0; 0 0 0 1 -1 0 0 0 0 0; 0 0 0 1 0 0
0 0 0 -1; 0 0 0 0 1 0 -1 0 0 0; 0 0 0 0 0 1 -1 0 0 0; 0 0 0 0 0 0 -1 1 0 0; 0 0 0 0 0
0 0 0 1 -1; 1 0 0 0 0 0 0 0 0 0; 0 1 0 0 0 0 0 0 0 0; 0 0 1 0 0 0 0 0 0 0; 0 0 0 1 0
0 0 0 0 0; 0 0 0 0 1 0 0 0 0 0; 0 0 0 0 0 1 0 0 0 0; 0 0 0 0 0 0 1 0 0 0; 0 0 0 0 0 0
0 1 0 0; 0 0 0 0 0 0 0 0 1 0; 0 0 0 0 0 0 0 0 0 1];
% A continuacion, se presenta la Ybus
Ybus=A.'*Yp*A;
Zbus = Ybus^-1;
%TOMAMOS LOS DATOS DE LA TABLA DE CARGA Y GENERACION DEL SISTEMAS, PARA PASAR LAS
POTENCIAS A P.U
%Recordemos que tenemos 100 MVA como potencia baase
%Potencias demandada
Sd=[0.4+0.15i;0;0.54+0.2i;0.4+0.35i;0.65+0.45i;0.3+0.15i;0.6+0.25i;0.1+0.08i;0.25+0.1
5i;0.3+0.2i];
%Potencias generadas
Sg= [0;0;0;0;0;0;0;1.1;0;0];
%Potenciaas especificadas
Ss= Sg-Sd;
%Voltajes inciales
V=[1.02; 1.0; 1.0; 1.0; 1.0; 1.0; 1.0; 1.01; 1.0; 1.0];
%Las potencias calculadas se encuentran con la ecuacion:
Sc = V.*conj(Ybus*V);
%Jacobiano y Potencias
er=1;
r=1;
W=[0; 0; 0; 0; 0; 0; 0; 0; 0];
D=[0; 0; 0; 0; 0; 0; 0; 0];
while er>0.00001
pág. 16
Sc=V(1:10,r).*conj(Ybus*V(1:10,r));
P=real(Ss(2:10,1))-real(Sc(2:10,1));
Q=[imag(Ss(2:7,1))-imag(Sc(2:7,1)); imag(Ss(9:10,1))-imag(Sc(9:10,1))];
er1=max(abs(P)-abs(W));
er2=max(abs(Q)-abs(E));
er=max(er1,er2);
J=zeros(17);
for k=1:9
for x=1:9
if k==x
J(k,x)=-imag(Sc(k+1,1))-
(abs(V(k+1,r))*abs(V(k+1,r))*imag(Ybus(k+1,x+1)));
else
J(k,x)=-
(abs(V(k+1,r))*abs(V(x+1,r))*abs(Ybus(k+1,x+1)))*(sin(angle(Ybus(k+1,x+1))-
imag(V(k+1,r))+imag(V(x+1,r))));
end
end
end
for k=10:15
for x=1:9
if k==(x+9)
J(k,x)=real(Sc(k-8,1))-(abs(V(k-8,r))*abs(V(k-8,r))*real(Ybus(k-
8,x+1)));
else
J(k,x)=-(abs(V(k-8,r)*abs(V(x+1,r))*(abs(Ybus(k-
8,x+1))))*cos(angle(Ybus(k-8,x+1))-imag(V(k-8,r))+imag(V(x+1,r))));
end
end
end
for k=16:17
for x=1:9
if k==16 && x==8
J(k,x)=-real(Sc(k-7,1))-(abs(V(k-7,r))*abs(V(k-7,r))*real(Ybus(k-
7,x+1)));
elseif k==17 && x==9
J(k,x)=-real(Sc(k-7,1))-(abs(V(k-7,r))*abs(V(k-7,r))*real(Ybus(k-
7,x+1)));
else
J(k,x)=-(abs(V(k-7,r)*abs(V(x+1,r))*(abs(Ybus(k-
7,x+1))))*cos(angle(Ybus(k-7,x+1))-imag(V(k-7,r))+imag(V(x+1,r))));
end
end
end
for k=1:9
for x=10:15
if x==(k+9)
pág. 17
J(k,x)=real(Sc(k+1,1))+(abs(V(k+1,r))*abs(V(k+1,r))*real(Ybus(k+1,x-
8)));
else
J(k,x)=abs(V(k+1,r)*abs(V(x-8,r))*(abs(Ybus(k+1,x-
8))))*(cos(angle(Ybus(k+1,x-8))-imag(V(k+1,r))+imag(V(x-8,r))));
end
end
for x=16:17
if x==(k+8)
J(k,x)=real(Sc(k+1,1))+(abs(V(k+1,r))*abs(V(k+1,r))*real(Ybus(k+1,x-
7)));
else
J(k,x)=abs(V(k+1,r)*abs(V(x-7,r))*(abs(Ybus(k+1,x-
7))))*(cos(angle(Ybus(k+1,x-7))-imag(V(k+1,r))+imag(V(x-7,r))));
end
end
end
for k=10:15
for x=10:15
if k==x
J(k,x)=imag(Sc(k-8,1))-(abs(V(k-8,r))*abs(V(k-8,r))*imag(Ybus(k-8,k-
8)));
else
J(k,x)=-(abs(V(k-8,r)*abs(V(x-8,r))*(abs(Ybus(k-8,x-
8))))*sin(angle(Ybus(k-8,x-8))-imag(V(k-8,r))+imag(V(x-8,r))));
end
end
for x=16:17
J(k,x)=-(abs(V(k-8,r)*abs(V(x-7,r))*(abs(Ybus(k-8,x-
7))))*(sin(angle(Ybus(k-8,x-7))-imag(V(k-8,r))+imag(V(x-7,r)))));
end
end
for k=16:17
for x=10:15
J(k,x)=-(abs(V(k-7,r)*abs(V(x-8,r))*(abs(Ybus(k-7,x-
8))))*(sin(angle(Ybus(k-7,x-8))-imag(V(k-7,r))+imag(V(x-8,r)))));
end
for x=16:17
if k==x
J(k,x)=imag(Sc(k-7,1))-(abs(V(k-7,r))*abs(V(k-7,r))*imag(Ybus(k-7,k-
7)));
else
J(k,x)=-(abs(V(k-7,r)*abs(V(x-7,r))*(abs(Ybus(k-7,x-
7))))*(sin(angle(Ybus(k-7,x-7))-imag(V(k-7,r))+imag(V(x-7,r)))));
end
end
end
PQ=[P;Q];
sV=(J^-1)*PQ;
V(1,r+1)=V(1,1);
si=[0; sV(1:9,1)];
for o=2:7
V(o,r+1)=V(o,r)+(sV(o+8,1));
end
pág. 18
for o=8
V(o,r+1)=V(o,r);
end
for o=9:10
V(o,r+1)=V(o,r)+(sV(o+7,1));
end
for u=2:10
V(u,r+1)=V(u,r+1)+1i*(si(u,1));
end
r=r+1;
W=P;
D=Q
if r>100
break;
end
end
S1=conj(V(1,r))*((V(1,1)*(Ybus(1,1)))+((Ybus(1,9)*V(9,r))+(Ybus(1,2)*V(2,r))))-
Ss(1,1);
I12=Ybus(1,2)*(V(1,r)-V(2,r));
I19=Ybus(1,9)*(V(1,r)-V(9,r));
I23=Ybus(2,3)*(V(2,r)-V(3,r));
I24=Ybus(2,4)*(V(2,r)-V(4,r));
I25=Ybus(2,5)*(V(2,r)-V(5,r));
I36=Ybus(3,6)*(V(3,r)-V(6,r));
I45=Ybus(4,5)*(V(4,r)-V(5,r));
I140=Ybus(4,10)*(V(4,r)-V(10,r));
I57=Ybus(5,7)*(V(5,r)-V(7,r));
I67=Ybus(6,7)*(V(6,r)-V(7,r));
I87=Ybus(8,7)*(V(8,r)-V(7,r));
I910=Ybus(9,10)*(V(9,r)-V(10,r));
fp12=(V(1,r)*conj(I12));
fp19=(V(1,r)*conj(I19));
fp23=(V(2,r)*conj(I23));
fp24=(V(2,r)*conj(I24));
fp5=(V(2,r)*conj(I25));
fp36=(V(3,r)*conj(I36));
fp45=(V(4,r)*conj(I45));
fp410=(V(4,r)*conj(I140));
fp57=(V(5,r)*conj(I57));
fp67=(V(6,r)*conj(I67));
fp87=(V(8,r)*conj(I87));
fp910=(V(9,r)*conj(I910));
fp=[fp12; fp19; fp23; fp24; fp25; fp36; fp45; fp410; fp57; fp67; fp87; fp910];
pl12=(V(1,r)*conj(I12))+(V(2,r)*conj(-I12));
pl19=(V(1,r)*conj(I19))+(V(9,r)*conj(-I19));
pl23=(V(2,r)*conj(I23))+(V(3,r)*conj(-I23));
pl24=(V(2,r)*conj(I24))+(V(4,r)*conj(-I24));
pl25=(V(2,r)*conj(I25))+(V(5,r)*conj(-I25));
pl36=(V(3,r)*conj(I36))+(V(6,r)*conj(-I36));
pl45=(V(4,r)*conj(I45))+(V(5,r)*conj(-I45));
pl410=(V(4,r)*conj(I140))+(V(10,r)*conj(-I140));
pág. 19
pl57=(V(5,r)*conj(I57))+(V(7,r)*conj(-I57));
pl67=(V(6,r)*conj(I67))+(V(7,r)*conj(-I67));
pl87=(V(8,r)*conj(I87))+(V(7,r)*conj(-I87));
pl910=(V(9,r)*conj(I910))+(V(10,r)*conj(-I910));
pl=[ pl12; pl19; pl23; pl24; pl25; pl36; pl45; pl410; pl57; pl67; pl87; pl910];
disp(V);
disp(S1);
disp(fp);
disp(pl);
SIMULACION
1. Simular el caso de análisis y comparar sus resultados utilizando el software
PowerWorld.
2. Con la finalidad de mejorar el perfil de voltaje en el sistema, se desea instalar un
banco de
capacitores de 30 MVAR a 230 kV nominal. Simule la instalación del banco en los
nodos 2, 3, 4 5 y
6. Observe cómo cambian los voltajes en los nodos en función de la ubicación del
banco de
capacitores.
3. Graficar los perfiles de voltaje obtenidos en los casos analizados.
4. Realice un análisis del sistema en base a los resultados obtenidos tanto de forma
analítica como
con el uso del software propuesto
pág. 20
Caso de análisis del presente proyecto
Figura 8. Simulación del Sistema de Potencia
A continuación, se presentarán las demás características del sistema simulado.
La Ybus que se obtuvo fue la siguiente:
Figura 9. Ybus del sistema simulado
pág. 21
El Jacobiano del sistema, nos dio de la siguiente manera:
Figura 10. Jacobiano del Sistema de Potencia simulado Parte 1
Figura 11. Jacobiano del Sistema de Potencia simulado Parte 2
Las características de los dos generadores que había en el sistema, fue de:
Figura 12. Características de los generadores del sistema Parte 1
Figura 13. Características de los generadores del sistema Parte 2
pág. 22
Los voltajes obtenidos y en sí, características generales del sistema:
Figura 14. Características del sistema
Ya que tenemos el sistema y conocimos sus características, podemos comenzar con los
capacitores, que se proponen con la finalidad de mejorar el perfil de voltaje en el sistema,
se desea instalar un banco de capacitores de 30 MVAR a 230 kV nominal. Simule la
instalación del banco en los nodos 2, 3, 4 5 y 6. Observe cómo cambian los voltajes en
los nodos en función de la ubicación del banco de capacitores.
Banco de capacitores en el nodo 2:
Figura 15. Banco de capacitores en el nodo 2
Banco de capacitores en el nodo 3:
Figura 16. Banco de capacitores en el nodo 3
pág. 23
Banco de capacitores en el nodo 4:
Figura 17. Banco de capacitores en el nodo 4
Banco de capacitores en el nodo 5:
Figura 18. Banco de capacitores en el nodo 5
Banco de capacitores en el nodo 6:
Figura 19. Banco de capacitores en el nodo 6
pág. 24
El reporte nos pide que grafiquemos los resultados de cada vez que se ubicó el banco
de capacitores, para analizar en donde es la mejor ubicación. Graficamos el estado
inicial, o sea sin banco de capacitores.
Grafica 1. Voltaje respecto a nodos iniciales
Voltaje (p.u)
Nodo
Empezamos con este grafica que tiene una forma muy similar entre todos los nodos.
Se espera que haya incrementos de voltaje en el nodo donde se colocara el banco de
capacitores y que cambien un poco también el voltaje de los demás nodos.
Grafica 2. Banco de capacitores Nodo 2 Grafica 3. Banco de capacitores Nodo 3
Voltaje (p.u)
Voltaje (p.u)
Nodo Nodo
pág. 25
Grafica 4. Banco de capacitores Nodo 4 Grafica 5. Banco de capacitores Nodo 5
Voltaje (p.u)
Voltaje (p.u)
Nodo Grafica 6. Nodo 6 Nodo
Banco de capacitores en el nodo 6
Nodo
Observando las graficas y sus cambios, nos damos cuenta de que el nodo 10 siempre
es el de menor voltaje siempre, los demás siempre están sobre valores mas constantes
casi siempre, mientras el 10 siempre es el menor, así que se puede utilizar este nodo
como referencia para tomar la decisión de en donde colocar el banco de capacitores, y
podemos ver que la mejor opción es instalar este banco en el nodo 4, ya que no solo
aumenta el voltaje en el nodo 10, si no es el que eleva el voltaje en mas nodos de las
demás ubicaciones para la instalación de dicho banco de capacitores.
pág. 26
ANALISIS DE RESULTADOS DE FORMA ANALITICA Y SIMULADA
Cabe recalcar que primero se realizó la simulación y después la solución analítica, por
lo que nos basamos en la simulación para considerar como las respuestas correctas
los valores del software.
• Podemos decir que la parte inicial, la Ybus si quedo bastante parecida, de
hecho, hay pocos valores que cambian por una décima, así que parece que con
la Ybus se cumplió.
• Con la parte del Jacobiano, la parte analítica con respecto a la simulada, hay
unos valores que están muy cercanos, pero no coinciden exactamente y hay
otros valores que si están muy lejos de coincidir.
• Para el Voltaje podemos observar que en la gran mayoría de nodos no se
acerca, quizás por eso se necesita o se propone el banco de capacitores,
porque en la solución analítica también lo pide, pero por cuestiones de tiempo y
falta de experiencia para programar este tipo de cosas, no se realizó, por lo que
los voltajes finales difieren aún más.
Falto mucha más exactitud en general, sabemos que obtener la Ybus es sencillo, pero
falto analizar más o igual incluso practicar un poco, ya que mi proceso de programación
fue bastante tardado, además de que la fecha ya se aplazo bastante.
pág. 27
CONCLUSION
El ejercicio demuestra de forma muy buena lo que es un sistema de potencia, la solución
analítica y la solución por simulador se complementan bastante bien, es una pequeña
demostración de como fluye la potencia a través de las líneas, como se controla y en si
todo lo que tiene que ver con un sistema de potencia, ya que es bastante complejo y hay
muchas opciones que no se han utilizado, al menos personalmente. Así que el simulador
es de gran ayuda, pues no da un panorama de como es esto en realidad, y la parte
analítica, nos sirve mucho para entender como se soluciona el problema de flujos de
potencia, prácticamente nos obliga el mismo problema a entrar de lleno a las ecuaciones
de potencia y a su solución. En resumen, es bastante interesante todo este tema del
análisis de flujos de potencia, como se compensan las perdidas, la complejidad de
coordinar la generación del suministrador y la carga que se va a consumir, toda la
seguridad que hay detrás de esto y, sobre todo, la importancia que tiene en nuestras
vidas a todas horas, que hay demasiada investigación en el tema, desconocía todo eso.
Y como yo como estudiante de ingeniería en Energías Renovables, es completamente
nuevo todo esto, si considero que es bastante importante este apartado y que mis
compañeros de carrera deben llevar esta materia, quizás no es necesario que sea tan
profundo como mis compañeros eléctricos, pero un parque eólico, una granja solar, la
transmisión de energía que viene de una central geotérmica y en cualquier generación
de electricidad que sea renovable puede aplicar esto, así que es importante. Y
personalmente, fue todo un desafío realizar este problema, quizás no este bien, pero
estuvo interesante.
pág. 28
Bibliografía
[1] I. Matulic, «Scielo,» RevActaNova, 2 6 2003. [En línea]. Available:
[Link] [Último
acceso: 3 6 2023].
[2] A. P. Sandoval, «Expansion de la transmision de un sistema de potencia de 57 nodos,» IPN, Estado
de Mexico, 2008.
[3] L. C. Cisneros, «[Link].,» ITM-IDE, Buenos Aires, 2012.
[4] J. Greinger y W. Stevenson, Analisis de Sistemas de Potencia. U.S: [Link]-Hill, 1996
pág. 29
