Method and system for evaluating inertia of power system and storage medium

Abstract

A method and a system for evaluating inertia of a power system and a storage medium. The method includes: injecting a cosine active power disturbance into the power system by small-disturbance injection, and obtaining frequency response at a node where the disturbance is injected, where the active power disturbance can be an energy storage, wind power, or photovoltaic power; acquiring an evaluation framework of inertia and frequency regulation capability of the power system according to relative characteristics of a frequency response function; and constructing a mathematical relationship between the impedance and frequency response characteristics according to a relationship among active power disturbance, frequency fluctuation and impedance.

Claims

1. A method for evaluating inertia of a power system, comprising: (S1) injecting active power disturbance signals varying in frequency into the power system at a node in a power grid, wherein at time t, the active power disturbance signals are each expressed as:
P(t)=P.sub.0+P.sub.m cos(2πf.sub.rt+φ.sub.p); Wherein P.sub.0 is an output power of a non-synchronous power supply under normal operation; P.sub.m is a magnitude of a cosine term P.sub.m cos(2πf.sub.rt+φ.sub.p) of an active power disturbance; φ.sub.p is a phase angle of the cosine term P.sub.m cos(2πf.sub.rt+φ.sub.p) of the active power disturbance; and f.sub.r is a disturbance frequency of the active power disturbance; and acquiring a first frequency response function FR(s) of the power system at a point of common coupling (PCC) where the active power disturbance is injected according to Fourier transform, or equivalently characterizing a second frequency response function FR(s) of the power system according to impedance response; wherein the first frequency response function FR(s) is expressed as: FR(s)≡F[f.sub.r]/P[f.sub.r]; wherein s=j2πf.sub.r; F[f.sub.r]=0.5f.sub.m exp(jφ.sub.f); P[f.sub.r]=0.5P.sub.m exp(jφ.sub.p); F[f.sub.r] is a Fourier component corresponding to a frequency f(t) at the PCC under the disturbance frequency f.sub.r; P[f.sub.r] is a Fourier component corresponding to an active power disturbance P(t) at the PCC under the disturbance frequency f.sub.r; and f(t)=f.sub.0+f.sub.m cos(2πf.sub.rt+φ.sub.f), wherein f.sub.0 is a rated frequency of the power system under normal operation, f.sub.m is a magnitude of a frequency cosine term f.sub.m cos(2πf.sub.r+φ.sub.f) at the PCC, φ.sub.f is a phase angle of the frequency cosine term f.sub.m cos(2πf.sub.rt+φ.sub.f) at the PCC, and j is an imaginary unit in a complex field; or the second frequency response function FR(s) is expressed as: $\mathrm{FR}\left(s\right)\equiv \frac{s}{2\pi u_{\mathrm{sd}}}{T}_{\mathrm{PLL}}\left(s\right)G_i\left(s\right)\frac{1}{u_{\mathrm{sd}}}{Z}_{\mathrm{qd}}\left(s\right);R\left(s\right)\approx \frac{3s}{2\pi u_{\mathrm{sd}}^2}{Z}_{\mathrm{qd}}\left(s\right),f_r<f_{\mathrm{PLL\_BW}};$ wherein $T_{\mathrm{PLL}}\left(s\right)=\frac{U_1{H}_{\mathrm{PLL}}\left(s\right)/s}{1+U_1{H}_{\mathrm{PLL}}\left(s\right)/s};Z_{\mathrm{qd}}\left(s\right)$  is a coupling impedance between a q-axis and a d-axis of the power system in a dq reference frame at the PCC; G.sub.i(s) is the closed-loop gain of the current control; H.sub.PLL(s) is a transfer function of a phase-locked loop; s=j2π.sub.r; u.sub.sd is a d-axis component of a grid voltage at the PCC; f.sub.PLL_BW indicates phase-locked loop bandwidth; and U.sub.1 is the grid voltage at the PCC; and (S2) acquiring a Bode diagram of the first frequency response function FR(s) or the second frequency response function FR(s) to obtain a resonance point and a direct current (DC) gain K.sub.DC_gain; and determining an inertia level H.sub.sys of the power system according to the resonance point and the DC gain K.sub.DC_gain through the following formulas: $C_{\mathrm{eq}}^H.Math.\frac{\mathrm{df}}{\mathrm{dt}}=\Delta P\mathrm{and}\frac{2H_{\mathrm{sys}}{S}_{\mathrm{sys}}}{f_0}.\frac{\mathrm{df}}{\mathrm{dt}}=\Delta P;$ wherein C.sub.eq.sup.H is an equivalent capacitance characterized by an inertia of the power system, and is expressed as $C_{\mathrm{eq}}^H=\frac{1}{4{\pi }^2{f}_x^2{L}_{\mathrm{eq}}};L_{\mathrm{eq}}$  is an equivalent inductance of power lines of the power system; f.sub.x is a resonance frequency; ΔP is an unbalanced active power of the power system; and S.sub.sys is a rated capacity of the power system.

2. The method of claim 1, wherein in step (S1), the active power disturbance signals are injected into the power system by using a non-synchronous power supply, wherein the active power disturbance has a double-loop control structure; an outer loop of the double-loop control structure is configured to perform a power control, wherein a power control model is expressed as $i_{\mathrm{dref}}=\frac{2P_{\mathrm{ref}}}{3u_{\mathrm{sd}}},$  and i.sub.qref=0; an inner loop of the double-loop control structure is configured to perform a current control, wherein a current control model is expressed as i.sub.d(s)=G.sub.i(s)i.sub.dref the phase-locked loop is achieved by q-axis phase-locking, wherein a relation between an output angle of the phase-locked loop and an input q-axis voltage is expressed as $\Delta {\theta }_{\mathrm{PLL}}\left(s\right)=\frac{H_{\mathrm{PLL}}\left(s\right)/s}{1+U_1{H}_{\mathrm{PLL}}\left(s\right)/s}{u}_{sq};$ wherein P.sub.ref is a reference value of an input active power disturbance; i.sub.dref is a d-axis current reference value, and i.sub.qref is a q-axis current reference value; G.sub.i(s) is the closed-loop gain of the current control; P.sub.PLL(s) is the transfer function of the phase-locked loop; s=j2πf.sub.r; U.sub.1 is the grid voltage at the PCC; u.sub.sd is the d-axis component of the grid voltage at the PCC; u.sub.sq is a q-axis component of the grid voltage at the PCC; and Δθ.sub.PLL(s) is an angle variation of the grid voltage at the PCC.

3. The method of claim 1, wherein the step (S2) further comprises: evaluating frequency regulation capability of the power system according to the DC gain K.sub.DC_gain, wherein a quasi-stable-state frequency deviation Δf.sub.ss of the power system is expressed as Δf.sub.ss=K.sub.DC_gain*ΔP, wherein ΔP is the unbalanced active power of the power system.

4. A system for evaluating inertia of a power system, comprising: a computer device; wherein the computer device is configured or programmed to implement the method of claim 1.

5. A system for evaluating inertia of a power system, comprising: a computer device; wherein the computer device is configured or programmed to implement the method of claim 3.

6. A non-transitory computer-readable storage medium, comprising: a program runnable on a computer device; wherein the program is configured or programmed to implement the method of claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flowchart of a method for evaluating inertia of a power system according to an embodiment of this application;

(2) FIG. 2 schematically illustrates a framework for evaluating inertia and frequency regulation capability of a power system based on frequency response according to an embodiment of this application;

(3) FIG. 3 structurally shows a system for evaluating inertia of a power system in simulation according to an embodiment of this application;

(4) FIG. 4a is a schematic diagram of a synchronous power generation model according to an embodiment of this application;

(5) FIG. 4b is a schematic diagram of a voltage tracking converter model according to an embodiment of this application;

(6) FIG. 5 is a block diagram depicting the injection of an active power disturbance according to an embodiment of this application;

(7) FIG. 6 is a block diagram of the active power disturbance phase-locked loop (PLL) according to an embodiment of this application;

(8) FIG. 7 shows frequency response FR(s) at bus 1 where the active power disturbance is injected;

(9) FIG. 8 shows frequency response FR(s) at bus 2 where the active power disturbance is injected under different droop control coefficients of a converter;

(10) FIG. 9 illustrates frequency deviation simulation results under the active power disturbance of 0.05 pu; and

(11) FIG. 10 shows Z.sub.qd(s) response of the system at the bus 1 under the active power disturbance.

DETAILED DESCRIPTION OF EMBODIMENTS

(12) This application provides a method for evaluating inertia of a power system, which does not rely on large-disturbance events, and merely needs impedance response to equivalently characterize frequency response characteristics of the power system, so as to evaluate the inertia level of the power system and the frequency regulation capability of various kinds of power-supply units.

(13) The method is specifically performed as follows.

(14) 1) An active power disturbance is injected into the power system using a non-synchronous power supply, such as energy storage, and is expressed as:
P(t)=P.sub.0+P.sub.m cos(2πf.sub.rt+φ.sub.p)

(15) Frequency at a node where the active power disturbance is injected is expressed as: f(t)=t.sub.0+f.sub.m cos(2πf.sub.rt+φ.sub.f), where P.sub.0 is an output power of the non-synchronous power supply under normal operation, f.sub.0 is a rated frequency of the non-synchronous power supply in the power system under normal operation, P.sub.m is a magnitude of a cosine term P.sub.m cos(2πf.sub.rt+φ.sub.p) of an active power disturbance, φ.sub.p is a phase angle of the cosine term P.sub.m cos(2πf.sub.rt+φ.sub.p) of the active power disturbance, f.sub.in is a magnitude of a frequency cosine term f.sub.m cos(2πf.sub.rt+φ.sub.f), φ.sub.f is a phase angle of a frequency cosine term f.sub.m cos(2πf.sub.rt+φ.sub.f), and f.sub.r is a disturbance frequency of the active power disturbance.

(16) 2) A first frequency response function FR(s) is acquired according to Fourier transform, and the first frequency response function FR(s) is expressed as FR(s) F[f.sub.r]/P[f.sub.r]; where s=j2πf.sub.r; F[f.sub.r]=0.5f.sub.m exp(jφ.sub.f); F[f.sub.r] is a Fourier component corresponding to a frequency f(t) under the disturbance frequency f.sub.r, P[f.sub.r] is a Fourier component corresponding to the active power disturbance P(t) under the disturbance frequency f.sub.r. A Bode diagram of the first frequency response function FR(s) of the power system under the active power disturbance is acquired to obtain a resonance point and a DC gain K.sub.DC_gain. An inertia level of the power system is determined according to the resonance point (the lowest point of the magnitude response), and the frequency regulation capability of the power system is determined according to the DC gain K.sub.DC_gain.

(17) 3) Combining with a relationship among active power disturbance, frequency fluctuation and impedance, a non-intrusive frequency response transfer function equivalent model based on impedance is established, that is, a second frequency response function FR (s) of the power system is equivalently characterized by the impedance response, so as to evaluate the inertia and frequency regulation capability of the non-intrusive power system based on impedance according to the evaluation framework in step 2).

(18) In this embodiment, the inertia level of the power system and the frequency regulation capability of various kinds of power-supply units are evaluated through the equivalent impedance response of the node where the power supply is connected, which are obtained by measuring or equivalently calculating, without relying on the large disturbance events in the power system.

(19) In step 1), an active power disturbance is injected into the power system using a non-synchronous power supply such as energy storage, which is taken as an active power disturbance signal. The active power disturbance signal is configured to use a grid-following converter, which fails to simulate or control virtual inertia, and enables to flexibly adjust an output power. The non-synchronous power supply is configured to adopt a double-loop control structure. An outer loop of the double-loop control structure is configured to perform a constant-power control, where a power control model is expressed as

(20) $i_{\mathrm{dref}}=\frac{2P_{\mathrm{ref}}}{3u_{sd}},$
and i.sub.qref=0. An inner loop of the double-loop control structure is configured to perform a current control, where a current control model is expressed as i.sub.d(s)=G.sub.i(s)i.sub.dref. The phase-locked loop is achieved by q-axis phase-locking, where a relation between an output angle of the phase-locked loop and an input q-axis voltage is expressed as:

(21) $\Delta {\theta }_{\mathrm{PLL}}\left(s\right)=\frac{H_{\mathrm{PLL}}\left(s\right)/s}{1+U_1{H}_{\mathrm{PLL}}\left(s\right)/s}{u}_{sq},$
where P.sub.ref is a reference value of an input active power disturbance, i.sub.dref is a d-axis current reference value, and i.sub.qref is a q-axis current reference value, G.sub.i(s) is the closed-loop gain of the current control, and is expressed as

(22) $G_i\left(s\right)=\frac{\left(s*K_p^c+K_i^c\right)/\left(s*\left(R_c+s*L_c\right)\right)}{1+\left(s*K_p^c+K_i^c\right)/\left(s*\left(R_c+s*L_c\right)\right)},$
and H.sub.PLL(s) is the transfer function of the phase-locked loop, and is expressed as

(23) $H_{\mathrm{PLL}}\left(s\right)=K_p^{\mathrm{pll}}+\frac{K_i^{\mathrm{pll}}}{s},s=j2\pi f_r,$
U.sub.1 is grid voltage at the PCC, u.sub.sd is the d-axis component of the grid voltage at the PCC, u.sub.sq is a q-axis component of the grid voltage at the PCC, and Δθ.sub.PLL(s) is an angle variation of the grid voltage at the PCC, K.sub.p.sup.c is a proportion coefficient of the current-control loop, K.sub.i.sup.c is an integration coefficient of the current-control loop, K.sub.p.sup.pll is a proportion coefficient of the phase-locked loop, K.sub.i.sup.pll is an integration coefficient of the phase-locked loop, R.sub.c is a resistance at the PCC, and L.sub.c is an inductance at the PCC.

(24) The step 2) is specifically performed as follows.

(25) 4.1) An input active power disturbance signal and an output frequency signal corresponding to it are subjected to Fourier transformed to obtain a first frequency response function calculation model of the power system, which is expressed as FR(s)≡F[f.sub.r]/P[f.sub.r], where s=j2πf.sub.r, F[f.sub.r]=0.5f.sub.m exp(jφ.sub.f); and P[f.sub.r]=0.5P.sub.m exp(jφ.sub.p). F[f.sub.r] is a Fourier component corresponding to a frequency f(t) at the PCC under the disturbance frequency f.sub.r, P[f.sub.r] is a Fourier component corresponding to an active power disturbance P(t) at the PCC under the disturbance frequency f.sub.r. The magnitude and phase angle of the frequency response of the power system in a certain frequency range is obtained via frequency sweeping, so as to obtain a corresponding Bode diagram of the frequency response of the power system, and obtain a resonance point.

(26) 4.2) The inertia of the power system is evaluated according to the resonance point on the Bode diagram and the equivalent inductance of the power grid, and the frequency regulation capability of the power system is evaluated according to the frequency response DC gain K.sub.DC_gain. The resonance point can be understood as a resonance occurred between the inertia (capacitive) of the power system and an equivalent inductance (inductive) of the power line, such that C.sub.eq.sup.H is an equivalent capacitance characterized by the inertia of the power system, and is expressed as

(27) 0$C_{\mathrm{eq}}^H=\frac{1}{4{\pi }^2{f}_x^2{L}_{\mathrm{eq}}},$
where L.sub.eq is the equivalent inductance of the of power lines of the power system; and f.sub.x is a resonance frequency, such that the inertia level H.sub.sys of the power system is determined through the following formulas:

(28) $C_{\mathrm{eq}}^H.Math.\frac{\mathrm{df}}{\mathrm{dt}}=\Delta P\mathrm{and}\frac{2H_{\mathrm{sys}}{S}_{\mathrm{sys}}}{f_0}.Math.\frac{\mathrm{df}}{\mathrm{dt}}=\Delta P.$
Frequency regulation capability of the power system is indirectly demonstrated through K.sub.DC_gain (mHz/MW), where a quasi-stable-state frequency deviation Δf.sub.ss of the power system is expressed as Δf.sub.ss=K.sub.DC_gain*ΔP.

(29) The step 3) is specifically performed as follows.

(30) 5.1) By means of harmonic linearization method combined with the harmonic signal flowchart, a second frequency response function of the power system is expressed as follows:

(31) $\mathrm{FR}\left(s\right)=-\frac{\mathrm{sj}}{2\pi u_{\mathrm{sd}}}{T}_{\mathrm{PLL}}\left(s\right)G_i\left(s\right)\frac{2}{3u_{\mathrm{sd}}}.Math.\frac{1}{2}\left[Z_{\mathrm{pp}}\left(s\right)+Z_{\mathrm{pn}}\left(s\right)-Z_{\mathrm{np}}\left(s\right)-Z_{\mathrm{nn}}\left(s\right)\right],$
where

(32) $T_{\mathrm{PLL}}\left(s\right)=\frac{U_1{H}_{\mathrm{PLL}}\left(s\right)/s}{1+U_1{H}_{\mathrm{PLL}}\left(s\right)/s};T_{\mathrm{PLL}}\left(s\right)$
is a gain of the phase-locked closed-loop, Z.sub.pp(s) is a positive sequence impedance in the sequence domain of the power system at the node where the active power disturbance is injected, Z.sub.nn(s) is a negative sequence impedance in the sequence domain of the power system at the node where the active power disturbance is injected, Z.sub.pn(s) is a positive sequence versus negative sequence coupling impedance in the sequence domain of the power system at the node where the active power disturbance is injected, and Z.sub.np(s) is a negative sequence versus positive sequence coupling impedance in the sequence domain of the power system at the node where the active power disturbance is injected.

(33) 5.2) According to a conversion relationship between the sequence domain and the dq domain is expressed as

(34) $\left[\begin{array}{cc}{Z}_{\mathrm{dd}}\left(s\right)& Z_{\mathrm{dq}}\left(s\right)\\ Z_{\mathrm{qd}}\left(s\right)& Z_{\mathrm{qq}}\left(s\right)\end{array}\right]=T^{-1}\left[\begin{array}{cc}{Z}_{\mathrm{pp}}\left(s\right)& Z_{\mathrm{pn}}\left(s\right)\\ Z_{\mathrm{np}}\left(s\right)& Z_{\mathrm{nn}}\left(s\right)\end{array}\right]T,$
where Z.sub.dd(s) is a coupling impedance between the d-axis and the d-axis of the power system in the dq domain at the node where the active power disturbance is injected, Z.sub.dq(s) is a coupling impedance between the q-axis and the q-axis of the power system in the dq domain at the node where the active power disturbance is injected, Z.sub.qd(s) is a coupling impedance between the q-axis and the d-axis of the power system in the dq domain at the node where the active power disturbance is injected, Z.sub.qq(s) is a coupling impedance between the q-axis and the q-axis of the power system in the dq domain at the node where the active power disturbance is injected, such that the second frequency response function in step 5.1) is simplified into:

(35) $\mathrm{FR}\left(s\right)\equiv \frac{T_{\mathrm{PLL}}\left(s\right)G_i\left(s\right)Z_{\mathrm{qd}}\left(s\right)}{2\pi u_{\mathrm{sd}}^2},$
where a conversion matrix is expressed as

(36) $T=\left[\begin{array}{cc}1& j\\ 1& -j\end{array}\right].$
If it is further assumed that an input frequency of the active power disturbance is less than the bandwidth of the phase-locked loop, the second frequency response function FR (s) can be further expressed as:

(37) $\mathrm{FR}\left(s\right)\approx \frac{3s}{4\pi u_{\mathrm{sd}}^2}{Z}_{\mathrm{qd}}\left(s\right),f_r<f_{\mathrm{PLL\_BW}}.$
In engineering practice, the coupling impedance of the q-axis and the d-axis can be obtained through measurement or small signal modeling, which is supported in an article (Yunjie Gu, Yitong Li, YueZhu, Timothy C. Green. Impedance-Based Whole-System Modeling for a Composite Grid via Embedding of Frame Dynamics[J]. IEEE TRANSACTIONS ON POWER SYSTEMS, 2020, 1-10.)

(38) 5.3) As illustrated in the frequency response equivalent function FR(s) of the power system obtained in 5.2), the frequency response function of the power system under the injection of the disturbance is characterized using the coupling impedance between the q-axis and d-axis of the power system, so as to evaluate the contribution of various power generation units on the inertia and frequency regulation capability of the power system using the impedance response model.

(39) As shown in FIG. 1, the method provided herein uses non-synchronized power sources such as energy storage or wind power/photovoltaic power, and adopts appropriate control methods to inject cosine active power disturbance into the power system to obtain the second frequency response function FR(s) at the node where the disturbance is injected in the power system. According to the analysis block diagram of the frequency response shown in FIG. 2, the inertia and of the power system and the frequency regulation capability of the non-synchronous power supply are evaluated.

(40) This embodiment is specifically performed as follows.

(41) As shown in FIG. 3, a two-machine power system consisting of a traditional synchronous power generator, a voltage structural converter and corresponding loads is established in a MATrix LABoratory (Matlab)/Simulink software. The synchronous power generator is shown in FIG. 4(a), and the main parameters of the synchronous power generator is illustrated in Table 1. The converter model is shown in FIG. 4(b), and the main parameters of the converter model is illustrated in Table 2. In this embodiment, the active power disturbance is a grid-following converter based on PLL, whose block diagram is shown in FIG. 5. In addition, the phase-locked loop is designed to track the q-axis voltage, and obtain the frequency of the power system via a proportional integral controller. The block frame of controlling the active power disturbance phase-locked loop (PLL) is shown in FIG. 6.

(42) TABLE-US-00001 TABLE 1 Main parameters of synchronous power generator Parameter Value Inertia time constant (H) 10 Damping coefficient (D) 5 Time coefficient of speed controller (T.sub.g) 0.5 Droop coefficient of speed controller (R.sub.g) 0.02 Electrical impedance of power generator (Z.sub.St) 0.01 + j0.02

(43) TABLE-US-00002 TABLE 2 Main parameters of voltage structural converter Parameter Value Rated angular frequency of power system (ω.sub.0) 314.159 Output filter inductance (L.sub.f) 0.05/314.159 Output filter resistance (R.sub.f) 0.01 Output filter capacitance (C.sub.f) 0.02/314.159 Impedance of connection between converter and power grid 0.02 + j0.01 (Z.sub.C) Droop control coefficient under active power disturbance 0.02 (K.sub.Dp) Bandwidth under alternating voltage control 250 Bandwidth under alternating current control 500

(44) It is demonstrated in FIG. 5, if the active power disturbance, which is expressed as: f(t)=f.sub.0+f.sub.m cos(2πf.sub.rt+φ.sub.f), with the injection frequency f.sub.r is injected into the power system at bus 1, the injected active power disturbance generates a frequency response with an oscillation frequency f.sub.r, which is expressed as: f(t)=f.sub.0+f.sub.m cos(2πf.sub.rt+φ.sub.f), and then a Bode diagram of the frequency response transfer function FR (s) at the node where the active power disturbance is injected is obtained combined with Fourier transform, which is shown in FIG. 7.

(45) A verification of evaluating the inertia in the power system is shown as follows.

(46) It can be illustrated from FIG. 7 that the frequency response magnitude of the power system is the smallest at a frequency of 4.463 Hz, which is resulted from a resonance occurred between the capacitive of the inertia of the power system and equivalent inductance of the power line. According to the evaluation framework of the inertia of the power system in FIG. 2, it is assumed that only the influence of the inertia of the power system is considered, C.sub.eq.sup.H (Capacitive impedance characterized by the inertia of the power system) in the system provided herein is obtained as 0.789 by combining with the relevant parameters of the power system, such that the time constant of the inertia of the power system is obtained as 9.8625 s, and the error between the time constant of the inertia of the power system and the true value of the time constant of the inertia of the power system is only 1.375%. Thus, the inertia level of the power system evaluated by the method provided herein can accurately evaluate the inertia level of the power system.

(47) A verification of the evaluation of the frequency regulation capability in the power system is shown as follows.

(48) As shown in FIG. 8, when the active power disturbance is injected at bus 2, the frequency response characteristics of the power system under different droop gains of the grid-forming converter are illustrated. It can be illustrated from FIG. 8 that when the droop gain of the grid-forming converter K.sub.Dp is 0.02, the corresponding DC gain of the frequency response FR (s) of the power system is 56.975 dB, such that when the load of the power system suddenly increases by 0.05 pu, the quasi-stable-state frequency deviation of the power system is calculated as 35.295 mHz according to the DC gain K.sub.DC_gain. And when the droop gain of the grid-forming converter K.sub.Dp is 0.04, the corresponding DC gain of frequency response of the power system is 60.749 dB, such that when the load of the power system suddenly increases by 0.05 pu, the quasi-stable-state frequency deviation of the power system is calculated as 57.5 mHz according to the DC gain K.sub.DC_gain. Compared with the simulation results of the Simulink software shown in FIG. 9, when the droop gain of the grid-forming converter K.sub.Dp is 0.02, the quasi-stable-state frequency deviation of the power system is 37.12 mHz under time-domain simulation, and when the droop gain of the grid-forming converter K.sub.Dp is 0.04, the quasi-stable-state frequency deviation of the power system is 57.34 mHz under time-domain simulation. Therefore, the DC gain of the frequency response FR (s) is capable of accurately evaluating the frequency regulation capability of the system/equipment.

(49) An equivalent verification of the frequency response of the power system based on impedance is shown as follows.

(50) By comparing with FIG. 10 and FIG. 7, it can be illustrated that the curve in FIG. 10 and a curve in FIG. 7 are almost the same in the low frequency range, and have almost the same resonance frequency. The DC gain on the two curves is also approximately equal. It can be seen that the frequency response equivalent model of the power system based on impedance can well characterize the actual frequency response process of the power system, and thus the inertia and the frequency regulation capability of the power system can be evaluated by obtaining the impedance response under the disturbance injection. In engineering practice, the impedance response is easier to be obtained, such that the method provided herein is extremely recommended in practical engineering with strong practicability.