Frequency modulation equipment capability evaluation method and device based on energy visual angle

Abstract

The invention discloses a frequency modulation equipment capability evaluation method and device based on an energy visual angle. The method comprises the following steps: acquiring equipment parameters, network parameters and preset fault information of a system, and initializing equivalent speed regulator model parameters by aggregating the equipment parameters; based on the equivalent speed regulator model parameters, the network parameters and preset fault information, establishing a vibration mechanical model of the power system; solving the vibration mechanical model according to a modal analysis method, and obtaining a common-mode component and an oscillation component of the frequency variation of the power system; determining the frequency variation of each unit in the power system according to the common-mode component and the oscillation component; according to the frequency variation and the equivalent speed regulator model parameters, calculating the total unit frequency modulation energy accumulated at the maximum frequency deviation moment of each unit; and calculating the frequency modulation capability contribution degree of each unit based on the total frequency modulation energy of each unit, and outputting a frequency modulation capability evaluation conclusion of each unit according to the frequency modulation capability contribution degree of each unit.

Description

Technical Field

[0001] This invention relates to the field of power system operation technology, and more specifically, to a method and apparatus for evaluating the capabilities of frequency regulation equipment based on an energy perspective. Background Technology

[0002] Currently, the mainstream indicators for evaluating the frequency regulation capability of synchronous generators, new energy generator sets, and other frequency regulation equipment are reserve capacity and frequency regulation response speed. However, these evaluation methods neglect the reflection of frequency space differences and are difficult to comprehensively consider the coupling effects of dynamic response and various limiting mechanisms. As the response modes of frequency regulation resources become more diversified and the differences in frequency regulation response speed and depth become increasingly prominent, the accuracy of evaluating the frequency regulation capability of multiple types of equipment using existing methods will decrease. The aforementioned technical bottleneck lies in the lack of a quantitative correlation between the transient frequency security level of the large power grid and various frequency regulation resources, making it difficult to quickly assess the contribution of each frequency regulation resource during the transient frequency response process of the large power grid.

[0003] With the construction of new power systems, the proportion of renewable energy units in the system is constantly increasing, squeezing the operating space of conventional units. On the one hand, because a large number of renewable energy units cannot provide inertia support, the system exhibits multi-inertia center characteristics and significant frequency spatial distribution characteristics; on the other hand, renewable energy reserves are costly and have limited short-term available frequency regulation energy. With the continuous increase in the proportion of renewable energy and the profound changes in grid frequency characteristics, existing technologies cannot fully meet the needs of analyzing and evaluating the frequency characteristics of large power grids. There is a lack of quantitative correlation between the transient frequency security level of a large power grid and various frequency regulation resources, making it difficult to quickly assess the contribution of each frequency regulation resource during the transient frequency response of a large power grid. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this invention provides a method and apparatus for evaluating the capabilities of frequency modulation equipment based on an energy perspective.

[0005] According to one aspect of the present invention, a method for evaluating the capability of frequency modulation equipment based on an energy perspective is provided, comprising:

[0006] The system acquires equipment parameters, network parameters, and preset fault information, and initializes the equivalent speed governor model parameters by aggregating the equipment parameters.

[0007] Based on the equivalent speed governor model parameters, network parameters, and preset fault information, a vibration dynamics model of the power system is established.

[0008] The common-mode component and oscillation component of the frequency variation of the power system are obtained by solving the vibration dynamics model using modal analysis.

[0009] The frequency variation of each unit in the power system is determined based on the common-mode component and the oscillation component.

[0010] Based on the frequency change and the parameters of the equivalent governor model, the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation is calculated; where the total frequency regulation energy is the sum of rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

[0011] Based on the total frequency regulation energy of each unit, the contribution of each unit to the frequency regulation capability is calculated, and the evaluation conclusion of the frequency regulation capability of each unit is output according to the contribution of each unit to the frequency regulation capability.

[0012] Optionally, initializing the equivalent governor model parameters includes:

[0013]

[0014] In the formula, S N1 and S N2 The equipment capacities of the first and second speed governors are respectively; R1 and R2 are the droop coefficients of the first and second speed governors respectively; and F is respectively... H1 and F H2 The high-pressure cylinder work ratio of the first speed governor and the second speed governor, respectively, T R1 and T R2 The equivalent first-order inertial time constants of the first and second speed controllers are respectively, and the aggregated equipment and control parameters are as follows: S N R is the aggregated equipment capacity, R is the aggregated adjustment coefficient, and F is the aggregated droop coefficient. H The power ratio of the high-pressure cylinder after polymerization, T R The equivalent first-order inertial time constant after aggregation. The maximum value of the equivalent first-order inertial time constant of the governor of a conventional unit, λ1, λ2, M, k M a, b, c, and d are all intermediate variables.

[0015] Alternatively, the expression for the vibration dynamics model is:

[0016] T J s 2 δ(s)+Dsδ(s)+B K δ(s)=-B F ΔP L (s)+ΔP mSYN (s)

[0017] =-B F ΔP L (s)+G SYN (s)sδ(s)

[0018] In the formula, δ is the column vector composed of the internal potential phase angle, and T is the internal potential phase angle. JD represents a diagonal matrix composed of the inertial time constant and damping coefficient of each unit, and B represents the diagonal matrix composed of the inertial time constant and damping coefficient of each unit. K To represent the stiffness matrix of the interaction between the rotors of each unit, B F The distribution coefficient matrix for allocating system disturbance power to the rotors of each unit, ΔP mSYN G is a column vector consisting of the changes in the unit's mechanical power in the original coordinate system. SYN (s) is a diagonal matrix composed of the primary frequency regulation transfer functions of each governor of the synchronous generator under the original coordinates, ΔP L It is a column vector composed of the disturbance power of each node in the system.

[0019] Optionally, the common-mode components and oscillation components of the frequency variation of the power system are obtained by solving the vibration dynamics model using modal analysis, including:

[0020] By solving |T J Ω 2 -B K | = 0 yields the natural oscillation frequency Ω, where Ω = [Ω1, Ω2, ..., Ω n ];

[0021] Solve for each oscillation frequency eigenvectors Φ (i) The vibration dynamics model in modal coordinates can be obtained, where Φ (i) The matrix Φ, composed of column vectors, is called the mode shape matrix;

[0022] Solving the vibration dynamics model in modal coordinates yields the frequency change ω of mode 1. p1 And the frequency change ω of the i-th mode with a modal order greater than 1 pi The frequency change ω of mode 1 p1 For common-mode components, the frequency change ω of the i-th mode with a modal order greater than 1 pi This is an oscillating component.

[0023] Optionally, the expression for the vibration dynamics model in modal coordinates is:

[0024]

[0025] In the formula, M p Let D be the mass matrix in modal coordinates. p K is the damping matrix in modal coordinates. p Let F be the elasticity matrix in modal coordinates. p Let ΔP be the excitation vector in modal coordinates. mSYNp A column vector representing the changes in the unit's mechanical power in modal coordinates;

[0026] Mode 1 frequency change ω p1The time-domain expression is:

[0027]

[0028] in,

[0029]

[0030] In the formula, F p1 For F p The first element, M p1 For M p The first diagonal element, D p1 D p The first diagonal element, k D k G A m μ, η All are intermediate variables;

[0031] The frequency change ω of the i-th mode with a modal order greater than 1 pi The time-domain expression is:

[0032]

[0033] in,

[0034]

[0035] In the formula, t is time, F pi For F p The i-th element, M pi For M p The i-th diagonal element, D pi D p The i-th diagonal element, ξ' i Ω' bi All of these are intermediate variables.

[0036] Optionally, the frequency change ω of the i-th unit i The time-domain expression is:

[0037]

[0038] In the formula, n is the number of synchronous generator sets in the system, t is the time, and ω is the time. pj Φ represents the frequency change of the j-th mode. (i,j) is the element in the i-th row and j-th column of the mode shape matrix.

[0039] Optionally, based on the frequency change and the parameters of the equivalent governor model, the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation is calculated, including:

[0040] The change in mechanical power of each unit is obtained based on the change in frequency.

[0041] Based on the mechanical power change and the parameters of the equivalent governor model, the rotor kinetic energy release / absorption, damping energy and primary frequency regulation energy of each unit are calculated.

[0042] The total frequency regulation energy of each unit is obtained by summing the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

[0043] Optionally, the change in mechanical power ΔP of each unit mi The calculation expression is:

[0044]

[0045] Where, ΔP mi_max Let ΔP be the upper limit of the mechanical power of the i-th unit. mtempi (t) is an intermediate variable, and its calculation expression is:

[0046]

[0047] Among them, R i Let T be the droop coefficient of the i-th generating unit. Ri Let be the equivalent first-order inertial time constant of the i-th unit;

[0048] The kinetic energy release / absorption of each unit's rotor Rotori The calculation expression is:

[0049]

[0050] in,

[0051]

[0052] In the formula, ω i Let S be the frequency change of the i-th unit. Ni Let T be the rated capacity of the i-th unit. Ji Let t be the inertial time constant of the i-th unit. nadir μ, η All are intermediate variables;

[0053] Damping energy E of each unit Di :

[0054]

[0055] Calculate the primary frequency regulation energy of each unit, where the frequency regulation energy E of the i-th unit is... PFRi for:

[0056]

[0057] In the formula, ΔP mi Let be the change in mechanical power of the i-th unit.

[0058] Optionally, the frequency regulation capability contribution C of each unit i The calculation expression is:

[0059]

[0060] In the formula, E totali Let be the total frequency regulation capacity of the i-th unit; n is the number of synchronous generator units in the system.

[0061] According to another aspect of the present invention, an energy-based frequency modulation device capability assessment apparatus is provided, comprising:

[0062] The aggregation module is used to obtain the system's equipment parameters, network parameters, and preset fault information, and to initialize the equivalent speed governor model parameters by aggregating the equipment parameters.

[0063] A module is established to build a vibration dynamics model of the power system based on the parameters of the equivalent speed governor model, network parameters, and preset fault information.

[0064] The solver module is used to solve the vibration dynamics model according to the modal analysis method to obtain the common-mode components and oscillation components of the frequency variation of the power system.

[0065] The determination module is used to determine the frequency variation of each unit in the power system based on the common-mode component and the oscillation component.

[0066] The first calculation module is used to calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation based on the frequency change and the parameters of the equivalent governor model; wherein, the total frequency regulation energy of the unit is the sum of the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy;

[0067] The second calculation module is used to calculate the frequency regulation capability contribution of each unit based on the total frequency regulation energy of each unit, and output the frequency regulation capability evaluation conclusion of each unit according to the frequency regulation capability contribution of each unit.

[0068] According to another aspect of the present invention, a computer-readable storage medium is provided, the storage medium storing a computer program for performing the methods described in any of the above aspects of the present invention.

[0069] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising: a processor; a memory for storing executable instructions of the processor; the processor being configured to read the executable instructions from the memory and execute the instructions to implement the method described in any of the preceding aspects of the present invention.

[0070] Therefore, this invention proposes an energy-based method for evaluating the capability of frequency regulation equipment. Through six modules—acquiring new energy and system parameters, an equivalent governor model, calculating frequency common-mode components and oscillation components, calculating the mechanical power variation of each unit, calculating the total frequency regulation energy of each unit and its proportion, and outputting the frequency regulation capability analysis results of each unit—it achieves rapid evaluation of the capability of frequency regulation equipment from an energy perspective. This patent constructs an energy-based method for evaluating the capability of frequency regulation equipment such as synchronous generators and new energy generator sets under the spatiotemporal distribution characteristics and complex constraints of frequency regulation. Attached Figure Description

[0071] Exemplary embodiments of the present invention can be more fully understood by referring to the following figures:

[0072] Figure 1 This is a flowchart illustrating an energy-based method for evaluating the capabilities of frequency modulation equipment, provided in an exemplary embodiment of the present invention.

[0073] Figure 2 This is a schematic diagram of the New England 68-node system topology provided in an exemplary embodiment of the present invention;

[0074] Figure 3 This is a schematic diagram of the structure of a frequency modulation device capability assessment device based on an energy perspective provided in an exemplary embodiment of the present invention;

[0075] Figure 4 This is the structure of an electronic device provided in an exemplary embodiment of the present invention. Detailed Implementation

[0076] Hereinafter, exemplary embodiments according to the present invention will be described in detail with reference to the accompanying drawings. Obviously, the described embodiments are merely some embodiments of the present invention, and not all embodiments of the present invention. It should be understood that the present invention is not limited to the exemplary embodiments described herein.

[0077] It should be noted that, unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps described in these embodiments do not limit the scope of the invention.

[0078] Those skilled in the art will understand that the terms "first," "second," etc., in the embodiments of the present invention are only used to distinguish different steps, devices, or modules, and do not represent any specific technical meaning, nor do they indicate a necessary logical order between them.

[0079] It should also be understood that in the embodiments of the present invention, "multiple" can refer to two or more, and "at least one" can refer to one, two or more.

[0080] It should also be understood that any component, data or structure mentioned in the embodiments of the present invention can generally be understood as one or more unless explicitly defined or given contrary instructions in the context.

[0081] Furthermore, the term "and / or" in this invention is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this invention generally indicates that the preceding and following related objects have an "or" relationship.

[0082] It should also be understood that the description of the various embodiments in this invention emphasizes the differences between the various embodiments, and the similarities or similarities can be referred to each other. For the sake of brevity, they will not be described in detail.

[0083] At the same time, it should be understood that, for ease of description, the dimensions of the various parts shown in the accompanying drawings are not drawn according to actual scale.

[0084] The following description of at least one exemplary embodiment is merely illustrative and is in no way intended to limit the invention or its application or use.

[0085] Techniques, methods, and equipment known to those skilled in the art may not be discussed in detail, but where appropriate, they should be considered part of the specification.

[0086] It should be noted that similar labels and letters in the following figures indicate similar items; therefore, once an item is defined in one figure, it does not need to be discussed further in subsequent figures.

[0087] The embodiments of this invention can be applied to electronic devices such as terminal devices, computer systems, and servers, and can operate together with a wide range of other general-purpose or special-purpose computing system environments or configurations. Well-known examples of terminal devices, computing systems, environments, and / or configurations suitable for use with electronic devices such as terminal devices, computer systems, and servers include, but are not limited to: personal computer systems, server computer systems, thin clients, thick clients, handheld or laptop devices, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputer systems, mainframe computer systems, and distributed cloud computing environments including any of the above systems, etc.

[0088] Electronic devices such as terminal devices, computer systems, and servers can be described in the general context of computer system executable instructions (such as program modules) executed by a computer system. Typically, program modules can include routines, programs, object programs, components, logic, data structures, etc., which perform specific tasks or implement specific abstract data types. Computer systems / servers can be implemented in distributed cloud computing environments, where tasks are executed by remote processing devices linked through communication networks. In distributed cloud computing environments, program modules can reside on local or remote computing system storage media, including storage devices.

[0089] Exemplary method

[0090] Figure 1 This is a schematic flowchart of an energy-based capability assessment method for frequency modulation devices provided in an exemplary embodiment of the present invention. This embodiment can be applied to electronic devices, such as… Figure 1 As shown, the energy-based frequency modulation equipment capability assessment method 100 includes the following steps:

[0091] Step 101: Obtain the system's equipment parameters, network parameters, and preset fault information, and initialize the equivalent speed governor model parameters by aggregating the equipment parameters.

[0092] Step 102: Based on the equivalent speed governor model parameters, network parameters, and preset fault information, establish a vibration dynamics model of the power system;

[0093] Step 103: Solve the vibration dynamics model using modal analysis to obtain the common-mode component and oscillation component of the frequency variation of the power system;

[0094] Step 104: Determine the frequency variation of each unit in the power system based on the common-mode component and the oscillation component;

[0095] Step 105: Calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation based on the frequency change and the parameters of the equivalent governor model; whereby the total frequency regulation energy of the unit is the sum of the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

[0096] Step 106: Based on the total frequency regulation energy of each unit, calculate the frequency regulation capability contribution of each unit, and output the frequency regulation capability evaluation conclusion of each unit according to the frequency regulation capability contribution of each unit.

[0097] Specifically, this invention constructs an energy-based method for evaluating the capability of frequency regulation equipment, addressing the challenge of assessing the frequency regulation capability of synchronous generators, new energy generator sets, and other frequency regulation equipment under complex constraints and the characteristics of frequency spatiotemporal distribution. The method includes the following steps:

[0098] Acquire the system's equipment parameters, network parameters, and preset fault information; aggregate the equipment parameters to obtain the equivalent speed governor model parameters;

[0099] Based on the equivalent governor model parameters, network parameters, and preset fault information, a vibration dynamics model is established; the common-mode components and oscillation components of the system are obtained by solving the vibration dynamics model using modal analysis; and the frequency of the unit is determined based on the common-mode components and oscillation components.

[0100] The contribution of each unit is evaluated by calculating the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation and its proportion; where the frequency regulation energy is the sum of rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

[0101] Output the evaluation conclusions of the frequency regulation capability of each unit.

[0102] Furthermore, by aggregating the device parameters, the equivalent speed governor model parameters are obtained, specifically:

[0103]

[0104] The parameters of the two speed controllers are denoted as follows: the equipment capacities are S respectively. N1 and S N2 The adjustment coefficients are R1 and R2, and the power ratio of the high-pressure cylinder is F. H1 and F H2 The equivalent first-order inertial time constant is T. R1 and T R2 The aggregated equipment and control parameters are as follows: Equipment capacity S N The adjustment coefficient is R, and the power ratio of the high-pressure cylinder is F. H The equivalent first-order inertial time constant is T. R , It is the maximum value of the equivalent first-order inertial time constant of the conventional unit speed governor, λ1, λ2, M, k M a, b, c, and d are all intermediate variables.

[0105] Furthermore, based on the equivalent governor model parameters, network parameters, and preset fault information, a vibration dynamics model is established, including:

[0106] Based on the equivalent governor model parameters, network parameters, and preset fault information, the multi-machine system, under the original coordinates, has a frequency-dynamic vibration model, specifically as follows:

[0107] T J s 2 δ(s)+Dsδ(s)+B K δ(s)=-B F ΔP L (s)+ΔP mSYN (s)

[0108] =-B F ΔP L (s)+G SYN (s)sδ(s)

[0109] Among them, the column vector composed of the internal potential phase angle δ, T J D represents a diagonal matrix composed of the inertial time constant and damping coefficient of each unit, and B represents the diagonal matrix composed of the inertial time constant and damping coefficient of each unit. K To represent the stiffness matrix of the interaction between the rotors of each unit, B F The distribution coefficient matrix for allocating system disturbance power to the rotors of each unit, ΔP mSYN G is a column vector consisting of the changes in the unit's mechanical power in the original coordinate system. SYN (s) is a diagonal matrix composed of the primary frequency regulation transfer functions of each governor of the synchronous generator under the original coordinates, ΔP L It is a column vector composed of the disturbance power of each node in the system.

[0110] Furthermore, by solving the above vibration dynamics model using modal analysis, the common-mode components and oscillatory components of the system frequency variation are obtained, including:

[0111] By solving |T J Ω 2 -B K | = 0 yields the natural oscillation frequency Ω, where Ω = [Ω1, Ω2, ..., Ω n Solve for each eigenvectors Φ (i) , with Φ (i) The matrix Φ, composed of column vectors, is called the mode shape matrix. From this, the vibration dynamics model in modal coordinates can be obtained as follows:

[0112]

[0113] Where M p Let D be the mass matrix in modal coordinates. p K is the damping matrix in modal coordinates. p Let F be the elasticity matrix in modal coordinates. p Let ΔP be the excitation vector in modal coordinates. mSYNp A column vector representing the changes in the unit's mechanical power in modal coordinates;

[0114] Solving the above equation yields the frequency change ω of mode 1. p1 The time-domain expression is:

[0115]

[0116] in,

[0117]

[0118] F p1 For F p The first element, M p1 For M p The first diagonal element, D p1 D p The first diagonal element.

[0119] The frequency change ω of the i-th mode with a modal order greater than 1 pi The time-domain expression is:

[0120]

[0121] in,

[0122]

[0123] F pi For F p The i-th element, M pi For M p The i-th diagonal element, D pi D p The i-th diagonal element. Further, based on the common-mode component and the oscillation component, the frequency variation of the unit is determined, including:

[0124] Transform the modal coordinates back to the original coordinates to obtain the frequency change ω of the i-th unit. i The time-domain representation is:

[0125]

[0126] Where n is the number of synchronous generator sets in the system.

[0127] Furthermore, based on the parameters of the equivalent governor model and the frequency variation of the units, the mechanical power variation ΔP of each unit is calculated. mi The calculation expression is:

[0128]

[0129] Where, ΔP mi_max Let ΔP be the upper limit of the mechanical power of the i-th unit. mtempi (t) is an intermediate variable, and its calculation expression is:

[0130]

[0131] Among them, R i Let T be the droop coefficient of the i-th generating unit. Ri Let be the equivalent first-order inertial time constant of the i-th unit;

[0132] Furthermore, the sum of the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy for each unit is calculated, including:

[0133] Calculate the energy released / absorbed by the rotor kinetic energy of each unit, E. Rotori :

[0134]

[0135] in,

[0136]

[0137] Calculate the damping energy E of each unit. Di :

[0138]

[0139] Calculate the primary frequency regulation energy of each unit, where the frequency regulation energy E of the i-th unit is... PFRi for:

[0140]

[0141] Furthermore, the total frequency regulation energy E of each unit is calculated. totali :

[0142] E totali =E Rotori +E Di +E PFRi

[0143] Furthermore, the contribution C of the frequency regulation capability of each unit is calculated. i :

[0144]

[0145] In a specific embodiment of the present invention, based on Figure 2 The New England 68-node example shown was used for application verification. The inertia and governor parameters of the 16 units are shown in Table 1. The loads at buses 21, 28, and 49 were modified to HVDCs with a rated power of 1000MW to simulate faults such as HVDC commutation failure. The disturbance fault was a DC blocking fault at bus 21, resulting in a power deficit of 1000MW in the system.

[0146] Table 1 shows a comparison of the frequency regulation energy and its proportion of each unit calculated based on the method of this invention.

[0147] Table 1 Calculation results of main parameters and indicators of the unit

[0148]

[0149] Therefore, this invention proposes an energy-based method for evaluating the capability of frequency regulation equipment. Through six modules—acquiring new energy and system parameters, an equivalent governor model, calculating frequency common-mode components and oscillation components, calculating the mechanical power variation of each unit, calculating the total frequency regulation energy of each unit and its proportion, and outputting the frequency regulation capability analysis results of each unit—it achieves rapid evaluation of the capability of frequency regulation equipment from an energy perspective. This patent constructs an energy-based method for evaluating the capability of frequency regulation equipment such as synchronous generators and new energy generator sets under the spatiotemporal distribution characteristics and complex constraints of frequency regulation.

[0150] Exemplary apparatus

[0151] Figure 3 This is a schematic diagram of the structure of a frequency modulation device capability assessment device based on an energy perspective, provided in an exemplary embodiment of the present invention. Figure 3 As shown, the device 300 includes:

[0152] The aggregation module 310 is used to acquire the system's equipment parameters, network parameters, and preset fault information, and to initialize the equivalent speed governor model parameters by aggregating the equipment parameters.

[0153] Module 320 is established to build a vibration dynamics model of the power system based on the parameters of the equivalent speed governor model, network parameters, and preset fault information.

[0154] Solver module 330 is used to solve the vibration dynamics model according to the modal analysis method to obtain the common mode component and oscillation component of the frequency change of the power system;

[0155] The determination module 340 is used to determine the frequency variation of each unit in the power system based on the common-mode component and the oscillation component.

[0156] The first calculation module 350 is used to calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation based on the frequency change and the parameters of the equivalent governor model; wherein, the total frequency regulation energy of the unit is the sum of the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy;

[0157] The second calculation module 360 ​​is used to calculate the frequency regulation capability contribution of each unit based on the total frequency regulation energy of each unit, and output the frequency regulation capability evaluation conclusion of each unit according to the frequency regulation capability contribution of each unit.

[0158] Optionally, initializing the equivalent governor model parameters includes:

[0159]

[0160] In the formula, S N1 and S N2 The equipment capacities of the first and second speed governors are respectively; R1 and R2 are the droop coefficients of the first and second speed governors respectively; and F is respectively... H1 and F H2 The high-pressure cylinder work ratio of the first speed governor and the second speed governor, respectively, T R1 and T R2 The equivalent first-order inertial time constants of the first and second speed controllers are respectively, and the aggregated equipment and control parameters are as follows: S N R is the aggregated equipment capacity, R is the aggregated adjustment coefficient, and F is the aggregated droop coefficient. H The power ratio of the high-pressure cylinder after polymerization, T R The equivalent first-order inertial time constant after aggregation. The maximum value of the equivalent first-order inertial time constant of the governor of a conventional unit, λ1, λ2, M, k M a, b, c, and d are all intermediate variables.

[0161] Alternatively, the expression for the vibration dynamics model is:

[0162] T J s 2 δ(s)+Dsδ(s)+B K δ(s)=-B F ΔP L (s)+ΔP mSYN (s)

[0163] =-B F ΔP L (s)+G SYN (s)sδ(s)

[0164] In the formula, δ is the column vector composed of the internal potential phase angle, and T is the internal potential phase angle. J D represents a diagonal matrix composed of the inertial time constant and damping coefficient of each unit, and B represents the diagonal matrix composed of the inertial time constant and damping coefficient of each unit. K To represent the stiffness matrix of the interaction between the rotors of each unit, B F The distribution coefficient matrix for allocating system disturbance power to the rotors of each unit, ΔP mSYN G is a column vector consisting of the changes in the unit's mechanical power in the original coordinate system. SYN (s) is a diagonal matrix composed of the primary frequency regulation transfer functions of each governor of the synchronous generator under the original coordinates, ΔP L It is a column vector composed of the disturbance power of each node in the system.

[0165] Optionally, the common-mode components and oscillation components of the frequency variation of the power system are obtained by solving the vibration dynamics model using modal analysis, including:

[0166] By solving |T J Ω 2 -B K | = 0 yields the natural oscillation frequency Ω, where Ω = [Ω1, Ω2, ..., Ω n ];

[0167] Solve for each oscillation frequency eigenvectors Φ (i) The vibration dynamics model in modal coordinates can be obtained, where Φ (i) The matrix Φ, composed of column vectors, is called the mode shape matrix;

[0168] Solving the vibration dynamics model in modal coordinates yields the frequency change ω of mode 1. p1 And the frequency change ω of the i-th mode with a modal order greater than 1 pi The frequency change ω of mode 1 p1 For common-mode components, the frequency change ω of the i-th mode with a modal order greater than 1 pi This is an oscillating component.

[0169] Optionally, the expression for the vibration dynamics model in modal coordinates is:

[0170]

[0171] In the formula, M p Let D be the mass matrix in modal coordinates. p K is the damping matrix in modal coordinates. p Let F be the elasticity matrix in modal coordinates. p Let ΔP be the excitation vector in modal coordinates. mSYNpA column vector representing the changes in the unit's mechanical power in modal coordinates;

[0172] Mode 1 frequency change ω p1 The time-domain expression is:

[0173]

[0174] in,

[0175]

[0176] In the formula, F p1 For F p The first element, M p1 For M p The first diagonal element, D p1 D p The first diagonal element, k D k G A m μ, η All are intermediate variables;

[0177] The frequency change ω of the i-th mode with a modal order greater than 1 pi The time-domain expression is:

[0178]

[0179] in,

[0180]

[0181] In the formula, t is time, F pi For F p The i-th element, M pi For M p The i-th diagonal element, D pi D p The i-th diagonal element, ξ' i Ω' bi All of these are intermediate variables.

[0182] Optionally, the frequency change ω of the i-th unit i The time-domain expression is:

[0183]

[0184] In the formula, n is the number of synchronous generator sets in the system, t is the time, and ω is the time. pj Φ represents the frequency change of the j-th mode. (i,j) is the element in the i-th row and j-th column of the mode shape matrix.

[0185] Optionally, the first computing module 350 includes:

[0186] The change in mechanical power of each unit is obtained based on the change in frequency.

[0187] Based on the mechanical power change and the parameters of the equivalent governor model, the rotor kinetic energy release / absorption, damping energy and primary frequency regulation energy of each unit are calculated.

[0188] The total frequency regulation energy of each unit is obtained by summing the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

[0189] Optionally, the change in mechanical power ΔP of each unit mi The calculation expression is:

[0190]

[0191] Where, ΔP mi_max Let ΔP be the upper limit of the mechanical power of the i-th unit. mtempi (t) is an intermediate variable, and its calculation expression is:

[0192]

[0193] Among them, R i Let T be the droop coefficient of the i-th generating unit. Ri Let be the equivalent first-order inertial time constant of the i-th unit;

[0194] The kinetic energy release / absorption of each unit's rotor Rotori The calculation expression is:

[0195]

[0196] in,

[0197]

[0198] In the formula, ω i Let S be the frequency change of the i-th unit. Ni Let T be the rated capacity of the i-th unit. Ji Let t be the inertial time constant of the i-th unit. nadir μ, η All are intermediate variables;

[0199] Damping energy E of each unit Di :

[0200]

[0201] Calculate the primary frequency regulation energy of each unit, where the frequency regulation energy E of the i-th unit is...PFRi for:

[0202]

[0203] In the formula, ΔP mi Let be the change in mechanical power of the i-th unit.

[0204] Optionally, the frequency regulation capability contribution C of each unit i The calculation expression is:

[0205]

[0206] In the formula, E totali Let be the total frequency regulation capacity of the i-th unit; n is the number of synchronous generator units in the system.

[0207] Exemplary electronic device

[0208] Figure 4 This is the structure of an electronic device provided in an exemplary embodiment of the present invention. For example... Figure 4 As shown, the electronic device 40 includes one or more processors 41 and a memory 42.

[0209] The processor 41 may be a central processing unit (CPU) or other form of processing unit with data processing and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions.

[0210] The memory 42 may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, and the processor 41 may execute the program instructions to implement the methods of the software programs of the various embodiments of the present invention described above, and / or other desired functions. In one example, the electronic device may also include an input device 43 and an output device 44, these components being interconnected via a bus system and / or other forms of connection mechanisms (not shown).

[0211] In addition, the input device 43 may also include, for example, a keyboard, a mouse, etc.

[0212] The output device 44 can output various information to the outside. The output device 44 may include, for example, a display, a speaker, a printer, and a communication network and its connected remote output devices, etc.

[0213] Of course, for the sake of simplicity, Figure 4 Only some of the components of this electronic device relevant to the present invention are shown, omitting components such as buses, input / output interfaces, etc. In addition, the electronic device may include any other suitable components depending on the specific application.

[0214] Exemplary computer program product and computer readable storage medium

[0215] In addition to the methods and apparatus described above, embodiments of the present invention may also be computer program products, which include computer program instructions that, when executed by a processor, cause the processor to perform the steps in the methods according to various embodiments of the present invention described in the "Exemplary Methods" section above.

[0216] The computer program product can be written in any combination of one or more programming languages ​​to perform the operations of the embodiments of the present invention. The programming languages ​​include object-oriented programming languages ​​such as Java and C++, as well as conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on the user's computing device, partially on the user's computing device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0217] Furthermore, embodiments of the present invention may also be computer-readable storage media storing computer program instructions thereon, which, when executed by a processor, cause the processor to perform the steps of the methods according to various embodiments of the present invention described in the "Exemplary Methods" section above.

[0218] The computer-readable storage medium may be any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or any combination thereof. More specific examples (a non-exhaustive list) of readable storage media include: an electrical connection having one or more wires, a portable disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof.

[0219] The basic principles of the present invention have been described above with reference to specific embodiments. However, it should be noted that the advantages, benefits, and effects mentioned in the present invention are merely examples and not limitations, and should not be considered as essential features of each embodiment of the present invention. Furthermore, the specific details disclosed above are for illustrative and facilitative purposes only, and are not limitations. These details do not limit the present invention to the necessity of employing the aforementioned specific details.

[0220] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For system embodiments, since they largely correspond to method embodiments, the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments.

[0221] The block diagrams of devices, systems, devices, and systems involved in this invention are merely illustrative examples and are not intended to require or imply that they must be connected, arranged, or configured in the manner shown in the block diagrams. As those skilled in the art will recognize, these devices, systems, devices, and systems can be connected, arranged, and configured in any manner. Words such as “comprising,” “including,” “having,” etc., are open-ended terms meaning “including but not limited to,” and are used interchangeably with them. The terms “or” and “and” as used herein refer to the terms “and / or,” and are used interchangeably with them unless the context clearly indicates otherwise. The term “such as” as used herein refers to the phrase “such as but not limited to,” and is used interchangeably with it.

[0222] The methods and systems of the present invention may be implemented in many ways. For example, they may be implemented by software, hardware, firmware, or any combination of software, hardware, and firmware. The above-described order of steps for the methods is for illustrative purposes only, and the steps of the methods of the present invention are not limited to the order specifically described above unless otherwise specifically stated. Furthermore, in some embodiments, the present invention may also be implemented as a program recorded on a recording medium, the program comprising machine-readable instructions for implementing the methods according to the present invention. Thus, the present invention also covers recording media storing programs for performing the methods according to the present invention.

[0223] It should also be noted that in the systems, apparatus, and methods of the present invention, the components or steps can be disassembled and / or recombined. These disassemblies and / or recombinations should be considered equivalents of the present invention. The above description of the disclosed aspects is provided to enable any person skilled in the art to make or use the invention. Various modifications to these aspects will be readily apparent to those skilled in the art, and the general principles defined herein can be applied to other aspects without departing from the scope of the invention. Therefore, the invention is not intended to be limited to the aspects shown herein, but rather to be carried out within the widest scope consistent with the principles and novel features disclosed herein.

[0224] The above description has been given for purposes of illustration and description. Furthermore, this description is not intended to limit the embodiments of the invention to the forms disclosed herein. Although numerous exemplary aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations thereof.

Claims

1. A method for evaluating the capability of frequency modulation equipment based on an energy perspective, characterized in that, include: The system acquires equipment parameters, network parameters, and preset fault information, and initializes the equivalent speed governor model parameters by aggregating the equipment parameters. Based on the equivalent speed governor model parameters, the network parameters, and the preset fault information, a vibration dynamics model of the power system is established. The vibration dynamics model is solved using modal analysis to obtain the common-mode component and oscillation component of the frequency variation of the power system. The frequency variation of each unit in the power system is determined based on the common-mode component and the oscillation component. Based on the frequency change and the parameters of the equivalent governor model, calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation; wherein, the total frequency regulation energy of the unit is the sum of rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy; Based on the total frequency regulation energy of each unit, the frequency regulation capability contribution of each unit is calculated, and the frequency regulation capability evaluation conclusion of each unit is output according to the frequency regulation capability contribution of each unit.

2. The method according to claim 1, characterized in that, The initialization parameters of the equivalent speed governor model include: In the formula, S N1 and S N2 The equipment capacities of the first and second speed governors are respectively; R1 and R2 are the droop coefficients of the first and second speed governors respectively; and F is respectively... H1 and F H2 The high-pressure cylinder work ratio of the first speed governor and the second speed governor, respectively, T R1 and T R2 The equivalent first-order inertial time constants of the first and second speed controllers are respectively, and the aggregated equipment and control parameters are as follows: S N R is the aggregated equipment capacity, R is the aggregated adjustment coefficient, and F is the aggregated droop coefficient. H The power ratio of the high-pressure cylinder after polymerization, T R The equivalent first-order inertial time constant after aggregation. The maximum value of the equivalent first-order inertial time constant of the governor of a conventional unit, λ1, λ2, M, k M a, b, c, and d are all intermediate variables.

3. The method according to claim 1, characterized in that, The expression for the vibration dynamics model is: T J s 2 δ(s)+Dsδ(s)+B K δ(s)=-B F ΔP L (s)+ΔP mSYN (s) =-B F ΔP L (s)+G SYN (s)sδ(s) In the formula, δ is the column vector composed of the internal potential phase angle, and T is the internal potential phase angle. J D represents a diagonal matrix composed of the inertial time constant and damping coefficient of each unit, and B represents the diagonal matrix composed of the inertial time constant and damping coefficient of each unit. K To represent the stiffness matrix of the interaction between the rotors of each unit, B F The distribution coefficient matrix for allocating system disturbance power to the rotors of each unit, ΔP mSYN G is a column vector consisting of the changes in the unit's mechanical power in the original coordinate system. SYN (s) is a diagonal matrix composed of the primary frequency regulation transfer functions of each governor of the synchronous generator under the original coordinates, ΔP L It is a column vector composed of the disturbance power of each node in the system.

4. The method according to claim 1, characterized in that, The vibration dynamics model is solved using modal analysis to obtain the common-mode components and oscillation components of the power system frequency variation, including: By solving |T J Ω 2 -B K | = 0 yields the natural oscillation frequency Ω, where Ω = [Ω1, Ω2, ..., Ω n ]; Solve for each oscillation frequency eigenvectors Φ (i) The vibration dynamics model in modal coordinates can be obtained, where Φ (i) The matrix Φ, composed of column vectors, is called the mode shape matrix; Solving the vibration dynamics model in modal coordinates yields the frequency change ω of mode 1. p1 And the frequency change ω of the i-th mode with a modal order greater than 1 pi The frequency change ω of mode 1 p1 For the common-mode component, the frequency change ω of the i-th mode with a modal order greater than 1 pi This is an oscillating component.

5. The method according to claim 4, characterized in that, The expression for the vibration dynamics model in modal coordinates is: In the formula, M p Let D be the mass matrix in modal coordinates. p K is the damping matrix in modal coordinates. p Let F be the elasticity matrix in modal coordinates. p Let ΔP be the excitation vector in modal coordinates. mSYNp A column vector representing the changes in the unit's mechanical power in modal coordinates; The frequency change of mode 1 ω p1 The time-domain expression is: in, In the formula, F p1 For F p The first element, M p1 For M p The first diagonal element, D p1 D p The first diagonal element, k D k G A m μ, η All are intermediate variables; The frequency change ω of the i-th mode with a modal order greater than 1 pi The time-domain expression is: in, In the formula, t is time, F pi For F p The i-th element, M pi For M p The i-th diagonal element, D pi D p The i-th diagonal element, ξ' i Ω' bi All of these are intermediate variables.

6. The method according to claim 1, characterized in that, The frequency change ω of the i-th unit i The time-domain expression is: In the formula, n is the number of synchronous generator sets in the system, t is the time, and ω is the time. pj Φ represents the frequency change of the j-th mode. (i,j) is the element in the i-th row and j-th column of the mode shape matrix.

7. The method according to claim 1, characterized in that, Based on the frequency change and the parameters of the equivalent governor model, calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation, including: The mechanical power change of each unit is obtained based on the frequency change. Based on the mechanical power change and the equivalent governor model parameters, calculate the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy of each unit. The total frequency regulation energy of each unit is obtained by summing the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy.

8. The method according to claim 7, characterized in that, Mechanical power change ΔP of each unit mi The calculation expression is: Where, ΔP mi_max Let ΔP be the upper limit of the mechanical power of the i-th unit. mtempi (t) is an intermediate variable, and its calculation expression is: Among them, R i Let T be the droop coefficient of the i-th generating unit. Ri Let be the equivalent first-order inertial time constant of the i-th unit; The kinetic energy release / absorption of each unit's rotor Rotori The calculation expression is: in, In the formula, ω i Let S be the frequency change of the i-th unit. Ni Let T be the rated capacity of the i-th unit. Ji Let t be the inertial time constant of the i-th unit. nadir μ, η All are intermediate variables; Damping energy E of each unit Di : Calculate the primary frequency regulation energy of each unit, where the frequency regulation energy E of the i-th unit is... PFRi for: In the formula, ΔP mi Let be the change in mechanical power of the i-th unit.

9. The method according to claim 1, characterized in that, The contribution of each unit to the frequency regulation capability C i The calculation expression is: In the formula, E totali Let be the total frequency regulation capacity of the i-th unit; n is the number of synchronous generator units in the system.

10. A frequency modulation equipment capability assessment device based on an energy perspective, characterized in that, include: The aggregation module is used to acquire the system's equipment parameters, network parameters, and preset fault information, and to initialize the equivalent speed governor model parameters by aggregating the equipment parameters. A module is established to establish a vibration dynamics model of the power system based on the equivalent speed governor model parameters, the network parameters, and the preset fault information. The solution module is used to solve the vibration dynamics model according to the modal analysis method to obtain the common-mode component and oscillation component of the frequency change of the power system; The determination module is used to determine the frequency variation of each unit in the power system based on the common-mode component and the oscillation component. The first calculation module is used to calculate the total frequency regulation energy accumulated by each unit at the moment of maximum frequency deviation based on the frequency change and the parameters of the equivalent speed governor model; wherein, the total frequency regulation energy of the unit is the sum of the rotor kinetic energy release / absorption, damping energy, and primary frequency regulation energy; The second calculation module is used to calculate the frequency regulation capability contribution of each unit based on the total frequency regulation energy of each unit, and output the frequency regulation capability evaluation conclusion of each unit according to the frequency regulation capability contribution of each unit.

11. A computer-readable storage medium, characterized in that, The storage medium stores a computer program for performing the method described in any one of claims 1-11.

12. An electronic device, characterized in that, The electronic device includes: processor; Memory used to store the processor's executable instructions; The processor is configured to read the executable instructions from the memory and execute the instructions to implement the method described in any one of claims 1-11.

Images