skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Complex-LASSO Approach for Localizing Forced Oscillations in Power Systems
We study the problem of localizing multiple sources of forced oscillations (FOs) and estimating their characteristics, such as frequency, phase, and amplitude, using noisy PMU measurements. For each source location, we model the input oscillation as a sum of unknown sinusoidal terms. This allows us to obtain a linear relationship between measurements and the inputs at the unknown sinusoids’ frequencies in the frequency domain. We determine these frequencies by thresholding the em- pirical spectrum of the noisy measurements. Assuming sparsity in the number of FOs’ locations and the number of sinusoids at each location, we cast the location recovery problem as an 1-norm regularized least squares problem in the complex domain—i.e., complex-LASSO (linear shrinkage and selection operator). We numerically solve this optimization problem using the complex- valued coordinate descent method, and show its efficiency on the IEEE 68-bus, 16 machine and WECC 179-bus, 29-machine systems.  more » « less
Award ID(s):
1934766
PAR ID:
10353956
Author(s) / Creator(s):
Date Published:
Journal Name:
2022 IEEE Power & Energy Society General Meeting (PESGM)
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper explores the use of changepoint detection (CPD) for an improved time-localization of forced oscillations (FOs) in measured power system data. In order for the autoregressive moving average plus sinusoids (ARMA+S) class of electromechanical mode meters to successfully estimate modal frequency and damping from data that contains a FO, accurate estimates of where the FO exists in time series are needed. Compared to the existing correlation-based method, the proposed CPD method is based on upon a maximum likelihood estimator (MLE) for the detection of an unknown number changes in signal mean to unknown levels at unknown times. Using the pruned exact linear time (PELT) dynamic programming algorithm along with a novel refinement technique, the proposed approach is shown to provide a dramatic improvement in FO start/stop time estimation accuracy while being robust to intermittent FOs. These findings were supported though simulations with the minniWECC model. 
    more » « less
  2. We consider the problem of finding the lowest order stable rational transfer function that interpolates a set of given noisy time and frequency domain data points. Our main result shows that exploiting results from rational interpolation theory allows for recasting this problem as minimizing the rank of a matrix constructed from the frequency domain data (the Loewner matrix) along with the Hankel matrix of time domain data, subject to a semidefinite constraint that enforces stability and consistency between the time and frequency domain data. These results are applied to a practical problem: identifying a system from noisy measurements of its time and frequency responses. The proposed method is able to obtain stable low order models using substantially smaller matrices than those reported earlier and consequently in a fraction of the computation time. 
    more » « less
  3. Abstract The inverse problem in gravimetry is to find a domain 𝐷 inside the reference domain Ω from boundary measurements of gravitational force outside Ω.We found that about five parameters of the unknown 𝐷 can be stably determined given data noise in practical situations.An ellipse is uniquely determined by five parameters.We prove uniqueness and stability of recovering an ellipse for the inverse problem from minimal amount of data which are the gravitational force at three boundary points.In the proofs, we derive and use simple systems of linear and nonlinear algebraic equations for natural parameters of an ellipse.To illustrate the technique, we use these equations in numerical examples with various location of measurements points on ∂ ⁡ Ω \partial\Omega .Similarly, a rectangular 𝐷 is considered.We consider the problem in the plane as a model for the three-dimensional problem due to simplicity. 
    more » « less
  4. State estimation (SE) is an important energy management system application for power system operations. Linear state estimation (LSE) is a variant of SE based on linear relationships between state variables and measurements. LSE estimates system state variables, including bus voltage magnitudes and angles in an electric power transmission network, using a network model derived from the topology processor and measurements. Phasor measurement units (PMUs) enable the implementation of LSE by providing synchronized high-speed measurements. However, as the size of the power system increases, the computational overhead of the state-of-the-art (SOTA) LSE grows exponentially, where the practical implementation of LSE is challenged. This paper presents a distributed linear state estimation (D-LSE) at the substation and area levels using a hierarchical transmission network topology processor (H-TNTP). The proposed substation-level and area-level D-LSE can efficiently and accurately estimate system state variables at the PMU rate, thus enhancing the estimation reliability and efficiency of modern power systems. Network-level LSE has been integrated with H-TNTP based on PMU measurements, thus enhancing the SOTA LSE and providing redundancy to substation-level and area-level D-LSE. The implementations of D-LSE and enhanced LSE have been investigated for two benchmark power systems, a modified two-area four-machine power system and the IEEE 68 bus power system, on a real-time digital simulator. The typical results indicate that the proposed multilevel D-LSE is efficient, resilient, and robust for topology changes, bad data, and noisy measurements compared to the SOTA LSE. 
    more » « less
  5. State estimation (SE) is an important energy management system application for power system operations. Linear state estimation (LSE) is a variant of SE based on linear relationships between state variables and measurements. LSE estimates system state variables, including bus voltage magnitudes and angles in an electric power transmission network, using a network model derived from the topology processor and measurements. Phasor measurement units (PMUs) enable the implementation of LSE by providing synchronized high-speed measurements. However, as the size of the power system increases, the computational overhead of the state-of-the-art (SOTA) LSE grows exponentially, where the practical implementation of LSE is challenged. This paper presents a distributed linear state estimation (D-LSE) at the substation and area levels using a hierarchical transmission network topology processor (H-TNTP). The proposed substation-level and area-level D-LSE can efficiently and accurately estimate system state variables at the PMU rate, thus enhancing the estimation reliability and efficiency of modern power systems. Network-level LSE has been integrated with H-TNTP based on PMU measurements, thus enhancing the SOTA LSE and providing redundancy to substation-level and area-level D-LSE. The implementations of D-LSE and enhanced LSE have been investigated for two benchmark power systems, a modified two-area four-machine power system and the IEEE 68 bus power system, on a real-time digital simulator. The typical results indicate that the proposed multilevel D-LSE is efficient, resilient, and robust for topology changes, bad data, and noisy measurements compared to the SOTA LSE. 
    more » « less