Modelado de relés de Distancia

en ATP-EMTP para Estudios de

Protecciones en Sistemas de
Potencia
MSc. Ing. Rubén Guamán Medina
ENDE TRANSMISION S. A.
Cbba, 30-Octubre-2024
Escenarios de simulación en SEP
ESCENARIOS DE SIMULACION

ESTADO ESTABILIDAD TRANSITORIOS

ESTACIONARIO TRANSITORIA ELECTROMAGNETICOS

• ESTABILIDAD ANGULAR
• FLUJOS DE CARGA
• ESTABILIDAD VOLTAJE • SOBRETENSIONES POR
• CORTOCIRCUITOS DESCARGAS ATMOSFERICAS
• OSCILACION DE POTENCIA
• PROTECCIONES • TRV DE INTERRUPTOPTORES
• PROTECCIONES
• COMPORTAMIENTO MAGNETICO DE
EQUIPOS
• SATURACIÓN DE NUCLEOS DE
TRANSFORMADORES
• FERRORESONANCIA

Programas basados en fasores •

•
PROTECCIONES
ETC
Escenarios de simulación en SEP
TRANSITORIOS
ELECTROMAGNETICOS

PROGRAMAS DE PROGRAMA EMT PARA LA

EMTP MODELANDO RED Y
PROPOSITO RED Y PROGRAMA
RELE
GENERAL EXTERNO PARA EL RELE

• MODELADO • MODELADO
MODELADO DE LA RED DEL RELE • MODELA LA
PROTECCION Y LA PROTECCION
RED • MODELA LA RED

• ATP/EMTP • MATLAB
• PSCAD • MATHCAD
• MICROTRAN
• ATP/EMTP
• MATLAB • EMTP-RV
• PSCAD
• MATHCAD • Otros
• ETC
Objetivos de la Investigación
El objetivo de este trabajo de tesis fue
proponer un modelo genérico de relé
distancia, y establecer una metodología válida
incluso para modelar relés numéricos
detallados en el ambiente de simulación de
transitorios electromagnéticos, para análisis de
desempeño de los mismos ante eventos de
contingencia eléctricos.
Aplicar el modelo desarrollado al estudio de
fallas en una red reducida de ENDE
Transmisión.
Principio de operación
Z3

Z2
t
Z1
tz2 Z1 tz3
tz1
0 1 2 3 4 5
Margen de
Coordinación
USC V=0
RTC Z1
3
ISC
21
RTP ഥ𝑆𝐶
𝑈
𝑍𝐹ҧ = No es un relé de distancia
ҧ
𝐼𝑆𝐶
Conexión de Transformadores de
corriente y tensión
a
Barra b
c

TP

Ua Ub Uc
Relé 21
Ia Ib Ic

Ia

Ic
TC Ib

52
Línea

Conexión universal de un relé numérico digital

Impedancias de secuencia de una
línea de transmisión
En la figura muestra una línea perfectamente
simétrica, es decir que tiene una transposición
completa y todas sus impedancias propias son iguales
y también sus impedancias mutuas.
Acoplamiento Magnético
𝑰𝒂 ഥ𝒔
𝒁
a
𝑰𝒃 ഥ𝒔
𝒁 ഥ𝒎
𝒁 ഥ𝒎
𝒁
b
ഥ𝒔
𝒁 ഥ𝒎
𝒁
𝑰𝒄
c
ഥ𝒂 𝑽
𝑽 ഥ𝒃 𝑽
ഥ𝒄 ഥ ′𝒂 𝑽
𝑽 ഥ ′𝒃 ഥ ′𝒄
𝑽
Impedancia de secuencia positiva y
cero
Se puede demostrar que:

𝑍1ҧ = 𝑍𝑆ҧ − 𝑍𝑚
ҧ

𝑍0ҧ = 𝑍𝑆ҧ + 2𝑍𝑚

ҧ
Dónde:
◼ 𝑍1ҧ = Impedancia de secuencia positiva de la línea
◼ 𝑍𝑆ҧ = Impedancia propia de fase
◼ 𝑍𝑚 ҧ = Impedancia mutua
Para el caso de la impedancia de secuencia negativa la
deducción es similar y es conocido que para líneas de
transmisión 𝑍1ҧ = 𝑍2ҧ .
Medición de la Impedancia para una
falla trifásica a tierra
Para una falla trifásica franca a una distancia "𝑚"
desde la ubicación del relé, tenemos:
𝑉ത𝐹𝑎 = 𝑉ത𝐹𝑏 = 𝑉ത𝐹𝑐 = 0

𝑚 ∙ 𝑍𝐿ҧ (1 − 𝑚) ∙ 𝑍𝐿ҧ
𝐼𝑎ҧ
a
𝐼𝑏ҧ
b
𝐼𝑐ҧ
c

𝑉ത𝑎 𝑉ത𝑏 𝑉ത𝑐 𝑉ത𝐹𝑎 𝑉ത𝐹𝑏 𝑉ത𝐹𝑐

Falla Trifásica - Tierra

Medición de la Impedancia para una
falla trifásica a tierra

Luego en el punto del relé, y considerando que:

𝐼𝑏ҧ + 𝐼𝑐ҧ = −𝐼𝑎ҧ .

𝑉ത𝑎 = 𝑚 𝑍𝑆ҧ 𝐼𝑎ҧ + 𝑍𝑀

ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ + 0 = 𝑚 𝑍𝑆ҧ − 𝑍𝑀
ҧ 𝐼𝑎ҧ

𝑉ത𝑏 = 𝑚 𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑆ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ + 0 = 𝑚 𝑍𝑆ҧ − 𝑍𝑀
ҧ 𝐼𝑏ҧ
𝑉ത𝑐 = 𝑚(𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑀
ҧ 𝐼𝑏ҧ + 𝑍𝑆ҧ 𝐼𝑐ҧ ) + 0 = 𝑚(𝑍𝑆ҧ − 𝑍𝑀
ҧ )𝐼𝑐ҧ
Medición de la Impedancia para una
falla trifásica a tierra
Reemplazando 𝑍1ҧ = 𝑍𝑆ҧ − 𝑍𝑚
ҧ , tenemos

𝑉ത𝑎 = 𝑚 𝑍𝑆ҧ − 𝑍𝑀
ҧ 𝐼𝑎ҧ = 𝑚𝑍1𝐿
ҧ 𝐼𝑎ҧ

𝑉ത𝑏 = 𝑚(𝑍𝑆ҧ − 𝑍𝑀
ҧ )𝐼𝑏ҧ = 𝑚𝑍1𝐿
ҧ 𝐼𝑏ҧ

𝑉ത𝑐 = 𝑚(𝑍𝑆ҧ − 𝑍𝑀
ҧ )𝐼𝑐ҧ = 𝑚𝑍1𝐿
ҧ 𝐼𝑐ҧ
Donde:
◼ 𝑍1𝐿ҧ = Impedancia de secuencia positiva de la línea
◼ 𝑍𝑆ҧ = Impedancia propia de fase
◼ 𝑍𝑚 ҧ = Impedancia mutua
◼ 𝑉ത𝑎 , 𝑉ത𝑏 , 𝑉ത𝑐 = Voltajes Fase-Neutro en el punto de ubicación del relé
◼ 𝐼𝑎ҧ , 𝐼𝑏ҧ , 𝐼𝑐ҧ = Corrientes de fase en el punto de ubicación del relé
◼ 𝑚 = Distancia a la falla en “por unidad” (pu)
Medición de la Impedancia para una
falla trifásica a tierra
La impedancia medida por el relé será:
𝑉ത𝑎 𝑉ത𝑏 𝑉ത𝑐
ҧ ;
= 𝑚𝑍1𝐿 ҧ ;
= 𝑚𝑍1𝐿 ҧ
= 𝑚𝑍1𝐿
𝐼𝑎ҧ 𝐼𝑏ҧ 𝐼𝑐ҧ

El desarrollo algebraico muestra que un único

elemento de distancia es suficiente para detectar fallas
trifásicas equilibradas.
Para una línea simétricamente transpuesta, es decir
las impedancias propias de cada fase iguales, y las
impedancias mutuas entre fases también iguales.
Medición de la Impedancias para una
falla bifásica
Para este tipo de falla, la condición de contorno
relevante es que las tensiones de la fase A y B son
iguales.
𝑉ത𝐹𝑏 = 𝑉ത𝐹𝑐

𝑚 ∙ 𝑍𝐿ҧ (1 − 𝑚) ∙ 𝑍𝐿ҧ
𝐼𝑎ҧ
a
𝐼𝑏ҧ
b
𝐼𝑐ҧ
c

𝑉ത𝑎 𝑉ത𝑏 𝑉ത𝑐 𝑉ത𝐹𝑏 𝑉ത𝐹𝑐 Falla Bifásica B-C

Medición de la Impedancias para una
falla bifásica
Sumando las caídas de tensión para las tres fases.

𝑉ത𝑎 = 𝑚 𝑍𝑆ҧ 𝐼𝑎ҧ + 𝑍𝑀

ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ + 𝑉ത𝐹𝑎

𝑉ത𝑏 = 𝑚 𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑆ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ + 𝑉ത𝐹𝑏
Restando
𝑉ത𝑐 = 𝑚(𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑀
ҧ 𝐼𝑏ҧ + 𝑍𝑆ҧ 𝐼𝑐ҧ ) + 𝑉ത𝐹𝑐
Medición de la Impedancias para una
falla bifásica
Haciendo lleva a.
𝑉ത𝑏 − 𝑉ത𝑐 = 𝑚 𝑍𝑆ҧ − 𝑍𝑀
ҧ (𝐼𝑏ҧ − 𝐼𝑐ҧ )

𝑉ത𝑏 − 𝑉ത𝑐 = 𝑚𝑍1𝐿

ҧ (𝐼𝑏ҧ − 𝐼𝑐ҧ )
𝑉ത𝑏 − 𝑉ത𝑐
ҧ
= 𝑚𝑍1𝐿
𝐼𝑏ҧ − 𝐼𝑐ҧ

Un desarrollo similar puede ser efectuado para fallas

A-B y C-A.
Como conclusión, para medir adecuadamente la
distancia a la falla para una falla bifásica B-C se debe
hacer el cociente entre la tensión Fase-Fase con la
diferencia de corrientes de las fases en falla.
Medición de la Impedancias para una
falla monofásica a tierra
Asumiendo las siguiente condición de contorno.

𝑉ത𝐹𝑎 = 0 (𝑠𝑖 𝑅𝐹 = 0)

𝑚 ∙ 𝑍𝐿ҧ (1 − 𝑚) ∙ 𝑍𝐿ҧ
𝐼𝑎ҧ
a
𝐼𝑏ҧ
b
𝐼𝑐ҧ
c

Falla Monofásica
𝑉ത𝑎 𝑉ത𝑏 𝑉ത𝑐 𝑉ത𝐹𝑎
a Tierra A-T
Medición de la Impedancias para una
falla monofásica a tierra
De similar forma que las anteriores desarrollamos las
ecuaciones.

𝑉ത𝑎 = 𝑚(𝑍𝑆ҧ 𝐼𝑎ҧ + 𝑍𝑀

ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ ) + 0

𝑉ത𝑏 = 𝑚(𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑆ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ ) + 𝑉ത𝐹𝑏

𝑉ത𝑐 = 𝑚(𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑀
ҧ 𝐼𝑏ҧ + 𝑍𝑆ҧ 𝐼𝑐ҧ ) + 𝑉ത𝐹𝑐
Medición de la Impedancias para una
falla monofásica a tierra

Realizamos las transformaciones algebraicas para

buscar la impedancia medida.
Agregamos

𝑉ത𝑎 = 𝑚(𝑍𝑆ҧ 𝐼𝑎ҧ − 𝑍𝑀

ҧ 𝐼𝑎ҧ + 𝑍𝑀
ҧ 𝐼𝑎ҧ + 𝑍𝑀
ҧ 𝐼𝑏ҧ + 𝑍𝑀
ҧ 𝐼𝑐ҧ )

= 𝑚((𝑍𝑆ҧ −𝑍𝑀
ҧ )𝐼𝑎ҧ + 𝑍𝑀
ҧ (𝐼𝑎ҧ + 𝐼𝑏ҧ + 𝐼𝑐ҧ ))

ҧ − 𝑍1𝐿
𝑍0𝐿 ҧ ҧ − 𝑍1𝐿
𝑍0𝐿 ҧ
ҧ 𝐼𝑎ҧ +
= 𝑚 𝑍1𝐿 𝐼𝑟ҧ = 𝑚𝑍1𝐿
ҧ 𝐼𝑎ҧ + 𝐼𝑟ҧ
3 ҧ
3𝑍1𝐿

ҧ 𝐼𝑎ҧ + 𝑘ത0 𝐼𝑟ҧ ;

𝑉ത𝑎 = 𝑚𝑍1𝐿
𝑍0ҧ = 𝑍𝑆ҧ + 2𝑍𝑚
ҧ

𝑍1ҧ = 𝑍𝑆ҧ − 𝑍𝑚
ҧ
Medición de la Impedancias para una
falla monofásica a tierra
Donde: ҧ − 𝑍1𝐿
𝑍0𝐿 ҧ
𝑘ത0 =
ҧ
3𝑍1𝐿
Finalmente
𝑉ത𝑎
ҧ
= 𝑚𝑍1𝐿
𝐼𝑎ҧ + 𝑘ത0 𝐼𝑟ҧ
De acuerdo a la última ecuación, para medir
adecuadamente la distancia a la falla, es necesario
dividir la tensión de la fase fallada sobre una corriente
compensada que resulta de la suma de la corriente de
fase de la fase en falla más la corriente residual
multiplicada por el factor de compensación de secuencia
cero 𝑘0 .
Medición de la Impedancias para fallas
sólidas ത
𝑘 =
𝑍ҧ − 𝑍ҧ
ҧ 0
0𝐿 1𝐿
3𝑍1𝐿
Tabla 2.1 Tensiones y corrientes para medición de impedancias de falla, fallas sólidas

Elementos de Distancia ഥ
𝑽 ത𝑰 ഥ 𝒎𝒆𝒅𝒊𝒅𝒐
𝒁
𝑉ത𝑎 − 𝑉ത𝑏
AB 𝑉ത𝑎 − 𝑉ത𝑏 𝐼𝑎ҧ − 𝐼𝑏ҧ ҧ
= 𝑚𝑍1𝐿
ҧ𝐼𝑎 − 𝐼𝑏ҧ

Elementos 𝑉ത𝑏 − 𝑉ത𝑐

BC 𝑉ത𝑏 − 𝑉ത𝑐 𝐼𝑏ҧ − 𝐼𝑐ҧ ҧ
= 𝑚𝑍1𝐿
de Fase 𝐼𝑏ҧ − 𝐼𝑐ҧ

𝑉ത𝑐 − 𝑉ത𝑎
CA 𝑉ത𝑐 − 𝑉ത𝑎 𝐼𝑐ҧ − 𝐼𝑎ҧ ҧ
= 𝑚𝑍1𝐿
𝐼𝑐ҧ − 𝐼𝑎ҧ

𝑉ത𝑎
A 𝑉ത𝑎 𝐼𝑎ҧ + 𝑘ത 0 𝐼𝑟ҧ ҧ
= 𝑚𝑍1𝐿
𝐼𝑎ҧ + 𝑘ത 0 𝐼𝑟ҧ

Elementos 𝑉ത𝑏
B 𝑉ത𝑏 𝐼𝑏ҧ + 𝑘ത 0 𝐼𝑟ҧ ҧ
= 𝑚𝑍1𝐿
de Tierra 𝐼𝑏ҧ + 𝑘ത 0 𝐼𝑟ҧ

𝑉ത𝑐
C 𝑉ത𝑐 𝐼𝑐ҧ + 𝑘ത 0 𝐼𝑟ҧ ҧ
= 𝑚𝑍1𝐿
𝐼𝑐ҧ + 𝑘ത 0 𝐼𝑟ҧ
Medición de la impedancia aparente
con resistencia de falla
La Figura 2.14 muestra un sistema de doble fuente
con una falla en el extremo de la línea (solo por
simplicidad).
Para este sistema veremos cómo resultan las
ecuaciones de la medición de la impedancia aparente
tomando en cuenta la resistencia de falla.
𝐸∠𝛿 Barra S Barra R 𝐸∠0
𝑍𝑆 𝑍𝐿 𝑍𝑅

𝑅𝐹
21

Figura 2.14 Unifilar de un Sistema con doble fuente

Fuente: Elaboración propia
Impedancia medida con resistencia de
falla
Para el análisis de líneas con doble fuente es muy útil
el uso de los factores de distribución de corriente, que
para el caso del circuito de la anterior diapositiva se
puede representar de forma esquemática de la forma:

𝑅𝐹
′
𝑍𝑆 𝑍 𝑍𝑅

𝐼 𝐼𝑇 𝐼𝑅
Impedancia medida con resistencia de
falla
Aplicando la conocida regla de divisor de corriente se
puede encontrar que en general
𝑍𝑅
𝐼= ′
𝐼𝑇 = 𝐶𝐼𝑇
𝑍𝑆 + 𝑍 + 𝑍𝑅
Donde
𝑍𝑅
𝐶=
𝑍𝑆 + 𝑍 ′ + 𝑍𝑅
𝐶 se denomina el factor de distribución de corriente.
Podemos notar que la resistencia de falla no interviene
en los factores de distribución de corriente, por lo que
estos factores no dependen del hecho si el
cortocircuito es sólido o resistivo.
Falla monofásica con resistencia de falla
𝑍𝑅
𝐶=
Tabla 2.2 Impedancia aparente medida para fallas con resistencia de falla 𝑍𝑆 + 𝑍𝐿 + 𝑍𝑅

Tipo Falla
Bucle en ഥ
𝑽 𝑰ത ഥ 𝒎𝒆𝒅𝒊𝒅𝒐
𝒁
Falla
ഥ𝑎
A 𝑉ത𝑎 𝐼𝑎ҧ 𝑉
= 𝑍1′ +
𝑅𝐹
𝐼𝑎ҧ 𝐶1
Falla ഥ𝑏
B 𝑉ത𝑏 𝐼𝑏ҧ 𝑉
= 𝑍1′ +
𝑅𝐹
Trifásica 𝐼𝑏ҧ 𝐶1
ഥ𝑐
C 𝑉ത𝑐 𝐼𝑐ҧ 𝑉
ҧ𝐼𝑐
= 𝑍1′ +
𝑅𝐹
𝐶1
ഥ𝑎 − 𝑉ഥ𝑏
AB 𝑉ത𝑎 − 𝑉ത𝑏 𝐼𝑎ҧ − 𝐼𝑏ҧ 𝑉
ҧ𝐼𝑎 − 𝐼𝑏ҧ
= 𝑍1′ +
𝑅𝐹
𝐶1 + 𝐶2
Falla ഥ𝑏 − 𝑉ഥ𝑐
BC 𝑉ത𝑏 − 𝑉ത𝑐 𝐼𝑏ҧ − 𝐼𝑐ҧ 𝑉
= 𝑍1′ +
𝑅𝐹
Bifásica 𝐼𝑏ҧ − 𝐼𝑐ҧ 𝐶1 + 𝐶2
ഥ𝑐 − 𝑉ഥ𝑎
CA 𝑉ത𝑐 − 𝑉ത𝑎 𝐼𝑐ҧ − 𝐼𝑎ҧ 𝑉
= 𝑍1′ +
𝑅𝐹
𝐼𝑐ҧ − 𝐼𝑎ҧ 𝐶1 + 𝐶2
Vഥa 3R F
A ഥa
V Iaҧ + kത 0 Irҧ ҧIa + kത 0 Iҧr
= Z1′ +
C1 + C2 + (1 + 3k 0 )C0

Falla 𝑉ഥ𝑏 3R F
= Z1′ +
Monofásica B 𝑉ത𝑏 𝐼𝑏ҧ + 𝑘ത 0 𝐼𝑟ҧ ҧ𝐼𝑏 + 𝑘ത 0 𝐼𝑟ҧ C1 + C2 + (1 + 3k 0 )C0

a tierra
ഥ𝑐
𝑉 3R F
= Z1′ +
C 𝑉ത𝑐 𝐼𝑐ҧ + 𝑘ത 0 𝐼𝑟ҧ 𝐼𝑐ҧ + 𝑘ത 0 𝐼𝑟ҧ C1 + C2 + (1 + 3k 0 )C0
Comparadores del relé de distancia
Los ángulos limitadores del comparador de fases son
conocidos como ángulos coincidentes, y se obtiene
una característica frontera simétrica si 𝛽1 = 𝛽2 .
La Figura muestra una representación del comparador
de amplitud y fase.
𝑆0 𝑆1
Comparador Comparador
𝑆𝑅 de Amplitud Disparo 𝑆2 de Fase Disparo

Operación de Operación de
disparo cuando disparo cuando
𝑆0 ≥ 𝑆𝑅 −𝛽1 ≤ 𝛽 ≤ 𝛽2
𝛽 = 𝐴𝑟𝑔(𝑆1 Τ𝑆2 )
Comparador de fase
𝑆1ҧ
Formación
de la señal
𝐴ҧ
Comparador Disparo
𝐵ത de Fase
Formación
𝑆2ҧ
de la señal

𝑆1ҧ = 𝑘1 𝐴ҧ + 𝑘2 𝐵ത

𝑆2ҧ = 𝑘3 𝐴ҧ + 𝑘4 𝐵ത
Donde:
◼ 𝑆1ҧ y 𝑆2ҧ son las señales construidas para comparar.

◼ 𝑘1 , 𝑘2 , 𝑘3 y 𝑘4 constantes definidas escalares.

Comparador coseno
Tomando
𝑆1ҧ = 𝑆1 𝑒 𝑗𝜃1
𝑆2ҧ = 𝑆2 𝑒 𝑗𝜃2
Se usa para comparar si el ángulo de dos señales
fasoriales 𝑆1ҧ y 𝑆2ҧ , están comprendidas en el rango
entre [−90°, +90°].
Según esta restricción, si el ángulo está entre −90°y
+ 90° la función debe dar algún permiso.
Una función matemática que más se adecua para
esto, es la función COSENO.
La función coseno es positiva desde −90°y +90° .
Comparador Coseno
Esto será satisfecho siempre y cuando
−𝜋Τ2 ≤ 𝜃1 − 𝜃2 ≤ + 𝜋Τ2 (2.95)

Siendo 𝜃1 y 𝜃2 los ángulos de fase de las cantidades

construidas 𝑆1 y 𝑆2 .
Es posible demostrar que la característica de este
comparador también es un círculo en un plano alfa o
beta
Comparador de fase
Entonces con
𝑆1ҧ = 𝑆1 𝑒 𝑗𝜃1
(2.92)
𝑆2ҧ = 𝑆2 𝑒 𝑗𝜃2

La condición de operación ser resuelve con:

𝑅𝑒 𝑆1ҧ 𝑆2ҧ ∗ = 𝑅𝑒 𝑆1 𝑆2 𝑒 𝑗(𝜃1 −𝜃2 ≥ 0 (2.93)

Esta expresión debe contener 𝑆1 y 𝑆2 diferentes de

cero o
𝑅𝑒 𝑒 𝑗(𝜃1 −𝜃2 ≥ 0 (2.94)

𝐶𝑜𝑠 𝜃1 − 𝜃1 ≥ 0
Comparador coseno
NO 90º
OPERA OPERA
𝑆1ҧ
Cos(φ) 𝜑 < 90°
CONDICION
DE 180º
REF
OPERACION: NO 𝑆2ҧ
OPERA OPERA
- 90º
90º
NO
-90º 90º φ CONDICION OPERA OPERA
DE 𝑆2ҧ
RESTRICCION: 180º REF
𝜑 > −90°

OPERA
NO
𝑆1ҧ OPERA
- 90º

Se usa para comparar si el ángulo de dos señales

fasoriales 𝑆1ҧ y 𝑆2ҧ , están comprendidas en el rango
entre [−90°, +90°].
Según esta restricción, si el ángulo está entre −90°y
+ 90° la función debe dar algún permiso.
Comparador coseno
Para la implementación digital de este comparador se
hace uso de la siguiente relación matemática.
Siendo 𝑆1ҧ y 𝑆2ҧ dos fasores y 𝑆2ҧ ∗ el conjugado del fasor
𝑆2ҧ .
𝑆1ҧ ∙ 𝑆2ҧ ∗ = 𝑆1 𝑆2 𝑐𝑜𝑠𝜃 + 𝑗 𝑆1 𝑆2 𝑠𝑒𝑛𝜃 (2.106)

Donde: 𝑆1ҧ = 𝑆1 ∠𝜃1 ; 𝑆2ҧ = 𝑆2 ∠𝜃2 ; 𝜃 = 𝜃1 − 𝜃2

La componente real de (2.106) tiene el máximo valor
cuando 𝑆1ҧ y 𝑆2ҧ están en fase y cero cuando 𝑆1ҧ y 𝑆2ҧ
están en cuadratura. Comparando este valor si es
positivo o no, obtenemos el comparador de fase tipo
coseno.
𝑅𝑒 𝑆1ҧ ∙ 𝑆2ҧ ∗ = 𝑆1 𝑆2ҧ 𝑐𝑜𝑠𝜃
Comparador coseno
Y para obtener el comparador coseno solo hacemos

𝑅𝑒 𝑆1ҧ ∙ 𝑆2ҧ ∗ > 0 (2.108)

Como referencia los comparadores coseno analógicos

eran:
Relés Electromecánicos
◼ Bobina móvil.
◼ Disco y de cilindro de inducción (Principio de Ferraris
modificado por medio de desfasamiento).

Relés Estáticos
◼ Puente rectificador de onda completa con circuito integrador.
◼ Circuito lógico combinacional con circuito integrador.
Comparadores del relé mho
𝑍ҧ𝑅
Punto de Alcance
𝑍ҧ𝑚𝑒𝑑 del Relé (Punto de
Balance)
a)

C
A X
B
𝑍𝑅ҧ − 𝑍𝑚𝑒𝑑
ҧ
B

𝑍𝑅ҧ D ത 1
𝑆"
90°
b)
o C
ҧ
𝑍𝑚𝑒𝑑 +90º
ത 2
𝑆"
c)

R
A

Figura 2.27 Principio de la protección de distancia, medición de la impedancia de falla

Fuente: Elaboración propia
Comparadores del relé mho
𝑆"ҧ 1 = 𝑍𝑅ҧ − 𝑍𝑚𝑒𝑑
ҧ

𝑆"ҧ 2 = 𝑍𝑚𝑒𝑑
ҧ

Multiplicando ambas señales por la corriente del bucle

de falla 𝐼𝐿ҧ .
𝑆′ҧ 1 = 𝐼𝐿ҧ ∙ 𝑍𝑅ҧ − 𝐼𝐿ҧ ∙ 𝑍𝑚𝑒𝑑
ҧ

𝑆′ҧ 2 = 𝐼𝐿ҧ ∙ 𝑍𝑚𝑒𝑑

ҧ

Reemplazando ILҧ ∙ 𝑍𝑚𝑒𝑑

ҧ por la tensión del bucle de falla
𝑉ത𝐿 tenemos las señales efectivamente comparadas.
𝑆1ҧ = 𝐼𝐿ҧ ∙ 𝑍𝑅ҧ − 𝑉ത𝐿 ∗
𝑅𝑒 𝑆1ҧ ∙ 𝑆2ҧ > 0
𝑆2ҧ = 𝑉ത𝐿
El modelo de relé genérico
El modelo propuesto
Toma de Tratamiento
Simulación ATP Filtro digital muestras por Estimación
Anti-Aliasing de fasores algoritmos Salida a ATP
ciclo

Estimador Algoritmos de
Muestreo Salida a ATP
de fasores protección
TCs, TPs

El modelo inicia con las señales de TCs y TPs modelados

en ATP como parte del modelo del sistema de potencia.
Estas señales son de naturaleza digital, debido a que las
señales de corriente y tensión son ondas simuladas de
acuerdo al paso de integración de la simulación de
EMTP.
El efecto aliasing
t
n= 1 2 3 4 5 6 7 8 9
0
Frecuencia

10
𝑓
9 𝑠
1
𝑓
9 𝑠

Frecuencia

1
𝑓
3 𝑠

∆𝑡 = intérvalo de muestreo
1
∆𝑡
= 𝑓𝑠 = frecuencia de muestreo
Puntos de
muestreo

Figura 3.7 Efecto aliasing en muestreo de señales.

Fuente: IEE [13]
Teorema del muestreo
El teorema del muestreo establece que una señal
analógica, conteniendo componentes con una
máxima frecuencia de 𝑓𝑐 , puede ser completamente
representada por muestras siempre que la frecuencia
de muestreo 𝑓𝑠 sea al menos dos veces 𝑓𝑐 .
Es decir, la frecuencia de muestreo debe ser al
menos el doble de la frecuencia de la señal de
interés.
𝑓𝑠 ≥ 2𝑓𝑐
Filtro Butterworth
La respuesta en amplitud 𝐻(𝜔) de un filtro pasa
bajo Butterworth de orden N está dado por
1
𝐻(𝜔) =
𝜔 2𝑁
1+( )
𝜔𝑐

Donde:
N es el orden del filtro
𝜔0 es la frecuencia de corte
Filtro Butterworth
y la función normalizada es
1
෡
𝐻(𝑠) = (3.20)
𝑠 − 𝑠1 𝑠 − 𝑠2 … 𝑠 − 𝑠𝑛

La función de transferencia final es obtenida

reemplazando
𝑠
෡
𝐻(𝑠) = 𝐻( )
𝜔𝑐

(3.21)
3.2.1 Filtro Butterworth
En la Figura 3.3 se muestra la respuesta en
frecuencia de un filtro Butterworth de 1er a 10mo
orden.

Figura 3.3 Respuesta filtro Butterworth.

Fuente: Elaboración propia en MATLAB
Filtro Butterworth
En el modelo generalizado se decidió aplicar un filtro digital
Butterworth de tercer orden con una frecuencia de corte de
180 Hz, que emula al filtro analógico; diseñado para obtener
correctamente la frecuencia fundamental para propósitos de
la protección de distancia.

1
෡ 𝑠 =
𝐻
𝑠 3 + 2𝑠 2 + 2𝑠 + 1

Para un filtro con una frecuencia de corte de 180 Hz

tenemos que:
𝑟𝑎𝑑
𝜔𝑐 = 2𝜋𝑓𝑐 = 2 ∗ 𝜋 ∗ 180 = 1130.9734 [ ]
𝑠
Filtro Butterworth
Reemplazando 𝑠 por 𝑠/𝜔𝑐 en la ecuación del filtro
normalizado, y reordenando, se tiene la función de
transferencia de un filtro analógico Butterworth de
tercer orden

1446628844.7961
𝐻 𝑠 = 3 (3.25)
𝑠 + 2261.9467𝑠 2 + 2558201.4608𝑠 + 1446628844.7961
Filtrado Digital
Filtros de respuesta de Impulso Finito
(FIR)
Un filtro FIR es aquel que tiene una respuesta finita
al impulso, y se caracterizan por ser sistemas no
recursivos, es decir se basan en obtener la salida a
partir exclusivamente de las entradas actuales y
anteriores.
Así para un filtro de longitud N, tendremos:
𝑁−1

𝑦 𝑛 = 𝑏0 𝑥 𝑛 + 𝑏1 𝑥 𝑛 − 1 + ⋯ + 𝑏𝑁−1 𝑥 𝑛 − 𝑁 + 1 = ෍ 𝑏𝑘 𝑥 𝑛 − 𝑘 (3.61)
𝑘=0

Donde 𝑏𝑘 , son los coeficientes del filtro.

Filtros de respuesta de Impulso Finito
(FIR)
La salida puede escribirse como la convolución de la
entrada 𝑥[𝑛] con la respuesta impulsional ℎ[𝑘]:
+∞ 0, 𝑘<0
𝑦 𝑛 = ෍ ℎ[𝑘] ∙ 𝑥 𝑛 − 𝑘 𝑐𝑜𝑛 ℎ 𝑘 = ቐℎ 𝑘 , 0≤𝑘<𝑁 (3.62)
𝑘=−∞ 0, 𝑘≥𝑁

Luego, la expresión puede reescribirse como:

𝑁−1

𝑦 𝑛 = ෍ ℎ[𝑘] ∙ 𝑥 𝑛 − 𝑘 (3.63)
𝑘=0

Donde identificando (3.61) y (3.62 𝑏𝑘 = ℎ[𝑘], es

decir los coeficientes del filtro es la respuesta al
impulso [22].
Filtros de respuesta de Impulso Infinito
(IIR)
En este tipo de filtros la salida depende además de
salidas anteriores y que, estando en reposo, al ser
estimulados con una entrada de impulso, su salida no
vuelve al reposo, de ahí el nombre de filtros de
respuesta de impulso infinita (IIR).
Su ecuación en diferencias generales es de la forma,
[22]:
𝑦 𝑛 = 𝑏0 𝑥 𝑛 + 𝑏1 𝑥 𝑛 − 1 + ⋯ + 𝑏𝑀 𝑥 𝑛 − 𝑀 − 𝑎1 𝑦 𝑛 − 1 − 𝑎2 𝑦 𝑛 − 2 − ⋯ − 𝑎𝑁 𝑦 𝑛 − 𝑁
𝑀 𝑁

= ෍ 𝑏𝑘 ∙ 𝑥 𝑛 − 𝑘 − ෍ 𝑎𝑘 ∙ 𝑦 𝑛 − 𝑘 (3.64)
𝑘=0 𝑘=1

Donde el orden es el mayor de M y N.

Filtros de respuesta de Impulso Infinito
(IIR)
La función de transferencia en Z (transformada Z) de
este tipo de filtro es:[22]:
𝑀
෌𝑘=0 𝑏𝑘 ∙ 𝑧 −𝑘
𝐻𝑧 = (3.65)
σ𝑁
𝑘=0 𝑎𝑘 ∙ 𝑧 −𝑘

Existen dos filosofías de diseños de filtros IIR.

Filtros de respuesta de Impulso Infinito
(IIR)
◼ INDIRECTA. Se basa en aplicar filtros análogos diseñados
previamente, transformaciones que los conviertan en
digitales con las mismas características. Para esto hay tres
métodos fundamentales:
 Diseño por impulso invariante
 Diseño por analogía y aproximación de derivadas y
 Diseño por transformación bilineal

◼ DIRECTA. Se propone el diseño de filtros digitales

imponiendo una serie de condiciones a la respuesta para
determinar los coeficientes. Hay dos métodos simples que
son:
 Diseño por la aproximación de Padé
 Diseño por la aproximación de mínimos cuadrados
Digitalización del Filtro
En el modelo propuesto, se utilizó el método de la transformada
bilineal como procedimiento para digitalizar el filtro analógico
Butterworth de 3er orden normalizado
El apéndice D muestra los cálculos realizados para encontrar los
coeficientes del filtro digital con los siguientes datos:

fc≔180 Hz frecuencia de corte

fs≔28800 frecuencia de muestreo [Hz]
Wc≔2∙π∙180=Wc=1130.9734 [rad/seg]

ts≔1/fs=0.00003472 [seg] Período de muestreo corte

𝑏0∙𝑧+𝑏1∙𝑧 −1 +𝑏2∙𝑧 −2 +𝑏3∙𝑧 −3

𝐻 𝑧 =
1+𝑎1∙𝑧 −1 +𝑎2∙𝑧 −2 +𝑎3∙𝑧 −3
Digitalización del Filtro
Los coeficientes calculados son:

𝑏0 = 0.000007281172475 𝑎0 = 1
𝑏1 = 0.000021843517425 𝑎1 = − 2.92146522319489
𝑏2 = 0.000021843517425 𝑎2 = 2.84598405576111
𝑏3 = 0.000007281172475 𝑎3 = −0.924460583186421

Finalmente, por inspección podemos expresar la ecuación en

dominio del tiempo del filtro IIR en la forma general.

𝑦 𝑛 =𝑏 0 𝑥 𝑛 +𝑏 1 𝑥 𝑛−1 +𝑏 2 𝑥 𝑛−2 +𝑏 3 𝑥 𝑛−3

+𝑎 1 𝑦 𝑛 − 1 + 𝑎 2 𝑦 𝑛 − 2 + 𝑎 3 𝑦(𝑛 − 3)
Modelo de la etapa de muestreo
Para el modelo del relé se propuso utilizar un
muestreo de sus señales analógicas de 24 muestras
por ciclo (𝑓𝑠 = 1200 𝐻𝑧, a una frecuencia de sistema
50 Hz.), y se ha elegido una frecuencia de simulación
de EMT de 28800 𝐻𝑧 (múltiplo de 24 de la frecuencia
de muestreo del relé).
El paso de integración ajustado para la simulación en
EMT, debe ser el mismo valor que se utiliza para el
cálculo de los coeficientes del filtro digital.
Las muestras por ciclo se puede calcular de
𝑇 𝑓𝑚
𝑚𝑝𝑐 = 𝑚𝑝𝑐 =
∆𝑡 𝑓𝑛

Dando para 28800 Hz 576 mpc.

Algoritmos de estimación de fasores

pu
Ciclo completo
Medio Ciclo
1

−1

Figura 3.14 Ventanas de datos de medición de los algoritmos.

Fuente: Reproducción de [23].
Algoritmo Fourier de ciclo completo
De acuerdo con la teoría de análisis de Fourier,
cualquier función periódica 𝑓(𝑡) puede ser
representada por las siguientes series de Fourier de
armónicos discretos [32]:
∞ ∞
𝑎0
𝑓 𝑡 = + ෍ 𝑎𝑛 𝑐𝑜𝑠 𝑛𝜔0 𝑡 + ෍ 𝑏𝑛 𝑠𝑒𝑛(𝑛𝜔0 𝑡) (3.102)
2
𝑛=1 𝑛=1
Donde:
◼ 𝜔0 es la frecuencia angular fundamental = 2𝜋𝑓0 = 2𝜋Τ𝑇
◼ 𝑇 es el periodo de la frecuencia fundamental
◼ 𝑛𝜔0 es la frecuencia angular del enésimo armónico
◼ 𝑡 es la variable independiente del tiempo
Algoritmo Fourier de ciclo completo
Si tomamos por ejemplo una forma de onda de
voltaje 𝑣(𝑡), tenemos
∞ ∞
𝑎0
𝑣 𝑡 = + ෍ 𝑎𝑛 𝑐𝑜𝑠 𝑛𝜔0 𝑡 + ෍ 𝑏𝑛 𝑠𝑒𝑛(𝑛𝜔0 𝑡) (3.106)
2
𝑛=1 𝑛=1

y
𝑡0 +𝑇
2 (3.107)
𝑎𝑛 = න 𝑣 𝑡 𝑐𝑜𝑠(𝑛𝜔0 𝑡) 𝑑𝑡 , 𝑛 = 0, 1, 2, . . .
𝑇
𝑡0

𝑡0 +𝑇
2
𝑏𝑛 = න 𝑣 𝑡 𝑠𝑒𝑛(𝑛𝜔0 𝑡) 𝑑𝑡 , 𝑛 = 1, 2, . . . (3.108)
𝑇
𝑡0
Algoritmo Fourier de ciclo completo
Si el tiempo bajo consideración es 𝑡0 , entonces 𝑉𝑥
puede ser encontrado de la ecuación (3.107):
𝑡0 +𝑇
2
𝑉𝑥 = 𝑎1 = න 𝑣 𝑡 𝑐𝑜𝑠(𝜔0 𝑡) 𝑑𝑡 (3.109)
𝑇
𝑡0

Supóngase que 𝑁 el número de muestras por ciclo de

la componente fundamental, ∆𝑡 el intervalo de
muestreo, 𝑡𝑖 = 𝑖∆𝑡 el tiempo de la i-ésima muestra y
𝑇 = 𝑁∆𝑡 el periodo del componente fundamental.
Entonces, la integral (3.109) puede ser evaluada
usando el método rectangular de integración
Algoritmo Fourier de ciclo completo
Dando
2
𝑉𝑥 ≈ [𝑣 𝑡0 𝑐𝑜𝑠 𝜔0 𝑡0 + 𝑣 𝑡1 𝑐𝑜𝑠 𝜔0 𝑡1 + ⋯ 𝑣 𝑡𝑖 𝑐𝑜𝑠 𝜔0 𝑡𝑖 + ⋯
𝑁∆𝑡
⋯ + 𝑣 𝑡𝑁−1 𝑐𝑜𝑠 𝜔0 𝑡𝑁−1 + 𝑣 𝑡𝑁 𝑐𝑜𝑠 𝜔0 𝑡𝑁 ]∆𝑡

𝑁
2 2𝜋𝑖
𝑉𝑥 ≈ ෍ 𝑣𝑖 𝑐𝑜𝑠 (3.110)
𝑁 𝑁
𝑖=0

donde, 𝑣𝑖 = 𝑣(𝑡𝑖 ) es la i-ésima muestra de la forma

de onda de voltaje. Similarmente, para 𝑉𝑦 tenemos:

𝑁
2 2𝜋𝑖
𝑉𝑦 ≈ ෍ 𝑣𝑖 𝑠𝑒𝑛 (3.113)
𝑁 𝑁
𝑖=0
Algoritmo Fourier de ciclo completo
La Figura 3.18 muestra la respuesta en frecuencia del
algoritmo de onda completa, donde los componentes
armónicos son filtrados.
output
input 1.0

0.75

0.5

0.25

0
0 1 2 3 4 5 6
𝑤
𝑤0

Figura 3.18 Respuesta de frecuencia del algoritmo de Fourier de onda completa

Fuente: Johns A. T., Salman S. K. [2]
Implementacion en ATP y
Resultados de Simulación
ATP y ATPDraw
◼ Relación entre ATP/EMTP y ATPdraw

*.ACP Editor de texto

project File ASCII

*.ALC
Line data

USP *.ATP
Library Input File

*.PCH ATP
library
(tpbig.exe)

MATLAB GTPPLOT
MATHCAD
*.PL4
Output File
PLOTXY COMTRADE
*.DAT, *.CFG
OMICRON
Archivo
generado
por ATPDraw
TACS y MODELS de ATP

◼ - TACS – Transient Analysis of Control Circuits

◼ - MODELS – Programmed modeling of control and other circuits
MODELS
MODELS es un lenguaje de descripción técnica general
soportado por el solucionador de simulación de ATP. Es
un lenguaje de programación dirigido a la simulación en
el dominio del tiempo, y se utiliza como herramienta para
desarrollar modelos de componentes de circuito y control
que no se puedan construir fácilmente con el conjunto
de componentes existentes disponibles en ATP y TACS.
También se puede utilizar para generar señales o para
analizar mediciones del circuito. Finalmente, se puede
usar como interfaz entre ATP y programas externos.
Uno de los propósitos de tener MODELS disponibles en
ATP es:
◼ Tener la flexibilidad de un lenguaje de programación completo
sin tener que interactuar con ATP a nivel de programación;
MODELS BEGIN NEW DATA CASE
...
.... datos misceláneos de líneas
MODELS
..... código
MODEL
......código
1 MODEL ............
............
ENDMODEL
USE Grupo MODELS
......código
ENDUSE
RECORD
.... código
ENDMODELS
C ....
....Sigue archivo de datos de EMTP
 MODEL nombre_model
-- comentarios
MODELS DATA --declaraciones
.... entradas
CONST -- declaraciones
.....entradas
HISTORY –- directiva
..... entradas
DELAY CELLS -– directiva
..... entradas
-- comentarios
VAR -- declaraciones
..... entradas
INIT -- procedimiento
..... entradas
ENDINIT
-- comentarios
EXEC –- procedimiento
.... algoritmos de simulación del model
ENDEXEC
ENDMODEL
FOREIGN MODELS
MODELS proporciona una interfaz predefinida que puede
ser utilizada por el usuario para conectar programas
externos "FOREIGN" a una simulación.
La interfaz se define con cuatro matrices que transportan
los valores de datos, entrada, salida y variables entre
MODELS y el procedimiento externo.
Cada modelo externo debe proporcionar tanto un
procedimiento de ejecución como un procedimiento de
inicialización, con una función correspondiente al
procedimiento EXEC y al procedimiento INIT de un
modelo. Por ejemplo, un modelo FORTRAN se escribiría
como:
FOREIGN MODEL EN FORTRAN
construir fácilmente con el conjunto de componentes
existentes disponibles en ATP y TACS.

SUBROUTINE execsub(xdata, xin, xout, xvar)

DIMENSION xdata(*), xin(*), xout(*), xvar(*)
... execution procedure ...
RETURN
ENTRY inisub(xdata, xin, xout, xvar)
... initialization procedure ...
END
Integración del modelo en ATP
◼ Subrutina conteniendo el modelo del relé (encabezamiento)
Compilado de un nuevo tpbig.exe
El compilador utilizado por ATP-EMTP es el G77.exe
de Mingw de GNU, de distribución gratuita.
Se requiere los archivos fuente de ATP que son
distribuidos a usuarios licenciados.
Es necesario el uso de un programa utilitario
denominado make.exe, también de distribución
gratuita de GNU. Sirve para compilar y enlazar varios
archivos y librerías que vienen de forma separada
Para utilizar make.exe es necesario preparar un
archivo makefile donde se especifican los archivos a
compilar, el compilador a utilizar y las opciones de
compilación.
Es de gran ayuda para el desarrollador del modelo
porque ahorra tiempo y comandos de compilación
Archivo
makefile
Compilado de un nuevo tpbig.exe
 MODELS
INPUT va {v(VSECA)}
vb {v(VSECB)}
Invocación al modelo con MODELS vc {v(VSECC)}
ia {i(CTSA)}
ib {i(CTSB)}
ic {i(CTSC)}
OUTPUT vam, vap, vbm, vbp, vcm, vcp
iam, iap, ibm, ibp, icm, icp
adf, rts, zap, trp, pi2, mp2
zan, ires, zag, zaa

MODEL UmssAtpRel FOREIGN UMSS_ATP_REL {ixdata:9, ixin:7, ixout:40, ixvar:181}

USE umssAtpRel AS umssAtpRel

DATA xdata[1..9] := [50.0 28800.0 1200.0 600.0 2000.0 44.95 1.41 0.76 0.17]
INPUT xin[1..7] := [va vb vc ia ib ic t]
HISTORY xvar[1..176] := 0.0
xvar[177] := 1.0
xvar[178..181] := 0.0
xout[37..40] := 0.0
OUTPUT vam := xout[4], vap := xout[5]
zan :=xout[37], ires :=xout[38]
zag := xout[39], zaa :=xout[40]
ENDUSE
RECORD vafil AS vafil
vbfil AS vbfil
vcfil AS vcfil
zaa AS zaa
ENDMODELS
C ------------------------------------------------------------------------------
Archivo *.bat para corer el modelo en
ATP
◼ Este archivo es equivalente a hacer correr ATP por línea de
comando

◼ Como resultado de una simulación en ATP, se generan archivos

gráficos con extensión .pl4, donde se registran las señales
instantáneas de corrientes y tensiones, así como también señales
digitales de estados y otras magnitudes que han sido
programadas para su registro. Este archivo .pl4 se lo puede
visualizar haciendo uso del programa PlotXY.
Red de prueba
Se usó el modelo en ATP de la red reducida a ocho
barras de 230 kV y 115 kV del sistema de transmisión
de Bolivia, operado por ENDE Transmisión.
MAZOCRUZ PALCA

S1 VINTO115 VINTO230 CUMBRE

C1 36 km
193 km 31 km
36 km S2
Compensación
2TZ-1 PAL115 3TZ-2
Equiv. de red Serie Equiv. de red
S3

Equiv. de red

25 km 244 km 244 km

SANTIVAÑEZ

117 km
PAGADOR

S4

Equiv. de red
Aplicación del modelo a una línea
Los TCs t TPs están modelados de forma ideal en
ATPDraw
Estudios de simulación
Las fallas simuladas fueron:
◼ Evaluación de las diferentes etapas internas del modelo:
filtrado analógico, muestreo de las señales analógicas para
digitalización, estimación de fasores y evaluación del
algoritmo de operación del relé ante una falla.
◼ Protección local de una línea de transmisión para fallas
bifásicas y monofásicas a tierra localizadas a 10, 50, 75 y
85% de la longitud de la línea protegida.
◼ Protección de respaldo de una línea de transmisión para
fallas bifásicas y monofásicas a tierra localizadas a 10, 50,
75 y 100% de la longitud de la línea de transmisión remota.
◼ Fallas en zona reversa de la protección de distancia.
◼ Falla bifásica y monofásica a tierra en bornes del banco de
compensación serie.
 Evaluación del procesamiento de señales
El filtro anti-aliasing

a) a)

b)
b)

Figura 5.6 Tensiones y corrientes de falla, Figura 5.7 Señal de tensión Vb en

valores secundarios, a) Señales de voltaje, diferentes etapas, a) Señal entregada por el
b) Señales de corriente TP, b) Señal filtrada anti-Aliasing
Evaluación del procesamiento de señales
1.81 ms
El filtro anti-aliasing

Figura 5.8 Retraso de grupo de la

señal de tensión Vb

El desfase medido en
tiempo es de 1.81 ms,
debido al retraso de
grupo

Figura 5.9 Respuesta del retraso de grupo

filtro Butterworth
Evaluación del procesamiento de señales
Etapa de muestreo

Figura 5.10 Señal de tensión Vb Señal muestreada

La señal de simulación de ATP es 28800 Hz, que corresponde a

576 muestras por ciclo. Para uso en el modelo de relé, es
necesario re-muestrear la señal a una frecuencia de 24
muestras por ciclo. Estas 24 muestras por cada ciclo son las
graficadas en la Figura 5.10
Evaluación del procesamiento de señales
Etapa de estimación de fasores

a) a)

25ms 45ms

b) b)

Figura 5.11 Señal de tensión Vb, etapa de Figura 5.12 Señal de corriente Ib, etapa de
Estimación de fasores. a) Magnitud, b) Fase Estimación de fasores. a) Magnitud, b) Fase
 Ajustes unidad de distancia LT Santiváñez-
Pagador

Dato Valor T. Operación Característica Criterio utilizado

Longitud LT 116.71 km - - -
TC_Ratio 600:1 - - -
TP_Ratio 2000:1 - - -
Z1_LT 37.45∠80.95° Ωprim - - -
Z0_LT 123.01∠73.77° Ωprim - - -
Z1_LT 11.24∠80.95° Ωsec - - -
Z0_LT 36.90∠73.77° Ωsec - - -
ZONA 1 8.98 Ωsec 0 Mho 80%Z1_LT
ZONA 2 13.49 Ωsec 420 ms Mho 120%Z1_LT
ZONA 3 16.38 Ωsec 840 ms Mho 120%(Z1_LT+Z_LTady)
Operación Zona 1 de fases
El propósito de este estudio fue estudiar el
desempeño de la protección de distancia y evaluar
los alcances de zona cuando ocurren diferentes fallas
en la línea de transmisión, donde el relé de distancia
proporciona protección local primaria.
La condición de operación del comparador de fases
está dado por

π π
ത ത
− ≤ arg(S1 ΤS2 ) ≤
2 2
Operación Zona 1 de fases
Se simularon fallas fase-fase A-B al 10%, 50%, 75%
y 85% de la longitud de la línea desde el extremo
Santiváñez. La falla introducida a los 0.1 s de la
simulación es una falla bifásica A-B sin resistencia de
falla (R_f=0).
Operación Zona 1 de fases
Figura 5.14 Señal de corriente y tensión, fallas al a) 10% de SAN-PAG,
b) 50% de SAN-PAG, c) 75% de SAN-PAG y d) 85% de SAN-PAG

a)

b)

c)

d)
Operación Zona 1 de fases
Figura 5.15 Diferencia angular entre S1 y S2, fallas al a) 10% de SAN-
PAG, b) 50% de SAN-PAG, c) 75% de SAN-PAG y d) 85% de SAN-PAG

a)
c)

b) d)
Operación Zona 1 de fases
Figura 5.16 Señal de disparo de Z1 fallas al a) 10% de SAN-PAG, b)
50% de SAN-PAG, c) 75% de SAN-PAG y d) 85% de SAN-PAG

a) c)

b) d)
Operación Zona 1 de fases
Figura 5.17 Falla en el plano R-X, falla A-B al a) 10% de SAN-PAG, b)
50% de SAN-PAG, c) 75% de SAN-PAG y d) 85% de SAN-PAG

a) b)

c) d)
Operación Zona 1 de fases
Figura 5.18 Impedancia medida loop A-B al a) 10% de SAN-PAG, b) 50% de
SAN-PAG, c) 75% de SAN-PAG y d) 85% de SAN-PAG

a) b)

c) d)

Los valores coinciden con el punto de falla simulado, ej. al 50% de la longitud de la
línea la impedancia debería ser 11.24/2=5.62 Ohm_ sec, y el modelo arroja un
valor medido en la fig b) de 5.33 Ohms_sec.
Conclusiones
Resumen conclusiones
El modelado de relés para su simulación en el ambiente
de transitorios electromagnéticos puede proporcionar
una importante fuente de información tanto para
fabricantes, planificadores, operadores y estudiantes.
El trabajo reportado en esta tesis ha realizado las
siguientes aportaciones:
◼ El proyecto de investigación ha presentado una
metodología clara y detallada para poder modelar
relés en ATP/EMTP con un modelo de relé genérico.
◼ Se ha desarrollado un relé de distancia digital
numérico genérico integrado a motor solucionador
de ATP.
Resumen conclusiones
En todas estas simulaciones se ha encontrado el
desempeño esperado del modelo, tomando en cuenta
las características de la red.
Estos resultados aportarán información relevante para
entender el desempeño esperado de un relé de
distancia en la red de ENDE transmisión.
Los resultados permitieron examinar las variables
internas del relé, para evaluar su desempeño.
Se mostro una metodología detallada de integración del
modelo del relé dentro el motor de ATP, aplicable a
cualquier otro modelo de relé.
 Gracias…
Preguntas?
Link para descargar las diapositivas