Dfcarriong (talk | contribs) |
Dfcarriong (talk | contribs) |
||
Line 12: | Line 12: | ||
==1 Introducción== | ==1 Introducción== | ||
− | El estimador de estado tiene como objetivo principal detectar, identificar y eliminar las mediciones y/o parámetros erróneos de los Sistemas Eléctricos de Potencia (SEP) y proporcionar el estado más probable del mismo | + | El estimador de estado tiene como objetivo principal detectar, identificar y eliminar las mediciones y/o parámetros erróneos de los Sistemas Eléctricos de Potencia (SEP) y proporcionar el estado más probable del mismo [1], mediante un proceso que toma como entrada las mediciones del SEP y obtiene como salida las variables de estado (voltaje) del SEP. Hoy en día, la estimación del estado dinámico es una tarea crucial para el monitoreo y control de los sistemas eléctricos, particularmente después de una desconexión o pérdida de una línea de trasmisión, para realizar una operación restaurativa adecuada del sistema eléctrico [Kumar2018]. En la actualidad se pueden encontrar varias investigaciones enfocadas en aspectos y metodologías de diferentes estimadores de estado, que van desde su formulación matemática, hasta los problemas para su implementación [Shabaninia2012]. |
Los Estimadores de Estado (EE) requieren de toda la topología de red, parámetros de las líneas de transmisión, así como las mediciones provenientes de los sistemas de medición convencional SCADA, y de la Unidad de Medición Fasorial (PMU) [Habiballah2018] y se proponen la condición de estimar todos los fasores de voltaje de cada nodo del SEP, estas funciones se realizan en los centros de control modernos [Singh2018]. | Los Estimadores de Estado (EE) requieren de toda la topología de red, parámetros de las líneas de transmisión, así como las mediciones provenientes de los sistemas de medición convencional SCADA, y de la Unidad de Medición Fasorial (PMU) [Habiballah2018] y se proponen la condición de estimar todos los fasores de voltaje de cada nodo del SEP, estas funciones se realizan en los centros de control modernos [Singh2018]. | ||
Line 651: | Line 651: | ||
===BIBLIOGRAPHY=== | ===BIBLIOGRAPHY=== | ||
− | + | [1] S. Venugopal and R. Sunitha, “Robust Dynamic State EstimationDuring Missing of PMU Data,”2018 4th International Conferencefor Convergence in Technology (I2CT) SDMIT Ujire, Mangalore,India. Oct 27-28, 2018, pp. 1–6, 2020. [Online]. Available:doi:10.1109/i2ct42659.2018.9058140. |
El presente trabajo propone diseñar y aplicar un algoritmo de Estimador de Estado en los Sistemas Eléctricos, aplicando la metodología de los Mínimos Cuadrados Ponderados, mediante un proceso denominado estimador hibrido no lineal el mismo que utiliza las mediciones provenientes de las PMU y las mediciones convencionales, además, de encontrar las ubicaciones óptimas de PMU respetando las restricciones de observabilidad y redundancia para los Sistemas Eléctricos de Potencia, a través de la programación lineal entera mixta, considerando las contingencias N-1 en el SEP. El estimador de estado propuesto se ajusta a las mediciones obtenidas, ante la presencia de contingencias en el SEP. Las simulaciones que se obtienen como resultado de la implementación del estimador de estado hibrido, muestran que la propuesta mejoran la precisión del estimador y la velocidad de convergencia, considerando el mínimo error de mediciones del sistema y número de PMU colocadas en el Sistema de Potencia, luego de presentarse una contingencia, que pueden ser la salida de una PMU o la pérdida de una de las línea de transmisión del SEP.
Palabras Clave: Estimador de estado, PMU, Mínimos Cuadrados, Ubicación Óptima de PMU, Sistemas eléctricos de Potencia.
The present research proposes to design and apply a state estimator algorithm in electrical power systems, applying the Weighted Least Squares methodology, through a process called non-linear hybrid estimator, the same one that uses measurements from PMUs and conventional measurements. in addition, to find the optimal PMU locations respecting the observability and redundancy restrictions for the Electric Power Systems, through mixed integer linear programming, considering the N-1 contingencies in the SEP. The proposed state estimator is adjusted to the measurements obtained, in the presence of contingencies in the SEP. The simulations obtained as a result of the implementation of the hybrid state estimator show that the proposal improves the precision of the estimator and the speed of convergence, considering the minimum error of system measurements and the number of PMUs placed in the Power System, after a contingency occurs, which may be the departure of a PMU or the loss of one of the transmission lines of the SEP.
Keywords: State estimator, PMU, Least squares, Optimal PMU location, Electrical power system.
El estimador de estado tiene como objetivo principal detectar, identificar y eliminar las mediciones y/o parámetros erróneos de los Sistemas Eléctricos de Potencia (SEP) y proporcionar el estado más probable del mismo [1], mediante un proceso que toma como entrada las mediciones del SEP y obtiene como salida las variables de estado (voltaje) del SEP. Hoy en día, la estimación del estado dinámico es una tarea crucial para el monitoreo y control de los sistemas eléctricos, particularmente después de una desconexión o pérdida de una línea de trasmisión, para realizar una operación restaurativa adecuada del sistema eléctrico [Kumar2018]. En la actualidad se pueden encontrar varias investigaciones enfocadas en aspectos y metodologías de diferentes estimadores de estado, que van desde su formulación matemática, hasta los problemas para su implementación [Shabaninia2012].
Los Estimadores de Estado (EE) requieren de toda la topología de red, parámetros de las líneas de transmisión, así como las mediciones provenientes de los sistemas de medición convencional SCADA, y de la Unidad de Medición Fasorial (PMU) [Habiballah2018] y se proponen la condición de estimar todos los fasores de voltaje de cada nodo del SEP, estas funciones se realizan en los centros de control modernos [Singh2018].
Los EE tradicionales utilizan las mediciones provenientes de los flujos de potencia e inyecciones de potencia real y reactiva, tomando en cuenta las mediciones de los fasores de las corrientes proveniente de las PMUs en la red; se puede plantear un esquema de Estimador de Estado Híbrido (EEH) con mediciones provenientes de las PMU y SCADA, debido a que no hay la cantidad suficiente de PMU instaladas en el sistema eléctrico [Noopura2017].
Una de las formulaciones efectivas para solución de los EE tradicionales y EEH [Zhang2016], es la metodología de los Mínimos Cuadrados Ponderados (WLS), el EEH considera las mediciones convencionales y fasoriales, los resultados demuestran que las simulaciones del algoritmo de EEH mejora en gran medida la precisión y obtienen un mejor comportamiento en comparación con el estimador de estado tradicional [Lin2018].
La determinación del comportamiento del SEP en tiempo real, basados en la medición de los diferentes parámetros eléctricos es uno de los retos actuales que tiene los operadores del sistema de potencia [Ganjkhani2019]; para lo cual se han desarrollado diferentes técnicas que usan la información obtenida por las PMUs [Nanda2017]. La presencia de contingencias genera dificultades en la velocidad de respuesta de los algoritmos estimadores de estado y errores por el cambio en la topología del sistema de eléctrico [Raposo2017].
El trabajo de la PMU permite realizar, con gran precisión y velocidad, el análisis dinámico del estado del SEP, utilizando la técnica de estimación de fasores, este mecanismo es esencial para análisis e interpretación de los datos recopilados de manera correcta [IsraelAkingeneye]. Por lo tanto, un proceso de estimación de fasores tiene mucha importancia en diseño de mecanismos de control y diversos sistemas de protección [Shrivastava2018].
La tecnología de medición de sincrofasores tiende a convertirse en una tecnología madura con el despliegue de una gran cantidad de PMU en muchos sistemas eléctricos importantes en todo el mundo, en las últimas dos décadas [Bundoc2019] con las adecuadas redes de comunicaciones las mediciones recopiladas por las PMUs, en diferentes ubicaciones del sistema, se sincronizan y se ajustan en el tiempo, para después combinarse y proporcionar una fotografía completa y precisa del SEP, en diferentes puntos y a grandes distancias logrando que se interconecten en tiempo real. Lo cual permite reaccionar, determinar y tomar acciones a los operadores del SEP, para mantener la continuidad y calidad del servicio [Almasabi2018], existen algunos obstáculos para utilizar los datos de la PMU en aplicaciones en tiempo real; de los cuales podemos mencionar: la necesidad de un software compatible con el Concentrador de datos Phasor Data Concentrator (PDC), para organizar los datos de las PMUs [Bi2015], además, la comunicación de ancho de banda para transmitir una gran cantidad de datos generados por las PMU ubicadas en diferentes sitios del sistema eléctrico. Las PMUs no se limitan solo a recolectar datos, además, tienen las funciones para la desconexión automática de elementos del SEP, considerando ciertas condiciones. Así mismo, los PDC permiten al Sistema de Monitores de Control y Protección (WAMPC), acoplarse a diferentes sistemas de control para mejorar la detección y resolución de fallas [Patil2018].
La PMU mide directamente los fasores de voltaje y corriente del nodo donde está instalado y los fasores actuales de las líneas que inciden en ese nodo. Por lo tanto, no es exigente instalar PMU en todos los nodos del SEP, para logar que el sistema sea completamente observable. Existen varios algoritmos para la ubicación óptima de PMU [Aminifar2013,Aminifar2014], para su selección de ubicación e instalación estratégicas y para que se realicen de manera inteligente, con un número máximo de mediciones disponibles que aumentan la observabilidad del sistema [Ghosh2018].
Desde los últimos años, algunas herramientas matemáticas directas como la programación lineal entera mixta (MILP), el algoritmo codicioso (voraz), así como varias técnicas no deterministas como algoritmo genético (GA), optimización de enjambre de partículas (PSO), etc.; han sido utilizados ampliamente, para obtener una solución al problema de ubicación óptima de PMU [Tawfik2018].
La observabilidad de una o varias PMUs se refiere a la región en la cual se puede estimar confiablemente la capacidad de observación del sistema eléctrico [Hyacinth2021,Shahriar2018], la cual, puede reducirse si algunas PMUs dejan de estar disponibles, por lo cual se requiere maximizar la utilización de diferentes PMUs para mantener la observabilidad del sistema y la precisión bajo cualquier escenario de medición incompleto [Zhang2017,Tarif2019]. Debido a que los costos asociados (gestión, compilación de gran cantidad de datos de los sincrofasores) para la adquisición e instalación de PMU es elevado, la ubicación óptima de las PMUs es un papel crítico, la programación lineal entera mixta, es la más utilizada, para encontrar ubicaciones óptimas de PMU, para la completa observabilidad del sistema ante las contingencias que pueden presentarse en los sistemas eléctricos [Ahmed2019, CarrionGalarza2019].
La presente investigación propone un modelo de Estimador de Estado, que utilice las mediciones tanto de los PMU y SCADA, con el cual, se pueda determinar el comportamiento del Sistema de Potencia con el menor error luego de una contingencia N-1 que pueden ser la salida de PMUs o pérdida de líneas de transmisión del SEP y mantener el número adecuado de PMU en el SEP. Con el objeto de verificar la precisión del EE ante la pérdida de las mediciones provenientes de las líneas de transmisión y las PMUs del SEP, se utilizó los modelos de prueba de la IEEE de 14, 30 y 118 barras.
El trabajo consta de V secciones, que se detallan a continuación: sección I introducción, sección II metodología para la ubicación óptima de PMU; sección III Estimador de Estado, tomando en cuenta las mediciones convencionales y las provenientes de los PMU; sección IV presenta las gráficas obtenidas de las simulaciones, ante las contingencias que pueden presentarse en el SEP y los beneficion de la implementación del estimador de estado híbrido y para finalizar en la sección V se detallan las conclusiones obtenidas del presente trabajo.
La observabilidad de los SEP se fundamentan en las Leyes de Ohm y Kirchhoff, mediante tres reglas, directa, indirecta y la calculada, con las cuales se puede obtener los fasores de corriente y voltaje en todos los nodos del sistema [CarrionGalarza2019].
Para el despliegue óptimo de PMUs en el SEP, se debe construir como paso inicial, la matriz cuadrada nxn de conectividad (CM), la misma que varía según la topología de la red, para cada SEP [Ahmed2019], [CarrionGalarza2019].
|
(1) |
Donde:
i=nodos de salida
j=nodos de llegada
Con el fin de minimizar el número de PMU a ser ubicadas se describe la función objetivo en (2); tomando en cuenta las restricciones de redundancia y observabilidad, descritas en las restricciones: (3) indica que al menos una PMU se encuentra instalada en el SEP, (4) mantiene al máximo la redundancia del SEP y (5) restringe al 100 % la observabilidad del sistema.
FO:
|
(2) |
Restricciones:
|
(3) |
|
(4) |
|
(5) |
Donde:
=números de PMU
=nodos de ubicación de PMU
n=nodos del sistema
= es una variable de observabilidad (1 si es observable, caso contrario 0)
X_i=\left 1 si requiere de PMU en el nodo i
0 no requiere PMU \right.
Ubicación óptima de PMU considerando las contingencias del SEP Las contingencias se definen como eventos muy frecuentes presentes en los sistemas eléctricos, por lo cual, los instrumentos de medición deben mantener los registros de los parámetros eléctricos, aún cuando las contingencias N-1 continúen presentes en el sistema. Para ubicar las PMUs de manera óptima, se debe considerar las contingencias como la pérdida de una PMU o la salida de una línea de transmisión. [CarrionGalarza2019].
Ante la salida de una PMU, se pretende que el sistema sea observable completamente, por lo cual, cada nodo debería ser observado por 2 PMU, la función objetivo se define en (6); restringida por: (7) indica que al menos una PMU se encuentra instalada en el SEP, (8) mantiene al máximo la redundancia del SEP y (9) restringe al 100 % la observabilidad del sistema.
FO: Restricciones:
|
(7) |
|
(8) |
|
(9) |
Cuando se presenta en el sistema la salida de una línea de transmisión, la misma debe estar asociada a una PMU o a su vez cada nodo debe ser observado por 2 PMU; con la seguridad de no perder la observabilidad del sistema, se aplica la función objetivo (10); restringida por: (11) que indica que al menos una PMU se encuentra instalada en el SEP, (12) mantiene al máximo la redundancia del sistema y (13) & (14) restringe al 100 % la observabilidad del sistema.
FO:
|
(10) |
Restricciones:
|
(11) |
|
(12) |
|
(13) |
|
(14) |
Esta sección describe la metodología propuesta para la elaboración del EEH que integra las mediciones fasoriales ( y ) al proceso de EE convencional (, , y ), obteniendo como resultado un EE que genera la óptima estimación de las variables de estado (voltaje en magnitud y ángulo) del sistema, con el menor número de iteraciones, manteniendo la redundancia de las mediciones, cuando se presentan contingencias.
El EE convencional se basa en la metodología de WLS cuya función objetivo es minimizar la suma ponderada de los cuadrados de las desviaciones, resulado de la diferencia de las mediciones estimadas y las mediciones reales z [Lin2018] y utiliza el grupo de mediciones, convencionales (, , y ), además, de las magnitudes de tensiones nodales, la ecuación no lineal que relaciona las medidas y el vector del estado se define:
|
(15) |
Donde:
z=vector de mediciones (, , y )
h(x)=vector de las funciones no lineales
e=vector de error de medición y número de mediciones
También podemos escribir (15) de la siguiente manera:
|
(16) |
Las cantidades con sombrero son los valores estimados de las cantidades sin sombrero correspondientes, los errores estimados se lo pueden obtener mediante:
|
(17) |
La matriz jacobiana del sistema se obtiene a partir de las derivadas parciales con respeto al vector de estado x, la cual se desea minimizar y se define en (18).
|
(18) |
El vector h (x) relaciona a las mediciones con las variables de estado, se encuentra integrado por funciones no lineales, las cuales, calculan el valor verdadero de las mediciones del sistema, permitiendo linealizar el problema de estimación alrededor de un valor dado para todos los componentes del vector de estado, tenemos la matriz de ganancia G definida para los valores de k= 0, 1,2... hasta una apropiada convergencia, definida de la siguiente manera:
|
(19) |
Como se desea minimizar la función, mediante la deducción de WLS, remplazando (19) en (18), resulta:
|
(20) |
Para resolver la ecuación de estimador de estado , se sigue el mismo procedimiento que el de flujos de potencia por Newton Raphson, los valores en módulo y ángulo de las variables de estado (fasor voltaje) ingresan como datos iniciales de entrada, para la segunda iteración y el proceso continua hasta obtener el menor error al cual converge el sistema. De forma alternativa podemos escribir (20), de la siguiente manera:
|
(21) |
La matriz de errores medidos se denomina, matriz de covarianza (R), los valores de cada elemento fuera de la diagonal son cero porque los errores se han supuesto independientes, los valores esperados de la diagonal corresponden a las varianzas de , obtenemos una matriz a partir de:
|
(22) |
Además, definimos los términos de las derivadas parciales dependientes de x, de la siguiente manera:
|
(23) |
Para realizar la estimación de estado utilizando las corrientes medidas por las PMUs se procede de manera similar para el caso del EE que utiliza las mediciones convencionales (, , y ), al vector de estado se lo expresa en parte real e imaginaria, las funciones se las obtienen a partir del modelo generalizado, estas funciones relacionan al vector de estado con las mediciones de corriente fasoriales [Noopura2017,Lin2018], que se expresan de la siguiente manera:
|
(24) |
Donde:
= parte real de corriente
= parte imaginaria de corriente
= enlace del elemento conectado entre i y j
= magnitud de voltajes de los nodos i y j
= ángulo de voltajes de los nodos i y j
La matriz jacobiana la podemos definir de la siguiente manera:
|
(25) |
El modelo híbrido de EE, que considera las mediciones de , , , y , , es definido mediante:
|
(26) |
La matriz de covarianzas híbrida está constituida de las matrices de covarianza del EE convencional e Híbrido, la matriz se define como la de covarianza de las mediciones provenientes de las mediciones real e imaginaria de las corrientes, la matriz de covarianzas híbrida se define de la siguiente manera:
|
(27) |
Una vez definida la matriz de covarianzas y la matriz jacobiana para el EEH, podemos rescribir la ecuación (20), con lo cual obtenemos:
|
(28) |
Donde: = mediciones provenientes del EE tradicional y la provenientes de los PMU parte real e imaginaria
= vector de las funciones no lineales, que calcula los fujos de potencia, inyecciones de potencia activa y reactiva; corrientes real e imaginaria del SEP.
Diagrama de Flujo
En el siguiente diagrama de flujo se muestra el procedimiento para la resolución del EEH, mediante el método WLS, como primer paso se ingresa los datos de las mediciones del SEP tanto las convencionales (, , , ) como las provenientes de las PMUs ( y ) y las características de las líneas de transmisión (,) y topología del SEP. Si se excede el número de iteraciones permitidas, el algoritmo despliega un mensaje mencionando que el sistema no converge, si el es menor al fijado en las condiciones iniciales y la solución encontrada se encuentra dentro del número adecuado de iteraciones, el algoritmo despliega los valores de los vectores de estado.
![]() |
.
La Tabla 1, muestra el contraste de número de PMU en el SEP, ante la presencia de contingencias y sin contingencias, cuando se consideran las contingencias aumenta el número de PMU a ser colocadas en el SEP, con la ventaja que ante la pérdida de algún elemento del SEP se conserva la observabilidad del todo el sistema.
Ubicación de PMU | ||
Modelo IEEE | Sin contingencia | Con contingencia |
14 nodos | 2, 6, 7, 9 | 1, 3, 6, 8, 9, 10, 13 |
30 nodos | 1, 2, 6, 10, 11, 12, 15, 19, 25, 29 | 2, 3, 7, 8, 10, 11, 12, 13, 15, 16, 19, 21, 24, 26, 27, 29 |
118 nodos | 1, 5, 9, 12, 15, 17, 21, 25, 28, 34, 37, 40, 45, 49, 52, 56, 62, 64, 68, 70, 71, 76, 77, 80, 85, 87, 91, 94, 101, 105, 110, 114 | 1, 5, 6, 8, 10, 11, 12, 15, 17, 19, 21, 23, 26, 27, 29, 32, 34, 36, 37, 38, 40, 42, 44, 46, 49, 50, 51, 53, 56, 59, 61, 64, 67, 68, 72, 73, 74, 75, 76, 78, 80, 83, 85, 87, 89, 91, 92, 94, 92, 100, 102, 105, 106, 109, 111, 112, 114, 116, 117 |
El EEH muestra un mejor desempeño, debido a que toma en cuenta mediciones de , , , , y , que logran que el número de iteraciones disminuya y aumente su precisión, en contraste con el EEC y el EE PMU. Las mediciones del EE PMU son incluidas en el proceso de EEC, para logar una matriz H como se muestra en (26), que ante la pérdida de las mediciones provenientes de las líneas de transmisión en el SEP o de las PMUs mantengan la precisión y calidad del vector de estado resultante, como se puede corroborar en los siguientes escenarios de contingencias analizados.
![]() |
Figure 1: Variación de voltaje en p.u. debido a la disminución de PMU, para el EEH |
![]() |
Figure 2: Variación de ángulo en grados, debido a la disminución de PMU, para el EEH |
![]() |
Figure 3: Variación de voltaje en p.u. debido a la disminución de PMU, para el EE-PMU |
![]() |
Figure 4: Variación de ángulo en grados, debido a la disminución de PMU, para el EE-PMU |
Para el EE-PMU, el las figuras 3 y 4, se puede observar que, ante la salida de la PMU 1, se presenta el peor escenario para el sistema eléctrico, debido a que la pérdida de la PMUs causa distorsiones y la pérdida de la precisión el estimador debido a que la salida de las PMUs disminuye la cantidad de mediciones del estimador de estado. Siendo la pérdida de la PMU 1, la más severa para el sistema, obteniendo el valor máximo de voltaje en p.u. de 1.269 y la máxima diferencia angular de 9.060 grados, entre los valores reales y los estimados.
![]() |
Figure 5: Variación de voltaje en p.u., debido a la pérdida de líneas de transmisión, para el EEH |
![]() |
Figure 6: Variación del ángulo, debido a la pérdida de líneas de transmisión, para el EEH |
![]() |
Figure 7: Variación de voltaje en p.u., debido a la pérdida de líneas de transmisión, para el EE PMU |
![]() |
Figure 8: Variación del ángulo, debido a la pérdida de líneas de transmisión, para el EE PMU |
![]() |
Figure 9: Variación de voltaje en p.u., debido a la pérdida de líneas de transmisión, pare el EEC |
![]() |
Figure 10: Variación del ángulo, debido a la pérdida de líneas de transmisión, para el EEC |
Para el EEH, en los gráficos 5 y 6 se puede observar que, ante la salida de las mediciones asociadas a las líneas de transmisión, el estimador de estado mantiene su precisión, ya que las mediciones convencionales mantienen la precisión del estimador, el peor escenario para el sistema eléctrico, se presenta ante la pérdida de la línea de Transmision 4-7, obteniendo el valor máximo de voltaje en p.u. de 1.121 y la máxima diferencia angular de 3.947 grados, entre los valores reales y los estimados.
Para el EE PMU, en los gráficos 7 y 8 se puede observar que, ante la salida de las PMUs, el estimador de estado presenta valores excesivos en los valores estimados, ya que las pérdidas de mediciones disminuye la precisión del estimador, el peor escenario para el sistema eléctrico, se presenta ante la pérdida de la línea de Transmision 2-3, obteniendo el valor máximo de voltaje en p.u. de 1.178 y la máxima diferencia angular de 8.386 grados, entre los valores reales y los estimados.
Para el EEC, en los gráficos 9 y 10 se puede observar que, ante la salida de las mediciones asociadas a las líneas de transmisión, el estimador de estado presenta valores excesivos en los valores estimados, ya que la pérdida de mediciones disminuye la precisión del estimador, el peor escenario para el sistema eléctrico, se presenta ante la pérdida de la línea de Transmision 5-2, obteniendo el valor máximo de voltaje en p.u. de 1.43 y la máxima diferencia angular de 39.14 grados, entre los valores reales y los estimados. Además ante la pérdida de la línea de Transmision 4-2 el sistema no converge, es decir no tiene una solución óptima.
Ante la pérdida de las líneas de transmisión del SEP, se evidencia la disminución de precisión del vector de estado, resultado de los estimadores de estado, como se puede apreciar en las figuras. Obteniendo los peores escenarios para el EEC, ya que debido a la pérdida de las mediciones relacionada con las líneas de transmisión, el estimador de estado pierde su precisión y calidad en los valores estimados; el EE PMU presenta mejor precisión en comparación del EEC, ya que la pérdida de mediciones asociadas a las líneas de Transmision son compensadas con el resto de las mediciones provenientes de las PMU y el EEH mantiene su precisión y robustez ante la perdida de mediciones, siendo para el EEH, el que presenta los mejores resultados para el vector de estado estimado y la mejor opción de estimación de estado para el SEP, debido a que existen mediciones convencionales y fasoriales, que se complementan para mantener la precisión del mismo.
PMU | |||
Tipo de Estimador | Modelo IEEE | Tiempo (seg) | No. de iteraciones |
EEH | 14 nodos | 0.0044, 0.0264, 0.0047 | 9, 31, 9 |
30 nodos | 0.0313, 0.1201, 0.0213 | 13, 15, 14 | |
118 nodos | 0.2381, 0.3012, 0.2582 | 14, 15, 14 | |
EE PMU | 14 nodos | 0.001, 0.0328, 0.0025 | 8, 6, 7 |
30 nodos | 0.0327, 0.0457, 0.0311 | 12, 11, 13 | |
118 nodos | 0.2471, 0.2842, 0.3584 | 15, 15, 18 |
En la tabla 2 se visualiza, el número de iteraciones y el tiempo en segundos, que tarda el código para los EEH y EE PMU, debido a la pérdida de PMU en el sistema eléctrico, que generan los peores escenarios del vector de estado. Para los sistemas de la IEEE de prueba de 14, 30 y 118 nodos. El EEH presenta los mejores tiempo de respuesta y número de iteraciones en contraste con el EE PMU.
Es eminente el aumento de números de PMU en el sistema eléctrico de potencia que considera contingencias en contraste con el sistema que no considera contingencias, como la salida de una PMU o la pérdida de una línea de transmisión, sin embargo, esto ayuda a mejorar la precisión y calidad del estimador de estado. Las mediciones provenientes de las PMUs mantienen la observabilidad y redundancia del SEP.
Ante la pérdida de PMUs en el sistema eléctrico, se evidencia una reducción en la precisión de los EEH y EE PMU, siendo el EEH, el cual presenta una mejor respuesta con respecto al EE PMU, ya que el aporte de las mediciones convencionales mantiene la calidad del vector de estado tanto en magnitud y ángulo. Las mediciones convencionales, equivalen aproximadamente al número de mediciones fasoriales de las PMUs que se pierden ante una contingencia, motivo por el cual no deben ser desechadas en el proceso de estimación de estado híbrido. Ante este escenario de pérdida de PMUs en el SEP, el EEH presenta la mejor respuesta y calidad del vector de estado.
Con la perdida de las líneas de Transmision del SEP, se evidencia disminución considerable en la precisión de los EEH, EEPMU y EEC, siendo este último el que mayor número de iteraciones y tiempos de solución del algoritmo presenta, seguido del EE PMU, el EEH la mejor opción de estimación de estado para el SEP debido a que existen mediciones convencionales y fasoriales, que se complementan y mantienen la calidad y precisión del vector de estado estimado, además, su implementación en los actuales estimadores de estado convencionales, no requiere grandes cambios en el software existente.
El objetivo de los SEP a largo plazo es contar con el mayor número posible de PMUs, sin embargo, los costos asociados (gestión, compilación de gran cantidad de datos de los sincrofasores) para la instalación de PMU son elevados, motivo por el cual ante las contingencias que pueden presentarse en los SEP, como la pérdida de una PMU o la salida de una línea de transmisión, el papel de las mediciones convencionales es de gran importancia para el proceso de estimación de estado, ya que las mismas mantienen la precisión del estimador sin incrementar el número de iteraciones y el tiempo de respuesta del estimador de estado.
Al grupo de docentes de la Universidad Politécnica Salesiana, Maestría en Electricidad por compartir experiencias personales adquiridas a lo largo de su vida profesional y conocimientos, de manera especial, a Diego, tutor del proyecto de tesis, quien ha supervisado con paciencia y rectitud, este trabajo.
[1] S. Venugopal and R. Sunitha, “Robust Dynamic State EstimationDuring Missing of PMU Data,”2018 4th International Conferencefor Convergence in Technology (I2CT) SDMIT Ujire, Mangalore,India. Oct 27-28, 2018, pp. 1–6, 2020. [Online]. Available:doi:10.1109/i2ct42659.2018.9058140.
Published on 29/03/22
Accepted on 14/03/22
Submitted on 09/12/21
Volume 38, Issue 1, 2022
DOI: 10.23967/j.rimni.2022.03.010
Licence: CC BY-NC-SA license
Are you one of the authors of this document?