Instituto Federal de Goiás Departamento de Indústria Curso de Engenharia Elétrica Caíque Torres Santos Thalys Eduardo Ferreira Rezende Comparação do método simplificado de regulação de tensão em relação ao fluxo de potência radial de Céspedes. Jataí – GO
Comparação entre o método simplificado de regulação de tensão e o método de cespedes para fluxo de potência radial para determinar a tensão nos barramentos de um sistema elétrico de potência
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Instituto Federal de Goiás
Departamento de Indústria
Curso de Engenharia Elétrica
Caíque Torres Santos
Thalys Eduardo Ferreira Rezende
Comparação do método simplificado de regulação de tensão em relação ao fluxo de
potência radial de Céspedes.
Jataí – GO
2
Caíque Torres Santos
Thalys Eduardo Ferreira Rezende
Comparação do método simplificado de regulação de tensão em relação ao fluxo de
potência radial de Céspedes.
Trabalho apresentado para avaliação na
disciplina de Qualidade de Energia do curso de
Engenharia Elétrica do Instituto Federal de
Goiás
Prof. Luciano de Souza da Costa e Silva
IFG
Jataí – GO
2014
3
Lista de figuras
Figura 1:Linha de distribuição com oito alimentadores ................................................................ 5
Figura 2: Circuito equivalente entre subestação, linha de distribuição e carga. ............................ 6
Figura 3: Representação esquemática para demonstração do método. ......................................... 6
Figura 4: Queda de tensão real. ..................................................................................................... 6
Figura 5: Tabela coleta de dados calculados e de fácil visualização ............................................ 7
Figura 6 ......................................................................................................................................... 7
Figura 7 ......................................................................................................................................... 7
Figura 8: Equacionamento ............................................................................................................ 8
Figura 9: Equação (1) .................................................................................................................... 8
Figura 10 ....................................................................................................................................... 9
Figura 11: Perfil de tensão para demanda máxima. ................................................................... 12
Figura 12: Perfil de tensão para demanda mínima. .................................................................... 12
); else ativaacumulada(n)=0; reativaacumulada(n)=0; perdasativas(n)=0; perdasreativas(n)=0; end P(n)=P(n)+perdasativas(n)+ativaacumulada(n); Q(n)=Q(n)+perdasreativas(n)+reativaacumulada(n); end v=D(:,7); for n=2:9 eq=[1 0 (2*(D(n,3)*P(n)+D(n,4)*Q(n))-(D(D(n,2),7)^2)) 0
(P(n)^2+Q(n)^2)*(D(n,3)^2+D(n,4)^2)]; c=roots(eq); D(n,7)= max(abs(c)); end erro=max(abs(D(:,7)-v)); h=h+1; end vdmin=D(:,7); formatSpec='Docente: Luciano de Souza da Costa e Silva\nDiscentes:
Caique Torres & Thalys Rezende \nDisciplina:Qualidade de energia
elétrica\n\n'; fprintf(formatSpec,h,erro); fprintf('\n\nPerfil de tensão para demanda máxima através do método
de soma de potências de cespédes:\n\n'); for n=2:9 formatSpec='O módulo da tensão na barra %d é: %4.4f V.\n'; fprintf(formatSpec,n-1,vdmax(n)); end fprintf('\n\nPerfil de tensão para demanda mínima através do método
de soma de potências de cespédes:\n\n'); for n=2:9 formatSpec='O módulo da tensão na barra %d é: %4.4f V.\n'; fprintf(formatSpec,n-1,vdmin(n)); end
k(n)=(real(zcabo)*fp(n)+imag(zcabo)*sin(acos(fp(n))))*100/(13800^2); end for(n=8:-1:1) deltavpercentualmax(n)=k(n)*H(n,1)*S3facmax(n); deltavpercentualmin(n)=k(n)*H(n,1)*S3facmin(n); end for(n=1:8) if(n==1) deltavperacmax(n)=deltavpercentualmax(n); deltavperacmin(n)=deltavpercentualmin(n); else deltavperacmax(n)=deltavpercentualmax(n)+deltavperacmax(n-1); deltavperacmin(n)=deltavpercentualmin(n)+deltavperacmin(n-1); end end for(n=1:8) if(n==1) vbmsmax(n)=13800-(deltavpercentualmax(n))*138; vbmsmin(n)=13800-(deltavpercentualmin(n))*138; else vbmsmax(n)=vbmsmax(n-1)-(deltavpercentualmax(n))*138; vbmsmin(n)=vbmsmin(n-1)-(deltavpercentualmin(n))*138; end end fprintf('\n\nPerfil de tensão para demanda máxima através do método
simplificado:\n\n'); for n=1:8 formatSpec='O módulo da tensão na barra %d é: %4.4f V.\n'; fprintf(formatSpec,n,vbmsmax(n)); end fprintf('\n\nPerfil de tensão para demanda mínima através do método
simplificado:\n\n'); for n=1:8 formatSpec='O módulo da tensão na barra %d é: %4.4f V.\n'; fprintf(formatSpec,n,vbmsmin(n)); end vbmax(1)=13800; vbmin(1)=13800; for(n=1:8)
12
vbmax(n+1)=vbmsmax(n); vbmin(n+1)=vbmsmin(n); end
%Comparação dos gráficos figure(1); plot(distancia,vdmax,distancia,vbmax,distancia,vadli,distancia,vadls,d