pith. sign in

arxiv: 2605.14099 · v1 · pith:3VOSNDOWnew · submitted 2026-05-13 · 📡 eess.SY · cs.SY

Frequency Nadir-Constrained Power System Restoration Planning with Energy Storage

Xiangyu Zou , Amir Reza Nikzad , Ilyas Farhat , John W. Simpson-Porco This is my paper

Pith reviewed 2026-05-15 05:15 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords power system restorationfrequency nadirenergy storage systemsMILP optimizationblack-start planningfrequency stabilitytransmission systems
0
0 comments X

The pith

A frequency-constrained MILP framework plans black-start restoration while using energy storage to keep frequency deviations safe and shorten recovery time.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Power system operators must restore electricity after blackouts without triggering dangerous frequency swings that could damage equipment or cause further collapse. This paper builds an optimization model that sequences the reconnection of generators, loads, and energy storage units while predicting the lowest frequency each action will produce. A new prediction method tailored to energy storage integration supplies the constraints that keep those frequency drops inside safe bounds. The resulting plans are tested on a modified IEEE 9-bus system in both MATLAB and PSS/E, confirming that frequency security holds while total restoration time decreases. Readers care because shorter, more secure recovery directly reduces economic losses and public inconvenience after major outages.

Core claim

The paper establishes a multiperiod mixed-integer linear programming model for transmission-system black-start restoration that embeds a frequency nadir prediction method developed specifically for systems containing energy storage, thereby enforcing safe frequency limits on every restorative action and demonstrating faster overall recovery through storage coordination.

What carries the argument

Frequency nadir prediction method for ESS-integrated systems, embedded as constraints inside the multiperiod MILP restoration planning framework.

If this is right

  • Restoration sequences can be computed that coordinate energy storage to reduce total recovery duration.
  • Frequency deviations remain inside prescribed safe limits throughout the plan, as confirmed by MATLAB and PSS/E simulations.
  • The formulation applies directly to systems that combine synchronous machines with energy storage units.
  • Black-start plans become feasible even when early restorative actions would otherwise risk excessive frequency drops.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Extending the same nadir-prediction approach to include other dynamic limits such as voltage or rotor-angle stability could produce more complete restoration models.
  • The MILP structure may support real-time re-optimization if updated system measurements become available during an actual event.
  • Testing the method on larger transmission networks would reveal whether solution times remain practical for operational use.

Load-bearing premise

The frequency nadir prediction method accurately estimates the minimum frequency reached after each restorative action when energy storage is present.

What would settle it

A full dynamic simulation of the computed restoration plan on the modified IEEE 9-bus system that produces a frequency nadir below the safe limit enforced by the optimization.

Figures

Figures reproduced from arXiv: 2605.14099 by Amir Reza Nikzad, Ilyas Farhat, John W. Simpson-Porco, Xiangyu Zou.

Figure 1
Figure 1. Figure 1: Four phases of NBSU start-up. To capture the start-up phase of each generator in the MILP framework, auxiliary binary variables bgc, bgr, bgo ∈ {0, 1} G×T are introduced, where b i gc[k] = 1 if generator i is in the cranking phase at step k, b i gr[k] = 1 if it is ramping, and b i go[k] = 1 if it is online; all three are zero when the generator is offline. Since the cranking and ramping durations are known… view at source ↗
Figure 2
Figure 2. Figure 2: ESS grid interconnection. The state of charge (energy stored) is denoted by Es ∈ R S×T , with initial value Es,0 ∈ R S, and is bounded by the maximum storage capacity Es ∈ R S, i.e., 0S ≤ Es[k] ≤ Es, k ∈ {1, . . . , T}. (11) The power flow to and from the battery is decom￾posed into nonnegative charging and discharging components P in s , P out s ∈ R S×T , where P in s [k], P out s [k] ≥ 0S. Accounting for… view at source ↗
Figure 4
Figure 4. Figure 4: The model consists of a governor that detects frequency [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: Response of frequency response following a disturbance. [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Closed-loop system model with PFR and ESS support. [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Open-loop system model under the ramp approximation. [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Rolling-horizon optimization restoration planning. At each iteration G2 [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Modified IEEE 9-bus system. Numbers in brackets denote the indices [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 12
Figure 12. Figure 12: ESS power setpoints and state of charge during restoration under [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Center-of-inertia frequency profiles (PSS/E simulation) for proposed [PITH_FULL_IMAGE:figures/full_fig_p009_13.png] view at source ↗
read the original abstract

Power system restoration following blackouts must ensure frequency stability throughout the recovery process. This paper proposes a frequency-constrained mixed-integer linear programming (MILP) framework for black-start restoration planning in transmission systems with synchronous machines and energy storage systems. To prevent excessive frequency deviations caused by restorative actions, a frequency nadir prediction method is developed for power systems with energy storage system (ESS) integration and incorporated into a multiperiod optimization framework. The formulation ensures that frequency deviations resulting from restorative actions remain within prescribed safe limits. Furthermore, the presented framework leverages ESSs to enhance frequency security and recovery speed. Case studies on a modified IEEE 9-bus system demonstrate that the computed restoration plan maintains frequency security, as validated through MATLAB and PSS/E simulations, while reducing restoration time through ESS coordination.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a frequency-nadir-constrained MILP framework for transmission-system black-start restoration planning that incorporates a custom nadir-prediction method for systems containing both synchronous machines and energy storage. The method is embedded as linear constraints in a multi-period optimization model whose objective is to minimize restoration time while keeping frequency deviations within safe limits; the resulting schedule is validated on a modified IEEE 9-bus system via MATLAB and PSS/E time-domain simulations.

Significance. If the embedded nadir predictor is shown to be sufficiently accurate and non-optimistic, the framework would provide a practical, optimization-based tool for operators to coordinate ESS during restoration, potentially shortening recovery times without compromising frequency security—an increasingly relevant capability as inverter-based resources and storage proliferate.

major comments (2)
  1. [Frequency Nadir Prediction Method] §4 (or the section presenting the nadir-prediction method): the linearization/approximation of the swing-equation dynamics plus ESS droop/inertia response is load-bearing for the central claim, yet no explicit error bound, validation against nonlinear simulation at each restorative step, or sensitivity analysis to load-pickup timing and multi-machine synchronization is provided. The IEEE 9-bus PSS/E check only confirms the final schedule, not whether the predictor was tight or systematically optimistic on the critical nadir events.
  2. [Case Studies] §5 (Case Studies): the reported restoration-time reduction and frequency-security claim rest on a single modified IEEE 9-bus topology; without additional test systems (e.g., larger transmission networks or systems with higher ESS penetration) or a direct comparison of predicted versus simulated nadir at every restorative action, it is unclear whether the MILP constraints remain valid under realistic sequential dynamics.
minor comments (2)
  1. [Notation] Notation for ESS state-of-charge limits and droop coefficients should be defined once and used consistently across the prediction equations and the MILP constraints.
  2. [Figures] Figure captions and axis labels in the simulation results should explicitly state whether the plotted frequency traces are from the embedded predictor or from the full nonlinear PSS/E model.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major point below, indicating where revisions will be made to strengthen the validation of the nadir predictor and the case-study evidence.

read point-by-point responses

  1. Referee: [Frequency Nadir Prediction Method] §4 (or the section presenting the nadir-prediction method): the linearization/approximation of the swing-equation dynamics plus ESS droop/inertia response is load-bearing for the central claim, yet no explicit error bound, validation against nonlinear simulation at each restorative step, or sensitivity analysis to load-pickup timing and multi-machine synchronization is provided. The IEEE 9-bus PSS/E check only confirms the final schedule, not whether the predictor was tight or systematically optimistic on the critical nadir events.

    Authors: We agree that explicit error quantification and step-wise validation would strengthen the central claim. In the revised manuscript we will add (i) an error-bound analysis comparing the linear nadir predictor against full nonlinear swing-equation solutions over a range of load-pickup magnitudes and timings, (ii) a sensitivity study with respect to synchronization delays and multi-machine interactions, and (iii) a table reporting predicted versus simulated nadir values at every restorative action. These additions will demonstrate that the embedded constraints are neither systematically optimistic nor overly conservative. revision: yes

  2. Referee: [Case Studies] §5 (Case Studies): the reported restoration-time reduction and frequency-security claim rest on a single modified IEEE 9-bus topology; without additional test systems (e.g., larger transmission networks or systems with higher ESS penetration) or a direct comparison of predicted versus simulated nadir at every restorative action, it is unclear whether the MILP constraints remain valid under realistic sequential dynamics.

    Authors: The IEEE 9-bus system is a standard, well-documented benchmark that permits detailed PSS/E validation; the single-system study is therefore sufficient to illustrate the framework’s core mechanics and frequency-security guarantees. Nevertheless, we acknowledge the value of per-step nadir comparisons. In revision we will insert these comparisons (as noted in the response to the first comment) and will add a brief discussion of scalability to larger networks in the conclusions. We do not add new test systems at this stage because the computational and modeling effort would exceed the scope of a major revision, but the requested per-action validation will be provided. revision: partial

Circularity Check

0 steps flagged

No circularity: nadir predictor derived independently and validated externally

full rationale

The paper develops a frequency nadir prediction method from system dynamics (swing equation plus ESS droop/inertia) and embeds the resulting linear constraints into the MILP restoration planner. This is a standard forward modeling step, not a self-referential fit. The IEEE 9-bus case studies then validate the full schedule via independent MATLAB and PSS/E time-domain simulations, providing external falsifiability. No self-citation chain, no fitted parameter renamed as prediction, and no reduction of the central claim to its own inputs by construction appear in the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the framework appears to build on standard power system models and MILP techniques.

pith-pipeline@v0.9.0 · 5438 in / 935 out tokens · 31275 ms · 2026-05-15T05:15:35.023400+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    Electricity system resilience in a world of increased climate change and cybersecurity risk,

    E. L. Ratnam, K. G. Baldwin, P. Mancarella, M. Howden, and L. See- beck, “Electricity system resilience in a world of increased climate change and cybersecurity risk,”The Electricity Journal, vol. 33, no. 9, p. 106833, 2020

    work page 2020

  2. [2]

    Shedding light on the economic costs of long-duration power outages: A review of resilience assessment methods and strategies,

    M. Macmillan, K. Wilson, S. Baik, J. P. Carvallo, A. Dubey, and C. A. Holland, “Shedding light on the economic costs of long-duration power outages: A review of resilience assessment methods and strategies,” Energy Research & Social Science, vol. 99, p. 103055, 2023

    work page 2023

  3. [3]

    Technical analysis of the august 14, 2003, blackout,

    U.S.-Canada Power System Outage Task Force, “Technical analysis of the august 14, 2003, blackout,” NERC, Tech. Rep., 2004

    work page 2003

  4. [4]

    System restoration from blackstart resources,

    North American Electric Reliability Corporation (NERC), “System restoration from blackstart resources,” NERC, Atlanta, GA, USA, Tech. Rep. EOP-005-3, 2018

    work page 2018

  5. [5]

    U. G. Knight,Power Systems in Emergencies: From Contingency Plan- ning to Crisis Management. Wiley, 2013

    work page 2013

  6. [6]

    An integrated approach for power system restoration planning,

    F. Qiu and P. Li, “An integrated approach for power system restoration planning,”Proceedings of the IEEE, vol. 105, no. 7, pp. 1234–1252, 2017

    work page 2017

  7. [7]

    Part 7.8: Ontario power system restoration plan,

    Independent Electricity System Operator (IESO), “Part 7.8: Ontario power system restoration plan,” IESO, Toronto, ON, Canada, Tech. Rep., 2022

    work page 2022

  8. [8]

    Anderson and A

    P. Anderson and A. Fouad,Power System Control and Stability. Wiley- IEEE Press, 2003

    work page 2003

  9. [9]

    Inertia and the power grid: A guide without the spin,

    P. Denholm, T. Mai, R. W. Kenyon, B. Kroposki, and M. O’Malley, “Inertia and the power grid: A guide without the spin,” National Renewable Energy Lab.(NREL), Golden, CO, US, Tech. Rep., 2020

    work page 2020

  10. [10]

    MISO power system restoration plan,

    Midcontinent Independent System Operator (MISO), “MISO power system restoration plan,” MISO, Carmel, IN, United States, Tech. Rep., 2025

    work page 2025

  11. [11]

    Power system restoration - a task force report,

    M. Adibi, P. Clelland, L. Fink, H. Happ, R. Kafka, J. Raine, D. Scheurer, and F. Trefny, “Power system restoration - a task force report,”IEEE Transactions on Power Systems, vol. 2, no. 2, pp. 271–277, 1987

    work page 1987

  12. [12]

    New approaches in power system restoration,

    M. Adibi, L. Fink, J. Giri, D. Kirschen, S. Shahidehpour, and J. Zaborszky, “New approaches in power system restoration,”IEEE Transactions on Power Systems, vol. 7, no. 4, pp. 1428–1434, 1992

    work page 1992

  13. [13]

    Power system restoration: a literature review from 2006 to 2016,

    Y . Liu, R. Fan, and V . Terzija, “Power system restoration: a literature review from 2006 to 2016,”Journal of Modern Power Systems and Clean Energy, vol. 4, no. 3, pp. 332–341, 2016

    work page 2006

  14. [14]

    Transmission system restoration with co-optimization of repairs, load pickups, and generation dispatch,

    C. Coffrin and P. Van Hentenryck, “Transmission system restoration with co-optimization of repairs, load pickups, and generation dispatch,” International Journal of Electrical Power & Energy Systems, vol. 72, pp. 144–154, 2015

    work page 2015

  15. [15]

    A milp-based power system parallel restoration model with the support of mobile energy storage systems,

    Y . Xie, S. Cai, J. Wang, Y . Chen, and C. Qiu, “A milp-based power system parallel restoration model with the support of mobile energy storage systems,”Electric Power Systems Research, vol. 234, p. 110592, 2024

    work page 2024

  16. [16]

    Transmission system repair and restoration,

    P. Van Hentenryck and C. Coffrin, “Transmission system repair and restoration,”Mathematical Programming, vol. 151, no. 1, pp. 347–373, Jun. 2015

    work page 2015

  17. [17]

    Recursive restoration refinement: A fast heuristic for near-optimal restoration prioritization in power systems,

    N. Rhodes, C. Coffrin, and L. Roald, “Recursive restoration refinement: A fast heuristic for near-optimal restoration prioritization in power systems,”Electric Power Systems Research, vol. 212, p. 108454, 2022

    work page 2022

  18. [18]

    Feedback control approaches for restoration of power grids from blackouts,

    J. M. Miller, H. N. V . Pico, I. Dobson, A. Bernstein, and B. Cui, “Feedback control approaches for restoration of power grids from blackouts,”Electric Power Systems Research, vol. 211, p. 108414, 2022

    work page 2022

  19. [19]

    Service restoration using deep reinforcement learning and dynamic microgrid formation in distribution networks,

    M. A. Igder and X. Liang, “Service restoration using deep reinforcement learning and dynamic microgrid formation in distribution networks,” IEEE Transactions on Industry Applications, vol. 59, no. 5, pp. 5453– 5472, 2023

    work page 2023

  20. [20]

    Automatic power system restoration with inrush current estimation for industrial facility,

    A. Papasani, K. Zia, and W.-J. Lee, “Automatic power system restoration with inrush current estimation for industrial facility,”IEEE Transactions on Industry Applications, vol. 57, no. 6, pp. 5772–5781, 2021

    work page 2021

  21. [21]

    Review of restoration technology for renewable-dominated electric power systems,

    C. Chen, H. Liang, X. Zhai, J. Zhang, S. Liu, Z. Lin, and L. Yang, “Review of restoration technology for renewable-dominated electric power systems,”Energy Conversion and Economics, vol. 3, no. 5, pp. 287–303, 2022

    work page 2022

  22. [22]

    The future of power system restoration: Using distributed energy resources as a force to get back online,

    M. Braun, J. Brombach, C. Hachmann, D. Lafferte, A. Klingmann, W. Heckmann, F. Welck, D. Lohmeier, and H. Becker, “The future of power system restoration: Using distributed energy resources as a force to get back online,”IEEE Power and Energy Magazine, vol. 16, no. 6, pp. 30–41, 2018

    work page 2018

  23. [23]

    Advanced techniques of power system restoration and practical applications in transmission grids,

    D. Sharma, C. Lin, X. Luo, D. Wu, K. Thulasiraman, and J. N. Jiang, “Advanced techniques of power system restoration and practical applications in transmission grids,”Electric Power Systems Research, vol. 182, p. 106238, 2020

    work page 2020

  24. [24]

    A review of restoration experience from historical blackouts and a decision support framework for parallel restoration with a case study,

    H. De Silva, M. Nadarajah, M. Haque, M. Islam, and R. Madurai Elavarasan, “A review of restoration experience from historical blackouts and a decision support framework for parallel restoration with a case study,”Electric Power Systems Research, vol. 248, p. 111933, 2025

    work page 2025

  25. [25]

    Incorporating wind energy in power system restoration planning,

    A. Golshani, W. Sun, Q. Zhou, Q. P. Zheng, and Y . Hou, “Incorporating wind energy in power system restoration planning,”IEEE Transactions on Smart Grid, vol. 10, no. 1, pp. 16–28, 2019

    work page 2019

  26. [26]

    Governor rate-constrained opf for primary frequency control adequacy,

    H. Ch ´avez, R. Baldick, and S. Sharma, “Governor rate-constrained opf for primary frequency control adequacy,”IEEE Transactions on Power Systems, vol. 29, no. 3, pp. 1473–1480, 2014

    work page 2014

  27. [27]

    Analytical method to aggregate multi-machine sfr model with applications in power system dynamic studies,

    Q. Shi, F. Li, and H. Cui, “Analytical method to aggregate multi-machine sfr model with applications in power system dynamic studies,”IEEE Transactions on Power Systems, vol. 33, no. 6, pp. 6355–6367, 2018

    work page 2018

  28. [28]

    An analytical model for frequency nadir prediction following a major disturbance,

    L. Liu, W. Li, Y . Ba, J. Shen, C. Jin, and K. Wen, “An analytical model for frequency nadir prediction following a major disturbance,”IEEE Transactions on Power Systems, vol. 35, no. 4, pp. 2527–2536, 2020

    work page 2020

  29. [29]

    A method for incorporating frequency nadir limits in power system restoration planning,

    X. Zou, I. Farhat, and J. W. Simpson-Porco, “A method for incorporating frequency nadir limits in power system restoration planning,” inIEEE Electrical Power and Energy Conference (EPEC), 2025, pp. 254–259

    work page 2025

  30. [30]

    Graph decomposition for constructing blackstart restoration strategies in benchmark cases,

    A. B. Birchfield, “Graph decomposition for constructing blackstart restoration strategies in benchmark cases,”Electric Power Systems Research, vol. 212, p. 108402, 2022

    work page 2022

  31. [31]

    Bullo,Lectures on Network Systems

    F. Bullo,Lectures on Network Systems. CreateSpace Independent Publishing Platform, 2018

    work page 2018

  32. [32]

    DC power flow revisited,

    B. Stott, J. Jardim, and O. Alsac, “DC power flow revisited,”IEEE Transactions on Power Systems, vol. 24, no. 3, pp. 1290–1300, 2009

    work page 2009

  33. [33]

    A frequency-constrained MILP framework for power system restoration planning,

    X. Zou, “A frequency-constrained MILP framework for power system restoration planning,” Master of Applied Science Thesis, University of Toronto, Toronto, Canada, 2025

    work page 2025

  34. [34]

    Dynamic equivalents for average system frequency behavior following major disturbances,

    M. L. Chan, R. D. Dunlop, and F. Schweppe, “Dynamic equivalents for average system frequency behavior following major disturbances,” IEEE Transactions on Power Apparatus and Systems, vol. 91, no. 4, pp. 1637–1642, 1972

    work page 1972

  35. [35]

    J. H. Chow and J. J. Sanchez-Gasca,Power System Modeling, Compu- tation, and Control. John Wiley & Sons, 2020

    work page 2020

  36. [36]

    Neplan AG,Turbine-Governor Models: Standard Dynamic Turbine- Governor Systems in NEPLAN Power System Analysis Tool, 2015

    work page 2015

  37. [37]

    Y ALMIP: A toolbox for modeling and optimization in MATLAB,

    J. L ¨ofberg, “Y ALMIP: A toolbox for modeling and optimization in MATLAB,” inProceedings of the CACSD Conference, Taipei, Taiwan, 2004

    work page 2004

  38. [38]

    Gurobi Optimizer Reference Manual,

    Gurobi Optimization, LLC, “Gurobi Optimizer Reference Manual,”

    work page

  39. [39]

    Available: https://www.gurobi.com

    [Online]. Available: https://www.gurobi.com

    work page