NUEVO FACTOR DE DISTORSIÓN TOTAL POR ARMÓNICAS (THD) PARA

SEÑALES APERIODICAS
DISTORSIÓN ARMONICA DE VENTANA MOVIL (STHD)
Leonidas Sayas P. Msc, Ing.
[email protected]

RESUMEN: denominadas no lineales, siendo el responsable

principal de estas deformaciones, el desarrollo
En la evaluación del sistema eléctrico de potencia espectacular de los semiconductores y la electrónica
(SEP) ante la presencia de dispositivos no lineales, de potencia, desde principios de los sesenta, que ha
se supone que las tensiones y corrientes son generalizado el uso de rectificadores, tanto en
periódicas pero no sinusoidales, y se analizan aparatos de potencia como en otros pequeños pero
utilizando herramientas, como la transformada de muy numerosos tales como ordenadores, televisores,
Fourier. VHS, DVD, iluminación controlada etc. A causa del
Lamentablemente, existen cargas industriales que incremento de estas cargas no lineales, la polución
consumen una intensidad de corriente que no es armónica de las tensiones y corrientes han ido
periódica. aumentando produciendo un mal funcionamiento de
Así en estas situaciones donde las señales no son los distintos elementos ya sea aguas abajo y aguas
periódicas, y por tanto no existe su desarrollo en arriba del punto de acoplamiento común (PAC) del
serie de Fourier, herramientas como la transformada SEP[1].
de Fourier de ventana móvil o la transformada La problemática de la existencia de armónicas en el
wavelet permiten analizar el espectro de frecuencia SEP ha derivado en el estudio de la calidad en la
de la señal aperiodica. onda de tensión y onda de corriente, cuyo objetivo es
definir índices o factores y normativas que
En este trabajo se desarrolla un nuevo índice a partir cuantifiquen y limiten la distorsión de las ondas a
de la transformada de Fourier de ventana móvil niveles tolerables para el sistema y los receptores.
(STHD Shot-Term Harmonic Distortion), y con el El control de esta calidad supone el análisis previo
objetivo de mejorar la caracterización de la calidad del espectro de las mismas para posteriormente
del producto en régimen no periódico, se propone comparar los resultados obtenidos con los índices o
una extensión de los factores habitualmente factores que describen dicha calidad.
utilizados en condiciones periódicas. Este nuevo
conjunto de indicadores para analizar señales no
periódicas convergen a los índices habituales
utilizados en condiciones periódicas y además (h)
J NL = PI( h ) + PG( h ) + PLh ......(2)
proporcionan resultados consistentes a partir de
algoritmos estándar y permiten, una interpretación 1.2 Flujo de carga armónica :
física de los resultados. Los métodos de análisis armónicos de redes
eléctricas sirven para determinar como afecta la
presencia de cargas no lineales a los valores de las
1. INTRODUCCIÓN magnitudes del sistema (tensiones, corrientes,
1.1 La polución armónica: potencia etc.)
,A la frecuencia fundamental el circuito se comporta
Si bien las corrientes y tensiones en el sistema con su equivalente thevenin alimentando a una carga
eléctrico de potencia nunca han sido exactamente no lineal y a frecuencia armónica el sistema se
senoidales, hasta hace algunos años no era representa a la carga no lineal como una fuente de
frecuente encontrarse en la practica con algún corriente.
problema originado por la deformación de dichas
ondas o lo que es lo mismo, por la presencia de
armónicas, que son originados por las cargas
 VG ZL VL

2. CALIDAD DEL PRODUCTO

PI 2.1 Consideraciones:
El problema del incremento de la distorsión, ha
NL
VG derivado, entre otras cosas, en la definición de
índices y normativas cuya finalidad es cuantificar y
PNL limitar la distorsión de las ondas eléctricas a nivel del
SEP y que los receptores puedan tolerar. La
normativa existente supone el análisis previo del
espectro frecuencial de las señales del sistema, con
el fin de determinar su distorsión. Es decir, si la señal
Fig.1.1 Circuito a frecuencia fundamental
es periódica, la Transformada de Fourier (FT) para
señales continuas, la transformada discreta de
(1)
PNL = PU(1) + PP(1) + PH(1) .......(1) Fourier (DFT) para señales discretas o la
transformada rápida de Fourier (FFT) para señales
A frecuencia fundamental, la PG generada se discretas con un numero de muestras potencia de
absorbe principalmente en la carga no lineal PNL y dos, serían las herramientas a utilizar. En este caso
una pequeña parte se consume en al línea PI. La de señales aperiodicas, existen herramientas como
tensión fundamental en los bornes del generador VG la transformada de Fourier de Ventana Móvil
es mayor a la tensión en la línea VL. (STDFT, Short Time Discrete Fourier Transform) o el
Una parte de la potencia absorbida por la carga no análisis mediante la transformada wavelet que
lineal (PNL) se convierte en potencia útil (PU) y en permite analizar el espectro frecuencial de la señal.
pérdidas de potencia activa en la misma (PP), Para el caso de señales periódicas existe numerosos
índices de calidad las mismas que se encuentran en
mientras que el resto (PH) se convierte en corrientes
las normas IEC, IEEE, CIGRE, CENELEC, NTCSE
armónicas que son inyectadas por la carga no lineal etc. Estos índices están basados en el desarrollo de
a la red. Fourier de la señal, por lo que dejan de ser útiles
VGh ZL VLh para señales a periódicas.
Lamentablemente, existen muchas cargas que no
consumen corrientes periódicas. Entre ellas se
PIh encuentran los hornos de arco eléctricos, las
maquinas de soldar eléctrica, laminadoras de acero,
PGh y ciertos variadores de velocidad para esquemas de
ZG
Jh control especializados. Existen otros casos en las
cuales la corriente de carga es no periódica: la carga
puede variar, puede estar bajo controles que causen
PLh una variación de ella, o la tensión que produce puede
estar sufriendo variaciones no periódicas. Debido a
que pocos procesos son realmente estacionarios, las
variaciones en las corrientes de carga y la aparición
Fig.1.2 Circuito a frecuencia armónica de ondas a periódicas son sucesos corrientes.
A una de las frecuencias armónicas, se puede considerar 2.2 Transformadas frecuenciales e índices de
que la carga no lineal se comporta como una fuente de calidad:
corriente de valor Jh cuyo valor depende de la tensión en
sus propios bornes, mientras que el generador aparece La aplicación de los distintos índices de calidad
solo con su impedancia interna. La intensidad armónica
2π
inyectada por la carga no lineal provoca caídas de tensión 1 N −1
∑
jk Tsn
en la línea y en el generador. g ´n = G´ k e T
(n = 0 ÷ N - 1)...(9)
N k=0
Por lo tanto las tensiones del generador y de la carga
están distorsionadas. Además, que las corrientes supones, la transformación de las señales
temporales a señales en el dominio de la frecuencia
armónicas Jh inyectadas por la carga no lineal distinguiendo en ambos caso si la señal es periódica
dependen del valor y de la forma de onda de la o no. Es preciso abordar estos conceptos antes de
tensión en bornes de dicha carga; esta dependencia iniciar con el desarrollo del tema.
se conoce como INTERACCIÓN ARMÓNICA [5]
2.2.1 Transformadas en el dominio de la
N −1 2π
frecuencia − jk
∑
Tsn
G ´k = g´ n e
T
(k = 0 ÷ N - 1)...(10)
En general, cualquier función continua, g(t), se puede n=0

expresar a partir de la ecuación de la Transformada A diferencia lo que ocurría con la serie de Fourier en
de Fourier (FT), tiempo continuo, en la serie de Fourier en tiempo
discreto sólo existen N funciones sinusoidales de
1 ∞ frecuencia múltiplo de la frecuencia fundamental,
g (t ) =
2π ∫
−∞
G ( w)e jwt dw (w ∈ ℜ) ...(3)
f1=1/T, dado que son periódicas. Esta frecuencia
fundamental sería la resolución en frecuencia.
Donde la obtención de la función compleja G(w) se
realiza a partir de la FT de la función continua, Debido a la naturaleza finita y discreta de los
∞ ordenadores, cualquier señal que sea procesada
G ( w) = FT {g (t )} = ∫ g (t )e − jwt dw (w ∈ ℜ) ...(4) tendrá que ser discreta y de función finita, al igual
−∞
que las transformadas correspondientes. Así
cualquier señal de tiempo continuo, g(t), podrá ser
Si la función continua estudiada es periódica, g´(t),
muestreada para su procesamiento de forma que se
de periodo T, se puede expresar a partir de su serie obtenga la secuencia discreta gn=g(nTs) con n=0,..,N-
de Fourier, 1, siendo N el numero de muestras, Ts el intervalo de
Donde w1=2πf1=2π/T es la pulsación fundamental de muestreo y T=NTs la duración de la señal.
la señal y f1 su frecuencia fundamental. La secuencia discreta así obtenida puede ser
La serie de Fourier no es mas que la representación procesada para obtener su respuesta frecuencial a
de una función continua y periódica, g´(t), como partir de la ecuación (10);esta respuesta será la
suma ponderada de senos de frecuencia múltiplo de transformada discreta de Fourier, DFT, de la
la frecuencia fundamental, f1, de la señal. secuencia gn, y corresponde a la respuesta
1 ∞ frecuencial de la secuencia g´n que se formaría por
g ´( t ) = ∑ G k e jkw 1 t ...(5)
T k =-∞ repetición periódica de gn,
La obtención de los coeficientes espectrales de N −1 2π
− jk
G k = DFT {g n }∑ g n e
Ts n
T
ponderación Gk se realiza con la ecuación de ...(11)
análisis, n =0
para el mismo intervalo, Se observa que los términos
∫ Gk así obtenidos corresponden a una frecuencia
− jkw
G k = g ´( t ) e 1t
dt ...(6)
T fk=k/T=kf1, donde f1=1/T es la resolución en el
En general cualquier función discreta, gn, se puede espectro frecuencial, o frecuencia mínima del calculo
expresar a partir de la ecuación de síntesis de la realizado.
transformada de Fourier, El cálculo de DTF considera el tratamiento de
señales discretas de larga duración y por lo tanto
1 ∞
2π ∫2π
gn = G (Ω)e jΩn dΩ ...(7) requiere procedimientos computacionales eficientes.
En este sentido si el número de muestras de la señal
donde la obtención de la función compleja G(Ω) se discreta tiene un número de muestras potencias de
realiza a partir de la transformada de Fourier en el dos se emplea la transformada rápida de Fourier
tiempo discreto, FFT.
∞
Si la señal de tiempo continuo, g(t),presenta un
G (Ω ) = ∑
n =-∞
g ne − jΩ n
...(8)
contenido armónico variable con el tiempo, se recurre
a la transformada de Fourier de ventana móvil, STFT,
siendo Ω=wTs Ts es el intervalo de muestreo de la
si se trabaja en el tiempo continuo, o a la
señal continua. A diferencia de la transformada de
tiempo continuo, G(w), la transformada en el tiempo transformada discreta de Fourier de ventana móvil,
STDFT, si se trabaja en tiempo discreto a partir del
discreto, G(Ω), es periódica de periodo 2π.
muestreo de la señal g(t).
Si la función discreta estudiada es periódica, g´n, de
periodo N ó T=NTs donde Ts es el intervalo de
muestreo, se puede expresar a partir de su serie de
Fourier,
La obtención de los coeficientes espectrales de
ponderación G´k se realiza con la ecuación de
análisis,
 La función ventana utilizada W(t-t´) es la función
ventana rectangular de amplitud unitaria mostrada en
la Figura 2.1, siendo ΔT el ancho de dicha ventana.
Fuera del intervalo correspondiente el ancho ΔT la
función ventana tiene un valor igual a cero. Como
consecuencia, el producto de la señal y la función
ventana en el intervalo donde el valor de esta ultima
no es nulo, es igual a la señal original.
n-0 .... n -N-1
Fig. 2.2Muestreo de señales para su procesamiento.
Aplicación de una función ventana rectangular En el caso de señales discretas, se trabaja a partir de
una ventana de ancho N muestras que se mueve en
el tiempo ,según Figura 2.2 De este modo, la ventana
La idea básica de la transformada STFT consiste en contiene una secuencia de muestras obtenidas a
multiplicar partir de la señal g(t), que para el caso de la ventana
de orden m se podría expresar como gn,m = g((m-
1)∆T+nTs).w(nTs)con n=0 .., N-1m= 0, 1, ..., siendo
la señal a analizar g(t) por una función ventana cuya Ts el tiempo de muestreo y ∆T = N Ts la duración de
posición varia en el tiempo, es decir, la señal se la ventana.
divide en pequeños segmentos de corta duración, Así se puede aplicar la DFT al conjunto de N
Δ T, Fig. 2.1 y aplicar la transformada de Fourier en muestras de la ventana m-ésima y la definición de la
cada uno de dichos segmentos. De esta manera, se STDFT resulta,
obtiene el contenido frecuencial de la señal durante
2π
un periodo corto de tiempo. La totalidad del espectro N −1 − jk
∑g
t .n
resultante contiene la evolución en el dominio de la Gk,m = STDFT{gn,m} = e
n, m
T
(K = 0 ÷ N − 1)...(13)
frecuencia de la señal, y por tanto, la STFT tiene en n =0

cuenta la evolución temporal de la señal, reflejada en

dicho espectro [7, 8]. Al igual que sucedía en la expresión (11) los
términos Gk,m corresponden a una frecuencia fk =
g(t) k/∆T =Kfmin donde fmin= 1/∆T seria la resolución en
el espectro frecuencial, o frecuencia mínima, del
SEÑAL
calculo realizado.
t
Tal como se ha comentado al caracterizar la STFT, la
función ventana utilizada es la función ventana
w ( t - t´) rectangular de amplitud unitaria, es decir:
1.0
VENTANA
t´ t ⎛1n = 0,.., N − 1 ⎞
ΔΤ wn = w( nT s ) = ⎜⎜ ⎟⎟
⎝ 0casocontra rio ⎠
g(t) . w ( t - t´)
La figura 2.3 muestra una manera de representar el
FRACCION calculo realizado por la transformada STDF de una
DE SEÑAL señal g(t).
t
La señal se divide en ventanas de duración ∆T
(resolución en el dominio del tiempo) tomándose en
cada ventana N muestras de dicha señal. Sobre
Figura 2.1.Aplicación de una función ventana a una señal esas, muestras, y para cada ventana, se realiza la
DFT obteniéndose el espectro de la señal para unas
Así, la definición de la STFT es, determinadas frecuencias. El espectro así obtenido
esta constituido por N valores y la resolución en el
G ( w, t´) = STFT {g (t )} = dominio de la frecuencia de dicho espectro, o
∞ frecuencia mínima, es fmin, el cual viene fijado por la
∫−∞
g (t ).w(t − t´).e − jwt dt ( w ∈ R ).......(12) duración temporal de la ventana según fmin = 1/∆T.
 - El primer valor del espectro frecuencial, de
Ventana m-ésima
f V1 Vm cualquier ventana, Go,m, corresponde al valor
de continua fo=0fmin=0.
(N-1)f min
- La posición de la frecuencia fundamental de
la señal, f1, en el eje de las frecuencias se
Valor G k,m
obtendría según k1 = int( f 1 / f min) + 1 ,
K f min f min
siendo la función int( f 1 / f min) el redondeo
Frecuencia k del cociente f 1 / f min al entero mas
próximo.
0

La función ventana rectangular se eligió por su

sencillez y por su correlación con la transformada
g (t)
de Fourier. Sin embargo, otras funciones
N muestras ventana, como la función de Gauss o la función
de Hanning, podrían ser utilizadas en la
realización de otros estudios, tomando como
partida lo expuesto.

0 ........ t 2.2.2 Índices de distorsión total por armónicas

m T
T
Fig.2.3 Representación de STDFT de una señal g(t)
Existe un amplio conjunto de índices y medidas para
describir la calidad de potencia eléctrica de señales
periódicas, pero ninguno describe de forma univoca a
Así, las señales por lo que no son aplicables en todos los
- Para una ventana m se obtendrían Gk,m casos y, en muchas ocasiones, estos índices pueden
valores, con k = 0 ÷ N − 1 , y la intersección ocultar información en lugar de mostraría.
marcada en dibujo captaría la componente
Gk,m que corresponde a la frecuencia Para tensiones y corrientes periódicas se ha
desarrollado un amplio conjunto de índices de
fk=kfmin y a un instante de tiempo t m =m∆T. calidad de potencia, los más comunes están
- Los resultados de DFT, se decir los valores definidos para señales que se pueden descomponer
Gk,m, serán distintos para cada ventana. según series de Fourier, y su definición se expresa a
- El tamaño de la ventana, ∆T, influirá en los partir de los términos fundamentales y armónicos de
resultados obtenidos ya que fija la frecuencia la señal. En referencia a los índices de señales
mínima fmin y por lo tanto las frecuencias del periódicas, ya han sido históricamente estudiados y
espectro estudiado, fk=kfmin . por lo tanto efectuaremos un breve resumen, tanto
- Una ventana de ancho pequeño proporciona para señales periódicas y aperiódicas.
una frecuencia mínima, fmin, de valor
elevado por lo que, aunque la resolución en 2.2.2.1 Señales periódicas
e tiempo se buena, la resolución en el Tal como se ha visto, toda onda distorsionada que
dominio de la frecuencia no lo es, y sea periódica, de periódico T, puede descomponerse
viceversa. en suma de ondas sinusoidales de frecuencia
- Si se desea obtener la máxima información múltiplo de la frecuencia fundamental f1=1/T según la
de la señal se deberá trabajar con ec. (5) a partir de ec. (6) si la señal es continua o
frecuencias mínimas, fmin, pequeñas, de según (9) a partir de (10) si la señal es discreta. Se
forma que se estudien la mayor parte de las denomina componente fundamental a la componente
frecuencias, y esto significa trabajar con de frecuencia f1 y armónicos a las frecuencias que
ventanas, ∆T, grandes. Sin embargo, el son múltiplos de este. Así, los armónicos describen
tamaño de las ventanas vendrá limitado por en forma matemática, en régimen permanente, la
la capacidad de cálculo que se disponga y distorsión de las ondas de tensión y corriente que
por la longitud de la muestra con la que se permite abordar el estudio de los sistemas de
pueda trabajar. potencia en presencia de cargas no lineales.
Para las señales periódicas, los índices mas
extendidos para cuantificar la distorsión armónica de
la señal son el índice de distorsión armónica total aplicados. Entre estos índices se encuentran los
(TDH, Total Harmonic Distortion ) y el índice de siguientes:
distorsión armónica individual (HD, Harmonic - Factor de cresta (Crest Factor) tensión o
Distortion), que se define, según el IEEE, a partir de corriente de pico dividido por 2 veces su
los componentes fundamental y armónicas de la
valor eficaz. Se utiliza en la evaluación de
señala a estudiar, por ejemplo la tensión v(t) cuyos
fatiga de los dieléctricos, y en la evaluación
armónicos se caracterizarían con V1 , Vh (h
de la fatiga de componentes electrónicos de
=2,...., ∞ )
potencia.

∑ (V )
∞
h 2 - Valor de pico: valor máximo positivo mínimo
V h negativo de la tensión o la corriente. Se
h =.2
THD = 1
...(14) y HD h = ...(15) utiliza en los mismos ámbitos que el anterior.
V V1 - Valor eficaz de la onda: Aplicable al calculo
de la amplitud total de señales periódicas y
La IEC llama a la expresión (14) factor de distorsión
(DF; distorsión Factor). Ambos índices presentan la algunas aperiódicas.
ventaja de que son fáciles y rápidos de calcular. Sin - Ancho de banda: Ancho (en Hz. o rad/s) del
espectro para el cual la magnitud de dicho
embargo, en su definición se puede perder
información del espectro total de la señal ya que no espectro cae 3 dB por debajo del valor de
describe de forma unívoca la señal. pico del espectro. Tiene aplicaciones en
filtros, canales de comunicación,....
También la NTCSE Peruana, denomina como factor
de distorsión total por armónicos y tampoco describe - Factor de potencia: Potencia media
la señal, expresado por: instantánea dividida por el producto del valor
eficaz de la tensión por el valor eficaz de la
corriente. Puede ser aplicado en casos
THD = ( ∑ i = 2..40
(Vi2 / V N2 ) × 100% ...(16) periódicos debido a que la potencia media y
los valores eficaces de tensión y corriente
son calculables.
Una valor comúnmente utilizado como frontera entre
alta y baja distorsión es el correspondiente a un THD
3. DISTORSIÓN ARMÓNICA DE VENTANA MÓVIL:
igual al 5%. Sin embargo, este valor es demasiado
alto para ser permitido en sistemas de transmisión,
pero para algunos sistemas de distribución es un El punto de partida para la definición de un nuevo
valor adecuado. índice que pueda ser aplicable para la medida de la
distorsión de una señal aperiódica es el concepto de
La tabla 2.1 muestra los limites de la distorsión de la la STDFT, (13). Esta herramienta matemática se
tensión dados por la normativa IEEE Stándard 519- utiliza para extender el uso de índices de calidad
1992 [9 existentes al caso de tensiones y corrientes
Distorsión de Tensión Distorsión de la tensión
aperiódicas. Por lo tanto, al realizar este estudio se
Tensión considero que el nuevo índice debería cumplir lo
individual (%) total
siguiente:
Hasta 69 KV 3.0 5.0 • La definición del nuevo índice debe
69.001 KV a 161 KV 1.5 2.5 converger a la definición habitual en el caso
Más de 161 KV 1.0 1.5 de que la señal sea periódica.
Li i d l di ió d l ió ú l i id d IEEE • Su calculo debe ser realizado mediante un
Tabla 2.1 Limites de la distorsión de la tensión según IEEE-519 algoritmo estándar que proporciones
resultados consistentes.
La NTCSE fija como limites de 3% para alta y muy • La definición del nuevo índice debe tener
alta tensión y 8% para sistemas de media y baja correspondientes, en la medida de lo posible,
tensión.[10] con procesos físicos como son las perdidas
Ninguno de los índices citados proporciona del sistema, fatiga, o desviación de la
información acerca de la amplitud de la onda. operación sinusoidal.
El objetivo de esta definición es, como se ha
comentado previamente, obtener un índice que, para
el caso de señales sinusoidales en la que se
2.2.2.2 Señales aperiódicas encuentran presentes componentes de frecuencia
En el caso de tensiones y corrientes aperiódicas, múltiplos enteros de la frecuencia fundamental f1,
existen una serie de índices que pueden ser converja hacia el índice de distorsión armónica THD.
Además, el nuevo índice debería ser consistente con componente continua debería iniciarse en
el concepto del teorema de Parseval [6], el cual K=0.
indica que la energía contenida en una señal que se
puede descomponer en una serie de Fourier es
función de los términos de dicha serie.
La Figura 3.1, tal como lo hacia la Figura 2.3,
3.1 Caracterización: muestra una manera de representar el calculo del
STHD.
La primera idea es utilizar la transformada STDFT de
un modo similar a como se utiliza la DFT para Ventana m-ésima
señales periódicas al definir el THD, expresión (15 o f V1 Vm
16). De este modo, el nuevo índice contendrá
información de las diferentes componentes (N-1)fmin
frecuenciales presentes en la señal, no solo de los
componentes frecuenciales presentes en la seña, no .
.
Gkm = STDFTkm
solo de los componentes armónicos, y podrá ser .

usado tanto para caracterizar señales periódicas K fmin f min

como aperiódicas.
.
Frecuencia k
.
Dada una señal, se define distorsión armónica de .

Ventana Móvil (STHD, Short Term Harmonic

Distortion), para cada una de las m ventanas para las 0
que se ha calculado la STDFT según (13), como
STHDo .............. STHDm ................

N −1 Fig. 3.1 Representación del cálculo del STHD

∑
k =1
STDFT 2
k ,m

k ≠ k1 A partir de ella se puede concluir

STDH m = ...( 17 ) - una vez realizada la STDFT sobre cada
STDFT K 1 ,m ventana, los elementos obtenidos como
resultado, por ejemplo los Gk.m con K= 1 ÷
N-1 para la ventana m-ésima, se
Donde, combinaran siguiendo la definición del TDH,
- N es el numero de frecuencias para las que (14 ó 16), para obtener el correspondiente
se ha calculado la STDFT. índice STDHm asociado a dicha ventana.
- K 1 = int( f 1 / f min) es la posición de la - Así, para cada ventana se obtendrá una
frecuencia fundamental de la señal, f1, en el valor de STDHm de forma que es
espectro de frecuencias calculado. conveniente representar el STHD de una
- fmin es la frecuencia mínima del espectro en señal como una función del tiempo (este
el dominio de la frecuencia. Dicha frecuencia aspecto, por ejemplo, la variación del THD
caracterizará la resolución del estudio con el tiempo, es un concepto familiar para
frecuencial. los ingenieros en el área de la calidad de
- Se ha utilizado la expresión STDFTk,m para suministro eléctrico).
renombrar a los Gk.m componentes del - La precisión en el cálculo del índice STHD
espectro. depende del tamaño de la ventana. Así, a
- Comparando las expresiones (14 o 16) y (17) medida que se incrementa su tamaño, ∆T,
se observa que la expresión (17) es una disminuye la frecuencia mínima, fmin por lo
generalización de la (14 o 16) que engloba que se captura el diverso contenido
todos los posibles componentes frecuencia frecuencial de la señal de una manera más
les que pueden salir al estudiar una señal, en precisa y se aporta mayor información al
este caso aperiódica. Por este motivo, la cálculo del STHD: esta idea se muestra en la
sumatoria de la expresión (17) recorre todas Figura 3.2. A partir de un determinado
las componentes frecuenciales (K= 1, ..., N- umbral, el STHD deja de incrementarse ya
1) a excepción de la fundamental (k= K1) y la que toda la información de las frecuencias
componente continua (K=0) de forma similar presentes en la señal han sido capturadas y
a la definición del TDH (14 o 16). Si en la una resolución, o frecuencia mínima, mas
expresión se hubiese considerado la
 fina no proporciona información adicional de obtenido, es decir 0.2, pues en este caso los índices
la señal. En la figura 3.2 (fmin)lim puede STHD y THD deberían coincidir.
estar limitado por la capacidad de calculo de
la que se disponga y por la longitud de la
muestra de la señal con la que se trabaje. f (Hz) V1 V2 V3 V4 V5 V6

- La transformada STDFT ofrece la posibilidad 1500 0.226 0.227 0.226 0.226 0.227 0.226

(a partir del trabajo con una fmin suficiente ... ...

pequeña) de capturar componentes de la 1080 0.095 0.076 0.095 0.095 0.076 0.095

señal que no son “verdaderos” armónicos y 1020 16.881 16.888 16.881 16.888 16.881 16.888

aportarlos al nuevo índice de calidad 960 0.487 0.491 0.487 0.487 0.491 0.487

definido. Esto es importante ya que los ... ...

múltiplos no enteros de la frecuencia 120 0.222 0.224 0.222 0.222 0.224 0.222

fundamental, interarmónicos y fracciones de 60 83.7 83.701 83.7 83.7 83.701 83.7

armónicos, han dado históricamente a 0 0.4 0.401 0.4 0.401 0.401 0.4
0.2027 0.2028 0.2027 0.2028 0.2028 0.2027
ingenieros e investigadores problemas en el STHDm

uso de los índices de calidad tradicionales. Tabla 3.1 STDFT de una señal

STHD I/(fmin)lim

4. CONCLUSIONES
Las conclusiones mas relevantes al termino de este
trabajo son las siguientes:

ƒ Definición de un nuevo índice de calidad, a

STHD se incrementa partir del concepto de transformada de
con el ancho de la
ventana
Fourier de ventana móvil, con el objetivo de
mejorar la caracterización de la calidad del
servicio eléctrico en régimen no periódico. La
Ancho de ventana T definición se realizó a partir de la extensión
de la distorsión armónica total.
Fig. 3.2 Variación del STHD en función de ΔT para una señal
aperiodica
ƒ Se verifico que estos nuevos índices
convergen a los índices en regímenes
periódicos.
3.2 Aplicación:
ƒ Una extensión de este análisis se puede
Estudio de una señal como un caso ilustrativo, efectuar para los índices conocidos como
definida por: son, Interferencias en la señal telefónica,
TIF, Perdidas de potencia activa, Perdida de
g(t) = cos (2 π 60t) +0.2.cos (2 π I020t) vida útil de componentes eléctricos o la fatiga
de dieléctricos, los que están directamente
Donde la frecuencia fundamental se ha tomado relacionados con la interferencia de la señal,
f1=60 Hz. así mismo el factor de potencia, el factor K
La señal considerada presenta un único armónico a que relaciona las perdidas en
1020 Hz (h=17) de valor 0.2 por lo que se podría transformadores en el caso de corrientes de
calcular el THD cuyo valor seria 0.2. carga no sinusoidales.
Para esta señal se calculo la STDFT, que se muestra
en la Tabla 3.1, siguiendo el patrón de la Figura 3.1.
Las filas de la Tabla 3.1 corresponden a diferentes ƒ Del mismo modo esta extensión se ha verificado
frecuencias espaciadas con una frecuencia fmin= en la evaluación del conductor neutro en función
de la distorsión de la onda de corriente frente a
60Hz. las columnas corresponden a ventanas una señal a periódica, presentando una variación
espaciadas un intervalo de tiempo ∆T = 16.67 ms. El hasta del 10% respecto a la evaluación
valor de STHD es evaluado en cada ventana y considerando el THD, según la expresión.
mostrado en la última fila de la Tabla. In 3STHD
Se observa como los resultados se aproximan =
bastante al valor teórico que se debería haber IL 1 + STHD 2
 His areas of interests include power system analysis,
harmonics, protection of power systems, Power system
5. REFERENCIAS reliability and power quality; he has also published papers
in the fields of power quality, protection systems, power
systems and systems grounding
[1] SAYAS POMA, L., "SEMINARIO Optimización de las Mr. Sayas is a registered professional engineer in the city
Instalaciones Eléctricas", PROCOBRE PERÚ, noviembre of Lima
1997.
[2] SAYAS POMA, L., "Armónicos en Sistemas Eléctricos
de Potencia", Curso para fiscalizadores I OSINERG, julio
2002.
[3] SAYAS POMA, L., "Calidad de la Energía Eléctrica",
Curso para Profesionales de Electro Puno S.A., abril 2001.
[4] SAYAS POMA, L., "Calidad de la Energía Eléctrica en
Instalaciones de Baja Tensión", Conferencia Técnica, XIV
CONIMERA, octubre del 2001.
[5] SAYAS POMA, L., "Soluciones Electrotécnicas de
armónicos en Sistemas Eléctricos de Potencia", Curso para
fiscalizadores de las Unidades de Transmisión y
Distribución I OSINERG, octubre 2003.
[6] HWEI P. HSU, "Analysis de Fourier “, Texto ADDISON-
WESLEY IBEROAMERICANA. Wayne State University,
Michigan, 1987
[7] Herraiz, G. T. Heydt, E. O´Neill-Carrillo, "Power Quality
for Aperiodic Voltages and Currents.” IEEE Trans. On
Power Delyvery, Vol. 15, No 2, April 2000, pp.784-790
[8]S. Herraiz, S. J. Tesis Doctoral “Aportaciones al estudio
del flujo armónicos de cargas” , Universidad Politécnica de
Cataluña, Barcelona, mayo de 2002.
[9]IEEE Recommended Practices and Requirementes for
Harmonic Control in Electrical Power Systems. IEEE Std
519-1992.
[10] Norma Técnica de Calidad de los Servicios Eléctricos-
1997, NTCSE-Peruana., ministerio de Energía y Minas
Perú.

BIOGRAPHY:

L. Sayas P., was born in Huancayo, Peru, on April 22,

1962. He receiving his degree with honors in Electrical
Engineering from the University of Center in 1984 and he
has a Masters Degree in Power Systems in 1995, from the
University Engineering of Peru. He has over 20 years of
generation, transmission and distribution systems
experience. He was an Professor of Power Systems at The
Superior Institute of Technology (1993-2003), also he was
head Assistant of the Electrical Department from 2001-
2003.
He was an advisory of PROCOBRE, TECSUR and He was
a General Manager of PT Power Systems Technologies.
Now, he is Supervisor of electrical power system
interruption at OSINERG.