## Comparative Study of Data-driven Area Inertia Estimation Approaches on WECC Power Systems

Bendong Tan, Jiangkai Peng, Ningchao Gao, Junbo Zhao, Jin Tan

Abstract -With the increasing integration of inverter-based resources into the power grid, there has been a notable reduction in system inertia, potentially compromising frequency stability. To assess the suitability of existing area inertia estimation techniques for real-world power systems, this paper presents a rigorous comparative analysis of system identification, measurement reconstruction, and electromechanical oscillation-based area inertia estimation methodologies, specifically applied to the large-scale and multi-area WECC 240-bus power system. Comprehensive results show that the system identification-based approach exhibits superior robustness and accuracy relative to its counterparts.

power quantification and insufficient PMU placements. Furthermore, in the context of systems with a significant presence of IBRs, inertia might be unevenly distributed throughout the system [11]. This introduces spatial variances, nuances that system-level inertia estimation fails to capture. While the overall system inertia might seem satisfactory, specific regions may experience an inertia deficit, potentially giving rise to stability concerns.

Index Terms -Area inertia estimation, WECC 240-bus power system, system identification, measurement reconstruction, electromechanical oscillation.

## I. INTRODUCTION

The role of inertia in ensuring the frequency stability of power systems is pivotal [1]. Nonetheless, the substantial integration of inverter-based resources (IBRs), such as wind energy and solar, has led to a marked decline in system inertia. Consequently, obtaining an accurate understanding of inertia distribution becomes imperative for the secure operation of power systems [2].

Extensive research has been dedicated to inertia estimation at the system level, focusing on quantifying the cumulative inertia within an entire power system. This approach conceptualizes all generators as a singular equivalent entity, the dynamics of which are articulated through the aggregated swing equation, framed under the center of inertia (COI) coordinate system [1]. Regarding the system-level inertia estimation, both the total disturbed active power and the COI frequency are requisite [3]. System inertia can be estimated via the dynamics of the aggregated swing equation or the definition of inertia. The COI frequency can be depicted as the weighted average frequency derived from phasor measurement units (PMUs) positioned at various locations [4]-[6]. On the other hand, the total disturbed power can be obtained via recorded large disturbance [5] or the total power variation at each generator [7].

However, acquiring the total disturbed power poses challenges, primarily attributed to the complexities in disturbance

This work was authored in part by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. This material is based upon work supported by the U.S. Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE) under the Solar Energy Technologies Office Award Number 37772. B. Tan and J. Zhao are with University of Connecticut, Storrs, CT 06269. This work was based on Bendong Tan's summer internship at NREL. J. Peng, N. Gao and J. Tan are with National Renewable Energy Laboratory, Golden, CO 80401. (e-mail: jin.tan@nrel.gov)

Hence, if system operators are inclined towards understanding inertia distribution, area inertia estimation emerges as a more favorable approach. Analogous to system-level inertia estimation, area inertia determination relies on the total disturbed power within the area and the corresponding area COI frequency. However, the cumulative disturbed power within a specific area can be ascertained by aggregating the changes in tie-line power [10]. Approaches for area inertia estimation can be categorized into three distinct categories: (1) System identification-based methods; (2) Measurement reconstructionbased methods; (3) Electromechanical oscillation-based methods. Within the realm of system identification-based methodologies, both [9] and [11] employ a reduced-order model, i.e., the transfer function, to emulate the dynamic input-output relationship between the total disturbed area power and the area COI frequency. This results in a straightforward computation of area inertia, anchored on the inertia's fundamental definition. Transitioning to the measurement reconstructionbased techniques, [12] deploys polynomial expansion to reconstruct both power deviation and frequency. Concurrently, [13] harnesses dynamic mode decomposition (DMD), yielding an analytical framework for deducing area inertia. Delving into the electromechanical oscillation-based strategies, both [14] and [15] establish the relationship between the electromechanical oscillation mode and area inertia, facilitating the quantification of the latter.

The previously discussed inertia estimation methodologies, tested on benchmark systems with coherent generator dynamics, have yet to be comprehensively compared in the same benchmark system or real-world systems under practical testing scenarios where areas are not strictly defined by generator coherency. The goal of this paper is to understand the pros and cons of each method, i.e, on the WECC 240-bus system, a reduced model of the US western grid. Here, frequency stability is maintained through dispatch in balanced areas, which may not always correspond to generator coherency [16], [17].

## II. PROBLEM FORMULATION

In a specific regional power system, as shown in Fig. 1, the dynamics between the area's tie-line power and its center

Fig. 1. Diagram of the multi-area power system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-a6a0a0b3-93ec-4693-be16-f51717dd1303.jpg)

of inertia (COI) frequency can be effectively characterized by representing the system with an equivalent synchronous generator. This representation yields an aggregated swing equation given by:

$$\dot { \omega } _ { \text{COI}, i } = \frac { 1 } { 2 H _ { i } } \left ( P _ { m, i } - P _ { e, i } - D _ { i } \Delta \omega _ { \text{COI}, i } \right ) \quad ( 1 ) \quad \text{where} \\ \text{there} \, \dots \, \text{andnotac the COI} \, \text{fromcan} \, \text{for} \, \text{the} \, i \, \text{th} \, a \, \text{or}.$$

where ω COI ,i denotes the COI frequency for the i -th area:

$$\omega \ w c o l, i \text{ onwos uv } \times \omega \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ w p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p\ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u v \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \ u p \$$

where ω m and S m are the generator m's rotor speed and rated power respectively. P m,i and P e,i are representative of the tieline total mechanical and electrical power for the i -th area, respectively. Additionally, H i signifies the equivalent inertia and D i the damping factor pertinent to the i -th area. Since primary frequency control can be formulated as K i ∆ ω COI ,i , where K i is the primary frequency control coefficient, the perturbations effects associated with P m,i can be included into term D i ∆ ω COI ,i in (1). Therefore, (1) can be further reformulated in an incremental form:

where ∆ P e,i signifies the deviation in electrical power from its steady-state value. When P m,i undergoes dynamic changes, as observed during the initiation of primary frequency control, (3) yields a higher estimate for the equivalent damping factor.

$$\text{Immature in an element from} \,. \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \   \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \
 \Delta \dot { \omega } _ { \text{COL}, i } = \frac { 1 } { 2 H _ { i } } \left ( - \Delta P _ { e, i } - D _ { i } \Delta \omega _ { \text{COL}, i } \right ) \quad \quad \quad$$

To estimate H i , the primary challenge lies in accurately determining ∆˙ ω COI ,i . Relying solely on the finite difference method for frequency measurements can introduce numerical errors. Therefore, existing state-of-art methods try to estimate inertia without directly calculating ∆˙ ω COI ,i .

## III. AREA INERTIA ESTIMATION APPROACHES

In the following section, we benchmark cutting-edge methodologies from each category of area inertia estimation on a large-scale practical power system. Specifically, we adopt the techniques proposed in [11], [13], and [15] which represent the system identification-based, measurement reconstructionbased, and electromechanical oscillation-based methods, respectively. It is worth mentioning that enhancements have been incorporated into the techniques outlined in [11] and [13] to address and mitigate numerical challenges.

## A. System Identification-based Method

The Laplace transformation applied to both sides of (3) results in the relation ∆ ω COI ,i ( s ) ∆ P e,i ( s ) ≈ -1 2 H s i + D i . Acknowledging the influences of additional control components linked with the generator, the central premise of [11] is to employ a higherorder transfer function, denoted as G s i ( ) , to act as a surrogate model elucidating the relationship between ω COI ,i and ∆ P e,i . The form of G s i ( ) is given by:

$$\dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ G _ { i } ( s ) = \frac { 1 } { a _ { n } s ^ { n } + \cdots + a _ { 1 } s + a _ { 0 } }, \quad ( 4 ) \\ \cdot a \ \text{cimifies the coefficient of the nth order term } \text{Given}$$

where a n signifies the coefficient of the n th-order term. Given the presence of other control modules linked to generators, n = 2 is deemed sufficient to approximate the relationship between ω COI ,i and ∆ P e,i . The parameters intrinsic to G s i ( ) can be discerned through the N4SID (Numerical Algorithm for Subspace State Space System Identification) method, leveraging the input ∆ P e,i and output ∆ ω COI ,i [11].

Subjecting G s i ( ) to a step function (emulating a scenario where ∆ P e,i = 1 pu) yields the step response g step ( ) t . This step response can be approximated using an N p -order polynomial:

$$g _ { \text{step} } ( t ) = c _ { N _ { p } } t ^ { N _ { p } } + \dots + c _ { 1 } t + c _ { 0 } \text{ \quad \ \ } ( 5 )$$

where t is the temporal variable and c N p is the coefficient of the N p -order term of g step ( ) t . Given the small magnitude of D i , which can be approximated as negligible, we propose to estimate H i with (6):

$$\text{ited} _ { \text{tie} ^ { - } } \left ( \begin{array} { c } \text{cumm} \, \infty \, \infty. \\ \widetilde { H } _ { i } = - \, \frac { \Delta P _ { e, i } ( t ) } { 2 \dot { \omega } _ { \text{COI,i} } ( t ) } \right | _ { t = 0 } = - \, \frac { 1 } { 2 \dot { g } _ { \text{step} } ( t ) } \right | _ { t = 0 } = - \frac { 1 } { 2 c _ { 1 } } \end{array} \right ) ( 6 )$$

where ˜ H i is estimated equivalent inertia for the i -th area. In the swing equation, inertia values of synchronous generator range from 1.75 s to approximately 10 s, and the damping factor varies from the magnitude of 10 -2 to 10 -1 [2], [18]. Hence, the damping effect is small enough to be neglected. The exclusion of D i in expression (6) is expected to result in a slightly biased estimate of inertia. Generators with similar rotor angle response following a disturbance are called coherent generators. The coherency between generators, particularly over large electrical distances, may not be significant, potentially leading to substantial errors in ω COI as an aggregate frequency representation. However, since the system identification-based method solely utilizes dynamics at the initial stage ( t = 0 ), the coherency of area generators exerts minimal influence on its performance.

## B. Measurement Reconstruction-based Method

A typical approach for measurement reconstructionbased method is using DMD to reconstruct the measurements of X ∈ R 2 N a × M , where X = [∆ ω COI , 1 , ∆ ω COI , 2 , · · · , ∆ ω COI ,N a , ∆ P e, 1 ∆ P e, 2 · · · , ∆ P e,N a ] T [13]. N a is the number of areas in the power system, while M is the number of measurements.

$$s \text{ the number of measurements.} \\ A _ { c } \left [ \begin{matrix} \varphi _ { k } ^ { \Delta \omega _ { c o l } } \\ - - \\ \varphi _ { k } ^ { \Delta P _ { c } } \end{matrix} \right ] = \lambda _ { k } \left [ \begin{matrix} \varphi _ { k } ^ { \Delta \omega _ { c o l } } \\ - - \\ \varphi _ { k } ^ { \Delta P _ { e } } \end{matrix} \right ] \\ \text{re $A_\L=\ln(A_\L)\/\L\tcol:\L\t$ is the samnlin} \, \text{time interval} \colon A_\L.$$

where A c = ln ( A d ) / ∆ t ; ∆ t is the sampling time interval; A d is the state transition matrix that leads to X 2: M = A d X 1: M -1 ; X 2: M denotes the measurements from time instants 1 to M -1 ; φ ∆ P e k and φ ∆ P e k denote the the k -th DMD modes with respect to ∆ ω COI and ∆ P e ; λ k is the k -th eigenvalue. As a result, measurements can be reconstructed as:

$$\left [ \begin{array} { c } \Delta \omega _ { \text{COI} } \\ \Delta P _ { e } \end{array} \right ] = \sum _ { k = 1 } ^ { 2 N _ { a } } e ^ { \lambda _ { k } t } b _ { k } \left [ \begin{array} { c } \varphi _ { k } ^ { \Delta \omega _ { \text{COI} } } \\ \varphi _ { k } ^ { \Delta P _ { e } } \end{array} \right ] \quad \quad ( 8 )$$

where b k is the initial value coefficient corresponding to k -th eigenvalues. By applying (8) to (3), (9) can be derived:

$$\left [ \begin{matrix} \mathcal { R } \left ( \varphi ^ { \Delta \omega _ { \text color} } \lambda b \right ) & \mathcal { R } \left ( \varphi ^ { \Delta \omega _ { \text color} } b \right ) \\ \mathbb { m } \left ( \varphi ^ { \Delta \omega _ { \text color} } \lambda b \right ) & \mathbb { m } \left ( \varphi ^ { \Delta \omega _ { \text color} } b \right ) \right ] \left [ \begin{array} { c } 2 \widetilde { H } \\ D \end{array} \right ] = - \left [ \begin{array} { c } \mathcal { R } \left ( \varphi ^ { \Delta P _ { \text } } b \right ) \\ \mathbb { m } \left ( \varphi ^ { \Delta P _ { \text } } b \right ) \end{array} \right ] \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        }                                                                                                                                                                                                        \end{matrix} } \right ] \text{ } \text{ } \text{ } \text{                                                                                                                                                                                                     \text{                                                                                                                                                                                                    \text{  inerti                                                                                                                                                                                                       }                                                                                                                                                                                                        |                                                                                                                                                                                                        -                                                                                                                                                                                                        +                                                                                                                                                                                                        <                                                                                                                                                                                                        </                                                                                                                                                                                                        >                                                                                                                                                                                                       .                                                                                                                                                                                                        )                                                                                                                                                                                                        ]                                                                                                                                                                                                        {$$

In the work of [13], inertia is inferred through the resolution of (9). Nonetheless, our observations suggest potential inaccuracies in scenarios where the aforementioned equation admits multiple solutions. In an attempt to constrain and refine the solution space, we introduce additional constraints:

$$0 \leq D \leq 1 \quad \quad \ \ ( 1 0 ) \ | \, \frac { 1 } { 2 }$$

By jointly solving (9) and (10), we posit that the inertia estimation converges to a more pragmatic and precise value.

## C. Electromechanical Oscillation-based Method

The study by [15] establishes a connection between area inertia and electromechanical oscillations, facilitating the inference of area inertia directly from these oscillations. Taking into account the frequency bandwidth B of the electromechanical oscillation, the Fourier transform of (3) can be derived as:

$$2 H _ { i } \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta \dot { \omega } _ { \text{COI}, i } ( \gamma ) + D _ { i } \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta \omega _ { \text{COI}, i } ( \gamma ) + \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta P _ { e, i } ( \gamma ) = 0 \left | \int _ { \substack { \stackrel { \sim } { \stackrel { \sim } { \stackrel { \sim } { \sim } } } } } ^ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \SSigma } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \SSigma } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrelrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \longrightarrow \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim } { \sim } \int _ { \stackrel { \sim } { \sim }$$

(11)

where γ is the frequency in the oscillation band. According to (11), H i can be estimated via (12).

## D. Comparisons among Various Methods

From the comprehensive analysis of the methodologies presented, it is evident that none necessitate the computation of COI frequency derivatives, a process known to encounter numerical challenges, particularly in the presence of noise or outliers.

- (1) The pros of each method are listed as follows:
- · System Identification: This method offers ease of use due to its limited hyperparameters, which minimally impact its performance.
- · Measurement Reconstruction: Beyond inertia estimation, this method also furnishes damping estimation, thereby yielding additional dynamic information.
- · Electromechanical Oscillation: Similar to measurement reconstruction-based method, this method can also provide damping estimation.

(2) The cons of each method and corresponding measures are summarized as follows:

- · System Identification: Despite yielding a biased estimation due to the omission of the damping factor, the resultant error is deemed acceptable, considering the relatively minor magnitude of the damping factor.

]

the final inertia estimation. Consequently, noise processing and careful choice of the start point of the response measurements is essential.

- · Electromechanical Oscillation: As delineated in (12), inertia estimation accuracy is contingent upon the bandwidth B , underscoring its sensitivity to B . Thus, precise estimation of B is crucial.

Fig. 2. One-line diagram of WECC 240-bus power system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-340ef1cf-a91b-4603-beb7-123c0eb1bc06.jpg)

## IV. NUMERICAL RESULTS

To comprehensively compare the pros and cons of each method, extensive validations have been conducted on the WECC 240-bus power system, where a total of 13 regional power grids have been delineated according to balanced areas, as shown in Fig. 2. Every scenario was examined in the presence of a disturbance, characterized by a three-phase shortcircuit fault at Bus 1002, initiated at t = 1 s and cleared by t = 1 1 . s. The entirety of the simulation spanned 10 seconds.

## A. Parameter Sensitivity Analysis for Different Methods

- · Measurement Reconstruction: Given that b , the initial value coefficient associated with eigenvalues in measurement reconstruction, is determined by the first data point in response measurements, noise considerably influences

Upon reviewing various methodologies, it becomes evident that the polynomial order N p , the initial value b , and the oscillation bandwidth B serve as hyperparameters for system identification, measurement reconstruction, and the electromechanical oscillation-based method, respectively. For b , it is imperative to designate the specific instant at which at which

$$\widetilde { H } _ { i } = \frac { 1 } { 2 } \frac { \text{Re} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta \omega _ { \text{col}, i } ( \gamma ) \right ] \text{Im} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta P _ { e, i } ( \gamma ) \right ] - \text{Re} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta P _ { e, i } ( \gamma ) \right ] \text{Im} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta \omega _ { \text{col}, i } ( \gamma ) \right ] } { \text{Re} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } j \gamma \Delta \omega _ { \text{col}, i } ( \gamma ) \right ] \text{Im} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } \Delta \omega _ { \text{col}, i } ( \gamma ) \right ] \text{Im} \left [ \sum _ { \gamma = 0 } ^ { \gamma = B } j \gamma \Delta \omega _ { \text{col}, i } \right ] } \quad ( \text{12)}$$

Fig. 3. Maximum area inertia estimation error under different hyperparameters for various methods. (a) Measurement reconstruction-based method; (b) System identification-based method; (c) Electromechanical oscillation-based method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-53c2e184-d645-418a-9f02-f1cbc0c21d53.jpg)

Fig. 4. Area inertia estimation error for various methods. (a) System identification-based method; (b) Measurement reconstruction-based method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-5453422a-d3e7-469b-b273-7158502c66ba.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-cdcd540d-7b4d-4b56-bbb7-4411a4871a8a.jpg)

the starting point of the measurement response is selected. To evaluate the efficacy of each method, we introduce two indices: the maximum estimation error (MEE) and the estimation error (EE):

$$\text{(EE):} \\ \text{MEE} = \max _ { 1 \leq i \leq N _ { a } } \left ( \left | \frac { \widetilde { H } _ { i } - H _ { i } } { H _ { i } } \right | \times 1 0 0 \% \right ), \text{EE} = \left | \frac { \widetilde { H } _ { i } - H _ { i } } { H _ { i } } \left | \times 1 0 0 \% \text{oisy measure} \\ \text{From Eio $2$ it is availat that on increased $N$ another $\alpha$ on each $\alpha$ and $\colon$}$$

From Fig. 3, it is evident that an increased N p enhances the fitting quality for g step ( ) t , thereby improving the accuracy of area inertia estimation. Specifically, N p = 4 emerges as an optimal choice, yielding the minimal MEE for the system identification-based method. For the measurement reconstruction-based method, the initial values b are chosen at several time steps post-disturbance clearing. As depicted in Fig. 3(a), selecting b shortly after disturbance clearing results in a MEE exceeding 400%. This can be attributed to the fact that while DMD offers a linear approximation of the dynamic system, measurements taken shortly post-disturbance exhibit pronounced nonlinear characteristics. Consequently, b is determined using the measurement from the 14-th time step post-disturbance to minimize the estimation error. Regarding the electromechanical oscillation-based method, its elevated MEE, reaching up to 100%, stems from the formulation of (12) through the continuous domain's Fourier transform. This necessitates approximation via the discrete Fourier transform, given that measurements are discretely acquired [15]. Owing to the substantial estimation error associated with the electromechanical oscillation-based method, it is excluded from subsequent comparative analyses.

to follow a zero-mean Gaussian distribution with a standard deviation of σ = 4 × 10 -3 [19]. Notably, this is greater than the realistic noise encountered in measurement data, which has an upper limit of 1% total vector error (TVE). To mimic the data preprocessing typically encountered in practical operations, noisy measurements undergo smoothing using a low-pass filter with a cut-off frequency of 5Hz [11]. Fig. 4 illustrates the inertia estimation errors for each area as determined by both the system identification and measurement reconstructionbased methods. In both noise-free and noisy conditions, the system identification-based method consistently yields estimation errors below 5% for all areas, with the exception of AZ. This underscores the method's robustness. In contrast, the measurement reconstruction-based method exhibits area inertia estimation errors exceeding 100%. Such discrepancies might arise from noise-induced alterations in the initial values b , which in turn precipitate significant estimation errors.

## C. Estimation Results under Different Operation Conditions

## B. Area Inertia Estimation Results Under Noise

In this section, the performance of system identification and measurement reconstruction-based methods are validated under measurement noise. The measurement noise is postulated

This study extends the validation of both the system identification and measurement reconstruction-based methods across diverse operational conditions. The two-week hourly operational data, spanning January 1-14, 2019, was obtained using NREL's Multi-timescale Integrated Dynamic and Scheduling Model (MIDAS) [20]. For each operational condition, a consistent disturbance and measurement noise, as previously described, is introduced. Owing to page constraints, we present the estimation results exclusively for states CO, OR, WA, CA1, and CA2, as these exhibit notable inertia fluctuations over time, as shown in Figs. 5 and 6. From the data, it is evident that the total system load undergoes periodic fluctuations, leading to the disconnection or connection of generators. As a result, the total inertia in each area exhibits time-varying characteristics. The system identification-based method showcases consistent estimation performance, with the maximum

Fig. 5. Inertia estimation results of system identification-based method for each area from January 1st to January 14th, 2019.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-dbdf4473-6fac-48be-94d8-26eb54b661c3.jpg)

Fig. 6. Inertia estimation results of measurement reconstruction-based method for each area from January 1st to January 14th, 2019.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00511v1-figure-6618bb0b-2962-43fc-a484-e505c90906c1.jpg)

estimation error remaining below 15%, as depicted in Fig. 5. While its estimation outcomes tend to be lower than the actual values, indicating a biased approach, its performance is nonetheless reliable. This further illustrates its suitability for regional power systems where generator coherency might not necessarily be maintained. Conversely, due to the influence of operational conditions on the initial value b , the measurement reconstruction-based method registers errors exceeding 100% for states OR, WA, CA1, and CA2, even though its estimations for CO appear accurate, as shown in Fig. 6. This underscores the limited robustness of the measurement reconstructionbased method in estimating area inertia. In summary, the system identification-based method emerges as the more favorable approach for practical power system area inertia estimation.

## V. CONCLUSIONS

In this paper, we undertook an exhaustive comparative analysis to assess the efficacy of system identification, measurement reconstruction, and electromechanical oscillation-based methods for area inertia estimation in a large-scale practical power system, specifically the WECC 240-bus power system. Numerical evaluations under conditions of noise and timevarying operations underscore that the system identificationbased method offers superior robustness and accuracy in area inertia estimation, despite its slight estimation bias, compared to the other methods.

## REFERENCES

- [1] P. Kundur, Power System Stability and Control . McGraw-Hill, 1994.
- [2] B. Tan, et al. , 'Power system inertia estimation: Review of methods and the impacts of converter-interfaced generations,' Int. J. Elect. Power Energy Syst. , vol. 134, p. 107362, 2022.
- [3] P. Wall and V. Terzija, 'Simultaneous estimation of the time of disturbance and inertia in power systems,' IEEE Trans. Power Syst. , vol. 29, no. 4, pp. 2018-2031, 2014.
- [4] J. Schiffer, et al. , 'Online estimation of power system inertia using dynamic regressor extension and mixing,' IEEE Trans. Power Syst. , vol. 34, no. 6, pp. 4993-5001, 2019.
- [5] B. Wang, et al. , 'Power system inertia estimation method based on maximum frequency deviation,' IET Renew. Power Gener. , vol. 16, no. 3, pp. 622-633, 2022.
- [6] C. Phurailatpam, et al. , 'Measurement based estimation of inertia in ac microgrids,' IEEE Trans. Sustain. Energy , pp. 1-1, 2019.
- [7] P. Makolo, et al. , 'Online inertia estimation for power systems with high penetration of res using recursive parameters estimation,' IET Renew. Power Gener. , vol. 15, no. 12, pp. 2571-2585, 2021.
- [8] F. Zeng, et al. , 'Online identification of inertia distribution in normal operating power system,' IEEE Trans. Power Syst. , vol. 35, no. 4, pp. 3301-3304, 2020.
- [9] K. Tuttelberg, et al. , 'Estimation of power system inertia from ambient wide area measurements,' IEEE Trans. Power Syst. , vol. 33, no. 6, pp. 7249-7257, 2018.
- [10] D. Li, et al. , 'Area inertia estimation of power system containing wind power considering dispersion of frequency response based on measured area frequency,' IET Gener. Transm. Distrib. , vol. 16, no. 22, pp. 4640-4651, 2022.
- [11] F. Zeng, et al. , 'Online estimation of power system inertia constant under normal operating conditions,' IEEE Access , vol. 8, pp. 101426-101436, 2020.
- [12] C. Phurailatpam, et al. , 'Measurement-based estimation of inertia in AC microgrids,' IEEE Trans. Sustain. Energy , vol. 11, no. 3, pp. 1975-1984, July 2020.
- [13] D. Yang, et al. , 'Data-driven estimation of inertia for multiarea interconnected power systems using dynamic mode decomposition,' IEEE Trans. Ind. Inform. , vol. 17, no. 4, pp. 2686-2695, April 2021.
- [14] G. Cai, et al. , 'Inertia estimation based on observed electromechanicaloscillation response for power systems,' IEEE Trans. Power Syst. , vol. 34, no. 6, pp. 4291-4299, Nov. 2019..
- [15] B. Wang, et al. , 'An improved electromechanical oscillation-based inertia estimation method,' IEEE Trans. Power Syst. , vol. 37, no. 3, pp. 2479-2482, May 2022.
- [16] H. Yuan, et al. , 'Developing a reduced 240- busWECC dynamic model for frequency response study of high renewable integration,' in Proc. IEEE PES Transmiss. Distrib. Conf. Expo. , 2020, pp. 1-5.
- [17] National Renewable Energy Laboratory, 'Test case repository for high renewable study,' 2020. [Online]. Available: https://www.nrel.gov/grid/test-case-repository.html.
- [18] R. K. Panda, et al. , 'Online estimation of system inertia in a power network utilizing synchrophasor Measurements,' IEEE Trans. Power Syst. , vol. 35, no. 4, pp. 3122-3132, July 2020.
- [19] R. Singh, et al. , 'Choice of estimator for distribution system state estimation'. IET Gener. Transm. Distrib. , vol. 3, no. 7, pp. 666, 2009.
- [20] J. Tan, et al. , 'Multi-timescale integrated dynamics and scheduling for solar (MIDAS-solar) (final technical report),' United States, 2023. [Online]. Available: https://doi.org/10.2172/1972321.
Bendong Tan Jiangkai Peng Ningchao Gao Junbo Zhao Jin Tan
Comparative Study of Data-driven Area Inertia Estimation Approaches on WECC Power Systems
With the increasing integration of inverter-based resources into the power
grid, there has been a notable reduction in system inertia, potentially
compromising frequency stability. To assess the suitability of existing area
inertia estimation techniques for real-world power systems, this paper presents
a rigorous comparative analysis of system identification, measurement
reconstruction, and electromechanical oscillation-based area inertia estimation
methodologies, specifically applied to the large-scale and multi-area WECC
240-bus power system. Comprehensive results show that the system
identification-based approach exhibits superior robustness and accuracy
relative to its counterparts.