Memorial do Curso de Engenharia Elétrica

__________________________________________________________

FLUXO DE POTÊNCIA
NO POWER WORLD

1
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Fluxo de potência no Power World

Cleber J. Rodrigues1
Michel Mussoline2
Matheus Dias
Wesley Vilmar

1 Acadêmico do VII Período do Curso de Engenharia Elétrica da Faculdade Comunitária União

das Américas – UNIAMÉRICA, Foz do Iguaçu-PR. E-mail: [email protected]
2 Acadêmico do VII Período do Curso de Engenharia Elétrica da Faculdade Comunitária União

das Américas – UNIAMÉRICA, Foz do Iguaçu-PR. E-mail: [email protected]

2
Artigo produzido no Componente Curricular Metodologia Científica no II Semestre de 2018 da
Centro Universitário União das Américas – UNIAMÉRICA.

2
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Como o propósito de calcular o fluxo de potência, entre duas barras tanto manual quanto
automático.

1. Considere uma linha de transmissão cuja tensão nominal é de 220 kV e tem uma
impedância de 𝑍=5+ 𝑗40 𝛺 por fase. Considere que os módulos das tensões nas
barras extremas possam ser controlados. Para os casos seguintes calcule a potência
ativa e a potência reativa saindo de cada extremo da linha e as perdas ativa e reativa
na linha, utilizando como base 220 kV e 100MVA:

1 |𝑽𝟏 | = 220kV |𝑽𝟐 | = 220kV 𝜹𝟏𝟐 = 𝟕º

2 |𝑽𝟏 | = 220kV |𝑽𝟐 | = 230kV 𝜹𝟏𝟐 = 𝟕º

3 |𝑽𝟏 | = 220kV |𝑽𝟐 | = 220kV 𝜹𝟏𝟐 = 𝟖º

 𝑉1 = 𝑉𝑘
 𝑉2 = 𝑉𝑚

 Quais fatos podem ser percebidos em relação ao sentido em que

fluem as potências levando em consideração os módulos e ângulos das
tensões?

Na primeira linha corresponde como referência para os demais, na segunda linha a 𝑄𝑝𝑒𝑟𝑑𝑎𝑠
é mais vulnerável quando existe alteração no módulo da tensão e no caso 3 a 𝑃𝑝𝑒𝑟𝑑𝑎𝑠 ela é
mais sensível a alteração do angulo.

3
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

 O que se pode perceber em relação às variações das potências

comparando os três casos?

Nas três situações foi possui evidenciar pouca variação da potência com relação a
alteração das tensões ao ângulo modificado

𝑷𝒌𝒎=𝒈𝒌𝒎|𝑽𝒌|𝟐 − 𝒈𝒌𝒎|𝑽𝒌||𝑽𝒎|𝒄𝒐𝒔𝜽𝒌𝒎 − 𝒃𝒌𝒎|𝑽𝒌||𝑽𝒎|𝒔𝒆𝒏𝜽𝒌𝒎 𝑸𝒌𝒎

=−𝒃𝒌𝒎|𝑽𝒌|𝟐 −𝒈𝒌𝒎|𝑽𝒌||𝑽𝒎|𝒔𝒆𝒏𝜽𝒌𝒎 + 𝒃𝒌𝒎|𝑽𝒌||𝑽𝒎|𝒄𝒐𝒔𝜽𝒌𝒎

Tabela 1 - Resultado obtidos no Matlab

𝑃𝐾𝑚 𝑃𝑚𝐾 𝑃𝑝𝑒𝑟𝑑𝑎𝑠 𝑄𝐾𝑚 𝑄𝑚𝐾 𝑄𝑝𝑒𝑟𝑑𝑎𝑠

1 -1,463 0,0222 -0,0926 0,0926 0,2703 0,1777
2 1,462 -1,436 0,0263 -0,6384 0,8487 0,2103
3 1,673 -1,644 0,029 -0,0913 0,3232 0,232

Para calcular a potência ativa e a potência reativa saindo de cada extremo da linha
e as perdas ativa e reativa na linha, temos que primeiro calcular a impedância
base para obtermos a impedância em por unidade (pu) para que assim acharmos
a admitancia em por unidade (pu).

𝑉 2 (220. 103 )²
𝑍𝑔 = = = 484
𝑃 100. 106
o 𝑍𝑔= impedância base
o 𝑉= tensão base
o 𝑃= potencia base

𝑍𝑓 5 + 40𝑖 5 10
𝑍𝑘𝑚𝑝𝑢 = = = +𝑗
𝑍𝑔 484 484 121
𝑟𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑟𝑒𝑎𝑡â𝑛𝑐𝑖𝑎
𝑝𝑎𝑟𝑡𝑒 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡𝑒 𝑖𝑚𝑎𝑔𝑖𝑛á𝑟𝑖𝑎
o 𝑍𝑝𝑢= impedância por fase em pu (por unidade)
o 𝑍𝑓= impedância de fase

4
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

1 1
𝑍𝑘𝑚𝑝𝑢 = = = 1,4892 − 𝑗11,9138
𝑍𝑘𝑚𝑝𝑢 5 10
484 + 𝑗 121

𝑦𝑝𝑢= admitância por fase em pu (por unidade)

Depois de terminado esses cálculos podemos seguir em frente e calcular a

potência ativa e a potência reativa e as perdas ativa e reativa na linha.

𝑘𝑚= analise do lado 1 ao lado 2 da linha

𝑚𝑘= analise do lado 2 ao lado 1 da fase
o 𝑃𝑘𝑚 𝑜𝑢 𝑃𝑚𝑘 = potência ativa
o 𝑄𝑘𝑚 𝑜𝑢 𝑄𝑚𝑘 = resistência (parte real)
o 𝑏𝑘𝑚 ou 𝑏𝑚𝑘 = reatância (parte imaginária)
o 𝑄𝑘𝑚 ou 𝑄𝑚𝑘 = potência reativa
o 𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = perdas ativas
o 𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = perdas reativas

Situação 1:

𝑃𝑘𝑚 = 𝑔𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛𝜃𝑘𝑚

𝑃𝑘𝑚 = 1,4892|1|2 − 1,4892|1||1| cos(7°) − (−11,9138)|1||1|𝑠𝑒𝑛(7°)
𝑃𝑘𝑚 = 1,463
𝑄𝑘𝑚 = −𝑏𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛𝜃𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘
𝑄𝑘𝑚 = −(−11,9138)|1|2 − 1,4892|1||1|𝑠𝑒𝑛(7°) + (−11,9138)|1||1| cos(7°)
𝑄𝑘𝑚 = −0,0926
𝑃𝑚𝑘 = 𝑔𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚| cos(−𝜃) 𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚
𝑃𝑚𝑘 = 1,4892|1|2 − 1,4892|1||1| cos(−7°) − (−11,9138)|1||1|𝑠𝑒𝑛(−7°)
𝑃𝑚𝑘 = −1,4408
𝑄𝑚𝑘 = −𝑏𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚| cos(−𝜃) 𝑘𝑚
𝑄𝑚𝑘 = −(−11,9138)|1|2 − 1,4892|1||1|𝑠𝑒𝑛(−7°) + (−11,9138)|1||1| cos(−7°)
𝑄𝑚𝑘 = 0,2703
𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑃𝑘𝑚 + 𝑃𝑚𝑘

𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 1,463 + (−1,4408) = 0,0222

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑄𝑘𝑚 + 𝑄𝑚𝑘

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = −0,0926 + 0,2703 = 0,1777

5
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Situação 2

𝑃𝑘𝑚 = 𝑔𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛𝜃𝑘𝑚

𝑃𝑘𝑚 = 1,4892|1|2 − 1,4892|1||1,045| cos(7°) − (−11,9138)|1||1,045|𝑠𝑒𝑛(7°)
𝑃𝑘𝑚 = 1,462
𝑄𝑘𝑚 = −𝑏𝑘𝑚|𝑉𝑘|² − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛𝜃𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘𝑚
𝑄𝑘𝑚 = −(−11,9138)|1|2 − 1,4892|1||1,045|𝑠𝑒𝑛(7°) + (−11,9138)|1||1,045| cos(7°)
𝑄𝑘𝑚 = −0,6384
𝑃𝑚𝑘 = 𝑔𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚| cos(−𝜃) 𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚
𝑃𝑚𝑘 = 1,4892|1|2 − 1,4892|1||1,045| cos(−7°) − (−11,9138)|1||1,045|𝑠𝑒𝑛(−7°)
𝑃𝑚𝑘 = −1,436
𝑄𝑚𝑘 = −𝑏𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚| cos(−𝜃) 𝑘𝑚
𝑄𝑚𝑘 = −(−11,9138)|1|2 − 1,4892|1||1,045|𝑠𝑒𝑛(−7°)
+ (−11,9138)|1||1,045| cos(−7°)
𝑄𝑚𝑘 = 0,8487
𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑃𝑘𝑚 + 𝑃𝑚𝑘

𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 1,462 + (−1,436) = 0,0263

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑄𝑘𝑚 + 𝑄𝑚𝑘

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = −0,6384 + 0,8487 = 0,2103

6
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Stuação 3:

𝑃𝑘𝑚 = 𝑔𝑘𝑚|𝑉𝑘|² − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛𝜃

𝑃𝑘𝑚 = 1,4892|1|2 − 1,4892|1||1| cos(7°) − (−11,9138)|1||1|𝑠𝑒𝑛(8°)

𝑃𝑘𝑚 = 1,673

𝑄𝑘𝑚 = −𝑏𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛𝜃𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑐𝑜𝑠𝜃𝑘𝑚

𝑄𝑘𝑚 = −(−11,9138)|1|2 − 1,4892|1||1|𝑠𝑒𝑛(7°) + (−11,9138)|1||1| cos(8°)

𝑄𝑘𝑚 = −0,0913

𝑃𝑚𝑘 = 𝑔𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚| cos(−𝜃) 𝑘𝑚 − 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚

𝑃𝑚𝑘 = 1,4892|1|2 − 1,4892|1||1| cos(−7°) − (−11,9138)|1||1|𝑠𝑒𝑛(−8°)

𝑃𝑚𝑘 = −1,644

𝑄𝑚𝑘 = −𝑏𝑘𝑚|𝑉𝑘|2 − 𝑔𝑘𝑚|𝑉𝑘||𝑉𝑚|𝑠𝑒𝑛(−𝜃)𝑘𝑚 + 𝑏𝑘𝑚|𝑉𝑘𝑚||𝑉𝑚| cos(−𝜃) 𝑘𝑚

𝑄𝑚𝑘 = −(−11,9138)|1|2 − 1,4892|1||1|𝑠𝑒𝑛(−7°) + (−11,9138)|1||1| cos(−8°)

𝑄𝑚𝑘 = 0,3232

𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑃𝑘𝑚 + 𝑃𝑚𝑘

𝑃𝑝𝑒𝑟𝑑𝑎𝑠 = 1,673 + (−1,644) = 0,029

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = 𝑄𝑘𝑚 + 𝑄𝑚𝑘

𝑄𝑝𝑒𝑟𝑑𝑎𝑠 = −0,0913 + 0,3232 = 0,232

7
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Código usado no Matlab:

clear all
clc
Vb = 220000;
Pb = 100000000;
Zbase = Vb/(Pb/Vb);
Zkm = 5+i*40;
Ykm = 1/Zkm;
Zpu = Zkm/Zbase;
Ypu = 1/Zpu;
gkm = real(Ypu);
bkm = imag(Ypu);
%%situação 1
Vk1 = 220000/Vb;
Vm1 = 220000/Vb;
angulokm1 = 7;
Pkm1 = (gkm*(abs(Vk1)^2)) - (gkm*abs(Vk1)*abs(Vm1)*cosd(angulokm1)) -
(bkm*abs(Vk1)*abs(Vm1)*sind(angulokm1));
Qkm1 = (-bkm*(abs(Vk1)^2)) - (gkm*abs(Vk1)*abs(Vm1)*sind(angulokm1)) +
bkm*abs(Vk1)*abs(Vm1)*cosd(angulokm1);
angulomk1 = angulokm1*(-1);
Pmk1 = (gkm*(abs(Vm1)^2)) - (gkm*abs(Vm1)*abs(Vk1)*cosd(angulomk1)) -
(bkm*abs(Vm1)*abs(Vk1)*sind(angulomk1));
Qmk1 = (-bkm*(abs(Vm1)^2)) - (gkm*abs(Vm1)*abs(Vk1)*sind(angulomk1)) +
bkm*abs(Vm1)*abs(Vk1)*cosd(angulomk1);
Pperdas1 = Pkm1+Pmk1;
Qperdas1 = Qkm1+Qmk1;
%%situação 2
Vk2 = 220000/Vb;
Vm2 = 230000/Vb;
angulokm2 = 7;
Pkm2 = (gkm*(abs(Vk2)^2)) - (gkm*abs(Vk2)*abs(Vm2)*cosd(angulokm2)) -
(bkm*abs(Vk2)*abs(Vm2)*sind(angulokm2));
Qkm2 = (-bkm*(abs(Vk2)^2)) - (gkm*abs(Vk2)*abs(Vm2)*sind(angulokm2)) +
bkm*abs(Vk2)*abs(Vm2)*cosd(angulokm2);
angulomk2 = angulokm2*(-1);
Pmk2 = (gkm*(abs(Vm2)^2)) - (gkm*abs(Vm2)*abs(Vk2)*cosd(angulomk2)) -
(bkm*abs(Vm2)*abs(Vk2)*sind(angulomk2));
Qmk2 = (-bkm*(abs(Vm2)^2)) - (gkm*abs(Vm2)*abs(Vk2)*sind(angulomk2)) +
bkm*abs(Vm2)*abs(Vk2)*cosd(angulomk2);
Pperdas2 = Pkm2+Pmk2;
Qperdas2 = Qkm2+Qmk2;
%%situação 3
Vk3 = 220000/Vb;
Vm3 = 220000/Vb;
angulokm3 = 8;
Pkm3 = (gkm*(abs(Vk3)^2)) - (gkm*abs(Vk3)*abs(Vm3)*cosd(angulokm3)) -
(bkm*abs(Vk3)*abs(Vm3)*sind(angulokm3));
Qkm3 = (-bkm*(abs(Vk3)^2)) - (gkm*abs(Vk3)*abs(Vm3)*sind(angulokm3)) +
bkm*abs(Vk3)*abs(Vm3)*cosd(angulokm3);
angulomk3 = angulokm3*(-1);

8
 Memorial do Curso de Engenharia Elétrica
__________________________________________________________

Pmk3 = (gkm*(abs(Vm3)^2)) - (gkm*abs(Vm3)*abs(Vk3)*cosd(angulomk3)) -

(bkm*abs(Vm3)*abs(Vk3)*sind(angulomk3));
Qmk3 = (-bkm*(abs(Vm3)^2)) - (gkm*abs(Vm3)*abs(Vk3)*sind(angulomk3)) +
bkm*abs(Vm3)*abs(Vk3)*cosd(angulomk3);
Pperdas3 = Pkm3+Pmk3;
Qperdas3 = Qkm3+Qmk3;
%%resultado
disp([' Fluxo de Potência ']);
disp([' Pkm | Pmk | Pperdas | Qkm | Qmk | Qperdas |']);
disp(['Situação 1: ',num2str(Pkm1,4),' | ',num2str(Pmk1,4),' |
',num2str(Pperdas1,4),' | ',num2str(Qkm1,3),' | ',num2str(Qmk1,4),' |
',num2str(Qperdas1,4),' |']);
disp(['Situação 2: ',num2str(Pkm2,4),' | ',num2str(Pmk2,4),' |
',num2str(Pperdas2,4),' | ',num2str(Qkm2,4),' | ',num2str(Qmk2,4),' |
',num2str(Qperdas2,4),' |']);
disp(['Situação 3: ',num2str(Pkm3,4),' | ',num2str(Pmk3,4),' |
',num2str(Pperdas3,4),' | ',num2str(Qkm3,3),' | ',num2str(Qmk3,4),' |
',num2str(Qperdas3,3),' |']);
disp([' ']);