Centro Federal de Educação Tecnológica de Minas

Gerais
Departamento de Engenharia Elétrica
Engenharia Elétrica

E STIMAÇÃO DE P ARÂM ETROS DE L INHAS DE

T RANSM ISSÃO – E STUDO DE C ASO

Luis Felipe Soares de Souza

Julho/2015
Centro Federal de Educação Tecnológica de Minas Gerais
Departamento de Engenharia Elétrica
Amazonas 7675 – Nova Gameleira – Belo Horizonte, Minas Gerais – 30510-
000
(31) – 3319-6834. [email protected]

Luis Felipe Soares de Souza

E STIMAÇÃO DE P ARÂM ETROS DE L INHAS DE

T RANSM ISSÃO – E STUDO DE C ASO

Trabalho de Conclusão de Curso submetido

à banca examinadora designada pelo
Colegiado do Departamento de Engenharia
Elétrica do Centro Federal de Educação
Tecnológica, como parte dos requisitos
necessários à obtenção do grau de bacharel
em Engenharia Elétrica.

Área de Concentração: Sistemas Elétricos de

Potência
Orientador (a): Patrícia Romeiro Silva Jota
CEFET-MG

Belo Horizonte
Departamento de Engenharia Elétrica
Julho/2015
Folha de Aprovação a ser anexada
 Agradecimentos

Agradeço a Deus acima de todas as coisas;

Aos meus pais, Luiz e Denise, meu braço forte;
Irmãos, queridos e companheiros;
Ana, minha inspiração;
E amigos de todas as horas.
À Orientadora Patrícia e demais professores do CEFET meus sinceros
agradecimentos.
Aos colegas da CEMIG: obrigado pela disponibilidade.
 Resumo

As linhas de transmissão, parte integrante dos sistemas elétricos de potência,

e responsáveis pelo transporte de energia das usinas geradoras às unidades
consumidoras podem ser modeladas matematicamente, de forma que os parâmetros
deste modelo podem ser estimados por meio de métodos os quais serão
apresentados neste trabalho.
O modelo utilizado, conhecido como modelo 𝜋, apresenta quatro parâmetros:
resistência e reatância em série, susceptância e condutância em paralelo, sendo que
este apresenta valores desprezíveis.
Problemas usuais nos estudos de linhas de transmissão assumem que estes
parâmetros já são conhecidos, bem como os dados de entrada ou saída, sendo
possível assim calcular as variáveis de saída e entrada respectivamente. Este trabalho
apresenta um problema inverso no qual as medidas de entrada e saída são
conhecidas e desta forma, pode-se estimar os parâmetros da linha (INDULKAR;
RAMALINGAM,2007).
As medidas de entrada e saída citadas tratam-se de tensões, correntes, e
portanto, potências nos terminais das linhas. Estas medidas estão disponíveis nas
unidades de monitoramento fasorial (PMUs), os quais processam os dados que são
adquiridos por medidores digitais nas subestações. Na ausência de dados dos PMUs,
os Sistemas Supervisórios de Aquisição de Dados (SCADA) fornecem as medidas
necessárias.
Este estudo de caso utiliza as medidas provenientes do Sistema Supervisório
da Companhia Energética de Minas Gerais, CEMIG, da linha de transmissão
(Mesquita – Baguari), com tensão nominal de 230kV e extensão de 67,51 km. Foram
utilizados dados colhidos nos dias 29 e 30 de junho de 2015 permitindo a estimação
dos parâmetros da linha além de estudos da variação destes ao longo do dia.
Um dos métodos utilizados permite além da obtenção dos três parâmetros
(resistência e reatância em série, e susceptância shunt), a obtenção do ângulo entre
as tensões dos dois terminais da linha. O ângulo entre essas tensões possibilita
estudos do fluxo de potência ativa em uma linha, além da verificação de que esta
opera em regiões de estabilidade.
i
 Para que esses parâmetros sejam estimados, torna-se necessário a
manipulação de histórico de dados de medidas elétricas: tensões, correntes e
potência nos terminais da linha. Manipulação que se dá por meio de métodos
matemáticos como o método dos mínimos quadrados, e o método de Newton-
Raphson.
Este trabalho então, consiste em um estudo de caso, no qual se estima os
parâmetros de linhas de transmissão através do histórico de dados do Sistema de
Medição Fasorial da Companhia Energética de Minas Gerais (CEMIG).

ii
 Abstract

The transmission lines, one of the components of the Electric Power System,
and responsible for moving energy from the power station to the customers of
electricity, can be mathematically modeled in such way that the parameters of this
model can be estimated.
The 𝜋 model, has four parameters: series resistance and reactance, and shunt
susceptance and conductance which is a negligible value.
Usual problems of transmission line studies, consider that these values area
already known, as well as the input and output data, making possible to calculate
output and input respectively. This paper presents an inverse problem in which
measurements from both ends of the lines are known so it is possible to calculate the
transmission line parameters. (INDULKAR; RAMALINGAM, 2007).
Input and output measurements are voltages, currents and power on the ends
of the lines. These measurements are available on Phasor Measurement Units
(PMUs), which process the data measured on the substations. In the absence of
PMUs, Supervisory Control and Data Acquisition (SCADA) are responsible for the
measurements.
This case study uses measurements from the Phasor Measurement System
from CEMIG, of the line (Mesquita – Baguari), with nominal voltage of 230 kV and
67,51 km long distance. The data from June 29 and 30 of 2015 allows the estimation
of parameters and studies about it.
One of the methods used allows also the estimation of the angle between the
voltage phasors, which tells the active power flux.
The estimation is possible by manipulating historical data of electric
measurements through math methods, such as Least Square Method and Newton-
Raphson Method.
So this paper is a case study in which the parameters of transmission lines from
CEMIG are estimated.

iii
 Sumário

Capítulo 1 ......................................................................................................................................................... 10
Introdução .............................................................................................................................................................................. 10

1.1. Relevância do tema em questão ...................................................................................................... 10

1.2. Objetivos do trabalho ........................................................................................................................ 11

1.3. Metodologia............................................................................................................................................ 12

1.4. Organização do trabalho .................................................................................................................. 12

Capítulo 2 ......................................................................................................................................................... 13
Métodos Numéricos aplicados na estimação de parâmetros de linhas de transmissão .................... 13

2.1. Estimação de Parâmetros em Sistemas Elétricos de Potência .............................................. 13

2.2. Métodos dos mínimos quadrados lineares e o processo de estimação ............................. 13

2.3. Método iterativo para resolução de sistemas não lineares ................................................... 16

2.4. Método de Newton-Raphson .......................................................................................................... 16

Capítulo 3 ......................................................................................................................................................... 21
Linhas de transmissão: modelos, parâmetros e fluxo de potência ativa ................................................... 21

3.1. Primeiro estudo – módulos de tensão, corrente e potência disponíveis .......................... 22

3.2. Segundo estudo – medidas fasoriais de corrente e tensão disponíveis ............................ 27

Capítulo 4 ......................................................................................................................................................... 30
Medidas disponíveis – Sistemas de medição em Sistemas Elétricos ............................................................ 30

4.1. Sistema de Medição Fasorial.......................................................................................................... 30

4.2. Utilização do SISMEF no Brasil ......................................................................................................... 32

4.3. Localização de Faltas do tipo Curto-Circuito e a Estimação de Parâmetros .................... 34

4.4. Utilização de Unidades Terminais Remotas na aquisição de dados do Sistema Elétrico

.............................................................................................................................................................................. 36

4.5. Conclusão acerca do uso das medidas para estimação de parâmetros ............................. 37

Capítulo 5 ......................................................................................................................................................... 38
Estudo de Caso – Linha de Transmissão da CEMIG: Mesquita - Baguari .................................................. 38

iv
5.1. Newton-Raphson e a estimação do ângulo da corrente 𝑰𝒓 ..................................................... 40

5.2. Método dos mínimos quadrados e análise dos parâmetros .................................................. 44

Capítulo 6 ......................................................................................................................................................... 46
Análise dos Resultados – Variações dos Parâmetros da linha ....................................................................... 46

6.1. Resistência em Série.......................................................................................................................... 46

6.2. Reatância em série .............................................................................................................................. 50

6.3. Susceptância Shunt............................................................................................................................ 53

6.4. Métodos numéricos ............................................................................................................................ 54

6.5. Sincronização dos PMUs ................................................................................................................. 54

Capítulo 7 ......................................................................................................................................................... 56
Conclusão ................................................................................................................................................................................ 56

v
 Lista de Figuras
Figura 1. Ajuste realizado por uma reta ................................................................................................................................... 14

Figura 2. Método de Newton Raphson ..................................................................................................................................... 17

Figura 3. Aproximação - Newton Raphson............................................................................................................................. 18

Figura 4. Três iterações para encontro da solução ........................................................................................................... 19

Figura 5. Circuito π - Modelo das linhas de transmissão ............................................................................................... 21

Figura 6. Diagrama fasorial ............................................................................................................................................................. 23

Figura 7. Fluxo de potência ativa e o ângulo 𝛿 .................................................................................................................... 27

Figura 8. Elementos constituintes de um PMU .................................................................................................................... 30

Figura 9. Parâmetros em tempo real ......................................................................................................................................... 33

Figura 10. Tela de monitoração em tempo real do Projeto MedFasee ................................................................... 34

Figura 11. Fluxograma do Método Proposto por Silveira ............................................................................................... 35

Figura 12. Sistema de Gerenciamento de Energia ............................................................................................................ 36

Figura 13. Linha de transmissão: MESQUITA - BAGUARI ........................................................................................... 39

Figura 14. Barra – Baguari .............................................................................................................................................................. 41

Figura 15. Barra – Mesquita ........................................................................................................................................................... 41

Figura 16. Seção DLIN - Arquivo ONS ..................................................................................................................................... 41

Figura 17. Variação da Resistência ao longo do dia ......................................................................................................... 44

Figura 18. Variação da Reatância ao longo do dia ............................................................................................................ 45

Figura 19. Variação da Susceptância ao longo do dia..................................................................................................... 45

Figura 20. Variação da resistência ............................................................................................................................................. 47

Figura 21. Dispersão da resistência em função da corrente ........................................................................................ 48

Figura 22. Histograma – Resistência......................................................................................................................................... 48

Figura 23. Variação da Corrente .................................................................................................................................................. 49

Figura 24. Variação da resistência ............................................................................................................................................. 49

Figura 25. Gráfico de dispersão da resistência .................................................................................................................... 50

Figura 26. Condutores na linha de transmissão .................................................................................................................. 50

Figura 27. Variação da reatância ................................................................................................................................................. 51

vi
Figura 28. Gráfico de dispersão da reatância ....................................................................................................................... 52

Figura 29. Histograma - reatância ............................................................................................................................................... 52

Figura 30. Variação da susceptância ........................................................................................................................................ 53

Figura 31. Mudanças de ângulo no PMU .............................................................................................................................. 55

vii
 Lista de Tabelas

Tabela 1. Processo iterativo - Newton Raphson ...................................................................................................................... 20

Tabela 2. Medidas - Linha de Transmissão 1 ............................................................................................................................. 25

Tabela 3. Parâmetros obtidos com o software MATLAB ...................................................................................................... 26

Tabela 4. Medidas - Linha Bauru – Cabreúva ............................................................................................................................ 29

Tabela 5. Parâmetros encontrados - Linha Bauru - Cabreúva ........................................................................................... 29

Tabela 6. Parâmetros – ONS......................................................................................................................................................... 39

Tabela 7. Parâmetros – CEMIG ................................................................................................................................................... 40

Tabela 8. Parâmetros obtidos no método de Newton Raphson ........................................................................................ 43

Tabela 9. Cinco primeiras medições realizadas pelo PMU .................................................................................................. 43

Tabela 10. Medição após 1 minuto ............................................................................................................................................. 44

viii
 Lista de Abreviações

CEMIG - Companhia Energética de Mians Gerais

EESEP - Estimação de Estados em Sistemas Elétricos de Potência
GPS - Sistema de Posicionamento Global
PDC - Phasor Data Concentrator
PPS - Sistema de Posicionamento Preciso
PMUs - Unidade de Medição Fasorial
RTU - Unidade Terminal Remota
SCADA - Sistema de Supervisão e Aquisição de Dados
SISMEF - Sistema de Medição Fasorial
ONS - Operador Nacional do Sistema

ix
 Capítulo 1

Introdução

Este capítulo tem como objetivo apresentar a relevância da estimação de

parâmetros nos Sistemas Elétricos de Potência, os objetivos deste estudo de caso
bem como a metodologia utilizada para realização deste.

1.1. Relevância do tema em questão

As operações dos sistemas elétricos são extremamente importantes pelo nível

de instantaneidade requerido. Nos centros de operação e distribuição das companhias
energéticas, operadores acompanham as grandezas elétricas em tempo real e de
acordo com seus comportamentos agem sob o sistema através de telecomandos.
A modelagem da rede elétrica em tempo real, tem como o Estimador de Estado
sua ferramenta computacional mais importante. Um Estimador de Estado em
Sistemas Elétricos de Potência (EESEP), processa informações redundantes
(medidas analógicas e lógicas), e através destas fornece as tensões complexas do
sistema. Medidas analógicas referem-se a medidas de fluxo de potência e as lógicas
correspondem à estados de chaves e disjuntores. Elas são provenientes do Sistema
de Supervisão e Aquisição de Dados (SCADA). Através do histórico de medidas deste
pode - se chegar aos parâmetros de linhas de transmissão.
Em uma Estimação de Estado supõe-se que os valores do modelo 𝜋 de linhas
de transmissão não apresentam incerteza, porém, o que se sabe, é que essa
suposição não é verdadeira. Estes parâmetros podem estar incorretos por causa de
estimativas grosseiras, acarretando na degradação de resultados da EESEP. Por
outro lado, a utilização de parâmetros confiáveis pode aumentar a margem de
transferência de potência e assim, trazer benefícios econômicos ao sistema. Ou seja,
os parâmetros de elementos da rede como as linhas de transmissão podem tornar o
processo de Estimação de Estado mais eficaz. A indutância da linha, por exemplo,
introduz os maiores erros sobre os valores estimados.

10
 Erros em parâmetros podem não só influenciar na Estimação de Estado como
na análise de contingências, curto-circuito, planejamento de transmissão, etc.
(ALBERTINI, 2010).
Este trabalho consiste na apresentação de métodos matemáticos para
estimação destes parâmetros e pelo fato destes contribuírem no processo de
Estimação de Estado é que este projeto se torna relevante. Por meio de um
considerável número de medidas foi possível a análise da influência das condições de
tempo nos parâmetros da linha e como resultado das equações utilizadas pode-se
obter também o ângulo entre as tensões dos terminais do modelo de linha utilizado.
Ou seja, foi possível analisar o fluxo de potência ativa e verificação de operação
estável da linha.

1.2. Objetivos do trabalho

Diante da relevância do trabalho apresentada no item anterior propõe-se neste

trabalho, principalmente, estimar os parâmetros de linhas de transmissão através de
histórico de dados reais, provenientes de Sistemas Supervisórios reais, no caso, o
Sistema Supervisório da Cemig. Diante das possibilidades de manipulação desse
histórico, dois métodos matemáticos são utilizados: um para resolução de sistemas
lineares, mínimos quadrados lineares, e outro para resolução de equações não
lineares, Newton Raphson. Assim, será possível estimar os parâmetros da linha, para
posterior comparação com os parâmetros apresentados no sistema real. As principais
análises, como já citado na seção 1.1 são da influência das condições climáticas nos
parâmetros, do fluxo de potência ativa na linha e verificação de operação desta em
níveis estáveis.
O trabalho tem, também, como objetivo, mostrar a relevância do histórico de
dados dos sistemas supervisórios das companhias energéticas. Ou seja, permitir
entender que, quando processadas, as medidas acumuladas no tempo podem ser
traduzidas em informações.

11
1.3. Metodologia

Para que fosse possível a estimação de parâmetros, os métodos

implementados em estimação em Sistemas Elétricos de Potência foram estudados.
Trabalhos foram replicados e estudados para compreensão das possibilidades de
estudo em uma linha de transmissão. Após estudo de conceitos de estimação,
implementação de métodos numéricos e trabalhos já realizados, buscou-se os dados
reais da Companhia Energética de Minas Gerais (CEMIG), afim de estudá-los e assim
implementar os métodos necessários. Com medidas colhidas em intervalos a cada
segundo nos dias 29 e 30 de Junho de 2015 e através dos trabalhos já realizados na
área, entendeu-se que são várias as possibilidades de estimação nos estudos de linha
de transmissão, inclusive estudos do fluxo de potência ativa.

1.4. Organização do trabalho

Primeiramente, apresenta-se a relevância deste trabalho. Conceitos de

estimação e um breve histórico desta nos Sistemas Elétricos de Potência são
apresentados no capítulo 2. Os métodos matemáticos utilizados no trabalho são
apresentados no capítulo 3. O Capítulo 4 mostra o modelo 𝜋 de linha de transmissão
e suas equações. O próximo passo, neste trabalho, também no capítulo 4, foi a
escolha de um método apropriado para utilizar os dados colhidos das linhas da CEMIG
e através do software MATLAB os métodos foram utilizados. Os parâmetros obtidos
são apresentados e posteriormente analisados no capítulo 5.

12
 Capítulo 2

Métodos Numéricos aplicados na estimação de

parâmetros de linhas de transmissão

2.1. Estimação de Parâmetros em Sistemas Elétricos de Potência

Um Sistema Elétrico de Potência se divide, basicamente, em três processos:

Geração, Transmissão e Distribuição de Energia. Dentro destes processos, os
principais elementos podem ser citados:
 Linha de Transmissão;
 Transformador de Potência;
 Gerador;
 Carga
Para cada um desses elementos têm-se um modelo, e cada modelo possui
parâmetros a serem estimados. O processo de estimação envolve a utilização de
métodos numéricos, seja para o tratamento estatístico de uma série de parâmetros
calculados, ou até mesmo para o cálculo destes diante de equações não lineares de
um modelo. Diante das medidas disponíveis, o cálculo dos parâmetros pode ser direto
ou pode envolver métodos mais rebuscados de resolução de equações. Este trabalho
apresenta dois métodos para estimação de parâmetros de linhas de transmissão:
método dos mínimos quadrados, para resolução de sistemas lineares e o método de
Newton-Raphson, para resolução de sistemas não lineares.

2.2. Métodos dos mínimos quadrados lineares e o processo de

estimação

Um processo de estimação de parâmetros se baseia em conhecimentos sobre

a teoria de amostragem e sobre inferência estatística. Considerando uma distribuição
de pares de medição (𝑥 𝑒 𝑓(𝑥)) a estimação de parâmetros consistiria em encontrar
os parâmetros da melhor função que se ajusta a esta distribuição. Considerando-se

13
os dados experimentais do caso da figura 1 percebe-se claramente que o melhor
ajuste seria realizado por uma reta.

Figura 1. Ajuste realizado por uma reta

Fonte: Elaborada pelo autor

Ao observar o gráfico da figura 1, claramente se observa que a reta teórica,

função de ajuste, não passa por todos os pontos de medição. Alguns pontos se
encontram com ordenadas superiores a reta teórica e outros com ordenadas
inferiores. A essa diferença entre ordenada das medições e ordenada do ajuste se dá
o nome de resíduo. O método dos mínimos quadrados minimiza a soma dos
quadrados dos resíduos afim de determinar parâmetros desconhecidos (FILHO,
2013).
Seja o resíduo determinado por ∆𝑦 = 𝑦𝑚𝑒𝑑 − 𝑓(𝑥𝑖 ) , onde 𝑦𝑚𝑒𝑑 equivale a
ordenada de medição e 𝑓(𝑥𝑖 ) a função de ajuste, a soma de todos os resíduos é dada
pela equação 2.1:

𝑁
𝑄= ∑ 𝑦𝑚𝑒𝑑− 𝑓(𝑥𝑖 )
𝑖=1

(2.1)
Calculando-se a soma dos quadrados encontra-se a equação 2.1 que de
acordo com o método dos mínimos quadrados lineares deve ser minimizada.

14
 𝑁 2
2
𝑄 = ∑ (𝑦𝑚𝑒𝑑 − 𝑓(𝑥𝑖 ))
𝑖=1

(2.2)
A função de ajuste 𝑓(𝑥) utilizada no caso da figura 1 foi a função de uma reta,
dada por: 𝑓(𝑥) = 𝑎𝑥 + 𝑏. Cada distribuição de dados permite o ajuste por diferentes
funções, podendo estas serem uma parábola, função cúbica, exponencial, etc.
Substituindo-se o ajuste por reta na equação 2.2, obtêm-se a equação abaixo:

𝑁
𝑄 2 (𝑎, 𝑏) = ∑ (𝑦𝑚𝑒𝑑 − 𝑎 − 𝑏𝑥𝑖 )2
𝑖=1

(2.3)
Para um conjunto de medidas tem que:

𝑥1
∆𝑦1 𝑦𝑚𝑒𝑑1 1 𝑥2
𝑎 ∆𝑦 𝑦 1 𝑥3
𝐴 = ( ) , ∆𝑦 = ( 2 ) , 𝑦𝑚𝑒𝑑 = ( 𝑚𝑒𝑑2 ) 𝑒 𝑋 =
𝑏 ⋮ ⋮ ⋮ ⋮
∆𝑦3 𝑦𝑚𝑒𝑑3 1
( 𝑥3 )

Apesar da função de ajuste não ser linear para 𝑥 observa-se que ela é linear
para os parâmetros, e desta forma:

𝑋 × 𝐴 = 𝑌𝑚𝑒𝑑
(2.4)
Aplicando-se a matriz transposta de 𝑋 nos dois lados da equação:

𝑋 𝑇 × 𝑋 × 𝐴 = 𝑋 𝑇 × 𝑌𝑚𝑒𝑑
(2.5)
Aplicando-se agora a inversa de (𝑋 𝑇 × 𝑋) dos dois lados da equação:

(𝑋 𝑇 × 𝑋)−1 × 𝑋 𝑇 × 𝑋 × 𝐴 = (𝑋 𝑇 × 𝑋)−1 𝑋 𝑇 × 𝑌𝑚𝑒𝑑

(2.6)
Como a matriz inversa de uma matriz multiplicada por ela mesmo é igual a
matriz identidade 𝐼 tem-se que:

15
 𝐼 × 𝐴 = (𝑋 𝑇 × 𝑋)−1 𝑋 𝑇 × 𝑌𝑚𝑒𝑑
(2.7)
e portanto, a matriz dos coeficientes 𝐴 é dada por:
𝐴 = (𝑋 𝑇 × 𝑋)−1 𝑋 𝑇 × 𝑌𝑚𝑒𝑑
(2.8)
Lembrando que o ajuste pode ser dado por qualquer tipo de função e
dependerá da distribuição de dados em questão.
Nesta visão, a estimação está relacionada a obtenção dos parâmetros diante
de várias amostras, as quais são as diversas entradas e saídas de um processo.
Alguns modelos proporcionam a existência de sistemas não lineares e dessa forma,
a estimação consiste na utilização de métodos iterativos para resolução de equações
não lineares. O método apresentado na seção 2.3, de Newton-Raphshon é um
exemplo de procedimento iterativo. Obviamente, os parâmetros obtidos para uma
amostra podem ser reunidos às outras amostras e assim estes podem ser tratados
também, estatisticamente.

2.3. Método iterativo para resolução de sistemas não lineares

Métodos iterativos tem como objetivo encontrar a raiz 𝑥 aproximada de uma

função 𝑓(𝑥). Basicamente, o procedimento destes métodos se inicia considerando-
se uma raiz inicial e a partir de iterações e com uma condição de convergência
encontra-se a raiz aproximada da função. A cada iteração têm-se um novo valor de
𝑥 e um novo valor de 𝑓(𝑥) , até que este novo valor possua erro dentro da tolerância
pré-estabelecida (VIANA, 2001). Um exemplo de método iterativo para resolução de
sistemas não lineares é o método de Newton Raphson, explicado na seção 2.4.

2.4. Método de Newton-Raphson

Considerado por muitos o melhor método para encontro das raízes de uma
função real este método é atribuído a Isaac Newton e Joseph Raphson (AZEVEDO,
2013). Como citado na seção 2.3 o objetivo é encontrar valores de 𝑥 para o qual a
função 𝑓(𝑥) se anula. Ou seja, na figura 2, a solução, caso a tolerância estabelecida
16
esteja em um intervalo no qual 𝑓(𝑥𝑛+2 ) está inserido, será 𝑥𝑛+2, após realizar-se
iterações que passassem por 𝑥𝑛 e 𝑥𝑛+1 .

Figura 2. Método de Newton Raphson

Fonte: http://fourier.eng.hmc.edu/e176/lectures/NM/node20.html(2015)

Considerando-se um ponto 𝑥0 suficientemente próximo de 𝑥𝑛+2, pode-se

estimar a distância entre esses pontos por meio da expansão da função 𝑓(𝑥),
conhecida como expansão em série de Taylor (em torno de 𝑥0 ) (HAFFNER, 2007).

1 𝜕𝑓(𝑥0 ) 1 𝜕 2 𝑓(𝑥0 )
𝑓(𝑥) = 𝑓(𝑥0 ) + (𝑥 − 𝑥0 ) + (𝑥 − 𝑥0 )2 + ⋯
1! 𝜕𝑥 2! 𝜕𝑥 2
(2.9)
Desprezando-se os termos de ordem superior a 2 tem-se:

1 𝜕𝑓(𝑥0 )
𝑓(𝑥) = 𝑓(𝑥0 ) + (𝑥 − 𝑥0 )
1! 𝜕𝑥
(2.10)
E sendo Δ𝑋 = 𝑥 − 𝑥0 a equação 3.17 fica da seguinte forma:

1 𝜕𝑓(𝑥0 )
𝑓(𝑥0 + Δ𝑋 ) = 𝑓(𝑥0 ) + . ( Δ𝑋) = 𝑓(𝑥0 ) + 𝑓 ′ (𝑥0 ). Δ𝑋
1! 𝜕𝑥
(2.11)

17
 Como Δ𝑋 = 𝑥 − 𝑥0 , tem-se que:

𝑓(𝑥0 + Δ𝑋 ) ≈ 𝑓(𝑥 ) = 0
(2.12)
Sendo assim:
𝑓(𝑥0 ) + 𝑓 ′ (𝑥0 ). Δ𝑋 ≈ 0
(2.13)
𝑓(𝑥0 )
Δ𝑥 ≈ −
𝑓 ′ (𝑥0 )
(2.14)
Observe que o resultado é aproximado e portanto 𝑥0 + Δ𝑋 é na verdade 𝑥1 , e
não a solução 𝑥𝑠 conforme figura 3:

Figura 3. Aproximação - Newton Raphson

Fonte: Notas de aula Professora Patrícia Romeiro da Silva Jota

𝑥1 só será a solução procurada se o valor da tolerância estabelecida estiver

antes de 𝑓(𝑥1 ) , imagem referida a 𝑥1 .
Os passos citados abaixo descrevem o algoritmo para encontro da raiz da
função 𝑓(𝑥). Logo após os passos um exemplo é apresentado (JOTA, 2015).

1) O contador de iterações 𝑣 é inicializado ->𝑣 = 0;

2) Escolhe-se um ponto inicial 𝑥 = 𝑥 𝑣 = 𝑥 0 ;

18
 3) Calcula-se o valor da função 𝑓(𝑥) no ponto 𝑥 𝑣 , obtendo-se 𝑓(𝑥 𝑣 );
4) Se o valor de 𝑓(𝑥 𝑣 ) for menor que a tolerância 𝜀 estabelecida: 𝑓(𝑥 𝑣 )<𝜀,
então 𝑥 𝑣 será a solução procurada. Caso contrário, deve-se continuar o
processo;
5) Calcula-se a derivada de 𝑓(𝑥 𝑣 ) e então obtêm-se o valor de Δ𝑥 a partir da
equação 3.21;
6) O novo ponto para nova iteração 𝑣 = 1 é calculado: 𝑥 𝑣+1 = 𝑥 𝑣 + ∆𝑥 𝑣 e
então o procedimento é repetido até que se encontre 𝑓(𝑥 𝑣 ) com valor abaixo
da iteração 𝜀 pré-estabelecida.
A figura 4 mostra três iterações para o encontro da solução. Repare que
𝑥 (3) não se equivale a 𝑥𝑠 porém já se encontra dentro do valor de tolerância 𝜀.

Figura 4. Três iterações para encontro da solução

Fonte: Jota(2015) – Notas de aula

O exemplo a seguir demonstra a aplicação do método (HAFFNER, 2007).

𝑓1 → 2𝑥 + 𝑦 = 2𝑥1 + 𝑥2 = 4
(2.15)
𝑓2 → 2𝑥 + 𝑦 2 = 2𝑥1 + 𝑥2 2 = 6
(2.16)

19
 Repare que, para efeitos didáticos, a variável 𝑦 foi substituída por 𝑥2 .
Diferentemente do exemplo literal, o problema agora é n-dimensional e a derivada
agora é feita por meio da matriz Jacobiana 𝐽 formada pelas derivadas parciais das
funções 𝑓1 e 𝑓2 em relação a 𝑥1 e 𝑥2 . Desta forma, para o exemplo em questão:
𝜕𝑓1 (𝑥) 𝜕𝑓1 (𝑥)
𝜕𝑓(𝑥) 𝜕𝑥1 𝜕𝑥2 2 1
𝐽 = = = [ ]
𝜕𝑥 𝜕𝑓2 (𝑥) 𝜕𝑓2 (𝑥) 2 2𝑥2
[ 𝜕𝑥1 𝜕𝑥2 ]
(2.17)
A Tabela 1 mostra os resultados parciais do processo iterativo considerando-
se os passos de 1 à 6 e uma tolerância 𝜀𝑥 = 𝜀𝑦 = 0,001: (HAFFNER, 2007)

𝒗 𝒙𝒗𝟏 𝒇𝟏 (𝒙𝒗𝟏 ) −[𝑱(𝒙𝒗 )]

−𝟏 ∆𝒙𝒗𝟏
𝒙𝒗𝟐 𝒇𝟐 (𝒙𝒗𝟐 ) ∆𝒙𝒗𝟐
𝟎 0 −1 −0,6 0,1 0,9
0,2 −0,2
3 3 −0,8
𝟏 0,9 0 −0,647 0,147 0,094
0,294 −0,294
2,2 0,64 −0,188
𝟐 0,994 0 −0,665 0,165 0,00586
0,331 −0,331
2,012 0,0353 −0,0117
𝟑 0,99997 0
2,000046 0,000137

Tabela 1. Processo iterativo - Newton Raphson

FONTE: http://slhaffner.phpnet.us/sistemas_de_energia_1/seI7a.pdf

Repare que em 𝑣 = 2 apesar de 𝑓1 (𝑥1𝑣 ) apresentar tolerância menor do que

estabelecida, o processo deve continuar já que a condição de parada é que o vetor
de funções formado por 𝑓1 (𝑥1𝑣 ) e 𝑓2 (𝑥2𝑣 ) tenha tolerância 𝜀 = 0,001. Condição satisfeita
em 𝑣 = 3.

20
 Capítulo 3

Linhas de transmissão: modelos, parâmetros e fluxo

de potência ativa

Parte fundamental de um sistema elétrico de potência, as linhas de

transmissão, responsável pelo transporte da energia gerada às unidades
consumidoras, são modeladas através do circuito π, representado na figura 5:

Figura 5. Circuito π - Modelo das linhas de transmissão

FONTE:INDULKAR; RAMALINGAM, 2007,p.2, com modificações

Onde 𝑅 e 𝑋 representam a resistência e a reatância em série , formando a

impedância em série 𝑍𝜋 . e 𝐺 e 𝐵 são a condutância e susceptância shunt formando a
admitância 𝑌𝜋 em derivação do modelo. Os valores de 𝐺 são desprezíveis e
geralmente são desconsiderados. Em casos de linhas curtas (< 80𝑘𝑚) 𝐵 é
desprezível e portanto sobram os parâmetros 𝑅 e 𝑋 para serem calculados. Linhas
longas requerem correções feitas em função do comprimento 𝑙 da linha e da constante
de propagação 𝛾 da linha, dada pela equação 3.1:

𝑍𝜋 .𝑌𝜋
𝑐𝑜𝑠ℎ−1 ( + 1)
2
𝛾=
𝑙
(3.1)

21
 Onde 𝑍𝜋 e 𝑌𝜋 são os valores das impedâncias e admitancias do modelo π
A partir do cálculo desta é possível calcular os valores reais de 𝑍 e 𝑌 da linha
por meio das equações 3.2 e 3.3 (FERNANDES, 2014):

𝑍𝜋
𝑍=
𝑠𝑒𝑛ℎ(𝛾. 𝑙)
(3.2)
Onde 𝑍𝜋 é o valor calculado da impedância em série antes da correção.

𝛾
𝑌=
𝑙
(3.3)
Os valores encontrados de 𝑍 e 𝑌 estão respectivamente em 𝛺/𝑘𝑚 e
𝑚ℎ𝑜/𝑘𝑚 ou 𝑆/𝑘𝑚 e portanto devem ser multiplicados pelo comprimento 𝑙 da linha
para cálculo de seus parâmetros totais.

3.1. Primeiro estudo – módulos de tensão, corrente e potência

disponíveis

Por meio de uma análise de circuitos, Lei de Kirchhoff das correntes, pode-se
encontrar tensão e corrente de entrada , 𝑉𝑠 e 𝐼𝑠 em função da tensão e corrente de
saída: 𝑉𝑠 e 𝐼𝑠 . Equações 3.4 à 3.6 mostram como deduzir a equação 3.7.
A corrente 𝐼𝑥 pode ser calculada pela razão entre a diferença das tensões 𝑉𝑠 e 𝑉𝑟 e a
impedância 𝑍:

𝑉𝑠 − 𝑉𝑟
𝐼𝑥 =
𝑍
(3.4)
Como 𝐼𝑟 , pela lei de Kirchhoff das correntes é dada por:

𝐼𝑟 = 𝐼𝑥 – 𝑉𝑟 . 𝑌
(3.5)

22
 Substituindo-se (3.4) em (3.5):

𝑉𝑠 − 𝑉𝑟
𝐼𝑟 + 𝑉𝑟 . 𝑌 =
𝑍
(3.6)
Desta forma, tem-se a tensão de entrada 𝑉𝑠 em função da tensão e corrente de
saída: 𝑉𝑟 (INDULKAR; RAMALINGAM ,2007):

𝑉𝑠 = 𝑍. (𝐼𝑟 + 𝑉𝑟 ). 𝑌 + 𝑉𝑟
(3.7)
A corrente 𝐼𝑥 também pode ser escrita em função de 𝑉𝑠 e 𝐼𝑠 , como representado
na equação 3.8:
𝐼𝑥 = 𝐼𝑠 – 𝑉𝑠 . 𝑌
(3.8)
E substituindo-se (3.4) e (3.7) em (3.8) pode-se obter a corrente de entrada 𝐼𝑠
em função de 𝑉𝑟 e 𝐼𝑟 .

𝐼𝑠 = (2. 𝑌 + 𝑍. 𝑌 2 ). 𝐼𝑟 + (1 + 𝑍. 𝑌). 𝐼𝑟
(3.9)
Lembrando-se que as medidas de tensão e corrente são fasoriais.
Substituindo-se os valores de 𝑍 e 𝑌 nas equações e considerando-se o diagrama
fasorial abaixo, com 𝑉𝑟 como tensão de referência, obtem-se as equações 3.10 e 3.11:

Figura 6. Diagrama fasorial

Fonte: INDULKAR;RAMALINGAM, 2007, p.2.

23
̅𝑠 = 𝑉𝑟 + 𝑗. 𝑅. 𝐵. 𝑉𝑟 − 𝐵. 𝑋. 𝑉𝑟 + 𝑅. 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 − 𝑗. 𝑅. 𝐼𝑟 , 𝑠𝑒𝑛𝜑𝑟 + 𝑗. 𝑋. 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 + 𝑋. 𝐼𝑟 . 𝑠𝑒𝑛𝜑𝑟
𝑉
(3.10)

𝐼̅𝑠 = 𝑗. 2. 𝐵. 𝑉𝑟 − 𝑅. 𝐵 2 . 𝑉𝑟 − 𝑗. 𝑋. 𝐵 2 . 𝑉𝑟 + 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 − 𝑗. 𝐼𝑟 , 𝑠𝑒𝑛𝜑𝑟 + 𝑗. 𝑅. 𝐵. 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 …

… − 𝑅. 𝐵. 𝐼𝑟 𝑠𝑒𝑛𝜑𝑟 − 𝐵. 𝑋. 𝐼𝑟 𝑐𝑜𝑠𝜑𝑟 + 𝑗. 𝐵. 𝑋. 𝐼𝑟 . 𝑠𝑒𝑛𝜑𝑟
(3.11)

̅𝑠 e 𝐼̅𝑠 foi usada para relembrar que as medidas são fasoriais e

A simbologia 𝑉
portanto as equações 3.10 e 3.11 podem ser divididas em duas equações cada.
Equações 3.12 e 3.13 mostram parte real e imaginária de 3.10. Equações 3.14 e 3.15
mostram parte real e imaginária de 3.11.

𝑉𝑠 . 𝑐𝑜𝑠𝜕 = 𝑉𝑟 − 𝐵. 𝑋. 𝑉𝑟 + 𝑅. 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 + 𝑋. 𝐼𝑟 . 𝑠𝑒𝑛𝜑𝑟

(3.12)
𝑉𝑠 . 𝑠𝑒𝑛𝜕 = 𝑅. 𝐵. 𝑉𝑟 − 𝑅. 𝐼𝑟 , 𝑠𝑒𝑛𝜑𝑟 + 𝑋. 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟
(3.13)
𝐼𝑠 . (cos(𝜕 − 𝜑𝑠 )) = −𝑅. 𝐵 2 . 𝑉𝑟 + 𝐼𝑟 . 𝑐𝑜𝑠𝜑𝑟 − 𝐵. 𝑋. 𝐼𝑟 𝑐𝑜𝑠𝜑𝑟 − 𝑅. 𝐵. 𝐼𝑟 . 𝑠𝑒𝑛𝜑𝑟
(3.14)
𝐼𝑠 . (sen(𝜕 − 𝜑𝑠 )) = 2. 𝐵. 𝑉𝑟 − 𝑋. 𝐵 2 . 𝑉𝑟 − 𝐼𝑟 𝑠𝑒𝑛𝜑𝑟 + 𝑅. 𝐵. 𝐼𝑟. 𝑐𝑜𝑠𝜑𝑟 + 𝐵. 𝑋. 𝐼𝑟 . 𝑠𝑒𝑛𝜑𝑟
(3.15)
Desta forma, têm-se 4 equações não lineares e 4 raízes a serem estimadas.
Um vetor de funções 𝐹 pode ser obtido com todos os termos das equações do lado
esquerdo igualando-os a zero: 𝐹 = [𝑓1 , 𝑓2, 𝑓3, 𝑓4, ]𝑇 e as raízes a serem estimadas são 𝑥 =
[𝑋, 𝑅, 𝐵, 𝛿]𝑇 .
Diante da não linearidade das equações, o método de Newton-Raphson foi
implementado e, como pode ser visto no exemplo da seção 2.4, o problema n-
dimensional, com um vetor formado por 4 funções e com 4 raízes a serem estimadas
a matriz Jacobiana 𝐽 é dada por:

24
 𝜕𝑓1 𝜕𝑓1 𝜕𝑓1 𝜕𝑓1
𝜕𝑋 𝜕𝑅 𝜕𝐵 𝜕𝛿
𝜕𝑓2 𝜕𝑓2 𝜕𝑓2 𝜕𝑓2
𝐽 = 𝜕𝑋 𝜕𝑅 𝜕𝐵 𝜕𝛿
𝜕𝑓3 𝜕𝑓3 𝜕𝑓3 𝜕𝑓3
𝜕𝑋 𝜕𝑅 𝜕𝐵 𝜕𝛿
𝜕𝑓4 𝜕𝑓4 𝜕𝑓4 𝜕𝑓4
[ 𝜕𝑋 𝜕𝑅 𝜕𝐵 𝜕𝛿 ]
(3.16)
Considerando os seguintes valores de medidas (INDULKAR; RAMALINGA,
2007):

𝑽𝒔 (𝑲𝑽) 127
𝑰𝒔 (𝑲𝑨) 0,416
𝑷𝒔 (𝑴𝑾) 150
𝑸𝒔 (𝑴𝑽𝑨𝒓) 50
𝑽𝒓 (𝑲𝑽) 102
𝑰𝒓 (𝑲𝑨) 0,44
𝑷𝒓 (𝑴𝑾) 135
𝑸𝒓 (𝑴𝑽𝒂𝒓) 5,7

Tabela 2. Medidas - Linha de Transmissão 1

Fonte: INDULKAR; RAMALINGAM, 2007,p.2.

Observe que por meio das medidas de potência (𝑀𝑊, 𝑀𝑉𝐴𝑟) é possível obter
os valores dos fatores de potência: 𝑐𝑜𝑠𝜑𝑠 e 𝑐𝑜𝑠𝜑𝑟 . Portanto é só substituir os valores
da Tabela 2 nas equações 3.12 à 3.15. Observe que nesse caso as medidas de ângulo
de tensão e corrente não são fornecidas. Isso acontece em situações o qual não se
tem disponível as Unidades de Medição Fasorial, (PMUs), do inglês Phasor
Measurement Units. Na ausência deste, os dados coletados pelo SCADA (Sistema de
Supervisão e Aquisição de Dados), como os da Tabela 2 devem ser usados
(KATHERINE, 2013). Com as medidas fasoriais disponíveis, o problema pode se
tornar mais simples e o método dos mínimos quadrados lineares pode ser
implementado. É o que será apresentado na próxima seção.
Após implementação do Método de Newton-Raphson com código
implementado em MATLAB e tolerância 𝜀 = 0,0001 os parâmetros 𝑥 = [𝑋, 𝑅, 𝐵, 𝛿]𝑇 ,

25
foram obtidos e foram compatíveis com os valores disponíveis (INDULKAR;
RAMALINGA, 2007), conforme Tabela 3:

𝐗(Ω) 146,4181
𝐑(𝛀) 26,7484
𝑩(𝑺) 0,0005
𝜹(°) 30,9

Tabela 3. Parâmetros obtidos com o software MATLAB

Fonte: Elaborada pelo autor

Com a tolerância estabelecida, o número de iterações realizadas foi igual a 13.

Daí a necessidade da utilização de software já que treze iterações demandam tempo
se realizadas manualmente.
Além da estimação dos parâmetros, este caso possibilita a estimação do
ângulo 𝛿 e observando-se o diagrama fasorial da Figura 7 este é exatamente o ângulo
entre as tensões 𝑉𝑠 e 𝑉𝑟 . O ângulo 𝛿 permite a análise do fluxo de potência ativa em
uma linha, como pode ser visto na Figura 7.

26
 Figura 7. Fluxo de potência ativa e o ângulo 𝛿
Fonte: http://slhaffner.phpnet.us/sistemas_de_energia_1/seI2.pdf

No caso apresentado, observa-se que com ângulo de 30,9° a transmissão de

potência ativa é da barra 𝑠 para a barra 𝑟 e ocorre de maneira estável.

3.2. Segundo estudo – medidas fasoriais de corrente e tensão

disponíveis

Diante das equações provenientes do modelo π, a estimação dependerá das

medidas que estão disponíveis. De posse de medidas fasoriais de corrente e tensão
o problema se reduz a duas equações e para utilizar o mesmo número de raízes dessa
equação utiliza-se os parâmetros 𝑍 e 𝑌. Como 𝑍 = 𝑅 + 𝑗. 𝑋 e 𝑌 = 𝐺 + 𝑗. 𝐵 , a parte real
de Z nos dá a resistência em série da linha, a parte imaginária de Z nos dá a reatância
da linha. Como 𝐺 é desprezível , o valor de 𝑌 encontrado equivale a susceptância
shunt da linha. Analisando-se a figura 5, as correntes 𝐼𝑠 e 𝐼𝑟 podem ser obtidas por
meio das equações 3.17 e 3.18:

27
 𝑉𝑠 − 𝑉𝑟
𝐼𝑠 = + 𝑉𝑠 . 𝑌
𝑍
(3.17)
𝑉𝑠 − 𝑉𝑟
𝐼𝑟 = + 𝑉𝑟 . 𝑌
𝑍
(3.18)

Equações 3.17 e 3.18 podem ser escritas na forma matricial

(FERNANDES,2014):

𝐼 𝑉 − 𝑉𝑟 𝑉𝑠 1
[ 𝑠] = [ 𝑠 ].[ ]
𝐼𝑟 𝑉𝑠 − 𝑉𝑟 𝑉𝑟 𝑍
𝑌
(3.19)

Lembre-se que os valores de 𝐼𝑠 , 𝐼𝑟 ,𝑉𝑠 e 𝑉𝑟 são medidas fasoriais.

Observe que a vetor de correntes pode ser entendida como o vetor de entrada
e o vetor de tensões como matriz de saída, buscando assim encontrar os valores de
1
e 𝑌. Claramente, 𝑍 será o inverso do valor encontrado. Com analogia a equação 2.4
𝑍

𝑌𝑚𝑒𝑑 é a matriz das correntes, 𝑋 a matriz das tensões e A contém os valores dos
parâmetros a serem estimados.
O estudo em questão utilizou as seguintes medidas fasoriais da linha de
BAURU-CABREUVA, 235km,440kV. (FERNANDES, 2014), apresentados na Tabela
4:

28
 Tabela 4. Medidas - Linha Bauru – Cabreúva
Fonte: http://www.eps.ee.kth.se/personal/luigiv/events/brazil-
workshop2014/III%20Workshop%20PMU_Fernandes.pdf , com modificações

Substituindo-se os valores de corrente e tensão da tabela na equação matricial

1
3.19, foi possível, através do método de mínimos quadrados obter a matriz: [ 𝑍 ], com
𝑌
0.0044 − 0.0117𝑖 1
valores iguais a: [ ]. Dessa forma 𝑍 = 0.0044−0.0117𝑖 = 27.8695 +
0.0000 + 0.0013𝑖
74.8670𝑖 e 𝑌 = 0.0013𝑖. Lembrando-se que esta linha é considerada uma linha longa:
𝑙 = 235𝑘𝑚 e portanto deve-se fazer a correção hiperbólica descrita nas equações 3.1
à 3.3. Dessa forma, encontra-se para os parâmetros da linha os seguintes valores,
conforme Tabela 5:

𝑹(Ω) 28.7648
𝑿(Ω) 75.8856
𝑩(𝑺) 0.0013

Tabela 5. Parâmetros encontrados - Linha Bauru - Cabreúva

Fonte: http://www.eps.ee.kth.se/personal/luigiv/events/brazil-
workshop2014/III%20Workshop%20PMU_Fernandes.pdf , com modificações

29
 Capítulo 4

Medidas disponíveis – Sistemas de medição em

Sistemas Elétricos

4.1. Sistema de Medição Fasorial

Diante das possíveis estimações apresentadas no capítulo 3 este capítulo

apresenta e explica os equipamentos utilizados na medição das grandezas do sistema
elétrico. Como nova tendência dos Sistemas de Gerenciamento de Energia (Energy
Management System) os Sistemas de Medição Fasorial (SISMEF) vem trazendo
novas possibilidades nos estudos do sistema elétrico.
Observe que na seção 3.2 as medidas disponíveis são Correntes e Tensões
Fasoriais e a aquisição dessas medidas é feita por meio do dispositivo presente no
SISMEF, denominado Unidade de Medição Fasorial (PMU).
O PMU é o principal elemento de um SISMEF. Primeiramente, esse elemento
adquire os valores de corrente e tensão e então estas são convertidas em fasores,
podendo assim ser enviadas a um dispositivo que concentra esses dados. A figura 8
mostra aquilo que constitui um PMU.

Figura 8. Elementos constituintes de um PMU

Fonte: BENETTI, 2009, p. 9.
30
 Por meio da figura 8 observa-se que um PMU é constituído por um sistema de
aquisição; filtros anti-aliasing, circuito oscilador e módulo conversor Analógico/Digital
(A/D); e por um microprocessador. É este que efetua o tratamento matemático das
amostras. Para fornecimento de uma referência de tempo ao conversor D/A, o PMU
utiliza um módulo o qual recebe e processa o sinal do GPS. Alguns PMUs possuem
um buffer de memória local. Através deste é possível a gravação digital de
perturbações.
Primeiramente, para que um processo de aquisição de dados seja
sincronizado, independentemente da localização das subestações nas quais estão
situadas os PMUs torna-se necessário a utilização de GPS (Sistema de
Posicionamento Global) para envio de um PPS (Serviço de Posicionamento Preciso).
A precisão é maior que que 1 𝜇𝑆 e o sinal é enviado com precisão igual para estações
distantes (BENETTI, 2009).
As aquisições simultâneas são feitas a partir dos secundários dos
Transformadores de Potencial e de Corrente. Após conversão A/D, emprega-se
diferentes métodos, com o objetivo de se efetuar o cálculo das grandezas fasoriais. A
ferramenta matemática de maior uso nessas aplicações é a Transformada Discreta
de Fourier. Segundo Benetti (2009), após processamento, os fasores são então
enviados ao Concentrador de Dados (Phasor Data Concentrator - PDC). Para que se
tenha um desempenho de qualidade do SISMEF, é necessário que o PDC apresente
uma série de características, dentre elas:
 Robustez para suporte das falhas de hardware e software evitando possíveis
perdas de informação;
 Por conta da presença de vários PMUs nas várias estações, prever
expansão do processamento de dados;
 Adquirir, processar e disponibilizar os dados fasoriais de forma simultânea;
A medição fasorial sincronizada é uma tecnologia desenvolvida para aumentar
a confiabilidade dos sistemas elétricos diante de um cenário que as redes elétricas ao
redor do mundo tem sido operadas em condições próximas aos limites de sua
capacidade. Pode ser considerada uma técnica recente quando os anos de início de
implantação do sistema nos Estados Unidos e na China por exemplo são analisados:
2003 e 2002 respectivamente (BENETTI, 2009). Com pesquisa desde 1994 e

31
implantação a mais de 10 anos os chineses apontam algumas funções básicas dos
PMUs:
 Monitoração em tempo real;
 Monitoração de distúrbios;
 Integração dos dados do sistema elétrico;
 Monitoração da operação de geradores;
 Utilização de medidas fasoriais na estimação de estados;
 Aprovação de modelos para simulações;
 Identificação e localização de faltas;
Na identificação e localização de faltas a estimação de parâmetros de linhas de
transmissão é utilizada e seção 4.3 explica como este processo acontece. Antes
seção 4.2 mostra a utilização do SISMEF no Brasil

4.2. Utilização do SISMEF no Brasil

A Figura 9 mostra a estimação dos parâmetros em tempo real exibida no

openPDC Manager, Concentrador de Dados da General Electric (GE) para a linha de
Cabreúva – Baurú, a mesma da tabela 3.

32
 Figura 9. Parâmetros em tempo real
Fonte: http://www.eps.ee.kth.se/personal/luigiv/events/brazil-
workshop2014/III%20Workshop%20PMU_Fernandes.pdf

O projeto pioneiro no Brasil é o Projeto MedFasee, o qual teve início em 2003

com parceria entre a Universidade de Santa Caratina (UFSC), a empresa Reason
Recnologia S.A. e a Financiadora de Estudos e Projetos (FINEP). Neste projeto ,3
PMUs foram instalados nas capitais do três estados do sul do Brasil afim de monitorar
grandezas do sistema elétrico brasileiro. Os fasores de tensão eram enviados através
da internel ao PDC instalado no Laboratório de Planejamento de Sistemas de Energia
Elétrica (LabPlan) em Florianópolos, Santa Catarina. A Figura 10 mostra os três locais
nos quais foram instalados os PMUs.

33
 Figura 10. Tela de monitoração em tempo real do Projeto MedFasee
Fonte: BENETTI, 2009, p. 21.

Com a instalação de todos os equipamentos, o Projeto MedFasee terminou

de fato em 2009. A partir deste projeto, alcançou-se desenvolvimentos acadêmicos
no que tange a melhoria da estabilidade angular dos sinais dos sistemas elétricos bem
como a análise da dinâmica destes.
Em Minas Gerais, a Companhia Energética de Minas Gerais (CEMIG) já tem
implantado um SISMEF desde Outubro de 2010. Por meio do sistema, a companhia
tem realizado a avaliação da estabilidade de tensão, permitindo o controle preventivo
e o restabelecimento de ilhas e linhas, para casos de desequilíbrio.

4.3. Localização de Faltas do tipo Curto-Circuito e a Estimação de

Parâmetros

A rápida localização de faltas em Sistemas de Transmissão para o

reestabelecimento da operação normal do sistema é extremamente importante. A
diminuição do número de linhas de transmissão pode levar um sistema à instabilidade
e agentes reguladores podem penalizar as companhias energéticas caso não
mantenham a continuidades desses serviços.

34
 Técnicas de localização de faltas utilizam a estimação de parâmetros das linhas
de transmissão em tempo real. Observe, na Figura 11, que a estimação está presente
no algoritmo de Silveira (BENETTI, 2009).

Figura 11. Fluxograma do Método Proposto por Silveira

Fonte: BENETTI, 2009, p.82.

Algoritmos de localização de falta como o de Silveira realizam um processo

de varredura da rede, onde se busca minimizar o erro entre o ponto estimado de falta
e o ponto real (SANTO; PEREIRA; ZANETTA, 2011). Como a cada estimação durante
a varredura novos parâmetros estão associados aos novos cálculos, tem-se a
necessidade de estimação de parâmetros de linhas de transmissão em tempo real. E
essa é uma das funcionalidades disponibilizadas pelos PMUs.

35
4.4. Utilização de Unidades Terminais Remotas na aquisição de dados
do Sistema Elétrico

Um Sistema de Supervisão e Aquisição de Dados (SCADA) funciona como um

grande banco de dados das medições do sistema elétrico. Dados que são colhidos
das Unidades Terminais Remotas (Remote Terminal Units – RTUs), relés de proteção
e transdutores análogos aos PMUs, apresentados na seção 4.3. Acontece que,
diferentemente dos PMUs, os dados nas RTUs são colhidos de forma assíncrona e
lenta o que não permite informações precisas como os ângulos de tensão e corrente
nos dois terminais de uma linha. A Figura 12 mostra um Sistema de Gerenciamento
de Energia com RTUs e aquisição dos dados realizada pelo SCADA:

Figura 12. Sistema de Gerenciamento de Energia

Fonte: ALSTOM

Diferentemente das amostras por segundo nos PMUs as medidas provenientes

das RTUs podem ter amostras de 15 em 15 minutos, não sendo possível a observação
de perturbações de curta duração na rede.

36
4.5. Conclusão acerca do uso das medidas para estimação de
parâmetros

Na seção 3.1, primeiro estudo, pode-se observar que as equações trabalham

com os módulos de corrente e tensão e os ângulos do modelo podem ser obtidos por
meio de medições de potência ativa e reativa. Claramente, o modelo também é
aplicável no caso de disponibilidade medições fasoriais, porém de posse somente de
módulos de corrente, tensão e potência o estudo da seção 3.2 não é possível.
O capítulo 5 mostra um estudo de caso em que se estima os parâmetros reais
de uma linha de transmissão por meio de dados colhidos no SISMEF da CEMIG.

37
 Capítulo 5

Estudo de Caso – Linha de Transmissão da CEMIG:

Mesquita - Baguari

Este estudo consiste em replicar as técnicas utilizada nas seções 4.1 e 4.2 de
forma a comparar os parâmetros calculados com os parâmetros reais de uma linha de
transmissão da CEMIG, e analisar o fluxo de potência nesta. Na seção 4.2 o método
foi implementado considerando-se um vetor de corrente e um vetor de tensão cada
um com 13 medidas. Nesta seção, os parâmetros são estimados para cada 1 minuto
no período entre dia 29 e 30 de junho de 2015 e dessa forma (com um total de 1440
amostras). A partir desta estimação, pretende-se analisar o comportamento dos
parâmetros da linha de transmissão. A linha Mesquita – Baguari possui comprimento
de 67,51 km e apesar de seu comprimento estar inserido no padrão de linhas curtas
a linha é considerada de média de distância. Os valores delimitados de 80 à 240 km
são teóricos e como pode ser visto nos casos de referência do Operador Nacional do
Sistema (ONS) Elétrico bem como nos dados disponibilizados pela CEMIG a linha
possui susceptância shunt. A linha Mesquita – Baguari possui nível de tensão igual à
230 kV e portanto é parte da chamada rede básica: sistema constituído por todas as
subestações e linhas de transmissão em tensões de 230 kV ou superior (Agência
Nacional de Energia Elétrica[ANEEL],1998). Os parâmetros disponibilizados pelo
ONS se encontram na Tabela 6. A linha de transmissão estudada é apresentada na
Figura 13:

38
 Figura 13. Linha de transmissão: MESQUITA - BAGUARI

R (%): X (%) B (%):

1.43 6.64 11.79
Tabela 6. Parâmetros – ONS
Fonte: http://www.ons.org.br/plano_ampliacao/casos_refer_regime_permanente.aspx

Observe que os valores dos parâmetros disponibilizados pela ONS se

encontram em porcentagem sendo que a tensão de base e potência de base são
respectivamente 230 kV e 100 MVA. Dessa forma, a impedância Z de base será:

𝑉𝑏𝑎𝑠𝑒 2
𝑍𝑏𝑎𝑠𝑒 = = 529𝛺
𝑆𝑏𝑎𝑠𝑒
(5.1)
E admitância base:

1
𝑌𝑏𝑎𝑠𝑒 = = 1,89𝑚𝑆
𝑍𝑏𝑎𝑠𝑒
(5.2)
Para obter os valores reais dos parâmetros, multiplica-se os valores em
porcentagem pelas bases utilizadas:

39
 1,43 × 529
𝑅= = 7.5647𝛺
100
(5.3)
6,64 × 529
𝑋= = 35,1256𝛺
100
(5.4)
15,89560 × 0,00189
𝐵= = 222,831µ𝑠
100
(5.5)
Os valores reais fornecidos pela CEMIG condizem com os do NOS, conforme
Tabela 7:

R (Ω): X (Ω) B/2(µS):

7,58100 35,14300 111,41000
Tabela 7. Parâmetros – CEMIG
Fonte: CEMIG

Os dados fornecidos pelo Sistema de Medição Fasorial (SISMEF) da CEMIG

consistem nas medições de tensão e corrente e seus respectivos ângulos. O projeto
do SISMEF começou em Outubro de 2010 e utiliza os PMUs instalados nas
subestações da CEMIG. Eles são sincronizados com o Sistema de Posicionamento
Global (GPS) os quais providenciam a posição e o tempo da informação permitindo
com que os operadores do sistema tenham atualização mais precisa deste
(ANDRADE, 2013).

5.1. Newton- Raphson e a estimação do ângulo da corrente 𝑰𝒓

Esta seção segue os procedimentos da seção 4.1. É importante ressaltar a

necessidade de adaptação das medidas da CEMIG e observância das similaridades
e divergências em relação ao modelo apresentado na seção 4.1. Primeiramente, para
identificação dos sentidos das correntes 𝐼𝑠 e 𝐼𝑟 ,representadas no modelo da Figura 5,
os casos de referência do ONS foram consultados e a partir do manual do ANAREDE,
programa de análise de redes mais utilizado em sistemas elétricos de potência, foi
possível a verificação dos sentidos. Os arquivos do caso de referência possuem dois
40
códigos principais: DBAR e DLIN, os quais representam todas as barras e linhas da
rede básica respectivamente. Cada barra possui uma identificação numérica e como
pode ser visto nas Figuras 14 e 15 as barras de Baguari e Mesquita puderam ser
identificadas na seção DBAR:

Figura 14. Barra – Baguari

Fonte: http://www.ons.org.br/plano_ampliacao/casos_refer_regime_permanente.aspx

Figura 15. Barra – Mesquita

Fonte: http://www.ons.org.br/plano_ampliacao/casos_refer_regime_permanente.aspx

As barras puderam ser identificadas através da numeração à direita, as quais

indicam o nível de tensão de 230 𝑘𝑉 Identificadas as barras verificou-se a presença
da linha na seção DLIN. Por meio desta seção foi possível a verificação do sentido da
corrente, já que, as colunas (De) e (Para) da seção DLIN indicam o sentido do fluxo.
A Figura 16 mostra a seção DLIN do arquivo do ONS.

Figura 16. Seção DLIN - Arquivo ONS

Fonte: http://www.ons.org.br/plano_ampliacao/casos_refer_regime_permanente.aspx

41
 Dessa forma, com Baguari no lado r, tem-se como referência a corrente saindo
da barra e indo para o lado s. Caso o ângulo entre 𝐼𝑟 e 𝑉𝑟 apresente comportamento
conforme situação da tabela 4 (𝜑𝑟 > 90°), a corrente tem sentido oposto ao sentido
de convenção. Para o caso da CEMIG, saindo da barra. Por esse motivo de inversão
do sentido de corrente equações 3.7 e 3.9 ficam da seguinte forma:
𝑉𝑠 = (1 + 𝑍𝑌). 𝑉𝑟 − 𝑍. 𝐼 𝑟
(5.6)
𝐼𝑠 = (2. 𝑌 + 𝑍. 𝑌 2 ). 𝑉𝑟 + (−1 − 𝑍. 𝑌). 𝐼𝑟
(5.7)
As medidas colhidas pelos PMUs mostram que entre às 18 horas do dia 29 e 2
da manhã do dia 30 𝜑𝑟 > 90° e portanto as equações 5.6 e 5.7 foram aplicadas.
Para as medições da CEMIG, ao invés da estimação do ângulo entre as
tensões preferiu-se a estimação do ângulo da corrente 𝐼𝑟 .

O ângulo entre as tensões dessa linha foi verificado como muito pequeno.
Dessa forma, para efeitos de eficiência do método escolheu-se um valor maior de
ângulo para que os parâmetros não tivessem mais de uma solução. Isso porque o
valor da Susceptância Shunt da linha é muito pequeno e no caso de duas soluções
muito pequenas as outras duas soluções poderiam assumir mais de um valor, o que
torna o estudo da seção 3.1 um pouco ineficiente. No exemplo dado, o ângulo 𝛿 entre
as tensões é quase que 10 vezes o ângulo entre as tensões da linha da CEMIG.
Tabela 8 mostra os conjuntos de soluções obtidos para as cinco primeiras medidas
representadas na tabela 9:

42
 𝑿(𝑶𝒉𝒎𝒔) 56,4302 56,5304 56,8461 57,2438 57,3895
𝑹(Ohms) 23,9981 24,0747 24,0909 24,0087 24,2958
𝑩(S) 0,0000 0,0000 0,0000 0,0000 0,0000
𝝋𝒓 (º) 197,48 197,38 197,21 197 197,05

Tabela 8. Parâmetros obtidos no método de Newton Raphson

Fonte: Elaborada pelo autor

BGR-VB1 BGR-IB1 MQT-VB1 MQT-IB1

BGR-VB1 BGR-IB1 MQT-VB1 MQT-IB1
Voltage Current Voltage Current
Voltage Current Voltage Current
Phase Phase Phase Phase
Magnitude Magnitude Magnitude Magnitude
Angle Angle Angle Angle
-
134735,10 160,2289 112,369 135745,90 2,9251 110,812 14,4732
178,6153
-
134725,50 158,4709 112,132 135745,90 4,6831 110,634 16,2312
176,8573
-
134701,70 155,6180 111,514 135719,90 7,5360 110,115 19,0841
174,0044
-
134678,10 151,6935 110,791 135683,60 11,4605 109,517 23,0086
170,0799
-
134673,50 147,4699 110,423 135698,70 15,6841 109,123 27,2322
165,8563
-
134749,40 143,5641 109,236 135752,20 19,5899 107,893 31,1380
161,9505

Tabela 9. Cinco primeiras medições realizadas pelo PMU

Fonte: Elaborada pelo autor

Observe que, apesar de distantes dos parâmetros de projeto, os cálculos por

Newton – Raphson mostram precisão. A solução, depois de um minuto, mesmo tempo
colhido na utilização de mínimos quadrados, possui valores próximos dos primeiros 5
segundos, conforme observado na Tabela 10:

43
 𝑿(𝑶𝒉𝒎𝒔) 58,9424
𝑹(Ohms) 25,1499
𝑩(S) 0,0000
𝝋𝒓 (º) 197,52

Tabela 10. Medição após 1 minuto

Fonte: Elaborada elo autor

Para validação do método, substituiu-se os valores dos parâmetros nas

equações 3.7 e 3.9 provando-se que os valores obtidos são as raízes das equações
do modelo. Maiores análises acerca do método serão feitas no capítulo 6.

5.2. Método dos mínimos quadrados e análise dos parâmetros

Nesta seção, o método dos mínimos quadrados, conforme feito no estudo da

seção 3.1 foi implementado para obtenção dos parâmetros da linha de transmissão
da CEMIG. A partir da estimação de minuto em minuto, feita com 12 amostras em
cada, foi possível traçar curvas de representação da variação dos parâmetros ao
longo do dia. (Figuras 17 à 19). Repare que algumas regiões apresentam variações
bruscas dos parâmetros chegando a valores de grande proporção. A análise e motivos
dessas variações bem como o detalhamento por faixas de tempo das Figuras 17 à 19
se encontram no capítulo 6.

Figura 17. Variação da Resistência ao longo do dia

44
 Figura 18. Variação da Reatância ao longo do dia

Figura 19. Variação da Susceptância ao longo do dia

45
 Capítulo 6

Análise dos Resultados – Variações dos Parâmetros da

linha

6.1. Resistência em Série

Após obtenção dos valores dos parâmetros da linha de transmissão:

MESQUITA – BAGUARI, primeiramente analisou-se as variações dos valores de
resistência em relação ao projeto da linha, presente na tabela 6. O principal motivo na
variação de resistência é a temperatura do condutor. Condutores de linha de
transmissão são geralmente de alumínio principalmente pelo custo e peso reduzido,
é o caso da linha em questão. Além disso, possui menor tendência à ionização do ar.
A equação 6.1 mostra a variação da resistência de um condutor em função de sua
temperatura:

𝑅 = 𝑅20 (1 + 𝛼20 (𝑇 − 20))

(6.1)
Onde:
 𝑅20 = Resistência à 20ºC;
 𝛼20 = Coeficiente de temperatura à 20ºC;
 𝑇 = Temperatura real
Admitindo-se que o projeto da linha considera o valor de resistência à 20ºC, e
que 𝑅 é o valor real do parâmetro, valor calculado pelo método dos mínimos
quadrados, pode-se determinar a temperatura do condutor para cada ponto isolando-
se T na equação 6.1 e dessa forma:

𝑅 + 𝑅20 × (20 × 𝛼20 − 1)

𝑇 =
𝑅20 × 𝛼20
(6.2)

Esse valor de temperatura obtido considera somente as características do

condutor. Diante disso, o que se esperaria em um gráfico de dispersão: Resistência x

46
Corrente, é a maior ocorrência de resistências altas em temperaturas altas, já que a
corrente alta aumenta a temperatura do cabo. Porém, o que se sabe, é que a
resistência calculada pode também ter seu valor associado a temperatura ambiente.
Os gráficos apresentados nas figuras 20 e 21 abaixo mostram respectivamente a
variação da resistência ao longo do dia no horário de 18:00, do dia 29, até 02:00, do
dia 30 de junho de 2015 e depois o gráfico de dispersão da resistência em função da
corrente. A correlação positiva entre maiores valores de resistência e entre maiores
valores de corrente não pode ser comprovada, porém em horários como 20:30 - 22:00
supõe-se que a resistência diminui em consequência da diminuição da temperatura
do dia. Pela variação apresentada nas Figuras 20 e 21 não se pode observar relações
diretas entre resistência e temperatura do dia.
As marcações de cores vermelha, verde e amarela associam o gráfico da figura
20 ao gráfico da figura 21. Pela marcação vermelha na Figura 20 observa-se que os
valores de resistência de 16Ω permanecem até o horário de aproximadamente 19:45
e nesse momento, de acordo com a figura 21, as correntes apresentam valores na
faixa de 100ª. Variação entre 12 e 14 Ω (marcação verde), aproximadamente, é
observada com corrente variando entre 200 A e 100 A respectivamente, o que não
era esperado. Isto porque, ressaltando novamente, altas correntes interferem na
temperatura do cabo, e consequentemente na alta resistência deste. A marcação
amarela mostra os maiores valores de resistência no período, na Figura 20, resultando
de dados estranhos às outras tendências, representados na Figura 21. Sugere-se,
diante das variações de resistência, estudos com detalhamento para investigação
desse fenômeno.

Figura 20. Variação da resistência

Período: 18:00 (29/06/2015) - 02:00 (30/07/2015)

47
 Figura 21. Dispersão da resistência em função da corrente
(18:00 - 02:00)

O histograma da Figura 22 mostra a maior frequência de ocorrência em

valores de, aproximadamente 16 Ω para a resistência. Figura 23 mostra a variação da
corrente no período de 18:00 às 2 da manhã. Repare que, conforme já observado a
resistência da linha não acompanha a corrente e entre 20 e 22 horas por exemplo,
nos valores de maior corrente é onde ocorre as menores resistências.

Figura 22. Histograma – Resistência

Período: 18:00 (29/06/2015) - 02:00 (30/07/2015)

48
 Figura 23. Variação da Corrente
Período:18:00 (29/06/2015) - 02:00 (30/07/2015)

Figuras 24 e 25 mostram a mesma situação e agora para o período de 12:00

até às 18:00 do dia 30 de junho de 2015. Regiões em vermelho, verde claro e verde
escuro apresentam as regiões de maior permanência do valor de corrente em 13 Ω,
15, 5 Ω e 15 Ω aproximadamente, respectivamente. Valores que estão associados a
120 A, 130A e 140 à 150 A aproximadamente, conforme observado na Figura 25. A
região destacada em amarelo, na Figura 24 mostra o salto do valor de resistência de
13 A para 14 A aproximadamente. Ou seja, região de pouca ocorrência no gráfico de
dispersão da figura, com valores entre 100 e 110 A. A partir do gráfico da Figura 25
pode-se observar a tendência da correlação positiva entre corrente e resistência por
conta do aquecimento do condutor. Neste momento, as possíveis temperaturas altas
podem estar contribuindo para o maior valor de resistência no período.

Figura 24. Variação da resistência

(12:00 - 18:00)
49
 Figura 25. Gráfico de dispersão da resistência
(12:00 - 18:00)

6.2. Reatância em série

Como pode ser observado, também para reatância, os valores encontrados são
maiores que os valores reais. A variação da reatância, em caso de condutores
compostos, que possuem dois ou mais fios, está relacionada a distância entre os fios
de um mesmo condutor, e a distância entre os fios de dois ou mais condutores
presentes em uma linha. Figura 26 mostra as distâncias citadas:

Figura 26. Condutores na linha de transmissão

Como a reatância indutiva 𝑋 é diretamente proporcional a indutância, 𝐿, para

uma frequência constante 𝑓, 𝑋 = 2𝜋𝑓𝐿 a análise do aumento da reatância pode ser

50
dado considerando as equações de indutância, apresentadas abaixo. As equações
consideram dois condutores inicialmente, os mesmos apresentados na Figura 26:

𝑀
𝜇 √𝐷𝑎𝐴 𝐷𝑎𝐵 𝐷𝑎𝐶 … 𝐷𝑎𝑀
𝐿𝑎 = 𝑛 𝑙𝑛 𝑛
2𝜋 √𝑟𝑎′ 𝐷𝑎𝑏 𝐷𝑎𝑐 … 𝐷𝑎𝑀
(6.3)
𝑀
𝜇 √𝐷𝑏𝐴 𝐷𝑏𝐵 𝐷𝑏𝐶 … 𝐷𝑎𝑀
𝐿𝑏 = 𝑛 𝑙𝑛 𝑛 ′
2𝜋 √𝑟𝑏 𝐷𝑏𝑎 𝐷𝑏𝑐 … 𝐷𝑏𝑛
(6.4)

A razão, para aumento da reatância em relação ao valor do projeto, é,

inicialmente e hipoteticamente a mesma apresentada para variação da resistência. Se
a linha, no dia apresentado, possui maior carregamento, existe a possibilidade das
distâncias apresentadas nas equações 6.3 e 6.4 se alterarem devido ao aquecimento
dos condutores. Se a reatância aumentou, sugere-se que a distância entre os
condutores aumentou mais do que a distância entre os fios de um mesmo condutor.
Além disso, com o aquecimento, o cabo pode ter seu comprimento aumentando e o
maior número de loops de corrente aumenta a indutância nos cabos.
Os gráficos apresentados nas Figuras 27 e 28 mostram a variação da reatância
e o gráfico de dispersão desta em função da corrente entre às 18 horas do dia 29 e
02 da manhã do dia 30.

Figura 27. Variação da reatância

Período: 18:00 (29/06/2015) - 02:00 (30/07/2015)

51
 Figura 28. Gráfico de dispersão da reatância

(18:00 - 02:00 )

Repare que, acompanhando-se a região de permanência dos valores da

resistência, as reatâncias apresentam, no mesmo horário e em mesma corrente
valores iguais a aproximadamente 65 Ω. Histograma na Figura 29 mostra a
distribuição normal em torno do valor de 65 Ω.

Figura 29. Histograma - reatância

(18:00 - 02:00)

52
 Aproveita-se aqui para observar que os valores médios encontrados no período
entre 18 horas do dia 29 e 2 horas do dia 30, para resistência: 16 Ω e reatância 65 Ω
equivalem a aproximadamente o dobro do projeto da linha: resistência de 7,58 Ω e
reatância de 35,14. Sugerindo-se assim a empresa prestadora dos dados deste
trabalho, CEMIG, a revisão dos valores presentes

6.3. Susceptância Shunt

O valor da susceptância shunt é, teoricamente, zero. A linha utilizada por

possuir distância menor que 80 km é uma linha curta. Acontece que não se pode
considerar esses limites de forma pragmática. Como pode ser observado, a própria
CEMIG considera um valor de susceptância shunt para a linha. No método de Newton-
Raphson, pouco inexato em relação aos parâmetros obtidos no método dos mínimos
quadrados, valores de susceptância shunt são iguais a zero o que mostra que a
solução inicial realmente partiu para iterações que chegassem a uma solução
plausível. A variação da susceptância shunt no período de 18 horas do dia 29 e 2
horas do dia 30 de junho é apresentada na Figura 30. Repare que o eixo das
ordenadas no gráfico apresenta unidade com ordem de grandeza de 10−5 :

Figura 30. Variação da susceptância

(18:00 - 02:00)

53
6.4. Métodos numéricos

O método dos mínimos quadrados apresentam realmente os parâmetros reais

da linha por se tratar, no caso, quase que um cálculo direto dos parâmetros. Já o
método de Newton Raphson se mostra dependente de um bom valor inicial para dar
prosseguimento as iterações. O estudo apresentado na seção 3.1 se mostrou
questionável uma vez que o ângulo 𝛿 entre as tensões neste exemplo é grande. O
ângulo entre as tensões na linha de transmissão da CEMIG é na maioria das vezes
10 vezes menor que 30 graus, ângulo do caso da seção 3.1 Dessa forma, para o
conjunto de parâmetros a serem estimados na seção 3.1, tem-se dois valores muito
pequenos, e por isso, entendeu-se que a resistência e a reatância poderiam assumir
mais de um valor. Por isso a escolha de estimação de um ângulo maior, no caso o
ângulo da corrente 𝐼𝑟 .

6.5. Sincronização dos PMUs

Os valores estranhos observados, principalmente no horário entre às 2 da

manhã e meio-dia do dia 30 de Junho estão atrelados a variação dos fasores das
medidas. Com comportamento conforme observado na Figura 31 os ângulos entre as
tensões e correntes da barra variam de tal forma que se sugere estudos para
compreensão destes de forma a possibilitar o entendimento da estimação de
parâmetros por meio de medições fasoriais. Como pode ser visto ao longo do capítulo
6, somente em algumas regiões do dia, nas quais era possível a implementação das
equações do circuito π foi possível a estimação dos parâmetros

54
Figura 31. Mudanças de ângulo no PMU

55
 Capítulo 7

Conclusão

Este trabalho possibilitou entender que estimação de parâmetros envolve muito

mais do que o simples cálculo destes. O processo de estimação envolve escolha e
aplicação de métodos numéricos e diante das possíveis ineficiências destes deve-se
compará-los de modo a compreender possíveis erros de estimação. Muitas vezes, os
erros não estarão na implementação dos métodos. O processo envolve avaliação de
causas da inexatidão ou imprecisão dos resultados.
A tarefa de trabalhar com dados reais instiga a compreensão também dos
processos de medição e os equipamentos e tecnologias utilizados para esta. O estudo
das PMUs e UTRs e seus métodos de funcionamento para estudo é extremamente
importante. Como pode ser visto, cada tipo de dado fornecido dirige a um modelo, a
cada tipo de modelo dirige a um método numérico.
Ao longo do trabalho, pode-se compreender a relevância deste tema e a
aplicação para localização de faltas é um trabalho sugerido que como pode ser visto,
gera resultados de confiabilidade do sistema bem como diminuição de prejuízos às
companhias. Encoraja-se também, o avanço dos estudos das possibilidades que o
Sistema de Medição Fasorial oferece. Observou-se que este apresenta uma série de
aplicações e como é algo relativamente novo, estudos sobre o sistema vem avançado
e gerando possibilidades ao sistema elétrico.
Este trabalho trouxe conhecimento e se mostrou interdisciplinar, com conceitos
que envolviam Sistemas Elétricos de Potência, Planejamento e Operação do Sistema,
Eletrotécnica, Matemática e Computação. Desta forma, o trabalho se mostra relevante
à sociedade científica.

56
 Referências Bibliográficas

Albertini, M. (2010). Metodologia para depuração off-line de parâmetros série e shunt de

linhas de transmissão através de diversas amostras de medidas. Área de Concentração
de Sistemas Elétricos de Potência. Escola de Engenharia de São Carlos da Universidade de
São Carlos.

Azevedo M.(2013). Método dos Mínimos Quadrados. Curso MAT01169. Departamento de

Matemática Pura e Aplicada.

Beiguelman B.(2011). A distribuição normal. Disponível em:

<http://www.ufpa.br/dicas/biome/bionor.htm>. Acesso em 29 de outubro de 2014.

Benetti, M. (2009). Utilização de Medição Fasorial Sincronizada em Metodologias de

Localização de Faltas em Linhas de Transmissão. Programa de Pós-Graduação em
Engenharia Elétrica. Universidade Federal de Santa Catarina.

Fernandes R.(2014). Estabilidade de tensão em regime permanente, oscilação

eletromecânica, estimação de parâmetros de linha com uso de PMU’s, em tempo real e
perspectvas da IEC/TR 61850-90-5. Disponível em: <http
://www.eps.ee.kth.se/personal/luigiv/events/brazilworkshop2014/III%20Workshop%20PMU_
Fernandes.pdf>. Acesso em 17 de Março de 2015.

FILHO, Képler O.(2013). Mínimos Quadrados. Departamento de Astronomia do Instituo de

Física. Universidade Federal do Rio Grande do Sul.

Haffner, S.(2007). Modelagem e Análise de Sistemas Elétricos em Regime Permanente:

A linha de transmissão. Disponível em:
<http://slhaffner.phpnet.us/sistemas_de_energia_1/seI3.pdf>. Acesso em 31 de março de
2015.

Indulkar, C.; Ramalingam, K. (2007). Estimation of Transmission Line Parameters from

Measurements. Departameto de Engenhaia Elétrica e Complexo ATS, India.

57
Katherine R., Tomas J., et al. (2013). Estimation of Line Parameters from Historical Data.
46th Hawaii Internation Conference On System Sciences.

Meister, D. (2006). Aplicação do Método dos Mínimos Quadrados na estimação de

Parâmetros do modelo de um transformador. Departamento de Engenharia Elétrica.
Faculdade de Tecnologia. Universidade de Brasília. Universidade Estadual do Ceará.

Mendonza, M. (2011). Identificação de Parâmetros de Geradores Síncronos.

Departamento de Sistemas e Controle de Energia. Faculdade de Engenharia Elétrica e de
Computação. Universidade Estadual de Campinas

Neves, M. (2008). Modelagem de Carga em Sistemas de Energia Elétrica: Modelo

matemática e sua validação com testes de campo. Faculdade de Engenharia.
Universidade Federal de Juiz de Fora.

Notas de aula: Professora Patrícia Romeiro da Silva Jota. – Slides de SEP.

Pereira, Fernando S.(2005). Análise de Observabilidade Para o Estimador de Estados e

Parâmetros. Escola de Engenharia São Carlos.

Santo, S.; Pereira, C.; Zanetta, L. (2011). Localização De Faltas Em Linhas De

Transmissão Com Múltiplas Derivações Utilizando Dados De Somente Dois Terminais.
Revista Controle & Automação/Vol.22 no.2/Março e Abril 2011.

Viana, G. (2001). Métodos iterativos para resolução de sistemas de equações não

lineares. Faculdade Lourenço Filho.

58