Exploring the impact of wind penetration on power system equilibrium using a numerical continuation approach

In this paper we investigate how the equilibrium characteristics of conventional power systems may change with an increase in wind penetration. We first derive a differentialalgebraic model of a power system network consisting of synchronous generators, loads and a wind power plant modeled by wind turbines and doubly-fed induction generators (DFIG). The models of these three components are coupled via nonlinear power flow equations. In contrast to the traditional approach for solving the power flows via iterative methods that often lead to a local solution, we apply a recently developed parameterhomotopy based numerical continuation algorithm to compute all possible solutions. The method obtains all the solutions of the power flow equations over multiple values of the wind penetration. We observe that depending on the penetration limit and the setpoint value for the magnitude of the wind bus voltage, the system may exhibit several undesired or even unstable equilibria. We illustrate these results through a detailed simulation of a 5-machine power system model with wind injection, and highlight how the solutions may be helpful for small-signal stability assessment.

[1]  Andrew J. Sommese,et al.  The numerical solution of systems of polynomials - arising in engineering and science , 2005 .

[2]  Anjan Bose,et al.  Smart Transmission Grid Applications and Their Supporting Infrastructure , 2010, IEEE Transactions on Smart Grid.

[3]  James S. Thorp,et al.  An efficient algorithm to locate all the load flow solutions , 1993 .

[4]  Konstantin Turitsyn,et al.  Numerical polynomial homotopy continuation method to locate all the power flow solutions , 2014, 1408.2732.

[5]  B. C. Lesieutre,et al.  Counterexample to a Continuation-Based Algorithm for Finding All Power Flow Solutions , 2013, IEEE Transactions on Power Systems.

[6]  A. Morgan,et al.  Coefficient-parameter polynomial continuation , 1989 .

[7]  Bruno Buchberger,et al.  Bruno Buchberger's PhD thesis 1965: An algorithm for finding the basis elements of the residue class ring of a zero dimensional polynomial ideal , 2006, J. Symb. Comput..

[8]  L. Ni,et al.  Parallel processing for the load flow of power systems: the approach and applications , 1989, Proceedings of the 28th IEEE Conference on Decision and Control,.

[9]  N. Jenkins,et al.  State-Space Modeling of Wind Turbine Generators for Power System Studies , 2013, IEEE Transactions on Industry Applications.

[10]  Thomas J. Overbye,et al.  Computation of a practical method to restore power flow solvability , 1995 .

[11]  I. Hiskens,et al.  Exploring the Power Flow Solution Space Boundary , 2001, IEEE Power Engineering Review.

[12]  J. Yorke,et al.  The cheater's homotopy: an efficient procedure for solving systems of polynomial equations , 1989 .

[13]  Konstantin S. Turitsyn,et al.  Appearance of multiple stable load flow solutions under power flow reversal conditions , 2014, 2014 IEEE PES General Meeting | Conference & Exposition.

[14]  Manfred Morari,et al.  Numerical Algebraic Geometry for Optimal Control Applications , 2011, SIAM J. Optim..

[15]  Hsiao-Dong Chiang,et al.  On the existence and uniqueness of load flow solution for radial distribution power networks , 1990 .

[16]  Antonio Montes,et al.  A New Algorithm for Discussing Gröbner Bases with Parameters , 2002, J. Symb. Comput..

[17]  David A. Cox,et al.  Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, 3/e (Undergraduate Texts in Mathematics) , 2007 .

[18]  Babu Narayanan,et al.  POWER SYSTEM STABILITY AND CONTROL , 2015 .