Sistema de Detección DPs Técnicas Microondas PDF

Sistema de detección de descargas
parciales en transformadores de potencia

basado en técnicas de microondas.
David Leonardo Galindo Huertas
Universidad Nacional de Colombia

Facultad de Ingenieria, Eléctrica y Electrónica, Maestría en Automatización Industrial
Bogotá D.C, Colombia
2019
Sistema de detección de descargas
parciales en transformadores de potencia
David Leonardo Galindo Huertas
Tesis presentada como requisito parcial para optar al título de:

Magister en Ingeniería Automatización Industrial
Director:
Ph.D Javier Leonardo Araque Quijano
Codirector:
MS Francisco Amortegui Gil
Línea de Investigación:
Detección y diagnostico de fallos
Grupo de Investigación:
CMUN
Universidad Nacional de Colombia

Facultad de Ingenieria, Eléctrica y Electrónica, Maestría en Automatización Industrial
Bogotá D.C, Colombia
2019
A mi sobrino quien me demuestra día a día que
no existen imposibles
Resumen
Esta investigación presenta los resultados de un sistema de detección de descargas parciales

de manera incipiente, siendo este el mayor problema de fallos en transformadores de potencia
ubicados en la parte activa principalmente en los devanados del transformador. La ruptura del
dieléctrico representa para la maquina eléctrica alto grado de degradación del sistema, para
resolver este problema se propone un sistema de detección basado en técnicas microondas
usando software definido por radio y una tarjeta de adquisición USRP B210, en una primera
fase se desarrolla simulación en el software CST del comportamiento del campo EM en
un transformador de potencia, después se presenta la selección de tecnologías y montajes
usados en el sistema, a continuación un algoritmo es desarrollado para adquirir las señales
de campo E y por ultimo se presentan los resultados del sistema de la señal capturada por
radiofrecuencia y una señal base de 60Hz capturada por la tarjeta de sonido del computador.
Lista de Figuras
2-1. Vector de fallo DPs. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2-2. Número de electrones por segundo impulsados por el campo eléctrico. . . . . 5
2-3. Cantidad de electrones ionizados por probabilidad de ionización . . . . . . . 5
2-4. Interacción de electrones y iones en el dipolo . . . . . . . . . . . . . . . . . . 6
2-5. Proceso de generación de múltiples avalanchas y corriente de la DPs . . . . . 7
2-6. Descargas superficiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2-7. Descarga interna . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2-8. Campo Ea distribuido y campo de la vacuola Ec . . . . . . . . . . . . . . . 9
2-9. Métodos de medida de Descargas parciales . . . . . . . . . . . . . . . . . . . 11
2-10.Circuito equivalente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2-11.Simulación en LTSPICE. Ciclo repetitivo de descargas parciales en una onda
sinusoidal de 60 Hz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2-12.Corriente de DPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2-13.Modelo de dipolo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2-14.Pulso Gaussiano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2-15.Dipolo eléctrico y punto de observación . . . . . . . . . . . . . . . . . . . . . 15
2-16.Campo eléctrico lejano producido por DPs . . . . . . . . . . . . . . . . . . . 16
2-17.Derivada pulso Gaussiano . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2-18.Métodos de detección de fallos según Isermann . . . . . . . . . . . . . . . . . 18
4-1. Diagrama de flujo para simulación . . . . . . . . . . . . . . . . . . . . . . . . 22

4-2. Modelo 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4-3. Ubicación antenas en el modelo de simulación . . . . . . . . . . . . . . . . . 24
4-4. Sonda S3 resultados en el dominio del tiempo . . . . . . . . . . . . . . . . . 25
4-5. Sonda S4 resultados en el dominio del tiempo . . . . . . . . . . . . . . . . . 26
4-7. Sonda S3 resultados en el dominio de la frecuencia . . . . . . . . . . . . . . . 26
4-6. Espectro en frecuencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
4-8. Sonda S4 resultados en el dominio de la frecuencia . . . . . . . . . . . . . . . 27
5-1. Ingreso de antenas al transformador [Siegel et al., 2017] . . . . . . . . . . . . 29

5-2. Antenas modelo CAD en CST, a) modelo monopolo, b)modelo gota, c) Acu-
mulación de cargas en el monopolo y d) Acumulación de cargas en la gota . . 30
5-3. Antenas utilizadas para la detección de DPs de izquierda a derecha: Antena
QUAD GSM, Antena Wifi, Antena JCG015 . . . . . . . . . . . . . . . . . . 31
x Lista de Figuras
5-4. USRP B210 y Arquitectura USRP B210 . . . . . . . . . . . . . . . . . . . . 32

5-5. Circuito adquisición y tarjeta de sonido . . . . . . . . . . . . . . . . . . . . 32
5-6. Tecnologías usadas en el proceso de análisis de DPs . . . . . . . . . . . . . . 34
5-7. Cavidad Resonante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
5-8. Montaje de laboratorio Alta tensión . . . . . . . . . . . . . . . . . . . . . . . 35
5-9. Clasificación de descargas parciales basado en electrodos. a) DPs internas, b)
DPs superficial, c) DPs corona [Kreuger et al., 1993] . . . . . . . . . . . . . . 36
5-10.Electrodos para simular descargas tipo corona, superficiales e internas . . . . 36
5-11.Tecnología utilizada para la detección de DPs . . . . . . . . . . . . . . . . . 37
6-1. Diagrama de flujo Visión general del procesamiento de datos . . . . . . . . . 39

6-2. Diagrama de flujo Gnu-Radio . . . . . . . . . . . . . . . . . . . . . . . . . . 40
6-3. Post procesamiento de las señales . . . . . . . . . . . . . . . . . . . . . . . . 41
6-4. Modelo para descomposición de Wavelet . . . . . . . . . . . . . . . . . . . . 42
6-5. Supresión de señal con técnica de Wavelet . . . . . . . . . . . . . . . . . . . 43
6-6. Algoritmo de detección de DPs . . . . . . . . . . . . . . . . . . . . . . . . . 44
7-1. Encendido eléctrico cocina . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

7-2. Medidas de ensayo controlado . . . . . . . . . . . . . . . . . . . . . . . . . . 46
7-3. Montaje generación de pulsos controlados . . . . . . . . . . . . . . . . . . . . 47
7-4. Limpieza de ruido de DPs superficiales detectadas con dos antenas . . . . . . 47
7-5. DPs Corona detectadas con dos antenas . . . . . . . . . . . . . . . . . . . . 48
7-6. Descargas tipo Corona . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
7-7. DPs Corona detectadas con dos antenas . . . . . . . . . . . . . . . . . . . . 51
7-8. Montaje generación DPs Superficiales . . . . . . . . . . . . . . . . . . . . . . 51
7-9. Descargas tipo Superficial . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
7-10.DPs Superficiales detectadas con dos antenas . . . . . . . . . . . . . . . . . . 53
7-11.Barrido en frecuencias 100 MHz - 900 MHz . . . . . . . . . . . . . . . . . . 54
A-1. Antena Wifi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

A-2. Antena JCG015 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
A-3. Antena QUAD GSM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
Índice de cuadros
4-1. Parámetros estructurales para el transformador de potencia simulado . . . . 23
4-2. Amplitud del pulso Vs Th . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4-3. Materiales usados en la simulación . . . . . . . . . . . . . . . . . . . . . . . 24
4-4. Ubicación de sondas de medición . . . . . . . . . . . . . . . . . . . . . . . . 24
Contenido
Resumen VII
Lista de figuras IX
Lista de tablas XI
1. Introducción 2
2. Marco teórico y estado del arte 4

2.1. Clasificación y espectro electromagnético de las DPs. . . . . . . . . . . . . . 4
2.1.1. Descargas parciales (DPs) . . . . . . . . . . . . . . . . . . . . . . . . 4
2.1.2. Proceso de generación de DPs . . . . . . . . . . . . . . . . . . . . . . 4
2.1.3. Clasificación de DPs según su origen . . . . . . . . . . . . . . . . . . 8
2.1.4. Métodos de medida de descargas parciales . . . . . . . . . . . . . . . 10
2.2. Modelos de DPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2.1. Modelo capacitivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2.2. Modelo Dipolo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3. Técnicas de detección de fallos . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.4. Radio definido por software . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.5. Estado del arte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3. Planteamiento del problema y objetivos 20

3.1. Descripción del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.2. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.2.1. Objetivo General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.2.2. Objetivos Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4. Espectro electromagnético producido por DPs en transformadores de potencia. 22

4.1. Proceso de simulación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4.1.1. Modelamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4.1.2. Preprocesamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.1.3. Post-procesamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.2. Resultados de la simulación . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.2.1. Dominio de tiempo . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.2.2. Dominio de frecuencia . . . . . . . . . . . . . . . . . . . . . . . . . . 26
Contenido 1
4.2.3. Conclusión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5. Tecnología análisis de espectro electromagnético y montaje para generación

de DPs 28
5.0.1. Selección antenas para medir DPs . . . . . . . . . . . . . . . . . . . . 28
5.0.2. Dispositivo para la captura del espectro EM basado en radio definido
por software . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
5.0.3. Captura banda base de 60 Hz . . . . . . . . . . . . . . . . . . . . . . 32
5.0.4. Procesamiento de la señal . . . . . . . . . . . . . . . . . . . . . . . . 33
5.0.5. Circuito para generación de DPs . . . . . . . . . . . . . . . . . . . . . 34
5.0.6. Conclusión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
6. Algoritmo de adquisición, procesamiento y detección de DPs 38

6.1. Estructura General del código . . . . . . . . . . . . . . . . . . . . . . . . . . 38
6.2. Gnu Radio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
6.3. Post procesamiento Python . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7. Pruebas experimentales y Resultados 45

7.1. Pruebas preliminares experimento controlado . . . . . . . . . . . . . . . . . . 45
7.2. Denoising . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
7.3. Detección de DPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
7.4. Descargas tipo Corona . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
7.5. Descargas tipo Superficial . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
8. Conclusiones y recomendaciones 55
8.1. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
8.2. Recomendaciones y trabajos futuros . . . . . . . . . . . . . . . . . . . . . . . 56
8.3. Aportes del trabajo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
A. Medición de antenas usadas (Rohde y Schwarz Spectrum Analyzer) 58
Bibliografía 61
1. Introducción
Según el CIGRE (Comité Andino del Consejo Internacional de Grandes Redes Eléctricas),
el mayor porcentaje de fallos en transformadores de potencia se localiza en la parte activa
principalmente en los devanados del transformador. La ruptura del dieléctrico y la generación
de DPs (descargas parciales) entre devanados es una de las fallas más comunes; son un
fenómeno estocástico que se puede describir en tres estados:
El primer estado se presenta como streamer creando corrientes de 10 a 100 mA, al segundo
estado se le conoce como townsend donde la corriente disminuye, por último, se presenta un
estado conocido como pitting. Estas DPs representan un alto grado de degradación del siste-
ma por sobrecalentamiento, ruptura del dieléctrico, envejecimiento mecánico, deformación,
oxidación, entre otros ocasionando una falla en su primera fase incipiente y con el tiempo en
una falla funcional.
Existen varias técnicas de detección y localización de DPs en transformadores, en la norma
IEC 60270 se especifican las técnicas de detección de DPs que se producen en equipos eléctri-
cos de tensión alterna y continua. La medida eléctrica convencional de DPs es la más común
de las técnicas de detección planteando tres circuitos que se diferencian por la disposición
de la impedancia de medición con respecto al objeto de medida, su base de funcionamien-
to es generar un camino de baja impedancia al pulso generado por la DPs. Sin embargo
esta técnica requiere desconexión del equipo de alto voltaje, también existen técnicas no
convencionales IEC 62478 que no requieren la desconexión y no invasivas.
Las DPs son eventos individuales y además están acompañados de emisiones de luz, sonido,
pulsos electromagnéticos, produciendo reacciones químicas y calor.
La norma IEC-62478 acopia las técnicas de detección no convencionales, una de estas técnicas
es la detección del espectro electromagnético generado por una DPs. La detección se puede
hacer en el dominio de la frecuencia o en el dominio del tiempo; se hace necesario el uso de
sensores adecuados con un amplio ancho de banda debido a que las DPs generan interferencia
electromagnética produciendo ondas desde HF hasta el rango de las UHF (3 MHz hasta 3
GHz).
Uno de los desafíos más grandes en detección por este método es la ubicación de las antenas
en el transformador, hasta el momento se usan dos métodos: en el primero son ingresadas
en una ventana en la que se toman muestras de aceite, sin embargo solo dos antenas pueden
ser ingresadas. El segundo método consiste en abrir agujeros en el tanque del transformador
por lo que cuenta con una gran restricción en el momento de la instalación.
El mayor problema de la ubicación de estas antenas es la variación que pueden generar en
3
el campo eléctrico debido que algunos puntos del transformador son altamente críticos tal
como alrededor de los devanados. Una metodología para identificar la mejor posición del
arreglo de antenas es expuesta en este trabajo.
Las principales ventajas de este método detección en comparación a los métodos clásicos
de baja frecuencia presentadas en la norma IEC 60270 son: mayor inmunidad a ruido y
perturbaciones, ubicación de origen de la DP. Además de la posibilidad de hacer localización
de la fuente de DPs
En Colombia el alto precio de tecnologías de detección y diagnóstico de fallos impide que
muchas empresas que distribuyen energía tengan acceso a estos sistemas. La solución a nivel
nacional es la detección por medio de una medida directa como lo establece la norma IEC
60270. Este método requiere la desconexión del equipo lo que genera costos por interrupción
del sistema eléctrico.
Las metodologías nombradas en la norma IEC-62478 no requieren la desconexión del com-
ponente, por lo que se ofrecen alternativas en disminución de costos de diagnóstico. La lo-
calización de un fallo incipiente evitando que los transformadores tengan una falla funcional
representa una disminución drástica en lo que refiere a costos de mantenimiento correctivo
y reemplazo de componentes.
La técnica planteada se basa en el análisis del espectro electromagnético producido por las
DP y se propone una técnica basada en detección de microondas, una práctica novedosa y
poco estudiada a nivel nacional en fallos en transformadores de potencia.
2. Marco teórico y estado del arte
2.1. Clasificación y espectro electromagnético de las

DPs.
2.1.1. Descargas parciales (DPs)
La norma IEC 60270 define una descarga parcial como una descarga eléctrica y localizada
que cortocircuita sólo parcialmente el aislamiento entre conductores y que puede ser o no,
adyacente a un conductor [Standard, 2001]. El proceso de generación de DPs está condicio-
nada al medio de generación, el cual debe tener aunque sea parcialmente moléculas de gas,
provocadas por la vaporización del liquido aislante.
2.1.2. Proceso de generación de DPs

Cuando una descarga parcial ocurre asume diferentes formas: (Spark-times-True Puseless
glow-Pseudoglow). Todas estas descargas son emisiones del cátodo hacia el ánodo conocidas
como emission-sustained discharges y son esencialmente townsend discharges. Estas descar-
gas se caracterizan por tener un débil grado de ionización [Bartnikas, 2002].
En la Figura 2-1 se presenta el proceso de generación de DPs en un dieléctrico.
Figura 2-1.: Vector de fallo DPs.

2.1 Clasificación y espectro electromagnético de las DPs. 5
Aparición de micro espacios en el material aislante(t=0).

En un primer momento t=0 se empieza a apreciar micro-espacios conocidos como va-
cuolas en el aislante, producidos por: envejecimiento del mismo, humedad, sobrecarga,
entre otros. Estos micro-espacios pueden ser representados como un dipolo en el que
un campo electromagnético es producido por la influencia de cargas a una distancia de
separación entre el ánodo y cátodo de la vacuola [Shibuya et al., 2010].
Determinación Electrones disponibles (t=1).

En t=1 el proceso de ionización inicia por photoemision lanzando electrones que viajan
impulsados por el campo eléctrico desde el cátodo hacia el ánodo creando n0 cargas
por segundo como se observa en la Figura 2-2. Cuando todas las cargas son recogidas
por el ánodo se genera una corriente de saturación, esta descarga es débil, difícilmente
visible y depende de una fuente externa.
Figura 2-2.: Número de electrones por segundo impulsados por el campo eléctrico.
Townsend Discharge generados por los coeficientes de ionización (t=2).

Este fenómeno es observado por Townsend en 1900 introduciendo un termino conocido
como “First ionization coefficent” representado por α (Número de ionizaciones por
electrón por unidad de longitud a lo largo del campo eléctrico E).
Figura 2-3.: Cantidad de electrones ionizados por probabilidad de ionización

6 2 Marco teórico y estado del arte
Por definición; si es la cantidad de electrones llegando a la posición x y siendo α la

probabilidad de ionización por unidad longitud x + dx
El incremento del número de electrones debido a las ionizaciones es igual al electrón
que arriba a una nueva posición dx por la probabilidad de que se genere una nueva
ionización ver ecuación (2-1):
dN = N · α · dx (2-1)
Integrando
Z N Z X
dN
= α · dx (2-2)
N0 N 0
Si α es independiente de x entonces la cantidad de electrones cruzando por la nueva

posición dx incrementa de manera exponencial ver Figura 2-3 a lo que se le conoce
como avalancha.
Figura 2-4.: Interacción de electrones y iones en el dipolo
Ṅ = N0 · eα·d (2-3)
N0 es la cantidad de electrones por segundo, en un tiempo t > 0 se alcanzara una

corriente de saturación, entonces la corriente que atraviesa la vacuola de gas es igual
a: ecuación 2-4 donde i0 es la corriente de saturación
Por cada electrón ionizado se crea un ion positivo que viaja en dirección contraria por
lo tanto la corriente que atraviesa el circuito es siempre la misma ver Figura 2-4.
i = i0 · eα·d (2-4)
Descarga auto-sostenida(t=3)
Ahora en la Figura 2-5 se observa como iones cargados positivamente arriban al cá-
todo y la avalancha continua hacia el ánodo. La mayoría de iones son neutralizados
en el cátodo pero algunos pueden causar un emisión secundaria provocando una nue-
va avalancha, entonces la corriente puede incrementar con cada nueva avalancha de
electrones.
El arribo de iones cargados positivamente al cátodo es condicionado a una probabili-
dad, conocida como segundo coeficiente γ. La suma de series de avalanchas da como
resultado una corriente caracterizada por la ecuación 2-5
(a) Proceso de generación de mul- (b) Estado de la descarga

tiples avalanchas
Figura 2-5.: Proceso de generación de múltiples avalanchas y corriente de la DPs
i = i0 · eα·d /1 − γ(eα·d − 1) (2-5)
Con el tiempo y debido al fenómeno repetitivo de las DPs se producen cambios en

la geometría de la vacuola, provocando un aumento en la corriente que genera una
descarga conocida como streamer discharge en la Figura 2-5.
Cambios químicos, aparición del fallo(t=4)

Por último, se produce un cambio químico debido al incremento de temperatura en el

medio dieléctrico y desgaste mecánico. Como consecuencia final, el árbol de la descarga
en una cavidad de mayor tamaño se produce ruptura del dieléctrico.
2.1.3. Clasificación de DPs según su origen

Las DPs pueden ser clasificar principalmente en tres tipos: internas, superficiales y tipo
corona.
DPs Superficiales
Las descargas superficiales ocurren entre dos materiales aislantes ver Figura. 2-6 La
diferencia en permitividades produce el aumento del campo en cualquiera de los ma-
teriales y esto puede ocasionar una ruptura. [Calva et al., 2010].
Figura 2-6.: Descargas superficiales
DPs Corona
Las descargas tipo corona ocurren normalmente entre dos electrodos de pequeño ra-
dio de curvatura rodeados por medio gaseoso. Los electrodos tienden a tener forma
circulares la descarga toma una forma circular o corona. Se producen en el aire y el
elemento que actúa como electrodo, estas descargas no tienen ruptura total en los dos
electrodos.
DPs Internas
Las DPs internas son producidas por la existencia de burbujas en aislantes líquidos o es-
pacios en superficies sólidas con baja rigidez dieléctrica [Díaz Gautier and Martínez Tarifa, 2002].
Se pueden considerar como descargas de dos niveles que comienzan con la descarga de
la cavidad debido a huecos en el dieléctrico, esta cavidad está rodeada por electrodos o
completamente rodeada por aislamiento ver Figura 2-7 [Negm et al., 2017]. Se carac-
teriza por ser la más difícil de diagnosticar ya que el problema no presenta síntomas
visibles ni audibles.
Figura 2-7.: Descarga interna
Los aislantes están sometidos a una diferencia de potencial que produce un campo eléctrico
homogéneo en todos sus puntos Ea , cuando el aislante tiene vacuolas en el caso de aislantes
líquidos o imperfecciones en aislantes sólidos el campo eléctrico aumenta en esos puntos Ec
ver Figura 2-8.
Figura 2-8.: Campo Ea distribuido y campo de la vacuola Ec
En cuanto mayor sea la permitividad dentro de la imperfección o vacuola ε mayor será el

campo Ec con respecto a Ea ver ecuación 2-6. La tensión aumenta hasta romper la rigidez
dieléctrica de la burbuja produciendo una avalancha de electrones.
D = ε0 · εr · Ea (2-6)
Las DPs aparecen en una señal sinusoidal por la corriente AC del transformador con una
tensión conocida como Vi (tensión de ignición). El campo eléctrico generado por Vi genera
una avalancha de electrones descritas al empezar este capítulo. Cuando la tensión disminuye
hasta un punto conocido como Ve (tensión de extinción) la descarga parcial desaparece, el
proceso se repite cada vez que la tensión Va aumenta provocando que Vc (tensión en la
vacuola) aumente debido a la pendiente en la onda de tensión alterna. Las DPs se presentan
en los ciclos positivos como en lo ciclos negativos ver Figura 2-11.
La tensión de ignición Vi varia dependiendo del tipo de geometría de la vacuola. También
influye el aumento del estrés eléctrico dentro de un gas causando la ruptura de la rigidez
eléctrica, ya que la tensión de ignición es una variable que depende directamente de este
principio descrito por la ley experimental de Parchen; la cual establece que para una distancia
constante, la tensión necesaria para crear DPs en una cavidad disminuye a medida que la
presión también lo hace.
2.1.4. Métodos de medida de descargas parciales

Existen varias técnicas de detección y localización de DPs en transformadores. En la norma
IEC 60270 [Standard, 2001] se especifican las técnicas de detección de DPs que se producen
en equipos eléctricos de tensión alterna y continua.
La medida eléctrica convencional de DPs es la más común de las técnicas de detección,
planteando tres circuitos que se diferencian por; la disposición de la impedancia de medición
con respecto al objeto de medida.
La base de funcionamiento es generar un camino de baja impedancia al pulso generado por
la DPs.
Sin embargo esta técnica requiere desconexión del equipo de alta tensión (este método será
explicado mas adelante).
También existen técnicas no convencionales IEC 62478 [Standard, 2012] que no requieren la
desconexión y no son invasivas.
Las DPs son eventos individuales y además están acompañados de emisiones de luz, sonido,
pulsos electromagnéticos, produciendo reacciones químicas y calor ver Figura 2-9.
Estas técnicas tienen ventajas, limitaciones, bandas de operación y dependen de diferentes
tipos de sensores.
La más común de las técnicas utilizadas es la técnica basada en medición directa o modelo
eléctrico (20 KHz a 20 MHz), la principal ventaja de esta técnica es la sensibilidad y calidad
de la señal adquirida sin embargo no puede ser usada con el transformador en servicio.
Otra técnica muy común está basada en la detección de las ondas mecánicas (150 KHz)
generadas por las DPs; la principal ventaja de esta técnica es que puede ser usada en sitio,
sin embargo tiene poca sensibilidad.
Por último la técnica basada en el análisis de espectro electromagnético generado por las
DPs es conocida como UHF (300 Mhz a 3 GHz); la principal ventaja de esta técnica es que
2.2 Modelos de DPs 11
puede ser usada con el transformador en servicio, la desventaja de esta técnica es que esta
expuesta a diferentes fuentes de ruido por la banda de radio, televisión, GSM, WI-FI, etc. Es
una técnica que hasta ahora se está reglamentando por lo que no hay método de calibración.
Optical Effects
Optical
(Light)
Acustic
Pressure Wave
Mechanical
(Sound) Opto-acustic
IEC 62478
Discharge Effects IEC 60270

PD Dielectric Losses High Electrical
Frecuency Wavez
HF/VHF/UHF
Chemical Effects Chemical
Heat
Figura 2-9.: Métodos de medida de Descargas parciales
2.2. Modelos de DPs

2.2.1. Modelo capacitivo
Un modelo capacitivo describe las DPs. Consiste en el uso de un circuito equivalente y
descrito por la norma IEC 60270, en el que se mide la carga aparente. German y Philip-
pof en 1993 [Achillides et al., 2008] desarrollaron el primer circuito equivalente basado en
capacitores para simular DPs ver Figura 2-10.
La descripción del modelo capacitivo resulta del paralelo entre un capacitor Ca que representa
la capacitancia del dieléctrico restante sano y dos capacitores en serie Cb y Cc que representa
la capacitancia del dieléctrico y la capacitancia de la imperfección respectivamente.
Simulación LTSPICE
Para comprender el concepto del circuito equivalente se lleva a cabo una simulación en
LTSPICE.
Las DPs pueden ser simuladas usando un circuito equivalente. El circuito comprende un
condensador de acople, una impedancia de entrada, un condensador de medición y el objeto
de medida representado por Ca, Cb, Cc ver Figura 2-10-c [Negm et al., 2017].
Como resultado de la simulación ver Figura 2-11 el circuito es sometido a una tensión de
entrada de 5KV y con un switch que representa el electrodo que descarga, el capacitor Cc es
Ca Cb Ca Cb
Va
Vb Va Vb
} Ct
Ck
Cc Cc h
Vc Vc
Zm MI
a b c
Figura 2-10.: Circuito equivalente
Figura 2-11.: Simulación en LTSPICE. Ciclo repetitivo de descargas parciales en una onda
sinusoidal de 60 Hz
la capacitancia de la imperfección. Cada vez que se supera la tensión de ignición se presenta

una descarga por ciclo y tiene una tensión en Cc y una corriente ver Figura 2-12 que atraviesa
Cc.
Corrient e de descarga
0.06
0.05
Corriente (I)
0.04
0.03
0.02
0.01
0.00
−2 −1 0 1 2
Tiem po Segundos 1e− 9+ 4.1665e− 2
Figura 2-12.: Corriente de DPs
2.2.2. Modelo Dipolo

En 1980 Pederson propuso un modelo basado en un dipolo y no en un capacitor, el modelo
de dipolo describe el evento DPs y la distribución del campo eléctrico entre los electrodos
de prueba ver Figura 2-13.
cc
i i0(t)
v(t)
d
cc
v(t)
Figura 2-13.: Modelo de dipolo
Las corrientes que atraviesan el dipolo son difíciles de observar, varios estudios revelan que
la forma de onda de la corriente de descarga puede ser asociada con una forma gaussiana
[Callender et al., 2017], [Boggs and Stone, 1982], [Li et al., 2012b], [Nanao et al., 2008]. Es-
te pulso es representado por dos parámetros: amplitud en mA y (FWHM) anchura a media

altura Th ver Figura 2-14.
Una DPs consiste en el flujo de electrones e iones que se mueven a través de un campo
eléctrico; la velocidad de los electrones a una distancia de separación (generalmente) muy
pequeña es mayor que la de los iones.
La medición de la corriente de descarga revelará un pulso (causado por electrones) de mayor
magnitud pero poca duración seguido de un pulso duración corta pero mas duradera.
La parte de la descarga que contiene el movimiento de electrones puede ser convenientemente
modelada como una función Gaussiana ver ecuación 2-8, donde Imax es la corriente máxima
de la DPs y Th es (FWHM) anchura a media altura. .
" #
2
t
i0 (t) = Imax exp − (2-7)
t0
T
T0 = √h (2-8)
2 log 2
Ip Tipico pulso gaussiano de DPs

0.0010 Th = 1e 09seg
Th = 2e 09seg
Th = 3e 09seg
0.0008 Th = 7.5e 10seg
Am plit ud en m A
Th = 1.5e 09seg
Th = 5.3e 09seg
0.0006
Th = 7e 09seg
0.0004
0.0002 Th
0.0000
0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00
Tiem po en segundos 1e 8
Figura 2-14.: Pulso Gaussiano
Se supone que la carga verdadera Q se mueve a una distancia d que corresponda a la longitud
de vacío ver Figura 2-13 . Este movimiento define la variación del tiempo dipolo eléctrico
ver ecuación 2-9.
" #
Z t 2
t ∆p t
ρ(t) = d I0 · exp − dt = erf (2-9)
0 t0 2 t0
La onda electromagnética generada en el tiempo y basada en el modelo de dipolo eléctrico

es analizada usando el vector de Hertz [Shibuya et al., 2010]
r
ρ· t−
Π= c (2-10)
4 · π · ε0 · r
En la ecuación ε0 y c son la permitividad y velocidad de la luz en el espacio libre res-
pectivamente ver ecuación 2-10. La variable r es la distancia entre el dipolo y el punto de
observación.
Los campos E y H de acuerdo con la teoría electromagnética se describen mediante las
siguientes ecuaciones ver ecuación 2-11 .

∂Π
E = ∇ × (∇ × Π), H = ε0 ∇ × (2-11)
∂t
Hz
y
Dipolo
Eléctrico
y Punto de
Ex
observación
Figura 2-15.: Dipolo eléctrico y punto de observación
Considerando que ρ se encuentra en el origen y alineado con el eje x ver Figura 2-15 y el punto
de observación es tomado en el eje y, ya que corresponde a la dirección de mayor radiación
del campo, se obtienen las siguientes expresiones de los campos E y H ver ecuaciones 2-12,
2-13 .

|y| 2  2 
   
∆p t − |y| |y|
c  t − c   ∆p  t − c  
Ex = exp −  − exp −  
2 · ε0 · c2 · π 3/2 · t30 · |y| t0 4 · ε0 · c · π 3/2 · y 2 t0
    
|y|
 
∆p t−
− erf
 c 
8 · ε0 · π · y 3 
t0

(2-12)

|y| 2  2 
   
∆p t − |y| |y|
c  t − c   ∆p  t − c  
Hz = − exp −  + exp −  
2 · c · π 3/2 · t30 · |y| t0  4 · π 3/2 · t0 · y 2 t0
   
(2-13)
El campo de radiación que se estudiará se define en el primer término de las ecuaciones 2-12
y 2-13.Teniendo en cuenta esto se va a cambiar el punto de observación, definido por D como
se muestra en la Figura 2-16, de allí el campo eléctrico se expresa como ver ecuación 2-14.
|y|
∆p · (τ ) t−
Ex = exp[−(τ )2 ] τ = c (2-14)
2 · ε0 · c2 · π 3/2 · t30 · D t0
De la ecuación 2-14 el primer termino de la ecuación corresponde al valor pico del campo
eléctrico y el segundo a la derivada de función Gaussiana. Por lo tanto la forma de la onda
se describe con la siguiente función ver ecuación 2-15.
∆p √
Ex = √ 2
· [ 2e · τ · exp[−(τ )2 ]] (2-15)
2· 2 3/2
2e · ε0 · c · π · t0 · D
La intensidad del campo radiado depende de la rapidez con la que cambie el dipolo fuente.
Por lo que se puede deducir que el campo eléctrico seguirá la derivada del pulso Gaussiano
ver Figura 2-17.
Hz
(x,y,z)
y
Onda Ex
Electromagnética
Dipolo
Eléctrico
Δp
D
x
Figura 2-16.: Campo eléctrico lejano producido por DPs

[Shibuya et al., 2010]
2.3 Técnicas de detección de fallos 17
Figura 2-17.: Derivada pulso Gaussiano
2.3. Técnicas de detección de fallos

Los diversos métodos y técnicas utilizados para detectar fallos pueden clasificarse de acuerdo
con la naturaleza del conocimiento disponible sobre el proceso [Aranguren Z. and Tarantino A., 2009].
Detección basada en métodos analíticos: Utiliza herramientas matemáticas o analíticas

(modelos matemáticos precisos, procesado de señales) para realizar sus funciones.
Detección basada en conocimiento: Utiliza el conocimiento de un experto en el fenó-

meno estudiado incluye herramientas de la Inteligencia Artificial.
Existen distintos tipos de detección de fallos basados en técnicas como: analíticas, de cono-
cimiento, modelos, señales y síntomas. Cada uno de estos tipos de detección de fallos tiene
sus fortalezas y debilidades.
Dependiendo de las características del sistema y del número de señales disponibles por va-
riable en este caso DPs ver Figura 2-18. Por la naturaleza repetitiva y probabilística de una
DPs un sistema de detección basado en un análisis estadístico puede ser implementado en
el sistema planteado. Este método sera descrito en la sección de post procesamiento de la
señal.
Figura 2-18.: Métodos de detección de fallos según Isermann
2.4. Radio definido por software

En 1992 Joseph Mitola introduce el Radio Definido por Software (SDR, por sus siglas en
inglés) diseñado para dispositivos basados en comunicaciones inalámbricas.
El objetivo de estos dispositivos es reprogramar y reconfigurar sus características de ad-
quisición y generación por medio de herramientas basadas en software, para incrementar la
cantidad de funciones de los mismos usándose para diferentes estándares de comunicación.
El SDR se compone por una antena, un convertidor analógico digital (ADC) en la parte de
recepción y un convertidor digital analógico (DAC) y otra antena en la etapa de transmisión
[Autores, 2018].
2.5. Estado del arte

Ricardo Albarracin, realiza una investigación en la que se realizan medidas de descarga
parciales usando diferentes tipos de antenas, usando un analizador de espectro para la cap-
tura de las señales producidas por tres tipos de electrodos (superficial, interna y corona)
[Albarracín et al., 2013].
A. J. Reid, M. D. Judd, R. A. Fouracre realizan un trabajo en el que se miden descargas
parciales usando el método planteado por la IEC60270 combinado con técnicas basadas en
la captura del espectro electromagnético [Judd, 2011].
2.5 Estado del arte 19
K. Raja, F. Devaux, and S. Lelaidier de igual manera realizan un estudio de descargas par-
ciales usando como método de generación los típicos electrodos orientando su investigación
al análisis de la señal adquirida en función de la onda base de 50Hz [Raja et al., 2002].
Tianyan Jiang y su equipo presentan simulaciones y experimentos para la detección de DPs
usando el método clásico de generación por medio de electrodos como objetos de ensayo las
señales de descarga parcial en ondas en radiofrecuencia están ligeramente influenciadas por
la reflexión de los componentes del transformador, mientras que son muy susceptibles a la
refracción de los componentes del transformador. [Jiang et al., 2013]
Una gran variedad de estudios que pueden ser consultados para conocer los detalles de cada
uno [Jahangir et al., 2017], [Sinaga et al., 2012], [Robles et al., 2013] entre otros, todos estas
investigaciones parten de los mismos factores realizando el típico montaje para la generación
y posterior adquisición en UHF.
Sin embargo en la búsqueda del estado del arte se empiezan a conocer en los últimos años
trabajos relacionados con SDR (Radio definido por software), en H. Mohamed y otros rea-
lizan trabajos relacionados con adquisición usando tarjetas de bajo costo basados en SDR
[Mohamed et al., 2016].
3. Planteamiento del problema y
objetivos
3.1. Descripción del problema

La confiabilidad de un sistema de energía depende del adecuado funcionamiento de cada uno
de los elementos que la componen. Los transformadores de potencia, el elemento más crítico
de una subestación, suele presentar fallos por motivos tales como: humedad, fugas, sobre
tensión, mala protección eléctrica, mala instalación, entre otras.
Actualmente se cuenta con varias técnicas de monitoreo en línea de parámetros clave en
transformadores como lo es el DGA (Dissolved gas analysis), análisis térmico entre otros.
El fallo más común en los transformadores de potencia son las descargas parciales que pueden
ser detectadas mediante varias técnicas como lo indica la norma IEC 60270. Es una falla
impredecible y catastróficas para el sistema [Fernández, 2013]. Generando un proceso de
degradación que evoluciona gradualmente hasta presentarse la falla, que en ocasiones puede
ser funcional.
El transformador es una máquina eléctrica y presenta una gran cantidad de ruido electromag-
nético de diferentes fuentes lo que implica un estudio detallado del tipo de transformador,
caracterización del fallo con el fin de calibrar correctamente la antena planteada y distinguir
descarga parcial en un ambiente de ruido.
La técnica planteada basada en técnicas de detección de microondas en la frecuenta UHF
plantea diferentes problemáticas como lo son; la ubicación de dichas antenas, la morfología y
el posterior procesamiento de la señal. Características de la técnica que aún están en estudio
por lo que se brinda una gran posibilidad de desarrollo de nuevas soluciones a cada una de
las fases de detección de DPs.
Las DPs s son las primeras manifestaciones de un posible fallo potencial. La detección de DPs
en la fase incipiente se refleja en una rápida toma de decisiones por parte del departamento
técnico de una subestación eléctrica, produce ahorro energético y por lo tanto disminución
en gastos de mantenimiento.
La detección oportuna de una degradación en su fase incipiente hace la diferencia entre some-
ter el transformador a una preparación mayor o llevar el activo a mantenimiento preventivo
basado en su condición, evitando un fallo potencial o en el peor de los casos un fallo funcional
lo que obligaría a realizar mantenimiento correctivo o sustituir el transformador.
3.2 Objetivos 21
3.2. Objetivos
3.2.1. Objetivo General
Desarrollar un sistema de detección de descargas parciales en transformadores de potencia
3.2.2. Objetivos Específicos

Estudiar el espectro electromagnético de las descargas parciales en transformadores de
potencia.
Seleccionar la tecnología para el análisis de espectro electromagnético.
Diseñar una antena apropiada para detección de descargas parciales en un transforma-

dor de potencia.
Desarrollar algoritmos de procesamiento para la detección de descargas parciales.
Caracterizar el desempeño del sistema en laboratorio.

4. Espectro electromagnético
producido por DPs en
transformadores de potencia.
Con el fin de comprender la interacción del campo EM generado por una descarga parcial
se propone realizar una simulación y resolver el objetivo propuesto ”Estudiar el espectro
electromagnético de las descargas parciales en transformadores de potencia ”.
4.1. Proceso de simulación

La simulación se realizó en el dominio del tiempo utilizando la solución CST Microwave Stu-
dio; la figura 4-1 muestra el proceso para este fin. Los pasos serán explicados a continuación
Solucionador
Modelamiento Preprocesamiento Posprocesamiento
electromagnético
• Interfaz gráfica • Selección del

• Exportar datos
• Importar CAD/3D método de
• Visualización del
• Pulso Gaussiano solución
campo
• Materiales • Enmallado
Figura 4-1.: Diagrama de flujo para simulación
4.1.1. Modelamiento
Importar modelo CAD/3D
Un transformador trifásico fue modelado en un software 3D-CAD e importado en CST ver
figura 4-2. Los parámetros estructurales son mostrados en la tabla 4-1
En esta simulación se usaron diferentes parámetros para el pulso Gaussiano representando
diferentes tipos de descarga [Shibuya et al., 2010] ver tabla 4-2
4.1 Proceso de simulación 23
Figura 4-2.: Modelo 3D
Tabla 4-1.: Parámetros estructurales para el transformador de potencia simulado

Parámetros en mm
Parte
Radio interior Radio exterior Largo Ancho Altura
Devanado H-V 265 300 - - 774
Devanado L-V 215 250 - - 774
Núcleo - 160,5 - - -
Tanque - - 2200 750 1500
Tabla 4-2.: Amplitud del pulso Vs Th

Th (ns) 0,75 1 1,5 2 2,4 3 5,3 7
Imax (mA) 1 1 1 1 27 1 20 0,8
Materiales
La selección de los parámetros involucrados en la simulación fueron tomados de [Xu et al., 2015],
donde la simulación se realizó en una sola fase del transformador modelado como devanado
sólido. Los parámetros del material se muestran en la tabla 4-3
24 4 Espectro electromagnético producido por DPs en transformadores de potencia.
Tabla 4-3.: Materiales usados en la simulación

Parámetros de materiales
Material
Categoría σ (s/m) µ r r
Núcleo Dieléctrico - 1.0 -
Devanado Conductor 2x10^4 2000 -
Tanque Conductor 5.8x10^7 1.0 -
Aceite Conductor 5x10^6 1.0 2.2
4.1.2. Preprocesamiento
Dado el uso del solucionador de dominio en el tiempo, la malla base es de tipo hexaédrico.
Se utilizó un enmallado multiescala para que cada detalle geométrico se mostrara de forma
correcta manteniendo para todos los casos un tamaño de malla no mayor a una décima
de longitud por onda. Antes de la simulación se debe definir las salidas deseadas. En este
estudio, parte del análisis consiste en evaluar las posibles ubicaciones para los sensores de
campo. Con este fin, las sondas de campos se ubicaron como se indica en la tabla 4-4.
Tabla 4-4.: Ubicación de sondas de medición

S1 S2 S3 S4
X -1022 1022 1022 -1022
Y 1350 1350 -50 -50
Z -200 450 -200 450
Las tabla muestra las sondas ubicadas en diferentes puntos del transformador, se seleccio-
naron las esquinas por ser los puntos mas lejanos y normalmente las válvulas de toma de
muestra del aceite DN50 Y DN80 suelen ser ubicadas cerca de las equinas 4-3.
Figura 4-3.: Ubicación antenas en el modelo de simulación

4.2 Resultados de la simulación 25
4.1.3. Post-procesamiento
Finalmente, los resultados obtenidos en el punto de observación seleccionado son visualizados
en el dominio de tiempo y en frecuencia. La representación en el dominio del tiempo permite
la apreciación del efecto del entorno geométrico complejo de propagación en la forma de
onda simple proporcionada en la entrada, también permite la evaluación del tiempo de
llegada para la localización de la DP. La representación en el dominio de frecuencia se utiliza
para identificar los componentes de frecuencia más relevantes y como estos se relacionan
con los parámetros de las DP y los detalles geométricos del transformador; este análisis es
útil para comprender los requerimientos para colocar un elemento de detección (ubicación,
polarización y banda) y evaluar las técnicas de procesamiento de señales de banda estrecha.
4.2. Resultados de la simulación

Por brevedad, solo se presenta un subconjunto de los resultados. Nos centraremos en los
puntos de observación S3 y S4 (respectivamente el más cercano y el más alejado de la DP) y
los pulsos con el valor Th de 1ns y 3ns. Los resultados se presentan por separado para cada
componente del campo: Ex , Ey y Ez .
4.2.1. Dominio de tiempo

Las figuras 4-4 y 4-5 presentan los resultados de tiempo en las sondas S3 y S4 respectivamen-
te. La forma de onda observada en S3 tiene un pico alto en la componente Y correspondiente
a la primera llegada, tiene también copias atenuadas-retrasadas debido a múltiples reflexiones
en el entorno complejo.
1e−4 Campo elect ico Th = 1ns en la posición S3 1e−5 Campo elect ico Th = 3ns en la posición S3
X 1.0 X
Y Y
Z Z
0.5 0.5
0.0
E (V/m)
E (V/m)
0.0
−0.5
−0.5 −1.0
−1.5
0 20 40 60 80 100 0 20 40 60 80 100
Tiempo en (ns) Tiempo en (ns)
(a) Sonda S3 1 ns (b) Sonda S3 3 ns
Figura 4-4.: Sonda S3 resultados en el dominio del tiempo

26 4 Espectro electromagnético producido por DPs en transformadores de potencia.
1e−4 Campo electrico Th = 1ns en la posición S4 1e−5 Campo electrico Th = 3n en la po ición S4

X X
Y Y
1 Z 2 Z
0 1
E (V/m)
E (V/m)
−1 0
−2 −1
0 20 40 60 80 100 −2 0 20 40 60 80 100
Tiempo en (ns) Tiempo en (n )
Figura 4-5.: Sonda S4 resultados en el dominio del tiempo
4.2.2. Dominio de frecuencia
La figura 4-6 muestra el espectro de los pulsos de corriente que presentas las DP para los
valores Th de 1ns y 3ns, mientras que las figuras 4-7 y 4-8 representan el espectro de los
componentes del campo eléctrico en las sondas S3 y S4 respectivamente. Se puede obser-
var como se esperaba que los componentes de frecuencia más relevantes están determinados
por la duración del pulso, mientras que las fluctuaciones rápidas de frecuencia se deben a
la respuesta de frecuencia de la cavidad del transformador y sus componentes. Dados los
valores típicos para Th , estos resultados sugieren centrarse en la banda alrededor de 1Ghz
y proporcionar pautas para aprovechar mejor el rango dinámico del sistema de sensores.
Otra observación es que no hay diferencia sustancial entre la información proporcionada por
cada uno de los tres componentes del campo, de modo que se puede elegir la polarización
del sensor para minimizar la perturbación del entorno del transformador, para facilitar la
instalación/reemplazo o favorecer cualquier criterio práctico.
Ca po eléctrico espectro en frec(encia S3 (Th=1 ns) Ca po eléctrico espectro en frec(encia S3 (Th=3 ns)
−250 −250
X X
−260 Y −260 Y
Z −270 Z
−270
−280
E dB(V/ )
E dB(V/ )
−280
−290
−290 −300
−300 −310
−310 −320
−320 −330
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
Frequency (GHz) Frequency (GHz)
Figura 4-7.: Sonda S3 resultados en el dominio de la frecuencia

4.2 Resultados de la simulación 27
Curren pulse frequency spec rum (Th=1 ns,Th=3 ns)

0 Th = 1e − 09seg
Th = 3e − 09seg
−20
|S_pd (f)|(dB)
−40
−60
−80
−100
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
frequency en Hz 1e9
Figura 4-6.: Espectro en frecuencia
Campo e −c(rico espec(ro en frec)encia S4 (Th=1 ns) −250

Campo e −c(rico espec(ro en frec)encia S4 (Th=3 ns)
X X
−250 Y −260 Y
−260 Z Z
−270
−270
E dB(V/m)
E dB(V/m)
−280
−280
−290
−290
−300
−300
−310
−310
−320
−320
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
Frequency (GHz) Frequency (GHz)
Figura 4-8.: Sonda S4 resultados en el dominio de la frecuencia
4.2.3. Conclusión
Los resultados de la simulación son útiles para obtener las pautas para el diseño de elementos
de detección, sus ubicaciones y los algoritmos de procesamiento posteriores. Se ha encontrado
que cualquiera de los componentes de campo proporciona esencialmente la misma informa-
ción para elegir la polarización del sensor teniendo en cuenta los requisitos básicos como
fácil instalación/remplazo, reduce el riesgo de nuevas descargas y el costo. También se ha
observado que la ubicación del sensor no tiene un efecto significativo en las formas de onda
observadas, por lo que este parámetro también puede seleccionarse según su conveniencia
y compatibilidad. Finalmente, con el transformador en operación, un sistema de sensores
múltiples puede proporcionar una ubicación aproximada de una DP en función del análisis
del tiempo de llegada.
5. Tecnología análisis de espectro
electromagnético y montaje para
generación de DPs
En este capítulo se selecciona la tecnología más adecuada para la adquisición y posterior
análisis de señales de radio frecuencia producidas por las DPs.
Primero se realizó un estudio del estado del arte para la medición de DPs usando técni-
cas UHF, basados en esto e identificando las limitaciones y necesidades de esta técnica se
seleccionaran los dispositivos para ser usados en este trabajo.
En el estudio de DPs y el uso de UHF se han planteado diferentes experimentos que se basan
en los tres tipos de DPs (interna, superficial, corona) y un circuito que emula la naturaleza
de la descarga, sin embargo se requiere un equipo especializado para realizar estos estudios
como analizadores de espectro, generadores de alta tensión y dispositivos de captura de la
banda base.
En este trabajo, se ha planteado hacer más accesible este procedimiento con equipo disponible
en la Universidad Nacional de Colombia. Se pretende usar dispositivos de menor precio y
promover el acceso de los dispositivos a la industria regional.
5.0.1. Selección antenas para medir DPs

Ubicación de antenas
El tanque del transformador tiene como objetivo principal blindar la parte activa del trans-
formador, por lo que existen pocas ubicaciones en las que se puede ingresar una antena para
la medición de DPs. La posición mas común para el ingreso de antenas al transformador es
por medio de las válvulas de drenaje estandarizadas DN50 Y DN80 ver Figura 5-1 ubicadas
según el diseñador, pero mas comúnmente en la parte baja del tanque del transformador.
Una antena monopolo encaja con la geometría requerida para el ingreso de antenas por
medio de las válvulas de drenaje. Las dimensiones varían entre 100 mm hasta 30 mm según
el rango de frecuencia para el que se quiera diseñar. Es necesario mantener la longitud en
un valor relativamente pequeño para evitar cualquier interferencia con la parte activa del
transformador. Una longitud de alrededor de 100 mm se considera un riesgo razonable para
la la integridad del transformador [Griffiths and Pilling, 2004].
29
(a) Ventanas de ingreso de sensores (b) Válvula

DN50
Figura 5-1.: Ingreso de antenas al transformador [Siegel et al., 2017]
Las dimensiones de la antena deben confirmarse siempre de acuerdo con el diseño del trans-
formador específico y el punto de inserción de la antena en el tanque.
Según la bibliografía consultada las antenas mas utilizadas para la detección de de DPs
son monopolos [Akbari et al., 2016] por la comodidad en su geometría para ser insertados
en el transformador [Jin et al., 2012], también son utilizadas antenas de banda ancha como
espirales y log-espirales [Li et al., 2012a]. Estas antenas tienen una banda mas ancha y son
ubicadas en ventanas dieléctricas del transformador, sin embargo son difíciles de ubicar por
la poca disponibilidad de ventanas dieléctricas en transformadores antiguos.
Este trabajo plantea un diseño de un monopolo ver Figura 5-2 con una punta en forma de
gota con el fin de disminuir la acumulación de cargas en la antena como se demuestra en la
Figura 5-2(d), el monopolo con punta cilíndrica acumula mas campo en por lo que puede
provocar nuevas descargas ver Figura 5-2(c).
Por disponibilidad en el laboratorio de comunicaciones de la Universidad Nacional tres an-

tenas fueron usadas para este proyecto ver Figura 5-3.Los parámetros S11 son presentados
en el anexo A y fueron medidos con el VNA (Rohde y Schwarz Spectrum Analyzer), modelo
FSH8 con un rango de frecuencia de 100KHz a 8GHz.
5.0.2. Dispositivo para la captura del espectro EM basado en

radio definido por software
Hardware para el uso de SDR
En los últimos años se han desarrollado diferentes dispositivos basados en SDR siendo los
mas comunes HackRF producido por Great Scott Gadgets, el bladeRF producido por nuand
y el USRP producido por Ettus.
La Tabla 5.0.2 compara las características de las diferentes opciones de Hardware disponi-
bles en el mercado actual, actualmente otros proyectos se han desarrollado pero no están
disponibles en el mercado o son muy recientes y difíciles de adquirir.
Basado en la Tabla 5.0.2 se selecciona como dispositivo de adquisición de datos la USRP
30 5 Tecnología análisis de espectro electromagnético y montaje para generación de DPs
(a) Monopolo (b) Gota
(c) Antena monopolo sometido a campo eléctrico
(d) Antena gota sometida a campo eléctrico
Figura 5-2.: Antenas modelo CAD en CST, a) modelo monopolo, b)modelo gota, c) Acu-
mulación de cargas en el monopolo y d) Acumulación de cargas en la gota
B210 por su versatilidad, disponibilidad y fiabilidad.
El código fuente UHD disponible para ser usado con GNU Radio está bien documentado y
la comunidad de desarrollo es amplia lo que facilita el uso de la misma.
Por último se pretende utilizar un dispositivo RTL a manera de prueba por su valor econó-
mico y con el fin de proyectar un sistema mas asequible para un desarrollo futuro.
31
Figura 5-3.: Antenas utilizadas para la detección de DPs de izquierda a derecha: Antena
QUAD GSM, Antena Wifi, Antena JCG015
Hack RF BladeRF RTL USRP B210

Rango de espectro 30 MHz a 6GHz 300 MHz a 3,8 GHz 500 kHz a 1766 MHz 60 MHz a 6 GHz
Bandwidth 20 MHz 28 MHz 2.4 MHz 61.44 MHz
Duplex Half Full no Full MIMO
ADC/DAC 8 bit 12 bit 8bit 12 bit
Frecuencia Maximo
20 Msps 40 Msps 2.048 Msps 61.44 Msps
muestreo(ADC/DAC)
Velocidad del puerto USB 2.0 USB 3.0 USB 2.0 USB 3.0
de comunicación (480 megabit) (5 Gigabit) (480 megabit) (5 Gigabit)
Microcontrolador LPC 43XX Cypress FX3 R820T Cypress FX3
HDL +
Open Source Disponible Disponible Host Code
Code schematics
# Entradas 1 RX 2 RX 1RX 2 RX
Disponibilidad no no si si
National
Fabricante Scott Gadgets nuand Realtek
Instruments
Costo $300,00 $650,00 $40,00 $1.100,00
USRP B210
La USRP B210 ver Figura 5-4 producida por Ettus linea de National Instruments es selec-
cionada por su rango de frecuencias 70 MHz a 6 GHz, fue diseñada para experimentación
de bajo costo. Las características son dos entradas y dos salidas con posibilidad de hacer
cambios en el frente de radio, un reloj, opcionalmente un dispositivo GPS y un conector de
USB 3.0.
El RFIC AD9361 proporciona hasta 56 MHz de ancho de banda en tiempo real, una FPGA
Spartan 6 abierto y re-programable y conectividad a través de una conexión USB 3.0.

La librería de software USRP Hardware Driver (UHD) permite el desarrollo con GNU Radio.
(a) (b)
Figura 5-4.: USRP B210 y Arquitectura USRP B210
5.0.3. Captura banda base de 60 Hz

Como se demostró en el capitulo anterior la señal de 60 Hz es un factor importante en el
análisis de DPs.
La clasificación de DPs suele hacerse con métodos basados en la ubicación de la descarga en
la onda de 60Hz, estos métodos son conocidos como diagramas PRPD o patrón resuelto en
fase [Kreuger, 1989].La banda en la que esta centrada está descrita en el estándar IEC60270
y se encuentra por debajo de los 500 KHz, representa la corriente y el voltaje producido por
DPs y se visualiza por medio de la onda seno de 60Hz en corriente alterna.
Las distribución de DPs superficiales en la onda seno aparecen cerca al máximo de cada
ciclo, las DPs corona sólo se reflejan en un semiciclo y las DPs internas aparecen en los pasos
por cero.
(a) (b) (c)
Figura 5-5.: Circuito adquisición y tarjeta de sonido
Para la captura de la señal sinusoidal de 60 Hz disponible en el laboratorio de alta tensión

de la Universidad Nacional de Colombia se utiliza una tarjeta de sonido Logitech ver Figura
5-5-c con una frecuencia de muestreo máxima de 48 KHz.
33
Antes de ingresar la señal por la tarjeta de sonido se adecua la señal en el rango de los mA que
también aísla la parte de alta tensión con la parte de baja tensión usando un opto acoplador
en medio ver Figura 5-5-a-b, es importante la realización de este circuito para evitar daños
en los equipos de baja tensión como la USRP, el computador y demás accesorios conectados
a los mismos.
Este circuito se compone de cuatro etapas de funcionamiento:
La primera etapa se encarga de reducir el voltaje al nivel de los mV.
La segunda etapa corresponde a el voltaje es montado en un nivel DC con el fin de

poder utilizar cualquier tipo de opto acoplador y no opto acopladores diseñados para
entradas AC.
La tercera etapa consiste es un filtro RC y amplificador de la señal que limpia algunas

componentes de ruido de la onda de 60 Hz.
En la cuarta etapa la señal es opto-acoplada y posteriormente transmitida por medio

de un cable R-C-A de audio hacia una tarjeta de audio.
5.0.4. Procesamiento de la señal
Para procesar la señal se usó software libre bajo el sistema operativo Linux.
La señal fue transmitida haciendo uso de la USRP o RTL por medio de el puerto USB 3.0,
es procesada usando GNU-Radio diseñado exclusivamente para SDR es de código abierto y
altamente usado en investigación académica, esta construido con una interfaz amigable en
la que el usuario puede hacer uso de bloques de programación con diferentes funciones como
sincronización con hardware SDR, filtros, decodificadores, visualización entre otros muchas
aplicaciones para procesamiento de señales.
Figura 5-6.: Tecnologías usadas en el proceso de análisis de DPs
Para esta aplicación en concreto los bloques de GNU-Radio no satisfacen la necesidad del
cambio en el frente de radio o cambio del punto de frecuencia de adquisición, sin embargo
este software ofrece al usuario la posibilidad de crear bloques propios escritos en Python o
c++.
El paquete UHD esta escrito en c++ por lo que Python se usa como interprete de este,
además sera usado para el pos procesamiento de la señal y posterior visualización de los
resultados.
En la Figura 5-6 se presenta en forma de infografía un resumen de como proceso la señal y
las diferentes herramientas tecnológicas usadas.
5.0.5. Circuito para generación de DPs

En el marco teórico se revisaron los modelos de DPs (capacitivo y dipolo), en el modelo
capacitivo se describió el uso de un circuito de apoyo para la medida de DPS que consistía
en un capacitor de acoplamiento, una impedancia de entrada, un capacitor de medida y una
resistencia de medida. El diagrama de este circuito se presenta en la Figura 5-8-c.
Este el circuito implementado en el laboratorio de alta tensión para detectar DPs, es común
utilizarlo ya que tiene variaciones en cuanto a la ubicación de la impedancia de medición
descritas en la norma IEC60270. El uso de diferentes tipos de condensadores para modificar
la forma de onda producida por la descarga y la variación de la tensión de entrada del circuito
nos permite generar una prueba de laboratorio completa.
35
Por último el laboratorio de alta tensión de la Universidad Nacional es un ambiente seguro

para las pruebas y cuenta con los equipos necesarios para la generación de DPs.
(a) Cavidad Resonante
Figura 5-7.: Cavidad Resonante
El transformador es una maquina eléctrica que está contenida por una caja metálica, para
simular el ambiente de un transformador se construyó un tanque metálico ver Figura 5-7
que se usará en algunas mediciones, pruebas y contendrá electrodos de diferentes tipos que
serán descritos en la siguiente sección.
El tanque metálico es importante integrarlo en el análisis debido a los cambios que puede
generar en la captura ya sea por atenuación de la señal o por las ondas reflejadas como lo
observamos en la simulación de CST presentada en capítulos anteriores.
El comportamiento de la una cavidad resonante es descrita por autores como [Tozzi et al., 2010]
[Albarracín et al., 2013] que especifican las frecuencias teóricas en las que el tanque resuena,
en este trabajo este análisis no es tenido en cuenta debido al enfoque, para este caso se
tendrá en cuenta captura de las señales detectadas en los alrededores del tanque que serian
las señales que pueden provenir de los electrodos que generan las DPs dentro de la caja
metálica.
(a) Montaje de laboratorio (b) Etapa de Adquisición (c) Circuito Generador DPs
Figura 5-8.: Montaje de laboratorio Alta tensión

Electrodos para generar DPs internas, superficiales y corona
Para producir DPs se pueden usar defectos de tipo industrial que son ingresados de manera
artificial en maquinas eléctricas o defectos estándar, estos son producidos a partir de la
geometría de diferentes electrodos o dieléctricos en medio de los mismos ver Figura 5-9 .
Figura 5-9.: Clasificación de descargas parciales basado en electrodos. a) DPs internas, b)

DPs superficial, c) DPs corona [Kreuger et al., 1993]
Los electrodos utilizados como objetos de ensayo en el laboratorio de alta tensión de la Uni-
versidad Nacional de Colombia permiten simular una descarga parcial y su comportamiento
repetitivo en la onda de 60 Hz.
Como fue estudiado en el primer capitulo de esta tesis se presentó la clasificación de descarga
mas comunes en maquinas eléctricas (internas, superficiales y corona) [Kreuger et al., 1993]
se propone manufacturar electrodos tipo punta-semiesfera, barra-plano, planos paralelos. La
geometría utilizada fue planteada por J. Vanegas de la Universidad del Valle [Vanegas, 2016].
Los electrodos son construidos con tecnicas de manufactura en torno ver Figura 5-10.
Según el tipo de descarga (Corna, interna, superficial) se observan en diferentes puntos de
la onda de 60 Hz.
(a) Corona (b) Superficiales (c) Internas
Figura 5-10.: Electrodos para simular descargas tipo corona, superficiales e internas
37
5.0.6. Conclusión
La tecnología seleccionada se resume en la Figura 5-11. Consta de tres etapas: la primera
corresponde a la etapa de generación de DPs usando un circuito descrito por la norma IEC
60270 compuesto por un capacitor de acople, una impedancia de entrada y tres electrodos,
una segunda etapa de adquisición de datos usando una USRP B200 y una Tarjeta de sonido
USB que están conectadas por medio de USB3.0 a un computador con Linux como siste-
ma operativo. Por último una etapa de procesamiento de datos basada en software libre y
explicada en el siguiente capitulo.
z
PD electodos
CK
oltaje
Antena
Tarjeta de sonido
Figura 5-11.: Tecnología utilizada para la detección de DPs

6. Algoritmo de adquisición,
procesamiento y detección de DPs
Para seleccionar la metodología a utilizar en el sistema de detección de DPs influye en gran
medida las características de los fallos que se van a detectar.
Una de las principales características de las DPs es su repetitividad, como es una falla incisiva
en el tiempo aun siendo un proceso estocástico cuando comienza continuará en el tiempo
hasta producir una ruptura del dieléctrico lo que conlleva a un fallo funcional, sin embargo
es una falla que evoluciona lentamente por lo que una detección en linea no es necesaria,
esto no quiere decir que una detección del fallo en su estado incipiente no sea importante.
Basado en estas características se pretende identificar un método de detección adecuado para
este trabajo.
Así las cosas por las características de la descarga se puede inferir que el comportamiento de
este fallo será muy similar en todas las ventanas de un vector de tiempo amplio, por lo que
un manejo estadístico calculando las distancias de separación entre un pulso y la media de la
señal muestreada nos da un indicio de existencia o no de descargas parciales. Es factible usar
una técnica de detección descrita por Isermann [Isermann, 2006] conocida como gráficos de
control estadístico en la que se calculen umbrales mínimos y maximos fijos ver Figura 2-18.
6.1. Estructura General del código

Una estructura general del código es descrita en esta sección ver Figura 6-1. La señal UHF
adquirida por medio de la USRP y la señal banda base de 60Hz con la tarjeta de sonido
envían a través de USB la información que sera procesada con GNU-Radio.
Se guardan los datos procesados con GNU-Radio haciendo uso de una Ramdisk que es una
memoria de alta velocidad, la cantidad de muestras por segundo se puede aumentar además
de evitar el cuello de botella entre el puerto USB y el guardado en Disco duro magnético en
el caso mas común.
Una vez guardados los archivos en el Ramdisk se genera un archivo que es leído usando
una función especial de conversión de datos complejos a binarios programado en Python, es
entonces cuando los datos leídos estructuran un dataframe generado con la librería Pandas
disponible para Python.
Cuando el dataframe es construido se procede a hacer un post-procesamiento a las señales
capturadas que consiste en limpiar la señal de ruido y por ultimo detectar si hay presencia
6.2 Gnu Radio 39
o no de DPs en la muestra capturada.

En las siguientes secciones se explica paso a paso el algoritmo de adquisición, post procesa-
miento, detección y visualización de las señales producidas en el entorno de laboratorio.
Figura 6-1.: Diagrama de flujo Visión general del procesamiento de datos
6.2. Gnu Radio

GNU-Radio es una herramienta exclusivamente diseñada para el desarrollo de sistemas SDR
(Software definido por radio), nos permite interactuar con hardware dedicado para este fin
en este caso una USRP B210. Este software permite programar en diagrama de bloques, se
uso Python como lenguaje de programación debido al propósito de mejorar la velocidad de
adquisición y generar algunas funciones que no están disponibles en los bloques de GNU-
Radio.
El flujo de programación esta representado en la Figura 6-2 y se describe a continuación.
Primero se realiza la configuración de la tarjeta de adquisición de datos USRP B210 o RTL
según el caso, se carga la librería UHD que es la base de código de USRP, desde esta librería
se establecen las diferentes configuraciones y variables como cantidad de entradas (antenas en
nuestro caso), ganancia, ancho de banda, frecuencia de muestreo y frecuencia de operación.
Cuando los valores iniciales son establecidos en la tarjeta se procede a establecer una serie
40 6 Algoritmo de adquisición, procesamiento y detección de DPs
de variables para el control de la adquisición de datos, estas variables son: Cantidad máxima
de datos, frecuencia inicial, frecuencia final y cantidad de datos por frecuencia.
Una vez realizado el proceso de establecimiento de todas las variables usadas en el programa
se comienza con la adquisición de los datos iniciando el flujo de datos, se van escribiendo
en un archivo y una condición de datos máximos permite que el archivo no se desborde
del tamaño aceptado por el ramdisk. Por otra parte un contador aumenta hasta alcanzar la
cantidad de datos por frecuencia, al alcanzar esta cantidad se suma el ancho de banda a la
frecuencia actual y se vuelve a sintonizar la frecuencia de operación.
Por ultimo el archivo es guardado en el ramdisk y posteriormente se mueve al disco duro
para el siguiente paso que es el post-procesamiento de la señal.
Figura 6-2.: Diagrama de flujo Gnu-Radio
6.3. Post procesamiento Python

Los datos guardados en el archivo son requeridos desde Python. GNU-Radio escribe estos
datos con un formato complejo por lo que se usa una función que convierte este formato
complejo a binario y se procede a hacer el tratamiento de los datos capturados ver Figura
6-3.
6.3 Post procesamiento Python 41
En primer lugar se usa la librería Pandas para generar una base de datos que esta compuesta
de datos capturados por la tarjeta USRP, la tarjeta de sonido, la banda a la que pertenece
cada dato y la etiqueta de tiempo, en el dataframe se sincroniza la señal proveniente de la
USRP con la tarjeta de sonido.
Cuando los datos están debidamente ubicados y sincronizados en el dataframe se procede a
limpiar el ruido proveniente de los mismos usando la transformada de Wavelet.
Graficos de control
estadistico
Figura 6-3.: Post procesamiento de las señales
Transformada de Wavelet
La transformada de Wavelet es una transformada comúnmente usada en el procesamiento
de señales como audio, imágenes entre otros. Es particularmente usada para señales no
periódicas. [Wu et al., 2015].
La transformada permite separar en coeficientes la mayor parte de ruido y la señal de in-
terés, lo que mitiga o elimina el ruido a través de una función de contracción de umbral
[Cunha et al., 2015]. La implementación práctica de la transformada de la onda en un banco
de filtros se muestra en la Figura 6-4.
En la primera etapa la señal pasa a través de un banco de filtros de dos canales compuesto
por un filtro pasa-bajo G(n) correspondiente a la función de escala, y un filtro pasa-alto
Figura 6-4.: Modelo para descomposición de Wavelet

[Cunha et al., 2015]
H(n) correspondiente a la función de Wavelet; seguido de un factor de dos muestreos de

bajada.
Las salidas de la primera etapa son dos señales con longitudes iguales a la mitad de la
señal de entrada. Las dos señales son los coeficientes a1 de aproximación y d1 de detalle.
En la siguiente etapa los coeficientes a1 son aplicados a los filtros pasa-bajo y pasa-alto
seguido de los muestreos de bajada para generar un nuevo coeficiente de aproximación a2
y un coeficiente de detalle d2 . Este proceso de descomposición se repite hasta obtener un
coeficiente de aproximación final aj y un coeficiente de detalle dj .
El banco de filtros son de tipo espejo en cuadratura (QMF), son diseñados de tal manera que
sea posible tener una reconstrucción perfecta de la señal original aplicando la transformada
de onda digital inversa (IDWT), que se puede implementar con un banco de filtros de síntesis
con factor de dos interpoladores en cada nivel. La supresión de la señal por Wavelet se resume
en la Figura 6-5.
La libreria pywavelet fue utilizada es una librería escrita en Python y de código abierto.
Gráficos de control estadístico

Se basa en el establecimiento de dos limites superior U CLx ver ecuación 6-1 e inferior LCLx
ver ecuación 6-2. El análisis de la media y la desviación estándar permite registrar las varia-
ciones que sobrepasan cualquiera de estos limites [Aranguren Z. and Tarantino A., 2009].
Los valores aceptados se establecen por medio del calculo de la media aceptable, si la media
esta dentro de los limites superior e inferior no hay existencia de fallo, pero si la medida
supera el valor del rango inferior o superior se puede considerar como presencia de fallo, en
este caso DPs.
U CLx = X + 3αx (6-1)

Donde U CLx es el limite superior, x es la media de la muestra y α es la desviación estándar
de la muestra.
LCLx = X − 3αx (6-2)
6.3 Post procesamiento Python 43
Determinar número de
niveles de descomposición
"j" y elegir los filtros
Wavelet para cada nivel j.
Obtener los valores del

coeficiente de umbral
Reconstruir la señal
aplicando IDWT a los
coeficientes de umbral
obteniendo la señal sin
ruido
Figura 6-5.: Supresión de señal con técnica de Wavelet

[Cunha et al., 2015]
Donde LCLx es el limite inferior, x es la media de la muestra y α es la desviación estándar

de la muestra.
El diagrama de flujo del algoritmo de gráficos de control estadístico presentado en la Figura
6-6, comienza determinando la frecuencia de muestreo, ventana observación y tamaño de
muestra, continuando con el calculo de la media y la desviación estándar, por ultimo se
determinan los limites inferior y superior.
Figura 6-6.: Algoritmo de detección de DPs

7. Pruebas experimentales y
Resultados
En este capitulo se describen las pruebas realizadas en laboratorio y resultados obtenidos.
7.1. Pruebas preliminares experimento controlado

Una prueba preliminar a las mediciones en el laboratorio de alta tensión fue desarrollada;
la finalidad de esta prueba fue evitar el daño que se podía producir por el manejo de alta
tensión en los equipos de adquisición en el laboratorio.
Para reducir el riesgo se pensó en un generador de pulsos de descarga por ruptura en un
material aislante el cual tuviera un periodo de descarga constante, fácil de instalar y adquirir.
La solución que se planteó para esta prueba fue usar el encendido electrónico de una estufa
ver Figura 7-1. Este elemento cumple con los requisitos propuestos en especial el periodo de
descarga constante. Es entonces que se empiezan a realizar pruebas, primero en un ambiente
abierto lo que permite el desarrollo del algoritmo de detección de una manera confiable y sin
necesidad del uso del laboratorio de alta tensión.
Figura 7-1.: Encendido eléctrico cocina
Así las cosas, se presentan los primeros resultados ver Figura 7-2. La Figura (a) es el resultado
de un barrido de frecuencias sin existencia de descargas, la Figura (b) es el resultado de un
barrido de frecuencias con existencia de descargas. Estas pruebas fueron realizadas usando
una tasa de adquisición de 10 Msps, ganancia de 20 dBs, frecuencia inicial en 70 MHz y
frecuencia final en 1 GHz. Como se puede observar en la Figura 7-2 (b) los datos muestran
46 7 Pruebas experimentales y Resultados
presencia de descargas y además se observa que en las frecuencias mas bajas la señal se
adquiere con mayor intensidad.
Ensayo controlado - barrido de frecuencias
Ensayo controlado - Barrido de frecuencias
0.8
Piso de ruido Pulsos controlados
Amplitud piso de ruido E (V/m)
Amplitud de pulsos controlados

0.004
0.6
0.4
0.002
Campo E (V/m)
0.2
0.000
0.0
0.2
0.002
0.4
0.004
0.6
0 1 2 3 4 5 0 1 2 3 4 5
1e7 1e7
Número de datos Número de datos
(a) Piso de ruido - barrido de frecuencias (b) Datos capturados prueba de pulsos con-
trolados - barrido de frecuencias
Figura 7-2.: Medidas de ensayo controlado
La parte activa de un transformador esta embebida en un tanque metálico por lo que se

planteó construir una cavidad resonante (Caja metálica) ver Figura 7-3 simulando el com-
portamiento de propagación de la descarga en el transformador.
Un electrodo es ingresado en la caja metálica y con este sistema se realizó la optimización del
algoritmo con el objetivo de alcanzar la máxima tasa de muestreo sin overflow o underflow,
el cambio de frecuencias sin producir cambios en los datos adquiridos, ganancia adecuada,
distancias seguras a la fuente de descarga, cantidad de datos por frecuencia y por ultimo
cantidad de antenas conectadas en el puerto de recepción. Esta prueba permitió variar las
frecuencias centrales por canal, las ganancias, las distancias y demás parámetros en el sistema
de adquisición.
El sistema de pruebas preliminar fue de gran importancia para esta investigación, con este
se logró la optimización del sistema de adquisición en un ambiente controlado, obteniendo
los siguientes logros.
Tasa de muestreo de hasta 25 Msps con solo una antena conectada.
Máximo 4 antenas conectadas (2 USRP B210) en un solo computador con una tasa de
muestreo de hasta 5MHz por canal.
Cambio del frente de radio sin generar perturbaciones en las muestras.
7.2. Denoising
Con el sistema de adquisición optimizado se procede a hacer la limpieza de ruido usando
coeficientes de Wavelet con los que se obtienen los siguientes resultados ver Figura 7-4, en la
Figura (a) se observan datos que pertenecen a DPs superficiales. Haciendo zoom a un pulso
7.3 Detección de DPs 47
Figura 7-3.: Montaje generación de pulsos controlados
de esta serie de DPs, Figura (b), se ve el rendimiento del algoritmo de limpieza de ruido
que sigue la forma de pulso eliminando el ruido a la izquierda y derecha pero sin remover
componentes importantes de los pulsos, además se observa el comportamiento descrito por
la simulación en CST en el capítulo 4 ver Figuras 4-4,4-5.
Pulso de descarga procesado tranformada wavelet

0.025
Procesamiento de señal transformada wavelet

0.08
0.020
Datos Procesados
0.06 Datos base
0.015
Amplitud de DPs
Amplitud de DPs
0.04
0.010
0.02
0.005
0.00
0.000
0.02
0.005
0.04
0.010
0.06
0 100000 200000 300000 400000 500000

− 60 − 40 − 20 0 20 40
+ 2.18e5
Número de muestras Número de muestras
(a) Limpieza de ruido usando Wavelet (b) Pulso de descarga procesado con Wavelet
Figura 7-4.: Limpieza de ruido de DPs superficiales detectadas con dos antenas
Con el sistema optimizado se procede a realizar pruebas en el laboratorio de alta tensión en

el que se producen DPs en un ambiente aproximado al comportamiento en un transformador
de potencia. Los resultados mostrados a continuación pertenecen a cada uno de los tipos de
descarga que se estudian.
7.3. Detección de DPs

El algoritmo de detección (Gráficos de control estadístico) ubica los limites UCLx y LCLx
determinando la existencia de DPs en todas las bandas de frecuencias medidas ver Figura
7-5.
Descarga superficial vs Tiempo

Dps_Canal_A
0.100
Banda base 60 Hz
UCLx
LCLx
0.075
0.050
Dps (Normalizada)
0.025
0.000
( 0.025
( 0.050
( 0.075
( 0.100
0.0 2.5 5.0 7.5 10.0 12.5 15.0 17.5

Tiem po (m s)
(a) Limpieza de ruido usando Wavelet
Descarga Superficia (s Tiempo

Dps_Cana _A
0.010
Banda base 60 Hz
UCL)
LCL)
0.005
0.000
Dps (Norma izada)
−0.005
−0.010
−0.015
−0.020
5.294 5.296 5.298 5.300 5.302 5.304 5.306 5.308

Tiempo (ms)
(b) Pulso de descarga procesado con Wavelet
Figura 7-5.: DPs Corona detectadas con dos antenas
7.4. Descargas tipo Corona

El montaje descrito en el capítulo 5 ver Figura 5-8 (c) basado en la norma IEC 60270
[Standard, 2001], se implementó con una tensión de entrada V in de 8 KV, una resistencia
7.4 Descargas tipo Corona 49
en serie a la fuente de 10 Mohm en paralelo dos capacitores de Ck 2nf y Ca 100 pf, una
resistencia de medición Zm de 500 ohm y un capacitor de medición Cm de 200 nf con un
objeto de prueba punta corona.
Los resultados de las pruebas realizadas se presentan en diferentes escenarios haciendo cam-
bios en los parámetros de adquisición ver Figura 7-6. La onda base de 60 Hz, en esta Figura
la señal esta recortada debido a un problema de saturación en el convertidor análogo-digital
de la tarjeta de sonido pues la tensión de entrada supera los límites del mismo.
Escenario (a): Tasa de muestreo de 10 MHz para la USRP y 48 KHz para la tarjeta
de audio, ganancia 2 dBs, 3 m de distancia entre la antena y el objeto de prueba, la
antena Wifi acoplada en 2,46 GHz con polarización vertical 90 con respecto al eje de
la tierra; se toman los datos en una frecuencia central de 100 MHz.
Escenario (b): Igual que el descrito escenario (a) variando la frecuencia de central a
300 MHz.
Escenario (c): Igual que el descrito escenario (a) variando la frecuencia de central a
500 MHz.
Escenario (d): Igual que el descrito escenario (a) variando la frecuencia de central a
800 MHz.
Se puede notar que las descargas tipo corona contienen una descarga por semi-ciclo ver Figura
7-6 (b), este comportamiento ya observado y descrito por diversos autores [Jahangir et al., 2017];
también se aprecia un incremento de la repetición de la descarga en el semi-ciclo positivo
Figura 7-6 (c), esto debido a la tensión de alimentación del objeto de prueba.
Con los resultados obtenidos por [Zhou et al., 2018] cuando la tensión aumenta se presentan
mas descargas en el semi-ciclo positivo.
Como se demostró en la simulación realizada en CST en frecuencias mas altas la amplitud
de la descarga disminuye es decir que la respuesta experimental obedece a los resultados de
la simulación.
En la Figura 7-7 se usan dos antenas con diferente polarización y distancia de la fuente de
generación para la adquisición de los datos, estos se tomaron en una frecuencia central de
300 Mhz.
La antena del canal A tiene una mejor respuesta debido a la polarización pues es perpendi-
cular a la descarga, la antena en el canal b esta polarizada en paralelo al objeto de prueba.
Respecto a las distancias las dos antenas están separadas por 1 m.
En detalle no se nota un cambio significativo en el tiempo de llegada entre los pulsos mas
importantes debido a la velocidad de adquisición de la tarjeta pues las descargas se generan
en tiempos de ns y la resolución de la tarjeta no permite una respuesta que pueda diferenciar
tiempos de llegada entre dos pulsos generados por la misma descarga.
Descarga Corona s Tiempo Descarga Corona vs Tiempo

0.100 Dps Dps
Banda base 60 Hz Banda base 60 H
0.4
0.075
0.050
0.2
Dps (Normalizada)
Dps (Normali ada)

0.025
0.000 0.0
−0.025
−0.2
−0.050
−0.075
−0.4
−0.100
0 10 20 30 40 50 0 10 20 30 40 50
Tiempo (ms) Tiempo (ms)
(a) Frecuencia central 100 MHz a 7.5 KV (b) Frecuencia central 300 MHz a 7.5 KV
Descarga Corona s Tiempo Descarga Corona s Tiempo

0.100 Dps 0.100 Dps
Banda base 60 Hz Banda base 60 Hz
0.075 0.075
0.050 0.050
Dps (Normalizada)
Dps (Normalizada)
0.025 0.025
0.000 0.000
−0.025 −0.025
−0.050 −0.050
−0.075 −0.075
−0.100 −0.100
0 10 20 30 40 50 0 10 20 30 40 50
(c) Frecuencia central 500 MHz a 7.5 KV (d) Frecuencia central 800 MHz a 7.5 KV
Figura 7-6.: Descargas tipo Corona
En las pruebas realizadas variando voltajes, frecuencia central, distancia entre antenas y
ganancia, se detectaron descargas tipo corona con amplitud normalizada de campo eléctrico
de 0.4 y con tensión de 8 KV por lo que las antenas se pueden ubicar hasta una distancia
de 5 m para obtener una detección satisfactoria.
7.5. Descargas tipo Superficial

El montaje descrito en el capítulo 5 ver Figura 5-8 (c) basado en la norma IEC 60270
[Standard, 2001], diferenciado de la descarga tipo corona con una variación que consiste en
retirar el capacitor de descarga Ca, la tensión de entrada V in es de 3 KV, una resistencia
en serie a la fuente de 10 Mohm en paralelo dos capacitores de Ck 2nf, una resistencia de
medición Zm de 500 ohm y un capacitor de medición Cm de 200 nf con un objeto de prueba
(punta - acrílico - plano) ver Figura 7-8 .
En los resultados presentados en la Figura 7-9 se detectaron DPs superficiales.
Los resultados varían respecto a las DPs tipo corona en la cantidad de DPs por ciclo, lo
primero que hay que denotar es la cantidad de pulsos por semi-ciclo debido a que las DPs
se producen en diferentes puntos de la superficie del aislante.
El ciclo repetitivo característico de DPs es una por cada vez que se supera el voltaje de
ignición, comportamiento explicado en el marco teórico ver Figura 2-11.
7.5 Descargas tipo Superficial 51
De carga Corona v Tiempo

0.03 Dp _Canal_A
Banda ba e 60 Hz
Dp _Canal_B
0.02
Dp (Normalizada)
0.01
0.00
−0.01
0 10 20 30 40 50
Tiempo (m )
(a) Frecuencia central 300 MHz a 7.5 KV dos antenas
Figura 7-7.: DPs Corona detectadas con dos antenas
Figura 7-8.: Montaje generación DPs Superficiales
Las descargas superficiales al igual que las descargas tipo corona como se demostró en la
simulación realizada en CST ver Figuras 4-7, 4-8 tienen la característica de disminuir la
amplitud con respecto a la frecuencia, por lo que en frecuencias de hasta 500 MHz se detectan
pulsos significativos de descarga y cuando aumenta la frecuencia de adquisición la amplitud
disminuye, no implica menos información pues en frecuencias mas bajas se encuentra también
ruido proveniente de radio, televisión, entre otros.
Descarga Superficial vs Tie po Descarga Superficial vs Tiempo

Dps_Canal_A Dps_Ca al_A
0.06
Banda base 60 H) Ba da base 60 H(
0.08 UCL( UCLx
LCL( LCLx
0.06 0.04
0.04
Dps (Nor ali)ada)
Dps (Normali(ada)
0.02
0.02
0.00
0.00
−0.02 )0.02
−0.04
)0.04
−0.06
)0.06
0 10 20 30 40 50 0 10 20 30 40 50
Tiempo ( s) Tiempo (ms)
(a) Frecuencia central 100 MHz a 3 KV (b) Frecuencia central 300 MHz a 3 KV
Descarga Superficial vs Tiempo Descarga Superficial vs Tiempo

Dps_Ca al_A Dps_Ca al_A
0.06 Ba da base 60 H( 0.06 Ba da base 60 H(
UCLx UCLx
LCLx LCLx
0.04 0.04
Dps (Normali(ada)
Dps (Normali(ada)
0.02 0.02
0.00 0.00
)0.02 )0.02
)0.04 )0.04
)0.06
0 10 20 30 40 50 0 10 20 30 40 50
(c) Frecuencia central 500 MHz a 3 KV (d) Frecuencia central 800 MHz a 3 KV
Figura 7-9.: Descargas tipo Superficial
En la Figura 7-10 a) se capturaron DPs tipo superficial con dos antenas y un frente de radio
sintonizado de 300 MHz y b) dos antenas y un frente de radio sintonizado en 500 MHz. Se
observa el cambio de amplitud en 300 MHz y 500 MHz, las señales provenientes de la antena
polarizada en paralelo a la fuente de descarga son de mayor amplitud a la antena polarizada
ortogonal a la fuente de descarga, sin embargo las dos antenas capturan la señal en el mismo
instante de tiempo debido a la poca separación entre las mismas.
Se recomienda hacer la detección en el rango de frecuencias UHF. En los 800 MHz se detectan
pulsos mas amplios debido a que la antena usada en las pruebas resuena en esta frecuencia
ver anexo a A-2.
Por último un barrido en frecuencia fue realizado para observar con detalle el comportamiento
del cambio de amplitud de campo capturado para cada banda de frecuencia aumentando de
10 MHz por banda, se realizo desde los 100 MHz hasta los 900 MHz ver Figura 7-11.
7.5 Descargas tipo Superficial 53
Descarga Superficial vs Tiempo

Dps_Ca al_A
Ba da base 60 H(
Dps_Ca al_B
0.04 UCLx
LCLx
0.02
Dps (Normali(ada)
0.00
)0.02
)0.04
)0.06
0 10 20 30 40 50
Tiempo (ms)
(a) Frecuencia central 300 MHz a 3 KV dos antenas
Descarga Superficial vs Tie po

Dps_Canal_A
Banda base 60 H)
0.075 Dps_Canal_B
UCL(
LCL(
0.050
Dps (Nor ali)ada)
0.025
0.000
−0.025
−0.050
−0.075
0 10 20 30 40 50
Tiempo ( s)
(b) Frecuencia central 500 MHz a 3 KV dos antenas
Figura 7-10.: DPs Superficiales detectadas con dos antenas

)Presencia de DPs superficiales barrido de frecuencias

100 MHz 200 Mhz 300 Mhz 400 Mhz 500 Mhz 600 Mhz 700 Mhz 800 Mhz 900 Mhz
Amplitud de DPs (Normalizada)
3Kv
Tiempo (ms)
Figura 7-11.: Barrido en frecuencias 100 MHz - 900 MHz

8. Conclusiones y recomendaciones
8.1. Conclusiones
Se desarrollo un sistema de detección de DPs basado en un circuito que emula las DPs en
transformadores de potencia con técnicas de microondas usando software definido por radio
y una tarjeta de adquisición USRP B210 con resultados satisfactorios relacionando el costo
del sistema con los sistemas actuales que usan el mismo principio.
El estudio de espectro electromagnético de las descargas parciales en transformadores

de potencia usando software de simulación CST arroja resultados que ayudan a se-
leccionar el tipo de antena siendo la mas adecuada un monopolo con una geometría
suave que no provoque rupturas y/o nuevas descargas en la interacción de campo elec-
tromagnético provocada por la parte activa del transformador con la geometría de la
antena; en cuanto a la dirección de la misma se concluye que en los tres ejes X, Y y
Z se obtiene el suficiente registro de campo para lograr la detección, con respecto a
las frecuencias de interés la banda UHF es ideal para la detección; por último estas
simulaciones aportan una idea de lo que pruebas experimentales pueden resultar.
La tecnología para el análisis de espectro electromagnético, tiene características simila-

res a los sistemas usados en el estado del arte, con diferencias en el precio del sistema, el
uso del montaje en laboratorio de alta tensión replica el comportamiento de descargas
superficiales,corona e internas.
Los algoritmos de procesamiento desarrollados para la detección de descargas parciales

en la fase adquisición y detección resuelven el problema de una manera eficiente, tenien-
do en cuenta la cantidad de muestras por segundo requerida para calcular una señal
de alta velocidad como lo son las DPs. Es entonces el software definido por radio una
herramienta suficiente para detectar DPs sin embargo se sugieren trabajos posteriores
para caracterizar, clasificar y localizar las fuentes de descarga.
Se caracterizó el sistema haciendo uso del laboratorio de alta tensión de la Universidad

Nacional con el que se logra reproducir el comportamiento de DPs en transformadores
de potencia. Los resultados obtenidos reflejan el comportamiento antes calculado por
la simulación en CST en cada uno de los tipos de descarga (superficial, corona, interna)
este sistema es comparable con los sistemas antes realizados por otros autores con la
56 8 Conclusiones y recomendaciones
ventaja de un precio mas bajo y con tecnologías mas fáciles de adquirir aportando en
el ámbito académico un banco de pruebas de descargas parciales económico, confiable
y con resultados satisfactorios.
8.2. Recomendaciones y trabajos futuros

Se recomienda realizar trabajos futuros en el ámbito de caracterización, clasificación
y localización de DPs con los que se pueda calcular aproximadamente el estado de
envejecimiento del transformador.
Se recomienda agregar una nueva etapa en el algoritmo en la parte de calibración de

las tarjetas utilizadas, la calibración se puede realizar con un instrumento de medición
como un analizador de espectro, esto con el fin de hacer cálculos de potencia y registros
de campo eléctrico.
Se recomienda hacer uso del sistema de detección para generar una base de datos con la
que se puedan realizar investigaciones posteriores además la construcción de diferentes
bancos de prueba para generar mas tipos de descarga.
Se recomienda hacer mediciones reales en un transformador en servicio para una em-

presa de construcción de transformadores o prestadora del servicio de energía eléctrica.
8.3. Aportes del trabajo

Modelo de simulación de DPs
Se desarrollo un modelo de simulación usando el software CST, en que se diseño una
geometría CAD de un transformador de potencia, los materiales de la parte activa y
aislantes del transformador se tuvieron en cuenta. Las respuestas de la simulación se
obtuvieron en el dominio de la frecuencia y tiempo.
El modelo de simulación fue presentado como articulo de congreso "Full wave simu-
lation of partial discharge in power TransformersÏCEAA International Conference on
Electromagnetics in Advanced Applications.
Circuito aislante de alta tensión con hardware de adquisición y procesamiento de datos

- Hardware de adquisición de banda base usando tarjeta de sonido común.
Se diseño y desarrollo un circuito que separa físicamente la parte de alta tensión con
el hardware de adquisición, esto con el fin de proteger la electrónica de adquisición y
procesamiento de datos.
En la etapa de adquisición de la banda de 60 Hz se uso una tarjeta de sonido de
computador lo que reduce los gastos del proyecto y plantea nuevas perspectivas de
8.3 Aportes del trabajo 57
procesamiento de señales generadas por circuitos de alta tensión usando radio definido
por software y una tarjeta de sonido común.
Diseño de una antena que puede ser situada en el interior del tanque del transformador.
Se diseño una antena monopolo con punta en forma de gota que reduce el riesgo de
nuevas descargas al ser ingresada en el tanque del transformador, esta antena tiene
características muy similares a las antenas presentes en el mercado que suelen ser
monopolos. En el proceso de diseño de la antena se tuvo en cuenta los puntos de
acceso de la misma y se pretende fabricar para ser probada en la industria.
Software de adquisición, almacenamiento, procesamiento y detección de DPs

Se desarrollo un software para la captura de DPs usando radio definido por software,
reduciendo los gastos del proyecto y facilitando el acceso de la información y a futuras
investigaciones.
Circuito equivalente de DPs y objetos de ensayo.

Se implemento un circuito emula DPs en el laboratorio de alta tensión de la universidad
Nacional de Colombia, futuras investigaciones pueden partir de este montaje y sus
respectivos objetos de prueba.
Sistema de adquisición de DPs

Un sistema de adquisición de DPs fue desarrollado usando hardware económico (USRP
y RTL) basado en radio definido por software, además del uso de herramientas de
software libre para el desarrollo del algoritmo de procesamiento detección de DPs. Lo
que ofrece nuevas oportunidades de desarrolla la industria nacional y avance en la
investigación en este tema con un sistema económico.
A. Medición de antenas usadas
(Rohde y Schwarz Spectrum
Analyzer)
Figura A-1.: Antena Wifi

59
Figura A-2.: Antena JCG015

60 A Medición de antenas usadas (Rohde y Schwarz Spectrum Analyzer)
Figura A-3.: Antena QUAD GSM

Bibliografía
[Achillides et al., 2008] Achillides, Z., Georghiou, G., and Kyriakides, E. (2008). Partial dis-
charges and associated transients: The induced charge concept versus capacitive modeling.
IEEE Transactions on Dielectrics and Electrical Insulation, 15(6):1507–1516.
[Akbari et al., 2016] Akbari, A., Werle, P., Azirani, M. A., and Mirzaei, H. R. (2016). Cha-
llenges in calibration of the measurement of partial discharges at ultrahigh frequencies in
power transformers. IEEE Electrical Insulation Magazine, 32(2):27–34.
[Albarracín et al., 2013] Albarracín, R., Rajo-iglesias, E., and Martínez-tarifa, J. M. (2013).
Antenna Parametrization for the Detection of Partial Discharges Antennas parametriza-
tion for the detection of partial discharges. IEEE Transactions on Dielectrics and Electrical
Insulation, 62(May):932–941.
[Aranguren Z. and Tarantino A., 2009] Aranguren Z., S. and Tarantino A., R. (2009). Meto-
dologías y tecnologías de detección y diagnóstico de fallas aplicadas a procesos Industriales.
Revista Colombiana de Tecnologías de Avanzada.
[Autores, 2018] Autores, V. (2018). gnuradio. url: https://www.gnuradio.org/about/. Ac-

cedido 20-05-2018.
[Bartnikas, 2002] Bartnikas, R. (2002). Partial discharges their mechanism, detection and
measurement. IEEE Transactions on Dielectrics and Electrical Insulation, 9(5):763–808.
[Boggs and Stone, 1982] Boggs, S. A. and Stone, G. C. (1982). Fundamental Limitations
In The Measurement Of Corona And Partial Discharge. IEEE Transactions on Electrical
Insulation, EI-17(2):143–150.
[Callender et al., 2017] Callender, G., Rapisarda, P., and Lewin, P. L. (2017). Improving
models of partial discharge activity using simulation. 2017 IEEE Electrical Insulation
Conference (EIC), 3(June):392–395.
[Calva et al., 2010] Calva, A., Torres, D., and Ramírez, C. (2010). Sistema de Localización
de Descargas Parciales en Línea. PhD thesis, Instituto Politecnico Nacional.
[Cunha et al., 2015] Cunha, C. F., Carvalho, A. T., Petraglia, M. R., and Lima, A. C. (2015).
A new wavelet selection method for partial discharge denoising. Electric Power Systems
Research, 125:184–195.
62 Bibliografía
[Díaz Gautier and Martínez Tarifa, 2002] Díaz Gautier, J. R. and Martínez Tarifa, J. M.
(2002). Caracterización de defectos en sistemas de aislamiento mediante detección de
descargas parciales. PhD thesis, Universidad Carlos III de Madrid Repositorio.
[Fernández, 2013] Fernández, S. (2013). Evaluación Por Riesgo Del Peso De Las Fallas.
pages 1–11.
[Griffiths and Pilling, 2004] Griffiths, H. and Pilling, N. (2004). Advances in High Voltage
Engineering. IET Power and Energy Series 40, pages 349–413.
[Isermann, 2006] Isermann, R. (2006). Fault-Diagnosis Systems, pages 49–58.
[Jahangir et al., 2017] Jahangir, H., Akbari, A., Werle, P., Akbari, M., and Szczechowski,
J. (2017). UHF characteristics of different types of PD sources in power transformers.
Electrical Engineering (ICEE), 2017 Iranian Conference on, 1(Icee):1242–1247.
[Jiang et al., 2013] Jiang, T., Li, J., He, Z., and Wang, P. (2013). Simulation and experiments
of UHF detection for partial discharge in power transformers. Annual Report - Conference
on Electrical Insulation and Dielectric Phenomena, CEIDP, pages 1262–1265.
[Jin et al., 2012] Jin, C.-h., Lee, J.-y., Park, D.-w., and Kil, G.-s. (2012). Detection of Partial
Discharges by a Monopole Antenna in Insulation Oil. 11-th Conference on Instrumenta-
tion, Measurement, Circiuts and Systems (IMCAS ’12), pages 27–30.
[Judd, 2011] Judd, M. D. (2011). Experience with UHF partial discharge detection and
location in power transformers. Electrical Insulation Conference (EIC), 3(June):201–205.
[Kreuger, 1989] Kreuger, F. (1989). Partial discharge detection in high-voltage equipment.

Butterworths.
[Kreuger et al., 1993] Kreuger, F. H., Gulski, E., and Krivda, A. (1993). Classification of
partial discharges. IEEE Transactions on Electrical Insulation, 28(6):917–931.
[Li et al., 2012a] Li, T., Rong, M., Zheng, C., and Wang, X. (2012a). Development simula-
tion and experiment study on UHF partial discharge sensor in GIS. IEEE Transactions
on Dielectrics and Electrical Insulation, 19(4):1421–1430.
[Li et al., 2012b] Li, X.-C., Niu, D.-Y., Xu, L.-F., Jia, P.-Y., and Chang, Y.-Y. (2012b).
Simulation of transition from Townsend mode to glow discharge mode in a helium dielectric
barrier discharge at atmospheric pressure. Chinese Physics B, 21(7):075204.
[Mohamed et al., 2016] Mohamed, H., Lazaridis, P., Upton, D., Khan, U., Saeed, B., Jaber,
A., Zhang, Y., Mather, P., Vieira, M. F., Barlee, K. W., Atkinson, D. S., Mihovska,
A., Gavrilovska, L., and Glover, I. A. (2016). Partial discharge detection using low cost
RTL-SDR model for wideband spectrum sensing. 2016 23rd International Conference on
Telecommunications, ICT 2016, 1(November 2017).
Bibliografía 63
[Nanao et al., 2008] Nanao, K., Murakami, Y., and Nagao, M. (2008). Analysis of internal
partial discharge based on PD current waveform. Annual Report - Conference on Electrical
Insulation and Dielectric Phenomena, CEIDP, pages 475–478.
[Negm et al., 2017] Negm, T. S., Refaey, M., and Hossam-Eldin, A. A. (2017). Modeling
and simulation of internal Partial Discharges in solid dielectrics under variable applied
frequencies. 2016 18th International Middle-East Power Systems Conference, MEPCON
2016 - Proceedings, pages 639–644.
[Raja et al., 2002] Raja, K., Devaux, F., and Lelaidier, S. (2002). Recognition of discharge
sources using UHF PD signatures. IEEE Electrical Insulation Magazine, 18(5):8–14.
[Robles et al., 2013] Robles, G., Sanchez-Fernandez, M., Albarracin Sanchez, R., Rojas-
Moreno, M. V., Rajo-Iglesias, E., and Martinez-Tarifa, J. M. (2013). Antenna parametri-
zation for the detection of partial discharges. IEEE Transactions on Instrumentation and
Measurement, 62(5):932–941.
[Shibuya et al., 2010] Shibuya, Y., Matsumoto, S., Tanaka, M., Muto, H., and Kaneda, Y.
(2010). Electromagnetic waves from partial discharges and their detection using patch
antenna. IEEE Transactions on Dielectrics and Electrical Insulation, 17(3):862–871.
[Siegel et al., 2017] Siegel, M., Beltle, M., and Tenbohlen, S. (2017). Application of UHF
Sensors for PD Measurement at Power Transformers. IEEE Transactions on Dielectrics
and Electrical Insulation Vol., 24(1):331–339.
[Sinaga et al., 2012] Sinaga, H. H., Phung, B. T., and Blackburn, T. R. (2012). Partial
discharge localization in transformers using UHF detection method. IEEE Transactions
on Dielectrics and Electrical Insulation, 19(6):1891–1900.
[Standard, 2001] Standard, B. (2001). High-voltage test techniques — Partial discharge

measurements.
[Standard, 2012] Standard, B. (2012). High Voltage Test Techniques - Measurement of Par-
tial Discharges by Electromagnetic and Acoustic Methods.
[Tozzi et al., 2010] Tozzi, M., Eng, P., and Carlo, G. (2010). Partial discharge in power
distribution electrical systems: pulse propagation models and detection optimization. PhD
thesis, University of Bologna.
[Vanegas, 2016] Vanegas, J. A. (2016). Desarrollo de un sistema prototipo de medición y

análisis de DP generadas en laboratorio. PhD thesis, Universidad del Valle.
[Wu et al., 2015] Wu, M., Cao, H., Cao, J., Nguyen, H. L., Gomes, J. B., and Krishnaswamy,
S. P. (2015). An overview of state-of-the-art partial discharge analysis techniques for
condition monitoring. IEEE Electrical Insulation Magazine, 31(6):22–35.
64 Bibliografía
[Xu et al., 2015] Xu, B., Li, Y., Zhao, J., Wang, X., and Zhang, R. (2015). Simulation of
propagation characteristic of ultra high frequency signals in power transformers using finite
difference time domain method. In 2015 5th International Conference on Electric Utility
Deregulation and Restructuring and Power Technologies (DRPT), pages 1591–1594. IEEE.
[Zhou et al., 2018] Zhou, Y., Zhou, Z., Zhang, L., Zhang, Y., Mo, Y., and Sun, J. (2018).
Characterization and comprehension of corona partial discharge in air under power fre-
quency to very low frequency voltage.

Sistema de Detección DPs Técnicas Microondas PDF

Cargado por

Copyright:

Formatos disponibles

Sistema de Detección DPs Técnicas Microondas PDF

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Sistema de Detección DPs Técnicas Microondas PDF

Cargado por

Copyright:

Formatos disponibles

Sistema de detección de descargas

parciales en transformadores de potencia

David Leonardo Galindo Huertas

Universidad Nacional de Colombia

David Leonardo Galindo Huertas

Tesis presentada como requisito parcial para optar al título de:

Universidad Nacional de Colombia

Esta investigación presenta los resultados de un sistema de detección de descargas parciales

4-1. Diagrama de flujo para simulación . . . . . . . . . . . . . . . . . . . . . . . . 22

5-1. Ingreso de antenas al transformador [Siegel et al., 2017] . . . . . . . . . . . . 29

5-4. USRP B210 y Arquitectura USRP B210 . . . . . . . . . . . . . . . . . . . . 32

6-1. Diagrama de flujo Visión general del procesamiento de datos . . . . . . . . . 39

7-1. Encendido eléctrico cocina . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

A-1. Antena Wifi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

2. Marco teórico y estado del arte 4

3. Planteamiento del problema y objetivos 20

4. Espectro electromagnético producido por DPs en transformadores de potencia. 22

5. Tecnología análisis de espectro electromagnético y montaje para generación

6. Algoritmo de adquisición, procesamiento y detección de DPs 38

7. Pruebas experimentales y Resultados 45

A. Medición de antenas usadas (Rohde y Schwarz Spectrum Analyzer) 58

2.1. Clasificación y espectro electromagnético de las

2.1.2. Proceso de generación de DPs

Figura 2-1.: Vector de fallo DPs.

Aparición de micro espacios en el material aislante(t=0).

Determinación Electrones disponibles (t=1).

Townsend Discharge generados por los coeficientes de ionización (t=2).

Figura 2-3.: Cantidad de electrones ionizados por probabilidad de ionización

Por definición; si es la cantidad de electrones llegando a la posición x y siendo α la

Si α es independiente de x entonces la cantidad de electrones cruzando por la nueva

Figura 2-4.: Interacción de electrones y iones en el dipolo

N0 es la cantidad de electrones por segundo, en un tiempo t > 0 se alcanzara una

(a) Proceso de generación de mul- (b) Estado de la descarga

Figura 2-5.: Proceso de generación de múltiples avalanchas y corriente de la DPs

i = i0 · eα·d /1 − γ(eα·d − 1) (2-5)

Con el tiempo y debido al fenómeno repetitivo de las DPs se producen cambios en

Cambios químicos, aparición del fallo(t=4)

Por último, se produce un cambio químico debido al incremento de temperatura en el

2.1.3. Clasificación de DPs según su origen

Figura 2-6.: Descargas superficiales

Figura 2-7.: Descarga interna

Figura 2-8.: Campo Ea distribuido y campo de la vacuola Ec

En cuanto mayor sea la permitividad dentro de la imperfección o vacuola ε mayor será el

2.1.4. Métodos de medida de descargas parciales

Discharge Effects IEC 60270

Chemical Effects Chemical

Figura 2-9.: Métodos de medida de Descargas parciales

2.2. Modelos de DPs

Figura 2-10.: Circuito equivalente

la capacitancia de la imperfección. Cada vez que se supera la tensión de ignición se presenta

Figura 2-12.: Corriente de DPs

2.2.2. Modelo Dipolo

Figura 2-13.: Modelo de dipolo

te pulso es representado por dos parámetros: amplitud en mA y (FWHM) anchura a media

Ip Tipico pulso gaussiano de DPs

Figura 2-14.: Pulso Gaussiano

La onda electromagnética generada en el tiempo y basada en el modelo de dipolo eléctrico

Figura 2-15.: Dipolo eléctrico y punto de observación