pith. sign in

arxiv: 2604.15686 · v1 · submitted 2026-04-17 · 📡 eess.SY · cs.SY

DAE-Aware Bayesian Inference for Joint Generator-Network Parameter Estimation

Pith reviewed 2026-05-10 08:56 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords Bayesian inferenceparameter estimationdifferential-algebraic equationspower systemsgenerator parametersnetwork parametersIEEE test systems
0
0 comments X

The pith

A joint Bayesian framework estimates generator and network parameters directly from DAE power system models.

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

The paper develops a Bayesian inference method that treats generator physics and network power-flow equations together in their native differential-algebraic form rather than approximating them as ordinary differential equations. This joint approach recovers values for generator inertia and damping along with branch resistances and reactances while producing well-calibrated uncertainty estimates. Demonstrations on the IEEE 9-bus system recover the parameters accurately, and a 39-bus case indicates the method continues to work on a larger joint-estimation task without requiring restrictive prior distributions.

Core claim

By embedding the known DAE structure into a physics-aware statistical model and using efficient posterior sampling, the framework performs simultaneous Bayesian calibration of generator and network parameters, yielding accurate point estimates and well-behaved posterior uncertainty on standard IEEE test systems without overly conservative priors.

What carries the argument

Joint Bayesian inference procedure that directly incorporates the differential-algebraic equation (DAE) model of coupled generator dynamics and network power flow to enable simultaneous parameter estimation.

If this is right

  • Generator inertia, damping, and network resistances and reactances can be estimated jointly from the same data set.
  • Posterior distributions remain well-behaved even when the number of unknown parameters grows with system size.
  • The method works on the IEEE 9-bus system and shows evidence of remaining practical on the IEEE 39-bus system.
  • Accurate recovery occurs without the need for overly restrictive prior distributions on the parameters.

Where Pith is reading between the lines

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

  • If measurement quality is high enough, the same framework could support online monitoring of parameter drift in operating grids.
  • The approach may extend to other engineering domains whose models are naturally expressed as differential-algebraic equations, such as chemical process networks or mechanical multibody systems.
  • Future work could test whether the sampling efficiency scales to systems with hundreds of buses or with partial observability.

Load-bearing premise

The DAE model structure is known and accurately represents the true system dynamics, while the available measurements supply enough information to identify both generator and network parameters together.

What would settle it

Running the procedure on the IEEE 9-bus system with known true parameter values and checking whether the recovered posterior means deviate substantially from those values or the uncertainty intervals fail to contain the true parameters.

Figures

Figures reproduced from arXiv: 2604.15686 by Abdallah Alalem Albustami, Ahmad F. Taha, Sankaran Mahadevan.

Figure 1
Figure 1. Figure 1: Overview of the DAE-aware Bayesian framework. PMU measure [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Normalized co-identifiability matrix IXY , showing the strength of cross-block coupling in θ = [M⊤, D⊤, r⊤, x⊤]⊤. and Stage 0 enforces power-flow feasibility. Initialization: a stagewise procedure first fits (M, D) to frequency residuals, then (r, x) to all channels, followed by a short joint polish; the resulting point initializes the chain and the local Gauss– Newton curvature used in (25). C. Co-Identif… view at source ↗
Figure 3
Figure 3. Figure 3: Prior-to-posterior update for representative parameters under the joint Bayesian model. Posterior marginals contract around the true values while [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Joint and decoupled posterior marginals for selected entries of [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
read the original abstract

This paper addresses the classic problem of parameter estimation (PE) in multimachine power system models. Such models are typically described by a set of nonlinear differential-algebraic equations (DAE), where generator physics and network power flow equations are coupled. DAE models are well established in classic power system textbooks, but parameter identification and estimation of generator inertia and damping together with network branch resistances and reactances for these models remain relatively underexplored. In contrast to prior approaches that rely on ODE approximations, this paper develops a joint Bayesian inference framework to perform PE of generator and network parameters while exploiting grid DAE models. It further combines physics-aware statistical modeling with computationally efficient posterior sampling to make joint Bayesian calibration practical. Results on the IEEE 9-bus system show accurate parameter recovery with well-behaved posterior uncertainty, while a short 39-bus study provides evidence that the framework remains effective on a materially larger joint-estimation problem. These results are obtained without requiring overly conservative priors.

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 develops a joint Bayesian inference framework for parameter estimation (PE) of generator inertia/damping and network resistances/reactances in multimachine power systems. It exploits the full nonlinear DAE structure (generator dynamics coupled to network power-flow equations) rather than ODE approximations, combines physics-aware statistical modeling with efficient posterior sampling, and reports accurate recovery with well-behaved posteriors on the IEEE 9-bus system plus a brief 39-bus demonstration, all without requiring overly conservative priors.

Significance. If the central claims hold, the work is significant for power-system model calibration: it provides a practical route to joint generator-network PE that respects the established DAE physics, yields uncertainty quantification, and scales at least to 39-bus cases. The explicit use of standard Bayesian sampling on a physics-based DAE model, together with the reported recovery accuracy on benchmark systems, supplies a reproducible template that can be extended to other DAE-governed networks.

major comments (2)
  1. [§3] §3 (Methodology): the claim that the framework is 'computationally efficient' for joint estimation rests on the integration of the DAE solver inside the posterior sampler; the manuscript must specify the numerical integration scheme, step-size control, and any approximations used when the algebraic variables are solved at each likelihood evaluation, because these choices directly determine both accuracy and wall-clock cost for the 39-bus case.
  2. [§4.2] §4.2 (9-bus results): the reported 'accurate parameter recovery' is quantified only at summary level; the manuscript should include the full posterior marginals or credible-interval coverage against the known ground-truth values for all estimated parameters (inertia, damping, R, X) so that readers can judge whether the posteriors are truly well-behaved or merely consistent with the chosen priors.
minor comments (2)
  1. Notation: the distinction between generator states, algebraic network variables, and the measurement vector should be made explicit in the first appearance of the DAE and likelihood equations to avoid ambiguity when the same symbols appear in both the physics model and the statistical model.
  2. The 39-bus study is described as 'short'; adding a brief table of wall-clock times or effective sample sizes per parameter would strengthen the scalability claim without lengthening the paper.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of our work and the recommendation for minor revision. The comments are constructive and help strengthen the reproducibility and transparency of the manuscript. We address each major comment below and will make the corresponding revisions.

read point-by-point responses
  1. Referee: [§3] §3 (Methodology): the claim that the framework is 'computationally efficient' for joint estimation rests on the integration of the DAE solver inside the posterior sampler; the manuscript must specify the numerical integration scheme, step-size control, and any approximations used when the algebraic variables are solved at each likelihood evaluation, because these choices directly determine both accuracy and wall-clock cost for the 39-bus case.

    Authors: We agree that additional detail on the numerical implementation is necessary for reproducibility and to substantiate the efficiency claim. In the revised manuscript we will expand §3 with a new subsection that specifies the DAE integration method (a variable-step-size implicit Runge-Kutta scheme with embedded error control), the absolute and relative tolerances employed, the Newton solver settings used for the algebraic power-flow equations at each likelihood evaluation, and any simplifications (e.g., sparse Jacobian handling). We will also report average wall-clock time per likelihood evaluation on the 39-bus system to quantify the computational cost. revision: yes

  2. Referee: [§4.2] §4.2 (9-bus results): the reported 'accurate parameter recovery' is quantified only at summary level; the manuscript should include the full posterior marginals or credible-interval coverage against the known ground-truth values for all estimated parameters (inertia, damping, R, X) so that readers can judge whether the posteriors are truly well-behaved or merely consistent with the chosen priors.

    Authors: We concur that summary statistics alone are insufficient for readers to fully evaluate posterior behavior. In the revised §4.2 we will add figures displaying the complete marginal posterior densities for every estimated parameter (generator inertias, dampings, and all network R and X values) on the 9-bus system. Each plot will include the 95 % credible interval and an explicit marker for the ground-truth value, allowing direct assessment of coverage and calibration independent of the prior. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper presents a joint Bayesian inference framework that applies standard posterior sampling to a known DAE model structure for generator and network parameter estimation. The likelihood is constructed from the DAE dynamics and available measurements, with priors chosen to avoid overly conservative assumptions; posterior sampling then yields parameter estimates and uncertainty. No load-bearing step reduces by the paper's own equations to a fitted input renamed as a prediction, nor does any central claim rest on a self-citation chain that itself lacks independent verification. The IEEE 9-bus and 39-bus results constitute empirical validation outside the inference procedure itself, rendering the derivation self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The framework rests on the standard assumption that the multimachine power system is accurately described by the established nonlinear DAE structure from power system theory, with Bayesian priors chosen to be non-overly-conservative but not further specified.

free parameters (1)
  • Bayesian prior distributions
    Priors on generator and network parameters are used in the inference but described only as not overly conservative; specific forms or hyperparameters not detailed in abstract.
axioms (1)
  • domain assumption The power system dynamics are correctly captured by the given set of nonlinear differential-algebraic equations coupling generators and network power flow
    Invoked as the basis for exploiting grid DAE models instead of ODE approximations.

pith-pipeline@v0.9.0 · 5478 in / 1367 out tokens · 50464 ms · 2026-05-10T08:56:54.271388+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Parameter estima- tion of multiterminal transmission lines using joint pmu and scada data,

    S. S. Mousavi-Seyedi, F. Aminifar, and S. Afsharnia, “Parameter estima- tion of multiterminal transmission lines using joint pmu and scada data,” IEEE Transactions on Power Delivery, vol. 30, no. 3, pp. 1077–1085, 2014

  2. [2]

    Estimation of transmission line parameters by using two least-squares methods,

    F. P. Albuquerque, E. C. M. Costa, L. H. Liboni, R. F. R. Pereira, and M. C. de Oliveira, “Estimation of transmission line parameters by using two least-squares methods,”IET Generation, Transmission & Distribution, vol. 15, no. 3, pp. 568–575, 2021

  3. [3]

    Extended kalman filtering based real-time dynamic state and parameter estimation using pmu data,

    L. Fan and Y . Wehbe, “Extended kalman filtering based real-time dynamic state and parameter estimation using pmu data,”Electric Power Systems Research, vol. 103, pp. 168–177, 2013

  4. [4]

    Unscented kalman filter for power system dynamic state estimation,

    G. Valverde and V . Terzija, “Unscented kalman filter for power system dynamic state estimation,”IET generation, transmission & distribution, vol. 5, no. 1, pp. 29–37, 2011

  5. [5]

    Armax-based transfer function model identification using wide-area measurement for adaptive and coordinated damping control,

    H. Liu, L. Zhu, Z. Pan, F. Bai, Y . Liu, Y . Liu, M. Patel, E. Farantatos, and N. Bhatt, “Armax-based transfer function model identification using wide-area measurement for adaptive and coordinated damping control,” IEEE Transactions on Smart Grid, vol. 8, no. 3, pp. 1105–1115, 2015

  6. [6]

    Synchronous machine parameter identification using particle swarm op- timization,

    G. I. Hutchison, B. Zahawi, K. Harmer, B. Stedall, and D. Giaouris, “Synchronous machine parameter identification using particle swarm op- timization,” in5th IET International Conference on Power Electronics, Machines and Drives (PEMD 2010). IET, 2010, pp. 1–4

  7. [7]

    Non-invasive identification of turbo-generator parameters from actual transient network data,

    G. Hutchison, B. Zahawi, K. Harmer, S. Gadoue, and D. Giaouris, “Non-invasive identification of turbo-generator parameters from actual transient network data,”IET Generation, Transmission & Distribution, vol. 9, no. 11, pp. 1129–1136, 2015

  8. [8]

    A bayesian approach for parameter estimation with uncertainty for dynamic power systems,

    N. Petra, C. G. Petra, Z. Zhang, E. M. Constantinescu, and M. Anitescu, “A bayesian approach for parameter estimation with uncertainty for dynamic power systems,”IEEE Transactions on Power Systems, vol. 32, no. 4, pp. 2735–2743, 2016

  9. [9]

    Bayesian inertia estimation via parallel mcmc hammer in power systems,

    W. Zhong, C. Li, M. Chu, Y . Che, S. Zhou, Z. Wu, and K. Liu, “Bayesian inertia estimation via parallel mcmc hammer in power systems,”Ener- gies, vol. 18, no. 15, p. 3905, 2025

  10. [10]

    Bayesian inference of parameters in power system dynamic models using trajectory sensitivities,

    R. Nagi, X. Huan, and Y . C. Chen, “Bayesian inference of parameters in power system dynamic models using trajectory sensitivities,”IEEE Transactions on Power Systems, vol. 37, no. 2, pp. 1253–1263, 2021

  11. [11]

    Continuous estimation of power system inertia using convolutional neural networks,

    D. Linaro, F. Bizzarri, D. Del Giudice, C. Pisani, G. M. Giannuzzi, S. Grillo, and A. M. Brambilla, “Continuous estimation of power system inertia using convolutional neural networks,”Nature Communications, vol. 14, no. 1, p. 4440, 2023

  12. [12]

    Beyond gcn: Accurate real- time inertia estimation with multi-head spatial-temporal graph attention network coupled with lstm,

    F. Albeladi, K. Basulaiman, and M. Barati, “Beyond gcn: Accurate real- time inertia estimation with multi-head spatial-temporal graph attention network coupled with lstm,”IEEE Transactions on Industry Applica- tions, 2025

  13. [13]

    Machine learning assisted inertia estimation using ambient measurements,

    M. Tuo and X. Li, “Machine learning assisted inertia estimation using ambient measurements,”IEEE Transactions on Industry Applications, vol. 59, no. 4, pp. 4893–4903, 2023

  14. [14]

    Application of physics- informed machine learning techniques for power grid parameter estima- tion,

    S. Lakshminarayana, S. Sthapit, and C. Maple, “Application of physics- informed machine learning techniques for power grid parameter estima- tion,”Sustainability, vol. 14, no. 4, p. 2051, 2022

  15. [15]

    Power system dynamics and stability,

    P. S. Kundur, “Power system dynamics and stability,”Power System Stability and Control, Third Edition, pp. II–1–II–3, 2017

  16. [16]

    Dynamic state estimation of nonlinear differential algebraic equation models of power networks,

    M. Nadeem, S. A. Nugroho, and A. F. Taha, “Dynamic state estimation of nonlinear differential algebraic equation models of power networks,” IEEE Transactions on Power Systems, vol. 38, no. 3, pp. 2539–2552, 2022

  17. [17]

    A. R. Bergen,Power systems analysis. Pearson Education India, 2009

  18. [18]

    Identifiability in linear models,

    J. Richmond, “Identifiability in linear models,”Econometrica: Journal of the Econometric Society, pp. 731–736, 1974

  19. [19]

    A new method for parameter identifiability of dae systems,

    L. Chen, X. Chen, J. Qian, and Z. Yao, “A new method for parameter identifiability of dae systems,” in2013 10th IEEE International Confer- ence on Control and Automation (ICCA). IEEE, 2013, pp. 280–283

  20. [20]

    Extended space method for parameter identifiability of dae systems,

    ——, “Extended space method for parameter identifiability of dae systems,”ISA transactions, vol. 53, no. 5, pp. 1476–1480, 2014

  21. [21]

    An analysis of structural and practical identifiability applied to a transient generator model,

    T. G. Landgraf, E. P. Cari, and L. F. Alberto, “An analysis of structural and practical identifiability applied to a transient generator model,” Electric Power Systems Research, vol. 206, p. 107817, 2022

  22. [22]

    Equation of state calculations by fast computing machines,

    N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, “Equation of state calculations by fast computing machines,” The journal of chemical physics, vol. 21, no. 6, pp. 1087–1092, 1953

  23. [23]

    Monte carlo sampling methods using markov chains and their applications,

    W. K. Hastings, “Monte carlo sampling methods using markov chains and their applications,” 1970

  24. [24]

    Markov chain monte carlo using an approximation,

    J. A. Christen and C. Fox, “Markov chain monte carlo using an approximation,”Journal of Computational and Graphical statistics, vol. 14, no. 4, pp. 795–810, 2005