Key Takeaways & Executive Findings
- •• • China's renewable energy installed capacity reached 1.84 billion kW by end-2025, with projections to become the primary electricity supply source by 2030; this scale renders conventional electromechanical transient models inadequate because they neglect fast electromagnetic dynamics, directly driving the need for EMT eigenvalue analysis in systems with thousands of state variables. • • Park transformation fails under asymmetric operating conditions by introducing negative-sequence double-frequency components and cannot extend to single-phase or multi-phase systems; double Park transformation resolves positive- and negative-sequence components but leaves zero-sequence components as AC quantities, fundamentally limiting its applicability to unbalanced systems. • • Shifted frequency analysis (SFA) converts AC state variables into DC envelope signals under steady-state, enabling equilibrium-point linearization; however, when implemented with real matrices via time-scale transformation, the state matrix dimension doubles, increasing eigenvalue computation burden by approximately 2x compared to complex-space formulations. • • Floquet-theory-based trajectory linearization eliminates time-periodicity by eigendecomposing the state transition matrix Φ(T,0), but the computational cost of this eigendecomposition scales as O(n³) for an n-dimensional system, making it impractical for large-scale MMC-based systems with thousands of states without sparse or partial eigenvalue techniques.
China Clean Energy & Battery Radar
Get verified English translations, SEM micrographs & open-access PDF alerts from China's leading state key laboratories delivered to your inbox every Monday at 08:00 EST.
Abstract
The escalating 'double-high' penetration of power electronics in new-type power systems has rendered conventional electromechanical transient small-signal stability analysis inadequate, necessitating electromagnetic transient (EMT) small-signal stability assessment. Eigenvalue analysis, grounded in rigorous theoretical foundations, is widely applied but faces two critical bottlenecks when scaled to large systems: the construction of EMT linearized state-space equations and the solution of high-order state-matrix eigenvalues. This review examines the urgent demand for EMT model eigenvalue analysis in large-scale new-type power systems. It systematically surveys research status and challenges across three domains: equilibrium-point modeling of EMT models, linearized state-space modeling, and efficient computation of critical eigenvalues. For equilibrium-point modeling, the paper evaluates Park transformation, double Park transformation, shifted frequency analysis (SFA), time-scale transformation, and Floquet-theory-based trajectory linearization. For linearized state-space modeling, it assesses component-level and network-level linearization strategies. For eigenvalue computation, it reviews partial eigenvalue algorithms, sparse matrix techniques, and model-order reduction methods. Key challenges include the inability of Park transformation to handle asymmetric or single-phase systems, the computational burden of Floquet transition matrix eigendecomposition, and the poor scalability of dense eigenvalue solvers for systems exceeding thousands of states. The paper concludes by identifying future research directions, including structure-preserving linearization, GPU-accelerated sparse eigenvalue algorithms, and data-driven model reduction, to enable practical EMT small-signal stability analysis for systems with 18.4 billion kW of installed renewable capacity by 2030.
1. Introduction
Conventional electromechanical transient small-signal stability analysis, built on quasi-steady-state phasor assumptions, has become structurally inadequate for new-type power systems dominated by power electronic converters. The 'double-high' characteristic—high penetration of renewables and high proportion of power electronics—introduces wideband oscillations spanning multiple frequency bands with time-varying characteristics and broad propagation ranges. These oscillations originate from small-signal instabilities that electromechanical models cannot capture because they neglect fast electromagnetic transients. While impedance analysis, open-loop modal analysis, complex torque coefficient methods, and time-domain simulation each offer partial insights, none provides the complete global state characterization and rigorous theoretical foundation of eigenvalue analysis. The bottleneck is not conceptual but computational: constructing EMT linearized state-space equations for large-scale systems is difficult because EMT models use instantaneous values, and solving eigenvalues of high-order state matrices is prohibitively expensive.
This review addresses the three-stage pipeline required for practical EMT eigenvalue analysis. First, equilibrium-point modeling must convert periodic AC state variables into DC quantities; existing approaches—Park transformation, double Park transformation, SFA, time-scale transformation, and Floquet-based trajectory linearization—each carry specific limitations under asymmetric, single-phase, or multi-frequency conditions. Second, linearized state-space construction must assemble component-level and network-level equations without introducing spurious dynamics or losing structural sparsity. Third, eigenvalue computation must extract critical modes from matrices with thousands to tens of thousands of states. The paper systematically evaluates research status and challenges across these three domains, then identifies future breakthroughs including structure-preserving formulations, sparse partial eigenvalue algorithms, and GPU-accelerated solvers to enable EMT small-signal stability analysis at the scale of China's 1.84 billion kW renewable fleet.
Loading authentic research manuscript (Pages 1–5)...
WANG Yuhong, SU Miaohong, GAO Shilin, YE Hua, CHEN Ying (2026). Eigenvalue Computation Techniques for Small-Signal Stability Analysis of Large-Scale New-Type Power Systems: A Review and Outlook. Power Automation Equipment. https://doi.org/10.16081/j.epae.202606020
Research & Educational Purpose Only: The translations, structured abstracts, analytical annotations, and data reports provided by SinoGreenTechare intended exclusively for academic research, internal corporate R&D, and educational benchmarking. They do not constitute formal engineering, chemical safety, legal, or professional advice.
Copyright & Intellectual Property Notice: Original copyright of the underlying source articles and experimental data remains with the respective authors, institutions, and original publishing journals. SinoGreenTech claims intellectual property only over its proprietary translations, analytical syntheses, and AEO structured enhancements in accordance with international fair use and academic citation principles.
Frequently Asked Questions
What specific failure mechanism prevents Park transformation from being used for equilibrium-point modeling in asymmetric or single-phase power systems?
Park transformation maps three-phase instantaneous quantities to a synchronously rotating dq0 reference frame, converting fundamental-frequency AC signals into DC quantities under ideal symmetric conditions. Under asymmetric operation, the transformation introduces negative-sequence double-frequency components that remain as AC quantities in the dq0 frame, preventing equilibrium-point linearization. Furthermore, Park transformation is mathematically defined only for three-phase systems; it cannot be extended to single-phase or multi-phase systems without modification. Double Park transformation partially resolves this by mapping positive- and negative-sequence components to separate dq+ and dq− frames, but zero-sequence components remain as AC quantities, leaving the equilibrium-point problem unsolved for unbalanced systems with zero-sequence paths.
What is the computational complexity bottleneck of Floquet-theory-based trajectory linearization for large-scale MMC systems, and how does it compare to SFA?
Floquet-based trajectory linearization requires eigendecomposition of the state transition matrix Φ(T,0), which scales as O(n³) for an n-dimensional system. For an MMC-based system with thousands of states, this becomes computationally prohibitive. Additionally, obtaining the periodic trajectory requires time-domain simulation over at least one fundamental period, adding simulation cost. In contrast, SFA converts AC variables to DC envelope signals through frequency shifting, enabling direct equilibrium-point linearization without eigendecomposition. However, SFA with real-matrix time-scale transformation doubles the state dimension, increasing eigenvalue computation cost by approximately 2x. The trade-off is between Floquet's O(n³) eigendecomposition and SFA's 2x dimension expansion, with neither approach providing a clear scalable solution for systems exceeding 10,000 states.
Why can't existing eigenvalue solvers handle high-order state matrices from large-scale EMT models, and what specific numerical failure modes occur?
Dense eigenvalue solvers such as QR iteration scale as O(n³) per iteration and require O(n²) memory, making them infeasible for state matrices with n > 10,000. For EMT models of large-scale systems with thousands of power electronic converters, the state matrix dimension can reach tens of thousands. Sparse eigenvalue algorithms such as Arnoldi and Lanczos iterations can extract a subset of critical eigenvalues, but they suffer from slow convergence when eigenvalues are clustered near the imaginary axis—precisely the region of interest for small-signal stability. Additionally, the state matrix may be ill-conditioned due to disparate time constants spanning from microsecond electromagnetic transients to second-level control dynamics, causing numerical instability in eigenvalue computation. Model-order reduction techniques can reduce dimension but may introduce spurious modes or miss critical eigenvalues if the reduction basis is poorly chosen.
What are the practical scalability bottlenecks for implementing EMT eigenvalue analysis on China's 1.84 billion kW renewable fleet by 2030?
Three bottlenecks dominate. First, model construction: each power electronic device (VSC, MMC, DFIG) requires detailed EMT modeling with switching-level or averaged models, and a system with thousands of converters generates state matrices with tens of thousands of states. Second, equilibrium-point computation: finding the periodic steady-state trajectory for a system with 18.4 billion kW of renewables requires time-domain simulation over multiple fundamental periods, which is computationally expensive for systems with high penetration of power electronics. Third, eigenvalue extraction: critical modes near the imaginary axis require high-resolution computation, but sparse solvers struggle with clustered eigenvalues. Current research directions include structure-preserving linearization to maintain sparsity, GPU-accelerated sparse eigenvalue algorithms to reduce computation time, and data-driven model reduction to decrease state dimension while preserving critical modes.
How does the time-scale transformation variant of SFA affect eigenvalue computation accuracy and cost compared to complex-space SFA?
Complex-space SFA represents state variables as complex envelopes, preserving the original state dimension n but requiring complex arithmetic. Time-scale transformation represents the same variables using real matrices, which doubles the state dimension to 2n but allows standard real-valued eigenvalue solvers. The 2n dimension increases computation cost by approximately 2x for dense solvers (O((2n)³) = 8O(n³) vs. O(n³) for complex arithmetic, though complex arithmetic itself costs roughly 4x real arithmetic, making the comparison nuanced). Accuracy is theoretically equivalent, but real-matrix formulations may exhibit different numerical conditioning. For large-scale systems where sparse solvers are necessary, the doubled dimension increases memory requirements and may slow convergence. The choice between complex-space and real-matrix SFA depends on available solver libraries and the specific sparsity structure of the state matrix.
Related Chinese Research & Cross-Citations
Impedance Modeling and Parameter Optimization Method for Wind Farms Considering Station-Level Control
Large-scale wind farms integrated into weak grids are susceptible to broadband oscillations, a problem that existing impedance models and control parameter optimization methods fail to address systematically because they neglect station-level control and frequency coupling effects. This paper proposes a station-level control strategy based on available power allocation and an adaptive compass search (ACS) algorithm for optimizing station control parameters. A sequence impedance method incorporating frequency coupling effects establishes the aggregated wind farm impedance, and the grid-connected multiple-input multiple-output (MIMO) system is decoupled into positive- and negative-sequence single-input single-output (SISO) impedance models. A high-precision wind farm impedance model that accounts for station-level control is constructed, and the Nyquist criterion evaluates the suppression effect of station control on broadband oscillations. The ACS algorithm optimizes station control parameters to enhance the adaptability of wind farm impedance to external grid impedance, thereby reducing oscillation risk. RT-LAB platform simulations validate the impedance modeling method and control parameter optimization. Results demonstrate that considering station-level control yields a high-precision wind farm impedance model. Compared with genetic algorithm (GA), particle swarm optimization (PSO), grey wolf optimizer (GWO), and sparrow search algorithm (SSA), ACS is more suitable for optimizing wind farm control parameters. Under specific weak-grid conditions, ACS-optimized station control parameters effectively improve grid-connected stability of the wind farm.
Active and Reactive Power Coordinated Voltage Support Control Method for High-Inertia Energy-Storage Synchronous Condenser
The high-inertia energy-storage synchronous condenser (SC-HI-ES) can provide reactive power support by adjusting excitation current and active power support by actively varying rotor speed. However, the absence of a prime mover makes its control capability difficult to quantify, preventing coordinated active and reactive power support for grid voltage. This paper analyzes the factors influencing SC-HI-ES power control capability, derives rotor current equations considering stator and rotor circuit constraints, and establishes the power controllable range (PCR) under rotor current constraints and speed variations. The relationship between the PCR and voltage support power demand is parsed, and a coordinated active-reactive power voltage support control method is proposed. Case studies validate the method's effectiveness. Results show that the proposed method enhances voltage support and avoids rotor current over-limit. Compared with methods ignoring speed variation or controlling only reactive power, the proposed method dynamically calculates PCR-based power demand under current speed, and through coordinated active and reactive power control, maximally satisfies the point of common coupling voltage target while ensuring rotor current does not exceed allowable values. This provides effective technical support for SC-HI-ES application in new-type power systems. Frequency control capability characterization and frequency-coordinated control are identified as future research directions.
Hardware-in-the-Loop Test Platform for Photovoltaic High-Frequency Controllers Based on UREP + FPGA
Real-time constraints prevent CPU-based electromagnetic transient simulators from accurately sampling high-frequency PWM signals at microsecond step sizes, leading to insufficient HIL test precision and even instability of the closed-loop system. This paper introduces a PWM averaging technique to construct a novel UREP (electromagnetic transient real-time simulator) + FPGA hardware-in-the-loop test platform. To achieve accurate PWM acquisition, the FPGA's nanosecond-level clock resolution increases the number of sampling points within a single switching period from a few points under microsecond simulation steps to several thousand points. The FPGA implements fixed-step PWM averaging, converting discrete switching states into continuous duty cycles that serve as control signals for the UREP-side inverter average model. This technique avoids frequent topology updates and state matrix reconstructions triggered by switching events, reduces high real-time computational resource consumption, improves numerical stability of the in-loop system, and enhances HIL test accuracy under large-step real-time constraints. Simulation cases and industrial-grade controller HIL tests verify the platform's feasibility and accuracy. The platform provides an efficient, reproducible verification scheme for photovoltaic controller control strategies and offers a feasible technical path for domestic electromagnetic transient real-time simulation platforms to conduct high-frequency controller HIL tests.
Design and Implementation of a Distributed Simulation Architecture for Modern Power Systems Based on Data Distribution Service
The integration of high-penetration renewable energy and power electronic converters has intensified the dynamic complexity of modern power systems, imposing stringent demands on simulation accuracy and computational efficiency. This paper proposes a distributed simulation architecture based on Data Distribution Service (DDS) that leverages dispersed computing resources to enhance scalability and efficiency while preserving numerical fidelity. The power system model is mathematically decoupled into multiple independently solvable basic subsystems, and a generic data transmission interface is designed using DDS. A time-consumption balancing scheme groups and distributes these subsystems across multiple devices, and a two-layer synchronization strategy enables efficient parallel simulation. The architecture is validated on a two-area four-machine system, the IEEE New England 10-machine 39-bus system, and the WECC 29-machine 179-bus system. Compared with centralized simulation, the average error on the two-area four-machine system remains below 1.05%. Simulation efficiency improvement reaches approximately 10% on the 10-machine 39-bus system and about 48% on the 29-machine 179-bus system. The results confirm high accuracy across different system scales and demonstrate that efficiency gains become more pronounced as system size increases, validating the architecture's scalability and compatibility. The proposed framework offers a promising pathway for large-scale power system simulation and supports future integration with edge computing and cross-platform deployment.
Heterogeneous Weighted Graph Partitioning and Decoupling Optimization Strategy for Electromagnetic Transient Parallel Simulation of AC/DC Distribution Networks
The increasing penetration of power electronic devices in AC/DC distribution networks imposes stringent computational demands on electromagnetic transient (EMT) parallel simulation. Conventional transmission-line delay decoupling methods are ill-suited to the strong electrical coupling characteristic of such networks. This paper proposes a non-delay decoupling parallel simulation acceleration framework based on heterogeneous weighted graph partitioning. An empirical computational cost evaluation model for each parallel decoupling stage is established, and a heterogeneous weighted graph model is constructed to precisely characterize the simulation computational complexity of AC/DC distribution network components, mapping matrix dimensions of device mathematical models to graph node weights. A multi-objective graph partitioning scheme is formulated that simultaneously balances partition computational overhead and minimizes the number of tie-line variables, complemented by an optimal partition number screening strategy. Simulation validation is conducted on three large-scale AC/DC distribution network composite test cases: IEEE 34-node, IEEE 123-node, and European Low Voltage (European LV) systems, all retrofitted with DC sections. Results demonstrate that the proposed empirical computational cost model achieves a simulation time fitting goodness-of-fit R² > 0.97, indicating high predictive accuracy. Under optimal partition configuration, the proposed method attains parallel speedup ratios of 12.11–15.77, significantly outperforming conventional natural partitioning schemes and effectively enhancing the EMT parallel simulation efficiency of AC/DC distribution networks.
Hierarchical Hybrid Non-Delay Decoupling Parallel Method for Electromagnetic Transient Simulation of DC-Collector Offshore Wind Farms
Real-time electromagnetic transient (EMT) simulation of DC-collector offshore wind farms is constrained by microsecond time steps, high model order from cascaded power electronic converters, and the inability of conventional decoupling methods to handle complex series-parallel topologies without introducing artificial delays. Existing non-delay decoupling methods, such as multi-area Thevenin equivalence (MATE) and compensation method, rely on branch tearing and are ill-suited for systems with numerous common-bus partitions, leading to excessive link variables and singular admittance matrices. This paper proposes a hierarchical hybrid non-delay decoupling parallel method that integrates MATE and its dual (node-tearing) approach through a layered architecture. The method constructs a mixed equivalent model tailored to DC-collector offshore wind farms, enabling flexible selection between unified and hierarchical solution modes for link variables based on operating conditions. A complete non-delay decoupling simulation workflow is established and validated on a DC series-parallel grid-connected offshore wind farm test case implemented in MATLAB. Results demonstrate that the proposed method reduces solution matrix dimensions and achieves significant speedup without compromising accuracy. For a 96-turbine wind farm, the hybrid decoupling model achieves an 11.99x speedup over the detailed model, compared to 8.77x for series-only decoupling and 6.22x for parallel-only decoupling. Mean absolute errors (MAE) for key variables remain below 0.1011, and root mean square errors (RMSE) below 0.6318, confirming high fidelity. The method enhances parallel simulation performance and offers a generalizable solution for real-time EMT simulation of large-scale offshore wind farms.