1. Introducción
En la actualidad, la energía eléctrica se ve como un producto y no solo como un servicio; como todo producto, se cuantifica por sus características como la calidad, el precio y la accesibilidad. La calidad del producto eléctrico se mide a partir de variables cuantificables, siendo una de ellas la distorsión armónica (THD), esta variable mide porcentualmente cuánto se deforma la forma de onda sinusoidal de voltaje o corriente respecto de su onda fundamental 1. Para calcular el THD es necesario encontrar las frecuencias armónicas inmersas en la onda eléctrica, que son frecuencias múltiplos de la frecuencia fundamental de trabajo del sistema; cuya amplitud va decreciendo conforme aumenta el múltiplo y la causa fundamental de la distorsión de las ondas de voltaje y corriente.
La interconexión de los sistemas eléctricos es cada vez más compleja. Una gran cantidad de cargas no lineales provocan la contaminación armónica grave y compleja en la red eléctrica; por ende, el control de los armónicos es de gran importancia. La señal en la aplicación práctica tiene la característica compleja de la interferencia del ruido, que incide directamente en el proceso de detección de armónicos 2. Por lo tanto, eliminar el ruido y hacer que la detección de armónicos sea estable lo más rápido posible es de gran importancia para evaluar la calidad de energía eléctrica que se está suministrando 3.
Hoy en día, las técnicas computacionales en paralelo permiten obtener una mejora en el tiempo de cálculo informático y el procesamiento de amplios conjuntos de datos 4. Como existe una gran cantidad de núcleos de computación, el tiempo de procesamiento se puede disminuir cuando una tarea compleja inicial se puede descomponer en subtareas independientes más simples asignando un núcleo de procesador a cada una, ejecutándolas simultáneamente 5-7.
La medición y el análisis de armónicos son partes clave de la evaluación de la calidad de la energía 8,9. Las razones para que se produzcan armónicos son varias, entre las más frecuentes están las originadas por los equipos electrónicos que consumen energía eléctrica de una forma no lineal, es decir, de una forma no continua en el tiempo; también está la variación de los parámetros de un circuito e impedancias, y la conexión y desconexión de equipos 6-10. Los armónicos en la demanda introducen fluctuaciones, caídas de tensión, desequilibrios en la red conduciendo a posibles daños y fallos de energía.
Para el análisis de la calidad de energía eléctrica ha sido pertinente el uso de técnicas de procesamiento digital de señales, técnicas que permiten la obtención de indicadores para evaluar la calidad del servicio, la literatura registra aplicaciones para el estudio del nivel de armónicos, medidas de precisión en sistemas eléctricos, estudio de potencia reactiva y análisis de transitorios entre otros indicadores.
El método clásico para la detección de distorsión armónica (THD) de una onda de voltaje o corriente presente en la red eléctrica es la transformada de Fourier (FT), de este método se derivan numerosas técnicas alternativas de procesamiento digital de señales para la localización de armónicos, a continuación se detallan las más relevantes 11:
Transformada de Fourier de corto tiempo (STFT).
Transformada de Wavelet (WT).
Transformada de Gabor (GT).
Transformada S (ST).
Filtros Digitales.
Filtros Kalman (KF).
Clase de Cohen.
Métodos Paramétricos.
Sensado Comprimido.
El presente trabajo analiza las ventajas y desventajas de cada una de estas técnicas alternativas, las cuales son de gran utilidad para la detección de armónicos. Centra su temática en la detección de armónicos mediante la técnica de sensado comprimido 12. En la Figura 1 se plasma la idea principal del presente artículo.
1.1. Trabajos relacionados
El crecimiento acelerado de la tecnología, por ende de los equipos electrónicos, hace que se aumente considerablemente las cargas no lineales, siendo estas fuentes de perturbaciones eléctricas que ocasionan anomalías en la tensión o corriente presentando variaciones de magnitud, o valores cuadráticos medios, con respecto a su valor nominal durante un intervalo de tiempo; estos problemas afectan directamente a los consumidores en tiempos y costos de producción 13. Existen varias técnicas y algoritmos para la detección de distorsión armónica en las señales eléctricas, mismas que tienen ventajas y desventajas. A continuación se citan las técnicas más usadas 12,13.
1.1.1. Transformada de Fourier (FT). - Es la técnica más conocida para el análisis en el dominio de la frecuencia del cual derivan las demás técnicas. Sin embargo, funciona bien para el caso de tiempo infinito de una señal estacionaria y no puede resolver ninguna información temporal asociada con las fluctuaciones, para resolver esto se usa la transformada de Fourier en tiempo corto (STFT) 12,14.
1.1.2. Transformada de Fourier de corto tiempo (STFT). - Divide la señal en pequeños segmentos, donde se puede suponer que estos segmentos de señal están estacionarios 1. El STFT funciona bien siempre que la ventana sea lo suficientemente corta en comparación con la tasa de fluctuación; las altas tasas de fluctuación pueden dar lugar a errores importantes 15.
1.1.3. Transformada de Wavelet (WT). - Es otra técnica para obtener características en el dominio tiempo-frecuencia. Se han propuesto muchos tipos de ondículas madre, incluida la ondícula compleja. Otros trabajos utilizan la combinación de WT y FFT para obtener ciertas características de la señal bajo análisis. La señal que se analiza se descompone primero en distintas representaciones: una rica en altas frecuencias y la otra en bajas frecuencias, procesando la señal se filtra en niveles sucesivos de detalle. Este proceso está fuertemente influenciado por el ruido superpuesto en la señal y la naturaleza iterativa de los algoritmos basados en el WT, también sufre de fugas espectrales y efectos de valla con impacto en la estimación 15-17.
1.1.4. Transformada de Gabor (GT). - Posee una buena resolución tiempo-frecuencia y buena relación señal-ruido. La constancia en el ancho de todos sus elementos limita su uso a altas frecuencias GT se realiza en cada muestra utilizando todos los valores de muestra en sus ventanas de tiempo, por lo tanto, la complejidad computacional de estos algoritmos está directamente asociada con la frecuencia de muestreo.
1.1.5. Transformada S (ST). - Esta transformada tiene una corrección de fase asegurando que las amplitudes máximas sean regiones de fase estacionarias, se puede referir la fase del espectro al origen del eje del tiempo o a un punto. Debido a que los anchos de las ventanas de frecuencia, las transformadas S son proporcionales a su frecuencia central, esto puede ocasionar una medición incorrecta de los armónicos 18. el ST es una extensión del GT y del WT, se basa en la forma de procesamiento de bloques y por lo tanto no satisface bien el requisito de tiempo real 12.
1.1.6. Filtros Digitales. - Brindan una menor complejidad computacional que wavelet y STFT. A pesar de que las frecuencias centrales de los filtros de paso de banda se establecen en las frecuencias armónicas, la salida no siempre corresponde a los componentes armónicos puros. Las bandas de frecuencia de los filtros de paso de banda pueden superponerse parcialmente entre sí y las propiedades de frecuencia de magnitud de estos filtros no son planas en sus respectivas bandas de frecuencia 12.
1.1.7. Filtros Kalman (KF). - Poseen buena estimación de la relación de cambio frecuencia - amplitud y contenido armónico, buena relación señal-ruido. No son muy resistentes al ruido aditivo, su complejidad de cálculo es alta la observación del estado no es lineal, podría hacer que el KF divergiera en algunas condiciones iniciales mal elegidas. La inestabilidad KF puede conducir a una estimación inexacta de los componentes armónicos y fundamentales de una sinusoide en ruido 19.
1.1.8. Clase de Cohen. - Reduce los problemas de interferencia dependiendo del núcleo, utilizado alta resolución en el dominio del tiempo y frecuencia. Las características de la función de distribución dependen de la elección del núcleo 17.
1.1.9. Métodos Paramétricos. Tiene una resolución excepcional de tiempo-frecuencia. Los parámetros de la señal de prueba y entrenamiento están permitidos en una amplia gama y de manera controlada. Gran consumo de recursos computacionales y tiempo de procesamiento requiere una buena suposición de la distribución estadística de los datos del patrón 12,14,19.
1.1.10. Sensado Comprimido (SC).- Es una técnica que ayuda a la recuperación de señales de manera eficiente respecto del costo computacional y tiempo de cálculo, permite representar señales mediante un número reducido de medidas, y a su vez, ofrece herramientas capaces de reconstruir las señales originales con muy poca pérdida de información 8,9. Debido a las características antes mencionadas esta técnica es la que admite la mayoría de las ventajas de las técnicas antes mencionadas, a más de su reconocimiento de armónicos, se puede reconstruir señales y crear filtros para la reducción de la distorsión armónica (THD) en señales eléctricas.
1.2. Formulación del problema
Las técnicas para localizar distorsiones armónicas, descritas en trabajos relacionados, parten del principio fundamental de la Transformada de Fourier (FT), la cual lleva la señal del dominio del tiempo al dominio de la frecuencia. Sin embargo, funciona bien para el caso de tiempo infinito de una señal estacionaria y no puede resolver ninguna información temporal asociada con las fluctuaciones.
Todas las técnicas concuerdan en que: en el dominio de la frecuencia es donde realizan su análisis para la detección de armónicos, encontrándose con los mayores desafíos: el ruido inmerso en la señal, grandes recursos computacionales, elevado tiempo de procesamiento de la señal, definición de ancho de ventana y la gran cantidad de datos que deben ser analizados continuamente. Lo que limita su desempeño en línea e implementación en plataformas tecnológicas como: procesadores de señales digitales, arreglos de puertas programables de campo, microprocesadores o microcontroladores.
2. Metodología
La metodología del sensado comprimido capta las características más relevantes de las técnicas analizadas y aporta con todas sus ventajas haciendo de esta técnica una potente herramienta para la detección de distorsión armónica, a más de su reconocimiento de armónicos, el sensado comprimido tiene la capacidad de estimar estados, localizar fallas eléctricas, en si muchas aplicaciones de estimación y reconstrucción de señales.
La detección de distorsión armónica en redes eléctricas de distribución por medio del sensado comprimido inicia adquiriendo una señal en el dominio del tiempo de tipo discreta. Esta señal es almacenada en un vector [¦]; de este vector se extrae una muestra aleatoria y reducida de la señal en estudio, dicha muestra es transformada al dominio de la frecuencia mediante una matriz diccionario [Y] para el caso de estudio se está usando la transformada discreta del coseno (DCT) como matriz diccionario.
La señal reducida en el dominio de la frecuencia es denominada vector de coeficientes totales [b], este vector es sometido a un proceso de optimización no lineal por medio de la norma ℓ1, dicho proceso clasifica a los coeficientes, extrae los más representativos y convierte en cero a los menos representativos, mediante el método de programación línea simplex revisado el cual la mayoría de los números que aparecen en el vector de coeficientes del método normal no se usan realmente en las iteraciones, por lo cual en el método revisado solo se calculan los valores necesarios para encontrar la solución óptima a través de matrices. El resultado del proceso de optimización no lineal es un vector de coeficientes reducidos o también llamado vector de mediciones [x], una combinación lineal del vector de mediciones [x] con una matriz de mediciones da como resultado una señal comprimida de similares características a la original con sus coeficientes reducidos, en la Figura 2 se ilustra el flujograma del proceso de detección de distorsión armónica mediante la técnica de sensado comprimido.
A continuación, se detalla la metodología matemática del sensado comprimido, con la ayuda del siguiente ejemplo, se cuenta con una señal teórica de voltaje de frecuencia fundamental de 60 (Hz) y una amplitud pico-pico de 200 (V) el tiempo de adquisición de datos de la señal es de 0.4 segundos, la señal es de tipo discreta con 4000 muestras esto implica que el dt es de 0.1 ms. La señal es una sinusoide pura desde t=0 hasta t=0.1 segundos.
De t=0.1 segundos a t= 0.2 segundos la señal fundamental se contamina con una señal de tercer armónico de amplitud 80 (V) pico-pico y frecuencia de 180 (Hz).
De t=0.2 segundos a t= 0.3 segundos la señal fundamental se contamina con una señal de quinto armónico de amplitud 100 (V) pico-pico y frecuencia de 300 (Hz).
De t=0.3 segundos a t= 0.4 segundos la señal fundamental se contamina con dos señales una de tercer armónico y una de quinto armónico, de amplitud 40 (V) pico-pico y frecuencia de 180 (Hz) y de amplitud 40 (V) pico-pico y frecuencia de 300 (Hz) respectivamente. La señal a ser estudiada con las técnicas de sensado comprimido se ilustran en la Figura 3.
Una señal rala o escasa no necesariamente debe cumplir con el teorema de Shannon/Nyquist que manifiesta que, para no perder la información en el proceso de reconstrucción de la señal a partir de la toma de muestras, la señal está limitada en banda y tasa de muestreo superior al doble de su ancho de banda. Según la técnica de sensado comprimido con unos pocos coeficientes se puede reconstruir la señal 6,20-22.
Esta metodología se basa en el tratamiento digital de señales, se inicia con la adquisición de una señal discreta la cual es almacenada en un vector [ ¦ ] para el presente ejemplo el vector [ ¦ ] es un vector fila de longitud m = (1x4000), esta señal es trasformada del dominio del tiempo al dominio de la frecuencia para extraer los coeficientes de la misma, por medio de una matriz denominada diccionario [ Y ]. La matriz diccionario [ Y ] o matriz subyacente puede ser modelada a partir de la trasformada rápida de Fourier (FFT), transformada Wavelet (DWT), transformada discreta del coseno (DCT), en sí cualquier trasformación lineal que lleve la señal sensada al dominio de la frecuencia, donde la señal se transforma en un vector de coeficientes [C], para el presente caso de estudio se usa como matriz diccionario [ Y ] la transformada del coseno (DCT), porque sus coeficientes son de tipo real 23-25.
La transformada como base en sensado comprimido, permite obtener representaciones dispersas de las señales, permitiendo usarlas como base en un sistema de sensado comprimido. La transformada discreta del coseno (DCT) es una transformada basada en la transformada de Fourier discreta (DTF), pero solo utilizan números reales 26. La DCT de una secuencia de datos 𝑋(𝑚),𝑚=0,1…,(𝑀―1) se define como la Ec.1 y 2:
Donde 𝐺𝑥(𝑘) es el k-ésimo coeficiente. La transformada inversa discreta coseno (IDCT) se define como la Ec.3:
Si la Ec.3 se escribe en su forma matricial A de tamaño M x donde: M describe la transformación coseno, entonces la propiedad de ortogonalidad puede ser expresada como la Ec.4:
Donde A T es la matriz transpuesta de A e I es la matriz identidad 26.
La matriz diccionario es una colección finita de coeficientes elementales de norma unitaria, la matriz diccionario [ Y ] permite capturar información de una señal dispersa para que esta pueda ser utilizada por los algoritmos para la recuperación de dicha señal 1. La Ec.5 muestra cómo obtener el vector de la señal discreta en dominio del tiempo en relación con [ Y ] y [ C ].
Donde:
¦ : Vector de la señal discreta en dominio del tiempo.
Y: Matriz diccionario (DCT)
C: Vector de coeficientes de la señal ¦ en el dominio de la frecuencia.
Una vez obtenido el vector de coeficientes [C] se analiza que efectivamente es un vector con la mayoría de los coeficientes tendiendo a cero por lo que se puede expresar que el vector [C] es realmente escaso. Ahora se debe obtener un vector aleatorio [t´] que es un subconjunto del vector [t], el vector aleatorio [t´] para el presente caso de estudio es 1/10 del vector [t], lo que significa n = 400 datos no repetitivos tomados aleatoriamente del vector [t], el vector [t´] es un vector de índices utilizado para obtener la muestra del vector columna [b] de la Ec.6 23-25,27.
Donde:
¦ : Vector de la señal discreta sensada en el dominio del tiempo.
F: Matriz de mediciones (subconjunto de filas del operador identidad).
b: Muestras aleatorias de la señal ¦.
La matriz [F] conocida como matriz de mediciones para el caso en estudio es una matriz de 400 filas por 4000 columnas, esto implica que tiene 400 vectores fila de longitud igual a 4000, cada vector fila de la matriz [ F ] tiene 3999 ceros y un 1 que corresponde donde se tomó la muestra de medición esto se expresa en la Figura 4.
En la construcción de la matriz de mediciones [F] es donde se puede analizar la potencialidad del sensado comprimido, si hablamos de un sistema eléctrico por ejemplo, la matriz [F] puede ser armada en función de la matriz de admitancias y se podría estimar estados en cada una de las barras, encontrar distorsiones armónicas en las barras, estimar las cargas de menor factor de potencia del sistema eléctrico, estimar pérdidas no técnicas del sistema, estimar despachos de carga, estimar curvas de demanda eso solo tomando como referencia a un sistema eléctrico, el sensado comprimido puede escalar a otras ramas de la ciencia y poder estimar parámetros en los que se tenga pocos datos con una buena aproximación y bajos recursos computacionales 13,28-35.
El vector [b] es un vector con (1/10) de muestras aleatorias de la señal ¦, por lo que la matriz de mediciones [Φ] es un subconjunto de las filas del operador identidad. En la Figura 5a se aprecia la señal de [¦ ] en el dominio del tiempo en forma de puntos, en la Figura 5b se observa los coeficientes de la señal en el dominio de la frecuencia.
Para obtener la matriz A se extrae la transformada discreta del coseno a la matriz de mediciones [F] de acuerdo a la Ec.7:
Para reconstruir la señal se minimiza el número de coeficientes del vector [C], obteniendo como respuesta un vector de coeficiente reducidos [x] que contiene los coeficientes optimizados de la señal [¦] para lograr este cometido se debe resolver la Ec.8 28,36:
Donde:
A: Matriz de coeficientes (m x n).
x: Vector de coeficientes reducidos.
b: Muestras aleatorias de la señal ¦
Como se trata de una compresión la matriz A es de tipo rectangular, con más columnas que filas (4000 columnas y 400 filas), para el cálculo del vector de coeficientes reducidos [x] implica resolver un sistema indeterminado de ecuaciones lineales con 10 veces más incógnitas que ecuaciones que es lo que se muestra en la Ec.8, para solucionar dichas indeterminaciones se usa las normas matriciales de aproximación, que se encuentran dadas en forma de vectores matriciales que ayudan a recuperar señales las cuales pueden llegar a tener seudonormas que al aproximar a los vectores se puede minimizar el error de aproximación o reconstrucción los tipos de normas pueden ser 𝑙1,𝑙2,𝑙∞.
Reconstrucción mínima de norma 𝑙2 Define la norma 𝑙𝑝 del vector de [x] con longitud N como la Ec.9:
Donde para obtener una recuperación aproximada de la señal analizada, es necesario emplear la Ec.1, donde la solución está contenida dentro del vector y empleando el concepto de aproximación para 𝑙2 se tiene la Ec.10:
Donde:
𝛹: es la matriz de reconstrucción de la señal.
𝑦: margen de error relativo.
La desventaja es que una minimización 𝑙2, generalmente no encuentra una solución adecuada, retornando un vector x disperso con varios elementos diferentes de cero. Reconstrucción mínima de norma 𝑙0. Debido a que la norma 𝑙2 hace una medición de la energía de la señal mas no su dispersión, es necesario considerar la norma que cuantifica el número de entradas diferentes de cero en el vector x. En consecuencia, un vector k-disperso contiene una norma 𝑙0 igual a k, dando como resultado un proceso de optimización modificado, como se observa en la Ec.11:
La reconstrucción en base a la norma 𝑙0 es capaz de recuperar la señal analizada con una alta probabilidad empleando únicamente M=k+1 mediciones gaussianas. La desventaja de este método radica que la solución de la Ec.11 es numéricamente inestable y del tipo NP-complete. Los problemas NP-complete son problemas muy complejos y no pueden ser resueltos por una máquina de Turing determinista. En consecuencia, son resueltos mediante máquinas no deterministas y un tiempo polinómico acotado. Los algoritmos basados en optimización convexa son otro enfoque de aproximación dispersa, los cuales remplazan las funciones combinatorias de la norma 𝑙0 con la función convexa 𝑙1 lo cual, convierte el problema combinatorio en un problema de optimización convexa. La norma 𝑙1 con la Ec.11 puede recuperar con una gran probabilidad k señales dispersas y comprimibles usando únicamente M ≥ cK log (N / K) iid mediciones gaussianas, como se indica en la Ec.12:
En la Figura 6 se puede observar la reconstrucción de las señales y los coeficientes del vector [x]. Este es un problema de optimización convexa que de forma práctica se reduce a un programa lineal conocido como seguimiento de bases, una vez encontrado el vector de mediciones [x] se puede reconstruir la señal mediante la Ec.13:
3. Resultados
De los problemas analizados en las diferentes técnicas de localización de distorsión armónicas en redes eléctricas, los más recurrentes son: los recursos computacionales, tiempo de procesamiento, ruido superpuesto en la señal y resolución de tiempo frecuencia de la señal. Todos estos problemas son superables con la técnica de sensado comprimido, es por esto que se propone esta técnica para la localización de distorsiones armónicas 1,37. En la Figura 7 se puede apreciar la comparación entre la señal original y la señal comprimida que se reconstruyó.
Como resultado de la señal teórica analizada se obtuvo la Figura 8, donde se puede analizar las frecuencias inmersas en la señal se observa 3 coeficientes más prominentes, en 60 Hz señal fundamental, en 180 Hz tercer armónico y en 300 Hz quinto armónico. Queda abierta la investigación para determinar el tiempo en el cual se producen los armónicos o alguna perturbación de tipo transitoria.
Una posible investigación o aplicación del algoritmo se sensado comprimido, es la detección de distorsiones armónicas en sistemas eléctricos, poder localizar las fuentes de THD sería de gran utilidad para las empresas distribuidoras de energía, para poder mitigar y mejorar su calidad de producto eléctrico, la detección se lograría con el uso de medidores sincrofasoriales (PMU) óptimamente ubicados en el sistema eléctrico, esta tecnología es costosa en la actualidad, pero a futuro sería una buena alternativa para la estimación de parámetros de calidad de energía, lo expuesto se ilustra en la Figura 9.
4. Conclusiones
Todas las técnicas para localizar distorsiones armónicas parten de la transformada de Fourier la cual lleva la señal del dominio del tiempo al dominio de la frecuencia, de esta transformada se derivan las demás transformadas, que son casos especiales de la transformada de Fourier. Los mayores desafíos a la hora de la detección de distorsión armónica en la red eléctrica son: el ruido inmerso en la señal en estudio, recursos computacionales, tiempo de procesamiento y ancho de ventana.
La metodología del sensado comprimido acoge las características más relevantes de las técnicas convencionales para la detección de distorsión armónica, con la diferencia que no cumple con el teorema de Nyquist, a más de su reconocimiento de distorsión armónica, se puede reconstruir señales y crear filtros para la reducción de la distorsión armónica (THD) en señales eléctricas.
Para el sensado comprimido siempre se requiere usar una matriz diccionario que no es más que usar alguna transformación lineal que lleve la señal del dominio del tiempo al dominio de la frecuencia, en el ejemplo del presente artículo se usó como matriz diccionario a la transformada del coseno (DCT), queda abierta la investigación probando con otras matrices diccionario que podrían encontrar una representación mejor de la señal utilizando un menor número de coeficientes.
La eficiencia del Sensado Comprimido (SC) depende del diccionario a utilizar mas no depende del ancho de banda de la señal bajo análisis, sino de la cantidad de información de la misma, y de la norma 𝑙1 de minimización no lineal que es de gran utilidad ya que recupera y aproxima las señales dispersas con errores muy bajos.
La técnica de SC es ideal para localizar frecuencias armónicas en estado estable en una señal eléctrica sea esta de voltaje o corriente, entendiendo por estado estable a una señal que posee un numero finito de muestras. Para detección de perturbaciones transitorias se requiere de un tiempo considerable de procesamiento computacional, debido a la cantidad de muestras que debe tomar para no perder detalles a la hora de reconstruir la señal.
En la construcción de la matriz de mediciones [F] es donde se puede analizar la potencialidad del sensado comprimido, esta matriz es armada en función de lo que se requiera analizar o estimar por ejemplo un sistema eléctrico de potencia, la matriz [F] puede ser armada en función de la matriz de admitancias y se podría estimar varios parámetros en un sistema eléctrico de potencia.